Intestinal microbiota of preterm infants differ over time and between hospitals
© Taft et al.; licensee BioMed Central Ltd. 2014
Received: 22 April 2014
Accepted: 22 August 2014
Published: 1 October 2014
Intestinal microbiota are implicated in risk of necrotizing enterocolitis (NEC) and sepsis, major diseases of preterm infants in neonatal intensive care units (NICUs). Rates of these diseases vary over time and between NICUs, but time and NICU comparisons of the intestinal microbiota of preterm infants are lacking.
We included 66 singleton infants <29 weeks gestational age with stool samples collected between postnatal days 3 to 21 of life who survived free of NEC and sepsis. Infants were enrolled during 2010 and 2011. Twenty-six infants were enrolled at hospital 1 in Cincinnati, OH, and 40 infants were enrolled at hospital 2 in Birmingham, AL. Samples collected from days 3–9 (“week 1”) and days 10–16 (“week 2”) were compared between years and hospitals. Microbial succession was compared between hospitals in 28 infants with samples from the first 3 weeks of life. DNA extracted from stool was used to sequence the 16S rRNA gene by Illumina MiSeq using universal primers. Resulting operational taxonomic unit tables were analyzed for differences between years and hospitals using linear discriminant analysis effect size algorithm (LEfSe; significance, p < 0.05).
Significant variation was observed in infant microbiota by year and hospital. Among hospital 1 infants, week 1 samples had more phylum Firmicutes (class Bacilli, families Clostridiaceae and Enterococcaceae) in 2010 and more phylum Proteobacteria (family Enterobacteriaceae) in 2011; week 2 samples did not significantly vary over time. However, among hospital 2 infants, the week 1 shift was nearly opposite, with more Proteobacteria (Enterobacteriaceae) in 2010 and more Firmicutes (Bacilli) in 2011; week 2 samples exhibited the same pattern. Regression analysis of clinical covariates found that antibiotic use had an important influence but did not explain these observed shifts in microbiota over time and between hospitals. Microbial succession also differed by hospital, with greater change in microbiota in hospital 1 than hospital 2 infants (p < 0.01, Jaccard distance).
Colonizing microbiota differ over time and between NICUs in ways that could be relevant to disease. Multi-site, longitudinal studies are needed to reliably define the impact of intestinal microbiota on adverse outcomes of preterm infants.
KeywordsInfants Premature Microbiome Geo-temporal analysis Microbial succession
John Snow’s investigation of cholera in 1854 pioneered the understanding of disease risk in relation to not only individual characteristics (person) but also geographic patterns (place) and temporal patterns (time). Since then, epidemiologic investigations have advanced our understanding of the carriage frequency, pathogenicity, and transmission dynamics of enteric organisms in relation to person, place, and time. With the advent of high-throughput sequencing methods, research on the human intestinal microbiome has led to important insights regarding the impact of personal characteristics such as age, gender, diet, and disease state on the colonizing microbial communities. By comparison, few studies have investigated geographic and temporal influences on human microbial communities. Understanding temporal and geographic influences on the colonizing human microbiota is a critical step towards understanding the extent to which studies undertaken in a population defined by time and location can be generalized to the same population at a later time or populations at other locations.
Geographic differences in the human intestinal microbiota have been more extensively studied than temporal differences within populations. For example, Bacteroidetes were shown to dominate the microbiomes of African children, while Firmicutes and Proteobacteria were enriched in the microbiota of European children. Similarly, Lee et al. reported differences in the microbiota of Korean and US twins, particularly in taxa belonging to Firmicutes and Bacteroidetes. In other comparative studies, the microbiota of built environments can differ significantly; such differences could be attributed to the broader local environment or to the inhabitants of the built environment. A study of the colonization patterns of preterm infants in three distinct Florida neonatal intensive care units (NICUs) demonstrated significant differences in the relative abundance of the phyla Bacteroidetes, Proteobacteria, and Firmicutes in the inhabitants of these NICUs. Variations in infant microbiota can lead to differences in environmental exposures and vice versa, as bacteria in the NICU environment are also found to colonize the preterm infant gut.
It is not known whether colonization patterns of intestinal microbial communities shift within human populations over calendar time, but temporal patterns are well documented with many enteric pathogens. Cholera outbreaks have one or two annual peaks, in spring and fall. Helicobacter pylori also shows peak colonization in children in spring and fall[8, 9]. Although it is unclear if commensals also exhibit such temporal patterns, one study found that the microflora of buildings varied over a single year of sampling, and in one of the studied buildings, there was a shift from Actinobacteria to Proteobacteria during the year.
To address the question of temporal and geographic differences in the intestinal microbiota, we conducted a prospective study of the intestinal microbiota of preterm infants at two large, geographically separated level III NICUs. We used 2 years of data from this cohort study to test the hypothesis that initial microbiota colonization of disease-free survivors changes over time and differs between birth hospitals.
Characteristics of infants with samples included in analysis, by week
Hospital 1, week 1
Hospital 2, week 1
Hospital 1, week 2
Hospital 2, week 2
n = 9
n = 13
n = 12
n = 17
n = 11
n = 13
n = 15
n = 21
Birthweight, mean ± sd (g)
963 ± 216
965 ± 261
938 ± 234
953 ± 223
985 ± 193
995 ± 261
914 ± 216
962 ± 207
Birth length, mean ± sd (cm)
35.8 ± 3.0
35.2 ± 2.9
34.0 ± 3.1
34.9 ± 2.9
35.8 ± 2.7
35.5 ± 2.9
33.7 ± 3.5
36.3 ± 6.9
GA at birth, weeks (25th–75th percentile)
Cesarean section (%)
Maternal abx at delivery (%)
Maternal age ± sd (years)a
27.4 ± 1.6
28.0 ± 1.4
22.5 ± 4.9
x27.8 ± 6.1
28.5 ± 6.7
29.1 ± 5.6
21.9 ± 5.1
27.7 ± 5.7
Maternal parity (25th–75th percentile)
Day of life of sample collection, median (25th–75th percentile)b
Days on abx in the first 7 days ± sd
4.0 ± 2.6
3.1 ± 2.1
4.7 ± 2.0
3.8 ± 1.7
Days on abx in the first 14 days ± sd
4.4 ± 3.1
3.6 ± 3.8
7.6 ± 4.4
5.4 ± 3.4
Days of life of first enteral feed, mean ± sd
4.6 ± 3.7
2.5 ± 0.66
3 ± 1.3
2.8 ± 1.2
3.6 ± 1.5
2.7 ± 0.75
3.6 ± 2.4
3 ± 1.6
Median percent formula feeds in the first 7 days (25th–75th percentile)
0% (0%, 0%)
0% (0%, 0%)
0% (0%, 35.7%)
0% (0%, 25%)
Median percent formula feeds in the first 14 days (25th–75th percentile)
0% (0%, 0%)
0% (0%, 0%)
8.3% (0%, 31.7%)
6.7% (0%, 18.2%)
In week 2 samples, comparing infants within hospital 1 by year revealed that colonization patterns were relatively stable between 2010 and 2011 demonstrating only OTU level differences (Figure 2C). In hospital 2 infants, colonization patterns changed between years with decreased Proteobacteria and increased Firmicutes in 2011 (Figure 2D). There were also differences at lower phylogenetic levels in hospital 2, with infants in 2010 enriched in Staphylococcaceae and infants in 2011 enriched in Enterobacteriaceae. Week 2 samples were different between the two hospitals, with hospital 2 enriched in Proteobacteria and hospital 1 enriched in Firmicutes. In 2011, the week 2 infant colonization patterns were similar between the hospitals, consistent with the observation of week 1 colonization patterns.
Results of GEE models of Proteobacteria relative abundance for weeks 1 and 2 sample within each hospital NICU
Year 2011 vs. 2010
Maternal antibiotic use
Infant gestational age
Year 2011 vs. 2010
Maternal age at delivery
Infant gestational age
Models were hospital specific. In both hospitals, models included the following covariates: birth year, maternal antibiotic use, infant antibiotic use in the first 14 days of life, infant gestational age, maternal age at delivery, delivery mode, and day of sample collection. Percentage of feeds that were formula was not included in the hospital 1 model because infants in that hospital received no formula prior to day of life 14. In hospital 1, whether or not a sample was stored in thioglycollate was included to account for any storage protocol differences, as all week 1 samples and 19 (80%) of 24 samples in week 2 were stored in thioglycollate. In hospital 2, we also analyzed the percentage of feeds that were formula; we did not analyze that variable in hospital 1, where all infants were human milk fed. Backwards elimination was used to determine the final model. Any variable with a p value greater than 0.10 that did not change the coefficient of the year variable by more than 10% was eliminated starting with the covariate with the highest p value. Both models at both hospitals found that birth year remained significant, but the other covariates dropped out of most models, as shown in Table 2.
Comparing the Jaccard index calculated from OTU tables confirmed this difference; the Jaccard index between week 1 and week 2, week 1 and week 3, and week 2 and week 3 samples was significantly lower for hospital 2 infants (Kruskal-Wallis, p < 0.01 for each of the three time points). A lower value for the Jaccard index indicates a greater degree of similarity between samples, indicating that infants at hospital 2 experienced less change in microbial colonization from week to week than did infants at hospital 1. As change in microbiota between weeks could be affected by duration of the interval of measurement, we examined the days of collection and time interval between samples for the 28 infants in this analysis. The mean day at the time of sample collection by hospital was the same for weeks 1 and 3 but differed in week 2 (week 1—6.4 for hospital 1 and 6.2 for hospital 2, week 2—12.9 for hospital 1 and 11.1 for hospital 2 [p < 0.05], and week 3—20.2 for hospital 1 and 20.4 for hospital 2). The mean number of days between sample collections was significantly different between proximal weeks but the same for week 1 to week 3 (the days between week 1 and week 2 samples were 6.5 days for hospital 1 and 4.9 days for hospital 2 [p < 0.05], between week 2 to week 3 samples were 7.7 days for hospital 1 and 9.2 days for hospital 2 [p < 0.05], and between week 1 to week 3 samples were 14.2 days for hospital 1 and 14.1 days for hospital 2). Since a greater rate of change in microbiota was observed between weeks in hospital 1 compared to hospital 2 for each time interval regardless of whether the time interval was shorter, longer, or the same, we conclude that hospital 1 infants experienced a more rapid turnover in bacteria and that this was not confounded by sample collection timing.
Consistent with previous reports, the microbiota of the preterm infants in our study were most frequently dominated by the phylum Proteobacteria; Firmicutes was the next most frequent phylum. Nevertheless, our observations support the hypothesis that temporal and geographic factors influence the intestinal microbial colonization of preterm infants. These temporal and geographic differences were not explained by factors known to influence the intestinal microbiome, such as antibiotics. Differences in Proteobacteria colonization between 2010 and 2011 remained even after modeling to adjust for the effect of antibiotic use on colonization patterns. Microbial shifts observed at the phylum level involved multiple taxa. Differences within NICUs over time of phylum Firmicutes largely included organisms of the class Bacilli, and within hospital 1, families Enterococcaceae and Clostridiaceae, while differences within NICUs over time of phylum Proteobacteria largely involved organisms of the family Enterobacteriaceae.
Differences in hospital environments over time could be due to changes in environmental management, clinical practices, ongoing quality improvement initiatives, patient population, or other factors. Hospital 1 has a necrotizing enterocolitis (NEC) task force which implemented changes to the standard of care in its NICU over the course of this study to reduce the NEC incidence rate. Changes included the introduction of a standard feeding protocol and use of donor milk when mother’s own milk is unavailable, both of which could potentially influence the developing microbiome of the infant. In February 2010, the hospital 2 NICU moved to a new location. Part of the change in infant colonization during this time could have resulted from the change in location and its evolving microflora. While we report these events as qualitative background information, our study was not designed to analyze the impact of specific events quantitatively.
We speculate that the observed patterns could be due to shifts in the microbiota of the NICU environment itself. Longer-term studies are needed to determine the extent to which temporal shifts in colonization recur in regular patterns, occur by the introduction of specific organisms, or are driven by changes in human behavior or the built environment. Studies with concurrent sampling of both infants and the hospital environment are needed to determine the extent to which hospital resident microbes influence infant colonization patterns and the changes in infant colonization patterns over time.
This report presents the microbial colonization patterns of relatively healthy infants and does not address disease outcomes per se. Nevertheless, the differences that we observed could have important clinical consequences. During the period of this study, the hospital NICUs varied in their rates of late-onset sepsis and NEC, two conditions in which intestinal microbiota are implicated. The NICU of hospital 2 had double the rate of late-onset sepsis observed at hospital 1, though sepsis rates remained constant between 2010 and 2011. Rates of NEC at hospital 1 were approximately double those at hospital 2 in 2010, but in 2011, the rates in the two hospitals were indistinguishable.
Nationally, 10% of preterm infants born <29 weeks gestational age develop NEC. This risk has been attributed in part to immature response to LPS-bearing organisms, specifically, excessive TLR4 signaling and hyper-inflammatory response to Proteobacteria. Investigators have noted that a surge or “bloom” in Proteobacteria (typically, Enterobacteriaceae) or a decline in Firmicutes precedes NEC[10–12]. The tendency towards surge in Proteobacteria and equivalent decline in Firmicutes may differ between NICUs. Indeed, in our study, in addition to observing differences between intestinal colonization between the hospitals by week of life, we also observed differences in microbial succession within infants over the first few weeks of postnatal life. Hospital 1 infants experienced a more rapid turnover of bacterial OTUs than did hospital 2 infants. This could have important implications for disease risk.
This study has several limitations: This study was a secondary analysis of a dataset collected for analysis of microbiota prior to onset of NEC. Data from both hospitals and years was available only for the first 2 weeks of postnatal life for most infants; thus, our comparison of hospitals is restricted to those weeks, and we examined microbial succession in the subset of infants with all 3 weeks of data available. Variation in storage protocol due to logistical challenges in the study resulted in use of thioglycollate in one hospital and not the other, which reduced our ability to confidently identify differences in Bacteroidetes, Propionibacteriaceae, and Leuconostocaceae, which are relatively of low abundance in preterm infants. While this methodologic variation is unfortunate, it did not impair the analysis of major taxa in this study. Furthermore, as in all epidemiologic studies, the findings of this study are susceptible to confounding over time or between hospitals. However, we found that the clinical characteristics of infants were generally well matched between hospitals and years, and modeling found that study years within hospitals remained the strongest predictor of the relative abundance of the Proteobacteria after controlling for potential confounding by clinical variables or methodologic factors.
Time of birth and hospital correlate with distinct premature infant intestinal colonization patterns. Identifying these differences in various institutions and over years may provide a biomarker for monitoring disease risk. Our data provide additional evidence for caution regarding the degree to which studies of the microbiome and diseases of prematurity can be generalized to other hospitals or even to the same hospital later in time. Indeed, our data supports the need for large, multi-site, multi-year epidemiological studies in preterm infants and other patient populations to understand the association between intestinal microbial colonization and disease risk. The cohort study from which infants were selected for this study contained too few subjects with NEC and sepsis to effectively assess temporal changes in case infants. This study also clearly indicates a need for careful matching strategies on time and site to examine disease risk in studies of the intestinal microbiota.
Data is available on NCBI SRA, Accession PRJNA63661. Additional metadata and the metadata dictionary are available as a supplementary file to this article (see Additional file1). This file contains all the metadata used in this article. The OTU table and taxa assignments after rarefaction are also available as supplementary file to this article (see Additional file2).
Study infants were a subset of the infants enrolled from one level III NICU in Cincinnati, OH (hospital 1), and one level III NICU in Birmingham, AL (hospital 2), in 2010 and 2011 as part of an ongoing cohort study of novel biomarkers for NEC. Study infants were ≤32 weeks gestational age. The average daily census of the NICU in hospital 1 was 46 with a capacity of 60, and the average daily census of the NICU at hospital 2 was 85 with a capacity of 120. From the overall cohort, this study was restricted to infants <29 weeks gestational age, singleton births, who survived free of proven NEC or proven sepsis until discharge, and were free of congenital anomalies. These inclusion criteria were applied to establish fundamental comparability between the two NICUs studied and remove the influence of disease outcomes per se or other potentially extraneous differences. In addition, all infants had at least one stool sample from week 1 or week 2 that was successfully extracted and sequenced using the methods described below. The Institutional Review Boards at Cincinnati Children’s Hospital and the two participating hospitals approved the study. Parent or guardian consent was obtained for all infants included in this study.
Sample collection and storage
Stool samples were collected from infants during the first weeks of life on postnatal days 5, 8, 11, 14, and 21 plus or minus 2 days. Due to intermittency of infant stooling, we categorized samples as having been collected in collection windows defined as 3–9 days of life (week 1), 10–16 days of life (week 2), and 17–23 days of life (week 3) in order to optimize the number of infant samples included for analysis at each time point. These time windows were selected because no study infant stooled prior to day of life 3, which is typical for premature infants. Samples were collected from soiled diapers, immediately refrigerated in the NICU, and transported to the laboratory where they remained in the refrigerator until processing for storage at -80°C. Cryogenic storage in hospital 1 utilized thioglycollate buffer except for collection day 11, which was biobanked without buffer to provide a sample for metabolomic analysis. However, when other samples from the same week were not available, day 11 sample was used for 16S analysis. Cryogenic storage in hospital 2 did not use thioglycollate or any other storage buffer.
To determine the impact of thioglycollate buffer on the microbiome, we compared samples from hospital 1 in 2009—seven stored without and seven stored with thioglycollate—matched on day of life of collection and infant clinical history. We found that samples stored without thioglycollate were enriched in Bacteroidetes, Leuconostocaceae, and Propionibacteriaceae. As samples from the hospital 2 NICU were stored without thioglycollate and samples from the hospital 1 NICU were stored with thioglycollate during the study period, identified between-site differences in these taxa were considered a potential artifact of the storage protocol, and analyses between sites were restricted to the taxa not influenced by storage protocol.
Stool extraction and 16S rDNA sequencing
If stored in thioglycollate, stool was thawed and centrifuged for 10 min at 4,000 × g and the supernatant was removed. For all samples, 100 μL of TE buffer with lysozyme and proteinase K was added to 0.24 g of thawed stool and vortexed for 10 min. An amount of 1.2 mL of buffer RLT with beta-mercaptoethanol was added to the sample and transferred to sterile bead beating tubes containing 0.3 g of 0.1-mm glass beads. Samples were homogenized for 3 min in a bead beater and centrifuged at 4,000 × g for 5 min to pellet debris. Supernatant was transferred to a clean microcentrifuge tube and spun at 4,000 × g for an additional 2 min to remove remaining debris. Supernatant was then transferred to a Qiagen AllPrep DNA spin column, and DNA was isolated using the Qiagen AllPrep DNA/RNA mini kit (Qiagen, Valencia, CA, USA).
Using extracted DNA, 180-nt paired-end reads were generated using established primers and protocols, with samples allocated across multiple Illumina MiSeq (Illumina, San Diego, CA, USA) runs. Read pairs were merged to create amplicon-spanning sequences that were then filtered to remove those with less than 70% identity to any read in the rRNA16S.gold.fasta reference set (http://drive5.com/uchime/uchime_download.html) using “usearch -usearch_global -id 0.70.” Using the UPARSE pipeline, software version usearch7.0.959_i86linux64, 79,076,883 sequences were processed. The following commands were used with default settings unless otherwise specified. Dereplication resulted in 35,605,130 sequences (-derep_fulllength); removal of singleton reads in 2,206,563 sequences (-sortbysize -minsize 2) and clustering yielded 7,249 OTU representative sequences (-cluster_otus). The OTU table was constructed by mapping reads to OTUs (-usearch_global -strand plus -id 0.97) and applying the python script uc2otutab.py (http://drive5.com/python/). Additional chimera filtering was not applied. QIIME version 1.6 was used to provide classifications of the OTU representative sequences using the gg_13_5 GreenGenes taxonomy and representative sequences constructed at 99% similarity. A phylogenetic tree was constructed within the QIIME package using FastTree and filtered PyNAST alignments of the OTU representative sequences. OTUs with a minimum count fraction of 0.0002 were removed from the OTU table in QIIME, resulting in 525 unique OTUs.
Comparisons between years and hospitals were restricted to week 1 and 2 samples, as hospital 2 had only a single sample from week 3 in 2011.
Differences in clinical characteristics among groups by week were tested using Fisher’s exact test for categorical variables and t test for continuous variables. To standardize comparisons of microbiota, we rarefied the OTU table to 2,000 reads per sample. Rarefaction randomly selects reads from the complete set obtained for each sample until the specified number of reads is obtained. This means that each sample has an equal chance of including rare OTUs so that the samples can be compared.
Alpha diversity was calculated for weeks 1 and 2 using two metrics: Simpson Diversity Index (1-D) and Chao1. Kruskal-Wallis (KW) was used to test for differences in alpha diversity by year and hospital. To examine beta diversity, we used non-metric dimensional scaling (NMDS) to ordinate the microbial communities based on both the unweighted and weighted UniFrac distance calculated in QIIME as described in Morrow et al.. The unweighted UniFrac examines presence/absence only while the weighted UniFrac accounts for abundance differences.
Significant differences in specific taxa between hospitals and by year were identified by linear discriminant analysis effect size (LEfSe). GEE models were used assuming linear and logistic relationships to test the association of taxa identified by LEfSe with hospitals after adjustment with other potential confounders. A backwards elimination approach was used to remove non-significant covariates from the models. Samples from the same infant in different weeks were included in the same model.
Differences in degree of succession between hospitals were tested using the Jaccard index. Values for the Jaccard index were calculated between the weeks 1 and 3 samples from 28 infants who had samples collected across the first 3 weeks of life. The Jaccard index is a distance metric used to show similarity over time. Identical communities will have a Jaccard index value of 0 while completely non-overlapping communities will have a Jaccard index value of 1. KW was used to test for differences in the intra-subject Jaccard index values by hospital and year.
We gratefully acknowledge Estelle Fischer and the Perinatal Institute clinical research core for clinical data collection and data entry; Barbara Davidson for sample collection planning and management oversight; Myra Johnson, Amy LeFevers, and Maria Hughes for sample collection and management; and Donna Wuest for assistance with manuscript preparation. This project has been funded in whole or in part with Federal funds from the National Institute of Child Health and Human Development, National Institutes of Health (grant numbers HD059140, HD13021, and HD27853); National Center for Research Resources, National Institutes of Health (grant numbers 5UL1RR026314-03, U54HG004969); National Human Genome Research Institute, National Institutes of Health (HG005969); Danone Research (PLF-5972-GD); National Institute of Environmental Health Sciences, National Institutes of Health (T32 ES 10957-11); and the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under Contract No. HHSN272200900018C.
- Gordis L: Epidemiology. 2009, Philadelphia, PA: Saunders Elsevier, 375-4Google Scholar
- De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, Collini S, Pieraccini G, Lionetti P: Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A. 2010, 107: 14691-14696. 10.1073/pnas.1005963107. published online EpubAug 17 (10.1073/pnas.1005963107)View ArticlePubMedPubMed CentralGoogle Scholar
- Lee S, Sung J, Lee J, Ko G: Comparison of the gut microbiotas of healthy adult twins living in South Korea and the United States. Appl Environ Microbiol. 2011, 77: 7433-7437. 10.1128/AEM.05490-11. published online EpubOct (10.1128/aem.05490-11)View ArticlePubMedPubMed CentralGoogle Scholar
- Rintala H, Pitkaranta M, Toivola M, Paulin L, Nevalainen A: Diversity and seasonal dynamics of bacterial community in indoor environment. BMC Microbiol. 2008, 8: 56-10.1186/1471-2180-8-56. doi:10.1186/1471-2180-8-56View ArticlePubMedPubMed CentralGoogle Scholar
- Torrazza RM, Ukhanova M, Wang X, Sharma R, Hudak ML, Neu J, Mai V: Intestinal microbial ecology and environmental factors affecting necrotizing enterocolitis. PLoS One. 2013, 8: e83304-10.1371/journal.pone.0083304. doi:10.1371/journal.pone.0083304View ArticlePubMedPubMed CentralGoogle Scholar
- Brooks B, Firek BA, Miller CS, Sharon I, Thomas BC, Baker R, Morowitz MJ, Banfield JF: Microbes in the neonatal intensive care unit resemble those found in the gut of premature infants. Microbiome. 2014, 2: 1-10.1186/2049-2618-2-1. doi:10.1186/2049-2618-2-1View ArticlePubMedPubMed CentralGoogle Scholar
- Altizer S, Dobson A, Hosseini P, Hudson P, Pascual M, Rohani P: Seasonality and the dynamics of infectious diseases. Ecol Lett. 2006, 9: 467-484. 10.1111/j.1461-0248.2005.00879.x. published online EpubApr (10.1111/j.1461-0248.2005.00879.x)View ArticlePubMedGoogle Scholar
- Bhuiyan TR, Qadri F, Saha A, Svennerholm AM: Infection by Helicobacter pylori in Bangladeshi children from birth to two years: relation to blood group, nutritional status, and seasonality. Pediatr Infect Dis J. 2009, 28: 79-85. 10.1097/INF.0b013e31818a5d9d. published online EpubFeb (10.1097/INF.0b013e31818a5d9d)View ArticlePubMedGoogle Scholar
- Pascual M, Bouma MJ, Dobson AP: Cholera and climate: revisiting the quantitative evidence. Microbes Infect/ Institut Pasteur. 2002, 4: 237-245. 10.1016/S1286-4579(01)01533-7. doi:237. published online EpubFebView ArticleGoogle Scholar
- Morrow AL, Lagomarcino AJ, Schibler KR, Taft DH, Yu Z, Wang B, Altaye M, Wagner M, Gevers D, Ward DV, Kennedy MA, Huttenhower C, Newburg DS: Early microbial and metabolomic signatures predict later onset of necrotizing enterocolitis in preterm infants. Microbiome. 2013, 1: 13-10.1186/2049-2618-1-13. doi:10.1186/2049-2618-1-13View ArticlePubMedPubMed CentralGoogle Scholar
- Mai V, Young CM, Ukhanova M, Wang X, Sun Y, Casella G, Theriaque D, Li N, Sharma R, Hudak M, Neu J: Fecal microbiota in premature infants prior to necrotizing enterocolitis. PLoS One. 2011, 6: e20647-10.1371/journal.pone.0020647. doi:10.1371/journal.pone.0020647View ArticlePubMedPubMed CentralGoogle Scholar
- Claud EC, Keegan KP, Brulc JM, Lu L, Bartels D, Glass E, Chang EB, Meyer F, Antonopoulos DA: Bacterial community structure and functional contributions to emergence of health or necrotizing enterocolitis in preterm infants. Microbiome. 2013, 1: 20-10.1186/2049-2618-1-20. doi:10.1186/2049-2618-1-20View ArticlePubMedPubMed CentralGoogle Scholar
- Weaver LT, Lucas A: Development of bowel habit in preterm infants. Arch Dis Child. 1993, 68: 317-320. 10.1136/adc.68.3_Spec_No.317. published online EpubMarView ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, Gormley N, Gilbert JA, Smith G, Knight R: Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6: 1621-1624. 10.1038/ismej.2012.8. published online EpubAug (10.1038/ismej.2012.8)View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC: UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013, 10: 996-998. 10.1038/nmeth.2604. published online EpubOct (10.1038/nmeth.2604)View ArticlePubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303. published online EpubMay (10.1038/nmeth.f.303)View ArticlePubMedPubMed CentralGoogle Scholar
- Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C: Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12: R60-10.1186/gb-2011-12-6-r60. doi:10.1186/gb-2011-12-6-r60)View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.