Skip to main content

Microbial signatures in amniotic fluid at preterm birth and association with bronchopulmonary dysplasia



Microbiome dysbiosis can have long-lasting effects on our health and induce the development of various diseases. Bronchopulmonary dysplasia (BPD) is a multifactorial disease with pre- and postnatal origins including intra-amniotic infection as main risk factor. Recently, postnatal pathologic lung microbiota colonization was associated with BPD. The objectives of this prospective observational cohort study were to describe differences in bacterial signatures in the amniotic fluid (AF) of intact pregnancies without clinical signs or risk of preterm delivery and AF samples obtained during preterm deliveries and their variations between different BPD disease severity stages.


AF samples were collected under sterile conditions during fetal intervention from intact pregnancies (n = 17) or immediately before preterm delivery < 32 weeks (n = 126). Metabarcoding based approaches were used for the molecular assessment of bacterial 16S rRNA genes to describe bacterial community structure.


The absolute amount of 16S rRNA genes was significantly increased in AF of preterm deliveries and detailed profiling revealed a reduced alpha diversity and a significant change in beta diversity with a reduced relative abundance of 16S rRNA genes indicative for Lactobacillus and Acetobacter while Fusobacterium, Pseudomonas, Ureaplasma and Staphylococcus 16S rRNA gene prevailed. Although classification of BPD by disease severity revealed equivalent absolute 16S rRNA gene abundance and alpha and beta diversity in no, mild and moderate/severe BPD groups, for some 16S rRNA genes differences were observed in AF samples. Bacterial signatures of infants with moderate/severe BPD showed predominance of 16S rRNA genes belonging to the Escherichia-Shigella cluster while Ureaplasma and Enterococcus species were enriched in AF samples of infants with mild BPD.


Our study identified distinct and diverse intrauterine 16S rRNA gene patterns in preterm infants immediately before birth, differing from the 16S rRNA gene signature of intact pregnancies. The distinct 16S rRNA gene signatures at birth derive from bacteria with varying pathogenicity to the immature lung and are suited to identify preterm infants at risk. Our results emphasize the prenatal impact to the origins of BPD.


Even more than 50 years after its original description, bronchopulmonary dysplasia (BPD) remains a highly prevalent and devastating morbidity following premature delivery [1]. The diagnosis of BPD poses a high risk of lifelong consequences for the overall health status and quality of life [2, 3]. Due to its pre- and postnatal multifactorial origins, therapeutic approaches to prevent BPD are still very restricted and of limited efficacy [4,5,6]. Prenatally, maternal preeclampsia/HELLP, intrauterine growth restriction and nicotine exposure constitute main risk factors. Postnatally, relevant contributors to the evolution of BPD are oxygen toxicity and mechanical ventilation. Postnatal infections aggravate the injury to the immature lung. Studies in experimental models demonstrated comparable pathomechanisms as observed during hyperoxia and mechanical ventilation [6]. Recently, microbiota colonization of the immature lung was coming into the focus of research and dysbiotic changes were observed in preterm infants after birth that later on induced BPD development [6]. However, the underlying clinical studies were of marked heterogeneity and restricted to preterm infants that required invasive mechanical ventilation where tracheal aspirates were available [7]. In these infants, a reduced alpha and beta diversity was observed, which was a result a reduced relative abundance of Lactobacilli and an increase of Enterobacteriaceae. In another study, bacterial colonization of the upper airways with Enterococcus Enterobacter Staphylococcus and Escherichia species was identified to be associated with BPD [8]. But the impact of microbiota on the immature lung and the risk for BPD remains controversial [9]. The two largest meta-analyses on this topic found an association between chorioamnionitis or the pulmonary colonization with bacteria of the genus Ureaplasma and BPD, but both reports point to the uncertainty due to potential biases in the published data despite the statistically significant impact on BPD [10, 11]. This uncertainty is strengthened by a recent large cohort analysis where the mode of birth, which highly impacts the microbial colonization at birth, had no influence on the incidence and severity of BPD [12].

This inconsistency in evidence for the role of microbiota in BPD development stands in contrast to the results from preclinical studies that unravelled the circumstances of attenuated and aggravated lung injury in connection to the time point and mass of bacterial colonization [9]. A single prenatal lipopolysaccharide (LPS) stimulus in rodents induced BPD typical pathologies [13]. However, in foetal lambs the intrauterine colonization with Ureaplasma did not alter air space and vascular development despite the induction of inflammation. The clinically relevant scenario of multifactorial exposures did not provoke a pulmonary inflammatory response when intra-amniotic Ureaplasma application preceded the exposure to LPS [14, 15]. Comparable results were reported when antenatal steroids were applied before LPS [16]. From these studies, the term immunotolerance induction in the lung was propagated.

While the shaping of the postnatal early life microbiome has attracted particular interest [17], it is an ongoing controversial discussion whether the amniotic cavity represents a sterile environment as long as no rupture of the amniotic membranes occurs [18]. Results from molecular barcoding-based studies on amniotic fluid, placenta and foetal samples questioned the original statement of a sterile environment and quantification by realtime PCR prevailed low levels of microbial DNA. But the presence of bacterial components including their 16S rRNA gene cannot separate DNA leakage into the amniotic cavity and real living bacteria translocation [19,20,21].

A better understanding of the consequences of prenatal pathologies is a prerogative due to their high impact on BPD and the association of in utero microbiota structures with the pulmonary outcome in preterm infants has not been fully elaborated. In our prospective cohort study, we aimed to specify bacterial community structures in the amniotic fluid at birth in premature deliveries in comparison with intact control pregnancies and their association with the severity of BPD.


Study design, data acquisition and parameter definitions

Preterm infants < 32 weeks of gestation were prospectively enrolled within the PROTECT-AIRR study at the perinatal centre of the university hospital of Giessen between July 2015 and May 2020.

Data were recorded from the electronic data management system and paper file records as executed before [22]. Prenatal antibiotic therapy was counted when applied within the last 7 days before delivery. Small for gestational age (SGA) was defined as birth weight < 10th percentile of the German perinatal registry [23]. Nosocomial infections were documented applying the criteria of the German NICU nosocomial surveillance system (NEO-KISS) [24]. Preterm deliveries were categorized as intraamniotic infection whenever one of the following criteria was fulfilled: histologic chorioamnionitis, elevated AF interleukin-6 > 3600 pg/ml, intractable premature contractions under tocolytic therapy or premature rupture of membranes preceding labour.

To determine the presence and severity of BPD, the NICHD consensus definition was applied: Infants on supplemental oxygen for ≥ 28 days of life were categorized as mild BPD if they were stably breathing room air at 36 weeks postmenstrual age. Infants with persistent oxygen requirements < 30% at 36 weeks were grouped as moderate BPD and ≥ 30% of oxygen and/or need for positive pressure support as severe BPD. Therapy with highflow nasal cannulas was graded as continuous positive airway pressure whenever a level of 3-cm H2O or higher was attained and the fraction of oxygen provided by low flow nasal cannulas was calculated as published [24,25,26].

Sample collection

N = 126 infants were included into the analyses whenever an amniotic fluid (AF) sample was available (Fig. 1). AF samples from preterm deliveries were mostly obtained during routine caesarean section. In n = 4 cases samples from amniotic puncture performed immediately before the decision to premature delivery were used. n = 18 patients were excluded from the analyses due to death before 36 weeks of gestation for non-pulmonary reasons (n = 7), necrotising enterocolitis or focal intestinal perforation (n = 8), cardiac surgery (n = 2) or milk aspiration (n = 1) due to the potential high impact on the pulmonary outcome independent of the natural course (Fig. 1). n = 17 control AF samples were obtained under sterile conditions from intact pregnancies without signs for premature delivery during prenatal interventions with a median gestational age of 25 weeks that were subjected to prenatal intervention for clinical indication (Fig. 1).

Fig. 1
figure 1

Study population flow chart. Flow chart of inclusion and exclusion of infants into the study population. Exclusion criteria included morbidities unrelated to the baseline pulmonary status but with high impact on the pulmonary outcome and AF samples with low reads bacterial signal

The highly probable microbial contamination of samples during spontaneous delivery of preterm infants with high loads of microbiota or their particles was excluded as all AF samples were obtained immediately before delivery by amniotic puncture or during caesarean section. AF samples were collected in sterile pyrogen-free protein low-bind tubes (Eppendorf SE, Hamburg, Germany) and stored at − 80° until analysis.

Nucleic acid extraction

1.5 mL of amniotic fluid was centrifuged for 15 min at 14,000×g. The obtained pellet was incubated with lysozyme (20 mg/ml) and Proteinase K (20 mg/ml) for 1 h at 37 °C. DNA was extracted following phenol–chloroform extraction protocol as described previously [27]. Additionally, blank extraction samples were performed to identify potential contamination during DNA extraction.

Quantitative real-time PCR

Quantitative real-time PCR of the 16S rRNA gene as proxy for bacterial load was performed on the ABI 7300 Cycler (Applied Biosystems, Germany) using the following reaction mixture: 12.5 µl 2 × Power SYBR Green master mix (Thermofisher Scientific, Germany), 5 pmol primers (FP 16S/RP 16S) [43], 0.5 µl 3% BSA and 2 µl DNA template in a total volume of 25 µl. PCR conditions were 10 min at 95 °C; 40 cycles of 45 s at 95 °C, 45 s at 58 °C, 45 s at 72 °C; 10 min 72 °C; 1 cycle of 15 s at 95 °C, 30 s at 60 °C, 15 s at 95 °C. All PCR products were checked on agarose gel. Amplification efficiency (calculated by Eff = [10^(− 1/slope) − 1]) was 90% with R2 of 0.994. The quantified gene copy numbers were normalized to 1 ml of amniotic fluid.

16S rRNA gene sequencing

Amplicon sequencing of the V3–V4 hypervariable region of the 16S rRNA gene was performed on a MiSeq Illumina instrument (MiSeq Reagent Kit v3 (600 Cycle); Illumina, San Diego, CA, USA) using the universal eubacterial primers 347F and 803R [28], extended with sequencing adapters to match Illumina indexing primers. To identify potential contamination during DNA extraction and amplification, both extraction and PCR no template control samples were performed. PCR was done using NEBNext high fidelity polymerase (New England Biolabs, Ipswich, USA) in a total volume of 25 µl (10 ng DNA template, 12.5 µl polymerase, 5 pmol of each primer). PCR conditions were 5 min at 98 °C; 35 cycles of 10 s at 98 °C, 30 s at 56 °C, 30 s at 72 °C; 10 min 72 °C. PCR products were purified using AMPure XP beads (Beckman Coulter Life Sciences, Indianapolis, USA) and quantified via PicoGreen assay. Subsequently, indexing PCR was performed using the Nextera XT Index Kit v2 (Illumina, Inc. San Diego, CA, USA) in a total volume of 25 µl (10 ng DNA template, 12.5 µl NEBNext high fidelity polymerase, 2.5 µl of each indexing primer) and the following PCR conditions: 30s at 98 °C; 8 cycles of 10 s at 98 °C, 30s at 55 °C, 30s at 72 °C; 5min 72 °C. Indexing PCR products were purified using AMPure XP beads, qualified and quantified via a Fragment Analyzer™ instrument (Advanced Analytical Technologies, Inc., Ankeny, USA) and pooled in an equimolar ratio of 4 nM.

Sequence data processing

Sequences were analysed using QIIME 2 software package release 2019.10.0 [29] with default parameters. FASTQ files were trimmed and merged with a minimum read length of 50 and minimum Phred score of 15 using AdapterRemoval [30]. Quality control was done using the QIIME 2 plugin DADA2 [31] by removing 20 bp n-terminally and truncating the reads at position 280 (forward) and 220 (reverse), respectively, resulting in unique amplicon sequence variants (ASV). For taxonomic analysis, Naive Bayes classifier was pre-trained on the SILVA_132_QIIME release 99% using the primer pair 347F / 803R. To exclude potential contamination, ASV occurring in extraction and PCR controls were removed from the dataset (33 and 2 ASV, respectively) (Additional file 1: Table S1). Furthermore, mitochondrial sequences and singletons were removed. Data normalization was done by subsampling to 5434 reads (lowest obtained read number). Rarefaction curves show saturation, indicating that this sampling depth was sufficient for further analysis. The sequence data obtained in this study are deposited in the short read archive of NCBI under accession number PRJNA817592.

Statistical analysis

All statistics were performed in R version 4.0.5 ( For qualitative parameters, Chi-Square Test and Fisher-Test were used as appropriate, with Benjamini–Hochberg p value adjustment for multiple pairwise comparison. Quantitative real-time PCR data were analyzed using Kruskal–Wallis test and Wilcoxon-rank sum test, respectively. Alpha diversity was calculated using species richness based on ASV number, Pielou’s evenness and Shannon diversity index. Beta diversity was analysed via unweighted and weighted UniFrac distance. For statistical purposes, Kruskal–Wallis test, Wilcoxon-rank sum test and PERMANOVA with Benjamini–Hochberg p value correction for multiple comparisons was used. To identify microbial taxa differing between the analysed groups, generalized linear models accounting for zero-inflation (R packages MASS and pscl) were used. For model validation, residual histograms, plots showing sample quantiles versus theoretical quantiles and plots showing residuals versus fitted values were checked for normal distribution and variance homogeneity of residuals [32]. Significant differences (p < 0.05) were calculated via ANOVA respective Wilcoxon and Kruskal–Wallis test using Benjamini–Hochberg p value adjustment for multiple pairwise comparisons. Additionally, methods for microbiome biomarker identification (LEfSe) and differential expression analysis (edgeR) were performed to verify the model results. Plots were created in R using ggplot2, ggpubr and metacoder.


Cohort characteristics

In line with the published literature, AF samples of n = 38 infants and n = 1 control rendered low reads bacterial 16S rRNA gene signal [33] leaving n = 70 preterm infants and n = 16 intact control pregnancies for the final analyses as presented in the flowchart in Fig. 1. Of importance, no systematic bias was introduced as the low read group and the cohort available for analysis did not differ for baseline characteristics (Table 1). In our analysis cohort, the median gestational age was 28 weeks (range 23–31) and the median birth weight 970 g (range 360–1800). Gender distribution was nearly balanced with 46% females and 76% received antenatal steroids. In line with the published literature, 56% of all preterm infants in our cohort fulfilled any BPD criterion and n = 14 still depended on respiratory support and/or supplemental oxygen at 36 weeks of gestation (Table 1 presenting median and quartiles) [34].

Table 1 Maternal and neonatal characteristics of the total study population

Microbial signatures in intact pregnancies and preterm deliveries

Comparing intact pregnancies and preterm deliveries, bacterial load in the AF was significantly higher whereas the alpha diversity was reduced in the preterm deliveries (Fig. 2A, B). Beta diversity showed clearly differences in 16S rRNA gene composition for preterm AF samples (Fig. 2C), with 16S rRNA gene sequences indicative for genera belonging to Proteobacteria, Firmicutes, Fusobacteria and Tenericutes differing significantly between AF samples of preterms and intact pregnancies. While ASV assigned to Acetobacter (Alphaproteobacteria) and Lactobacillus (Firmicutes) were significantly reduced in preterm AF samples, Pseudomonas (Gammaproteobacteria), Staphylococcus (Firmicutes), Fusobacterium (Fusobacteria) and Ureaplasma (Tenericutes) prevailed increased in relative 16S rRNA gene abundance (Fig. 2D). These results reveal for the first time a lower 16S rRNA gene load and higher richness in the AF of intact pregnancies compared to preterm deliveries.

Fig. 2
figure 2

Differences in 16S rRNA gene microbial abundance in amniotic fluid samples from preterm deliveries and intact pregnancies. A Box plots show the gene copy numbers of 16S rRNA gene per mL of AF samples of intact pregnancy (red) and preterm delivery (blue) and prevailed a significantly higher bacterial load in AF samples of preterm deliveries. B Box plots of different alpha diversity indices show that richness and Shannon diversity was significantly reduced in AF samples of preterm deliveries while evenness did not differ. Statistical analysis was performed using Wilcoxon Rank-Sum test. C NMDS plot of weighted Unifrac distances shows significant (p < 0.05) altered bacterial 16S rRNA gene community composition for AF samples of preterm deliveries. Statistical analysis was performed using PERMANOVA with Benjamini–Hochberg correction for multiple comparisons. D Heat tree including genera ≥ 0.1% of all reads. The tree shows the taxonomic information (domain to genus) and represents a comparison between preterm AF samples and samples of intact pregnancies. Coloured taxa are more abundant (based on log2-transformed ratio of median proportions) in the samples indicated (red: increased in preterm, blue: increased in intact). Significant changes (p < 0.05 based on Wilcoxon-Rank-Sum test with FDR correction) are marked with asterisks

Microbial signatures separated by BPD disease severity stages

To study the association of specific microbial 16S rRNA gene signatures in the AF at birth with the later development of BPD, infants were separated by the severity of BPD into no (n = 31), mild (n = 25) and moderate/severe (n = 14) BPD cases where the last category comprises the most severely affected infants with the highest risk for lifetime pulmonary sequelae. As expected, infants within the mild and moderate/severe BPD categories were more immature and had a lower gestational age at birth than infants within the no BPD category while they did not differ for maternal characteristics (Table 2). No significant differences in baseline parameters became evident between the two groups of infants fulfilling the BPD criterion besides lower birth weight (p = 0.045) (Table 2). Fitting to the other respiratory variables, the maximum fraction of oxygen required to meet the oxygen saturation target was significantly increased within the first 72 h in infants with moderate/severe BPD reflecting the more compromised gas exchange at birth (Table 2).

Table 2 Neonatal characteristics of patients included in the analyses separated for the severity of BPD

Total bacterial 16S rRNA gene load, gene alpha and beta diversity did not differ between the three BPD categories (except partly increased evenness of mild BPD) arguing for comparable absolute abundance and diversity pattern of bacterial taxa (Fig. 3A–C). However, at the level of single ASV, 16S rRNA sequences indicative for Enterococcus (Firmicutes) and Ureaplasma (Tenericutes) were present in higher relative abundance in the mild BPD compared to the no BPD group while 16S rRNA genes of Gardnerella (Actinobacteria) and members of the Escherichia-Shigella complex (Proteobacteria) displayed a significantly lower relative abundance (Fig. 3D). When infants with moderate/severe BPD were analysed in comparison to the mild BPD group, relative abundance of Proteobacteria within the Escherichia-Shigella group was significantly higher in the moderate/severe BPD cases while no significant variations were present for the further ASV that differed between infants with no and mild BPD (Fig. 3D). These species did as well not differ significantly in relative abundance between the no and the moderate/severe category. LEfSe analysis revealed LDA scores of 4.1–5.2 (p < 0.05), confirming those taxa being relevant biomarkers for the respective groups.

Fig. 3
figure 3

Disparities in relative but not absolute 16S rRNA gene microbial abundance in amniotic fluid samples from preterm deliveries depending on the BPD disease severity. A AF samples from preterm deliveries from Fig. 2 were separated for the BPD disease severity into no (red), mild (green) and moderate/severe (m/s, blue) BPD. Box plots show the gene copy numbers of 16S rRNA gene per ml of AF samples and prevailed no significant difference in total bacterial load. Statistical analysis was performed using Kruskal–Wallis and Wilcoxon Rank-Sum test with Benjamini–Hochberg correction for multiple comparisons, respectively. B Box plots of different alpha (diversity indices show that richness and Shannon diversity did not differ significantly between the samples of different disease severities while evenness of mild BPD was increased. Statistical analysis was performed using Kruskal–Wallis and Wilcoxon Rank-Sum test with Benjamini–Hochberg correction for multiple comparisons, respectively. C NMDS plot of weighted Unifrac distances shows no significant (p < 0.05) difference in bacterial 16S rRNA gene community composition for samples of different disease severities. Statistical analysis was performed using PERMANOVA with Benjamini–Hochberg correction for multiple comparisons. D Heat tree including genera ≥ 0.1% of all reads. The labelled tree on the lower left shows the taxonomic information (domain to genus) and is the key for the unlabelled smaller trees. Smaller trees represent a pairwise comparison between no, mild and moderate/severe (m/s) BPD. Coloured taxa are more abundant (based on log2-transformed ratio of median proportions) in the samples indicated in the respective tree subtitle colour. Significant changes (p < 0.05 based on Wilcoxon-Rank-Sum test with FDR correction) are marked with green rectangles


Distinct 16S rRNA genes composition in intact pregnancies and preterm deliveries

Our study demonstrates that loads as well as composition of 16S rRNA genes differ between prematurely born infants immediately before their deliveries and intact pregnancies of comparable gestational age. The absolute 16S rRNA gene abundance was increased in preterm deliveries while alpha diversity was decreased and beta diversity showed a significantly altered microbial community composition, arguing for important differences in microbial 16S rRNA gene abundance and composition. Our data are congruent with previous publications, where a higher vaginal bacterial load and a Lactobacillus poor bacterial community with the predominance of Ureaplasma and Fusobacterium were associated with preterm birth [35, 36]. The results of our AF samples indicate similar changes in AF bacterial DNA signatures immediately before premature delivery but it might be possible that further taxa including Acetobacter, Staphylococcus and Pseudomonas need to be considered. Therefore our association study might be suited to identify infants at risk for preterm delivery already before birth.

Separation of BPD disease severity stages by 16S rRNA genes at birth

Of highest relevance, our study provides for the first-time conclusive evidence that distinct AF bacterial 16S rRNA gene signatures immediately before delivery separate infants for different BPD disease severities. In contrast to the results for premature delivery, here neither the total 16S rRNA gene abundance nor the alpha and beta diversity segregated BPD categories but differences in ASV linked to the relative abundance of specific genera were observed. In AF samples of infants with mild BPD, ASV assigned to Ureaplasma and Enterococcus species were increased compared to no BPD samples. While comparable trends were observed in infants with moderate/severe BPD, ASV linked to the Escherichia-Shigella complex significantly differ between mild and moderate/severe BPD samples with Escherichia-Shigella being increased in the latter group. By reporting an association of prenatal AF 16S rRNA gene-based microbial community changes and BPD severity, the presented data are unique and in line with recent association studies on the consequences of aberrant microbial lung colonization after birth or postnatal infections [37, 38]. Two recent publications provided conclusive evidence, that changes in the vaginal microbiome and in the airway microbiome drive innate immune activation towards an inflammatory status of the metabolome [37, 39]. Although not investigated here, a similar connection between a dysbiotic structure of the AF microbiome and metabolomic alterations might be assumed and explain the prenatal inflammatory damage to the immature lung.

Results of the study in the context of published preclinical studies

Our clinical associations are solidified by a well-founded level of preclinical evidence where causal relations between microbial actions and BPD development were studied on an unprecedented level: Both prenatal and postnatal LPS exposure derived from Escherichia coli caused the typical severe hallmarks of BPD. The LPS effects on the immature lung were more pronounced with the increased dosage of LPS. Furthermore, lung injury could be worsened by synergistic effects of LPS and hyperoxia [13, 40]. Further elaboration of the model separated effects on the lung that were induced both by LPS or hyperoxia alone including exhaustion of antioxidative capacities while cell proliferation and vascular development and function were particularly affected during the combined exposure. Probably most important, the authors were able to ascribe adaptive immune responses to low dosages of LPS, and here subsequent hyperoxia did not aggravate lung injury. In contrast, high dosages of LPS plus hyperoxia caused the worst BPD phenotype that mimics the most frequent clinical scenario of excessive intraamniotic infection and subsequent need for postnatal supplemental oxygen [41]. Lastly and in line with the previous data, the germ-free conditions in gnotobiotic mice attenuated the pathologies of BPD provoked by hyperoxia, particularly inflammation that were traced back to the lack of preceding shaping of the host immune system [42]. Vice versa, prenatal antibiotic therapy disrupted the intestinal bacterial colonization of the offspring with commensals resulting in a more severe lung injury in the newborn hyperoxia mouse model [43]. When taking all these results into consideration, our results on the disparate microbial 16S rRNA gene signatures in infants with mild and moderate/severe BPD are not surprising and cannot be traced back solely to disparities in clinical baseline parameters including gestational age and birth weight. Here, the AF samples from intact pregnancies of equivalent gestational age as the moderate/severe BPD group substantiate our conclusions. Ureaplasma and most likely as well Enterococcus induce a restricted pro-inflammatory response or even an immunotolerance phenotype in animal studies [15, 40]. A similar phenomenon although not that pronounced was observed in the immune cells of premature infants [44]. In contrast, the overwhelming effects of bacteria of the Escherichia-Shigella complex with the release of LPS can induce a pronounced inflammatory response and lung injury during severe intraamniotic infection [41]. Interestingly, bacteria of the Escherichia-Shigella complex were also highly abundant in the group of infants without BPD. Here, the five weeks older gestational age at birth and the concomitant tremendous differences in maturity of the lung and the immune system may account for the absent predisposition for severe lung injury. It needs to be mentioned that Paracoccus was present in about half of samples studied. Although mainly present in diverse natural environments like soil, pristine or marine sediments, members of this genus could have been isolated from human skin or eye, too, and were implicated in opportunistic infections of humans [45, 46]. However, its relative proportion was < 0.5% except for 3 samples and is therefore not of relevance to our results.

Limitations of the study

We need to acknowledge several limitations of our study. Most importantly, we did not investigate living bacteria but executed 16S rRNA gene-based sequencing in our study. Parallel microbial cultivation was not feasible as all pregnant women with preterm labour and premature rupture of membranes received antibiotics before delivery. Repeated sampling was not feasible as decision to preterm delivery and birth intervals were mostly short and additional prenatal sample collection by amniotic puncture not justifiable for ethical reasons. Thereby, we are unable to allocate the origin of the 16S rRNA genes to transmigration of living bacteria that are rapidly killed or leakage of bacterial particles as described in septic patients and healthy controls before [47,48,49,50]. Particles might derive from the gut, oropharynx or other microbial niches of the body and transmission might be propagated by maternal local or systemic inflammation in our setting [35, 51, 52]. Taking all these considerations into account, the bacterial 16S rRNA gene transmission might be caused by maternal inflammation and accompany that of inflammatory particles into the amniotic cavity. Due to the nature of prematurity, some AF samples rendered low reads of bacterial 16S rRNA genes using appropriate sequencing techniques [50] and did not meet the threshold of > 5434 bacterial reads (which was taken to cover the microbial diversity within the sample based on rarefaction curve analysis) [33] preventing analysis of all samples collected. A systematic bias was excluded as cohort characteristics did not differ between excluded and included patients (Table 1). Furthermore, we cannot completely exclude sampling contamination. But such selective contamination of the sampling environment of the abdominal skin within one isolated surgery room used are unlikely. These limitations need to be carefully addressed in subsequent studies to clarify the presence of a livid intra-amniotic microbiome and its origin. The restriction in sample size derives from the single centre approach and is too low to further separate beyond the severity stage of BPD. Therefore, we were unable to respect the causes of preterm delivery, particularly for intraamniotic infection. This needs to be addressed within future larger scaled prospective multicentre studies. Furthermore, the number of samples available for analyses from intact pregnancies (n = 16) compared to the AF samples of preterm deliveries (n = 70) and the fraction of preterm infants developing moderate/severe BPD (n = 14 of n = 70) was limited.

Our cohort displays heterogenous baseline characteristics between infants without and with BPD. The high disparities for the gestational age at birth derive from the nature of BPD as immaturity of the lung constitutes the main risk factor. But it needs to be stated that nearly all very immature infants < 28 weeks of gestation develop more or less pronounced restrictions in lung function. Therefore, the focus on the comparison of mild and moderate/severe BPD cases is of higher relevance. The addition of AF samples from control pregnancies might expand our knowledge of our study on the differences in 16S rRNA genes [53]. Our data argue towards focussing on disease severities of BPD instead of comparing the yes/no criterion that compares infants of completely different prerequisites. From all these considerations it becomes clear that our data are not suited to specify the causality between Ureaplasma or Enterococcus colonization and the outcome of mild restrictions in lung function, but a more severe BPD disease severity in connection to Escherichia-Shigella species.


Overall, our results provide first-time evidence that microbial 16S rRNA gene signatures in the AF at birth differ between BPD disease categories. These data have the potential to guide future research directions towards the prenatal origins of BPD. The clear distinction of 16S rRNA gene microbiome signatures between intact pregnancies and preterm deliveries and the association with BPD severity categories in a limited number of participants argues towards a relevant impact on BPD. Future studies need to specify whether microbiota shaping therapies or immunomodulatory properties constitute appropriate therapeutic approaches or whether both aspects need to be considered as it is probably the case for the highly beneficial effects of breast milk to prevent BPD [54, 55]. A detailed knowledge on the mutual interference of microbiome structures and the immune system in the clinical arena can deliver the basis to bridge the insights from preclinical studies to the clinics and to direct future research strategies to prevent BPD. Our results argue that it might be necessary to initiate such a strategy before delivery. Finally, our results of distinct 16S rRNA gene-based microbiota signatures in different BPD disease categories demand such evaluation within further lung diseases and beyond and might give rise to a novel mechanistic concept of drivers of disease severity.

Availability of data and materials

The data that support the findings of this study are available on reasonable request from the corresponding author.



Bronchopulmonary dysplasia


Amniotic fluid




Small for gestational age


Amplicon sequence variant


  1. Northway WH, Rosan RC, Porter DY. Pulmonary disease following respirator therapy of hyaline-membrane disease. Bronchopulmonary dysplasia. N Engl J Med. 1967;276:357–68.

    Article  PubMed  Google Scholar 

  2. Bårdsen T, Røksund OD, Benestad MR, Hufthammer KO, Clemm HH, Mikalsen IB, et al. Tracking of lung function from 10 to 35 years after being born extremely preterm or with extremely low birth weight. Thorax. 2022;77:790–8.

    Article  PubMed  Google Scholar 

  3. Lee DMX, Tan AKS, Ng YPM, Amin Z. Quality of life of patients and caregivers affected by bronchopulmonary dysplasia: a systematic review. Qual Life Res. 2023;32:1859–69.

    Article  PubMed  Google Scholar 

  4. Laughon MM, Langer JC, Bose CL, Smith PB, Ambalavanan N, Kennedy KA, et al. Prediction of bronchopulmonary dysplasia by postnatal age in extremely premature infants. Am J Respir Crit Care Med. 2011;183:1715–22.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Schmidt B, Roberts R, Millar D, Kirpalani H. Evidence-based neonatal drug therapy for prevention of bronchopulmonary dysplasia in very-low-birth-weight infants. Neonatology. 2008;93:284–7.

    Article  PubMed  Google Scholar 

  6. Holzfurtner L, Shahzad T, Dong Y, Rekers L, Selting A, Staude B, et al. When inflammation meets lung development—an update on the pathogenesis of bronchopulmonary dysplasia. Mol Cell Pediatr. 2022;9:1–12.

    Article  Google Scholar 

  7. Pammi M, Lal CV, Wagner BD, Mourani PM, Lohmann P, Luna RA, et al. Airway microbiome and development of bronchopulmonary dysplasia in preterm infants: a systematic review. J Pediatr. 2019;204:126-133.e2.

    Article  PubMed  Google Scholar 

  8. Lauer T, Behnke J, Oehmke F, Baecker J, Gentil K, Chakraborty T, et al. Bacterial colonization within the first six weeks of life and pulmonary outcome in preterm infants <1000 g. J Clin Med. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Dong Y, Rivetti S, Lingampally A, Tacke S, Kojonazarov B, Bellusci S, Ehrhardt H. Insights into the black box of intra-amniotic infection and its impact on the premature lung: from clinical and preclinical perspectives. Int J Mol Sci. 2022.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Lowe J, Watkins WJ, Edwards MO, Spiller OB, Jacqz-Aigrain E, Kotecha SJ, Kotecha S. Association between pulmonary ureaplasma colonization and bronchopulmonary dysplasia in preterm infants: updated systematic review and meta-analysis. Pediatr Infect Dis J. 2014;33:697–702.

    Article  PubMed  Google Scholar 

  11. Hartling L, Liang Y, Lacaze-Masmonteil T. Chorioamnionitis as a risk factor for bronchopulmonary dysplasia: a systematic review and meta-analysis. Arch Dis Child Fetal Neonatal Ed. 2012;97:F8–17.

    Article  PubMed  Google Scholar 

  12. Ehrhardt H, Desplanches T, van Heijst A, Toome L, Fenton A, Torchin H, et al. Mode of delivery and incidence of bronchopulmonary dysplasia: results from the population-based EPICE cohort. Neonatology. 2022;119:464–73.

    Article  CAS  PubMed  Google Scholar 

  13. Stranik J, Kacerovsky M, Vescicik P, Faist T, Jacobsson B, Musilova I. A rodent model of intra-amniotic inflammation/infection, induced by the administration of inflammatory agent in a gestational sac, associated with preterm delivery: a systematic review. J Matern Fetal Neonatal Med. 2022;35:1592–600.

    Article  CAS  PubMed  Google Scholar 

  14. Polglase GR, Dalton RGB, Nitsos I, Knox CL, Pillow JJ, Jobe AH, et al. Pulmonary vascular and alveolar development in preterm lambs chronically colonized with Ureaplasma parvum. Am J Physiol Lung Cell Mol Physiol. 2010;299:L232–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kallapur SG, Kramer BW, Knox CL, Berry CA, Collins JJP, Kemp MW, et al. Chronic fetal exposure to Ureaplasma parvum suppresses innate immune responses in sheep. J Immunol. 2011;187:2688–95.

    Article  CAS  PubMed  Google Scholar 

  16. Kuypers E, Collins JJP, Kramer BW, Ofman G, Nitsos I, Pillow JJ, et al. Intra-amniotic LPS and antenatal betamethasone: inflammation and maturation in preterm lamb lungs. Am J Physiol Lung Cell Mol Physiol. 2012;302:L380–9.

    Article  CAS  PubMed  Google Scholar 

  17. Stiemsma LT, Michels KB. The role of the microbiome in the developmental origins of health and disease. Pediatrics. 2018.

    Article  PubMed  Google Scholar 

  18. Blaser MJ, Devkota S, McCoy KD, Relman DA, Yassour M, Young VB. Lessons learned from the prenatal microbiome controversy. Microbiome. 2021;9:8.

    Article  PubMed  PubMed Central  Google Scholar 

  19. de Goffau MC, Lager S, Sovio U, Gaccioli F, Cook E, Peacock SJ, et al. Human placenta has no microbiome but can contain potential pathogens. Nature. 2019;572:329–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Leiby JS, McCormick K, Sherrill-Mix S, Clarke EL, Kessler LR, Taylor LJ, et al. Lack of detection of a human placenta microbiome in samples from preterm and term deliveries. Microbiome. 2018;6:196.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Sterpu I, Fransson E, Hugerth LW, Du J, Pereira M, Cheng L, et al. No evidence for a placental microbiome in human pregnancies at term. Am J Obstet Gynecol. 2021;224:296.e1-296.e23.

    Article  CAS  PubMed  Google Scholar 

  22. Behnke J, Estreich V, Oehmke F, Zimmer K-P, Windhorst A, Ehrhardt H. Compatibility of rapid enteral feeding advances and noninvasive ventilation in preterm infants—an observational study. Pediatr Pulmonol. 2022;57:1117–26.

    Article  PubMed  Google Scholar 

  23. Voigt M, Schneider KT, Jährig K. Analyse des Geburtengutes des Jahrgangs 1992 der Bundesrepublik Deutschland. Teil 1: Neue Perzentilwerte für die Körpermasse von Neugeborenen. [Analysis of a 1992 birth sample in Germany. 1: New percentile values of the body weight of newborn infants]. Geburtshilfe Frauenheilkd. 1996;56:550–8.

    Article  CAS  PubMed  Google Scholar 

  24. Thiess T, Lauer T, Woesler A, Neusius J, Stehle S, Zimmer K-P, et al. Correlation of early nutritional supply and development of bronchopulmonary dysplasia in preterm infants <1,000 g. Front Pediatr. 2021;9: 741365.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wilkinson DJ, Andersen CC, Smith K, Holberton J. Pharyngeal pressure with high-flow nasal cannulae in premature infants. J Perinatol. 2008;28:42–7.

    Article  CAS  PubMed  Google Scholar 

  26. Walsh M, Engle W, Laptook A, Kazzi SNJ, Buchter S, Rasmussen M, Yao Q. Oxygen delivery through nasal cannulae to preterm infants: can practice be improved? Pediatrics. 2005;116:857–61.

    Article  PubMed  Google Scholar 

  27. Lueders T, Manefield M, Friedrich MW. Enhanced sensitivity of DNA- and rRNA-based stable isotope probing by fractionation and quantitative analysis of isopycnic centrifugation gradients. Environ Microbiol. 2004;6:73–8.

    Article  CAS  PubMed  Google Scholar 

  28. Nossa CW, Oberdorf WE, Yang L, Aas JA, Paster BJ, Desantis TZ, et al. Design of 16S rRNA gene primers for 454 pyrosequencing of the human foregut microbiome. World J Gastroenterol. 2010;16:4135–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Schubert M, Lindgreen S, Orlando L. AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res Notes. 2016;9:88.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Rosenblad AJJ. Faraway: extending the linear model with R: generalized linear, mixed effects and nonparametric regression models. Comput Stat. 2009;24:369–70.

    Article  Google Scholar 

  33. Gallacher D, Mitchell E, Alber D, Wach R, Klein N, Marchesi JR, Kotecha S. Dissimilarity of the gut–lung axis and dysbiosis of the lower airways in ventilated preterm infants. Eur Respir J. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Shah PS, Lui K, Sjörs G, Mirea L, Reichman B, Adams M, et al. Neonatal outcomes of very low birth weight and very preterm neonates: an international comparison. J Pediatr. 2016;177:144-152.e6.

    Article  PubMed  Google Scholar 

  35. Staude B, Oehmke F, Lauer T, Behnke J, Göpel W, Schloter M, et al. The microbiome and preterm birth: a change in paradigm with profound implications for pathophysiologic concepts and novel therapeutic strategies. Biomed Res Int. 2018;2018:7218187.

    Article  PubMed  PubMed Central  Google Scholar 

  36. DiGiulio DB, Callahan BJ, McMurdie PJ, Costello EK, Lyell DJ, Robaczewska A, et al. Temporal and spatial variation of the human microbiota during pregnancy. Proc Natl Acad Sci USA. 2015;112:11060–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lal CV, Kandasamy J, Dolma K, Ramani M, Kumar R, Wilson L, et al. Early airway microbial metagenomic and metabolomic signatures are associated with development of severe bronchopulmonary dysplasia. Am J Physiol Lung Cell Mol Physiol. 2018;315:L810–5.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Jung E, Lee BS. Late-onset sepsis as a risk factor for bronchopulmonary dysplasia in extremely low birth weight infants: a nationwide cohort study. Sci Rep. 2019;9:15448.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Pruski P, Correia GDS, Lewis HV, Capuccini K, Inglese P, Chan D, et al. Direct on-swab metabolic profiling of vaginal microbiome host interactions during pregnancy and preterm birth. Nat Commun. 2021;12:5967.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Shrestha AK, Bettini ML, Menon RT, Gopal VYN, Huang S, Edwards DP, et al. Consequences of early postnatal lipopolysaccharide exposure on developing lungs in mice. Am J Physiol Lung Cell Mol Physiol. 2019;316:L229–44.

    Article  CAS  PubMed  Google Scholar 

  41. Shrestha AK, Menon RT, El-Saie A, Barrios R, Reynolds C, Shivanna B. Interactive and independent effects of early lipopolysaccharide and hyperoxia exposure on developing murine lungs. Am J Physiol Lung Cell Mol Physiol. 2020;319:L981–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Dolma K, Freeman AE, Rezonzew G, Payne GA, Xu X, Jilling T, et al. Effects of hyperoxia on alveolar and pulmonary vascular development in germ-free mice. Am J Physiol Lung Cell Mol Physiol. 2019;318:L421–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Willis KA, Siefker DT, Aziz MM, White CT, Mussarat N, Gomes CK, et al. Perinatal maternal antibiotic exposure augments lung injury in offspring in experimental bronchopulmonary dysplasia. Am J Physiol Lung Cell Mol Physiol. 2020;318:L407–18.

    Article  CAS  PubMed  Google Scholar 

  44. Azizia M, Lloyd J, Allen M, Klein N, Peebles D. Immune status in very preterm neonates. Pediatrics. 2012;129:e967–74.

    Article  PubMed  Google Scholar 

  45. Dong K, Pu J, Yang J, Zhou G, Gao Y, Kang Z, et al. Whole-genome sequencing of Paracoccus species isolated from the healthy human eye and description of Paracoccus shanxieyensis sp. nov. Int J Syst Evol Microbiol. 2023;73:5838.

    Article  CAS  Google Scholar 

  46. Lasek R, Szuplewska M, Mitura M, Decewicz P, Chmielowska C, Pawłot A, et al. Genome structure of the opportunistic pathogen paracoccus yeei (Alphaproteobacteria) and identification of putative virulence factors. Front Microbiol. 2018;9:2553.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Gosiewski T, Ludwig-Galezowska AH, Huminska K, Sroka-Oleksiak A, Radkowski P, Salamon D, et al. Comprehensive detection and identification of bacterial DNA in the blood of patients with sepsis and healthy volunteers using next-generation sequencing method—the observation of DNAemia. Eur J Clin Microbiol Infect Dis. 2017;36:329–36.

    Article  CAS  PubMed  Google Scholar 

  48. Wang X, Buhimschi CS, Temoin S, Bhandari V, Han YW, Buhimschi IA. Comparative microbial analysis of paired amniotic fluid and cord blood from pregnancies complicated by preterm birth and early-onset neonatal sepsis. PLoS ONE. 2013;8: e56131.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Lim ES, Rodriguez C, Holtz LR. Amniotic fluid from healthy term pregnancies does not harbor a detectable microbial community. Microbiome. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Urushiyama D, Suda W, Ohnishi E, Araki R, Kiyoshima C, Kurakazu M, et al. Microbiome profile of the amniotic fluid as a predictive biomarker of perinatal outcome. Sci Rep. 2017;7:12171.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Cobo T, Vergara A, Collado MC, Casals-Pascual C, Herreros E, Bosch J, et al. Characterization of vaginal microbiota in women with preterm labor with intra-amniotic inflammation. Sci Rep. 2019;9:18963.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Fettweis JM, Serrano MG, Brooks JP, Edwards DJ, Girerd PH, Parikh HI, et al. The vaginal microbiome and preterm birth. Nat Med. 2019;25:1012–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Fawke J, Lum S, Kirkby J, Hennessy E, Marlow N, Rowell V, et al. Lung function and respiratory symptoms at 11 years in children born extremely preterm: the EPICure study. Am J Respir Crit Care Med. 2010;182:237–45.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Villamor-Martínez E, Pierro M, Cavallaro G, Mosca F, Kramer B, Villamor E. Probiotic supplementation in preterm infants does not affect the risk of bronchopulmonary dysplasia: a meta-analysis of randomized controlled trials. Nutrients. 2017.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Huang J, Zhang L, Tang J, Shi J, Qu Y, Xiong T, Mu D. Human milk as a protective factor for bronchopulmonary dysplasia: a systematic review and meta-analysis. Arch Dis Child Fetal Neonatal Ed. 2019;104:F128–36.

    Article  PubMed  Google Scholar 

Download references


We thank the parents of all infants included into the study for their consent to participate.


Open Access funding enabled and organized by Projekt DEAL. This study was supported by Clinical Research Unit KFO309/2. The sponsor had no role in the design of the study, data collection and analyses, execution, interpretation or writing of the study.

Author information

Authors and Affiliations



Conceptualization, BS, SG, MS and HE; Methodology, SG, MS and HE; Sample collection, BS, TF, FO, TK and HE; Sample processing, SG, SK and MS; Sequencing, SG, SK and MS; Collection of patient data, BS, TF, FO and TK; Visualization, SG; Writing of the original draft, BS, SG, MS and HE; Writing, reviewing and editing, all authors; all authors read and approved the final manuscript.

Corresponding author

Correspondence to Harald Ehrhardt.

Ethics declarations

Ethics approval and consent to participate

The study was conducted following the rules of the Declaration of Helsinki, was approved by the ethics committee of the Justus-Liebig-University of Gießen (Az 135/12) and registered at DRKS (DRKS00004600). Written informed consent was obtained from the parents of preterm infants and pregnant women intended for prenatal interventions after provision of oral and written information.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. List of Amplicon sequence variants (ASV) detected in blanks. Those ASV were considered as potential contaminants and removed from AF dataset before analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Staude, B., Gschwendtner, S., Frodermann, T. et al. Microbial signatures in amniotic fluid at preterm birth and association with bronchopulmonary dysplasia. Respir Res 24, 248 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: