Previous Article | Next Article 
Journal of Clinical Microbiology, November 2007, p. 3781-3784, Vol. 45, No. 11
0095-1137/07/$08.00+0 doi:10.1128/JCM.01603-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Multiple-Locus Variable-Number Tandem-Repeat Analysis for Rapid Typing of Candida glabrata
Frédéric Grenouillet,1,2*
Laurence Millon,1,2
Jean-Mathieu Bart,2
Sandrine Roussel,2
Isabelle Biot,1
Emeline Didier,1
Anne-Sophie Ong,1 and
Renaud Piarroux1,2
Mycology-Parasitology Department, University Hospital Jean Minjoz,1
SERF (Santé Environnement Rural) Team, EA 2276, Faculty of Medicine, University of Franche-Comté, 25030 Besançon, France2
Received 11 August 2007/
Returned for modification 27 August 2007/
Accepted 31 August 2007
 |
ABSTRACT
|
|---|
A multiple-locus variable-number tandem-repeat analysis (MLVA) using six microsatellite markers was assessed in 127 Candida glabrata isolates. Thirty-seven different genotypes, stable both in vitro and in vivo, were observed. The highest discriminatory power (D = 0.902) was reached by using only four markers. MLVA seems to be relevant for C. glabrata typing.
 |
TEXT
|
|---|
Candida glabrata has recently emerged as a major pathogen, causing mucosal and systemic infections (10, 17). Several methods, such as electrophoretic karyotyping, restriction enzyme analysis, Southern blotting with probes, randomly amplified polymorphic DNA, and multilocus sequence typing (MLST), have been used to distinguish and type C. glabrata isolates (2, 7, 8, 14, 21, 28). Microsatellite polymorphism analysis has been widely used for typing fungi (3, 4, 11, 15, 18, 25), and this could be an alternative, easy-to-perform, reproducible method suitable for large-scale studies of C. glabrata epidemiology. Recently, Foulet et al. described three polymorphic microsatellite markers to investigate the delineation of clinical C. glabrata isolates (12). The discriminatory power of this method was good but not optimal. The aim of our study was to assess a microsatellite-based multiple-locus variable-number tandem-repeat analysis (MLVA) using new markers for C. glabrata typing.
One hundred twenty-seven C. glabrata strains were analyzed, including four reference strains, 98 independent clinical isolates, and 25 epidemiologically related isolates. These 25 were from the blood cultures and peripheral site isolates of eight patients with C. glabrata candidemia. Genomic DNA was extracted by boiling with Chelex resin as previously described (6, 23). Six microsatellite markers were selected from the C. glabrata DNA sequences available in the GenBank database (29). Primer sequences were designed with Primer3 software (24), and locations in the C. glabrata genome were determined with the Genolevures database (http://cbi.labri.fr/Genolevures/) (Table 1). For each primer set, PCRs were performed with a 20-µl final volume containing 1 µl of DNA, each deoxynucleoside triphosphate at 200 µM, a forward primer and a 5'-dye-labeled reverse primer at 0.25 µM each, and 1 U of Taq DNA polymerase (Promega, Madison, WI). The amplification conditions were 5 min at 95°C, followed by 35 cycles of 95°C for 30 s, 52°C for 30 s, and 72°C for 45 s and then a final step of 72°C for 10 min. For any given isolate, amplicons of each PCR were pooled before multiplex fragment sizing with a Ceq 8000 Genetic Analyzer (Beckman Coulter, Fullerton, CA). Strain IHEM 9670 was run as a control in each experiment. As allele sizing depends on dyes and the analyzer used for electrophoresis, results were expressed as the exact size of the sequence (determined by sequencing of each representative allele for each locus) to allow further interlaboratory comparisons of MLVA results (13, 20). The reproducibility and stability of the method were assessed as described elsewhere (25). Primer specificity was checked by studying the 11 non-C. glabrata reference strains Candida albicans IHEM 9559, Candida dubliniensis IHEM 14280, Candida tropicalis CBS 1920, Candida parapsilosis IHEM 9557, Candida krusei IHEM 9560, Candida norvegensis IHEM 5575, Candida lusitaniae IHEM 10293, Candida nivariensis CBS 9983 and CBS 9984, Candida bracarensis CBS 10154, and Saccharomyces cerevisiae IHEM 6036. Primer specificity was also tested against the sequences in the GenBank database by using BLAST searches.
Sequencing of microsatellites was performed for 30 strains, including reference strains and at least 1 strain for each allele of a given locus by using the CEQ DTSC Quick Start Master Kit (Beckman Coulter, Fullerton, CA).
The discriminatory power (D) of MLVA was calculated with Hunter and Gaston's formula (16). To group the unrelated isolates according to their genetic distance, hierarchical clustering analysis was performed with R software (http://www.r-project.org) and the pvclust package (27). Potential relationships between the genotypes and origins of isolates (clinical data, sex, ward, and anatomical sites) were assessed by hierarchical clustering analysis with canonical discriminant analysis with Tanagra software (http://eric.univ-lyon2.fr/
ricco/tanagra/en/tanagra.html).
Thirty-seven different MLVA profiles of C. glabrata strains were observed (Fig. 1). Three clusters of 23, 20, and 10 genotypically similar isolates were identified (respectively, 22.5%, 19.6%, and 9.8% of the isolates). Twenty-six isolates each gave their own unique MLVA patterns. The discriminatory power (D) of each locus alone varied from 0.64 for Cg7 to 0.79 for Cg10. A D value of 0.902 was reached by combined use of the six markers, but MLVA data obtained by using only four of the six (Cg4, Cg5, Cg6, and Cg10) achieved the same results.

View larger version (32K):
[in this window]
[in a new window]
|
FIG. 1. MLVA-based dendrogram and genotype scores derived from results of the six markers for the 102 unrelated isolates. Hierarchical analysis by the Euclidean distance and the Ward clustering method was performed to classify the 37 genotypes. Allele size is expressed as the exact size of the sequence. The approximately unbiased P values (values on nodes, in percent) were calculated with a multiscale bootstrap (B = 1,000) and are shown only for P values of 80% and genetic distances of 0.2. Reference strain data are identified by the following superscript letters: a, IHEM 19154 and IHEM 19221; b, CBS 138; c, IHEM 9670. Abbreviations: Dig., digestive tract; Gen.Ur., genitourinary tract; Pul., pulmonary tract; Oro., oropharynx; Blood, blood cultures and heart valves; Oth., other sites (skin, nails).
|
|
Sequencing data showed that fragment size polymorphism was explainable by the variation of microsatellite repeats, with the exception of additional polymorphisms observed in the flanking regions for locus Cg6 (poly-T [8 or 9 T residues] and a 16-bp deletion for isolate B473) and Cg4 (a 2-bp deletion at allele 238). Homoplasy (i.e., identical lengths but different DNA sequences) was identified for three loci, Cg7, Cg10 (differences in alternations of CAA and CAG tags in the microsatellite sequence), and Cg6 (mutations in the flanking region). In our investigation, sequencing did not, however, increase the discriminatory power of our MLVA method based on fragment size analysis and thus no further studies were undertaken. Except for loci Cg5 (700-bp amplicon with C. lusitaniae) and Cg11 (120-bp amplicon for C. albicans, multiple bands for C. parapsilosis and C. norvegensis), primers failed to amplify the DNA of non-C. glabrata yeast species. No other potential cross-amplifications were detected with the BLAST software. MLVA markers were reproducible, and none showed any variations or microevolutions even after 300 generations in vitro. This stability was also confirmed after at least 40 days of in vivo maintenance in two patients with C. glabrata candidemia (patients 5 and 6, Table 2). Statistical analysis failed to show any correlations between the epidemiological characteristics of strains and MLVA data. Lastly, for each of the eight patients with C. glabrata candidemia, isolates from peripheral sites and blood culture showed the same MLVA genotype, confirming the endogenous origin of C. glabrata candidemia (Table 2).
In the present study, we investigated six new microsatellite markers for C. glabrata. By combining four of them, a discriminatory power (D) of 0.902, higher than those previously published, was achieved (12). As D values higher than 0.90 are needed for the accurate typing of epidemiologically related isolates (12, 16), microsatellite-based MLVA seems to be an attractive method for large epidemiological surveys and appears to be as discriminant as MLST (D = 0.898 computed from data obtained by Dodgson et al. [9]). Because of their specificity and typeability, microsatellite markers can also be used to distinguish phenotypically and genetically related yeast species, i.e., C. parapsilosis, C. metapsilosis, and C. orthopsilosis in the study by Lasker et al. and C. glabrata, C. nivariensis, and C. bracarensis in our study (1, 5, 18).
At a given locus, most of the possible numbers of repeats were present, demonstrating a continuum for increasing or decreasing the numbers of repeats. However, the distribution of alleles was irregular and several alleles were prominent, since three multilocus genotypes represented 52% of the isolates studied. Other studies using microsatellite or MLST data have highlighted the fact that distinct genetic clades of C. glabrata prevail in different geographical regions (8, 9, 12). Our unrelated isolates were collected from patients from a restricted geographical area, which could partially explain the unequal distribution of genotypes in the population studied. The other hypothesis is that some genotypes could have an ecological advantage (12).
Another finding of our study is that no microsatellite genotypes were associated with any of the clinical data recorded, as previously reported (8, 12, 19). Additionally, in the eight infected patients, we have shown that C. glabrata candidemia originated from colonizing isolates, as previously described for both C. glabrata (12) and C. albicans (4, 22).
The MLVA method, based on fragment size analysis of four microsatellite markers (Cg4, Cg5, Cg6, and Cg10), is easy to perform, discriminatory, and highly reproducible. It appears to be a powerful method for distinguishing epidemiologically related isolates and could be especially useful for studying nosocomial cross-transmission and the kinetics of the colonization-to-infection process. Use of additional microsatellite markers, such as those described elsewhere (12), could perhaps further improve the discriminatory power of our MLVA primer set.
 |
ACKNOWLEDGMENTS
|
|---|
This study was supported by a grant from the French Ministry of Health (PHRC Régional).
We are grateful to Lois Rose for editorial assistance.
 |
FOOTNOTES
|
|---|
* Corresponding author. Mailing address: Mycology-Parasitology Department, CHU Jean Minjoz, 25030 Besançon, France. Phone: (33) 381 66 91 64. Fax: (33) 381 66 8910. E-mail: fgrenouillet{at}chu-besancon.fr 
Published ahead of print on 12 September 2007. 
 |
REFERENCES
|
|---|
- Alcoba-Flórez, J., S. Mendez-Alvarez, J. Cano, J. Guarro, E. Perez-Roth, and M. P. Arevalo. 2005. Phenotypic and molecular characterization of Candida nivariensis sp. nov., a possible new opportunistic fungus. J. Clin. Microbiol. 43:4107-4111.[Abstract/Free Full Text]
- Arif, S., T. Barkham, E. G. Power, and S. A. Howell. 1996. Techniques for investigation of an apparent outbreak of infections with Candida glabrata. J. Clin. Microbiol. 34:2205-2209.[Abstract]
- Bart-Delabesse, E., J. Sarfati, J. P. Debeaupuis, W. van Leeuwen, A. van Belkum, S. Bretagne, and J. P. Latge. 2001. Comparison of restriction fragment length polymorphism, microsatellite length polymorphism, and random amplification of polymorphic DNA analyses for fingerprinting Aspergillus fumigatus isolates. J. Clin. Microbiol. 39:2683-2686.[Abstract/Free Full Text]
- Botterel, F., C. Desterke, C. Costa, and S. Bretagne. 2001. Analysis of microsatellite markers of Candida albicans used for rapid typing. J. Clin. Microbiol. 39:4076-4081.[Abstract/Free Full Text]
- Correia, A., P. Sampaio, S. James, and C. Pais. 2006. Candida bracarensis sp. nov., a novel anamorphic yeast species phenotypically similar to Candida glabrata. Int. J. Syst. Evol. Microbiol. 56:313-317.[Abstract/Free Full Text]
- de Lamballerie, X., C. Zandotti, C. Vignoli, C. Bollet, and P. de Micco. 1992. A one-step microbial DNA extraction method using "Chelex 100" suitable for gene amplification. Res. Microbiol. 143:785-790.[Medline]
- Di Francesco, L. F., F. Barchiesi, F. Caselli, O. Cirioni, and G. Scalise. 1999. Comparison of four methods for DNA typing of clinical isolates of Candida glabrata. J. Med. Microbiol. 48:955-963.[Abstract]
- Dodgson, A. R., C. Pujol, D. W. Denning, D. R. Soll, and A. J. Fox. 2003. Multilocus sequence typing of Candida glabrata reveals geographically enriched clades. J. Clin. Microbiol. 41:5709-5717.[Abstract/Free Full Text]
- Dodgson, A. R., C. Pujol, M. A. Pfaller, D. W. Denning, and D. R. Soll. 2005. Evidence for recombination in Candida glabrata. Fungal Genet. Biol. 42:233-243.[CrossRef][Medline]
- Fidel, P. L., Jr., J. A. Vazquez, and J. D. Sobel. 1999. Candida glabrata: review of epidemiology, pathogenesis, and clinical disease with comparison to C. albicans. Clin. Microbiol. Rev. 12:80-96.[Abstract/Free Full Text]
- Fisher, M., D. Aanensen, S. de Hoog, and N. Vanittanakom. 2004. Multilocus microsatellite typing system for Penicillium marneffei reveals spatially structured populations. J. Clin. Microbiol. 42:5065-5069.[Abstract/Free Full Text]
- Foulet, F., N. Nicolas, O. Eloy, F. Botterel, J. C. Gantier, J. M. Costa, and S. Bretagne. 2005. Microsatellite marker analysis as a typing system for Candida glabrata. J. Clin. Microbiol. 43:4574-4579.[Abstract/Free Full Text]
- Hahn, M., J. Wilhelm, and A. Pingoud. 2001. Influence of fluorophor dye labels on the migration behavior of polymerase chain reaction-amplified short tandem repeats during denaturing capillary electrophoresis. Electrophoresis 22:2691-2700.[CrossRef][Medline]
- Hamal, P., R. Kappe, and D. Rimek. 2001. Rate of transmission and endogenous origin of Candida albicans and Candida glabrata on adult intensive care units studied by pulsed field gel electrophoresis. J. Hosp. Infect. 49:37-42.[CrossRef][Medline]
- Hennequin, C., A. Thierry, G. F. Richard, G. Lecointre, H. V. Nguyen, C. Gaillardin, and B. Dujon. 2001. Microsatellite typing as a new tool for identification of Saccharomyces cerevisiae strains. J. Clin. Microbiol. 39:551-559.[Abstract/Free Full Text]
- Hunter, P., and M. Gaston. 1988. Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J. Clin. Microbiol. 26:2465-2466.[Abstract/Free Full Text]
- Krcmery, V., and A. J. Barnes. 2002. Non-albicans Candida spp. causing fungaemia: pathogenicity and antifungal resistance. J. Hosp. Infect. 50:243-260.[CrossRef][Medline]
- Lasker, B. A., G. Butler, and T. J. Lott. 2006. Molecular genotyping of Candida parapsilosis group I clinical isolates by analysis of polymorphic microsatellite markers. J. Clin. Microbiol. 44:750-759.[Abstract/Free Full Text]
- Lin, C. Y., Y. H. Chen, H. J. Lo, K. W. Chen, and S. Y. Li. 2007. Assessment of Candida glabrata strain relatedness by pulsed-field gel electrophoresis and multilocus sequence typing. J. Clin. Microbiol. 45:2452-2459.[Abstract/Free Full Text]
- Lindstedt, B. 2005. Multiple-locus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria. Electrophoresis 26:2567-2582.[CrossRef][Medline]
- Lockhart, S. R., S. Joly, C. Pujol, J. D. Sobel, M. A. Pfaller, and D. R. Soll. 1997. Development and verification of fingerprinting probes for Candida glabrata. Microbiology 143(Pt. 12):3733-3746.[Abstract]
- Marco, F., S. R. Lockhart, M. A. Pfaller, C. Pujol, M. S. Rangel-Frausto, T. Wiblin, H. M. Blumberg, J. E. Edwards, W. Jarvis, L. Saiman, J. E. Patterson, M. G. Rinaldi, R. P. Wenzel, and D. R. Soll. 1999. Elucidating the origins of nosocomial infections with Candida albicans by DNA fingerprinting with the complex probe Ca3. J. Clin. Microbiol. 37:2817-2828.[Abstract/Free Full Text]
- Provost, F., F. Laurent, R. Camacho Uzcategui, and P. Boiron. 1997. Molecular study of persistence of Nocardia asteroides and Nocardia otitidiscaviarum strains in patients with long-term nocardiosis. J. Clin. Microbiol. 35:1157-1160.[Abstract]
- Rozen, S., and H. Skaletsky. 2000. Primer3 on the WWW for general users and for biologist programmers. Methods Mol. Biol. 132:365-386.[Medline]
- Sampaio, P., L. Gusmao, C. Alves, C. Pina-Vaz, A. Amorim, and C. Pais. 2003. Highly polymorphic microsatellite for identification of Candida albicans strains. J. Clin. Microbiol. 41:552-557.[Abstract/Free Full Text]
- Shemer, R., Z. Weissman, N. Hashman, and D. Kornitzer. 2001. A highly polymorphic degenerate microsatellite for molecular strain typing of Candida krusei. Microbiology 147:2021-2028.[Abstract/Free Full Text]
- Suzuki, R., and H. Shimodaira. 2006. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22:1540-1542.[Abstract/Free Full Text]
- Vazquez, J. A., L. M. Dembry, V. Sanchez, M. A. Vazquez, J. D. Sobel, C. Dmuchowski, and M. J. Zervos. 1998. Nosocomial Candida glabrata colonization: an epidemiologic study. J. Clin. Microbiol. 36:421-426.[Abstract/Free Full Text]
- Wong, S., M. A. Fares, W. Zimmermann, G. Butler, and K. H. Wolfe. 2003. Evidence from comparative genomics for a complete sexual cycle in the asexual pathogenic yeast Candida glabrata. Genome Biol. 4:R10.[CrossRef][Medline]
Journal of Clinical Microbiology, November 2007, p. 3781-3784, Vol. 45, No. 11
0095-1137/07/$08.00+0 doi:10.1128/JCM.01603-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.