Single-cell detection of copy number changes reveals dynamic mechanisms of adaptation to antifungals in Candida albicans.

Genomic copy number changes are associated with antifungal drug resistance and virulence across diverse fungal pathogens, but the rate and dynamics of these genomic changes in the presence of antifungal drugs are unknown. Here we optimized a dual-fluorescent reporter system in the diploid pathogen Candida albicans to quantify haplotype-specific copy number variation (CNV) and loss of heterozygosity (LOH) at the single-cell level with flow cytometry. We followed the frequency and dynamics of CNV and LOH at two distinct genomic locations in the presence and absence of antifungal drugs in vitro and in a murine model of candidiasis. Copy number changes were rapid and dynamic during adaptation to fluconazole and frequently involved competing subpopulations with distinct genotypes. This study provides quantitative evidence for the rapid speed at which diverse genotypes arise and undergo dynamic population-level fluctuations during adaptation to antifungal drugs in vitro and in vivo.

Step-wise evolution of azole resistance through copy number variation followed by KSR1 loss of heterozygosity in Candida albicans.

Antimicrobial drug resistance poses a global health threat, requiring a deeper understanding of the evolutionary processes that lead to its emergence in pathogens. Complex evolutionary dynamics involve multiple mutations that can result in cooperative or competitive (clonal interference) effects. Candida albicans, a major fungal pathogen, displays high rates of copy number variation (CNV) and loss of heterozygosity (LOH). CNV and LOH events involve large numbers of genes and could synergize during evolutionary adaptation. Understanding the contributions of CNV and LOH to antifungal drug adaptation is challenging, especially in the context of whole-population genome sequencing. Here, we document the sequential evolution of fluconazole tolerance and then resistance in a C. albicans isolate involving an initial CNV on chromosome 4, followed by an LOH on chromosome R that involves KSR1. Similar LOH events involving KSR1, which encodes a reductase in the sphingolipid biosynthesis pathway, were also detected in independently evolved fluconazole resistant isolates. We dissect the specific KSR1 codons that affect fluconazole resistance and tolerance. The combination of the chromosome 4 CNV and KSR1 LOH results in a >500-fold decrease in azole susceptibility relative to the progenitor, illustrating a compelling example of rapid, yet step-wise, interplay between CNV and LOH in drug resistance evolution.

Erg251 has complex and pleiotropic effects on sterol composition, 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. Homozygous deletions of ERG251 resulted in accumulation of ergosterol intermediates consistent with the fitness defect in rich medium. Dysfunction of ERG251, together with FLC exposure, resulted in decreased accumulation of the toxic sterol (14-ɑ-methylergosta-8,24(28)-dien-3ß,6α-diol) and increased accumulation of non-toxic alternative sterols. The altered sterol composition of the ERG251 mutants had pleiotropic effects on transcription, filamentation, and stress responses including cell membrane, osmotic and oxidative stress. Interestingly, while dysfunction of ERG251 resulted in azole tolerance, it also led to transcriptional upregulation of ZRT2, a membrane-bound Zinc transporter, in the presence of FLC, and overexpression of ZRT2 is sufficient to increase azole tolerance in wild-type C. albicans. 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 demonstrates that single allele dysfunction of ERG251 is a recurrent and effective mechanism of acquired azole tolerance. We propose that altered sterol composition resulting from ERG251 dysfunction mediates azole tolerance as well as pleiotropic effects on stress response, filamentation and virulence.

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.

Aneuploidy Can Be an Evolutionary Diversion on the Path to Adaptation.

Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.

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.

Rapid Evolution of Multidrug Resistance in a Candida lusitaniae Infection during Micafungin Monotherapy.

Candida (Clavispora) lusitaniae is a rare, emerging non-albicans Candida species that can cause life-threatening invasive infections, spread within hospital settings, and rapidly acquire antifungal drug resistance, including multidrug resistance. The frequency and spectrum of mutations causing antifungal drug resistance in C. lusitaniae are poorly understood. Analyses of serial clinical isolates of any Candida species are uncommon and often analyze a limited number of samples collected over months of antifungal therapy with multiple drug classes, limiting the ability to understand relationships between drug classes and specific mutations. Here, we performed comparative genomic and phenotypic analysis of 20 serial C. lusitaniae bloodstream isolates collected daily from an individual patient treated with micafungin monotherapy during a single 11-day hospital admission. We identified isolates with decreased micafungin susceptibility 4 days after initiation of antifungal therapy and a single isolate with increased cross-resistance to micafungin and fluconazole, despite no history of azole therapy in this patient. Only 14 unique single nucleotide polymorphisms (SNPs) were identified between all 20 samples, including three different FKS1 alleles among isolates with decreased micafungin susceptibility and an ERG3 missense mutation found only in the isolate with increased cross-resistance to both micafungin and fluconazole. This is the first clinical evidence of an ERG3 mutation in C. lusitaniae that occurred during echinocandin monotherapy and is associated with cross-resistance to multiple drug classes. Overall, the evolution of multidrug resistance in C. lusitaniae is rapid and can emerge during treatment with only first-line antifungal therapy.

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.

Copy Number Variation and Allele Ratio Analysis in Candida albicans Using Whole Genome Sequencing Data.

Whole genome sequencing of human fungal pathogens has revolutionized the speed and accuracy in which sequence variants that cause antifungal resistance can be identified. Genome rearrangements resulting in copy number variation (CNV) are a significant source of acquired antifungal drug resistance across diverse fungal species. Some CNVs are transient in nature, while other CNVs are stable and well tolerated even in the absence of antifungal drugs. By visualizing whole genome sequencing read depth as a function of genomic location, CNVs and CNV breakpoints (genomic positions where the copy number changes occur relative to the rest of the genome) are rapidly identified. A similar analysis can be used to visualize allele ratio changes that occur across the genomes of heterozygous fungal species, both in the presence and absence of CNVs. This protocol walks through the bioinformatic analysis of CNVs and allele ratios utilizing free, open-source visualization tools. We provide code to use with an example dataset (matched antifungal drug-sensitive and drug-resistant Candida albicans isolates) and notes on how to expand this protocol to other fungal genomes.

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.

Stress combined with loss of the Candida albicans SUMO protease Ulp2 triggers selection of aneuploidy via a two-step process.

A delicate balance between genome stability and instability ensures genome integrity while generating genetic diversity, a critical step for evolution. Indeed, while excessive genome instability is harmful, moderated genome instability can drive adaptation to novel environments by maximising genetic variation. Candida albicans, a human fungal pathogen that colonises different parts of the human body, adapts rapidly and frequently to different hostile host microenvironments. In this organism, the ability to generate large-scale genomic variation is a key adaptative mechanism triggering dangerous infections even in the presence of antifungal drugs. Understanding how fitter novel karyotypes are selected is key to determining how C. albicans and other microbial pathogens establish infections. Here, we identified the SUMO protease Ulp2 as a regulator of C. albicans genome integrity through genetic screening. Deletion of ULP2 leads to increased genome instability, enhanced genome variation and reduced fitness in the absence of additional stress. The combined stress caused by the lack of ULP2 and antifungal drug treatment leads to the selection of adaptive segmental aneuploidies that partially rescue the fitness defects of ulp2Δ/Δ cells. Short and long-read genomic sequencing demonstrates that these novel genotypes are selected via a two-step process leading to the formation of novel chromosomal fragments with breakpoints at microhomology regions and DNA repeats.

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

Antifungal drug resistance and tolerance pose a serious threat to global public health. In the human fungal pathogen, Candida auris, resistance to triazole, polyene, and echinocandin antifungals is rising, resulting in multidrug resistant isolates. Here, we use genome analysis and in vitro evolution of 17 new clinical isolates of C. auris from clades I and IV to determine how quickly resistance mutations arise, the stability of resistance in the absence of drug, and the impact of genetic background on evolutionary trajectories. We evolved each isolate in the absence of drug as well as in low and high concentrations of fluconazole. In just three passages, we observed genomic and phenotypic changes including karyotype alterations, aneuploidy, acquisition of point mutations, and increases in MIC values within the populations. Fluconazole resistance was stable in the absence of drug, indicating little to no fitness cost associated with resistance. Importantly, two isolates substantially increased resistance to ≥256 µg/mL fluconazole. Multiple evolutionary pathways and mutations associated with increased fluconazole resistance occurred simultaneously within the same population. Strikingly, the subtelomeric regions of C. auris were highly dynamic as deletion of multiple genes near the subtelomeres occurred during the three passages in several populations. Finally, we discovered a mutator phenotype in a clinical isolate of C. auris. This isolate had elevated mutation rates compared to other isolates and acquired substantial resistance during evolution in vitro and in vivo supporting that the genetic background of clinical isolates can have a significant effect on evolutionary potential. IMPORTANCE Drug resistant Candida auris infections are recognized by the CDC as an urgent threat. Here, we obtained and characterized a set of clinical isolates of C. auris including multiple isolates from the same patient. To understand how drug resistance arises, we evolved these isolates and found that resistance to fluconazole, the most commonly prescribed antifungal, can occur rapidly and that there are multiple pathways to resistance. During our experiment, resistance was gained, but it was not lost, even in the absence of drug. We also found that some C. auris isolates have higher mutation rates than others and are primed to acquire antifungal resistance mutations. Furthermore, we found that multidrug resistance can evolve within a single patient. Overall, our results highlight the high stability and high rates of acquisition of antifungal resistance of C. auris that allow evolution of pan-resistant, transmissible isolates in the clinic.

A small molecule produced by Lactobacillus species blocks Candida albicans filamentation by inhibiting a DYRK1-family kinase.

The fungus Candida albicans is an opportunistic pathogen that can exploit imbalances in microbiome composition to invade its human host, causing pathologies ranging from vaginal candidiasis to fungal sepsis. Bacteria of the genus Lactobacillus are colonizers of human mucosa and can produce compounds with bioactivity against C. albicans. Here, we show that some Lactobacillus species produce a small molecule under laboratory conditions that blocks the C. albicans yeast-to-filament transition, an important virulence trait. It remains unexplored whether the compound is produced in the context of the human host. Bioassay-guided fractionation of Lactobacillus-conditioned medium linked this activity to 1-acetyl-ß-carboline (1-ABC). We use genetic approaches to show that filamentation inhibition by 1-ABC requires Yak1, a DYRK1-family kinase. Additional biochemical characterization of structurally related 1-ethoxycarbonyl-ß-carboline confirms that it inhibits Yak1 and blocks C. albicans biofilm formation. Thus, our findings reveal Lactobacillus-produced 1-ABC can prevent the yeast-to-filament transition in C. albicans through inhibition of Yak1.

The fitness costs and benefits of trisomy of each Candida albicans chromosome.

Candida albicans is a prevalent human fungal pathogen. Rapid genomic change, due to aneuploidy, is a common mechanism that facilitates survival from multiple types of stresses including the few classes of available antifungal drugs. The stress survival of aneuploids occurs despite the fitness costs attributed to most aneuploids growing under idealized lab conditions. Systematic study of the aneuploid state in C. albicans has been hindered by the lack of a comprehensive collection of aneuploid strains. Here, we describe a collection of diploid C. albicans aneuploid strains, each carrying one extra copy of each chromosome, all from the same genetic background. We tested the fitness of this collection under several physiological conditions including shifts in pH, low glucose, oxidative stress, temperature, high osmolarity, membrane stress, and cell wall stress. We found that most aneuploids, under most conditions, were less fit than their euploid parent, yet there were specific conditions under which specific aneuploid isolates provided a fitness benefit relative to the euploid parent strain. Importantly, this fitness benefit was attributable to the change in the copy number of specific chromosomes. Thus, C. albicans can tolerate aneuploidy of each chromosome and some aneuploids confer improved growth under conditions that the yeast encounters in its host niches.

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs.

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

Genome plasticity in Candida albicans is driven by long repeat sequences.

Genome rearrangements resulting in copy number variation (CNV) and loss of heterozygosity (LOH) are frequently observed during the somatic evolution of cancer and promote rapid adaptation of fungi to novel environments. In the human fungal pathogen Candida albicans, CNV and LOH confer increased virulence and antifungal drug resistance, yet the mechanisms driving these rearrangements are not completely understood. Here, we unveil an extensive array of long repeat sequences (65-6499 bp) that are associated with CNV, LOH, and chromosomal inversions. Many of these long repeat sequences are uncharacterized and encompass one or more coding sequences that are actively transcribed. Repeats associated with genome rearrangements are predominantly inverted and separated by up to ~1.6 Mb, an extraordinary distance for homology-based DNA repair/recombination in yeast. These repeat sequences are a significant source of genome plasticity across diverse strain backgrounds including clinical, environmental, and experimentally evolved isolates, and represent previously uncharacterized variation in the reference genome.

Selection of Candida albicans trisomy during oropharyngeal infection results in a commensal-like phenotype.

When the fungus Candida albicans proliferates in the oropharyngeal cavity during experimental oropharyngeal candidiasis (OPC), it undergoes large-scale genome changes at a much higher frequency than when it grows in vitro. Previously, we identified a specific whole chromosome amplification, trisomy of Chr6 (Chr6x3), that was highly overrepresented among strains recovered from the tongues of mice with OPC. To determine the functional significance of this trisomy, we assessed the virulence of two Chr6 trisomic strains and a Chr5 trisomic strain in the mouse model of OPC. We also analyzed the expression of virulence-associated traits in vitro. All three trisomic strains exhibited characteristics of a commensal during OPC in mice. They achieved the same oral fungal burden as the diploid progenitor strain but caused significantly less weight loss and elicited a significantly lower inflammatory host response. In vitro, all three trisomic strains had reduced capacity to adhere to and invade oral epithelial cells and increased susceptibility to neutrophil killing. Whole genome sequencing of pre- and post-infection isolates found that the trisomies were usually maintained. Most post-infection isolates also contained de novo point mutations, but these were not conserved. While in vitro growth assays did not reveal phenotypes specific to de novo point mutations, they did reveal novel phenotypes specific to each lineage. These data reveal that during OPC, clones that are trisomic for Chr5 or Chr6 are selected and they facilitate a commensal-like phenotype.

Draft Genome Assemblies of Clinical Isolates of Klebsiella pneumoniae V9011662 and Enterobacter hormaechei Entb306.

Enterobacter hormaechei and Klebsiella pneumoniae are pathogenic Enterobacteriaceae that have been associated with the spread of antibiotic resistance. Here, we report draft genome assemblies of an Enterobacter hormaechei clinical isolate and a multidrug-resistant clinical isolate of Klebsiella pneumoniae.