The Candida Genome Database (CGD): incorporation of Assembly 22, systematic identifiers and visualization of high throughput sequencing data

Hotspot gene conversion between FKS1 and FKS2 in echinocandin resistant Candida glabrata serial isolates

Candida glabrata (Nakaseomyces glabratus) is the most common cause of drug-resistant candidemia and is associated with a high mortality rate. Only a few mechanisms of drug resistance are known in C. glabrata, predominantly involving recurrent single nucleotide polymorphisms. The importance of structural variation in acquired drug resistance is not understood. We performed comparative phenotypic and genomic analyses of six serial bloodstream isolates of C. glabrata and identified novel mutations associated with resistance to echinocandins. Critically, we identified a novel gene conversion event between the hotspot 2 regions of FKS1 and FKS2 that was associated with increased resistance to micafungin. We further analyzed 621 publicly available C. glabrata genomes and found three additional examples of structural variation involving FKS1 and FKS2. Ultimately, drug resistance in C. glabrata involves structural variants that are missed with current diagnostic methods and need to be considered when designing and implementing more effective antifungal management strategies.

Transcriptional control of C. albicans white-opaque switching and modulation by environmental cues and strain background

The opportunistic fungal pathogen Candida albicans can undergo cellular transitions in response to environmental cues that impact its lifestyle and its interactions with the human host. This is exemplified by the white-opaque switch, which is a heritable transition between two phenotypic states that is regulated by a highly interconnected network of transcription factors (TFs). To obtain greater understanding of the transcriptional regulation of the switch, we generated a genome-wide, tetracycline-inducible TF library in the WO-1 strain background and identified those TFs whose forced expression induces white cells to switch to the opaque state. This set of opaque-inducing TFs was also evaluated for their ability to induce switching in a second strain background, that of the standard reference strain SC5314, as well as during growth on different laboratory media. These experiments identify 14 TFs that can drive white-to-opaque switching when overexpressed but that do so in a highly strain- and media-specific manner. In particular, changes in pH, amino acids, and zinc concentrations had marked effects on the ability of TFs to drive phenotypic switching. These results provide insights into the complex transcriptional regulation of switching in C. albicans and reveal that an interplay between genetic and environmental factors determines TF function and cell fate.IMPORTANCEThe white-opaque switch in Candida albicans represents a model system for understanding an epigenetic switch in a eukaryotic pathogen. Here, we generated an inducible library of the set of transcription factors (TFs) present in C. albicans and identify 14 TFs that can drive the white-to-opaque transition when ectopically expressed. We demonstrate that several of these TFs induce the switch in a highly strain- and media-specific manner. This highlights that both strain background and changes in experimental conditions (including different water sources) can profoundly impact the phenotypic consequences of TF overexpression. Moreover, the inducible TF library provides an invaluable tool for the further analysis of TF function in this important human pathogen.

Aneuploidy confers a unique transcriptional and phenotypic profile to Candida albicans

Inaccurate chromosome segregation can lead to the formation of aneuploid cells that harbor an imbalanced complement of chromosomes. Several fungal species are not only able to tolerate the detrimental effects of aneuploidy but can use it to adapt to environmental pressures. The fungal pathobiont Candida albicans frequently acquires supernumerary chromosomes that enable growth in the presence of antifungal drugs or in specific host niches, yet the transcriptional changes associated with aneuploidy are not globally defined. Here, a karyotypically diverse set of C. albicans strains revealed that expression generally correlated with gene copy number regardless of the strain karyotype. Unexpectedly, aneuploid strains shared a characteristic transcriptional profile that was distinct from a generalized environmental stress response previously defined in aneuploid yeast cells. This aneuploid transcriptional response led to altered growth and oxidative balances relative to euploid control strains. The increased expression of reactive oxygen species (ROS) mitigating enzymes in aneuploid cells reduced the levels of ROS but caused an acute sensitivity to both internal and external sources of oxidative stress. Taken together, our work demonstrates common transcriptional and phenotypic features of aneuploid C. albicans cells with consequences for infection of different host niches and susceptibility to environmental stimuli.

Protocol to identify genes important for Candida albicans fitness in diverse environmental conditions using pooled bar-seq screening approach

Identifying genes important for fitness in Candida albicans advances our understanding of this important pathogen of humans. Here, we present a functional genomics approach for assessing fitness through the quantification of strain-specific barcodes. We describe steps for library preparation, propagation of strains, genomic DNA extraction, amplification of barcodes, and sequencing. We then detail the computational analysis of data to determine effect size and statistical significance. For complete details on the use and execution of this protocol, please refer to Xiong et al.1.

The Candida Genome Database: annotation and visualization updates

The Candida Genome Database (CGD; www.candidagenome.org) is unique in being both a model organism database and a fungal pathogen database. As a fungal pathogen database, CGD hosts locus pages for 5 species of the best-studied pathogenic fungi in the Candida group. As a model organism database, the species Candida albicans serves as a model both for other Candida spp. and for non-Candida fungi that form biofilms and undergo routine morphogenic switching from the planktonic form to the filamentous form, which is not done by other model yeasts. As pathogenic Candida species have become increasingly drug resistant, the high lethality of invasive candidiasis in immunocompromised people is increasingly alarming. There is a pressing need for additional research into basic Candida biology, epidemiology and phylogeny, and potential new antifungals. CGD serves the needs of this diverse research community by curating the entire gene-based Candida experimental literature as it is published, extracting, organizing, and standardizing gene annotations. Gene pages were added for the species Candida auris, recent outbreaks of which have been labeled an "urgent" threat. Most recently, we have begun linking clinical data on disease to relevant Literature Topics to improve searchability for clinical researchers. Because CGD curates for multiple species and most research focuses on aspects related to pathogenicity, we focus our curation efforts on assigning Literature Topic tags, collecting detailed mutant phenotype data, and assigning controlled Gene Ontology terms with accompanying evidence codes. Our Summary pages for each feature include the primary name and all aliases for that locus, a description of the gene and/or gene product, detailed ortholog information with links, a JBrowse window with a visual view of the gene on its chromosome, summarized phenotype, Gene Ontology, and sequence information, references cited on the summary page itself, and any locus notes. The database serves as a community hub, where we link to various types of reference material of relevance to Candida researchers, including colleague information, news, and notice of upcoming meetings. We routinely survey the community to learn how the field is evolving and how needs may have changed. For example, we asked our users which species we should next add to CGD, and the clear answer was Candida tropicalis. A key future challenge is management of the flood of high-throughput expression data to make it as useful as possible to as many researchers as possible. The central challenge for any community database is to turn data into knowledge, which the community can access, use, and build upon.

PCR-based CRISPR/Cas9 system for fluorescent tagging: A tool for studying Candida parapsilosis virulence

Candida parapsilosis is persistent in a hospital environment hence it is often associated with nosocomial infections especially amongst low-birth weight neonates. Genetic modification is therefore important to characterise the physiological and virulence related properties of this fungus. A PCR-based CRISPR/Cas9 system has been adopted to facilitate the generation of fluorescent tagged prototroph isolates. We examined a total of eight fluorescent protein coding genes, out of which three were found to be applicable for simultaneous utilisation. We investigated three clinical isolates of C. parapsilosis in terms of their adherence to silicone and their uptake by J774.2 murine macrophages in competition assays. Interestingly, we found significant differences between them in both experiments where GA1 isolate was significantly less resistant to macrophage uptake and CDC317 was significantly more adherent to silicone material. In silico analysis of the agglutinin-like sequences (Als) exposed remarkable diversity in this protein family and additionally, the thorough analysis of the ALS genes revealed evidence of formation of a new gene by intrachromosomal recombination in the GA1 isolate. Finally, we provide a step by step protocol for the application of the PCR-based CRISPR/Cas9 system for fluorescently labelling C. parapsilosis isolates.

In vitro long-term exposure to chlorhexidine or triclosan induces cross-resistance against azoles in Nakaseomyces glabratus

BACKGROUND

Topical antiseptics are crucial for preventing infections and reducing transmission of pathogens. However, commonly used antiseptic agents have been reported to cause cross-resistance to other antimicrobials in bacteria, which has not yet been described in yeasts. This study aims to assess the in vitro efficacy of antiseptics against clinical and reference isolates of Candida albicans and Nakaseomyces glabratus, and whether prolonged exposure to antiseptics promotes the development of antifungal (cross)resistance.

METHODS

A high-throughput approach for in vitro resistance development was established to simultaneously expose 96 C. albicans and N. glabratus isolates to increasing concentrations of a given antiseptic - chlorhexidine, triclosan or octenidine. Susceptibility testing and whole genome sequencing of yeast isolates pre- and post-exposure were performed.

RESULTS

Long-term exposure to antiseptics does not result in the development of stable resistance to the antiseptics themselves. However, 50 N. glabratus isolates acquired resistance to azole antifungals after long-term exposure to triclosan or chlorhexidine, revealing newly acquired mutations in the PDR1 and PMA1 genes.

CONCLUSIONS

Chlorhexidine as well as triclosan, but not octenidine, were able to introduce selective pressure promoting resistance to azole antifungals. Although we assessed this phenomenon only in vitro, these findings warrant critical monitoring in clinical settings.

Draft genome sequence of a non-human primate-derived isolate of Candida parapsilosis

Candida parapsilosis is a common human commensal and opportunistic fungal pathogen that is also found in non-human primates (NHPs). Here, we report the first draft sequence of C. parapsilosis NCYC 4418, a fecal isolate from an adult cynomolgus macaque.

Structural insights into mechanisms of zinc scavenging by the Candida albicans zincophore Pra1

Candida albicans causes more than 400,000 life-threatening, and half a billion mucosal infections annually. In response to infection, the host limits availability of essential micronutrients, including zinc, to restrict growth of the invading pathogen. As assimilation of zinc is essential for C. albicans pathogenicity, its limitation induces the secretion of the zincophore protein Pra1 to scavenge zinc from the host. Pra1 also plays a number of important roles in host-pathogen interactions and is conserved in most fungi. However, the structure of fungal zincophores is not known. Here, we present the first cryogenic electron microscopy structures of C. albicans Pra1 in its apo- and zinc-bound states, at 2.8 and 2.5 Å resolution respectively. Our work reveals a hexameric ring-like assembly with multiple zinc binding sites. Through genetic studies, we show that one of these zinc binding sites is essential for C. albicans growth under zinc restriction but does not affect the inflammatory properties of Pra1. These data provide a foundation for future work to explore the structural basis of Pra1-mediated host-pathogen interactions, C. albicans zinc uptake, as well therapeutics development.

Mitochondrial targeting of Candida albicans SPFH proteins and requirement of stomatins for SDS-induced stress tolerance

The SPFH (stomatin, prohibitin, flotillin, and HflK/HflC) protein superfamily is conserved across all domains of life. Fungal SPFH proteins are required for respiration, stress adaptation, and membrane scaffolding. In the yeast Candida albicans, stomatin-like protein 3 (Slp3) forms punctate foci at the plasma membrane, and SLP3 overexpression causes cell death following exposure to the surfactant, SDS, and the oxidative stressor, H2O2. Here, we sought to determine the cellular localization and functionally characterize stomatin-like protein 2 (Slp2), prohibitin-1 (Phb1), prohibitin-2 (Phb2), and prohibitin-12 (Phb12) in C. albicans. Cytological and western blotting results showed that Slp2-Gfp/Rfp and prohibitin-Gfp fusion proteins localize to the mitochondrion in yeast cells. Growth assay results did not identify any respiration defects in a panel of stomatin and prohibitin mutant strains, suggesting that SPFH respiratory function has diverged in C. albicans from other model eukaryotes. However, a slp2Δ/Δ/slp3Δ/Δ double mutant strain grew poorly in the presence of 0.08% SDS, accumulated intracellular reactive oxidative species, and displayed aberrant ergosterol distribution in the plasma membrane. These phenotypes were not observed in slp2Δ/Δ or slp3Δ/Δ single mutants, indicating a possible indirect genetic interaction between SLP2 and SLP3. In addition, slp2Δ/Δ and slp2Δ/Δ/slp3Δ/Δ mutant strains were slightly resistant to the antifungal drug, fluconazole. Collectively, these findings reveal the cellular localization of Slp2, Phb1, Phb2, and Phb12, highlight the significance of stomatins in C. albicans SDS stress tolerance, and, for the first time, associate stomatins with antifungal resistance.

IMPORTANCE

Stomatins and prohibitins coordinate respiration and stress adaptation in fungi. Invasive mycoses caused by Candida albicans are a significant cause of morbidity, and candidemia patients show high mortality rates worldwide. Mitochondria are essential for C. albicans commensalism and virulence, and mitochondrial proteins are targets for antifungal interventions. C. albicans encodes five SPFH proteins: two stomatin-like proteins and three prohibitins. We have previously shown that Slp3 is important for C. albicans adaptation to various types of environmental stress. Moreover, synthetic compounds that bind to mammalian prohibitins inhibit C. albicans filamentation and are fungicidal. However, there is limited information available regarding the remaining SPFH proteins. Our findings show that mitochondrial localization of SPFH proteins is conserved in C. albicans. In addition, we demonstrate the importance of stomatins in plasma membrane and mitochondrial stress tolerance.

Alpha1-antitrypsin impacts innate host-pathogen interactions with Candida albicans by stimulating fungal filamentation

Our immune system possesses sophisticated mechanisms to cope with invading microorganisms, while pathogens evolve strategies to deal with threats imposed by host immunity. Human plasma protein α1-antitrypsin (AAT) exhibits pleiotropic immune-modulating properties by both preventing immunopathology and improving antimicrobial host defence. Genetic associations suggested a role for AAT in candidemia, the most frequent fungal blood stream infection in intensive care units, yet little is known about how AAT influences interactions between Candida albicans and the immune system. Here, we show that AAT differentially impacts fungal killing by innate phagocytes. We observed that AAT induces fungal transcriptional reprogramming, associated with cell wall remodelling and downregulation of filamentation repressors. At low concentrations, the cell-wall remodelling induced by AAT increased immunogenic β-glucan exposure and consequently improved fungal clearance by monocytes. Contrastingly, higher AAT concentrations led to excessive C. albicans filamentation and thus promoted fungal immune escape from monocytes and macrophages. This underscores that fungal adaptations to the host protein AAT can differentially define the outcome of encounters with innate immune cells, either contributing to improved immune recognition or fungal immune escape.

Acquired amphotericin B resistance leads to fitness trade-offs that can be mitigated by compensatory evolution in Candida auris

Candida auris is a growing concern due to its resistance to antifungal drugs, particularly amphotericin B (AMB), detected in 30 to 60% of clinical isolates. However, the mechanisms of AMB resistance remain poorly understood. Here we investigated 441 in vitro- and in vivo-evolved C. auris lineages from 4 AMB-susceptible clinical strains of different clades. Genetic and sterol analyses revealed four major types of sterol alterations as a result of clinically rare variations in sterol biosynthesis genes ERG6, NCP1, ERG11, ERG3, HMG1, ERG10 and ERG12. In addition, aneuploidies in chromosomes 4 and 6 emerged during resistance evolution. Fitness trade-off phenotyping and mathematical modelling identified diverse strain- and mechanism-dependent fitness trade-offs. Variation in CDC25 rescued fitness trade-offs, thereby increasing the infection capacity. This possibly contributed to therapy-induced acquired AMB resistance in the clinic. Our findings highlight sterol-modulating mechanisms and fitness trade-off compensation as risks for AMB treatment failure in clinical settings.

CRISPR-GRIT: Guide RNAs with Integrated Repair Templates Enable Precise Multiplexed Genome Editing in the Diploid Fungal Pathogen Candida albicans

Candida albicans, an opportunistic fungal pathogen, causes severe infections in immunocompromised individuals. Limited classes and overuse of current antifungals have led to the rapid emergence of antifungal resistance. Thus, there is an urgent need to understand fungal pathogen genetics to develop new antifungal strategies. Genetic manipulation of C. albicans is encumbered by its diploid chromosomes requiring editing both alleles to elucidate gene function. Although the recent development of CRISPR-Cas systems has facilitated genome editing in C. albicans, large-scale and multiplexed functional genomic studies are still hindered by the necessity of cotransforming repair templates for homozygous knockouts. Here, we present CRISPR-GRIT (Guide RNAs with Integrated Repair Templates), a repair template-integrated guide RNA design for expedited gene knockouts and multiplexed gene editing in C. albicans. We envision that this method can be used for high-throughput library screens and identification of synthetic lethal pairs in both C. albicans and other diploid organisms with strong homologous recombination machinery.

SUMO-targeted Ubiquitin Ligases as crucial mediators of protein homeostasis in Candida glabrata

Candida glabrata is an opportunistic human pathogen, capable of causing severe systemic infections that are often resistant to standard antifungal treatments. To understand the importance of protein SUMOylation in the physiology and pathogenesis of C. glabrata, we earlier identified the components of SUMOylation pathway and demonstrated that the deSUMOylase CgUlp2 is essential for pathogenesis. In this work we show that the CgUlp2 is essential to maintain protein homeostasis via the SUMO-targeted ubiquitin ligase pathway. The dual loss of deSUMOylase and specific ubiquitin ligase, CgSlx8, results in heightened protein degradation, rendering the cells vulnerable to various stressors. This degradation affects crucial processes such as purine biosynthesis and compromises mitochondrial function in the mutants. Importantly, the absence of these ubiquitin ligases impedes the proliferation of C. glabrata in macrophages. These findings underscore the significance of SUMOylation and SUMO-mediated protein homeostasis as pivotal regulators of C. glabrata physiology and capacity to survive in host cells. Understanding these mechanisms could pave the way for the development of effective antifungal treatments.

Strain variation in Candida albicans glycolytic gene regulation

Central carbon metabolism is vital for the proliferation of Candida albicans, a fungus that is prominent as a commensal and pathogen. Glycolytic genes are activated by overlapping activities of the transcription factors Tye7 and Gal4, as shown by studies in the SC5314 genetic background. However, regulatory relationships can vary among C. albicans isolates. Here, we analyzed Tye7- and Gal4-related phenotypes in five diverse clinical isolates of C. albicans. We tested growth properties and gene expression impact through Nanostring profiling and, for the two strains SC5314 and P87, RNA sequencing. Our results lead to three main conclusions. First, the functional redundancy of Tye7 and Gal4 for glycolytic gene activation is preserved among all strains tested. Second, at the gene expression level, strain P87 is an outlier with regard to tye7Δ/Δ impact, and strain SC5314 is an outlier with regard to gal4Δ/Δ impact. Third, while Gal4 is well known to be dispensable for induction of the GAL1, GAL7, and GAL10 galactose-specific metabolic genes, we find that gal4Δ/Δ mutants of several strains have a mild galactose fermentation defect, as assayed by growth on galactose with the respiration inhibitor antimycin A. Our findings indicate that even a central metabolic regulatory network is subject to strain variation and illustrates an unexpected genotype-phenotype relationship.The fungal commensal and pathogen Candida albicans rely upon metabolic flexibility to colonize and infect host niches. Central carbon metabolism is governed by two regulators, Tye7 and Gal4, as defined in the reference strain SC5314. Here, we have explored the impact of Tye7 and Gal4 on carbon utilization and gene expression across five diverse C. albicans clinical isolates. Novel aspects of this study are the finding that even a central metabolic regulatory network is subject to strain variation and the observation of an unexpected mutant phenotype.

Targeting Candida albicans O-acetyl-L-homoserine sulfhydrylase (Met15p) in antifungal treatment

Fungal infections are a serious threat to public health as they are becoming increasingly frequent. A major problem stems also from a rising fungal resistance to currently available antifungal therapies, therefore novel molecular targets are highly desirable. Exploration of enzymes participating in the biosynthesis pathways of essential amino acids such as L-methionine (L-Met) may provide new insights into pharmaceutical development. The MET15 gene from Candida albicans, encoding O-acetyl-L-homoserine sulfhydrylase (Met15p), an enzyme catalyzing the second step in that pathway, was cloned and expressed in two versions: as N and C-terminal oligo-His-tagged fusion proteins. The recombinant enzymes revealed appropriate activity, and catalyzed conversion of O-acetyl-L-homoserine and a sulfide ion to produce L-homocysteine. A new RP-HPLC-DAD method, using the enzymatic reaction product pre-column derivatization with 5,5'-dithio-bis-(2-nitrobenzoic acid) was developed and used by us to determine Met15p activity. Newly synthesized compounds as well as two commercially available exhibited a Met15p inhibitory effect which was related to antifungal activity. Fungal cells' sensitivity to inhibitors depending on the presence or absence of L-Met in the medium clearly indicated Met15p targeting. Moreover, the synergistic effect of the first methionine biosynthetic enzyme affecting inhibitor and Met15p inhibitors indicate that methionine biosynthesis pathway enzymes are promising molecular targets.

The Candida auris Hog1 MAP kinase is essential for the colonization of murine skin and intradermal persistence

Candida auris, a multidrug-resistant human fungal pathogen, was first identified in 2009 in Japan. Since then, systemic C. auris infections have now been reported in more than 50 countries, with mortality rates of 30%-60%. A major contributing factor to its high inter- and intrahospital clonal transmission is that C. auris, unlike most Candida species, displays unique skin tropism and can stay on human skin for a prolonged period. However, the molecular mechanisms responsible for C. auris skin colonization, intradermal persistence, and systemic virulence are poorly understood. Here, we report that C. auris Hog1 mitogen-activated protein kinase is essential for efficient skin colonization, intradermal persistence as well as systemic virulence. RNA-seq analysis of wild-type parental and hog1Δ mutant strains revealed marked downregulation of genes involved in processes such as cell adhesion, cell wall rearrangement, and pathogenesis in hog1Δ mutant compared to the wild-type parent. Consistent with these data, we found a prominent role for Hog1 in maintaining cell wall architecture, as the hog1Δ mutant demonstrated a significant increase in cell-surface β-glucan exposure and a concomitant reduction in chitin content. Additionally, we observed that Hog1 was required for biofilm formation in vitro and fungal survival when challenged with primary murine macrophages and neutrophils ex vivo. Collectively, these findings have important implications for understanding the C. auris skin adherence mechanisms and penetration of skin epithelial layers preceding bloodstream infections.

IMPORTANCE

Candida auris is a World Health Organization fungal priority pathogen and an urgent public health threat recognized by the Centers for Disease Control and Prevention. C. auris has a unique ability to colonize human skin. It also persists on abiotic surfaces in healthcare environments for an extended period of time. These attributes facilitate the inter- and intrahospital clonal transmission of C. auris. Therefore, understanding C. auris skin colonization mechanisms is critical for infection control, especially in hospitals and nursing homes. However, despite its profound clinical relevance, the molecular and genetic basis of C. auris skin colonization mechanisms are poorly understood. Herein, we present data on the identification of the Hog1 MAP kinase as a key regulator of C. auris skin colonization. These findings lay the foundation for further characterization of unique mechanisms that promote fungal persistence on human skin.

Biofilm matrix regulation by Candida glabrata Zap1 under acidic conditions: transcriptomic and proteomic analyses

The vaginal acidic environment potentiates the formation of Candida glabrata biofilms, leading to complicated and recurrent infections. Importantly, the production of matrix is known to contribute to the recalcitrant features of Candida biofilms. In this study, we reveal that Zap1 regulates the matrix of C. glabrata acidic biofilms and analyzed the modulation of their transcriptome (by microarrays) and matrix proteome (by LC-MS/MS) by Zap1. For that, the deletion mutant zap1Δ and its complemented strain zap1Δ::ZAP1 were constructed, and their biofilms were developed at pH 4 (adjusted with lactic acid). The results revealed that Zap1 is a negative regulator of the total amount of protein and carbohydrate in the biofilm matrix. Accordingly, various genes and matrix proteins with predicted functions in the regulation of carbohydrate metabolism, sugar binding, sugar transport, and adhesion (including Epa family) were repressed by Zap1. Nevertheless, the results also suggested that Zap1 is essential to the delivery and organization of some matrix components. Indeed, Zap1 was required for the secretion of 122 proteins to the matrix and induced the expression of 557 genes, including various targets involved in glucan metabolism. Additionally, Zap1 induced targets with roles in virulence, resistance to antifungals, and host immunity evasion, including yapsins, ERG family, and moonlighting proteins. Zap1 was also required for the secretion of acidic-specific matrix proteins, indicating a contribution to the response to the acidic environment. Overall, this study demonstrates that Zap1 is a relevant regulator of the biofilm matrix, contributing to a better understanding of C. glabrata acidic biofilms.IMPORTANCEThe rising prevalence of vulvovaginal candidiasis (VVC) and the increasing presence of Candida spp. with aggressive virulence features and low susceptibility to common antifungals, particularly Candida glabrata, have resulted in more severe, prolonged, and recurrent cases of VVC, with significant implications for patients. This research offers valuable insights into the molecular changes that contribute to the formation of C. glabrata biofilms in the acidic vaginal environment, representing a significant advancement in the understanding of C. glabrata's virulence. Notably, this study identified Zap1 as a critical regulator of C. glabrata biofilm matrix, with additional potential roles in adhesion, antifungal resistance, evasion of host immunity, and response to acidic conditions, making it a promising target for new therapeutic approaches. Importantly, Zap1 is the first regulator of the biofilm matrix to be identified in C. glabrata, and the elucidation of its targets (including genes and matrix proteins) lays a strong foundation for future research.

Collateral sensitivity counteracts the evolution of antifungal drug resistance in Candida auris

Antifungal drug resistance represents a serious global health threat, necessitating new treatment strategies. Here we investigated collateral sensitivity (CS), in which resistance to one drug increases sensitivity to another, and cross-resistance (XR), in which one drug resistance mechanism reduces susceptibility to multiple drugs, since CS and XR dynamics can guide treatment design to impede resistance development, but have not been systematically explored in pathogenic fungi. We used experimental evolution and mathematical modelling of Candida auris population dynamics during cyclic and combined drug exposures and found that especially CS-based drug cycling can effectively prevent the emergence of drug resistance. In addition, we found that a CS-based treatment switch can actively select against or eradicate resistant sub-populations, highlighting the potential to consider CS in therapeutic decision-making upon resistance detection. Furthermore, we show that some CS trends are robust among different strains and resistance mechanisms. Overall, these findings provide a promising direction for improved antifungal treatment approaches.

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.

Functional redundancy in Candida auris cell surface adhesins crucial for cell-cell interaction and aggregation

Candida auris is an emerging nosocomial fungal pathogen associated with life-threatening invasive disease due to its persistent colonization, high level of transmissibility and multi-drug resistance. Aggregative and non-aggregative growth phenotypes for C. auris strains with different biofilm forming abilities, drug susceptibilities and virulence characteristics have been described. Using comprehensive transcriptional analysis we identified key cell surface adhesins that were highly upregulated in the aggregative phenotype during in vitro and in vivo grown biofilms using a mouse model of catheter infection. Phenotypic and functional evaluations of generated null mutants demonstrated crucial roles for the adhesins Als4112 and Scf1 in mediating cell-cell adherence, coaggregation and biofilm formation. While individual mutants were largely non-aggregative, in combination cells were able to co-adhere and aggregate, as directly demonstrated by measuring cell adhesion forces using single-cell atomic force spectroscopy. This co-adherence indicates their role as complementary adhesins, which despite their limited similarity, may function redundantly to promote cell-cell interaction and biofilm formation. Functional diversity of cell wall proteins may be a form of regulation that provides the aggregative phenotype of C. auris with flexibility and rapid adaptation to the environment, potentially impacting persistence and virulence.

Nosocomial transmission of fluconazole-resistant Candida glabrata bloodstream isolates revealed by whole-genome sequencing

The clonal transmission of fluconazole-resistant Candida glabrata isolates within hospitals has seldom been analyzed by whole-genome sequencing (WGS). We performed WGS on 79 C. glabrata isolates, comprising 31 isolates from three premature infants with persistent C. glabrata bloodstream infection despite antifungal treatment in the same neonatal intensive care unit (NICU) in 2022 and 48 (27 fluconazole-resistant and 21 fluconazole-susceptible dose-dependent) bloodstream isolates from 48 patients in 15 South Korean hospitals from 2010 to 2022. Phylogenetic analysis based on WGS single-nucleotide polymorphisms (SNPs) distinguished the 79 isolates according to multilocus sequence typing (MLST) (17 sequence type [ST]3, 13 ST7, two ST22, 41 ST26, four ST55, and two ST59 isolates) and unveiled two possible clusters of nosocomial transmission among ST26 isolates. One cluster from two premature infants with overlapping NICU hospitalizations in 2022 encompassed 15 fluconazole-resistant isolates harboring pleiotropic drug-resistance transcription factor (Pdr1p) P258L (13 isolates) or N1086I (two isolates), together with 10 fluconazole-susceptible dose-dependent isolates lacking Pdr1p SNPs. The other cluster indicated unforeseen clonal transmission of fluconazole-resistant bloodstream isolates among five patients (four post-lung transplantation and one with diffuse interstitial lung disease) in the same hospital over 8 months. Among these five isolates, four obtained after exposure to azole antifungals harbored distinct Pdr1p SNPs (N1091D, E388Q, K365E, and R376Q). The findings reveal the transmission patterns of clonal bloodstream isolates of C. glabrata among patients undergoing antifungal treatment, exhibiting different levels of fluconazole susceptibility or distinct Pdr1p SNP profiles.

IMPORTANCE

The prevalence of fluconazole-resistant bloodstream infections caused by Candida glabrata is increasing globally, but the transmission of these resistant strains within hospitals has rarely been documented. Through whole-genome sequencing and epidemiological analyses, this study identified two potential clusters of C. glabrata bloodstream infections within the same hospital, revealing the transmission of clonal C. glabrata strains with different levels of fluconazole susceptibility or distinct transcription factor pleiotropic drug resistance protein 1 (Pdr1p) single-nucleotide polymorphism profiles among patients receiving antifungal therapy.

Temporal transcriptional response of Candida glabrata during macrophage infection reveals a multifaceted transcriptional regulator CgXbp1 important for macrophage response and fluconazole resistance

Candida glabrata can thrive inside macrophages and tolerate high levels of azole antifungals. These innate abilities render infections by this human pathogen a clinical challenge. How C. glabrata reacts inside macrophages and what is the molecular basis of its drug tolerance are not well understood. Here, we mapped genome-wide RNA polymerase II (RNAPII) occupancy in C. glabrata to delineate its transcriptional responses during macrophage infection in high temporal resolution. RNAPII profiles revealed dynamic C. glabrata responses to macrophages with genes of specialized pathways activated chronologically at different times of infection. We identified an uncharacterized transcription factor (CgXbp1) important for the chronological macrophage response, survival in macrophages, and virulence. Genome-wide mapping of CgXbp1 direct targets further revealed its multi-faceted functions, regulating not only virulence-related genes but also genes associated with drug resistance. Finally, we showed that CgXbp1 indeed also affects fluconazole resistance. Overall, this work presents a powerful approach for examining host-pathogen interaction and uncovers a novel transcription factor important for C. glabrata's survival in macrophages and drug tolerance.

Identification of Homeobox Transcription Factors in a Dimorphic Fungus Talaromyces marneffei and Protein-Protein Interaction Prediction of RfeB

Talaromyces marneffei is a thermally dimorphic fungus that can cause life-threatening systemic mycoses, particularly in immunocompromised individuals. Fungal homeobox transcription factors control various developmental processes, including the regulation of sexual reproduction, morphology, metabolism, and virulence. However, the function of homeobox proteins in T. marneffei has not been fully explored. Here, we searched the T. marneffei genome for the total homeobox transcription factors and predicted their biological relevance by performing gene expression analysis in different cell types, including conidia, mycelia, yeasts, and during phase transition. RfeB is selected for further computational analysis since (i) its transcripts were differentially expressed in different phases of T. marneffei, and (ii) this protein contains the highly conserved protein-protein interaction region (IR), which could be important for pathobiology and have therapeutic application. To assess the structure-function of the IR region, in silico alanine substitutions were performed at three-conserved IR residues (Asp276, Glu279, and Gln282) of RfeB, generating a triple RfeB mutated protein. Using 3D modeling and molecular dynamics simulations, we compared the protein complex formation of wild-type and mutated RfeB proteins with the putative partner candidate TmSwi5. Our results demonstrated that the mutated RfeB protein exhibited increased free binding energy, elevated protein compactness, and a reduced number of atomic contacts, suggesting disrupted protein stability and interaction. Notably, our model revealed that the IR residues primarily stabilized the RfeB binding sites located in the central region (CR). This computational approach for protein mutagenesis could provide a foundation for future experimental studies on the functional characterization of RfeB and other homeodomain-containing proteins in T. marneffei.

A Brg1-Rme1 circuit in Candida albicans hyphal gene regulation

Major Candida albicans virulence traits include its ability to make hyphae, to produce a biofilm, and to damage host cells. These traits depend upon expression of hypha-associated genes. A gene expression comparison among clinical isolates suggested that transcription factor Rme1, established by previous studies to be a positive regulator of chlamydospore formation, may also be a negative regulator of hypha-associated genes. Engineered RME1 overexpression supported this hypothesis, but no relevant rme1Δ/Δ mutant phenotype was detected. We reasoned that Rme1 may function within a specific regulatory pathway. This idea was supported by our finding that an rme1Δ/Δ mutation relieves the need for biofilm regulator Brg1 in biofilm formation. The impact of the rme1Δ/Δ mutation is most prominent under static or "biofilm-like" growth conditions. RNA sequencing (RNA-seq) of cells grown under biofilm-like conditions indicates that Brg1 activates hypha-associated genes indirectly via repression of RME1: hypha-associated gene expression levels are substantially reduced in a brg1Δ/Δ mutant and partially restored in a brg1Δ/Δ rme1Δ/Δ double mutant. An rme1Δ/Δ mutation does not simply bypass Brg1, because iron homeostasis genes depend upon Brg1 regardless of Rme1. Rme1 thus connects Brg1 to the targets relevant to hypha and biofilm formation under biofilm growth conditions.IMPORTANCECandida albicans is a major fungal pathogen of humans, and its ability to grow as a surface-associated biofilm on implanted devices is a common cause of infection. Here, we describe a new regulator of biofilm formation, RME1, whose activity is most prominent under biofilm-like growth conditions.

A human commensal-pathogenic fungus suppresses host immunity via targeting TBK1

Candida albicans stably colonizes humans but is the leading cause of hospital-acquired fungemia. Traditionally, masking immunogenic moieties has been viewed as a tactic for immune evasion. Here, we demonstrate that C. albicans blocks type I interferon (IFN-I) signaling via translocating an effector protein Cmi1 into host cells. Mechanistically, Cmi1 binds and inhibits TANK-binding kinase 1 (TBK1) to abrogate IFN-regulatory factor 3 (IRF3) phosphorylation, thereby suppressing the IFN-I cascade. Murine infection with a cmi1 mutant displays an exaggerated IFN-I response in both kidneys and bone-marrow-derived macrophages, leading to rapid fungal clearance and host survival. Remarkably, the lack of CMI1 compromises gut commensalism and increases IFN-I response in mouse colonic cells. These phenotypes of cmi1 are rescued by the depletion of IFN-I receptor. This work establishes the importance of TBK1 inhibition in fungal pathogenesis and reveals that a human commensal-pathogenic fungus significantly impacts host immunity during gut colonization and infection via delivering effector proteins into host cells.

Sinefungin, a natural nucleoside analog of S-adenosyl methionine, impairs the pathogenicity of Candida albicans

Candida albicans, an opportunistic fungal pathogen, causes life-threatening infections in immunocompromised patients. Current antifungals are limited by toxicity, drug-drug interactions, and emerging resistance, underscoring the importance of identifying novel treatment approaches. Here, we elucidate the impact of sinefungin, an analog of S-adenosyl methionine, on the virulence of C. albicans strain SC5314 and clinical isolates. Our data indicate that sinefungin impairs pathogenic traits of C. albicans including hyphal morphogenesis, biofilm formation, adhesion to epithelial cells, and virulence towards Galleria mellonella, highlighting sinefungin as an avenue for therapeutic intervention. We determine that sinefungin particularly disturbs N6-methyladenosine (m6A) formation. Transcriptome analysis of C. albicans hyphae upon sinefungin treatment reveals an increase in transcripts related to the yeast form and decrease in those associated with hyphae formation and virulence. Collectively, our data propose sinefungin as a potent molecule against C. albicans and emphasize further exploration of post-transcriptional control mechanisms of pathogenicity for antifungal design.

A Proteogenomic Pipeline for the Analysis of Protein Biosynthesis Errors in the Human Pathogen Candida albicans

Candida albicans is a diploid pathogen known for its ability to live as a commensal fungus in healthy individuals but causing both superficial infections and disseminated candidiasis in immunocompromised patients where it is associated with high morbidity and mortality. Its success in colonizing the human host is attributed to a wide range of virulence traits that modulate interactions between the host and the pathogen, such as optimal growth rate at 37 °C, the ability to switch between yeast and hyphal forms, and a remarkable genomic and phenotypic plasticity. A fascinating aspect of its biology is a prominent heterogeneous proteome that arises from frequent genomic rearrangements, high allelic variation, and high levels of amino acid misincorporations in proteins. This leads to increased morphological and physiological phenotypic diversity of high adaptive potential, but the scope of such protein mistranslation is poorly understood due to technical difficulties in detecting and quantifying amino acid misincorporation events in complex protein samples. We have developed and optimized mass spectrometry and bioinformatics pipelines capable of identifying rare amino acid misincorporation events at the proteome level. We have also analyzed the proteomic profile of an engineered C. albicans strain that exhibits high level of leucine misincorporation at protein CUG sites and employed an in vivo quantitative gain-of-function fluorescence reporter system to validate our LC-MS/MS data. C. albicans misincorporates amino acids above the background level at protein sites of diverse codons, particularly at CUG, confirming our previous data on the quantification of leucine incorporation at single CUG sites of recombinant reporter proteins, but increasing misincorporation of Leucine at these sites does not alter the translational fidelity of the other codons. These findings indicate that the C. albicans statistical proteome exceeds prior estimates, suggesting that its highly plastic phenome may also be modulated by environmental factors due to translational ambiguity.

Functional genomic analysis of genes important for Candida albicans fitness in diverse environmental conditions

Fungal pathogens such as Candida albicans pose a significant threat to human health with limited treatment options available. One strategy to expand the therapeutic target space is to identify genes important for pathogen growth in host-relevant environments. Here, we leverage a pooled functional genomic screening strategy to identify genes important for fitness of C. albicans in diverse conditions. We identify an essential gene with no known Saccharomyces cerevisiae homolog, C1_09670C, and demonstrate that it encodes subunit 3 of replication factor A (Rfa3). Furthermore, we apply computational analyses to identify functionally coherent gene clusters and predict gene function. Through this approach, we predict the cell-cycle-associated function of C3_06880W, a previously uncharacterized gene required for fitness specifically at elevated temperatures, and follow-up assays confirm that C3_06880W encodes Iml3, a component of the C. albicans kinetochore with roles in virulence in vivo. Overall, this work reveals insights into the vulnerabilities of C. albicans.

The spliceosome impacts morphogenesis in the human fungal pathogen Candida albicans

At human body temperature, the fungal pathogen Candida albicans can transition from yeast to filamentous morphologies in response to host-relevant cues. Additionally, elevated temperatures encountered during febrile episodes can independently induce C. albicans filamentation. However, the underlying genetic pathways governing this developmental transition in response to elevated temperatures remain largely unexplored. Here, we conducted a functional genomic screen to unravel the genetic mechanisms orchestrating C. albicans filamentation specifically in response to elevated temperature, implicating 45% of genes associated with the spliceosome or pre-mRNA splicing in this process. Employing RNA-Seq to elucidate the relationship between mRNA splicing and filamentation, we identified greater levels of intron retention in filaments compared to yeast, which correlated with reduced expression of the affected genes. Intriguingly, homozygous deletion of a gene encoding a spliceosome component important for filamentation (PRP19) caused even greater levels of intron retention compared with wild type and displayed globally dysregulated gene expression. This suggests that intron retention is a mechanism for fine-tuning gene expression during filamentation, with perturbations of the spliceosome exacerbating this process and blocking filamentation. Overall, this study unveils a novel biological process governing C. albicans filamentation, providing new insights into the complex regulation of this key virulence trait.IMPORTANCEFungal pathogens such as Candida albicans can cause serious infections with high mortality rates in immunocompromised individuals. When C. albicans is grown at temperatures encountered during human febrile episodes, yeast cells undergo a transition to filamentous cells, and this process is key to its virulence. Here, we expanded our understanding of how C. albicans undergoes filamentation in response to elevated temperature and identified many genes involved in mRNA splicing that positively regulate filamentation. Through transcriptome analyses, we found that intron retention is a mechanism for fine-tuning gene expression in filaments, and perturbation of the spliceosome exacerbates intron retention and alters gene expression substantially, causing a block in filamentation. This work adds to the growing body of knowledge on the role of introns in fungi and provides new insights into the cellular processes that regulate a key virulence trait in C. albicans.

Variations in candidalysin amino acid sequence influence toxicity and host responses

Candida albicans causes millions of mucosal infections in humans annually. Hyphal overgrowth on mucosal surfaces is frequently associated with tissue damage caused by candidalysin, a secreted peptide toxin that destabilizes the plasma membrane of host cells thereby promoting disease and immunopathology. Candidalysin was first identified in C. albicans strain SC5314, but recent investigations have revealed candidalysin "variants" of differing amino acid sequence in isolates of C. albicans, and the related species C. dubliniensis, and C tropicalis, suggesting that sequence variation among candidalysins may be widespread in natural populations of these Candida species. Here, we analyzed ECE1 gene sequences from 182 C. albicans isolates, 10 C. dubliniensis isolates, and 78 C. tropicalis isolates and identified 10, 3, and 2 candidalysin variants in these species, respectively. Application of candidalysin variants to epithelial cells revealed differences in the ability to cause cellular damage, changes in metabolic activity, calcium influx, MAPK signalling, and cytokine secretion, while biophysical analyses indicated that variants exhibited differences in their ability to interact with and permeabilize a membrane. This study identifies candidalysin variants with differences in biological activity that are present in medically relevant Candida species.

IMPORTANCE

Fungal infections are a significant burden to health. Candidalysin is a toxin produced by Candida albicans that damages host tissues, facilitating infection. Previously, we demonstrated that candidalysins exist in the related species C. dubliniensis and C. tropicalis, thereby identifying these molecules as a toxin family. Recent genomic analyses have highlighted the presence of a small number of candidalysin "variant" toxins, which have different amino acid sequences to those originally identified. Here, we screened genome sequences of isolates of C. albicans, C. dubliniensis, and C. tropicalis and identified candidalysin variants in all three species. When applied to epithelial cells, candidalysin variants differed in their ability to cause damage, activate intracellular signaling pathways, and induce innate immune responses, while biophysical analysis revealed differences in the ability of candidalysin variants to interact with lipid bilayers. These findings suggest that intraspecies variation in candidalysin amino acid sequence may influence fungal pathogenicity.

Deletion of the Candida albicans TLO gene family results in alterations in membrane sterol composition and fluconazole tolerance

Development of resistance and tolerance to antifungal drugs in Candida albicans can compromise treatment of infections caused by this pathogenic yeast species. The uniquely expanded C. albicans TLO gene family is comprised of 14 paralogous genes which encode Med2, a subunit of the multiprotein Mediator complex which is involved in the global control of transcription. This study investigates the acquisition of fluconazole tolerance in a mutant in which the entire TLO gene family has been deleted. This phenotype was reversed to varying degrees upon reintroduction of representative members of the alpha- and beta-TLO clades (i.e. TLO1 and TLO2), but not by TLO11, a gamma-clade representative. Comparative RNA sequencing analysis revealed changes in the expression of genes involved in a range of cellular functions, including ergosterol biosynthesis, mitochondrial function, and redox homeostasis. This was supported by the results of mass spectrometry analysis, which revealed alterations in sterol composition of the mutant cell membrane. Our data suggest that members of the C. albicans TLO gene family are involved in the control of ergosterol biosynthesis and mitochondrial function and may play a role in the responses of C. albicans to azole antifungal agents.

The putative prenyltransferase Nus1 is required for filamentation in the human fungal pathogen Candida albicans

Candida albicans is a major fungal pathogen of humans that can cause serious systemic infections in vulnerable immunocompromised populations. One of its virulence attributes is its capacity to transition between yeast and filamentous morphologies, but our understanding of this process remains incomplete. Here, we analyzed data from a functional genomic screen performed with the C. albicans Gene Replacement And Conditional Expression collection to identify genes crucial for morphogenesis in host-relevant conditions. Through manual scoring of microscopy images coupled with analysis of each image using a deep learning-based method termed Candescence, we identified 307 genes important for filamentation in tissue culture medium at 37°C with 5% CO2. One such factor was orf19.5963, which is predicted to encode the prenyltransferase Nus1 based on sequence homology to Saccharomyces cerevisiae. We further showed that Nus1 and its predicted interacting partner Rer2 are important for filamentation in multiple liquid filament-inducing conditions as well as for wrinkly colony formation on solid agar. Finally, we highlight that Nus1 and Rer2 likely govern C. albicans morphogenesis due to their importance in intracellular trafficking, as well as maintaining lipid homeostasis. Overall, this work identifies Nus1 and Rer2 as important regulators of C. albicans filamentation and highlights the power of functional genomic screens in advancing our understanding of gene function in human fungal pathogens.

Candida albicans natural diversity: a resource to dissect fungal commensalism and pathogenesis

Candida albicans is a ubiquitous fungus of humans. It is not only a component of the oral and intestinal microbiota of most healthy adults but also a major cause of mucosal disorders and life-threatening disseminated infections. Until recently, research on the biology and pathogenesis of the fungus was largely based on a single clinical isolate. We review investigations that have started to dissect a diverse set of C. albicans strains. Using different approaches to leverage the species' phenotypic and/or genetic diversity, these studies illuminate the wide range of interactions between fungus and host. While connecting genetic variants to phenotypes of interest remains challenging, research on C. albicans' natural diversity is central to understand fungal commensalism and pathogenesis.

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, oxidative stress resistance 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. Initial, large-scale genomic alterations like aneuploidies and loss of heterozygosity eventually resolved, as populations 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 severity of others like membrane stress responses diverged from Pi scarcity fitness. Among preliminary candidate genes for contributors to fitness recovery, those with links to TORC1 were overrepresented. Since Pi homeostasis differs substantially between fungi and humans, adaptive processes to Pi deprivation may harbor small-molecule targets that impact fungal growth, stress resistance and virulence.

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.

The zinc cluster transcription factor Znc1 regulates Rta3-dependent miltefosine resistance in Candida albicans

Zinc cluster transcription factors (ZCFs) are a family of transcription regulators that are almost exclusively found in the fungal kingdom. Activating mutations in the ZCFs Mrr1, Tac1, and Upc2 frequently cause acquired resistance to the widely used antifungal drug fluconazole in the pathogenic yeast Candida albicans. Similar to a hyperactive Tac1, a constitutively active form of the ZCF Znc1 causes increased fluconazole resistance by upregulating the multidrug efflux pump-encoding gene CDR1. Hyperactive forms of both Tac1 and Znc1 also cause overexpression of RTA3, which encodes a seven-transmembrane receptor protein involved in the regulation of asymmetric lipid distribution in the plasma membrane. RTA3 expression is also upregulated by miltefosine, an antiparasitic drug that is active against fungal pathogens and considered for treatment of invasive candidiasis, and rta3Δ mutants are hypersensitive to miltefosine. We found that activated forms of both Tac1 and Znc1 confer increased miltefosine resistance, which was dependent on RTA3 whereas CDR1 was dispensable. Intriguingly, the induction of RTA3 expression by miltefosine depended on Znc1, but not Tac1, in contrast to the known Tac1-dependent RTA3 upregulation by fluphenazine. In line with this observation, znc1Δ mutants were hypersensitive to miltefosine, whereas tac1Δ mutants showed wild-type tolerance. Forced expression of RTA3 reverted the hypersensitivity of znc1Δ mutants, demonstrating that the hypersensitivity was caused by the inability of the mutants to upregulate RTA3 in response to the drug. These findings establish Znc1 as a key regulator of miltefosine-induced RTA3 expression that is important for wild-type miltefosine tolerance.

IMPORTANCE

Transcription factors are central regulators of gene expression, and knowledge about which transcription factor regulates specific genes in response to a certain signal is important to understand the behavior of organisms. In the pathogenic yeast Candida albicans, the RTA3 gene is required for wild-type tolerance of miltefosine, an antiparasitic drug that is considered for treatment of invasive candidiasis. Activated forms of the transcription factors Tac1 and Znc1 cause constitutive overexpression of RTA3 and thereby increased miltefosine resistance, but only Tac1 mediates upregulation of RTA3 in response to the known inducer fluphenazine. RTA3 expression is also induced by miltefosine, and we found that this response depends on Znc1, whereas Tac1 is dispensable. Consequently, znc1Δ mutants were hypersensitive to miltefosine, whereas tac1Δ mutants showed wild-type tolerance. These findings demonstrate that Znc1 is the key regulator of RTA3 expression in response to miltefosine that is important for wild-type miltefosine tolerance.

Strain variation in the Candida albicans iron limitation response

Iron acquisition is critical for pathogens to proliferate during invasive infection, and the human fungal pathogen Candida albicans is no exception. The iron regulatory network, established in reference strain SC5314 and derivatives, includes the central player Sef1, a transcription factor that activates iron acquisition genes in response to iron limitation. Here, we explored potential variation in this network among five diverse C. albicans strains through mutant analysis, Nanostring gene expression profiling, and, for two strains, RNA-Seq. Our findings highlight four features that may inform future studies of natural variation and iron acquisition in this species. (i) Conformity: In all strains, major iron acquisition genes are upregulated during iron limitation, and a sef1Δ/Δ mutation impairs that response and growth during iron limitation. (ii) Response variation: Some aspects of the iron limitation response vary among strains, notably the activation of hypha-associated genes. As this gene set is tied to tissue damage and virulence, variation may impact the progression of infection. (iii) Genotype-phenotype variation: The impact of a sef1Δ/Δ mutation on cell wall integrity varies, and for the two strains examined the phenotype correlated with sef1Δ/Δ impact on several cell wall integrity genes. (iv) Phenotype discovery: DNA repair genes were induced modestly by iron limitation in sef1Δ/Δ mutants, with fold changes we would usually ignore. However, the response occurred in both strains tested and was reminiscent of a much stronger response described in Cryptococcus neoformans, a suggestion that it may have biological meaning. In fact, we observed that the iron limitation of a sef1Δ/Δ mutant caused recessive phenotypes to emerge at two heterozygous loci. Overall, our results show that a network that is critical for pathogen proliferation presents variation outside of its core functions.IMPORTANCEA key virulence factor of Candida albicans is the ability to maintain iron homeostasis in the host where iron is scarce. We focused on a central iron regulator, SEF1. We found that iron regulator Sef1 is required for growth, cell wall integrity, and genome integrity during iron limitation. The novel aspect of this work is the characterization of strain variation in a circuit that is required for survival in the host and the connection of iron acquisition to genome integrity in C. albicans.

Molecular Evaluation of the mRNA Expression of the ERG11, ERG3, CgCDR1, and CgSNQ2 Genes Linked to Fluconazole Resistance in Candida glabrata in a Colombian Population

INTRODUCTION

The study of Candida glabrata genes associated with fluconazole resistance, from a molecular perspective, increases the understanding of the phenomenon with a view to its clinical applicability.

OBJECTIVE

We sought to establish the predictive molecular profile of fluconazole resistance in Candida glabrata by analyzing the ERG11, ERG3, CgCDR1, and CgSNQ2 genes.

METHOD

Expression was quantified using RT-qPCR. Metrics were obtained through molecular docking and Fisher discriminant functions. Additionally, a predictive classification was made against the susceptibility of C. glabrata to fluconazole.

RESULTS

The relative expression of the ERG3, CgCDR1, and CgSNQ2 genes was higher in the fluconazole-resistant strains than in the fluconazole-susceptible, dose-dependent strains. The gene with the highest relative expression in the fluconazole-exposed strains was CgCDR1, and in both the resistant and susceptible, dose-dependent strains exposed to fluconazole, this was also the case. The molecular docking model generated a median number of contacts between fluconazole and ERG11 that was lower than the median number of contacts between fluconazole and ERG3, -CgCDR1, and -CgSNQ2. The predicted classification through the multivariate model for fluconazole susceptibility achieved an accuracy of 73.5%.

CONCLUSION

The resistant strains had significant expression levels of genes encoding efflux pumps and the ERG3 gene. Molecular analysis makes the identification of a low affinity between fluconazole and its pharmacological target possible, which may explain the lower intrinsic susceptibility of the fungus to fluconazole.

Long-term stability of acquired drug resistance and resistance associated mutations in the fungal pathogen Nakaseomyces glabratus (Candida glabrata)

The limited number of available antifungal drugs and the increasing number of fungal isolates that show drug or multidrug resistance pose a serious medical threat. Several yeast pathogens, such as Nakaseomyces glabratus (Candida glabrata), show a remarkable ability to develop drug resistance during treatment through the acquisition of genetic mutations. However, how stable this resistance and the underlying mutations are in non-selective conditions remains poorly characterized. The stability of acquired drug resistance has fundamental implications for our understanding of the appearance and spread of drug-resistant outbreaks and for defining efficient strategies to combat them. Here, we used an in vitro evolution approach to assess the stability under optimal growth conditions of resistance phenotypes and resistance-associated mutations that were previously acquired under exposure to antifungals. Our results reveal a remarkable stability of the resistant phenotype and the underlying mutations in a significant number of evolved populations, which conserved their phenotype for at least two months in the absence of drug-selective pressure. We observed a higher stability of anidulafungin resistance over fluconazole resistance, and of resistance-conferring point mutations as compared with aneuploidies. In addition, we detected accumulation of novel mutations in previously altered resistance-associated genes in non-selective conditions, which suggest a possible compensatory role. We conclude that acquired resistance, particularly to anidulafungin, is a long-lasting phenotype, which has important implications for the persistence and propagation of drug-resistant clinical outbreaks.

Development of a Shuttle Vector That Transforms at High Frequency for the Emerging Human Fungal Pathogen: Candida auris

Routine molecular manipulation of any organism is inefficient and difficult without the existence of a plasmid. Although transformation is possible in C. auris, no plasmids are available that can serve as cloning or shuttle vectors. C. auris centromeres have been well characterized but have not been explored further as molecular tools. We tested C. auris centromeric sequences to identify which, if any, could be used to create a plasmid that was stably maintained after transformation. We cloned all seven C. auris centromeric sequences and tested them for transformation frequency and stability. Transformation frequency varied significantly; however, one was found to transform at a very high frequency. A 1.7 Kb subclone of this sequence was used to construct a shuttle vector. The vector was stable with selection and maintained at ~1 copy per cell but could be easily lost when selection was removed, which suggested that the properties of the centromeric sequence were more Autonomously Replicating Sequence (ARS)-like than centromere-like when part of a plasmid. Rescue of this plasmid from transformed C. auris cells into E. coli revealed that it remained intact after the initial C. auris transformation, even when carrying large inserts. The plasmid was found to be able to transform all four clades of C. auris, with varying frequencies. This plasmid is an important new reagent in the C. auris molecular toolbox, which will enhance the investigation of this human fungal pathogen.

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.

Molecular Diagnostics for Invasive Fungal Diseases: Current and Future Approaches

Invasive fungal diseases (IFDs) comprise a growing healthcare burden, especially given the expanding population of immunocompromised hosts. Early diagnosis of IFDs is required to optimise therapy with antifungals, especially in the setting of rising rates of antifungal resistance. Molecular techniques including nucleic acid amplification tests and whole genome sequencing have potential to offer utility in overcoming limitations with traditional phenotypic testing. However, standardisation of methodology and interpretations of these assays is an ongoing undertaking. The utility of targeted Aspergillus detection has been well-defined, with progress in investigations into the role of targeted assays for Candida, Pneumocystis, Cryptococcus, the Mucorales and endemic mycoses. Likewise, whilst broad-range polymerase chain reaction assays have been in use for some time, pathology stewardship and optimising diagnostic yield is a continuing exercise. As costs decrease, there is also now increased access and experience with whole genome sequencing, including metagenomic sequencing, which offers unparalleled resolution especially in the investigations of potential outbreaks. However, their role in routine diagnostic use remains uncommon and standardisation of techniques and workflow are required for wider implementation.

Candida glabrata maintains two Hap1 homologs, Zcf27 and Zcf4, for distinct roles in ergosterol gene regulation to mediate sterol homeostasis under azole and hypoxic conditions

Candida glabrata exhibits innate resistance to azole antifungal drugs but also has the propensity to rapidly develop clinical drug resistance. Azole drugs, which target Erg11, is one of the three major classes of antifungals used to treat Candida infections. Despite their widespread use, the mechanism controlling azole-induced ERG gene expression and drug resistance in C. glabrata has primarily revolved around Upc2 and/or Pdr1. In this study, we determined the function of two zinc cluster transcription factors, Zcf27 and Zcf4, as direct but distinct regulators of ERG genes. Our phylogenetic analysis revealed C. glabrata Zcf27 and Zcf4 as the closest homologs to Saccharomyces cerevisiae Hap1. Hap1 is a known zinc cluster transcription factor in S. cerevisiae in controlling ERG gene expression under aerobic and hypoxic conditions. Interestingly, when we deleted HAP1 or ZCF27 in either S. cerevisiae or C. glabrata, respectively, both deletion strains showed altered susceptibility to azole drugs, whereas the strain deleted for ZCF4 did not exhibit azole susceptibility. We also determined that the increased azole susceptibility in a zcf27Δ strain is attributed to decreased azole-induced expression of ERG genes, resulting in decreased levels of total ergosterol. Surprisingly, Zcf4 protein expression is barely detected under aerobic conditions but is specifically induced under hypoxic conditions. However, under hypoxic conditions, Zcf4 but not Zcf27 was directly required for the repression of ERG genes. This study provides the first demonstration that Zcf27 and Zcf4 have evolved to serve distinct roles allowing C. glabrata to adapt to specific host and environmental conditions.

Metabolic reprogramming during Candida albicans planktonic-biofilm transition is modulated by the transcription factors Zcf15 and Zcf26

Candida albicans is a commensal of the human microbiota that can form biofilms on implanted medical devices. These biofilms are tolerant to antifungals and to the host immune system. To identify novel genes modulating C. albicans biofilm formation, we performed a large-scale screen with 2,454 C. albicans doxycycline-dependent overexpression strains and identified 16 genes whose overexpression significantly hampered biofilm formation. Among those, overexpression of the ZCF15 and ZCF26 paralogs that encode transcription factors and have orthologs only in biofilm-forming species of the Candida clade, caused impaired biofilm formation both in vitro and in vivo. Interestingly, overexpression of ZCF15 impeded biofilm formation without any defect in hyphal growth. Transcript profiling, transcription factor binding, and phenotypic microarray analyses conducted upon overexpression of ZCF15 and ZCF26 demonstrated their role in reprogramming cellular metabolism by regulating central metabolism including glyoxylate and tricarboxylic acid cycle genes. Taken together, this study has identified a new set of biofilm regulators, including ZCF15 and ZCF26, that appear to control biofilm development through their specific role in metabolic remodeling.

EpiCandIn: An open online resource for epidemiology of Candida infections in India

Background & objectives Candida spp. cause candidiasis in humans under conditions disrupting the host defence. While Candida albicans is the most reported cause of candidiasis, there is a surge in the incidence of infections by non-albicans Candida species (NACs), such as C. tropicalis, C. glabrata and C. auris. These species can infect all organs of the human body. To effectively manage these outbreaks, it is important to track the epidemiology of candidiasis. A consolidated resource describing the landscape of candidiasis in India is absent. Methods To address this gap, we have developed an online resource named Epidemiology of Candida Infections in India (EpiCandIn) by manually curating published literature on Candida infections in the Indian population obtained from PubMed and ScienceDirect databases. Results EpiCandIn contains data available since 1972 from 51 sites across 16 States and four Union Territories of India. It provides information on geographical location, Candida species, niche affected, disease characteristics and drug therapy details extracted from the publications. This resource is integrated with visualization tools. Interpretation & conclusions EpiCandIn will be useful for public health researchers and policymakers as it will help them gain insights into the emerging trends and management of Candida infections in India. It can be accessed at epicandin.bicnirrh.res.in.

Structural and Biophysical Dynamics of Fungal Plasma Membrane Proteins and Implications for Echinocandin Action in Candida glabrata

Fungal plasma membrane proteins represent key therapeutic targets for antifungal agents, yet their structure and spatial distribution in the native context remain poorly characterized. Herein, we employ an integrative multimodal approach to elucidate the structural and functional organization of plasma membrane protein complexes in Candida glabrata , focusing on prominent and essential membrane proteins, the polysaccharide synthase β-(1,3)-glucan synthase (GS) and the proton pump Pma1. Cryo-electron tomography (cryo-ET) and live cell imaging reveal that GS and Pma1 are heterogeneously distributed into distinct plasma membrane microdomains. Treatment with caspofungin, an echinocandin antifungal that targets GS, alters the plasma membrane and disrupts the native distribution of GS and Pma1. Based on these findings, we propose a model for echinocandin action that considers how drug interactions with the plasma membrane environment lead to inhibition of GS. Our work underscores the importance of interrogating the structural and dynamic characteristics of fungal plasma membrane proteins in situ to understand function and facilitate precisely targeted development of novel antifungal therapies.

In silico genome wide identification of long non-coding RNAs differentially expressed during Candida auris host pathogenesis

Candida auris is an invasive fungal pathogen of high concern due to acquired drug tolerance against antifungals used in clinics. The prolonged persistence on biotic and abiotic surfaces can result in onset of hospital outbreaks causing serious health threat. An in depth understanding of pathology of C. auris is highly desirable for development of efficient therapeutics. Non-coding RNAs play crucial role in fungal pathology. However, the information about ncRNAs is scanty to be utilized. Herein our aim is to identify long noncoding RNAs with potent role in pathobiology of C. auris. Thereby, we analyzed the transcriptomics data of C. auris infection in blood for identification of potential lncRNAs with regulatory role in determining invasion, survival or drug tolerance under infection conditions. Interestingly, we found 275 lncRNAs, out of which 253 matched with lncRNAs reported in Candidamine, corroborating for our accurate data analysis pipeline. Nevertheless, we obtained 23 novel lncRNAs not reported earlier. Three lncRNAs were found to be under expressed throughout the course of infection, in the transcriptomics data. 16 of potent lncRNAs were found to be coexpressed with coding genes, emphasizing for their functional role. Noteworthy, these ncRNAs are expressed from intergenic regions of the genes associated with transporters, metabolism, cell wall biogenesis. This study recommends for possible association between lncRNA expression and C. auris pathogenesis.