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

AutoESDCas: A Web-Based Tool for the Whole-Workflow Editing Sequence Design for Microbial Genome Editing Based on the CRISPR/Cas System

Genome editing is the basis for the modification of engineered microbes. In the process of genome editing, the design of editing sequences, such as primers and sgRNA, is very important for the accurate positioning of editing sites and efficient sequence editing. The whole process of genome editing involves multiple rounds and types of editing sequence design, while the development of related whole-workflow design tools for high-throughput experimental requirements lags. Here, we propose AutoESDCas, an online tool for the end-to-end editing sequence design for microbial genome editing based on the CRISPR/Cas system. This tool facilitates all types of genetic manipulation covering diverse experimental requirements and design scenarios, enables biologists to quickly and efficiently obtain all editing sequences needed for the entire genome editing process, and empowers high-throughput strain modification. Notably, with its off-target risk assessment function for editing sequences, the usability of the design results is significantly improved. AutoESDCas is freely available at https://autoesdcas.biodesign.ac.cn/with the source code at https://github.com/tibbdc/AutoESDCas/.

Hybrids unleashed: exploring the emergence and genomic insights of pathogenic yeast hybrids

Hybridisation is the crossing of two divergent lineages that give rise to offspring carrying an admixture of both parental genomes. Genome sequencing has revealed that this process is common in the Saccharomycotina, where a growing number of hybrid strains or species, including many pathogenic ones, have been recently described. Hybrids can display unique traits that may drive adaptation to new niches, and some pathogenic hybrids have been shown to have higher prevalence over their parents in human and environmental niches, suggesting a higher fitness and potential to colonise humans. Here, we discuss how hybridisation and its genomic and phenotypic outcomes can shape the evolution of fungal species and may play a role in the emergence of new pathogens.

Engineering of Yarrowia lipolytica as a platform strain for producing adipic acid from renewable resource

There is an increasing demand for bio-based dicarboxylic acids (DCA) as an eco-friendly alternatives to chemically synthesized DCA. Adipic acid, which is not naturally produced by microorganisms, is an essential DCA with significant industrial importance. This study aimed to develop a platform strain using Yarrowia lipolytica for efficient bioconversion of renewable resources into adipic acid. To prevent the complete oxidation of adipic acid, peroxisomal β-oxidation was engineered by selectively disrupting acyl-CoA oxidases. Furthermore, ω-oxidation activity was improved via introducing an additional copy of cytochrome P450 monooxygenase (ALK5) and reductase (CPR1) with fatty alcohol oxidase (FAO1). The production phase used SP92D medium in a two-stage bioconversion process, during which the engineered strain exhibited the highest production level, achieving a remarkable 9.7-fold increase compared to that of the parental strain. To our knowledge, this is the first report demonstrating that engineered Y. lipolytica can produce adipic acid from fatty acid methyl esters.

Recent gene selection and drug resistance underscore clinical adaptation across Candida species

Understanding how microbial pathogens adapt to treatments, humans and clinical environments is key to infer mechanisms of virulence, transmission and drug resistance. This may help improve therapies and diagnostics for infections with a poor prognosis, such as those caused by fungal pathogens, including Candida. Here we analysed genomic variants across approximately 2,000 isolates from six Candida species (C. glabrata, C. auris, C. albicans, C. tropicalis, C. parapsilosis and C. orthopsilosis) and identified genes under recent selection, suggesting a highly complex clinical adaptation. These involve species-specific and convergently affected adaptive mechanisms, such as adhesion. Using convergence-based genome-wide association studies we identified known drivers of drug resistance alongside potentially novel players. Finally, our analyses reveal an important role of structural variants and suggest an unexpected involvement of (para)sexual recombination in the spread of resistance. Our results provide insights on how opportunistic pathogens adapt to human-related environments and unearth candidate genes that deserve future attention.

Tandem gene duplications contributed to high-level azole resistance in a rapidly expanding Candida tropicalis population

Invasive diseases caused by the globally distributed commensal yeast Candida tropicalis are associated with mortality rates of greater than 50%. Notable increases of azole resistance have been observed in this species, particularly within Asia-Pacific regions. Here, we carried out a genetic population study on 1571 global C. tropicalis isolates using multilocus sequence typing (MLST). In addition, whole-genome sequencing (WGS) analysis was conducted on 629 of these strains, comprising 448 clinical invasive strains obtained in this study and 181 genomes sourced from public databases. We found that MLST clade 4 is the predominant azole-resistant clone. WGS analyses demonstrated that dramatically increasing rates of azole resistance are associated with a rapid expansion of cluster AZR, a sublineage of clade 4. Cluster AZR isolates exhibited a distinct high-level azole resistance, which was induced by tandem duplications of the ERG11A395T gene allele. Ty3/gypsy-like retrotransposons were found to be highly enriched in this population. The alarming expansion of C. tropicalis cluster AZR population underscores the urgent need for strategies against growing threats of antifungal resistance.

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.

Origin of fungal hybrids with pathogenic potential from warm seawater environments

Hybridisation is a common event in yeasts often leading to genomic variability and adaptation. The yeast Candida orthopsilosis is a human-associated opportunistic pathogen belonging to the Candida parapsilosis species complex. Most C. orthopsilosis clinical isolates are hybrids resulting from at least four independent crosses between two parental lineages, of which only one has been identified. The rare presence or total absence of parentals amongst clinical isolates is hypothesised to be a consequence of a reduced pathogenicity with respect to their hybrids. Here, we sequence and analyse the genomes of environmental C. orthopsilosis strains isolated from warm marine ecosystems. We find that a majority of environmental isolates are hybrids, phylogenetically closely related to hybrid clinical isolates. Furthermore, we identify the missing parental lineage, thus providing a more complete overview of the genomic evolution of this species. Additionally, we discover phenotypic differences between the two parental lineages, as well as between parents and hybrids, under conditions relevant for pathogenesis. Our results suggest a marine origin of C. orthopsilosis hybrids, with intrinsic pathogenic potential, and pave the way to identify pre-existing environmental adaptations that rendered hybrids more prone than parental lineages to colonise and infect the mammalian host.

Enhancing the Antifungal Efficacy of Fluconazole with a Diterpene: Abietic Acid as a Promising Adjuvant to Combat Antifungal Resistance in Candida spp

The increasing antifungal resistance rates against conventional drugs reveal the urgent need to search for new therapeutic alternatives. In this context, natural bioactive compounds have a critical role in antifungal drug development. Since evidence demonstrates that abietic acid, a diterpene found in Pinus species, has significant antimicrobial properties, this study aimed to evaluate the antifungal activity of abietic acid against Candida spp and its ability to potentiate the activity of fluconazole. Abietic acid was tested both individually and in combination with fluconazole against Candida albicans (CA INCQS 40006), Candida krusei (CK INCQS 40095), and Candida tropicalis (CT INCQS 40042). The microdilution method was used to determine the IC50 and the cell viability curve. Minimum Fungicidal Concentration (MFC) was determined by subculture in a solid medium. The plasma membrane permeability was measured using a fluorescent SYTOX Green probe. While the IC50 of the drugs alone ranged between 1065 and 3255 μg/mL, the IC50 resulting from the combination of abietic acid and fluconazole ranged between 7563 and 160.1 μg/mL. Whether used in combination with fluconazole or isolated, abietic acid exhibited Minimum Fungicidal Concentration (MFC) values exceeding 1024 μg/mL against Candida albicans, Candida krusei and Candida tropicalis. However, it was observed that the antifungal effect of fluconazole was enhanced when used in combination with abietic acid against Candida albicans and Candida tropicalis. These findings suggest that while abietic acid alone has limited inherent antifungal activity, it can enhance the effectiveness of fluconazole, thereby reducing antifungal resistance.

The evolution of anaerobic growth in Saccharomycotina yeasts

Humans rely on the ability of budding yeasts to grow without oxygen in industrial scale fermentations that produce beverages, foods, and biofuels. Oxygen is deeply woven into the energy metabolism and biosynthetic capabilities of budding yeasts. While diverse ecological habitats may provide wide varieties of different carbon and nitrogen sources for yeasts to utilize, there is no direct substitute for molecular oxygen, only a range of availability. Understanding how a small subset of budding yeasts evolved the ability to grow without oxygen could expand the set of useful species in industrial scale fermentations as well as provide insight into the cryptic field of yeast ecology. However, we still do not yet appreciate the full breadth of species that can growth without oxygen, what genes underlie this adaptation, and how these genes have evolved.

Genotyping and susceptibility testing uncovers large azole-resistant Candida tropicalis clade in Alexandria, Egypt

OBJECTIVES

Candida tropicalis is an emerging medically relevant Candida species. The yeast primarily causes opportunistic infections in intensive care units and is highly prevalent in tropical countries. The genetic diversity within this species is high, and nosocomial transmission has been reported. C. tropicalis genotyping of isolates from low- and middle-income countries is underrepresented when compared with that from high-income countries. Also, in Egypt, only limited genotyping has been conducted for C. tropicalis isolates, while antifungal resistance seems to increase, especially against azoles.

METHODS

Antifungal susceptibility testing was performed on 64 C. tropicalis isolates from ICU patients collected from multiple hospitals in Alexandria, Egypt. Genotyping by means of short tandem repeat (STR) and whole genome sequencing (WGS) single nucleotide polymorphism (SNP) analysis was performed.

RESULTS

Using antifungal susceptibility testing, fluconazole resistance was observed in 24 isolates (38%), of which 23 harboured an ERG11 G464S substitution, previously shown to cause resistance in Candida albicans. STR genotyping showed that these 23 isolates were related, forming a distinct resistant clade. WGS SNP analysis subsequently confirmed this genetic relationship, although isolates within this clade differed in at least 429 SNPs, suggesting that these were independently introduced.

CONCLUSION

Overall, STR and WGS SNP analysis of this collection indicates limited C. tropicalis nosocomial transmission in Alexandria, while the presence of this large azole-resistant C. tropicalis clade within this city hampers the treatment of intensive care unit patients.

Impact of pathogen genetics on clinical phenotypes in a population of Talaromyces marneffei from Vietnam

Talaromycosis, a severe and invasive fungal infection caused by Talaromyces marneffei, is difficult to treat and impacts those living in endemic regions of Southeast Asia, India, and China. While 30% of infections result in mortality, our understanding of the genetic basis of pathogenesis for this fungus is limited. To address this, we apply population genomics and genome-wide association study approaches to a cohort of 336 T. marneffei isolates collected from patients who enrolled in the Itraconazole vs Amphotericin B for Talaromycosis trial in Vietnam. We find that isolates from northern and southern Vietnam form two distinct geographical clades, with isolates from southern Vietnam associated with increased disease severity. Leveraging longitudinal isolates, we identify multiple instances of disease relapse linked to unrelated strains, highlighting the potential for multistrain infections. In more frequent cases of persistent talaromycosis caused by the same strain, we identify variants arising over the course of patient infections that impact genes predicted to function in the regulation of gene expression and secondary metabolite production. By combining genetic variant data with patient metadata for all 336 isolates, we identify pathogen variants significantly associated with multiple clinical phenotypes. In addition, we identify genes and genomic regions under selection across both clades, highlighting loci undergoing rapid evolution, potentially in response to external pressures. With this combination of approaches, we identify links between pathogen genetics and patient outcomes and identify genomic regions that are altered during T. marneffei infection, providing an initial view of how pathogen genetics affects disease outcomes.

Atmospheric humidity regulates same-sex mating in Candida albicans through the trehalose and osmotic signaling pathways

Sexual reproduction is prevalent in eukaryotic organisms and plays a critical role in the evolution of new traits and in the generation of genetic diversity. Environmental factors often have a direct impact on the occurrence and frequency of sexual reproduction in fungi. The regulatory effects of atmospheric relative humidity (RH) on sexual reproduction and pathogenesis in plant fungal pathogens and in soil fungi have been extensively investigated. However, the knowledge of how RH regulates the lifecycles of human fungal pathogens is limited. In this study, we report that low atmospheric RH promotes the development of mating projections and same-sex (homothallic) mating in the human fungal pathogen Candida albicans. Low RH causes water loss in C. albicans cells, which results in osmotic stress and the generation of intracellular reactive oxygen species (ROS) and trehalose. The water transporting aquaporin Aqy1, and the G-protein coupled receptor Gpr1 function as cell surface sensors of changes in atmospheric humidity. Perturbation of the trehalose metabolic pathway by inactivating trehalose synthase or trehalase promotes same-sex mating in C. albicans by increasing osmotic or ROS stresses, respectively. Intracellular trehalose and ROS signal the Hog1-osmotic and Hsf1-Hsp90 signaling pathways to regulate the mating response. We, therefore, propose that the cell surface sensors Aqy1 and Gpr1, intracellular trehalose and ROS, and the Hog1-osmotic and Hsf1-Hsp90 signaling pathways function coordinately to regulate sexual mating in response to low atmospheric RH conditions in C. albicans.

Genomic epidemiology of human candidaemia isolates in a tertiary hospital

Invasive candida infections are significant infections that may occur in vulnerable patients with high rates of mortality or morbidity. Drug-resistance rates also appear to be on the rise which further complicate treatment options and outcomes. The aims of this study were to describe the prevalence, molecular epidemiology, and genetic features of Candida bloodstream isolates in a hospital setting. The resistance mechanisms towards the two most commonly administered antifungals, fluconazole and anidulafungin, were determined. Blood culture isolates between 1 January 2018 and 30 June 2021 positive for Candida spp. were included. Susceptibility testing was performed using Etest. Whole-genome-sequencing was performed using Illumina NovaSeq with bioinformatics analysis performed. A total of 203 isolates were sequenced: 56 C. glabrata, 53 C. tropicalis, 44 C. albicans, 36 C. parapsilosis complex (consisting of C. parapsilosis, C. orthopsilosis, and C. metapsilosis), six C. krusei, five C. dubliniensis, and three C. auris. A single cluster of azole-resistant C. tropicalis, and four clusters of C. parapsilosis isolates were observed, suggesting possible transmission occurring over several years. We found 11.3%, and 52.7 % of C. tropicalis and C. parapsilosis, respectively, clustered with other isolates, suggesting exogenous sources may play a significant role of transmission, particularly for C. parapsilosis. The clusters spanned over several years suggesting the possibility of environmental reservoirs contributing to the spread. Limited clonality was seen for C. albicans. Several sequence types appeared to be dominant for C. glabrata, however the SNP differences varied widely, indicating absence of sustained transmission.

Evolution of loss of heterozygosity patterns in hybrid genomes of Candida yeast pathogens

BACKGROUND

Hybrids are chimeric organisms with highly plastic heterozygous genomes that may confer unique traits enabling the adaptation to new environments. However, most evolutionary theory frameworks predict that the high levels of genetic heterozygosity present in hybrids from divergent parents are likely to result in numerous deleterious epistatic interactions. Under this scenario, selection is expected to favor recombination events resulting in loss of heterozygosity (LOH) affecting genes involved in such negative interactions. Nevertheless, it is so far unknown whether this phenomenon actually drives genomic evolution in natural populations of hybrids. To determine the balance between selection and drift in the evolution of LOH patterns in natural yeast hybrids, we analyzed the genomic sequences from fifty-five hybrid strains of the pathogenic yeasts Candida orthopsilosis and Candida metapsilosis, which derived from at least six distinct natural hybridization events.

RESULTS

We found that, although LOH patterns in independent hybrid clades share some level of convergence that would not be expected from random occurrence, there is an apparent lack of strong functional selection. Moreover, while mitosis is associated with a limited number of inter-homeologous chromosome recombinations in these genomes, induced DNA breaks seem to increase the LOH rate. We also found that LOH does not accumulate linearly with time in these hybrids. Furthermore, some C. orthopsilosis hybrids present LOH patterns compatible with footprints of meiotic recombination. These meiotic-like patterns are at odds with a lack of evidence of sexual recombination and with our inability to experimentally induce sporulation in these hybrids.

CONCLUSIONS

Our results suggest that genetic drift is the prevailing force shaping LOH patterns in these hybrid genomes. Moreover, the observed LOH patterns suggest that these are likely not the result of continuous accumulation of sporadic events-as expected by mitotic repair of rare chromosomal breaks-but rather of acute episodes involving many LOH events in a short period of time.

Impact of pathogen genetics on clinical phenotypes in a population of Talaromyces marneffei from Vietnam

Talaromycosis, a severe and invasive fungal infection caused by Talaromyces marneffei , is difficult to treat and impacts those living in endemic regions of southeast Asia, India, and China. While 30% of infections result in mortality, our understanding of the genetic basis of pathogenesis for this fungus is limited. To address this, we apply population genomics and genome wide association study approaches to a cohort of 336 T. marneffei isolates collected from patients who enrolled in the Itraconazole versus Amphotericin B for Talaromycosis (IVAP) trial in Vietnam. We find that isolates from northern and southern Vietnam form two distinct geographical clades, with isolates from southern Vietnam associated with increased disease severity. Leveraging longitudinal isolates, we identify multiple instances of disease relapse linked to unrelated strains, highlighting the potential for multi-strain infections. In more frequent cases of persistent talaromycosis caused by the same strain, we identify variants arising over the course of patient infections that impact genes predicted to function in the regulation of gene expression and secondary metabolite production. By combining genetic variant data with patient metadata for all 336 isolates, we identify pathogen variants significantly associated with multiple clinical phenotypes. In addition, we identify genes and genomic regions under selection across both clades, highlighting loci undergoing rapid evolution, potentially in response to external pressures. With this combination of approaches, we identify links between pathogen genetics and patient outcomes and identify genomic regions that are altered during T. marneffei infection, providing an initial view of how pathogen genetics affects disease outcomes.

Short Tandem Repeat Genotyping and Antifungal Susceptibility Testing of Latin American Candida tropicalis Isolates

Candida tropicalis is emerging as one of the most common Candida species causing opportunistic infections in Latin America. Outbreak events caused by C. tropicalis were reported, and antifungal resistant isolates are on the rise. In order to investigate population genomics and look into antifungal resistance, we applied a short tandem repeat (STR) genotyping scheme and antifungal susceptibility testing (AFST) to 230 clinical and environmental C. tropicalis isolates from Latin American countries. STR genotyping identified 164 genotypes, including 11 clusters comprised of three to seven isolates, indicating outbreak events. AFST identified one isolate as anidulafungin-resistant and harboring a FKS1 S659P substitution. Moreover, we identified 24 clinical and environmental isolates with intermediate susceptibility or resistance to one or more azoles. ERG11 sequencing revealed each of these isolates harboring a Y132F and/or Y257H/N substitution. All of these isolates, except one, were clustered together in two groups of closely related STR genotypes, with each group harboring distinct ERG11 substitutions. The ancestral C. tropicalis strain of these isolates likely acquired the azole resistance-associated substitutions and subsequently spread across vast distances within Brazil. Altogether, this STR genotyping scheme for C. tropicalis proved to be useful for identifying unrecognized outbreak events and better understanding population genomics, including the spread of antifungal-resistant isolates.

Development and Application of a Short Tandem Repeat Multiplex Typing Assay for Candida tropicalis

Candida tropicalis is a clinically important yeast that causes candidemia in humans with a high mortality rate. The yeast primarily infects immunocompromised patients, and causes outbreaks in health care facilities. Antifungal resistant isolates have been reported. We developed a short tandem repeat (STR) typing scheme for C. tropicalis to enable fast, cost-effective, and high-resolution genotyping. For the development of the typing scheme, 6 novel STR markers were selected, combined into 2 multiplex PCRs. In total, 117 C. tropicalis isolates were typed, resulting in the identification of 104 different genotypes. Subsequently, the outcome of STR typing of 10 isolates was compared to single nucleotide polymorphism (SNP) calling from whole-genome sequencing (WGS). Isolates with more than 111 SNPs were differentiated by the typing assay. Two isolates, which were identical according to SNP analysis, were separated by STR typing in 1 marker. To test specificity, the STR typing was applied to 15 related yeast species, and we found no amplification of these targets. For reproducibility testing, 2 isolates were independently typed five times, which showed identical results in each experiment. In summary, we developed a reliable and multiplex STR genotyping for C. tropicalis, which was found to correlate well to SNP calling by WGS. WGS analysis from and extensive collection of isolates is required to establish the precise resolution of this STR assay. IMPORTANCE Candida tropicalis frequently causes candidemia in immunocompromised patients. C. tropicalis infections have a high mortality rate, and the yeast is able to cause outbreaks in health care facilities. Further, antifungal resistant isolates are on the rise. Genotyping is necessary to investigate potential outbreaks. Here, we developed and applied a STR genotyping scheme in order to rapidly genotype isolates with a high-resolution. WGS SNP outcomes were highly comparable with STR typing results. Altogether, we developed a rapid, high-resolution, and specific STR genotyping scheme for C. tropicalis.

Multilocus sequence typing (MLST) analysis reveals many novel genotypes and a high level of genetic diversity in Candida tropicalis isolates from Italy and Africa

BACKGROUND

Candida tropicalis is a human pathogenic yeast frequently isolated in Latin America and Asian-Pacific regions, although recent studies showed that it is also becoming increasingly widespread throughout several African and south-European countries. Nevertheless, relatively little is known about its global patterns of genetic variation as most of existing multilocus sequence typing (MLST) data come from Asia and there are no genotyped African isolates.

OBJECTIVES

We report detailed genotyping data from a large set of C. tropicalis isolates recovered from different clinical sources in Italy, Egypt and Cameroon in order to expand the allele/genotype library of MLST database (https://pubmlst.org/ctropicalis), and to explore the genetic diversity in this species.

METHODS

A total of 103 C. tropicalis isolates were genotyped using the MLST scheme developed for this species. All isolates were also tested for in vitro susceptibility to various antifungals to assess whether certain genotypes were associated with drug-resistance.

RESULTS AND CONCLUSIONS

A total of 104 different alleles were detected across the MLST-loci investigated. The allelic diversity found at these loci resulted in 51 unique MLST genotypes of which 36 (70.6%) were novel. Global optimal eBURST analysis identified 18 clonal complexes (CCs) and confirm the existence of a specific Italian-cluster (CC36). Three CCs were also statistically associated with fluconazole resistance, which was elevated in Cameroon and Egypt. Our data show high genetic diversity in our isolates suggesting that the global population structure of C. tropicalis is still poorly understood. Moreover, its clinical impact in Italy, Egypt and Cameroon appears to be relevant and should be carefully considered.

Clonality, inbreeding, and hybridization in two extremotolerant black yeasts

BACKGROUND

The great diversity of lifestyles and survival strategies observed in fungi is reflected in the many ways in which they reproduce and recombine. Although a complete absence of recombination is rare, it has been reported for some species, among them 2 extremotolerant black yeasts from Dothideomycetes: Hortaea werneckii and Aureobasidium melanogenum. Therefore, the presence of diploid strains in these species cannot be explained as the product of conventional sexual reproduction.

RESULTS

Genome sequencing revealed that the ratio of diploid to haploid strains in both H. werneckii and A. melanogenum is about 2:1. Linkage disequilibrium between pairs of polymorphic loci and a high degree of concordance between the phylogenies of different genomic regions confirmed that both species are clonal. Heterozygosity of diploid strains is high, with several hybridizing genome pairs reaching the intergenomic distances typically seen between different fungal species. The origin of diploid strains collected worldwide can be traced to a handful of hybridization events that produced diploids, which were stable over long periods of time and distributed over large geographic areas.

CONCLUSIONS

Our results, based on the genomes of over 100 strains of 2 black yeasts, show that although they are clonal, they occasionally form stable and highly heterozygous diploid intraspecific hybrids. The mechanism of these apparently rare hybridization events, which are not followed by meiosis or haploidization, remains unknown. Both extremotolerant yeasts, H. werneckii and even more so A. melanogenum, a close relative of the intensely recombining and biotechnologically relevant Aureobasidium pullulans, provide an attractive model for studying the role of clonality and ploidy in extremotolerant fungi.

Evolution and molecular bases of reproductive isolation

The most challenging problem in speciation research is disentangling the relative strength and order in which different reproductive barriers evolve. Here, we review recent developments in the study of reproductive isolation in yeasts. With over a thousand genome-sequenced isolates readily available for testing the viability, sterility, and fitness of both intraspecies and interspecies hybrid crosses, Saccharomyces yeasts are an ideal model to study such fundamental questions. Our survey demonstrates that, while chromosomal-level mutations are widespread at the intraspecific level, anti-recombination-driven chromosome missegregation is the primary reproductive barrier between species. Finally, despite their strength, all of these postzygotic barriers can be resolved through the asexual life history of hybrids.

Whole Genome Sequencing Shows Genetic Diversity, as Well as Clonal Complex and Gene Polymorphisms Associated with Fluconazole Non-Susceptible Isolates of Candida tropicalis

Resistance to azoles in Candida tropicalis is increasing and may be mediated by genetic characteristics. Using whole genome sequencing (WGS), we examined the genetic diversity of 82 bloodstream C. tropicalis isolates from two countries and one ATCC strain in a global context. Multilocus sequence typing (MLST) and single nucleotide polymorphism (SNP)-based phylogenies were generated. Minimum inhibitory concentrations (MIC) for antifungal agents were determined using Sensititre YeastOne YO10. Eleven (13.2%) isolates were fluconazole-resistant and 17 (20.5%) were classified as fluconazole-non susceptible (FNS). Together with four Canadian isolates, the genomes of 12 fluconazole-resistant (18 FNS) and 69 fluconazole-susceptible strains were examined for gene mutations associated with drug resistance. Fluconazole-resistant isolates contained a mean of 56 non-synonymous SNPs per isolate in contrast to 36 SNPs in fluconazole-susceptible isolates (interquartile range [IQR] 46−59 vs. 31−48 respectively; p < 0.001). Ten of 18 FNS isolates contained missense ERG11 mutations (amino acid substitutions S154F, Y132F, Y257H). Two echinocandin-non susceptible isolates had homozygous FKS1 mutations (S30P). MLST identified high genetic diversity with 61 diploid sequence types (DSTs), including 53 new DSTs. All four isolates in DST 773 were fluconazole-resistant within clonal complex 2. WGS showed high genetic variation in invasive C. tropicalis; azole resistance was distributed across different lineages but with DST 773 associated with in vitro fluconazole resistance.

The emerging threat antifungal-resistant Candida tropicalis in humans, animals, and environment

Antifungal resistance in humans, animals, and the environment is an emerging problem. Among the different fungal species that can develop resistance, Candida tropicalis is ubiquitous and causes infections in animals and humans. In Asia and some Latin American countries, C. tropicalis is among the most common species related to candidemia, and mortality rates are usually above 40%. Fluconazole resistance is especially reported in Asian countries and clonal spread in humans and the environment has been investigated in some studies. In Brazil, high rates of azole resistance have been found in animals and the environment. Multidrug resistance is still rare, but recent reports of clinical multidrug-resistant isolates are worrisome. The molecular apparatus of antifungal resistance has been majorly investigated in clinical C. tropicalis isolates, revealing that this species can develop resistance through the conjunction of different adaptative mechanisms. In this review article, we summarize the main findings regarding antifungal resistance and Candida tropicalis through an "One Health" approach.

Experimental evolution of drug resistance in human fungal pathogens

Experimental evolution in vitro is a powerful tool to uncover the factors that contribute to resistance evolution and understand the genetic basis of adaptation. Here, we discuss recent experimental evolution studies from human fungal pathogens. We synthesize the results to highlight the common threads that influence resistance acquisition. The picture that emerges is that drug resistance consistently appears readily and rapidly. Mutations are often found in an overlapping set of genes and genetic pathways known to be involved in drug resistance, including whole or partial chromosomal aneuploidy. The likelihood of acquiring resistance and cross-resistance between drugs seems to be influenced by the specific drug (not just drug class), level of drug, and strain genetic background. We discuss open questions, such as the potential for increases in drug tolerance to evolve in static drugs. We highlight opportunities to use this framework to probe how different factors influence the rate and nature of adaptation to antifungal drugs in fungal microbes through a call for increased reporting on all replicates that were evolved, not just those that acquired resistance.

The teenage years of yeast population genomics trace history, admixing and getting wilder

Population genomics studies the evolutionary processes that shape intraspecies genetic variations. In this review, I explore the insights into yeast-population genomics that have emerged from recent advances in sequencing. Genomes of the model Saccharomyces cerevisiae and many new yeast species from around the world are being used to address various aspects of population biology, including geographical origin, the level of introgression, domestication signatures, and outcrossing frequency. New long-read sequencing has enabled a greater capacity to quantify these variations at a finer resolution from complete de novo genomes at the population scale to phasing subgenomes of different origins. These resources provide a platform to dissect the relationship between phenotypes across environmental niches.

Genome plasticity in Candida albicans: A cutting-edge strategy for evolution, adaptation, and survival

Candida albicans is the most implicated fungal species that grows as a commensal or opportunistic pathogen in the human host. It is associated with many life-threatening infections, especially in immunocompromised persons. The genome of Candida albicans is very flexible and can withstand a wide assortment of variations in a continuously changing environment. Thus, genome plasticity is central to its adaptation and has long been of considerable interest. C. albicans has a diploid heterozygous genome that is highly dynamic and can display variation from small to large scale chromosomal rearrangement and aneuploidy, which have implications in drug resistance, virulence, and pathogenicity. This review presents an up-to-date overview of recent genomic studies involving C. albicans. It discusses the accumulating evidence that shows how mitotic recombination events, ploidy dynamics, aneuploidy, and loss of heterozygosity (LOH) influence evolution, adaptation, and survival in C. albicans. Understanding the factors that affect the genome is crucial for a proper understanding of species and rapid development and adjustment of therapeutic strategies to mitigate their spread.

Genome analysis of five recently described species of the CUG-Ser clade uncovers Candida theae as a new hybrid lineage with pathogenic potential in the Candida parapsilosis species complex

Candida parapsilosis species complex comprises three important pathogenic species: Candida parapsilosis sensu stricto, Candida orthopsilosis and Candida metapsilosis. The majority of C. orthopsilosis and all C. metapsilosis isolates sequenced thus far are hybrids, and most of the parental lineages remain unidentified. This led to the hypothesis that hybrids with pathogenic potential were formed by the hybridization of non-pathogenic lineages that thrive in the environment. In a search for the missing hybrid parentals, and aiming to get a better understanding of the evolution of the species complex, we sequenced, assembled and analysed the genome of five close relatives isolated from the environment: Candida jiufengensis, Candida pseudojiufengensis, Candida oxycetoniae, Candida margitis and Candida theae. We found that the linear conformation of mitochondrial genomes in Candida species emerged multiple times independently. Furthermore, our analyses discarded the possible involvement of these species in the mentioned hybridizations, but identified C. theae as an additional hybrid in the species complex. Importantly, C. theae was recently associated with a case of infection, and we also uncovered the hybrid nature of this clinical isolate. Altogether, our results reinforce the hypothesis that hybridization is widespread among Candida species, and potentially contributes to the emergence of lineages with opportunistic pathogenic behaviour.

Yeasts from temperate forests

Yeasts are ubiquitous in temperate forests. While this broad habitat is well-defined, the yeasts inhabiting it and their life cycles, niches, and contributions to ecosystem functioning are less understood. Yeasts are present on nearly all sampled substrates in temperate forests worldwide. They associate with soils, macroorganisms, and other habitats, and no doubt contribute to broader ecosystem-wide processes. Researchers have gathered information leading to hypotheses about yeasts' niches and their life cycles based on physiological observations in the laboratory as well as genomic analyses, but the challenge remains to test these hypotheses in the forests themselves. Here we summarize the habitat and global patterns of yeast diversity, give some information on a handful of well-studied temperate forest yeast genera, discuss the various strategies to isolate forest yeasts, and explain temperate forest yeasts' contributions to biotechnology. We close with a summary of the many future directions and outstanding questions facing researchers in temperate forest yeast ecology. Yeasts present an exciting opportunity to better understand the hidden world of microbial ecology in this threatened and global habitat.

Identification of a novel Candida metapsilosis isolate reveals multiple hybridization events

Candida metapsilosis is a member of the Candida parapsilosis species complex, a group of opportunistic human pathogens. Of all the members of this complex, C. metapsilosis is the least virulent, and accounts for a small proportion of invasive Candida infections. Previous studies established that all C. metapsilosis isolates are hybrids, originating from a single hybridization event between two lineages, parent A and parent B. Here, we use MinION and Illumina sequencing to characterize a C. metapsilosis isolate that originated from a separate hybridization. One of the parents of the new isolate is very closely related to parent A. However, the other parent (parent C) is not the same as parent B. Unlike C. metapsilosis AB isolates, the C. metapsilosis AC isolate has not undergone introgression at the mating type-like locus. In addition, the A and C haplotypes are not fully collinear. The C. metapsilosis AC isolate has undergone loss of heterozygosity with a preference for haplotype A, indicating that this isolate is in the early stages of genome stabilization.

Plasmid-Based CRISPR-Cas9 Editing in Multiple Candida Species

CRISPR-Cas9 technology radically changed the approach to genetic manipulation of both medically and industrially relevant Candida species, as attested by the ever-increasing number of applications to the study of pathogenesis, drug resistance, gene expression, and host pathogen interaction and drug discovery. Here, we describe the use of plasmid-based systems for high efficiency CRISPR-Cas9 gene editing into any strain of four non-albicans Candida species, namely, Candida parapsilosis, Candida orthopsilosis, Candida metapsilosis, and Candida tropicalis. The plasmids pCP-tRNA and pCT-tRNA contain all the elements necessary for expressing the CRISPR-Cas9 machinery, and they can be used in combination with a repair template for disrupting gene function by insertion of a premature stop codon or by gene deletion. The plasmids are easily lost in the absence of selection, allowing scarless gene editing and minimizing detrimental effects of prolonged Cas9 expression.

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.

Fungi as mediators linking organisms and ecosystems

Fungi form a major and diverse component of most ecosystems on Earth. They are both micro and macroorganisms with high and varying functional diversity as well as great variation in dispersal modes. With our growing knowledge of microbial biogeography, it has become increasingly clear that fungal assembly patterns and processes differ from other microorganisms such as bacteria, but also from macroorganisms such as plants. The success of fungi as organisms and their influence on the environment lies in their ability to span multiple dimensions of time, space, and biological interactions, that is not rivalled by other organism groups. There is also growing evidence that fungi mediate links between different organisms and ecosystems, with the potential to affect the macroecology and evolution of those organisms. This suggests that fungal interactions are an ecological driving force, interconnecting different levels of biological and ecological organisation of their hosts, competitors, and antagonists with the environment and ecosystem functioning. Here we review these emerging lines of evidence by focusing on the dynamics of fungal interactions with other organism groups across various ecosystems. We conclude that the mediating role of fungi through their complex and dynamic ecological interactions underlie their importance and ubiquity across Earth's ecosystems.

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.

Aneuploidy Underlies Tolerance and Cross-Tolerance to Drugs in Candida parapsilosis

Candida species are the most common human fungal pathogens worldwide. Although C. albicans remains the predominant cause of candidiasis, infections caused by non-albicans Candida species, including C. parapsilosis, are increasing. In C. albicans, genome plasticity has been shown to be a prevalent strategy of adaptation to stresses. However, the role of aneuploidy in C. parapsilosis is largely unknown. In this study, we found that six different aneuploid karyotypes conferred adaptation to the endoplasmic reticulum stress inducer tunicamycin (TUN) in C. parapsilosis. Interestingly, a specific aneuploidy including trisomy of chromosome 6 (Chr6x3) also enabled cross-tolerance to aureobasidin A (AbA), a sphingolipid biosynthesis inhibitor. Consistent with this, selection on AbA identified adaptors with three different aneuploid karyotypes, including Chr6x3, which also enabled cross-tolerance to both AbA and TUN. Therefore, as in other Candida species, recurrent aneuploid karyotypes enable the adaptation of C. parapsilosis to specific stresses, and specific aneuploidies enable cross-adaptation to different stresses.

IMPORTANCECandida parapsilosis is an emerging human fungal pathogen, especially prevalent in neonates. Aneuploidy, having uneven numbers of chromosomes, is a well-known mechanism for adapting to stress in Candida albicans, the most common human fungal pathogen. In this study, we exposed C. parapsilosis to two very different drugs and selected for rare cells that grew in one of the drugs. We found that the majority of isolates that grew in the drugs had acquired an extra copy of one of several aneuploid chromosomes and that specific aneuploid chromosomes appeared in several independent cell clones. Importantly, an extra copy of chromosome 6 was detected following selection in either one of the drugs, and this extra chromosome conferred the ability to grow in both drugs, a property called cross-adaptation, or cross-tolerance. Thus, this study highlights the genome plasticity of C. parapsilosis and the ability of an extra copy of a single chromosome to promote cell growth in the presence of more than one drug.

Continual Decline in Azole Susceptibility Rates in Candida tropicalis Over a 9-Year Period in China

Background

There have been reports of increasing azole resistance in Candida tropicalis, especially in the Asia-Pacific region. Here we report on the epidemiology and antifungal susceptibility of C. tropicalis causing invasive candidiasis in China, from a 9-year surveillance study.

Methods

From August 2009 to July 2018, C. tropicalis isolates (n = 3702) were collected from 87 hospitals across China. Species identification was carried out by mass spectrometry or rDNA sequencing. Antifungal susceptibility was determined by Clinical and Laboratory Standards Institute disk diffusion (CHIF-NET10-14, n = 1510) or Sensititre YeastOne (CHIF-NET15-18, n = 2192) methods.

Results

Overall, 22.2% (823/3702) of the isolates were resistant to fluconazole, with 90.4% (744/823) being cross-resistant to voriconazole. In addition, 16.9 (370/2192) and 71.7% (1572/2192) of the isolates were of non-wild-type phenotype to itraconazole and posaconazole, respectively. Over the 9 years of surveillance, the fluconazole resistance rate continued to increase, rising from 5.7 (7/122) to 31.8% (236/741), while that for voriconazole was almost the same, rising from 5.7 (7/122) to 29.1% (216/741), with no significant statistical differences across the geographic regions. However, significant difference in fluconazole resistance rate was noted between isolates cultured from blood (27.2%, 489/1799) and those from non-blood (17.6%, 334/1903) specimens (P-value < 0.05), and amongst isolates collected from medical wards (28.1%, 312/1110) versus intensive care units (19.6%, 214/1092) and surgical wards (17.9%, 194/1086) (Bonferroni adjusted P-value < 0.05). Although echinocandin resistance remained low (0.8%, 18/2192) during the surveillance period, it was observed in most administrative regions, and one-third (6/18) of these isolates were simultaneously resistant to fluconazole.

Conclusion

The continual decrease in the rate of azole susceptibility among C. tropicalis strains has become a nationwide challenge in China, and the emergence of multi-drug resistance could pose further threats. These phenomena call for effective efforts in future interventions.