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

The Airway Microbiota in Cystic Fibrosis: A Complex Fungal and Bacterial Community—Implications for Therapeutic Management

  • Laurence Delhaes ,

    Laurence.delhaes@pasteur-lille.fr

    Affiliations Center for Infection and Immunity of Lille (CIIL), Institut Pasteur de Lille, Biology and Diversity of Emerging Eukaryotic Pathogens (BDEEP), BP 245, Lille, France, INSERM U1019, Lille, France, UMR CNRS 8402, Lille, France, Department of Parasitology-Mycology, Faculty of Pharmacy, University Lille Nord de France, EA4547, Lille, France, Department of Microbiology, Lille Hospital, Faculty of Medicine, Lille, France

  • Sébastien Monchy,

    Affiliations LMGE, Laboratoire Microorganismes: Génome et Environnement, UMR CNRS 6023, Clermont Université, Blaise Pascal, BP 80026, Aubière, France, Université Lille Nord de France, Université du Littoral Côte d'Opale, ULCO, Laboratoire d'Océanologie et de Géoscience (LOG), UMR CNRS 8187, Wimereux, France

  • Emilie Fréalle,

    Affiliations Center for Infection and Immunity of Lille (CIIL), Institut Pasteur de Lille, Biology and Diversity of Emerging Eukaryotic Pathogens (BDEEP), BP 245, Lille, France, INSERM U1019, Lille, France, UMR CNRS 8402, Lille, France, Department of Parasitology-Mycology, Faculty of Pharmacy, University Lille Nord de France, EA4547, Lille, France, Department of Microbiology, Lille Hospital, Faculty of Medicine, Lille, France

  • Christine Hubans,

    Affiliation Genoscreen, Institut Pasteur of Lille, Lille, France

  • Julia Salleron,

    Affiliation Department of Biostatistics, Lille Hospital, Faculty of Medicine, Lille, France

  • Sylvie Leroy,

    Affiliation Department of Pneumology and Immuno-Allergology, CRCM adulte, Calmette Hospital, Lille, France

  • Anne Prevotat,

    Affiliation Department of Pneumology and Immuno-Allergology, CRCM adulte, Calmette Hospital, Lille, France

  • Frédérick Wallet,

    Affiliation Department of Microbiology, Lille Hospital, Faculty of Medicine, Lille, France

  • Benoit Wallaert,

    Affiliation Department of Pneumology and Immuno-Allergology, CRCM adulte, Calmette Hospital, Lille, France

  • Eduardo Dei-Cas,

    Affiliations Center for Infection and Immunity of Lille (CIIL), Institut Pasteur de Lille, Biology and Diversity of Emerging Eukaryotic Pathogens (BDEEP), BP 245, Lille, France, INSERM U1019, Lille, France, UMR CNRS 8402, Lille, France, Department of Parasitology-Mycology, Faculty of Pharmacy, University Lille Nord de France, EA4547, Lille, France, Department of Microbiology, Lille Hospital, Faculty of Medicine, Lille, France

  • Telesphore Sime-Ngando,

    Affiliation LMGE, Laboratoire Microorganismes: Génome et Environnement, UMR CNRS 6023, Clermont Université, Blaise Pascal, BP 80026, Aubière, France

  • Magali Chabé,

    Affiliations Center for Infection and Immunity of Lille (CIIL), Institut Pasteur de Lille, Biology and Diversity of Emerging Eukaryotic Pathogens (BDEEP), BP 245, Lille, France, INSERM U1019, Lille, France, UMR CNRS 8402, Lille, France, Department of Parasitology-Mycology, Faculty of Pharmacy, University Lille Nord de France, EA4547, Lille, France

  • Eric Viscogliosi

    Affiliations Center for Infection and Immunity of Lille (CIIL), Institut Pasteur de Lille, Biology and Diversity of Emerging Eukaryotic Pathogens (BDEEP), BP 245, Lille, France, INSERM U1019, Lille, France, UMR CNRS 8402, Lille, France, Department of Parasitology-Mycology, Faculty of Pharmacy, University Lille Nord de France, EA4547, Lille, France

Abstract

Background

Given the polymicrobial nature of pulmonary infections in patients with cystic fibrosis (CF), it is essential to enhance our knowledge on the composition of the microbial community to improve patient management. In this study, we developed a pyrosequencing approach to extensively explore the diversity and dynamics of fungal and prokaryotic populations in CF lower airways.

Methodology and Principal Findings

Fungi and bacteria diversity in eight sputum samples collected from four adult CF patients was investigated using conventional microbiological culturing and high-throughput pyrosequencing approach targeting the ITS2 locus and the 16S rDNA gene. The unveiled microbial community structure was compared to the clinical profile of the CF patients. Pyrosequencing confirmed recently reported bacterial diversity and observed complex fungal communities, in which more than 60% of the species or genera were not detected by cultures. Strikingly, the diversity and species richness of fungal and bacterial communities was significantly lower in patients with decreased lung function and poor clinical status. Values of Chao1 richness estimator were statistically correlated with values of the Shwachman-Kulczycki score, body mass index, forced vital capacity, and forced expiratory volume in 1 s (p = 0.046, 0.047, 0.004, and 0.001, respectively for fungal Chao1 indices, and p = 0.010, 0.047, 0.002, and 0.0003, respectively for bacterial Chao1 values). Phylogenetic analysis showed high molecular diversities at the sub-species level for the main fungal and bacterial taxa identified in the present study. Anaerobes were isolated with Pseudomonas aeruginosa, which was more likely to be observed in association with Candida albicans than with Aspergillus fumigatus.

Conclusions

In light of the recent concept of CF lung microbiota, we viewed the microbial community as a unique pathogenic entity. We thus interpreted our results to highlight the potential interactions between microorganisms and the role of fungi in the context of improving survival in CF.

Introduction

The human respiratory tract represents the major portal of entry for numerous microorganisms, primarily those occurring as airborne particles such as viral and bacterial entities, or fungal spores. Microorganism characteristics coupled with the local host immune response will determine whether they will be cleared or adhere and colonize the airways leading to acute or chronic pulmonary disease.

In cystic fibrosis (CF), mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene result in defective mucociliary clearance and, as a consequence, lead to the production of thick and sticky bronchial mucus, which facilitates the entrapment of airborne viruses, bacteria and fungal spores and provides a suitable environment for the growth of these microorganisms. In addition to bacteria, which are well known to cause recurrent exacerbations of CF-associated pulmonary disease and often determine the vital prognosis of patients [1], many fungi also colonize the respiratory tract of CF patients [2][5], although their involvement in respiratory infections remains controversial and largely unsolved [6], [7]. As an alternative to conventional methods (direct examination and microbiological cultures), new molecular techniques such as pyrosequencing, have been developed to offer a more complete view of the microbiota. In human samples, these molecular methods can distinguish microorganisms difficult to identify and/or those that are refractory to culture (such as Pneumocystis jirovecii, Scedosporium apiospermum, atypical mycobacteria, etc.), as well as new or as yet unknown pathogens [8][10]. The metagenomic approach has been recently used for the identification of human bacterial populations in the gut as well as in the mouth saliva and skin of patients [11][15]. In addition, metagenomic studies have been successfully used for providing an overview of community composition with semi-quantitative information [13], [14], [16][19]. Some studies have been published on the human respiratory tract, but only few have specifically focused on microbial diversity in CF [16][18], [20][25].

In the present study, we applied a molecular approach by pyrosequencing variable regions of bacterial 16S rDNA and fungal ITS2 genes in sputum samples from CF patients. Our aims were to explore the fungal and bacterial assemblages in CF patients to achieve a better understanding of species/taxon diversity and population dynamics of the microbiota, and their relevance for the clinical course of pulmonary disease in CF.

Results

Samples and patients

We prospectively collected eight sputum samples from four CF adult patients (median age of 29.5 years; Q1, 24.5; Q3, 34) who were all part of a long-term follow-up program at Lille's Adult CF center. Two temporal sputum samples were collected from each clinically stable patient with a sampling interval of 1 year; for Patient 4 the sampling interval was only three months (Tables 1 and 2). Three out of the four patients were homozygous or heterozygous for the ΔF508 mutation.

thumbnail
Table 1. Clinical data and treatment from CF patients included in the study.

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

thumbnail
Table 2. Microbiological data from CF patients included in the study.

https://doi.org/10.1371/journal.pone.0036313.t002

Results of the Shwachman-Kulczycki score (S-K score), body mass index (BMI), forced vital capacity (% of predicted FVC), and forced expiratory volume in 1 s (% of predicted FEV1) expressed as medians (Q1, Q3) were 55.0 (50.0; 90.0), 20.4 (17.1; 22.85), 75.0 (48.0; 111.5), and 53.0 (41.0; 90.0), respectively. None of the CF patients had pancreatic alterations. All patients had bacteria present in their sputa, as determined by culturing, but none were known to be chronically infected with Staphylococcus aureus or Burkholderia cepacia complex. Using conventional methods, we identified mucoid and non-mucoid Pseudomonas aeruginosa, meticillin-sensitive Staphylococcus aureus, Haemophilus influenzae, and Alcaligenes xylosoxidans (Table 2). No oropharyngeal flora was detected, and no mycobacteria were isolated. Two patients (Patients 1 and 4) were colonized by P. aeruginosa and treated with azithromycin, as recommended for chronic P. aeruginosa infections. Both patients had severe airway disease as assessed by the S-K score (S-K score ≤50%), BMI (under 16 for Patient 4) and standard spirometry (FEV1<50% of predicted FEV1) (Table 1). Patients 2, 3 and 4 were treated with inhaled corticosteroids, and Patient 4 received systemic corticoids (Table 1). Only Patient 1 received long-term itraconazole treatment (600 mg per day for 6 months) for allergic bronchopulmonary aspergillosis (ABPA).

Regarding fungi, Candida albicans and Geotrichum sp., and two filamentous species, Aspergillus fumigatus and Aspergillus flavus, were isolated from sample cultures. Aspergillus nidulans, Aspergillus terreus, S. apiospermum, Scedosporium prolificans, or Exophiala dermatitidis were not isolated. In addition, P. jirovecii colonization was retrospectively diagnosed in three out of four patients. Both sputum samples of Patient 2, as well as one of Patient 1 (sample 2) and Patient 4 (sample 1) were nested PCR-positive for P. jirovecii (Table 2) [26]. Aspergillus DNA was detected using an ultrasensitive real-time PCR assay [27] in five of the eight sputum samples (Table 2).

Overall richness and diversity of microbial community evaluated from pyrosequences

We obtained a total of 326,277 sequences from samples 1 and 2 of Patients 1, 2, 3 and 4 using primers for the prokaryote 16S rDNA gene, a result in agreement with recent published data [17] (Figure 1A). Using the fungus-specific ITS2 primers, we obtained a total of 133,317 sequences from these samples (Figure 1B). Once primer, tag and key fragments were removed, 93% and 85% of the sequences had lengths greater than 450–500 bp and 300–450 bp for the 16S rDNA and ITS2 loci, respectively.

thumbnail
Figure 1. Rarefaction curves.

These curves are representing the numbers of OTUs with respect to the number of pyrosequence reads obtained from each patient at different sampling times and using the two set of primers targeting prokaryotic 16S rDNA (A) and fungal ITS2 (B) loci.

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

The pyrosequences that presented similarities with sequences available in databases but that could not be classified to at least the level of kingdom using BLASTN and MEGAN software were designated as “not assigned” and excluded from subsequent diversity analyses. For each sputum sample, these sequences represented less than 5% of the 16S rDNA or ITS2 sequences included in analyses, except for Patient 1-sample 2 and Patient 2-sample 2, which showed 9.7% and 8.4% of non-assigned 16S rDNA and ITS2 sequences, and 29.4% of non-assigned ITS2 sequences, respectively. Pyrosequences without any similarity with sequences available in databases were designated as “no hits”, and may represent species not yet represented in databases. Unsurprisingly, there were more ‘no hits’ for ITS2 pyrosequences than in 16S rDNA pryosequences (Tables 3 and 4, Figures S1, S2, S3, S4), due to the massive amount of data available in the Silva SSU rDNA database compared to the ITS2dbScreen database created expressly for the present analysis (see Materials and Methods). Un-represented organisms in sequence databases have already been described as a limitation in the ability to placing reads in the phylogeny [28].

thumbnail
Table 3. Number of 16S-pyrosequencing reads assigned to each taxonomic group of Bacteria.

https://doi.org/10.1371/journal.pone.0036313.t003

thumbnail
Table 4. Number of ITS-pyrosequencing reads assigned to each taxonomic group of Fungi.

https://doi.org/10.1371/journal.pone.0036313.t004

For all patients and samples except one (i.e. Patient 1-sample 1 for the ITS2 locus), the rarefaction curves for the number of OTUs per pyrosequence reads reached a plateau, indicating that almost all OTUs present in each sample were detected. The apparent observed diversity was higher for the prokaryote 16S rDNA locus (Figure 1A) in comparison to the fungus-specific ITS2 locus (Figure 1B).

Calculated to analyze microbial diversity, Chao1 richness estimator values corroborated rarefaction curves, confirming high bacterial diversity (Figure 2A). Bacterial diversity was higher in samples from Patients 2 and 3 than in samples from Patients 1 and 4. Fungal diversity showed a similar pattern (Figure 2B).

thumbnail
Figure 2. Relation between species richness and clinical status (A) or lung function (B).

Total richness of prokaryotic and fungal communities from each patient-sample was expressed using the Chao1 richness estimator; each spot size is proportional to the corresponding Chao1 value. The clinical status is expressed as S-K score and BMI in Figure 2A, while lung function is expressed as FEV1 and FVC values in Figure 2B. Given to the absence of S-K score value from Patient 2-sample 2 (Table 1), this spot is missing in Figure 2A. Bacterial and fungal Chao1 values corresponding to Patient 1, Patient 2, Patient 3, and Patient 4 are represented in blue-, green-, red- and yellow-edged spots, respectively. Dark and light colour intensity is corresponding to the first and second sampling dates of each patient, respectively. Dark grey and light grey are corresponding to fungal and bacterial Chao1 richness values, respectively.

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

Comparison of the culture and pyrosequencing results

Our results confirmed the new genera recently identified in CF patients [8], [16], [17], [20], [23], [24], [29][31], with Gemella sp. being found in sputum samples of 3 out of 4 patients (Table 3, Figures S1AS4A). The most represented genera identified in the present study were Pseudomonas, Streptococcus, Haemophilus, and anaerobes, in agreement with published data [1], [16], [17], [20], [23], [24], [30], [31]. Bacteria belonging to Pseudomonas, Streptococcus, Prevotella, Fusobacterium, Haemophilus, Veillonella, and Porphyromonas genera were isolated as recently reported in either sputum or BAL samples from CF patients (Table 3) [16], [18], [20], [21], [23], [30].

High fungal diversity was also observed in samples, with more than 60% of the fungal species or genera obtained in pyrosequencing not identified by mycological cultures (Tables 1, 2 and 3, Figures S1BS4B). Among the 24 species or genera of micromycetes identified by pyrosequencing, only four were also isolated in culture. Using the metagenomic approach, we identified additional species, especially within the genera Candida and Aspergillus, which are microorganisms known to be involved in pulmonary diseases or infectious diseases [3], [4], [32][45]. Geotrichum sp., which represents an important pathogenic genus with teleomorphs assigned to Dipodascus [35], was identified to the family level using the pyrosequencing method (Table 4, Figure S1B), due to the stringent parameters chosen to assign ITS2 sequences. Phylogenetically distinct from the A. fumigatus cluster, non-fumigatus Aspergillus isolates were molecularly identified in Patient 2-sample 1 (Figure S5), in agreement with RT-PCR results (Table 2), which detects the mitochondrial DNA of A. fumigatus as well as other species such as A. flavus [27]. Aspergillus lentulus, which represents a species difficult to differentiate from A. fumigatus solely based on phenotype criteria but has decreased susceptibility to azoles [45][47], was isolated from Patient 3-sample 2 (Figure S3B). The major expected advance from the high-throughput sequencing method was its ability to identify difficult-to-culture micromycetes, such as P. jirovecii or Malassezia sp. Although nested-PCR targeting P. jirovecii was positive in 4 sputum samples (Table 2), high-throughput sequencing did not identify this fungus, probably because there is only one copy of the ITS2 locus in the Pneumocystis genome [48]. Malassezia restricta was identified in all patients (Table 4, Figures S1BS4B), and Malassezia globosa and Malassezia sympodialis were molecularly identified in Patient 2 samples (Figure S2B). These results are consistent with the lipophilic nature that characterizes these yeasts and prevents their growth in standard culture media because they require an exogenous source of fatty acids [35]. Malassezia spp are frequently found in the skin of warm-blooded vertebrates, and they are currently recognized as emerging infectious pathogens [4], [35]. Recently, Malassezia has been identified in sputum samples from CF patients [25].

Since results from the conventional and high-throughput sequencing techniques concurred, the pyrosequencing method was used to identify dominant taxa, estimate their diversity, and analyze their temporal distribution, based on data obtained from both bacterial and fungal primers.

Fungal diversity and associated patterns of bacterial flora

The relative amounts of each species were estimated from the number of assigned pyrosequences, and were represented by pie charts whose diameters are proportional to the number of assigned sequences (Figures S1, S2, S3, S4). According to recent publications [19], [22], the number of pyrosequences obtained corresponds to the number of genome copies present in the sputum sample. The median (Q1, Q3) number of microorganism genera per sputum sample was 3.5 (3; 7.5) micromycetes and 6.5 (5; 13.5) bacteria; these results were comparable to those of previous studies [1], [8], [9], [29], [31]. We observed bacterial diversity similar to that recently reported in CF patients using molecular methods [1], [8], [17], [20], [23], [24], [29][31], with anaerobic bacteria representing a large proportion of the detected species (ranging from 2% to 50% of total pyrosequences for Patient 1-sample 1 and Patient 3-sample 2, respectively). For the kingdom Fungi, the 133,317 pyrosequences corresponded to 30 species or genera, including 24 micromycetes and 6 basidiomycetous macroscopic fungi. Among them, filamentous fungi belonging to the genera Aspergillus (in particular Aspergillus fumigatus), and Penicillium have already been described as pathogens in CF patients [2][5], [35], [45]. Candida albicans and species from the Candida parapsilosis complex have been recently recognized as medically important organisms colonizing CF patients [2], [4], [40], [42], [43]. Although their clinical relevance is still matter of debate, long-term persistence of Candida strains have been described in CF respiratory tracts [4], [40], [42], [43], [49]. Clavispora is a yeast genus that includes Clavispora lusitaniae (teleomorph of Candida lusitaniae); this ascomycete has already been isolated from sputa [36], [40].

A significant proportion of other species were either fungi reported in asthma (Didymella exitialis, Penicillium camemberti), allergy diseases (Aspergillus penicilloides and Eurotium halophilicum) [32], [33], [37], [39], [50], or infectious diseases (Kluyveromyces lactis, Malassezia sp., non-neoformans Cryptococci, Chalara sp.) [34][36], [41]. The other species or genera represented environmental taxa, either described as wood-inhabiting fungi common in temperate regions of the Northern Hemisphere, such as cereal pathogens associated with air pollution or food preparation processes [51][53]. In addition, macromycetes living on the wood of tree species common in Europe including in northern France, were identified in sputum samples from Patients 1, 2, and 4 (Figures S1B, S2B, S4B), probably corresponding to the signature of the outdoor environment that the patients are exposed to [54][56].

A growing number of studies has revealed that bacterial [1], [8], [17], [20], [23], [24], [29], [30], [57], [58] and fungal [9], [25] community compositions vary greatly among patients. Diversity at sub-species levels has also been described in CF, mainly for bacteria such as P. aeruginosa [17], [57], [58], and to a lesser degree for fungi [9], [25] or viruses [18]. Therefore, the microbial community was currently considered to be a unique pathogenic entity with potential interactions between microorganisms [17], [59][61]. From the perspective of this microbiota concept, we phylogenetically analyzed the diversity of the main fungi and bacteria identified by pyrosequencing, considered the taxon composition of each sample with potential interactions between fungi and bacteria, and investigated its clinical significance.

Population dynamics of the microbial communities in CF airways and clinical relevance

Although we observed lower diversity in CF airways than in other communities such as human skin, gut, or water microbiomes [12], [14], [19], reduced diversity and richness of fungal and bacterial communities were associated with poor clinical status, as evaluated from the S-K score and BMI values, and decreased lung function, as measured by FEV1 and FVC values, in CF (Figure 2). Chao1 values were statistically correlated with values of the S-K score, BMI, FVC, and FEV1 (p = 0.046, 0.047, 0.004, and 0.001, respectively, for Chao1 indices of fungal species, and p = 0.010, 0.047, 0.002, and 0.0003, respectively, for Chao1 values of bacterial species). Moreover, fewer fungus species were detected in sputum samples with lower FEV1 values; the correlation trended toward significance (p = 0.062). In parallel, a significant correlation between the number of bacteria species detected per sputum sample and values of S-K scores, BMI, FVC and FEV1 was observed (p = 0.0005, 0.03, 0.0003, 0.016 respectively) in agreement with published data [23], [58].

As previously observed [1], [30], anaerobes were significantly isolated in association with Pseudomonas, when comparing the relative amount of reads in each sample (p = 0.0003). Using a phylogenetic method, most Pseudomonas pyrosequences proved to be highly similar and clustered with sequences of P. aeruginosa strains isolated from CF patients or endotracheal tube biofilms (Figure S6) [62]. They also exhibited high infraspecific diversity, in agreement with previous results [57], [58]. The next most common bacterial genus was Streptococcus, of which the Streptococcus milleri group (SMG) has been isolated in CF [1], [20], [30], [63], linked to pulmonary exacerbations [1], [63], and demonstrated to produce quorum-sensing signal molecules [63]. SMG-related Streptococcus were identified in Patient 2-sample 2 and Patient 3-sample 2 (Figure S7). These phylogenetically identified SMG members (their sequences clustered with the SMG sequences of Streptococcus anginosus, S. intermedius, and S. constellatus in figure S7, using Neighbor-joining approach) were not numerically dominant compared to other clades, in agreement with the clinically reported absence of pulmonary exacerbation. Phylogenetic analysis of pyrosequences corresponding to the genera Haemophilus and Malassezia did not provide any new information compared to the pyrosequencing analysis using BLASTN and MEGAN software.

Using the same phylogenetic method, we observed diversity among genotypes of C. albicans, C. parapsilosis and A. fumigatus, with the same genotypes shared between patients, and/or genotypes that persisted over time within patients (Figures S5, S8, S9), in agreement with published data [40], [64]. Candida albicans and C. parapsilosis represented typical dominant yeasts isolated from CF sputa [2], [4], [5], [40], [65] for which we observed diversity similar to that already reported (mainly a single predominant C. albicans genotype) [40]. Regarding the aspergilli, samples were mainly composed of A. fumigatus as shown in the phylogenetic analysis (Figure S5), except for Patient 2-sample 1 in which the Aspergillus genus showed a high diversity, including non-fumigatus Aspergillus (Figure S2B, and sequences in dark green in Figure S5). Among A. fumigatus pyrosequences of Patient 3, one genotype was predominant in both samples of the patient, in agreement with previous studies that have demonstrated the emergence of a single genotype from a multiple-genotype population when chronic infection has been established [64], [66], [67].

Several recent taxonomic studies have identified cryptic species within key clinical morpho-species of both yeast and molds, including the C. parapsilosis complex, the A. fumigatus species complex and the S. apiospermum complex, which are particularly involved in CF lung colonization [46], [47], [68][70]. Here, we were able to differentiate C. metapsilosis genotypes from C. parapsilosis genotypes (Figure S7), as well as A. lentulus from A. fumigatus (Figure S5). This may have therapeutic implications given the different antifungal susceptibility profiles of these species [40], [45], [71].

The relative amounts (expressed as percentage of reads in each sample) of C. albicans or A. fumigatus were not statistically correlated with any bacterial taxon, neither anaerobic bacteria, nor Pseudomonas, nor Streptococcus. Nevertheless C. albicans was frequently associated with P. aeruginosa (80% of cases), which may be related to its recently proposed core status [21] and the bidirectional signalling pathway observed [for review60], [72][75]. Patient 3-sample 2 had a high number of A. fumigatus pyrosequences (23.6%) and this was associated with a predominance of Streptococcus (44.4%), which is a genus known to produce quorum-sensing molecules and to induce interactions between microorganisms, particularly among SMG members isolated from CF patients [63]. Regarding the temporal changes in the microbiota in each patient, we observed similar patterns, namely a disappearance of or major decrease in some bacterial genera recently described as members of the “core” pulmonary microbiome [76] and known to be a part of the oral bacterial community coupled with the emergence of more pathogenic bacteria (loss of Prevotella in Patients 1, 2, 3, and 4, Gemella in Patients 2, 3, and 4, Veillonella in Patient 3 associated with the emergence of P. aeruginosa in Patients 1 and 3, and of H. influenzae in Patient 2). Similarly, C. albicans, A. fumigatus or A. lentulus were detected in the second samples of Patient 1 and Patients 2, 3, and 4, respectively, while fungal species or genera known to be poorly pathogenic disappeared.

On the whole, our study highlights the correlation between richness and diversity of fungal and bacterial microbiota (Figure 2). We therefore suggest that “colonization resistance” occurs in CF lower airways, similar to what has been proposed to explain the exclusion of pathogenic species from the gut and the mouth by the presence of a specific microbiota [59], [77][79]. This phenomenon may be due to a range of factors and microbe-microbe interactions, including the presence of “synergens” described as enhancing the pathogenicity of the whole microbiota [78], [79], that will be discussed below.

Discussion

Given the recent evidence that fungi may be of clinical relevance in the decline of CF lung function, associated with co-colonization of fungi and bacteria [5], [22], [49], [80], we coupled fungal analysis to the characterization of bacterial flora in sputum samples from CF adults using the pyrosequencing technique. We acknowledge that the present CF cohort is small but comparable to sample size recently published (from 4 to 14 sputum samples [17], [23], [76]), lacks a specific control group — which is difficult to choose [7], e.g. there can be extensive overlap of bacterial membership between the pulmonary microbiome of healthy subjects and patients with or without COPD [76] —, and probably is not completely representative of the full spectrum of CF pulmonary pathology. However, this pyrosequencing-based study of fungal and bacterial communities in the human airway confirmed the recently reported bacterial diversity (including anaerobes) in CF patients [8], [17], [20], [23], [24] as well as in COPD patients using BAL [76], and revealed complex fungal biota in sputum samples, with a majority of the fungal species or genera obtained by pyrosequencing not identified in cultures, most of them known to be pathogens. Using phylogenetic tools, we also found infraspecific diversity in C. albicans, C. parapsilosis and A. fumigatus similar to previous published data [40], [42], [64], [66], [67]. In parallel, cryptic and new unculturable (or difficult to grow in vitro) species have also been identified, most of them described as human pathogens. In agreement with a recent oligonucleotide array analysis [9], we showed that fungal microbiota colonizing the lower airways of CF patients is more diverse and complex than previously estimated with culture methods. Therefore, culture methods are probably inadequate for assessing CF respiratory fungal microbiota, although culture methods can be improved with increased standardization [3], [81] and are still required to determine drug susceptibility. Moreover, we have evidence that poor clinical status is associated with lower taxon diversity and richness in fungal and bacterial communities (decrease in S-K scores, BMI, FVC, and FEV1 values significantly associated with low Chao1 indices).

Our findings add support to (i) the pathogenicity of species derived from the oral cavity and usually considered as clinically insignificant such as anaerobes and SMG members, even if their role in infection and inflammation needs to be further elucidated [1], [8], [17], [20], [23], [29], [31], [63], [79], and (ii) the complex interaction between typical pathogens and microbiota, such as the association between P. aeruginosa and anaerobes [20], [30], [58], [59]. Since C. albicans and C. parapsilosis can also be part of oral flora, these yeasts can migrate from the oral environment, colonize and persist within the lower airways of CF patients [40], as proposed for bacteria [23]. Although the implication of C. albicans in the decline of CF lung function has been recently suggested [49], the clinical relevance of yeasts is still matter of debate, and remains to be confirmed. Given the airborne transmission of molds such as A. fumigatus, opportunistic molds represent the most common agents of fungal colonization and/or infection of the CF airways. Among them, A. fumigatus has been reported more and more frequently since the 2000s [3][5], [9], and is associated with clinical significance in CF [80] and modification in the population of genotypes during chronic colonization [64], [66], [67]. Fungal colonization (especially repeated or chronic colonization) may have a substantial impact on the development of CF pulmonary disease [43], [49], [80], but more studies are required to determine this fungal risk, especially in light of the concomitant bacterial biota.

Given the relationship between decreased microbiota diversity and poor clinical status, we hypothesize that the composition of the microbial community in CF airways is the result of dynamics that take into account the different microorganisms present as an entity with interactions at the intra-species level as well as at the inter-species level. This is somewhat similar to the constitution of oral microbial consortia for which the potential for infection or co-infection is realized when potential pathogens find suitable community partners and local conditions (host response, adhesion receptors, biofilm formation) [59]. It is well known that the heterogeneity of mucus composition in CF provides suitable conditions for chronic infection by a wide range of microorganisms. In particular, recent data indicate that reduced oxygen tension in CF lung promotes the growth of P. aeruginosa [82], [83], as well as other anaerobic bacteria [1], [30]. Candida albicans can also grow under anaerobic conditions, showing mating type modifications that may promote yeast development [84], [85].

In addition to local physical conditions, both bacteria and fungi possess the ability to form biofilm consortia [63], [82], [83], [86][88]. In this context, direct and indirect microbe-microbe interactions have been well documented, particularly those involving the major prokaryotic CF pathogen: P. aeruginosa (for review, see [89]). Pseudomonas aeruginosa can produce substances that modulate growth of other microorganisms, in particular fungi. Pseudomonas aeruginosa and C. albicans can coexist, or have an antagonistic influence as recently proposed between P. aeruginosa and A. fumigatus [60], [72], [90], [91]. Moreover, C albicans produces farnesol, which in addition to its quorum-sensing function regulating yeast morphogenesis and its ability to modify P. aeruginosa growth, also reduces competition from other fungi such as A. fumigatus [92]. Because a large proportion of bacteria have been shown to synergistically affect CF disease outcome by modifying the expression of virulence genes [79], it may not be surprising to find evidence of such synergistic interactions within the fungus community.

Thus, analyzing microbial diversity in polymicrobial samples such as CF sputa is the type of study for which metagenomic methods have been recently proposed [16], [18], [21], [78]. Our results, along with others [16][18], [24], [76], demonstrate the utility of high-throughput sequencing in identifying microorganisms to investigate the microbiome associated with chronic pulmonary diseases, such as CF or COPD. These results now need to be confirmed by further pyrosequencing studies, especially in large multicenter studies that will lead to a better understanding of the dynamics of such CF microbiota.

In the near future, microbiota complexity should be taken into account to analyze host-microbe interactions, which are bi-directional and probably not limited to the direct contact lung area (as proposed in ref. [61]). The analysis of CF pulmonary disease and its management should be reassessed in light of these interactions. This concept of CF lung microbiota has emerged recently from the scientific community working on the microbiology of the CF respiratory tract [61], [78], and entails coupling environmental microbiological approaches with community ecology analyses (i.e. analyzing species richness and relative species abundance in terms of either spatial or temporal distribution and dividing species into core and satellite groups, at an ecologically relevant spatial scale) [23], [31], [61], [76]. Furthermore, these molecular results should be combined with biological models, such as biofilm models or in vivo planktonic cultures as recently proposed, in order to elucidate the possible interaction between bacteria and fungi detected here [59], [63], [82], [83], [86][93].

Few culture-independent strategies have been developed to evaluate bacterial [1], [8], [16], [17], [20], [23], [24], [29], [31], fungal [9], and viral [18] diversities in sputum samples from CF patients. Thus, new high-throughput sequencing approaches offer more exhaustive coverage of the sequences present in PCR products, in particular when the new generations of automatic sequencers, such as the GS FLX Titanium System, are used. Compared to terminal restriction fragment length polymorphism (T-RFLP) analysis, high-throughput sequencing methods more accurately identify pathogens, because they are based on sequences instead of amplicon sizes that can be shared between two or more species. For example, S. sanguinis, S. parasanguinis and S. salivarius all generate a 576 bp T-RFLP fragment [1].

Nevertheless, these molecular strategies can have some confounding factors. One important drawback due to the basic PCR approach is the incapacity to reflect the viability of the microorganisms detected by DNA amplification, unless samples are pre-treated (with, for example, propidium monoazide, [61]). Furthermore, DNA extraction from clinical samples is the first crucial step in ensuring faithful molecular detection. Non-homogenous lysis of bacterial and fungal cells, which are known to require strong lysis in order to extract DNA, may introduce biases as in any method based on DNA amplification [94], [95]. In addition to DNA extraction efficiency that can vary between microorganisms, the choice of the PCR protocol, from primer design to the number of PCR cycles, can affect the results. In contrast to specific PCR targeting a specific pathogen, high-throughput methods as well as cloning/sequencing techniques, are based on amplification with primers targeting conservative regions of microorganism DNA. These techniques can thus identify any microorganism that is reasonably abundant within the sample without the need for prior prediction of which species may be present. This universal-primer approach leads to the preferential amplification of the most prevalent flora. This bias may explain the negative pyrosequencing results for P. jirovecii, which may be present in small numbers since only nested-PCR was positive (not detected upon direct examination). The clonal Sanger-sequencing approach would be more suitable than pyrosequencing methods for identifying microorganisms in relatively low abundance [8]. Improvements in amplicon length with the next generation of sequencers will determine the capacity to analyze amplicon diversity and to assign amplicons to species instead of genera. Additionally, the prominent advantage of pyrosequencing is its automation, which leads to increased standardization, from DNA extraction to sequencing analysis, allowing multicenter studies to be carried out at without compromising reproducibility.

Conclusion

The aim of microbiological diagnosis from CF patients is to provide data with which clinicians can make rational and effective therapeutic decisions. Given the currently acknowledged polymicrobial nature of CF sputa [1], [8], [9], [17], [29], [31], better knowledge of sputum microbiota would represent a major advance in our understanding of the disease. In light of this concept of CF lung microbiota [61], [78], [96], high-throughput sequencing, due to its potential for massive direct sequencing after a single run of DNA amplification and automation, appears to be the most promising approach. The present study should stimulate a debate over the best way to set up new studies with the aim of combing (i) new technology (deep-sequencing), (ii) ecological tools (to analyze dynamics, diversity and relative species abundance, as species distribution is ecologically important in terms of community interactions [31], [78]), and (iii) clinically relevant information (e.g. pulmonary exacerbation in which SMG bacteria have been implicated when chronic colonization by P. aeruginosa develops a loss of virulence [1]) as well as the impact of therapeutics (long-term antibiotics cause a decline in bacterial diversity and inadvertently allow P. aeruginosa to flourish [58]; little is known about the impact of azole on fungal biota in CF).

Clearly, further metagenomic research, for which a scientific framework is needed as are well-designed translational studies, is now warranted to enhance knowledge of the process that drives the progression of CF respiratory disease. A comprehensive view of bacterial plus fungal microbiota present in CF lower airways has the potential to dramatically improve survival in CF patients. Moreover, it will pave the way for developing personalized drug therapy strategies based on the manipulation of complex microflora (i.e. controlling growth of less desirable microorganisms or controlling biofilm-associated infections as recently proposed [59], [89]).

Materials and Methods

Sample collection and DNA extraction

Patients were eligible if they could be classified as clinically stable (i.e., being followed-up during their annual check-up without exacerbation status). All volunteers with CF were required to have a well-documented diagnosis, with either the two mutations identified in the CFTR gene or an abnormally high sweat chloride test (Table 1). The four CF individuals selected for the study consisted of two males and two females, with an age range of 19 to 39 years. All clinical, therapeutic, radiological, and biological data were collected by clinical staff at the time of the visit (Tables 1 and 2). Human sputum samples (two samples collected for each patient at two visits) were collected by expectoration into a sterile cup after a water rinse to prevent excessive salivary contamination, [17], [23], [31]. Sputa were homogenized for 30 min at 37°C with Digest-EUR® (Eurobio, France) in 1∶1 (v∶v) ratio (final volume of approximately 10 ml), and mycological cultures were performed after direct examination, as previously described [4]. Briefly, 20 µl aliquots of the digested sample were inoculated onto three growth media: CandiSelectTM4 (Bio-Rad; incubation at 37°C for 3 weeks), Sabouraud glucose peptone agar with 0.5 g/L amikacin (incubation at 25°C for 3 weeks), and 1∶2 diluted Sabouraud glucose agar with 0.5 g/L amikacin (incubation at 25°C for 3 weeks). All sputa were inoculated in parallel onto five agar plates including chocolate Poly ViteX agar, Columbia colistin-nalidixic acid agar, Bromo Cresol Purple agar, blood agar (all purchased from bioMérieux, France) and incubated at 37°C for 48 h) and Cepacia agar (purchased from AES Laboratory, France), and incubated at 30°C for 5 days). Colonies growing on these media were identified using conventional methods or spectrometry. Then, 200 µL of each digested sample were frozen at −20°C until use. Samples were first ground in liquid nitrogen with a mortar and pestle. DNA was then extracted using the High Pure PCR Template Preparation kit (Roche Applied Science, Germany) according to manufacturer's protocol, except for the proteinase K digestion step, which was performed for 1 h at 70°C rather than 10 min. Total DNA concentrations ranged from 50 to 75 ng/µL, using NanoDrop® ND-1000 spectrophotometer. A nested PCR targeting Pneumocystis jirovecii, a difficult-to-culture micromycete, and a real-time PCR targeting Aspergillus fumigatus, were retrospectively done as described previously [26], [27]. No significant PCR inhibitions were observed when DNA samples were diluted in 1/10.

Ethics Statement

Sputa from four CF patients who volunteered for the study were collected at the Lille Adult CF center, in accordance with the ethical guidelines of Lille University Hospital. This study was part of the “MucoFong” protocol and was approved by the Institutional Human Care and Use Committee of the Lille University Hospital (Comité de Protection des Personnes Nord Ouest IV - reference Number CPP 06/84; assurance number: SHAM 127795). Written informed consents were provided by study participants.

Pyrosequencing analysis

Two sets of primers were used to amplify the 16S rDNA and ITS2 loci from prokaryotes and fungi, respectively. The first set of primers, 3271-16S-F (TACGGRAGGCAGCAG) and 3271-16S-R (GGACTACCAGGGTATCTAAT), was designed to amplify a 465 bp region containing the complete V3 domain of all prokaryotic 16S rDNA genes [97]. The second set, composed of primers 3271-ITS2F (CARCAAYGGATCTCTTGG) and 3271-ITS2R (GATATGCTTAAGTTCAGCGGGT) was designed to amplify a 340–360 bp fragment of the ITS2 region from all major phyla of fungi, according to the use for reconstructing phylogenies at a higher taxonomical level of this region [98]. A 10 bp tag specific to each of the eight samples, a 4 bp TCAG key, and a 21 bp adapter for the GS FLX system, were added to the sequences of both primers sets. PCRs were carried out using standard conditions for Taq DNA polymerase with 10 ng of DNA as template. After the denaturation step at 95°C for 5 min, 35 cycles of amplification were performed with a GeneAmp PCR System cycler (Applied Biosystems) as follows: 30 s at 95°C, 30 s at 50°C and 1 min at 72°C. Each DNA sample was analyzed in duplicate. The Genoscreen company (Pasteur Institute of Lille, France) carried out the pyrosequencing. The library and the 454 GS FLX Titanium (Roche) pyrosequencing runs were prepared according to manufacturer's recommendations. We obtained 326,277 and 133,317 sequences with the first (16S prokaryotes) and second (ITS2 fungi) set of primers, respectively. The sequences or reads were classified according to the presence of the tag corresponding to each of the eight samples of interest. Primers, tag and key fragments were not included in sequence analysis.

For identification, the 16S rDNA gene sequences were compared to the Silva SSU rRNA database (http://www.arb-silva.de/) release 102 (updated on February 15, 2010) comprising 1,246,462 SSU rRNA sequences using BLASTN software [99]. For ITS2 sequence identification, we constructed a fungal ITS2 database, based on the following steps: (i) a search through the complete nucleotide database of GenBank for potential ITS2 sequences, (ii) selection of ITS2 sequences that included the sequences of the primers designed in the present study, and (iii) inclusion of human genome sequences that were 500 bp long with at least one of the two primers to filter sequences belonging to host human cells (indicated as “Homo sapiens” in the final taxonomic assignment of the pyrosequencing ITS2 reads). This ITS2 database, named ITS2dbScreen, is available on request via the web site of the Genoscreen company (www.genoscreen.fr).

BLAST results (with a 10−5 E-value threshold) were visualized using the metagenomic software MEGAN [100]. Based on NCBI taxonomy, this software explores the taxonomic content of the samples with the option “import BLASTN”. The program uses several thresholds to generate sequence-taxon matches. The “min-score” filter, corresponding to a bit score cutoff value, was set at 35 for 16S rDNA amplicons as previously described [19], and at 200 for ITS2 amplicons to obtain an alignment with a minimum of 100 nucleotides. The “top-percent” filter used to select hits whose scores lay within a given percentage of the highest bit score, was set at 10 and at 5 for 16S rDNA and the ITS2 loci, respectively. The “min-support core” filter, used to set a threshold for the minimum number of sequences that must be assigned to a taxon, was set to 5. These stringent parameters should result in a “conservative” assignment of many sequences to internal branches (i.e. with less precision) of the taxonomic tree. Distribution of the sequences was schematically represented by Neighbor-Joining (NJ) tree diagrams (Figures S1, S2, S3, S4).

Rarefaction curves and richness estimator

The quality of the sampling effort was assessed through the calculation of rarefaction curves, i.e. the number of operational taxonomic units (OTUs) with respect to the number of reads [101]. The prokaryote 16S rDNA and fungus ITS2 pyrosequences were sorted by decreasing length and clustered with 97% similarity using Uclust (http://www.drive5.com/usearch/) [102]. Rarefaction curves were calculated according to Uclust results using a perl script. The total richness of a community from a patient-sample was estimated using the Chao1 richness estimator [103]. This non-parametric estimator allows cross-sample comparison of species diversity. The Chao1 index was calculated from Uclust results using the formula: Chao1 = n+(n1*(n1−1))/(2*(n2+1)), where n is total number of OTUs, n1, the number of OTUs composed of one read, and n2, the number of OTUs composed of two reads. These diversity indices and richness estimators were then used to compare the relative complexities of communities and to estimate the completeness of sampling.

Phylogenetic analysis

The phylogenetic trees inferred from 16S rDNA and ITS2 pyrosequences were used to compare biodiversity of specific taxa or within genera between samples for the same patient and/or between patients. The bacterial 16S rDNA sequences corresponding to the genera Streptococcus, Haemophilus and Pseudomonas, and the fungal ITS2 rDNA reads corresponding to those of Aspergillus, Candida and Malassezia were extracted from the pyrosequencing database using MEGAN, individually sorted by size, and clustered by homology (with a 97% identity threshold) using Uclust [102]. The longest read (>400 bp) from each cluster was selected as the representative sequence and submitted to a BLAST search [99] on the non-redundant nucleotide database (NCBI) to determine an approximate phylogenetic affiliation. The representative sequences and reference sequences were aligned using Muscle [102] as implemented in the SeaView4 program [104]. The resulting alignments were manually curated with the BioEdit software (http://www.mbio.ncsu.edu/bioedit/bioedit.html), and phylogenetic trees were constructed using both the NJ method from the SeaView4 package [104] and the Bayesian method implemented in MrBayes3 software (http://mrbayes.csit.fsu.edu/index.php) [105]. Since topologies of the phylogenetic trees generated by the two methods were similar, only the NJ trees are shown. The reliability of internal branches was assessed using the bootstrap method implemented in SeaView4 with 1000 replicates; only probabilities of more than 50% are shown at the tree nodes. Phylogenetic trees were edited using Dendroscope [106].

The pyrosequences were deposited in GenBank-SRA under the accession number SRA049426.2.

Statistical analysis

Numerical variables were described as medians and interquartile ranges (Q1, Q3). To study the relationship between clinical data, taxon richness, and community composition of sputum samples, Spearman's correlation coefficient were calculated. P-values of less than 0.05 were considered as significant. All statistical analyses were performed using SAS software (SAS Institute, Cary, NC, USA; version 9.2).

Supporting Information

Figure S1.

Taxonomic assignment of the 16S rDNA (A) and ITS2 (B) reads obtained from Patient 1.

https://doi.org/10.1371/journal.pone.0036313.s001

(TIF)

Figure S2.

Taxonomic assignment of the 16S rDNA (a) and ITS2 (b) reads obtained from Patient 2.

https://doi.org/10.1371/journal.pone.0036313.s002

(TIF)

Figure S3.

Taxonomic assignment of the 16S rDNA (a) and ITS2 (b) reads obtained from Patient 3.

https://doi.org/10.1371/journal.pone.0036313.s003

(TIF)

Figure S4.

Taxonomic assignment of the 16S rDNA (a) and ITS2 (b) reads obtained from Patient 4. Footnotes for figures S1 to S4. Reads obtained from sputum samples of Patients 1–4 were analyzed using the software MEGAN, after BLASTN search against databases (see Material and Methods section). The MEGAN software plots on schematic trees represent the number of pyrosequence reads matching a particular taxonomical group. The tree displays all taxonomic groups identified from the assignment of reads obtained either with prokaryotic primers (Figures S1AS4A), or fungus-designed primers (Figures S1BS4B).

https://doi.org/10.1371/journal.pone.0036313.s004

(TIF)

Figure S5.

NJ-phylogenetic tree of ITS2 sequences from the genus Aspergillus .

https://doi.org/10.1371/journal.pone.0036313.s005

(TIF)

Figure S6.

NJ-phylogenetic tree of 16S rRNA sequences from the genus Pseudomonas .

https://doi.org/10.1371/journal.pone.0036313.s006

(TIF)

Figure S7.

NJ-phylogenetic tree of 16S rRNA sequences from the genus Streptococcus .

https://doi.org/10.1371/journal.pone.0036313.s007

(TIF)

Figure S8.

NJ-phylogenetic tree of ITS2 sequences from Candida albicans .

https://doi.org/10.1371/journal.pone.0036313.s008

(TIF)

Figure S9.

NJ-phylogenetic tree of ITS2 sequences from the Candida parapsilosis complex. Footnotes for figures S5 to S9. Neighbor-joining trees of the ITS2 or 16SrRNA sequences from the genus Aspergillus (Figure S5), Pseudomonas (Figure S6), and Streptococcus (Figure S7), and the species C. albicans (Figure S8) and the C. parapsilosis complex (Figure S9). The representative sequences corresponding to Patient 1 in blue, Patient 2 in green, Patient 3 in red and Patient 4 in yellow, while dark and light colour intensity were corresponding to the first and second sampling dates, respectively. Numbers in brackets indicate the number of reads composing each cluster. Clusters composed of reads that are at least 50% greater than the number of reads composing the most dominant cluster are in bold. Bootstrap values (threshold >50) are indicated at the nodes.

https://doi.org/10.1371/journal.pone.0036313.s009

(TIF)

Acknowledgments

The authors would like to thank Carolyn Engel-Gautier for English editing. This study was presented as late-breaker abstract for a poster presentation at the 51st ICAAC congress (M-1523a).

Author Contributions

Conceived and designed the experiments: LD ED-C TS-N EV. Performed the experiments: LD SM EF CH FW MC. Wrote the paper: LD SM EV MC. Performed the molecular analysis: LD MC SM. Collaborated on the molecular analysis: SL AP BW. Physicians in charge of the CF patients: SL AP BW. Involved in the statistical analysis: JS. In charge of the ITS2 database constitution: CH.

References

  1. 1. Sibley CD, Parkins MD, Rabin HR, Duan K, Norgaard JC, et al. (2008) A polymicrobial perspective of pulmonary infections exposes an enigmatic pathogen in cystic fibrosis patient. Proc Natl Acad Sci USA 105: 15070–15075.
  2. 2. Nagano Y, Millar BC, Goldsmith CE, Walker JM, Elborn JS, et al. (2008) Development of selective media for the isolation of yeasts and filamentous fungi from the sputum of adult patients with cystic fibrosis (CF). J Cyst Fibros 7: 566–572.
  3. 3. Pihet M, Carrère J, Cimon B, Chabasse D, Delhaes L, et al. (2009) Occurrence and relevance of filamentous fungi in respiratory secretions of patients with cystic fibrosis - a review. Med Mycol 47: 387–397.
  4. 4. Borman AM, Palmer MD, Delhaes L, Carrère J, Favennec L, et al. (2010) Lack of standardization in the procedures for mycological examination of sputum samples from CF patients: a possible cause for variations in the prevalence of filamentous fungi. Med Mycol 48: Suppl 1S88–97.
  5. 5. Sudfeld CR, Dasenbrook EC, Merz WG, Carroll KC, Boyle MP (2010) Prevalence and risk factors for recovery of filamentous fungi in individuals with cystic fibrosis. J Cyst Fibros 9: 110–116.
  6. 6. Horré R, Symoens F, Delhaes L, Bouchara JP (2010) Fungal respiratory infections in cystic fibrosis: a growing problem. Med Mycol 48: Suppl 1S1–3.
  7. 7. Hauser AR, Jain M, Bar-Meir M, McColley SA (2011) Clinical significance of microbial infection and adaptation in cystic fibrosis. Clin Microbiol Rev 24: 29–70.
  8. 8. Bittar F, Richet H, Dubus JC, Reynaud-Gaubert M, Stremler N, et al. (2008) Molecular detection of multiple emerging pathogens in sputa from cystic fibrosis patients. PLoS One 3: e2908.
  9. 9. Bouchara JP, Hsieh HY, Croquefer S, Barton R, Marchais V, et al. (2009) Development of an oligonucleotide array for direct detection of fungi in sputum samples from patients with cystic fibrosis. J Clin Microbiol 47: 142–152.
  10. 10. Chabé M, Aliouat-Denis CM, Delhaes L, Aliouat el M, Viscogliosi E, et al. (2011) Pneumocystis: from a doubtful unique entity to a group of highly diversified fungal species. FEMS Yeast Res 11: 2–17.
  11. 11. Colombo AP, Boches SK, Cotton SL, Goodson JM, Kent R, et al. (2009) Comparisons of subgingival microbial profiles of refractory periodontitis, severe periodontitis, and periodontal health using the human oral microbe identification microarray. J Periodontol 80: 1421–1432.
  12. 12. Grice EA, Kong HH, Conlan S, Deming CB, Davis J, et al. (2009) Topographical and temporal diversity of the human skin microbiome. Science 324: 1190–1192.
  13. 13. Wolcott RD, Gontcharova V, Sun Y, Dowd SE (2009) Evaluation of the bacterial diversity among and within individual venous leg ulcers using bacterial tag-encoded FLX and Titanium amplicon pyrosequencing and metagenomic approaches. BMC Microbiol 9: 226.
  14. 14. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, et al. (2010) A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464: 59–65.
  15. 15. Fava F, Danese S (2011) Intestinal microbiota in inflammatory bowel disease: Friend of foe? World J Gastroenterol 17: 557–566.
  16. 16. Armougom F, Bittar F, Stremler N, Rolain JM, Robert C, et al. (2009) Microbial diversity in the sputum of a cystic fibrosis patient studied with 16S rDNA pyrosequencing Eur J Clin Microbiol Infect Dis 28: 1151–1154.
  17. 17. Guss AM, Roeselers G, Newton IL, Young CR, Klepac-Ceraj V, et al. (2011) Phylogenetic and metabolic diversity of bacteria associated with cystic fibrosis. ISME J 5: 20–29.
  18. 18. Willner D, Furlan M, Haynes M, Schmieder R, Angly FE, et al. (2009) Metagenomic analysis of respiratory tract DNA viral communities in cystic fibrosis and non-cystic fibrosis individuals. PLoS One 4: e7370.
  19. 19. Monchy S, Sanciu G, Jobard M, Rasconi S, Gerphagnon M, et al. (2011) Exploring and quantifying fungal diversity in freshwater lake ecosystems using rDNA cloning/sequencing and SSU tag pyrosequencing. Environ Microbiol 13: 1433–1453.
  20. 20. Harris JK, De Groote MA, Sagel SD, Zemanick ET, Kapsner R, et al. (2007) Molecular identification of bacteria in bronchoalveolar lavage fluid from children with cystic fibrosis. Proc Natl Acad Sci U S A 104: 20529–20533.
  21. 21. Friaza V, la Horra C, Rodríguez-Domínguez MJ, Martín-Juan J, Cantón R, et al. (2010) Metagenomic analysis of bronchoalveolar lavage samples from patients with idiopathic interstitial pneumonia and its antagonic relation with Pneumocystis jirovecii colonization. J Microbiol Methods 82: 98–101.
  22. 22. Leclair LW, Hogan DA (2010) Mixed bacterial-fungal infections in the CF respiratory tract. Med Mycol 48: Suppl 1S125–32.
  23. 23. van der Gast CJ, Walker AW, Stressmann FA, Rogers GB, Scott P, et al. (2011) Partitioning core and satellite taxa from within cystic fibrosis lung bacterial communities. ISME J 5: 780–791.
  24. 24. Willner D, Haynes MR, Furlan M, Schmieder R, Lim YW, et al. (2011) Spatial distribution of microbial communities in the cystic fibrosis lung. ISME J.
  25. 25. Nagano Y, Elborn JS, Millar BC, Walker JM, Goldsmith CE, et al. (2010) Comparison of techniques to examine the diversity of fungi in adult patients with cystic fibrosis. Med Mycol 48: 166–176.
  26. 26. Montes-Cano MA, de la Horra C, Dapena FJ, Mateos I, Friaza V, et al. (2007) Dynamic colonisation by different Pneumocystis jirovecii genotypes in cystic fibrosis patients. Clin Microbiol Infect 13: 1008–1011.
  27. 27. Fréalle E, Decrucq K, Botterel F, Bouchindhomme B, Camus D, et al. (2009) Diagnosis of invasive aspergillosis using bronchoalveolar lavage in haematology patients: influence of bronchoalveolar lavage human DNA content on real-time PCR performance. Eur J Clin Microbiol Infect Dis 28: 223–232.
  28. 28. Sundquist A, Bigdeli S, Jalili R, Druzin ML, Waller S, et al. (2007) Bacterial flora-typing with targeted, chip-based Pyrosequencing. BMC Microbiol 7: 108.
  29. 29. Rogers GB, Carroll MP, Serisier DJ, Hockey PM, Jones G, et al. (2004) Characterization of bacterial community diversity in cystic fibrosis lung infections by use of 16s ribosomal DNA terminal restriction fragment length polymorphism profiling. J Clin Microbiol 42: 5176–5183.
  30. 30. Tunney MM, Field TR, Moriarty TF, Patrick S, Doering G, et al. (2008) Detection of anaerobic bacteria in high numbers in sputum from patients with cystic fibrosis. Am J Respir Crit Care Med 177: 995–1001.
  31. 31. Rogers GB, Skelton S, Serisier DJ, van der Gast CJ, Bruce KD (2010) Determining cystic fibrosis-affected lung microbiology: comparison of spontaneous and serially induced sputum samples by use of terminal restriction fragment length polymorphism profiling. J Clin Microbiol 48: 78–86.
  32. 32. Samson RA, van der Lustgraaf B (1978) Aspergillus penicilloides and Eurotium halophilicum in association with house-dust mites. Mycopathologia 64: 13–16.
  33. 33. Harries MG, Lacey J, Tee RD, Cayley GR, Taylor AJ (1985) Didymella exitialis and late summer asthma. Lancet 1: 1063–1066.
  34. 34. Cimerman M, Gunde-Cimerman N, Zalar P, Perkovic T (1999) Femur osteomyelitis due to a mixed fungal infection in a previously healthy man. J Clin Microbiol 37: 1532–1535.
  35. 35. Guarro J, Gené J, Stchigel AM (1999) Developments in fungal taxonomy. Clin Microbiol Rev 12: 454–500.
  36. 36. Lachance MA, Daniel HM, Meyer W, Prasad GS, Gautam SP, et al. (2003) The D1/D2 domain of the large-subunit rDNA of the yeast species Clavispora lusitaniae is unusually polymorphic. FEMS Yeast Res 4: 253–258.
  37. 37. Bearchell SJ, Fraaije BA, Shaw MW, Fitt BD (2005) Wheat archive links long-term fungal pathogen population dynamics to air pollution. Proc Natl Acad Sci U S A 102: 5438–5442.
  38. 38. Casagrande BF, Flückiger S, Linder MT, Johansson C, Scheynius A, et al. (2006) Sensitization to the yeast Malassezia sympodialis is specific for extrinsic and intrinsic atopic eczema. J Invest Dermatol 126: 2414–2421.
  39. 39. Merget R, Sander I, Rozynek P, Heinze E, Imoehl M, et al. (2008) Occupational immunoglobulin E-mediated asthma due to Penicillium camemberti in a dry-sausage packer. Respiration 76: 109–111.
  40. 40. Muthig M, Hebestreit A, Ziegler U, Seidler M, Müller FM (2010) Persistence of Candida species in the respiratory tract of cystic fibrosis patients. Med Mycol 48: 56–63.
  41. 41. Gomez-Lopez A, Pan D, Cuesta I, Alastruey-Izquierdo A, Rodriguez-Tudela JL, et al. (2010) Molecular identification and susceptibility profile in vitro of the emerging pathogen Candida kefyr. Diagn Microbiol Infect Dis 66: 116–119.
  42. 42. Gammelsrud KW, Sandven P, Høiby EA, Sandvik L, Brandtzaeg P, et al. (2011) Colonization by Candida in children with cancer, children with cystic fibrosis, and healthy controls. Clin Microbiol Infect. 10.1111/j.1469-0691.2011.03528.x.
  43. 43. Reihill JA, Moore JE, Elborn JS, Ennis M (2011) Effect of Aspergillus fumigatus and Candida albicans on pro-inflammatory response in cystic fibrosis epithelium. J Cyst Fibros 10: 401–406.
  44. 44. Miceli MH, Diaz JA, Lee SA (2011) Emerging opportunistic yeast infections. Lancet Infect Dis 11: 142–151.
  45. 45. Mortensen KL, Johansen HK, Fuursted K, Knudsen JD, Gahrn-Hansen B, et al. (2011) A prospective survey of Aspergillus spp. in respiratory tract samples: prevalence, clinical impact and antifungal susceptibility. Eur J Clin Microbiol Infect Dis 30: 1355–1363.
  46. 46. Balajee SA, Gribskov JL, Hanley E, Nickle D, Marr KA (2005) Aspergillus lentulus sp. nov., a new sibling species of A. fumigatus. Eukaryot Cell 4: 625–632.
  47. 47. Balajee SA, Nickle D, Varga J, Marr KA (2006) Molecular studies reveal frequent misidentification of Aspergillus fumigatus by morphotyping. Eukaryot Cell 5: 1705–1712.
  48. 48. A-Nahimana A, Francioli P, Blanc DS, Bille J, Wakefield AE, et al. (2000) Determination of the copy number of the nuclear rDNA and beta-tubulin genes of Pneumocystis carinii f. sp. hominis using PCR multicompetitors. J Eukaryot Microbiol 47: 368–372.
  49. 49. Chotirmall SH, O'Donoghue E, Bennett K, Gunaratnam C, O'Neill SJ, et al. (2010) Sputum Candida albicans presages FEV1‚ decline and hospital-treated exacerbations in cystic fibrosis. Chest 138: 1186–1195.
  50. 50. Peterson SW (2008) Phylogenetic analysis of Aspergillus species using DNA sequences from four loci. Mycologia 100: 205–226.
  51. 51. Moore JE, Xu J, Millar BC, Elshibly S (2002) Edible dates (Phoenix dactylifera), a potential source of Cladosporium cladosporioides and Sporobolomyces roseus: implications for public health. Mycopathologia 154: 25–28.
  52. 52. Vilela DM, Pereira GV, Silva CF, Batista LR, Schwan RF (2010) Molecular ecology and polyphasic characterization of the microbiota associated with semi-dry processed coffee (Coffea arabica L.). Food Microbiol 27: 1128–1135.
  53. 53. Hirooka Y, Rossman AY, Chaverri P (2011) A morphological and phylogenetic revision of the Nectria cinnabarina species complex. Stud Mycol 68: 35–56.
  54. 54. Inácio J, Portugal L, Spencer-Martins I, Fonseca A (2005) Phylloplane yeasts from Portugal: seven novel anamorphic species in the Tremellales lineage of the Hymenomycetes (Basidiomycota) producing orange-coloured colonies. FEMS Yeast Res 5(12): 1167–1183.
  55. 55. Sayers EW, Barrett T, Benson DA, Bryant SH, Canese K, et al. (2009) Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 37(Database issue): D5–15.
  56. 56. European Commission Health & Consumers Directorate-General Phlebiopsis gigantea. SANCO/163/08-rev.3 http://ec.europa.eu/food/plant/protection/evaluation/existactive/phlebiopsis_gigantea.pdf.
  57. 57. Fegan M, Francis P, Hayward AC, Fuerst JA (1991) Heterogeneity, persistence, and distribution of Pseudomonas aeruginosa genotypes in cystic fibrosis patients. Clin Microbiol 29: 2151–2157.
  58. 58. Klepac-Ceraj V, Lemon KP, Martin TR, Allgaier M, Kembel SW, et al. (2010) Relationship between cystic fibrosis respiratory tract bacterial communities and age, genotype, antibiotics and Pseudomonas aeruginosa. Environ Microbiol 12: 1293–1303.
  59. 59. Jenkinson HF, Lamont RJ (2005) Oral microbial communities in sickness and in health. Trends Microbiol 13: 589–595.
  60. 60. Hogan DA Talking to themselves: autoregulation and quorum sensing in fungi. (2006) Eukaryot Cell 5: 613–619.
  61. 61. Rogers GB, Carroll MP, Hoffman LR, Walker AW, Fine DA, et al. (2010) Comparing the microbiota of the cystic fibrosis lung and human gut. Gut Microbes 1: 85–93.
  62. 62. Perkins SD, Woeltje KF, Angenent LT (2010) Endotracheal tube biofilm inoculation of oral flora and subsequent colonization of opportunistic pathogens. Int J Med Microbiol 300: 503–511.
  63. 63. Grinwis ME, Sibley CD, Parkins MD, Eshaghurshan CS, Rabin HR, et al. (2009) Characterization of Streptococcus milleri group isolates from expectorated sputum of adult patients with cystic fibrosis. J Clin Microbiol 48: 395–401.
  64. 64. Vanhee LM, Symoens F, Bouchara JP, Nelis HJ, Coenye T (2008) High-resolution genotyping of Aspergillus fumigatus isolates recovered from chronically colonised patients with cystic fibrosis. Eur J Clin Microbiol Infect Dis 27: 1005–1007.
  65. 65. Gammelsrud KW, Sandven P, Høiby EA, Sandvik L, Brandtzaeg P, et al. (2011) Colonization by Candida in children with cancer, children with cystic fibrosis, and healthy controls. Clin Microbiol Infect 17: 1875–1881.
  66. 66. Neuvéglise C, Sarfati J, Debeaupuis JP, Vu Thien H, Just J, et al. (1997) Longitudinal study of Aspergillus fumigatus strains isolated from cystic fibrosis patients. Eur J Clin Microbiol Infect Dis 16: 747–750.
  67. 67. Cimon B, Symoens F, Zouhair R, Chabasse D, Nolard N, et al. (2001) Molecular epidemiology of airway colonisation by Aspergillus fumigatus in cystic fibrosis patients. J Med Microbiol 50: 367–374.
  68. 68. Tavanti A, Davidson AD, Gow NA, Maiden MC, Odds FC (2005) Candida orthopsilosis and Candida metapsilosis spp. nov. to replace Candida parapsilosis groups II and III. J Clin Microbiol 43: 284–292.
  69. 69. Gilgado F, Cano J, Gené J, Sutton DA, Guarro J (2008) Molecular and phenotypic data supporting distinct species statuses for Scedosporium apiospermum and Pseudallescheria boydii and the proposed new species Scedosporium dehoogii. J Clin Microbiol 46: 766–771.
  70. 70. Delhaes L, Harun A, Chen SC, Nguyen Q, Slavin M, et al. (2008) Molecular typing of Australian Scedosporium isolates showing genetic variability and numerous S. aurantiacum. Emerg Infect Dis 14: 282–290.
  71. 71. Lockhart SR, Messer SA, Pfaller MA, Diekema DJ (2008) Geographic distribution and antifungal susceptibility of the newly described species Candida orthopsilosis and Candida metapsilosis in comparison to the closely related species Candida parapsilosis. J Clin Microbiol 46: 2659–2664.
  72. 72. Hogan DA, Vik A, Kolter R (2004) A Pseudomonas aeruginosa quorum-sensing molecule influences Candida albicans morphology. Mol Microbiol 54: 1212–1223.
  73. 73. McAlester G, O'Gara F, Morrissey JP (2008) Signal-mediated interactions between Pseudomonas aeruginosa and Candida albicans. J Med Microbiol 57: 563–569.
  74. 74. Reen FJ, Mooij MJ, Holcombe LJ, McSweeney CM, McGlacken GP, et al. (2011) The Pseudomonas quinolone signal (PQS), and its precursor HHQ, modulate interspecies and interkingdom behaviour. FEMS Microbiol Ecol 77: 413–428.
  75. 75. Hall RA, Turner KJ, Chaloupka J, Cottier F, De Sordi L, et al. (2011) The quorum sensing molecules Farnesol/Homoserine lactone and Dodecanol operate via distinct modes of action in Candida albicans. Eukaryot Cell 10: 1034–1042.
  76. 76. Erb-Downward JR, Thompson DL, Han MK, Freeman CM, McCloskey L, et al. (2011) Analysis of the lung microbiome in the “healthy” smoker and in COPD. PLoS One 6: e16384.
  77. 77. Stecher B, Hardt WD (2008) The role of microbiota in infectious disease. Trends Microbiol 16: 107–114.
  78. 78. Sibley CD, Surette MG (2011) The polymicrobial nature of airway infections in cystic fibrosis: Cangene Gold Medal Lecture. Can J Microbiol 57: 69–77.
  79. 79. Sibley CD, Duan K, Fischer C, Parkins MD, Storey DG, et al. (2008) Discerning the complexity of community interactions using a Drosophila model of polymicrobial infections. PLoS Pathog 4: e1000184.
  80. 80. Amin R, Dupuis A, Aaron SD, Ratjen F (2011) The effect of chronic infection with Aspergillus fumigatus on lung function and hospitalization in patients with cystic fibrosis. Chest 137: 171–176.
  81. 81. European Cystic Fibrosis Society website (2012) Available: http://www.ecfs.eu/special-projects/mucofong-international-project, Updated 10th of February 2012.
  82. 82. Yoon SS, Hennigan RF, Hilliard GM, Ochsner UA, Parvatiyar K, et al. (2002) Pseudomonas aeruginosa anaerobic respiration in biofilms: relationships to cystic fibrosis pathogenesis. Dev Cell 3: 593–603.
  83. 83. Alvarez-Ortega C, Harwood CS (2007) Responses of Pseudomonas aeruginosa to low oxygen indicate that growth in the cystic fibrosis lung is by aerobic respiration. Mol Microbiol 65: 153–165.
  84. 84. Dumitru R, Hornby JM, Nickerson KW (2004) Defined anaerobic growth medium for studying Candida albicans basic biology and resistance to eight antifungal drugs. Antimicrob Agents Chemother 48: 2350–2354.
  85. 85. Dumitru R, Navarathna DH, Semighini CP, Elowsky CG, Dumitru RV, et al. (2007) In vivo and in vitro anaerobic mating in Candida albicans. Eukaryot Cell 6: 465–472.
  86. 86. Seidler MJ, Salvenmoser S, Müller FM (2008) Aspergillus fumigatus forms biofilms with reduced antifungal drug susceptibility on bronchial epithelial cells. Antimicrob Agents Chemother 52: 4130–4136.
  87. 87. Lattif AA, Mukherjee PK, Chandra J, Swindell K, Lockhart SR, et al. (2010) Characterization of biofilms formed by Candida parapsilosis, C. metapsilosis, and C. orthopsilosis. Int J Med Microbiol 300: 265–270.
  88. 88. Bjarnsholt T, Jensen PØ, Jakobsen TH, Phipps R, Nielsen AK, et al. (2010) Quorum sensing and virulence of Pseudomonas aeruginosa during lung infection of cystic fibrosis patients. PLoS One 5: e10115.
  89. 89. Rybtke MT, Jensen PØ, Høiby N, Givskov M, Tolker-Nielsen T, et al. (2011) The implication of Pseudomonas aeruginosa biofilms in infections. Inflamm Allergy Drug Targets 10: 141–157.
  90. 90. Mowat E, Rajendran R, Williams C, McCulloch E, Jones B, et al. (2010) Pseudomonas aeruginosa and their small diffusible extracellular molecules inhibit Aspergillus fumigatus biofilm formation. FEMS Microbiol Lett 313: 96–102.
  91. 91. Bandara HM, Yau JY, Watt RM, Jin LJ, Samaranayake LP (2010) Pseudomonas aeruginosa inhibits in-vitro Candida biofilm development. BMC Microbiol 10: 125.
  92. 92. Semighini CP, Hornby JM, Dumitru R, Nickerson KW, Harris SD (2006) Farnesol-induced apoptosis in Aspergillus nidulans reveals a possible mechanism for antagonistic interactions between fungi. Mol Microbiol 59: 753–764.
  93. 93. Tré-Hardy M, Nagant C, El Manssouri N, Vanderbist F, Traore H, et al. (2010) Efficacy of the combination of tobramycin and a macrolide in an in vitro Pseudomonas aeruginosa mature biofilm model. Antimicrob Agents Chemother 54: 4409–4415.
  94. 94. Fredricks DN, Smith C, Meier A (2005) Comparison of six DNA extraction methods for recovery of fungal DNA as assessed by quantitative PCR. J Clin Microbiol 43: 5122–5128.
  95. 95. Feinstein LM, Sul WJ, Blackwood CB (2009) Assessment of bias associated with incomplete extraction of microbial DNA from soil. Appl Environ Microbiol 75: 5428–5433.
  96. 96. Zemanick ET, Sagel SD, Harris JK (2011) The airway microbiome in cystic fibrosis and implications for treatment. Curr Opin Pediatr 23: 319–324.
  97. 97. Wang Y, Qian PY (2009) Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies. PLoS one 4: e7401.
  98. 98. Agüero-Chapin G, Sanchez-Rodriguez A, Hidalgo-Yanes PI, Pérez-Castillo Y, Molina-Ruiz R, et al. (2011) An alignment-free approach for eukaryotic ITS2 annotation and phylogenetic inference. PLoS One 6: e26638.
  99. 99. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403–410.
  100. 100. Huson DH, Auch AF, Qi J, Schuster SC (2007) MEGAN analysis of metagenomic data. Genome Res 17: 377–386.
  101. 101. Hughes J, Hellmann J, Ricketts T, Bohannan B (2001) Counting the uncountable: statistical approaches to estimating microbial diversity. Appl Environ Microbiol 67: 4399–4406.
  102. 102. Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26: 2460–2461.
  103. 103. Chao A (1984) Non-parametric estimation of the number of classes in a population. Scand J Stat 11: 265–270.
  104. 104. Gouy M, Guindon S, Gascuel O (2010) SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol 27: 221–224.
  105. 105. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.
  106. 106. Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, et al. (2007) Dendroscope- An interactive viewer for large phylogenetic trees. BMC Bioinformatics 8: 460.