Understand the genomic diversity and evolution of fungal pathogen Candida glabrata by genome-wide analysis of genetic variations

Genetic interaction analysis of Candida glabrata transcription factors CST6 and UPC2A in the regulation of respiration and fluconazole susceptibility

Candida glabrata is the second most common cause of invasive candidiasis and is widely known to have reduced susceptibility to fluconazole relative to many other Candida spp. Upc2A is a transcription factor that regulates ergosterol biosynthesis gene expression under conditions of sterol stress such as azole drug treatment or hypoxia. Through an in vitro microevolution experiment, we found that loss-of-function mutants of the ATF/CREB transcription factor CST6 suppresses the fluconazole hyper-susceptibility of the upc2A∆ mutant. Here, we confirm that the cst6∆ upc2A∆ mutants are resistant to fluconazole but not to hypoxia relative to the upc2A∆ mutant. Sterol analysis of these mutants indicates that this suppression phenotype is not due to restoration of ergosterol levels in the cst6∆ upc2A∆ mutant. Furthermore, increased expression of CDR1, the efflux pump implicated in the vast majority of azole-resistant C. glabrata strains, does not account for the suppression phenotype. Instead, our data suggest that this effect is due in part to increased expression of the adhesin EPA3, which has been shown by others to reduce fluconazole susceptibility in C. glabrata. In addition, we find that loss of both UPC2A and CST6 reduces the expression of mitochondrial and respiratory genes and that this also contributes to the suppression phenotype as well as to the resistance of cst6∆ to fluconazole. These latter data further emphasize the connection between mitochondrial function and azole susceptibility.

Transposon-sequencing (Tn-seq) of the Candida glabrata reference strain CBS138 reveals epigenetic plasticity, structural variation, and intrinsic mechanisms of resistance to micafungin

Candida glabrata (also called Nakaseomyces glabratus) is an opportunistic pathogen that can resist common antifungals and rapidly acquire multidrug resistance. A large amount of genetic variation exists between isolates, which complicates generalizations. Portable transposon-sequencing (Tn-seq) methods can efficiently provide genome-wide information on strain differences and genetic mechanisms. Using the Hermes transposon, the CBS138 reference strain and a commonly studied derivative termed 2001 were subjected to Tn-seq in control conditions and after exposure to varying doses of the clinical antifungal micafungin. The approach revealed large differences between these strains, including a 131-kb tandem duplication and a variety of fitness differences. Additionally, both strains exhibited up to 1,000-fold increased transposon accessibility in subtelomeric regions relative to the BG2 strain, indicative of open subtelomeric chromatin in these isolates and large epigenetic variation within the species. Unexpectedly, the Pdr1 transcription factor conferred resistance to micafungin through targets other than CDR1. Other micafungin resistance pathways were also revealed including mannosyltransferase activity and biosynthesis of the lipid precursor sphingosine, the inhibition of which by SDZ 90-215 and myriocin enhanced the potency of micafungin in vitro. These findings provide insights into the complexity of the C. glabrata species as well as strategies for improving antifungal efficacy.

Tn-seq of the Candida glabrata reference strain CBS138 reveals epigenetic plasticity, structural variation, and intrinsic mechanisms of resistance to micafungin

C. glabrata is an opportunistic pathogen that can resist common antifungals and rapidly acquire multidrug resistance. A large amount of genetic variation exists between isolates, which complicates generalizations. Portable Tn-seq methods can efficiently provide genome-wide information on strain differences and genetic mechanisms. Using the Hermes transposon, the CBS138 reference strain and a commonly studied derivative termed 2001 were subjected to Tn-seq in control conditions and after exposure to varying doses of the clinical antifungal micafungin. The approach revealed large differences between these strains, including a 131 kb tandem duplication and a variety of fitness differences. Additionally, both strains exhibited up to 1000-fold increased transposon accessibility in subtelomeric regions relative to the BG2 strain, indicative of open subtelomeric chromatin in these isolates and large epigenetic variation within the species. Unexpectedly, the Pdr1 transcription factor conferred resistance to micafungin through targets other than CDR1 . Other micafungin resistance pathways were also revealed including mannosyltransferase activity and biosynthesis of the lipid precursor sphingosine, the drugging of which by SDZ 90-215 or myriocin enhanced the potency of micafungin in vitro . These findings provide insights into complexity of the C. glabrata species as well as strategies for improving antifungal efficacy.

Summary

Candida glabrata is an emerging pathogen with large genetic diversity and genome plasticity. The type strain CBS138 and a laboratory derivative were mutagenized with the Hermes transposon and profiled using Tn-seq. Numerous genes that regulate innate and acquired resistance to an important clinical antifungal were uncovered, including a pleiotropic drug resistance gene (PDR1) and a duplication of part of one chromosome. Compounds that target PDR1 and other genes may augment the potency of existing antifungals.

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.

Recent gene selection and drug resistance underscore clinical adaptation across Candida species

Understanding how microbial pathogens adapt to treatments, humans and clinical environments is key to infer mechanisms of virulence, transmission and drug resistance. This may help improve therapies and diagnostics for infections with a poor prognosis, such as those caused by fungal pathogens, including Candida. Here we analysed genomic variants across approximately 2,000 isolates from six Candida species (C. glabrata, C. auris, C. albicans, C. tropicalis, C. parapsilosis and C. orthopsilosis) and identified genes under recent selection, suggesting a highly complex clinical adaptation. These involve species-specific and convergently affected adaptive mechanisms, such as adhesion. Using convergence-based genome-wide association studies we identified known drivers of drug resistance alongside potentially novel players. Finally, our analyses reveal an important role of structural variants and suggest an unexpected involvement of (para)sexual recombination in the spread of resistance. Our results provide insights on how opportunistic pathogens adapt to human-related environments and unearth candidate genes that deserve future attention.

Whole genome analysis of echinocandin non-susceptible Candida Glabrata clinical isolates: a multi-center study in China

BACKGROUND

Candida glabrata is an important cause of invasive candidiasis. Echinocandins are the first-line treatment of invasive candidiasis caused by C. glabrata. The epidemiological echinocandin sensitivity requires long-term surveillance and the understanding about whole genome characteristics of echinocandin non-susceptible isolates was limited.

RESULTS

The present study investigated the echinocandin susceptibility of 1650 C. glabrata clinical isolates in China from August 2014 to July 2019. The in vitro activity of micafungin was significantly better than those of caspofungin and anidulafungin (P < 0.001), assessed by MIC50/90 values. Whole genome sequencing was conducted on non-susceptible isolates and geography-matched susceptible isolates. Thirteen isolates (0.79%) were resistant to at least one echinocandin. Six isolates (0.36%) were solely intermediate to caspofungin. Common evolutionary analysis of echinocandin-resistant and echinocandin-intermediate isolates revealed genes related with reduced caspofungin sensitivity, including previously identified sphinganine hydroxylase encoding gene SUR2. Genome-wide association study identified SNPs at subtelometric regions that were associated with echinocandin non-susceptibility. In-host evolution of echinocandin resistance of serial isolates revealed an enrichment for non-synonymous mutations in adhesins genes and loss of subtelometric regions containing adhesin genes.

CONCLUSIONS

The echinocandins are highly active against C. glabrata in China with a resistant rate of 0.79%. Echinocandin non-susceptible isolates carried common evolved genes which are related with reduced caspofungin sensitivity. In-host evolution of C. glabrata accompanied intensive changing of adhesins profile.

Genomic Assembly of Clinical Candida glabrata (Nakaseomyces glabrata) Isolates Reveals within-Species Structural Plasticity and Association with In Vitro Antifungal Susceptibility

The opportunistic human pathogen Candida glabrata has become an increasingly important threat to human health, with infections globally characterized by high mortality rates and multidrug resistance. To face this threat, more efficient diagnostic and therapeutic approaches are required, underpinning research to help define the intraspecies epidemiology, genetic variability, and therefore, diagnostic and therapeutic target stability. Previous comparative genetics studies conducted on limited numbers of strains only revealed partial resolution of chromosomal settings. In this study, by combining short- and long-read genome sequencing, phenotypic characterization, and comparative genomics over a large set of strains, we detected strict relationships between large chromosomal rearrangements and phylogenetic clades, genes subjected to different selective pressures, and new sets of genes associated with resistance to antifungals. Overall, these results not only provide a fundamental contribution to our knowledge of C. glabrata evolution and epidemiology but may also lay the foundations for the future development of tailored therapeutic approaches. IMPORTANCE The human pathogen Candida glabrata has become a global threat to human health, with infections characterized by high mortality and multidrug resistance. We have obtained nine fully assembled genomes from clinical isolates through a combination of short- and long-read sequencing approaches. The quality and completeness of such genomes and their subsequent comparison to the broadest set of genomes so far allowed us to pinpoint chromosomal rearrangements in several genomes and detect phylogenetic clades that were not associated with geographic location or isolation source. We identified a new set of genes associated with resistance to antifungals coding for adhesin or adhesin-like proteins, suggesting C. glabrata resists antifungals by forming aggregates or adhering to the host tissue. These results, which provide a fundamental contribution to our knowledge of C. glabrata evolution and epidemiology, may initiate the development of precision medicine interventions for patients with suspected or proven invasive fungal infections.

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

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

Chromosome-level assemblies from diverse clades reveal limited structural and gene content variation in the genome of Candida glabrata

Background

Candida glabrata is an opportunistic yeast pathogen thought to have a large genetic and phenotypic diversity and a highly plastic genome. However, the lack of chromosome-level genome assemblies representing this diversity limits our ability to accurately establish how chromosomal structure and gene content vary across strains.

Results

Here, we expanded publicly available assemblies by using long-read sequencing technologies in twelve diverse strains, obtaining a final set of twenty-one chromosome-level genomes spanning the known C. glabrata diversity. Using comparative approaches, we inferred variation in chromosome structure and determined the pan-genome, including an analysis of the adhesin gene repertoire. Our analysis uncovered four new adhesin orthogroups and inferred a rich ancestral adhesion repertoire, which was subsequently shaped through a still ongoing process of gene loss, gene duplication, and gene conversion.

Conclusions

C. glabrata has a largely stable pan-genome except for a highly variable subset of genes encoding cell wall-associated functions. Adhesin repertoire was established for each strain and showed variability among clades.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01412-1.

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.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

A Pan-Cancer Analysis of Cystatin E/M Reveals Its Dual Functional Effects and Positive Regulation of Epithelial Cell in Human Tumors

Cystatin E/M (CST6), a representative cysteine protease inhibitor, plays both tumor-promoting and tumor-suppressing functions and is pursued as an epigenetically therapeutic target in special cancer types. However, a comprehensive and systematic analysis for CST6 in pan-cancer level is still lacking. In the present study, we explored the expression pattern of CST6 in multiple cancer types across ∼10,000 samples from TCGA (The Cancer Genome Atlas) and ∼8,000 samples from MMDs (Merged Microarray-acquired Datasets). We found that the dynamic expression alteration of CST6 was consistent with dual function in different types of cancer. In addition, we observed that the expression of CST6 was globally regulated by the DNA methylation in its promoter region. CST6 expression was positively correlated with the epithelial cell infiltration involved in epithelial-to-mesenchymal transition (EMT) and proliferation. The relationship between CST6 and tumor microenvironment was also explored. In particular, we found that CST6 serves a protective function in the process of melanoma metastasis. Finally, the clinical association analysis further revealed the dual function of CST6 in cancer, and a combination of the epithelial cell infiltration and CST6 expression could predict the prognosis for SKCM patients. In summary, this first CST6 pan-cancer study improves the understanding of the dual functional effects on CST6 in different types of human cancer.

Candida glabrata Antifungal Resistance and Virulence Factors, a Perfect Pathogenic Combination

In recent years, a progressive increase in the incidence of invasive fungal infections (IFIs) caused by Candida glabrata has been observed. The objective of this literature review was to study the epidemiology, drug resistance, and virulence factors associated with the C. glabrata complex. For this purpose, a systematic review (January 2001-February 2021) was conducted on the PubMed, Scielo, and Cochrane search engines with the following terms: "C. glabrata complex (C. glabrata sensu stricto, C. nivariensis, C. bracarensis)" associated with "pathogenicity" or "epidemiology" or "antibiotics resistance" or "virulence factors" with language restrictions of English and Spanish. One hundred and ninety-nine articles were found during the search. Various mechanisms of drug resistance to azoles, polyenes, and echinocandins were found for the C. glabrata complex, depending on the geographical region. Among the mechanisms found are the overexpression of drug transporters, gene mutations that alter thermotolerance, the generation of hypervirulence due to increased adhesion factors, and modifications in vital enzymes that produce cell wall proteins that prevent the activity of drugs designed for its inhibition. In addition, it was observed that the C. glabrata complex has virulence factors such as the production of proteases, phospholipases, and hemolysins, and the formation of biofilms that allows the complex to evade the host immune response and generate fungal resistance. Because of this, the C. glabrata complex possesses a perfect pathogenetic combination for the invasion of the immunocompromised host.

Genetic Susceptibility to Fungal Infections and Links to Human Ancestry

Over the ages, fungi have associated with different parts of the human body and established symbiotic associations with their host. They are mostly commensal unless there are certain not so well-defined factors that trigger the conversion to a pathogenic state. Some of the factors that induce such transition can be dependent on the fungal species, environment, immunological status of the individual, and most importantly host genetics. In this review, we discuss the different aspects of how host genetics play a role in fungal infection since mutations in several genes make hosts susceptible to such infections. We evaluate how mutations modulate the key recognition between the pathogen associated molecular patterns (PAMP) and the host pattern recognition receptor (PRR) molecules. We discuss the polymorphisms in the genes of the immune system, the way it contributes toward some common fungal infections, and highlight how the immunological status of the host determines fungal recognition and cross-reactivity of some fungal antigens against human proteins that mimic them. We highlight the importance of single nucleotide polymorphisms (SNPs) that are associated with several of the receptor coding genes and discuss how it affects the signaling cascade post-infection, immune evasion, and autoimmune disorders. As part of personalized medicine, we need the application of next-generation techniques as a feasible option to incorporate an individual’s susceptibility toward invasive fungal infections based on predisposing factors. Finally, we discuss the importance of studying genomic ancestry and reveal how genetic differences between the human race are linked to variation in fungal disease susceptibility.

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.

Phytol-Loaded Solid Lipid Nanoparticles as a Novel Anticandidal Nanobiotechnological Approach

Phytol is a diterpene alcohol and can be found as a product of the metabolism of chlorophyll in plants. This compound has been explored as a potential antimicrobial agent, but it is insoluble in water. In this study, we describe a novel approach for an interesting anticandidal drug delivery system containing phytol. Different formulations of phytol-loaded solid lipid nanoparticles (SLN) were designed and tested using a natural lipid, 1,3-distearyl-2-oleyl-glycerol (TG1). Different compositions were considered to obtain three formulations with 1:10, 1:5, and 1:3 w/w phytol/TG1 ratios. All the formulations were prepared by emulsification solvent evaporation method and had their physicochemical properties assessed. The biocompatibility assay was performed in the HEK-293 cell line and the antifungal efficacy was demonstrated in different strains of Candida ssp., including different clinical isolates. Spherical and uniform SLN (<300 nm, PdI < 0.2) with phytol-loading efficiency >65% were achieved. Phytol-loaded SLN showed a dose-dependent cytotoxic effect in the HEK-293 cell line. The three tested formulations of phytol-loaded SLN considerably enhanced the minimal inhibitory concentration of phytol against 15 strains of Candida spp. Considering the clinical isolates, the formulations containing the highest phytol/TG1 ratios showed MICs at 100%. Thus, the feasibility and potential of phytol-loaded SLN was demonstrated in vitro, being a promising nanocarrier for phytol delivery from an anticandidal approach.