Previous Article | Next Article ![]()
Journal of Clinical Microbiology, September 2004, p. 4253-4260, Vol. 42, No. 9
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.9.4253-4260.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Division of Dermatology, Department of Medicine, Sunnybrook and Women's College Health Science Center, the University of Toronto, Toronto,1 Mediprobe Laboratories, London, Ontario, Canada,2 Centraalbureau voor Schimmelcultures, Utrecht, The Netherlands3
Received 8 January 2004/ Returned for modification 3 March 2004/ Accepted 28 March 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The genus Malassezia, until recently, was characterized on the basis of rRNA sequences as consisting of seven species, including the lipid-dependent species M. furfur, M. sympodialis, M. globosa, M. obtusa, M. restricta, and M. slooffiae and the lipophilic species M. pachydermatis (12, 13). Recently, two new species have been identified: M. dermatis (39), which was isolated from atopic dermatitis patients, and M. equi (unpublished), a species that was isolated from the skin of horses (32). The latter species has only tentatively been named and still awaits formal description; the present study deals only with the eight described species.
Guillot et al. (17) introduced a physiological system based on lipid assimilation and other phenotypic characteristics for identifying the various Malassezia species. This phenotypic system has been used as the conventional method of identification, though in practice, the test results are not always easy to read. Therefore, several research groups have explored the use of molecular techniques, such as pulsed-field gel electrophoresis (3, 4, 36), randomly amplified polymorphic DNA analysis (1, 3), amplified fragment length polymorphism (AFLP) analysis (40), denaturing gradient gel electrophoresis (40), multilocus enzyme electrophoresis (30), sequencing analysis (37), restriction analysis of PCR amplicons of ribosomal sequences (8, 15, 16, 20, 28), and chitin synthase gene sequence analysis (1, 4, 8, 37), for the identification of Malassezia species. Recently, Gemmer et al. (9) devised a remarkably efficient, novel technique, terminal fragment length polymorphism analysis, for the rapid and reliable identification of Malassezia species. It eliminates the need for strain cultivation in direct investigations of Malassezia populations on skin samples, as well as the need for restriction enzyme digestion of nucleic acids. The technique, however, is not suitable for epidemiological typing, as its ability to show heterogeneity within a species is limited.
The objective of the present study was to facilitate an improved understanding of the epidemiology of Malassezia infections by determining which molecular characterization methodologies could most effectively be used for the identification and strain typing of Malassezia species. It is important to identify the different species correctly, as some studies have given preliminary indications that different Malassezia species may occupy well-defined niches on the human body and may play differential roles in causing various diseases (11, 19, 21). However, these studies have not all been in agreement about which species is prevalent or etiologically important with regard to a specific disease, suggesting that technical improvements in identification methods is critical in producing definitive results.
AFLP analysis is a universally applicable technique that has proven to be useful for the identification and strain typing of microorganisms (2, 5, 34). In this study, we investigated the use of AFLP analysis as well as rapid sequencing of the internal transcribed spacer (ITS) and large-subunit (LSU) regions of the nuclear ribosomal DNA (rDNA) for the identification of the eight currently recognized Malassezia species. We compared the two molecular methods with the conventional method of identification to see if any discrepancies could be observed. The molecular techniques were used to investigate the genetic diversity of Malassezia isolates obtained from a variety of patients in Canada. These patients were mostly from the province of Ontario, a region with a notably diverse population developed through immigration from many areas of the world. In order to understand the extent to which the Ontario isolates might reflect the overall genetic biodiversity of Malassezia spp. from human sources and to ensure that the techniques used were maximally developed as tools for identifying genetic types within Malassezia species, a diverse comparison sample of Malassezia strains from different worldwide geographic regions was chosen from the collection of the Centraalbureau voor Schimmelcultures (CBS; Utrecht, The Netherlands).
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
AFLP analysis. AFLP analysis was performed according to the manufacturer's instructions in the AFLP microbial fingerprinting protocol (PE Biosystems, Nieuwerkerk aan de IJssel, The Netherlands), with some modifications. Restriction and ligation were performed simultaneously on 10 ng of genomic DNA by using 1 U of MseI, 5 U of EcoRI, and 3 U of T4 DNA ligase (Biolabs, Westburg, The Netherlands). The sequences of the primers EcoRI and MseI were 5'-GACTGCGTACCAATTCAC-3' and 5'-GATGAGTCCTGAGTAAC-3', respectively. The adaptors used were EcoRI (5'-CTCGTAGACTGCGTACC-3', forward; 3'-CATCTGACGCATGGTTAA-5', reverse) and MseI (5'-GACGATGAGTCCTGAG-3', forward; 3'-CTACTCAGGACTCAT-5', reverse). The reaction was allowed to take place in a total volume of 5.5 µl with the following constituents: a 0.36 µM concentration of the EcoRI adaptor and a 3.64 µM concentration of the MseI adaptor from the AFLP microbial fingerprinting kit (PE Biosystems), 0.1 M NaCl, 0.91 mM Tris-HCl (pH 7.8), 0.18 mM MgCl2, 0.18 mM dithiothreitol, 18 µM ATP, and 91.36 µg of bovine serum albumin ml1. The restriction ligation mixture was incubated for 2 h at 37°C and later diluted by adding 25 µl of sterile double-distilled water. The first PCR was performed with two preselective primers (EcoRI core sequence and MseI core sequence) and the AFLP amplification core mix from the AFLP microbial fingerprinting kit, according to the manufacturer's manual, under the following conditions: 2 min at 72°C, followed by 20 cycles of 20 s at 94°C, 30 s at 56°C, and 2 min at 72°C each. The PCR product was diluted by adding 25 µl of sterile double-distilled water. The second PCR used more-selective primers, EcoRI-A FAM and MseI-G. The conditions were 2 min at 94°C; 10 cycles consisting of 20 s at 94°C, 30 s at 66°C (decreasing 1°C every step of the cycle), and 2 min at 72°C; and then 25 cycles consisting of 20 s at 94°C, 30 s at 56°C, and 2 min at 72°C. The samples were prepared for acrylamide capillary electrophoresis with the following loading mix: 2.0 µl of selective amplification product, 24 µl of deionized formamide, and 1 µl of GeneScan-500 labeled with 6-carboxy-X-rhodamine (Applied Biosystems, Foster City, Calif.) as an internal size standard. After incubation for 5 min at 95°C, the samples were run on an ABI 310 genetic analyzer for 30 min each. Data were analyzed with the Bionumerics software package (version 2.5; Applied Maths, Kortrijk, Belgium), by using (i) Pearson correlation based on similarities of the densitometric curves and (ii) the unweighted pair group method with arithmetic means analysis. The statistical reliability of the clusters was confirmed by using the cophenetic values, which calculate the correlation between the calculated similarities and the dendrogram-derived similarities.
PCR amplification and sequencing of ITS and LSU regions. The ITS I and ITS II regions and the D1 and D2 domains of the LSU region were amplified with primers V9 (5'-TGCGTTGATTACGTCCCTGC) and RLR3R (5'-GGTCCGTGTTTCAAGAC). The PCR conditions were 5 min at 94°C; 35 cycles of 45 s at 94°C, 40 s at 56°C, and 2 min at 72°C; followed by chilling at 4°C. The PCR products were purified by using GFX columns (Amersham Pharmacia Biotech Inc., Roosendaal, The Netherlands) and visualized on electrophoresis gel after ethidium bromide staining. The rDNA was sequenced with the BigDye terminator cycle sequencing kit (Applied Biosystems) and analyzed on an ABI Prism 3700 sequencer (Applied Biosystems) by using the standard conditions recommended by the vendor. The primers used in the sequence reaction were ITS4 (5'-TCCTCCGCTTATTGATATGC) and ITS5 (5'-GGAAGTAAAAGTCGTAACAAGG) for analysis of the ITS region and NL1 (5'-GCATATCAATAAGCGGAGGAAAAG) and RLR3R (5'-GGTCCGTGTTTCAAGAC) for analysis of the D1 and D2 domains of 26S rDNA. The PCR conditions were as follows: 25 PCR cycles of 96°C for 10 s (denaturation), 50°C for 5 s (annealing), and 60°C for 4 min (extension). The sequencing products were purified with Sephadex (Amersham Pharmacia). Sequences were assembled and edited with Seqman II software (DNAStar Inc., Madison, Wis.) and aligned with Megalign (DNAStar). The sequences were visually corrected. Phylogenetic analysis was performed by using PAUP (version 4) parsimony analysis, random stepwise addition, and tree bisection-reconnection.
Statistical analysis. Likelihood ratios were computed by using SPSS for Macintosh (release 11.0.2, 2003; SPSS, Chicago, Ill.).
Nucleotide sequence accession numbers. The nucleotide sequences determined in this study have been deposited in GenBank. The GenBank numbers for D1 and D2 domain and ITS sequences are listed in Table 1.
| RESULTS |
|---|
|
|
|---|
|
|
Because many of the subtype 1 and 6 isolates were from the trunk (back and chest), we tested whether these genotypes might be significantly associated with this body site. Isolate numbers permitted statistical testing for the significance of differences in M. furfur subtypes among isolates from the trunk, the head (including the scalp and ear, whether neonatal or not), and internal sites or sites of apparent systemic invasion (including catheter sites). In this isolate breakdown, differences among AFLP types were not significant (likelihood ratio, 13.039; P = 0.221).
The delineation of possible host differences among the subtypes was rendered problematic by the superimposed geographic differences among the isolates. Subtypes 1 and 6 contained mainly isolates from Ontario (11 of 13 total isolates), but the majority of isolates from adult skin (11 of 17) also came from the Ontario sample. The subtype 4 isolates from internal sites as well as the internal and neonatal isolates of subtype 5 were from Western Europe (Germany, France, The United Kingdom, and The Netherlands), but that region was the source of only one of the nonneonatal human skin isolates available for study. That isolate belongs to subtype 5. Statistical testing of AFLP type distributions between North America and Europe showed that the differences seen were significant (P < 0.001; likelihood ratio, 23.689); however, it was not possible to concomitantly test whether the differences observed were truly attributable to geographic factors rather than to the accompanying differences in host body sites. One M. furfur isolate, CBS 6000, clustered separately and was designated type 7. It was from human dandruff sampled in India. The only other Asian isolate in the sample was from adult skin sampled in Hong Kong. It belonged to subtype 6.
M. sympodialis formed a coherent AFLP cluster. The majority of strains were from the backs or chests of healthy subjects from Ontario. Some of the M. sympodialis strains studied were from patients who had skin lesions. The lesional and nonlesional isolates showed no genetic differences based on their AFLP types.
Tight clustering was seen among M. sympodialis and M. dermatis isolates; but in contrast, the clusters of M. pachydermatis, M. globosa, M. restricta, and M. slooffiae were relatively loosely structured. The low internal similarity values (Table 1) seen within each of these four species support the impression that they are genetically diverse species.
(ii) Putative new species. There were some strains that fell into a cluster that was not associated with any recognized Malassezia taxon (Fig. 1). These strains appear to belong to a species that is as yet undescribed.
LSU rDNA D1 and D2 domain and ITS sequence analysis. Sequence analyses of the LSU and ITS regions resulted in eight well-separated, distinct groups representing the eight different Malassezia spp. The D1 and D2 domain and ITS sequences were in the range of 550 and 600 bp, respectively. The nucleotide sequences determined in this study have been deposited in the CBS and GenBank databases, which may be further consulted for identification purposes (www.cbs.knaw.nl/databases/index.htm and www.ncbi.nlm.nih.gov). The phylogenetic trees based on the ITS and LSU sequences were in complete agreement with each other. Some species turned out to be genetically heterogenous with respect to the D1 and D2 domain and ITS sequences. This was true in particular for M. furfur (D1, D2, and ITS), M. sympodialis (ITS), M. pachydermatis (D1, D2, and ITS), M. globosa (D1, D2, and ITS), and M. slooffiae (D1, D2, and ITS). Phylogenetically, M. furfur and M. obtusa appear to be sister species. A sister species relationship is also found with M. restricta and M. globosa. In both trees, the phylogenetic origin of M. slooffiae was unclear.
Comparison of physiological and molecular methods. The results of the two molecular methods used in our study were concordant. However, some discrepancies were observed when the AFLP analysis and sequencing results were compared with the results of physiologically based identification. One isolate (98F) that had previously been identified as M. furfur was reidentified as M. sympodialis, while three isolates (VMB3, GMB1, and WF39) that had previously been identified as M. globosa were reidentified as M. sympodialis. One isolate (JSB2) that was formerly identified as M. slooffiae was reidentified as M. sympodialis, and one (SF2) that was formerly identified as M. sympodialis was recognized as M. slooffiae. Two isolates (GM420 and 4SJCMC1) which had previously been identified as M. pachydermatis and M. globosa, respectively, were reidentified as M. furfur. Overall, a misidentification rate of 13.8% was observed.
| DISCUSSION |
|---|
|
|
|---|
Theelen et al. (40) established AFLP analysis as a useful discriminatory technique for Malassezia species identification. AFLP analysis has the capacity to assay a much greater number of loci for polymorphisms than are surveyed by other PCR-based techniques (2) because it is based on the ligation of known sequences (adaptors) to a wide range of restriction fragments; these adaptors then function as targets for PCR primers (2, 34). The use of an internal size standard with every sample for normalization purposes greatly enhances the reproducibility of the results (5). Unlike results from most other PCR-based methods, which analyze only a part of the genome, the banding patterns in AFLP analysis illustrate a broad-ranging subsample of the whole genome and are much easier to analyze than those resulting from restriction fragment length polymorphism and random amplified polymorphic DNA analysis. Another advantage of AFLP analysis is that the patterns can be stored in an accessible general database for future comparison and species identification of additional Malassezia isolates.
Sequencing of the two most varied domains, D1 and D2, of the 26S rDNA has proved useful at both the phylogenetic and taxonomic levels, as the region is sufficiently varied to allow distinction between species, in particular basidiomycetes (27, 35). On the other hand, sequencing of the ITS region has also been found to be useful in discriminating closely related species (10, 23). In the case of Malassezia spp., Gupta et al. (20) observed that PCR-restriction fragment length polymorphism analysis of the ITS region was sufficient to resolve the differences between the physiologically similar species M. furfur, M. sympodialis, and M. slooffiae. Further, sequence diversity within various species has been observed, which suggests the presence of several genotypes within the species (11, 37). In keeping with these observations, we decided to sequence both of these regions and also to add AFLP analysis to further refine the precision of genotypic clustering in this group and to facilitate accurate species identification.
Following up on the studies of Theelen et al. (40), and testing numerous new clinical isolates, we confirmed that the AFLP analysis yields high-quality fingerprints and is species specific in typing Malassezia isolates, though some species contain multiple subtypes. We were able to characterize the major genotypes, including several novel ones, that were present in the Malassezia isolates from Ontario. AFLP clusters seen in M. furfur in our study suggested the existence of some previously undisclosed specificities. For example, our M. furfur genotype 1 isolates came mainly from sites on the chest and back, a remarkable finding given that many of the other Malassezia isolates screened came from body sites, such as the exposed areas of the face, that are only subtly physiologically differentiated from the skin of the trunk. In addition, subtypes 1 and 6 appeared to be especially associated with nonneonatal skin, a habitat in which populations of M. furfur isolates might be expected, to some extent, to have attained an ecological balance among themselves and with other species (i.e., to be living in a way approximating the ecological k-strategy or "climax community") (7). On the other hand, subtypes 4 and 5 were obtained mostly from sites of internal invasion in immunocompromised patients or from neonatal skin, i.e., from classic newly available sites potentially allowing aggressive, ruderal-type (ecological r-strategy) invasion. Human skin may support such strongly ruderal genotypes in relatively small numbers except in sites of skin disturbance, where they may have an advantage. (M. furfur, as a whole, with its relatively rapid growth and minimally fastidious character in vitro, may already constitute one of the more ruderal components of the spectrum of skin Malassezia species, as would be predicted from its regular involvement in nosocomial infection.)
We could not rule out the possibility that the distributional differences between the subtype 1 and 6 group and the subtype 4 and 5 group reflected not ecology but rather a difference in geographic distribution or a sampling artifact simulating such a difference. Subtypes 1 and 6 were mainly from the Ontario sample, consisting mainly of isolates from adult skin, and subtypes 4 and 5 were mainly from Western Europe, where available M. furfur isolates came mostly from in-hospital studies. It appears highly unlikely that Ontario could have strongly geographically specific Malassezia populations in contradistinction to those of Western Europe, since most Ontarians are of western European origin within their own lifetime or within a small number of generations. Without extensive follow-up sampling, however, the underlying causality of the distributional differences established among subtypes 1, 4, 5, and 6 cannot be discerned. In particular, a sample of hospital-associated isolates from Ontario and skin isolates from adult western Europeans would be of value. A relatively large sample of Malassezia isolates of all species, however, would have to be processed in order to capture sufficient M. furfur isolates to resolve the questions raised by the present study. Development of selective isolation media for individual Malassezia species would be of value.
The disclosure of the intriguing differences among human-associated subtypes, as well as the finding of two genotypes strongly associated with animal sources, dramatically illustrates the power of AFLP analysis in allowing us to discover potentially important epidemiological patterns within Malassezia species. AFLP analysis was also able to discriminate Malassezia strains that could not be identified by physiological tests (Fig. 1). As mentioned above, it is probable that these strains belong to an as-yet-undescribed species. Further investigation is in progress.
The discrepancies observed in our study between phenotypic and molecular methods of identification were comparable to those observed by Makimura et al. (28), who used ITS1 sequencing as a molecular identification standard. They observed that out of 46 clinical isolates that had formerly been identified as M. furfur by employing the conventional phenotypic approach, 22 were M. sympodialis and 5 were M. slooffiae. It is worth mentioning, however, that misidentification may not be the only potential source of such findings. Malassezia, as a lipophile, is relatively hydrophobic, and reliable single-cell cultures may be relatively hard to make. Clumping occurs in water, especially if oil from the medium is present, unless relatively vigorous efforts are made. It appears that cultures of slow-growing Malassezia spp. may conceal relatively inactive inocula of certain faster-growing species and that when the cultures senesce, the contaminating inocula are activated in a way that causes the slow-growing species to be overgrown. Nakabayashi et al. (31) observed that M. globosa, which had been detected in a primary culture, disappeared in several experiments and that only M. sympodialis remained in the culture. More than one Malassezia species may be present in a clinical sample or even in an apparent single colony on a contact plate, necessitating considerable effort to ensure that each Malassezia colony submitted for identification is in a pure state. In some situations, pharmacological selection pressure derived from topical drugs or hygienic materials may also influence the isolation of Malassezia spp. from clinical samples. Trace quantities of inhibitors carried over into primary cultures may result in the growth of some species being partially repressed until later subcultures attenuate the inhibitor levels. More research is required to find out whether the discrepancies in identification are due to technological error in the physiological methods or to sampling artifacts as mentioned above.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|