- Open Access
Intestinal dysbiosis in preterm infants preceding necrotizing enterocolitis: a systematic review and meta-analysis
Microbiomevolume 5, Article number: 31 (2017)
Necrotizing enterocolitis (NEC) is a catastrophic disease of preterm infants, and microbial dysbiosis has been implicated in its pathogenesis. Studies evaluating the microbiome in NEC and preterm infants lack power and have reported inconsistent results.
Methods and results
Our objectives were to perform a systematic review and meta-analyses of stool microbiome profiles in preterm infants to discern and describe microbial dysbiosis prior to the onset of NEC and to explore heterogeneity among studies. We searched MEDLINE, PubMed, CINAHL, and conference abstracts from the proceedings of Pediatric Academic Societies and reference lists of relevant identified articles in April 2016. Studies comparing the intestinal microbiome in preterm infants who developed NEC to those of controls, using culture-independent molecular techniques and reported α and β-diversity metrics, and microbial profiles were included. In addition, 16S ribosomal ribonucleic acid (rRNA) sequence data with clinical meta-data were requested from the authors of included studies or searched in public data repositories. We reprocessed the 16S rRNA sequence data through a uniform analysis pipeline, which were then synthesized by meta-analysis.
We included 14 studies in this review, and data from eight studies were available for quantitative synthesis (106 NEC cases, 278 controls, 2944 samples). The age of NEC onset was at a mean ± SD of 30.1 ± 2.4 weeks post-conception (n = 61). Fecal microbiome from preterm infants with NEC had increased relative abundances of Proteobacteria and decreased relative abundances of Firmicutes and Bacteroidetes prior to NEC onset. Alpha- or beta-diversity indices in preterm infants with NEC were not consistently different from controls, but we found differences in taxonomic profiles related to antibiotic exposure, formula feeding, and mode of delivery. Exploring heterogeneity revealed differences in microbial profiles by study and the target region of the 16S rRNA gene (V1-V3 or V3-V5).
Microbial dysbiosis preceding NEC in preterm infants is characterized by increased relative abundances of Proteobacteria and decreased relative abundances of Firmicutes and Bacteroidetes. Microbiome optimization may provide a novel strategy for preventing NEC.
Necrotizing enterocolitis (NEC) is a catastrophic disease that is a major cause of mortality in preterm infants who survive the first few days after birth . NEC occurs in 7% of infants born at less than 1500 g and up to 5% of admissions to the neonatal intensive care unit [2–5]. NEC is associated with a high mortality (15–30%) and long-term neurodevelopmental morbidity [2, 6].
The pathogenesis of NEC is not clear, and a unifying concept is lacking but microbial dysbiosis, formula feeding, and excessive inflammation have all been implicated [2, 7–9]. Compared to term infants, the intestinal microbiota of preterm infants has fewer bacterial species, less diversity, and increased proportions of potential pathogens [10, 11]. The microbial dysbiosis hypothesis of NEC is supported by the fact that NEC cannot be produced in germ free animals [12, 13] and by an association between early antibiotic use and NEC [14, 15]. Immune dysregulation in association with microbial dysbiosis including excessive toll-like receptor 4 (TLR4) signaling in response to lipopolysaccharide (LPS) [12, 16–18] and an exaggerated inflammatory response [8, 9] have been reported, but these are largely animal data.
Technological advances and availability of new molecular and analytic techniques such as those used in Human Microbiome Project [19, 20] have provided greater resolution in the evaluation of the neonatal intestinal microbiome. Most studies on the human microbiome do not have sufficient power to detect clinically important differences . Previously reported microbiome studies on NEC in preterm infants may have been underpowered to detect differences and have reported inconsistent results. We systematically reviewed studies that reported the intestinal microbiome in preterm infants who developed NEC in relation to controls and performed a meta-analysis of 16S ribosomal ribonucleic acid (rRNA) gene sequence data from eight of the included studies. To our knowledge, this is the first systematic review and meta-analysis of microbiome studies on NEC in preterm infants.
The primary objective of this review is to determine and describe intestinal microbial dysbiosis patterns in preterm infants preceding NEC (defined as stage II or stage III of Bell’s classification ). The secondary objective is to explore heterogeneity among studies that might explain inconsistency in the reported results.
We performed our systematic review according to the recommended “Preferred Reporting Items for Systematic Reviews and Meta-analyses” (PRISMA) guidelines  and “Meta-analysis of Observational Studies in Epidemiology” (MOOSE) consensus statement . We followed previously published methods for the meta-analyses of microbiome data .
Prospective or retrospective, case-control or cohort studies were included if they evaluated the neonatal intestinal microbiome in preterm infants with NEC compared to those infants without NEC, using culture-independent molecular techniques and reported α and β-diversity metrics and microbial profiles.
We searched for eligible studies using the Cochrane Neonatal Review Group’s (CNRG) search strategy (http://neonatal.cochrane.org/) without language restriction in April 2016. Search strategy and databases searched are outlined in Additional file 1.
Data collection and analyses
All titles and abstracts identified by our search strategy were screened for relevance by MP, and those deemed relevant were retrieved in full and evaluated for inclusion eligibility by MP and JN independently. All results were compared, and disagreements were resolved by mutual discussion. Relevant data were extracted from included studies and additional information to clarify the study design, and data was sought from the authors via email, for at least three attempts. The data extracted by the author were discussed for any discrepancies by input from a second author JN, and any conflicts were resolved by mutual discussion. In addition, 16S rRNA gene sequence data with clinical meta-data were requested from the principal investigators of the published studies. The methodological quality of each study was assessed using relevant items of the checklist proposed for observational studies by Viswanathan et al. .
Meta-analysis of microbiome data
We compiled all sequence files and subject phenotype data received from the investigators or downloaded from the National Center for Biotechnology Information Short Read Archive (NCBI SRA), Database of Genotypes and Phenotypes (NCBI dbGaP), and/or the European Nucleotide Archive (ENA) (accessed 06/27/2015). All sequences were processed using Quantitative Insights into Microbial Ecology package (QIIME 1.8.0) [25, 27]. Quality parameters of ≥25 average quality score, 200–1000 base length, no mismatches to barcode and primer, primer removal, no ambiguous bases, and maximum homopolymer length of 8 were used. As there were several 16S variable regions targeted, operational taxonomic units (OTUs) were assigned by closed-reference picking against the GreenGenes v13_8 database with the UCLUST algorithm in QIIME [28, 29]. We removed all samples that failed to yield at least 1200 sequences, resulting in the exclusion of one complete study  due to shallow sequencing depth. All analyses were performed on the rarefied data.
We derived the corrected gestational age (CGA) at diagnosis of NEC by adding gestational age at birth to the day of diagnosis of NEC for 61 individual cases of NEC. The distribution of NEC diagnoses in CGA (weeks) was tested for a continuous normal distribution with the Shapiro-Wilk method.
QIIME was used to calculate alpha diversity metrics, including observed OTUs, and the Shannon and Simpson diversity indices. Differences in alpha diversity metrics were tested with Mann-Whitney or Kruskal-Wallis tests, depending on the number of groups being compared, and Dunn’s multiple comparison testing was implemented for those with significant differences. We also utilized a negative binomial regression model to evaluate the degree to which number of species (our dependent variable) could be explained by gestational age, group (NEC vs. control), and a gestational age by group interaction term.
The weighted and unweighted UniFrac metrics were used for beta diversity comparisons. Differences in taxonomic relative abundance were evaluated with Mann-Whitney or Kruskal-Wallis tests and Benjamini-Hochberg false discovery rate correction.
Our search identified 6812 records. After removing duplicate reports, we screened 4217 records for inclusion and excluded 4193 records (PRISMA flow diagram, Fig. 1). Of the 24 full text articles reviewed for eligibility, 14 studies met our inclusion criteria. Ten excluded studies and reasons for exclusion are summarized in Table 1. Methodological quality of the 14 included studies was assessed using the item checklist and responses as “yes”, “no,” or “unclear” are reported in Table 2. The microbiota characteristics reported by the 14 included studies including alpha and beta diversity indices and microbial profiles are summarized in Additional file 2. It should be noted that when evaluating diversity in a non-complex microbial population, shifts may be attributed to one or few bacterial taxa. As such, differences in diversity may serve as a proxy for overrepresentation or underrepresentation of small number of organisms and are discussed later in the paper.
Results of the meta-analyses
Nine of 14 studies included in the systematic review provided sequence and clinical metadata for both the NEC cases and control patients. Data from one study  was excluded due to shallow sequencing depth, and hence, data from eight studies were incorporated in our quantitative synthesis (106 NEC cases, 278 controls, 2944 samples) [31–38] (Additional file 3). The results of the systematic review and meta-analysis are reported as recommened by PRISMA guidelines (Additional file 4).
Microbiome profile differences between NEC and controls
Corrected gestational age CGA at NEC diagnosis in preterm infants forms a normal distribution (n = 61, a subset of the population, Shapiro-Wilk test, p = 0.26) with a mean ± SD of 30.1 ± 2.4 weeks (median, 30 weeks) post-conception (Fig. 2). We analyzed data from samples collected at all CGAs and compared alpha-diversity metrics between infants with NEC and controls (Additional file 5). None of the alpha-diversity metrics calculated differed significantly between infants with NEC and controls (Fig. 3a–c). When alpha-diversity metrics (means of OTU richness and SDI) were plotted over CGA, we found that in control infants, OTU richness and SDI decreased from birth to their lowest level at about a CGA of 27 weeks and gradually increased to 36 weeks CGA (Additional file 6: Fig. S1a, b). The results of our negative binomial regression modeling support this pattern, indicating that the number of species increased significantly with gestational age (p < 0.001) and that NEC patients tended to have fewer species than controls after controlling for gestational age (p = 0.053). Principal coordinates analysis (PCoA) with the UniFrac metrics did not show clustering of samples as a function of NEC vs. control status (Fig. 3d, e).
Taxonomic abundances grouped by CGA showed consistent trends toward increased relative abundances of Proteobacteria and decreased relative abundances of Firmicutes and Bacteroidetes in infants with NEC, which contrasted with decreased relative abundances of Proteobacteria and increased abundances of Firmicutes observed in control infants (Fig. 4). Of note, this deviation in taxonomic abundances in preterm infants with NEC vs. controls happens around 30 weeks CGA (Fig. 4a), which was the mean and median age of onset of NEC in preterm infants in this pooled cohort. When examined at all CGAs, significant differences were observed with respect to the relative abundances of Proteobacteria, Firmicutes, and Bacteroidetes (Fig. 4b–d). Differences in the relative abundance of bacterial phyla and genera in NEC infants vs. controls are depicted in Fig. 4e, f, respectively.
Comparison of microbiota profiles between infants who had NEC compared to controls in relation to antibiotic exposure, diet, and mode of delivery
Investigating antibiotic exposure, which was defined inconsistently among the different studies (samples analyzed; controls with antibiotics n = 1822, controls without antibiotics, n = 107, NEC infants with antibiotics n = 273, and NEC infants without antibiotics, n = 29), OTU richness and SDI differed significantly between control infants who did not receive antibiotics compared to NEC infants who received antibiotics. OTU richness was significantly different between control infants with and without antibiotics, as was SDI between control infants with antibiotics and NEC infants with antibiotics (Additional file 7: Fig. S2a, b). PCoA by antibiotic exposure demonstrated clustering by the weighted and unweighted UniFrac metrics (Additional file 7: Fig. S2c, d). In conjunction with antibiotic exposure, we observed in control and NEC infants, increased relative abundances of Proteobacteria, and decreased relative abundances of Firmicutes, Actinobacteria, and Bacteroidetes. Among control infants, antibiotic exposure was associated with increased relative abundances of Klebsiella, unclassified Enterobacteriaceae, Proteus, Paenibacillus, Epulopiscium, and Pseudomonas at the genus level. Control infants not exposed to antibiotics had increased abundances of the genus Clostridium and unclassified Clostridiaceae (Additional file 7: Fig. S2e, f).
Investigating the role of diet (i.e., breast milk, formula, or both) (samples analyzed, breast milk-fed controls n = 967, formula-fed controls n = 55, both breast milk- and formula-fed controls n = 1300; breast milk-fed NEC n = 251, formula-fed NEC n = 15, both breast milk- and formula-fed NEC n = 241), significant differences in OTU richness but not SDI were observed between control breast milk-fed infants vs. NEC breast milk-fed infants and control subjects receiving breast milk vs. both breast milk and formula (Additional file 8: Fig. 3a, b). Meaningful clustering as a function of diet was not observed on PCoA plots (Additional file 8: Fig. S3c, d). Among control infants, we observed increased relative abundances of the phylum Firmicutes and decreased abundances of Proteobacteria in formula-fed infants compared to breast milk fed-infants (Additional file 8: Fig. S3e, f). Formula-fed babies who developed NEC had more Proteobacteria and less Firmicutes compared to breast milk-fed controls.
Investigating the role of mode of delivery (samples analyzed, vaginal delivery; control n = 683, NEC n = 178; C-section; control n = 1618, NEC n = 350), we found no differences in OTU richness or SDI betweeen NEC and controls but observed significant differences between control infants born by vaginal delivery compared to cesarean-section delivery (Additional file 9: Fig. S4a, b). No clustering was observed by mode of delivery on PCoA plots (Additional file 9: Fig. S4c, d). In control infants, we found increased relative abundance of the phylum Firmicutes in infants born by cesarean section and increased abundance of Bacteroidetes in infants born via vaginal delivery (Additional file 9: Fig. S4e, d).
Heterogeneity among included studies
Three studies in our meta-analysis targeted the V1-V3 region of the 16S rRNA gene, and five studies targeted the V3-V5 region. There were no differences in alpha-diversity metrics (OTU richness or SDI), nor was meaningful clustering observed within PCoA plots of UniFrac distances, as a function of variable region (Fig. 5a–d). In addition, studies targeting the V3-V5 region compared to the V1-V3 region reported increased relative abundances of Proteobacteria and decreased relative abundances of Firmicutes (Fig. 5e).
We did not find significant differences in alpha diversity among studies although some variations were noted (Fig. 6a, b), and some clustering as a function of study was observed on PCoA plots of UniFrac distances (Fig. 6c, d). NEC samples from studies by Normann et al. (targeting V3-V5) and Torraza et al. (targeting V1-V3) had significantly increased relative abundance of the phylum Firmicutes and significantly decreased relative abundance of the phylum Proteobacteria compared to other studies (Fig. 6e).
Our review identified 14 eligible studies of which data from eight studies were synthesized by meta-analysis. We assessed the methodological quality and risk of bias in the 14 included studies by the key components of study design, which has been deemed to be more important than assigning quality scores in observational studies . The methodological quality of the included observational studies was adequate . All studies had appropriate comparison groups, reported all our pre-specified outcomes, and results were plausible given the study limitations. Few studies evaluated outcomes in a subset of participants. Eleven out of 14 included studies addressed the issue of confounding by matching or stratification to avoid known confounders.
Alpha-diversity metrics reported by the included studies were inconsistent. Three studies reported a decrease in alpha diversity in preterm infants with NEC compared to controls [11, 32, 37], one study reported decreasing trends  and others did not report differences [31, 34, 36]. One study  included mixed models predicting Shannon diversity and controlling for gestational age at birth, route of delivery, and birth weight. Significant time-associated trends were identified among control subjects, which increased with respect to diversity over time, while NEC cases did not. Significant differences in beta-diversity metrics between NEC cases, and controls were reported by six studies, two of which were fingerprinting-based [39, 40] and four were based on sequencing small pools of PCR amplification products [11, 31, 33, 36].
Reported microbial profiles preceding NEC in preterm infants were variable across the 14 included studies. Eight of the 14 studies report an increase in the relative abundance of the phylum Proteobacteria (class Gammaproteobacteria or family Enterobacteriaceae) in the stools of infants who developed NEC [11, 30, 31, 33–37]. Various studies reported an increased relative abundance of Gammaproteobacteria and a decrease in other bacterial species , increased relative abundances of Enterococcus spp. and Citrobacter-like sequences , increased relative abundances of Proteobacteria, or a decreased relative abundance of Firmicutes in stools before NEC diagnosis . Increased relative abundances of Proteobacteria two weeks before and Actinobacteria 1 week before NEC diagnosis as well as lower counts of Bifidobacteria and Bacteriodetes have also been observed . Morrow et al. reported two types of intestinal dysbiosis associated with NEC ; one was dominated by Firmicutes and the other by Proteobacteria, specifically Enterobacteriaceae. All NEC cases in the fore-mentioned study lacked Propionibacterium. Some studies reported the abundance of Clostridium perfringens [35, 37, 39] in the stools of neonates with NEC, but one study reported the absence of Clostridia in association with NEC . Three studies reported unique OTUs in preterm infants with NEC: a unique bacterial OTU belonging to the Enterobacteriaceae family  and in the first stool samples, the presence of a novel sequence closest to Klebsiella pneumoniae . In our pooled meta-analysis, we found that decreased abundances of Firmicutes and Bacteroidetes and increased abundances of Proteobacteria precede the diagnosis of NEC in preterm infants.
We synthesized high-quality data available from eight studies in a sequence-based meta-analysis. The mean CGA at NEC diagnosis (n = 61) was 30.1 ± 2.4 weeks (median of 30 weeks). The Canadian network study suggested a bimodal distribution for the onset of NEC in preterm infants with an early onset and a late onset NEC . NEC occurred later after birth in infants born at younger gestational ages and earlier after birth in infants born at later gestational ages. Other investigators report a distribution for the onset of NEC in preterm infants similar to our findings [38, 42, 43]. We observed patterns of increased OTU richness over time where the number of species increased significantly with gestational age. This may represent the transition beyond the perinatally acquired microbiome which may be influenced by mode of delivery, the NICU environment, feeding, antibiotic exposure, or other factors and warrants further study.
A patterned progression of the intestinal microbiome has been observed in preterm infants . In a study of 922 16S rRNA sequence libraries from 58 infants, a patterned progression of bacterial classes from Bacilli to Gammaproteobacteria and then to Clostridia was described, and by CGA 33–36 weeks, the communities were well populated by anaerobes . We were unable to demonstrate this patterned progression in our meta-analysis due to lack of associated meta-data in some studies. In our meta-analyses, we found gradual shifts in the relative abundances of multiple phyla at about 27 weeks CGA, where decreased abundances of Firmicutes and Bacteroidetes and increased abundances of Proteobacteria precede the diagnosis of NEC in preterm infants. Proteobacteria are recognized by the innate immune system by TLR4 in the intestine, which may play a significant role in intestinal inflammation, enterocyte injury that may lead to the development of NEC in preterm infants . TLR4 expression is increased in the intestinal tract of preterm neonates and may regulate the balance between intestinal repair and injury . It is interesting to note that TLR4 mutant mice and TLR4 knock out mice are protected from NEC [17, 45–47], and inhibition of the TLR4 pathway may provide novel strategies in the treatment or prevention of NEC.
Our sequence-based analysis portrayed inconsistent differences in SDI, OTU richness, and Simpson diversity index between NEC cases and controls. Comparison of species richness over time found that when controlling for corrected gestational, NEC patients tended to have fewer species than controls (p = 0.053, negative binomial regression model). The absence of consistent differences with respect to alpha diversity may be due to the low starting point of microbial diversity that is observed in preterm infants, relatively small sample sizes in each of the included studies, or variability with respect to methodology for generating microbiome data. Of particular concern are the limited numbers of taxa present in preterm stool. This limitation places constraints on interpreting diversity changes as diversity in a non-complex population could reflect changes in only one taxon.
We did not observe distinct clustering of NEC and control samples by unweighted or weighted UniFrac metrics. The lack of consistency in distinct clustering of cases and controls may be due to the methodological, clinical, or study heterogeneity of the included studies. We found a consistent trend toward increased relative abundances of Proteobacteria and decreased in the relative abundances of Firmicutes and Bacteroidetes preceding NEC. This is consistent with the predominance of Proteobacteria in NEC infants reported in some studies [11, 31, 33, 36]. Although the pooling of data across studies increased sample size to give statistical power to detect differences, this benefit may have been outweighed by variation with respect to sampling and lab-based methodology (e.g., sample handling, storage, or DNA extraction method) associated with the different studies. Although no single genus or bacterium appears to cause or precede NEC, increased relative abundances of Proteobacteria appear to be consistent with NEC-associated dysbiosis. A recent study that studied fecal microbiome in preterm and term infants by shotgun metagenomic sequencing has reported association with uropathogenic Escherichia coli colonization and necrotizing enterocolitis , though this study was unusual in the predominance of E. coli among NEC cases. The clinical significance of the differences in phyla and genera between cases and controls of the pooled data clearly needs further scrutiny. If dysbiosis with increased Proteobacteria and decreased Firmicutes and Bacteroidetes is associated with NEC in preterm infants, then measures to balance these taxa by probiotic therapy or luminal antibiotics such as aminoglycosides may be beneficial.
We report our bacterial data at the taxonomic depth of phylum. However, we cannot state with certainty if phylum is the appropriate level to conduct these analyses. If the inciting mechanism is a shared phenotype across all genera within the class, then attempts to further refine an etiology (or protective group of bacteria) to a specific genus, species, clade, or genotype would obscure associations. Alternatively, identification of important strain level differences have been reported, which analyses focusing on higher taxonomic levels would overlook [33, 49]. Indeed, Warner et al. have been able to identify risk-associated populations at the class level (Gammaproteobacteria) and protection-associated populations at the genus level (the members of the putatively protective Negativicutes class were overwhelmingly Veillonella). We look forward to continuing analysis of these data sets and, in additional cohorts, in future attempts to confirm or refute the associations we identify and to identify injurious or protective microbes more narrowly, if possible.
We examined clinical determinants that are known to affect the developing intestinal microbiome, namely, antibiotic exposure preceding NEC, mode of delivery, and diet. We observed differences in the microbial profiles in the preterm infants who had antibiotics (any antibiotic) prior to the onset of NEC and found increased relative abundances of Proteobacteria and decreased relative abundances of Firmicutes and Actinobacteria. This may coincide with the proteobacterial bloom preceding NEC [11, 31, 33, 36]. Our observations support the results of studies that report association of prior antibiotic use with increased incidence of NEC [14, 15, 50], though we were not able to take into account various classes of antibiotics, which have been reported to be determinants of intra-gut bacterial community structure . We also sought taxonomic differences by mode of delivery (cesarean section or vaginal delivery) and found no differences in NEC infants but observed increased relative abundances of Firmicutes after cesarean section and increased relative abundances of Bacteroides after vaginal delivery in control infants. We also assessed taxonomic profiles of fecal microbiota based on type of feeding (breast milk, formula, or both). We found no differences in NEC infants but observed an increased relative abundance of Firmicutes in infant-fed formula compared to breast milk and an increased relative abundance of Proteobacteria in those fed breast milk compared to formula (controls). We identified no significant associations of mode of delivery or type of feeding with the microbial profiles of preterm infants who developed NEC.
We also investigated for potential sources of heterogeneity in the studies with respect to the variable region of the 16S rRNA gene targeted (V1-V3 vs. V3-V5) and by study and did observe differences in alpha-diversity metrics and UniFrac distances with respect to study and 16S rRNA target regions. Proteobacteria were more abundant and Firmicutes less abundant in studies targeting V3-V5 compared to V1-V3. NEC samples from studies by Normann et al. (targeting V3-V5) and Torraza et al. (targeting V1-V3) had significantly greater relative abundances of the phylum Firmicutes and significantly decreased relative abundances of the phylum Proteobacteria compared to other studies. It is widely recognized that the use of partial 16S rRNA gene sequences instead of the whole gene, or whole genome, can give rise to inconsistencies in evaluating microbiota composition. Primer bias may influence which members of the microbial community are amplified and detected, and OTU counts arising from different 16S rRNA gene regions can be inconsistent, challenging interpretation and, to some degree, limiting consensus with respect to which region best reflects the gut microbial community. Indeed, investigators have reported somewhat differing composition and varying phylogenetic divergence [52, 53]when multiple hypervariable regions were sequenced on the same samples simultaneously, leading to differential detection of members of Prevotella, Fusobacterium, Eubacterium, Enterococcus, Streptococcus, Granulicatella, Bacteroides, Porphyromonas, and Treponema [54, 55].
These studies were performed in different locations and hospitals, and recent studies have highlighted the role of geographical differences in the intestinal microbiota of preterm infants . We were unable to explore heterogeneity resulting from differences in DNA extraction protocols and sample processing methods due to considerable variation in included studies. The methodological variations were too great to provide statistical power for reasonable conclusions.
Strengths and limitations of the review
We adhered to the standard methods of conducting a systematic review as per recommended guidelines [23, 24]. We searched comprehensively for all eligible studies using clinically relevant inclusion criteria and strove to explain sources of heterogeneity by careful analysis of the study methods and reporting subgroup analyses. Unlike meta-analyses of randomized control trials, heterogeneity is a well-recognized problem in reviews of observational studies . Publication bias in studies reporting negative results may have excluded eligible unpublished studies. Poor reporting of study design, method of enrollment, and patient characteristics may hamper methodological assessment and external validity of the studies. The inconsistency in the microbiome profiles in the stools preceding NEC may be due to considerable clinical and methodological heterogeneity among the studies. In addition, while metadata was derived from both NCBI dbGaP as well as from investigators individually, not all metadata could be standardized and/or be available for incorporation into this analysis. The other reasons for inconsistency of results among studies may be differences in location emphasizing the importance of the unique NICU environment and associated infection control practices  and differences in time periods of study. We were not able to determine the time to event (NEC) analysis due to unavailability of day of NEC diagnosis from the meta-data. The time to event analysis may have enabled us to describe pre-NEC changes in microbiome more accurately.
Although there are many narrative reviews, to our knowledge, this is the first systematic review and meta-analysis of microbial dysbiosis preceding NEC. Our review supports the hypothesis that intestinal dysbiosis with an increase in microbes found primarily within the phylum Proteobacteria and a decrease in Firmicutes and Bacteroidetes precedes NEC in preterm infants. Furthermore, antibiotic usage is associated with increased abundance of Proteobacteria, which is also associated with higher incidence of NEC. In contrast, increased Proteobacteria were associated with intake of breast milk but type of feeding did not associate with development of NEC. The relationship between Proteobacteria and NEC seen in our review supports a strong association but should be viewed in the context of considerable variations in clinical and methodological characteristics of included studies. This certainly highlights the need for standardization of methods and transparent reporting of microbiome studies to advance this field. Our data with its limitations supports the microbial dysbiosis theory for the development of NEC but also highlights the lack of host response data to complete the story. Studies supporting a biologic plausibility will need to include metagenomics, transcriptomics, and additional data on inflammatory mediators. Proof of causality will need to further fulfill Koch’s postulates by recapitulating certain aspects of the disease in model systems. 16S rRNA gene sequencing often limits resolution to the level of the genus, whereas metagenomic data from whole genome shotgun sequencing may provide resolution at the strain level. Metagenomic analysis and metabolomics may reveal the differences in the functionality of the microbiome in patients who develop NEC. Optimizing the microbiome and correcting the microbial profile perturbations by microbial biotherapy may prevent NEC and improve clinical outcomes.
Neu J, Walker WA. Necrotizing enterocolitis. N Engl J Med. 2011;364:255–64.
Lin PW, Stoll BJ. Necrotising enterocolitis. Lancet. 2006;368:1271–83.
Fanaroff AA, Stoll BJ, Wright LL, Carlo WA, Ehrenkranz RA, Stark AR, Bauer CR, Donovan EF, Korones SB, Laptook AR, et al. Trends in neonatal morbidity and mortality for very low birthweight infants. Am J Obstet Gynecol. 2007;196:147. e141-148.
Lemons JA, Bauer CR, Oh W, Korones SB, Papile LA, Stoll BJ, Verter J, Temprosa M, Wright LL, Ehrenkranz RA, et al. Very low birth weight outcomes of the National Institute of Child health and human development neonatal research network, January 1995 through December 1996. NICHD Neonatal Research Network. Pediatrics. 2001;107:E1.
Holman RC, Stoll BJ, Curns AT, Yorita KL, Steiner CA, Schonberger LB. Necrotising enterocolitis hospitalisations among neonates in the United States. Paediatr Perinat Epidemiol. 2006;20:498–506.
Stoll BJ, Hansen NI, Adams-Chapman I, Fanaroff AA, Hintz SR, Vohr B, Higgins RD. Neurodevelopmental and growth impairment among extremely low-birth-weight infants with neonatal infection. JAMA. 2004;292:2357–65.
Claud EC, Walker WA. Bacterial colonization, probiotics, and necrotizing enterocolitis. J Clin Gastroenterol. 2008;42 Suppl 2:S46–52.
Nanthakumar NN, Fusunyan RD, Sanderson I, Walker WA. Inflammation in the developing human intestine: a possible pathophysiologic contribution to necrotizing enterocolitis. Proc Natl Acad Sci U S A. 2000;97:6043–8.
Nanthakumar N, Meng D, Goldstein AM, Zhu W, Lu L, Uauy R, Llanos A, Claud EC, Walker WA. The mechanism of excessive intestinal inflammation in necrotizing enterocolitis: an immature innate immune response. PLoS One. 2011;6:e17776.
Carlisle EM, Morowitz MJ. The intestinal microbiome and necrotizing enterocolitis. Curr Opin Pediatr. 2013;25:382–7.
Wang Y, Hoenig JD, Malin KJ, Qamar S, Petrof EO, Sun J, Antonopoulos DA, Chang EB, Claud EC. 16S rRNA gene-based analysis of fecal microbiota from preterm infants with and without necrotizing enterocolitis. ISME J. 2009;3:944–54.
Afrazi A, Sodhi CP, Richardson W, Neal M, Good M, Siggers R, Hackam DJ. New insights into the pathogenesis and treatment of necrotizing enterocolitis: toll-like receptors and beyond. Pediatr Res. 2011;69:183–8.
Musemeche CA, Kosloske AM, Bartow SA, Umland ET. Comparative effects of ischemia, bacteria, and substrate on the pathogenesis of intestinal necrosis. J Pediatr Surg. 1986;21:536–8.
Cotten CM, Taylor S, Stoll B, Goldberg RN, Hansen NI, Sanchez PJ, Ambalavanan N, Benjamin Jr DK. Prolonged duration of initial empirical antibiotic treatment is associated with increased rates of necrotizing enterocolitis and death for extremely low birth weight infants. Pediatrics. 2009;123:58–66.
Alexander VN, Northrup V, Bizzarro MJ. Antibiotic exposure in the newborn intensive care unit and the risk of necrotizing enterocolitis. J Pediatr. 2011;159:392–7.
Morowitz MJ, Poroyko V, Caplan M, Alverdy J, Liu DC. Redefining the role of intestinal microbes in the pathogenesis of necrotizing enterocolitis. Pediatrics. 2010;125:777–85.
Sodhi CP, Neal MD, Siggers R, Sho S, Ma C, Branca MF, Prindle Jr T, Russo AM, Afrazi A, Good M, et al. Intestinal epithelial toll-like receptor 4 regulates goblet cell development and is required for necrotizing enterocolitis in mice. Gastroenterology. 2012;143:708–18. e701-705.
Lu P, Sodhi CP, Hackam DJ. Toll-like receptor regulation of intestinal development and inflammation in the pathogenesis of necrotizing enterocolitis. Pathophysiology. 2014;21:81–93.
Gevers D, Knight R, Petrosino JF, Huang K, McGuire AL, Birren BW, Nelson KE, White O, Methe BA, Huttenhower C. The Human Microbiome Project: a community resource for the healthy human microbiome. PLoS Biol. 2012;10:e1001377.
Human Microbiome Project C. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.
Sze MA, Schloss PD. Looking for a signal in the noise: revisiting obesity and the microbiome. mBio. 2016;7(4):e01018–16. doi:10.1128/mBio.01018-16.
Bell MJ, Ternberg JL, Feigin RD, Keating JP, Marshall R, Barton L, Brotherton T. Neonatal necrotizing enterocolitis, therapeutic decisions based upon clinical staging. Ann Surg. 1978;187:1–7.
Moher D, Liberati A, Tetzlaff J, Altman DG, Group P. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. BMJ. 2009;339:b2535.
Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, Rennie D, Moher D, Becker BJ, Sipe TA, Thacker SB. Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis of observational studies in epidemiology (MOOSE) group. JAMA. 2000;283:2008–12.
Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vazquez-Baeza Y, Jansson JK, Gordon JI, Knight R. Meta-analyses of studies of the human microbiota. Genome Res. 2013;23:1704–14.
Viswanathan M, Berkman ND, Dryden DM, Hartling L. Assessing Risk of Bias and Confounding in Observational Studies of Interventions or Exposures: Further Development of the RTI Item Bank. Methods Research Report. (Prepared by RTI–UNC Evidence-based Practice Center under Contract No. 290-2007-10056-I). AHRQ Publication No. 13-EHC106-EF. Rockville: Agency for Healthcare Research and Quality; 2013. www.effectivehealthcare.ahrq.gov/reports/final.cfm. Accessed 27 June 2015.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Mshvildadze M, Neu J, Shuster J, Theriaque D, Li N, Mai V. Intestinal microbial ecology in premature infants assessed with non-culture-based techniques. J Pediatr. 2010;156:20–5.
Mai V, Young CM, Ukhanova M, Wang X, Sun Y, Casella G, Theriaque D, Li N, Sharma R, Hudak M, Neu J. Fecal microbiota in premature infants prior to necrotizing enterocolitis. PLoS One. 2011;6:e20647.
McMurtry VE, Gupta RW, Tran L, Blanchard EE, Penn D, Taylor CM, Ferris MJ. Bacterial diversity and Clostridia abundance decrease with increasing severity of necrotizing enterocolitis. Microbiome. 2015;3:11.
Morrow AL, Lagomarcino AJ, Schibler KR, Taft DH, Yu Z, Wang B, Altaye M, Wagner M, Gevers D, Ward DV, et al. Early microbial and metabolomic signatures predict later onset of necrotizing enterocolitis in preterm infants. Microbiome. 2013;1:13.
Normann E, Fahlen A, Engstrand L, Lilja HE. Intestinal microbial profiles in extremely preterm infants with and without necrotizing enterocolitis. Acta Paediatr. 2013;102:129–36.
Sim K, Shaw AG, Randell P, Cox MJ, McClure ZE, Li MS, Haddad M, Langford PR, Cookson WO, Moffatt MF, Kroll JS. Dysbiosis anticipating necrotizing enterocolitis in very premature infants. Clin Infect Dis. 2015;60:389–97.
Torrazza RM, Ukhanova M, Wang X, Sharma R, Hudak ML, Neu J, Mai V. Intestinal microbial ecology and environmental factors affecting necrotizing enterocolitis. PLoS One. 2013;8:e83304.
Zhou Y, Shan G, Sodergren E, Weinstock G, Walker WA, Gregory KE. Longitudinal analysis of the premature infant intestinal microbiome prior to necrotizing enterocolitis: a case-control study. PLoS One. 2015;10:e0118632.
Warner BB, Deych E, Zhou Y, Hall-Moore C, Weinstock GM, Sodergren E, Shaikh N, Hoffmann JA, Linneman LA, Hamvas A, et al. Gut bacteria dysbiosis and necrotising enterocolitis in very low birthweight infants: a prospective case–control study. Lancet. 2016;387(10031):1928–36.
de la Cochetiere MF, Piloquet H, des Robert C, Darmaun D, Galmiche JP, Roze JC. Early intestinal bacterial colonization and necrotizing enterocolitis in premature infants: the putative role of Clostridium. Pediatr Res. 2004;56:366–70.
Stewart CJ, Marrs EC, Magorrian S, Nelson A, Lanyon C, Perry JD, Embleton ND, Cummings SP, Berrington JE. The preterm gut microbiota: changes associated with necrotizing enterocolitis and infection. Acta Paediatr. 2012;101:1121–7.
Yee WH, Soraisham AS, Shah VS, Aziz K, Yoon W, Lee SK. Incidence and timing of presentation of necrotizing enterocolitis in preterm infants. Pediatrics. 2012;129:e298–304.
Gonzalez-Rivera R, Culverhouse RC, Hamvas A, Tarr PI, Warner BB. The age of necrotizing enterocolitis onset: an application of Sartwell’s incubation period model. J Perinatol. 2011;31:519–23.
Llanos AR, Moss ME, Pinzon MC, Dye T, Sinkin RA, Kendig JW. Epidemiology of neonatal necrotising enterocolitis: a population-based study. Paediatr Perinat Epidemiol. 2002;16:342–9.
La Rosa PS, Warner BB, Zhou Y, Weinstock GM, Sodergren E, Hall-Moore CM, Stevens HJ, Bennett Jr WE, Shaikh N, Linneman LA, et al. Patterned progression of bacterial populations in the premature infant gut. Proc Natl Acad Sci U S A. 2014;111:12522–7.
Sodhi C, Richardson W, Gribar S, Hackam DJ. The development of animal models for the study of necrotizing enterocolitis. Dis Model Mech. 2008;1:94–8.
Sodhi CP, Shi XH, Richardson WM, Grant ZS, Shapiro RA, Prindle Jr T, Branca M, Russo A, Gribar SC, Ma C, Hackam DJ. Toll-like receptor-4 inhibits enterocyte proliferation via impaired beta-catenin signaling in necrotizing enterocolitis. Gastroenterology. 2010;138:185–96.
Jilling T, Simon D, Lu J, Meng FJ, Li D, Schy R, Thomson RB, Soliman A, Arditi M, Caplan MS. The roles of bacteria and TLR4 in rat and murine models of necrotizing enterocolitis. J Immunol. 2006;177:3273–82.
Ward DV, Scholz M, Zolfo M, Taft DH, Schibler KR, Tett A, Segata N, Morrow AL. Metagenomic sequencing with strain-level resolution implicates uropathogenic E. coli in necrotizing enterocolitis and mortality in preterm infants. Cell Rep. 2016;14:2912–24.
Morowitz MJ, Denef VJ, Costello EK, Thomas BC, Poroyko V, Relman DA, Banfield JF. Strain-resolved community genomic analysis of gut microbial colonization in a premature infant. Proc Natl Acad Sci U S A. 2011;108:1128–33.
Zeissig S, Blumberg RS. Life at the beginning: perturbation of the microbiota by antibiotics in early life and its role in health and disease. Nat Immunol. 2014;15:307–10.
Gibson MK, Wang B, Ahmadi S, Burnham CA, Tarr PI, Warner BB, Dantas G. Developmental dynamics of the preterm infant gut microbiota and antibiotic resistome. Nat Microbiol. 2016;1:16024.
Barb JJ, Oler AJ, Kim HS, Chalmers N, Wallen GR, Cashion A, Munson PJ, Ames NJ. Development of an analysis pipeline characterizing multiple hypervariable regions of 16S rRNA using mock samples. PLoS One. 2016;11:e0148047.
Zheng W, Tsompana M, Ruscitto A, Sharma A, Genco R, Sun Y, Buck MJ. An accurate and efficient experimental approach for characterization of the complex oral microbiota. Microbiome. 2015;3:48.
Kumar PS, Brooker MR, Dowd SE, Camerlengo T. Target region selection is a critical determinant of community fingerprints generated by 16S pyrosequencing. PLoS One. 2011;6:e20956.
Taft DH, Ambalavanan N, Schibler KR, Yu Z, Newburg DS, Ward DV, Morrow AL. Intestinal microbiota of preterm infants differ over time and between hospitals. Microbiome. 2014;2:36.
Torrazza RM, Neu J. The altered gut microbiome and necrotizing enterocolitis. Clin Perinatol. 2013;40:93–108.
Schwiertz A, Gruhl B, Lobnitz M, Michel P, Radke M, Blaut M. Development of the intestinal bacterial composition in hospitalized preterm infants in comparison with breast-fed, full-term infants. Pediatr Res. 2003;54:393–9.
Bjorkstrom MV, Hall L, Soderlund S, Hakansson EG, Hakansson S, Domellof M. Intestinal flora in very low-birth weight infants. Acta Paediatr. 2009;98:1762–7.
Latuga MS, Stuebe A, Seed PC. A review of the source and function of microbiota in breast milk. Semin Reprod Med. 2014;32:68–73.
Sharon I, Morowitz MJ, Thomas BC, Costello EK, Relman DA, Banfield JF. Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization. Genome Res. 2013;23:111–20.
Grishin A, Papillon S, Bell B, Wang J, Ford HR. The role of the intestinal microbiota in the pathogenesis of necrotizing enterocolitis. Semin Pediatr Surg. 2013;22:69–75.
Raveh-Sadka T, Thomas BC, Singh A, Firek B, Brooks B, Castelle CJ, Sharon I, Baker R, Good M, Morowitz MJ, Banfield JF. Gut bacteria are rarely shared by co-hospitalized premature infants, regardless of necrotizing enterocolitis development. Elife. 2015; 4.
Millar MR, Linton CJ, Cade A, Glancy D, Hall M, Jalal H. Application of 16S rRNA gene PCR to study bowel flora of preterm infants with and without necrotizing enterocolitis. J Clin Microbiol. 1996;34(10):2506–10.
Smith B, Bodé S, Skov TH, Mirsepasi H, Greisen G, Krogfelt KA. Investigation of the early intestinal microflora in premature infants with/without necrotizing enterocolitis using two different methods. Pediatr Res. 2012;71(1):115–20.
We acknowledge the help of Toni-Ann Mistretta PhD in the initial stages of the study and for statistical advice. We also acknowledge statistical support from Joseph L. Hagan ScD.
Drs Tarr and Warner are funded by the NIH (NIAID; UH3AI083265) and the NIDDK-supported Digestive Diseases Research Core Center (P30DK052574) and the Washington University Institute of Clinical and Translational Sciences (UL1RR024992). Katherine Gregory is funded by the K23 NR011320-01 and Pilot/Feasibility Grant from Harvard Catalyst CTSA Grant UL1 RR 025758–01. AM is funded by the NIH grant NICHD R01 HD059140. Work in the Kroll laboratory was supported by funding from the Winnicott Foundation and the National Institute for Health Research (NIHR) Biomedical Research Centre based at Imperial Healthcare NHS Trust and Imperial College London.
Availability of data and materials
MP conceived the study, performed the systematic review, and wrote the manuscript. JC performed the quantitative synthesis of data and generation of figures. JN assisted in reviewing literature, guided the analyses, and assisted in writing the manuscript. EH and JV provided critical intellectual input in the study and the manuscript. VoM, PIT, BW, AM, KG, SK, VaM, MF, LE, and HEL provided data and intellectual input in the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This is a systematic review of studies and hence not applicable.
Search strategy for the systematic review in word document. (DOCX 14 kb)
Characteristics of included studies and results reported. This excel document describes the details the 14 included studies and results of alpha-, beta-diversity metrics and microbial profiles. (XLSX 16 kb)
Characteristics of studies included in the meta-analysis. This excel document describes the details of the data from 8 studies that were included in the meta-analysis. (XLSX 12 kb)
PRISMA checklist. This word document profiles the details of the systematic review that are recommended to be reported by the PRISMA statement. (DOC 63 kb)
Alpha diversity trends over corrected gestational age alpha-diversity metrics (mean with SD) by corrected gestational age (weeks). Fig. A depicts observed species (in operational taxonomic units, OTUs) and Fig. B depicts Shannon diversity index (SDI) against corrected gestational age (CGA) in controls and in infants with necrotizing enterocolitis (NEC). Data is represented in box and whisker plots with median and whiskers representing 10–90th centiles. The results of a fitted negative binomial regression model showed that the number of species increased significantly with gestational age (p < 0.001), and NEC patients tended to have fewer species than controls after controlling for gestational age (p = 0.053). (TIF 39 kb)
Comparison of microbiome profiles in NEC and control infants subgrouped by antibiotic usage before the onset of NEC. Observed species (operational taxonomic unit, OTU) richness (a) and Shannon diversity index (SDI) values (b) in cases and controls are subgrouped by antibiotic usage. Data is represented in box and whisker plots with median and whiskers representing 10–90th centiles. Significant differences were observed in OTU richness between controls who were exposed to antibiotics compared to controls who were not (*p < 0.05). No other significant differences were observed. Figs. c and d depict weighted and unweighted UniFrac distances in PCoA plots of NEC and controls subgrouped by antibiotic usage, and no notable clustering was observed. The proportion of variation explained along each axis is listed in parenthesis with the axis labels. Figs. e and f represent mean relative abundances of bacterial taxa in NEC cases and controls stratified by antibiotic usage at the phylum and genus levels, respectively. (TIF 712 kb)
Comparison of microbiota profiles in NEC and control infants subgrouped by diet. Observed species (operational taxonomic unit, OTU) richness (A) and Shannon diversity index (SDI) values (B) in cases and controls are subgrouped by diet (breast milk, formula, or both breastmilk and formula). Data is represented in box and whisker plots with median and whiskers representing 10–90th centiles. We observed significant differences in OTUs between these groups; breast milk-fed controls vs. breast milk-fed NEC infants (*p < 0.05) and both milk-fed controls vs. breast milk-fed controls (**p < 0.05). No other significant differences were observed. Figs. c and d depict weighted and unweighted PCoA plots of UniFrac distances of NEC cases and controls stratified by diet, and no notable clustering was observed. The proportion of variation explained along each axis is listed in parenthesis with the axis labels. Figs. e and f represent mean relative abundance distributions by NEC cases and controls subgrouped by diet at the phylum and genus levels, respectively. (TIF 931 kb)
Comparison of microbiota profiles in NEC and control infants subgrouped by mode of delivery. Observed species (operational taxonomic unit, OTU) richness (a) and Shannon diversity index (SDI) values (b) in cases and controls are subgrouped by mode of delivery: cesarean section (C-sec) or vaginal (Vag) delivery. Data is represented in box and whisker plots with median and whiskers representing 10–90th centiles. We observed significant differences in OTUs and SDI between controls delivered by C-section vs. vaginal delivery (*p < 0.05) and controls delivered by C-section vs. NEC infants delivered by vaginal delivery (**p <0.05). No other significant differences were observed. Figs. c and d depict weighted and unweighted PCoA plots of UniFrac distances of NEC cases and controls (all CGAs) stratified by mode of delivery, respectively, and no notable clustering was observed. The proportion of variation explained along each axis is listed in parenthesis with the axis labels. Figs. e and f represent mean relative abundance distributions by NEC cases and controls subgrouped by mode of delivery at the phylum and genus levels, respectively. (TIF 1002 kb)