Characterization and quantification of the fungal microbiome in serial samples from individuals with cystic fibrosis
© Willger et al.; licensee BioMed Central Ltd. 2014
Received: 28 April 2014
Accepted: 25 September 2014
Published: 3 November 2014
Human-associated microbial communities include fungi, but we understand little about which fungal species are present, their relative and absolute abundances, and how antimicrobial therapy impacts fungal communities. The disease cystic fibrosis (CF) often involves chronic airway colonization by bacteria and fungi, and these infections cause irreversible lung damage. Fungi are detected more frequently in CF sputum samples upon initiation of antimicrobial therapy, and several studies have implicated the detection of fungi in sputum with worse outcomes. Thus, a more complete understanding of fungi in CF is required.
We characterized the fungi and bacteria in expectorated sputa from six CF subjects. Samples were collected upon admission for systemic antibacterial therapy and upon the completion of treatment and analyzed using a pyrosequencing-based analysis of fungal internal transcribed spacer 1 (ITS1) and bacterial 16S rDNA sequences. A mixture of Candida species and Malassezia dominated the mycobiome in all samples (74%–99% of fungal reads). There was not a striking trend correlating fungal and bacterial richness, and richness showed a decline after antibiotic therapy particularly for the bacteria. The fungal communities within a sputum sample resembled other samples from that subject despite the aggressive antibacterial therapy. Quantitative PCR analysis of fungal 18S rDNA sequences to assess fungal burden showed variation in fungal density in sputum before and after antibacterial therapy but no consistent directional trend. Analysis of Candida ITS1 sequences amplified from sputum or pure culture-derived genomic DNA from individual Candida species found little (<0.5%) or no variation in ITS1 sequences within or between strains, thereby validating this locus for the purpose of Candida species identification. We also report the enhancement of the publically available Visualization and Analysis of Microbial Population Structures (VAMPS) tool for the analysis of fungal communities in clinical samples.
Fungi are present in CF respiratory sputum. In CF, the use of intravenous antibiotic therapy often does not profoundly impact bacterial community structure, and we observed a similar stability in fungal species composition. Further studies are required to predict the effects of antibacterials on fungal burden in CF and fungal community stability in non-CF populations.
Subjects with the genetic disease cystic fibrosis (CF) often have complex, chronic mixed bacterial-fungal lung infections that cause severe lung damage that can progress to respiratory failure. Molecular profiling studies strongly support a model in which CF lung infections contain multiple bacteria[1–3], and recent work has revealed that these communities are often surprisingly stable within subjects over time[2, 4–7]. Fungal species, such as Candida albicans and Aspergillus fumigatus, are also commonly detected in sputum cultures. Culture-based studies have detected C. albicans in ~40%–75% CF subjects[8–21] and described its association with reduced lung function[16, 22]. A large European study found that the presence of fungi by laboratory culture correlated with a 5%–10% lower lung function in CF. Two prospective studies have linked the detection of either C. albicans or A. fumigatus to lower lung function and increased frequency of disease exacerbations[16, 22]. Other fungi have also been detected in the CF lung[19, 23–27]. Challenges associated with the culturing of fungi, due to slow growth, specific medium requirements, and the lack of quantitative methods due in part to filamentous morphologies, have made it difficult to describe fungal species and their relative burdens in the CF lung. Consequently, despite data that suggest that the mycobiome poses a health threat to individuals with CF[16, 22], there is not yet a comprehensive standard of care upon the detection of fungi in sputum.
Individuals with CF-associated lung infections usually experience periods of stable lung function interspersed with episodes of disease exacerbation, which is loosely defined as a period of impaired lung function or the worsening of other health metrics that requires hospitalization for the administration of one or more intravenous antibiotics. The effects of antibiotic treatment on bacterial numbers and community structure are not fully understood and variable across subjects. One recent study found that antibiotic therapy does not cause reproducible shifts in bacterial community structure or bacterial density within sputum, while another study reported that antibiotic use was associated with decreased diversity and short-term shifts in community structure. The effects of antibiotic therapy on fungal loads or fungal species composition in the lung are not well known. Several studies have indicated that antibiotic usage precedes the increased detection of fungi in sputum[31–36]. Furthermore, antibiotic therapy in other systems correlates with an increased incidence of fungal infections, such as bloodstream infections[37–40], oral candidosis, and vaginal candidiasis, and increased fungal loads in the gastrointestinal tract[43–45]. A full characterization of the fungal species within sputum and other host-associated microbial communities will require the community-wide implementation of lysis methods that are effective for all bacterial and fungal species and the use of a combination of molecular methods that provide information on all of the species present. Most molecular analyses of host-associated fungi have relied upon amplification and sequencing of the internal transcribed spacer 1 (ITS1) region within the ribosomal RNA operon[46–49]. Due to a high degree of variation between even closely related species[50–52], this region can serve as a unique taxonomic identifier in both traditional capillary and next generation DNA sequencing studies[46–49].
In this report, we assessed the fungi within the CF lung in subjects hospitalized for treatment with antibacterial antibiotics and advanced the development of tools for the analysis of the fungal component of the microbiome (mycobiome) within clinical samples. Using primers that amplify the ITS1 region, we characterized the mycobiome and bacterial microbiome in expectorated sputum samples from six subjects before and after antibiotic therapy. Candida species dominated the mycobiomes, and their relative abundances were stable over time in contrast to less stable bacterial community structures. Using these samples, we also quantified the numbers of fungal genomes within each sample and found variability in how fungal densities in sputum changed over time and after treatment with antibacterial drugs. Lastly, towards facilitating future mycobiome studies, we assess intra-strain variability in Candida ITS1 sequences and describe improvement of the fungal ITS1 reference database associated with VAMPS (Visualization and Analysis of Microbial Population Structures) and its use in analyses of fungal ITS1 sequences amplified from host-associated microbial communities.
Adapting a publicly available tool for mycobiome analysis
To aid the analysis of fungal ITS1 sequences, particularly those amplified from clinical samples, we improved the ITS reference database for Global Alignment for Sequence Taxonomy (GAST) analyses within VAMPS, a widely used web-based tool that resolves phylogenetic affinity of marker gene sequences from diverse bacterial communities ([2, 55–57] as examples). The reference database for assigning the taxonomic affinity of fungal ITS1 amplicon sequences derives from the UNITE database. Our modifications included a number of corrections to the taxonomic information within records for medically relevant fungi, and the removal of some sequences that were exact matches to well-characterized plant sequences that had been incorrectly annotated as being of fungal origin. We also increased representation of the human-associated fungi Candida and Malassezia in the database by adding all Genbank database entries that contained the keywords “ITS”, “Candida”, and “Malassezia” and had complete annotations. This step increased the number of Candida species from 251 to 284 and the number of unique accession IDs containing the keyword “Candida” from 930 to 2021. For Malassezia, 72 unique accession IDs containing the keyword “Malassezia” from Genbank with complete genus and species information were added.
Overall assessment of the fungi identified in sputum from inpatient subjects
Using VAMPS, we analyzed the fungal community in CF sputum from six hospitalized subjects that received subject-specific antibacterial treatments (Additional file1: Table S1). Using primers that match conserved flanking sequences, we PCR-amplified the fungal ITS1 sequences from sputum DNA samples for pyrosequencing analyses. We determined the reproducibility of the fungal ITS1 amplification and deep sequencing protocol by analyzing technical replicate libraries created from the same DNA that had been extracted from a de-identified sputum sample collected for the purposes of the method development. The results from the two technical replicates were very similar in terms of the fraction of reads assigned to each fungus (<1% difference in percent abundance per sample) (Additional file2: Figure S1).
To validate the results obtained by GAST analyses, we employed BLAST to assign fungal ITS1 sequences to species and genus level categories (Additional file6: Table S5). For the BLAST searches, we batch blasted examples of all unique reads and assigned a genus or species to each set of sequences. Comparison of GAST and BLAST results initially highlighted the need for database improvement (described above), and once the database had been modified, the two methods yielded a similar number of taxa. Because GAST reports the consensus annotation for all of the sequences that match the query, the developmental state of the ITS1 database may limit the assignment of genus or species level information. In this case, sequences clustered by GAST in VAMPS can be easily retrieved for further analysis.
Because the characterization of fungi within host-associated microbial communities is at an early stage, we sought to determine if there were any ITS1 sequences that have not been previously described. In initial analyses, up to 13.3% of sequences in a sample could not be identified as fungi (BLAST scores of greater than e - 15 for the top hit). On average, samples had 3.6% of reads in this category. Further analysis querying Genbank’s Nucleotide collection (nr/nt) database revealed that most of these unassigned sequences corresponded to ITS1 sequences in the genomes of food-associated plants such as Pisum sativum (peas), Allium cepa (onions), Allium sativum (garlic), Solanum tuberosum (potato), and Glycine max (soybean) suggesting contamination of the sputum sample from subjects’ diets. The primers used to amplify fungal ITS1 sequences (Additional file7: Table S2) displayed high levels of similarity to comparable regions in plants. VAMPS allows one to limit analyses to specific taxonomic levels or to specific groups (for instance, fungi or Ascomycetes), and this eliminated the problem of contaminating sequences. A very small number of unassignable sequences (13 out of 300,072 reads) were also identified, but they were most often single reads that were only found within single samples and were not characterized further.
Comparison of intra- and inter-species ITS1 sequence differences in Candida spp
Some fungal genomes may contain up to 600 copies of the ribosomal RNA genes (rDNA) and intraspecific ITS variability can exceed 24.2%[59–62]. In contrast, C. albicans typically has 150–200 copies of the ITS1 locus, with one report estimating less than 0.2% ITS variability within a strain. Since Candida represented the predominant fungal genus in our CF subjects, we used Illumina MiSeq sequencing to investigate levels of ITS1 sequence variation within and among clinical C. albicans isolates. When the ITS1 region was amplified from genomic DNA from pure cultures of C. albicans wild type SC5314 and seven clinical isolates, then analyzed using Illumina MiSeq, the predominant sequence (>99% of reads) for each strain was the same (100% identical across all isolates). The second most abundant sequence accounted for only ~0.02% of reads. Similarly, we observed no evidence for intra-strain variation in the ITS1 sequences from the analysis of genomic DNAs isolated from pure cultures of three other Candida species (C. dubliniensis, C. parapsilosis, and C. kruseii) (data not shown).
When we compared the most abundant C. albicans ITS1 sequence in the sputum-derived datasets to the dominant sequence from the clinical isolates, we found that they were 100% identical (data not shown). Alignments of the three most common ITS1 sequences that resolved to C. albicans from each sample, with the ITS1 sequences from the genome information for C. albicans strains SC5314 and WO-1, displayed 99% identity. The second most abundant C. albicans ITS1 454 sequence amplified from sputum DNA differed from the most abundant ITS1 sequence by one missing “T” at position 60, which most likely reflects a 454 homopolymer sequencing error. By way of comparison, the closely related Candida species C. albicans and C. dubliniensis have ITS1 sequences that have 94.6% identity.
Description of the mycobiome and bacterial microbiome across serial sputum samples
Consistent with published work[4, 6, 65], the bacterial community also displayed stability in the major community members in many subjects (Figure 3), though some changes in the bacteria with the highest relative abundance were observed in subjects #2 and #6 (Figure 3). The Bray-Curtis dissimilarity of the bacterial communities showed a less clear picture than for the fungi, but it still was apparent that the samples from the same subject showed some degree of relatedness. Bacterial community profiles from the same subject were marginally more similar (SPT) but the difference was not significant by a Tukey’s honest significant test (Figure 2B); like for the fungi, the STX samples were not more similar than the UNR samples. There was less of a clear trend in sample clustering in the PCA scores plot, and the larger number of bacterial genera in each sample may have contributed to this complexity (Additional file8: Figure S3).
Analysis of minor members of the mycobiome
The remaining fungal genera represented relatively low abundance populations that occurred in one or more samples (Figure 1B and Additional file5: Table S4 for identities). When using the BLAST-based mycobiome analysis (see Additional file11), the number of taxa identified within a sample ranged from 5 taxa (resolved to 3 species, 1 genus, and 1 cluster of unidentified sequences in the post-treatment sample from subject #2) to 60 taxa (resolved to 43 species, 16 genera, and 1 cluster of unidentified sequences in the pre-treatment sample of subject #1) (Additional file12: Table S5). Minor members (<0.05% of annotated ITS1 sequences) of the mycobiome included other Candida species, Trichosporon (an emerging pathogen in CF subjects[27, 66]), and the human pathogens Rhodotorula, Fusarium, and Penicillium (Additional file5: Table S4).
Analysis of the changes in community structure and absolute abundance before and after antibacterial therapy
We determined total bacterial burden without pre-amplification by targeting the 16S rRNA locus (Figure 6B). On average, we detected ~9.2 × 108 bacterial genome copies/mg dry weight sputum with a maximum of 1.1 × 1010 copies in the pre-treatment sample of subject #2 and minimum of 4.0 × 105 copies in the post-treatment of subject #2. As with the fungal quantification, we did not observe a specific trend in bacterial density in the pre- and post-treatment samples. The bacterial sputum density dropped in four of the patients, most significantly in subject #2 in which the bacterial burden within sputum dropped ~27,500-fold to a level close to the detection limit. In two subjects, the bacterial level went up approximately five-fold. There was no correlation between the changes in bacterial density and changes in fungal density (Figure 6).
Dominant members of fungal communities, which we defined as those that account for more than 80% of all assigned reads per sample (Candida and Malassezia species) (Figure 1A) persisted throughout treatment without significant change in relative abundance suggesting that there is stability within the mycobiome. Many minor members of the fungal community appeared to be ephemeral or sporadic, and thus did not show changes in prevalence or abundance that could be attributed to antibacterial treatment (Figure 5 and Additional file13: Table S7). Based on the depth of sequencing and the rarefaction analyses, the apparent instability of minor members is unlikely an artifact of sequencing depth, but rather likely a reflection of their transient nature. The fungal species we detected within CF sputum mirrors results in other recently published work. For example, denaturing high-performance liquid chromatography analysis of the ITS1 region amplified from fungi within CF sputum also found C. albicans, M. restricta and M. globosa. In a study of four stable CF subjects by Delhaes et al., these fungi were also well represented. While C. albicans is commonly observed in CF sputum, the presence of Malassezia is less well understood in part because Malassezia cannot be cultured without the use of specialized medium. Investigation of the fungal communities on the human skin found M. restricta and M. globosa to be very common, and future work will be required to determine if Malassezia are living in the lung or if the signal is coming from dead cells or DNA. Mournier et al. also detected more fungi using molecular methods when compared to culture-based analyses, and Candida and Malassezia were among them. In addition, they found that A. fumigatus was detected by culture in some samples but not by molecular techniques. While A. fumigatus was not seen in clinical microbiological analyses of the sputa in this study, we were surprised by its complete absence from our datasets. We have detected A. fumigatus in other samples indicating that our methods can amplify its ITS1 region. A recent shotgun metagenomic analysis of CF sputum DNA also found significant levels of C. albicans and A. fumigatus by culture but found no DNA evidence for A. fumigatus and a small number of other Candida sequences. In contrast to what was observed in CF sputum, Charlson et al. analyzed the fungi in the lavages of healthy individuals, and they found that healthy volunteers’ bronchoalveolar lavage (BAL) yielded scant fungal amplification. Although Candida was identified in several healthy subjects’ oropharyngeal washes, it was absent in their BAL.
In this study, we explored whether inpatient treatment, which included multiple antibacterials, impacted the fungal community in a consistent way, perhaps indicating changes in biological niches or environmental conditions upon perturbation of the bacterial community. Our data suggest that the use of multiple antibiotics did not necessarily have strong effects on either the bacterial or fungal communities. While the number of bacterial genera (richness) declined slightly but significantly (Additional file10: Figure S4), there was not a consistently change to the structure of the community with similar relative abundances of dominant species in four out of six subjects (Figure 2 and Additional file12: Figure S5). In these same samples, we did not observe a consistent increase in the number of fungal genera (Additional file10: Figure S4) and the fungal species that dominated the samples did not change markedly (Figure 1A). This lack of difference between pre- and post-samples may well be due to the fact that these subjects have been exposed to years of inhaled and intravenous antibacterial therapy, and the bacteria in their lungs have adapted or evolved in this context. Thus, future studies in subjects without routine exposure to antibiotics are needed in order to assess the effects of antibacterials on the resident fungi. We and others have characterized the interactions between Pseudomonas aeruginosa and C. albicans, including effects on fungal morphology[75, 76], viability[77, 78], metabolism[76, 79], and virulence factor regulation[80, 81], and the high relative abundance of P. aeruginosa in the bacterial communities and C. albicans in the fungal communities supports the hypothesis that these species have the potential to impact one another in the lung.
The therapy used in the treatment of disease exacerbation also had no consistent impact on fungal density in patient sputum samples. Measurements of fungal and bacterial genomes per milligram of sputum showed an inconsistent trend across serial samples collected before and after antibacterial treatment (Figure 6) suggesting the need for either a larger number of subjects or alternative methods for the measurement of live microbes including bacteria and fungi. Previous reports in other body sites have linked antibiotics to increased fungal burden[82–88]. In a recent study, Dolville et al. investigated the gut microbiome of mice over 76 days of treatment with the antibiotics vancomycin, ampicillin, neomycin, and metronidazole and subsequent recovery. They observed that bacterial numbers dropped in abundance more than 3-fold while fungi increased ~40-fold in abundance upon treatment but found that both numbers went back to pre-treatment levels upon recovery. Interestingly, the bacterial microbiome was detectably different from Candida-free controls and Candida became more abundant than at the time point prior to treatment. While we observed that fungal density in sputum had the potential to vary (either increase or decrease), we found no common trend. It is important to note that all of these patients have a significant history of prior antibiotic exposure and thus some adaptation may have occurred. The effects of antibiotic therapy in fungal populations may differ from those obtained upon first exposure to a drug.
Fungal load has been associated with lung disease severity[16, 20, 22] suggesting that it is important to measure fungal density in sputum in addition to the mere presence of specific species. We found a range of fungal loads across samples (Figure 6), and this range is consistent with previous reports. For example, Bauernfeind et al. found between 104 and 107 colony forming units per gram of sputum of Candida using culture methods. Charlson et al. highlight the importance of measuring bacterial and fungal loads and to compare it to background contamination by upper respiratory tract flora to identify lung-enriched taxa. Here, we applied qPCR quantification of the 16S rDNA locus for total bacteria, the 18S rDNA locus for total fungi, and species-specific primers to quantify three Candida species. One limitation of DNA-based quantification protocols is the inability to differentiate between nucleic acids from live and dead species. Methods to inactivate DNA from dead or lysed cells would enhance our understanding of the active microbiota in the lung.
Finally, we report the use of VAMPS for analysis of the ITS1 sequencing results. VAMPS, which has provided a user-friendly tool for the analysis of bacterial communities[2, 91], now has an improved database that will expand the ability to analyze fungal ITS1 sequences from clinical samples. The VAMPS database, derived from the UNITE database, contains most of the available ITS1 sequences. However, further curation will improve taxonomic information in existing entries, and the inclusion of additional high-quality, taxonomically resolved sequences that further aid in the analysis and taxon assignment in future human mycobiome studies.
The data presented in this report found stable fungal communities in the sputum of six individuals with CF. Though to a lesser extent, stability was also mirrored in the bacterial communities. We hypothesize that the resident bacteria have developed a tolerance or resistance to antimicrobials due to repeated prior exposure, and thus, little perturbation of co-colonizing fungi is observed. In addition to numerically dominant fungal populations that are predominantly members of the genus Candida, there are less abundant fungi, many of which were less stable across serial samples from the same subject. To enhance future analyses of the mycobiome, we developed and modified a fungal ITS1 database for use in VAMPS, a publically available tool for the analyses of visualization of data on microbial communities. Future studies will continue to investigate the effects of antibacterial drugs on fungal and bacterial organisms and their densities in sputum.
Study participants and ethics statement
This study was performed in accordance with the guidelines set by the Committee for the Protection of Human Subjects (CPHS) and approved by the Institutional Review Board (IRB) at Dartmouth College. Sputum samples were obtained from six subjects over the age of 18 with a confirmed diagnosis of CF and who had previous cultures positive for P. aeruginosa after written, informed consent. Information about the study participants is provided in Additional file1: Table S1.
Serial sputum samples were collected from six CF subjects at the time of hospitalization for treatment of a respiratory disease exacerbation and after a complete course of intravenous antibacterial therapy, which often included the use of multiple antibiotics (Additional file1: Table S1). For some subjects, an intermediate sample during therapy was also collected. Study participants expectorated sputum into a sterile specimen cup, and the sputum sample was stored at -80°C. In preparation for DNA extraction, the sample was thawed and homogenized by passing the sputum multiple times through needles of increasing gauges. Sputum (500 μl) was dispensed into a 2-ml screw cap tube and frozen at -80°C for at least 1 h prior to lyophilization for at least 12 h. The dry weight of the sputum was determined before DNA extraction.
Cells within the dried sputum sample were disrupted using a mixture sterile glass beads (equivalent amounts of 0.1-, 0.5-, and 1-mm diameter beads) in a bead beater (Biospec, Bartlesville, OK, USA); the tubes were agitated four times for 30 s each time with a cooling step in between pulses. The ground samples were resuspended in 300 μl of TE + DTT (TE amended with DTT at a final concentration of 0.08% added from a 2% stock solution) containing lysozyme (3 mg/ml) and lyticase (10 U/ml) and incubated for 30 min at 37°C. Cell lysis buffer (500 μl) (Qiagen Puregene Core Kit B, QIAGEN Inc., Valencia, CA, USA) was added, and the mixture was incubated for 15 min at 80°C. To remove RNA, RNase (1.5 μl) (QIAGEN Inc.) was added and the samples were incubated for 30 min at 37°C. Lysates were chilled on ice for 1 min, 200 μl of Protein Precipitation Solution (Qiagen Puregene Core Kit B, QIAGEN Inc.) was added, and the solutions were mixed vigorously for 20 s. Cell debris was sedimented by centrifugation at 13,000 rpm for 3 min, and the supernatant was transferred to a new 1.5-ml tube prior to addition of 600 μl of 100% isopropanol. After mixing by inversion, the DNA was precipitated by centrifugation at 13,000 rpm for 20 min. The DNA pellet was washed with 300 μl of 70% ethanol and air dried before resuspension in 100–200 μl of DNA hydration solution (Qiagen Puregene Core Kit B, QIAGEN Inc.). The DNA concentrations were measured using a Nanodrop.
Fungal ITS1 amplification and sequencing
The ITS1 region was amplified from sputum template DNA using the primers ITS1F and ITS1R (Additional file7: Table S2) in two rounds of PCR. The reaction components were 1× Qiagen PCR buffer, 0.75 mM MgCl2, 7% (v/v) dimethyl sulfoxide (DMSO), 0.2 mM deoxynucleotide triphosphates (dNTPs), 0.4 μM primers, 1.25 units Taq polymerase (QIAGEN Inc.), 100 ng of template DNA, and water in a 25-μl reaction volume. Cycling conditions were 95°C 5 min; 95°C 30 s, 50°C 30 s, 72°C 45 s for 10 cycles; and 72°C 2 min. PCR products were purified using the MinElute PCR Purification Kit (QIAGEN Inc.) and resuspended in 10 μl of EB buffer. The entire 10 μl volume was used as template for a second round of PCR using primers containing adapter and barcode sequences for Roche GS-FLX amplicon pyrosequencing (Additional file2: Table S2). Reaction components and cycling conditions were as described above, but with 30 cycles. Amplified DNA was purified with 0.7 volumes of Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA) and quantified on a Bioanalyzer 2100 (Agilent, Palo Alto, CA, USA). Purified amplicon libraries were sequenced using the Titanium amplicon protocol on a GS-FLX (Roche Applied Science, Indianapolis, IN, USA) following the methods outlined in. Sequence data were quality filtered as described in[54, 92]. We identified and removed chimeric reads by querying against the curated reference ITS1 database, as well as using the sequences as seeds in UCHIME de novo. The raw sequencing results of the ITS1 analysis are deposited in the Sequence Read Archive (NCBI BioProject PRJNA246028; National Center for Biotechnology Information (NCBI), US National Library of Medicine, Bethesda, MD, USA) and are summarized in Additional file5: Table S4.
For the MiSeq analysis of the ITS1 region amplified from clinical isolates of Candida spp., template DNA was isolated using the Epicentre Masterpure Yeast DNA Purification Kit. Primers are listed in Additional file7: Table S2. The amplification and sequencing methods used are comparable to those reported.
Bacterial 16S v6 amplification and sequencing
The bacterial communities from the pre- and post-antibacterial treated CF subjects were characterized by Illumina HiSeq 16S rRNA V6 amplicon sequencing (Illumina Inc., San Diego, CA, USA) as previously described by Eren et al.. Briefly, we amplified the bacterial V6 region from each sample in triplicate 33 μL reactions, pooled the triplicates, and verified successful amplification on a Caliper LabChip GX (Perkin Elmer, Waltham, MA, USA). We cleaned amplicon products with a QIAQuick clean up assay (QIAGEN Inc.), and used fusion V6 primers to barcode each sample in a multiplexing strategy. Final barcoded products were quantitated with a Picogreen assay (Invitrogen/Life Technologies, Grand Island, NY, USA) and pooled in equimolar concentrations with respect to the target size. All PCR reactions contained 1× HiFi Buffer, 2 mM MgSO4, 0.02 U/μL Platinum Taq polymerase (Invitrogen), 0.2 mM each dNTPs (ThermoFisher Scientific, Milwaukee, WI, USA), and up to 5 ng template. The final pool was size selected (Pippin Prep, SageScience, Beverly, MA, USA), quantitated (Kapa Biosystems, Woburn, MA, USA) and sequenced on a HiSeq 1000 in a 2 × 100 nt sequencing run. Quality filtering and chimera checks were performed as described above. The raw sequencing results of the bacterial 16S v6 analysis are deposited in the Sequence Read Archive (NCBI BioProject PRJNA256169; National Center for Biotechnology Information (NCBI), US National Library of Medicine, Bethesda, MD, USA) and are summarized in Additional file9: Table S6.
Microbiome and mycobiome clustering and taxonomic assignment
The website http://vamps.mbl.edu provides access to individual reads, taxon assignments, and descriptions of individual clusters. To create the databases that are accessed through VAMPS researchers at the Marine Biology Laboratories at Woods Hole merged paired-end 100 bp reads into single consensus reads that completely covered the V6 region, and they assigned taxonomy with GAST against a curated SILVA database. To compare operational taxonomic unit (OTU) clustering performance, they employed the default UCLUST method at a 97% similarity threshold with minimum cluster size of 2 using quantitative insights into microbial ecology (QIIME) (v1.5).
To verify adequate sequencing depth of the mycobiome and the microbiome, we performed rarefaction analysis of the two datasets. For both rarefaction analyses, we used raw counts of all subject samples generated by VAMPS with the taxonomy depth selector set to “species” to create a read matrix. This table was rarefied using an R-based protocol described on the webpage “Individual Based Rarefaction using R-package” (http://www.jennajacobs.org/R/rarefaction.html). The resulting table was transferred into GraphPad Prism 6 to plot the rarefaction curves.
Euler and Venn diagrams
To visualize the impact of antibacterial treatment on the prevalence of fungal and bacterial taxa detected in CF subjects either immediately after an exacerbation and before starting an antibacterial therapy (pre), approximately 2 weeks afterwards (post) or while hospitalized (inter), we used the Additional file5: Tables S4 and Additional file9: Table S6 to generate the diagrams. First, we copied the data for pre, post and if applicable inter samples into the list-fields of an online Venn diagram generator (http://www.bioinformatics.lu/venn.php). Using the R package “RVennDiagram” (http://cran.r-project.org/web/packages/VennDiagram/index.html), we generated Euler (“draw.pairwise.venn”) and Venn (“draw.triple.venn”) diagrams following the instructions in the manual (http://cran.r-project.org/web/packages/VennDiagram/VennDiagram.pdf).
Quantification of fungi, Candida species, and bacteria
In order to quantify minor members of the mycobiome, we performed analyses using either sputum DNA or pre-amplified DNA derived from ten rounds of PCR starting with 100 ng of sputum DNA using primers specific to a conserved region of the 18S rDNA (FungiQuant primers) or Candida species-specific primers (Additional file7: Table S2). We confirmed that the pre-amplification step was linear with respect to the initial sample; we performed multiple control experiments using genomic DNA isolated from different fungal species at concentrations ranging from 30 pg/μl to 300 ng/μl. The preamplification allowed for the analysis of minor members of the mycobiota. The reaction components of the first PCR were 1× Qiagen PCR buffer, 0.75 mM MgCl2, 7% (v/v) DMSO, 0.2 mM dNTPs, 0.4 μM primers, 1.25 units Taq polymerase (QIAGEN Inc.), 100 ng sputum DNA, and water in a 25 μl reaction volume. Cycling conditions were 95°C 5 min; 95°C 30 s, 50°C 30 s, 72°C 45 sec for 10 cycles; and 72°C 2 min. These PCR products were purified using the MinElute PCR Purification Kit (QIAGEN Inc.), and resuspended in 10 μl of EB buffer. Quantitative PCR (qPCR) was conducted in 20 μl reaction volumes with the iQ SYBR Green Supermix (Bio-Rad Laboratories, Hercules, CA, USA). Reactions contained 2 μl of the amplicon suspension from the first round PCR and 0.2 μM of the appropriate primers. To quantify total fungal burden, we used the primers FungiQuant_RT_F and FungiQuant_RT_R. The aforementioned species-specific qPCR primers were designed to assess Candida burdens based on gene sequences available from the Candida Genome Database (http://www.candidagenome.org/). Genes used for primer design were CAWG_05066 (C. albicans), CPAR2_301290 (C. parapsilosis), Cd36_16280 (C. dubliniensis) and CTRG_03824 (C. tropicalis) and were selected based on conservation within species and a lack of cross reactivity with bacterial, fungal, and human sequences. qPCR was performed utilizing a CFX96 real-time PCR detection system combined with a C1000 thermal cycler (Bio-Rad Laboratories). All PCRs were done in duplicate, and data were analyzed with the CFX96 System gene expression software. Melt curve analysis was performed after the PCR was complete to confirm the absence of nonspecific amplification products. Standard curves containing a known numbers of genome equivalents were used to calculate genome numbers per mg dry weight.
For total bacterial quantification, 100 ng of sputum DNA was used as the PCR reaction template. To generate a standard curve for quantification purposes we prepared 10-fold dilutions of P. aeruginosa genomic DNA isolated from 102 to 106 cells. The 16S rDNA genes were amplified by qPCR, using the universal bacterial primers “total bacteria_F” and “total bacteria_R” at a final concentration of 0.2 μM[101, 102] (Additional file7: Table S2). qPCR was conducted in 20 μl reaction volumes with the iQ SYBR Green Supermix (Bio-Rad Laboratories). qPCRs for each sputum sample were performed and analyzed on the CFX96 PCR system as described above.
Bray-Curtis dissimilarity was measured using the ecodist and vegan packages in R statistical software. Significant differences between groups were established using analysis of variance with Tukey’s honest significant difference post-test. Tests producing a P value less than 0.05 were deemed significant. The program Prism 6 (GraphPad, San Diego, CA, USA) was used for remaining statistical tests.
Global Alignment for Sequence Taxonomy
internal transcribed spacer 1
Marine Biology Laboratories at Woods Hole
operational taxonomic unit
polymerase chain reaction
quantitative insights into microbial ecology
ribosomal RNA genes
visualization and analysis of microbial population structures.
The authors would like to thank Dr. Jason Stajich and Dr. Greg (Dai) Gu for their support creating the BLAST-based mycobiome analysis tool and for their constructive comments. Furthermore, we would like to thank all other researchers in the Hogan lab as well as the anonymous reviewers for their support and useful suggestions. Research reported in this publication was supported by grants from the National Institutes of Health to DAH (R01GM108492), an award from the Hitchcock Foundation (DAH), and a pilot SYNERGY grant (AA and DAH) and a grant from the Alfred P. Sloan Foundation’s Microbiology of the Built Environment (MLS). Support was also provided by 2R37AI83256-06 (G.A.O.). Support from the Cystic Fibrosis Foundation Research Development Program (STANTO07R0) and an Institutional Development Award (IDeA) from the National Institute of General Medical Sciences of the National Institutes of Health under grant number P30GM106394 facilitated these studies through support of the Dartmouth Lung Biology Translational Research Core. The work presented is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
- Tunney MM, Klem ER, Fodor AA, Gilpin DF, Moriarty TF, McGrath SJ, Muhlebach MS, Boucher RC, Cardwell C, Doering G, Elborn JS, Wolfgang MC: Use of culture and molecular analysis to determine the effect of antibiotic treatment on microbial community diversity and abundance during exacerbation in patients with cystic fibrosis. Thorax. 2011, 66 (7): 579-584. 10.1136/thx.2010.137281.PubMedGoogle Scholar
- Filkins LM, Hampton TH, Gifford AH, Gross MJ, Hogan DA, Sogin ML, Morrison HG, Paster BJ, O’Toole GA: Prevalence of streptococci and increased polymicrobial diversity associated with cystic fibrosis patient stability. J Bacteriol. 2012, 194 (17): 4709-4717. 10.1128/JB.00566-12.PubMedPubMed CentralGoogle Scholar
- Spilker T, Vandamme P, Lipuma JJ: A multilocus sequence typing scheme implies population structure and reveals several putative novel Achromobacter species. J Clin Microbiol. 2012, 50 (9): 3010-3015. 10.1128/JCM.00814-12.PubMedPubMed CentralGoogle Scholar
- Zhao J, Schloss PD, Kalikin LM, Carmody LA, Foster BK, Petrosino JF, Cavalcoli JD, VanDevanter DR, Murray S, Li JZ, Young VB, LiPuma JJ: Decade-long bacterial community dynamics in cystic fibrosis airways. Proc Natl Acad Sci U S A. 2012, 109 (15): 5809-5814. 10.1073/pnas.1120577109.PubMedPubMed CentralGoogle Scholar
- Sibley CD, Parkins MD, Rabin HR, Duan K, Norgaard JC, Surette MG: A polymicrobial perspective of pulmonary infections exposes an enigmatic pathogen in cystic fibrosis patients. Proc Natl Acad Sci U S A. 2008, 105 (39): 15070-15075. 10.1073/pnas.0804326105.PubMedPubMed CentralGoogle Scholar
- Fodor AA, Klem ER, Gilpin DF, Elborn JS, Boucher RC, Tunney MM, Wolfgang MC: The adult cystic fibrosis airway microbiota is stable over time and infection type, and highly resilient to antibiotic treatment of exacerbations. PLoS One. 2012, 7 (9): e45001-10.1371/journal.pone.0045001.PubMedPubMed CentralGoogle Scholar
- Stressmann FA, Rogers GB, van der Gast CJ, Marsh P, Vermeer LS, Carroll MP, Hoffman L, Daniels TW, Patel N, Forbes B, Bruce KD: Long-term cultivation-independent microbial diversity analysis demonstrates that bacterial communities infecting the adult cystic fibrosis lung show stability and resilience. Thorax. 2012, 67 (10): 867-873. 10.1136/thoraxjnl-2011-200932.PubMedGoogle Scholar
- Leclair LW, Hogan DA: Mixed bacterial-fungal infections in the CF respiratory tract. Med Mycol. 2010, 48 (O1): S125-S132. 10.3109/13693786.2010.521522.PubMedGoogle Scholar
- Morales DK, Hogan DA: Candida albicans interactions with bacteria in the context of human health and disease. PLoS Pathog. 2010, 6 (4): e1000886-10.1371/journal.ppat.1000886.PubMedPubMed CentralGoogle Scholar
- Peleg AY, Hogan DA, Mylonakis E: Medically important bacterial-fungal interactions. Nat Rev Microbiol. 2010, 8 (5): 340-349. 10.1038/nrmicro2313.PubMedGoogle Scholar
- Hoppe JE, Theurer-Mainka U, Stern M: Comparison of three methods for culturing throat swabs from cystic fibrosis patients. J Clin Microbiol. 1995, 33 (7): 1896-1898.PubMedPubMed CentralGoogle Scholar
- Bakare N, Rickerts V, Bargon J, Just-Nubling G: Prevalence of Aspergillus fumigatus and other fungal species in the sputum of adult patients with cystic fibrosis. Mycoses. 2003, 46 (1–2): 19-23.PubMedGoogle Scholar
- Haase G, Skopnik H, Groten T, Kusenbach G, Posselt HG: Long-term fungal cultures from sputum of patients with cystic fibrosis. Mycoses. 1991, 34 (9–10): 373-376.PubMedGoogle Scholar
- Hughes WT, Kim HK: Mycoflora in cystic fibrosis: some ecologic aspects of Pseudomonas aeruginosa and Candida albicans. Mycopathol Mycol Appl. 1973, 50 (3): 261-269. 10.1007/BF02053377.PubMedGoogle Scholar
- Bauernfeind A, Bertele RM, Harms K, Horl G, Jungwirth R, Petermuller C, Przyklenk B, Weisslein-Pfister C: Qualitative and quantitative microbiological analysis of sputa of 102 patients with cystic fibrosis. Infection. 1987, 15 (4): 270-277. 10.1007/BF01644137.PubMedGoogle Scholar
- Chotirmall SH, O’Donoghue E, Bennett K, Gunaratnam C, O’Neill SJ, McElvaney NG: Sputum Candida albicans presages FEV1 decline and hospitalized exacerbations in cystic fibrosis. Chest. 2010, 138 (5): 1186-1195. 10.1378/chest.09-2996.PubMedGoogle Scholar
- Doern GV, Brogden-Torres B: Optimum use of selective plated media in primary processing of respiratory tract specimens from patients with cystic fibrosis. J Clin Microbiol. 1992, 30 (10): 2740-2742.PubMedPubMed CentralGoogle Scholar
- Güngör O, Tamay Z, Güler N, Erturan Z: Frequency of fungi in respiratory samples from Turkish cystic fibrosis patients. Mycoses. 2012, 56 (2): 123-129.PubMedGoogle Scholar
- Delhaes L, Monchy S, Frealle E, Hubans C, Salleron J, Leroy S, Prevotat A, Wallet F, Wallaert B, Dei-Cas E, Sime-Ngando T, Chabe M, Viscogliosi E: The airway microbiota in cystic fibrosis: a complex fungal and bacterial community–implications for therapeutic management. PLoS One. 2012, 7 (4): e36313-10.1371/journal.pone.0036313.PubMedPubMed CentralGoogle Scholar
- Navarro J, Rainisio M, Harms HK, Hodson ME, Koch C, Mastella G, Strandvik B, McKenzie SG: Factors associated with poor pulmonary function: cross-sectional analysis of data from the ERCF. European Epidemiologic Registry of Cystic Fibrosis. Eur Respir J. 2001, 18 (2): 298-305. 10.1183/09031936.01.00068901.PubMedGoogle Scholar
- Storey DG, Ujack EE, Rabin HR, Mitchell I: Pseudomonas aeruginosa lasR transcription correlates with the transcription of lasA, lasB, and toxA in chronic lung infections associated with cystic fibrosis. Infect Immun. 1998, 66 (6): 2521-2528.PubMedPubMed CentralGoogle Scholar
- Amin R, Dupuis A, Aaron SD, Ratjen F: The effect of chronic infection with Aspergillus fumigatus on lung function and hospitalization in patients with cystic fibrosis. Chest. 2010, 137 (1): 171-176. 10.1378/chest.09-1103.PubMedGoogle Scholar
- Masoud-Landgraf L, Badura A, Eber E, Feierl G, Marth E, Buzina W: Modified culture method detects a high diversity of fungal species in cystic fibrosis patients. Med Mycol. 2013, 52 (2): 179-186.Google Scholar
- Parize P, Billaud S, Bienvenu AL, Bourdy S, le Pogam MA, Reix P, Picot S, Robert R, Lortholary O, Bouchara JP, Durieu I: Impact of Scedosporium apiospermum complex seroprevalence in patients with cystic fibrosis. J Cyst Fibros. 2014, doi:10.1016/j.jcf.2014.01.011Google Scholar
- Bernhardt A, Sedlacek L, Wagner S, Schwarz C, Wurstl B, Tintelnot K: Multilocus sequence typing of Scedosporium apiospermum and Pseudallescheria boydii isolates from cystic fibrosis patients. J Cyst Fibros. 2013, 12 (6): 592-598. 10.1016/j.jcf.2013.05.007.PubMedGoogle Scholar
- Packeu A, Lebecque P, Rodriguez-Villalobos H, Boeras A, Hendrickx M, Bouchara JP, Symoens F: Molecular typing and antifungal susceptibility of Exophiala isolates from patients with cystic fibrosis. J Med Microbiol. 2012, 61 (Pt 9): 1226-1233.PubMedGoogle Scholar
- Hickey PW, Sutton DA, Fothergill AW, Rinaldi MG, Wickes BL, Schmidt HJ, Walsh TJ: Trichosporon mycotoxinivorans, a novel respiratory pathogen in patients with cystic fibrosis. J Clin Microbiol. 2009, 47 (10): 3091-3097. 10.1128/JCM.00460-09.PubMedPubMed CentralGoogle Scholar
- Nagano Y, Elborn JS, Millar BC, Walker JM, Goldsmith CE, Rendall J, Moore JE: Comparison of techniques to examine the diversity of fungi in adult patients with cystic fibrosis. Med Mycol. 2010, 48 (1): 166-76 e1. 10.3109/13693780903127506.PubMedGoogle Scholar
- Respiratory cultures from Cystic Fibrosis patients, in Clinical Microbiology Procedures Handbook. Edited by: Garcia L. 2010, ASM PressGoogle Scholar
- Price KE, Hampton TH, Gifford AH, Dolben EL, Hogan DA, Morrison HG, Sogin ML, O’Toole GA: Unique microbial communities persist in individual cystic fibrosis patients throughout a clinical exacerbation. Microbiome. 2013, 1 (1): 27-10.1186/2049-2618-1-27.PubMedPubMed CentralGoogle Scholar
- Jubin V, Ranque S, Stremler N, Le B, Sarles J, Dubus JC: Risk factors for Aspergillus colonization and allergic bronchopulmonary aspergillosis in children with cystic fibrosis. Pediatr Pulmonol. 2010, 45 (8): 764-771. 10.1002/ppul.21240.PubMedGoogle Scholar
- Burns J, Van Dalfsen J, Shawar R, Otto K, Garber R, Quan J, Montgomery A, Albers G, Ramsey B, Smith A: Effect of intermittent administration of inhaled tobramycin on respiratory microbial flora in patients with cystic fibrosis. J Infect Dis. 1999, 179: 1190-1196. 10.1086/314727.PubMedGoogle Scholar
- Cheer SM, Waugh J, Noble S: Inhaled tobramycin (TOBI): a review of its use in the management of Pseudomonas aeruginosa infections in patients with cystic fibrosis. Drugs. 2003, 63 (22): 2501-2520. 10.2165/00003495-200363220-00015.PubMedGoogle Scholar
- Sudfeld CR, Dasenbrook EC, Merz WG, Carroll KC, Boyle MP: Prevalence and risk factors for recovery of filamentous fungi in individuals with cystic fibrosis. J Cyst Fibros. 2010, 9 (2): 110-116. 10.1016/j.jcf.2009.11.010.PubMedPubMed CentralGoogle Scholar
- Kerr J: Inhibition of fungal growth by Pseudomonas aeruginosa and Pseudomonas cepacia isolated from patients with cystic fibrosis. J Infect. 1994, 28: 305-310. 10.1016/S0163-4453(94)91943-7.PubMedGoogle Scholar
- Mastella G, Rainisio M, Harms HK, Hodson ME, Koch C, Navarro J, Strandvik B, McKenzie SG: Allergic bronchopulmonary aspergillosis in cystic fibrosis. A European epidemiological study. Epidemiologic Registry of Cystic Fibrosis. Eur Respir J. 2000, 16 (3): 464-471. 10.1034/j.1399-3003.2000.016003464.x.PubMedGoogle Scholar
- Krcmery V, Matejicka F, Pichnova E, Jurga L, Sulcova M, Kunova A, West D: Documented fungal infections after prophylaxis or therapy with wide spectrum antibiotics: relationship between certain fungal pathogens and particular antimicrobials?. J Chemother. 1999, 11 (5): 385-390. 10.1179/joc.19126.96.36.1995.PubMedGoogle Scholar
- Hebert C, Villaran R, Tolentino J, Best L, Boonlayangoor S, Pitrak D, Lin M, Weber SG: Prior antimicrobial exposure and the risk for bloodstream infection with fluconazole-non-susceptible Candida strains. Scand J Infect Dis. 2010, 42 (6–7): 506-509.PubMedGoogle Scholar
- Wey SB, Mori M, Pfaller MA, Woolson RF, Wenzel RP: Risk factors for hospital-acquired candidemia. A matched case–control study. Arch Intern Med. 1989, 149 (10): 2349-2353. 10.1001/archinte.1989.00390100145030.PubMedGoogle Scholar
- Charles PE, Dalle F, Aube H, Doise JM, Quenot JP, Aho LS, Chavanet P, Blettery B: Candida spp. colonization significance in critically ill medical patients: a prospective study. Intensive Care Med. 2005, 31 (3): 393-400. 10.1007/s00134-005-2571-y.PubMedGoogle Scholar
- Soysa NS, Samaranayake LP, Ellepola AN: Antimicrobials as a contributory factor in oral candidosis–a brief overview. Oral Dis. 2008, 14 (2): 138-143. 10.1111/j.1601-0825.2006.01357.x.PubMedGoogle Scholar
- Xu J, Schwartz K, Bartoces M, Monsur J, Severson RK, Sobel JD: Effect of antibiotics on vulvovaginal candidiasis: a MetroNet study. J Am Board Fam Med. 2008, 21 (4): 261-268. 10.3122/jabfm.2008.04.070169.PubMedGoogle Scholar
- Samonis G, Anastassiadou H, Dassiou M, Tselentis Y, Bodey GP: Effects of broad-spectrum antibiotics on colonization of gastrointestinal tracts of mice by Candida albicans. Antimicrob Agents Chemother. 1994, 38 (3): 602-603. 10.1128/AAC.38.3.602.PubMedPubMed CentralGoogle Scholar
- Samonis G, Gikas A, Toloudis P, Maraki S, Vrentzos G, Tselentis Y, Tsaparas N, Bodey G: Prospective study of the impact of broad-spectrum antibiotics on the yeast flora of the human gut. Eur J Clin Microbiol Infect Dis. 1994, 13 (8): 665-667. 10.1007/BF01973996.PubMedGoogle Scholar
- Giuliano M, Barza M, Jacobus NV, Gorbach SL: Effect of broad-spectrum parenteral antibiotics on composition of intestinal microflora of humans. Antimicrob Agents Chemother. 1987, 31 (2): 202-206. 10.1128/AAC.31.2.202.PubMedPubMed CentralGoogle Scholar
- Leaw SN, Chang HC, Sun HF, Barton R, Bouchara JP, Chang TC: Identification of medically important yeast species by sequence analysis of the internal transcribed spacer regions. J Clin Microbiol. 2006, 44 (3): 693-699. 10.1128/JCM.44.3.693-699.2006.PubMedPubMed CentralGoogle Scholar
- Pryce TM, Palladino S, Kay ID, Coombs GW: Rapid identification of fungi by sequencing the ITS1 and ITS2 regions using an automated capillary electrophoresis system. Med Mycol. 2003, 41 (5): 369-381. 10.1080/13693780310001600435.PubMedGoogle Scholar
- Ghannoum MA, Jurevic RJ, Mukherjee PK, Cui F, Sikaroodi M, Naqvi A, Gillevet PM: Characterization of the oral fungal microbiome (mycobiome) in healthy individuals. PLoS Pathog. 2010, 6 (1): e1000713-10.1371/journal.ppat.1000713.PubMedPubMed CentralGoogle Scholar
- Ferrer C, Colom F, Frases S, Mulet E, Abad JL, Alio JL: Detection and identification of fungal pathogens by PCR and by ITS2 and 5.8S ribosomal DNA typing in ocular infections. J Clin Microbiol. 2001, 39 (8): 2873-2879. 10.1128/JCM.39.8.2873-2879.2001.PubMedPubMed CentralGoogle Scholar
- Borman AM, Linton CJ, Miles SJ, Johnson EM: Molecular identification of pathogenic fungi. J Antimicrob Chemother. 2008, 61 (Suppl 1): i7-i12.PubMedGoogle Scholar
- Spiess B, Seifarth W, Hummel M, Frank O, Fabarius A, Zheng C, Morz H, Hehlmann R, Buchheidt D: DNA microarray-based detection and identification of fungal pathogens in clinical samples from neutropenic patients. J Clin Microbiol. 2007, 45 (11): 3743-3753. 10.1128/JCM.00942-07.PubMedPubMed CentralGoogle Scholar
- Kittelmann S, Naylor GE, Koolaard JP, Janssen PH: A proposed taxonomy of anaerobic fungi (class neocallimastigomycetes) suitable for large-scale sequence-based community structure analysis. PLoS One. 2012, 7 (5): e36866-10.1371/journal.pone.0036866.PubMedPubMed CentralGoogle Scholar
- Huse SM, Mark DB, Welch A, Voorhis A, Shipunova HG, Morrison AM, Eren , Sogin ML: VAMPS: a website for visualization and analysis of microbial population structures. BMC Bioinformatics. 2014, 15 (1): 41-10.1186/1471-2105-15-41.PubMedPubMed CentralGoogle Scholar
- Huse SM, Dethlefsen L, Huber JA, Mark D, Welch DA, Relman , Sogin ML: Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing. PLoS Genet. 2008, 4 (11): e1000255-10.1371/journal.pgen.1000255.PubMedPubMed CentralGoogle Scholar
- Huse SM, Welch DM, Morrison HG, Sogin ML: Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ Microbiol. 2010, 12 (7): 1889-1898. 10.1111/j.1462-2920.2010.02193.x.PubMedPubMed CentralGoogle Scholar
- Huber JA, Mark Welch DB, Morrison HG, Huse SM, Neal PR, Butterfield DA, Sogin ML: Microbial population structures in the deep marine biosphere. Science. 2007, 318 (5847): 97-100. 10.1126/science.1146689.PubMedGoogle Scholar
- Madan JC, Koestler DC, Stanton BA, Davidson L, Moulton LA, Housman ML, Moore JH, Guill MF, Morrison HG, Sogin ML, Hampton TH, Karagas MR, Palumbo PE, Foster JA, Hibberd PL, O’Toole GA: Serial analysis of the gut and respiratory microbiome in cystic fibrosis in infancy: interaction between intestinal and respiratory tracts and impact of nutritional exposures. MBio. 2012, 3 (4): e00251-12.PubMedPubMed CentralGoogle Scholar
- Koljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AF, Bahram M, Bates ST, Bruns TD, Bengtsson-Palme J, Callaghan TM, Douglas B, Drenkhan T, Eberhardt U, Duenas M, Grebenc T, Griffith GW, Hartmann M, Kirk PM, Kohout P, Larsson E, Lindahl BD, Lucking R, Martin MP, Matheny PB, Nguyen NH, Niskanen T, Oja J, Peay KG, Peintner U, Peterson M: Towards a unified paradigm for sequence-based identification of fungi. Mol Ecol. 2013, 22 (21): 5271-5277. 10.1111/mec.12481.PubMedGoogle Scholar
- Bruns T, White T, Taylor J: Fungal molecular systematics. Annu Rev Ecol Syst. 1991, 22: 525-564. 10.1146/annurev.es.22.110191.002521.Google Scholar
- Yao C, Frederiksen RA, Magill CW: Length heterogeneity in ITS 2 and the methylation status of CCGG and GCGC sites in the rRNA genes of the genus Peronosclerospora. Curr Genet. 1992, 22 (5): 415-420. 10.1007/BF00352443.PubMedGoogle Scholar
- Longo AV, Rodriguez D, Da Silva D, Leite LF, Toledo C, Mendoza A, Burrowes PA, Zamudio KR: ITS1 copy number varies among Batrachochytrium dendrobatidis strains: implications for qPCR estimates of infection intensity from field-collected amphibian skin swabs. PLoS One. 2013, 8 (3): e59499-10.1371/journal.pone.0059499.PubMedPubMed CentralGoogle Scholar
- Nilsson RH, Kristiansson E, Ryberg M, Hallenberg N, Larsson KH: Intraspecific ITS variability in the kingdom fungi as expressed in the international sequence databases and its implications for molecular species identification. Evol Bioinform Online. 2008, 4: 193-201.PubMedPubMed CentralGoogle Scholar
- Bittar F, Rolain JM: Detection and accurate identification of new or emerging bacteria in cystic fibrosis patients. Clin Microbiol Infect. 2010, 16 (7): 809-820. 10.1111/j.1469-0691.2010.03236.x.PubMedGoogle Scholar
- Sibley CD, Rabin H, Surette MG: Cystic fibrosis: a polymicrobial infectious disease. Future Microbiol. 2006, 1 (1): 53-61. 10.2217/174609188.8.131.52.PubMedGoogle Scholar
- Zemanick ET, Harris JK, Wagner BD, Robertson CE, Sagel SD, Stevens MJ, Accurso FJ, Laguna TA: Inflammation and airway microbiota during cystic fibrosis pulmonary exacerbations. PLoS One. 2013, 8 (4): e62917-10.1371/journal.pone.0062917.PubMedPubMed CentralGoogle Scholar
- Hirschi S, Letscher-Bru V, Pottecher J, Lannes B, Jeung MY, Degot T, Santelmo N, Sabou AM, Herbrecht R, Kessler R: Disseminated Trichosporon mycotoxinivorans, Aspergillus fumigatus, and Scedosporium apiospermum coinfection after lung and liver transplantation in a cystic fibrosis patient. J Clin Microbiol. 2012, 50 (12): 4168-4170. 10.1128/JCM.01928-12.PubMedPubMed CentralGoogle Scholar
- Kothavade RJ, Kura MM, Valand AG, Panthaki MH: Candida tropicalis: its prevalence, pathogenicity and increasing resistance to fluconazole. J Med Microbiol. 2010, 59 (Pt 8): 873-880.PubMedGoogle Scholar
- Davies AN, Brailsford S, Broadley K, Beighton D: Oral yeast carriage in patients with advanced cancer. Oral Microbiol Immunol. 2002, 17 (2): 79-84. 10.1046/j.0902-0055.2001.00095.x.PubMedGoogle Scholar
- Xu J, Mitchell TG: Geographical differences in human oral yeast flora. Clin Infect Dis. 2003, 36 (2): 221-224. 10.1086/345672.PubMedGoogle Scholar
- Silva S, Negri M, Henriques M, Oliveira R, Williams DW, Azeredo J: Candida glabrata, Candida parapsilosis and Candida tropicalis: biology, epidemiology, pathogenicity and antifungal resistance. FEMS Microbiol Rev. 2012, 36 (2): 288-305. 10.1111/j.1574-6976.2011.00278.x.PubMedGoogle Scholar
- Mounier J, Gouello A, Keravec M, Le Gal S, Pacini G, Debaets S, Nevez G, Rault G, Barbier G, Hery-Arnaud G: Use of denaturing high-performance liquid chromatography (DHPLC) to characterize the bacterial and fungal airway microbiota of cystic fibrosis patients. J Microbiol. 2014, 52 (4): 307-314. 10.1007/s12275-014-3425-5.PubMedGoogle Scholar
- Findley K, Oh J, Yang J, Conlan S, Deming C, Meyer JA, Schoenfeld D, Nomicos E, Park M, Kong HH, Segre JA: Topographic diversity of fungal and bacterial communities in human skin. Nature. 2013, 498 (7454): 367-370. 10.1038/nature12171.PubMedPubMed CentralGoogle Scholar
- Hauser PM, Bernard T, Greub G, Jaton K, Pagni M, Hafen GM: Microbiota present in cystic fibrosis lungs as revealed by whole genome sequencing. PLoS One. 2014, 9 (3): e90934-10.1371/journal.pone.0090934.PubMedPubMed CentralGoogle Scholar
- Charlson ES, Diamond JM, Bittinger K, Fitzgerald AS, Yadav A, Haas AR, Bushman FD, Collman RG: Lung-enriched organisms and aberrant bacterial and fungal respiratory microbiota after lung transplant. Am J Respir Crit Care Med. 2012, 186 (6): 536-545. 10.1164/rccm.201204-0693OC.PubMedPubMed CentralGoogle Scholar
- Hogan DA, Vik A, Kolter R: A Pseudomonas aeruginosa quorum-sensing molecule influences Candida albicans morphology. Mol Microbiol. 2004, 54 (5): 1212-1223. 10.1111/j.1365-2958.2004.04349.x.PubMedGoogle Scholar
- Morales DK, Grahl N, Okegbe C, Dietrich LE, Jacobs NJ, Hogan DA: Control of Candida albicans metabolism and biofilm formation by Pseudomonas aeruginosa phenazines. MBio. 2013, 4 (1): e00526-12.PubMedPubMed CentralGoogle Scholar
- Hogan DA, Kolter R: Pseudomonas-Candida interactions: an ecological role for virulence factors. Science. 2002, 296 (5576): 2229-2232. 10.1126/science.1070784.PubMedGoogle Scholar
- Gibson J, Sood A, Hogan DA: Pseudomonas aeruginosa-Candida albicans interactions: localization and fungal toxicity of a phenazine derivative. Appl Environ Microbiol. 2009, 75 (2): 504-513. 10.1128/AEM.01037-08.PubMedPubMed CentralGoogle Scholar
- Lindsay AK, Morales DK, Liu Z, Grahl N, Zhang A, Willger SD, Myers LC, Hogan DA: Analysis of Candida albicans mutants defective in the Cdk8 module of mediator reveal links between metabolism and biofilm formation. PLoS Genet. 2014, 10 (10): e1004567-10.1371/journal.pgen.1004567.PubMedPubMed CentralGoogle Scholar
- Cugini C, Morales DK, Hogan DA: Candida albicans-produced farnesol stimulates Pseudomonas quinolone signal production in LasR-defective Pseudomonas aeruginosa strains. Microbiology. 2010, 156 (Pt 10): 3096-3107.PubMedPubMed CentralGoogle Scholar
- McAlester G, O’Gara F, Morrissey JP: Signal-mediated interactions between Pseudomonas aeruginosa and Candida albicans. J Med Microbiol. 2008, 57 (Pt 5): 563-569.PubMedGoogle Scholar
- Samonis G, Gikas A, Anaissie EJ, Vrenzos G, Maraki S, Tselentis Y, Bodey GP: Prospective evaluation of effects of broad-spectrum antibiotics on gastrointestinal yeast colonization of humans. Antimicrob Agents Chemother. 1993, 37 (1): 51-53. 10.1128/AAC.37.1.51.PubMedPubMed CentralGoogle Scholar
- Mulligan ME, Citron DM, McNamara BT, Finegold SM: Impact of cefoperazone therapy on fecal flora. Antimicrob Agents Chemother. 1982, 22 (2): 226-230. 10.1128/AAC.22.2.226.PubMedPubMed CentralGoogle Scholar
- Abbott J: Clinical and microscopic diagnosis of vaginal yeast infection: a prospective analysis. Ann Emerg Med. 1995, 25 (5): 587-591. 10.1016/S0196-0644(95)70168-0.PubMedGoogle Scholar
- Oksala E: Factors predisposing to oral yeast infections. Acta Odontol Scand. 1990, 48 (1): 71-74. 10.3109/00016359009012736.PubMedGoogle Scholar
- Leyden JJ, Marples RR: Ecologic principles and antibiotic therapy in chronic dermatoses. Arch Dermatol. 1973, 107 (2): 208-211. 10.1001/archderm.1973.01620170020006.PubMedGoogle Scholar
- Ben-Ami R, Olshtain-Pops K, Krieger M, Oren I, Bishara J, Dan M, Wiener-Well Y, Weinberger M, Zimhony O, Chowers M, Weber G, Potasman I, Chazan B, Kassis I, Shalit I, Block C, Keller N, Kontoyiannis DP, Giladi M: Antibiotic exposure as a risk factor for fluconazole-resistant Candida bloodstream infection. Antimicrob Agents Chemother. 2012, 56 (5): 2518-2523. 10.1128/AAC.05947-11.PubMedPubMed CentralGoogle Scholar
- Bow EJ, Louie TJ: Changes in endogenous microflora among febrile granulocytopenic patients receiving empiric antibiotic therapy: implications for fungal superinfection. CMAJ. 1987, 137 (5): 397-403.PubMedPubMed CentralGoogle Scholar
- Dollive S, Chen YY, Grunberg S, Bittinger K, Hoffmann C, Vandivier L, Cuff C, Lewis JD, Wu GD, Bushman FD: Fungi of the murine gut: episodic variation and proliferation during antibiotic treatment. PLoS One. 2013, 8 (8): e71806-10.1371/journal.pone.0071806.PubMedPubMed CentralGoogle Scholar
- Rogers GB, Cuthbertson L, Hoffman LR, Wing PA, Pope C, Hooftman DA, Lilley AK, Oliver A, Carroll MP, Bruce KD, van der Gast CJ: Reducing bias in bacterial community analysis of lower respiratory infections. ISME J. 2013, 7 (4): 697-706. 10.1038/ismej.2012.145.PubMedPubMed CentralGoogle Scholar
- Young VB, Raffals LH, Huse SM, Vital M, Dai D, Schloss PD, Brulc JM, Antonopoulos DA, Arrieta RL, Kwon JH, Reddy KG, Hubert NA, Grim SL, Vineis JH, Dalal S, Morrison HG, Eren AM, Meyer F, Schmidt TM, Tiedje JM, Chang EB, Sogin ML: Multiphasic analysis of the temporal development of the distal gut microbiota in patients following ileal pouch anal anastomosis. Microbiome. 2013, 1 (1): 9-10.1186/2049-2618-1-9.PubMedPubMed CentralGoogle Scholar
- Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM: Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 2007, 8 (7): R143-10.1186/gb-2007-8-7-r143.PubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27 (16): 2194-2200. 10.1093/bioinformatics/btr381.PubMedPubMed CentralGoogle Scholar
- Huse SM, Young VB, Morrison HG, Antonopoulos DA, Kwon J, Dalal S, Arrieta R, Hubert NA, Shen L, Vineis JH, Koval JC, Sogin ML, Chang EB, Raffals LE: Comparison of brush and biopsy sampling methods of the ileal pouch for assessment of mucosa-associated microbiota of human subjects. Microbiome. 2014, 2 (1): 5-10.1186/2049-2618-2-5.PubMedPubMed CentralGoogle Scholar
- Eren AM, Vineis JH, Morrison HG, Sogin ML: A filtering method to generate high quality short reads using illumina paired-end technology. PLoS One. 2013, 8 (6): e66643-10.1371/journal.pone.0066643.PubMedPubMed CentralGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41 (Database issue): D590-D596.PubMedPubMed CentralGoogle Scholar
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26 (19): 2460-2461. 10.1093/bioinformatics/btq461.PubMedGoogle Scholar
- Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R: PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010, 26 (2): 266-267. 10.1093/bioinformatics/btp636.PubMedPubMed CentralGoogle Scholar
- Jacobs J: Individual Based Rarefaction using R-package. 2011Google Scholar
- Liu CM, Kachur S, Dwan MG, Abraham AG, Aziz M, Hsueh PR, Huang YT, Busch JD, Lamit LJ, Gehring CA, Keim P, Price LB: FungiQuant: a broad-coverage fungal quantitative real-time PCR assay. BMC Microbiol. 2012, 12: 255-10.1186/1471-2180-12-255.PubMedPubMed CentralGoogle Scholar
- Horz HP, Vianna ME, Gomes BP, Conrads G: Evaluation of universal probes and primer sets for assessing total bacterial load in clinical samples: general implications and practical use in endodontic antimicrobial therapy. J Clin Microbiol. 2005, 43 (10): 5332-5337. 10.1128/JCM.43.10.5332-5337.2005.PubMedPubMed CentralGoogle Scholar
- Maeda H, Fujimoto C, Haruki Y, Maeda T, Kokeguchi S, Petelin M, Arai H, Tanimoto I, Nishimura F, Takashiba S: Quantitative real-time PCR using TaqMan and SYBR Green for Actinobacillus actinomycetemcomitans, Porphyromonas gingivalis, Prevotella intermedia, tetQ gene and total bacteria. FEMS Immunol Med Microbiol. 2003, 39 (1): 81-86. 10.1016/S0928-8244(03)00224-4.PubMedGoogle Scholar
- Goslee SC, Urban DL: The ecodist package for dissimilarity-based analysis of ecological data. J Stat Softw. 2007, 22: 1-19.Google Scholar
- Dixon P: VEGAN, a package of R functions for community ecology. J Veg Sci. 2003, 14 (6): 927-930. 10.1111/j.1654-1103.2003.tb02228.x.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.