Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Data Mining of Lung Microbiota in Cystic Fibrosis Patients

  • Jianguo Li ,

    Contributed equally to this work with: Jianguo Li, Chunyan Hao, Lili Ren

    Affiliation Modern Research Center for Traditional Chinese Medicine, Shanxi University, Taiyuan 030006, China

  • Chunyan Hao ,

    Contributed equally to this work with: Jianguo Li, Chunyan Hao, Lili Ren

    Affiliation College of Chemical & Biological Engineering, Taiyuan University of Science & Technology, Taiyuan 030021, China

  • Lili Ren ,

    Contributed equally to this work with: Jianguo Li, Chunyan Hao, Lili Ren

    Affiliation MOH Key Laboratory of System Pathogen Biology and Christophe Mérieux Laboratory, IPB, CAMS-Foundation Mérieux, Institute of Pathogen Biology (IPB), Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100730, China

  • Yan Xiao,

    Affiliation MOH Key Laboratory of System Pathogen Biology and Christophe Mérieux Laboratory, IPB, CAMS-Foundation Mérieux, Institute of Pathogen Biology (IPB), Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100730, China

  • Jianwei Wang,

    Affiliation MOH Key Laboratory of System Pathogen Biology and Christophe Mérieux Laboratory, IPB, CAMS-Foundation Mérieux, Institute of Pathogen Biology (IPB), Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100730, China

  • Xuemei Qin

    qinxm001@126.com

    Affiliation Modern Research Center for Traditional Chinese Medicine, Shanxi University, Taiyuan 030006, China

Abstract

The major therapeutic strategy used to treat exacerbated cystic fibrosis (CF) is antibiotic treatment. As this approach easily generates antibiotic-resistant strains of opportunistic bacteria, optimized antibiotic therapies are required to effectively control chronic and recurrent bacterial infections in CF patients. A promising future for the proper use of antibiotics is the management of lung microbiota. However, the impact of antibiotic treatments on CF microbiota and vice versa is not fully understood. This study analyzed 718 sputum samples from 18 previous studies to identify differences between CF and uninfected lung microbiota and to evaluate the effects of antibiotic treatments on exacerbated CF microbiota. A reference-based OTU (operational taxonomic unit) picking method was used to combine analyses of data generated using different protocols and platforms. Findings show that CF microbiota had greater richness and lower diversity in the community structure than uninfected control (NIC) microbiota. Specifically, CF microbiota showed higher levels of opportunistic bacteria and dramatically lower levels of commensal bacteria. Antibiotic treatment affected exacerbated CF microbiota notably but only transiently during the treatment period. Limited decrease of the dominant opportunistic bacteria and a dramatic decrease of commensal bacteria were observed during the antibiotic treatment for CF exacerbation. Simultaneously, low abundance opportunistic bacteria were thriving after the antibiotic treatment. The inefficiency of the current antibiotic treatment against major opportunistic bacteria and the detrimental effects on commensal bacteria indicate that the current empiric antibiotic treatment on CF exacerbation should be reevaluated and optimized.

Introduction

Cystic fibrosis (CF), as an inherited disease of the secretory glands, is caused by defective mutations of the CF transmembrane conductance regulator (CFTR) gene [13]. Thickened mucus secretions resulting from the defective CFTR enable the development and persistence of pulmonary bacterial infections [4, 5]. These chronic pulmonary infections, which are mainly caused by opportunistic bacteria such as Pseudomonas aeruginosa and Staphylococcus aureus, induce intense bronchial neutrophilic inflammation. As a result of chronic polymicrobial infection and inflammation, CF is a devastating life-threatening disease with the majority of mortalities due to respiratory failure [6, 7]. Over the course of CF lung disease two categories have been identified according to disease severity: exacerbated and clinically stable. Pulmonary exacerbation is the hallmark of CF and is characterized by aggravated pulmonary symptoms and decreased pulmonary function. A clinically stable condition is the stationary phase before or after pulmonary exacerbation.

Unlike the traditional view of a sterile lung, accumulating evidence shows that a microbial community resides in the lung tissue; this is referred to as the lung microbiota [811]. Healthy lung microbiota predominately consist of commensal bacteria, including Streptococcus, Veillonella, Prevotella, and Actinomyces. [8, 1214]. Maintenance of a healthy lung microbiota is important for reducing the risk of respiratory tract diseases [12, 15, 16]. Decreased diversity and increased amounts of opportunistic bacteria are notable features of CF lung microbiota compared to uninfected control (NIC) microbiota [1720]. A substantial reduction in the diversity of airway microbiota paired with stabilized microbial density is associated with CF lung disease progression [8, 17, 2123].

Chronic and recurrent colonization of opportunistic bacteria causes devastating lung damage and respiratory decline. Intravenous (IV) injection of broad-spectrum antibiotics is the leading therapy for treating pulmonary exacerbation. A typical IV-antibiotic treatment usually requires a duration of two weeks and is inefficient at controlling infections because of the easily generated antibiotic-resistant bacterial strains [2429]. Thus, effective control mechanisms against opportunistic bacterial infections are needed without generating antibiotic-resistant strains to optimize therapies [22, 3034].

Management of lung microbiota provides a promising alternative for the treatment of CF with pulmonary exacerbation [8, 10]. Probiotic commensal bacteria of NIC microbiota have substantial and continuous effects on human health and physiological development, including in maturing the immune system and in preventing pathogen invasion [35]. By regulating pathogen activity and enhancing the natural immune system using probiotic supplementary therapies and antibiotic treatments that precisely target opportunistic bacteria, the growth of antibiotic-resistant bacteria can be prevented and the toxicity of these bacteria to NIC microbiota can be reduced [36]. However, several key questions must be answered before microbiota profiling can inform clinical practice. First, how many opportunistic bacteria are there in CF lung microbiota? Second, aside from opportunistic bacteria what are the differences between CF and NIC microbiota? Third, how do antibiotic treatments influence opportunistic and commensal bacteria? Numerous studies have attempted to answer these questions, but little instructive information has been obtained. This is partially due to the small sample size enrolled, which has limited the ability to investigate potential correlations between lung microbiota and CF [12].

This study analyzed the largest number of airway microbiota from CF patients and uninfected controls from published studies to evaluate the differences between uninfected and CF-derived microbiota. The efficacy of the antibiotic treatments in regulating the composition of lung microbiota in CF patients was also evaluated. In general, intravenous antibiotics treatment is the most popular strategies for CF exacerbation, while oral intake of antibiotics is commonly used to control seizures. The antibiotic treatment mentioned in this study refers to the large overdose intravenous therapy for CF exacerbation. Overall, our findings show that CF microbiota were richer and had a lower diversity in the community structure than NIC microbiota The obtained data provide insights necessary to evaluate the current therapeutic strategy.

Materials and methods

1. Definitions

The items used in this study were the same as those defined by Rabin et al [8] with some modifications. Briefly, enrolled CF patients were divided into two groups according to disease stages: baseline CF and CF with pulmonary exacerbation (defined as a hospitalization for respiratory symptoms requiring large-dose antibiotic therapy.) The pulmonary exacerbation stage was divided into three subgroups according to the stages of antibiotic administration: BeforeTreat, DuringTreat, and AfterTreat. NIC (uninfected control) was defined an individual without acute or chronic infections such as viral infections, tuberculosis, chronic obstructive pulmonary disease, bronchiectasis, or idiopathic pulmonary fibrosis. Antibiotic therapy or antibiotic treatment in this study specifically referred to the administration of antibiotics necessary to treat CF with pulmonary exacerbation, excluding the use of antibiotics for chronic maintenance.

2. Data collection and preprocessing

High-throughput sequencing data of 16s rRNA were collected from public databases (Table 1). Bacterial 16s rRNA genes were PCR amplified and sequenced on different sequencing platforms across the studies (Roche 454 or Illumina Hiseq/Miseq). Primers and targeted regions are listed in Table 1. Raw sequences were filtered and processed using the default parameters of the NGS QC Toolkit version 2.3 [37]. Briefly, for data generated by the Roche 454 (pyrosequencing) platform, sequences that were outside of 200–1000 nucleotides in length, with more than six ambiguous bases, contained homopolymer regions (> 6 bp), had at least two mismatches in the primer, or could not be assigned to a sample by the corresponding barcode were excluded. For data generated by the Illumina (Hiseq and Miseq) platform, paired-end reads were first merged if possible. Reads were then truncated if more than three consecutive low-quality base calls were found, and reads were excluded if ambiguous bases were found after quality trimming.

thumbnail
Table 1. Information of the data exploited in this study.

https://doi.org/10.1371/journal.pone.0164510.t001

3. OTU picking

To compare sequences from different regions of the 16s rRNA gene, a reference mapping protocol [3839] was used to group the sequences into operational taxonomic units (OTUs). Sequences derived from the same bacteria were assumed to match the same reference sequence regardless of the region or length of the targeted 16s rRNA gene. Specifically, the Greengenes database (May 2013 version) [40] of near full-length 16s rRNA was used as the reference database. Close reference OTU picking module in QIIME (version 1.8.0) was used to assign the OTUs [41], and a 97% minimum pairwise nucleotide sequence identity threshold was employed to align sequences in UCLUST [42]. Sequences with identity percentages lower than 97% to any of the reference sequences in the Greengenes database were excluded from further analysis. Moreover, samples with a proportion of assignable sequences less than 85% were excluded. Samples with less than 1000 sequences after quality filtering and OTU assignment were also discarded.

4. Alpha diversity measurements

Shared OTUs among different sample groups were plotted in Venn diagrams using Venny software (http://bioinfogp.cnb.csic.es/tools/venny/index.html). Alpha diversity measurements of community richness and diversity were analyzed by estimators of Chao1 and Shannon indices, respectively. Default parameters were applied and boxplots were created using the “Phyloseq” in R package [43]. OTUs of each sample across all studies were then grouped by different taxonomic levels using the workflow script of QIIME according to the 454 overview tutorial (http://qiime.org/tutorials/tutorial.html).

5. Beta diversity measurements

To compare the bacterial composition between different groups, unweighted UniFrac analyses of bray-curtis dissimilarity were performed. For this purpose, one thousand sequences were randomly selected from each sample. All analyses were carried out using QIIME.

The 16s rRNA sequencing data and corresponding metadata used in this study are available in publically accessible databases (NCBI SRA [http://www.ncbi.nlm.nih.gov/sra], MG-RAST [https://metagenomics.anl.gov], and VAMPs [https://vamps.mbl.edu/]) with accession numbers listed in Table 1.

6. Statistical analyses

Analysis of variance (ANOVA) tests were performed to assess OTUs that significantly differed between sample categories. T-tests were used to determine α-diversity significance. Analysis of similarities (ANOSIM) tests were performed to determine β-diversity significance.

Results

Differential community richness and diversity among CF and NIC microbiota

Alpha diversity measurements were used to assess the variation in community structure between the microbiota of CF patients and NICs. Specifically, Chao1 and Shannon indices were employed to analyze community richness and diversity, respectively (Fig 1). These results revealed greater richness of CF Baseline (Fig 1A) and lower diversity of Exacerbated CF (Fig 1B) in CF microbiota than the microbiota of NIC. Comparable levels of bacterial abundance were observed between CF and NIC microbiota, and these levels were remarkably lower than those of CF-Baseline (CF-Base) microbiota (p < 0.01, Fig 1A). The microbiota of CF patients with pulmonary exacerbation (CF-Exa) had the lowest community diversity. These results suggest that the microbial community of the CF-Exa microbiota has lower community diversity and richness than NIC microbiota.

thumbnail
Fig 1. Boxplots for alpha diversity measurements of the microbiota of CF subgroups.

Chao1 and Shannon indices were used to determine the richness and diversity of the bacterial community, respectively. Chao1 (a) and Shannon (b) measurements of the microbiota of Uninfected Control (NIC), CF baseline, and Exacerbated CF, Chao1 (c) and Shannon (d) measurements of the microbiota of Before, During and AfterTreatment of exacerbated CF. The upper edge of each box represents the 75% line, the lower edge represents the 25% line, and the inner line of each box represents the 50% line. ** refers to statistically highly significant as p < 0.01.

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

As the microbiota of the two CF stages and NICs exhibited divergence in the alpha diversity measurements, beta diversity indices were then calculated to depict this divergence. The unweighted UniFrac-PCoA analysis demonstrated that the microbiota of NICs can be easily distinguished from that of CF patients (Fig 2A and 2B). No marked differences were observed between the microbiota of CF-Base and CF-Exa (Fig 2F). This analysis suggests that there is a differential divergence between the community structure of the NIC microbiota and CF-Base/CF-Exa microbiota, whereas no remarkable divergence is observed between the CF-Base and CF-Exa microbiota.

thumbnail
Fig 2. Unweighted UniFrac PCoA plot of healthy and CF microbiota.

Distance comparisons were made for microbiota of: (a) Uninfected control (NIC) and CF baseline, (b) NIC and Exacerbated CF, (c) NIC and After Treatment of Exacerbated CF, (d) NIC and On Treatment of Exacerbated CF, (e) NIC and Before Treatment of Exacerbated CF, (f) CF baseline and Exacerbated CF, (g) Before, During and After Treatment of Exacerbated CF, (h) CF baseline and Before, During, After Treatment of Exacerbated CF. Each point represents a sample from the studies listed in Table 1. The color for each sample group can be found on the figure.

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

CF microbiota was featured with dominant opportunistic bacteria and dramatically decreased commensal bacteria

Venn diagrams were then generated to analyze the differences in the bacterial composition between the CF and NIC microbiota (Fig 3). The NIC microbiota possessed far fewer bacterial species (1470) than the CF-Base (3764) and the CF-Exa microbiota (2942) (Fig 3A). About 76.5% (1125/1470) of the taxonomies found in the NIC microbiota were also observed in CF microbiota. More specifically, 73.6% (1082/1470) of these taxonomies were present in CF-Base microbiota, 55.9% (822/1470) were present in CF-Exa microbiota, and 53.0% (779/1470) were found in both types of CF microbiota (Fig 3A). About 46.4% (682/1470) of the taxonomies in NIC microbiota were distinct from the treatment-related CF-Exa microbiota (Fig 3C). In contrast, only 29.4% (852/2942) and 28.7% (1082/3764) of the taxonomies identified in CF-Base and CF-Exa microbiota, respectively, were found in NIC microbiota. Furthermore, all of the taxonomies found in CF-Base microbiota were also present in treatment-related CF-Exa microbiota (Fig 3D). These results demonstrate that CF microbiota consist of much higher abundances of taxonomies than NIC microbiota, whereas CF-Exa microbiota harbor the largest number of distinct taxonomies that differ from those of CF-Base and NIC microbiota.

thumbnail
Fig 3. Venn diagrams for the microbiota of healthy individuals and patients in the various CF subgroups.

Comparisons were made for: (a) Uninfected control (NIC), CF baseline, and Exacerbated CF; (b) Before, During and After Treatment of Exacerbated CF; (c) NIC and Before, During, After Treatment of Exacerbated CF; (d) CF baseline and Before, During, After Treatment of Exacerbated CF. The number inside each circle represents the amounts of OTUs observed.

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

To probe the differences in taxonomies present in CF and NIC microbiota further, we compared the rank abundance of OTUs at the genus level (Fig 4). Only OTUs with relative abundances larger than 3% were shown. Pseudomonas was the major component of CF-Base and CF-Exa microbiota with a relative abundance of approximately 30% (Fig 4); this predominance of Pseudomonas was not observed in NIC microbiota (Fig 4). Similarly, Staphylococcus and Alcaligenaceae were only observed in CF microbiota with high relative abundance (Fig 4) but not in the NIC microbiota (Fig 4). Streptococcus, Veillonella, and Prevotella were the dominant genera in NIC microbiota (Fig 4). These results demonstrate the distinct architectures of CF and NIC microbiota.

thumbnail
Fig 4. Rank abundance of the dominant bacteria in the microbiota of healthy individuals and patients in the various CF subgroups.

Relative abundance of each genus in the three groups (Uninfected control (NIC), CF baseline, and Exacerbated CF) was shown. Only the bacterial genus with relative abundances higher than 3% were shown. Symbol representation of each genus was listed at the right side of the figure.

https://doi.org/10.1371/journal.pone.0164510.g004

Notable but transient effects of antibiotic treatments were observed in the microbiota of CF patients with pulmonary exacerbation

To evaluate the effects of antibiotic treatment on CF-Exa microbiota, community structures of treatment-related microbiota were compared (Fig 1C and 1D). Consistent with antibiotic treatments, the Exa-DuringTreat microbiota exhibited remarkably reduced community richness (Fig 1C; p < 0.01) and slightly increased community diversity (Fig 1D; p = 0.24) compared to the microbiota of Exa-BeforeTreat. These findings demonstrate that antibiotic treatment plays an important role on the community structure of the microbiota in patients of CF with pulmonary exacerbation. However, the community richness and diversity of Exa-AfterTreat microbiota returned to levels similar to those of Exa-BeforeTreat microbiota (Fig 1C and 1D), indicating that the effects of antibiotic treatment on the structure of microbiota are transient. These results demonstrate that the notable decrease in community richness in CF-Exa microbiota due to antibiotic treatment was transient, as the architecture of the community returned to Exa-Before after the withdrawal of antibiotics.

To depict the divergence of microbiota caused by antibiotic treatment further, beta diversity indices were analyzed (Fig 2G and 2H). The Unweighted UniFrac-PCoA showed that CF-Exa microbiota could not be discriminated from one another in the community structure (Fig 2G). Moreover, no visible differences were observed between CF-Base microbiota and treatment-related CF-Exa microbiota (Fig 2H). These results demonstrate that the notable effects on community richness caused by antibiotic treatment were transient and did not result in the visible divergence in community structures observed in treatment-related microbiota.

Antibiotic treatment causes limited effects on the dominant opportunistic bacteria but significant effects on the commensal bacteria

Next, OTUs of CF microbiota were compared to further evaluate the impact of antibiotic treatment on the bacterial composition of CF-Exa microbiota (Fig 3B). The Exa-AfterTreat microbiota covered all of the OTUs of the Exa-BeforeTreat and the Exa -DuringTreat microbiota, 59.9% (3394/5663) of which were distinct from the other two groups. Of the 3394 distinct taxonomies in Exa-AfterTreat microbiota, 1.7% (58/3394) and 66.7% (2263/3394) were shared with NIC (Fig 3C) and CF-Base (Fig 3D) microbiota. The rest (1671/3394, 49.2%) of the Exa-AfterTreat microbiota were unique. These results indicate that the withdrawal of antibiotics in CF patients with pulmonary exacerbation correlated with the proliferation of various taxonomies in CF-Exa microbiota.

To further assess the differences in bacterial composition among CF exacerbation microbiota, OTUs of the microbiota from treatment-related stages of CF with pulmonary exacerbation were compared (Fig 4). Only the OTUs with relative abundances more than 0.5% were shown. Antibiotic treatment dramatically reduced the number of OTUs with high relative abundance (Fig 4). The relative abundances of the Staphylococcus, Gemella, Actinomyces, Moraxellaceae, and Fusobacterium genera were notably reduced, whereas those of, the Prevotella and Streptococcus genera were increased, and those of the Pseudomonas and Alcaligenaceae genera were not affected. After the withdrawal of antibiotics, the richness of the OTUs with high relative abundances was either restored or greater than that of Exa-BeforeTreat microbiota (Fig 4). Specifically, the relative abundances of Staphylococcus, Gemella, Actinomyces, and Moraxellaceae were restored (Fig 4). These results demonstrate that antibiotic treatments had no notable effect on the relative abundance of the dominant opportunistic bacteria, but they significantly reduced the relative abundance of commensal bacteria.

Consistency of the results across studies despite variables

To verify that the results above were not affected by the cross-sectional data enrolled in this study, CF-Exa microbiota from individual studies were analyzed separately using the same close reference OTU picking strategy described above (Fig 5). Analyses of the microbiota data from SRP025173 demonstrated that the relative abundances of the opportunistic bacteria, Pseudomonas and Staphylococcus, were slightly affected by antibiotic treatments but were restored after the withdrawal of antibiotics, whereas the relative abundances of the commensal bacteria, Streptococcus and Prevotella, were significantly reduced and were not restored after the withdrawal of antibiotic treatment (Fig 5A). Analyses of the microbiota data from PC02 demonstrated that the relative abundances of the opportunistic bacteria Pseudomonas and Burkholderia were higher in microbiota after antibiotic treatments compared to those before treatment (Fig 5B). Although the presence of the opportunistic bacteria Staphylococcus was diminished following antibiotic treatments, new opportunistic bacteria including Burkholderia appeared with higher relative abundance. Moreover, the relative abundances of the commensal bacteria were significantly reduced following treatment (Fig 5C). These results confirm that antibiotic treatments significantly affect the levels of commensal bacteria but do not impact the dominant opportunistic bacteria following the withdrawal of antibiotics.

thumbnail
Fig 5. Dynamic variance of the relative abundances for the various bacteria in microbiota of exacerbated CF generated from individual studies.

The relative abundances of the top two opportunistic and commensal bacteria are shown. The blank, zebra crossing, and solid column represent the Before, During and After Treatment of Exacerbated CF, respectively. (a) Data from SRP025173, (b) PC02, and (c) SRP044880. Detailed information for these three studies is listed in Table 1. * refers to statistically significant as p < 0.05, ** refers to statistically highly significant as p < 0.01.

https://doi.org/10.1371/journal.pone.0164510.g005

To verify the consistency of results across studies despite inherent variables, beta diversity measurements of Clinically Stable, Exacerbated, and NIC samples were compared using query of study assessments of the 16s rRNA region (Fig 6). No significant differences were observed across these studies (Fig 6A) or in the 16s rRNA regions (Fig 6B) in the Clinically Stable (Fig 6A and 6B), Exacerbated (Fig 6C and 6D), and NIC (Fig 6E and 6F) samples. These results highlight the consistency of the results in this and other studies.

thumbnail
Fig 6. Beta-diversity consistency across studies and 16s rRNA regions.

Samples of (a-b) CF baseline, (c-d) Exacerbated CF, and (e-f) Uninfected control (NIC) were compared with queries of studies (a, c, e) and 16s rRNA regions (b, d, f).

https://doi.org/10.1371/journal.pone.0164510.g006

Discussion

This study analyzed a large cohort of lung microbiota samples from CF patients and healthy individuals from cross-sectional studies to compare the community structures and to evaluate the impact of antibiotic treatment on lung microbiota. A sophisticated methodology for handling cross-sectional microbiota studies was used to ensure the reliability of the acquired results [3839]. Overall our results demonstrate that antibiotic treatments for CF exacerbation have limited effects on opportunistic bacteria but dramatic effects on commensal bacteria. These findings highlight the critical need to reevaluate and optimize the current strategy of antibiotic treatments used to treat CF patients with pulmonary exacerbations.

Healthy microbiota in the airway mucosa is essential to properly shape the immune response. Changing a well-balanced “healthy” microbiota to an unhealthy restricted microbiota renders the airway mucosa increasingly susceptible to pathogens such as Pseudomonas and consequent lung injury [44, 45, 46]. Thus, comparing the differences between CF and NIC microbiota to evaluate the dynamic changes in CF lung microbiota is important for understanding the resilience of the polymicrobial ecosystem to intense antibiotic treatments. Unlike previous studies (Table 1), this study enrolled uninfected healthy lung microbiota to serve as a reference for comparison between CF and NIC microbiota. We also evaluated the effects of antibiotic treatments on CF lung microbiota. Consistent with the results of other chronic lung diseases (i.e., chronic obstructive pulmonary disease [16], asthma [12], interstitial pneumonia [47], etc.), the structure of the bacterial community in the NIC microbiota was completely different from that found in the CF microbiota. CF microbiota had greater community richness and lower diversity compared to NIC microbiota, indicating that the colonization of opportunistic bacteria had important effects on the community structure of the lung microbiota. These differences were further supported by beta diversity analyses and were consistent with previous reports [48], demonstrating that CF microbiota is distinct from NIC microbiota.

Many studies have long searched for a correlation between certain lung microbes and the severity of CF disease. To date, only the loss of microbial diversity has been associated with severe lung disease in CF [49]; this notion was further supported by the results of the present study. However, the driving force for this decreased microbial diversity remains unclear because there are many complexities involved in the progression of CF disease including lung function, patient age, and antibiotic treatment. Antibiotic treatment, rather than lung function or patient age, was reported to be the primary cause of the decrease in microbiota diversity noted in CF patients who declined from clinically stable to exacerbated states [50]. Nevertheless, indistinguishable bacterial community structures were observed between CF patients deemed as clinically stable and those categorized as exacerbated [18, 51]. However, heavy doses of antibiotics are commonly used to control exacerbated CF states rather than maintain clinically stable states [33]. Thus, the correlation between antibiotic treatments and CF disease progression is weak.

Although intense empiric antibiotic treatment is commonly used to control CF with pulmonary exacerbation, its effects on CF-Exa microbiota must be evaluated to further optimize this strategy [27, 33, 52]. Here, we evaluated healthy and CF clinically stable microbiota as references, which revealed that the most obvious effect of antibiotic treatments on CF-Exa microbiota was the significant reduction in community richness (Fig 1). These results are consistent with previous reports [33, 50]. However, there was no reduction in the opportunistic bacteria targeted by the antibiotics. The relative abundance of Pseudomonas was not significantly affected by antibiotic treatment (Fig 4), which is also consistent with previous reports [18, 51], Thus, since Pseudomonas is the major opportunistic bacteria found in CF patients with pulmonary exacerbation, more effective antibiotics targeting Pseudomonas species are required to control these bacteria. Interestingly, the relative abundance of Staphylococcus, another important opportunistic bacteria, was dramatically reduced by antibiotic treatment but these abundance levels returned to their original levels after the withdrawal of antibiotics. This noteworthy phenomenon has also been reported in a previous study [22] and demonstrates the limitations of current antibiotic therapies. These drawbacks highlight the need to reevaluate and optimize current antibiotic treatments for CF patients with pulmonary exacerbation. Because large overdoses of antibiotics have been commonly prescribed for treating CF to ensure a curative effect, it is reasonable to conclude that the varied outcomes of antibiotics should not be the key factor responsible for the differences observed in this study. Furthermore, the concordance of results between individual data and combined data enrolled in this study (Fig 5) also support this speculation. That our results are consistent with previous studies (Fig 5) also suggests that the methods and results of this study are highly reliable.

In summary, this study employed microbiota data from a collection studies to analyze the differences between the microbiota of uninfected healthy individuals and CF patients and to evaluate the impact of antibiotic treatments for CF. CF-Exa microbiota had comparable richness but lower diversity compared to NIC microbiota. Antibiotic treatments for CF patients with pulmonary exacerbation slightly reduced the presence of opportunistic bacteria, but markedly reduced the abundance of commensal bacteria and allowed the proliferation of untargeted opportunistic bacteria. Thus, this study highlights the deleterious influences of antibiotic treatment to microbiota of CF exacerbation, which emphasizes the need to optimize empiric antibiotic treatments.

Acknowledgments

We thank Dr. Mingkun Li from the Fondation Merieux for his critical discussions and helpful suggestions.

Author Contributions

  1. Conceptualization: JL JW XQ.
  2. Data curation: JL CH.
  3. Formal analysis: JL CH LR YX.
  4. Funding acquisition: JL.
  5. Investigation: JL CH LR.
  6. Methodology: JL CH LR.
  7. Project administration: JL LR JW.
  8. Resources: XQ JW.
  9. Software: JL CH YX.
  10. Supervision: XQ JW.
  11. Validation: LR JW.
  12. Visualization: JL CH LR.
  13. Writing – original draft: JL.
  14. Writing – review & editing: LR XQ JW.

References

  1. 1. Sloane PA, Rowe SM. Cystic fibrosis transmembrane conductance regulator protein repair as a therapeutic strategy in cystic fibrosis. Curr Opin Pulm Med. 2010; 16(6): 591–597. pmid:20829696
  2. 2. Milla CE, Moss RB. Recent advances in cystic fibrosis. Clin in Chest Med. 2012; 33(2): 307–328.
  3. 3. Stenbit AE, Flume PA. Pulmonary exacerbations in cystic fibrosis. Curr Opin Pulm Med. 2011; 17(6): 442–447. pmid:21881509
  4. 4. Tildy BE, Rogers DF. Therapeutic Options for Hydrating Airway Mucus in Cystic Fibrosis. Pharmacology.2015; 95(3–4): 117–132. pmid:25823699
  5. 5. Kreda SM, Davis CW, Rose MC. CFTR, mucins, and mucus obstruction in cystic fibrosis. Cold Spring Harb Perspect Med. 2012; 2(9): a009589. pmid:22951447
  6. 6. Pressler T. Targeting airway inflammation in cystic fibrosis in children: past, present, and future. Paediatr Drugs. 2011; 13(3): 141–147. pmid:21500868
  7. 7. Cantin AM, Hartl D, Konstan MW, Chmiel JF. Inflammation in cystic fibrosis lung disease: Pathogenesis and therapy. J Cyst Fibros. 2015; 14(4): 419–430. pmid:25814049
  8. 8. Rabin HR, Surette MG. The cystic fibrosis airway microbiome. Curr Opin Pulm Med. 2012; 18(6): 622–627. pmid:22965275
  9. 9. Beck JM. ABCs of the lung microbiome. Ann Am Thorac Soc. 2014; 11 Suppl 1(11):S3–6. pmid:24437402
  10. 10. Zemanick ET, Sagel SD, Harris JK. The airway microbiome in cystic fibrosis and implications for treatment. Curr Opin Pediatr. 2011; 23(3): 319–324. pmid:21494150
  11. 11. Cernadas M. It takes a microbiome: commensals, immune regulation, and allergy. Am J Respir Crit Care Med. 2011; 184(2): 149–150. pmid:21765026
  12. 12. Beck JM, Young VB, Huffnagle GB. The microbiome of the lung. Transl Res. 2012; 160(4): 258–266. pmid:22683412
  13. 13. Dickson RP, Erb-Downward JR, Huffnagle GB. The role of the bacterial microbiome in lung disease. Expert Rev Respir Med. 2014; 7(3): 245–257.
  14. 14. Field TR, Sibley CD, Parkins MD, Rabin HR, Surette MG. The genus Prevotella in cystic fibrosis airways. Anaerobe. 2010; 16(6): 337–344.
  15. 15. Morris A, Beck JM, Schloss PD, Campbell TB, Crothers K, Curtis JL, et al. Comparison of the respiratory microbiome in healthy nonsmokers and smokers. Am J Respir Crit Care Med. 2013; 187(10): 1067–1075. pmid:23491408
  16. 16. Erb-Downward JR, Thompson DL, Han MK, Freeman CM, McCloskey L, Schmidt LA, et al. Analysis of the lung microbiome in the "healthy" smoker and in COPD. PLoS One. 2011; 6(2): e16384. pmid:21364979
  17. 17. Nguyen LD, Viscogliosi E, Delhaes L. The lung mycobiome: an emerging field of the human respiratory microbiome. Front Microbiol. 2015; 6: 89. pmid:25762987
  18. 18. Fodor AA, Klem ER, Gilpin DF, Elborn JS, Boucher RC, Tunney MM, et al. The adult cystic fibrosis airway microbiota is stable over time and infection type, and highly resilient to antibiotic treatment of exacerbations. PLoS One. 2012; 7(9): e45001. pmid:23049765
  19. 19. Hampton TH, Green DM, Cutting GR, Morrison HG, Sogin ML, Gifford AH, et al. The microbiome in pediatric cystic fibrosis patients: the role of shared environment suggests a window of intervention. Microbiome. 2014; 2(1): 1–5.
  20. 20. Nelson A, De Soyza A, Perry JD, Sutcliffe IC, Cummings SP. Polymicrobial challenges to Koch's postulates: ecological lessons from the bacterial vaginosis and cystic fibrosis microbiomes. Innate Immun. 2012; 18(5): 774–783. pmid:22377802
  21. 21. Lynch SV, Bruce KD. The cystic fibrosis airway microbiome. Curr Opin Pulmon Med. 2012; 18 (6): 622–627.
  22. 22. Willger SD, Grim SL, Dolben EL, Shipunova A, Hampton TH, Morrison HG, et al. Characterization and quantification of the fungal microbiome in serial samples from individuals with cystic fibrosis. Microbiome. 2014; 2(1): 1–15.
  23. 23. Zemanick ET, Wagner BD, Robertson CE, Stevens MJ, Szefler SJ, Accurso FJ, et al. Assessment of airway microbiota and inflammation in cystic fibrosis using multiple sampling methods. Ann Am Thorac Soc. 2014; 12(2): 221–229.
  24. 24. Saiman L. Infection prevention and control in cystic fibrosis. Curr Opin Infect Dis. 2011; 24(4): 390–395. pmid:21543979
  25. 25. Ashish A, Shaw M, Winstanley C, Ledson MJ, Walshaw MJ. Increasing resistance of the Liverpool Epidemic Strain (LES) of Pseudomonas aeruginosa (Psa) to antibiotics in cystic fibrosis (CF)—a cause for concern? J Cyst Fibros. 2012; 11(3): 173–179. pmid:22146482
  26. 26. Jain K, Smyth AR. Current dilemmas in antimicrobial therapy in cystic fibrosis. Expert Rev Respir Med. 2012; 6(4): 407–422. pmid:22971066
  27. 27. Hurley MN, Forrester DL, Smyth AR. Antibiotic adjuvant therapy for pulmonary infection in cystic fibrosis. Paediatr Respir Rev. 2013; 13(3): 160–161.
  28. 28. Sherrard LJ, Tunney MM, Elborn JS. Antimicrobial resistance in the respiratory microbiota of people with cystic fibrosis. Lancet. 2014; 384(9944): 703–713. pmid:25152272
  29. 29. Cuthbertson L, Rogers GB, Walker AW, Oliver A, Hafiz T, Hoffman LR, et al. Time between collection and storage significantly influences bacterial sequence composition in sputum samples from cystic fibrosis respiratory infections. J Clin Microbiol. 2014; 52(8): 3011–3016. pmid:24920767
  30. 30. Hassett DJ, Korfhagen TR, Irvin RT, Schurr MJ, Sauer K, Lau GW, et al. Pseudomonas aeruginosa biofilm infections in cystic fibrosis: insights into pathogenic processes and treatment strategies. Expert Opin Ther Targets. 2010; 14(14): 117–130.
  31. 31. Narasimhan M, Cohen R. New and investigational treatments in cystic fibrosis. Ther Adv Respir Dis. 2011; 5(4): 275–282. pmid:21372122
  32. 32. Hurley MN, Forrester DL, Smyth AR. Antibiotic adjuvant therapy for pulmonary infection in cystic fibrosis. Paediatr Respir Rev. 2013; 13(3): 160–161.
  33. 33. Smith DJ, Badrick AC, Zakrzewski M, Krause L, Bell SC, Anderson GJ, et al. Pyrosequencing reveals transient cystic fibrosis lung microbiome changes with intravenous antibiotics. Eur Respir J. 2014; 44: 922–930. pmid:25034564
  34. 34. Kramer R, Sauer-Heilborn A, Welte T, Jauregui R, Brettar I, Guzman CA, et al. High individuality of respiratory bacterial communities in a large cohort of adult cystic fibrosis patients under continuous antibiotic treatment. PLoS One. 2015; 10(2): e0117436. pmid:25671713
  35. 35. Kosiewicz MM, Zirnheld AL, Alard P. Gut microbiota, immunity, and disease: A complex relationship. Front Microbiol. 2011; 2(1): 180.
  36. 36. Sekirov I, Shannon L, Russell L, Antunes CM, Finlay BB. Gut microbiota in health and disease. Physiol Rev. 2010; 90 (3): 859–904. pmid:20664075
  37. 37. Patel RK, Jain M. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS One. 2012; 7(2): e30619. pmid:22312429
  38. 38. Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vazquez-Baeza Y, et al. Meta-analyses of studies of the human microbiota. Genome Res. 2013; 23(10): 1704–1714. pmid:23861384
  39. 39. Walters WA, Xu Z, Knight R. Meta-analyses of human gut microbes associated with obesity and IBD. FEBS Lett. 2014; 588(22): 4223–4233. pmid:25307765
  40. 40. McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012; 6(3): 610–618. pmid:22134646
  41. 41. Kuczynski J, Stombaugh J, Walters WA, Gonzalez A, Caporaso JG, Knight R. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr Protoc Microbiol. 2011; Chapter 10: Unit 10.7.
  42. 42. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010; 26(19): 2460–2461. pmid:20709691
  43. 43. McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013; 8(4): e61217. pmid:23630581
  44. 44. Lemon KP, Armitage GC, Relman DA, Fischbach MA. Microbiota-targeted therapies: an ecological perspective. Sci Transl Med. 2012; 4(137): 4113–4122.
  45. 45. Segal LN, Rom WN, Weiden MD. Lung microbiome for clinicians. New discoveries about bugs in healthy and diseased lungs. Ann Am Thorac Soc. 2014; 11(1): 108–116. pmid:24460444
  46. 46. Huang YJ, Charlson ES, Collman RG, Colombini-Hatch S, Martinez FD, Senior RM. The role of the lung microbiome in health and disease. A National Heart, Lung, and Blood Institute workshop report. Am J Respir Crit Care Med. 2013; 187(12): 1382–1387. pmid:23614695
  47. 47. Molyneaux PL, Cox MJ, Willis-Owen SA, Mallia P, Russell KE, Russell AM, et al. The role of bacteria in the pathogenesis and progression of idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2014; 190: 906–913. pmid:25184687
  48. 48. Blainey PC, Milla CE, Cornfield DN, Quake SR. Quantitative analysis of the human airway microbial ecology reveals a pervasive signature for cystic fibrosis. Sci Transl Med. 2012; 4(153): 2569–2575.
  49. 49. Flight WG, Smith A, Paisey C, Marchesi JR, Bull MJ, Norville PJ, et al. Rapid detection of emerging pathogens and the loss of microbial diversity associated with severe lung disease in cystic fibrosis. J Clin Microbiol. 2015; 53(7):2022–2029. pmid:25878338
  50. 50. Zhao J, Schloss PD, Kalikin LM, Carmody LA, Foster BK, Petrosino JF, et al. Decade-long bacterial community dynamics in cystic fibrosis airways. Proc Natl Acad Sci U S A. 2012; 109(15): 5809–5814. pmid:22451929
  51. 51. Price KE, Hampton TH, Gifford AH, Dolben EL, Hogan DA, Morrison HG, et al. Unique microbial communities persist in individual cystic fibrosis patients throughout a clinical exacerbation. Microbiome. 2013; 1: 27. pmid:24451123
  52. 52. Conrad D, Haynes M, Salamon P, Rainey PB, Youle M, Rohwer F. Cystic fibrosis therapy: a community ecology perspective. Am J Respir Cell Mol Biol. 2012; 48(2): 150–156. pmid:23103995
  53. 53. Zemanick ET, Harris JK, Wagner BD, Robertson CE, Sagel SD, Stevens MJ, et al. Inflammation and Airway Microbiota during Cystic Fibrosis Pulmonary Exacerbations. PLoS one. 2013; 8(4): e62917. pmid:23646159
  54. 54. Rutebemberwa A, Stevens MJ, Perez MJ, Smith LP, Sanders L, Cosgrove G, et al. Novosphingobium and Its Potential Role in Chronic Obstructive Pulmonary Disease: Insights from Microbiome Studies. PLoS one. 2014; 9(10): e111150 pmid:25340840
  55. 55. Cuthbertson L, Rogers GB, Walker AW, Oliver A, Hoffman LR, Carroll MP, et al. Implications of multiple freeze-thawing on respiratory samples for culture-independent analysis. J Cyst Fibros. 2015; 14(4): 464–467. pmid:25459563
  56. 56. Filkins LM, Hampton TH, Gifford AH, Gross AH. Gross MJ, Hogan DAet al. Prevalence of Streptococci and Increased Polymicrobial Diversity Associated with Cystic Fibrosis Patient Stability. J Bacterology. 2012; 194(17):4709–4717.
  57. 57. Boutin S, Graeber SY, Weitnauer M, Panitz J, Stahl M, Clausznitzer D, et al. Comparison of Microbiomes from Different Niches of Upper and Lower Airways in Children and Adolescents with Cysitc Firbosis. PLoS one. 2015; 10(1): e0116029. pmid:25629612
  58. 58. Maughan H, Wang PW, CAballero JD, Fung P, Gong Y, Donaldson SL, et al. Analysis of the Cystic Fibrosis Lung Microbiota via Serial Illumina Sequencing of Bacterial 16s rRNA Hypervaribale Regions. PLoS one. 2012; 7(10): e45791. pmid:23056217
  59. 59. Bassis CM, Erb-Downward JRE, Dickson RP. Freeman CM, Schmidt TM, Young VB, et al. Analysis of the Upper Respiratory Tract Microbiotas as the Source of the Lung and Gastric Microbiotas in Healthy Individuals. mBio. 2015; 6(2): e00037–15. pmid:25736890
  60. 60. Pragman AA, Kim HB, Reilly CS, Wendt C, Isaacson RE. The Lung Microbiome in Moderate and Severe Chronic Obstructive Pulmonary Disease. PLoS one. 2012; 7(10): e47305. pmid:23071781