Genomic Diversity across Candida auris Clinical Isolates Shapes Rapid Development of Antifungal Resistance In Vitro and In Vivo

Genome-wide analysis of in vivo-evolved Candida auris reveals multidrug-resistance mechanisms

Candida auris, an emerging and multidrug-resistant fungal pathogen, has led to numerous outbreaks in China. While the resistance mechanisms against azole and amphotericin B have been studied, the development of drug resistance in this pathogen remains poorly understood, particularly in in vivo-generated drug-resistant strains. This study employed pathogen whole-genome sequencing to investigate the epidemiology and drug-resistance mutations of C. auris using 16 strains isolated from two patients. Identification was conducted through Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, and antimicrobial susceptibilities were assessed using broth microdilution and Sensititre YeastOne YO10. Whole-genome sequencing revealed that all isolates belonged to the South Asian lineage, displaying genetic heterogeneity. Despite low genetic variability among patient isolates, notable mutations were identified, including Y132F in ERG11 and A585S in TAC1b, likely linked to increased fluconazole resistance. Strains from patient B also carried F214L in TAC1b, resulting in a consistent voriconazole minimum inhibitory concentration of 4 µg/mL across all isolates. Furthermore, a novel frameshift mutation in the SNG1 gene was observed in amphotericin B-resistant isolates compared to susceptible ones. Our findings suggest the potential transmission of C. auris and emphasize the need to explore variations related to antifungal resistance. This involves analyzing genomic mutations and karyotypes, especially in vivo, to compare sensitive and resistant strains. Further monitoring and validation efforts are crucial for a comprehensive understanding of the mechanisms of drug resistance in C. auris.

Beyond resistance: antifungal heteroresistance and antifungal tolerance in fungal pathogens

Fungal infections are increasing globally, causing alarmingly high mortality and economic burden. In addition to antifungal resistance, other more subtle drug responses appear to increase the likelihood of treatment failures. These responses include heteroresistance and tolerance, terms that are more well-defined for antibacterial drugs, but are also evident in pathogenic fungi. Here, we compare these antifungal responses with similarly named antibacterial responses, and we review recent advances in how we understand the routes by which antifungal heteroresistance and tolerance emerge.

Emerging Antifungal Resistance in Fungal Pathogens

Purpose of Review

Over recent decades, the number of outbreaks caused by fungi has increased for humans, plants (including important crop species) and animals. Yet this problem is compounded by emerging antifungal drug resistance in pathogenic species. Resistance develops over time when fungi are exposed to drugs either in the patient or in the environment.

Recent Findings

Novel resistant variants of fungal pathogens that were previously susceptible are evolving (such as Aspergillus fumigatus) as well as newly emerging fungal species that are displaying antifungal resistance profiles (e.g. Candida auris and Trichophyton indotineae).

Summary

This review highlights the important topic of emerging antifungal resistance in fungal pathogens and how it evolved, as well as how this relates to a growing public health burden.

Erg251 has complex and pleiotropic effects on azole susceptibility, filamentation, and stress response phenotypes

Ergosterol is essential for fungal cell membrane integrity and growth, and numerous antifungal drugs target ergosterol. Inactivation or modification of ergosterol biosynthetic genes can lead to changes in antifungal drug susceptibility, filamentation and stress response. Here, we found that the ergosterol biosynthesis gene ERG251 is a hotspot for point mutations during adaptation to antifungal drug stress within two distinct genetic backgrounds of Candida albicans. Heterozygous point mutations led to single allele dysfunction of ERG251 and resulted in azole tolerance in both genetic backgrounds. This is the first known example of point mutations causing azole tolerance in C. albicans. Importantly, single allele dysfunction of ERG251 in combination with recurrent chromosome aneuploidies resulted in bona fide azole resistance. Homozygous deletions of ERG251 caused increased fitness in low concentrations of fluconazole and decreased fitness in rich medium, especially at low initial cell density. Dysfunction of ERG251 resulted in transcriptional upregulation of the alternate sterol biosynthesis pathway and ZRT2, a Zinc transporter. Notably, we determined that overexpression of ZRT2 is sufficient to increase azole tolerance in C. albicans. Our combined transcriptional and phenotypic analyses revealed the pleiotropic effects of ERG251 on stress responses including cell wall, osmotic and oxidative stress. Interestingly, while loss of either allele of ERG251 resulted in similar antifungal drug responses, we observed functional divergence in filamentation regulation between the two alleles of ERG251 (ERG251-A and ERG251-B) with ERG251-A exhibiting a dominant role in the SC5314 genetic background. Finally, in a murine model of systemic infection, homozygous deletion of ERG251 resulted in decreased virulence while the heterozygous deletion mutants maintain their pathogenicity. Overall, this study provides extensive genetic, transcriptional and phenotypic analysis for the effects of ERG251 on drug susceptibility, fitness, filamentation and stress responses.

Whole genome sequencing analysis demonstrates therapy-induced echinocandin resistance in Candida auris isolates

Candida auris is an emerging, multidrug-resistant yeast, causing outbreaks in healthcare facilities. Echinocandins are the antifungal drugs of choice to treat candidiasis, as they cause few side effects and resistance is rarely found. Previously, immunocompromised patients from Kuwait with C. auris colonisation or infection were treated with echinocandins, and within days to months, resistance was reported in urine isolates. To determine whether the development of echinocandin resistance was due to independent introductions of resistant strains or resulted from intra-patient resistance development, whole genome sequencing (WGS) single-nucleotide polymorphism (SNP) analysis was performed on susceptible (n = 26) and echinocandin-resistant (n = 6) isolates from seven patients. WGS SNP analysis identified three distinct clusters differing 17-127 SNPs from two patients, and the remaining isolates from five patients, respectively. Sequential isolates within patients had a maximum of 11 SNP differences over a time period of 1-10 months. The majority of isolates with reduced susceptibility displayed unique FKS1 substitutions including a novel FKS1M690V substitution, and nearly all were genetically related, ranging from only three to six SNP differences compared to susceptible isolates from the same patient. Resistant isolates from three patients shared the common FKS1S639F substitution; however, WGS analysis did not suggest a common source. These findings strongly indicate that echinocandin resistance is induced during antifungal treatment. Future studies should determine whether such echinocandin-resistant strains are capable of long-term colonisation, cause subsequent breakthrough candidiasis, have a propensity to cross-infect other patients, or remain viable for longer time periods in the hospital environment.

Acquired Triazole Resistance Alters Pathogenicity-Associated Features in Candida auris in an Isolate-Dependent Manner

Fluconazole resistance is commonly encountered in Candida auris, and the yeast frequently displays resistance to other standard drugs, which severely limits the number of effective therapeutic agents against this emerging pathogen. In this study, we aimed to investigate the effect of acquired azole resistance on the viability, stress response, and virulence of this species. Fluconazole-, posaconazole-, and voriconazole- resistant strains were generated from two susceptible C. auris clinical isolates (0381, 0387) and compared under various conditions. Several evolved strains became pan-azole-resistant, as well as echinocandin-cross-resistant. While being pan-azole-resistant, the 0381-derived posaconazole-evolved strain colonized brain tissue more efficiently than any other strain, suggesting that fitness cost is not necessarily a consequence of resistance development in C. auris. All 0387-derived evolved strains carried a loss of function mutation (R160S) in BCY1, an inhibitor of the PKA pathway. Sequencing data also revealed that posaconazole treatment can result in ERG3 mutation in C. auris. Despite using the same mechanisms to generate the evolved strains, both genotype and phenotype analysis highlighted that the development of resistance was unique for each strain. Our data suggest that C. auris triazole resistance development is a highly complex process, initiated by several pleiotropic factors.

The Microevolution of Antifungal Drug Resistance in Pathogenic Fungi

The mortality rates of invasive fungal infections remain high because of the limited number of antifungal drugs available and antifungal drug resistance, which can rapidly evolve during treatment. Mutations in key resistance genes such as ERG11 were postulated to be the predominant cause of antifungal drug resistance in the clinic. However, recent advances in whole genome sequencing have revealed that there are multiple mechanisms leading to the microevolution of resistance. In many fungal species, resistance can emerge through ERG11-independent mechanisms and through the accumulation of mutations in many genes to generate a polygenic resistance phenotype. In addition, genome sequencing has revealed that full or partial aneuploidy commonly occurs in clinical or microevolved in vitro isolates to confer antifungal resistance. This review will provide an overview of the mutations known to be selected during the adaptive microevolution of antifungal drug resistance and focus on how recent advances in genome sequencing technology have enhanced our understanding of this process.

Recurrent copy number variations in the human fungal pathogen Candida parapsilosis

Candida parapsilosis is an opportunistic fungal pathogen with increasing incidence in hospital settings worldwide; however, we lack a comprehensive understanding of the mechanisms promoting its virulence and drug resistance. Bergin et al. systematically quantify the frequency and effect of copy number variation (CNV) across 170 diverse clinical and environmental isolates of C. parapsilosis (Bergin SA, Zhao F, Ryan AP, Müller CA, Nieduszynski CA, Zhai B, Rolling T, Hohl TM, Morio F, Scully J, Wolfe KH, Butler G, 2022, mBio, https://doi.org/10.1128/mbio.01777-22). Using a combination of both short- and long-read whole genome sequencing techniques, they determine the structure and copy number of two CNVs that arose recurrently throughout the evolution of these isolates. Each CNV predominantly amplifies one coding sequence (ARR3 or RTA3); however, the amplitude and recombination breakpoints are variable across the isolates. Amplification of RTA3 correlates with drug resistance and deletion causes drug susceptibility. This study highlights the need for further research into the mechanisms and dynamics of CNV formation and the impact of these CNVs on virulence and drug resistance across diverse fungal pathogens.

Nothing makes sense in drug resistance except in the light of evolution

Our ability to fight infectious diseases is being increasingly compromised due to the emergence and spread of pathogens that become resistant to one or several drugs. This phenomenon is ubiquitous among pathogens and has parallels in cancer treatment. Given the urgency of the problem, there is a need for a paradigm shift in drug therapy toward one in which the objective to prevent the evolution of drug resistance is considered alongside the main objective of eliminating the infection or tumor. Here, I stress the importance of considering an evolutionary perspective to achieve this goal, and review recent advances in this direction, including therapies that exploit the fitness trade-offs of resistance.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization

Introduction

The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens.

Methods

TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens.

Results

A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances.

Discussion

TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.

Intra-clade Heterogeneity in Candida auris: Risk of Management

Candida auris has emerged as a significant nosocomial fungal pathogen with a high risk of pathogenicity. Since the initial detection of C. auris in 2009, it gained lots of attention with a recent alert by the Centers for Disease Control and Prevention (CDC) due to its high infectivity and drug resistance. Several studies showed the capability of C. auris to secrete lytic enzymes, germinate, and form a biofilm that eventually results in interactions with the host cells, leading to serious infections. Other studies demonstrated a decrease in susceptibility of C. auris strains to available antifungals, which may be caused by mutations within the target genes, or the drug efflux pumps. However, the contribution of C. auris heterogeneity in pathogenicity and drug resistance is not well studied. Here, we shed light on the factors contributing to the development of heterogeneity in C. auris. These include phenotypic changes, biofilm formation, mechanisms of drug resistance, host invasion, mode of transmission, and expression of virulence factors. C. auris exhibits different phenotypes, particularly aggregative, and non-aggregative forms that play an important role in fungal heterogeneity, which significantly affects drug resistance and pathogenicity. Collectively, heterogeneity in C. auris significantly contributes to ineffective treatment, which in turn affects the fungal pathogenicity and drug resistance. Therefore, understanding the underlying reasons for C. auris heterogeneity and applying effective antifungal stewardship could play a major role in controlling this pathogen.

One population, multiple lifestyles: Commensalism and pathogenesis in the human mycobiome

Candida auris and Candida albicans can result in invasive fungal diseases. And yet, these species can stably and asymptomatically colonize human skin and gastrointestinal tracts. To consider these disparate microbial lifestyles, we first review factors shown to influence the underlying microbiome. Structured by the damage response framework, we then consider the molecular mechanisms deployed by C. albicans to switch between commensal and pathogenic lifestyles. Next, we explore this framework with C. auris to highlight how host physiology, immunity, and/or antibiotic receipt are associated with progression from colonization to infection. While treatment with antibiotics increases the risk that an individual will succumb to invasive candidiasis, the underlying mechanisms remain unclear. Here, we describe several hypotheses that may explain this phenomenon. We conclude by highlighting future directions integrating genomics with immunology to advance our understanding of invasive candidiasis and human fungal disease.

Antifungal Drug Concentration Impacts the Spectrum of Adaptive Mutations in Candida albicans

Invasive fungal infections are a leading global cause of human mortality. Only three major classes of antifungal drugs are widely used, and resistance to all three classes can arise rapidly. The most widely prescribed antifungal drug, fluconazole, disseminates rapidly and reaches a wide range of concentrations throughout the body. The impact of drug concentration on the spectrum and effect of mutations acquired during adaptation is not known for any fungal pathogen, and how the specific level of a given stress influences the distribution of beneficial mutations has been poorly explored in general. We evolved 144 lineages from three genetically distinct clinical isolates of Candida albicans to four concentrations of fluconazole (0, 1, 8, and 64 μg/ml) and performed comprehensive phenotypic and genomic comparisons of ancestral and evolved populations. Adaptation to different fluconazole concentrations resulted in distinct adaptive trajectories. In general, lineages evolved to drug concentrations close to their MIC50 (the level of drug that reduces growth by 50% in the ancestor) tended to rapidly evolve an increased MIC50 and acquired distinct segmental aneuploidies and copy number variations. By contrast, lineages evolved to drug concentrations above their ancestral MIC50 tended to acquire a different suite of mutational changes and increased in drug tolerance (the ability of a subpopulation of cells to grow slowly above their MIC50). This is the first evidence that different concentrations of drug can select for different genotypic and phenotypic outcomes in vitro and may explain observed in vivo drug response variation.

Strategies to Prevent Transmission of Candida auris in Healthcare Settings

Purpose of Review

Candida auris, a recently recognized yeast pathogen, has become a major public health threat due to the problems associated with its accurate identification, intrinsic and acquired resistance to antifungal drugs, and its potential to easily contaminate the environment causing clonal outbreaks in healthcare facilities. These outbreaks are associated with high mortality rates particularly among older patients with multiple comorbidities under intensive care settings. The purpose of this review is to highlight strategies that are being adapted to prevent transmission of C. auris in healthcare settings.

Recent Findings

Colonized patients shed C. auris into their environment which contaminates surrounding equipment. It resists elimination even by robust decontamination procedures and is easily transmitted to new patients during close contact resulting in outbreaks. Efforts are being made to rapidly identify C. auris-infected/C. auris-colonized patients, to determine its susceptibility to antifungals, and to perform effective cleaning and decontamination of the environment and isolation of colonized patients to prevent further transmission.

Summary

Rapid and accurate identification of hospitalized patients infected/colonized with C. auris, rapid detection of its susceptibility patterns, and appropriate use of infection control measures can help to contain the spread of this highly pathogenic yeast in healthcare settings and prevent/control outbreaks.

Directed Evolution Detects Supernumerary Centric Chromosomes Conferring Resistance to Azoles in Candida auris

Candida auris exhibits resistance to multiple antifungal drug classes and sterilization agents, posing threats to the immunocompromised worldwide. Among the four major geographical clades, the East Asian clade 2 isolates of C. auris are mostly drug susceptible. In this study, we experimentally evolved one such drug-susceptible isolate for multiple generations in the presence of the antifungal compound fluconazole and analyzed changes in the karyotype, DNA sequence, and gene expression profiles in three evolved drug-resistant isolates. Next-generation sequencing and electrophoretic karyotyping confirm the presence of segmental aneuploidy as supernumerary chromosomes originating from centromere-inclusive chromosomal duplication events in two such cases. A 638-kb region and a 675-kb region, both of which originated from chromosome 5 and contained its centromere region, are instances of supernumerary chromosome formation identified in two evolved fluconazole-resistant isolates. Loss of the supernumerary chromosomes from the drug-resistant isolates results in a complete reversal of fluconazole susceptibility. Transcriptome analysis of the third isolate identified overexpression of drug efflux pumps as a possible non-aneuploidy-driven mechanism of drug resistance. Together, this study reveals how both aneuploidy-driven and aneuploidy-independent mechanisms may operate in parallel in an evolving population of C. auris in the presence of an antifungal drug, in spite of starting from the same strain grown under similar conditions, to attain various levels of fluconazole resistance. IMPORTANCE Fungal pathogens develop drug resistance through multiple pathways by acquiring gene mutations, increasing the copy number of genes, or altering gene expression. In this study, we attempt to understand the mechanisms of drug resistance in the recently emerged superbug, C. auris. One approach to studying this aspect is identifying various mechanisms operating in drug-resistant clinical isolates. An alternative approach is to evolve a drug-susceptible isolate in the presence of an antifungal compound and trace the changes that result in drug resistance. Here, we evolve a drug-susceptible isolate of C. auris in the laboratory in the presence of a widely used antifungal compound, fluconazole. In addition to the already known changes like overexpression of drug efflux pumps, this study identifies a novel mechanism of azole resistance by the emergence of additional chromosomes through segmental duplication of chromosomal regions, including centromeres. The centric supernumerary chromosome helps stable amplification of a set of genes with an extra copy to confer fluconazole resistance.

Drug Resistance and Evolvability in an Emerging Human Fungal Pathogen

Over the past decade, Candida auris has emerged as a highly transmissible human fungal pathogen. Because of its ability to transmit between patients in hospitals and its ability to rapidly develop drug resistance, C. auris presents unique challenges. However, at a genetic and genomic level we still understand relatively little about how drug resistance develops in this pathogen. Burrack et al. use experimental evolution and whole-genome sequencing to identify mutations correlated with fluconazole resistance in C. auris. They identify interesting genomic features, including highly plastic subtelomeric regions and whole chromosomal and segmental aneuploidies. Excitingly, they also identify the first example of a hypermutator strain in C. auris. In comparison with the model human fungal pathogen Candida albicans, C. auris is more likely to undergo mutation and less likely to undergo copy number variation in response to drug selection, which may be linked to differences in base ploidy level.