Abstract
Lung infections play a critical role in cystic fibrosis (CF) pathogenesis. CF respiratory tract is now considered to be a polymicrobial niche and advances in high-throughput sequencing allowed to analyze its microbiota and mycobiota. However, no NGS studies until now have characterized both communities during CF pulmonary exacerbation (CFPE). Thirty-three sputa isolated from patients with and without CFPE were used for metagenomic high-throughput sequencing targeting 16S and ITS2 regions of bacterial and fungal rRNA. We built inter-kingdom network and adapted Phy-Lasso method to highlight correlations in compositional data. The decline in respiratory function was associated with a decrease in bacterial diversity. The inter-kingdom network revealed three main clusters organized around Aspergillus, Candida, and Scedosporium genera. Using Phy-Lasso method, we identified Aspergillus and Malassezia as relevantly associated with CFPE, and Scedosporium plus Pseudomonas with a decline in lung function. We corroborated in vitro the cross-domain interactions between Aspergillus and Streptococcus predicted by the correlation network. For the first time, we included documented mycobiome data into a version of the ecological Climax/Attack model that opens new lines of thoughts about the physiopathology of CF lung disease and future perspectives to improve its therapeutic management.
Similar content being viewed by others
Introduction
Lung infections play a critical role in cystic fibrosis (CF) pathogenesis where they can lead to significant acute decrease of lung function, known as CF pulmonary exacerbation (CFPE). Developments of next-generation sequencing (NGS) approaches allowed us to understand microbiome composition that can contribute to lung physiopathology in both healthy individuals and patients with lung disease. More recently, NGS together with advances into statistical network inference tools allowed to analyze the microbial airway communities, appreciate the inter-kingdom equilibrium of respiratory floras, and therefore develop understanding of microbial communities as a whole1,2,3,4,5,6,7.
Acute CFPEs represent major clinical events that significantly decline the lung function, contribute to disease progression and require adapted, prompt anti-infectious treatment8,9,10,11. Omics approaches confirmed associations between bacterial community and exacerbation12,13,14,15,16,17,18,19,20. Apart from bacteria that are well known agents causing dramatic recurrent CFPEs, respiratory viruses have been recently found to be associated with CFPE, independently of bacteria or as trigger for bacterial infection21,22,23. Fungal species such as Candida albicans and Aspergillus fumigatus have been associated with decreased lung function as well as increased frequency of exacerbations24,25. However, only limited number of studies evaluated mycobiome during CFPE12,26, and never in association with bacterial community to address inter-kingdom crosstalk.
Although much has been learned about CFPEs, our current knowledge about its pathophysiology and usefulness of diagnostic criteria remains limited. Since CFPEs are difficult to predict, new biomarkers that may derive from current NGS approaches are needed to complete existing assessments9,13,14. Given the prevalence of fungal colonization, the predominance of allergic bronchopulmonary aspergillosis (ABPA), the growing role of fungal infection in CF and the stabilizing function that fungi appear to play in lungs1,27,28,29,30, we hypothesized that fungi are involved in CF pulmonary disease evolution and in CFPEs. In the present study we considered both fungal and bacterial communities as a unique pathogenic entity, and interpreted our results according to a co-occurrence/co-exclusion pattern approach that highlights relationships between microorganisms sharing the same ecological niche, i.e. the CF respiratory tract. According to our inter-kingdom network results, we selected species belonging to Aspergillus and Streptococcus genera, and documented experimentally the predicted positive growth of A. fumigatus when it is co-cultured with species from Streptococcus mitis group. By adapting Phy-Lasso method we pointed out the fungi relevantly associated with CFPE and/or a decline in lung function. Collectively, these data provided an opportunity to reflect on a revised version of the CF Climax/Attack model (CAM) derived from ecology research field31. In CAM model, CF lung infection is no longer interpreted as an invasion of a single pathogen such as Pseudomonas aeruginosa but as dynamic changes within the whole microbial community composed of two types of microbial populations (an Attack transient but virulent population associated with CFPE and a Climax chronic population driving the long-term prognostic of CF lung disease)4,7,13,14,31,32. While it would be still early to draw definitive conclusions, our results pave the way for future studies on respiratory mycobiome that should help us to begin deciphering the physiopathology of CF lung disease and to improve its therapeutic management. They also suggest the complexity of all respiratory microbiome componentsâ connections, which might yield important insights within the CF disease.
Results
Patients and sample characteristics
Among the 37 patients initially screened, samples of patients 2, 7, 28 and 35 were excluded because their corresponding rarefaction curves and/or numbers of sequences were inadequate to allow biodiversity comparison. Samples from 33 patients, corresponding to 246,228 and 304,914 bacterial and fungal reads respectively, were used to evaluate microbial composition, diversity, and perform network analysis. Patientsâ characteristics, spirometry data, presence of conventionally cultured pathogens in sputum, and therapeutic management at sampling time are recorded in Table 1. Seventeen patients were classified as having CFPE while 16 were not. According to body mass index (BMI), Shwachman-Kulczycki score (S-K score), and spirometry data (in particular Forced Expiratory Volume in 1âs (FEV1) expressed as percent of predicted value), patients were in a medium to advanced stage of CF disease, with several co-morbidities (exocrine pancreatic dysfunction in 90% of cases, and diabetes mellitus or impaired glucose tolerance in 33% of cases). Colonization by Pseudomonas aeruginosa was highly prevalent (96.7%, Table 1); all excepted one patient (a 23-year-old woman, patient 14) had positive mycological cultures. Candida albicans, other Candida sp., and A. fumigatus were observed in 66.6%, 18.2%, and 66.6% of patients respectively, as determined by conventional culturing. Most of the patients received rhDNAse and/or inhaled steroids (100%), low dose of azithromycin (69.7%) and antibiotics (90.1% of cases administered orally). Antifungal therapy was prescribed for one patient in agreement with ABPA diagnosis33. While 11 (33%) patients were treated with systemic steroids, no other systemic immunosuppressive therapy was used. No significant difference for those variables was observed between patients with and without CFPE, except for the median age, which was significantly lower in patients with CFPE.
FEV1 was associated with changes in alpha diversity, and CFPE with shifts in minority populations
Rarefied data was used to summarize microbial composition and its changes among groups (Fig. 1). These microbial communities inferred from NGS were consistent with results obtained by culturing (Table 1); however greater diversity was evident via NGS approach. We observed a significant positive linear relationship between FEV1 values and Chao1 indexes of bacterial communities (Fig. 1a). This linear relationship was also observed when using Shannon diversity indexes (pâ=â0.04), but not using Simpson diversity indexes (pâ=â0.10). Different diversity metrics capture, indeed, different information. Microbial populations were divided into majority and minority ones12,34. Microbiota and mycobiota of patients with CFPE exhibited a significant decrease in Haemophilus, Neisseria, and Cladosporium genera among majority populations, as well as changes in relative abundances of bacteria significantly associated with a loss in alpha diversity of fungal and bacterial minority populations (Table 2). To evaluate better inter-kingdom equilibrium between the microbial communities, we applied a ratio of fungal to bacterial diversity previously reported in IBD35. Fungi-to-bacteria diversity ratios showed a dominance of bacteria on fungus with a median at 0.64 (IQRâ=â[0.43; 4.25]), without any significant disequilibrium between bacterial and fungal loads among patients with and without CFPE. Median of the Basidiomycota/Ascomycota relative abundance ratio was 0.03 (IQRâ=â[0.01; 0.19]).
Inter-kingdom correlation and network analysis
We analyzed co-presence or absence patterns between bacteria and fungi, by building a correlation matrix at genus levels of the whole set of NGS data and representing it as a graphical network (Figs. 2 and 3). We inferred an inter-kingdom network by plotting bacterial genera significantly correlated with at least one fungal genus and vice versa. The network was composed of 26 nodes and 23 edges and it revealed three main clusters organized around Aspergillus, Candida and Scedosporium genera plus three disconnected pairs of microorganisms (Fig. 2). Among these pairs, two contained environmental fungi: Ramularia (a phytopathogen) and Penicillium connected with the class of Gammaproteobacteria and the phylum of Proteobacteria, respectively. These fungi likely represent environmental transitient fungal elements, as previously proposed12,36,37. Among the six negatively correlated OTUs identified in our model, Fusarium sp. appeared to display a negative correlation with Gemella sp., the latter identified as keystone genus predicting CFPE4,13. The other positively correlated OTUs predict interactions between Aspergillus, Candida, or Scedosporium and bacteria belonging to Capnocytophaga, Parvimonas, Streptococcus, or Veillonella. Apart from Scedosporium species that were described as belonging to Climax population31, our network results preferentially associate fungi to the bacterial population identified during CFPE and referred as Attack population in the ecological Climax/Attack Model (CAM) recently proposed4,13,31.
Assessment of network interactions using in vitro co-cultures
To document inter-kingdom interactions that we identified in correlation matrix analysis, we performed co-culture experiments between A. fumigatus and Streptococcus mitis or Streptococcus oralis. We selected these microorganisms because of their well-established medical importance in CF (at the species level), high quantity in our cohort associated with a robust positive co-variation (at the genus level), and known conditions for in vitro cultures. As predicted, A. fumigatus growth appeared to be enhanced by both S. mitis and S. oralis (Fig. 4). As our experimental method was focused on A. fumigatus growth that has a duplication time clearly different from Streptococcus sp. (24â48âh for Aspergillus versus about 1âh for Streptococcus), it could not inform on bacterial growth.
Microbiome and mycobiome markers of CF lung dysfunction
To fully predict the bacterial and fungal communities associated with altered lung function, we applied a feature selection method (bootstrap-enhanced Phy-Lasso) to identify genera relevantly associated with CFPE (Fig. 5a) and/or FEV1 decline (Fig. 5b). In patientsâ microbiomes Haemophilus, Pseudomonas, Staphylococcus, Streptococcus, Candida, Fusarium, Penicillium, and Scedosporium were selected 150 times out of 200 bootstrap replicates and were negatively associated with the presence of CFPE. Malassezia and Aspergillus were positively associated with CFPE (Fig. 5a). Pseudomonas, Penicillium, and Scedosporium were negatively associated with FEV1, an increase in their relative abundance being associated with a decrease in FEV1 measure (Fig. 5b). Altogether, these results highlighted the complexity of the microbial community interacting within the CF respiratory tract, and suggested the suitability of developing ecological models such as CAM. Collectively, our findings (inter-kingdom network analysis, in vitro co-culture results, and feature selection based on bootstrap-enhanced Phy-Lasso) pave the way for deciphering the role of fungi in CF lung disease at the ecological level by proposing a new version of the recently-described CAM model31 (Fig. 6).
Discussion
CFPEs represent key intermittent events in CF disease that are poorly defined but clearly associated with the decline in lung function and accelerated disease progression. They require specifically adapted anti-microbial treatments8,9,10,11 but those treatments fail to recover lung function of 25% of patients38. Therefore omics approaches represent promising tools to identify and characterize relevant biomarkers able to predict CFPE and to improve therapeutic monitoring13. Considering that CF respiratory tract is polymicrobial, i.e. composed of bacteria, viruses and fungi, and since all these microorganisms are able to induce CFPEs, the large data set generated by omics methods should take into account all resident lung microorganisms. Here, we performed a combined analysis of CF airway microbiome and mycobiome, from the microbial diversity to the inter-kingdom network investigation in relation with patient clinical status.
We reported a decrease in alpha diversity of bacteria associated with FEV1 decline and, in agreement with published data12,14,15,16,17,19,20,26,36,37,39,40,41,42,43,44. We identified Pseudomonas, Streptococcus, Haemophilus, Candida, Fusarium, Penicillium, and Scedosporium as dominant genera (Fig. 1) and numerous shifts in both fungal and bacterial minority populations during CFPE (Table 2). Similarly, Cuthbertson et al.42 evidenced a greater turnover rate within the rare species compared to the microbial core. Moreover, the simultaneous analysis of mycobiota and microbiota allowed us to appreciate the inter-kingdom interactions. The observed ratios of fungus-to-bacteria diversity were 10 times higher compared to the published digestive ratios35, supporting the notion that lung mycobiome is derived from inhaled spores and transient species12,36,37,45. Our estimation of the Basidiomycota/Ascomycota relative abundance ratios was lower than digestive ones35 revealing the excess of Ascomycota, the phylum to which fungal pathogens belong. Together, these results suggest a colonization of airways during CF following an exposure to fungal spores.
Our study design has some notable limitations: (i) First, even though the clinical features of our CF population (Table 1) are congruent with published data21,26,27,37,39,44, the limited cohort we analyzed retrospectively was based on samples stored at â20â°C, which did not allowed us to realize RNA extraction and meta-transcriptomic analysis. (ii) Secondly, we used sputum samples, which could be contaminated with a buccal flora. However, this sample type is easy to assess, not invasive compared to bronchial alveolar lavage and sufficiently robust in predicting microbiome as previously demonstrated7,39,46. (iii) The small sample size (as for the majority of studies about the association between respiratory microbiota or mycobiota and CF12,26,44,47) decreased the statistical power to recognize small differences if present. (iv) Finally, the complex structure and the high-dimensional data generated by NGS analysis required specific statistical methods and computational tools that have to be carefully evaluated. We acknowledge that although correlation does not indicate causation; it can provide a reasonable starting point to assess inter-kingdom network, especially when supported by clinical knowledge about pathogenic species (Figs. 2 and 5) and experimental in-vitro corroboration (Fig. 4). For the first time, the mycobiome aspect of the CAM model recently adapted from the ecology field4,31 is conjectured from CF fungal data (Fig. 6).
Niche-microbe interaction models have been developed to better interpret the dynamics, composition and role of lung microbial communities4,7,31,48,49,50. These models were focused on bacterial communities and suggested that relative abundances of bacteria were not driven by a neutral model of dispersion in CF50. In CF respiratory tract, microbiome assembly is resulting from dynamic selective pressures that are changing with each stage of lung disease. In different layers characterized by local nutrient availability, pH, oxygen pressure, host immune response and antibiotic treatment, a selection of microbiome members occurred and consequently a shift in microbial composition. So far, only CAM model has attempted to formally describe the lung mycobiome31.
Even if the role of anaerobes is still matter of debate in terms of community composition, metabolome, and resistome51,52, two types of bacterial populations have been highlighted in CAM model. These two populations have distinct metabolic potential defined by carbon source: an Attack transient but virulent population associated with CFPE and a Climax chronic population driving the long-term prognostic4,7,13,14,32. While anaerobes have been first affiliated to Climax population31, the Attack community appeared to be composed of anaerobes (Prevotella, Parvimonas, Veillonella, Porphyromonas, and Fusobacterium), facultative anaerobes (Streptococcus, Granulicatella, and Rothia), and the obligatory fermentative (Gemella) that have been associated with fermentation pathway producing organics acids and decrease in pH4,7,13,14,32. The Climax community is composed of Pseudomonas, Staphylococcus, Stenotrophomonas, and Achromobacter that were strongly associated with ornithine biosynthesis, transport of putrescine, and that were producing ammonia and an increase in pH4,13. Taking into account this model4,7,13,14,32 and our results, data allowed us to make small steps toward achieving a revised version of CAM model in the CF context (Fig. 6), in which the mycobiome role is completed as follow.
The decline in lung function in CF is driven by clinical changes especially by CFPEs, which require active treatment to recover from CFPE and go back to a steady state that may include chronic colonization4,19,20. The triggering elements of state shifts are not well known, but it has been proposed that microbiome changes acquired through environmental exposure or aspiration participate to CFPE occurrence leading to the Attack population ((a) in Fig. 6)31,53. Unexpectedly, both Staphylococcus and Pseudomonas appeared to be negatively correlated with CFPE (Fig. 5). This result may be explained by a CF population composed of adult patients with a rather mild to moderate lung diseases than a severe one (Table 1). Besides being involved into CFPE, a recent interest research has emerged in determining whether Staphylococcus aureus can have both damaging and protecting effects for the CF host according to the differences phases in the evolution of the CF airways infection54. Our results (Fig. 5b) agree with earlier studies in which colonization with P. aeruginosa and methicillin-resistant S. aureus was strongly associated with a loss of lung function (a loss of FEV1) in CF55 while microbiota communities appeared to not differ based on pulmonary exacerbation status of CF patients39,56. After CFPE occurred, the perturbed microbial ecosystem will return to its original stable state (Resilience (b) in Fig. 6) or if the perturbation is not controlled, will evolve to a new ecosystem considered as stable but composed of a novel Climax population that differs from the previous community (Adaptation (bâ) in Fig. 6). According to the CFPE severity, different filters (changes in nutrient sources, level of microorganism growth and virulence, strong host innate immune and inflammatory response, antimicrobial treatment pressure (c) in Fig. 6) will select this new population well adapted to the airway remodeling, in a circular relationship. As microbial community composition and disease progression proceed hand-in-hand, assembly of CF microbiome and mycobiome will result from dynamic events, especially from competitive exclusion or co-occurrence events observed through OTU network and related functional analysis4,13. In agreement with the bacterial community partitioned by carbon source4, our results suggest fungal species preferentially as Attack members (Fig. 5), which is consistent with clinical and mycological data28,29,30. Among them, Aspergillus and Malassezia were significantly associated with CFPE (Figs. 5a and 6). This result is in agreement with the ability of Aspergillus to be responsible for ABPA and related exacerbations in CF28,29. It goes in the same direction of previous results demonstrating an increased risk of CFPE requiring hospitalization when patients are colonized with A. fumigatus25. As Aspergillus was previously assigned to Climax population31, our findings have to be confirmed in a larger CF population. Malassezia are lipophilic yeasts, and require specific culture conditions that are not routinely used when a pulmonary exacerbation is monitored. Therefore, Malassezia have not been associated with CFPE yet. However, they have been identified as biomarker of chronic inflammatory diseases such as asthma57, or inflammatory bowel disease35. It may reflects a clear advantage (cross - feeding) that these yeasts get in sharing ecological niche with anaerobes and Streptococcus, since they produce multiple secreted lipases that may help in recycling lipids produced by fermenting bacteria13,58. In addition, Malassezia -which are part of the normal skin flora- are able to develop at acid pH (skin surface pH is below 5).
On the other side, Scedosporium were significantly associated with a decline in FEV1 (Fig. 5b), a result congruent with their ability to induce chronic detrimental colonization of CF lung during Climax phases12,14,15,17,19,20,42,59,60,61,62,63. As Pseudomonas do, Scedosporium is able to produce secondary metabolites with antimicrobial activity64. It secretes a polyketide boydone A with an anti-staphylococcus activity against methicillin-resistant isolates, which fits well with the course of microbial airway colonization clinically observed in CF. Moreover, the recent whole genome sequencing of Scedosporium species (https://genome.jgi.doe.gov/programs/fungi/index.jsf) revealed a number of genes encoding putative ankyrin motif-containing proteins, methyltransferases and oxidoreductases larger than the one of A. fumigatus genome65. Whether this larger magnitude of activities allows the mold to use a wider range of nutritive substrates or to be able to synthesize specialized metabolites including ammonia fermentation remains to be elucidated64,65,66, but it may participate in its successful establishment within the respiratory tract of CF patients. Over time, the interchange of these states will lead to the emergence of new Climax communities that are more and more adapted to the airway remodeling, leading to a less diverse community, composed of microorganisms highly resistant to antibiotic or antifungal agents such as Pseudomonas or Scedosporium.
In summary, microbial interactions are clearly complex in natural communities such as lungs, due to the impact of a multitude of microbes (bacteria, fungi, but also phages and viruses) able to interact67,68 and to host immune response that appears dual68,69. We show that NGS approach combined with inference network and ecological model analysis are useful in helping to decipher physiopathology of CF lung disease, and can be considered as a promising tool for improving our therapeutic protocols. As we clearly improved our knowledge on sequenced genomes of the most frequent molds in CF (https://genome.jgi.doe.gov/programs/fungi/1000fungalgenomes.jsf), such inter-kingdom analysis should be easily reproduced in a larger cohort, completed to include all significant members of the CF lung microbial community, and extended to analyze the meta-transcriptome and metabolome, in order to confirm or to correct the CAM model efficiency and to define accurate biomarkers of CFPEs. Additional in vitro (cell cultures) and/or in vivo (animal) models are required to support our future findings and to fully understand the complex microbial interactions that drive lung microbial community evolution in CF.
Materials and Methods
Patients and sputum samples
Among CF patients enrolled into a 3-year multicenter national prospective study (PHRC number 06/1902, acronym: MucoFong), 37 sputa consecutively collected at hospital centers of Lille (nâ=â19), Dunkerque (nâ=â3), and Grenoble (nâ=â15) were sequenced using NGS. They were alternatively selected from the first 64 sample collection21 according to CFPE presence or absence, if spirometry, therapeutic, radiological and biological data were collected, and if residual sputum sample was sufficient to perform DNA extraction. The clinical status (pulmonary exacerbation or stable period) of each patient was defined by the physician as follow. Recent changes in clinical parameters and/or modifications of the pulmonary function provide criteria of exacerbation according to European Respiratory Society statement10.
Sputa were collected and analyzed according to a standardized protocol previously described21,70. DNA extraction was performed as previously described21.
Ethics statement
MucoFong study was approved by the institutional ethics committees of Lille Hospital (Reference Number: CPP-06/84), and a written informed consent was provided by all participants. All methods were performed in accordance with the relevant guidelines.
Targeted metagenomic library preparation, sequencing and taxonomic assignation
V3-V4 and ITS2 regions of bacterial and fungal rDNA were amplified and sequenced as previously described37,44. All reactions were prepared in a sterile PCR hood, and negative control reactions were performed and examined by electrophoresis on 1.5% agarose gel. We obtained 429,565 and 443,270 pyrosequences with prokaryotic and eukaryotic loci, respectively. Sequencing data processing was performed using the Metabiote© pipeline where potential chimeras, read with poor quality (â¤20 quality scores), short read (â¤150 nucleotides for 16S rRNA gene, â¤100 nucleotides for ITS2 rRNA gene), homopolymers, singletons, and doubletons were removed before the creation of Operational Taxonomic Unit (OTU). OTUs were created applying complete-linkage clustering with 97% similarity criteria (UClust v1.2.22q); OTUsâ annotation was performed using Ribosomal Database Project classifier against Greengenes (v13.8 containing 203,452 16âS rRNA sequences) and UNITE (v12.11 containing 97,868 ITS1-ITS2 sequences) databases (Supplementary Table 1). Rarefaction curves were calculated to determine whether deep sequencing was sufficient to accurately characterize the bacterial and fungal community diversity, as previously described44.
Microbial community composition analysis
Bacterial and fungal data were analyzed at phylum to species levels. Clinical status was assessed using CFPE (patient acutely exacerbated or stable), S-K score, BMI, and the respiratory capacity measured by FEV1 value. The later was also expressed as categorical by using cut-off values: âmildâ (FEV1â>â70%), âmoderateâ (FEV1 between 40â70%) andâsevereâ ventilatory deficit (FEV1â<â40%)59,71. Microbial composition was analyzed according to clinical status measured as both continuous and categorical variables.
Alpha diversity was measured though Simpson, Shannon, and Chao1 indexes. Richness and evenness differences between clinical status groups (S-K score, BMI, and FEV1) were tested using Wilcoxon signed-rank test. Since these alpha diversity indexes do not reflect accurately community composition (i.e. two communities could have the same diversity index but different community structures), principal coordinate analysis (PCoA) of beta diversity analysis was performed using Bray-Curtis distance. Structural differences between groups were tested using analysis of dissimilarities (ANOSIM) test. Bacterial and fungal community diversities were analyzed through the vegan and fossil R packages (v 3.4.3). Their structures were analyzed using QIIME (v 1.9.1) software and vegan package. To assess inter-kingdom equilibrium between the microbial communities, we computed a fungus-to-bacteria diversity ratio and a Basidiomycota-to-Ascomycota relative abundance ratio35 that were compared between patient groups using Wilcoxon signed-rank test. Relative abundances of OTUs were compared between groups as follow: for each feature (be OTU, genus, etc.) present in at least 10% of the patients, we first tested the distributional adequacy of the Poisson, the negative binomial and the Gaussian to the normalized non-zero data. Secondly, we considered the model presenting the best fitted distribution and we tested the differences between groups. When none of the three distributions fitted, Wilcoxon signed-rank test was used. Goodness-of-fit tests for discrete distributions were implemented through the vcd package (v 1.4.4). The Shapiro-Wilk test was used to test for normality as well as graphical check. Bioconductor packages DESeq. 2 (v 1.22.1) and edgeR (v 3.16.1) implemented differential analysis for the model with discrete distributions, while metagenomeSeq package (v 1.18.0) implemented differential analysis for the zero inflated Gaussian model72. Finally, each sample was divided into majority and minority populations (abundant population referring to an OTU relative abundance greater than or equal to 1%, and minority population as having a relative abundance of less than 1%), since their contribution may differ12,34. According to Carmody et al.15, we defined a dominant OTU as the most abundant OTU having a relative abundance higher or equal than twice the relative abundance of the second most abundant OTU for each sputum sample. We compared both abundant and minority populations as well as dominant OTUs between clinical status groups.
Inter-kingdom network analysis
To assess the relationships among microorganisms and to determine their underlying physio-pathological senses, joint analysis of fungal and bacterial relative abundance data was performed by estimating their correlation matrix.
According to Gevers et al.73, to test the correlation while accounting for the compositional nature of data (i.e. proportions that sum to 1), the permutation-renormalization bootstrap method (ReBoot) was used. We built both the null distribution of correlations from renormalized permuted data which represents the correlation structure arising purely from the compositionality, and the distribution of correlations from bootstrap sampled data referring to the confidence interval of the observed correlation. The compositional null distribution and the bootstrap distribution were compared by Z-test with the variance pooled from both distributions. The ccrepe Bioconductor package (v 1.10.0) implemented this ReBoot method, and was used here; it also allows applying network analysis up from 20 samples. The significance of association between features was based on the Spearman correlation as similarity measure (at the 0.05 level). We generated 1,000 bootstrap replicates.
In vitro co-cultures
A. fumigatus strain CBS144.89, and two clinical strains of S. oralis and S. mitis isolated from two CF patients and kindly provided by Pr. Lehours (Bacteriology Department, Bordeaux hospital) were used for in vitro experiments. These two Streptococcus isolates were selected according to their prevalence in CF upper airways and their role in CF inflammation and disease74,75,76,77. Co-cultures (10 million bacterial cells plus 1 million Aspergillus conidia) were conducted in brain-heart infusion (BHI) medium, with moderate shaking, at 37â°C, under aerobic conditions as recently described1. Growth of A. fumigatus alone or in co-culture with S. oralis or S. mitis was measured every 24â hours for five consecutive days. At each time point, a diluted aliquot of microbial BHI-mixture was plated on chocolate agar; fungal colony forming units (CFUs) were counted after a 24-h incubation period. The experiment was done in triplicate.
Identification of taxonomic traits independently associated with the clinical outcomes
Among a high number of OTUs (greater than the number of patient that we assimilated to a high-dimensional dataset), we determined which ones were significantly associated with clinical status by applying Least Absolute Shrinkage and Selection Operator (Lasso) regression78, after adjusting for the contribution of the other OTUs. Lasso linear or logistic regression models were used when the clinical status was measured, respectively, as a continuous (by FEV1) or binary (by CFPE) variable. In Lasso regression, a penalty on absolute coefficient size is added to the usual loss functions used for regression problems (mean squared error in the linear model and negative logistic log likelihood in the logistic model). In one hand, the introduction of a penalty often improves the prediction accuracy due to the bias-variance trade-off. On the other hand, if the amount of penalty is sufficiently large, some coefficients are shrunk to exact zero, thus estimation and OTU selection are simultaneously achieved. The Lasso technique has been adapted to hierarchical data. Rush et al.3 recently proposed a Phy-Lasso method that incorporates the phylogenetic tree structure characteristic of microbiome OTUs. This Phy-Lasso method applied a hierarchical model for each taxonomic level; the corresponding code being available from the author79. We added Phy-Lasso linear regression to the toolbox, and estimated the optimal amount of penalty from the data as follow. Because of the small sample size of the present study, we used leave-one-out cross-validation (LOO-CV) in which the model is fitted for n-1 patients and the predicted clinical status for the left-out patient is compared with his/her actual clinical status. This procedure is repeated n times, and the average agreement of the predicted and observed clinical status computed. Among values varying between high and low amounts of penalties, the good amount of penalty is chosen such that the predicted error is minimized. In linear regression (such as FEV1), we considered the predicted quadratic error. In logistic regression (such as CFPE), we considered the predicted classification error determined from a cutoff probability of 0.5. While the Lasso has excellent properties in dimensional reduction and estimation, it over-selects OTUs to reduce the prediction errors when a cross-validation method is used. To address this point, we applied the bootstrap-enhanced Lasso (Bolasso, based on intersecting bootstrapped Lasso estimations80), given its appealing asymptotic consistency properties and its simple implementation. In Bolasso, only OTUs frequently chosen by Phy-Lasso over bootstrap samples are selected, which improve results stability. We generated 200 bootstrap replicates. The frequency thresh-old of Bolasso was specified to be conservative (70%), to not miss any OTU associated with clinical status. PhyLasso associated function as well as all the codes used in this study are available at: https://github.com/psBiostat/MucoFong-study.git, and summarized in the Supplementary Table 2.
Revision of the climax-attack model (CAM) previously proposed in CF
Based on our results of inter-kingdom network analysis, in vitro co-cultures and feature selection based on bootstrap-enhanced Phy-Lasso, we proposed a CAM version enriched with the present mycobiome data (Fig. 6). In this new CAM version, we confirmed affiliations of bacteria previously proposed and documented4,7,13,14,31,32. Regarding mycobiome, we confirmed the previous Scedosporium genus assignation to Climax population, as it was negatively associated with both CFPE and elevated value of FEV1 (see feature selection in Fig. 5). Based on feature selection (Fig. 5), we proposed the same affiliation for Penicillium.
As Malassezia and Aspergillus are associated with CFPE (Fig. 5a), we proposed to affiliate these genera to the Attack population, which fit well with the positive network interaction (Fig. 2), and the co-cultures results between Aspergillus and Streptococcus (Fig. 4).
References
Tipton, L. et al. Fungi stabilize connectivity in the lung and skin microbial ecosystems. Microbiome 6, 12 (2018).
OâBrien, S. & Fothergill, J. L. The role of multispecies social interactions in shaping Pseudomonas aeruginosa pathogenicity in the cystic fibrosis lung. FEMS Microbiol. Lett. 364 (2017).
Rush, S. T., Lee, C. H., Mio, W. & Kim, P. T. The Phylogenetic LASSO and the Microbiome. ArXiv160708877 Q-Bio Stat (2016).
Quinn, R. A. et al. Ecological networking of cystic fibrosis lung infections. NPJ Biofilms Microbiomes 2, 4 (2016).
Kurtz, Z. D. et al. Sparse and compositionally robust inference of microbial ecological networks. Plos Comput. Biol. 11, e1004226 (2015).
Berry, D. & Widder, S. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front. Microbiol. 5, 219 (2014).
Whiteson, K. L. et al. The upper respiratory tract as a microbial source for pulmonary infections in cystic fibrosis. Parallels from island biogeography. Am. J. Respir. Crit. Care Med. 189, 1309â1315 (2014).
Bhatt, J. M. Treatment of pulmonary exacerbations in cystic fibrosis. Eur. Respir. Rev. Off. J. Eur. Respir. Soc. 22, 205â216 (2013).
Stenbit, A. E. & Flume, P. A. Pulmonary exacerbations in cystic fibrosis. Curr. Opin. Pulm. Med. 17, 442â447 (2011).
Bilton, D. et al. Pulmonary exacerbation: towards a definition for use in clinical trials. Report from the EuroCareCF Working Group on outcome parameters in clinical trials. J. Cyst. Fibros. Off. J. Eur. Cyst. Fibros. Soc. 10(Suppl 2), S79â81 (2011).
Goss, C. H. & Burns, J. L. Exacerbations in cystic fibrosis. 1: Epidemiology and pathogenesis. Thorax 62, 360â367 (2007).
Nguyen, L. D. N. et al. Effects of Propidium Monoazide (PMA) Treatment on Mycobiome and Bacteriome Analysis of Cystic Fibrosis Airways during Exacerbation. Plos One 11, e0168860 (2016).
Quinn, R. A. et al. Metabolomics of pulmonary exacerbations reveals the personalized nature of cystic fibrosis disease. PeerJ 4, e2174 (2016).
Quinn, R. A. et al. A Winogradsky-based culture system shows an association between microbial fermentation and cystic fibrosis exacerbation. ISME J. 9, 1024â1038 (2015).
Carmody, L. A. et al. Changes in cystic fibrosis airway microbiota at pulmonary exacerbation. Ann. Am. Thorac. Soc. 10, 179â187 (2013).
Tunney, M. M. et al. Lung microbiota and bacterial abundance in patients with bronchiectasis when clinically stable and during exacerbation. Am. J. Respir. Crit. Care Med. 187, 1118â1126 (2013).
Zemanick, E. T. et al. Inflammation and airway microbiota during cystic fibrosis pulmonary exacerbations. Plos One 8, e62917 (2013).
Filkins, L. M. et al. Prevalence of streptococci and increased polymicrobial diversity associated with cystic fibrosis patient stability. J. Bacteriol. 194, 4709â4717 (2012).
Fodor, A. A. 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 7, e45001 (2012).
Zhao, J. et al. Decade-long bacterial community dynamics in cystic fibrosis airways. Proc. Natl. Acad. Sci. USA 109, 5809â5814 (2012).
Goffard, A. et al. Virus and cystic fibrosis: rhinoviruses are associated with exacerbations in adult patients. J. Clin. Virol. Off. Publ. Pan Am. Soc. Clin. Virol. 60, 147â153 (2014).
Willner, D. et al. Case studies of the spatial heterogeneity of DNA viruses in the cystic fibrosis lung. Am. J. Respir. Cell Mol. Biol. 46, 127â131 (2012).
Lysholm, F. et al. Characterization of the viral microbiome in patients with severe lower respiratory tract infections, using metagenomic sequencing. Plos One 7, e30875 (2012).
Chotirmall, S. H. et al. Sputum Candida albicans presages FEV1 decline and hospital-treated exacerbations in cystic fibrosis. Chest 138, 1186â1195 (2010).
Amin, R., Dupuis, A., Aaron, S. D. & Ratjen, F. The effect of chronic infection with Aspergillus fumigatus on lung function and hospitalization in patients with cystic fibrosis. Chest 137, 171â176 (2010).
Willger, S. D. et al. Characterization and quantification of the fungal microbiome in serial samples from individuals with cystic fibrosis. Microbiome 2, 40 (2014).
Delhaes, L. et al. Prevalence, geographic risk factor, and development of a standardized protocol for fungal isolation in cystic fibrosis: Results from the international prospective study âMFIPâ. J. Cyst. Fibros. 18, 212â220 (2018).
Armstead, J., Morris, J. & Denning, D. W. Multi-country estimate of different manifestations of aspergillosis in cystic fibrosis. Plos One 9, e98502 (2014).
Middleton, P. G., Chen, S. C.-A. & Meyer, W. Fungal infections and treatment in cystic fibrosis. Curr. Opin. Pulm. Med. 19, 670â675 (2013).
Speirs, J. J., van der Ent, C. K. & Beekman, J. M. Effects of Aspergillus fumigatus colonization on lung function in cystic fibrosis. Curr. Opin. Pulm. Med. 18, 632â638 (2012).
Conrad, D. et al. Cystic fibrosis therapy: a community ecology perspective. Am. J. Respir. Cell Mol. Biol. 48, 150â156 (2013).
Lim, Y. W. et al. Clinical insights from metagenomic analysis of sputum samples from patients with cystic fibrosis. J. Clin. Microbiol. 52, 425â437 (2014).
Stevens, D. A. et al. Allergic bronchopulmonary aspergillosis in cystic fibrosisâstate of the art: Cystic Fibrosis Foundation Consensus Conference. Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 37(Suppl 3), S225â264 (2003).
Galand, P. E., Casamayor, E. O., Kirchman, D. L. & Lovejoy, C. Ecology of the rare microbial biosphere of the Arctic Ocean. Proc. Natl. Acad. Sci. USA 106, 22427â22432 (2009).
Sokol, H. et al. Fungal microbiota dysbiosis in IBD. Gut 66, 1039â1048 (2017).
Kramer, R. et al. Cohort Study of Airway Mycobiome in Adult Cystic Fibrosis Patients: Differences in Community Structure between Fungi and Bacteria Reveal Predominance of Transient Fungal Elements. J. Clin. Microbiol. 53, 2900â2907 (2015).
Delhaes, L. et al. The airway microbiota in cystic fibrosis: a complex fungal and bacterial communityâimplications for therapeutic management. Plos One 7, e36313 (2012).
Sanders, D. B. et al. Failure to recover to baseline pulmonary function after cystic fibrosis pulmonary exacerbation. Am. J. Respir. Crit. Care Med. 182, 627â632 (2010).
Heirali, A. A. et al. The effects of inhaled aztreonam on the cystic fibrosis lung microbiome. Microbiome 5, 51 (2017).
Bacci, G. et al. Pyrosequencing Unveils Cystic Fibrosis Lung Microbiome Differences Associated with a Severe Lung Function Decline. PloS One 11, e0156807 (2016).
Hogan, D. A. et al. Analysis of Lung Microbiota in Bronchoalveolar Lavage, Protected Brush and Sputum Samples from Subjects with Mild-To-Moderate Cystic Fibrosis Lung Disease. Plos One 11, e0149998 (2016).
Cuthbertson, L. et al. Respiratory microbiota resistance and resilience to pulmonary exacerbation and subsequent antimicrobial intervention. ISME J. 10, 1081â1091 (2016).
Bos, L. D. J., Meinardi, S., Blake, D. & Whiteson, K. Bacteria in the airways of patients with cystic fibrosis are genetically capable of producing VOCs in breath. J. Breath Res. 10, 047103 (2016).
Botterel, F. et al. Fungal and Bacterial Diversity of Airway Microbiota in Adults with Cystic Fibrosis: Concordance Between Conventional Methods and Ultra-Deep Sequencing, and Their Practical use in the Clinical Laboratory. Mycopathologia 183, 171â183 (2018).
Richardson, M., Bowyer, P. & Sabino, R. The human lung and Aspergillus: You are what you breathe in? Med. Mycol. 57, S145âS154 (2019).
Feigelman, R. et al. Sputum DNA sequencing in cystic fibrosis: non-invasive access to the lung microbiome and to pathogen details. Microbiome 5, 20 (2017).
Kim, S. H. et al. Global Analysis of the Fungal Microbiome in Cystic Fibrosis Patients Reveals Loss of Function of the Transcriptional Repressor Nrg1 as a Mechanism of Pathogen Adaptation. Plos Pathog. 11, e1005308 (2015).
Boutin, S. & Dalpke, A. H. Acquisition and adaptation of the airway microbiota in the early life of cystic fibrosis patients. Mol. Cell. Pediatr. 4, 1 (2017).
Dickson, R. P. et al. Spatial Variation in the Healthy Human Lung Microbiome and the Adapted Island Model of Lung Biogeography. Ann. Am. Thorac. Soc. 12, 821â830 (2015).
Venkataraman, A. et al. Application of a neutral community model to assess structuring of the human lung microbiome. mBio 6 (2015).
Lamoureux, C., Guilloux, C.-A., Beauruelle, C., Jolivet-Gougeon, A. & Héry-Arnaud, G. Anaerobes in cystic fibrosis patientsâ airways. Crit. Rev. Microbiol. 1â15, https://doi.org/10.1080/1040841X.2018.1549019 (2019).
Caverly, L. J. & LiPuma, J. J. Good cop, bad cop: anaerobes in cystic fibrosis airways. Eur. Respir. J. 52 (2018).
Rogers, G. B., Hoffman, L. R., Carroll, M. P. & Bruce, K. D. Interpreting infective microbiota: the importance of an ecological perspective. Trends Microbiol. 21, 271â276 (2013).
Hurley, M. N. & Smyth, A. R. Staphylococcus aureus in cystic fibrosis: pivotal role or bit part actor? Curr. Opin. Pulm. Med. 24, 586â591 (2018).
Hector, A. et al. Microbial colonization and lung function in adolescents with cystic fibrosis. J. Cyst. Fibros. Off. J. Eur. Cyst. Fibros. Soc. 15, 340â349 (2016).
Zemanick, E. T. et al. Airway microbiota across age and disease spectrum in cystic fibrosis. Eur. Respir. J. 50 (2017).
van Woerden, H. C. et al. Differences in fungi present in induced sputum samples from asthma patients and non-atopic controls: a community based case control study. BMC Infect. Dis. 13, 69 (2013).
Guillot, J., Hadina, S. & Guého, E. The genus Malassezia: old facts and new concepts. Parassitologia 50, 77â79 (2008).
Kerem, E. et al. Factors associated with FEV1 decline in cystic fibrosis: analysis of the ECFS patient registry. Eur. Respir. J. 43, 125â133 (2014).
Pages-Monteiro, L. et al. Strong incidence of Pseudomonas aeruginosa on bacterial rrs and ITS genetic structures of cystic fibrosis sputa. Plos One 12, e0173022 (2017).
Whelan, F. J. et al. Longitudinal sampling of the lung microbiota in individuals with cystic fibrosis. Plos One 12, e0172811 (2017).
Somayaji, R. et al. Long-term clinical outcomes of âPrairie Epidemic Strainâ Pseudomonas aeruginosa infection in adults with cystic fibrosis. Thorax 72, 333â339 (2017).
Russell, G. K., Gadhok, R. & Simmonds, N. J. The destructive combination of Scediosporium apiosperum lung disease and exuberant inflammation in cystic fibrosis. Paediatr. Respir. Rev. 14(Suppl 1), 22â25 (2013).
Staerck, C. et al. The secreted polyketide boydone A is responsible for the anti-Staphylococcus aureus activity of Scedosporium boydii. FEMS Microbiol. Lett. 364 (2017).
Vandeputte, P. et al. Draft Genome Sequence of the Pathogenic Fungus Scedosporium apiospermum. Genome Announc. 2 (2014).
Han, Z., Kautto, L. & Nevalainen, H. Secretion of Proteases by an Opportunistic Fungal Pathogen Scedosporium aurantiacum. Plos One 12, e0169403 (2017).
Krüger, W., Vielreicher, S., Kapitan, M., Jacobsen, I. D. & Niemiec, M. J. Fungal-Bacterial Interactions in Health and Disease. Pathog. Basel Switz. 8 (2019).
Budden, K. F. et al. Functional effects of the microbiota in chronic respiratory disease. Lancet Respir. Med., https://doi.org/10.1016/S2213-2600(18)30510-1 (2019).
Chiu, L. et al. Protective Microbiota: From Localized to Long-Reaching Co-Immunity. Front. Immunol. 8 1678 (2017).
Coron, N. et al. Toward the Standardization of Mycological Examination of Sputum Samples in Cystic Fibrosis: Results from a French Multicenter Prospective Study. Mycopathologia 183, 101â117 (2018).
Abbott, J. et al. What defines a pulmonary exacerbation? The perceptions of adults with cystic fibrosis. J. Cyst. Fibros. Off. J. Eur. Cyst. Fibros. Soc. 8, 356â359 (2009).
Paulson, J. N., Stine, O. C., Bravo, H. C. & Pop, M. Differential abundance analysis for microbial marker-gene surveys. Nat. Methods 10, 1200â1202 (2013).
Gevers, D. et al. The treatment-naive microbiome in new-onset Crohnâs disease. Cell Host Microbe 15, 382â392 (2014).
Enaud, R. et al. Intestinal Inflammation in Children with Cystic Fibrosis Is Associated with Crohnâs-Like Microbiota Disturbances. J. Clin. Med. 8 (2019).
Prevaes, S. M. P. J. et al. Development of the Nasopharyngeal Microbiota in Infants with Cystic Fibrosis. Am. J. Respir. Crit. Care Med. 193, 504â515 (2016).
Bernarde, C. et al. Impact of the CFTR-potentiator ivacaftor on airway microbiota in cystic fibrosis patients carrying a G551D mutation. Plos One 10, e0124124 (2015).
Maeda, Y. et al. Population structure and characterization of viridans group streptococci (VGS) including Streptococcus pneumoniae isolated from adult patients with cystic fibrosis (CF). J. Cyst. Fibros. Off. J. Eur. Cyst. Fibros. Soc. 10, 133â139 (2011).
Tibshirani, R. The lasso method for variable selection in the Cox model. Stat. Med. 16, 385â395 (1997).
Rush, S. T. A. The Phylogenetic LASSO and the Microbiome: Metagenomic Modeling in Fecal Microbiota Transplantation. (2017).
Bach, F. R. Bolasso: model consistent Lasso estimation through the bootstrap. in Proceedings of the 25th international conference on Machine learning - ICMLâ08 33â40 (ACM Press), https://doi.org/10.1145/1390156.1390161 (2008).
Acknowledgements
Laurence Delhaes has received a research grant from âVaincre La Mucoviscidoseâ (RF20160501626); she is member of ESGHAMI and NGS-MycA working groups of ESCMID and ECCM respectively. Genoscreen Company provided support in the form of Louise-Eva Vandenborghtâs salary (AST CT 2016-394). All the funders had no role in article design, decision to publish, or preparation of the manuscript. The authors would like to thank Katarzyna Hooks (postdoctoral researcher at University of Bordeaux) for having edited the English language of the text, and Binbin Xu (postdoctoral researcher at Inserm BPH U1219) for his help on testing the R and Fortran codes.
Author information
Authors and Affiliations
Consortia
Contributions
Mucofong Investigation group, M.F., L.E.V. and L.D. collected sputum samples and clinical data. L.E.V., N.C., R.E. and L.D. participated to microbial analysis including NGS P.S., L.E.V. realized the data analysis and figures/tables conception. P.S., L.E.V., F.F., M.A. and R.T. performed the statistical analyses including that of NGS R.T. and L.D. contributed to study design. P.S., L.E.V. and L.D. wrote the manuscript with significant contributions of other co-authors, especially N.C., R.E., M.A., T.S., P.B. and R.T. All authors read and gave approval on the final manuscript.
Corresponding author
Ethics declarations
Competing interests
Pr L. Delhaes is member of the scientific advisory board of âVaincre la Mucoviscidoseâ, and of the executive committee of ESGHAMI (ESCMID working group). She also did punctual consulting for pharmaceutical companies (Pfizer, Gilead, and MSD) and was invited to participate to scientific congresses as compensations. The other authors declare no potential conflict of interest.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
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 license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Soret, P., Vandenborght, LE., Francis, F. et al. Respiratory mycobiome and suggestion of inter-kingdom network during acute pulmonary exacerbation in cystic fibrosis. Sci Rep 10, 3589 (2020). https://doi.org/10.1038/s41598-020-60015-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-60015-4
This article is cited by
-
Affinity of Malassezia and Other Yeasts for Pulmonary Lipids
Mycopathologia (2025)
-
Metagenomics Applied to the Respiratory Mycobiome in Cystic Fibrosis
Mycopathologia (2024)
-
The respiratory microbiota alpha-diversity in chronic lung diseases: first systematic review and meta-analysis
Respiratory Research (2022)
-
Gut bacteriobiota and mycobiota are both associated with Day-28 mortality among critically ill patients
Critical Care (2022)
-
The treatment of Qibai Pingfei Capsule on chronic obstructive pulmonary disease may be mediated by Th17/Treg balance and gut-lung axis microbiota
Journal of Translational Medicine (2022)