Microbiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Microbiology 152 (2006), 1539-1549; DOI  10.1099/mic.0.28626-0
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lan, L.
Right arrow Articles by Xu, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lan, L.
Right arrow Articles by Xu, J.
Agricola
Right arrow Articles by Lan, L.
Right arrow Articles by Xu, J.
Microbiology 152 (2006), 1539-1549; DOI  10.1099/mic.0.28626-0
© 2006 Society for General Microbiology

Multiple gene genealogical analyses suggest divergence and recent clonal dispersal in the opportunistic human pathogen Candida guilliermondii

Lisa Lan and Jianping Xu

Department of Biology, McMaster University, 1280 Main St West, Hamilton, ON L8S 4K1, Canada

Correspondence
Jianping Xu
jpxu{at}mcmaster.ca


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Candida guilliermondii is a haploid opportunistic pathogen accounting for about 2 % of human blood yeast infections. Recent analyses using multilocus enzyme electrophoresis and karyotyping suggest that strains from human sources traditionally designated C. guilliermondii in fact include at least two species, C. guilliermondii and Candida fermentati. However, the patterns of molecular variation within and between these two species remain largely unknown. In this study, DNA fragments were sequenced from five genes for each of 37 strains collected from Canada, China, the Philippines and Tanzania. The analyses identified significant sequence differences between C. guilliermondii and C. fermentati. The five gene genealogies showed no apparent incongruence, suggesting a predominantly clonal reproductive structure for both species in nature. Indeed, two large clones of C. guilliermondii were identified, with one from Ontario, Canada, and the other from China. Interestingly, the results indicate that strains currently designated C. guilliermondii may contain additional divergent lineages. On the practical side, the results revealed several diagnostic molecular markers that can be used in clinical microbiology laboratories to distinguish C. guilliermondii and C. fermentati. The multiple gene genealogical analyses conducted here revealed significant divergence and clonal dispersal in this important pathogenic yeast complex.


Abbreviations: IA, index of association; K2P, Kimura two-parameter; MP, maximum parsimony; PH, partition homogeneity; T-PTP, topology-dependent permutation tail probability


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Emerging infectious diseases are increasingly recognized as significant threats to human health and social and economic wellbeing. Most emerging infections are believed to result from anthropogenic activities such as increased contact between previously disjunct human and wildlife populations, the disturbance of human and animal ecological niches, and the rapid dispersal of infectious agents accompanying human travel. A prominent feature of emerging human infectious agents is the widespread geographical distribution of clones and genotypes. For example, in the model human fungal pathogens Candida albicans and Cryptococcus neoformans, strains with similar or identical multilocus genotypes are commonly found among diverse geographical locations in South and North America, Australia, Eurasia and Africa (Kidd et al., 2005Down; McEwen et al., 2000Down; Xu, 2004aDown; Xu et al., 1999bDown, cDown, 2000bDown). While most of our attention has focused on model pathogens, the population biology of other opportunistic fungal pathogens such as Candida guilliermondii remains poorly understood.

C. guilliermondii is a haploid ascomycete yeast (Doi et al., 1992Down; Imshenetskii et al., 1979Down). Like most other species in the genus Candida, C. guilliermondii is a common constituent of the normal human microflora (Kam & Xu, 2002Down; Odds, 1988Down; Tey et al., 2003Down; Xu & Mitchell, 2003aDown). However, when the immune system of the host is compromised due to organ transplantation, chemotherapy, or infections by other agents such as the human immunodeficiency virus, C. guilliermondii and other Candida species can cause severe, sometimes life-threatening infections (Xu et al., 1999aDown). Whereas Candida albicans is the most commonly isolated pathogenic yeast in humans, other species such as C. guilliermondii are commonly found (Dos Santos & Soares, 2005Down; Krcmery & Barnes, 2002Down; Odds, 1988Down; Sandven, 2000Down; Tey et al., 2003Down). Globally, C. guilliermondii accounts for about 2 % of all candidaemia (Krcmery & Barnes 2002Down; Sandven, 2000Down; Tey et al., 2003Down). However, in certain geographical areas such as Italy, India and Brazil, it may account for over 10 % of all candidaemia (Sandven, 2000Down).

Unlike the more common C. albicans, whose primary ecological niche is on or in mammals and for which a complete sexual life cycle has not been found, C. guilliermondii can be found in a variety of ecological niches and has an identifiable sexual life cycle in the laboratory (Odds, 1988Down; Barnett et al., 2000Down). However, whether such distinct lifestyles contribute to different population genetic patterns remains unknown. Before the mid-1990s, the identification of human pathogenic (as well as non-pathogenic) fungi was mostly based on morphological and physiological features. The use of molecular markers in the last decade has uncovered several previously unknown species such as Candida dubliniensis, a close relative of C. albicans (Sullivan et al., 1995Down), and two close relatives of C. parapsilosis, C. orthopsilosis and C. metapsilosis (Tavanti et al., 2005Down). In clinical samples of C. guilliermondii, multilocus enzyme electrophoresis (MLEE) and karyotyping using pulsed-field gel electrophoresis (PFGE) revealed that this species in fact contained two different species, C. fermentati and C. guilliermondii (Bai et al., 2000Down; San Millan et al., 1997Down). However, readily applicable diagnostic markers to distinguish these two species remain elusive. Furthermore, the patterns of DNA sequence variation within and between these two closely related species remain unknown.

The objective of this study was to use the multiple gene genealogy (also called multilocus sequence typing or MLST; Maiden et al., 1998Down) approach to understand the evolution and population genetic structure of C. guilliermondii and C. fermentati. Gene genealogical analyses have been used to address a variety of ecological and evolutionary genetic issues, including the rates of dispersal, hybridization and divergence, and the role of clonality and recombination in the structure of fungal populations (Xu, 2006Down; Xu et al., 2002Down, 2005Down; Xu & Mitchell, 2003bDown). Given a group of clonally reproducing organisms, the genealogies of multiple genes are expected to be the same, since evolutionary changes arise mainly through independent mutations. In contrast, in a group of organisms in which sexual recombination has occurred to some degree, the gene genealogies may be expected to differ, since sexual mating and meiosis result in the shuffling of genes from different parts of a genome (Xu, 2004bDown).

This study was designed to address the following specific questions. First, how much divergence is there between C. guilliermondii and C. fermentati at the DNA sequence level? Will the qualitative differences between the two species as revealed by MLEE and karyotyping be substantiated at the DNA sequence level? Second, are genealogies from different genes congruent? Congruent gene genealogies suggest a predominantly clonal mode of reproduction in natural populations of the species. Because C. guilliermondii is capable of sexual reproduction in the laboratory, we hypothesize that there should be evidence for gene genealogy incongruence and recombination in natural populations of C. guilliermondii. Third, what is the geographical pattern of molecular variation? Specifically, are strains from the same area more similar to each other than to those from other areas? And fourth, are there widespread clones or clonal lineages in C. guilliermondii as has been found for several other human fungal pathogens (Pujol et al., 2005Down; Xu, 2004aDown)? To address these questions, we analysed a collection of 37 strains using DNA sequences from fragments of five genes.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Yeast isolates.
Thirty-seven isolates were obtained from diverse geographical areas in four countries. Three strains were isolated from two healthy hosts in an aboriginal community at the First Nations Reserve just west of Hamilton, Ontario, Canada. Twelve strains were obtained from the regional Clinical Microbiology Laboratory in Toronto, Ontario. All 12 strains were from unrelated patients with yeast infections by C. guilliermondii. Six strains were collected from separate patients in the Philippines. Two strains were isolated from the oral cavities of patients infected with HIV in Tanzania. Fourteen strains were collected from three provinces in China: eight strains from a small farming village (Shitou Village) in Jiangxi Province in southeast China; one strain from Xi'an, the capital of Shaanxi Province in central China; and five strains from Chengdu, the capital of Sichuan Province in southwest China.

These 37 strains were initially identified as C. guilliermondii based on standard clinical mycological procedures using growth phenotypes on agar plates, microscopic features and API32 substrate utilization profiles (Warren & Hazen, 1999Down; Kam & Xu, 2002Down). Their identifications were confirmed and further distinctions were made using sequence analysis of the internal transcribed spacer (ITS) regions 1 and 2 using the PCR primers ITS1 and ITS4. These 37 strains were found to belong to two different species, 32 being C. guilliermondii and 5 being C. fermentati (see below). Details of these strains and their species designations are presented in Table 1Down. Except for strains LM85 and LM86, which were obtained from the same host but different body sites, all strains were isolated from a distinct host/source (Table 1Down).


View this table:
[in this window]
[in a new window]
 
Table 1. Strains of Candida guilliermondii and Candida fermentati used in this study

 
DNA manipulation.
High-molecular-mass genomic DNA was isolated using a previously described technique (Xu et al., 2000aDown). Fragments of five genes were studied: ACT1 (actin 1 gene), EF1 (translation elongation factor 1), ITS of the rRNA gene cluster, RIBO (riboflavin synthetase) and TOPOII (DNA topoisomerase II). These genes were chosen because they were conserved, housekeeping genes, to ensure successful amplification and sequencing of gene fragments from all strains. The chromosomal locations, primer sequences and original GenBank accession numbers for these genes in C. guilliermondii are given in Table 2Down. These GenBank sequences were used to design specific primers for the amplification and sequencing of each gene from each of the 37 strains. The primers and the expected PCR product size for all five gene fragments are presented in Table 2Down. PCR amplification was performed according to standard protocols (Xu et al., 2002Down, 2000bDown; Xu & Mitchell, 2003bDown). PCR products were cleaned using the microCLEAN purification kit (Qiagen) and sequenced using both the forward and reverse primers using an ABI3100 automated DNA sequencer (MoBix Laboratory, McMaster University). For each of the five gene fragments, strains that showed DNA sequences differing by only one or a few nucleotides from the dominant alleles were sequenced again to confirm the observed base substitutions. In all cases, the second round of sequencing produced sequences identical to those obtained in the first round.


View this table:
[in this window]
[in a new window]
 
Table 2. Gene and primer information for the five genes analysed in this study

 
Data analyses.
Phylogenetic analyses of the five individual gene fragments as well as the combined data were performed using the PAUP* software (Swofford, 2004Down). The most parsimonious trees were obtained by heuristic searches based on 500 random sequence additions. Mid-point rooting was used for all phylogenetic trees.

To test for clonality and recombination in samples of the two species, we implemented three tests. In the first, we calculated the overall index of association (IA). In this test, the observed data were compared against the null hypothesis that alleles from different loci were randomly associating with each other (i.e. a recombining population structure). A clonal population structure would be suggested if the null hypothesis were rejected; i.e. if alleles showed significant association with each other (Xu, 2005Down). In the second test, the proportion of pairwise loci that were phylogenetically incompatible was calculated. A phylogenetic incompatibility occurs between two loci with two alleles each when all four possible genotypes are found in the population. To perform the above two tests, each unique haplotype of each gene was counted as a unique allele and the data were then imported to the program MULTILOCUS version 1.3 (Agapow & Burt, 2001Down) for analyses. The underlying assumptions, formulae, and the inferences of statistical significance of these two tests can be found on the program homepage (Agapow & Burt, 2001Down). In the third test, the congruence of gene genealogies was assessed using the partition homogeneity (PH) test (Farris et al., 1995Down). In the PH test, the null hypothesis is congruence and clonality, the opposite of that for the index of association test (Xu, 2005Down). These three complementary tests were applied to samples of both C. guilliermondii and C. fermentati individually.

To test whether there were significant phylogenetic patterns based on geographical origins, we used the topology-dependent permutation tail probability (T-PTP) test (Faith, 1991Down). This test compares the lengths of maximum-parsimony (MP) trees with and without the monophyletic constraint defined by the described trait, geography. If the constrained trees are significantly longer than the unconstrained MP tree, the results suggest a lack of phylogenetic pattern based on the specific trait. The statistical significance of the T-PTP test was derived from permutation of the sequence data under the assumption of non-monophyly to generate a null distribution of tree lengths. Statistical support for non-monophyly is achieved when over 95 % of all permuted datasets have tree lengths shorter than the MP tree generated with the constraint of monophyly (Faith, 1991Down). The T-PTP test was implemented in PAUP* (Swofford, 2004Down).

Based on the derived phylogenetic pattern, the existence of potentially divergent lineages within each of the two species was further analysed using evolutionary distances within and between groups of strains. In this study, the common Kimura two-parameter (K2P) distance implemented in PAUP* was used. The K2P model treats transitions and transversions differently and uses observed empirical substitution patterns to derive the optimum weighing schemes among various types of transitions and transversions (Swofford, 2004Down).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We successfully obtained DNA sequences from fragments of all five genes for all 37 strains. For each of the 185 (37x5) DNA sequences, we observed no evidence of heterozygosity (i.e. double peaks in the sequencing chromatograms). This observation is consistent with the haploid nature of the C. guilliermondii and C. fermentati genomes. In the following analyses, we therefore treat the strains as haploid. The composite multilocus genotype for each strain is presented in Table 1Up. The patterns of DNA sequence variation are described below.

Significant molecular divergence between C. guilliermondii and C. fermentati
The BLASTn analyses of the ITS gene sequences identified that the 37 strains analysed here belonged to two different species. Thirty-two strains had ITS sequences 99·7–100 % identical to that of the type sequence of C. guilliermondii, and five strains matched 100 % of the type sequence of C. fermentati. Within the ITS1 and ITS2 regions, four polymorphic nucleotide sites distinguish these two species (Fig. 1Down).


Figure 1
Figure 1
View larger version (17K):
[in this window]
[in a new window]
 
Fig. 1. Maximum-parsimony trees for 37 strains of C. guilliermondii and C. fermentati from each of five gene fragments sequenced. (a) ACT1, tree length=14, consistency index (CI)=1·0; (b) EF1, tree length=20, CI=1·0; (c) ITS, tree length=7, CI=1·0; (d) RIBO, tree length=96, CI=0·938; (e) TOPOII, tree length=63, CI=0·984. The branch lengths are proportional to the amount of sequence divergence within each gene tree. The number of nucleotide changes for each branch is presented. Note the differences among genes in the scale bars: 0·5 change for (a) and (b), 0·1 change for (c), 1 change for (d) and (e).

 
The other four gene fragments analysed here also showed consistent divergence between C. guilliermondii and C. fermentati (Fig. 1Up, Table 3Down). These results contrast with those from traditional identification methods using morphological characteristics and physiological features, which failed to distinguish strains of these two species. Among the five gene fragments, RIBO showed the greatest divergence between the two species, followed by TOPOII, EF1, ACT1 and the least divergent ITS gene. The mean K2P distance between pairs of strains for RIBO (0·1413) was about twice that of TOPOII and almost 20 times that of the ITS gene fragment (Table 3Down). When all five genes are considered together, the mean K2P distance between C. guilliermondii and C. fermentati is 0·0511. This divergence is over 20 times the mean divergence between strains within both species. The results here thus support independent species status for C. guilliermondii and C. fermentati.


View this table:
[in this window]
[in a new window]
 
Table 3. Mean pairwise Kimura two-parameter distances (x10–2) within and between C. guilliermondii and C. fermentati

 
Limited nucleotide variation within populations of C. guilliermondii and C. fermentati
Except the ITS gene fragment, which showed no variation among the five strains of C. fermentati, all other gene fragments analysed here showed various degrees of nucleotide variation within populations of these two species. However, compared to the extensive variation between the two species, the variation within samples of both species was very limited. In addition, the rank orders among the five genes in the degree of nucleotide variation within both species were different from that between the two species. Specifically, within C. guilliermondii, the most divergent was RIBO, followed by EF1, TOPOII, ITS and ACT1. Within C. fermentati, the most divergent was TOPOII, followed by RIBO, ACT1, EF1 and ITS. These results indicated that the most divergent genes within these two species were not necessarily the most divergent between the species and vice versa.

The number of haplotypes among genes also varied in populations of the two species and was not directly correlated with the degree of nucleotide variation within or between the two examined species. Specifically, in C. guilliermondii, two haplotypes were found for ACT1, three for EF1, four each for TOPOII and ITS, and six for RIBO (Table 1Up). In C. fermentati, one haplotype was found for ITS, two for EF1, and three each for ACT1, RIBO and TOPOII (Table 1Up). The combined nucleotide data identified a total of 10 multilocus genotypes among the 32 strains of C. guilliermondii and 4 multilocus genotypes among the five strains of C. fermentati.

Evidence for clonal reproduction in C. guilliermondii
Several lines of evidence suggest that clonal reproduction is the main mode of reproduction in natural populations of C. guilliermondii. First, we observed that several multilocus genotypes were shared by multiple strains, all from unrelated hosts. Specifically, among 10 multilocus genotypes identified for the 32 strains of C. guilliermondii, three were shared by a total of 25 strains: one genotype was shared by three strains from the Philippines, one by 12 strains from Ontario, Canada, and the third by 10 strains from diverse geographical locations in China (Table 1Up, Fig. 2Down). The second line of evidence came from gene genealogical comparisons. The PH test revealed that genealogies from all five genes are congruent among each other in the sample of C. guilliermondii (P=0·81). The third line of evidence came from linkage disequilibrium analyses. The alleles from the five analysed loci showed no evidence of random association (P<0·001, Table 4Down). Indeed, we found no evidence for phylogenetic incompatibility in this sample of C. guilliermondii (P<0·001, Table 4Down).


Figure 2
View larger version (8K):
[in this window]
[in a new window]
 
Fig. 2. One of 28 maximum-parsimony trees from the combined DNA sequences of the ACT1, EF1, ITS, RIBO and TOPOII genes. Tree length=204, CI=0·946; scale bar, 5 changes.

 

View this table:
[in this window]
[in a new window]
 
Table 4. Analyses of association among haplotypes in populations of C. guilliermondii and C.fermentati

 
Similar evidence exists for clonal reproduction in the natural population of C. fermentati. The PH test identified that genealogies from the five genes are statistically congruent among each other in this sample of C. fermentati (P=0·48). No evidence for phylogenetic incompatibility was observed. Furthermore, despite the small sample size, linkage disequilibrium analysis rejected the null hypothesis of random mating (Table 4Up). Among the five strains analysed here, two shared identical multilocus genotypes. These two strains were isolated from the same host, one from the oral cavity and the other from a fingernail. However, it should be pointed out that the small sample sizes analysed here might not be representative of the global population of C. fermentati. Additional samples will be needed before a more robust conclusion can be drawn about the predominant mode of reproduction in natural populations of C. fermentati.

Geographical pattern of molecular variation in C. guilliermondii
Because of the small sample size for C. fermentati, our analysis of the geographical pattern of molecular variation was restricted to C. guilliermondii. As indicated above and shown in Fig. 2Up, three clones were found among the 32 strains of C. guilliermondii. All three clones were restricted to a specific geographical region: one was found only in the Philippines, one only in Ontario, Canada, and the third only in China. These results suggest a geographical pattern for clonal distribution in C. guilliermondii. Interestingly, these three clones were genetically very similar to each other, with only 1–2 nucleotide difference(s) (out of a total 3208 nucleotides) separating them. In addition, five other genotypes with one strain each (four from China and one from Toronto, Canada) differed at only one nucleotide position from the large clone from China. The small number of nucleotide differences suggests a recent common origin for these 30 strains. However, when all 32 strains of C. guilliermondii were analysed together using the T-PTP test, there was no statistically robust geographical pattern (P<0·05). The lack of a strict geographical pattern is evident in the combined phylogenetic tree (Fig. 2Up). Specifically, strains from none of the three regions (Ontario, the Philippines and China) were exclusively clustered together with each clade containing only strains from one specific region but not from other regions. Taken together, the evidence from the above analyses suggests the existence of dispersion among geographical populations of C. guilliermondii.

Are there additional divergent lineages in C. guilliermondii?
While strains PH144 and PH145 from the Philippines had ITS sequences indistinguishable from the majority of strains from China, their sequences in the ACT1, EF1, RIBO and TOPOII genes were consistently different from the other 30 C. guilliermondii strains (Table 1Up, Figs 1 and 2UpUp). The combined analysis of all five genes revealed that the mean K2P distance of the two strains from the other 30 strains of C. guilliermondii (K2P distance 0·01274±0·00035; mean±SD, n=60 pairwise comparisons) was about 30 times the divergence among strains within the main C. guilliermondii clade (0·00042±0·00031, n=435 pairwise comparisons). These results suggest that the current species C. guilliermondii may contain additional evolutionarily divergent lineages.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
This is believed to be the first multiple gene genealogical analysis and comparison of natural strains of C. guilliermondii and C. fermentati. The lack of heterozygosity within strains for any of the five genes in any of the 37 strains is consistent with the haploid nature of the genomes in C. guilliermondii and C. fermentati. Our analyses revealed significant molecular divergence between populations of these two species. While there was no strict geography-based pattern of molecular variation, we identified several region-specific clones of C. guilliermondii. Samples of both species exhibited evidence of clonal reproduction and limited or no evidence for recombination. Interestingly, our results also suggested that there might be additional divergent lineages within C. guilliermondii.

The observed significant molecular divergence between C. guilliermondii and C. fermentati in all five genes analysed suggested that potentially useful molecular markers could be developed to distinguish strains of the two species. Such a marker system could be highly valuable in clinical microbiology laboratories. Based on the consensus DNA sequences for each of the two species, a list of diagnostic restriction polymorphic sites distinguishing C. guilliermondii and C. fermentati was generated (Table 5Down). These polymorphic sites can be assayed through simple PCR reactions using primers listed in Table 2Up, followed by restriction enzyme digests of the PCR products and gel electrophoresis of the digested DNA. Such assays are unambiguous, rapid, economical, and ideal for clinical application.


View this table:
[in this window]
[in a new window]
 
Table 5. Examples of restriction site polymorphisms in fragments of five genes that distinguish C. guilliermondii and C. fermentati

For each gene, the restriction site polymorphisms were derived from consensus sequences of 32 strains of C. guilliermondii and 5 strains of C. fermentati. Only examples of restriction enzymes that cut only once or not at all are given here for each gene.

 
Because C. guilliermondii is capable of sexual reproduction in the laboratory, we hypothesized that there should be evidence for genealogical incongruence and recombination in natural populations of C. guilliermondii. Contrary to our expectation, we found no evidence for recombination in either C. guilliermondii or C. fermentati populations analysed here. These results contrast with those from the common yeast C. albicans, for which evidence for recombination in natural populations is often found despite our failure to identify its complete sexual reproductive cycle in the laboratory (Pujol et al., 2005Down; Xu et al., 1999cDown). It should be pointed out that the lack of evidence for recombination in the samples of C. guilliermondii and C. fermentati was not due to small sample sizes per se. In our analyses of clonality and recombination, three tests with two opposite null models, strict clonality (implemented in the incompatibility ratio test and the PH test) and random mating (implemented in the linkage equilibrium or IA test) were conducted and compared. In population genetic analyses, small sample size typically results in the lack of power to reject the null hypothesis. The fact that the null hypothesis of random mating in the index of association test was rejected in samples of both species suggests unambiguous evidence for clonality in the analysed samples. However, it should be emphasized that although the analysed samples are consistent with clonal reproduction in these two species, for several reasons, our data cannot exclude the possibility of recombination in either species.

First, although our samples were collected from diverse geographical regions without any a priori set of criteria, the samples analysed here were relatively small and may not be representative of the global populations of the two species. It is possible that samples from other geographical areas or natural environments may exhibit evidence of recombination. Analysis of additional samples is needed to identify the potential existence of such geographical or ecological-niche-specific population structures in these two species.

Second, the overall level of molecular variation is relatively limited within each of the two species. As a result, the two tests that used strict clonality as the null hypotheses may have lacked the power to reject the null hypothesis. Additional analyses using more polymorphic markers such as microsatellite loci may offer greater power to reject the hypothesis of strict clonality (Xu, 2005Down). Such markers should also allow a more critical examination of potentially divergent lineages within C. guilliermondii. However, it should be emphasized that, despite the low level of nucleotide variation within the two samples analysed here, the null hypothesis of random mating was rejected for both in the test of linkage equilibrium.

Third, the five marker genes may be functionally linked and thus significantly associated. Physically, these five genes are located on four different chromosomes (or supercontigs; Table 2Up) as shown in the recent assembly of a 10x coverage genome sequencing of a C. guilliermondii strain by the Broad Institute http://www.broad.mit.edu/annotation/fungi/candida_guilliermondii/background.html). Only two genes, ITS and TOPOII, are located on the same chromosome (supercontig 5). This supercontig is about 1·2 million base pairs long and ITS and TOPOII are located about 860 kb away from each other, almost at opposite ends of the chromosome. At present, the physical locations of the five genes are unknown in the C. fermentati genome. However, although physically these genes are far apart from each other in the C. guilliermondii genome, it is possible that they are functionally linked and interact with each other epistatically to influence fitness of strains. Epistatic interactions could allow the build-up of linkage disequilibrium in natural populations (Xu, 2004bDown). Analysing different marker genes should allow a direct test of this possibility.

Lastly, because of the limited genetic variation within each of the two samples, it is highly likely that even if sexual reproduction were common in populations of C. guilliermondii and C. fermentati, the analyses conducted here might not have been able to discriminate it from clonal reproduction. For example, two forms of sexual reproduction, crossing between close relatives (i.e. inbreeding) and self-fertilization, can both generate progeny with identical genotypes that would be recognized as belonging to the same asexual clones (Xu, 2005Down; Xu et al., 2005Down).

We found geographically specific clones of C. guilliermondii in three locations (China, the Philippines, and Ontario, Canada). This result suggests that there might be large locally adapted clones in this species. The Ontario clone contained strains from the metropolitan Toronto area as well as from a small native community about 70 km west of Toronto, two economically and ecologically very different areas. The clone from China contained strains from three locations separated by hundreds to thousands of kilometres: a village in southeast China (Shitou Village, Jiangxi Province), a large city in central China (Xi'an, Shaanxi Province), and another large city in southwest China (Chengdu, Sichuan Province). Furthermore, the close similarity in the DNA sequences among the three clones as well as the other five strains (TOR13, ZG57, ZG70, ZG248, ZG379; Fig. 2Up) suggests a recent common origin of the 30 C. guilliermondii strains (except strains PH144 and PH145).

Like most micro-organisms that show wide geographical distribution of certain genotypes but where the mechanisms for such patterns remain unknown (e.g. Kidd et al., 2005Down; Xu, 2004aDown; Xu et al., 1999bDown, 2005Down), the reasons for the wide distribution of some of the genotypes in C. guilliermondii are also unknown at present. However, there are several possibilities. First, human travelling in the past several decades might have contributed to the spread of certain genotypes. Second, the populations of C. guilliermondii from Asia and North America might have been the result of a severe population bottleneck associated with prehistoric humans when they migrated from Africa to other parts of the world a few hundred thousand years ago. Third, unlike C. albicans, for which mammals are the primary host, C. guilliermondii can be found in diverse environments such as water and soil, as well as a variety of animals (Odds, 1988Down; Barnett et al., 2000Down). Therefore, there could be dispersal through these environmental sources. All three possibilities could have contributed to clonal dispersal in C. guilliermondii.

The identification of additional divergent lineages within the current species of C. guilliermondii further supports the notion that clinical strains previously identified as C. guilliermondii consist of a group of closely related phylogenetic species. At present, the medical significance of the species and divergent clones is unknown. However, using the information obtained here, and with the availability of a genome sequence (http://www.broad.mit.edu/annotation/fungi/candida_guilliermondii/background.html), the analysis of additional clinical and environmental strains of C. guilliermondii using these and additional genetic markers should establish a solid evolutionary framework to allow a comprehensive examination of the epidemiology, pathogenesis, and control and treatment of C. guilliermondii infections.


    ACKNOWLEDGEMENTS
 
We thank many people who have contributed to the strain collections: Dr Glen Bulmer for the six strains from the Philippines, Laura Montour for the three strains from the First Nations Reserve in Ontario, Dr Debby Yamamura for the 12 strains from Toronto, Ontario, and Dr Wiley Schell of Duke University for the two strains from Tanzania. We thank Heather Yoell for comments. This work was funded by the Premier's Research Excellence Award, Genome Canada, and the Natural Science and Engineering Research Council (NSERC) of Canada.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Agapow, P. M. & Burt, A. (2001). Indices of multilocus linkage disequilibrium. Mol Ecol Notes 1, 101–102.[CrossRef]

Bai, F. Y., Liang, H. Y. & Jia, J. H. (2000). Taxonomic relationships among the taxa in the Candida guilliermondii complex, as revealed by comparative electrophoretic karyotyping. Int J Syst Evol Microbiol 50, 417–422.[Abstract]

Barnett, J. A., Payne, R. W. & Yarrow, D. (2000). Yeast Characteristics and Identification, 3rd edn. Cambridge: Cambridge University Press.

Doi, M., Homma, M., Chindamporn, A. & Tanaka, K. (1992). Estimation of chromosome number and size by pulsed-field gel electrophoresis (PFGE) in medically important Candida species. J Gen Microbiol 138, 2243–2251.[Medline]

Dos Santos, A. L. & Soares, R. M. (2005). Candida guilliermondii isolated from HIV-infected human secretes a 50 kDa serine proteinase that cleaves a broad spectrum of proteinaceous substrates. FEMS Immunol Med Microbiol 43, 13–20.[Medline]

Faith, D. L. (1991). Cladistic permutation tests for monophyly and nonmonophyly. Syst Zool 40, 366–375.[CrossRef]

Farris, J. S., Kallersjo, M., Kluge, A. G. & Bult, C. (1995). Testing significance of incongruence. Cladistics 10, 315–319.[CrossRef]

Imshenetskii, A. A., Solntseva, L. I. & Afanas'eva, L. N. (1979). Biomass accumulation by polyploid forms of Candida guilliermondii yeasts. Mikrobiologiia 48, 418–422.[Medline]

Kam, A. & Xu, J. (2002). Diversity of commensal yeasts from within and between healthy hosts. Diag Microbiol Infect Dis 43, 19–28.[CrossRef][Medline]

Kidd, S. E., Guo, H., Bartlett, K., Xu, J. & Kronstad, J. W. (2005). Comparative gene genealogies indicate that two clonal lineages of Cryptococcus gattii in British Columbia resemble strains from other geographical areas. Eukaryot Cell 4, 1629–1638.[Abstract/Free Full Text]

Krcmery, V. & Barnes, A. J. (2002). Non-Candida albicans spp. causing fungaemia: pathogenicity and antifungal resistance. J Hosp Infect 50, 243–260.[CrossRef][Medline]

Maiden, M. C., Bygraves, J. A., Feil, E. & 10 other authors (1998). Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95, 3140–3145.[Abstract/Free Full Text]

McEwen, J. G., Taylor, J. W., Carter, D. & 8 other authors (2000). Molecular typing of pathogenic fungi. Med Mycol 38, 189–197.

Odds, F. C. (1988). Candida and Candidiasis. London: Baillière Tindall.

Pujol, C., Dodgson, A. & Soll, D. R. (2005). Population genetics of ascomycetes pathogenic to humans and animals. In Evolutionary Genetics of Fungi, pp. 149–188. Edited by J. Xu. Norwich, UK: Horizon Scientific Press.

Sandven, P. (2000). Epidemiology of candidemia. Rev Iberoam Micol 17, 73–81.[Medline]

San Millan, R. M., Wu, L. C., Salkin, I. F. & Lehmann, P. F. (1997). Clinical isolates of Candida guilliermondii include Candida fermentati. Int J Syst Bacteriol 47, 385–393.[Abstract/Free Full Text]

Sullivan, D. J., Westerneng, T. J., Haynes, K. A., Bennett, D. E. & Coleman, D. C. (1995). Candida dubliniensis sp. nov.: phenotypic and molecular characterization of a novel species associated with oral candidosis in HIV-infected individuals. Microbiology 141, 1507–1521.[Abstract]

Swofford, D. L. (2004). PAUP* – Phylogenetic Analysis Using Parsimony*and Other Methods, version 4. Sunderland, MA: Sinauer Associates.

Tavanti, A., Davidson, A. D., Gow, N. A., Maiden, M. C. & Odds, F. C. (2005). Candida orthopsilosis and Candida metapsilosis spp. nov. to replace Candida parapsilosis groups II and III. J Clin Microbiol 43, 284–292.[Abstract/Free Full Text]

Tey, R., Yan, Z., Han, S., Li, X., Lazazzera, K., Sun, S. & Xu, J. (2003). Ecology and epidemiology of drug resistance in human pathogenic yeasts. Res Adv Microbiol 3, 67–85.

Warren, N. G. & Hazen, K. C. (1999). Candida, Cryptococcus, and other yeasts of medical significance. In Manual of Clinical Microbiology, 7th edn, pp. 1184–1199. Edited by P. R. Murray, E. J. Baron, M. A. Pfaller, F. C. Tenover & R. H. Yolken. Washington, DC: American Society for Microbiology.

Xu, J. (2004a). Gene genealogical analyses of human fungal pathogens. In Pathogenic Fungi: Structural Biology and Taxonomy, pp. 299–318. Edited by G. San-Blas & R. Calderone. Norwich: Horizon Scientific Press.

Xu, J. (2004b). The prevalence and evolution of sex in microorganisms. Genome 47, 775–780.[Medline]

Xu, J. (2005). Fundamentals of fungal molecular population genetic analyses. Evolutionary Genetics of Fungi, pp. 87–116. Edited by J. Xu. Norwich: Horizon Scientific Press.

Xu, J. (2006). Microbial ecology in the age of genomics and metagenomics: concepts, tools and recent advances. Mol Ecol (in press).

Xu, J. & Mitchell, T. G. (2003a). Geographic differences in human oral yeast flora. Clin Infect Dis 36, 221–224.[Medline]

Xu, J. & Mitchell, T. G. (2003b). Comparative gene genealogical analyses of strains of serotype AD identify recombination in populations of serotypes A and D in the human pathogenic yeast Cryptococcus neoformans. Microbiology 149, 2147–2154.[Abstract/Free Full Text]

Xu, J., Boyd, C. M., Livingstone, E., Meyer, W., Madden, J. F. & Mitchell, T. G. (1999a). Species and genotypic diversities and similarities of pathogenic yeasts colonizing women. J Clin Microbiol 37, 3835–3843.[Abstract/Free Full Text]

Xu, J., Vilgalys, R. & Mitchell, T. G. (1999b). Lack of genetic differentiation between two geographic samples of Candida albicans isolated from patients infected with the human immunodeficiency virus. J Bacteriol 181, 1369–1373.[Abstract/Free Full Text]

Xu, J., Mitchell, T. G. & Vilgalys, R. (1999c). PCR-RFLP analyses reveal both extensive clonality and local genetic differentiation in Candida albicans. Mol Ecol 8, 59–74.[CrossRef][Medline]

Xu, J., Ramos, A. R., Vilgalys, R. & Mitchell, T. G. (2000a). Clonal and spontaneous origins of fluconazole resistance in Candida albicans. J Clin Microbiol 38, 1214–1220.[Abstract/Free Full Text]

Xu, J., Vilgalys, R. & Mitchell, T. G. (2000b). Multiple gene genealogies reveal recent dispersion and hybridization in the human pathogenic fungus Cryptococcus neoformans. Mol Ecol 9, 1471–1481.[CrossRef][Medline]

Xu, J., Luo, G., Vilgalys, R. J., Brandt, M. E. & Mitchell, T. G. (2002). Multiple origins of hybrid strains of Cryptococcus neoformans with serotype AD. Microbiology 148, 203–212.[Abstract/Free Full Text]

Xu, J., Cheng, M., Tan, Q. & Pan, Y. (2005). Molecular population genetics of basidiomycete fungi. In Evolutionary Genetics of Fungi, pp. 221–252. Edited by J. Xu. Norwich, UK: Horizon Scientific Press.

Received 24 October 2005; revised 17 January 2006; accepted 18 January 2006.


This article has been cited by other articles:


Home page
Eukaryot CellHome page
F. C. Odds and M. D. Jacobsen
Multilocus Sequence Typing of Pathogenic Candida Species
Eukaryot. Cell, July 1, 2008; 7(7): 1075 - 1084.
[Full Text] [PDF]


Home page
GeneticsHome page
S. Sun and J. Xu
Genetic Analyses of a Hybrid Cross Between Serotypes A and D Strains of the Human Pathogenic Fungus Cryptococcus neoformans
Genetics, November 1, 2007; 177(3): 1475 - 1486.
[Abstract] [Full Text] [PDF]


Home page
MicrobiologyHome page
J. Xu, H. Guo, and Z.-L. Yang
Single nucleotide polymorphisms in the ectomycorrhizal mushroom Tricholoma matsutake
Microbiology, July 1, 2007; 153(7): 2002 - 2012.
[Abstract] [Full Text] [PDF]


Home page
MicrobiologyHome page
S. Sun, H. Guo, and J. Xu
Multiple gene genealogical analyses reveal both common and distinct population genetic patterns among replicons in the nitrogen-fixing bacterium Sinorhizobium meliloti.
Microbiology, November 1, 2006; 152(Pt 11): 3245 - 3259.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lan, L.
Right arrow Articles by Xu, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lan, L.
Right arrow Articles by Xu, J.
Agricola
Right arrow Articles by Lan, L.
Right arrow Articles by Xu, J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS
Copyright © 2006 Society for General Microbiology.