Previous Article | Next Article ![]()
Journal of Clinical Microbiology, January 2005, p. 277-283, Vol. 43, No. 1
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.1.277-283.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Grup de Micologia Veterinària, Departament de Sanitat i d'Anatomia Animals, Universitat Autònoma de Barcelona, Bellaterra, Spain
Received 9 July 2004/ Returned for modification 5 September 2004/ Accepted 12 September 2004
|
|
|---|
|
|
|---|
Since the genus Malassezia was created by Baillon in 1889, its taxonomy has been a matter of controversy. A few years ago, the genus was reclassified (14, 16) on the basis of studies of morphological, ultrastructural, physiological, and genetic characteristics, and seven species, which are now widely accepted, were proposed. They are the previously described Malassezia furfur, Malassezia pachydermatis, and Malassezia sympodialis and four new species, Malassezia obtusa, Malassezia globosa, Malassezia restricta, and Malassezia slooffiae.
Recently, four new lipid-dependent species belonging to the genus Malassezia have been described or proposed (19, 25, 32, 33). From these four species, Malassezia dermatis (32), Malassezia nana (19), and the tentatively named "Malassezia equi" (25) are close genetically to the type strain of M. sympodialis, and most of them also have some common morphological and physiological characteristics. M. sympodialis Simmons et Guého 1990 (30) was initially isolated from the auditory tract of a healthy human male and from the scalp of an AIDS patient suffering from tinea capitis. Although its pathogenic role is not yet elucidated, it is commonly isolated from healthy as well as diseased skin (1). M. dermatis was isolated from human patients with atopic dermatitis (32), M. nana was isolated from animals (cat and cows) (19), and "M. equi" was isolated from normal equine skin (25).
In our laboratory, Crespo et al. (3) described for the first time lipid-dependent yeasts associated with otitis externa in cats with similar morphological characteristics and some identical physiological characteristics to the M. sympodialis type species. Lipid-dependent yeasts were also isolated from the skin of horses and different domestic ruminants, as they are the major population of the lipophilic mycobiota in these animals. Difficulty in obtaining a high level of certainty in the identification of some of these lipid-dependent strains by means of these physiological tests was reported (7). The speciation of lipid-dependent isolates from animals by means of physiological tests presents some difficulties, and some of them cannot even be identified (3, 7, 8, 9, 19, 25).
Nucleic acid sequencing has already become an important tool that is useful for the identification of many microorganisms, including yeast species. Molecular systematics of yeasts has emphasized either coding and noncoding regions of the ribosomal DNA. Several studies determined that sequencing of the D1/D2 region could be used for the identification of most basidiomycetous yeast species and that closely related species or strains required sequencing of the internal transcribed spacer (ITS) region (10, 29).
The aim of this work was the DNA characterization of lipid-dependent strains from various domestic animal species close to the M. sympodialis type species by using D1/D2 26S and ITS-5.8S rRNA gene sequencing analysis to understand their phylogenetic relationships and to analyze their specific genetic variation. Phylogenetic relationships with the recently described or proposed M. sympodialis-related species are also analyzed and discussed.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Type strains and lipid-dependent isolates from animals used in this study
|
PCR and DNA sequencing of rRNA genes of the isolates. The variable D1 and D2 regions of the 26S rRNA gene (at the 5' end of the nuclear large subunit of the rRNA gene) were amplified by PCR by using the conserved fungal oligonucleotide primers NL1 and NL4 (26). The amplification process consisted of a predenaturation step at 94°C for 5 min, followed by 30 cycles of denaturation at 94°C for 45 s, annealing at 51°C for 1 min, and extension at 72°C for 3 min plus a final extension of 10 min at 72°C.
ITS and 5.8S rRNA genes were amplified by PCR with the conserved fungal primer pair ITS5 and ITS4 (35). The amplification process consisted of a predenaturation step at 94°C for 5 min, followed by 35 cycles of denaturation at 95°C for 30 s, annealing at 50°C for 1 min, and extension at 72°C for 1 min plus a final extension of 7 min at 72°C.
The amplified PCR products were purified by using a GFX PCR DNA and gel band purification kit (Amersham Pharmacia Biotech, Uppsala, Sweden) by following the supplier's protocol. The purified PCR products were sequenced directly by using a BigDye Terminator, version 3.1, cycle sequencing kit (Applied Biosystems, Gouda, The Netherlands) according to the manufacturer's recommendations. The 26S rRNA D1/D2 and ITS1-5.8S-ITS2 rRNA gene sequences were determined by using NL1 and NL4 and ITS5 and ITS4 as the sequencing primers, respectively. An Applied Biosystems 3100 sequencer was used to obtain the DNA sequences.
Sequence similarity and phylogenetic analyses. The sequences were aligned by using Clustal X software (34). Phylogenetic and molecular evolutionary analyses were performed by using MEGA, version 2.1 (21). For neighbor-joining analysis (28), the distances between sequences were calculated by using Kimura's two-parameter model (20). A bootstrap analysis was conducted with 1,000 replications (11).
Nucleotide sequence accession number. The nucleotide sequences of the D1/D2 26S rRNA gene and ITS-5.8S determined in this study have been deposited in GenBank, and their accession numbers are given in Table 1.
|
|
|---|
![]() View larger version (26K): [in a new window] |
FIG. 1. Molecular phylogenetic tree constructed by using the sequences of D1/D2 26S rRNA gene of members of the genus Malassezia and related Ustilaginomycetes species. The numbers at the branch points are the percentages of 1,000 bootstrapped data sets that supported the specific internal branches. Species with GenBank accession numbers represent sequences obtained from GenBank.
|
![]() View larger version (24K): [in a new window] |
FIG. 2. Molecular phylogenetic tree constructed by using the ITS1-5.8S rRNA gene-ITS2 sequences of members of the genus Malassezia. The numbers at the branch points are the percentages of 1,000 bootstrapped data sets that supported the specific internal branches. Species with GenBank accession numbers represent sequences obtained from GenBank.
|
![]() View larger version (37K): [in a new window] |
FIG. 3. Alignment of D1/D2 sequences of M. sympodialis strains. M. sympodialis CBS 7222 is used as the model. Hyphens designate gaps added to permit alignment. Identical nucleotides are indicated by dots. Lanes: 1, M. sympodialis CBS 7222, MA73, and MA182; 2, MA477; 3, MA88, MA461, and MA470; 4, MA250; 5, MA146; 6, "M. equi" AJ305330; 7, M. dermatis AB070361; 8, MA80, MA125, MA333, MA383, and MA400; 9, M. nana AB075224; 10, MA231, MA237, MA261, MA262, MA314, MA419, and MA424.
|
The sequences of cluster I were close to the M. dermatis AB070361 sequence, which included 579 bp. These sequences differed from the sequence of M. dermatis at seven positions (dissimilarity, 1.21%).
Cluster II included isolates from horses grouping together with the "M. equi" AJ305330 sequence. The nucleotide analysis of the 26S rRNA gene of these strains showed two substitutions over 578 bp among them. The sequences of "M. equi" AJ305330 and MA146 were identical and differed from the sequence of M. sympodialis CBS 7222 at seven positions (dissimilarity, 1.21%). The strains MA88, MA461, and MA470 showed identical sequences and differed from the sequence of M. sympodialis CBS 7222 at nine positions (dissimilarity, 1.56%). The sequence of strain MA250 differed from the sequence of M. sympodialis CBS 7222 at eight positions (dissimilarity, 1.38%).
Cluster III comprised isolates mainly from goats. The section of DNA sequenced included 578 bp, and the sequences of all of the strains were identical, differing from the sequence of M. sympodialis CBS 7222 at eight positions (dissimilarity, 1.38%).
Cluster IV included isolates mainly from cats grouping together with the M. nana AB075224 sequence. The section of DNA sequenced included 577 bp. The sequences of all strains were identical and differed from the sequence of M. nana at 2 positions and from the sequence of M. sympodialis CBS 7222 at 16 positions (dissimilarity, 2.77%).
Sequences of ITS regions and 5.8S rRNA gene. The aligned sequences of the 5.8S rRNA gene with the two intergenic spacers ITS1 and ITS2 are shown in Fig. 4.
![]() View larger version (37K): [in a new window] |
FIG. 4. Alignment of ITS1-5.8S rRNA gene-ITS2 sequences of M. sympodialis strains. M. sympodialis CBS 7222 is used as the model. Hyphens designate gaps added to permit alignment. Identical nucleotides are indicated by dots. Lanes: 1, M. sympodialis CBS 7222; 2, MA182; 3, MA477; 4, MA73; 5, M. dermatis AB070356; 6, MA80, MA125, MA333, MA383, and MA400; 7, MA88; 8, MA146 and MA250; 9, MA470; 10, MA231, MA237, MA261, MA262, MA314, MA419, and MA424.
|
The sequences of cluster I were similar. The length of the section of DNA sequenced ranged from 619 bp in MA477 to 623 bp in MA73. The sequence of MA182 included 621 bp and showed one base pair substitution from M. sympodialis CBS 7222. The nucleotide sequences analysis of these strains showed five substitutions among them and differences of less than four substitutions from those of M. sympodialis CBS 7222.
The sequences of cluster II showed identical ITS2 sequences and differed from the ITS2 sequence of M. sympodialis CBS 7222 at 31 positions. The sequences of strains MA470 and MA461 differed from the sequences of M. sympodialis CBS 7222 at 33 positions, MA146 and MA250 at 32 positions, and MA88 at 30 positions at the ITS1 region. Nucleotide sequence analysis of these isolates showed 3 substitutions among them at the ITS1 region and differences in more than 61 substitutions from those of M. sympodialis CBS 7222 (dissimilarity, 9.8 to 10.3%). Unfortunately, we were not able to analyze the "M. equi" ITS sequences because there is no sequence deposited in GenBank and we could not find any "M. equi" type strain in culture collections.
The sequences of cluster III were identical (620 bp) and differed from the sequence of M. sympodialis CBS 7222 at 39 positions: 18 in the ITS1 region and 21 in the ITS2 region (dissimilarity, 6.28%). The sequences of cluster III were close to the M. dermatis AB070356 sequence.
The sequences of cluster IV were identical (646 bp) and differed from the sequence of M. sympodialis CBS 7222 at 95 positions: 60 in the ITS1 region, 1 in the 5.8S region, and 34 in the ITS2 region (dissimilarity, 15.29%). The ITS1 sequence of M. nana AB075225 was identical to the ITS1 sequences of our isolates from cluster IV.
|
|
|---|
In our study, three isolates from different animal species (horse, lamb, and pig) clustered with the type species M. sympodialis CBS 7222. Only one strain had a single different nucleotide at the D1/D2 sequence. We observed ITS variation, but the strains demonstrated >99% sequence similarity. As previously reported, the gene sequences of ITS1 of Malassezia species are highly diverse (22).
At present, D1 and D2 26S rRNA gene sequences from almost all yeasts are used for species identification or phylogenetic analysis. Peterson and Kurtzman (27) correlated the biological species concept with the phylogenetic species concept by the extent of nucleotide substitutions in 26S rRNA gene sequences. Their study demonstrated that strains of a single biological species show less than 1% substitution in this region. It has been reported (2, 31) that conspecific strains generally demonstrated
99% sequence similarity in the ITS1 and ITS2 regions. According to this species concept, two new lipid-dependent species genetically close to M. sympodialis have been described: M. dermatis and M. nana. These species have similar phenotypes (some common physiological characteristics) to M. sympodialis (Table 2).
|
View this table: [in a new window] |
TABLE 2. Some physiological characteristics of lipid-dependent Malassezia species
|
Cluster II included isolates from horses grouping close to the "M. equi" AJ305330 sequence, showing nearly identical D1/D2 and ITS sequences and being conspecific strains. Dissimilarities between these isolates and the M. sympodialis strain in their D1/D2 regions of 26S and ITS-5.8S were more than 1 and 9%, respectively. At the moment, there is not a valid description yet for "M. equi" and no physiological characteristics have been reported (25).
Cluster III comprised isolates mainly from goats. All strains were completely identical in both 26S rRNA gene and ITS regions. Dissimilarities between these isolates and the M. sympodialis strain in their D1/D2 regions of 26S and ITS-5.8S were 1.3 and 6.28%, respectively. They were conspecific and clearly distinct from M. sympodialis.
M. nana was isolated from animals in Japan and Brazil. It morphologically and physiologically resembles M. dermatis and M. sympodialis, but the phylogenetic trees based on the D1/D2 and ITS1 regions were clearly distinct from these two species (19). The type strain isolated from a cat in Japan was slightly different from the Brazilian strains isolated from cows. Our isolates from cluster IV, mainly from cats, showed the same ITS1 sequence as the M. nana type strain, and the D1/D2 sequences were nearly identical. All strains were completely identical in both 26S rRNA gene and ITS regions, and the dissimilarities between these isolates and the M. sympodialis strain in their D1/D2 regions of 26S and ITS-5.8S were 2.77 and 15.29%, respectively. This group includes two strains associated with otitis externa in cats (MA262 and MA314) (3) with identical physiological characteristics to M. sympodialis (Table 2). The type strain of M. nana (NUSV 1003) also showed physiological characteristics identical to those of M. sympodialis (18) (Table 2).
Special micromorphological characteristics, such as small size of the cells in comparison to other Malassezia spp. (3, 19), buds formed in a narrow base (3, 14, 19, 30), and occasional sympodial budding (3, 14, 30), have been cited for M. sympodialis-related species. However, the separation of Malassezia species based on morphological characteristics is unreliable (16), and their routine identification is currently based on the use of some physiological properties (14, 17), avoiding the use of morphological criteria, which are considered to be more subjective.
On the other hand, the difficulty in obtaining a high level of certainty in the identification of some lipid-dependent strains by means of physiological tests has been reported (3, 7, 8, 9, 18, 25). In addition, the isolation and identification of these strains continue to be difficult due to the low viability associated with some isolate types and the lack of suitable methods for its isolation and preservation (5).
Malassezia species are opportunistic yeasts and belong to the normal cutaneous microbiota. Some of them may act as pathogens when exposed to certain changes in the skin microclimate. In this study, we included lipid-dependent isolates from healthy skin and associated with diseased skin. The association of lipid-dependent isolates with human skin disease is well known. For animals, this association has been recently described (3, 4, 7, 19), and unfortunately, there are few lipid-dependent strains available in culture collections. Using D1/D2 and ITS sequencing analysis, strains associated with skin disease do not differ genetically from isolates from healthy skin. Genetic type IV included isolates from the ears of healthy cats and those from cats and a dog with otitis externa. All strains were completely identical in both the 26S rRNA gene and ITS regions. Usually, this technique does not allow differentiation between the two kinds of isolates. Other molecular techniques, such as randomly amplified polymorphic DNA or amplified fragment length polymorphism, could be more powerful tools to find clinical correlations.
Formal descriptions of possible new taxa will require a standard phenotypic study, which is not a part of this report. Some strains have enough variation to define possible new species or subtypes. Furthermore, some of the sequences included for comparison in the phylogenetic analysis belong to the recently described M. sympodialis-related species and were grouped close or even together with our strains in the proposed clusters. Scorzetti et al. (29), working with other kinds of basidiomycetous yeasts, considered that trees provided by ITS, D1/D2, and other genes, provide focal points for biological and systematic analyses of these yeasts. A sequence, per se, does not describe a species, and descriptions of new species should not rely solely on nucleotide data.
Finally, it is not clear whether the value of these genetic differences is to define species or whether they only demonstrate genetic variation among strains from different origins within M. sympodialis. Further studies are required to investigate whether the divergence between M. sympodialis and our strains is sufficient to resolve them as individual species or whether it only indicates that this species is in the course of differentiation and adaptation to specific animal hosts.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»