Global analysis of mutations driving microevolution of a heterozygous diploid fungal pathogen

The dynamics of loss of heterozygosity events in genomes

Genomic instability is a hallmark of tumorigenesis, yet it also plays an essential role in evolution. Large-scale population genomics studies have highlighted the importance of loss of heterozygosity (LOH) events, which have long been overlooked in the context of genetic diversity and instability. Among various types of genomic mutations, LOH events are the most common and affect a larger portion of the genome. They typically arise from recombination-mediated repair of double-strand breaks (DSBs) or from lesions that are processed into DSBs. LOH events are critical drivers of genetic diversity, enabling rapid phenotypic variation and contributing to tumorigenesis. Understanding the accumulation of LOH, along with its underlying mechanisms, distribution, and phenotypic consequences, is therefore crucial. In this review, we explore the spectrum of LOH events, their mechanisms, and their impact on fitness and phenotype, drawing insights from Saccharomyces cerevisiae to cancer. We also emphasize the role of LOH in genomic instability, disease, and genome evolution.

Elevated mutation rates in multi-azole resistant Aspergillus fumigatus drive rapid evolution of antifungal resistance

The environmental use of azole fungicides has led to selective sweeps across multiple loci in the Aspergillus fumigatus genome causing the rapid global expansion of a genetically distinct cluster of resistant genotypes. Isolates within this cluster are also more likely to be resistant to agricultural antifungals with unrelated modes of action. Here we show that this cluster is not only multi-azole resistant but has increased propensity to develop resistance to next generation antifungals because of variants in the DNA mismatch repair system. A variant in msh6-G233A is found almost exclusively within azole resistant isolates harbouring the canonical cyp51A azole resistance allelic variant TR34/L98H. Naturally occurring isolates with this msh6 variant display up to 5-times higher rate of mutation, leading to an increased likelihood of evolving resistance to other antifungals. Furthermore, unlike hypermutator strains, the G233A variant conveys no measurable fitness cost and has become globally distributed. Our findings further suggest that resistance to next-generation antifungals is more likely to emerge within organisms that are already multi-azole resistant due to close linkage between TR34/L98H and msh6-G233A, posing a major problem due to the prospect of dual use of novel antifungals in clinical and agricultural settings.

Candida albicans strains adapted to the mouse gut are resistant to bile salts via a Flo8-dependent mechanism

Candidaalbicans normally colonizes the human gastrointestinal tract as a commensal. Studying fungal factors involved in colonizing the mammalian gastrointestinal tract requires mouse models with altered microbiota. We have obtained strains of C.albicans through microevolution in the mouse gut for a prolonged period (one year) that display a substantial increase in fitness in this niche. These strains show resistance to bile salts, an increase in their adhesion to the intestinal mucosa, and are unable to filament in response to serum. Genetic analysis revealed some alterations, mainly a triploidy of chr7, a whole chr6 homozygosis, and an SNP in the FLO8 gene (located in the chr6), resulting in a truncated protein version. A wild type FLO8 gene complemented filamentation and bile salt sensitivity but showed an intermediate fitness phenotype in colonization. Alterations in bile salt sensitivity were also evident in bmt mutants, defective in β-mannosylation, and transcriptional targets of Flo8, suggesting a link between the fungal cell wall and mammalian gut colonization via the Flo8 transcriptional regulator.

Blood cultures contain populations of genetically diverse Candida albicans strains that may differ in echinocandin tolerance and virulence

It is unknown whether within-patient Candida albicans diversity is common during bloodstream infections (BSIs). We determined whole genome sequences of 10 C. albicans strains from blood cultures (BCs) in each of 4 patients. BCs in 3 patients contained mixed populations of strains that differed by large-scale genetic variants, including chromosome (Chr) 5 or 7 aneuploidy (n=2) and Chr1 loss of heterozygosity (n=1). Chr7 trisomy (Tri7) strains from patient MN were attenuated for hyphal and biofilm formation in vitro compared to euploid strains, due at least in part to NRG1 over-expression. Nevertheless, representative Tri7 strain M1 underwent filamentation during disseminated candidiasis (DC) in mice. M1 was more fit than euploid strain M2 during DC and mouse gastrointestinal colonization, and in blood ex vivo. M1 and M2 exhibited identical echinocandin minimum inhibitory concentrations, but M2 was more tolerant to micafungin in vitro. Furthermore, M2 was more competitive with M1 in mouse kidneys following micafungin treatment than it was in absence of micafungin. Tri7 strains represented 74% of patient MN's baseline BC population, but after 1d and 3d of echinocandin treatment, euploid strains were 93% and 98% of the BC population, respectively. Findings suggest that echinocandin tolerant, euploid strains were a subpopulation to more virulent Tri7 strains at baseline and then were selected upon echinocandin exposure. In conclusion, BCs in at least some patients are comprised of diverse C. albicans populations not recognized by the clinical lab, rather than single strains. Clinical relevance of C. albicans diversity and antifungal tolerance merits further investigation.

The pathobiology of human fungal infections

Human fungal infections are a historically neglected area of disease research, yet they cause more than 1.5 million deaths every year. Our understanding of the pathophysiology of these infections has increased considerably over the past decade, through major insights into both the host and pathogen factors that contribute to the phenotype and severity of these diseases. Recent studies are revealing multiple mechanisms by which fungi modify and manipulate the host, escape immune surveillance and generate complex comorbidities. Although the emergence of fungal strains that are less susceptible to antifungal drugs or that rapidly evolve drug resistance is posing new threats, greater understanding of immune mechanisms and host susceptibility factors is beginning to offer novel immunotherapeutic options for the future. In this Review, we provide a broad and comprehensive overview of the pathobiology of human fungal infections, focusing specifically on pathogens that can cause invasive life-threatening infections, highlighting recent discoveries from the pathogen, host and clinical perspectives. We conclude by discussing key future challenges including antifungal drug resistance, the emergence of new pathogens and new developments in modern medicine that are promoting susceptibility to infection.

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.

Advances in the design and delivery of RNA vaccines for infectious diseases

RNA medicines represent a paradigm shift in treatment and prevention of critical diseases of global significance, e.g., infectious diseases. The highly successful messenger RNA (mRNA) vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were developed at record speed during the coronavirus disease 2019 pandemic. A consequence of this is exceptionally shortened vaccine development times, which in combination with adaptability makes the RNA vaccine technology highly attractive against infectious diseases and for pandemic preparedness. Here, we review state of the art in the design and delivery of RNA vaccines for infectious diseases based on different RNA modalities, including linear mRNA, self-amplifying RNA, trans-amplifying RNA, and circular RNA. We provide an overview of the clinical pipeline of RNA vaccines for infectious diseases, and present analytical procedures, which are paramount for characterizing quality attributes and guaranteeing their quality, and we discuss future perspectives for using RNA vaccines to combat pathogens beyond SARS-CoV-2.

Controlling Candida: immune regulation of commensal fungi in the gut

The intestinal microbiome harbors fungi that pose a significant risk to human health as opportunistic pathogens and drivers of inflammation. Inflammatory and autoimmune diseases are associated with dysbiotic fungal communities and the expansion of potentially pathogenic fungi. The gut is also the main reservoir for disseminated fungal infections. Immune interactions are critical for preventing commensal fungi from becoming pathogenic. Significant strides have been made in defining innate and adaptive immune pathways that regulate intestinal fungi, and these discoveries have coincided with advancements in our understanding of the fungal molecular pathways and effectors involved in both commensal colonization and pathogenesis within the gut. In this review, we will discuss immune interactions important for regulating commensal fungi, with a focus on how specific cell types and effectors interact with fungi to limit their colonization or pathogenic potential. This will include how innate and adaptive immune pathways target fungi and orchestrate antifungal immune responses, in addition to how secreted immune effectors, such as mucus and antimicrobial peptides, regulate fungal colonization and inhibit pathogenic potential. These immune interactions will be framed around our current understanding of the fungal effectors and pathways regulating colonization and pathogenesis within this niche. Finally, we highlight important unexplored mechanisms by which the immune system regulates commensal fungi in the gut.

Evolution and strain diversity advance exploration of Candida albicans biology

Fungi were some of the earliest organismal systems used to explore mutational processes and its phenotypic consequences on members of a species. Yeasts that cause significant human disease were quickly incorporated into these investigations to define the genetic and phenotypic drivers of virulence. Among Candida species, Candida albicans has emerged as a model for studying genomic processes of evolution because of its clinical relevance, relatively small genome, and ability to tolerate complex chromosomal changes. Here, we describe major recent findings that used evolution of strains from defined genetic backgrounds to delineate mutational and adaptative processes and include how nascent exploration into naturally occurring variation is contributing to these conceptual frameworks. Ultimately, efforts to discern adaptive mechanisms used by C. albicans will continue to divulge new biology and can better inform treatment regimens for the increasing prevalence of fungal disease.

Assessment of Genetic Diversity and the Population Structure of Species from the Fusarium fujikuroi Species Complex Causing Fusarium Stalk Rot of Maize

Fusarium stalk rot (FSR), caused by the Fusarium species complex, is an economic threat to maize cultivation all over the world. We investigated the population structure and genetic diversity of Fusarium species obtained from five major maize-growing regions in India. The Tef-1α locus was used for phylogenetic analysis of geographically distinct isolates of Fusarium verticillioides, F. andiyazi, F. proliferatum, F. nygamai, and F. acutatum causing FSR. A phylogenetic tree showed monophyletic, polyphyletic, and paraphyletic groupings reflecting the complex evolutionary history and genetic diversity within the genus. Monophyletic groupings depicting strong bootstrap support were shown to have a single common ancestor and genetic coherence with limited genetic divergence among sequences. Polyphyletic groupings also presented significant genetic differentiation within the F. verticillioides sequences from diverse ecological niches. Nucleotide diversity of moderate level 0.02471 reflected genetic variations within populations that were attributed to factors such as mutation, genetic drift, or varying selection pressures. The Fst value of 0.98205 is particularly indicative of high genetic differentiation, implying that most of the genetic variance is due to differences between populations rather than within them. F. verticillioides, with 57 sequences, showed low genetic diversity with three segregating sites and a low haplotype diversity of 0.19486, suggesting the founder effect, where a reduced population expands from a limited genetic pool. The total data estimates across all populations for haplotype analysis showed 72 sequences, 44 segregating sites, and 9 haplotypes with a haplotype diversity of 0.48513. The evolutionary dynamics showed genetic differentiations among Fusarium species causing FSR. AMOVA indicated high within-population variations, depicting a substantial genetic diversity within individual populations. The results offer a comprehensive framework for discussing the implications of genetic diversity in pathogen management and the evolutionary dynamics of the Fusarium species causing FSR in maize in the Indian subcontinent.

Approaches for identifying and measuring heteroresistance in azole-susceptible Candida isolates

Heteroresistance to antifungal agents poses a significant challenge in the treatment of fungal infections. Currently, the absence of established methods for detecting and measuring heteroresistance impedes progress in understanding this phenomenon in fungal pathogens. In response to this gap, we present a comprehensive set of new and optimized methods designed to detect and quantify azole heteroresistance in Candida albicans. Here, we define two primary assays for measuring heteroresistance: population analysis profiling, based on growth on solid medium, and single-cell assays, based on growth in liquid culture. We observe good correlations between the measurements obtained with liquid and solid assays, validating their utility for studying azole heteroresistance. We also highlight that disk diffusion assays could serve as an additional tool for the rapid detection of heteroresistance. These methods collectively provide a versatile toolkit for researchers seeking to assess heteroresistance in C. albicans. They also serve as a critical step forward in the characterization of antifungal heteroresistance, providing a framework for investigating this phenomenon in diverse fungal species and in the context of other antifungal agents. Ultimately, these advancements will enhance our ability to effectively measure antifungal drug responses and combat fungal infections.IMPORTANCEHeteroresistance involves varying antimicrobial susceptibility within a clonal population. This phenomenon allows the survival of rare resistant subpopulations during drug treatment, significantly complicating the effective management of infections. However, the absence of established detection methods hampers progress in understanding this phenomenon in human fungal pathogens. We propose a comprehensive toolkit to address this gap in the yeast Candida albicans, encompassing population analysis profiling, single-cell assays, and disk diffusion assays. By providing robust and correlated measurements through both solid and liquid assays, this work will provide a framework for broader applications across clinically relevant Candida species. These methods will enhance our ability to understand this phenomenon and the failure of antifungal therapy.

Rapid evolution of an adaptive multicellular morphology of Candida auris during systemic infection

Candida auris has become a serious threat to public health. The mechanisms of how this fungal pathogen adapts to the mammalian host are poorly understood. Here we report the rapid evolution of an adaptive C. auris multicellular aggregative morphology in the murine host during systemic infection. C. auris aggregative cells accumulate in the brain and exhibit obvious advantages over the single-celled yeast-form cells during systemic infection. Genetic mutations, specifically de novo point mutations in genes associated with cell division or budding processes, underlie the rapid evolution of this aggregative phenotype. Most mutated C. auris genes are associated with the regulation of cell wall integrity, cytokinesis, cytoskeletal properties, and cellular polarization. Moreover, the multicellular aggregates are notably more recalcitrant to the host antimicrobial peptides LL-37 and PACAP relative to the single-celled yeast-form cells. Overall, to survive in the host, C. auris can rapidly evolve a multicellular aggregative morphology via genetic mutations.

Adaptation to an amoeba host drives selection of virulence-associated traits and genetic variation in saprotrophic Candida albicans

Amoebae are micropredators that play an important role in controlling fungal populations in ecosystems. However, the interaction between fungi and their amoebic predators suggests that the pressure from predatory selection can significantly influence the development of fungal virulence and evolutionary processes. Thus, the purpose of this study was to investigate the adaptation of saprotrophic Candida albicans strains during their interactions with Acanthamoeba castellanii. We conducted a comprehensive analysis of survival after co-culture by colony counting of the yeast cells and examining yeast cell phenotypic and genetic characteristics. Our results indicated that exposure to amoebae enhanced the survival capacity of environmental C. albicans and induced visible morphological alterations in C. albicans, particularly by an increase in filamentation. These observed phenotypic changes were closely related to concurrent genetic variations. Notably, mutations in genes encoding transcriptional repressors (TUP1 and SSN6), recognized for their negative regulation of filamentous growth, were exclusively identified in amoeba-passaged isolates, and absent in unexposed isolates. Furthermore, these adaptations increased the exposed isolates' fitness against various stressors, simultaneously enhancing virulence factors and demonstrating an increased ability to invade A549 lung human epithelial cells. These observations indicate that the sustained survival of C. albicans under ongoing amoebic predation involved a key role of mutation events in microevolution to modulate the ability of these isolates to change phenotype and increase their virulence factors, demonstrating an enhanced potential to survive in diverse environmental niches.

Candida albicans' inorganic phosphate transport and evolutionary adaptation to phosphate scarcity

Phosphorus is essential in all cells' structural, metabolic and regulatory functions. For fungal cells that import inorganic phosphate (Pi) up a steep concentration gradient, surface Pi transporters are critical capacitators of growth. Fungi must deploy Pi transporters that enable optimal Pi uptake in pH and Pi concentration ranges prevalent in their environments. Single, triple and quadruple mutants were used to characterize the four Pi transporters we identified for the human fungal pathogen Candida albicans, which must adapt to alkaline conditions during invasion of the host bloodstream and deep organs. A high-affinity Pi transporter, Pho84, was most efficient across the widest pH range while another, Pho89, showed high-affinity characteristics only within one pH unit of neutral. Two low-affinity Pi transporters, Pho87 and Fgr2, were active only in acidic conditions. Only Pho84 among the Pi transporters was clearly required in previously identified Pi-related functions including Target of Rapamycin Complex 1 signaling and hyphal growth. We used in vitro evolution and whole genome sequencing as an unbiased forward genetic approach to probe adaptation to prolonged Pi scarcity of two quadruple mutant lineages lacking all 4 Pi transporters. Lineage-specific genomic changes corresponded to divergent success of the two lineages in fitness recovery during Pi limitation. In this process, initial, large-scale genomic alterations like aneuploidies and loss of heterozygosity were eventually lost as populations presumably gained small-scale mutations. Severity of some phenotypes linked to Pi starvation, like cell wall stress hypersensitivity, decreased in parallel to evolving populations' fitness recovery in Pi scarcity, while that of others like membrane stress responses diverged from these fitness phenotypes. C. albicans therefore has diverse options to reconfigure Pi management during prolonged scarcity. Since Pi homeostasis differs substantially between fungi and humans, adaptive processes to Pi deprivation may harbor small-molecule targets that impact fungal growth and virulence.

Candida albicans exhibits heterogeneous and adaptive cytoprotective responses to antifungal compounds

Candida albicans, an opportunistic human pathogen, poses a significant threat to human health and is associated with significant socio-economic burden. Current antifungal treatments fail, at least in part, because C. albicans can initiate a strong drug tolerance response that allows some cells to grow at drug concentrations above their minimal inhibitory concentration. To better characterize this cytoprotective tolerance program at the molecular single-cell level, we used a nanoliter droplet-based transcriptomics platform to profile thousands of individual fungal cells and establish their subpopulation characteristics in the absence and presence of antifungal drugs. Profiles of untreated cells exhibit heterogeneous expression that correlates with cell cycle stage with distinct metabolic and stress responses. At 2 days post-fluconazole exposure (a time when tolerance is measurable), surviving cells bifurcate into two major subpopulations: one characterized by the upregulation of genes encoding ribosomal proteins, rRNA processing machinery, and mitochondrial cellular respiration capacity, termed the Ribo-dominant (Rd) state; and the other enriched for genes encoding stress responses and related processes, termed the Stress-dominant (Sd) state. This bifurcation persists at 3 and 6 days post-treatment. We provide evidence that the ribosome assembly stress response (RASTR) is activated in these subpopulations and may facilitate cell survival.

Genotypic diversity and unrecognized antifungal resistance among populations of Candida glabrata from positive blood cultures

The longstanding model is that most bloodstream infections (BSIs) are caused by a single organism. We perform whole genome sequencing of five-to-ten strains from blood culture (BC) bottles in each of ten patients with Candida glabrata BSI. We demonstrate that BCs contain mixed populations of clonal but genetically diverse strains. Genetically distinct strains from two patients exhibit phenotypes that are potentially important during BSIs, including differences in susceptibility to antifungal agents and phagocytosis. In both patients, the clinical microbiology lab recovered a fluconazole-susceptible index strain, but we identify mixed fluconazole-susceptible and -resistant populations. Diversity in drug susceptibility is likely clinically relevant, as fluconazole-resistant strains were subsequently recovered by the clinical laboratory during persistent or relapsing infections. In one patient, unrecognized respiration-deficient small colony variants are fluconazole-resistant and significantly attenuated for virulence during murine candidiasis. Our data suggest a population-based model of C. glabrata genotypic and phenotypic diversity during BSIs.

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.

The use of whole-genome sequencing and development of bioinformatics to monitor overlapping outbreaks of Candida auris in southern Nevada

A Candida auris outbreak has been ongoing in Southern Nevada since August 2021. In this manuscript we describe the sequencing of over 200 C. auris isolates from patients at several facilities. Genetically distinct subgroups of C. auris were detected from Clade I (3 distinct lineages) and III (1 lineage). Open-source bioinformatic tools were developed and implemented to aid in the epidemiological investigation. The work herein compares three methods for C. auris whole genome analysis: Nullarbor, MycoSNP and a new pipeline TheiaEuk. We also describe a novel analysis method focused on elucidating phylogenetic linkages between isolates within an ongoing outbreak. Moreover, this study places the ongoing outbreaks in a global context utilizing existing sequences provided worldwide. Lastly, we describe how the generated results were communicated to the epidemiologists and infection control to generate public health interventions.

Loss of Heterozygosity and Its Importance in Evolution

Loss of heterozygosity (LOH) is a mitotic recombination event that converts heterozygous loci to homozygous loci. This mutation event is widespread in organisms that have asexual reproduction like budding yeasts, and is also an important and frequent mutation event in tumorigenesis. Mutation accumulation studies have demonstrated that LOH occurs at a rate higher than the point mutation rate, and can impact large portions of the genome. Laboratory evolution experiments of heterozygous yeasts have revealed that LOH often unmasks beneficial recessive alleles that can confer large fitness advantages. Here, I highlight advances in understanding dominance, fitness, and phenotypes in laboratory evolved heterozygous yeast strains. I discuss best practices for detecting LOH in intraspecific and interspecific evolved clones and populations. Utilizing heterozygous strain backgrounds in laboratory evolution experiments offers an opportunity to advance our understanding of this important mutation type in shaping adaptation and genome evolution in wild, domesticated, and clinical populations.

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.

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.

Genotypic diversity and unrecognized antifungal resistance among populations of Candida glabrata from positive blood cultures

The longstanding paradigm is that most bloodstream infections (BSIs) are caused by a single organism. We performed whole genome sequencing of five-to-ten strains from blood culture (BC) bottles in each of ten patients with Candida glabrata BSI. We demonstrated that BCs contained mixed populations of clonal but genetically diverse strains. Genetically distinct strains from two patients exhibited phenotypes that were potentially important during BSIs, including differences in susceptibility to antifungal agents and phagocytosis. In both patients, the clinical microbiology lab recovered a fluconazole-susceptible index strain, but we identified mixed fluconazole-susceptible and â€"resistant populations. Diversity in drug susceptibility was likely clinically relevant, as fluconazole-resistant strains were subsequently recovered by the clinical laboratory during persistent or relapsing infections. In one patient, unrecognized respiration-deficient small colony variants were fluconazole-resistant and significantly attenuated for virulence during murine candidiasis. Our data suggest a new population-based paradigm of C. glabrata genotypic and phenotypic diversity during BSIs.

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.

Aneuploidy and gene dosage regulate filamentation and host colonization by Candida albicans

Aneuploidy is a frequent occurrence in fungal species where it can alter gene expression and promote adaptation to a variety of environmental cues. Multiple forms of aneuploidy have been observed in the opportunistic fungal pathogen Candida albicans, which is a common component of the human gut mycobiome but can escape this niche and cause life-threatening systemic disease. Using a barcode sequencing (Bar-seq) approach, we evaluated a set of diploid C. albicans strains and found that a strain carrying a third copy of chromosome (Chr) 7 was associated with increased fitness during both gastrointestinal (GI) colonization and systemic infection. Our analysis revealed that the presence of a Chr 7 trisomy resulted in decreased filamentation, both in vitro and during GI colonization, relative to isogenic euploid controls. A target gene approach demonstrated that NRG1, encoding a negative regulator of filamentation located on Chr 7, contributes to increased fitness of the aneuploid strain due to inhibition of filamentation in a gene dosage-dependent fashion. Together, these experiments establish how aneuploidy enables the reversible adaptation of C. albicans to its host via gene dosage-dependent regulation of morphology.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

Fungal resilience and host-pathogen interactions: Future perspectives and opportunities

We are constantly exposed to the threat of fungal infection. The outcome-clearance, commensalism or infection-depends largely on the ability of our innate immune defences to clear infecting fungal cells versus the success of the fungus in mounting compensatory adaptive responses. As each seeks to gain advantage during these skirmishes, the interactions between host and fungal pathogen are complex and dynamic. Nevertheless, simply compromising the physiological robustness of fungal pathogens reduces their ability to evade antifungal immunity, their virulence, and their tolerance against antifungal therapy. In this article I argue that this physiological robustness is based on a 'Resilience Network' which mechanistically links and controls fungal growth, metabolism, stress resistance and drug tolerance. The elasticity of this network probably underlies the phenotypic variability of fungal isolates and the heterogeneity of individual cells within clonal populations. Consequently, I suggest that the definition of the fungal Resilience Network represents an important goal for the future which offers the clear potential to reveal drug targets that compromise drug tolerance and synergise with current antifungal therapies.

Genome-wide analysis of heat stress-stimulated transposon mobility in the human fungal pathogen Cryptococcus deneoformans

We recently reported transposon mutagenesis as a significant driver of spontaneous mutations in the human fungal pathogen Cryptococcus deneoformans during murine infection. Mutations caused by transposable element (TE) insertion into reporter genes were dramatically elevated at high temperatures (37° vs. 30°) in vitro, suggesting that heat stress stimulates TE mobility in the Cryptococcus genome. To explore the genome-wide impact of TE mobilization, we generated transposon accumulation lines by in vitro passage of C. deneoformans strain XL280α for multiple generations at both 30° and at the host-relevant temperature of 37°. Utilizing whole-genome sequencing, we identified native TE copies and mapped multiple de novo TE insertions in these lines. Movements of the T1 DNA transposon occurred at both temperatures with a strong bias for insertion between gene-coding regions. By contrast, the Tcn12 retrotransposon integrated primarily within genes and movement occurred exclusively at 37°. In addition, we observed a dramatic amplification in copy number of the Cnl1 (Cryptococcus neoformans LINE-1) retrotransposon in subtelomeric regions under heat-stress conditions. Comparing TE mutations to other sequence variations detected in passaged lines, the increase in genomic changes at elevated temperatures was primarily due to mobilization of the retroelements Tcn12 and Cnl1. Finally, we found multiple TE movements (T1, Tcn12, and Cnl1) in the genomes of single C. deneoformans isolates recovered from infected mice, providing evidence that mobile elements are likely to facilitate microevolution and rapid adaptation during infection.

Similar evolutionary trajectories in an environmental Cryptococcus neoformans isolate after human and murine infection

A pet cockatoo was the suspected source of Cryptococcus neoformans recovered from an immunocompromised patient with cryptococcosis based on molecular analyses available in 2000. Here, we report whole genome sequence analysis of the clinical and cockatoo strains. Both are closely related MATα strains belonging to the VNII lineage, confirming that the human infection likely originated from pet bird exposure. The two strains differ by 61 single nucleotide polymorphisms, including eight nonsynonymous changes involving seven genes. To ascertain whether changes in these genes are selected for during mammalian infection, we passaged the cockatoo strain in mice. Remarkably, isolates obtained from mouse tissue possess a frameshift mutation in one of the seven genes altered in the human sample (LQVO5_000317), a gene predicted to encode an SWI-SNF chromatin-remodeling complex protein. In addition, both cockatoo and patient strains as well as mouse-passaged isolates obtained from brain tissue had a premature stop codon in a homologue of ZFC3 (LQVO5_004463), a predicted single-zinc finger containing protein, which is associated with larger capsules when deleted and reverted to a full-length protein in the mouse-passaged isolates obtained from lung tissue. The patient strain and mouse-passaged isolates show variability in virulence factors, with differences in capsule size, melanization, rates of nonlytic expulsion from macrophages, and amoeba predation resistance. Our results establish that environmental strains undergo genomic and phenotypic changes during mammalian passage, suggesting that animal virulence can be a mechanism for genetic change and that the genomes of clinical isolates may provide a readout of mutations acquired during infection.

Pol32, an accessory subunit of DNA polymerase delta, plays an essential role in genome stability and pathogenesis of Candida albicans

Candida albicans is a pathobiont that inflicts serious bloodstream fungal infections in individuals with compromised immunity and gut dysbiosis. Genomic diversity in the form of copy number alteration, ploidy variation, and loss of heterozygosity as an adaptive mechanism to adverse environments is frequently observed in C. albicans. Such genomic variations also confer a varied degree of fungal virulence and drug resistance, yet the factors propelling these are not completely understood. DNA polymerase delta (Polδ) is an essential replicative DNA polymerase in the eukaryotic cell and is yet to be characterized in C. albicans. Therefore, this study was designed to gain insights into the role of Polδ, especially its non-essential subunit Pol32, in the genome plasticity and life cycle of C. albicans. PCNA, the DNA clamp, recruits Polδ to the replication fork for processive DNA replication. Unlike in Saccharomyces cerevisiae, the PCNA interaction protein (PIP) motif of CaPol32 is critical for Polδ's activity during DNA replication. Our comparative genetic analyses and whole-genome sequencing of POL32 proficient and deficient C. albicans cells revealed a critical role of Pol32 in DNA replication, cell cycle progression, and genome stability as SNPs, indels, and repeat variations were largely accumulated in pol32 null strain. The loss of pol32 in C. albicans conferred cell wall deformity; Hsp90 mediated azoles resistance, biofilm development, and a complete attenuation of virulence in an animal model of systemic candidiasis. Thus, although Pol32 is dispensable for cell survival, its function is essential for C. albicans pathogenesis; and we discuss its translational implications in antifungal drugs and whole-cell vaccine development.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Impact of changes at the Candida albicans cell surface upon immunogenicity and colonisation in the gastrointestinal tract

The immunogenicity of Candida albicans cells is influenced by changes in the exposure of microbe-associated molecular patterns (MAMPs) on the fungal cell surface. Previously, the degree of exposure on the C. albicans cell surface of the immunoinflammatory MAMP β-(1,3)-glucan was shown to correlate inversely with colonisation levels in the gastrointestinal (GI) tract. This is important because life-threatening systemic candidiasis in critically ill patients often arises from translocation of C. albicans strains present in the patient's GI tract. Therefore, using a murine model, we have examined the impact of gut-related factors upon β-glucan exposure and colonisation levels in the GI tract. The degree of β-glucan exposure was examined by imaging flow cytometry of C. albicans cells taken directly from GI compartments, and compared with colonisation levels. Fungal β-glucan exposure was lower in the cecum than the small intestine, and fungal burdens were correspondingly higher in the cecum. This inverse correlation did not hold for the large intestine. The gut fermentation acid, lactate, triggers β-glucan masking in vitro, leading to attenuated anti-Candida immune responses. Additional fermentation acids are present in the GI tract, including acetate, propionate, and butyrate. We show that these acids also influence β-glucan exposure on C. albicans cells in vitro and, like lactate, they influence β-glucan exposure via Gpr1/Gpa2-mediated signalling. Significantly, C. albicans gpr1Δ gpa2Δ cells displayed elevated β-glucan exposure in the large intestine and a corresponding decrease in fungal burden, consistent with the idea that Gpr1/Gpa2-mediated β-glucan masking influences colonisation of this GI compartment. Finally, extracts from the murine gut and culture supernatants from the mannan grazing gut anaerobe Bacteroides thetaiotaomicron promote β-glucan exposure at the C. albicans cell surface. Therefore, the local microbiota influences β-glucan exposure levels directly (via mannan grazing) and indirectly (via fermentation acids), whilst β-glucan masking appears to promote C. albicans colonisation of the murine large intestine.

Recent insights into the evolution of mutation rates in yeast

Mutation is the origin of all genetic variation, good and bad. The mutation process can evolve in response to mutations, positive or negative selection, and genetic drift, but how these forces contribute to mutation-rate variation is an unsolved problem at the heart of genetics research. Mutations can be challenging to measure, but genome sequencing and other tools have allowed for the collection of larger and more detailed datasets, particularly in the yeast-model system. We review key hypotheses for the evolution of mutation rates and describe recent advances in understanding variation in mutational properties within and among yeast species. The multidimensional spectrum of mutations is increasingly recognized as holding valuable clues about how this important process evolves.

Acquisition of cross-azole tolerance and aneuploidy in Candida albicans strains evolved to posaconazole

A number of in vitro studies have examined the acquisition of drug resistance to the triazole fluconazole, a first-line treatment for many Candida infections. Much less is known about posaconazole, a newer triazole. We conducted the first in vitro experimental evolution of replicates from 8 diverse strains of Candida albicans in a high level of the fungistatic drug posaconazole. Approximately half of the 132 evolved replicates survived 50 generations of evolution, biased toward some of the strain backgrounds. We found that although increases in drug resistance were rare, increases in drug tolerance (the slow growth of a subpopulation of cells in a level of drug above the resistance level) were common across strains. We also found that adaptation to posaconazole resulted in widespread cross-tolerance to other azole drugs. Widespread aneuploidy was observed in evolved replicates from some strain backgrounds. Trisomy of at least one of chromosomes 3, 6, and R was identified in 11 of 12 whole-genome sequenced evolved SC5314 replicates. These findings document rampant evolved cross-tolerance among triazoles and highlight that increases in drug tolerance can evolve independently of drug resistance in a diversity of C. albicans strain backgrounds.

Host's Immunity and Candida Species Associated with Denture Stomatitis: A Narrative Review

Denture-related Candida stomatitis, which has been described clinically in the literature, is either localized or generalized inflammation of the oral mucosa in connection with a removable prosthesis. During this inflammatory process, the mycobacterial biofilm and the host’s immune response play an essential role. Among microorganisms of this mixed biofilm, the Candida species proliferates easily and changes from a commensal to an opportunistic pathogen. In this situation, the relationship between the Candida spp. and the host is influenced by the presence of the denture and conditioned both by the immune response and the oral microbiota. Specifically, this fungus is able to hijack the innate immune system of its host to cause infection. Additionally, older edentulous wearers of dentures may experience an imbalanced and decreased oral microbiome diversity. Under these conditions, the immune deficiency of these aging patients often promotes the spread of commensals and pathogens. The present narrative review aimed to analyze the innate and adaptive immune responses of patients with denture stomatitis and more particularly the involvement of Candida albicans sp. associated with this pathology.

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.

Increased Virulence and Large-Scale Reduction in Genome Size of Tetraploid Candida albicans Evolved in Nematode Hosts

Candida albicans is an opportunistic fungal pathogen of humans, yet the within-host dynamics of C. albicans infection are not clear. While C. albicans is commonly diploid, it exhibits a range of ploidies, including tetraploidy. Previous work found that tetraploid C. albicans populations exhibited rapid adaptation and significant genome instability when evolved in vitro. Host immune function alters the rate and magnitude of C. albicans virulence evolution, but the effects of the host immunity on tetraploid C. albicans populations are unclear. Here, we tested the effects of the host immunity on genome stability and virulence evolution of tetraploid C. albicans using experimental evolution. We selected for C. albicans increased virulence within either immunocompetent or immunocompromised Caenorhabditis elegans hosts. After nine passages we observed a response to selection for increased virulence. Both populations exposed to either immunocompetent or immunocompromised hosts increased virulence after passage through C. elegans hosts. However, the C. albicans populations passaged through immunocompetent hosts under selection exhibited unique temporal dynamics, a rapid increase in virulence and then subsequent loss of virulence. Most C. albicans populations exhibited genome size reduction within six passages, however populations exposed to immunocompetent hosts exhibited the most rapid transition to ~diploid. Therefore, we found that tetraploids rapidly increase in virulence and decrease genome size within host environments. Further, the combination of selection for greater virulence in the presence of immunocompetent hosts results in major virulence fluctuations and genome size changes. Thus, host immunity significantly impacts the evolutionary trajectories of tetraploid C. albicans.

Multiple Stochastic Parameters Influence Genome Dynamics in a Heterozygous Diploid Eukaryotic Model

The heterozygous diploid genome of Candida albicans displays frequent genomic rearrangements, in particular loss-of-heterozygosity (LOH) events, which can be seen on all eight chromosomes and affect both laboratory and clinical strains. LOHs, which are often the consequence of DNA damage repair, can be observed upon stresses reminiscent of the host environment, and result in homozygous regions of various sizes depending on the molecular mechanisms at their origins. Recent studies have shed light on the biological importance of these frequent and ubiquitous LOH events in C. albicans. In diploid Saccharomyces cerevisiae, LOH facilitates the passage of recessive beneficial mutations through Haldane’s sieve, allowing rapid evolutionary adaptation. This also appears to be true in C. albicans, where the full potential of an adaptive mutation is often only observed upon LOH, as illustrated in the case of antifungal resistance and niche adaptation. To understand the genome-wide dynamics of LOH events in C. albicans, we constructed a collection of 15 strains, each one carrying a LOH reporter system on a different chromosome arm. This system involves the insertion of two fluorescent marker genes in a neutral genomic region on both homologs, allowing spontaneous LOH events to be detected by monitoring the loss of one of the fluorescent markers using flow cytometry. Using this collection, we observed significant LOH frequency differences between genomic loci in standard laboratory growth conditions; however, we further demonstrated that comparable heterogeneity was also observed for a given genomic locus between independent strains. Additionally, upon exposure to stress, three outcomes could be observed in C. albicans, where individual strains displayed increases, decreases, or no effect of stress in terms of LOH frequency. Our results argue against a general stress response triggering overall genome instability. Indeed, we showed that the heterogeneity of LOH frequency in C. albicans is present at various levels, inter-strain, intra-strain, and inter-chromosomes, suggesting that LOH events may occur stochastically within a cell, though the genetic background potentially impacts genome stability in terms of LOH throughout the genome in both basal and stress conditions. This heterogeneity in terms of genome stability may serve as an important adaptive strategy for the predominantly clonal human opportunistic pathogen C. albicans, by quickly generating a wide spectrum of genetic variation combinations potentially permitting subsistence in a rapidly evolving environment.

Evolution of the human pathogenic lifestyle in fungi

Fungal pathogens cause more than a billion human infections every year, resulting in more than 1.6 million deaths annually. Understanding the natural history and evolutionary ecology of fungi is helping us understand how disease-relevant traits have repeatedly evolved. Different types and mechanisms of genetic variation have contributed to the evolution of fungal pathogenicity and specific genetic differences distinguish pathogens from non-pathogens. Insights into the traits, genetic elements, and genetic and ecological mechanisms that contribute to the evolution of fungal pathogenicity are crucial for developing strategies to both predict emergence of fungal pathogens and develop drugs to combat them.

A chromosomal-level reference genome of the widely utilized Coccidioides posadasii laboratory strain "Silveira"

Coccidioidomycosis is a common fungal disease that is endemic to arid and semi-arid regions of both American continents. Coccidioides immitis and Coccidioides posadasii are the etiological agents of the disease, also known as Valley Fever. For several decades, the C. posadasii strain Silveira has been used widely in vaccine studies, is the source strain for production of diagnostic antigens, and is a widely used experimental strain for functional studies. In 2009, the genome was sequenced using Sanger sequencing technology, and a draft assembly and annotation were made available. In this study, the genome of the Silveira strain was sequenced using single molecule real-time sequencing PacBio technology, assembled into chromosomal-level contigs, genotyped, and the genome was reannotated using sophisticated and curated in silico tools. This high-quality genome sequencing effort has improved our understanding of chromosomal structure, gene set annotation, and lays the groundwork for identification of structural variants (e.g. transversions, translocations, and copy number variants), assessment of gene gain and loss, and comparison of transposable elements in future phylogenetic and population genomics studies.

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.

Clonal evolution in serially passaged Cryptococcus neoformans × deneoformans hybrids reveals a heterogenous landscape of genomic change

Cryptococcus neoformans × deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high (∼10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in sub-telomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments.

Candida albicans Isolates 529L and CHN1 Exhibit Stable Colonization of the Murine Gastrointestinal Tract

Candida albicans is a pathobiont that colonizes multiple niches in the body including the gastrointestinal (GI) tract but is also responsible for both mucosal and systemic infections. Despite its prevalence as a human commensal, the murine GI tract is generally refractory to colonization with the C. albicans reference isolate SC5314. Here, we identify two C. albicans isolates, 529L and CHN1, that stably colonize the murine GI tract in three different animal facilities under conditions where SC5314 is lost from this niche. Analysis of the bacterial microbiota did not show notable differences among mice colonized with the three C. albicans strains. We compared the genotypes and phenotypes of these three strains and identified thousands of single nucleotide polymorphisms (SNPs) and multiple phenotypic differences, including their ability to grow and filament in response to nutritional cues. Despite striking filamentation differences under laboratory conditions, however, analysis of cell morphology in the GI tract revealed that the three isolates exhibited similar filamentation properties in this in vivo niche. Notably, we found that SC5314 is more sensitive to the antimicrobial peptide CRAMP, and the use of CRAMP-deficient mice modestly increased the ability of SC5314 to colonize the GI tract relative to CHN1 and 529L. These studies provide new insights into how strain-specific differences impact C. albicans traits in the host and advance CHN1 and 529L as relevant strains to study C. albicans pathobiology in its natural host niche. IMPORTANCE Understanding how fungi colonize the GI tract is increasingly recognized as highly relevant to human health. The animal models used to study Candida albicans commensalism commonly rely on altering the host microbiome (via antibiotic treatment or defined diets) to establish successful GI colonization by the C. albicans reference isolate SC5314. Here, we characterize two C. albicans isolates that can colonize the murine GI tract without antibiotic treatment and can therefore be used as tools for studying fungal commensalism. Importantly, experiments were replicated in three different animal facilities and utilized three different mouse strains. Differential colonization between fungal isolates was not associated with alterations in the bacterial microbiome but rather with distinct responses to CRAMP, a host antimicrobial peptide. This work emphasizes the importance of C. albicans intraspecies variation as well as host antimicrobial defense mechanisms in defining the outcome of commensal interactions.

Host defense mechanisms induce genome instability leading to rapid evolution in an opportunistic fungal pathogen

The ability to generate genetic variation facilitates rapid adaptation in stressful environments. The opportunistic fungal pathogen Candida albicans frequently undergoes large-scale genomic changes, including aneuploidy and loss-of heterozygosity (LOH), following exposure to host environments. However, the specific host factors inducing C. albicans genome instability remain largely unknown. Here, we leveraged the genetic tractability of nematode hosts to investigate whether innate immune components, including antimicrobial peptides (AMPs) and reactive oxygen species (ROS), induced host-associated C. albicans genome instability. C. albicans associated with immunocompetent hosts carried multiple large-scale genomic changes including LOH, whole chromosome, and segmental aneuploidies. In contrast, C. albicans associated with immunocompromised hosts deficient in AMPs or ROS production had reduced LOH frequencies and fewer, if any, additional genomic changes. To evaluate if extensive host-induced genomic changes had long-term consequences for C. albicans adaptation, we experimentally evolved C. albicans in either immunocompetent or immunocompromised hosts and selected for increased virulence. C. albicans evolved in immunocompetent hosts rapidly increased virulence, but not in immunocompromised hosts. Taken together, this work suggests that host-produced ROS and AMPs induces genotypic plasticity in C. albicans which facilitates rapid evolution.

Heterogeneous Mutation Rates and Spectra in Yeast Hybrids

Mutation rates and spectra vary between species and among populations. Hybridization can contribute to this variation, but its role remains poorly understood. Estimating mutation rates requires controlled conditions where the effect of natural selection can be minimized. One way to achieve this is through mutation accumulation experiments coupled with genome sequencing. Here, we investigate 400 mutation accumulation lines initiated from 11 genotypes spanning intralineage, interlineage, and interspecific crosses of the yeasts Saccharomyces paradoxus and S. cerevisiae and propagated for 770 generations. We find significant differences in mutation rates and spectra among crosses, which are not related to the level of divergence of parental strains but are specific to some genotype combinations. Differences in number of generations and departures from neutrality play a minor role, whereas polyploidy and loss of heterozygosity impact mutation rates in some of the hybrid crosses in an opposite way.

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. IMPORTANCE Candida 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.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

Extreme diversification driven by parallel events of massive loss of heterozygosity in the hybrid lineage of Candida albicans

Candida albicans is the most commonly reported species causing candidiasis. The taxonomic classification of C. albicans and related lineages is controversial, with Candida africana (syn. C. albicans var. africana) and Candida stellatoidea (syn. C. albicans var. stellatoidea) being considered different species or C. albicans varieties depending on the authors. Moreover, recent genomic analyses have suggested a shared hybrid origin of C. albicans and C. africana, but the potential parental lineages remain unidentified. Although the genomes of C. albicans and C. africana have been extensively studied, the genome of C. stellatoidea has not been sequenced so far. In order to get a better understanding of the evolution of the C. albicans clade, and to assess whether C. stellatoidea could represent one of the unknown C. albicans parental lineages, we sequenced C. stellatoidea type strain (CBS 1905). This genome was compared to that of C. albicans and of the closely related lineage C. africana. Our results show that, similarly to C. africana, C. stellatoidea descends from the same hybrid ancestor as other C. albicans strains and that it has undergone a parallel massive loss of heterozygosity.

Tunicamycin Potentiates Antifungal Drug Tolerance via Aneuploidy in Candida albicans

How cells exposed to one stress are later able to better survive other types of stress is not well understood. In eukaryotic organisms, physiological and pathological stresses can disturb endoplasmic reticulum (ER) function, resulting in "ER stress." Here, we found that exposure to tunicamycin, an inducer of ER stress, resulted in the acquisition of a specific aneuploidy, chromosome 2 trisomy (Chr2x3), in Candida albicans. Importantly, the resulting aneuploidy also conferred cross-tolerance to caspofungin, a first-line echinocandin antifungal, as well as to hydroxyurea, a common chemotherapeutic agent. Exposure to a range of tunicamycin concentrations induced similar ER stress responses. Extra copies of one Chr2 gene, MKK2, affected both tunicamycin and caspofungin tolerance, while at least 3 genes on chromosome 2 (ALG7, RTA2, and RTA3) affected only tunicamycin and not caspofungin responses. Other Chr2 genes (RNR1 and RNR21) affected hydroxyurea tolerance but neither tunicamycin nor caspofungin tolerance. Deletion of components of the protein kinase C (PKC) or calcineurin pathways affected tolerance to both tunicamycin and caspofungin, supporting the idea that the ER stress response and echinocandin tolerance are regulated by overlapping stress response pathways. Thus, antifungal drug tolerance can arise rapidly via ER stress-induced aneuploidy. IMPORTANCE Candida albicans is a prevalent human fungal commensal and also a pathogen that causes life-threatening systemic infections. Treatment failures are frequent because few therapeutic antifungal drug classes are available and because drug resistance and tolerance limit drug efficacy. We found that C. albicans rapidly overcomes the cellular stress induced by the drug tunicamycin by duplicating chromosome 2. Also, chromosome 2 duplication confers tolerance not only to tunicamycin but also to the following two unrelated drugs: caspofungin, an antifungal drug, and hydroxyurea, a chemotherapeutic. Cross tolerance to the three drugs involves different sets of genes, although some genetic pathways affect the tolerance to two of these three drugs. This work highlights a serious concern, namely, that changes in whole chromosome copy number can occur in response to one type of stress, and yet, they may facilitate the emergence of tolerance to multiple drugs, including the few antifungal drug classes available to treat Candida infections.

Accounting for the Biological Complexity of Pathogenic Fungi in Phylogenetic Dating

In the study of pathogen evolution, temporal dating of phylogenies provides information on when species and lineages may have diverged in the past. When combined with spatial and epidemiological data in phylodynamic models, these dated phylogenies can also help infer where and when outbreaks occurred, how pathogens may have spread to new geographic locations and/or niches, and how virulence or drug resistance has developed over time. Although widely applied to viruses and, increasingly, to bacterial pathogen outbreaks, phylogenetic dating is yet to be widely used in the study of pathogenic fungi. Fungi are complex organisms with several biological processes that could present issues with appropriate inference of phylogenies, clock rates, and divergence times, including high levels of recombination and slower mutation rates although with potentially high levels of mutation rate variation. Here, we discuss some of the key methodological challenges in accurate phylogeny reconstruction for fungi in the context of the temporal analyses conducted to date and make recommendations for future dating studies to aid development of a best practices roadmap in light of the increasing threat of fungal outbreaks and antifungal drug resistance worldwide.

Functional Roles of Homologous Recombination and Non-Homologous End Joining in DNA Damage Response and Microevolution in Cryptococcus neoformans

DNA double-strand breaks (DSBs) are the most deleterious type of DNA lesions because they cause loss of genetic information if not properly repaired. In eukaryotes, homologous recombination (HR) and non-homologous end joining (NHEJ) are required for DSB repair. However, the relationship of HR and NHEJ in DNA damage stress is unknown in the radiation-resistant fungus Cryptococcus neoformans. In this study, we found that the expression levels of HR- and NHEJ-related genes were highly induced in a Rad53-Bdr1 pathway-dependent manner under genotoxic stress. Deletion of RAD51, which is one of the main components in the HR, resulted in growth under diverse types of DNA damage stress, whereas perturbations of KU70 and KU80, which belong to the NHEJ system, did not affect the genotoxic stresses except when bleomycin was used for treatment. Furthermore, deletion of both RAD51 and KU70/80 renders cells susceptible to oxidative stress. Notably, we found that deletion of RAD51 induced a hypermutator phenotype in the fluctuation assay. In contrast to the fluctuation assay, perturbation of KU70 or KU80 induced rapid microevolution similar to that induced by the deletion of RAD51. Collectively, Rad51-mediated HR and Ku70/Ku80-mediated NHEJ regulate the DNA damage response and maintain genome stability.