Temporal bacterial and metabolic development of the preterm gut reveals specific signatures in health and disease
© The Author(s). 2016
Received: 25 June 2016
Accepted: 1 December 2016
Published: 29 December 2016
The preterm microbiome is crucial to gut health and may contribute to necrotising enterocolitis (NEC), which represents the most significant pathology affecting preterm infants. From a cohort of 318 infants, <32 weeks gestation, we selected 7 infants who developed NEC (defined rigorously) and 28 matched controls. We performed detailed temporal bacterial (n = 641) and metabolomic (n = 75) profiling of the gut microbiome throughout the disease.
A core community of Klebsiella, Escherichia, Staphyloccocus, and Enterococcus was present in all samples. Gut microbiota profiles grouped into six distinct clusters, termed preterm gut community types (PGCTs). Each PGCT reflected dominance by the core operational taxonomic units (OTUs), except of PGCT 6, which had high diversity and was dominant in bifidobacteria. While PGCTs 1–5 were present in infants prior to NEC diagnosis, PGCT 6 was comprised exclusively of healthy samples. NEC infants had significantly more PGCT transitions prior to diagnosis. Metabolomic profiling identified significant pathways associated with NEC onset, with metabolites involved in linoleate metabolism significantly associated with NEC diagnosis. Notably, metabolites associated with NEC were the lowest in PGCT 6.
This is the first study to integrate sequence and metabolomic stool analysis in preterm neonates, demonstrating that NEC does not have a uniform microbial signature. However, a diverse gut microbiome with a high abundance of bifidobacteria may protect preterm infants from disease. These results may inform biomarker development and improve understanding of gut-mediated mechanisms of NEC.
KeywordsPreterm infant Gut microbiome Necrotising enterocolitis Metabolomics
Survival after preterm birth is increasing, but this is associated with increased numbers of infants developing necrotising enterocolitis (NEC) and late onset sepsis (LOS), with ~10% and ~20% of very low birth weight (VLBW) infants being affected, respectively . NEC and LOS are now responsible for more deaths after the first week of life in extremely preterm infants than any other single pathology , and both may result in significantly impaired outcomes in survivors. Antecedent risk factors for NEC vary, and the pathophysiology remains poorly defined. Currently accepted staging, e.g. Bell’s criteria, is deemed outdated  and there is increasing acceptance that a clinical diagnosis of “NEC” may simply be the final common pathway of multiple disease processes and LOS from gut bacterial translocation sharing similar processes .
Bacterial colonisation is necessary for the development of NEC, and specific bacterial taxa have been associated with onset, particularly those from the Proteobacteria phylum [5, 6] such as Enterobacter , Escherichia , Sphingomonas , and Klebsiella sp . Studies from different research groups (and hence different neonatal intensive care units) identify a variety of signatures potentially associated with disease states. The two largest studies exploring NEC have shown conflicting findings, with no differences between NEC and control infants (n = 163 patients (21 with NEC), 482 samples) in a denaturing gradient gel electrophoresis (DGGE) and culture exploration  and more recently an increased relative abundance of Gammaproteobacteria and reduced Negativicutes (n = 166 patients, 3587 stools) . Moving beyond 16S rRNA bacterial profiling, a metagenomic (shotgun sequencing) study of NEC samples during an apparent outbreak found no distinct microbial strain in NEC infants at the time of diagnosis . Studies have also reported an absence of “pathogenic” bacteria in tissue resections from NEC infants  and no difference in the total bacterial load .
Previous studies utilising metabolomics to understand NEC pathogenesis and identify potential biomarkers are limited, with no published data on stool metabolomics. However, Morrow et al. (2003) coupled bacterial profiling of stool to metabolomic analysis of matched urine, showing early dysbiosis of the bacterial community and a high alanine to histidine ratio that was associated with microbial characteristics and good overall prediction of NEC . Other metabolomic studies have utilised serum/plasma, employing either LCMS or gas chromatograph MS (GCMS), showing metabolomic profiles differ in disease compared to controls, but robust biomarkers were not found [15, 16].
To further progress understanding of the bacterial and host processes involved in NEC, we recruited a large cohort of preterm infants sampling stool daily where possible (n = >300 patients/>3000 samples) capturing key health related issues with precise definitions. Using strict classification for NEC and only including patients with robust temporal sampling before and after disease diagnosis, we present comprehensive longitudinal gut microbiome data on 641 samples from those with disease (n = 7) and well-matched non-diseased controls (n = 28). A subset of 75 samples (n = 16) also underwent ultra performance liquid chromatography mass spectrometry (UPLC-MS). We hypothesised that there would be detectable and coherent differences in bacterial and metabolomic profiles between healthy and diseased infants before the diagnosis of NEC, with associated bacterial signatures potentially predictive of health or disease onset.
Enterobacteriaceae, Escherichia, Enterococcus, and Staphylococcus were present in every sample (core microbiome). The unclassified Enterobacteriaceae was cultured and identified by Matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF-MS) and full-length 16S sequencing as Klebsiella (K. oxytoca and K. pneumoniae). The developing gut microbiome in patients diagnosed with NEC was highly dynamic and individual, with no clear causative organism identified in the patients diagnosed with NEC (Additional file 1: Figure S1 and Additional file 2: Table S1). Comparably, the bacterial profiles from each control patient showed a high degree of intra- and inter-patient variability, with large temporal shifts in the core and satellite operational taxonomic units (OTUs) (Additional file 3: Figure S2).
Interestingly, the relative abundance of Staphylococcus decreased in all groups following birth, with few samples containing a high relative abundance following day 25 of life. This was reflected in PGCT 3 (Staphylococcus dominated), which was the only PGCT associated with age. Additionally, no statistical difference in the relative abundance of any OTU was found between vaginal and caesarean section delivery from samples collected in the initial week of life. Similarly for the alpha diversity, while the Shannon diversity of vaginally delivered infants was higher than those delivered by caesarean section, this increase was not significant (Additional file 8: Figure S7). No significant difference was also found in the alpha diversity of the gut microbiome between infants born at different gestational ages (P = 0.16).
The extensive longitudinal sampling employed in this study demonstrated the extent of intra- and inter-individual variation in the development of the preterm gut microbiota in both health and disease. To facilitate comparisons between samples prior to disease diagnosis and from healthy (control) infants, we clustered samples into one of six specific PGCT based on the bacterial profiles. The PGCT analysis showed the majority of patients exhibited multiple PGCTs during the study, with only one patient (control 234) having the same PGCT (PGCT 1) throughout the entirety of sampling. Indeed, controls generally showed significantly increased stability compared to NEC infants prior to disease diagnosis.
No specific microbial signature was found prior to the diagnosis of NEC, compared to matched healthy controls. Although no microorganism has been reproducibly associated with NEC, many studies have reported potential associations with specific taxa [7, 9, 10, 20, 21], including a recent large multicentre study which found increased Gammaproteobacteria and Negativicutes prior to NEC . The latter study conducted analysis at the class level, whereas the current study used genus level taxonomic classification. Notably, analysing the current dataset at class level showed no significant difference in the relative abundance of Gammaproteobacteria (class) between the preNEC and control samples (Additional file 9: Figure S8). It should also be noted that by nature of “association” between specific taxa and NEC, the reported taxa are not ubiquitously increased in NEC and will also be present in comparable abundance in some healthy controls, which is in accordance with our data. Thus, the presence of such taxa alone may not improve our understanding of NEC pathogenesis. Our study is also consistent with another large preterm cohort (163 infants, 21 with NEC, analysing 482 samples) who report no associations using molecular DGGE profiling  and a recent powerful metagenomic study with analysis to the strain level [12, 22]. The latter metagenomic study highlights the importance of lower taxonomic level analysis in disease association studies. Further characterisation of Escherichia coli isolates from NEC and control infants in this study support the finding that each infant harboured unique strains (Additional file 10: Figure S9). While no specific microbial signature or taxa was evident in disease, the instability of the developing gut microbiome in NEC infants may be a “trigger” initiating and promoting a constant exaggerated inflammatory response, ultimately leading to tissue damage. Such findings need to be validated and quantified in independent cohorts, using rigorous longitudinal sampling prior to disease and similarly in controls.
Although no specific taxa was elucidated in bacterial profiling, metabolomic profiling showed that the five most discriminatory metabolites in the NEC samples at the time of diagnosis were significantly elevated compared to all the control samples. While the metabolites were most significant at diagnosis, the metabolites increased in intensity prior to diagnosis and reduced following NEC diagnosis. The temporal increase in metabolites prior to diagnosis, whilst levels in the control samples remained consistently low, suggests such metabolites may be predictive of disease onset. This finding offers potential in biomarker development for detecting the risk NEC, which precedes the clinical diagnosis. Based on the findings here, the progression toward NEC may be detected between 1–2 weeks prior to current clinical diagnosis, which is consistent with the previous data . Such findings require validation in large multicentre studies, but the finding that 4/5 metabolites most associated with NEC were involved in linoleate metabolism is intriguing given that trans-10, cis-12-conjugated linoleic acid (t10c12-CLA) negatively impacts on lipid metabolism, mediated by alterations in gut microbiome . As more metabolomic investigations are performed, the mechanism of linoleate metabolism in an inflammatory-mediated damage may be better elucidated, but the results of this study support the concept that it is, at least in part, mediated by the microbiome.
Analysing a large number of samples from healthy control infants revealed a novel healthy microbial signature (PGCT 6), characterised by dominance of bifidobacteria and a significantly higher overall diversity. In accordance with existing data, we found that control infants have increasing diversity whereas patients diagnosed with NEC have significantly reduced diversity while on the NICU [10, 19, 20], although the microbiome of NEC infants converge to that of controls within 1 year following NICU discharge . As well as bifidobacteria, PGCT 6 also had significantly higher relative abundance of lactobacilli. Interestingly, because of the associated health benefits, these genera are the most common bacteria used in manufactured probiotic products . The metabolomic profiling showed high concordance with the 16S bacterial profiling results. Strikingly, the metabolites most associated with NEC were significantly reduced in PGST 6 and the metabolites most associated with control were increased in PGST 6. This suggests that changes in the microbial community have a direct impact on the metabolic functioning on the host, which may ultimately contribute to disease pathophysiology.
Despite the recent trial showing no effect of supplemental Bifidobacterium breve , probiotics containing Bifidobacterium spp. have generally shown reduced incidences of NEC [27, 28] and bifidobacteria are known to play fundamental roles in the development and maintenance of a healthy gut microbiome [29, 30]. Thus, the high abundance of bifidobacteria may contribute to the improved health status by aiding gut maturation and preventing over-expression of potentially damaging pathways, although it is also possible that a higher abundance of bifidobacteria simply acts as a marker for a more healthy gut, rather than being directly causative in preterm infants. This is supported by the finding in this study that prostaglandin formation is significantly increased at NEC diagnosis and reduced in PGCT 6, with this pathway known to mediate inflammation [31, 32]. Caesarean or vaginal delivery did not significantly alter the preterm gut microbiome, even in the initial weeks of life, although PGCT 6 did occur more frequently in vaginally delivered infants compared to caesarean section (42 vs. 25, respectively).
Despite the collection of >3000 samples from >300 infants, the number of infants diagnosed with a disease from this cohort remains small due to both application of a strict definition of disease states and the inclusion criteria requiring robust temporal sampling prior to and following disease diagnosis. The extent of inter-individual microbial variability in both healthy and diseased infants suggests that no specific bacterial signature is likely to cause NEC in our unit, and this finding is unlikely to be influenced by additional cases of NEC.
No consistent microbial signature or taxa was associated with NEC, but the number of PGCT transitions, and thus microbiome instability, was significantly increased in NEC and may contribute to, or be a marker of, the inflammatory-mediated damage in NEC pathogenesis. A novel diverse bacterial community with high relative abundance of bifidobacteria was found only in healthy samples. The metabolomics and 16S data was highly concordant, where metabolites associated with NEC were significantly reduced in the healthy PGCT. Strategies aiming to improve outcomes in preterm infants which focus on attempts to promote healthy PGCT profiles, rather than focusing on PGCTs potentially associated with later disease, are warranted.
Participants and study design
Summary of patient samples and demographic per group
Control (n = 28)
NEC (n = 7)
Number of stool samples
Birth weight (g)a
Birth mode (CS/vaginal)
Antibiotic prediagnosis (days)a
Antibiotic total (days)a
Seven well-sampled cases of NEC were selected and each matched to four well-sampled controls based on gestational age (GA; +/− 1 week), birth weight, and delivery mode. The 35 infants contributed a total 641 analysed stool samples. Samples were stored immediately at −20 °C and long-term at −80 °C, with all samples having comparable collection/storage. 16S rRNA gene bacterial profiling was performed on all the samples in the study. Metabolomic profiling was performed on a subset of 16 infants (75 stools), 6 NEC and 10 matched controls. Samples were selected for analysis relative to disease diagnosis at day of life (DOL) −14 (time point 1; TP1), −7 (TP2), 0 (TP3), +7 (TP4), and +14 (TP5), and matched by DOL alone in controls.
16S rRNA gene Bacterial Profiling
Nucleic acid extraction of stool was carried out on 100 mg of sample using the PowerLyzer™ PowerSoil® DNA Isolation Kit (MoBio, CA, USA) in accordance with the manufacturer’s instructions. Bacterial profiling utilised the 16S rRNA gene targeting variable region four based on the Schloss wet-lab MiSeq SOP and resulting raw fastq data were processed using Mothur (version 1.31.2) as described previously . Briefly, combined reads were trimmed to 275 reads with 0 ambiguous bases. Chimeric sequences were detected by Chimera.uchime and removed from downstream analysis. Alignment was generated via the Silva v4 database  and Chloroplast, Mitochondria, unknown, Archaea, and Eukaryota linages were removed from the analysis. In total, 44,515,418 reads were obtained and raw sequences were deposited in MG-RAST under the accession numbers 4516545.3-4516585.3.
UPLC-MS metabolomic profiling
Metabolomic profiling was performed as previously described . Briefly, 100 mg stool was homogenised (80% methanol), vortexed for 15 min, centrifuged (10,000 × g), and lyophilised. Reverse-phase ultra-performance LCMS tandem mass spectrometry (UPLC-MS/MS) was performed using an Accucore C18 column (2.6 μm, 150 × 2.1 mm) at 40 °C, 3.0 μl injection, and 300 μl/min flow rate. Gradients increased from 5% acetronitrile (ACN) to 95% ACN over 22 min, followed by 8 min wash and re-equilibration. Samples were run randomly in triplicate on a Q-Exactive (Thermo) using HESI with high resolution (70,000) positive and negative switching. The mass range was set from 100–1000 m/z. SIEVE (Version 2.2 beta) was used to process the Thermo RAW files by component extraction detecting a total of 11,612 components (8343 positive and 3269 negative).
Bioinformatic and statistical analysis
Analysis and visualization of microbiome communities was conducted in R , utilizing the phyloseq package  to import sample data and calculate alpha- and beta-diversity metrics. Each sample was rarefied to 4397 reads. Significance of categorical variables were determined using the non-parametric Mann-Whitney test for two category comparisons or the Kruskal-Wallis test when comparing three or more categories. All p values were adjusted for multiple comparisons with the FDR algorithm.
An R script was implemented and made publicly available by DiGiulio et al.  was employed for linear mixed-effects modelling, partition around medoid based clustering, and Markov chain modelling to cluster all the samples within the study into PGCTs and determine temporal development in health and disease. Bray-Curtis was used to calculate the distance between all the samples, and this was denoised by extraction of the most significant principal coordinates analysis (PCoA) eigenvectors before applying the PAM algorithm. Gap statistic was used to determine the number of clusters. Ten days prior to diagnosis of disease, network analysis was used to classify PreNEC, based on previous data . To determine if the stability was potentially contributing to the disease pathophysiology, the total number of unique PGCTs and the number of transitions between PGCTs were compared for NEC patients, to matched controls. Only the pre diagnosis samples and control samples in the first 28 days of life were included (last predisease sample was day of life 28).
UPLC-MS data was filtered to include only m/z features that occurred in >20% of samples. Metabolite annotation and pathway enrichment was performed using Mummichog . Mummichog was used to determine significant pathways at TP3 by comparing NEC vs. control. Only pathways identified in Mummichog were included in subsequent analysis. SIMCA 13.0 (Umetrics, Stockholm, Sweden) was used to model the metabolomics dataset by first plotting an unsupervised principal coordinate analysis (PCA) and if clear grouping was apparent a orthogonal partial least squared-discriminatory analysis (OPLS-DA) was performed. Only OPLS-DA with R2 > 0.8 and Q2 > 0.5 were considered robust and predictive, respectively. The variable importance plot (VIP) was used to determine the five metabolites most associated with each class . MiniTab 17 was used to generate box plots and determine the significance by Tukey’s method for multiple comparisons.
Late onset sepsis
Neonatal intensive care unit
Preterm gut community type
Ultra-performance liquid chromatography mass spectrometry tandem mass spectrometry
This work was supported by funding from the Tiny Lives charity (Newcastle upon Tyne, UK), Newcastle upon Tyne Hospitals NHS Charity, and in part by an unrestricted educational grant from Nestle UK. The content is solely the responsibility of the authors.
Availability of data and materials
The data set supporting the results of this article is available in the MG-RAST repository [accession numbers 4516545.3-4516585.3].
CJS, JEB, NDE, and SPC were involved in the study concept and design. CJS, AJP, BA, TS, AN, ECLM, DPS, JFP, and JP were all involved in the collection and sorting of the samples and data acquisition and analysis. CJS, JEB, NDE, and SPC were primarily involved in drafting the manuscript, and all authors contributed to critical revisions and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Ethical approval was obtained from the County Durham and Tees Valley Research Ethics Committee. Parental informed consent was given.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- The SIFT Investigators Group. Early enteral feeding strategies for very preterm infants: current evidence from Cochrane reviews. Arch Dis Child Fetal Neonatal Ed. 2013;98:F470–2.View ArticleGoogle Scholar
- Neu J, Walker WA. Necrotizing enterocolitis. NEJM. 2011;364:255–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Jakaitis BM, Denning PW. Commensal and probiotic bacteria may prevent NEC by maturing intestinal host defenses. Pathophysiology. 2014;21:47–54.View ArticlePubMedGoogle Scholar
- Gordon P, Christensen R. Mapping the new world of necrotizing enterocolitis (NEC): review and opinion. EJ Neonatol Res. 2012;2:145–72.PubMedPubMed CentralGoogle Scholar
- Morrow AL, Lagomarcino AJ, Schibler KR, Taft DH, Yu Z, Wang B, et al. Early microbial and metabolomic signatures predict later onset of necrotizing enterocolitis in preterm infants. Microbiome. 2013;1:13.View ArticlePubMedPubMed CentralGoogle Scholar
- Warner BB, Deych E, Zhou Y, Hall-Moore C, Weinstock GM, Sodergren E, et al. Gut bacteria dysbiosis and necrotising enterocolitis in very low birthweight infants: a prospective case-control study. Lancet. 2016;387:1928–36. Elsevier.View ArticlePubMedGoogle Scholar
- Stewart C, Marrs E, Magorrian S, Nelson A, Lanyon C, Perry J, et al. The preterm gut microbiota: changes associated with necrotizing enterocolitis and infection. Acta Paediatr. 2012;101:1121–7.View ArticlePubMedGoogle Scholar
- Jenke AC, Postberg J, Mariel B, Hensel K, Foell D, Däbritz J, et al. S100A12 and hBD2 Correlate with the composition of the fecal microflora in ELBW infants and expansion of E. coli is associated with NEC. Biomed Res Int. 2013;2013:150372.View ArticlePubMedPubMed CentralGoogle Scholar
- Stewart CJ, Nelson A, Scribbins D, Marrs ECL, Perry JD, Embleton ND, et al. Bacterial and fungal viability in the preterm gut: NEC and sepsis. Arch Dis Child Fetal Neonatal Ed. 2013;98:F298–303.View ArticlePubMedGoogle Scholar
- Sim K, Shaw AG, Randell P, Cox MJ, McClure ZE, Li M-S, et al. Dysbiosis anticipating necrotizing enterocolitis in very premature infants. Clin Infect Dis. 2014;60:389–97.
- 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:115–20. Nature Publishing Group.View ArticlePubMedGoogle Scholar
- Raveh-Sadka T, Thomas BC, Singh A, Firek B, Brooks B, Castelle CJ, et al. Gut bacteria are rarely shared by co-hospitalized premature infants, regardless of necrotizing enterocolitis development. Elife. 2015;4:1–25.View ArticleGoogle Scholar
- Smith B, Bodé S, Petersen BL, Jensen TK, Pipper C, Kloppenborg J, et al. Community analysis of bacteria colonizing intestinal tissue of neonates with necrotizing enterocolitis. BMC Microbiol. 2011;11:73.View ArticlePubMedPubMed CentralGoogle Scholar
- Abdulkadir B, Nelson A, Skeath T, Marrs ECL, Perry JD, Cummings SP, et al. Stool bacterial load in preterm infants with necrotising enterocolitis. Early Hum Dev. 2016;95:1–2.View ArticlePubMedGoogle Scholar
- Wilcock A, Begley P, Stevens A, Whatmore A, Victor S. The metabolomics of necrotising enterocolitis in preterm babies: an exploratory study. J Matern Fetal Neonatal Med. 2016;29:758–62.View ArticlePubMedGoogle Scholar
- Stewart CJ, Nelson A, Treumann A, Skeath T, Cummings SP, Embleton ND, et al. Metabolomic and proteomic analysis of serum from preterm infants with necrotising entercolitis and late onset sepsis. Pediatr Res. 2015;79:425–31.View ArticlePubMedPubMed CentralGoogle Scholar
- DiGiulio DB, Callahan BJ, McMurdie PJ, Costello EK, Lyell DJ, Robaczewska A, et al. Temporal and spatial variation of the human microbiota during pregnancy. Proc Natl Acad Sci. 2015;112:11060–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SSK, McCulle SL, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci. 2011;108:4680–7.View ArticlePubMedGoogle Scholar
- Stewart CJ, Marrs ECL, Nelson A, Lanyon C, Perry JD, Embleton ND, et al. Development of the preterm gut microbiome in twins at risk of necrotising enterocolitis and sepsis. PLoS One. 2013;8:e73465.View ArticlePubMedPubMed CentralGoogle Scholar
- Mai V, Young CM, Ukhanova M, Wang X, Sun Y, Casella G, et al. Fecal microbiota in premature infants prior to necrotizing enterocolitis. PLoS One. 2011;6:e20647.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Stewart CJ, Cummings SP. Gut bacteria and necrotizing enterocolitis: cause or effect? Trends Microbiol. 2015;23:332–3.View ArticlePubMedGoogle Scholar
- Marques TM, Wall R, O’Sullivan O, Fitzgerald GF, Shanahan F, Quigley EM, et al. Dietary trans-10, cis-12-conjugated linoleic acid alters fatty acid metabolism and microbiota composition in mice. Br J Nutr. 2015;113:728–38.View ArticlePubMedGoogle Scholar
- Stewart C, Skeath T, Nelson A, Fernstad S, Marrs E, Perry J, et al. Preterm gut microbiota and metabolome following discharge from intensive care. Sci Rep. 2015;5:17141.View ArticlePubMedPubMed CentralGoogle Scholar
- AlFaleh K, Anabrees J. Probiotics for prevention of necrotizing enterocolitis in preterm infants. Cochrane Database Syst Rev. 2014;4:CD005496.Google Scholar
- Costeloe K, Hardy P, Juszczak E, Wilks M, Millar MR. Bifidobacterium breve BBG-001 in very preterm infants: a randomised controlled phase 3 trial. Lancet. 2016;387:649–60.View ArticlePubMedGoogle Scholar
- Repa A, Thanhaeuser M, Endress D, Weber M, Kreissl A, Binder C, et al. Probiotics (Lactobacillus acidophilus and Bifidobacterium bifidum) prevent NEC in VLBW infants fed breast milk but not formula. Pediatr Res. 2015;77:381–8.View ArticlePubMedGoogle Scholar
- Embleton ND, Berrington JE, McGuire W, Stewart CJ, Cummings SP. Lactoferrin: Antimicrobial activity and therapeutic potential. Semin Fetal Neonatal Med. 2013. doi:10.1016/j.siny.2013.02.001. [Epub ahead of print].
- Butel M, Suau A, Campeotto F. Conditions of bifidobacterial colonization in preterm infants: a prospective analysis. J Pediatr Gastroenterol Nutr. 2007;44:577–82.View ArticlePubMedGoogle Scholar
- Abdulkadir B, Nelson A, Skeath T, Marrs ECL, Perry JD, Cummings SP, et al. Routine use of probiotics in preterm infants: longitudinal impact on the microbiome and metabolome. Neonatology. 2016;109:239–47.View ArticlePubMedGoogle Scholar
- Kuehl FA, Egan RW. Prostaglandins, arachidonic acid, and inflammation. Science. 1980;210:978–84.View ArticlePubMedGoogle Scholar
- Ricciotti E, FitzGerald GA. Prostaglandins and inflammation. Arterioscler Thromb Vasc Biol. 2011;31:986–1000.View ArticlePubMedPubMed CentralGoogle Scholar
- Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS One. 2011;6:e27310.View ArticlePubMedPubMed CentralGoogle Scholar
- Team RC. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014. Available from: http://www.r-project.org/.Google Scholar
- McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.View ArticlePubMedPubMed CentralGoogle Scholar
- Li S, Park Y, Duraisingham S, Strobel FH, Khan N, Soltow QA, et al. Predicting network activity from high throughput metabolomics. PLoS Comput Biol. 2013;9:e1003123.View ArticlePubMedPubMed CentralGoogle Scholar
- Trivedi DK, Iles RK. The application of SIMCA P+ in shotgun metabolomics analysis of ZICⓇHILIC-MS spectra of human urine—experience with the Shimadzu IT-T of and profiling solutions data extraction software. J Chromatogr Sep Tech. 2012;03:145.Google Scholar