Previous Article | Next Article ![]()
Journal of Clinical Microbiology, June 2005, p. 2816-2823, Vol. 43, No. 6
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.6.2816-2823.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Laboratório de Pesquisas Clínicas, Centro de Pesquisas René Rachou, FIOCRUZ, Brasil, and Muséum National d'Histoire Naturelle, Paris, France,1 Departamento de Física, Universidade Federal de Minas Gerais, Minas Gerais, Brazil,2 Muséum National d'Histoire Naturelle, Paris, France,3 Centraalbureau voor Schimmelcultures, Utrecht, The Netherlands, and Institute of Biodiversity and Ecosystem Dynamics, University of Amsterdam, Amsterdam, The Netherlands4
Received 18 November 2004/ Returned for modification 14 January 2005/ Accepted 23 February 2005
|
|
|---|
|
|
|---|
Several comparative molecular techniques have been applied to display polymorphism among strains of black yeasts. Mitochondrial DNA (13) or mitochondrial cytochrome oxidase have been used for diversity at the population level. Chitin synthase (12), in contrast, has insufficient resolution to recognize species but rather displays phylogenetic relationships between species aggregates. The newly described species mentioned above were invariably supported by sequence data of the internal transcribed spacer (ITS) rRNA gene region. Marked ITS differences can sometimes be found in groups that are otherwise monomorphic, such as Ochroconis (G. S. de Hoog and H.-J. Choi, unpublished data). However, the taxonomic value of rRNA gene ITS has been questioned, as it does not provide a sufficient level of resolution in Trichoderma (15), Alternaria (11, 19), and other genera. It is apparent that ITS divergence rates cannot be used as a gold standard for taxonomic differences all over the fungal kingdom, but particular levels of divergence have been reached in some groups, among which are the black yeasts and their allies. In a number of cases (4, 5, 27) the entities found were found to coincide with phenetic characters that had been neglected or overlooked in earlier taxonomic systems, underlining the taxonomic predictivity of ITS characters.
ITS divergence between black yeast species is marked, species being clearly separated by at least 1% sequence diversity (3). A major drawback of this ITS diversity in black yeasts is the large degree of divergence when the entire ascomycete family Herpotrichiellaceae, to which they belong through Capronia teleomorphs (24), is compared. Such comparisons may be necessary, because morphology does not provide a priori clues to phylogenetic positions (9). Over larger phylogenetic distances, large parts of the ITS region may not be alignable with confidence (20) and therefore have been excluded from the analysis.
Recently, the algorithm DNA-Walk Divergence (DNAWD) has become available (14). The algorithm's robustness relies on global alignment of complete sequences, which resolves the arbitrariness often associated with the precise identification of indel mutation sites, while the integral nature of the DNA-walk still allows for the sensitivity in their occurrence. This allows for comparison of entire ITS sequences including a broad taxonomic distance range. It can be shown that a properly defined DNA-walk is the exact equivalent to the sequence composition, and therefore mutations and dislocations produce divergences in the walk geometry which can be measured. In this paper the DNAWD program is applied to the phylogenetic reconstruction of black yeast strains, exploring total ITS domains of the rRNA gene. A small-subunit (SSU) rRNA gene tree is used to confirm grouping as given by ITS trees.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Strains analyzed with collection reference and GenBank accession numbersa (SSU numbers of those species that were also analyzed for ITS data are also mentioned)
|
ITS sequences were aligned using the BioNumerics package (Applied Maths, Kortrijk, Belgium). Due to gaps necessary for alignment, the ITS 1 domain spanned 258 positions (real lengths, 197 to 215 bp), the 5.8S gene 157 positions, and the ITS2 domain 224 positions (real lengths, 171 to 211 bp). Positions 97 to 155 and 184 to 208 (ITS1) and 456 to 484 and 541 to 631 (ITS2) could not confidently be aligned and were excluded from most of the comparisons. The same alignment was used for a comparison with different packages applying different algorithms: Phylip package (v.3.572c) (6), according to maximum likelihood; BioNumerics according to Ward's averaging (27), without Indels, using Kimura's two-parameter model, distance estimation through neighbor joining using the Treecon package with Kimura's two-parameter model (25), and Paup (23) according to parsimony. In all comparisons, insertions and deletions (indels) were not taken into account. Exophiala pisciphila (AF050273) was taken as the outgroup.
Homology analysis by DNA-walk divergence. DNA-walks are defined by incrementing walk steps for each nucleotide in the sequence (for example, a positive step for purines and negative for pyrimidines). The DNA-walk divergence method (14) makes simultaneous three-dimensional walk comparisons (representing three composition skews), AG-TC, AC-TG, and AT-CG for each pair of sequences. One sequence slides against the other until the minimum squared walk difference is found, corresponding to a global alignment. This is then taken as a measure of their distance since statistically independent mutations and indels increase the mean square walk differences linearly. The resulting distance matrices are then fed to the Kitsch program of the Phylip package (v. 3.572c) (6), which generates trees with contemporary leaves.
The analyses were carried out for both the ITS1 and ITS2 regions separately, defined within the consensus limits ATCATTA to CAACAAC and TGTTGGA to AGGTTGA, respectively (Fig. 1a). The choice of the intervals was based on apparent homogeneity of sequence dissimilarities among species and can be appreciated in the DNA-walk graph of Fig. 1b. The lengths and average dissimilarities of the sections analyzed are also given in Fig. 1a.
![]() View larger version (19K): [in a new window] |
FIG. 1. a. Diagram showing ITS regions defined from strict consensus sequences used in DNAWD analysis. The mean size in base pairs and mean percent variability among species for each region are also noted. b. AT-CG DNA-walk coordinates for seven different species (P. verrucosa, P. americana, E. dermatitidis, F. pedrosoi, E. pisciphila, C. carrionii and C. bantiana). The ITS domains used in DNAWD analysis are also shown.
|
|
|
|---|
![]() View larger version (55K): [in a new window] |
FIG. 2. Consensus tree of 40 strains of black yeasts and related fungi based on confidently aligned rRNA gene ITS sequences, i.e., excluding positions 97 to 155 and 184 to 208 (ITS1), and 456 to 484 and 541 to 631 (ITS2), using the neighbor joining algorithm in the Treecon package with Kimura-2 correction. Bootstrap values of >90 from 100 resampled data sets are shown: solid bar, branches also recognized with Ward's averaging; dotted bar, branches also recognized with UPGMA; dashed bars, branches also recognized with parsimony.
|
![]() View larger version (52K): [in a new window] |
FIG. 3. Neighbor-joining tree based on 1,660 positions of 70 SSU rRNA gene sequences generated with the ARB package. Phaeomoniella chlamydospora CBS 101359 was used as the outgroup. Strains with numbers 1 to 5 are proven to belong to SSU groups with the same numbering indicated by vertical bars, the cluster delimitation derived from groups recognized by Haase et al. ITS groups (A to G) are superimposed, indicated with arrows.
|
![]() View larger version (16K): [in a new window] |
FIG. 4. Unrooted tree made with the Kitsch program in the Phylip package for construction of the DNAWD distance matrix of the entire ITS region of 40 black yeasts and related fungi. Group indications are those of the rRNA gene ITS in Fig. 2 (letters) and the SSU rRNA gene recognized by Haase et al. in Fig. 3 (numbers).
|
![]() View larger version (22K): [in a new window] |
FIG. 5. Unrooted tree made with the Kitsch program in the Phylip package for construction of the DNAWD distance matrix of ITS2 of 40 black yeasts and related fungi. The radial scale measures twice the percent dissimilarity in composition. Branches crossing the dotted circle (phenon circle) group about 13 species-aggregates with dissimilarities below 6.3%. The inner circle estimates large-distance dispersion and covers the most unreliable branching far away from the leaves.
|
|
|
|---|
Our selection of species and strains analyzed is more or less representative for the hyper-variability of the ITS domains of black yeasts and allies. With growing distances between strains, alignments quickly can no longer be done with confidence. Meaningful comparison to reveal evolutionary relationships between some of the species is then only possible when small selections of the ITS domains are taken into account, as mentioned in Materials and Methods. Phylogenetic interpretation is strongly hampered, as trees lack any resolved substructure, being based on too low number of mutations in the small portions that are still comparable. Due to the large distance between most of the groups, the clusters were robust and remained strictly identical with any algorithm applied, particularly in the lack of a hierarchical structure of main branches (Fig. 2).
A substructure, i.e., with separate species in a more or less hierarchical ordering, was revealed with ITS in the Phialophora/Cladophialophora carrionii cluster (group A) only. When the neighbor joining algorithm was applied to the complete ITS domain, which could be confidently aligned in the mentioned cluster, C. carrionii was found paraphyletically to Phialophora. This matches phenetic data because members of the clade share potential production of phialides and are agents of human chromoblastomycosis. These characters are also found in Fonsecaea (group C). In the SSU tree of the Herpotrichiellaceae published by Haase et al. (9) Fonsecaea has members of Cladophialophora other than C. carrionii as sister group species (cluster 2 of Haase et al.) (9). Cladophialophora species other than C. carrionii cause phaeohyphomycoses rather than chromoblastomycoses.
In our updated SSU tree of the same family (Fig. 3) members of Fonsecaea, Phialophora, and Cladophialophora (A to C) all are united in a single clade (2), showing phylogenetic coherence of the Herpotrichiellaceae with high degrees of virulence. In the DNA-walk tree based on ITS2 (Fig. 5) Fonsecaea, Phialophora, and Cladophialophora are all found in a main subdivision of the tree (II) where distances between species are smaller than the isolated, more distantly related taxa opposite of the dotted line (I) in Fig. 5. The radial scale measures twice the percent dissimilarity in composition. Branches crossing the dotted circle (phenon circle) group about 13 species-aggregates with dissimilarities below 6.3%.
With DNA-walking, the analysis of ITS domain data resulted in slightly different trees when either the complete ITS region (Fig. 4) or ITS2 alone (Fig. 5) was used. All trees recognized all species correctly. The question of whether the structure is phylogenetically interpretable can be verified by a comparison with SSU rRNA gene data of the same fungi. The most comprehensive phylogenetic overview of the family is that of Haase et al. (9). These authors noted that the tree was poorly resolved but recognized five approximate clades, with Exophiala dermatitidis, Cladophialophora bantiana, E. spinifera, E. nigra, and Coniosporium perforans as core species, as listed above. Ecologically these groups seem to be meaningful. Clade 1 is a thermophilic yeast group around the neurotrope E. dermatitidis (10), clade 2 is a group with virulent species with accent on species causing chromoblastomycosis and brain disease (7), clade 3 is the E. spinifera/E. jeanselmei complex (5), clade 4 is a meso- to psychrophilic group of Exophiala species in showers and ocean waters (9), and clade 5 contains meristematic species inhabiting rocks (22). This ecological support for the groups suggests that this reflects the phylogeny of these groups optimally.
In the SSU tree based on our expanded and updated data set (Fig. 3) some further, ecologically defined clades can be located. ITS groups F and G are psychrophilic species similar to the fish pathogen E. pisciphila, although unexpectedly E. salmonis is found in E, which is a cluster of mutually distantly related environmental species. In the SSU tree these taxa are mutually remote (Fig. 3). Exophiala oligosperma is a close relative of E. jeanselmei (5) and indeed is found together with this species in ITS cluster C (Fig. 5), matching the E. spinifera complex SSU group 3 (Fig. 3). Some dispersed strains previously listed as "E. spinifera" (Fig. 3) are now known to represent individual taxa, such as IFM 14855, which has been described as E. nishimurae (26).
With the DNA-walk of the entire ITS domain, no coherent substructure was recognizable. SSU-clade 2/ITS-clades A to C were found at the longest mutual distances (Fig. 4). In contrast, a major subdivision (I and II) was consistently recognized based on ITS2 alone, with E. pisciphila CBS 464.81, representing a distantly related, undescribed species, in an external position (G). Main branches I and II contained five to nine subclusters representing the individual species analyzed (Fig. 5). C. carrionii and Phialophora (2/A) were found to be interrelated, with Fonsecaea (2/C) branching off at the base; all are agents of chromoblastomycosis. Cladophialophora bantiana (2/B) was found as well-individualized clusters paraphyletically to psychrophilic Exophiala species (F).
In conclusion, the DNA-walk showed consistent results in the recognition of individual species, whether these consisted of several strains or a single strain. This was achieved with the complete ITS domain as well as with rRNA gene spacers 1 and 2 separately. In addition, ITS2 alone also enabled phylogenetic reconstruction to some extent, which was not possible with any other algorithm using the same data set. DNA-walk divergence proved to be a powerful tool for the analysis of ribosomal genes and their evolutionary interpretation. The gold standard for phylogenetic interpretations is the small or large rRNA gene subunits, which in general contain more information. However, in the case of the black yeasts and their relatives, these genes often contain degrees of variability too small to resolve individual species (9). DNA-walk divergence bridges the gap between current methods of phylogenetic reconstruction. Another straightforward application of this technique could be the molecular identification of the agents of phaeohyphomycosis and chromoblastomycosis isolated in the clinical laboratory.
Bert Gerrits van den Ende and Adile Cegcel are thanked for technical assistance.
|
|
|---|
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»