Multiple genome analysis of Candida glabrata clinical isolates renders new insights into genetic diversity and drug resistance determinants.

The emergence of drug resistance significantly hampers the treatment of human infections, including those caused by fungal pathogens such as Candida species. Candida glabrata ranks as the second most common cause of candidiasis worldwide, supported by rapid acquisition of resistance to azole and echinocandin antifungals frequently prompted by single nucleotide polymorphisms (SNPs) in resistance associated genes, such as PDR1 (azole resistance) or FKS1/2 (echinocandin resistance). To determine the frequency of polymorphisms and genome rearrangements as the possible genetic basis of C. glabrata drug resistance, we assessed genomic variation across 94 globally distributed isolates with distinct resistance phenotypes, whose sequence is deposited in GenBank. The genomes of three additional clinical isolates were sequenced, in this study, including two azole resistant strains that did not display Gain-Of-Function (GOF) mutations in the transcription factor encoding gene PDR1. Genomic variations in susceptible isolates were used to screen out variants arising from genome diversity and to identify variants exclusive to resistant isolates. More than half of the azole or echinocandin resistant isolates do not possess exclusive polymorphisms in PDR1 or FKS1/2, respectively, providing evidence of alternative genetic basis of antifungal resistance. We also identified copy number variations consistently affecting a subset of chromosomes. Overall, our analysis of the genomic and phenotypic variation across isolates allowed to pinpoint, in a genome-wide scale, genetic changes enriched specifically in antifungal resistant strains, which provides a first step to identify additional determinants of antifungal resistance. Specifically, regarding the newly sequenced strains, a set of mutations/genes are proposed to underlie the observed unconventional azole resistance phenotype.

Genomic evolution towards azole resistance in Candida glabrata clinical isolates unveils the importance of CgHxt4/6/7 in azole accumulation.

The increasing prevalence of candidosis caused by Candida glabrata is related to its ability to acquire azole resistance. Although azole resistance mechanisms are well known, the mechanisms for azole import into fungal cells have remained obscure. In this work, we have characterized two hexose transporters in C. glabrata and further investigate their role as potential azole importers. Three azole susceptible C. glabrata clinical isolates were evolved towards azole resistance and the acquired resistance phenotype was found to be independent of CgPDR1 or CgERG11 mutations. Through whole-genome sequencing, CgHXT4/6/7 was found to be mutated in the three evolved strains, when compared to their susceptible parents. CgHxt4/6/7 and the 96% identical CgHxt6/7 were found to confer azole susceptibility and increase azole accumulation in C. glabrata cells, strikingly rescuing the susceptibility phenotype imposed by CgPDR1 deletion, while the identified loss-of-function mutation in CgHXT4/6/7, leads to increased azole resistance. In silico docking analysis shows that azoles display a strong predicted affinity for the glucose binding site of CgHxt4/6/7. Altogether, we hypothesize that hexose transporters, such as CgHxt4/6/7 and CgHxt6/7, may constitute a family of azole importers, involved in clinical drug resistance in fungal pathogens, and constituting promising targets for improved antifungal therapy.

Prediction of Gene and Genomic Regulation in Candida Species, Using the PathoYeastract Database: A Comparative Genomics Approach.

The ability of living organisms to survive changing environmental conditions is dependent on the implementation of gene expression programs underlying adaptation and fitness. Transcriptional networks can be exceptionally complex: a single transcription factor (TF) may regulate hundreds of genes, and multiple TFs may regulate a single gene-depending on the environmental conditions. Moreover, the same TF may act as an activator or repressor in distinct conditions. In turn, the activity of regulators themselves may be dependent on other TFs, as well as posttranscriptional and posttranslational regulation. These traits greatly contribute to the intricate networks governing gene expression programs.In this chapter, a step-by-step guide of how to use PathoYeastract, one of several interconnecting databases within the YEASTRACT+ portal, to predict gene and genomic regulation in Candida spp. is provided. PathoYeastract contains a set of analysis tools to study regulatory associations in human pathogenic yeasts, enabling: (1) the prediction and ranking of TFs that contribute to the regulation of individual genes; (2) the prediction of the genes regulated by a given TF; and (3) the prediction and ranking of TFs that regulate a genome-wide transcriptional response. These capabilities are illustrated, respectively, with the analysis of: (1) the TF network controlling the C. glabrata QDR2 gene; (2) the regulon controlled by the C. glabrata TF Rpn4; and (3) the regulatory network controlling the C. glabrata transcriptome-wide changes induced upon exposure to the antifungal drug fluconazole. The newest potentialities of this information system are explored, including cross-species network comparison. The results are discussed considering the performed queries and integrated with the current knowledge on the biological data for each case-study.

Erg25 Controls Host-Cholesterol Uptake Mediated by Aus1p-Associated Sterol-Rich Membrane Domains in Candida glabrata.

The uptake of cholesterol from the host is closely linked to the proliferation of pathogenic fungi and protozoa during infection. For some pathogenic fungi, cholesterol uptake is an important strategy for decreasing susceptibility to antifungals that inhibit ergosterol biosynthesis. In this study, we show that Candida glabrata ERG25, which encodes an enzyme that demethylates 4,4-dimethylzymosterol, is required for cholesterol uptake from host serum. Based on the screening of C. glabrata conditional knockdown mutants for each gene involved in ergosterol biosynthesis, ERG25 knockdown was found to decrease lethality of infected mice. ERG25 knockdown impairs the plasma membrane localization of the sterol importer Aus1p, suggesting that the accumulated 4,4-dimethylzymosterol destabilizes the lipid domain with which Aus1p functionally associates. ERG25 knockdown further influences the structure of the membrane compartment of Can1p (MCC)/eisosomes (ergosterol-rich lipid domains), but not the localization of the membrane proteins Pma1p and Hxt1p, which localize to sterol-poor domains. In the sterol-rich lipid domain, Aus1p-contining domain was mostly independent of MCC/eisosomes, and the nature of these domains was also different: Ausp1-contining domain was a dynamic network-like domain, whereas the MCC/eisosomes was a static dot-like domain. However, deletion of MCC/eisosomes was observed to influence the localization of Aus1p after Aus1p was transported from the endoplasmic reticulum (ER) through the Golgi apparatus to the plasma membrane. These findings suggest that ERG25 plays a key role in stabilizing sterol-rich lipid domains, constituting a promising candidate target for antifungal therapy.

Characterization of the Candida glabrata Transcription Factor CgMar1: Role in Azole Susceptibility.

The prevalence of antifungal resistance in Candida glabrata, especially against azole drugs, results in difficult-to-treat and potentially life-threatening infections. Understanding the molecular basis of azole resistance in C. glabrata is crucial to designing more suitable therapeutic strategies. In this study, the role of the transcription factor encoded by ORF CAGL0B03421g, here denominated as CgMar1 (Multiple Azole Resistance 1), in azole susceptibility was explored. Using RNA-sequencing, CgMar1 was found to regulate 337 genes under fluconazole stress, including several related to lipid biosynthesis pathways. In this context, CgMar1 and its target CgRSB1, encoding a predicted sphingoid long-chain base efflux transporter, were found to contribute to plasma membrane sphingolipid incorporation and membrane permeability, decreasing fluconazole accumulation. CgMar1 was found to associate with the promoter of CgRSB1, which contains two instances of the CCCCTCC consensus, found to be required for CgRSB1 activation during fluconazole stress. Altogether, a regulatory pathway modulating azole susceptibility in C. glabrata is proposed, resulting from what appears to be a neofunctionalization of a Hap1-like transcription factor.

Marine Sponge and Octocoral-Associated Bacteria Show Versatile Secondary Metabolite Biosynthesis Potential and Antimicrobial Activities against Human Pathogens.

Marine microbiomes are prolific sources of bioactive natural products of potential pharmaceutical value. This study inspected two culture collections comprising 919 host-associated marine bacteria belonging to 55 genera and several thus-far unclassified lineages to identify isolates with potentially rich secondary metabolism and antimicrobial activities. Seventy representative isolates had their genomes mined for secondary metabolite biosynthetic gene clusters (SM-BGCs) and were screened for antimicrobial activities against four pathogenic bacteria and five pathogenic Candida strains. In total, 466 SM-BGCs were identified, with antimicrobial peptide- and polyketide synthase-related SM-BGCs being frequently detected. Only 38 SM-BGCs had similarities greater than 70% to SM-BGCs encoding known compounds, highlighting the potential biosynthetic novelty encoded by these genomes. Cross-streak assays showed that 33 of the 70 genome-sequenced isolates were active against at least one Candida species, while 44 isolates showed activity against at least one bacterial pathogen. Taxon-specific differences in antimicrobial activity among isolates suggested distinct molecules involved in antagonism against bacterial versus Candida pathogens. The here reported culture collections and genome-sequenced isolates constitute a valuable resource of understudied marine bacteria displaying antimicrobial activities and potential for the biosynthesis of novel secondary metabolites, holding promise for a future sustainable production of marine drug leads.

The N.C.Yeastract and CommunityYeastract databases to study gene and genomic transcription regulation in non-conventional yeasts.

Responding to the recent interest of the yeast research community in non-Saccharomyces cerevisiae species of biotechnological relevance, the N.C.Yeastract (http://yeastract-plus.org/ncyeastract/) was associated to YEASTRACT + (http://yeastract-plus.org/). The YEASTRACT + portal is a curated repository of known regulatory associations between transcription factors (TFs) and target genes in yeasts. N.C.Yeastract gathers all published regulatory associations and TF-binding sites for Komagataellaphaffii (formerly Pichia pastoris), the oleaginous yeast Yarrowia lipolytica, the lactose fermenting species Kluyveromyces lactis and Kluyveromyces marxianus, and the remarkably weak acid-tolerant food spoilage yeast Zygosaccharomyces bailii. The objective of this review paper is to advertise the update of the existing information since the release of N.C.Yeastract in 2019, and to raise awareness in the community about its potential to help the day-to-day work on these species, exploring all the information available in the global YEASTRACT + portal. Using simple and widely used examples, a guided exploitation is offered for several tools: (i) inference of orthologous genes; (ii) search for putative TF binding sites and (iii) inter-species comparison of transcription regulatory networks and prediction of TF-regulated networks based on documented regulatory associations available in YEASTRACT + for well-studied species. The usage potentialities of the new CommunityYeastract platform by the yeast community are also discussed.

From a genome assembly to full regulatory network prediction: the case study of Rhodotorula toruloides putative Haa1-regulon.

Numerous genomes are sequenced and made available to the community through the NCBI portal. However, and, unlike what happens for gene function annotation, annotation of promoter sequences and the underlying prediction of regulatory associations is mostly unavailable, severely limiting the ability to interpret genome sequences in a functional genomics perspective. Here we present an approach where one can download a genome of interest from NCBI in the GenBank Flat File (.gbff) format and, with a minimum set of commands, have all the information parsed, organized and made available through the platform web interface. Also, the new genomes are compared with a given genome of reference in search of homologous genes, shared regulatory elements and predicted transcription associations. We present this approach within the context of Community YEASTRACT of the YEASTRACT + portal, thus benefiting from immediate access to all the comparative genomics queries offered in the YEASTRACT + portal. Besides the yeast community, other communities can install the platform independently, without any constraints. In this work, we exemplify the usefulness of the presented tool, within Community YEASTRACT, in constructing a dedicated database and analysing the genome of the highly promising oleaginous red yeast species Rhodotorula toruloides currently poorly studied at the genome and transcriptome levels and with limited genome editing tools. Regulatory prediction is based on the conservation of promoter sequences and available regulatory networks. The case-study examined is focused on the Haa1 transcription factor-a key regulator of yeast resistance to acetic acid, an important inhibitor of industrial bioconversion of lignocellulosic hydrolysates. The new tool described here led to the prediction of a RtHaa1 regulon with expected impact in the optimization of R. toruloides robustness for lignocellulosic and pectin-rich residue biorefinery processes.

From the first touch to biofilm establishment by the human pathogen Candida glabrata: a genome-wide to nanoscale view.

Candida glabrata is an opportunistic pathogen that adheres to human epithelial mucosa and forms biofilm to cause persistent infections. In this work, Single-cell Force Spectroscopy (SCFS) was used to glimpse at the adhesive properties of C. glabrata as it interacts with clinically relevant surfaces, the first step towards biofilm formation. Following a genetic screening, RNA-sequencing revealed that half of the entire transcriptome of C. glabrata is remodeled upon biofilm formation, around 40% of which under the control of the transcription factors CgEfg1 and CgTec1. Using SCFS, it was possible to observe that CgEfg1, but not CgTec1, is necessary for the initial interaction of C. glabrata cells with both abiotic surfaces and epithelial cells, while both transcription factors orchestrate biofilm maturation. Overall, this study characterizes the network of transcription factors controlling massive transcriptional remodelling occurring from the initial cell-surface interaction to mature biofilm formation.

Role of CgTpo4 in Polyamine and Antimicrobial Peptide Resistance: Determining Virulence in Candida glabrata.

Candida glabrata is an emerging fungal pathogen whose success depends on its ability to resist antifungal drugs but also to thrive against host defenses. In this study, the predicted multidrug transporter CgTpo4 (encoded by ORF CAGL0L10912g) is described as a new determinant of virulence in C. glabrata, using the infection model Galleria mellonella. The CgTPO4 gene was found to be required for the C. glabrata ability to kill G. mellonella. The transporter encoded by this gene is also necessary for antimicrobial peptide (AMP) resistance, specifically against histatin-5. Interestingly, G. mellonella's AMP expression was found to be strongly activated in response to C. glabrata infection, suggesting AMPs are a key antifungal defense. CgTpo4 was also found to be a plasma membrane exporter of polyamines, especially spermidine, suggesting that CgTpo4 is able to export polyamines and AMPs, thus conferring resistance to both stress agents. Altogether, this study presents the polyamine exporter CgTpo4 as a determinant of C. glabrata virulence, which acts by protecting the yeast cells from the overexpression of AMPs, deployed as a host defense mechanism.

Transcriptome-wide differences between Saccharomyces cerevisiae and Saccharomyces cerevisiae var. boulardii: Clues on host survival and probiotic activity based on promoter sequence variability.

Although Saccharomyces cerevisiae and S. cerevisiae var. boulardii share more than 95% genome sequence homology, only S. cerevisiae var. boulardii displays probiotic activity. In this study, the transcriptomic differences exhibited by S. cerevisiae and S. cerevisiae var. boulardii in intestinal like medium were evaluated. S. cerevisiae was found to display stress response overexpression, consistent with higher ability of S. cerevisiae var. boulardii to survive within the human host, while S. cerevisiae var. boulardii exhibited transcriptional patterns associated with probiotic activity, suggesting increased acetate biosynthesis. Resorting to the creation of a S. cerevisiae var. boulardii genomic database within Yeastract+, a possible correlation between loss or gain of transcription factor binding sites in S. cerevisiae var. boulardii promoters and the transcriptomic pattern is discussed. This study suggests that S. cerevisiae var. boulardii probiotic activity, when compared to S. cerevisiae, relies, at least partially, on differential expression regulation, based on promoter variability.

A new regulator in the crossroads of oxidative stress resistance and virulence in Candida glabrata: The transcription factor CgTog1.

Candida glabrata is a prominent pathogenic yeast which exhibits a unique ability to survive the harsh environment of host immune cells. In this study, we describe the role of the transcription factor encoded by the gene CAGL0F09229g, here named CgTog1 after its Saccharomyces cerevisiae ortholog, as a new determinant of C. glabrata virulence. Interestingly, Tog1 is absent in the other clinically relevant Candida species (C. albicans, C. parapsilosis, C. tropicalis, C. auris), being exclusive to C. glabrata. CgTog1 was found to be required for oxidative stress resistance and for the modulation of reactive oxygen species inside C. glabrata cells. Also, CgTog1 was observed to be a nuclear protein, whose activity up-regulates the expression of 147 genes and represses 112 genes in C. glabrata cells exposed to H2O2, as revealed through RNA-seq-based transcriptomics analysis. Given the importance of oxidative stress response in the resistance to host immune cells, the effect of CgTOG1 expression in yeast survival upon phagocytosis by Galleria mellonella hemocytes was evaluated, leading to the identification of CgTog1 as a determinant of yeast survival upon phagocytosis. Interestingly, CgTog1 targets include many whose expression changes in C. glabrata cells after engulfment by macrophages, including those involved in reprogrammed carbon metabolism, glyoxylate cycle and fatty acid degradation. In summary, CgTog1 is a new and specific regulator of virulence in C. glabrata, contributing to oxidative stress resistance and survival upon phagocytosis by host immune cells.

Assessing regulatory features of the current transcriptional network of Saccharomyces cerevisiae.

The capacity of living cells to adapt to different environmental, sometimes adverse, conditions is achieved through differential gene expression, which in turn is controlled by a highly complex transcriptional network. We recovered the full network of transcriptional regulatory associations currently known for Saccharomyces cerevisiae, as gathered in the latest release of the YEASTRACT database. We assessed topological features of this network filtered by the kind of supporting evidence and of previously published networks. It appears that in-degree distribution, as well as motif enrichment evolve as the yeast transcriptional network is being completed. Overall, our analyses challenged some results previously published and confirmed others. These analyses further pointed towards the paucity of experimental evidence to support theories and, more generally, towards the partial knowledge of the complete network.

Saccharomyces boulardii: What Makes It Tick as Successful Probiotic?

Saccharomyces boulardii is a probiotic yeast often used for the treatment of GI tract disorders such as diarrhea symptoms. It is genetically close to the model yeast Saccharomyces cerevisiae and its classification as a distinct species or a S. cerevisiae variant has long been discussed. Here, we review the main genetic divergencies between S. boulardii and S. cerevisiae as a strategy to uncover the ability to adapt to the host physiological conditions by the probiotic. S. boulardii does possess discernible phenotypic traits and physiological properties that underlie its success as probiotic, such as optimal growth temperature, resistance to the gastric environment and viability at low pH. Its probiotic activity has been elucidated as a conjunction of multiple pathways, ranging from improvement of gut barrier function, pathogen competitive exclusion, production of antimicrobial peptides, immune modulation, and trophic effects. This review summarizes the participation of S. boulardii in these mechanisms and the multifactorial nature by which this yeast modulates the host microbiome and intestinal function.

Candida glabrata Transcription Factor Rpn4 Mediates Fluconazole Resistance through Regulation of Ergosterol Biosynthesis and Plasma Membrane Permeability.

The ability to acquire azole resistance is an emblematic trait of the fungal pathogen Candida glabrata Understanding the molecular basis of azole resistance in this pathogen is crucial for designing more suitable therapeutic strategies. This study shows that the C. glabrata transcription factor (TF) CgRpn4 is a determinant of azole drug resistance. RNA sequencing during fluconazole exposure revealed that CgRpn4 regulates the expression of 212 genes, activating 80 genes and repressing, likely in an indirect fashion, 132 genes. Targets comprise several proteasome and ergosterol biosynthesis genes, including ERG1, ERG2, ERG3, and ERG11 The localization of CgRpn4 to the nucleus increases upon fluconazole stress. Consistent with a role in ergosterol and plasma membrane homeostasis, CgRpn4 is required for the maintenance of ergosterol levels upon fluconazole stress, which is associated with a role in the upkeep of cell permeability and decreased intracellular fluconazole accumulation. We provide evidence that CgRpn4 directly regulates ERG11 expression through the TTGCAAA binding motif, reinforcing the relevance of this regulatory network in azole resistance. In summary, CgRpn4 is a new regulator of the ergosterol biosynthesis pathway in C. glabrata, contributing to plasma membrane homeostasis and, thus, decreasing azole drug accumulation.

Screening the Drug:H+ Antiporter Family for a Role in Biofilm Formation in Candida glabrata.

Biofilm formation and drug resistance are two key pathogenesis traits exhibited by Candida glabrata as a human pathogen. Interestingly, specific pathways appear to be in the crossroad between the two phenomena, making them promising targets for drug development. In this study, the 10 multidrug resistance transporters of the Drug:H+ Antiporter family of C. glabrata were screened for a role in biofilm formation. Besides previously identified players in this process, namely CgTpo1_2 and CgQdr2, two others are shown to contribute to biofilm formation: CgDtr1 and CgTpo4. The deletion of each of these genes was found to lead to lower biofilm formation, in both SDB and RPMI media, while their expression was found to increase during biofilm development and to be controlled by the transcription factor CgTec1, a predicted key regulator of biofilm formation. Additionally, the deletion of CgDTR1, CgTPO4, or even CgQDR2 was found to increase plasma membrane potential and lead to decreased expression of adhesin encoding genes, particularly CgALS1 and CgEPA1, during biofilm formation. Although the exact role of these drug transporters in biofilm formation remains elusive, our current model suggests that their control over membrane potential by the transport of charged molecules, may affect the perception of nutrient availability, which in turn may delay the triggering of adhesion and biofilm formation.

YEASTRACT+: a portal for cross-species comparative genomics of transcription regulation in yeasts.

The YEASTRACT+ information system (http://YEASTRACT-PLUS.org/) is a wide-scope tool for the analysis and prediction of transcription regulatory associations at the gene and genomic levels in yeasts of biotechnological or human health relevance. YEASTRACT+ is a new portal that integrates the previously existing YEASTRACT (http://www.yeastract.com/) and PathoYeastract (http://pathoyeastract.org/) databases and introduces the NCYeastract (Non-Conventional Yeastract) database (http://ncyeastract.org/), focused on the so-called non-conventional yeasts. The information in the YEASTRACT database, focused on Saccharomyces cerevisiae, was updated. PathoYeastract was extended to include two additional pathogenic yeast species: Candida parapsilosis and Candida tropicalis. Furthermore, the NCYeastract database was created, including five biotechnologically relevant yeast species: Zygosaccharomyces baillii, Kluyveromyces lactis, Kluyveromyces marxianus, Yarrowia lipolytica and Komagataella phaffii. The YEASTRACT+ portal gathers 289 706 unique documented regulatory associations between transcription factors (TF) and target genes and 420 DNA binding sites, considering 247 TFs from 10 yeast species. YEASTRACT+ continues to make available tools for the prediction of the TFs involved in the regulation of gene/genomic expression. In this release, these tools were upgraded to enable predictions based on orthologous regulatory associations described for other yeast species, including two new tools for cross-species transcription regulation comparison, based on multi-species promoter and TF regulatory network analyses.

Draft Genome Sequences of Three Clinical Isolates of the Pathogenic Yeast Candida glabrata.

Here, we report the draft genome sequences of three Candida glabrata clinical isolates, 040, 044, and OL152. The isolates were recovered from patients admitted to Centro Hospitalar de S. João (CHSJ) in Porto, Portugal. Isolates 040 and 044 were taken from blood samples, while isolate OL152 was collected from urine.

Divergent Approaches to Virulence in C. albicans and C. glabrata: Two Sides of the Same Coin.

Candida albicans and Candida glabrata are the two most prevalent etiologic agents of candidiasis worldwide. Although both are recognized as pathogenic, their choice of virulence traits is highly divergent. Indeed, it appears that these different approaches to fungal virulence may be equally successful in causing human candidiasis. In this review, the virulence mechanisms employed by C. albicans and C. glabrata are analyzed, with emphasis on the differences between the two systems. Pathogenesis features considered in this paper include dimorphic growth, secreted enzymes and signaling molecules, and stress resistance mechanisms. The consequences of these traits in tissue invasion, biofilm formation, immune system evasion, and macrophage escape, in a species dependent manner, are discussed. This review highlights the observation that C. albicans and C. glabrata follow different paths leading to a similar outcome. It also highlights the lack of knowledge on some of the specific mechanisms underlying C. glabrata pathogenesis, which deserve future scrutiny.