Genome Comparisons of Candida glabrata Serial Clinical Isolates Reveal Patterns of Genetic Variation in Infecting Clonal Populations

Implication of genotypes for prognosis of Candida glabrata bloodstream infections

BACKGROUND

Genotyping isolates of a specific pathogen may demonstrate unique patterns of antimicrobial resistance, virulence or outcomes. However, evidence for genotype-outcome association in Candida glabrata is scarce. We aimed to characterize the mycological and clinical relevance of genotypes on C. glabrata bloodstream infections (BSIs).

METHODS

Non-duplicated C. glabrata blood isolates from hospitalized adults were genotyped by MLST, and further clustered by the unweighted pair group method with arithmetic averages (UPGMA). A clonal complex (CC) was defined by UPGMA similarities of >90%. Antifungal susceptibility testing was performed by a colorimetric microdilution method and interpreted following CLSI criteria.

RESULTS

Of 48 blood isolates evaluated, 13 STs were identified. CC7 was the leading CC (n = 14; 29.2%), including 13 ST7. The overall fluconazole and echinocandin resistance rates were 6.6% and 0%, respectively. No specific resistance patterns were associated with CC7 or other CCs. Charlson comorbidity index (adjusted OR, 1.49; 95% CI, 1.05-3.11) was the only predictor for CC7. By multivariable Cox regression analyses, CC7 was independently associated with 28 day mortality [adjusted HR (aHR), 3.28; 95% CI, 1.31-8.23], even after considering potential interaction with neutropenia (aHR, 3.41; 95% CI, 1.23-9.42; P for interaction, 0.24) or limited to 34 patients with monomicrobial BSIs (aHR, 2.85; 95% CI, 1.15-7.08). Also, the Kaplan-Meier estimate showed greater mortality with CC7 (P = 0.003). Fluconazole resistance or echinocandin therapy had no significant impact on mortality.

CONCLUSIONS

Our data suggested comorbid patients were at risk of developing CC7 BSIs. Further, CC7 was independently associated with worse outcomes.

Candida glabrata hospital isolate from Lebanon reveals micafungin resistance associated with increased chitin and resistance to a cell-surface-disrupting agent

OBJECTIVES

This study aimed to identify the resistance mechanisms to micafungin and fluconazole in a clinical isolate of Candida glabrata.

METHODS

The isolate was whole-genome sequenced to identify amino acid changes in key proteins involved in antifungal resistance, and the isolate was further characterised by pathogenicity-related phenotypic assays that supported the sequencing results.

RESULTS

Amino acid substitutions were detected in 8 of 17 protein candidates. Many of these substitutions were novel, including in CHS3, CHS3B, and KRE5, which are involved in the development of micafungin resistance. Regarding fluconazole resistance, overexpression of efflux pumps was observed. Our isolate did not exhibit an increased virulence potential compared with the control strain; however, a significant increase in chitin content and potential to resist the cell surface disruptant sodium dodecyl sulphate was observed.

CONCLUSIONS

This clinical Candida glabrata isolate experienced a change in cell wall architecture, which correlates with the development of micafungin resistance.

Decoding the role of oxidative stress resistance and alternative carbon substrate assimilation in the mature biofilm growth mode of Candida glabrata

BACKGROUND

Biofilm formation is viewed as a vital mechanism in C. glabrata pathogenesis. Although, it plays a significant role in virulence but transcriptomic architecture and metabolic pathways governing the biofilm growth mode of C. glabrata remain elusive. The present study intended to investigate the genes implicated in biofilm growth phase of C. glabrata through global transcriptomic approach.

RESULTS

Functional analysis of Differentially expressed genes (DEGs) using gene ontology and pathways analysis revealed that upregulated genes are involved in the glyoxylate cycle, carbon-carbon lyase activity, pre-autophagosomal structure membrane and vacuolar parts whereas, down- regulated genes appear to be associated with glycolysis, ribonucleoside biosynthetic process, ribosomal and translation process in the biofilm growth condition. The RNA-Seq expression of eight selected DEGs (CgICL1, CgMLS1, CgPEP1, and CgNTH1, CgERG9, CgERG11, CgTEF3, and CgCOF1) was performed with quantitative real-time PCR (RT-qPCR). The gene expression profile of selected DEGs with RT-qPCR displayed a similar pattern of expression as observed in RNA-Seq. Phenotype screening of mutant strains generated for genes CgPCK1 and CgPEP1, showed that Cgpck1∆ failed to grow on alternative carbon substrate (Glycerol, Ethanol, Oleic acid) and similarly, Cgpep1∆ unable to grow on YPD medium supplemented with hydrogen peroxide. Our results suggest that in the absence of glucose, C. glabrata assimilate glycerol, oleic acid and generate acetyl coenzyme-A (acetyl-CoA) which is a central and connecting metabolite between catabolic and anabolic pathways (glyoxylate and gluconeogenesis) to produce glucose and fulfil energy requirements.

CONCLUSIONS

The study was executed using various approaches (transcriptomics, functional genomics and gene deletion) and it revealed that metabolic plasticity of C. glabrata (NCCPF-100,037) in biofilm stage modulates its virulence and survival ability to counter the stress and may promote its transition from commensal to opportunistic pathogen. The observations deduced from the present study along with future work on characterization of the proteins involved in this intricate process may prove to be beneficial for designing novel antifungal strategies.

Genomic description of acquired fluconazole- and echinocandin-resistance in patients with serial Candida glabrata isolates

Candida glabrata is one of the most common causes of systemic candidiasis, often resistant to antifungal medications. To describe the genomic context of emerging resistance, we conducted a retrospective analysis of 82 serially collected isolates from 33 patients from population-based candidemia surveillance in the United States. We used whole-genome sequencing to determine the genetic relationships between isolates obtained from the same patient. Phylogenetic analysis demonstrated that isolates from 29 patients were clustered by patient. The median SNPs between isolates from the same patient was 30 (range: 7-96 SNPs), while unrelated strains infected four patients. Twenty-one isolates were resistant to echinocandins, and 24 were resistant to fluconazole. All echinocandin-resistant isolates carried a mutation either in the FKS1 or FKS2 HS1 region. Of the 24 fluconazole-resistant isolates, 17 (71%) had non-synonymous polymorphisms in the PDR1 gene, which were absent in susceptible isolates. In 11 patients, a genetically related resistant isolate was collected after recovering susceptible isolates, indicating in vivo acquisition of resistance. These findings allowed us to estimate the intra-host diversity of C. glabrata and propose an upper boundary of 96 SNPs for defining genetically related isolates, which can be used to assess donor-to-host transmission, nosocomial transmission, or acquired resistance. IMPORTANCE In our study, mutations associated to azole resistance and echinocandin resistance were detected in Candida glabrata isolates using a whole-genome sequence. C. glabrata is the second most common cause of candidemia in the United States, which rapidly acquires resistance to antifungals, in vitro and in vivo.

Molecular epidemiology and antimicrobial resistance of vaginal Candida glabrata isolates in Namibia

Candida glabrata is the most common non-albicans Candida species that causes vulvovaginal candidiasis (VVC). Given the intrinsically low susceptibility of C. glabrata to azole drugs, investigations into C. glabrata prevalence, fungal susceptibility profile, and molecular epidemiology are necessary to optimise the treatment of VVC. This molecular epidemiological study was conducted to determine antifungal drug profile, single nucleotide polymorphisms (SNPs) associated with phenotypic antifungal resistance and epidemic diversity of C. glabrata isolates from women with VVC in Namibia. Candida glabrata isolates were identified using phenotypic and molecular methods. Antifungal susceptibility of strains was determined for fluconazole, itraconazole, amphotericin B, and anidulafungin. Whole genome sequencing was used to determine SNPs in antifungal resistance genes and sequence type (ST) allocation. Among C. glabrata isolates, all (20/20; 100%) exhibited phenotypic resistance to the azole class antifungal drug, (fluconazole), and phenotypic susceptibility to the polyene class (amphotericin B), and the echinocandins (anidulafungin). Non-synonymous SNPs were identified in antifungal resistance genes of all fluconazole-resistant C. glabrata isolates including ERG6 (15%), ERG7 (15%), CgCDR1 (25%), CgPDR1 (60%), SNQ2 (10%), FKS1 (5.0%), FKS2 (5.0%), CgFPS1 (5.0%), and MSH2 (15%). ST15 (n = 8/20, 40%) was predominant. This study provides important insight into phenotypic and genotypic antifungal resistance across C. glabrata isolates from women with VVC in Namibia. In this study, azole resistance is determined by an extensive range of SNPs, while the observed polyene and echinocandin resistance-associated SNPs despite phenotypic susceptibility require further investigation.

Overlooked Candida glabrata petites are echinocandin tolerant, induce host inflammatory responses, and display poor in vivo fitness

IMPORTANCE

Candida glabrata is a major fungal pathogen, which is able to lose mitochondria and form small and slow-growing colonies, called "petite." This attenuated growth rate has created controversies and questioned the clinical importance of petiteness. Herein, we have employed multiple omics technologies and in vivo mouse models to critically assess the clinical importance of petite phenotype. Our WGS identifies multiple genes potentially underpinning petite phenotype. Interestingly, petite C. glabrata cells engulfed by macrophages are dormant and, therefore, are not killed by the frontline antifungal drugs. Interestingly, macrophages infected with petite cells mount distinct transcriptomic responses. Consistent with our ex vivo observations, mitochondrial-proficient parental strains outcompete petites during systemic and gut colonization. Retrospective examination of C. glabrata isolates identified petite prevalence a rare entity, which can significantly vary from country to country. Collectively, our study overcomes the existing controversies and provides novel insights regarding the clinical relevance of petite C. glabrata isolates.

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.

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

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

Overlooked Candida glabrata petites are echinocandin tolerant, induce host inflammatory responses, and display poor in vivo fitness

Small colony variants (SCVs) are relatively common among some bacterial species and are associated with poor prognosis and recalcitrant infections. Similarly, Candida glabrata - a major intracellular fungal pathogen - produces small and slow-growing respiratory-deficient colonies, termed "petite." Despite reports of clinical petite C . glabrata strains, our understanding of petite behavior in the host remains obscure. Moreover, controversies exist regarding in-host petite fitness and its clinical relevance. Herein, we employed whole-genome sequencing (WGS), dual-RNAseq, and extensive ex vivo and in vivo studies to fill this knowledge gap. WGS identified multiple petite-specific mutations in nuclear and mitochondrially-encoded genes. Consistent with dual-RNAseq data, petite C . glabrata cells did not replicate inside host macrophages and were outcompeted by their non-petite parents in macrophages and in gut colonization and systemic infection mouse models. The intracellular petites showed hallmarks of drug tolerance and were relatively insensitive to the fungicidal activity of echinocandin drugs. Petite-infected macrophages exhibited a pro-inflammatory and type I IFN-skewed transcriptional program. Interrogation of international C . glabrata blood isolates ( n =1000) showed that petite prevalence varies by country, albeit at an overall low prevalence (0-3.5%). Collectively, our study sheds new light on the genetic basis, drug susceptibility, clinical prevalence, and host-pathogen responses of a clinically overlooked phenotype in a major fungal pathogen.

Importance

Candida glabrata is a major fungal pathogen, which is able to lose mitochondria and form small and slow-growing colonies, called "petite". This attenuated growth rate has created controversies and questioned the clinical importance of petiteness. Herein, we have employed multiple omicstechnologies and in vivo mouse models to critically assess the clinical importance of petite phenotype. Our WGS identifies multiple genes potentially underpinning petite phenotype. Interestingly, petite C. glabrata cells engulfed by macrophages are dormant and therefore are not killed by the frontline antifungal drugs. Interestingly, macrophages infected with petite cells mount distinct transcriptomic responses. Consistent with our ex-vivo observations, mitochondrial-proficient parental strains outcompete petites during systemic and gut colonization. Retrospective examination of C. glabrata isolates identified petite prevalence a rare entity, can significantly vary from country to country. Collectively, our study overcomes the existing controversies and provides novel insights regarding the clinical relevance of petite C. glabrata isolates.

Whole-Genome Sequence Analysis of Candida glabrata Isolates from a Patient with Persistent Fungemia and Determination of the Molecular Mechanisms of Multidrug Resistance

Whole-genome sequencing (WGS) was used to determine the molecular mechanisms of multidrug resistance for 10 serial Candida glabrata bloodstream isolates obtained from a neutropenic patient during 82 days of amphotericin B (AMB) or echinocandin therapy. For WGS, a library was prepared and sequenced using a Nextera DNA Flex Kit (Illumina) and the MiseqDx (Illumina) instrument. All isolates harbored the same Msh2p substitution, V239L, associated with multilocus sequence type 7 and a Pdr1p substitution, L825P, that caused azole resistance. Of six isolates with increased AMB MICs (≥2 mg/L), three harboring the Erg6p A158fs mutation had AMB MICs ≥ 8 mg/L, and three harboring the Erg6p R314K, Erg3p G236D, or Erg3p F226fs mutation had AMB MICs of 2-3 mg/L. Four isolates harboring the Erg6p A158fs or R314K mutation had fluconazole MICs of 4-8 mg/L while the remaining six had fluconazole MICs ≥ 256 mg/L. Two isolates with micafungin MICs > 8 mg/L harbored Fks2p (I661_L662insF) and Fks1p (C499fs) mutations, while six isolates with micafungin MICs of 0.25-2 mg/L harbored an Fks2p K1357E substitution. Using WGS, we detected novel mechanisms of AMB and echinocandin resistance; we explored mechanisms that may explain the complex relationship between AMB and azole resistance.

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.

Macrophage internalization creates a multidrug-tolerant fungal persister reservoir and facilitates the emergence of drug resistance

Candida glabrata is a major fungal pathogen notable for causing recalcitrant infections, rapid emergence of drug-resistant strains, and its ability to survive and proliferate within macrophages. Resembling bacterial persisters, a subset of genetically drug-susceptible C. glabrata cells can survive lethal exposure to the fungicidal echinocandin drugs. Herein, we show that macrophage internalization induces cidal drug tolerance in C. glabrata, expanding the persister reservoir from which echinocandin-resistant mutants emerge. We show that this drug tolerance is associated with non-proliferation and is triggered by macrophage-induced oxidative stress, and that deletion of genes involved in reactive oxygen species detoxification significantly increases the emergence of echinocandin-resistant mutants. Finally, we show that the fungicidal drug amphotericin B can kill intracellular C. glabrata echinocandin persisters, reducing emergence of resistance. Our study supports the hypothesis that intra-macrophage C. glabrata is a reservoir of recalcitrant/drug-resistant infections, and that drug alternating strategies can be developed to eliminate this reservoir.

The Candida glabrata Parent Strain Trap: How Phenotypic Diversity Affects Metabolic Fitness and Host Interactions

Reference strains improve reproducibility by standardizing observations and methodology, which has ultimately led to important insights into fungal pathogenesis. However, recent investigations have highlighted significant genotypic and phenotypic heterogeneity across isolates that influence genetic circuitry and virulence within a species. Candida glabrata is the second leading cause of candidiasis, a life-threatening infection, and undergoes extensive karyotype and phenotypic changes in response to stress. Much of the work conducted on this pathogen has focused on two sequenced strains, CBS138 (ATCC 2001) and BG2. Few studies have compared these strains in detail, but key differences include mating type and altered patterns of expression of EPA adhesins. In fact, most C. glabrata isolates and BG2 are MATa, while CBS138 is MATα. However, it is not known if other phenotypic differences between these strains play a role in our understanding of C. glabrata pathogenesis. Thus, we set out to characterize metabolic, cell wall, and host-interaction attributes for CBS138 and BG2. We found that BG2 utilized a broader range of nitrogen sources and had reduced cell wall size and carbohydrate exposure than CBS138, which we hypothesized results in differences in innate immune interactions and virulence. We observed that, although both strains were phagocytosed to a similar extent, BG2 replicated to higher numbers in macrophages and was more virulent during Galleria mellonella infection than CBS138 in a dose-dependent manner. Interestingly, deletion of SNF3, a major nutrient sensor, did not affect virulence in G. mellonella for BG2, but significantly enhanced larval killing in the CBS138 background compared to the parent strain. Understanding these fundamental differences in metabolism and host interactions will allow more robust conclusions to be drawn in future studies of C. glabrata pathogenesis. IMPORTANCE Reference strains provide essential insights into the mechanisms underlying virulence in fungal pathogens. However, recent studies in Candida albicans and other species have revealed significant genotypic and phenotypic diversity within clinical isolates that are challenging paradigms regarding key virulence factors and their regulation. Candida glabrata is the second leading cause of candidiasis, and many studies use BG2 or CBS138 for their investigations. Therefore, we aimed to characterize important virulence-related phenotypes for both strains that might alter conclusions about C. glabrata pathogenesis. Our study provides context for metabolic and cell wall changes and how these may influence host interaction phenotypes. Understanding these differences is necessary to support robust conclusions about how virulence factors may function in these and other very different strain backgrounds.

Disclosing azole resistance mechanisms in resistant Candida glabrata strains encoding wild-type or gain-of-function CgPDR1 alleles through comparative genomics and transcriptomics

The pathogenic yeast Candida glabrata is intrinsically resilient to azoles and rapidly acquires resistance to these antifungals, in vitro and in vivo. In most cases azole-resistant C. glabrata clinical strains encode hyperactive CgPdr1 variants, however, resistant strains encoding wild-type CgPDR1 alleles have also been isolated, although remaining to be disclosed the underlying resistance mechanism. In this study, we scrutinized the mechanisms underlying resistance to azoles of 8 resistant clinical C. glabrata strains, identified along the course of epidemiological surveys undertaken in Portugal. Seven of the strains were found to encode CgPdr1 gain-of-function variants (I392M, E555K, G558C, and I803T) with the substitutions I392M and I803T being herein characterized as hyper-activating mutations for the first time. While cells expressing the wild-type CgPDR1 allele required the mediator subunit Gal11A to enhance tolerance to fluconazole, this was dispensable for cells expressing the I803T variant indicating that the CgPdr1 interactome is shaped by different gain-of-function substitutions. Genomic and transcriptomic profiling of the sole azole-resistant C. glabrata isolate encoding a wild-type CgPDR1 allele (ISTB218) revealed that under fluconazole stress this strain over-expresses various genes described to provide protection against this antifungal, while also showing reduced expression of genes described to increase sensitivity to these drugs. The overall role in driving the azole-resistance phenotype of the ISTB218 C. glabrata isolate played by these changes in the transcriptome and genome of the ISTB218 isolate are discussed shedding light into mechanisms of resistance that go beyond the CgPdr1-signalling pathway and that may alone, or in combination, pave the way for the acquisition of resistance to azoles in vivo.

Using genomics to understand the mechanisms of virulence and drug resistance in fungal pathogens

Fungal pathogens pose an increasingly worrying threat to human health, food security and ecosystem diversity. To tackle fungal infections and improve current diagnostic and therapeutic tools it is necessary to understand virulence and antifungal drug resistance mechanisms in diverse species. Recent advances in genomics approaches have provided a suitable framework to understand these phenotypes, which ultimately depend on genetically encoded determinants. In this work, we review how the study of genome sequences has been key to ascertain the bases of virulence and drug resistance traits. We focus on the contribution of comparative genomics, population genomics and directed evolution studies. In addition, we discuss how different types of genomic mutations (small or structural variants) contribute to intraspecific differences in virulence or drug resistance. Finally, we review current challenges in the field and anticipate future directions to solve them. In summary, this work provides a short overview of how genomics can be used to understand virulence and drug resistance in fungal pathogens.

Population genetics and microevolution of clinical Candida glabrata reveals recombinant sequence types and hyper-variation within mitochondrial genomes, virulence genes, and drug targets

Candida glabrata is the second most common etiological cause of worldwide systemic candidiasis in adult patients. Genome analysis of 68 isolates from 8 hospitals across Scotland, together with 83 global isolates, revealed insights into the population genetics and evolution of C. glabrata. Clinical isolates of C. glabrata from across Scotland are highly genetically diverse, including at least 19 separate sequence types that have been recovered previously in globally diverse locations, and 1 newly discovered sequence type. Several sequence types had evidence for ancestral recombination, suggesting transmission between distinct geographical regions has coincided with genetic exchange arising in new clades. Three isolates were missing MATα1, potentially representing a second mating type. Signatures of positive selection were identified in every sequence type including enrichment for epithelial adhesins thought to facilitate fungal adhesin to human epithelial cells. In patent microevolution was identified from 7 sets of recurrent cases of candidiasis, revealing an enrichment for nonsynonymous and frameshift indels in cell surface proteins. Microevolution within patients also affected epithelial adhesins genes, and several genes involved in drug resistance including the ergosterol synthesis gene ERG4 and the echinocandin target FKS1/2, the latter coinciding with a marked drop in fluconazole minimum inhibitory concentration. In addition to nuclear genome diversity, the C. glabrata mitochondrial genome was particularly diverse, with reduced conserved sequence and conserved protein-encoding genes in all nonreference ST15 isolates. Together, this study highlights the genetic diversity within the C. glabrata population that may impact virulence and drug resistance, and 2 major mechanisms generating this diversity: microevolution and genetic exchange/recombination.

Functional variability in adhesion and flocculation of yeast megasatellite genes

Megasatellites are large tandem repeats found in all fungal genomes but especially abundant in the opportunistic pathogen Candida glabrata. They are encoded in genes involved in cell-cell interactions, either between yeasts or between yeast and human cells. In the present work, we have been using an iterative genetic system to delete several Candida glabrata megasatellite-containing genes and found that 2 of them were positively involved in adhesion to epithelial cells, whereas 3 genes negatively controlled adhesion. Two of the latter, CAGL0B05061g or CAGL0A04851g, were also negative regulators of yeast-to-yeast adhesion, making them central players in controlling Candida glabrata adherence properties. Using a series of synthetic Saccharomyces cerevisiae strains in which the FLO1 megasatellite was replaced by other tandem repeats of similar length but different sequences, we showed that the capacity of a strain to flocculate in liquid culture was unrelated to its capacity to adhere to epithelial cells or to invade agar. Finally, to understand how megasatellites were initially created and subsequently expanded, an experimental evolution system was set up, in which modified yeast strains containing different megasatellite seeds were grown in bioreactors for more than 200 generations and selected for their ability to sediment at the bottom of the culture tube. Several flocculation-positive mutants were isolated. Functionally relevant mutations included general transcription factors as well as a 230-kbp segmental duplication.

The effect of antifungal resistance development on the virulence of Candida species

In recent years, the relevance of diseases associated with fungal pathogens increased worldwide. Members of the Candida genus are responsible for the greatest number of fungal bloodstream infections every year. Epidemiological data consistently indicate a modest shift toward non-albicans species, albeit Candidaalbicans is still the most recognizable species within the genus. As a result, the number of clinically relevant pathogens has increased, and, despite their distinct pathogenicity features, the applicable antifungal agents remained the same. For bloodstream infections, only three classes of drugs are routinely used, namely polyenes, azoles and echinocandins. Antifungal resistance toward all three antifungal drug classes frequently occurs in clinical settings. Compared with the broad range of literature on virulence and antifungal resistance of Candida species separately, only a small portion of studies examined the effect of resistance on virulence. These studies found that resistance to polyenes and echinocandins concluded in significant decrease in the virulence in different Candida species. Meanwhile, in some cases, resistance to azole type antifungals resulted in increased virulence depending on the species and isolates. These findings underline the importance of studies aiming to dissect the connections of virulence and resistance in Candida species.

Antifungal drug resistance: Deciphering the mechanisms governing multidrug resistance in the fungal pathogen Candida glabrata

Fungal infections pose a significant health burden. In a new study, Candida glabrata isolates resistant to mechanistically distinct antifungals were evolved, unveiling mutations in ERG3 as a novel mechanism by which this pathogen evolves multidrug resistance.

Narrow mutational signatures drive acquisition of multidrug resistance in the fungal pathogen Candida glabrata

Fungal infections are a growing medical concern, in part due to increased resistance to one or multiple antifungal drugs. However, the evolutionary processes underpinning the acquisition of antifungal drug resistance are poorly understood. Here, we used experimental microevolution to study the adaptation of the yeast pathogen Candida glabrata to fluconazole and anidulafungin, two widely used antifungal drugs with different modes of action. Our results show widespread ability of rapid adaptation to one or both drugs. Resistance, including multidrug resistance, is often acquired at moderate fitness costs and mediated by mutations in a limited set of genes that are recurrently and specifically mutated in strains adapted to each of the drugs. Importantly, we uncover a dual role of ERG3 mutations in resistance to anidulafungin and cross-resistance to fluconazole in a subset of anidulafungin-adapted strains. Our results shed light on the mutational paths leading to resistance and cross-resistance to antifungal drugs.

Multicenter evaluation of attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy-based method for rapid identification of clinically relevant yeasts

Fourier transform infrared (FTIR) spectroscopy has demonstrated applicability as a reagent-free whole-organism fingerprinting technique for both microbial identification and strain typing. For routine application of this technique in microbiology laboratories, acquisition of FTIR spectra in the attenuated total reflectance (ATR) mode simplifies the FTIR spectroscopy workflow, providing results within minutes after initial culture without prior sample preparation. In our previous central work, 99.7% correct species identification of clinically relevant yeasts was achieved by employing an ATR-FTIR-based method and spectral database developed by our group. In this study, ATR-FTIR spectrometers were placed in 6 clinical microbiology laboratories over a 16-month period and were used to collect spectra of routine yeast isolates for on-site identification to the species level. The identification results were compared to those obtained from conventional biochemical tests and/or matrix-assisted laser desorption/ionization time of flight mass spectrometry. Isolates producing discordant results were reanalyzed by routine identification methods, ATR-FTIR spectroscopy and PCR gene sequencing of the D1/D2 and ITS regions. Among the 573 routine clinical yeast isolates collected and identified by the ATR-FTIR-based method, 564 isolates (98.4%) were correctly identified at the species level while the remaining isolates were inconclusive with no misidentifications. Due to the low prevalence of Candida auris in routine isolates, additional randomly selected C. auris (n = 24) isolates were obtained for evaluation and resulted in 100% correct identification. Overall, the data obtained in our multicenter evaluation study using multiple spectrometers and system operators indicate that ATR-FTIR spectroscopy is a reliable, cost-effective yeast identification technique that provides accurate and timely (∼3 minutes/sample) species identification promptly after the initial culture.

Comparative genomic analysis of clinical Candida glabrata isolates identifies multiple polymorphic loci that can improve existing multilocus sequence typing strategy

Candida glabrata is the second leading cause of candidemia in many countries and is one of the most concerning yeast species of nosocomial importance due to its increasing rate of antifungal drug resistance and emerging multidrug-resistant isolates. Application of multilocus sequence typing (MLST) to clinical C. glabrata isolates revealed an association of certain sequence types (STs) with drug resistance and mortality. The current C. glabrata MLST scheme is based on single nucleotide polymorphisms (SNPs) at six loci and is therefore relatively laborious and costly. Furthermore, only a few high-quality C. glabrata reference genomes are available, limiting rapid analysis of clinical isolates by whole genome sequencing. In this study we provide long-read based assemblies for seven additional clinical strains belonging to three different STs and use this information to simplify the C. glabrata MLST scheme. Specifically, a comparison of these genomes identified highly polymorphic loci (HPL) defined by frequent insertions and deletions (indels), two of which proved to be highly resolutive for ST. When challenged with 53 additional isolates, a combination of TRP1 (a component of the current MLST scheme) with either of the two HPL fully recapitulated ST identification. Therefore, our comparative genomic analysis identified a new typing approach combining SNPs and indels and based on only two loci, thus significantly simplifying ST identification in C. glabrata. Because typing tools are instrumental in addressing numerous clinical and biological questions, our new MLST scheme can be used for high throughput typing of C. glabrata in clinical and research settings.

Genome and transcriptome of a pathogenic yeast, Candida nivariensis

We present a highly contiguous genome and transcriptome of the pathogenic yeast, Candida nivariensis. We sequenced both the DNA and RNA of this species using both the Oxford Nanopore Technologies and Illumina platforms. We assembled the genome into an 11.8 Mb draft composed of 16 contigs with an N50 of 886 Kb, including a circular mitochondrial sequence of 28 Kb. Using direct RNA nanopore sequencing and Illumina cDNA sequencing, we constructed an annotation of our new assembly, supplemented by lifting over genes from Saccharomyces cerevisiae and Candida glabrata.

Diversity of GPI-anchored fungal adhesins

Selective adhesion of fungal cells to one another and to foreign surfaces is fundamental for the development of multicellular growth forms and the successful colonization of substrates and host organisms. Accordingly, fungi possess diverse cell wall-associated adhesins, mostly large glycoproteins, which present N-terminal adhesion domains at the cell surface for ligand recognition and binding. In order to function as robust adhesins, these glycoproteins must be covalently linkedto the cell wall via C-terminal glycosylphosphatidylinositol (GPI) anchors by transglycosylation. In this review, we summarize the current knowledge on the structural and functional diversity of so far characterized protein families of adhesion domains and set it into a broad context by an in-depth bioinformatics analysis using sequence similarity networks. In addition, we discuss possible mechanisms for the membrane-to-cell wall transfer of fungal adhesins by membrane-anchored Dfg5 transglycosidases.

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Evolutionary Divergence in DNA Damage Responses among Fungi

Cell cycle checkpoints and DNA repair pathways contribute to maintaining genome integrity and are thought to be evolutionarily ancient and broadly conserved. For example, in the yeast Saccharomyces cerevisiae and humans, DNA damage induces activation of a checkpoint effector kinase, Rad53p (human homolog Chk2), to promote cell cycle arrest and transcription of DNA repair genes.

Cell cycle checkpoints and DNA repair pathways contribute to maintaining genome integrity and are thought to be evolutionarily ancient and broadly conserved. For example, in the yeast Saccharomyces cerevisiae and humans, DNA damage induces activation of a checkpoint effector kinase, Rad53p (human homolog Chk2), to promote cell cycle arrest and transcription of DNA repair genes. However, recent studies have revealed variation in the DNA damage response networks of some fungi. For example, Shor et al. (mBio 11:e03044-20, 2020, https://doi.org/10.1128/mBio.03044-20) demonstrate that in comparison to S. cerevisiae, the fungal pathogen Candida glabrata has reduced activation of Rad53p in response to DNA damage. Consequently, some downstream targets that contribute to S. cerevisiae genome maintenance, such as DNA polymerases, are transcriptionally downregulated in C. glabrata. Downregulation of genome maintenance genes likely contributes to higher rates of mitotic failure and cell death in C. glabrata. This and other recent findings highlight evolutionary diversity in eukaryotic DNA damage responses.

Clonal evolution of Candida albicans, Candida glabrata and Candida dubliniensis at oral niche level in health and disease

Background:Candida species have long been recognised as aetiological agents of opportunistic infections of the oral mucosa, and more recently, as players of polymicrobial interactions driving caries, periodontitis and oral carcinogenesis.

Methods: We studied the clonal structure of Candida spp. at oral niche resolution in patients (n = 20) with a range of oral health profiles over 22 months. Colonies from oral micro-environments were examined with multilocus sequencing typing.

Results:Candida spp. identified were C. albicans, C. glabrata and C. dubliniensis. Increased propensity for micro-variations giving rise to multiple diploid strain types (DST), as a result of loss of heterozygosity, was observed among C. albicans clade 1 isolates compared to other clades. Micro-variations among isolates were also observed in C. dubliniensis contra to expectations of stable population structures for this species. Multiple sequence types were retrieved from patients without clinical evidence of oral candidosis, while single sequence types were isolated from oral candidosis patients.

Conclusion: This is the first study to describe the clonal population structure, persistence and stability of Candida spp. at oral niche level. Future research investigating links between Candida spp. clonality and oral disease should recognise the propensity to micro-variations amongst oral niches in C. albicans and C. dubliniensis identified here.

Human Pathogenic Candida Species Respond Distinctively to Lactic Acid Stress

Several Candida species are opportunistic human fungal pathogens and thrive in various environmental niches in and on the human body. In this study we focus on the conditions of the vaginal tract, which is acidic, hypoxic, glucose-deprived, and contains lactic acid. We quantitatively analyze the lactic acid tolerance in glucose-rich and glucose-deprived environment of five Candida species: Candidaalbicans, Candida glabrata, Candida parapsilosis, Candida krusei and Candida tropicalis. To characterize the phenotypic space, we analyzed 40–100 clinical isolates of each species. Each Candida species had a very distinct response pattern to lactic acid stress and characteristic phenotypic variability. C. glabrata and C. parapsilosis were best to withstand high concentrations of lactic acid with glucose as carbon source. A glucose-deprived environment induced lactic acid stress tolerance in all species. With lactate as carbon source the growth rate of C. krusei is even higher compared to glucose, whereas the other species grow slower. C. krusei may use lactic acid as carbon source in the vaginal tract. Stress resistance variability was highest among C. parapsilosis strains. In conclusion, each Candida spp. is adapted differently to cope with lactic acid stress and resistant to physiological concentrations.

The Quiet and Underappreciated Rise of Drug-Resistant Invasive Fungal Pathogens

Human fungal pathogens are attributable to a significant economic burden and mortality worldwide. Antifungal treatments, although limited in number, play a pivotal role in decreasing mortality and morbidities posed by invasive fungal infections (IFIs). However, the recent emergence of multidrug-resistant Candida auris and Candida glabrata and acquiring invasive infections due to azole-resistant C. parapsilosis, C. tropicalis, and Aspergillus spp. in azole-naïve patients pose a serious health threat considering the limited number of systemic antifungals available to treat IFIs. Although advancing for major fungal pathogens, the understanding of fungal attributes contributing to antifungal resistance is just emerging for several clinically important MDR fungal pathogens. Further complicating the matter are the distinct differences in antifungal resistance mechanisms among various fungal species in which one or more mechanisms may contribute to the resistance phenotype. In this review, we attempt to summarize the burden of antifungal resistance for selected non-albicansCandida and clinically important Aspergillus species together with their phylogenetic placement on the tree of life. Moreover, we highlight the different molecular mechanisms between antifungal tolerance and resistance, and comprehensively discuss the molecular mechanisms of antifungal resistance in a species level.

Candida and Candidiasis-Opportunism Versus Pathogenicity: A Review of the Virulence Traits

One of the most important questions in microbiology nowadays, is how apparently harmless, commensal yeasts like Candida spp. can cause a rising number of infections. The occurrence of the disease requires firstly the attachment to the host cells, followed by the invasion of the tissue. The adaptability translates into a rapid ability to respond to stress factors, to take up nutrients or to multiply under different conditions. By forming complex intracellular networks such as biofilms, Candida spp. become not only more refractive to antifungal therapies but also more prone to cause disease. The inter-microbial interactions can enhance the virulence of a strain. In vivo, the fungal cells face a multitude of challenges and, as a result, they develop complex strategies serving one ultimate goal: survival. This review presents the virulence factors of the most important Candida spp., contributing to a better understanding of the onset of candidiasis and raising awareness of the highly complex interspecies interactions that can change the outcome of the disease.

Evaluation of Molecular Epidemiology, Clinical Characteristics, Antifungal Susceptibility Profiles, and Molecular Mechanisms of Antifungal Resistance of Iranian Candida parapsilosis Species Complex Blood Isolates

Clonal expansion of fluconazole resistant (FLZ-R) Candida parapsilosis isolates is increasingly being identified in many countries, while there is no study exploring the antifungal susceptibility pattern, genetic diversity, and clinical information for Iranian C. parapsilosis blood isolates. Candida parapsilosis species complex blood isolates (n = 98) were recovered from nine hospitals located in three major cities, identified by MALDI-TOF MS, and their genetic relatedness was examined by AFLP fingerprinting. Antifungal susceptibility testing followed CLSI-M27-A3 and ERG11, MRR1 and hotspots 1/2 (HS1/2) of FKS1 were sequenced to assess the azole and echinocandin resistance mechanisms, respectively. Ninety-four C. parapsilosis and four Candida orthopsilosis isolates were identified from 90 patients. Only 43 patients received systemic antifungal drugs with fluconazole as the main antifungal used. The overall mortality rate was 46.6% (42/90) and death mostly occurred for those receiving systemic antifungals (25/43) relative to those not treated (17/47). Although, antifungal-resistance was rare, one isolate was multidrug-resistant (FLZ = 16 μg/ml and micafungin = 8 μg/ml) and the infected patient showed therapeutic failure to FLZ prophylaxis. Mutations causing azole and echinocandin resistance were not found in the genes studied. AFLP revealed five genotypes (G) and G1 was the main one (59/94; 62.7%). Clinical outcome was significantly associated with city (P = 0.02, α <0.05) and Mashhad was significantly associated with mortality (P = 0.03, α <0.05). Overall, we found a low level of antifungal resistance for Iranian C. parapsilosis blood isolates, but the noted MDR strain can potentially become the source of future infections and challenge the antifungal therapy in antifungal-naïve patients. AFLP typing results warrants confirmation using other resolutive typing methods.

Functional networks of co-expressed genes to explore iron homeostasis processes in the pathogenic yeast Candida glabrata

Candida glabrata is a cause of life-threatening invasive infections especially in elderly and immunocompromised patients. Part of human digestive and urogenital microbiota, C. glabrata faces varying iron availability, low during infection or high in digestive and urogenital tracts. To maintain its homeostasis, C. glabrata must get enough iron for essential cellular processes and resist toxic iron excess. The response of this pathogen to both depletion and lethal excess of iron at 30°C have been described in the literature using different strains and iron sources. However, adaptation to iron variations at 37°C, the human body temperature and to gentle overload, is poorly known. In this study, we performed transcriptomic experiments at 30°C and 37°C with low and high but sub-lethal ferrous concentrations. We identified iron responsive genes and clarified the potential effect of temperature on iron homeostasis. Our exploration of the datasets was facilitated by the inference of functional networks of co-expressed genes, which can be accessed through a web interface. Relying on stringent selection and independently of existing knowledge, we characterized a list of 214 genes as key elements of C. glabrata iron homeostasis and interesting candidates for medical applications.

The role of Lactobacillus species in the control of Candida via biotrophic interactions

Microbial communities have an important role in health and disease. Candida spp. are ubiquitous commensals and sometimes opportunistic fungal pathogens of humans, colonizing mucosal surfaces of the genital, urinary, respiratory and gastrointestinal tracts and the oral cavity. They mainly cause local mucosal infections in immune competent individuals. However, in the case of an ineffective immune defense, Candida infections may become a serious threat. Lactobacillus spp. are part of the human microbiome and are natural competitors of Candida in the vaginal environment. Lactic acid, low pH and other secreted metabolites are environmental signals sensed by fungal species present in the microbiome. This review briefly discusses the ternary interaction between host, Lactobacillus species and Candida with regard to fungal infections and the potential antifungal and fungistatic effect of Lactobacillus species. Our understanding of these interactions is incomplete due to the variability of the involved species and isolates and the complexity of the human host.

Molecular Typing of Candida glabrata

The yeast Candida glabrata has emerged, second only to Candida albicans, to be one of the most frequently isolated fungi in clinical specimen from human. Its frequent resistance towards azole antifungal drugs and the high capacity to form biofilms on indwelling catheters of individual isolates render it an often difficult to treat pathogen. Hence, there is a notably increasing scientific and clinical interest in this species. This has led to the development of a variety of molecular tools for genetic modification, strain collections, and last but not least different approaches to analyse the population structure among isolates of different geographical and clinical contexts. Often, these are used to study correlations (or the absence thereof) with different pathogenicity, virulence, or drug resistance traits. Three molecular methods have been used to type within the C. glabrata population on a genetic level by multiple studies: multi-locus sequence typing, microsatellite length polymorphisms, and clustering of whole-genome sequencing data, and these are subject of this review.

Clinical and Laboratory Development of Echinocandin Resistance in Candida glabrata: Molecular Characterization

The pathogenic yeast Candida glabrata has become a public health issue due to the increasing number of echinocandin resistant clinical strains reported. In this study, acquisition and development of resistance to this antifungal class were studied in serial C. glabrata isolates from five patients admitted in two Spanish hospitals with a resistant profile against echinocandins associated with different mutations in hot-spot 1 of FKS2 gene. For two of these patients susceptible FKS wild-type isolates obtained prior to resistant ones were also investigated. Isolates were genotyped using multilocus sequence typing and microsatellite length polymorphism techniques, which yielded comparable results. Susceptible and resistant isolates from the same patient had the same genotype, being sequence type (ST) 3 the most prevalent among them. Isolates with different FKS mutations but the same ST were present in the same patient. MSH2 gene alterations were also studied to investigate their correlation with antifungal resistance acquisition but no association was found with antifungal resistance nor with specific genotypes. In vitro exposure to increasing concentrations of micafungin to susceptible isolates developed colonies carrying FKS mutations in agar plates containing a minimum concentration of 0.06 mg/L of micafungin after less than 48 h of exposure. We investigated the correlation between development of resistance and genotype in a set of susceptible strains after being in vitro exposed to micafungin and anidulafungin but no correlation was found. Mutant prevention concentration values and spontaneous growth frequencies after selection with both echinocandins were statistically similar, although FKS mutant colonies were more abundant after micafungin exposure (p < 0.001). Mutation S663P and F659 deletion were the most common ones found after selection with both echinocandins.