- Research article
- Open Access
- Open Peer Review
Genetic diversity of medically important and emerging Candida species causing invasive infection
BMC Infectious Diseases volume 15, Article number: 57 (2015)
Genetic variation in the ribosomal DNA (rDNA) internal transcribed spacer (ITS) region has been studied among fungi. However, the numbers of ITS sequence polymorphisms in the various Candida species and their associations with sources of invasive fungal infections remain poorly investigated. Here, we characterized the intraspecific and interspecific ITS diversity of Candida spp. strains collected from patients with bloodstream or oroesophageal candidiasis.
We selected cultures of representative medically important species of Candida as well as some rare and emerging pathogens. Identification was performed by micromorphology and by biochemical testing using an ID32C® system, as well as by the sequencing of rDNA ITS. The presence of intraspecific ITS polymorphisms was characterized based on haplotype networks, and interspecific diversity was characterized based on Bayesian phylogenetic analysis.
Among 300 Candida strains, we identified 76 C. albicans, 14 C. dubliniensis, 40 C. tropicalis, 47 C. glabrata, 34 C. parapsilosis (sensu stricto), 31 C. orthopsilosis, 3 C. metapsilosis, 21 Meyerozyma guilliermondii (C. guilliermondii), 12 Pichia kudriavzevii (C. krusei), 6 Clavispora lusitaniae (C. lusitaniae), 3 C. intermedia, 6 Wickerhamomyces anomalus (C. pelliculosa), and 2 C. haemulonii strains, and 1 C. duobushaemulonii, 1 Kluyveromyces marxianus (C. kefyr), 1 Meyerozyma caribbica (C. fermentati), 1 Pichia norvegensis (C. norvegensis), and 1 Lodderomyces elongisporus strain. Out of a total of seven isolates with inconsistent ID32C® profiles, ITS sequencing identified one C. lusitaniae strain, three C. intermedia strains, two C. haemulonii strains and one C. duobushaemulonii strain. Analysis of ITS variability revealed a greater number of haplotypes among C. albicans, C. tropicalis, C. glabrata and C. lusitaniae, which are predominantly related to endogenous sources of acquisition. Bayesian analysis confirmed the major phylogenetic relationships among the isolates and the molecular identification of the different Candida spp.
Molecular studies based on ITS sequencing are necessary to identify closely related and emerging species. Polymorphism analysis of the ITS rDNA region demonstrated its utility as a genetic marker for species identification and phylogenetic relationships as well as for drawing inferences concerning the natural history of hematogenous infections caused by medically important and emerging Candida species.
Invasive candidiasis is recognized as a cause of morbidity and mortality in tertiary care hospitals worldwide [1,2]. Although the majority of cases of invasive yeast infection are attributed to Candida albicans, there are increasing rates of infection by non-C. albicans species in various parts of the world [2,3]. Conventional methods used by reference centers for the identification of medically important yeasts have been progressively replaced by PCR-based methods and proteomics [4,5]. However, in resource-limited settings, commercial biochemical tests still represent the cornerstone for the identification of human yeast pathogens . These methods are usually time-consuming and have potential limitations with respect to the accurate identification of rare pathogens and of cryptic species of Candida .
The establishment of DNA sequencing as a gold standard method for yeast identification by clinical laboratories has been hindered by several factors, including its cost, the lack of well-trained professionals , the limitations of the currently accepted DNA barcode system for fungi , and poor standardization of quality controls to ensure the accuracy of molecular methods . Other important issues include the quality of reference nucleotide sequences derived from well-characterized yeast collections deposited in public genomic databases and the limited number of fungal species, especially those related to human infections, for which data are present in sequence databases [10,11]. Nevertheless, progress in microorganism genomics and its application to the taxonomy of fungal pathogens has enabled an extensive review of several genera and the recognition of cryptic species within formerly recognized taxons, for example, the C. parapsilosis species complex, the C. guilliermondii complex, the C. haemulonii complex, C. rugosa and others [12-15]. Consequently, there is a need to expand public nucleotide databases to include sequences of emerging fungal pathogens [10,11].
Despite the limitations cited above, the utilization of the ribosomal DNA (rDNA) internal transcribed spacer (ITS) region for sequence analysis appears to be the most reliable strategy for the accurate and rapid molecular identification of fungal pathogens that infect humans [16,17]. Furthermore, polymorphisms in the ITS region have been extensively addressed in phylogenetic, taxonomic and population dynamics studies and are particularly useful for the delineation of Candida species and strains [16,18].
Haplotype analysis through the detection of single nucleotide polymorphisms (SNPs) found in particular target DNA sequences has been shown to be useful in the estimation of the intraspecies genetic diversity of fungal species, and this technique has also been used in population genetics studies of fungi that cause human and animal infections, permitting their evaluation from phylogenetic, biogeographic and epidemiologic perspectives [19-21].
Several molecular investigations characterizing the frequency of outbreaks and nosocomial clusters in patients with candidemia have been performed [22-25]. The methods used in these investigations have included pulsed-field gel electrophoresis (PFGE), random amplification of polymorphic DNA (RAPD), analysis of restriction fragment length polymorphisms (RFLP), PCR fingerprinting and multilocus sequence typing (MLST) [26-28]. In this context, analysis of nosocomial clustering of candidemia may be useful for determining the true frequency of exogenously acquired Candida infections transmitted to patients by the hands of caregivers and by contamination associated with invasive medical procedures.
In the present study, we aimed to determine the potential use of the rDNA ITS region as a molecular marker for evaluating genetic diversity within and among clinically important and emerging Candida species from a large Brazilian yeast collection characterized by conventional and molecular methods. The presence of intraspecific ITS variability was characterized based on haplotype networks, and Bayesian analysis was used to develop phylogenetic inferences. The incorporation of ITS sequences of human pathogenic Candida species into public nucleotide sequence databases and the reliability of ITS sequence data were also addressed.
Selection of microorganisms
For this study, 300 strains of Candida spp. were selected from the large yeast stock culture collection of the Laboratório Especial de Micologia, Escola Paulista de Medicina, Universidade Federal de São Paulo, Brazil. All fungal isolates were collected between 1997 and 2011 during multicenter surveillance studies conducted at Brazilian medical centers [1,29-31]. We selected cultures of representative medically important species of Candida, as well as cultures of some rare or emerging pathogens. With the exception of C. dubliniensis strains (n = 14), which were isolated from patients with oroesophageal infection, all species (n = 286) were obtained from blood cultures of patients with fungemia. In addition, the following 11 reference/type strains were included: C. albicans SC5314, C. albicans ATCC 24433, C. glabrata ATCC 2001, C. lusitaniae ATCC 66035, C. krusei ATCC 6258, C. tropicalis ATCC 13803, C. parapsilosis ATCC 22019, C. orthopsilosis ATCC 96141, C. metapsilosis ATCC 96143, C. guilliermondii CBS 566 and C. dubliniensis CBS 7987. The clinical isolates and the reference/type strains were identified simultaneously using phenotypic and molecular methods. In addition to the nucleotide sequences of the reference/type strains obtained and described above, the sequences of the following eight strains were used for haplotype and phylogenetic analyses: Lodderomyces elongisporus ATCC 11503 (GenBank accession number: NR_111593.1), C. kefyr ATCC 60480 (GenBank accession number: GU256755.1), C. pelliculosa CBS 606 (MycoBank accession number: 346023), C. fermentati CBS 2022 (GenBank accession number: EU568913.1), C. norvegensis CBS 2128 (GenBank accession number: AB278167.1), C. intermedia WM 811 (GenBank accession number: EF568011.1), C. haemulonii CBS 10970 (GenBank accession number: JX459674.1) and C. duobushaemulonii CBS 7798 (GenBank accession number: JX459666.1). This work was approved by the institutional board on ethics in research of the Universidade Federal de São Paulo, Brazil (CEP008/11).
Conventional identification of Candida species
Cultures of Candida spp. stored at -80°C were plated on CHROMagarTMCandida (CHROMagar Microbiology, Paris, France) prepared according to the manufacturer’s instructions for 48 h at 37°C to obtain pure colonies and for presumptive identification of Candida spp. Slide cultures on cornmeal agar medium with 1.2% Tween 80 were prepared to evaluate the presence of chlamydoconidia, blastoconidia and pseudohyphae. The ability to grow on Difco™ Sabouraud Dextrose Agar (SDA) (Becton Dickinson & Co. Sparks, MD, USA) plates at 42°C after 48 h of culture  or on hypertonic Sabouraud broth at 37°C for 96 h was used to discriminate between C. albicans and C. dubliniensis [32,33]. Reference strains of C. albicans (ATCC 24433) and C. dubliniensis (CBS 7987) were used as controls. With the exception of C. albicans and C. dubliniensis, the biochemical profiles of all Candida species were evaluated using a commercial ID32C® system according to the manufacturer’s instructions (bioMérieux, Marcy-l’Étoile, France).
Molecular identification of Candida species by sequencing of rDNA ITS
Total genomic DNA was extracted from the Candida isolates using PrepMan® Ultra Sample Preparation Reagent (Applied Biosystems, Inc., Foster City, CA, USA) according to the manufacturer’s instructions. PCR for the amplification of the ITS region was performed using the forward primer V9G (5′-TTACGTCCCTGCCCTTTGTA-3′) and the reverse primer LS266 (5′-GCATTCCCAAACAACTCGACTC-3′) . The total length of the amplified product was approximately 924 base pairs for C. albicans. A total reaction volume of 25 μl containing 40 ng/ml of genomic DNA, 10 pmol/μl of each primer, and PCR Master Mix with 50 units/ml of Taq DNA polymerase, 3 mM MgCl2, 400 μM dNTPs (Promega, Madison, WI, USA), and sterile water were used for the PCR reactions, which were performed in a Veriti 96-well Thermal Cycler (Applied Biosystems, Inc., Foster City, CA, USA) under the following conditions: an initial denaturation step at 94°C for 5 min, 35 cycles of denaturation at 94°C for 1 min, annealing at 56°C for 30 s, and extension at 72°C for 2 min, and a final extension step at 72°C for 10 min. Positive (DNA from C. albicans ATCC 24433) and negative (sample lacking DNA) controls were included in all assays. The amplicons were verified by electrophoresis at 90 volts in 1% agarose/SYBR® Safe DNA Stain (Invitrogen, Carlsbad, CA, USA) gels and photographed using a UV transilluminator.
PCR products were subjected to dideoxynucleotide sequencing with a Big Dye Terminator Reaction Kit v3.1 (Applied Biosystems, Inc., Foster City, CA, USA), using the forward primers V9G and ITS1 (5′-TCCGTAGGTGAACCTGCGG-3′) and the reverse primers LS266 and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) [34,35] according to the manufacturer’s instructions. After purification and denaturation, the samples were run on an automated ABI 3130 genetic analyzer (Applied Biosystems, Inc., Foster City, CA, USA). For the sequencing reactions, a total of six sequences were used, including three forward strands (one strand sequenced with V9G and two with ITS1) and three reverse strands (one strand sequenced with LS266 and two with ITS4) for each strain to increase confidence in the sequencing data for the detection of nucleotide polymorphisms and to avoid experimental artifacts.
Consensus sequence assembly and editing were performed using the programs Phred/Phrap and the sequence editor Consed [36-38]. The error probability for each called base was assessed, considering a Phred score > 40, which was associated with a base call accuracy of 99.99%. High-quality consensus sequences were obtained for analysis with assembly errors of less than one per 100 base pairs after editing. The consensus sequences obtained in our study were aligned and compared with sequences deposited in public genomic databases (GenBank, NCBI, USA and CBS database, the Netherlands). To ensure the high accuracy of the results obtained using the nucleotide sequence alignment tools, an e-value of less than 10-5 and a maximum identity of equal to or higher than 98% were considered for the correct identification of Candida at the species level.
Analysis of intraspecific and interspecific diversities of ITS sequences of Candida spp.
Haplotype analysis was conducted to assess the ITS intraspecific variability of the Candida species. In the present study, a haplotype was defined as a unique combination of SNPs along a sequence, i.e., each different sequence in an alignment. The ITS sequences of the Candida spp. were aligned and edited using the muscle algorithm implemented in SEAVIEW program 4.2.12 , excluding 18S and 28S rDNA. The complete ITS sequences, including ITS1 and ITS2 and the 5.8S region of rDNA were subjected to sequence polymorphism analysis using DnaSP version 5.10 software . The analysis was based on the number of haplotypes (Hap), variable sites, haplotype diversity (Hd), and nucleotide diversity (Pi) . In brief, haplotype diversity is a measure of the occurrence of a single haplotype in a given species, considering the number of sequences analyzed and the total number of haplotypes found. Values range from 0 to 1, with those closer to 1 indicating higher variability. Nucleotide diversity is a measure of the average number of nucleotide differences per site between two sequences. Haplotype 1 corresponded with the reference or type strain of each Candida spp.
For the ITS haplotype network, a total of 319 DNA sequences of clinical (n = 300) and reference/type strains (n = 19) comprising 17 Candida species and 1 non-Candida species (Lodderomyces elongisporus) were aligned. The haplotype network file (Roehl data file) was created using DnaSP v.5.10 software, considering gaps. The network was generated by the median-joining method  using Network v4.612 software (http://www.fluxus-engineering.com/).
To analyze the conservation of phylogenetic relationships among the different Candida species, sequences representative of each haplotype (Additional file 1: Table S1) and the reference sequence found in the comparisons with the public genomic databases were aligned and edited using the muscle algorithm implemented in SEAVIEW program 4.2.12  and considering only the ITS1-ITS2 and 5.8S-rDNA regions. Phylogenetic inference was performed using MrBayes 3.02  with the default priors as input and was run twice with four chains. The number of generations was 2 million, and data were saved every 100 generations along with the GTR (general time reversible) model, the shape of the gamma distribution parameters, and the proportions of invariant sites estimated during the run. Bootstrap analysis was conducted by evaluating 1,000 pseudoreplicates of the alignment with SEAVIEW program 4.2.12  using the neighbor-joining method.
Nucleotide sequence accession numbers
The rDNA ITS sequences of the clinical strains of Candida spp. and Lodderomyces elongisporus were deposited in the GenBank database with the accession numbers KC408939 to KC408999. For complete information including strain name, species name and accession numbers, see Additional file 1: Table S1.
In the present study, one hundred percent concordance between the phenotypic and molecular methods was obtained for a collection of clinically important yeast isolates containing representative numbers of strains (≥ 40 strains) as follows: C. albicans (n = 76), C. tropicalis (n = 40) and C. glabrata (n = 47) (Additional file 2: Table S2). Growth tests in hypertonic broth and at 42°C allowed for the presumptive identification of 14 C. dubliniensis and 76 C. albicans isolates, and ITS sequencing confirmed the identification of these two species, indicating 100% concordance (Additional file 2: Table S2). As expected, among 69 isolates phenotypically identified as C. parapsilosis (sensu lato), ITS sequencing distinguished 34 C. parapsilosis (sensu stricto), 31 C. orthopsilosis and three C. metapsilosis isolates and one Lodderomyces elongisporus isolate. A total of 22 isolates were identified as C. guilliermondii by biochemical tests, while DNA sequence analysis detected 21 Meyerozyma guilliermondii (teleomorph of C. guilliermondii) isolates and one Meyerozyma caribbica (teleomorph of C. fermentati) isolate.
Seven isolates that exhibited inconclusive results based on the ID32C® system were identified by ITS sequencing as C. lusitaniae (n = 1), C. intermedia (n = 3), C. haemulonii (n = 2) and C. duobushaemulonii (n = 1). Of six isolates identified as C. lusitaniae by molecular methods, one showed an inconsistent profile with the biochemical test. ITS sequencing distinguished species within the C. haemulonii complex, including two C. haemulonii isolates and one C. duobushaemulonii isolate (see Additional file 2: Table S2).
The presence of sequence polymorphisms along the total fragment length of the ITS rDNA (ITS1 and ITS2, including the 5.8S region) was analyzed in clinically important Candida species with representative numbers of strains (> 20 strains), as illustrated in Figure 1. A comparison of the ITS sequences of the clinical strains to the sequence of the reference or type strain of each Candida species revealed that nucleotide sequence variability per site was increased in ITS1 compared with ITS2. Five polymorphic sites in ITS1 and three in ITS2 were identified for C. albicans in addition to five variations in ITS1 and five in ITS2 for C. tropicalis, 13 in ITS1 and seven in ITS2 for C. glabrata, one in ITS1 for C. parapsilosis (sensu stricto), six in ITS1 and one in ITS2 for C. orthopsilosis, and one in ITS2 for M. guilliermondii. Sequence variations in the 5.8S region were detected only for C. albicans (one site) and C. tropicalis (one site).
Based on the results of rDNA ITS sequence analysis, we also determined the presence of intraspecific variability in the Candida species, of which three or more strains were analyzed and compared to type/reference strains (Table 1). High intraspecific variation was found for C. albicans, C. tropicalis, C. glabrata, C. metapsilosis, P. kudriavzevii, C. lusitaniae and C. intermedia, as evidenced by the large numbers of haplotypes and variable sites and haplotype (Hd = 0.6179 to 0.8571) and nucleotide (Pi = 0.00164 to 0.02815) diversity. It is worth mentioning that C. lusitaniae showed the highest haplotype diversity (Hd = 0.8571) despite the low number of isolates of this species (n = 6). The measurement of nucleotide diversity also revealed more nucleotide differences per site in C. intermedia (Pi = 0.02815) and C. lusitaniae (Pi = 0.02605) than in the other species tested.
On the other hand, minor intraspecific variations were observed among isolates of C. parapsilosis (sensu stricto) (Hd = 0.2924 and Pi = 0.00068), C. orthopsilosis (Hd = 0.4456 and Pi = 0.00083), and M. guilliermondii (Hd = 0.4848 and Pi = 0.00094). Another interesting finding was the presence of only one ITS haplotype for each of Wickerhamomyces anomalus (n = 7) and C. haemulonii (n = 3), as demonstrated by the lack of polymorphic sites and haplotype and nucleotide diversities in these species.
The ITS haplotype network constructed by the alignment of 319 sequences (317 sequences of Candida spp. and 2 L. elongisporus sequences) showed the presence of 67 haplotypes (Figure 2A). At a total of 873 sites analyzed, we found 716 variable positions (Hd = 0.9556), considering sites with gaps. Notable genetic diversity was observed, as demonstrated by the large numbers of haplotypes in clinically relevant species and in some rare species of Candida, including 12 C. albicans (n = 78), 9 C. tropicalis (n = 41), 11 C. glabrata (n = 48), 5 P. kudriavzevii (n = 13), 5 C. lusitaniae (n = 7) and 3 C. intermedia haplotypes (n = 4). In contrast, only 2 haplotypes each were found for C. parapsilosis (sensu stricto) (n = 35), C. dubliniensis (n = 15) and M. guilliermondii (n = 22), and only 1 haplotype each was found for W. anomalus (n = 7) and C. haemulonii (n = 3). With respect to acquisition sources of Candida infection (Figure 2B), high ITS diversity was found among species predominantly associated with endogenous sources of infection (C. albicans, C. tropicalis, C. glabrata and C. lusitaniae), while low genetic diversity was observed in species predominantly related to exogenous infection sources (C. parapsilosis species complex, M. guilliermondii and W. anomalus).
Bayesian analysis (Figure 3) confirmed the molecular identification of and the major phylogenetic relationships among the various Candida spp. isolates, with high posterior probabilities and bootstrap values in the branches. Intraspecies variation was also observed within the clades, demonstrating the conservation of the haplotypes described herein.
In our study, large numbers of ITS sequences were generated based on quality control criteria designed to ensure the accuracy of the molecular data used to identify the Candida isolates at the species level. The quality control criteria included PCR controls, the use of triplicate reactions with both DNA strands for sequencing, the estimation of assembly errors for high-quality consensus sequences, and BLAST search parameters. We deposited 60 ITS sequences of the most important fungal pathogens and rare or emerging pathogens, including 17 Candida species and 1 non-Candida species (Lodderomyces elongisporus), in the NCBI database (GenBank, USA). Fifty-eight ITS sequences, comprising the sequences of all Candida species and L. elongisporus submitted to GenBank, were deposited in the ISHAM-ITS reference database, which can be accessed at http://www.isham.org or at http://its.mycologylab.org.
The accurate identification of Candida spp. is essential for selecting the most effective therapeutic strategies to control invasive fungal infections caused by these species . Here, the phenotypic methods used performed well in identifying clinically important Candida species, including C. albicans, C. tropicalis, the C. parapsilosis complex and C. glabrata. However, the ID32C® system failed to identify one C. lusitaniae isolate, three C. intermedia isolates and three isolates belonging to the C. haemulonii species complex. Some limitations of conventional methods have been described for discrimination between C. lusitaniae and the closely related species C. pulcherrima . Conventional identification of the C. haemulonii species complex is still limited, however, because it is not included in the database of current commercial biochemical systems, such as ID32C® . Nevertheless, the reliable identification of rare and emerging Candida species that cause hematogenous infection, such as C. lusitaniae, C. intermedia and the C. haemulonii complex, was possible using DNA sequencing, corroborating with the results presented by other authors [6,44,45].
For epidemiological or diagnostic purposes, DNA-based methods have been used to differentiate the species forming the C. parapsilosis complex . Out of a total of 69 isolates phenotypically identified as C. parapsilosis (sensu lato), 34 were determined to be C. parapsilosis (sensu stricto), 31 were found to be C. orthopsilosis, three were found to be C. metapsilosis and one was determined to be Lodderomyces elongisporus by ITS sequencing. Lodderomyces elongisporus was initially described as a teleomorph of C. parapsilosis . Although L. elongisporus is considered to be an uncommon cause of human disease, isolates have been identified by DNA sequence analysis of some clinical sources [6,48]. M. guilliermondii is also currently considered to comprise a complex formed by Debaryomyces hansenii (teleomorph of C. famata), M. caribbica, C. carpophila and C. xestobii , which can be differentiated by molecular methods, including sequencing of the ITS and D1/D2 rDNA regions . Here, we found one isolate of M. caribbica by sequence analysis, whereas only M. guilliermondii isolates were detected by biochemical testing.
We also investigated intraspecific rDNA ITS polymorphisms of isolates from 17 Candida species and 1 non-Candida species (Lodderomyces elongisporus). Some authors have addressed the degree of variability within the rDNA ITS region among fungi and have also observed higher intraspecific diversity within the ITS1 region than within the ITS2 region [50,51]. Based on haplotype and network analysis, a high intraspecific variability of ITS sequences was observed in species that primarily represent endogenous sources of infection, including C. albicans, C. tropicalis, C. glabrata and C. lusitaniae . Species with high genetic diversity are most frequently human commensals, and this finding could explain the existence of additional genetic adaptation within normal microbiota with older evolutionary origins. Although the natural mode of reproduction of C. albicans is known to be clonal, other mechanisms that increase the genetic variability of this species could occur, including recombination [52,53]. Pfaller et al.  have tested 47 C. lusitaniae isolates, obtaining 28 different karyotype profiles and 25 different types of restriction endonuclease analysis of genomic DNA (REAG) profiles. Our data confirm the great diversity among C. lusitaniae haplotypes, possibly indicating the existence of a non-clonal form of propagation. In contrast, less intraspecific variation was observed in the isolates of C. parapsilosis and C. orthopsilosis, in which the primary mode of infection is thought to be exogenous . Previous studies have shown a lower sequence variability of C. parapsilosis (sensu stricto) compared with C. orthopsilosis and C. metapsilosis isolates [56,57]. Furthermore, we detected no ITS polymorphisms in six W. anomalus isolates. Barchiesi and colleagues  have studied 46 clinical isolates of W. anomalus using RFLP and RAPD and have found that all of these isolates produce similar band patterns. The lack of ITS variability of W. anomalus strains could be explained by the small number of sequenced isolates (n = 6), by possible clonal origin or even by a low mutation rate of the gene selected in our present study. In contrast with species that are more prevalent as human commensal organisms, other Candida species, such as the C. parapsilosis complex, M. guilliermondii and W. anomalus, have been reported to be exogenous organisms that have been isolated from environmental sources (plants, soil, insects, and food), medical devices (central venous catheter and parenteral nutrition), or the hands of health care workers [2,55,59,60]. Although they are considered to be rare causative agents of fungemia, W. anomalus and M. guilliermondii have been associated with the occurrence of nosocomial outbreaks and pseudo-outbreaks, respectively, especially in pediatric intensive care units [61,62]. Species with low genetic variability that are more often associated with exogenous transmission could be less well adapted to human hosts, may predominantly undergo clonal reproduction, and might have recently diverged during their evolutionary histories [12,58,63].
In the present study, despite the limitation conferred by the use of only one genetic marker to estimate intraspecific diversity, it was possible to use the ITS haplotype network to identify remarkable genetic polymorphisms among clinical isolates comprising 17 Candida species. Of course, the accuracy of haplotype analysis is influenced by the molecular markers selected as well as by the number of strains in the database and the network methods used . Importantly, all of the ITS sequences used for the evaluation of intraspecific genetic variation were generated by the direct sequencing of PCR products. The provision of reliable ITS sequences was based on quality control criteria that were applied from PCR to sequence analysis, for example, the amplification of products with a total length that included the entire ITS region, the use of a cut-off number of sequences to generate the consensus sequence for each strain, and the use of well-established computational tools for sequence assembly and editing.
Recent molecular epidemiological surveys aiming to assess genetic relationships among strains and temporal and geographic distributions of clustered isolates as well as to identify outbreaks and the recurrence of BSIs investigated the occurrence of the nosocomial clustering of candidemia caused by the most prevalent species of Candida [23,25]. Using PCR fingerprinting, a population-based study conducted in Iceland has shown that 18.7% to 39.9% of all cases of candidemia are nosocomial clusters primarily caused by C. albicans, C. tropicalis and C. parapsilosis . In a recent epidemiologic study conducted by Maganti and colleagues , nosocomial clusters have been shown to represent 33% of total isolates causing candidemia in Canadian hospitals. Using MLST analysis, the genetic relatedness of Candida isolates has been assessed in two Brazilian studies of candidemia, in which different clusters have been found among isolates of C. albicans  and C. tropicalis . Our present study reinforces the utility of sequence polymorphism analysis for evaluating the relationships among Candida species and its application for examining the molecular epidemiology of fungal diseases.
Comparative analysis of our Candida ITS sequences with those deposited in public genomic databases allowed for the definition of appropriate quality parameters for use with the sequence alignment search tools of the NCBI and CBS databases, as demonstrated by the high identities of our sequences at the species level to the sequences in the nucleotide repositories. However, some caution must be exercised with respect to species identification using DNA sequencing because several factors may influence the accuracy of this molecular assay, including the methodologies and programs used for sequence analysis and the limited number of representative pathogenic fungal species deposited in genomic databases, which may include incomplete sequences as well as errors in species nomenclature. Analysis of the ITS sequence regions of different taxonomic groups revealed that up to 20% of sequences deposited in GenBank may contain errors in species identification and/or outdated nomenclature and may lack descriptive and updated annotations .
In conclusion, although conventional methods can be used to reliably identify the most common Candida species, molecular studies based on ITS sequencing are necessary for the identification of closely related and emerging species. DNA sequence polymorphisms, especially those observed in ITS1, were more likely to be found in Candida species primarily involved in endogenously acquired infections. Thus, C. albicans, C. tropicalis, C. glabrata and C. lusitaniae showed high genetic diversity, while species predominantly associated with exogenously acquired infections, such as C. parapsilosis, M. guilliermondii and W. anomalus, showed low intraspecific variability. In addition, our findings indicate the importance of generating accurate ITS sequence regions based on quantitative parameters, such as those used in this study. The use of this criterion may help to increase the number of good-quality DNA sequences of Candida species deposited in public genomic databases, ensuring reliable information for future studies involving the epidemiological, clinical and molecular characterization of opportunistic fungi.
American Type Culture Collection
Nucleotide Basic Local Alignment Search Tool
Centraalbureau voor Schimmelcultures
Internal transcribed spacer
National Center for Biotechnology Information
Polymerase chain reaction
Random amplification of polymorphic DNA
Restriction fragment length polymorphisms
Multilocus sequence typing
Pulsed-field gel electrophoresis
Restriction endonuclease analysis of genomic DNA
- MALDI-TOF MS:
Matrix-assisted laser desorption/ionization-time of flight
General time reversible
Colombo AL, Nucci M, Park BJ, Nouér SA, Arthington-Skaggs B, da Matta DA, et al. Epidemiology of candidemia in Brazil: a nationwide sentinel surveillance of candidemia in eleven medical centers. J Clin Microbiol. 2006;44:2816–23.
Pfaller MA, Diekema DJ. Epidemiology of invasive candidiasis: a persistent public health problem. Clin Microbiol Rev. 2007;20:133–63.
Nishikaku AS, Melo ASA, Colombo AL. Geographic trends in invasive candidiasis. Curr Fungal Infect Rep. 2010;4:210–8.
Pincus DH, Orenga S, Chatellier S. Yeast identification–past, present, and future methods. Med Mycol. 2007;45:97–121.
Posteraro B, De Carolis E, Vella A, Sanguinetti M. MALDI-TOF mass spectrometry in the clinical mycology laboratory: identification of fungi and beyond. Expert Rev Proteomics. 2013;10:151–64.
Cendejas-Bueno E, Gomez-Lopez A, Mellado E, Rodriguez-Tudela JL, Cuenca-Estrella M. Identification of pathogenic rare yeast species in clinical samples: comparison between phenotypical and molecular methods. J Clin Microbiol. 2010;48:1895–9.
Ellepola ANB, Morrison CJ. Laboratory diagnosis of invasive candidiasis. J Microbiol. 2005;43 Spec:65–84.
Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for fungi. Proc Natl Acad Sci U S A. 2012;109:6241–6.
Clinical and Laboratory Standards Institute. Interpretive criteria for identification of bacteria and fungi by DNA target sequencing: approved guideline. CLSI document MM18-A. Wayne, PA: Clinical and Laboratory Standards Institute; 2008.
Bridge PD, Roberts PJ, Spooner BM, Panchal G. On the unreliability of published DNA sequences. New Phytol. 2003;160:43–8.
Balajee SA, Borman AM, Brandt ME, Cano J, Cuenca-Estrella M, Dannaoui E, et al. Sequence-based identification of Aspergillus, Fusarium, and Mucorales species in the clinical mycology laboratory: where are we and where should we go from here? J Clin Microbiol. 2009;47:877–84.
Tavanti A, Davidson AD, Gow NAR, Maiden MCJ, Odds FC. Candida orthopsilosis and Candida metapsilosis spp. nov. to replace Candida parapsilosis groups II and III. J Clin Microbiol. 2005;43:284–92.
Vaughan-Martini A, Kurtzman CP, Meyer SA, O’Neill EB. Two new species in the Pichia guilliermondii clade: Pichia caribbica sp. nov., the ascosporic state of Candida fermentati, and Candida carpophila comb. nov. FEMS Yeast Res. 2005;5:463–9.
Cendejas-Bueno E, Kolecka A, Alastruey-Izquierdo A, Theelen B, Groenewald M, Kostrzewa M, et al. Reclassification of the Candida haemulonii complex as Candida haemulonii (C. haemulonii group I), C. duobushaemulonii sp. nov. (C. haemulonii group II), and C. haemulonii var. vulnera var. nov.: three multiresistant human pathogenic yeasts. J Clin Microbiol. 2012;50:3641–5361.
Chaves GM, Terçarioli GR, Padovan ACB, Rosas RC, Ferreira RC, Melo ASA, et al. Candida mesorugosa sp. nov., a novel yeast species similar to Candida rugosa, isolated from a tertiary hospital in Brazil. Med Mycol. 2013;51:231–42.
Iwen PC, Hinrichs SH, Rupp ME. Utilization of the internal transcribed spacer regions as molecular targets to detect and identify human fungal pathogens. Med Mycol. 2002;40:87–109.
Ciardo DE, Lucke K, Imhof A, Bloemberg GV, Böttger EC. Systematic internal transcribed spacer sequence analysis for identification of clinical mold isolates in diagnostic mycology: a 5-year study. J Clin Microbiol. 2010;48:2809–13.
Chen YC, Eisner JD, Kattar MM, Rassoulian-Barrett SL, Lafe K, Yarfitz SL, et al. Identification of medically important yeasts using PCR-based detection of DNA sequence polymorphisms in the internal transcribed spacer 2 region of the rRNA genes. J Clin Microbiol. 2000;38:2302–10.
Bartelli TF, Ferreira RC, Colombo AL, Briones MRS. Intraspecific comparative genomics of Candida albicans mitochondria reveals non-coding regions under neutral evolution. Infect Genet Evol. 2013;14:302–12.
Cogliati M, Zamfirova RR, Tortorano AM, Viviani MA. Molecular epidemiology of Italian clinical Cryptococcus neoformans var. grubii isolates. Med Mycol. 2013;51:499–506.
Rodrigues AM, de Hoog G, Zhang Y, de Camargo ZP. Emerging sporotrichosis is driven by clonal and recombinant Sporothrix species. Emerg Microbes Infect. 2014;3:e32.
Clark TA, Slavinski SA, Morgan J, Lott T, Arthington-Skaggs BA, Brandt ME, et al. Epidemiologic and molecular characterization of an outbreak of Candida parapsilosis bloodstream infections in a community hospital. J Clin Microbiol. 2004;42:4468–72.
Asmundsdóttir LR, Erlendsdóttir H, Haraldsson G, Guo H, Xu J, Gottfredsson M. Molecular epidemiology of candidemia: evidence of clusters of smoldering nosocomial infections. Clin Infect Dis. 2008;47:e17–24.
Da Matta DA, Melo AS, Colombo AL, Frade JP, Nucci M, Lott TJ. Candidemia surveillance in Brazil: evidence for a geographical boundary defining an area exhibiting an abatement of infections by Candida albicans group 2 strains. J Clin Microbiol. 2010;48:3062–7.
Maganti H, Yamamura D, Xu J. Prevalent nosocomial clusters among causative agents for candidemia in Hamilton, Canada. Med Mycol. 2011;49:530–8.
Robles JC, Koreen L, Park S, Perlin DS. Multilocus sequence typing is a reliable alternative method to DNA fingerprinting for discrimination among strains of Candida albicans. J Clin Microbiol. 2004;42:2480–8.
Lockhart SR, Pujol C, Dodgson AR, Soll DR. Deoxyribonucleic acid fingerprinting methods for Candida species. Methods Mol Med. 2005;118:15–25.
Pujol C, Soll DR. DNA fingerprinting Candida species. Methods Mol Biol. 2009;499:117–29.
da Matta DA, de Almeida LP, Machado AM, Azevedo AC, Kusano EJU, Travassos NF, et al. Antifungal susceptibility of 1000 Candida bloodstream isolates to 5 antifungal drugs: results of a multicenter study conducted in São Paulo, Brazil, 1995-2003. Diagn Microbiol Infect Dis. 2007;57:399–404.
Nucci M, Queiroz-Telles F, Alvarado-Matute T, Tiraboschi IN, Cortes J, Zurita J, et al. Epidemiology of candidemia in Latin America: a laboratory-based survey. PLoS One. 2013;8:e59373.
Milan EP, de Laet Sant’ Ana P, de Azevedo Melo AS, Sullivan DJ, Coleman DC, Lewi D, et al. Multicenter prospective surveillance of oral Candida dubliniensis among adult Brazilian human immunodeficiency virus-positive and AIDS patients. Diagn Microbiol Infect Dis. 2001;41:29–35.
Sullivan DJ, Westerneng TJ, Haynes KA, Bennett DE, Coleman DC. Candida dubliniensis sp. nov.: phenotypic and molecular characterization of a novel species associated with oral candidosis in HIV-infected individuals. Microbiology. 1995;141:1507–21.
Alves SH, Pipolo Milan E, De Laet Sant’Ana P, Oliveira LO, Santurio JM, Colombo AL. Hypertonic sabouraud broth as a simple and powerful test for Candida dubliniensis screening. Diagn Microbiol Infect Dis. 2002;43:85–6.
van den Ende AHG G, de Hoog GS. Variability and molecular diagnostics of the neurotropic species Cladophialophora bantiana. Stud Mycol. 1999;43:151–62.
White TJ, Bruns S, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR protocols: a guide to methods applications. New York: Academic Press Inc; 1990. p. 315–22.
Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998;8:186–94.
Ewing B, Hillier L, Wendl M, Green P. Base-calling of automated sequencer traces using Phred. I. Accuracy assessment. Genome Res. 1998;8:175–85.
Gordon D, Abajian C, Green P. Consed: a graphical tool for sequence finishing. Genome Res. 1998;8:195–202.
Gouy M, Guindon S, Gascuel O. SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221–4.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Nei M. Molecular evolutionary genetics. New York: Columbia University Press; 1987.
Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.
Ronquist F, Huelsenbeck JP. MrBayes 3: bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Noël T, Favel A, Michel-Nguyen A, Goumar A, Fallague K, Chastin C, et al. Differentiation between atypical isolates of Candida lusitaniae and Candida pulcherrima by determination of mating type. J Clin Microbiol. 2005;43:1430–2.
Ruan S-Y, Chien J-Y, Hou Y-C, Hsueh P-R. Catheter-related fungemia caused by Candida intermedia. Int J Infect Dis. 2010;14:e147–9.
Souza ACR, Ferreira RC, Goncalves SS, Quindos G, Eraso E, Bizerra FC, et al. Accurate identification of Candida parapsilosis (sensu lato) by use of mitochondrial DNA and Real-Time PCR. J Clin Microbiol. 2012;50:2310–4.
Hamajima K, Nishikawa A, Shinoda T, Fukazawa Y. Deoxyribonucleic acid base composition and its homology between two forms of Candida parapsilosis and Lodderomyces elongisporus. J Gen Appl Microbiol. 1987;33:299–302.
Lockhart SR, Messer SA, Pfaller MA, Diekema DJ. 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. 2008;46:2659–64.
Desnos-Ollivier M, Ragon M, Robert V, Raoux D, Gantier J-C, Dromer F. Debaryomyces hansenii (Candida famata), a rare human fungal pathogen often misidentified as Pichia guilliermondii (Candida guilliermondii). J Clin Microbiol. 2008;46:3237–42.
Leaw SN, Chang HC, Sun HF, Barton R, Bouchara J-P, Chang TC. Identification of medically important yeast species by sequence analysis of the internal transcribed spacer regions. J Clin Microbiol. 2006;44:693–9.
Nilsson RH, Kristiansson E, Ryberg M, Hallenberg N, Larsson K-H. Intraspecific ITS variability in the kingdom fungi as expressed in the international sequence databases and its implications for molecular species identification. Evol Bioinform. 2008;4:193–201.
Tavanti A, Gow NAR, Maiden MCJ, Odds FC, Shaw DJ. Genetic evidence for recombination in Candida albicans based on haplotype analysis. Fungal Genet Biol. 2004;41:553–62.
Odds FC, Jacobsen MD. Multilocus sequence typing of pathogenic Candida species. Eukaryot Cell. 2008;7:1075–84.
Pfaller MA, Messer SA, Hollis RJ. Strain delineation and antifungal susceptibilities of epidemiologically related and unrelated isolates of Candida lusitaniae. Diagn Microbiol Infect Dis. 1994;20:127–33.
Trofa D, Gácser A, Nosanchuk JD. Candida parapsilosis, an emerging fungal pathogen. Clin Microbiol Rev. 2008;21:606–25.
Tavanti A, Hensgens LAM, Ghelardi E, Campa M, Senesi S. Genotyping of Candida orthopsilosis clinical isolates by amplification fragment length polymorphism reveals genetic diversity among independent isolates and strain maintenance within patients. J Clin Microbiol. 2007;45:1455–62.
Tavanti A, Hensgens LAM, Mogavero S, Majoros L, Senesi S, Campa M. Genotypic and phenotypic properties of Candida parapsilosis sensu strictu strains isolated from different geographic regions and body sites. BMC Microbiol. 2010;10:203.
Barchiesi F, Tortorano AM, Di Francesco LF, Rigoni A, Giacometti A, Spreghini E, et al. Genotypic variation and antifungal susceptibilities of Candida pelliculosa clinical isolates. J Med Microbiol. 2005;54:279–85.
Calderone RA. Candida and Candidiasis. Washington, DC: ASM Press; 1998.
Kam AP, Xu J. Diversity of commensal yeasts within and among healthy hosts. Diagn Microbiol Infect Dis. 2002;43:19–28.
Pasqualotto AC, Sukiennik TCT, Severo LC, de Amorim CS, Colombo AL. An outbreak of Pichia anomala fungemia in a Brazilian pediatric intensive care unit. Infect Control Hosp Epidemiol. 2005;26:553–8.
Medeiros EAS, Lott TJ, Colombo AL, Godoy P, Coutinho AP, Braga MS, et al. Evidence for a pseudo-outbreak of Candida guilliermondii fungemia in a university hospital in Brazil. J Clin Microbiol. 2007;45:942–7.
Lan L, Xu J. Multiple gene genealogical analyses suggest divergence and recent clonal dispersal in the opportunistic human pathogen Candida guilliermondii. Microbiology. 2006;152:1539–49.
Woolley SM, Posada D, Crandall KA. A comparison of phylogenetic network methods using computer simulation. PLoS One. 2008;3:e1913.
Magri MMC, Gomes-Gouvêa MS, de Freitas VLT, Motta AL, Moretti ML, Shikanai-Yasuda MA. Multilocus sequence typing of Candida tropicalis shows the presence of different clonal clusters and fluconazole susceptibility profiles in sequential isolates from candidemia patients in Sao Paulo, Brazil. J Clin Microbiol. 2013;51:268–77.
We are grateful to Dr. Paulo B. Paiva for the bioinformatics support for the sequence assembly and editing. We also thank Ricardo A. Siqueira for the help with the phenotypic identification. This study was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo - FAPESP (grant # 2007/08575-1) and Conselho Nacional de Desenvolvimento Científico e Tecnológico – CNPq.
The authors declare that they have no competing interest. The authors alone are responsible for the content and writing of this paper.
KBM performed the phenotypic and molecular identifications and haplotype analysis. ASN participated in the molecular identification and sequence and haplotype analyses. AMR performed haplotype network analysis. ACP assisted with sequence analysis and phylogenetic inferences. RCF contributed to the sequence assembly and editing and to haplotype analysis. ASAM contributed to the study design and revision of the manuscript. MRSB assisted in the establishment of the sequence analysis criteria. ALC designed the research and supervised the study. KBM, ASN, AMR, ACP, ASAM and ALC wrote the manuscript. All authors approved and contributed to the final manuscript.
Karina Bellinghausen Merseguel and Angela Satie Nishikaku contributed equally to this work.
Ribosomal DNA ITS sequences of Candida spp. and Lodderomyces elongisporus deposited in public sequence databases. Species names, strains, haplotypes and GenBank accession numbers.
Comparative analysis of the phenotypic and molecular methods used to identify Candida species. The yeast strains were isolated from bloodstream (n = 286) and oroesophageal (n = 14) infections.
About this article
Cite this article
Merseguel, K.B., Nishikaku, A.S., Rodrigues, A.M. et al. Genetic diversity of medically important and emerging Candida species causing invasive infection. BMC Infect Dis 15, 57 (2015) doi:10.1186/s12879-015-0793-3
- Candida spp.
- DNA sequencing
- Internal transcribed spacer region
- Haplotype network
- Nosocomial acquisition of candidemia