Next Article in Journal
Yeast Bloodstream Infections in the COVID-19 Patient: A Multicenter Italian Study (FiCoV Study)
Next Article in Special Issue
Innate Pulmonary Phagocytes and Their Interactions with Pathogenic Cryptococcus Species
Previous Article in Journal
The Antidepressant Sertraline Affects Cell Signaling and Metabolism in Trichophyton rubrum
Previous Article in Special Issue
Cryptococcus neoformans Causing Meningoencephalitis in Adults and a Child from Lima, Peru: Genotypic Diversity and Antifungal Susceptibility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Analyses of Multi-Locus Sequence Typing Data Reveal Geographic Differentiation, Hybridization, and Recombination in the Cryptococcus gattii Species Complex

Department of Biology, McMaster University, Hamilton, ON L8S 4K1, Canada
*
Author to whom correspondence should be addressed.
J. Fungi 2023, 9(2), 276; https://doi.org/10.3390/jof9020276
Submission received: 31 January 2023 / Revised: 16 February 2023 / Accepted: 18 February 2023 / Published: 20 February 2023

Abstract

:
Cryptococcus gattii species complex (CGSC) is a basidiomycete haploid yeast and globally distributed mammalian pathogen. CGSC is comprised of six distinct lineages (VGI, VGII, VGIII, VGIV, VGV, and VGVI); however, the geographical distribution and population structure of these lineages is incompletely described. In this study, we analyze published multi-locus sequence data at seven loci for 566 previously recorded sequence types (STs) encompassing four distinct lineages (VGI, VGII, VGIII, and VGIV) within the CGSC. We investigate indicators of both clonal dispersal and recombination. Population genetic analyses of the 375 STs representing 1202 isolates with geographic information and 188 STs representing 788 isolates with ecological source data suggested historically differentiated geographic populations with infrequent long-distance gene flow. Phylogenetic analyses of sequences at the individual locus and of the concatenated sequences at all seven loci among all 566 STs revealed distinct clusters largely congruent with four major distinct lineages. However, 23 of the 566 STs (4%) each contained alleles at the seven loci belonging to two or more lineages, consistent with their hybrid origins among lineages. Within each of the four major lineages, phylogenetic incompatibility analyses revealed evidence for recombination. However, linkage disequilibrium analyses rejected the hypothesis of random recombination across all samples. Together, our results suggest evidence for historical geographical differentiation, sexual recombination, hybridization, and both long-distance and localized clonal expansion in the global CGSC population.

1. Introduction

Cryptococcosis is an infection caused by the inhalation of the human pathogenic Cryptococcus (HPC) from the environment. The opportunistic haploid yeast initially establishes within the lungs, and without treatment or in the case of long, asymptomatic latent periods, it can enter the bloodstream and into the central nervous system, leading to meningoencephalitis and potentially death [1]. Due to the enormous global burden, especially among patients with HIV infection [2], HPC was recently classified as a priority fungal pathogen by the World Health Organization (WHO). HPC is commonly found in nature, including in bird droppings, trees, and soils, all serving as environmental reservoirs of infectious propagules for humans [3,4,5,6]. HPC consists of two sister species complexes: the Cryptococcus gattii species complex (CGSC) and the Cryptococcus neoformans species complex (CNSC) [7].
For over 100 years, CGSC was referred to as C. neoformans var. gattii and treated as a variety of C. neoformans [8]. Because of this, C. gattii has often been overlooked and is most likely underrepresented in the available historical literature. It was Dr. June Kwon-Chung who emphasized the difference between isolates of C. neoformans and C. gattii and presented that C. gattii was predominantly found in tropical and subtropical climates, while C. neoformans was shown to be distributed globally [9,10]. Consistent with her observations, most of the early data regarding C. gattii came from regions of Australia and Papua New Guinea [11,12,13,14,15]. Early findings also suggested that C. gattii favored immunocompetent hosts, contrasting C. neoformans, which predominantly infects immunocompromised hosts [10].
C. gattii began attracting broader attention after an unusual outbreak started in the early 2000s on Vancouver Island, a large island off the coast of British Columbia, Canada. The outbreak has since spread to mainland BC and the Pacific Northwest in the U.S. [16,17]. Indeed, the frequency of C. gattii infection and colonization has increased dramatically in humans, animals, and the environment across these temperate regions [17,18,19,20]. Previously linked to eucalyptus trees [20,21,22], C. gattii has been isolated from soil, air, water, and >10 tree species native to Canada’s west coast [19]. In addition, the infection rate on Vancouver Island was found to be ~9 times higher than in the endemic region of Australia [23]. The proliferation of C. gattii in the temperate region as well as an increase in human populations have generated several hypotheses about its emergence and mechanism of dispersal. Surveys of air samples have found particles indictive of C. gattii spores [23], likely explaining its increasingly broad dispersals and high virulence, as spores are known to be more infectious than yeast cells [24]. In addition, basidiospores are produced from sexual reproduction, suggesting high levels of fertility among Vancouver Island isolates [25,26]. However, sampling also revealed isolates were nearly exclusively of mating type α [26], suggesting that opposite-sex mating for C. gattii in nature is likely uncommon [25]. In contrast, haploid fruiting and same-sex mating have been observed in C. gattii, both of which could have contributed to basidiospore production [27,28,29]. Indeed, evidence for recombination in nature has been reported for molecular type VGII within C. gattii [30]. This ability for mating and potential for meiotic recombination could also contribute to its adaptability, range expansion, and virulence. Together, these features of C. gattii call for greater efforts to understanding the epidemiology, pathogenesis, and evolution.
Although both CNSC and CGSC have a well-defined sexual phase [25], historically, it has been viewed as a predominantly clonal pathogen, with little evidence for sexual recombination in nature [31,32]. However, after fertile populations were reported on Vancouver Island, multiple studies were published regarding the mode of reproduction, where evidence for both clonal reproduction [33,34,35] and sexually recombination were present among defined populations in nature [36,37,38]. At present, the extent of recombination in the global populations of C. gattii remains unknown.
Since the 1990s, a variety of molecular methods has been used to identify the genotypes of strains of the human pathogenic Cryptococcus, including PCR fingerprinting [8], amplified fragment length polymorphism (AFLP) [39], multi-locus sequence typing [40], and whole-genome sequencing [41]. These studies have revealed the diversity and heterogeneity among C. gattii isolates and established that C. gattii is a species complex containing multiple distinct lineages. These lineages are frequently referred to as molecular types and increasingly more as species [42]. Formally, the C. gattii species complex is divided into six distinct lineages (VGI, VGII, VGIII, VGIV, VGV, and VGVI) [3,43,44]. In this report, we focus on the four lineages with isolates represented in the International Fungal Multi Locus Sequence Typing (IFMLST) database. Evidence for hybridization has been reported between some of the lineages within CGSC [30] as well as between CGSC and CNSC [45]. The extent of hybridization among the lineages and recombination within individual lineages in the global population of CGSC remains unknown.
To standardize genotyping and enhance data sharing and communication, in 2009, the International Society for Human and Animal Mycology (ISHAM) implemented the multi-locus sequence typing (MLST) scheme for global genotyping of CNSC and CGSC isolates. The proposed scheme recommends the use of partial DNA sequences from seven unlinked nuclear loci (CAP59, GPD1, LAC1, PLB1, SOD1, URA5, and IGS1) [46]. For each locus, all unique DNA sequences are assigned an allele type (AT), which in combination creates an allelic profile that is assigned a sequence type (ST). A publicly accessible database (The International Fungal Multi Locus Sequence Typing, IFMLST) stores the allele and sequence data for CNSC and CGSC [47]. This allows consistent and efficient comparisons across labs. Over the last 14 years, thousands of HPC isolates have been genotyped using the IFMLST scheme and with the data published and stored in the centralized website. The availability of this standardized genotyping method and genotype dataset enables global population genetic analyses of HPC. Indeed, our recent study of the Cryptococcus neoformans species complex revealed evidence for recombination and statistically significant geographic structuring of the two species within CNSC: C. neoformans and C. deneoformans [48].
In this study, we examined the global population structure of both the total CGSC samples and the individual lineages within CGSC by analyzing the published MLST data of this species complex. We extracted information from a total of 1202 non-redundant isolates from 31 published studies. We evaluated the overall geographical pattern of genetic variation as well as investigated potential evidence for recombination within individual lineages and hybridization among the lineages.

2. Materials and Methods

2.1. Strains, Genotypes and Metadata

We carried out a literature search through PubMed for the term “Cryptococcus MLST”. All papers were examined for relevant MLST data associated with isolates of CGSC, including C. bacillisporus, C. gattii, C. deuterogattii, and C. tetragattii, and all reported data regarding CGSC were retrieved and consolidated. All MLST data for strains of CGSC published by November 2022 were collected from the publicly available International Fungal Multi Locus Sequence Typing database. For each unique sequence type, the database stores its allelic profile and DNA sequences across seven loci. A total of 1202 isolates from 31 studies were included in our meta-analyses of population geographic structures. Protocols for DNA extraction, amplification, sequencing, and allelic profiling as well as detailed descriptions of these isolates and their genotype data can be found in the original reports [34,35,37,38,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75].

2.2. Phylogenetic Distribution of Strains and Genotypes

To determine the relationships among STs within CGSC, DNA sequences at the seven MLST loci were concatenated for each sequence type. The concatenated sequences were then imported into MEGA X and aligned through Muscle [76,77] to reconstruct a phylogenetic tree using a neighbor-joining method based on the K2P distance model (Kimura 2-parameter model) [77]. Using the iTOL software [78], previously described lineage classifications from the IFMLST database were added to their associated ST to determine the distributions of the observed lineages in the overall phylogeny. Additionally, geographical and ecological source data for each ST found in the published literature were added to the specific ST to demonstrate their distributions on the phylogenetic tree.

2.3. Phylogenetic Distribution of Alleles

To investigate the partitioning of ATs among lineages, single-locus phylogenetic trees were reconstructed for all seven MLST loci using the DNA sequences for each AT in the IFMLST database. Three of the ATs associated with STs were excluded from the analysis because they did not have corresponding DNA sequences on the MLST database (AT91 at the CAP59 locus) and (AT43 and AT44 at the PLB1 locus). MEGA X software was used to reconstruct phylogenetic trees for the seven loci with neighbor-joining method based on muscle alignment and the K2P distance model (Kimura 2-parameter model) [76,77]. The trees were amended with metadata associated with each ST using the iTOL software [76]. In addition, the lineage data associated with the ST containing each AT was uploaded to the circular single locus trees to identify ATs shared across lineages.

2.4. Population Genetic Analyses

Population genetic analyses were performed to determine if the total CGSC and individual lineages within CGSC were genetically subdivided at the continental and country levels using GenAlEx V.6.5. [79]. Strains were subdivided based on the collected metadata and evaluated within five taxonomy-based samples: the complete CGSC sample, the VGI sample, the VGII sample, the VGIII sample, and the VGIV sample. For each taxonomic sample, both the clone-corrected (CC) and non-clone-corrected (NCC) datasets were analyzed. The CC dataset included only one isolate representing each ST, and the NCC dataset included all isolates retrieved from the literature.
An analysis of molecular variance (AMOVA) was conducted for each of the ten datasets (five taxonomic CC datasets and five taxonomic NCC datasets) at the country and continental levels. Additionally, Wright’s FST values between pairs of geographical populations were determined. To prevent potential biases due to small sample sizes, only subpopulations with >five isolates were included in these analyses. Statistical significance was determined by comparing the observed results with 1000 permutated datasets established with a null hypothesis of no genetic differentiations among geographic populations within each analyzed dataset.
Phylogenetic compatibility and linkage disequilibrium analyses were then conducted using the clone-corrected allelic profiles at the seven loci across 566 STs to test for evidence of recombination within and among individual lineages. Tests were performed across 17 population samples: the total CGSC sample, the four CGSC lineages (VGI, VGII, VGIII, and VGIV) at the global level, the CGSC sample divided by each continent (Africa, Asia, Europe, Oceania, South America, and North America), as well as six population samples of individual lineage by continent combinations for samples with more than 20 recorded STs (Asia VGI, Europe VGI, Oceania VGI, North America VGI, North America VGII, and South America VGII). Information regarding the underlying concepts of these analyses and how they were conducted can be found in the MultiLocus V1.3 manual [80].

3. Results

As of November 2022, there were 566 sequence types (STs) with affiliated DNA sequence data for all seven loci stored on the Cryptococcus MLST database for CGSC. Geographical data were available for 375 STs (66%) representing 1202 isolates collected from 31 published reports along with associated metadata. However, the remaining 191 STs (33%) were not associated with any geographic information. Our population genetic analyses focused on the 375 STs with geographic information. However, for recombination and hybridization analyses, all STs were included. Below, we outline the collected data regarding the 1202 isolates and illustrate the results of our analyses.

3.1. Geographical and Ecological Distributions

The 375 STs with geographic information representing 1202 isolates were subdivided into four previously proposed lineages (VGI, VGII, VGIII, and VGIV) as indicated by original authors. Among these four lineages, VGII was the most prevalent (721 isolates; 60%), followed by VGI, VGIII, and VGIV, containing 225 (19%), 221 (18%), and 34 (3%) isolates, respectively. The CGSC isolates as well as lineages VGI and VGII were distributed across all six continents. However, among the isolates in the MLST database, VGIV was not represented in North America, and VGIII was not recorded in Africa or Asia. In total, the CGSC isolates were distributed across 45 countries. Among the four lineages, VGI was recorded from the most countries (31), followed by VGII (25), VGIII (12), and VGIV (9). Out of the 1202 isolates with geographical data, only 788 isolates (65.5%) included ecological source data, while the remaining 414 isolates (34.5%) were not associated with any source information. Here, we broadly categorize the isolates into three ecological sources: clinical, environmental, and veterinary. Among the 788 isolates with ecological source information, 416 were collected from clinical sources (53%), followed by environmental (259 isolates; 33%) and veterinary (113 isolates; 14%) sources (Supplementary Table S1). The total CGSC, VGII, and VGIII were reported from all three sources, but VGI was not recorded from veterinary sources, and VGIV was not recorded from environmental sources. Further details pertaining to ecological source data regarding the four lineages is provided in Supplementary Table S2 and summarized in Table 1.
The geographical distribution of the non-redundant 1202 CGSC isolates as well as the isolates subdivided into the four lineages are shown in Table 2. These 1202 isolates were recorded across six continents and 45 countries. At the species complex level, most CGSC isolates in the database were from South America (517; 43%), followed by Asia (196; 16%), North America (179; 14.9%), Oceania (159; 13%), Africa (86; 7%), and Europe (65; 5%). Among the lineages, VGII was the most prevalent across South America, Asia, and Africa. However, VGIII was the most prevalent in North America and VGI the most prevalent in Europe. At the country/region level, most isolates were recorded from Brazil (349; 29%), followed by Australia (153; 12.7%), Colombia (145; 12%), the USA (141; 11.7%), and Taiwan (111; 9%). The remaining 38 countries each contained <100 isolates and together were represented by 303 of the isolates in the published literature.
Table 3 summarizes the geographic and ecological distributions of the 375 STs in the published MLST literature for CGSC with associated metadata. Geographically, among the 375 STs, none were reported across all six continents or a combination of five continents. Only three STs (representing 61 isolates) were each reported from four of the six continents, seven STs (representing 254 isolates) were each reported from three of the six continents, and twenty-four STs (representing 226 isolates) were each reported from two of the six continents. The remaining 341 STs (representing 658 isolates) were reported in only one of the continents. Among the continents, the smallest number of STs were reported from Africa (24 STs representing 32 isolates), and the largest was from South America (156 STs representing 331 isolates) (Table 3). Ecologically, among the 375 STs with geographic information, 188 had ecological niche data. Of these 188 STs, 5 (representing 237 isolates) were each reported from all three ecological niches, 17 (representing 205 isolates) were each reported from two of the three ecological niches, and 166 (representing 346 isolates) were each reported from only one ecological niche. The majority, 123 STs (representing 255 isolates), were reported from only clinical sources (Table 3). The detailed geographic associations of all 375 STs at the continental and country/region levels are presented in Supplementary Tables S3 and S4, respectively.

3.2. Geographic AMOVA

The overall objective of our AMOVA analyses was to determine the amount of genetic variation contributed by geographic separations at the continental and country/region levels. Below, we briefly summarize the results. At the continental level analyses, in the non-clone-corrected sample, genetic variations within continents contributed 83%, 79%, 84%, 88%, and 92% of the total observed genetic variations in the total CGSC population, VGI population, VGII population, VGIII population, and VGIV population, respectively. The remaining 17%, 21%, 16%, 12%, and 8% of the genetic variations were due to among continents. The within-continent and among-continent contributions for four of the populations (total CGSC, VGI, VGII, and VGIII) were statistically significant at the p < 0.001 level, and the VGIV population was statistically significant at the p < 0.05 level (Table 4). In the five clone-corrected samples, genetic variations within continents contributed 95%, 85%, 96%, 97%, and 97% of the total observed genetic variations in the entire CGSC population, the VGI population, VGII population, VGIII population, and VGIV population, respectively, all greater percentages than those without clone corrections. The remaining 5%, 15%, 4%, 3%, and 3% of the total genetic variations were attributed to among continents. The among-continent contributions for three of the five populations (total CGSC, VGI, and VGII) were statistically significant at p < 0.001, while for the VGIII and VGIV populations, they were not statistically significant, with p = 0.07 and p = 0.2, respectively (Table 4). The pairwise comparisons between continents for the five taxonomic samples are shown in Table 5. In addition, variable numbers of pairwise country/region level comparisons also showed statistically significant genetic differentiations. The detailed pairwise country/region level comparisons for the total CGSC, VGI, VGII, VGIII, and VGIV samples are presented in Supplementary Tables S5, S6, S7, S8, and S9, respectively.

3.3. Ecological Niche AMOVA

Further AMOVA analyses were conducted to determine the amount of genetic variation contributed by ecological niche separation within five countries/region each with a sample size greater than 100 (Australia, Brazil, Colombia, the USA, and Taiwan). Both the non-clone-corrected and clone-corrected populations were analyzed and divided into clinical, environmental, and veterinary sample populations. Only ecological populations with >9 isolates were considered. Because of the small sample size, the clone-corrected Australia and clone-corrected Taiwan populations were not analyzed. The remaining three clone-corrected sample populations (Brazil, Colombia, and the USA) failed to reject the null hypothesis that ecological niche separations contributed significantly to the total observed genetic variation within each country. However, all five non-clone-corrected population samples rejected the null hypothesis of no genetic differentiation among ecological samples at the p < 0.001 level, representing 22%, 11%, 10%, 13%, and 10% of genetic variation coming from among ecological niches in Australia, Brazil, Colombia, the USA, and Taiwan, respectively. These results are summarized in Table 6, and the pairwise comparisons are demonstrated in Supplementary Table S10.

3.4. DNA Sequence Variation

The allelic profiles of all 566 STs, including the allele type (AT) number at each of the seven MLST loci (CAP59, GPD1, LAC1, IGS1, PLB1, SOD1, and URA5) for each ST were retrieved from the IFMLST database. A summary of the total reported ATs for each of the seven loci as well as the proportion of shared ATs across the four VG lineages is shown in Table 7. The AT numbers among the seven loci ranged from 46 at the PLB1 locus to 136 at the SOD1 locus. In addition, the proportion of alleles shared among the four VG lineages ranged from 1 to 6, with the GPD1 locus containing the most shared alleles (6) and loci CAP59, PLB1, and URA5 each having only one shared allele. The differences in length (of base pair) among ATs at each locus ranged from 2 bp for CAP59 to 82 bp for IGS1. The range of occurrence of each AT within each of the four lineages is shown in Supplementary Table S11.

3.5. Phylogenetic Analysis

Phylogenetic investigations of the concatenated sequences at the seven MLST loci showed four clearly separated major ST groups, largely corresponding to four VG lineages. However, two STs (ST402 and ST351) were grouped inconsistently to the lineages that they were originally assigned in the published reports. In addition, two STs with previously undefined lineages (ST2 and ST71) were clearly grouped into VGI and VGIV, respectively (Figure 1).
Phylogenetic analyses of the single-locus trees were overall consistent with the phylogenetic pattern based on concatenated gene sequences. However, there was several interesting observations. First, 1–2 allele types at each of the seven loci showed great divergence from their expected clusters. These included AT96 at the CAP59 locus, AT63 and AT64 at the GPD1 locus, AT64 at the LAC1 locus, AT46 and AT47 at the PLB1 locus, AT113 at the IGS1 locus, AT133 and AT134 at the SOD1 locus, and AT59 at the URA5 locus. Single-locus phylogenetic trees including all ATs can be seen in Figure 2 (for locus GPD1) and Supplementary Figures S1–S6 (for other six loci). ATs across five of the seven single-locus trees (CAP59, LAC1, PLB1 IGS1, and SOD1) showed no inconsistencies between expected and observed VG lineage clustering. However, the GPD1 and URA5 single-locus trees contain one and two instances of inconsistent VG lineage grouping. At the GPD1 locus, AT52 (as part of ST400 and ST409) was expected to be grouped within the VGII cluster based on their original assignment in the database but instead was placed within the VGI cluster. At the URA5 locus, AT35 (as part of ST336) and AT34 (as part of ST419) were both grouped within the VGII cluster but were expected to be within the VGIII and VGIV clusters, respectively.

3.6. Inconsistent Allele Type Clustering

Interestingly, each of the seven single-locus trees included at least one occurrence of shared alleles among VG lineages representing potential signatures of recombination or loss of heterozygosity after hybridization. In total, of the 566 STs, 23 (4%) contained an unexpected allele shared with a VG lineage not associated with their originally proposed lineage. Interestingly, ST402 and ST351, which were inconsistently grouped in Figure 2, include six and four shared alleles, respectively. Out of the consistently grouped STs, two (ST431 and ST524) contained two shared alleles each, involving a total of three VG lineages. The remaining 19 STs each contained at least one allele shared between two VG lineages. The complete allelic profiles of these 23 STs are displayed in Table 8.

3.7. Recombination and Linkage Disequilibrium

We examined potential indicators of recombination by investigating two standard signatures: phylogenetic incompatibility and linkage equilibrium. Here, we analyzed a total of 17 population samples including the total CGSC population and the four major VG lineages (VGI, VGII, VGIII, and VGIV) at the global level, the total CGSC sample separated by continent (Africa, Asia, Europe, Oceania, North America, and South America), as well as individual lineages separated by continent for populations >20 STs (Asia VGI, Europe VGI, Oceania VGI, North America VGI, North America VGII, and South America VGII). In the phylogenetic incompatibility test, we found that none of the 17 population samples showed 100% phylogenetic compatibility (Table 9). Five of the seventeen population samples (Global CGSC, Global VGI, Global VGII, South America CGSC, and South America VGII) explicitly showed no phylogenetic compatibility among the seven loci. This was consistent with evidence of recombination at the global level among 21 pairs of loci within each of these three samples. VGIII showed 19 of the 21 pairwise loci to be phylogenetically incompatible, while VGIV revealed 16 of 21 pairwise loci to be phylogenetically incompatible. Together, the phylogenetic compatibility test showed that all four VG lineages contained evidence for recombination at the global level. The continental-level analyses displayed much lower frequencies of phylogenetic incompatibility for all continents except South America, which remained consistent for strong evidence for a highly recombing population.
At the global level, linkage disequilibrium analyses revealed that in all five samples, the null hypothesis of random recombination was rejected (Table 9). However, within each of the five samples, variable numbers of pairs of loci showed no significant deviation from those expected under the random recombination null hypothesis. For example, in the global VGI sample, 10 of the 21 loci pairs had observed genotype frequencies not significantly different from random recombination, all involving the CAP59 or IGS1 loci. Similarly, in the global VGIV sample, 5 of the 21 loci pairs had observed genotype frequencies not significantly different from random recombination, all involving the GPD1 locus. In comparison, no evidence for linkage equilibrium across all 21 loci pairs was observed within the total global CGSC sample, the global VGII sample, and the global VGIII sample. The complete allelic associations among loci for all 17 samples are shown in Supplementary Table S12.

4. Discussion

This study extracted MLST data of CGSC from 31 peer-reviewed publications and the IFMLST database. The IFMLST database currently consists of 566 established sequence types (STs) grouped into four VG lineages (VGI, VGII, VGIII, and VGIV) representing the recognized MLST genotypes of CGSC. The extracted MLST data were used to investigate the population structure of the total CGSC and each of the four lineages within CGSC. Below, we discuss the results of our analyses.

4.1. VG Lineage Distribution

As reported in the original literature [34,35,37,38,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75], each of the four CGSC lineages (VGI, VGII, VGIII, and VGIV) in the IFMLST database are broadly distributed across multiple continents. However, the four lineages showed different patterns: the VGI lineage was recorded in 30 of the 45 countries with CGSC, followed by VGII (in 25 countries), VGIII lineage (in 11 countries), and VGIV lineage (in 9 countries). This summarized distribution is consistent with past literature illustrating a greater worldwide distribution of the VGI and VGII lineages compared to the VGIII and VGIV lineages [55]. Among these four lineages at the global level, VGII was the most abundantly represented in the MLST database and was higher than VGI lineage, which was often considered as the most abundant lineage worldwide [55]. However, within Europe as well as in 15 of 45 countries/regions, VGI lineage was the most prominent. Congruent with our global results in the published MLST literature, the VGII lineage was the most prevalent within Africa, Asia, Oceania, and South America as well as in 15 of 45 countries. The VGIII lineage was the most prominent within North America as well as in 6 of 45 countries. This was unexpected, as the VGII lineage is credited with causing the outbreak across Vancouver Island and the Pacific Northeast [23]. Indeed, excluding isolates associated with this outbreak, the VGIII lineage was the most abundant in North America. Such a pattern was used as evidence for an environmental shift of the VGII lineage beyond tropical and sub-tropical regions into temperate North America [81]. Our meta-analyses identified no isolate from the VGIII lineage within Africa or Asia, which is consistent with previous studies that found no isolate of the VGIII lineage in China [56,61]. The VGIV lineage was the most abundant in four countries, including a clonal expansion of ST173 representing 43 isolates in Ivory Coast [81,82,83].

4.2. Geographic and Ecological Structuring

Our global population genetic analyses included 1202 isolates representing 375 STs from 45 countries across six continents. Several STs were broadly distributed, consistent with long-distance clonal expansion of these STs across countries and continents. However, of the 375 STs, 341 STs (90%) were sampled from a single continent, and none were recorded from all six or a combination of five continents. This pattern is different from what was observed in the sister-species complex CNSC, where several STs (ST5 and ST93) were broadly and abundantly distributed [48]. While a lack of CGSC sampling might have contributed to the observed pattern, our population analyses revealed statistically significant differentiations among continental and national populations of CGSC, suggesting relatively limited gene flow across the globe. In the total CGSC sample, the observed genetic differentiations were mostly due to variations in the distribution of CGSC lineages as discussed in the last section as well as localized clonal expansion of certain STs. Indeed, clonal expansion of specific STs was found for all lineages, most of which were localized to individual or subset of countries and continents. For instance, the most abundantly collected ST in VGI was ST106, represented by 37 isolates; in VGII, it was ST7, represented by 167 isolates; in VGIII, it was ST9, represented by 27 isolates; and in VGIV, it was ST336, represented by 5 isolates, and all were predominantly found within one to a few countries. Localized clonal expansion can potentially bias the allele and ST frequencies, contributing to the observed genetic difference among geographic populations [84]. Because of this, we further analyzed clone-corrected samples where only one representative from each geographical region was included. Using clone-corrected samples, the total genetic variance among continents was reduced by 12%, 6%, 12%, 9%, and 5% across the total CGSC population, VGI lineage, VGII lineage, VGIII lineage, and VGIV lineage, respectively. In addition, using clone-corrected samples, the total genetic variance among countries was reduced by 8%, 17%, 23%, 10%, and 12% across the total CGSC population, VGI lineage, VGII lineage, VGIII lineage, and VGIV lineage, respectively. While the contributions by geographic separation were reduced in the clone-corrected samples, in both the total CGSC as well as the VGI and VGII lineages, statistically significant genetic differentiations were still observed and were consistent with the existence of endemic genetic variations within most regional populations. Overall, these results suggest greater evidence for geographical structuring and less evidence for global gene flow in CGSC than in C. neoformans species complex [48].
Ecological data were available for 788 isolates represented by 188 STs. Across the whole CGSC population, VGI lineage, VGII lineage, VGIII lineage, and VGIV lineage, most isolates (416, 52, 284, 70, and 9, respectively, for the five sample types) were sampled from a clinical source, followed by (259, 51, 172, 36, and 0, respectively) from environmental sources, and the remaining (113, 0, 74, 38, and 1, respectively) from veterinary sources. The overall prevalence of isolates recorded from clinical settings is likely due to a greater effort in obtaining and analyzing clinical specimens than from environmental and veterinary sources. Because patient infection by CGSC is from environmental sources, it is important that we increase our efforts to understand environmental populations of this pathogen and identify the potential factors that may contribute to the differences between clinical and environmental samples and how factors in the environment such as fungicide usage may impact clinical antifungal resistance in CGSC [85]. Indeed, our analyses of ecological samples within each of five countries/regions with >100 isolates of CGSC showed statistically significant genetic differentiation in the non-clone-corrected samples but limited differences in the clone-corrected samples. Such results are consistent with the overall similarities among ecological samples in their allelic and sequence types at the seven loci within each of the five countries/regions but that the different ecological niches had clonally expanded frequencies of different allelic and sequence types. Our results suggest that most environmental CGSC alleles and sequence types are capable of infecting humans and animals, but they likely differ in their virulence. Similarly, many human fungal pathogens capable of infecting animals may also serve as a reservoir for the evolution of virulent strains to humans. Understanding the relationships among the ecological populations could have significant implications in our understanding of the evolution of virulence and drug resistance in fungal pathogens [86]. Indeed, in CGSC, genotype sharing among strains from the three different sources were common [Table 3]. For example, CGSC was transferred to a new region within Australia by transporting infected koalas [75]. Similarly, the dispersal of CGSC related to the Vancouver Island outbreak was facilitated by tree cutting and transportation of trees and wood products [87].

4.3. Hybridization and Recombination

We performed phylogenetic analyses using the published 566 STs from the IFMLST database to examine the overall population structure of C. gattii species complex using distance-based models. Overall, the phylogenetic results were consistent with the four main lineages as revealed based on PCR fingerprinting, AFLP, and/or PCR-RFLP of the URA5 gene fragment. There were two cases of inconsistent observed grouping: ST351 and ST402 observed in Figure 1. ST351 was previously designated within the VGI lineage [47] and contains ATs at the CAP59, IGS1, and URA5 loci within VGI. However, the concatenated sequence of ST351 was grouped within the VGII branch and with ATs at the remaining four loci GPD1, LAC1, PLB1, and SOD1 clustered within the VGII lineage. ST402 was previously designated within the VGIV lineage [47] but contains only the AT at the URA5 locus within the VGIV cluster. Interestingly, ST402 contains AT at the PLB1 locus being grouped within the VGII lineage, while alleles at the remaining five loci were clustered within the VGI lineage. The dominance of VGI alleles (5/7) within ST402 led to the clustering of ST402 within the VGI branch. In most cases, lineage designation of strains of human pathogenic Cryptococcus was derived based on PCR-restriction fragment length polymorphism at the URA5 locus [88]. Aside from ST351 and ST402, our analyses identified 21 other STs with ATs from two or more lineages. Together, these results are consistent with hybridization involving all four lineages within CGSC. Based on sequence polymorphisms in the mitochondrial genomes, previous studies showed evidence of hybridization between VGI and VGII lineages [30] and between VGII and VGIII lineages [89]. In addition, both mitochondrial and nuclear genome recombination were suggested within the VGI and VGII lineages [30,33,37]. Though most isolates were classified correctly to individual lineages based on this method, our analyses here suggest that caution should be taken when using the PCR-RFLP results at the URA5 locus for lineage assignment.
In addition to hybridizations among lineages, we identified evidence for recombination within each of the four lineages. The results are consistent with sexual reproduction playing a role in natural populations of CGSC. However, some of the observed phylogenetic incompatibility could be due to reverse mutations, convergent evolution, and/or incomplete lineage sorting. Additional investigations are needed to determine whether some of the nucleotide substitutions at the seven MLST loci have been positively selected and potentially contributed to the observed phylogenetic incompatibilities (and hence recombination). Similar to those observed for CNSC [48], evidence for asexual reproductions were prevalent in CGSC. Asexual reproduction enables rapid colonization of a suitable habitat [90] and spread of an infectious clone while sexual recombination increases the likelihood of generating a diversity of genotypes and potentially facilitating population adaptation to changing environments [91]. Clonality and sexual reproduction are not mutually exclusive, and frequent studies have documented cases of both in many fungal species [52,63,69,92,93]. In CGSC, one report illustrated a prevalence for clonal dispersal in Cucuta, Colombia [59], while another study reported an exceptionally diverse population with evidence for sexual recombination across multiple regions of Colombia [37]. Similarly, a study that looked at STs associated with the Vancouver Island outbreak (ST20 and ST7) concluded there was significant evidence for both clonal expansion and sexual reproduction within the VGII lineage [69].
Several theories have been brought forward to illustrate the potential origin(s) of the Vancouver Island outbreak thought to involve ST20 and ST7 [35,72], with most involving a mixture of clonal expansion and sexual recombination [23,26]. ST7 is a globally distributed sequence type [47]. For instance, 11 out of 12 CGSC isolates identified in Thailand were identical to ST7 associated with the Vancouver outbreak [35], and another report credited clonal expansion of ST7 as responsible for an ongoing VGII outbreak throughout Australia [34]. Our study recorded 167 isolates of ST7, making it the most prevalent ST of CGSC globally, represented across 10 countries/regions (Australia, Brazil, Canada, China, Colombia, South Korea, Taiwan, Thailand, and the USA), demonstrating a global presence of a clonally dispersed ST. In contrast, ST20 is endemic to the Amazon region [72] and was shown to display a large, scattered distribution within the region [39]. Our study recorded 40 cases of ST20, 31 of which were sampled from Brazil. A recent study showed evidence for high rates of recombination within the Brazilian VGII population and suggested that the Vancouver Island outbreak might have emerged from the Brazilian population, followed by clonal dispersal [38,94].
C. gattii species complex has a known sexual phase with two mating types, MATa and MATα, defined at the mating type (MAT) locus [95,96]. In the case of opposite-sex mating, meiotic progeny is expected to have even proportions of the two mating types. However, that is not typically observed in nature [33]. Reports predominantly describe a MATα-dominant population within CGSC and the sister species CNSC. However, some populations, specifically in Australia, have been described as predominantly MATa [97]. Since the MAT locus is not included within the MLST scheme, this study did not account for the distribution of MATa versus MATα strains in our analyses. Nonetheless, it seems uncommon to observe an equal distribution of both mating types in nature. Evidence suggesting sexual recombination within populations with a bias mating type can be explained by same sex mating [23,27,96]. This has been extensively discussed in the Vancouver Island outbreak population [27].
Sexual mating and reproduction among lineages have been observed in the lab. However, most hybrid progeny in laboratory crosses were diploid/aneuploid, with evidence of heterozygosity across most chromosomes [98]. In contrast, the observed hybrids were homozygous at all seven loci. The lack of heterozygosity could be achieved through a series of loss of heterozygosity (LOH) events during mitotic divisions after the initial hybridization [99]. Hybridization has also been recorded between isolates of HPC sister species CGSC and CNSC as well as between distinct lineages (C. neoformans and C. deneoformans) of CNSC [45,100]. Due to the distinct chromosomal structures of different lineages within HPC, meiosis in hybrids among lineages is often incomplete, and diploid or aneuploid sterile offspring are often produced [100]. During subsequent clonal budding, the hybrids could lose heterozygosity, resulting in recombinant haploid yeast cells that are genetically different from both parental strains. Together, these results suggest that although the CGSC population contains historically differentiated lineages, there is ongoing genetic exchange within and among these lineages, potentially generating distinct genotypes that are highly adaptive and virulent to humans.
In conclusion, our meta-analyses of the published MLST data allowed us to determine the genotype distribution and population structure of the C. gattii species complex. Our analyses illustrated geographic structuring of the CGSC population. While evidence for clonal dispersal was found, there is evidence for significant genetic differentiations in all four lineages at both the continental and country levels, with localized clonal expansion contributing to the observed geographic differences within most lineages. However, evidence for 19 incidences of ATs shared across VG lineages and significant evidence for recombination within all lineages demonstrates that genetic exchanges within and among lineages are ongoing. Continued population monitoring of CGSC, especially from environmental sources, is needed in order to effectively cope with the threats caused by this deadly pathogen.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof9020276/s1, Figure S1: Phylogenetic tree at CAP59; Figure S2: Phylogenetic tree at LAC1; Figure S3: Phylogenetic tree at IGS1; Figure S4: Phylogenetic tree at PLB1; Figure S5: Phylogenetic tree at SOD1; Figure S6: Phylogenetic tree at URA5; Table S1: Ecological distribution of CGSC sequence types; Table S2: Ecological distribution of CGSC molecular types; Table S3: Continent distribution of CGSC sequence types; Table S4: Country distribution of CGSC sequence types; Table S5: Pairwise geographic population genetic differentiation of the total CGSC sample; Table S6: Pairwise geographic population genetic differentiation of the total VGI sample; Table S7: Pairwise geographic population genetic differentiation of the total VGII sample; Table S8: Pairwise geographic population genetic differentiation of the total VGIII sample; Table S9: Pairwise geographic population genetic differentiation of the total VGIV sample; Table S10: Pairwise comparison among ecological niches; Table S11: Occurrence of allele types at each of the seven loci among the 566 STs; Table S12: Linkage disequilibrium analyses results.

Author Contributions

Conceptualization, J.X.; methodology, M.H. and J.X.; software, M.H.; validation, M.H. and J.X.; formal analysis, M.H.; investigation, M.H.; resources, M.H. and J.X.; data curation, M.H.; writing—original draft preparation, M.H. and J.X.; writing—review and editing, J.X.; visualization, M.H.; supervision, J.X.; project administration, J.X.; funding acquisition, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Sciences and Engineering Research Council of Canada (grant number RGPIN-2020-05732).

Institutional Review Board Statement

Not applicable. This study analyzed published and publicly available data.

Informed Consent Statement

Not applicable. This study analyzed published and publicly available data.

Data Availability Statement

All data analyzed in this study are cited and are summarized in the manuscript.

Acknowledgments

We thank the ISHAM Cryptococcus MLST subcommittee for their work. We thank Wieland Meyer and his colleagues for maintaining the Cryptococcus MLST database.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Harris, J.; Lockhart, S.; Chiller, T. Cryptococcus gattii: Where do we go from here? Med. Mycol. 2012, 50, 113–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Gushiken, C.A.; Saharia, K.K.; Baddley, W.J. Cryptococcosis. Infect. Dis. Clin. N. Am. 2021, 35, 493–514. [Google Scholar] [CrossRef] [PubMed]
  3. Hagen, F.; Khayhan, K.; Theelen, B.; Kolecka, A.; Polacheck, I.; Sionov, E.; Falk, R.; Parnmen, S.; Lumbsch, H.T.; Boekhout, T. Recognition of seven species in the Cryptococcus gattii/Cryptococcus neoformans species complex. Fungal Genet. Biol. 2015, 78, 16–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lazera, M.S.; Cavalcanti, M.A.S.; Londero, A.T.; Trilles, L.; Nishikawa, M.M.; Wanke, B. Possible primary ecological niche of Cryptococcus neoformans. Med. Mycol. 2000, 38, 379–383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Xu, J.; Manosuthi, W.; Banerjee, U.; Zhu, L.-P.; Chen, J.H.; Kohno, S.; Izumikawa, K.; Chen, Y.H.; Sungkanuparph, S.; Harrison, T.S.; et al. Cryptococcosis in Asia. In Cryptococcus: From Human Pathogen to Model Organism; Heitman, J., Kwon-Chung, J., Perfect, J., Casadevall, A., Eds.; ASM Press: Washington, DC, USA, 2011; pp. 287–298. [Google Scholar]
  6. Randhawa, H.S.; Kowshik, T.; Chowdhary, A.; Preeti Sinha, K.; Khan, Z.U.; Sun, S.; Xu, J. The expanding host tree species spectrum of Cryptococcus gattii and Cryptococcus neoformans and their isolations from surrounding soil in India. Med. Mycol. 2008, 46, 823–833. [Google Scholar] [CrossRef] [Green Version]
  7. Maziarz, E.K.; Perfect, J.R. Cryptococcosis. Infect. Dis. Clin. N. Am. 2016, 30, 179–206. [Google Scholar] [CrossRef] [Green Version]
  8. Kwon-Chung, K.J.; Boekhout, T.; Fell, J.W.; Diaz, M. (1557) Proposal to conserve the name Cryptococcus gattii against C. hondurianus and C. bacillisporus (Basidiomycota, Hymenomycetes, Tremellomycetidae). Taxon 2002, 51, 804–806. [Google Scholar] [CrossRef]
  9. Kwon-Chung, K.J.; Bennett, J.E. High Prevalence of Cryptococcus neoformans var. gattii in Tropical and Subtropical Regions. Zent. Bakteriol. Mikrobiol. Hygiene. 1. Abt. Originale. A Med. Mikrobiol. Infekt. Parasitol. 1984, 257, 213–218. [Google Scholar] [CrossRef]
  10. Kwon-Chung, K.J.; Bennett, J.E. Epidemiologic differences between the two varieties of Cryptococcus neoformans. Am. J. Epidemiol. 1984, 120, 123–130. [Google Scholar] [CrossRef]
  11. Seaton, R.A.; Verma, N.; Naraqi, S.; Wembri, J.P.; Warrell, D.A. Visual loss in immunocompetent patients with Cryptococcus neoformans var. gattii meningitis. Trans. R. Soc. Trop. Med. Hyg. 1997, 91, 44–49. [Google Scholar] [CrossRef]
  12. Speed, B.; Dunt, D. Clinical and Host Differences between Infections with the Two Varieties of Cryptococcus neoformans. Clin. Infect. Dis. 1995, 21, 28–34. [Google Scholar] [CrossRef]
  13. Ellis, D.H. Cryptococcus neoformans var. gattii in Australia. J. Clin. Microbiol. 1987, 25, 430–431. [Google Scholar] [CrossRef] [Green Version]
  14. Mitchell, D.H.; Sorrell, T.C.; Allworth, A.M.; Heath, C.H.; McGregor, A.R.; Papanaoum, K.; Richards, M.J.; Gottlieb, T. Cryptococcal Disease of the CNS in Immunocompetent Hosts: Influence of Cryptococcal Variety on Clinical Manifestations and Outcome. Clin. Infect. Dis. 1995, 20, 611–616. [Google Scholar] [CrossRef] [PubMed]
  15. Lalloo, D.; Fisher, D.; Naraqi, S.; Laurenson, I.; Temu, P.; Sinha, A.; Saweri, A.; Mavo, B. Cryptococcal meningitis (C. neoformans var. gattii) leading to blindness in previously healthy melanesian adults in Papua New Guinea. QJM Int. J. Med. 1994, 87, 343–349. [Google Scholar] [CrossRef]
  16. Fyfe, M.; MacDougall, L.; Romney, M.; Starr, M.; Pearce, M.; Mak, S. Cryptococcus gattii infections on Vancouver Island, British Columbia, Canada: Emergence of a tropical fungus in a temperate climate. Can Commun. Dis. Rep. 2008, 34, 1–13. [Google Scholar] [PubMed]
  17. Stephen, C.; Lester, S.; Black, W.; Fyfe, M.; Raverty, S. Multispecies outbreak of cryptococcosis on southern Vancouver Island, British Columbia. Can. Vet. J. 2002, 43, 792–794. [Google Scholar]
  18. Galanis, E.; Macdougall, L.; Kidd, S.; Morshed, M. Epidemiology of Cryptococcus gattii, British Columbia, Canada, 1999–2007. Emerg. Infect. Dis. 2010, 16, 251–257. [Google Scholar] [CrossRef]
  19. Datta, K.; Bartlett, K.H.; Baer, R.; Byrnes, E.; Galanis, E.; Heitman, J.; Hoang, L.; Leslie, M.J.; MacDougall, L.; Magill, S.S.; et al. Spread of Cryptococcus gattii into Pacific Northwest region of the United States. Emerg. Infect. Dis. 2009, 15, 1185–1191. [Google Scholar] [CrossRef]
  20. Bartlett, K.H.; Kidd, S.E.; Kronstad, J.W. The emergence of Cryptococcus gattii in British Columbia and the Pacific Northwest. Curr. Infect. Dis. Rep. 2008, 10, 58–65. [Google Scholar] [CrossRef]
  21. Gugnani, H.C.; Mitchell, T.G.; Litvintseva, A.P.; Lengeler, K.B.; Heitman, J.; Kumar, A.; Basu, S.; Paliwal-Joshi, A. Isolation of Cryptococcus gattii and Cryptococcus neoformans var. grubii from the flowers and bark of Eucalyptus trees in India. Med. Mycol. 2005, 43, 565–569. [Google Scholar] [CrossRef] [Green Version]
  22. Campisi, E.; Mancianti, F.; Pini, G.; Faggi, E.; Gargani, G. Investigation in central Italy of the possible association between Cryptococcus neoformans var. Gattii and Eucalyptus camaldulensis. Eur. J. Epidemiol. 2003, 18, 357–362. [Google Scholar] [CrossRef]
  23. Kidd, S.E.; Hagen, F.; Tscharke, R.L.; Huynh, M.; Bartlett, K.H.; Fyfe, M.; MacDougall, L.; Boekhout, T.; Kwon-Chung, K.J.; Meyer, W. A rare genotype of Cryptococcus gattii caused the cryptococcosis outbreak on Vancouver Island (British Columbia, Canada). Proc. Natl. Acad. Sci. USA 2004, 101, 17258–17263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sukroongreung, S.; Kitiniyom, K.; Nilakul, C.; Tantimavanich, S. Pathogenicity of basidiospores of Filobasidiella neoformans var. neoformans. Med. Mycol. 1998, 36, 419–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hull, C.M.; Heitman, J. Genetics of Cryptococcus neoformans. Annu. Rev. Genet. 2002, 36, 557–615. [Google Scholar] [CrossRef] [PubMed]
  26. Fraser, J.A.; Subaran, R.L.; Nichols, C.B.; Heitman, J. Recapitulation of the Sexual Cycle of the Primary Fungal Pathogen Cryptococcus neoformans var. gattii: Implications for an Outbreak on Vancouver Island, Canada. Eukaryot. Cell 2003, 2, 1036–1045. [Google Scholar] [CrossRef] [Green Version]
  27. Fraser, J.A.; Giles, S.S.; Wenink, E.C.; Geunes-Boyer, S.G.; Wright, J.R.; Diezmann, S.; Allen, A.; Stajich, J.E.; Dietrich, F.S.; Perfect, J.R.; et al. Same-sex mating and the origin of the Vancouver Island Cryptococcus gattii outbreak. Nature 2005, 437, 1360–1364. [Google Scholar] [CrossRef] [PubMed]
  28. Lin, X.; Hull, C.M.; Heitman, J. Sexual reproduction between partners of the same mating type in Cryptococcus neoformans. Nature 2005, 434, 1017–1021. [Google Scholar] [CrossRef]
  29. Wickes, B.L.; Mayorga, M.E.; Edman, U.; Edman, J.C. Dimorphism and haploid fruiting in Cryptococcus neoformans: Association with the alpha-mating type. Proc. Natl. Acad. Sci. 1996, 93, 7327–7331. [Google Scholar] [CrossRef] [Green Version]
  30. Xu, J.; Yan, Z.; Guo, H. Divergence, hybridization, and recombination in the mitochondrial genome of the human pathogenic yeast Cryptococcus gattii. Mol Ecol 2009, 18, 2628–2642. [Google Scholar] [CrossRef]
  31. Franzot, S.P.; Hamdan, J.S.; Currie, B.P.; Casadevall, A. Molecular epidemiology of Cryptococcus neoformans in Brazil and the United States: Evidence for both local genetic differences and a global clonal population structure. J. Clin. Microbiol. 1997, 35, 2243–2251. [Google Scholar] [CrossRef] [Green Version]
  32. Kwon-Chung, K.J.; Bennett, J.E. Distribution of α and α mating types of Cryptococcus neoformans among natural and clinical isolates. Am. J. Epidemiol. 1978, 108, 337–340. [Google Scholar] [CrossRef]
  33. Chowdhary, A.; Hiremath, S.S.; Sun, S.; Kowshik, T.; Randhawa, H.S.; Xu, J. Genetic differentiation, recombination and clonal expansion in environmental populations of Cryptococcus gattii in India. Environ. Microbiol. 2011, 13, 1875–1888. [Google Scholar] [CrossRef] [PubMed]
  34. Carriconde, F.; Gilgado, F.; Arthur, I.; Ellis, D.; Malik, R.; van de Wiele, N.; Robert, V.; Currie, B.J.; Meyer, W. Clonality and α-a Recombination in the Australian Cryptococcus gattii VGII Population—An Emerging Outbreak in Australia. PLoS ONE 2011, 6, e16936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kaocharoen, S.; Ngamskulrungroj, P.; Firacative, C.; Trilles, L.; Piyabongkarn, D.; Banlunara, W.; Poonwan, N.; Chaiprasert, A.; Meyer, W.; Chindamporn, A. Molecular epidemiology reveals genetic diversity amongst isolates of the Cryptococcus neoformans/C. gattii species complex in Thailand. PLoS Negl. Trop. Dis. 2013, 7, e2297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Campbell, L.T.; Currie, B.J.; Krockenberger, M.; Malik, R.; Meyer, W.; Heitman, J.; Carter, D. Clonality and Recombination in Genetically Differentiated Subgroups of Cryptococcus gattii. Eukaryot Cell 2005, 4, 1403–1409. [Google Scholar] [CrossRef] [Green Version]
  37. Vélez, N.; Escandón, P. Multilocus sequence typing (MLST) of clinical and environmental isolates of Cryptococcus neoformans and Cryptococcus gattii in six departments of Colombia reveals high genetic diversity. Rev. Soc. Bras. Med. Trop. 2020, 53, e20190422. [Google Scholar] [CrossRef] [PubMed]
  38. Souto, A.C.; Bonfietti, L.X.; Ferreira-Paim, K.; Trilles, L.; Martins, M.; Ribeiro-Alves, M.; Pham, C.D.; Martins, L.; Dos Santos, W.; Chang, M.; et al. Population Genetic Analysis Reveals a High Genetic Diversity in the Brazilian Cryptococcus gattii VGII Population and Shifts the Global Origin from the Amazon Rainforest to the Semi-arid Desert in the Northeast of Brazil. PLoS Negl. Trop. Dis. 2016, 10, e0004885. [Google Scholar] [CrossRef] [Green Version]
  39. Meyer, W.; Castañeda, A.; Jackson, S.; Huynh, M.; Castañeda, E. Molecular Typing of IberoAmerican Cryptococcus neoformans Isolates. Emerg. Infect. Dis. 2003, 9, 189–195. [Google Scholar] [CrossRef]
  40. Boekhout, T.; Theelen, B.; Diaz, M.; Fell, J.W.; Hop, W.C.J.; Abeln, E.C.A.; Dromer, F.; Meyer, W. Hybrid genotypes in the pathogenic yeast Cryptococcus neoformans. Microbiology 2001, 147, 891–907. [Google Scholar] [CrossRef] [Green Version]
  41. Hong, N.; Chen, M.; Xu, J. Molecular Markers Reveal Epidemiological Patterns and Evolutionary Histories of the Human Pathogenic Cryptococcus. Front. Cell. Infect. Microbiol. 2021, 11, 683670. [Google Scholar] [CrossRef]
  42. Yamamura, D.; Xu, J. Update on Pulmonary Cryptococcosis. Mycopathologia 2021, 186, 717–728. [Google Scholar] [CrossRef] [PubMed]
  43. Farrer, R.A.; Chang, M.; Davis, M.J.; Van Dorp, L.; Yang, D.-H.; Shea, T.; Sewell, T.R.; Meyer, W.; Balloux, F.; Edwards, H.M.; et al. A New Lineage of Cryptococcus gattii (VGV) Discovered in the Central Zambezian Miombo Woodlands. mBio 2019, 10, e02306–e02319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Monroy-Nieto, J.; Bowers, J.R.; Montfort, P.; Adame, G.; Taverna, C.G.; Yaglom, H.; Sykes, J.E.; Brady, S.; Mochon, A.B.; Meyer, W.; et al. Phylogenomic Placement of American Southwest-Associated Clinical and Veterinary Isolates Expands Evidence for Distinct Cryptococcus gattii VGVI. Microorganisms 2022, 10, 1681. [Google Scholar] [CrossRef] [PubMed]
  45. Samarasinghe, H.; You, M.; Jenkinson, T.S.; Xu, J.; James, T.Y. Hybridization Facilitates Adaptive Evolution in Two Major Fungal Pathogens. Genes 2020, 11, 101. [Google Scholar] [CrossRef] [Green Version]
  46. Meyer, W.; Aanensen, D.M.; Boekhout, T.; Cogliati, M.; Diaz, M.R.; Esposto, M.C.; Fisher, M.; Gilgado, F.; Hagen, F.; Kaocharoen, S.; et al. Consensus multi-locus sequence typing scheme for Cryptococcus neoformans and Cryptococcus gattii. Med. Mycol. 2009, 47, 561–570. [Google Scholar] [CrossRef] [Green Version]
  47. Meyer, W. International Fungal Multi Locus Sequence Typing Database. Available online: https://mlst.mycologylab.org/page/Home1 (accessed on 15 January 2023).
  48. Hitchcock, M.; Xu, J. Analyses of the Global Multilocus Genotypes of the Human Pathogenic Yeast Cryptococcus neoformans Species Complex. Genes 2022, 13, 2045. [Google Scholar] [CrossRef]
  49. Xess, I.; Pandey, M.; Dabas, Y.; Agarwal, R.; Das, S.; Srivastava, P.M.V.; Thakur, R.; Sharma, S.; Mani, P.; Biswas, A.; et al. Multilocus Sequence Typing of Clinical Isolates of Cryptococcus from India. Mycopathologia 2021, 186, 199–211. [Google Scholar] [CrossRef]
  50. Rocha, D.F.S.; Cruz, K.S.; Santos, C.; Menescal, L.S.F.; Neto, J.; Pinheiro, S.B.; Silva, L.M.; Trilles, L.; Braga de Souza, J.V. MLST reveals a clonal population structure for Cryptococcus neoformans molecular type VNI isolates from clinical sources in Amazonas, Northern-Brazil. PLoS ONE 2018, 13, e0197841. [Google Scholar] [CrossRef] [Green Version]
  51. Moslem, M.; Fatahinia, M.; Kiasat, N.; Mahmoudabadi, A.Z. Genotypic diversity of Iranian Cryptococcus neoformans using multilocus sequence typing (MLST) and susceptibility to antifungals. Mol. Biol. Rep. 2021, 48, 4201–4208. [Google Scholar] [CrossRef]
  52. Huang, C.; Tsui, C.K.M.; Chen, M.; Pan, K.; Li, X.; Wang, L.; Chen, M.; Zheng, Y.; Zheng, D.; Chen, X.; et al. Emerging Cryptococcus gattii species complex infections in Guangxi, southern China. PLoS Negl. Trop. Dis. 2020, 14, e0008493. [Google Scholar] [CrossRef]
  53. Firacative, C.; Roe, C.C.; Malik, R.; Ferreira-Paim, K.; Escandón, P.; Sykes, J.E.; Castañón-Olivares, L.R.; Contreras-Peres, C.; Samayoa, B.; Sorrell, T.C.; et al. MLST and Whole-Genome-Based Population Analysis of Cryptococcus gattii VGIII Links Clinical, Veterinary and Environmental Strains, and Reveals Divergent Serotype Specific Sub-populations and Distant Ancestors. PLoS Negl. Trop. Dis. 2016, 10, e0004861. [Google Scholar] [CrossRef] [Green Version]
  54. Bellet, V.; Roger, F.; Krasteva, D.; Gouveia, T.; Drakulovski, P.; Pottier, C.; Bertout, S. Multilocus sequence typing of strains from the Cryptococcus gattii species complex from different continents. Mycoses 2022, 65, 88–96. [Google Scholar] [CrossRef]
  55. Meyer, W.; Trilles, L. Genotyping of the Cryptococcus neoformans/Cryptococcus gattii species complex. Aust. Biochem. 2010, 41, 12–16. [Google Scholar]
  56. Wu, S.-Y.; Kang, M.; Liu, Y.; Chen, Z.-X.; Xiao, Y.-L.; He, C.; Ma, Y. Molecular epidemiology and antifungal susceptibilities of Cryptococcus species isolates from HIV and non-HIV patients in Southwest China. Eur. J. Clin. Microbiol. Infect. Dis. 2021, 40, 287–295. [Google Scholar] [CrossRef] [PubMed]
  57. Jin, L.; Cao, J.R.; Xue, X.Y.; Wu, H.; Wang, L.F.; Guo, L.; Shen, D.X. Clinical and microbiological characteristics of Cryptococcus gattii isolated from 7 hospitals in China. BMC Microbiol. 2020, 20, 73. [Google Scholar] [CrossRef] [PubMed]
  58. Maruyama, F.H.; de Paula, D.A.J.; Menezes, I.D.G.; Favalessa, O.C.; Hahn, R.C.; de Almeida, A.D.B.P.F.; Sousa, V.R.F.; Nakazato, L.; Dutra, V. Genetic Diversity of the Cryptococcus gattii Species Complex in Mato Grosso State, Brazil. Mycopathologia 2019, 184, 45–51. [Google Scholar] [CrossRef]
  59. Firacative, C.; Torres, G.; Meyer, W.; Escandón, P. Clonal Dispersal of Cryptococcus gattii VGII in an Endemic Region of Cryptococcosis in Colombia. J. Fungi 2019, 5, 32. [Google Scholar] [CrossRef] [Green Version]
  60. Dou, H.T.; Xu, Y.C.; Wang, H.Z.; Li, T.S. Molecular epidemiology of Cryptococcus neoformans and Cryptococcus gattii in China between 2007 and 2013 using multilocus sequence typing and the DiversiLab system. Eur. J. Clin. Microbiol. Infect. Dis. 2015, 34, 753–762. [Google Scholar] [CrossRef]
  61. Wu, S.-Y.; Lei, Y.; Kang, M.; Xiao, Y.-L.; Chen, Z.-X. Molecular characterisation of clinical Cryptococcus neoformans and Cryptococcus gattii isolates from Sichuan province, China. Mycoses 2015, 58, 280–287. [Google Scholar] [CrossRef] [PubMed]
  62. Bellet, V.; Kassi, F.K.; Krasteva, D.; Roger, F.; Drakulovski, P.; Mossou, C.; Kouakou, G.A.; Doumbia, A.; Delaporte, E.; Menan, H.; et al. First report of cryptococcosis due to Cryptococcus gattii sensu stricto VGI in an Ivorian HIV negative patient. J. Med. Mycol. 2021, 31, 101113. [Google Scholar] [CrossRef]
  63. Lin, K.-H.; Lin, Y.-P.; Ho, M.-W.; Chen, Y.-C.; Chung, W.-H. Molecular epidemiology and phylogenetic analyses of environmental and clinical isolates of Cryptococcus gattii sensu lato in Taiwan. Mycoses 2021, 64, 324–335. [Google Scholar] [CrossRef] [PubMed]
  64. Fan, X.; Xiao, M.; Chen, S.; Kong, F.; Dou, H.-T.; Wang, H.; Xiao, Y.-L.; Kang, M.; Sun, Z.-Y.; Hu, Z.-D.; et al. Predominance of Cryptococcus neoformans var. grubii multilocus sequence type 5 and emergence of isolates with non-wild-type minimum inhibitory concentrations to fluconazole: A multi-centre study in China. Clin. Microbiol. Infect. 2016, 22, 887.e1–887.e9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Barcellos, V.A.; Martins, L.M.S.; Fontes, A.C.L.; Reuwsaat, J.C.V.; Squizani, E.D.; de Sousa Araújo, G.R.; Frases, S.; Staats, C.C.; Schrank, A.; Kmetzsch, L.; et al. Genotypic and Phenotypic Diversity of Cryptococcus gattii VGII Clinical Isolates and Its Impact on Virulence. Front. Microbiol. 2018, 9, 132. [Google Scholar] [CrossRef] [Green Version]
  66. Cogliati, M.; Desnos-Ollivier, M.; McCormick-Smith, I.; Rickerts, V.; Ferreira-Paim, K.; Meyer, W.; Boekhout, T.; Hagen, F.; Theelen, B.; Inácio, J.; et al. Genotypes and population genetics of Cryptococcus neoformans and Cryptococcus gattii species complexes in Europe and the Mediterranean area. Fungal Genet. Biol. 2019, 129, 16–29. [Google Scholar] [CrossRef]
  67. Pinheiro, S.B.; Sousa, E.S.; Cortez, A.C.A.; da Silva Rocha, D.F.; Menescal, L.S.F.; Chagas, V.S.; Gómez, A.S.P.; Cruz, K.S.; Santos, L.O.; Alves, M.J.; et al. Cryptococcal meningitis in non-HIV patients in the State of Amazonas, Northern Brazil. Braz. J. Microbiol. 2021, 52, 279–288. [Google Scholar] [CrossRef] [PubMed]
  68. Brito-Santos, F.; Reis, R.S.; Coelho, R.A.; Almeida-Paes, R.; Pereira, S.A.; Trilles, L.; Meyer, W.; Wanke, B.; Lazéra, M.D.S.; Gremião, I.D.F. Cryptococcosis due to Cryptococcus gattii VGII in southeast Brazil: The One Health approach revealing a possible role for domestic cats. Med. Mycol. Case Rep. 2019, 24, 61–64. [Google Scholar] [CrossRef]
  69. Brito-Santos, F.; Barbosa, G.G.; Trilles, L.; Nishikawa, M.M.; Wanke, B.; Meyer, W.; Carvalho-Costa, F.A.; Lazéra Mdos, S. Environmental isolation of Cryptococcus gattii VGII from indoor dust from typical wooden houses in the deep Amazonas of the Rio Negro basin. PLoS ONE 2015, 10, e0115866. [Google Scholar] [CrossRef] [Green Version]
  70. Lizarazo, J.; Escandón, P.; Agudelo, C.I.; Firacative, C.; Meyer, W.; Castañeda, E. Retrospective study of the epidemiology and clinical manifestations of Cryptococcus gattii infections in Colombia from 1997–2011. PLoS Negl. Trop. Dis. 2014, 8, e3272. [Google Scholar] [CrossRef] [Green Version]
  71. Fang, L.-f.; Zhang, P.-p.; Wang, J.; Yang, Q.; Qu, T.-t. Clinical and microbiological characteristics of cryptococcosis at an university hospital in China from 2013 to 2017. Braz. J. Infect. Dis. 2020, 24, 7–12. [Google Scholar] [CrossRef]
  72. Brito-Santos, F.; Trilles, L.; Firacative, C.; Wanke, B.; Carvalho-Costa, F.A.; Nishikawa, M.M.; Campos, J.P.; Junqueira, A.C.V.; Souza, A.C.; Lazéra, M.D.S.; et al. Indoor Dust as a Source of Virulent Strains of the Agents of Cryptococcosis in the Rio Negro Micro-Region of the Brazilian Amazon. Microorganisms 2020, 8, 682. [Google Scholar] [CrossRef]
  73. Fu, Y.; Xu, M.; Zhou, H.; Yao, Y.; Zhou, J.; Pan, Z. Microbiological and clinical characteristics of cryptococcemia: A retrospective analysis of 85 cases in a Chinese hospital. Med. Mycol. 2019, 58, 478–484. [Google Scholar] [CrossRef]
  74. Kassi, F.K.; Drakulovski, P.; Bellet, V.; Roger, F.; Chabrol, A.; Krasteva, D.; Doumbia, A.; Landman, R.; Kakou, A.; Reynes, J.; et al. Cryptococcus genetic diversity and mixed infections in Ivorian HIV patients: A follow up study. PLOS Negl. Trop. Dis. 2019, 13, e0007812. [Google Scholar] [CrossRef] [Green Version]
  75. Schmertmann, L.J.; Danesi, P.; Monroy-Nieto, J.; Bowers, J.; Engelthaler, D.M.; Malik, R.; Meyer, W.; Krockenberger, M.B. Jet-Setting Koalas Spread Cryptococcus gattii VGII in Australia. mSphere 2019, 4, e00216-19. [Google Scholar] [CrossRef] [Green Version]
  76. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  77. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  78. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  79. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Agapow, P.-M.; Burt, A. Indices of multilocus linkage disequilibrium. Mol. Ecol. Notes 2001, 1, 101–102. [Google Scholar] [CrossRef]
  81. Lockhart, S.R.; Iqbal, N.; Harris, J.R.; Grossman, N.T.; DeBess, E.; Wohrle, R.; Marsden-Haug, N.; Vugia, D.J. Cryptococcus gattii in the United States: Genotypic Diversity of Human and Veterinary Isolates. PLoS ONE 2013, 8, e74737. [Google Scholar] [CrossRef]
  82. Cogliati, M. Global Molecular Epidemiology of Cryptococcus neoformans and Cryptococcus gattii: An Atlas of the Molecular Types. Science 2013, 2013, 675213. [Google Scholar] [CrossRef] [Green Version]
  83. Litvintseva, A.P.; Thakur, R.; Reller, L.B.; Mitchell, T.G. Prevalence of Clinical Isolates of Cryptococcus gattii Serotype C among Patients with AIDS in Sub-Saharan Africa. J. Infect. Dis. 2005, 192, 888–892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Xu, J. Fundamentals of fungal molecular population genetic analyses. Curr. Issues Mol. Biol. 2006, 8, 75–89. [Google Scholar] [PubMed]
  85. Fisher, M.C.; Alastruey-Izquierdo, A.; Berman, J.; Bicanic, T.; Bignell, E.M.; Bowyer, P.; Bromley, M.; Brüggemann, R.; Garber, G.; Cornely, O.A.; et al. Tackling the emerging threat of antifungal resistance to human health. Nat. Rev. Microbiol. 2022, 20, 557–571. [Google Scholar] [CrossRef]
  86. Xu, J. Assessing global fungal threats to humans. mLife 2022, 1, 223–240. [Google Scholar] [CrossRef]
  87. Kidd, S.E.; Bach, P.J.; Hingston, A.O.; Mak, S.; Chow, Y.; MacDougall, L.; Kronstad, J.W.; Bartlett, K.H. Cryptococcus gattii Dispersal Mechanisms, British Columbia, Canada. Emerg. Infect. Dis. 2007, 13, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Velegraki, A.; Kiosses, V.G.; Kansouzidou, A.; Smilakou, S.; Mitroussia-Ziouva, A.; Legakis, N.J. Prospective use of RFLP analysis on amplified Cryptococcus neoformans URA5 gene sequences for rapid identification of varieties and serotypes in clinical samples. Med. Mycol. 2001, 39, 409–417. [Google Scholar] [CrossRef] [Green Version]
  89. Farrer, R.A.; Desjardins, C.A.; Sakthikumar, S.; Gujja, S.; Saif, S.; Zeng, Q.; Chen, Y.; Voelz, K.; Heitman, J.; May, R.C.; et al. Genome Evolution and Innovation across the Four Major Lineages of Cryptococcus gattii. mBio 2015, 6, e00868-15. [Google Scholar] [CrossRef] [Green Version]
  90. Fisher, M.C.; Hanage, W.P.; de Hoog, S.; Johnson, E.; Smith, M.D.; White, N.J.; Vanittanakom, N. Low Effective Dispersal of Asexual Genotypes in Heterogeneous Landscapes by the Endemic Pathogen Penicillium marneffei. PLOS Pathog. 2005, 1, e20. [Google Scholar] [CrossRef] [Green Version]
  91. Goddard, M.R.; Godfray, H.C.; Burt, A. Sex increases the efficacy of natural selection in experimental yeast populations. Nature 2005, 434, 636–640. [Google Scholar] [CrossRef]
  92. Billmyre, R.B.; Croll, D.; Li, W.; Mieczkowski, P.; Carter, D.A.; Cuomo, C.A.; Kronstad, J.W.; Heitman, J. Highly Recombinant VGII Cryptococcus gattii Population Develops Clonal Outbreak Clusters through both Sexual Macroevolution and Asexual Microevolution. mBio 2014, 5, e01494-14. [Google Scholar] [CrossRef] [Green Version]
  93. Hiremath, S.S.; Chowdhary, A.; Kowshik, T.; Randhawa, H.S.; Sun, S.; Xu, J. Long-distance dispersal and recombination in environmental populations of Cryptococcus neoformans var. grubii from India. Microbiology 2008, 154 Pt 5, 1513–1524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Hagen, F.; Ceresini, P.C.; Polacheck, I.; Ma, H.; van Nieuwerburgh, F.; Gabaldón, T.; Kagan, S.; Pursall, E.R.; Hoogveld, H.L.; van Iersel, L.J.J.; et al. Ancient Dispersal of the Human Fungal Pathogen Cryptococcus gattii from the Amazon Rainforest. PLoS ONE 2013, 8, e71148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Kwon-Chung, K.J. A New Genus, Filobasidiella, The Perfect State of Cryptococcus neoformans. Mycologia 1975, 67, 1197–1200. [Google Scholar] [CrossRef]
  96. Kwon-Chung, K.J. A New Species of Filobasidiella, The Sexual State of Cryptococcus neoformans B and C Serotypes. Mycologia 1976, 68, 942–946. [Google Scholar] [CrossRef]
  97. Halliday, C.L.; Bui, T.; Krockenberger, M.; Malik, R.; Ellis, D.H.; Carter, D.A. Presence of α and Mating Types in Environmental and Clinical Collections of Cryptococcus neoformans var. gattii Strains from Australia. J. Clin. Microbiol. 1999, 37, 2920–2926. [Google Scholar] [CrossRef] [Green Version]
  98. You, M.; Xu, J. What Are the Best Parents for Hybrid Progeny? An Investigation into the Human Pathogenic Fungus Cryptococcus. J. Fungi 2021, 7, 299. [Google Scholar] [CrossRef] [PubMed]
  99. You, M.; Lin, Y.M.; Dobrin, A.; Xu, J. Genetic and Phenotypic Diversities in Experimental Populations of Diploid Inter-Lineage Hybrids in the Human Pathogenic Cryptococcus. Microorganisms 2021, 9, 1579. [Google Scholar] [CrossRef]
  100. Samarasinghe, H.; Xu, J. Hybrids and hybridization in the Cryptococcus neoformans and Cryptococcus gattii species complexes. Infect. Genet. Evol. 2018, 66, 245–255. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic tree illustrating the relationship among 566 sequence types (ST) of the C. gattii species complex. The tree was constructed using the neighbor-joining method based on K2P distance in MEGA11. The branch lengths are proportional to sequence divergence among sequence types. Three large circles and ten small circles are used to display the association of metadata with each ST. The outermost large circle, containing three small circles, demonstrates the ecological niche each ST was isolated from, with each color representing a specific niche (clinical, environmental, and veterinary). The middle large circle contains only one circle, displaying four colors, each representing one of the four CGSC lineages (VGI, VGII, VGIII, and VGIV). The inner most circle, containing six smaller circles, displays the continental distribution of each ST, with each color representing one continent (Africa, Asia, Europe, Oceania, North America, and South America). Legends on the left indicate the correspondence between each color and the associate metadata for each ST. STs shaded in blue contain alleles shared between more than one lineage.
Figure 1. Phylogenetic tree illustrating the relationship among 566 sequence types (ST) of the C. gattii species complex. The tree was constructed using the neighbor-joining method based on K2P distance in MEGA11. The branch lengths are proportional to sequence divergence among sequence types. Three large circles and ten small circles are used to display the association of metadata with each ST. The outermost large circle, containing three small circles, demonstrates the ecological niche each ST was isolated from, with each color representing a specific niche (clinical, environmental, and veterinary). The middle large circle contains only one circle, displaying four colors, each representing one of the four CGSC lineages (VGI, VGII, VGIII, and VGIV). The inner most circle, containing six smaller circles, displays the continental distribution of each ST, with each color representing one continent (Africa, Asia, Europe, Oceania, North America, and South America). Legends on the left indicate the correspondence between each color and the associate metadata for each ST. STs shaded in blue contain alleles shared between more than one lineage.
Jof 09 00276 g001
Figure 2. Phylogenetic tree illustrating the relationship among 61 allele types (AT) found at the GPD1 locus of the C. gattii species complex. The tree was constructed using the neighbor-joining method based on K2P distance in MEGA11. The branch lengths are proportional to sequence divergence among allele types. Four colors are used to display in which lineage each AT is recorded. Legend on the left indicates which color corresponds to which lineage (VGI, VGII, VGIII, and VGIV). Six of the sixty-one ATs were associated with STs belonging to at least two different lineages.
Figure 2. Phylogenetic tree illustrating the relationship among 61 allele types (AT) found at the GPD1 locus of the C. gattii species complex. The tree was constructed using the neighbor-joining method based on K2P distance in MEGA11. The branch lengths are proportional to sequence divergence among allele types. Four colors are used to display in which lineage each AT is recorded. Legend on the left indicates which color corresponds to which lineage (VGI, VGII, VGIII, and VGIV). Six of the sixty-one ATs were associated with STs belonging to at least two different lineages.
Jof 09 00276 g002
Table 1. Summary lineage distribution of strains of the C. gattii species complex genotyped using the ISHAM multi-locus sequence typing scheme.
Table 1. Summary lineage distribution of strains of the C. gattii species complex genotyped using the ISHAM multi-locus sequence typing scheme.
Sample GroupNo. of Recorded STNo. of ST with Geographical DataNo. of Isolates with Geographical DataFound in No. of ContinentsFound in No. of CountriesNo. of ST with Ecological DataNo. of Isolates with Ecological DataFound in No. of Ecological Sources
Total CGSC56637512026451887883
VGI16097225631231032
VGII2181817216251055303
VGIII7572221412561443
VGIV1112434597102
ST previously unassigned to a lineage21111111
No., number; Ecological sources are defined as: clinical, environmental, and veterinary.
Table 2. Summary geographic distributions of isolates of the C. gattii species complex genotyped using the ISHAM multi-locus sequence typing scheme.
Table 2. Summary geographic distributions of isolates of the C. gattii species complex genotyped using the ISHAM multi-locus sequence typing scheme.
Regionn%VGIVGIIVGIIIVGIVRegionn%VGIVGIIVGIIIVGIV
Africa867.15%2045021North America17914.89%25151390
Ivory Coast453.74%24300USA14111.73%2071140
DRC121.00%10200Mexico252.08%10240
South Africa121.00%4008Canada100.83%3700
Zimbabwe121.00%00012Aruba10.08%0100
Rwanda30.25%2001Cuba10.08%1000
Kenya20.17%2000Guatemala10.08%0010
Asia19616.31%8714109Oceania15913.23%2412942
Taiwan1119.23%446601Australia15312.73%2012922
China433.58%281400New Zealand30.25%1020
Thailand171.41%31400Papa New Guinea30.25%3000
India100.83%6004
Malaysia100.83%3304
Singapore30.25%3000
Japan10.08%0100
South Korea10.08%0100
Europe655.41%441731South America51743.01%25416751
Spain141.16%12110Brazil34929.03%934000
Greece110.92%8300Colombia14512.06%1363681
France90.75%5400Argentina100.83%1720
Italy90.75%9000Venezuela50.42%0320
The Netherlands80.67%2600Paraguay20.17%0020
Germany50.42%2120Peru20.17%2000
Portugal30.25%3000Uruguay20.17%0200
Belgium20.17%2000Chile10.08%0010
Switzerland20.17%1100French Guiana10.08%0100
Denmark10.08%0100
Sweden10.08%0001
n, total isolates; %, percent of isolates of the total recorded in CGSC.
Table 3. Summary distribution of 375 sequence types of the C. gattii species complex among continents and ecological niches.
Table 3. Summary distribution of 375 sequence types of the C. gattii species complex among continents and ecological niches.
Distribution PatternsSpecific Continent(s)/Ecological Niche(s)Number of Sequence TypesNumber of Isolates
Geographic
In all six continentsAfrica + Asia + Europe + N America + S America + Oceania00
In five continents only 00
In four continents only
Asia + Europe + N America + S America250
Asia + Europe + N America + Oceania111
Other combinations of four continents00
In three continents only
Africa + Asia + Europe356
Asia + Europe + S America18
Asia + N America + S America2186
Asia + N America + Oceania14
Other combinations of three continents00
In two continents only
Africa + Asia12
Africa + Europe13
Africa + S America13
Asia + Europe25
Asia + N America341
Asia + S America342
Asia + Oceania14
Europe + N America24
Europe + S America28
N America + S America595
N America + Oceania26
S America + Oceania116
Other combinations of two continents00
In one continent only
Africa2432
Asia3074
Europe3145
N America75131
S America156331
Oceania2545
Ecological niches
In all three niches 5237
In two niches only
Clinical + Veterinary531
Clinical + Environmental10165
Veterinary + Environmental29
In one niche only
Clinical123255
Veterinary1922
Environmental2469
Table 4. Analysis of molecular variance at the continental level.
Table 4. Analysis of molecular variance at the continental level.
Non-Clone-Corrected Clone-Corrected
dfMSEst.Var%p-ValuedfMSEst.Var%p-Value
Total CGSC
Among continents598.3730.54017%0.001514.2830.1755%0.001
Within continents11962.6462.64983%0.0014173.1003.10095%0.001
Total12013.0453.189 4223.233.275
VGI
Among continents520.4560.53221%0.001510.3290.40815%0.001
Within continents2192.0502.05079%0.0011152.2642.26485%0.001
Total2242.462.582 1202.62.672
VGII
Among continents446.5160.41716%0.00144.7980.1024%0.001
Within continents7132.1782.17884%0.0011902.7232.72396%0.001
Total 7172.422.594 1942.7662.826
VGIII
Among continents133.5650.32112%0.00114.6130.0823%0.066
Within continents2122.3182.31888%0.001712.6572.65797%0.066
Total 2132.462.639 722.682.739
VGIV
Among continents14.2170.1738%0.02212.7720.0573%0.196
Within continents282.0362.03692%0.022192.2102.21097%
Total 292.112.209 202.22.267
df, degrees of freedom; MS, mean square; Est.Var, estimated variance; %, percentage of variance; p-value based on 1000 permutations.
Table 5. Pairwise population comparison at the continental level.
Table 5. Pairwise population comparison at the continental level.
Non-Clone-CorrectedClone-Corrected
Total CGSC
AfricaAsiaEuropeOceaniaN. AmericaAfricaAsiaEuropeOceaniaN. America
Asia0.119 *** 0.093 ***
Europe0.322 ***0.220 *** 0.184 ***0.043 ***
Oceania0.410 ***0.362 ***0.183 *** 0.053 ***0.035 ***0.056 ***
N. America0.307 ***0.255 ***0.110 ***0.061 *** 0.120 ***0.015 ***0.070 ***0.075 ***
S. America0.223 ***0.169 ***0.053 ***0.128 ***0.065 ***0.036 ***0.026 ***0.101 ***0.034 ***0.034 ***
VGI
Asia0.173 *** 0.051 ***
Europe0.095 ***0.210 *** 0.083 ***0.062 ***
Oceania0.453 ***0.307 ***0.491 *** 0.403 ***0.220 ***0.412 ***
N. America0.189 ***0.104 ***0.242 ***0.130 *** 0.153 ***0.037 ***0.165 ***0.123 ***
S. America0.145 ***0.093 ***0.161 ***0.237 ***0.045 ***0.089 ***0.010 ***0.102 ***0.262 ***0.045 ***
VGII
AsiaN/A N/A
EuropeN/A0.191 *** N/A0.009 ***
OceaniaN/A0.320 ***0.388 *** N/A0.015 ***0.061 ***
N. AmericaN/A0.229 ***0.078 ***0.168 *** N/A0.027 ***0.020 ***0.066 ***
S. AmericaN/A0.135 ***0.051 ***0.158 ***0.025 ***N/A0.00 ***0.049 ***0.047 ***0.051 ***
VGIII
AsiaN/A N/A
EuropeN/AN/A N/AN/A
OceaniaN/AN/AN/A N/AN/AN/A
N. AmericaN/AN/AN/AN/A N/AN/AN/AN/A
S. AmericaN/AN/AN/AN/A0.122 ***N/AN/AN/AN/A0.030
VGIV
Asia0.078 * 0.025
EuropeN/AN/A N/AN/A
OceaniaN/AN/AN/A N/AN/AN/A
N. AmericaN/AN/AN/AN/A N/AN/AN/AN/A
S. AmericaN/AN/AN/AN/AN/AN/AN/AN/AN/AN/A
* p < 0.05; *** p < 0.001; N/A, not analyzed due to small sample size.
Table 6. Analysis of molecular variance by ecological niche.
Table 6. Analysis of molecular variance by ecological niche.
Non-Clone-CorrectedClone-Corrected
dfMSEst.Var%p-ValuedfMSEst.Var%p-Value
Australia
Among ecological niches223.80.3222%0.001N/AN/AN/AN/AN/A
Within ecological niches2141.11.1478%0.001N/AN/AN/AN/AN/A
Total2161.361.46 N/AN/AN/A
Brazil
Among ecological niches232.10.2411%0.00122.990.021%0.318
Within ecological niches4001.991.9989%0.001972.662.6699%0.318
Total4022.122.24 992.672.67
Colombia
Among ecological niches228.50.2410%0.00112.730.00%0.502
Within ecological niches3322.22.1690%0.001293.033.03100%0.502
Total33432.31 302.403.03
Taiwan
Among ecological niches219.20.2413%0.001N/AN/AN/AN/AN/A
Within ecological niches2551.581.5887%0.001N/AN/AN/AN/AN/A
Total257 1.82 N/AN/AN/A
USA
Among ecological niches229.60.2610%0.00112.440.00%0.54
Within ecological niches3202.292.2990%0.001452.782.78100%0.54
Total3222.462.56 462.782.78
df, degrees of freedom; MS, mean square; Est.Var, estimated variance; %, percentage of variance; p-value based on 1000 permutations. N/A, not analyzed due to small sample size.
Table 7. Allelic variation among the seven loci used for the multi-locus sequence typing of the C. gattii species complex. The number of shared alleles at each locus between VG lineages.
Table 7. Allelic variation among the seven loci used for the multi-locus sequence typing of the C. gattii species complex. The number of shared alleles at each locus between VG lineages.
GeneGene NameChromosomeLength (bp)Total Allele Number in CGSCNumber of Shared Alleles between VG Lineages
CAP59Capsular-associated protein1556–558951
GPD1Glyceraldehyde-3-phosphate dehydrogenase7544–549616
IGS1Ribosomal RNA intergenic spacer2624–7061104
LAC1Laccase8472–475632
PLB1Phospholipase12532–535461
SOD1Cu, Zn superoxide dismutase5700–7131365
URA5Orotidine monophosphate pyrophosphorylase8637–639591
Table 8. Allelic profiles of sequence types with shared alleles among VG lineages.
Table 8. Allelic profiles of sequence types with shared alleles among VG lineages.
STNo. of Inconsistent ATsNo. of VG Lineages RepresentedOriginally Assigned VG LineageCAP59GPD1IGS1LAC1PLB1SOD1URA5
ST1112VGIV17102139311
ST5312VGI165401353212
ST16712VGII23314187
ST18112VGII832242873
ST22312VGIII42286141315117
ST33512VGIII2912119132822
ST35012VGII51125416166
ST35142VGI4416594131414
ST35512VGII3557418582
ST35612VGI1663554512
ST35712VGII162542457
ST38612VGIV63112136211
ST40012VGII25256211907
ST40263VGIV25114613137311
ST40812VGI7142906511338
ST40912VGII2525621187
ST41912VGIV23421773534
ST43123VGI4445913136815
ST45112VGII86144263
ST50812VGIV23410317443534
ST52423VGI44410213136815
ST52812VGI7156906511338
ST55012VGI1651255323
No., number; Pink represents the VGI lineage, light blue represents the VGII lineage, dark blue represents the VGIII lineage, yellow represents the VGIV lineage, and orange represents STs inconsistently grouped in Figure 1.
Table 9. Summary of genotypic diversity and phylogenetic incompatibility based on the sequence type and VG lineage information assigned by the MLST database. ** p < 0.01; *** p < 0.001.
Table 9. Summary of genotypic diversity and phylogenetic incompatibility based on the sequence type and VG lineage information assigned by the MLST database. ** p < 0.01; *** p < 0.001.
PopulationNumberPhylogenetic Compatibility
(% of 21 Pairs)
Index of Association
Global CGSC56600.82 ***
Global VGI16001.27 ***
Global VGII21800.28 ***
Global VGIII759.5%1.46 ***
Global VGIV11123.8%0.14 **
Africa CGSC3152.4%2.10 ***
Asia CGSC5071.4%1.99 ***
Europe CGSC4676.1%2.40 ***
Oceania CGSC3490.4%3.10 ***
North America CGSC9414.2%1.96 ***
South America CGSC17400.55 ***
Asia VGI2976.1%1.63 ***
Europe VGI2885.7%1.51 ***
Oceania VGI2190.4%1.05 ***
North America VGI2390.4%1.43 ***
North America VGII5919%1.37 ***
South America VGII14801.60 ***
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hitchcock, M.; Xu, J. Global Analyses of Multi-Locus Sequence Typing Data Reveal Geographic Differentiation, Hybridization, and Recombination in the Cryptococcus gattii Species Complex. J. Fungi 2023, 9, 276. https://doi.org/10.3390/jof9020276

AMA Style

Hitchcock M, Xu J. Global Analyses of Multi-Locus Sequence Typing Data Reveal Geographic Differentiation, Hybridization, and Recombination in the Cryptococcus gattii Species Complex. Journal of Fungi. 2023; 9(2):276. https://doi.org/10.3390/jof9020276

Chicago/Turabian Style

Hitchcock, Megan, and Jianping Xu. 2023. "Global Analyses of Multi-Locus Sequence Typing Data Reveal Geographic Differentiation, Hybridization, and Recombination in the Cryptococcus gattii Species Complex" Journal of Fungi 9, no. 2: 276. https://doi.org/10.3390/jof9020276

APA Style

Hitchcock, M., & Xu, J. (2023). Global Analyses of Multi-Locus Sequence Typing Data Reveal Geographic Differentiation, Hybridization, and Recombination in the Cryptococcus gattii Species Complex. Journal of Fungi, 9(2), 276. https://doi.org/10.3390/jof9020276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop