Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2011; 6(10): e25959.
Published online 2011 October 10. doi:  10.1371/journal.pone.0025959
PMCID: PMC3189942

Molecular Identification of Bacteria in Tracheal Aspirate Fluid from Mechanically Ventilated Preterm Infants

Josef Neu, Editor



Despite strong evidence linking infections to the pathogenesis of bronchopulmonary dysplasia (BPD), limitations of bacterial culture methods have precluded systematic studies of airway organisms relative to disease outcomes. Application of molecular bacterial identification strategies may provide new insight into the role of bacterial acquisition in the airways of preterm infants at risk for BPD.


Serial (within 72 hours, 7, 14, and 21 days of life) tracheal aspirate samples were collected from 10 preterm infants with gestational age ≤34 weeks at birth, and birth weight of 500–1250 g who required mechanical ventilation for at least 21 days. Samples were analyzed by quantitative real time PCR assays for total bacterial load and by pyrosequencing for bacterial identification.


Subjects were diagnosed with mild (1), moderate (3), or severe (5) BPD. One patient died prior to determination of disease severity. 107,487 sequences were analyzed, with mean of 3,359 (range 1,724–4,915) per sample. 2 of 10 samples collected <72 hours of life contained adequate bacterial DNA for successful sequence analysis, one of which was from a subject exposed to chorioamnionitis. All other samples exhibited bacterial loads >70copies/reaction. 72 organisms were observed in total. Seven organisms represented the dominant organism (>50% of total sequences) in 31/32 samples with positive sequences. A dominant organism represented>90% of total sequences in 13 samples. Staphylococcus, Ureaplasmaparvum, and Ureaplasmaurealyticum were the most frequently identified dominant organisms, but Pseudomonas, Enterococcus, and Escherichia were also identified.


Early bacterial colonization with diverse species occursafter the first 3 days of life in the airways of intubated preterm infants, and can be characterized by bacterial load and marked species diversity. Molecular identification of bacteria in the lower airways of preterm infants has the potential to yield further insight into the pathogenesis of BPD.


Bronchopulmonary dysplasia (BPD) was initially attributed to mechanical ventilation and oxygen induced lung injury in infants born preterm, as first described by Northway and colleagues in 1967 [1]. Inflammation of the developing lung from these interventions and other causes is a predominate theory underlying thepathogenesisof BPD [2], [3]. Both prenatal and postnatal sources of inflammation, including infection, may directly disrupt lung growth [4] and predispose infants to BPD [5], [6]. Because fetal development takes place in a normally bacteria-free environment [7], premature exposure and colonization of the neonatal respiratory tract by bacteria may impact the neonatal immune response, increase the risk for early and late infections, and contribute to the subsequent development of BPD. Recent studies have shown that airway colonization with specific organisms during the neonatal period was associated with persistent wheeze in childhood [8], suggesting that early specific bacterial exposure and/or the host response contribute to future disease. Yet, little is known about the timing of bacterial exposure and acquisition in the airways of preterm infants, mechanisms that promote or impede colonization of specific organisms, or the inflammatory responses to colonizing bacteria in relation to injury to the developing lung.

Recently developed molecular methods have suggested that the lungs of healthy older children and adults are rich with microbes, which has led to a new understanding of airway colonization and has challenged the traditional theory that human lower airways are sterile [9], [10]. These methods further suggest that traditional culture methods are not adequate to investigate the contribution of microbial communities to disease. In addition to common bacterial organisms that can be readily cultured, there has been clear recognition of the potential pathological impact of other bacteria, such as Ureaplasma, Mycoplasma, and Chlamydia, as causes of perinatal infection that may increase the risk for the subsequent development of BPD [11], [12], [13], [14]. Thus, traditional culture methods are likely not sensitive enough to provide meaningful information on many potential pathological organisms that can affect outcomes in preterm neonates. Culture-independent molecular techniques offer the opportunity to identify the full spectrum of bacteria in biological samples, and these methods may provide new insight into the role of bacterial acquisition in the respiratory tract of newborn preterm infants and the development of BPD. Molecular bacterial identification methods can detail the total bacterial load in a sample, the diversity of bacteria present, and relative abundance of specific bacteria based on the relative proportion of sequences in a sample.

Thus, we hypothesized that early colonization with diverse microbial agents occurs early after birth in ventilated preterm infants, and that application of molecular bacterial identification strategies to serial airway samples can define changes in the microbiome in infants at risk for chronic lung disease. In order to investigate the airway microbiome in patients at high risk for BPD, we performed molecular identification of bacteria in serial tracheal aspirate samples from intubated preterm infants by isolating rRNA genes harvested by PCR to determine the composition and the timing of acquisition of bacterial communities in the lower airways. We compared these results to those obtained by clinically indicated standard laboratory cultures and described the timing of antibiotic administration in relation to the identification of airway bacterial communities.


Study Population

All data and specimens were obtained as part of a prospective observational research protocol, approved by the Colorado Multiple Institutional Review Board. Written informed consent was received from the parents/guardians of all participants. Assent was not obtained because all participants were infants. Criteria for enrollment included a gestational ageof 34 weeks or less, birth within the previous 72 hours, and a birth weight of 500 to 1250 g. To be included in this study analysis, infants also needed to be intubated at enrollment and remain intubated for at least 21 days with successful collection of tracheal aspirate samples at each of 4 times points (within 72 hours of birth, and 7, 14, and 21 days of life). Exclusion criteria werelethal congenital anomalies, complex congenital heart disease (includingan atrial septal defect larger than 1 cm and a ventricular septaldefect larger than 2 mm), or anticipated death prior to hospital discharge. Maternal medical history including complications of pregnancy and birth, and subject data, including clinically obtained respiratory tract cultures and antibiotic administration from the NICU course were collected as part of the prospective study. BPD status was assessed at 36 weeks post conceptual age via head hood oxygen reduction test as previously described [15], anddisease severity categorieswere adjusted for altitude in Denver, Colorado (1600 m).

Specimen Collection and processing

Tracheal aspirates were collected by instillation of 0.5 mlof sterile 0.9% saline (without preservative) via the endotracheal tube (ETT) and suction applied to a sterile catheter inserted to 0.5 cm below the tip of the ETTas the catheter was slowly withdrawn. The suction catheter was cleared of retained secretions by aspirating an additional 1 ml of the sterile saline. Specimens were immediately spun at 250× g for 20 minutes at °C. The supernatant was removed and the remaining pellet was resuspended with 0.2 mlnormal saline and stored at −70°C.

DNA Extraction, PCR and Sequencing

DNA was prepared from each sample using the Qiagen EZ1 Advanced platform. 200 µl of clinical sample was extracted using the bacterial card and tissue kit per manufacturer's instructions. The bacterial load of each sample was estimated using a TaqMan quantitative PCR (qPCR) assay described by Nadkarni [16]. Each DNA sample was assayed in triplicate, and assays with coefficient of variation (standard deviation/mean) greater than 20% were repeated. We have validated the reproducibility of this assay in comparable airway samples from children with cystic fibrosis [17]. Copy number was established using a standard curve assembled from 10-fold serial dilutions of plasmid DNA containing a cloned rRNA gene.

DNA sequencing was performed using barcoded primers (27F–338R) compatible with the Roche 454 pyrosequenceras previously reported [18]. Each DNA was amplified in triplicate PCR assays, which were pooled prior to assessing the product by agarose gel electrophoresis. A negative control was run for each barcode in parallel. Any assays where the negative control was positive were repeated. Amplicons of the correct length were arrayed in 96-well format for normalization. Amplicon concentration was normalized using the SequalPrep normalization plate (Invitrogen) prior to mixing [19]. Equal volumes were mixed, concentrated by evaporation and gel purified using the Montage Gel Purification Kit (Millipore). The resulting amplicon pool was sequenced per manufacturer's instructions using the Roche Genome Sequencer FLX system at the Consortium for Comparative Genomics sequencing facility at the University of Colorado Denver.

Sequence analysis

Sequence data (fasta and quality files) were assigned to samples by bar code and screened for low level quality defects (short sequences <150 nt in length, >1 sequence ambiguity, best read with quality ≥20 over a 10 nt moving window) by the software program BARTAB [20]. Non-bacterial sequences were removed from the dataset by requiring a close match with a bacterial secondary structure rRNA model within Infernal [21]. Sequences identified as potential chimeras by ChimeraSlayer [22] were also removed from the dataset. Taxonomic assignment was done by the Ribosomal Database Project Classifier software [23]. We used a BLAST database containing bacterial isolates identified in the Silva ARB database (version 104) to assign species names [24]. The taxonomic information was used to construct sequence groups with identical taxonomic rank, which were used to calculate ecology statistics for each sample. Shannon diversity and Chao1 were computed, using rarefaction (100 replicates) to compensate for differences in number of sequences per sample, using the biodiv functionality contained within the XplorSeq Toolkit [25], [26].

Statistical analysis

Data were presented graphically, and descriptive statistics were calculated using SAS v 9.2 (Cary, NC). The percent of total sequence count was calculated for each individual organism. Multilevel plots were constructed, displaying the percent of total sequences, clinical culture results, and antibiotic courses for each infant. Bacterial load and Shannon diversity index were plotted individually for each patient. Formal statistical testing was not performed.


Of 49 subjects enrolled in the larger observational cohort at the time of analysis, 10 patients met the ventilation and sampling criteria for inclusion in this study. Clinical characteristics of the 10 patients are listed in Table 1. Patient 1 died at day of life 40 from methicillin-resistant Staphylococcus aureus sepsis; this patient did not have an autopsy performed. The remaining patients were all diagnosed with BPD, ranging from mild to severe. All patients received broad spectrum antibiotics, empirically or for confirmed infection, at some point during the first 21 days of life (Figure 1).

Figure 1
Bacterial Composition of Tracheal Aspirate Samples for Each Patient by Day of Collection.
Table 1
Patient Characteristics.

The average yield of DNA from each sample was 32.9 ng/µl (range: 0.4–143.6). A total of 107,487 sequences were analyzed, with an average number of 3,359 (range: 1,724–4,915) per sample. Only 4 of 10 samples collected <72 hours of life had detectable bacterial signal by qPCR to assess bacterial load. Further, only 2 of these samples contained adequate load for successful sequence analysis of the bacterial community. All samples collected from day 7–21 of life contained detectable bacterial load (>70 copies/reaction), and each sample was successfully amplified for bacterial identification. The distribution of bacterial loads by day of collection for each subject is shown in Figure 2A. Seventy-two organisms were observed in total. The average number of organisms identified per sample was 7 (range 2–23). Sequences observed ≥1.0% of total in a sample represent 19organisms (96.5% of total sequences; 1–9 organisms per sample). Chao1, a nonparametric estimate of actual species richness in a sample, predicted on average, 8 organisms in each sample (range 1–46). The diversity of organisms in samples, determined by the Shannon Index, is displayed by day of collection for each subject in Figure 2B.

Figure 2
Comparison of Bacterial Load and Shannon Diversity Index by Age.

A dominant organism, defined as an organism representing >50% of total sequences was identified in31 of 32 samples with positive sequences(7 organisms represented). A dominant organism represented >90% of total sequences in 13 samples (6 organisms). Distributions of bacterial sequences in the samples for each patient are displayed in Figure 1 along with results of clinically obtained cultures and the timing of antibiotic administration. Ureaplasmaparvum was the dominant organism in the two positive samples from <72 hours, one of which was from a subject exposed to chorioamnionitis (subject 4). Ureaplasmaparvum or Ureaplasmaurealyticum represented the dominant organism in 9 samples from 6 subjects, 5 samples of which were from subjects without evidence of Ureaplasma on the initial sample. Staphylococcus species represented the dominant organism in 19 samples (17 were coagulase negative Staphylococcus). Of these, Staphylococcus represented>90% of sequences in7 samples. Staphylococcus aureus was the dominant species in 2 samples (both collected on day 21). Pseudomonasaeruginosa (1 sample; day 14), Enterococcus faecalis (1 sample; day 21), and Escherichia coli (1 sample; day 21) were the other dominant organisms identified in samples.

Based on suspected lung infection by the physicians managing these patients during the course of the study, 19 additional tracheal aspirate samples were obtained for bacterial culture within at least 72 hours of research tracheal aspirate collections (Figure 1). Of these, 14 (74%) had positive bacterial growth, with specific organisms identified in 11 samples and “mixed upper respiratory flora (MURF)” reported in the remaining 3 samples. Sequencing identified all bacteria that grew in the clinically indicated cultures (n = 11 positive cultures). Six of these 11 cultures grew 2 organisms, however, the bacteria identified by culture did not always represent the organisms with the most sequences identified by molecular methods (n = 5 cultures). Non-cultivatable bacteria, such as Ureaplasma, were a significant source of discrepancies between the molecular strategies versus standard culture approaches.

Potential pathogens identified by sequencing that were not identified by culture include Ureaplasma (3 cultures), Pseudomonas (1 culture), and E. coli (1culture). Of the 5 samples that were obtained for clinical indications that had no bacterial growth, the paired research samples revealed Staphylococcus species as the dominant organism in 3 samples, Ureaplasma species as the dominant organism in one sample, and no bacterial amplification in the remaining sample. Systemic antibiotics were administered for some period of time during the first 21 days of life in all patients. The median duration of antibiotic administration was 12 days (range, 2–21). The median duration on 2 or more antibiotics was 10 days (range, 2–19). All samples collected after 72 hours of life that had Ureaplasma species as the dominant organism (n = 7) had been treated with 2 or more antibiotics in the interval between research samples.


We performed molecular identification of bacteria in serial tracheal aspirate samples from intubated preterm infants by isolating rRNA genes harvested by PCR to determine the composition of bacterial communities and the timing of acquisition in the airways of these infants. We found that during the first 72 hours after birth, the majority of tracheal aspirate samples had low or undetectable bacterial loads. However, within 7 days of life, sustained bacterial colonization with diverse species is present in the airways of intubated preterm infants. 96.9% of samples with bacterial detection exhibited a dominant organism, with Staphylococcus and Ureaplasma species most frequently identified as the dominant organisms. Molecular methods identified all species detected by traditional clinical cultures, but also identified potential pathogens, including Ureaplasma, that were not detected by standard bacterial culture methods. Thus, improved resolution of bacterial detection with regard to overall bacterial load, species diversity within the microbiome or the identification of specific organisms for directed therapy, has the potential to yield further insights into the pathogenesis of BPD and late respiratory outcomes of preterm newborns.

These findings are important because early respiratory infection and inflammation have been shown to impair alveolar development [27] and contribute to the development of BPD [28]. Inflammatory injury disrupts lung growth in part by production of pro-inflammatory cytokines [29]. Although pro-inflammatory cytokines have been associated with the development of BPD and adverse outcomes [30], [31], [32], variable assessments of the pro-inflammatory response, and lack of context specificity (identification of all aspects of the perinatal milieu contributing to disease) have precluded the use of cytokines as adequate biomarkers for BPD. While antenatal and postnatal infections, hyperoxia, and mechanical ventilation have garnered much attention as inflammation provoking stimuli, the possible impact of airway colonization on the development of BPD is less clear and has generated contradictory results [33]. The presence of Ureaplasma in the respiratory tract of preterm infants, even without signs of infection, has been associated with a pro-inflammatory response and an increased risk for BPD [34]. Although earlier trials of erythromycin treatment for the treatment of Ureaplasma in preterm infants failed to show benefit [35], [36], [37], [38], a recent trial of azithromycin treatment revealed a significant reduction in the endpoint of BPD or death in the subgroup of patients who were found to be colonized or infected with Ureaplasma [39]. Thus, Ureaplasma remains a plausible but still unconfirmed contributor to BPD. Ureaplasma was the dominant organism in at least one sample for 6 of the 10 infants in our study. Ureaplasma was found to be the dominant organism in late samples even when not detected in initial samples, suggesting the possibility of post-natal acquisition and/or the presence of environmental factors, such as antibiotic administration, that encourage selection of this genus. The association between antibiotic administration, Ureaplasma colonization, and risk for BPD should be explored further in future studies. Exposure to colonizing bacteria may impact the preterm infant's immune response and subsequent lung development. Improved resolution of bacterial detection techniques as demonstrated here in conjunction with assessments of inflammation could improve our understanding of the developing host immune response in premature infants and its relationship to pathogen exposure.

Molecular methods have contradicted the traditional theory that human lower airways are sterile by demonstrating that the lungs of healthy older children and adults are rich with microbial flora [9], [10], and that traditional culture methods are not adequate to investigate the contribution of microbial communities to disease. Our results, together with a previous report [40] suggest that microbial communities are present early in the lower airways of intubated preterm infants. Rapidly accumulating evidence demonstrates that microbial communities play an important role in human physiology [41], [42]. Several diseases have been associated with changes in composition and diversity of these microbial communities [43], [44], [45]. Some have shown that microbiota of limited diversity and relative absence of commensal organisms are associated with increased inflammation, barrier permeability, and disease status [46], [47]. Early evaluations of the respiratory tract microbiota in pulmonary diseases such as cystic fibrosis, asthma, and COPD [9], [10] also reveal shifts in microbial composition compared to healthy patients, suggesting the composition of the respiratory tract microbiota can play a direct role in pathophysiology of lung disease. Neonatal oropharyngeal colonization with Streptococcus pneumoniae, Moraxella catarrhalis, Haemophilusinfluenzae, or a combination of these organisms, but not colonization with Staphylococcusaureus, was significantly associated with persistent wheeze in childhood [8]. Children living on farms are exposed to a wider range of microbes than were children living in urban environments, and these children appear relatively protected from asthma and atopy [48]. Together, these studies suggest the respiratory tract microbiota early in life can also modulate the host immune response and future lung disease.

Only 2of 10 samples collected within 72 hours after birth exhibited sufficient bacterial load for sequence analysis, suggesting the lower airways of preterm infants typically have comparatively low bacterial loads. One of these two samples was obtained from an infant exposed to chorioamnionitis and revealed Ureaplasmaparvum, a known cause of chorioamnionitis, as the dominant organism. Although a diverse array of bacteria were detected within the samples of these 10 infants, the diversity of organisms within samples was far less than detected in respiratory samples of older children and adults [9], [10]. Further, bacterial diversity was not associated with age of the infant or BPD severity.

The use of antibiotics administered to treat confirmed or potential infections may have long term effects on the composition of the airway microbial community. Potential negative consequences include the elimination of commensal bacteria and selection for other pathogenic bacteria, especially those that are not readily detected by clinical culture techniques and that promote chronic airway inflammation. Further studies will need to investigate the potential role of antibiotics on the composition of the airway microbiota and the impact on early lung development.

There are several limitations to our study. First, the small sample size in this pilot study precludes rigorous evaluation of the presence of specific bacteria, bacterial loads, or species diversity to the risk of BPD. All of the infants included in this study were less than 28 weeks gestation and required mechanical ventilation for at least 21 days, which represents a very high risk population for severe BPD. Serial evaluation of the airway microbiome in preterm infants who do not develop BPD will be important to determine how the composition of the airway microbiome is related to risk of BPD. Since preterm infants who do not develop BPD often do not require long periods of mechanical ventilation support, noninvasive methods of airway sampling will be required. Second, all infants in this study were treated with systemic antibiotics which are likely to impact the evolution of microbial composition, limiting the ability to ascertain the natural acquisition of airway colonization after birth. The presence of endotracheal tubes also impacts the access of microbes to lower airways. Although we attempted to compare molecular detection methods to clinical culture results, the samples were not simultaneously obtained. Therefore, the discrepancies between the two results cannot be accounted for by methodology alone. In addition, sampling of tracheal aspirates may not adequately reflect bacterial exposure at the alveolar level. Simultaneous sampling of the upper airway, lower airways and alveolar samples are required to provide a more complete assessment of bacterial colonization. Finally, future studies are needed to more directly link the presence of specific airway organisms, airway bacterial load, and diversity of the microbiota with altered innate immunity and inflammatory cytokines to better understand how early colonization impacts lung inflammation.

In summary, early bacterial colonization with diverse species are present in the airways of intubated preterm infants, and can be characterized by bacterial load and species diversity. Molecular detection of Ureaplasma species in the airways of preterm infants late in the course of mechanical ventilation provides further support for the concept that Ureaplasma may be involved in the pathogenesis of BPD and requires further investigation. Molecular identification of bacteria in the lower airways of preterm infants has the potential to yield further insight into the pathogenesis of BPD and to identify novel therapeutic targets.


Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was funded by National Institutes of Health, National Center for Research Resources K23 RR021021, and Colorado Clinical and Translational Sciences Institute UL1 RR025780,; National Heart Lung Blood Institute R01 HL085703,; Thrasher Research Fund, The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Northway WH, Jr, Rosan RC, Porter DY. Pulmonary disease following respirator therapy of hyaline-membrane disease. Bronchopulmonary dysplasia. N Engl J Med. 1967;276:357–368. [PubMed]
2. Groneck P, Speer CP. Inflammatory mediators and bronchopulmonary dysplasia. Arch Dis Child Fetal Neonatal Ed. 1995;73:F1–3. [PMC free article] [PubMed]
3. Speer CP. Inflammation and bronchopulmonary dysplasia. Semin Neonatol. 2003;8:29–38. [PubMed]
4. Speer CP. Inflammatory Mechanisms in Bronchopulmonary Dysplasia. In: Abman SH, editor. Bronchopulmonary dysplasia. New York: Informa Healthcare; 2010. xvii, 499.
5. Van Marter LJ, Dammann O, Allred EN, Leviton A, Pagano M, et al. Chorioamnionitis, mechanical ventilation, and postnatal sepsis as modulators of chronic lung disease in preterm infants. J Pediatr. 2002;140:171–176. [PubMed]
6. Speer CP. Inflammation and bronchopulmonary dysplasia: a continuing story. Semin Fetal Neonatal Med. 2006;11:354–362. [PubMed]
7. DiGiulio DB, Romero R, Amogan HP, Kusanovic JP, Bik EM, et al. Microbial prevalence, diversity and abundance in amniotic fluid during preterm labor: a molecular and culture-based investigation. PLoS One. 2008;3:e3056. [PMC free article] [PubMed]
8. Bisgaard H, Hermansen MN, Buchvald F, Loland L, Halkjaer LB, et al. Childhood asthma after bacterial colonization of the airway in neonates. N Engl J Med. 2007;357:1487–1495. [PubMed]
9. Hilty M, Burke C, Pedro H, Cardenas P, Bush A, et al. Disordered microbial communities in asthmatic airways. PLoS One. 2010;5:e8578. [PMC free article] [PubMed]
10. Huang YJ, Nelson CE, Brodie EL, Desantis TZ, Baek MS, et al. Airway microbiota and bronchial hyperresponsiveness in patients with suboptimally controlled asthma. The Journal of allergy and clinical immunology. 2011;127:372–381 e373. [PMC free article] [PubMed]
11. Da Silva O, Gregson D, Hammerberg O. Role of Ureaplasma urealyticum and Chlamydia trachomatis in development of bronchopulmonary dysplasia in very low birth weight infants. Pediatr Infect Dis J. 1997;16:364–369. [PubMed]
12. Garland SM, Bowman ED. Role of Ureaplasma urealyticum and Chlamydia trachomatis in lung disease in low birth weight infants. Pathology. 1996;28:266–269. [PubMed]
13. Iles R, Lyon A, Ross P, McIntosh N. Infection with Ureaplasma urealyticum and Mycoplasma hominis and the development of chronic lung disease in preterm infants. Acta Paediatr. 1996;85:482–484. [PubMed]
14. van Waarde WM, Brus F, Okken A, Kimpen JL. Ureaplasma urealyticum colonization, prematurity and bronchopulmonary dysplasia. Eur Respir J. 1997;10:886–890. [PubMed]
15. Walsh MC, Wilson-Costello D, Zadell A, Newman N, Fanaroff A. Safety, reliability, and validity of a physiologic definition of bronchopulmonary dysplasia. J Perinatol. 2003;23:451–456. [PubMed]
16. Nadkarni MA, Martin FE, Jacques NA, Hunter N. Determination of bacterial load by real-time PCR using a broad-range (universal) probe and primers set. Microbiology. 2002;148:257–266. [PubMed]
17. Zemanick ET, Wagner BD, Sagel SD, Stevens MJ, Accurso FJ, et al. Reliability of quantitative real-time PCR for bacterial detection in cystic fibrosis airway specimens. PLoS One. 2010;5:e15101. [PMC free article] [PubMed]
18. Hamady M, Walker JJ, Harris JK, Gold NJ, Knight R. Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex. Nat Methods. 2008;5:235–237. [PMC free article] [PubMed]
19. Harris JK, Sahl JW, Castoe TA, Wagner BD, Pollock DD, et al. Comparison of normalization methods for construction of large, multiplex amplicon pools for next-generation sequencing. Appl Environ Microbiol. 76:3863–3868. [PMC free article] [PubMed]
20. Frank DN. BARCRAWL and BARTAB: software tools for the design and implementation of barcoded primers for highly multiplexed DNA sequencing. BMC Bioinformatics. 2009;10:362. [PMC free article] [PubMed]
21. Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009;25:1335–1337. [PMC free article] [PubMed]
22. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome research. 2011;21:494–504. [PubMed]
23. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and environmental microbiology. 2007;73:5261–5267. [PMC free article] [PubMed]
24. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic acids research. 2007;35:7188–7196. [PMC free article] [PubMed]
25. Frank DN. XplorSeq: A Software Environment for Integrated Management and Phylogenetic Analysis of Metagenomic Sequence Data. BMC Bioinformatics. 2008;9:420. [PMC free article] [PubMed]
26. Frank DN. XplorSeq tool kit: Tools for nucleic acids sequence analysis. Phyloware website. 2010. Available at: Accessed 2011 Sept 15.
27. Kallapur SG. Contribution of inflammation to lung injury and development. Archives of Disease in Childhood - Fetal and Neonatal Edition. 2006;91:F132–F135. [PMC free article] [PubMed]
28. Speer CP. Chorioamnionitis, Postnatal Factors and Proinflammatory Response in the Pathogenetic Sequence of Bronchopulmonary Dysplasia. Neonatology. 2009;95:353–361. [PubMed]
29. Bose CL, Dammann CE, Laughon MM. Bronchopulmonary dysplasia and inflammatory biomarkers in the premature neonate. Arch Dis Child Fetal Neonatal Ed. 2008;93:F455–461. [PubMed]
30. Jonsson B, Tullus K, Brauner A, Lu Y, Noack G. Early increase of TNF alpha and IL-6 in tracheobronchial aspirate fluid indicator of subsequent chronic lung disease in preterm infants. Arch Dis Child Fetal Neonatal Ed. 1997;77:F198–201. [PMC free article] [PubMed]
31. Cayabyab RG, Jones CA, Kwong KY, Hendershott C, Lecart C, et al. Interleukin-1beta in the bronchoalveolar lavage fluid of premature neonates: a marker for maternal chorioamnionitis and predictor of adverse neonatal outcome. J Matern Fetal Neonatal Med. 2003;14:205–211. [PubMed]
32. Kotecha S, Wilson L, Wangoo A, Silverman M, Shaw RJ. Increase in interleukin (IL)-1 beta and IL-6 in bronchoalveolar lavage fluid obtained from infants with chronic lung disease of prematurity. Pediatr Res. 1996;40:250–256. [PubMed]
33. Groneck P, Schmale J, Soditt V, Stutzer H, Gotze-Speer B, et al. Bronchoalveolar inflammation following airway infection in preterm infants with chronic lung disease. Pediatr Pulmonol. 2001;31:331–338. [PubMed]
34. Viscardi RM, Hasday JD. Role of Ureaplasma species in neonatal chronic lung disease: epidemiologic and experimental evidence. Pediatr Res. 2009;65:84R–90R. [PMC free article] [PubMed]
35. Baier RJ, Loggins J, Kruger TE. Failure of erythromycin to eliminate airway colonization with ureaplasma urealyticum in very low birth weight infants. BMC pediatrics. 2003;3:10. [PMC free article] [PubMed]
36. Jonsson B, Karell AC, Ringertz S, Rylander M, Faxelius G. Neonatal Ureaplasma urealyticum colonization and chronic lung disease. Acta Paediatrica. 1994;83:927–930. [PubMed]
37. Lyon AJ, McColm J, Middlemist L, Fergusson S, McIntosh N, et al. Randomised trial of erythromycin on the development of chronic lung disease in preterm infants. Archives of disease in childhood Fetal and neonatal edition. 1998;78:F10–14. [PMC free article] [PubMed]
38. Waites KB, Sims PJ, Crouse DT, Geerts MH, Shoup RE, et al. Serum concentrations of erythromycin after intravenous infusion in preterm neonates treated for Ureaplasma urealyticum infection. The Pediatric Infectious Disease Journal. 1994;13:287–293. [PubMed]
39. Ballard HO, Shook LA, Bernard P, Anstead MI, Kuhn R, et al. Use of azithromycin for the prevention of bronchopulmonary dysplasia in preterm infants: a randomized, double-blind, placebo controlled trial. Pediatric Pulmonology. 2011;46:111–118. [PubMed]
40. Stressmann FA, Connett GJ, Goss K, Kollamparambil TG, Patel N, et al. The use of culture-independent tools to characterize bacteria in endo-tracheal aspirates from pre-term infants at risk of bronchopulmonary dysplasia. J Perinat Med. 2010;38:333–337. [PubMed]
41. Suerbaum S. Microbiome analysis in the esophagus. Gastroenterology. 2009;137:419–421. [PubMed]
42. Li P, Hotamisligil GS. Metabolism: Host and microbes in a pickle. Nature. 2010;464:1287–1288. [PubMed]
43. de La Serre CB, Ellis CL, Lee J, Hartman AL, Rutledge JC, et al. Propensity to high-fat diet-induced obesity in rats is associated with changes in the gut microbiota and gut inflammation. Am J Physiol Gastrointest Liver Physiol. 2010;299:G440–448. [PubMed]
44. van Vliet MJ, Harmsen HJ, de Bont ES, Tissing WJ. The role of intestinal microbiota in the development and severity of chemotherapy-induced mucositis. PLoS Pathog. 2010;6:e1000879. [PMC free article] [PubMed]
45. Yang L, Lu X, Nossa CW, Francois F, Peek RM, et al. Inflammation and intestinal metaplasia of the distal esophagus are associated with alterations in the microbiome. Gastroenterology. 2009;137:588–597. [PMC free article] [PubMed]
46. Cani PD, Possemiers S, Van de Wiele T, Guiot Y, Everard A, et al. Changes in gut microbiota control inflammation in obese mice through a mechanism involving GLP-2-driven improvement of gut permeability. Gut. 2009;58:1091–1103. [PMC free article] [PubMed]
47. Dimmitt RA, Staley EM, Chuang G, Tanner SM, Soltau TD, et al. Role of postnatal acquisition of the intestinal microbiome in the early development of immune function. J Pediatr Gastroenterol Nutr. 2010;51:262–273. [PMC free article] [PubMed]
48. Ege MJ, Mayer M, Normand AC, Genuneit J, Cookson WO, et al. Exposure to environmental microorganisms and childhood asthma. The New England journal of medicine. 2011;364:701–709. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science