A coupled process of same- and opposite-sex mating generates polyploidy and genetic diversity in Candida tropicalis

Hotspot mutations and genomic expansion of ERG11 are major mechanisms of azole resistance in environmental and human commensal isolates of Candida tropicalis

OBJECTIVES

Infections caused by azole-resistant Candida tropicalis strains are increasing in clinical settings. The reason for this epidemical change and the mechanisms of C. tropicalis azole resistance are not fully understood.

METHODS

In this study, we performed biological and genomic analyses of 239 C. tropicalis strains, including 115 environmental and 124 human commensal isolates.

RESULTS

Most (99.2%) of the isolates had a baseline diploid genome. The strains from both environmental and human niches exhibit similar abilities to survive under stressful conditions and produce secreted aspartic proteases. However, the human commensal isolates exhibited a stronger ability to filament than the environmental strains. We found that 19 environmental isolates (16.5%) and 24 human commensal isolates (19.4%) were resistant to fluconazole. Of the fluconazole-resistant strains, 37 isolates (86.0%) also exhibited cross-resistance to voriconazole. Whole-genome sequencing and phylogenetic analyses revealed that both environmental and commensal isolates were widely distributed in a number of genetic clusters, but the two populations exhibited a close genetic association. The majority of fluconazole-resistant isolates were clustered within a single clade (X).

CONCLUSIONS

The combination of hotspot mutations (Y132F and S154F) and genomic expansion of ERG11, which encodes the azole target lanosterol 14-α-demethylase and represents a major target of azole drugs, was a major mechanism for the development of azole resistance. The isolates carrying both hotspot mutations and genomic expansion of ERG11 exhibited cross-resistance to fluconazole and voriconazole. Moreover, the azole-resistant isolates from both the environmental and human commensal niches showed similar genotypes.

In Vivo Microevolutionary Analysis of a Fatal Case of Rhinofacial and Disseminated Mycosis Due to Azole-Drug-Resistant Candida Species

Ten Candida species strains were isolated from the first known fatal case of rhinofacial and rhino-orbital-cerebral candidiasis. Among them, five strains of Candida parapsilosis complex were isolated during the early stage of hospitalization, while five strains of Candida tropicalis were isolated in the later stages of the disease. Using whole-genome sequencing, we distinguished the five strains of C. parapsilosis complex as four Candida metapsilosis strains and one Candida parapsilosis strain. Antifungal susceptibility testing showed that the five strains of C. parapsilosis complex were susceptible to all antifungal drugs, while five C. tropicalis strains had high minimum inhibitory concentrations to azoles, whereas antifungal-drug resistance gene analysis revealed the causes of azole resistance in such strains. For the first time, we analyzed the microevolutionary characteristics of pathogenic fungi in human hosts and inferred the infection time and parallel evolution of C. tropicalis strains. Molecular clock analysis revealed that azole-resistant C. tropicalis infection occurred during the first round of therapy, followed by divergence via parallel evolution in vivo. The presence/absence variations indicated a potential decrease in the virulence of genomes in strains isolated following antifungal drug treatment, despite the absence of observed clinical improvement in the conditions of the patient. These results suggest that genomic analysis could serve as an auxiliary tool in guiding clinical diagnosis and treatment.

Ploidy Variation and Spontaneous Haploid-Diploid Switching of Candida glabrata Clinical Isolates

The human fungal pathogen Candida glabrata is phylogenetically closely related to Saccharomyces cerevisiae, a model eukaryotic organism. Unlike S. cerevisiae, which has both haploid and diploid forms and a complete sexual cycle, C. glabrata has long been considered a haploid and asexual species. In this study, we analyzed the ploidy states of 500 clinical isolates of C. glabrata from four Chinese hospitals and found that approximately 4% of the isolates were in or able to spontaneously switch to an aneuploid (genomic DNA, 1N-2N), diploid (2N), or hyperdiploid (>2N) form under in vivo or in vitro conditions. Stable diploid cells were identified in 3% of the isolates (15/500). Of particular interest, one clinical strain existed only in the diploid form. Multilocus sequence typing (MLST) assays revealed two major genetic clusters (A and B) of C. glabrata isolates. Most of the isolates (70%) from China belonged to the A cluster, whereas most of the isolates from other countries (such as Iran, Japan, United States, and European countries) belonged to the B cluster. Further investigation indicated that C. glabrata cells of different ploidy forms differed in a number of respects, including morphologies, antifungal susceptibility, virulence, and global gene expression profiles. Additionally, C. glabrata could undergo spontaneous switching between the diploid and haploid forms under both in vitro and in vivo conditions. Given the absence of an apparent sexual phase, one would expect that the ploidy shifts could function as an alternative strategy that promotes genetic diversity and benefits the ability of the fungus to rapidly adapt to the changing environment. IMPORTANCE The human fungal pathogen Candida glabrata has long been thought to be a haploid organism. Here, we report the population structure and ploidy states of 500 clinical isolates of C. glabrata from China. To our surprise, we found that the ploidy of a subset of clinical isolates varied dramatically. Some isolates were in or able to switch to an aneuploid, diploid, or hyperdiploid form. C. glabrata cells with different ploidy differed in a number of biological respects, including morphology, antifungal susceptibility, virulence, and global gene expression profile. Given the absence of an apparent sexual phase in this fungus, we propose that ploidy switching could be a strategy for rapid adaptation to environmental changes and could function as an alternative to sexual reproduction.

Epistatic genetic interactions govern morphogenesis during sexual reproduction and infection in a global human fungal pathogen

Cellular development is orchestrated by evolutionarily conserved signaling pathways, which are often pleiotropic and involve intra- and interpathway epistatic interactions that form intricate, complex regulatory networks. Cryptococcus species are a group of closely related human fungal pathogens that grow as yeasts yet transition to hyphae during sexual reproduction. Additionally, during infection they can form large, polyploid titan cells that evade immunity and develop drug resistance. Multiple known signaling pathways regulate cellular development, yet how these are coordinated and interact with genetic variation is less well understood. Here, we conducted quantitative trait locus (QTL) analyses of a mapping population generated by sexual reproduction of two parents, only one of which is unisexually fertile. We observed transgressive segregation of the unisexual phenotype among progeny, as well as a large-cell phenotype under mating-inducing conditions. These large-cell progeny were found to produce titan cells both in vitro and in infected animals. Two major QTLs and corresponding quantitative trait genes (QTGs) were identified: RIC8 (encoding a guanine-exchange factor) and CNC06490 (encoding a putative Rho-GTPase activator), both involved in G protein signaling. The two QTGs interact epistatically with each other and with the mating-type locus in phenotypic determination. These findings provide insights into the complex genetics of morphogenesis during unisexual reproduction and pathogenic titan cell formation and illustrate how QTL analysis can be applied to identify epistasis between genes. This study shows that phenotypic outcomes are influenced by the genetic background upon which mutations arise, implicating dynamic, complex genotype-to-phenotype landscapes in fungal pathogens and beyond.

Parasexuality of Candida Species

While most fungi have the ability to reproduce sexually, multiple independent lineages have lost meiosis and developed parasexual cycles in its place. Emergence of parasexual cycles is particularly prominent in medically relevant fungi from the CUG paraphyletic group of Candida species. Since the discovery of parasex in C. albicans roughly two decades ago, it has served as the model for Candida species. Importantly, parasex in C. albicans retains hallmarks of meiosis including genetic recombination and chromosome segregation, making it a potential driver of genetic diversity. Furthermore, key meiotic genes play similar roles in C. albicans parasex and highlights parallels between these processes. Yet, the evolutionary role of parasex in Candida adaptation and the extent of resulting genotypic and phenotypic diversity remain as key knowledge gaps in this facultative reproductive program. Here, we present our current understanding of parasex, the mechanisms governing its regulation, and its relevance to Candida biology.

Phenotypic Switching and Filamentation in Candida haemulonii, an Emerging Opportunistic Pathogen of Humans

Phenotypic plasticity is a common strategy adopted by fungal pathogens to adapt to diverse host environments. Candida haemulonii is an emerging multidrug-resistant human pathogen that is closely related to Candida auris. Until recently, it was assumed that C. haemulonii is incapable of phenotypic switching or filamentous growth. In this study, we report the identification of three distinct phenotypes in C. haemulonii: white, pink, and filament. The white and pink phenotypes differ in cellular size, colony morphology, and coloration on phloxine B- or CuSO₄-containing agar. Switching between the white and pink cell types is heritable and reversible and is referred to as “the primary switching system.” The additional switch phenotype, filament, has been identified and exhibits obviously filamentous morphology when grown on glycerol-containing medium. Several unique characteristics of the filamentous phenotype suggest that switching from or to this phenotype poses as a second yeast-filament switching system. The yeast-filament switch is nonheritable and temperature-dependent. Low temperatures favor the filamentous phenotype, whereas high temperatures promote filament-yeast transition. We further demonstrated that numerous aspects of the distinct cell types differ in numerous biological aspects, including their high temperature response, specific gene expression, CuSO₄ tolerance, secreted aspartyl protease (SAP) activity, and virulence. Therefore, transition among the three phenotypes could enable C. haemulonii to rapidly adapt to, survive, and thrive in certain host niches, thereby contributing to its virulence.

IMPORTANCE The capacity to switch between distinct cell types, known as phenotypic switching, is a common strategy adopted by Candida species to adapt to diverse environments. Despite considerable studies on phenotypic plasticity of various Candida species, Candida haemulonii is considered to be incapable of phenotypic switching or filamentous growth. Here, we report and describe filamentation and three distinct phenotypes (white, pink, and filament) in C. haemulonii. The three cell types differ in cellular and colony appearance, gene expression profiles, CuSO₄ tolerance, and virulence. C. haemulonii cells switch heritably and reversibly between white and pink cell types, which is referred to as the “primary switching system.” Switching between pink and filamentous phenotypes is nonheritable and temperature-dependent, representing a second switching system. As in other Candida species, switching among distinct morphological types may provide C. haemulonii with phenotypic plasticity for rapid responses to the changing host environment, and may contribute to its virulence.

Is Natural Population of Candida tropicalis Sexual, Parasexual, and/or Asexual?

Candida tropicalis is one of the most common opportunistic yeast pathogens of humans, especially prevalent in tropical and subtropical regions. This yeast has broad ecological distributions, can be found in both terrestrial and aquatic ecosystems, including being associated with a diversity of trees, animals, and humans. Evolutionary theory predicts that organisms thriving in diverse ecological niches likely have efficient mechanisms to generate genetic diversity in nature. Indeed, abundant genetic variations have been reported in natural populations (both environmental and clinical) of C. tropicalis. However, at present, our understanding on how genetic diversity is generated in natural C. tropicalis population remains controversial. In this paper, I review the current understanding on the potential modes of reproduction in C. tropicalis. I describe expectations of the three modes of reproduction (sexual, parasexual, and asexual) and compare them with the observed genotypic variations in natural populations. Though sexual and parasexual reproduction cannot be excluded, the analyses suggest asexual reproduction alone could explain all the observations reported so far. The results here have implications for understanding the evolution and epidemiology of C. tropicalis and other related human fungal pathogens.

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

The population genetics of ploidy change in unicellular fungi

Changes in ploidy are a significant type of genetic variation, describing the number of chromosome sets per cell. Ploidy evolves in natural populations, clinical populations, and lab experiments, particularly in fungi. Despite a long history of theoretical work on this topic, predicting how ploidy will evolve has proven difficult, as it is often unclear why one ploidy state outperforms another. Here, we review what is known about contemporary ploidy evolution in diverse fungal species through the lens of population genetics. As with typical genetic variants, ploidy evolution depends on the rate that new ploidy states arise by mutation, natural selection on alternative ploidy states, and random genetic drift. However, ploidy variation also has unique impacts on evolution, with the potential to alter chromosomal stability, the rate and patterns of point mutation, and the nature of selection on all loci in the genome. We discuss how ploidy evolution depends on these general and unique factors and highlight areas where additional experimental evidence is required to comprehensively explain the ploidy transitions observed in the field and the lab.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.