The gut microbiota comprises thousands of bacterial species with a collective gene pool 150 times larger than the human genome with tremendous functional diversity (1). A gut microbiota, characteristic of Crohn’s disease (CD), has been repeatedly described (2, 3). Although much progress has been made in understanding interactions between the gut microbiota and intestinal immune system and their involvement CD etiology (4), it remains unclear whether changes in microbiota precede or follow CD onset. Thus far, all studies exploring the colonic microenvironment of CD patients have been of cross-sectional and hence unable to address reverse causality (2, 3). Assessing, during the course of disease modifying treatment, serial changes in the gut microbiota, particularly organisms associated with colonic inflammation, will offer better insights into the role of gut bacteria or their products in the etiology of CD. Such a study has not been performed before.
Exclusive enteral nutrition (EEN) is an established treatment in active pediatric CD (5). In contrast to biologics, the mode of EEN action remains unknown. We have recently shown that EEN treatment in children with active CD induced reduction in numbers of presumptively protective gut bacterial species (e.g., Faecalibacterium prausnitzii) and metabolites (e.g., butyrate) that were paradoxically associated with disease improvement (6). These intriguing observations in conjunction with the simple nutrient and food ingredients of EEN make it plausible to speculate that EEN works by modulating the gut microbiota. Assessing changes in the gut microbiota during a course of EEN can offer clues to which bacteria and/or their metabolites may be implicated in the etiology of CD.
Most previous studies have explored the gut microbiota in CD by targeting certain limited microbial species and bacterial metabolites for which primers/probes were available (6, 7). Recent advances in sequencing technologies allow a hypothesis-free assessment of the entire microbial community dynamics. Using massively parallel ultra-deep sequencing and state-of-the-art computational approaches, we have characterized microbial community structure and genetic functional capacity in fecal samples from children with active CD and healthy controls by both targeting the bacterial 16S rRNA genes and employing shotgun metagenomics. Using the same methodology, we studied, for the first time, the association of the pediatric CD microbiota with intestinal inflammatory markers and changes before, during, and after induction treatment with EEN.
Fecal samples were collected from 23 children (males: n=13; 6.9–14.7 years) with active CD (Montreal classification; Disease location: L2, n=3; L2+L4, n=4; L3, n=3; L3+L4, n=13; Disease behavior: B1, n=20; B2, n=2; B3, n=1) and for 15 of these (11 newly diagnosed), with at least two consecutive samples, during 8 weeks of EEN (Modulen, Nestle, Vevey, Switzerland) as described previously (6). A maximum of five serial samples were collected per patient. The first sample (A) was collected before or within 6 days of EEN initiation (89% collected within 4 days), two during EEN (B:˜16 & C:˜32 days), and one close to end of treatment (D:˜54 days). A final sample (E) was collected when patients returned to habitual diet (E:˜63 days after EEN). The number of participants with 1, 2, 3, 4, and 5 samples were 5, 3, 3, 2, and 10, respectively. Two fecal samples were collected, at least 2 months apart, from 21 healthy children (males: n=12; 4.6–16.9 years) with no known family history of inflammatory bowel disease as a control group. The healthy children were recruited from the local community via poster advertisement. No participant had received antibiotics for 3 months prior to the study. Concomitant treatment was prospectively recorded (Supplementary Table S1 online).
Serum albumin, CRP, erythrocyte sedimentation rate, and fecal calprotectin (8) were measured and Paediatric Crohn’s Disease Activity Index (PCDAI) (9) calculated at EEN initiation, end of EEN, and when patients returned to habitual diet. All these disease activity indices significantly decreased, and 62% of participants achieved clinical remission (PCDAI <10) at the end of EEN (Supplementary Table S1).
This study received ethics approval by the Yorkhill Research Ethics Committee (05/S0708/66). Carer and patients provided written consent.
Sequencing of 16S rRNA gene was performed to characterize microbial community composition. Detailed description of the 16S rRNA sequencing is available in Supplementary Methods. Briefly, bacterial DNA was isolated using the chaotropic method (6). 16S rRNA gene sequencing of the V4 region was performed on the MiSeq (Illumina, Essex, UK) platform using 2 × 250 bp paired-end reads. The V4 region was amplified using fusion Golay adapters barcoded on the reverse strand.
Shotgun metagenome sequencing
Shotgun metagenome samples, indicative of microbiota genetic functional capacity, were prepared for 69 of the samples (CD, n=56; Healthy, n=13) with the Nextera XT Prep Kit (Illumina, FC-131–1096, UK) and the Illumina dual-barcoding Nextera XT Index kit (Illumina, FC-131–1002, UK). Detailed description of the shotgun metagenome sequencing is available in Supplementary Methods. Clusters were generated on-board a HiSeq 2500 (Illumina) instrument and sequencing performed using TruSeq Rapid SBS Kit reagents (Illumina, FC-402–4001, FC-402–4002). Sequencing was performed following a paired-end 150 cycle recipe.
Detailed description of bioinformatics analysis is available as Supplementary File (Supplementary Methods). In brief, microbiota composition was characterized at multiple levels of resolution using genus, operational taxonomic unit (OTU), and ‘oligotyping’ assignments from the 16S rRNA sequencing, the latter two approximating the species level. Shotgun metagenomics reads were used for species taxonomic profiling using MetaPhlAn and assignments to functional modules, orthologs, and pathways through alignment to Kyoto Encyclopedia of Genes and Genomes (KEGG) using HUMAnN (10). The latter are indicative of the genetic functional capacity of the entire gut microbiota.
Statistical analysis was performed on log-transformed data, and Benjamini–Hochberg false-discovery rate was used to account for multiple comparisons. We used the Adonis function of the R vegan package to determine the overall amount and significance of differences between controls and CD before EEN. Each marker was tested independently using Kruskal–Wallis. For the 16S rRNA markers (genera, OTU, and oligotypes), only those with mean relative abundance >0.01% were tested. For the metagenome data (MetaPhlAn species and HUManN KEGG modules and orthologs), we used all markers as there were fewer of them. We used a similar procedure as for comparisons between controls and CD children to test for the changes between sample D (end of EEN) and E (habitual diet) for 16S rRNA and metagenome markers.
For CD subjects with two or more samples and who were still on or had just stopped EEN treatment (≤3 days of treatment end), individual marker data regressed against EEN days with a subject-specific intercept using the GLM function. Thus, any confounding effect of different sample collection times, particularly after EEN initiation and completion, was naturally accounted for within the analysis.
We correlated the log-transformed relative abundance of each marker against fecal calprotectin using Kendall’s rank correlation. For 16S rRNA markers, we tested only those with mean relative abundance >0.01% but for metagenome markers used all markers. All markers with P-value <0.01 were included in stepwise multivariate regression of calprotectin as a function of all these markers.
We characterized the microbial community using 16S rRNA gene sequencing in 117 samples (n=78 CD participants and n=39 healthy controls) and shotgun metagenomics in 65 (n=53 from CD participants and n=12 controls). Only samples with greater than 5,000 filtered and overlapped 16S rRNA or 100,000 shotgun metagenome reads were included.
Bacterial community structure (16SrRNAgene and shotgun metagenomics)
Controls vs. CD patients before EEN
The mean Shannon diversity, a measure of gut microbiota richness and evenness, was higher in controls than CD children prior to EEN (Controls vs. CD: 10.23 vs. 8.00; P=0.006; Figure 1). A higher diversity was observed using OTU than genus analysis. However, the reduction in diversity associated with CD was similar with both approaches (P=4.2 × 10−5, Supplementary Figure S1). In contrast, beta diversity, the variability of the OTU community structure between samples, was higher for CD children than healthy controls. This was apparent from Figure 1 but also using the betadisper function of the R vegan package. The dissimilarity of microbial community structure, among groups, expressed with the average Bray–Curtis distance for samples, was significantly higher in CD children (CD vs. Controls: 0.4913 vs. 0.4089; P=0.048). This effect was more significant using OTU analysis (CD vs. Controls: 0.5848 vs. 0.4868; P<0.001).
The microbial community structure was also different between controls and CD patients, with 9.9% of the variation in community structure explained by sample groups (P=0.001). An association between age and genus community structure was observed for controls (P=0.033) but not CD children. However, this effect was weak and swamped by the effect that grouping (Controls vs. CD) had on genus community structure. Similarly, no difference was observed (P-value=0.697) in genus community structure between samples of CD patients collected one day before EEN (n=7) compared with samples collected within 6 days of EEN (n=10). The effect of disease affecting the ileum on community diversity and structure was explored and compared with subjects without ileal disease, but the difference did not reach statistical significance (P=0.176).
Figure 2 displays the relative abundance of the top 20 genera present in both CD and healthy controls. In total, 36 genera out of 84 (43%) tested differed significantly between CD children and controls (Table 1). Genera belonging to the Lachnospiraceae were more abundant in controls, (e.g., Coprococcus, Pseudobutyrivibrio, and Ruminococcus) as were Subdoligranum and Faecalibacterium from Clostridiales. In contrast, other members of Clostridiales such as Peptostreptococcus were more abundant in CD children (Table 1). A number of genera were also more abundant in CD children from diverse phyla such as the Actinobacteria member Atopobium or the Proteobacteria member Escherichia–Shigella (Table 1). Similar observations resulted from OTU assignments (Figure 2). Overall, 9.6% of the variation (P<0.001) was explained by group (Controls vs. CD). The relative abundance of a total of 141 OTUs out of 351 (42%) tested differed significantly between controls and those with CD. The 20 most significant are presented in Supplementary Table S2. Many OTUs, more abundant in controls, were affiliated with Lachnospiraceae, (e.g., OTU88; OTU10 (Pseudobutyrivibrio); OTU18 (Ruminococcus)), whereas OTU69 (Streptococcus) was more abundant in samples from CD patients (Supplementary Table S2).
In-depth analysis of the Lachnospiraceae in 1,191,470 reads using oligotyping resulted in 93 high-resolution oligotypes, explaining even more of the variation associated with CD children compared with controls (12.4%; P<0.001). Forty-four oligotypes (47%) differed significantly between the two groups. The 20 most significant oligotypes are in Supplementary Table 3. The findings of this high-resolution analysis were consistent with those observed through the genera and OTU assignments. Oligotypes classified as E. rectale and Roseburia inulinivorans were more abundant in controls, whereas R. gnavus and Clostridium clostridioformes were associated with CD.
For classification of shotgun metagenome reads to species level with MetaPhlAn, a tendency for higher mean Shannon diversity was found for controls (CD vs. Controls: 18.49 vs. 14.30; P=0.230, Supplementary Figure S1). The two groups explained 15% of the variation in community structure (P<1.0 × 10−5), but at the whole community level no differences in species proportions were observed (R2=0.0733, P=0.476; Figure S2 online). In total, 36 species had a false discovery rate <10% (Supplementary Table S4). These results were consistent with other computational approaches using 16S rRNA sequencing. Species such as Bifidobacterium adolescentis, Ruminococcus bromii, Eubacterium spp., F. prausnitzii, Coprococcus eutactus, and Subdoligranulum variabile were more abundant in controls, whereas in CD children Streptococcus anginosus, Enterococcus faecalis, and R. gnavus were overrepresented along with other species (Supplementary Table S4).
Impact of EEN on CD
During EEN, the Shannon diversity, at the genus level, decreased (P=0.037; Figure 1). This became apparent even after 15 days on EEN (sample B) with a minimum diversity observed by ˜30 days (sample C), a slight recovery toward the end of EEN (sample D) and complete recovery to pre-treatment levels when patients returned to habitual diet (sample E) (Figure 1). When genus diversity was regressed against days on EEN, using subject-dependent linear regression, to overcome any confounding effects of the time of sample collection during EEN, we found (Supplementary Figure S2) that every 10 days on EEN 0.6 genus diversity equivalents are lost (coeff (SE): −0.06 (0.012); P=2.36 × 10−5). Similar results were observed for classification of shotgun metagenome reads to species with MetaPhlAn (P=0.0023) and a trend based on OTU analysis (P=0.140; Supplementary Figure S1).
At the community structure level, EEN made the CD microbial community even more dissimilar to that typical of healthy controls (Figure 1). Sample time explained 6.4% of the variation in the microbial community structure of CD children (P=0.094). At all points during EEN, there was no difference in microbial community structure between patients with and without ileal disease involvement. Changes in the top 20 taxa during EEN are displayed in Figure 2. Microbial community structure at each EEN sample time point was independent of other concomitant treatment or whether this was the first or a subsequent EEN course suggesting that the observed changes were associated with EEN only.
When relative abundance of each bacterial taxon, based on OTU and genus assignments, was regressed against days on EEN, allowing for a subject-dependent intercept, representing different starting abundances in different children, the vast majority of significant changes represented a reduction in relative abundance. Indeed, of 34 genera that significantly changed, only one, Lactococcus, increased with EEN (Table 2). Interestingly, some of the most negatively impacted genera were those whose abundance was already lower in CD children, before EEN initiation, than controls (e.g., Bifidobacterium, Ruminococcus, and Faecalibacterium). A similar impact was observed in OTU analysis (Supplementary Table S5), with a large number of OTUs decreasing during EEN. For example, OTU3, assigned to Bifidobacterium, decreased in some children by two logs of relative abundance (Figure 3), whereas OTU4, assigned to R. gnavus, increased in most subjects during EEN (Figure 3). Oligotyping of Lachnospiraceae confirmed these patterns, with all oligotypes significantly decreasing in abundance during EEN including oligotypes that were significantly higher or lower in CD (Supplementary Table S6).
Correlations with fecal calprotectin
We analyzed our data set by direct correlation to calprotectin from the same sample to allow proxy analysis of colonic inflammation. Twenty-six genera significantly correlated with calprotectin in CD children (Supplementary Table S7); both strong positive and negative associations were found. Performing stepwise regression analysis, minimizing the overall model Akaike Information Criterion, to determine which associations were most important accounting for correlations between genera, seven of 14 genera were selected (Supplementary Table S8). The relative abundance of three selected genera was negatively associated with calprotectin (Anaerostipes, Coprococcus, and unknown Erysipelotrichaceae). Four were positively associated (Bacteroides, Gemella, Haemophilus, and Atopobium). Overall, the final model explained a remarkable 69% of the variation in fecal calprotectin concentration (Supplementary Table S8). A higher number of 35 OTUs significantly correlated (Supplementary Table S9) with calprotectin, and 14 of these were selected by stepwise regression (Supplementary Table S10). The final model predicted 78% of the variation in calprotectin levels among CD children. Among those OTUs that were most significantly positively associated with fecal calprotectin were OTU122 (Atopobium), OTU1073 (Lachnospiraceae), OTU45 (Streptococcus), OTU1086 (Peptostreptococcus), and OTU41 (Bacteroides). There were a number of OTUs negatively associated with calprotectin including OTU3 (Bifidobacterium), OTU21 (Lachnospiraceae), and OTU106 (Anaerostipes). Using oligotyping, no significant correlations were found between oligotypes and calprotectin, after correcting for multiple comparisons.
Genetic functional capacity using shotgun metagenomics
Healthy children compared with CD patients before EEN
When reads were assigned to metabolic KEGG modules with HUMAnN, a significantly larger mean Shannon diversity, in richness equivalents, was observed in CD children compared with controls (CD vs. Controls: 48.1 vs. 43.2; P=0.05 Figure 4). Healthy and CD children differed at the level of KEGG modules, with marginal significance explaining 9.4% of the variation (P=0.059; Supplementary Figure S3). Several KEGG modules and orthologs (Supplementary Table S11) were different between the two groups, but their significance became marginal (P<0.100) when adjusted for multiple testing (Table 3). Modules more abundant in samples from CD patients included ubiquinone and lipopolysaccharide biosynthesis and the twin-arginine translocation (Tat) system, whereas key processes such as fatty acid biosynthesis, initiation, and sulfur reduction were overrepresented in controls (Table 3).
Impact of EEN on CD and correlations with fecal calprotectin
Four KEGG modules correlated with the number of days on EEN including pathways for biotin and thiamine biosynthesis (negative), the spermidine/putrescine transport system, and the shikimic acid pathway (positive; Figure 4 and Supplementary Table S12). Unlike relationships observed with taxa abundance, there were no significant correlations between calprotectin and relative abundance of KEGG modules.
Comparison between start, end of EEN, and return to habitual diet
In contrast to the impact of EEN on community structure and KEGG modules, there were no significant changes, adjusting for multiple comparisons, between sample A or end of EEN (D) and return to habitual diet (E) for all data types (genera, OTUs, oligotypes, shotgun metagenomics). Although a mean increase in relative abundance of several taxa was observed as CD children returned to habitual diet, these changes were not statistically significant. Those OTUs associated with fecal calprotectin (Supplementary Table S10) showed no significant change, but it is interesting to note that those negatively associated with calprotectin did increase in abundance, whereas those positively associated in general remained low at return to the habitual diet in unadjusted analysis (Supplementary Table S13).
In this study, we detailed a comprehensive sequential analysis of the entire fecal microbiota, its taxonomic profile, genetic functional capacity, and associations with fecal calprotectin, during a course of EEN in children with CD. Using a hypothesis-free approach, this study identified various bacterial taxa associated with CD and other that significantly changed during treatment with EEN. This study complements our previous publication of quantitative changes in concentration of selected bacterial metabolites and major bacterial groups previously implicated in the onset of CD.
Several significant differences in taxon community structure and relative abundance between CD children and controls were observed. The gut microbiota of CD children presented lower taxonomic diversity per individual but higher variation between individuals perhaps attributed to variations in disease behavior, location, or inflammation severity. Shotgun metagenome analysis revealed, for the first time, that despite lower taxonomic bacterial diversity a higher level of genetic functional diversity was associated with the CD microbiota prior to EEN treatment. This probably reflects a greater range of functional roles that can be exploited by the gut microbiota in inflamed gut, in active disease. This diverse microbial functional capacity tended to decrease during EEN to levels similar to healthy controls. These data suggest a higher degree of functional redundancy, with multiple species performing similar roles in healthy children, and may be important in maintaining gut health. In active CD, this functional redundancy may be lost and then restored during EEN.
Before treatment, we found bacterial taxa that were both more and less abundant in children with CD. Several of these findings concur with previous work (2, 3, 6), but here the methodology employed is not only hypothesis-free, unselective of certain taxa, but also more comprehensive, detailing bacterial taxonomy down to species, OTU, and oligotype levels. Numerous OTUs and species were less abundant in CD children, including classic commensals such as Faecalibacterium spp. and Bifidobacterium spp., but also a number of organisms that receive less attention, such as Eubacterium rectale and R. obeum. Conversely, some bacteria were more prevalent in CD, including E. coli/Shigella spp., Streptococcus spp., Peptostreptococcus spp., and Atopobium spp. More interestingly, several of these differences were not uniform or of the same direction within a genus, but varied between species and OTUs, which emphasizes the importance and need for an in-depth microbiota analysis when the role of the gut microbiota is explored. For example, R. gnavus was higher, but Ruminococcus obeum was lower in CD.
Although the results agree with previous research (2, 3, 6), possible reverse causality complicates interpretation of their role in the etiology of CD. Their importance and the relevance of these findings can be further unraveled only if these taxa are related to clinical outcomes, colonic inflammation, or changed during the course of disease modifying interventions. In the current study, correlations between fecal calprotectin, global fecal microbial composition, and genetic functional capacity were explored. Several taxa, which were different between CD and controls, did not associate with calprotectin and may have a less protagonist role in the etiology of CD. Instead, we found several other that did associate with calprotectin, negatively and positively. A subset of taxa was identified explaining a sizable amount of the variation in fecal calprotectin. Among them, OTU3 (Bifidobacterium spp.) had the strongest negative and OTU122 (Atopobium spp.) the strongest positive association with calprotectin in multivariate regression analysis.
Treatment with EEN induced a major reduction in the relative abundance of several species, some of which were already at lower abundance compared with controls. This is a paradox, as we would expect EEN treatment to normalize the perceived “dysbiotic” microbiota toward a healthier state. Instead, EEN treatment shifted the microbiota to a more “dysbiotic” state, with a greater NMDS distance from the healthy group compared with samples before EEN (Figure 1). Moreover, EEN caused reduction in the relative abundance of gut bacteria that were negatively associated with calprotectin, as well as those that correlated positively with it. Among these, OTU122 (Atopobium) was strongly negatively associated with the number of EEN days, but the commensal strain Bifidobacterium OTU3 was even more strongly negatively impacted. The particular Atopobium OTU identified, OTU122, has 99% identity in the V4 region with Atopobium parvulum, which is associated with halitosis (11), and has been associated with disease activity in pediatric CD (12) and hence is a good candidate for further research.
In agreement with our previous observations using qPCR (6), we observed a higher abundance of E. coli in CD. The E. coli population is of particular interest, as adherent invasive E. coli are overrepresented in CD patients (13). Moreover, the higher representation of KEGG modules coding for ubiquinone in CD children, before EEN, coincides with this observation of a higher relative abundance of Enterobacteriaceae. Ubiquinone is produced only by organisms capable of aerobic respiration, its main role is maintaining membrane stability by acting as a crucial membrane chain-terminating antioxidant. In nematodes (14), ubiquinone-deficient E. coli strains were more efficiently lysed in the pharynx or the intestine than their wild-type counterparts. Therefore, ubiquinone or bacterial respiration may act as virulence factors, allowing bacteria to colonize and harm the host. Overrepresentation of ubiquinone also fits with previous work (15). Likewise, other KEGG modules that were in higher abundance in CD, such as secretion and (Tat) systems, or module encoding for LPS synthesis, may be linked with overrepresentation of Enterobacteriaceae or other pathobionts in CD.
During EEN treatment, representation of KEGG modules relevant to the biosynthesis of biotin and thiamine, two B-complex vitamins, decreased. This may be indicative of the reduction in bacteria that bear genes encoding for these vitamins during EEN, such as bifidobacteria and E. coli spp., or changes in the synthesis of short or medium-chain fatty acids that requires these vitamins. Alternatively, supply of these vitamins through EEN may induce bacterial redundancy or indolence to produce them. In contrast, genes encoding for aromatic amino acid metabolism and transport of spermidine/putrescine, which have a major role in cell growth, and may indicate increased epithelial cell renewal and tissue healing (16), were overrepresented during EEN. Against expectations, no association between KEGG modules and calprotectin was found in CD. This probably reflects the lack of level in detail of currently available metabolic databases to resolve the key processes by which organisms impact the gut epithelium and that taxonomy serves as a better proxy to functional role.
The observational design of this study impedes us from suggesting a causative association between these findings and their role in CD pathogenesis or indeed in proposing a mechanism of therapeutic action of EEN based on the changes in the gut microbiota. Intervention studies are required with the aim of elucidating whether maintenance of such changes, following a successful course of EEN, can prolong disease remission and whether replication of such changes occur with other types of formula and dietary intervention. Once this is proven, a plausible mechanism of action for EEN, mediated by modulation of the gut microbiota, could be suggested. Conversely, these changes are likely to be direct consequences of a diet low in fiber, rich in protein, fat, simple sugars and micronutrients, modulating various microbial communities, bacterial pathways, and metabolite production. Whether the clinical efficacy of EEN could be enhanced by selectively tailoring its composition to stimulate the growth of beneficial microbes, suppressed during EEN, is also an intriguing concept to explore. Moreover, the very distinct CD microbiome we observed in this study may hold implications for improvement in diagnostic and prognostic classification of the disease, particularly in conjunction with other established disease biomarkers such as fecal calprotectin.
This study has some limitations. We performed analysis in feces, enabling a non-invasive collection of serial samples, but results may differ from mucosal adherent bacteria. Despite correcting for multiple testing, our participant’s sample size was modest, and this study may be underpowered to explore secondary outcomes such as associations between disease location and behavior with microbiota characteristics. Also, as a community expresses only a variable subset of its genome, at any time, metatranscriptomic and metabolomic data are required to complement and progress this work further. Some patients were on concomitant treatment, and there was a modest age difference between controls and CD children. However, any effect was minimal with the changes demonstrated driven predominantly by EEN or differences according to grouping (CD vs. Healthy) on community structure. This reconfirms evidence that diet is a strong determinant of the microbial microenvironment (17). Although PCDAI and CRP were highly correlated with fecal calprotectin (r=0.60; P=7.924 × 10−06 and r=0.42; P=4.644 × 10−03), we preferred to use and present data only for the latter, an observer-independent and reliable marker of colonic inflammation, as this study aimed specifically to associate a direct marker of colonic inflammation with microbiota composition and functional capacity measured in the same fecal sample. Moreover, preliminary analysis showed that fecal calprotectin was a better predictor of community structure explaining 10.5% of the variation in genus relative abundances across these samples compared with 9.3% and 3.5% for the PCDAI scores and CRP values, respectively.
Although there is still a degree of reverse causality to address even with the methodological approach applied here, this study successfully removed a great deal of “microbial signal noise”, potentially irrelevant to the microbial origins of CD, and thus focuses the field for more targeted research. Future interventions aiming to induce or maintain EEN-induced changes in the gut microbiota are needed to elucidate the significance of these observed changes in the etiology of CD and mechanism of action of EEN.
We thank participants and their families, the Department of Paediatric Gastroenterology, Hepatology, and Nutrition. The IBD team at Royal Hospital for Children, Glasgow, is supported by the Catherine McEwan Foundation and the Yorkhill IBD fund. The Open Access publication fee of this article was kindly sponsored by the Crohn’s in Childhood Research Association.