Development of the Preterm Gut Microbiome in Twins at Risk of Necrotising Enterocolitis and Sepsis

Abstract

The preterm gut microbiome is a complex dynamic community influenced by genetic and environmental factors and is implicated in the pathogenesis of necrotising enterocolitis (NEC) and sepsis. We aimed to explore the longitudinal development of the gut microbiome in preterm twins to determine how shared environmental and genetic factors may influence temporal changes and compared this to the expressed breast milk (EBM) microbiome. Stool samples (n = 173) from 27 infants (12 twin pairs and 1 triplet set) and EBM (n = 18) from 4 mothers were collected longitudinally. All samples underwent PCR-DGGE (denaturing gradient gel electrophoresis) analysis and a selected subset underwent 454 pyrosequencing. Stool and EBM shared a core microbiome dominated by Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae. The gut microbiome showed greater similarity between siblings compared to unrelated individuals. Pyrosequencing revealed a reduction in diversity and increasing dominance of Escherichia sp. preceding NEC that was not observed in the healthy twin. Antibiotic treatment had a substantial effect on the gut microbiome, reducing Escherichia sp. and increasing other Enterobacteriaceae.

This study demonstrates related preterm twins share similar gut microbiome development, even within the complex environment of neonatal intensive care. This is likely a result of shared genetic and immunomodulatory factors as well as exposure to the same maternal microbiome during birth, skin contact and exposure to EBM. Environmental factors including antibiotic exposure and feeding are additional significant determinants of community structure, regardless of host genetics.

About the Authors:

Christopher J. Stewart

Affiliation: Faculty of Health and Life Sciences, University of Northumbria, Newcastle upon Tyne, United Kingdom

Emma C. L. Marrs

Affiliation: Department of Microbiology, Freeman Hospital, Newcastle upon Tyne, United Kingdom

Andrew Nelson

Affiliation: Faculty of Health and Life Sciences, University of Northumbria, Newcastle upon Tyne, United Kingdom

Clare Lanyon

Affiliation: Faculty of Health and Life Sciences, University of Northumbria, Newcastle upon Tyne, United Kingdom

John D. Perry

Affiliation: Department of Microbiology, Freeman Hospital, Newcastle upon Tyne, United Kingdom

Nicholas D. Embleton

Affiliation: Newcastle Neonatal Service, Royal Victoria Infirmary, Newcastle upon Tyne, United Kingdom

Stephen P. Cummings

* E-mail: [email protected]

Affiliation: Department of Microbiology, Freeman Hospital, Newcastle upon Tyne, United Kingdom

Janet E. Berrington

Affiliation: Newcastle Neonatal Service, Royal Victoria Infirmary, Newcastle upon Tyne, United Kingdom

Introduction

The gut microbiome is crucial to both health (immunomodulation/protection, nutrition and metabolism) and disease (inflammation, diabetes, obesity and allergy) [1]. Due to the complexities of both the microbial community and factors that affect it, exploring individual variables (including diet, medical interventions/exposures and genetic components) is challenging. Studying twins or higher order multiples may provide unique insights. Healthy twins have been shown to develop a comparable gut microbiome after term birth [2], in childhood [3] and into adulthood [4], suggesting genetic or shared environmental factors shape the gut community. However, there is currently a lack of research exploring the temporal changes of the gut microbiome in preterm twins.

Preterm neonates provide a unique cohort to study the gut microbiome due to intensive care practises and the susceptibility of these infants to disease, such as necrotising enterocolitis (NEC) and sepsis, resulting in differential microbiome development compared to term infants [5]. NEC is an inflammatory disease process of the intestine with a pathogenesis strongly linked to immune dysregulation, enteral feeding, and abnormal bacterial colonisation [5]. Changes in the microbial community have been associated with NEC and sepsis in unrelated preterm neonates [6,7]. Studying twins with NEC may help elucidate the role of specific exposures that may be key to reducing incidence of these diseases. Despite previous studies implicating dysbiosis of the gut microbiome in disease development, no single causative agent has been convincingly identified as a biomarker of NEC. A previous study of 20 infants, including 4 sets of discordant twin pairs (one with NEC, one without) confirmed that the communities in the bowels of NEC infants were different to healthy controls [8]. However, sampling was inconsistent and was not longitudinal, thus any temporal changes in microbial community dynamics that may have been associated with the pathophysiology of NEC were not captured.

For term infants vaginal delivery and receipt of maternal breast milk are key factors that facilitate the development of a ‘healthy’ microbiome. Breast milk contains many immunomodulatory factors that support growth and prevent infection including lysozyme, lactoferrin, and oligosaccharides as well as live bacteria which regulate host-microbe interactions [9] and modify infant gut microbiome development [10]. Preterm infants are less likely to experience vaginal birth or only receive breast milk feeds and are more likely to experience many medical interventions that affect the microbiome. To what extent the infant gut microbiome reflects the maternal breast milk microbiome is currently unclear.

In a twin cohort we aimed to explore the longitudinal development of the gut bacterial community after preterm birth by analysis of stool and expressed breast milk (EBM). In addition we focused on the development of NEC in one set of twins with regular longitudinal sampling and where only one infant was diagnosed with NEC.

Materials and Methods

Patients and samples

Preterm infants (<32 weeks gestation) from multiple births cared for in the Royal Victoria Infirmary (Newcastle upon Tyne, UK) had stool collected between April 2010 and February 2012. Stool samples (n = 173) were salvaged directly from the nappy from 12 sets of twins and 1 set of triplets (n = 27) from birth until discharge. EBM (n = 18) from 4 mothers was also analysed as EBM samples were only salvaged toward to latter stages of the study and often insufficient volume was acquired. EBM samples were obtained from the remaining liquid in the syringe after continual feed to the neonate through a catheter tube. All samples were frozen at -20°C within 24 hours. This unit has standardised feeding practices, antibiotic and antifungal use. Fluconazole is used as antifungal prophylaxis if <26 weeks/<1kg and in babies developing NEC. Antibiotics, specifically Penicillin and Gentamycin, were administered to all infants following birth and stopped at 48 hours unless proven infection. Clinical information was obtained from patient notes: NEC was confirmed independently by two experienced neonatal clinicians (JEB/NDE): sepsis was defined by positive blood culture, along with antibiotic treatment for a minimum of 5 days and clinical signs suggestive of infection.

DNA Extraction

Nucleic acid extraction of stool was carried out on 100 mg of sample using the PowerLyzer™ PowerSoil® DNA Isolation Kit (MoBio, CA, USA). Nucleic acid extraction of EBM was performed on 1.8 mL of EBM using the Powerfood™ Microbial DNA Isolation Kit (MoBio, CA, USA). Both kits were used in accordance with the manufacturer’s instructions.

PCR-DGGE

PCR amplification targeted the V3 region within the bacterial 16S rRNA gene and was performed with universal primers for total bacterial amplification; V3F-GC (5'-CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGGCCTACGGGAGGCAGCAG-3') and V3R (5'-ATTACCGCGGCTGCTGG-3') as described by Muyzer et al. [11], according to the methods of Baxter and Cummings [12].

DGGE (denaturing gradient gel electrophoresis) analysis was performed using the D-Code DGGE system (Bio-Rad, CA, USA) and selected DGGE bands (n=17) underwent excision and sequencing as previously described [13] (Table S1). A ladder of amplicons derived from pure culture organisms was loaded on each DGGE gel for multiple gel alignment. The analysis of the DGGE banding patterns employed TotalLab Phoretix 1D and 1D Pro software (Nonlinear Dynamics, Newcastle upon Tyne, UK). Sequencing was carried out on an ABI 3730XL sequencer (Eurofins MWG Operon, London, UK) and BLASTn was used to determine the closest sequence homologies. DGGE sequence data was deposited to the National Center for Biotechnology Information GenBank under the accession numbers KC686617-KC686622. Analysis of sequence data was carried out at the genus level where possible. As the sequences generated were based on a ~200 bp fragment of the 16S rRNA gene, robust species level classification is not feasible.

454 Pyrosequencing

A subset of key samples (Stool = 36EBM = 5) was selected for metagenomic sequencing: primarily EBM and stool from the triplets (patients 145/147/148) and a discordant twin pair where one developed NEC (patients 139/140). Although other twin sets in the study were discordant for NEC and/or sepsis, patients 139/140 provided the most robust sampling pre, during, and post disease diagnosis. 454 pyrosequencing was carried out using the bifidobacteria-optimised primer set 357F (5'-CTCCTACGGGAGGCAGCAGAN-3') and 926Rb (5'-CCGTCAATTYMTTTRAGT-3') [14]. The DNA sequencing was performed using the 454 Genome Sequencer FLX Titanium platform (Roche, IN, USA) by the Research and Testing Laboratory (RTL, TX, USA) using previously described methods [15].

The raw sequencing reads were quality filtered in QIIME 1.6.0 [16] using the split-library.py script with the following criteria: 1) exact matches to bar code tags; 2) no ambiguous bases; 3) maximum of 5 mismatches to primer sequence; 4) read-lengths between 200–700 base pairs (bp); 5) average quality score of >25 in a sliding window of 50 bp. Remaining high quality sequences were clustered into operational taxonomic units (OTUs) at 97% similarity using UCLUST [17]. Representative sequences for each OTU were aligned using PyNAST [18] and taxonomic identities were assigned using RDP-classifier (version 2.2) [19] with 50% as confidence value threshold. Detection of potentially chimeric sequences was performed using ChimeraSlayer [20] and chimeric sequences were removed from downstream analysis prior to tree building using FastTree [21]. Sequences were deposited in MG-RAST under the accession numbers 4516545.3-4516585.3

Statistical analysis

The DGGE profiles obtained from both stool and EBM samples were analysed using multivariate partial least squares discriminant analysis (PLS-DA) using SIMCA 13.0 (Umetrics, Stockholm, Sweden) [22]. To check that data was adhering to multivariate normalities, Hotelling’s T2 tolerance limits were calculated and set at 0.95. Shannon diversity indices (H'), a diversity index giving a measure of both species richness and the evenness of their abundance, was calculated using the PAST package [23].

The sequence reads generated using 454 pyrosequencing were analysed using weighted UniFrac [24] and visualised using principal coordinate analysis (PCoA). Rarefaction curves were produced to 5000 sequences for all samples (Figure S1). Martin’s P-test of significance [25] was applied to determine if significant differences were occurring in the profiles between and within patients. For all statistical analysis, a probability (P) value of <0.05 was deemed to be statistically significant.

Ethics Statement

Initial collection was part of routine service and all samples were collected during the course of normal treatment. Ethical approval was obtained from the County Durham and Tees Valley Research Ethics Committee to include molecular techniques in August 2010. For all infant stool samples, prospective parental informed consent was documented at the point of donation to have the samples stored for research purpose. For EBM, informed consent was documented at the point of donation to have the samples stored for research purposes from March 2011.

Results

Patients

27 patients (12 twin pairs 1 triplet set) contributed to the study: five developed definite NEC and five developed sepsis (two of whom also had NEC). Demographic information is summarised in Table 1 and a detailed version based on single samples can be found in the supplementary information (Table S2). No infant received supplemental probiotics.

Patient Delivery Mode Gestation (weeks) Birth Weight (g) Sex Chorionicity EBM ever NEC NEC Diagnsis (DOLa) BCb+ BCb (DOLa) Organism Number of samples collected

22 Caesarean 27 870 M DDc Y N N 2

23 Caesarean 27 790 F DDc Y N N 6

28 Caesarean 28 1250 M MMd Y N N 6

29 Caesarean 28 1180 M MMd Y Y 17 N 1

39 Caesarean 25 780 M DDc N N Y 15 CoNSg 8

41 Caesarean 25 820 M DDc N N Y 15 S. aureus 8

46 Caesarean 26 830 M DDc Y N N 6

47 Vaginal 26 760 F DDc Y Y 45 Y 03/08 CoNSG / Pseudomonas sp. 5

51 Caesarean 27 1060 M DDc Y Y 16 Y 17/56 CoNSg / CoNSg 3

55 Caesarean 27 1100 M DDc Y N N 2

68 Vaginal 26 760 M DDc Y N N 4

70 Vaginal 26 860 M DDc Y N Y 40 CoNSg 2

92 Caesarean 25 740 M DDc Y N Y 8 K. pneumoniae 1

93 Caesarean 25 670 M DDc Y N N 5

100 Caesarean 27 1050 M MDd Y N N 6

101 Caesarean 27 910 M MDd Y N N 5

112 Vaginal 25 700 F DDc Y N N 12

113 Vaginal 25 680 M DDc Y N Y 57 CoNSg 7

135 Caesarean 29 910 F MMe (TTTS)c N N Y 54 S. aureus + K. pneumoniae 8

136 Caesarean 29 1275 F MMe (TTTS)c N N N 3

139 Caesarean 30 1470 M DDc Y Y 28 N 25

140 Caesarean 30 1455 F DDc Y N N 20

145 Caesarean 31 990 M MMe Y N N 10

148 Caesarean 31 1455 M MMe Y N N 4

147 Caesarean 31 1540 M DDc Y N N 5

151 Vaginal 27 1020 F DDc Y N N 4

154 Vaginal 27 1060 M DDc Y Y 21 N 5

Table 1. Summary of patient demographics.

a Day of life.

b Blood culture.

c Dichorionic diamniotic.

d Monochorionic diamniotic.

e Monochorionic monoamniotic.

f TTTS – Twin to Twin Transfusion Syndrome.

g Coagulase negative Staphylococcus.

CSV

Download CSV

Stool profiles

PLS-DA of all samples from all patients based on DGGE data demonstrated twin pairs had comparable profiles which were distinct from unrelated individuals (Figure 1A). Samples which fell outside the ellipse (indicating Hotellings T2 range, at 95% confidence) belonged to the triplets or were associated with late onset infection. Specifically, this comprised sample 139t collected one day post NEC diagnosis as well as samples 39c and 41b collected three days prior and on the day of sepsis diagnosis, respectively. Based on DGGE analysis of the cohort, a subset of longitudinal samples from the triplets (represented by red squares) and twin pair discordant for NEC (represented by orange circles) which showed significant variation in community development were selected for pyrosequencing. The pyrosequencing data was analysed at the genus level using weighted UniFrac. In accordance with the DGGE data, samples grouped with their related twin, showing high intra-sample similarities in the development of the gut microbiome (Figure 1B). This is reflected in the bar plots which show each set of twins and the triplets developed a distinct gut microbiome (Figure S2). Proteobacteria and Firmicutes dominated samples in both DGGE and pyrosequencing. From pyrosequencing, 5 OTUs were found in the core microbiome (present in over 85% of samples) in stools from the families Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae.

[Figure omitted. See PDF.]

Figure 1. Community profiles of gut (stool) microbiome from preterm multiples.

Subjects are symbolised based on related multiples. a) Partial least squares discriminant analysis (PLS-DA) based on DGGE data of all stool samples. The ellipse indicates Hotellings T2 range, at 95% confidence. Selected sample labels removed for clarity. b) Weighted UniFrac PCoA based on pyrosequencing data at genus level, generated from a subset of stool samples.

https://doi.org/10.1371/journal.pone.0073465.g001

Comparison of breast milk with respective stool samples

The DGGE profiles of all EBM samples were compared with the respective stool sample. EBM samples showed a relative lack of diversity making them cluster near the origin: despite this PLS-DA did reveal EBM samples clustered with the stool samples of the respective set of multiples (Figure 2A). There are cases where different EBM samples from the same mother cluster separately showing the EBM microbiome was not stable (e.g. BM148e / BM147b / BM145a / BM145i or BM139b / BM139f / BM139q).

[Figure omitted. See PDF.]

Figure 2. Comparison of breast milk with respective stool profiles.

Subjects are symbolised based on related multiples. a) Partial least squares discriminant analysis (PLS-DA) based on DGGE data of all breast milk (EBM) samples matched to respective stools. The ellipse indicates Hotellings T2 range, at 95% confidence. b) Weighted UniFrac PCoA based on pyrosequencing data at genus level, generated from a subset of EBM and respective stool samples.

https://doi.org/10.1371/journal.pone.0073465.g002

EBM from the triplet set and twin set (patients 139/140) also underwent pyrosequencing. EBM samples from each mother showed high intra-sample similarity (Figure 2B). For twins 139/140 the EBM profiles (BM139/BM140) clustered closely with the respective stool (139s/140s). This clustering was less robust in the triplet set although triplet EBM (BM145/BM147/BM148) was still more similar to triplet stool (145a/147b/148b) than stool of other infants. The connection of EBM samples with stool can be visualised in the bar plot (Figure S2). Three OTUs were found in the core microbiome in EBM and, like in the stool core microbiome, were from families Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae.

Bacterial community development in twin set 139/140

A reduction in diversity was observed in patient 139 at least 5 days before NEC was apparent clinically (Figure 3A), which was not shared by the twin. A less pronounced reduction in diversity did occur in the control twin (140) earlier in development (day 18) that coincided with antibiotic administration and diversity was quickly re-established when antibiotics ceased. A similar trend was observed when the H' of the pyrosequencing data was performed (data not shown).

[Figure omitted. See PDF.]

Figure 3. Development of gut microbiome in twin pair 139/140 mapped to life events.

P – Penicillin, G – Gentamicin, F- Fluctoxacillin, A - Amoxicillin, M – Metronidazole, s –Start of antibiotics, e – End of antibiotics, 72hr – full enteral feed (at least 150 ml/kg/day) sustained for 72 hours. a) Shannon Diversity indices (H') of twin pair based on DGGE data b) Turnover of the most prevalent bacterial OTUs throughout the first 36 days of life in twin 139 where antibiotics were prescribed for NEC. c) Turnover of most prevalent bacterial OTUs throughout the first 34 days of life in twin 140 where antibiotics were prescribed due to pyrexia (fever).

https://doi.org/10.1371/journal.pone.0073465.g003

Pyrosequencing of matched twin samples (full profiles in the bar plot; Figure S2) allowed the 12 most abundant OTUs to be incorporated into the analysis (Figure 3). These results are in accordance with the DGGE data and allow greater detail about specific contributors to the overall diversity to be observed. While community structure in the twins was initially comparable, twin 139 showed reducing diversity with an increasing abundance of Escherichia sp., before NEC diagnosis on day 28 (Figure 3B). Conversely, in the sibling there is an increase in diversity and Escherichia sp. was present in much lower abundance (Figure 3C). After antibiotic receipt (day 16 in twin 140 and day 28 in twin 139) both twins demonstrate reduced Escherichia sp. abundance and an increased abundance of other Enterobacteriaceae, rapidly reversing in twin 140 on antibiotic cessation. In twin 139 NEC and subsequent antibiotic treatment significantly (P=0.028) altered the bacterial community in comparison to its sibling, with a new notable bloom in Klebsiella sp. and a smaller increase in Parascardovia sp. (family – Bifidobacteriaceae).

Discussion

In this study, we hypothesised that the bacterial community in related twins and a set of triplets would be comparable and reflect maternal EBM community. In cases of twin pairs discordant for NEC we hypothesized that they would show differences in microbial community development prior to NEC onset.

We have demonstrated in a preterm population with multiple clinical exposures that the development of the gut microbiome is more similar between genetically related individuals than between other preterm infants. However, due to similar environmental exposures encompassed by related individuals, this may not be a direct result of host genetics. Interestingly, community structure was similar for all triplets even though two (145/148) were monochorionic monoamniotic and the other (147) was dichorionic. This suggests that shared factors (genetic or environmental) are important in determining the gut microbiome even in an environment with many complex variable factors that also affect community development [2,26].

Proteobacteria and Firmicutes dominated the gut microbiome over the initial weeks of life as previously reported [27]. Stool and EBM shared a core microbiome of the families Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae. Others describing the EBM microbiome noted the presence of Streptococcaceae which in was present in low abundance in our maternal cohort [28,29]. We also detected a low abundance of bifidobacteria and lactobacilli in stool and EBM despite their reported prevalence by others [10,28,30]. This may be attributable to differences among subjects, organisms residing in the unit, and detection by differing techniques [29], but occurred despite the use of optimised universal primers in pyrosequencing designed to facilitate the detection of Bifidobacteria [14]. Notably, the cohort in this study is significantly preterm which can result in delayed colonisation with bifidobacteria; geographical or demographic differences may also be attributable to the low prevalence of this organism [2,7]. The milk microbiome was not stable throughout lactation [31] and EBM appeared to be an on-going source of new bacteria contributing to the dynamic nature of the gut microbiome. Furthermore, genera which typically reside on adult skin including Staphylococcus, Corynebacterium, and Propionibacterium were found in high abundance in the gut microbiome suggesting skin contact may be an important source of bacterial acquisition, even within the nursery environment [32].

The importance of the gut microbiome in disease is increasingly recognised, despite a lack of consistent causative agent between studies. In a focused temporal exploration from a twin pair discordant for NEC, we showed clear changes attributable to antibiotic exposure and NEC development, with effects on the dominance of Escherichia sp. and the abundance of other Enterobacteriaceae sp. [33]. The significantly different community observed in sample 139t is likely attributable to a temporary bloom in Klebsiella sp. following NEC diagnosis, which was reduced in the subsequent sample following broad spectrum antibiotic administration.

A twin study methodology has been utilised in other inflammatory bowel diseases such as ulcerative colitis and Crohn’s disease, but there are few comparable studies in NEC. Interestingly, analogous to NEC, dysbiosis is a major factor in the pathogenesis of these diseases, consistent with a lack of a single causative agent [34]. Specifically, a decreased diversity in the gut microbiome of ulcerative colitis and Crohn’s disease patients compared to healthy controls has been noted [8,35], and multiple studies report increased abundances of Proteobacteria, particularly Escherichia sp. [36,37]. This increase in Proteobacteria is an emerging theme in NEC development [6,13], but is probably one of several factors needed for NEC development. Escherichia spp. are reported pathogens [38] and the association of this genus with inflammatory mediated disease warrants further investigation.

Our use of molecular approaches for community profiling circumvents the known limitations of culturing human gut species [39]. We utilised cost effective DGGE to educate sample selection for 454 pyrosequencing. Studies on preterm infants, especially of multiple births, are difficult due to the exclusivity of the cohort and thus only 27 patients could be included. With the exception of the twin set 139/140, twins discordant for disease often lacked informative longitudinal samples; in large due to the withdrawal of feeds resulting in reduced excrement. Although the number and timing of samples collected from each set of twins was generally comparable, varying numbers of samples were collected between twin sets which may bias some analysis. Overall, the data generated using DGGE and pyrosequencing were in agreement, perhaps due to primer sets encompassing the V3 hypervariable region of the 16S rRNA genome in both analyses. Pyrosequencing allows larger fragments to be amplified allowing the use of bifidobacteria-optimised universal primers. Despite this, bifidobacteria was not found to be a prevalent genus [2,7] which requires further study due to the potential use of this taxa as probiotic in therapeutic intervention [40].

In summary, this study represents a unique temporal analysis of the gut microbiome in preterm twins, cared for within the complex environment of neonatal intensive care. Although twins discordant for late onset infection showed differences in gut microbiome development, overall, related infants harboured bacterial communities more similar to each other than nonrelated infants. As well as shared genetic and immunomodulatory factors, this is likely a result of exposure to the same maternal microbiota during birth, skin contact and exposure to EBM. We have also shown that other environmental factors, particularly antibiotic exposure, have additional significant effects on the gut microbiome in genetically related infants. Antibiotics have been shown to alter the gut microbiome in term [2] and adult populations [41] and the exact role of individual antibiotics in altering the preterm gut microbiome warrants further investigation. We have further noted potential concurrence between community changes associated with other inflammatory mediated diseases, such as ulcerative colitis and Crohn’s disease, and those increasingly reported to occur in NEC.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Rarefactions curves produced in QIIME to 5000 sequences comparing all samples (stool and expressed breast milk).

https://doi.org/10.1371/journal.pone.0073465.s001

(TIF)

Figure S2.

Order level bar plot of all samples (stool and expressed breast milk) which underwent 454 pyrosequencing.

https://doi.org/10.1371/journal.pone.0073465.s002

(TIF)

Table S1.

DGGE band identities and prevalence.

https://doi.org/10.1371/journal.pone.0073465.s003

(XLSX)

Table S2.

Complete table of patient demographics.

https://doi.org/10.1371/journal.pone.0073465.s004

(XLSX)

Author Contributions

Conceived and designed the experiments: SPC CJS JEB. Performed the experiments: CJS AN ECLM. Analyzed the data: CJS CL SPC. Contributed reagents/materials/analysis tools: JDP NDE JEB. Wrote the manuscript: CJS SPC JEB NDE JDP.

Citation: Stewart CJ, Marrs ECL, Nelson A, Lanyon C, Perry JD, Embleton ND, et al. (2013) Development of the Preterm Gut Microbiome in Twins at Risk of Necrotising Enterocolitis and Sepsis. PLoS ONE 8(8): e73465. https://doi.org/10.1371/journal.pone.0073465

References

1. Sekirov I, Russell SL, Antunes LC, Finlay BB (2010) Gut microbiota in health and disease. Physiol Rev 90: 859–904. doi:10.1152/physrev.00045.2009. PubMed: 20664075.

2. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO (2007) Development of the human infant intestinal microbiota. PLOS Biol 5: 1556–1573. PubMed: 17594176.

3. Stewart JA, Chadwick VS, Murray A (2005) Investigations into the influence of host genetics on the predominant eubacteria in the faecal microflora of children. J Med Microbiol 54: 1239–1242. doi:10.1099/jmm.0.46189-0. PubMed: 16278440.

4. Zoetendal E, Akkermans A, Akkermans-van Vliet W, Visser J, Vos W (2001) The host genotype affects the bacterial community in the human gastronintestinal tract. Microb Ecol Health Dis 13: 129–134. doi:10.1080/089106001750462669.

5. Berrington J, Stewart C, Embleton N, Cummings S (2012) Gut microbiota in preterm infants: assessment and relevance to health and disease. Arch Dis Child Fetal Neonatal Ed. Available: http://fn.bmj.com/content/early/2012/09/24/archdischild-2012-302134.short. Accessed 2013 February 27.

6. Mai V, Young CM, Ukhanova M, Wang X, Sun Y et al. (2011) Fecal microbiota in premature infants prior to necrotizing enterocolitis. PLOS ONE 6: e20647. doi:10.1371/journal.pone.0020647. PubMed: 21674011.

7. Stewart CJ, Marrs EC, Magorrian S, Nelson A, Lanyon C et al. (2012) The preterm gut microbiota: changes associated with necrotizing enterocolitis and infection. Acta Paediatr 101: 1121–1127. doi:10.1111/j.1651-2227.2012.02801.x. PubMed: 22845166.

8. Wang Y, Hoenig JD, Malin KJ, Qamar S, Petrof EO et al. (2009) 16S rRNA gene-based analysis of fecal microbiota from preterm infants with and without necrotizing enterocolitis. ISME J 3: 944–954. doi:10.1038/ismej.2009.37. PubMed: 19369970.

9. Martín R, Langa S, Reviriego C, Jimínez E, Marín ML et al. (2003) Human milk is a source of lactic acid bacteria for the infant gut. J Pediatr 143: 754–758. Available: http://www.ncbi.nlm.nih.gov/pubmed/14657823. doi:10.1016/j.jpeds.2003.09.028. PubMed: 14657823.

10. Martín R, Jiménez E, Heilig H, Fernández L, Marín ML et al. (2009) Isolation of bifidobacteria from breast milk and assessment of the bifidobacterial population by PCR-denaturing gradient gel electrophoresis and quantitative real-time PCR. Appl Environ Microbiol 75: 965–969. doi:10.1128/AEM.02063-08. PubMed: 19088308.

11. Muyzer G, De Waal EC, Uitterlinden AG (1993) Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl Environ Microbiol 59: 695–700. PubMed: 7683183.

12. Baxter J, Cummings SP (2006) The current and future applications of microorganism in the bioremediation of cyanide contamination. Antonie Van Leeuwenhoek 90: 1–17. doi:10.1007/s10482-006-9057-y. PubMed: 16683094.

13. Stewart CJ, Nelson A, Scribbins D, Marrs EC, Lanyon C et al. (2013) Bacterial and fungal viability in the preterm gut: NEC and sepsis. Arch Dis Child Fetal Neonatal Ed, 98: F298–303. PubMed: 23426613.

14. Sim K, Cox MJ, Wopereis H, Martin R, Knol J et al. (2012) Improved Detection of Bifidobacteria with Optimised 16S rRNA-Gene Based Pyrosequencing. PLOS ONE 7: e32543. doi:10.1371/journal.pone.0032543. PubMed: 22470420.

15. Dowd SE, Sun Y, Secor PR, Rhoads DD (2008) Survey of bacterial diversity in chronic wounds using pyrosequencing, DGGE, and full ribosome shotgun sequencing. BMC Microbiol 8: 43. doi:10.1186/1471-2180-8-43. PubMed: 18325110.

16. Caporaso JG, Kuczynski J, Stombaugh J (2010) QIIME allows analysis of high-throughput community sequencing data. Nature 7: 335–336. PubMed: 20383131.

17. Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26: 2460–2461. doi:10.1093/bioinformatics/btq461. PubMed: 20709691.

18. Caporaso JG, Bittinger K (2010) PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26: 266–267. doi:10.1093/bioinformatics/btp636. PubMed: 19914921.

19. Wang Q, Garrity GM, Tiedje JM, Cole JR (2007) Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 73: 5261–5267. doi:10.1128/AEM.00062-07. PubMed: 17586664.

20. Haas BJ, Gevers D, Earl AM (2011) Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res 21: 494–504. doi:10.1101/gr.112730.110. PubMed: 21212162.

21. Price MN, Dehal PS, Arkin AP (2010) FastTree 2-approximately maximum-likelihood trees for large alignments. PLOS ONE 5: e9490. doi:10.1371/journal.pone.0009490. PubMed: 20224823.

22. Eriksson L, Johansson E, Kettaneh-Wold N (2006) Multi- and Megavariate Data Analysis, Part 2, Advanced Applications and Method Extensions. MKS Umetrics AB.

23. Hammer Ø, Harper D (2001) PAST: Paleontological statistics software package for education and data analysis. Palaeo Electronica 4: 1–9.

24. Hamady M, Lozupone C, Knight R (2009) Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J 4: 17–27. PubMed: 19710709.

25. Martin AP (2002) Phylogenetic approaches for describing and comparing the diversity of microbial communities. Appl Environ Microbiol 68: 3673-3682. doi:10.1128/AEM.68.8.3673-3682.2002. PubMed: 12147459.

26. Hamady M, Lozupone C, Knight R (2009) Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J 4: 17–27. PubMed: 19710709.

27. Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD et al. (2005) Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A 102: 11070–11075. doi:10.1073/pnas.0504978102. PubMed: 16033867.

28. Mshvildadze M, Neu J, Shuster J, Theriaque D, Li N et al. (2010) Intestinal microbial ecology in premature infants assessed with non-culture-based techniques. J Pediatr 156: 20–25. doi:10.1016/j.jpeds.2009.06.063. PubMed: 19783002.

29. Collado MC, Delgado S, Maldonado A, Rodríguez JM (2009) Assessment of the bacterial diversity of breast milk of healthy women by quantitative real-time PCR. Lett Appl Microbiol 48: 523–528. doi:10.1111/j.1472-765X.2009.02567.x. PubMed: 19228290.

30. Hunt KM, Foster JA, Forney LJ, Schütte UME, Beck DL et al. (2011) Characterization of the diversity and temporal stability of bacterial communities in human milk. PLOS ONE 6: e21313. doi:10.1371/journal.pone.0021313. PubMed: 21695057.

31. Penders J, Thijs C, Vink C, Stelma FF, Snijders B et al. (2006) Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics 118: 511–521. doi:10.1542/peds.2005-2824. PubMed: 16882802.

32. Cogen AL, Nizet V, Gallo RL (2008) Skin microbiota: a source of disease or defence? Br J Dermatol 158: 442–455. doi:10.1111/j.1365-2133.2008.08437.x. PubMed: 18275522.

33. Madan JC, Salari RC, Saxena D, Davidson L, O’Toole GA et al. (2012) Gut microbial colonisation in premature neonates predicts neonatal sepsis. Arch Dis Child Fetal Neonatal Ed 97: F456–F462. PubMed: 22562869.

34. Lepage P, Häsler R, Spehlmann ME, Rehman A, Zvirbliene A et al. (2011) Twin study indicates loss of interaction between microbiota and mucosa of patients with ulcerative colitis. Gastroenterology 141: 227–236. doi:10.1053/j.gastro.2011.04.011. PubMed: 21621540.

35. Gophna U, Sommerfeld K, Gophna S, Doolittle WF, Veldhuyzen van Zanten SJO (2006) Differences between tissue-associated intestinal microfloras of patients with Crohn’s disease and ulcerative colitis. J Clin Microbiol 44: 4136–4141. doi:10.1128/JCM.01004-06. PubMed: 16988016.

36. Darfeuille-Michaud A, Boudeau J, Bulois P, Neut C, Glasser AL et al. (2004) High prevalence of adherent-invasive Escherichia coli associated with ileal mucosa in Crohn’s disease. Gastroenterology 127: 412–421. doi:10.1053/j.gastro.2004.04.061. PubMed: 15300573.

37. Sasaki M, Sitaraman SV, Babbin BA, Gerner-Smidt P, Ribot EM et al. (2007) Invasive Escherichia coli are a feature of Crohn’s disease. Lab Invest 87: 1042–1054. doi:10.1038/labinvest.3700661. PubMed: 17660846.

38. Kaper JB, Nataro JP, Mobley HL (2004) Pathogenic Escherichia coli. Nat Rev Microbiol 2: 123–140. doi:10.1038/nrmicro818. PubMed: 15040260.

39. Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L et al. (2005) Diversity of the human intestinal microbial flora. Science 308: 1635–1638. doi:10.1126/science.1110591. PubMed: 15831718.

40. Embleton ND, Yates R (2008) Probiotics and other preventative strategies for necrotising enterocolitis. Semin Fetal Neonatal Med 13: 35–43. doi:10.1016/j.siny.2007.09.003. PubMed: 17974513.

41. Simões CD, Maukonen J (2013) Habitual Dietary Intake Is Associated with Stool Microbiota Composition in Monozygotic Twins. J Nutr, 143: 417–23. doi:10.3945/jn.112.166322. PubMed: 23343669. PubMed: 23343669.