Secondary Logo

Journal Logo

The constitution and functional prediction of the microbiota in necrotizing enterocolitis with a gestational age of over 28 weeks

Feng, JinXing MDa; He, Yu MDb; Liu, Dong MDc; Li, Luquan MDb; Chen, Jingyu MDd; Yu, Jialin MDe,∗

Section Editor(s): Kantarçeken., Bülent

doi: 10.1097/MD.0000000000017206
Research Article: Observational Study
Open

Background: To explore the features and function of gut microbiota in necrotizing enterocolitis patients over 28 gestational age weeks through a case-control study.

Methods: Fecal samples from patients with NEC over 28 gestational week age and matched control cases were collected. DNA of the fecal samples was extracted for 16 s rRNA sequencing to estimate the composition of the microbiota. Functional inference analyses were conducted through PICRUSt based on the sequencing raw data.

Results: There was no significant difference in the total diversity of microbiota between the fecal samples from the patients with NEC and the controls (P = .40). Propionibacterium was more abundant in the NEC cases than in the controls. Conversely, Lactobacillus, Phascolarctobacterium, and Streptococcus_salivarius were found to be more plentiful in the controls through LEfSe analysis. Functional inference analysis revealed that the xenobiotic biodegradation and metabolic activity was lower in the NEC cases than in the controls (P < .05).

Conclusion: The NEC cohort with a gestational age of over 28 weeks has a different pattern of microbiota compared with the controls. Functional inference analysis indicated that the potential function of the microbiota may also differ between these groups.

aDepartment of Neonatology, Shenzhen Children's Hospital, Shenzhen

bDepartment of Neonatology, Children's Hospital, Chongqing Medical University, Chongqing

cDepartment of Neonatology, Shenzhen People's Hospital, Shenzhen

dDepartment of Ultrasonography, Children's Hospital, Chongqing Medical University, Chongqing

eDepartment of Neonatology, Shenzhen University General Hospital, China.

∗Correspondence: Jialin Yu, Department of Neonatology, Shenzhen University General Hospital, 1098, Xueyuan Road, Nanshan District, Shenzhen 518055, China (e-mail: yujialin486@126.com).

Abbreviations: IL = interleukin, IQR = interquartile ranges, KEGG = Kyoto Encyclopedia of Genes and Genomes, KO = KEGG orthologs, LEfSe = Linear discriminant analysis effect size, NEC = necrotizing enterocolitis, OTUs = operational taxonomic units, PCoA = principal coordinate analysis, PICRUSt = Phylogenetic Investigation of Communities by Reconstruction of Unobserved States, PPAR = peroxisome proliferator-activated receptor, salivarius = Streptococcus salivarius, TLR = toll-like receptor.

How to cite this article: Feng J, He Y, Liu D, Li L, Chen J, Yu J. The constitution and functional prediction of the microbiota in necrotizing enterocolitis with a gestational age of over 28 weeks. Medicine. 2019;98:40(e17206).

JXF and YH contributed equally to this paper.

The authors have no conflicts of interest to disclose.

This is an open access article distributed under the terms of the Creative Commons Attribution-Non Commercial License 4.0 (CCBY-NC), where it is permissible to download, share, remix, transform, and buildup the work provided it is properly cited. The work cannot be used commercially without permission from the journal. http://creativecommons.org/licenses/by-nc/4.0

Back to Top | Article Outline

1 Introduction

Necrotizing enterocolitis (NEC) is among the most common life threatening diseases in neonates and occurs mostly in preterm babies.[1] NEC has become one of the most difficult diseases to eradicate with a high mortality and thus has been a priority for research.[2] Although the exact cause of NEC still remains unknown, it is widely accepted that the colonization of abnormal microbiota plays a dominant role in the pathogenesis of NEC. Early studies have been limited by the use of a conventional culture based technique which is unable to thoroughly describe the human gastrointestinal microbiota, as more than 80% of bacteria cannot be cultivated.[3] In the last decade, advances in the 16S small subunit bacterial ribosomal RNA (rRNA) gene based sequencing method has provided a more detailed picture of the composition of the human intestinal microbiota.[4] Preliminary studies compared the bacterial populations found in patients with NEC to those found in healthy individuals. However, there was variability in these results, and the association between an increase in Proteobacteria, Clostridia, Staphylococci and a decrease in Firmicutes have been reported in different studies.[5–11] Some studies have also emphasized the correlation between a decrease in the diversity of the microbiota and the development of NEC.[10,12] However, most of the references above focused on extremely preterm babies. NEC is still one of the most common diseases which cause deaths in late preterm babies or even in term babies. The relationship between microbiota and NEC in preterm babies of over 28 gestational weeks remains unclear.[13]

Moreover, most previous studies have only described the constitution of gut microbes in patients with NEC and little is known about the molecular mechanisms through which these microbes contribute to NEC pathogenesis and the function of these microbes. Metagenomics is the best way to describe the function of microbiota. However, the amount of DNA required for this technique limits its use in neonates as the volume of fecal sample is limited. Advances in the 16S analysis technology has allowed us to apply functional inference analysis, a computational approach to predict the functional composition of a metagenome using marker gene data and a database of reference genomes, as an alternative choice in the identification of the functional factors that differ between patients with NEC and healthy individuals.

In this study, we compared the intestinal microbiota in 16 premature infants with NEC to that of 16 matched and unaffected control infants over 28 weeks through the 16S rRNA sequencing method and conducted a functional inference analysis to determine the potential functional shift that is present in the intestinal microbiota between patients with NEC and unaffected infants.

Back to Top | Article Outline

2 Materials and methods

2.1 Patients and sample collection

This study was a part of microbiota research of preterm babies and was approved by the Institutional Review Board for Human Studies of the Shenzhen Children's Hospital. The study was carried out from May 2016 to February 2018. All infants with definite or advanced NEC, corresponding to the Bell stages II and III, and with gestational age over 28 weeks were included in the NEC cohort. Infants with congenital diseases were excluded from the study. Healthy infants who matched the infants with NEC in gestational age, birth weight, date of birth (+/−2 months), mode of delivery, and feeding patterns were enrolled as controls.

Fecal samples from the patients with NEC at 28 gestational weeks or over were collected by the medical staff once a diagnosis of NEC was confirmed, with an average collection time of 10 hours after the diagnosis of NEC. The samples of the control infants were collected at the same postnatal day. The freshly evacuated faeces were gathered from diapers into sterile tubes and were then transported to the laboratory immediately. All samples were stored at −80°C for further processing.

Back to Top | Article Outline

2.2 DNA extraction

DNA was extracted from 250 mg of the frozen faecal samples using the QIAamp FAST DNA Stool Mini-Kit (Qiagen, Germany) according to the manufacturer's instructions. The extracted DNA was eluted using the 50 μL ATE buffer from the kit.

Back to Top | Article Outline

2.3 PCR amplification and illumina MiSeq sequencing

DNA was amplified according to the V3-V4 region of bacterial 16S rRNA gene with universal bacteria primers: 338F (5’-ACTCCTACGG-GAGGCAGCA-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’), which contain an 8-base barcode sequence unique to each sample. The cleansed PCR products were pooled in equimolar amounts and submitted for paired-end sequencing (2 × 250) on an Illumina MiSeq platform according to the standard protocols.

Back to Top | Article Outline

2.4 Bioinformatic processing

Raw FASTq files were quality-filtered by QIIME according to the index sequence. Reads with over one nucleotide mismatch in the primer sequence, ambiguous characters, a length shortened by more than 50bp or those which could not be assembled were discarded. The sequences were binned into Operational Taxonomic Units (OTUs) using a similarity level threshold of 97%. The phylogenetic affiliation of each 16S rRNA gene sequence was analyzed by Ribosomal Database Project.

Next, we applied Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt), a technique that is based on the correlation between phylogeny and functions, to analyze the phylogenetic proximity between the public genomes and the 16S rRNA OTUs.[14,15] Additionally, the popular Kyoto Encyclopedia of Genes and Genomes (KEGG) was employed for annotation.[15] Functional predictions were exported as KEGG orthologs (KO). We obtained all KO annotations from the v3.5 of IMG25 to produce a table of 6909 abundances for all of the genomes that had identifiers in the Greengenes reference tree. Four levels of KO were used to compare the bacterial functional profiles.

Back to Top | Article Outline

2.5 Statistical analysis

Data were analyzed using SAS 9.0 (SAS Institute Inc, Cary, NC). Qualitative or categorical variables are expressed as frequencies and proportions while data with non-normal distribution are expressed as median and interquartile ranges (IQR). The differences between the two groups were tested using the Fisher Exact test for categorical variables and the nonparametric Wilcoxon paired test for continuous variables. The linear discriminant analysis effect size (LEfSe) was used to identify the phylogenetic features that differed significantly between all of the NEC cases and the controls[16] with an LDA score threshold of >3.0, based on the Kruskal-Wallis test. For all statistical analysis, a probability (P) value of <.05 was set as statistical significance.

Back to Top | Article Outline

3 Results

3.1 General information

During the study period a total of 21 patients were diagnosed with NEC. One case was excluded because of insufficient fecal sample collection, 2 because of congenital heart disease, and 2 patients died before sample collection could be completed. Finally, a total of 16 infants with NEC (NEC group) and 16 normal infants (control group) were enrolled in this study in the Shenzhen Children's Hospital, Shenzhen, China. The general demographic information is displayed in Table 1. Briefly, no differences were found between the NEC group and the control group for sex, gestational age, birth weight, delivery patterns, or feeding patterns. Among the samples, 5 cases progressed to stage III, 3 cases received surgery, and 2 infants died.

Table 1

Table 1

Back to Top | Article Outline

3.2 Microbiota analysis

A total of 32 fecal samples were analyzed through high-throughput sequencing. A total of 451,987 rRNA sequences with an average of 14,125 sequences per sample were obtained. These sequences were clustered at the 97.5% similarity level to obtain OTUs containing similar sequences for further microbiota analysis. The mean number of OTUs in each sample was 22.

As shown in Figure 1, the species richness, a diversity index that was determined by a Chao-based estimate of the total numbers of OTUs present, did not differ between the NEC and control cases (P = .60). Furthermore, we examined the samples using principal coordinate analysis (PCoA), which is based on UniFrac metric, to determine the microbiota structure. This analysis indicated that no clustering was detected according to infant status (Fig. 2). In the analysis of the overall microbiota structure for each group at the phylum and genus level, Proteobacteria and Firmicutes were the most abundant flora in both groups and these floras occupied over 90% of the total bacteria species composing the samples. However, no differences were found in the amount of Proteobacteria (P = .81) or Firmicutes (P = .99) between the 2 groups (Fig. 3). We, therefore, employed LEfSe analysis to investigate any difference in bacteria between the 2 groups. Interestingly, Propionibacteriales was more abundant in the NEC group than in the control group, while the abundance of Lactobacillus, Phascolarctobacterium and Streptococcus_salivarius was higher in the control group than in the NEC group (Fig. 4).

Figure 1

Figure 1

Figure 2

Figure 2

Figure 3

Figure 3

Figure 4

Figure 4

Back to Top | Article Outline

3.3 Functional inference analysis

Next, we performed functional inference analysis based on PICRUSt. By collapsing the data at KEEG level 2, we found that xenobiotic biodegradation and metabolism was significantly different between the groups (P < .05). Other methods of analysis, including environmental information processing, genetic information processing, and cellular processing did not show any differences between the NEC and control groups (Table 2).

Table 2

Table 2

Back to Top | Article Outline

4 Discussion

NEC has been identified as one of the most common causes of death in preterm babies at 26 to 28 gestational weeks.[17] As a result, most of the studies that have focused on the pathogenesis of NEC have studied extremely preterm babies. In contrast with the findings in developed countries, in developing countries, such as China, only 10% of babies that are born before 28 weeks survive. However, the survival rate of infants that are born at this stage is more than 90% in developed countries. Thus, the survival gap between countries is 10: 90.[18] Most patients with NEC in developing countries are born over 28 weeks. This study explored the constitution and function of the microbiota in patients with NEC over 28 weeks in China to provide the possible pathogenic factors.

In this study, our data were to some extent consistent with previous studies that identified no difference in the overall richness between patients with NEC (at 26–28 gestational weeks) and control infants.[7,9] However, at the phylum level, no changes were found in the Proteobacteria, Bacteroidetes, or Firmicutes between patients with NEC over 28 gestational weeks and the control group. This finding is different to those of previous studies.[7,9,11,19] This may be because the gut is well colonized in infants approaching 33 to 36 gestational weeks (the median gestational age was 34.8 for infants with NEC) in our study. Thus, the regular pattern for bacteria colonization in extremely premature babies may not apply to those babies at risk of NEC with a gestational age of over 28 weeks.[20]

Through IEfSe analysis, we noted that the level of Streptococcus (S.)_salivarius differed most significantly between the 2 groups. S. salivarius is one of the first colonizers of the gut after birth and, therefore, may contribute to the establishment of gut immune homeostasis and the modulation of host inflammatory responses. It has been reported that live S salivarius strains regulates innate intestinal immune inflammation through the prohibition of the NF-κB pathway and modulates peroxisome proliferator-activated receptor γ (PPARγ) transcriptional activity on intestinal epithelial cells.[21,22] The NF-κB pathway is a traditional signal involved in the promotion of transcription of proinflammatory factors and PPARγ has been shown to facilitate an anti-inflammatory antioxidant response that interacts with different enzymatic pathways, such as cyclooxygenase 2. Both of these factors have been reported to take part in the pathogenesis of NEC.[23–25] The other function of S._salivarius indicate that it may limit the expansion of other pathogens.[26] The above evidence could explain why we observed that S._salivarius was more abundant in the control group.

Meanwhile, it is not a surprise that we identified an increase in the abundance of Lactobacillus in the control group. Lactobacillus has been widely accepted to be a kind of probiotic. Many randomized clinic trials have demonstrated that Lactobacillus has a preventative action in NEC.[27–29] The mechanism though which Lactobacillus prevents NEC appears to involved the attenuation of interleukin (IL)-1β-induced IL-8 and IL-6 expression, decreased toll-like receptor (TLR)2, TLR4, and TLR9 expression, increased levels of specific negative regulators of inflammation (such as Single Ig IL-1-related receptor and Toll interacting protein), and the restoration of the ratio of CD4(+) Foxp3(+) regulatory T cells, which can inhibit the inflammatory activity of TH17 cells.[30–33] In addition, we discovered that Propionibacteriales was more abundant in the NEC group than in the control group, while Phascolarctobacterium was more abundant in the control group than in the NEC group. Both of these the bacteria can produce propionic acid that has been reported to maintain the intestinal immune homeostasis.[34–36] However, both of these bacteria occupied less than 0.1 percent of the total bacteria in our analysis. Thus, the influence of these bacteria may be very small.

Through function inference analysis, we also identified that NEC cases showed impaired xenobiotic biodegradation and metabolism. Xenobiotics contain drugs environmental pollutants, dietary supplements, and food additives.[37] Dysregulation of xenobiotic metabolism may contribute to some intestinal diseases, such as ulcerative colitis.[38] Thus, this may be the reason for the increased rate of the dysregulation of xenobiotic metabolism in patients with NEC.

The most valuable result of this study is the difference in the patterns of microbiota that were observed between patients with NEC and controls at a gestational age of over 28 weeks, in contrast to those of extremely preterm babies. One limitation for this research is that we did not explore the differences in microbiota between different gestational groups in the NEC cases below 28 gestational weeks. However, we could compare the composition of the microbiota that we observed with those that have been reported previously. Another limitation is our research is a small scale single-center study. The local geographic area may also affect the constitution of micobiota. Moreover, the fact that the incidence rate for NEC over 28 weeks’ gestation is low (previously reported as low as 0.2%) limits our sample size.[39] Thus, an independent multi-center validation cohort is necessary to confirm what we found in this study.

Back to Top | Article Outline

5 Conclusion

Lactobacillus and Streptococcus_salivarius are less abundant in NEC at over 28 gestational weeks than in controls. In addition, NEC cases presented with lower xenobiotic biodegradation and metabolism by functional prediction analysis.

Back to Top | Article Outline

Acknowledgment

The authors thank the members of the Department of Neonatology, Children's hospital, Chongqing medical university, for their help analyzing the data.

Back to Top | Article Outline

Author contributions

Conceptualization: JinXing Feng, Jialin Yu.

Formal analysis: Yu He.

Funding acquisition: JinXing Feng, Dong Liu.

Investigation: JinXing Feng.

Resources: Dong Liu, Luquan Li, Jingyu Chen.

Supervision: Luquan Li, Jingyu Chen.

Writing – original draft: Yu He.

Writing – review & editing: Jialin Yu.

Back to Top | Article Outline

References

[1]. He Y, Zhong Y, Yu J, et al. Ultrasonography and radiography findings predicted the need for surgery in patients with necrotising enterocolitis without pneumoperitoneum. Acta Paediatr 2016;105:e151–5.
[2]. Neu J, Walker WA. Necrotizing enterocolitis. N Engl J Med 2011;364:255–64.
[3]. Eckburg PB, Bik EM, Bernstein CN, et al. Diversity of the human intestinal microbial flora. Science (New York, NY) 2005;308:1635–8.
[4]. Relman DA. The search for unrecognized pathogens. Science (New York, NY) 1999;284:1308–10.
[5]. Morrow AL, Lagomarcino AJ, Schibler KR, et al. Early microbial and metabolomic signatures predict later onset of necrotizing enterocolitis in preterm infants. Microbiome 2013;1:13.
[6]. Grishin A, Bowling J, Bell B, et al. Roles of nitric oxide and intestinal microbiota in the pathogenesis of necrotizing enterocolitis. J Pediatr Surg 2016;51:13–7.
[7]. Torrazza RM, Ukhanova M, Wang X, et al. Intestinal microbial ecology and environmental factors affecting necrotizing enterocolitis. PLoS One 2013;8:e83304.
[8]. Stewart CJ, Marrs ECL, Magorrian S, et al. The preterm gut microbiota: changes associated with necrotizing enterocolitis and infection. Acta Paediatr 2012;101:1121–7.
[9]. Mai V, Young CM, Ukhanova M, et al. Fecal microbiota in premature infants prior to necrotizing enterocolitis. PLoS One 2011;6:e20647.
[10]. Wang Y, Hoenig JD, Malin KJ, et al. 16S rRNA gene-based analysis of fecal microbiota from preterm infants with and without necrotizing enterocolitis. ISME J 2009;3:944–54.
[11]. Warner BB, Deych E, Zhou Y, et al. Gut bacteria dysbiosis and necrotising enterocolitis in very low birthweight infants: a prospective case-control study. Lancet 2016;387:1928–36.
[12]. Sanz Y, Stewart CJ, Marrs ECL, et al. Development of the preterm gut microbiome in twins at risk of necrotising enterocolitis and sepsis. PLoS One 2013;8:e73465.
[13]. Tomashek KM, Shapiro-Mendoza CK, Davidoff MJ, et al. Differences in mortality between late-preterm and term singleton infants in the United States, 1995–2002. J Pediatr 2007;151456:e451450-456.
[14]. Langille MGI, Zaneveld J, Caporaso JG, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol 2013;31:814–21.
[15]. Arboleya S, Sánchez B, Solís G, et al. Impact of prematurity and perinatal antibiotics on the developing intestinal microbiota: a functional inference study. Int J Mol Sci 2016;17:649.
[16]. Segata N, Izard J, Waldron L, et al. Metagenomic biomarker discovery and explanation. Genome Biol 2011;12:R60.
[17]. Patel RM, Kandefer S, Walsh MC, et al. Causes and timing of death in extremely premature infants from 2000 through 2011. N Engl J Med 2015;372:331–40.
[18]. Blencowe H, Cousens S, Oestergaard MZ, et al. National, regional, and worldwide estimates of preterm birth rates in the year 2010 with time trends since 1990 for selected countries: a systematic analysis and implications. Lancet (London, England) 2012;379:2162–72.
[19]. Mshvildadze M, Neu J, Shuster J, et al. Intestinal microbial ecology in premature infants assessed with non–culture-based techniques. J Pediatr 2010;156:20–5.
[20]. La Rosa PS, Warner BB, Zhou Y, et al. Patterned progression of bacterial populations in the premature infant gut. Proc Natl Acad Sci U S A 2014;111:12522–7.
[21]. Kaci G, Goudercourt D, Dennin V, et al. Anti-Inflammatory Properties of Streptococcus salivarius, a commensal bacterium of the oral cavity and digestive tract. Appl Environ Microbiol 2014;80:928–34.
[22]. Schuch R, Couvigny B, de Wouters T, et al. Commensal Streptococcus salivarius Modulates PPARγ transcriptional activity in human intestinal epithelial cells. PLoS One 2015;10:e0125371.
[23]. Baregamian N, Mourot JM, Ballard AR, et al. PPAR-γ agonist protects against intestinal injury during necrotizing enterocolitis. Biochem Biophys Res Commun 2009;379:423–7.
[24]. Corsini I, Polvani S, Tarocchi M, et al. Peroxisome proliferator-activated receptor-γ agonist pioglitazone reduces the development of necrotizing enterocolitis in a neonatal preterm rat model. Pediatr Res 2016;81:364–8.
[25]. Hunter CJ, De Plaen IG. Inflammatory signaling in NEC: role of NF-κB, cytokines and other inflammatory mediators. Pathophysiology 2014;21:55–65.
[26]. Doyle H, Pierse N, Tiatia R, et al. The effect of the oral probiotic streptococcus salivarius (K12) on group a streptococcus pharyngitis. Pediatr Infect Dis J 2017;1.
[27]. Braga TD, da Silva GA, de Lira PI, et al. Efficacy of Bifidobacterium breve and Lactobacillus casei oral supplementation on necrotizing enterocolitis in very-low-birth-weight preterm infants: a double-blind, randomized, controlled trial. Am J Clin Nutr 2011;93:81–6.
[28]. Sari FN, Dizdar EA, Oguz S, et al. Oral probiotics: Lactobacillus sporogenes for prevention of necrotizing enterocolitis in very low-birth weight infants: a randomized, controlled trial. Eur J Clin Nutr 2011;65:434–9.
[29]. Fernandez-Carrocera LA, Solis-Herrera A, Cabanillas-Ayon M, et al. Double-blind, randomised clinical assay to evaluate the efficacy of probiotics in preterm newborns weighing less than 1500 g in the prevention of necrotising enterocolitis. Arch Dis Child Fetal Neonatal Ed 2013;98:F5–9.
[30]. Ganguli K, Meng D, Rautava S, et al. Probiotics prevent necrotizing enterocolitis by modulating enterocyte genes that regulate innate immune-mediated inflammation. Am J Physiol Gastrointest Liver Physiol 2013;304:G132–141.
[31]. Liu Y, Tran DQ, Fatheree NY, et al. Lactobacillus reuteri DSM 17938 differentially modulates effector memory T cells and Foxp3+ regulatory T cells in a mouse model of necrotizing enterocolitis. Am J Physiol Gastrointest Liver Physiol 2014;307:G177–186.
[32]. Liu Y, Fatheree NY, Dingle BM, et al. Lactobacillus reuteri DSM 17938 changes the frequency of Foxp3+ regulatory T cells in the intestine and mesenteric lymph node in experimental necrotizing enterocolitis. PLoS One 2013;8:e56547.
[33]. Good M, Sodhi CP, Ozolek JA, et al. Lactobacillus rhamnosus HN001 decreases the severity of necrotizing enterocolitis in neonatal mice and preterm piglets: evidence in mice for a role of TLR9. Am J Physiol Gastrointest Liver Physiol 306:G1021–32.
[34]. Maji A, Misra R, Dhakan DB, et al. Gut microbiome contributes to impairment of immunity in pulmonary tuberculosis patients by alteration of butyrate and propionate producers. Environ Microbiol 2018;20:402–19.
[35]. Tax G, Urbán E, Palotás Z, et al. Propionic acid produced by propionibacterium acnes strains contributes to their pathogenicity. Acta Derm Venereol 2016;96:43–9.
[36]. Corrêa-Oliveira R, Fachi JL, Vieira A, et al. Regulation of immune cell function by short-chain fatty acids. Clin Transl Immunology 2016;5:e73.
[37]. Johnson CH, Patterson AD, Idle JR, et al. Xenobiotic metabolomics: major impact on the metabolome. Annu Rev Pharmacol Toxicol 2012;52:37–56.
[38]. Langmann T, Moehle C, Mauerer R, et al. Loss of detoxification in inflammatory bowel disease: dysregulation of pregnane X receptor target genes. Gastroenterology 2004;127:26–40.
[39]. Battersby C, Modi N. Challenges in advancing necrotizing enterocolitis research. Clin Perinatol 2019;46:19–27.
Keywords:

functional inference; microbiota; necrotizing enterocolitis

Copyright © 2019 the Author(s). Published by Wolters Kluwer Health, Inc.