Comparative genomics of emerging pathogens in the Candida glabrata clade

Candida glabrata maintains two HAP1 ohnologs, HAP1A and HAP1B, for distinct roles in ergosterol gene regulation to mediate sterol homeostasis under azole and hypoxic conditions

Candida glabrata exhibits innate resistance to azole antifungal drugs but also has the propensity to rapidly develop clinical drug resistance. Azole drugs, which target Erg11, is one of the major classes of antifungals used to treat Candida infections. Despite their widespread use, the mechanism controlling azole-induced ERG gene expression and drug resistance in C. glabrata has primarily revolved around Upc2 and/or Pdr1. Phylogenetic and syntenic analyses revealed that C. glabrata, following a whole genome duplication event, maintained HAP1A and HAP1B, whereas Saccharomyces cerevisiae only retained the HAP1A ortholog, HAP1. In this study, we determined the function of two zinc cluster transcription factors, Hap1A and Hap1B, as direct regulators of ERG genes. In S. cerevisiae, Hap1, an ortholog of Hap1A, is a known transcription factor controlling ERG gene expression under aerobic and hypoxic conditions. Interestingly, deleting HAP1 or HAP1B in either S. cerevisiae or C. glabrata, respectively, showed altered susceptibility to azoles. In contrast, the strain deleted for HAP1A did not exhibit azole susceptibility. We also determined that the increased azole susceptibility in a hap1BΔ strain is attributed to decreased azole-induced expression of ERG genes, resulting in decreased levels of total ergosterol. Surprisingly, Hap1A protein expression is barely detected under aerobic conditions but is specifically induced under hypoxic conditions, where Hap1A is required for the repression of ERG genes. However, in the absence of Hap1A, Hap1B can compensate as a transcriptional repressor. Our study shows that Hap1A and Hap1B is utilized by C. glabrata to adapt to specific host and environmental conditions.

IMPORTANCE

Invasive and drug-resistant fungal infections pose a significant public health concern. Candida glabrata, a human fungal pathogen, is often difficult to treat due to its intrinsic resistance to azole antifungal drugs and its capacity to develop clinical drug resistance. Therefore, understanding the pathways that facilitate fungal growth and environmental adaptation may lead to novel drug targets and/or more efficacious antifungal therapies. While the mechanisms of azole resistance in Candida species have been extensively studied, the roles of zinc cluster transcription factors, such as Hap1A and Hap1B, in C. glabrata have remained largely unexplored until now. Our research shows that these factors play distinct yet crucial roles in regulating ergosterol homeostasis under azole drug treatment and oxygen-limiting growth conditions. These findings offer new insights into how this pathogen adapts to different environmental conditions and enhances our understanding of factors that alter drug susceptibility and/or resistance.

Chromosome level assemblies of Nakaseomyces (Candida) bracarensis uncover two distinct clades and define its adhesin repertoire

BACKGROUND

The Nakaseomyces clade is formed by at least nine described species among which three can be pathogenic to humans, namely Nakaseomyces glabratus (Candida glabrata), the second most-common cause of candidiasis worldwide, and two rarer emerging pathogens: Nakaseomyces (Candida) nivarensis and Nakaseomyces (Candida) bracarensis. Early comparative genomics analyses identified parallel expansions of subtelomeric adhesin genes in N. glabratus and N. nivarensis/bracarensis, and suggested possible links with the emergence of the virulence potential in these species. However, as shown for N. glabratus, the proper assessment of subtelomeric genes is hindered by the use of incomplete assemblies and reliance on a single isolate.

RESULTS

Here we sequenced seven N. bracarensis isolates and reconstructed chromosome level assemblies of two divergent strains. We show that N. bracarensis isolates belong to two diverging clades that have slightly different genomic structures. We identified the set of encoded adhesins in the two complete assemblies, and uncovered the presence of a novel adhesin motif, found mainly in N. bracarensis. Our analysis revealed a larger adhesin content in N. bracarensis than previously reported, and similar in size to that of N. glabratus. We confirm the independent adhesin expansion in these two species, which could relate to their different levels of virulence.

CONCLUSION

N. bracarensis clinical isolates belong to at least two differentiated clades. We describe a novel repeat motif found in N. bracarensis adhesins, which helps in their identification. Adhesins underwent independent expansions in N. glabratus and N. bracarensis, leading to repertoires that are qualitatively different but quantitatively similar. Given that adhesins are considered virulence factors, some of the observed differences could contribute to variations in virulence capabilities between N. glabratus and N. bracarensis.

Chromosome-level genome assembly of an auxotrophic strain of the pathogenic yeast Candida parapsilosis

We report the genome sequence of the pathogenic yeast Candida parapsilosis strain SR23 (CBS 7157) used in a number of experimental studies. The nuclear genome assembly consists of eight chromosome-sized contigs with a total size of 13.04 Mbp (N50 2.09 Mbp) and a G+C content of 38.7%.

In vivo CRISPR-Cas9 expression in Candida glabrata, Candida bracarensis, and Candida nivariensis: A versatile tool to study chromosomal break repair

The CRISPR-Cas9 system is extremely useful for genome editing in many species, including the model yeast Saccharomyces cerevisiae, and other yeast species. We have previously reported the use of an inducible CRISPR-Cas9 system in Candida glabrata, which allows genome editing but also the study of double-strand break (DSB) repair. We report, in this study, a comparable system for C. glabrata, relying on a new plasmid, which is more stable than the previous one. We also report the use of this plasmid to induce DSBs in two additional human pathogens, Candida bracarensis and Candida nivariensis. We examine lethality induced by an in vivo DSB in the three species and describe the different types of nonhomologous end-joining (NHEJ) events detected in these three pathogens.

Simplified J774A.1 macrophage assay for fungal pathogenicity demonstrates non-clinical Nakaseomyces glabratus strains survive better than lab strains

Nakaseomyces glabratus (formerly known as Candida glabrata ) is the second most common cause of candidiasis, whereas the closely related yeast, Saccharomyces cerevisiae, causes few infections. Macrophages can control N. glabratus infections through phagocytosis, but in cell culture, N. glabratus is able to persist in macrophages better than non-pathogenic yeast. Using J774A.1 macrophages, we simplified a standard persistence/survival assay by counting yeast cells with flow cytometry and incorporating an antifungal treatment. These improvements minimized wash steps and variation so fewer replicates were needed. Here, we demonstrate that loss of NgTUP11 does not lower pathogenicity, and that three non-clinical N. glabratus strains survive in macrophages better than a laboratory strain.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

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

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

The Enigma of NTH2 Gene in Yeasts

The enzymatic hydrolysis of the non-reducing disaccharide trehalose in yeasts is carried out by trehalase, a highly specific α-glucosidase. Two types of such trehalase activity are present in yeasts, and are referred to as neutral and acid enzymes. They are encoded by distinct genes (NTH1 and ATH1, respectively) and exhibit strong differences in their biochemical and physiological properties as well as different subcellular location and regulatory mechanisms. Whereas a single gene ATH1 codes for acid trehalase, the genome of some yeasts appears to predict the existence of a second redundant neutral trehalase, encoded by the NTH2 gene, a paralog of NTH1. In S. cerevisiae the corresponding two proteins share 77% amino acid identity, leading to the suggestion that NTH2 codes for a functional trehalase activity. However, Nth2p lacks any measurable neutral trehalase activity and disruption of NTH2 gene has no effect on this activity compared to a parental strain. Likewise, single nth1Δ and double nth1Δ/nth2Δ null mutants display no detectable neutral activity. Furthermore, disruption of NTH2 does not cause any apparent phenotype apart from a slight involvement in thermotolerance. To date, no evidence of a duplicated NTH gene has been recorded in other archetypical yeasts, like C. albicans or C. parapsilosis, and a possible regulatory mechanism of Nth2p remains unknown. Therefore, although genomic analysis points to the existence, in some yeasts, of two distinct genes encoding trehalase activities, the large body of biochemical and physiological evidence gathered from NTH2 gene does not support this proposal. Indeed, much more experimental evidence would be necessary to firmly validate this hypothesis.

Chromosome-level genome assembly of the yeast Lodderomyces beijingensis reveals the genetic nature of metabolic adaptations and identifies subtelomeres as hotspots for amplification of mating type loci

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterized its genome sequence.

The Hog1 MAPK substrate governs Candida glabrata-epithelial cell adhesion via the histone H2A variant

CgHog1, terminal kinase of the high-osmolarity glycerol signalling pathway, orchestrates cellular response to multiple external stimuli including surplus-environmental iron in the human fungal pathogen Candida glabrata (Cg). However, CgHog1 substrates remain unidentified. Here, we show that CgHog1 adversely affects Cg adherence to host stomach and kidney epithelial cells in vitro, but promotes Cg survival in the iron-rich gastrointestinal tract niche. Further, CgHog1 interactome and in vitro phosphorylation analysis revealed CgSub2 (putative RNA helicase) to be a CgHog1 substrate, with CgSub2 also governing iron homeostasis and host adhesion. CgSub2 positively regulated EPA1 (encodes a major adhesin) expression and host adherence via its interactor CgHtz1 (histone H2A variant). Notably, both CgHog1 and surplus environmental iron had a negative impact on CgSub2-CgHtz1 interaction, with CgHTZ1 or CgSUB2 deletion reversing the elevated adherence of Cghog1Δ to epithelial cells. Finally, the surplus-extracellular iron led to CgHog1 activation, increased CgSub2 phosphorylation, elevated CgSub2-CgHta (canonical histone H2A) interaction, and EPA1 transcriptional activation, thereby underscoring the iron-responsive, CgHog1-induced exchange of histone partners of CgSub2. Altogether, our work mechanistically defines how CgHog1 couples Epa1 adhesin expression with iron abundance, and point towards specific chromatin composition modification programs that probably aid fungal pathogens align their adherence to iron-rich (gut) and iron-poor (blood) host niches.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

SWI/SNF complex-mediated chromatin remodeling in Candida glabrata promotes immune evasion

Immune evasion is critical for fungal virulence. However, how the human opportunistic pathogen Candida glabrata (Cg) accomplishes this is unknown. Here, we present the first genome-wide nucleosome map of the macrophage-internalized Cg consisting of ∼12,000 dynamic and 70,000 total nucleosomes. We demonstrate that CgSnf2 (SWI/SNF chromatin remodeling complex-ATPase subunit)-mediated chromatin reorganization in macrophage-internalized Cg upregulates and downregulates the immunosuppressive seven-gene mannosyltransferase-cluster (CgMT-C) and immunostimulatory cell surface adhesin-encoding EPA1 gene, respectively. Consistently, EPA1 overexpression and CgMT-C deletion elevated IL-1β (pro-inflammatory cytokine) production and diminished Cg proliferation in macrophages. Further, Cgsnf2Δ had higher Epa1 surface expression, and evoked increased IL-1β secretion, and was killed in macrophages. Akt-, p38-, NF-κB- or NLRP3 inflammasome-inhibition partially reversed increased IL-1β secretion in Cgsnf2Δ-infected macrophages. Importantly, macrophages responded to multiple Candida pathogens via NF-κB-dependent IL-1β production, underscoring NF-κB signaling's role in fungal diseases. Altogether, our findings directly link the nucleosome positioning-based chromatin remodeling to fungal immunomodulatory molecule expression.

Genome-wide profiling of piggyBac transposon insertion mutants reveals loss of the F1F0 ATPase complex causes fluconazole resistance in Candida glabrata

Invasive candidiasis caused by non-albicans species has been on the rise, with Candida glabrata emerging as the second most common etiological agent. Candida glabrata possesses an intrinsically lower susceptibility to azoles and an alarming propensity to rapidly develop high-level azole resistance during treatment. In this study, we have developed an efficient piggyBac (PB) transposon-mediated mutagenesis system in C. glabrata to conduct genome-wide genetic screens and applied it to profile genes that contribute to azole resistance. When challenged with the antifungal drug fluconazole, PB insertion into 270 genes led to significant resistance. A large subset of these genes has a role in the mitochondria, including almost all genes encoding the subunits of the F1F0 ATPase complex. We show that deleting ATP3 or ATP22 results in increased azole resistance but does not affect susceptibility to polyenes and echinocandins. The increased azole resistance is due to increased expression of PDR1 that encodes a transcription factor known to promote drug efflux pump expression. Deleting PDR1 in the atp3Δ or atp22Δ mutant resulted in hypersensitivity to fluconazole. Our results shed light on the mechanisms contributing to azole resistance in C. glabrata. This PB transposon-mediated mutagenesis system can significantly facilitate future genome-wide genetic screens.

The oral microbiome: diversity, biogeography and human health

The human oral microbiota is highly diverse and has a complex ecology, comprising bacteria, microeukaryotes, archaea and viruses. These communities have elaborate and highly structured biogeography that shapes metabolic exchange on a local scale and results from the diverse microenvironments present in the oral cavity. The oral microbiota also interfaces with the immune system of the human host and has an important role in not only the health of the oral cavity but also systemic health. In this Review, we highlight recent advances including novel insights into the biogeography of several oral niches at the species level, as well as the ecological role of candidate phyla radiation bacteria and non-bacterial members of the oral microbiome. In addition, we summarize the relationship between the oral microbiota and the pathology of oral diseases and systemic diseases. Together, these advances move the field towards a more holistic understanding of the oral microbiota and its role in health, which in turn opens the door to the study of novel preventive and therapeutic strategies.

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.

RNA viruses, M satellites, chromosomal killer genes, and killer/nonkiller phenotypes in the 100-genomes S. cerevisiae strains

We characterized previously identified RNA viruses (L-A, L-BC, 20S, and 23S), L-A-dependent M satellites (M1, M2, M28, and Mlus), and M satellite-dependent killer phenotypes in the Saccharomyces cerevisiae 100-genomes genetic resource population. L-BC was present in all strains, albeit in 2 distinct levels, L-BChi and L-BClo; the L-BC level is associated with the L-BC genotype. L-BChi, L-A, 20S, 23S, M1, M2, and Mlus (M28 was absent) were in fewer strains than the similarly inherited 2µ plasmid. Novel L-A-dependent phenotypes were identified. Ten M+ strains exhibited M satellite-dependent killing (K+) of at least 1 of the naturally M0 and cured M0 derivatives of the 100-genomes strains; in these M0 strains, sensitivities to K1+, K2+, and K28+ strains varied. Finally, to complement our M satellite-encoded killer toxin analysis, we assembled the chromosomal KHS1 and KHR1 killer genes and used naturally M0 and cured M0 derivatives of the 100-genomes strains to assess and characterize the chromosomal killer phenotypes.

Evolution and functional diversification of yeast sugar transporters

While simple sugars such as monosaccharides and disaccharide are the typical carbon source for most yeasts, whether a species can grow on a particular sugar is generally a consequence of presence or absence of a suitable transporter to enable its uptake. The most common transporters that mediate sugar import in yeasts belong to the major facilitator superfamily (MFS). Some of these, for example the Saccharomyces cerevisiae Hxt proteins have been extensively studied, but detailed information on many others is sparce. In part, this is because there are many lineages of MFS transporters that are either absent from, or poorly represented in, the model S. cerevisiae, which actually has quite a restricted substrate range. It is important to address this knowledge gap to gain better understanding of the evolution of yeasts and to take advantage of sugar transporters to exploit or engineer yeasts for biotechnological applications. This article examines the full repertoire of MFS proteins in representative budding yeasts (Saccharomycotina). A comprehensive analysis of 139 putative sugar transporters retrieved from 10 complete genomes sheds new light on the diversity and evolution of this family. Using the phylogenetic lens, it is apparent that proteins have often been misassigned putative functions and this can now be corrected. It is also often seen that patterns of expansion of particular genes reflects the differential importance of transport of specific sugars (and related molecules) in different yeasts, and this knowledge also provides an improved resource for the selection or design of tailored transporters.

Epidemiology and susceptibility of Nakaseomyces (formerly Candida) glabrata bloodstream isolates from hospitalised adults in South Africa

During 2016-2017, Nakaseomyces glabrata (formerly Candida glabrata) caused 14% of cases of candidaemia in South Africa. We aimed to describe the clinical characteristics of adults with N. glabrata candidaemia at 20 sentinel hospitals (accounting for 20% (172/917) of cases) and the antifungal susceptibility of the corresponding isolates. A higher proportion of patients with N. glabrata candidaemia were older (median age: 55 years [interquartile range (IQR): 41-65 years] vs. 49 years [IQR: 35-63 years]; p = 0.04), female (87/164, 53% vs. 283/671, 42%; p = 0.01), admitted to a public-sector hospital (152/172, 88% vs. 470/745, 63%; p < 0.001), treated with fluconazole only (most with suboptimal doses) (51/95, 54% vs. 139/361, 39%; p < 0.001), and had surgery (47/172, 27% vs. 123/745, 17%; p = 0.001) and a shorter hospital stay (median 7 days [IQR: 2-20 days] vs. 13 days [IQR: 4-27 days]; p < 0.001) compared to patients with other causes of candidaemia. Eight N. glabrata isolates (6%, 8/131) had minimum inhibitory concentrations in the intermediate or resistant range for ≥ 1 echinocandin and a R1377K amino acid substitution encoded by the hotspot 2 region of the FKS2 gene. Only 11 isolates (8%, 11/131) were resistant to fluconazole. Patients with confirmed N. glabrata candidaemia are recommended to be treated with an echinocandin (or polyene), thus further guideline training is required.

Is There a Relationship Between Mating and Pathogenesis in Two Human Fungal Pathogens, Candida albicans and Candida glabrata?

Purpose of Review

Human fungal pathogens are rapidly increasing in incidence and readily able to evade the host immune responses. Our ability to study the genetic behind this has been limited due to the apparent lack of a sexual cycle and forward genetic tools. In this review, we discuss the evolution of mating, meiosis, and pathogenesis and if these processes are advantageous to pathogens.

Recent Findings

This review summarises what is currently known about the sexual cycles of two important human fungal pathogens, Candida albicans and Candida glabrata. This includes the identification of parasexual cycle in C. albicans and the observed low levels of recombination in C. glabrata populations.

Summary

In this review, we present what is currently known about the mating types and mating/sexual cycles of two clinically important human fungal pathogens, Candida albicans and Candida glabrata. We discuss the evolution of meiosis using the knowledge that has been amassed from the decades of studying Saccharomyces cerevisiae and how this can be applied to fungal pathogens. We further discuss how the evolution of pathogenesis has played a role in influencing mating processes in human fungal pathogens and compare sexual cycles between C. albicans and C. glabrata, highlighting knowledge gaps and suggesting how these two fungi have evolved distinct mating niches to allow the development of disease in a human host.

Parallel expansion and divergence of an adhesin family in pathogenic yeasts

Opportunistic yeast pathogens arose multiple times in the Saccharomycetes class, including the recently emerged, multidrug-resistant (MDR) Candida auris. We show that homologs of a known yeast adhesin family in Candida albicans, the Hyr/Iff-like (Hil) family, are enriched in distinct clades of Candida species as a result of multiple, independent expansions. Following gene duplication, the tandem repeat-rich region in these proteins diverged extremely rapidly and generated large variations in length and β-aggregation potential, both of which are known to directly affect adhesion. The conserved N-terminal effector domain was predicted to adopt a β-helical fold followed by an α-crystallin domain, making it structurally similar to a group of unrelated bacterial adhesins. Evolutionary analyses of the effector domain in C. auris revealed relaxed selective constraint combined with signatures of positive selection, suggesting functional diversification after gene duplication. Lastly, we found the Hil family genes to be enriched at chromosomal ends, which likely contributed to their expansion via ectopic recombination and break-induced replication. Combined, these results suggest that the expansion and diversification of adhesin families generate variation in adhesion and virulence within and between species and are a key step toward the emergence of fungal pathogens.

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.

A widespread inversion polymorphism conserved among Saccharomyces species is caused by recurrent homogenization of a sporulation gene family

Saccharomyces genomes are highly collinear and show relatively little structural variation, both within and between species of this yeast genus. We investigated the only common inversion polymorphism known in S. cerevisiae, which affects a 24-kb 'flip/flop' region containing 15 genes near the centromere of chromosome XIV. The region exists in two orientations, called reference (REF) and inverted (INV). Meiotic recombination in this region is suppressed in crosses between REF and INV orientation strains such as the BY x RM cross. We find that the inversion polymorphism is at least 17 million years old because it is conserved across the genus Saccharomyces. However, the REF and INV isomers are not ancient alleles but are continually being re-created by re-inversion of the region within each species. Inversion occurs due to continual homogenization of two almost identical 4-kb sequences that form an inverted repeat (IR) at the ends of the flip/flop region. The IR consists of two pairs of genes that are specifically and strongly expressed during the late stages of sporulation. We show that one of these gene pairs, YNL018C/YNL034W, codes for a protein that is essential for spore formation. YNL018C and YNL034W are the founder members of a gene family, Centroid, whose members in other Saccharomycetaceae species evolve fast, duplicate frequently, and are preferentially located close to centromeres. We tested the hypothesis that Centroid genes are a meiotic drive system, but found no support for this idea.

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.

Confirmation of fifth Candida auris clade by whole genome sequencing

Candida auris has emerged globally as a multidrug-resistant pathogen causing outbreaks in health care facilities. Whole genome sequencing (WGS) analysis has identified four major clades, while earlier WGS data from a single Iranian isolate suggested the existence of a potential fifth clade. Here, we confirm the existence of this fifth clade by providing WGS data of another four Iranian isolates. These clade V isolates differed less than 100 single-nucleotide polymorphisms (SNPs) between each other, while they were separated from the other clades by more than 200,000 SNPs. Two of these isolates were resistant to fluconazole and were found to harbour mutations in the TAC1b and ERG11 genes.

A novel Candida glabrata doxycycline-inducible system for in vitro/in vivo use

Candida glabrata is an important pathogen causing superficial to invasive disease in human. Conditional expression systems are helpful in addressing the function of genes and especially when they can be applied to in vivo studies. Tetracycline-dependent regulation systems have been used in diverse fungi to turn-on (Tet-on) or turn-off (Tet-off) gene expression either in vitro but also in vivo in animal models. Up to now, only a Tet-off expression has been constructed for gene expression in C. glabrata. Here, we report a Tet-on gene expression system which can be used in vitro and in vivo in any C. glabrata genetic background. This system was used in a mice model of systemic infection to demonstrate that the general amino acid permease Gap1 is important for C. glabrata virulence.

Using in vivo transcriptomics and RNA enrichment to identify genes involved in virulence of Candida glabrata

Candida species are the most commonly isolated opportunistic fungal pathogens in humans. Candida albicans causes most of the diagnosed infections, closely followed by Candida glabrata. C. albicans is well studied, and many genes have been shown to be important for infection and colonization of the host. It is however less clear how C. glabrata infects the host. With the help of fungal RNA enrichment, we here investigated for the first time the transcriptomic profile of C. glabrata during urinary tract infection (UTI) in mice. In the UTI model, bladders and kidneys are major target organs and therefore fungal transcriptomes were addressed in these organs. Our results showed that, next to adhesins and proteases, nitrogen metabolism and regulation play a vital role during C. glabrata UTI. Genes involved in nitrogen metabolism were upregulated and among them we show that DUR1,2 (urea amidolyase) and GAP1 (amino acid permease) were important for virulence. Furthermore, we confirmed the importance of the glyoxylate cycle in the host and identified MLS1 (malate synthase) as an important gene necessary for C. glabrata virulence. In conclusion, our study shows with the support of in vivo transcriptomics how C. glabrata adapts to host conditions.

Mating-Type Switching in Budding Yeasts, from Flip/Flop Inversion to Cassette Mechanisms

Mating-type switching is a natural but unusual genetic control process that regulates cell identity in ascomycete yeasts. It involves physically replacing one small piece of genomic DNA by another, resulting in replacement of the master regulatory genes in the mating pathway and hence a switch of cell type and mating behavior. In this review, we concentrate on recent progress that has been made on understanding the origins and evolution of mating-type switching systems in budding yeasts (subphylum Saccharomycotina). Because of the unusual nature and the complexity of the mechanism in Saccharomyces cerevisiae, mating-type switching was assumed until recently to have originated only once or twice during yeast evolution. However, comparative genomics analysis now shows that switching mechanisms arose many times independently-at least 11 times in budding yeasts and once in fission yeasts-a dramatic example of convergent evolution. Most of these lineages switch mating types by a flip/flop mechanism that inverts a section of a chromosome and is simpler than the well-characterized 3-locus cassette mechanism (MAT/HML/HMR) used by S. cerevisiae. Mating-type switching (secondary homothallism) is one of the two possible mechanisms by which a yeast species can become self-fertile. The other mechanism (primary homothallism) has also emerged independently in multiple evolutionary lineages of budding yeasts, indicating that homothallism has been favored strongly by natural selection. Recent work shows that HO endonuclease, which makes the double-strand DNA break that initiates switching at the S. cerevisiae MAT locus, evolved from an unusual mobile genetic element that originally targeted a glycolytic gene, FBA1.

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.

Adhesion as a Focus in Trichoderma-Root Interactions

Fungal spores, germlings, and mycelia adhere to substrates, including host tissues. The adhesive forces depend on the substrate and on the adhesins, the fungal cell surface proteins. Attachment is often a prerequisite for the invasion of the host, hence its importance. Adhesion visibly precedes colonization of root surfaces and outer cortex layers, but little is known about the molecular details. We propose that by starting from what is already known from other fungi, including yeast and other filamentous pathogens and symbionts, the mechanism and function of Trichoderma adhesion will become accessible. There is a sequence, and perhaps functional, homology to other rhizosphere-competent Sordariomycetes. Specifically, Verticillium dahliae is a soil-borne pathogen that establishes itself in the xylem and causes destructive wilt disease. Metarhizium species are best-known as insect pathogens with biocontrol potential, but they also colonize roots. Verticillium orthologs of the yeast Flo8 transcription factor, Som1, and several other relevant genes are already under study for their roles in adhesion. Metarhizium encodes relevant adhesins. Trichoderma virens encodes homologs of Som1, as well as adhesin candidates. These genes should provide exciting leads toward the first step in the establishment of beneficial interactions with roots in the rhizosphere.

Tup1 Paralog CgTUP11 Is a Stronger Repressor of Transcription than CgTUP1 in Candida glabrata

TUP1 is a well-characterized repressor of transcription in Saccharomyces cerevisiae and Candida albicans and is observed as a single-copy gene. We observe that most species that experienced a whole-genome duplication outside of the Saccharomyces genus have two copies of TUP1 in the Saccharomycotina yeast clade. We focused on Candida glabrata and demonstrated that the uncharacterized TUP1 homolog, C. glabrata TUP11 (CgTUP11), is most like the S. cerevisiae TUP1 (ScTUP1) gene through phenotypic assays and transcriptome sequencing (RNA-seq). Whereas CgTUP1 plays a role in gene repression, it is much less repressive in standard growth media. Through RNA-seq and reverse transcription-quantitative PCR (RT-qPCR), we observed that genes associated with pathogenicity (YPS2, YPS4, and HBN1) are upregulated upon deletion of either paralog, and loss of both paralogs is synergistic. Loss of the corepressor CgCYC8 mimics the loss of both paralogs, but not to the same extent as the Cgtup1Δ Cgtup11Δ mutant for these pathogenesis-related genes. In contrast, genes involved in energy metabolism (CgHXT2, CgADY2, and CgFBP1) exhibit similar behavior (dependence on both paralogs), but deletion of CgCYC8 is very similar to the Cgtup1Δ Cgtup11Δ mutant. Finally, some genes (CgMFG1 and CgRIE1) appear to only be dependent on CgTUP11 and CgCYC8 and not CgTUP1. These data indicate separable and overlapping roles for the two TUP1 paralogs and that other genes may function as the CgCyc8 corepressor. Through a comparison by RNA-seq of Sctup1Δ, it was found that TUP1 homologs regulate similar genes in the two species. This work highlights that studies focused only on Saccharomyces may miss important biological processes because of paralog loss after genome duplication. IMPORTANCE Due to a whole-genome duplication, many yeast species related to C. glabrata have two copies of the well-characterized TUP1 gene, unlike most Saccharomyces species. This work identifies roles for the paralogs in C. glabrata, highlights the importance of the uncharacterized paralog, called TUP11, and suggests that the two paralogs have both overlapping and unique functions. The TUP1 paralogs likely influence pathogenicity based on tup mutants upregulating genes that are associated with pathogenicity.

From Genome Variation to Molecular Mechanisms: What we Have Learned From Yeast Mitochondrial Genomes?

Analysis of genome variation provides insights into mechanisms in genome evolution. This is increasingly appreciated with the rapid growth of genomic data. Mitochondrial genomes (mitogenomes) are well known to vary substantially in many genomic aspects, such as genome size, sequence context, nucleotide base composition and substitution rate. Such substantial variation makes mitogenomes an excellent model system to study the mechanisms dictating mitogenome variation. Recent sequencing efforts have not only covered a rich number of yeast species but also generated genomes from abundant strains within the same species. The rich yeast genomic data have enabled detailed investigation from genome variation into molecular mechanisms in genome evolution. This mini-review highlights some recent progresses in yeast mitogenome studies.

The intestinal microbiota and metabolites in patients with anorexia nervosa

Brain-gut microbiota interactions are intensively studied in connection with various neurological and psychiatric diseases. While anorexia nervosa (AN) pathophysiology is not entirely clear, it is presumably linked to microbiome dysbiosis. We aimed to elucidate the gut microbiota contribution in AN disease pathophysiology. We analyzed the composition and diversity of the gut microbiome of patients with AN (bacteriome and mycobiome) from stool samples before and after renourishment, and compared them to healthy controls. Further, levels of assorted neurotransmitters and short-chain fatty acids (SCFA) were analyzed in stool samples by MS and NMR, respectively. Biochemical, anthropometric, and psychometric profiles were assessed. The bacterial alpha-diversity parameter analyses revealed only increased Chao 1 index in patients with AN before the realimentation, reflecting their interindividual variation. Subsequently, core microbiota depletion signs were observed in patients with AN. Overrepresented OTUs (operation taxonomic units) in patients with AN taxonomically belonged to Alistipes, Clostridiales, Christensenellaceae, and Ruminococcaceae. Underrepresented OTUs in patients with AN were Faecalibacterium, Agathobacter, Bacteroides, Blautia, and Lachnospira. Patients exhibited greater interindividual variation in the gut bacteriome, as well as in metagenome content compared to controls, suggesting altered bacteriome functions. Patients had decreased levels of serotonin, GABA, dopamine, butyrate, and acetate in their stool samples compared to controls. Mycobiome analysis did not reveal significant differences in alpha diversity and fungal profile composition between patients with AN and healthy controls, nor any correlation of the fungal composition with the bacterial profile. Our results show the changed profile of the gut microbiome and its metabolites in patients with severe AN. Although therapeutic partial renourishment led to increased body mass index and improved psychometric parameters, SCFA, and neurotransmitter profiles, as well as microbial community compositions, did not change substantially during the hospitalization period, which can be potentially caused by only partial weight recovery.

Analyzing the human gut mycobiome – a short guide for beginners

The human body is a dynamic ecosystem consisting of millions of microbes which are often comprised under the term microbiome. Compared to bacteria, which count for the overwhelming majority of the microbiome, the number of human-associated fungi is small and often underestimated. Nonetheless, they can be found in different host niches such as the gut, the oral cavity and the skin. The fungal community has several potential roles in health and disease of the human host. In this review we will focus on intestinal fungi and their interaction with the host as well as bacteria. We also summarize technical challenges and possible biases researchers must be aware of when conducting mycobiome analysis.

Candida glabrata Hst1-Rfm1-Sum1 complex evolved to control virulence-related genes

C. glabrata is an opportunistic fungal pathogen and the second most common cause of opportunistic fungal infections in humans, that has evolved virulence factors to become a successful pathogen: strong resistance to oxidative stress, capable to adhere and form biofilms in human epithelial cells as well as to abiotic surfaces and high resistance to xenobiotics. Hst1 (a NAD⁺-dependent histone deacetylase), Sum1 (putative DNA binding protein) and Rfm1 (connector protein) form a complex (HRS-C) and control the resistance to oxidative stress, to xenobiotics (the antifungal fluconazole), and adherence to epithelial cells. Hst1 is functionally conserved within the Saccharomycetaceae family, Rfm1 shows a close phylogenetic relation within the Saccharomycetaceae family while Sum1 displays a distant phylogenetic relation with members of the family and is not conserved functionally. CDR1 encodes for an ABC transporter (resistance to fluconazole) negatively controlled by HRS-C, for which its binding site is located within 223 bp upstream from the ATG of CDR1. The absence of Hst1 and Sum1 renders the cells hyper-adherent, possibly due to the overexpression of AED1, EPA1, EPA22 and EPA6, all encoding for adhesins. Finally, in a neutrophil survival assay, HST1 and SUM1, are not required for survival. We propose that Sum1 in the HRS-C diverged functionally to control a set of genes implicated in virulence: adherence, resistance to xenobiotics and oxidative stress.

The involvement of the Candida glabrata trehalase enzymes in stress resistance and gut colonization

Candida glabrata is an opportunistic human fungal pathogen and is frequently present in the human microbiome. It has a high relative resistance to environmental stresses and several antifungal drugs. An important component involved in microbial stress tolerance is trehalose. In this work, we characterized the three C. glabrata trehalase enzymes Ath1, Nth1 and Nth2. Single, double and triple deletion strains were constructed and characterized both in vitro and in vivo to determine the role of these enzymes in virulence. Ath1 was found to be located in the periplasm and was essential for growth on trehalose as sole carbon source, while Nth1 on the other hand was important for oxidative stress resistance, an observation which was consistent by the lower survival rate of the NTH1 deletion strain in human macrophages. No significant phenotype was observed for Nth2. The triple deletion strain was unable to establish a stable colonization of the gastrointestinal (GI) tract in mice indicating the importance of having trehalase activity for colonization in the gut.

The Flo Adhesin Family

The first step in the infection of fungal pathogens in humans is the adhesion of the pathogen to host tissue cells or abiotic surfaces such as catheters and implants. One of the main players involved in this are the expressed cell wall adhesins. Here, we review the Flo adhesin family and their involvement in the adhesion of these yeasts during human infections. Firstly, we redefined the Flo adhesin family based on the domain architectures that are present in the Flo adhesins and their functions, and set up a new classification of Flo adhesins. Next, the structure, function, and adhesion mechanisms of the Flo adhesins whose structure has been solved are discussed in detail. Finally, we identified from Pfam database datamining yeasts that could express Flo adhesins and are encountered in human infections and their adhesin architectures. These yeasts are discussed in relation to their adhesion characteristics and involvement in infections.

The Formation of Neochromosomes during Experimental Evolution in the Yeast Saccharomyces cerevisiae

Novel, large-scale structural mutations were previously discovered during the cultivation of engineered Saccharomyces cerevisiae strains in which essential tRNA synthetase genes were replaced by their orthologs from the distantly related yeast Yarrowia lipolytica. Among those were internal segmental amplifications forming giant chromosomes as well as complex segmental rearrangements associated with massive amplifications at an unselected short locus. The formation of such novel structures, whose stability is high enough to propagate over multiple generations, involved short repeated sequences dispersed in the genome (as expected), but also novel junctions between unrelated sequences likely triggered by accidental template switching within replication forks. Using the same evolutionary protocol, we now describe yet another type of major structural mutation in the yeast genome, the formation of neochromosomes, with functional centromeres and telomeres, made of extra copies of very long chromosomal segments ligated together in novel arrangements. The novel junctions occurred between short repeated sequences dispersed in the genome. They first resulted in the formation of an instable neochromosome present in a single copy in the diploid cells, followed by its replacement by a shorter, partially palindromic neochromosome present in two copies, whose stability eventually increased the chromosome number of the diploid strains harboring it.

Point mutations in Candida glabrata 3-hydroxy-3-methylglutaryl-coenzyme A reductase (CgHMGR) decrease enzymatic activity and substrate/inhibitor affinity

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) is a crucial enzyme in the ergosterol biosynthesis pathway. The aim of this study was to obtain, purify, characterize, and overexpress five point mutations in highly conserved regions of the catalytic domain of Candida glabrata HMGR (CgHMGR) to explore the function of key amino acid residues in enzymatic activity. Glutamic acid (Glu) was substituted by glutamine in the E680Q mutant (at the dimerization site), Glu by glutamine in E711Q (at the substrate binding site), aspartic acid by alanine in D805A, and methionine by arginine in M807R (the latter two at the cofactor binding site). A double mutation, E680Q-M807R, was included. Regarding recombinant and wild-type CgHMGR, in vitro enzymatic activity was significantly lower for the former, as was the in silico binding energy of simvastatin, alpha-asarone and the HMG-CoA substrate. E711Q displayed the lowest enzymatic activity and binding energy, suggesting the importance of Glu711 (in the substrate binding site). The double mutant CgHMGR E680Q-M807R exhibited the second lowest enzymatic activity. Based on the values of the kinetic parameters KM and Vmax, the mutated amino acids appear to participate in binding. The current findings provide insights into the role of residues in the catalytic site of CgHMGR.

Investigating Candida glabrata Urinary Tract Infections (UTIs) in Mice Using Bioluminescence Imaging

Urinary tract infections (UTIs) are quite common and mainly caused by bacteria such as Escherichia coli. However, when patients have urinary catheters, fungal infections comprise up to 15% of these types of infections. Moreover, fungal UTIs have a high mortality, due to rapid spreading of the fungi to the kidneys. Most fungal UTIs are caused by Candida species, among which Candida albicans and Candida glabrata are the most common. C. glabrata is an opportunistic pathogenic yeast, phylogenetically quite close to Saccharomyces cerevisiae. Even though it is commonly isolated from the urinary tract and rapidly acquires resistance to antifungals, its pathogenesis has not been studied extensively in vivo. In vivo studies require high numbers of animals, which can be overcome by the use of non-invasive imaging tools. One such tool, bioluminescence imaging, has been used successfully to study different types of C. albicans infections. For C. glabrata, only biofilms on subcutaneously implanted catheters have been imaged using this tool. In this work, we investigated the progression of C. glabrata UTIs from the bladder to the kidneys and the spleen. Furthermore, we optimized expression of a red-shifted firefly luciferase in C. glabrata for in vivo use. We propose the first animal model using bioluminescence imaging to visualize C. glabrata in mouse tissues. Additionally, this UTI model can be used to monitor antifungal activity in vivo over time.

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.

The fungivorous amoeba Protostelium aurantium targets redox homeostasis and cell wall integrity during intracellular killing of Candida parapsilosis

Predatory interactions among microbes are major evolutionary driving forces for biodiversity. The fungivorous amoeba Protostelium aurantium has a wide fungal food spectrum including foremost pathogenic members of the genus Candida. Here we show that upon phagocytic ingestion by the amoeba, Candida parapsilosis is confronted with an oxidative burst and undergoes lysis within minutes of processing in acidified phagolysosomes. On the fungal side, a functional genomic approach identified copper and redox homeostasis as primary targets of amoeba predation, with the highly expressed copper exporter gene CRP1 and the peroxiredoxin gene PRX1 contributing to survival when encountered with P. aurantium. The fungicidal activity was largely retained in intracellular vesicles of the amoebae. Following their isolation, the content of these vesicles induced immediate killing and lysis of C. parapsilosis in vitro. Proteomic analysis identified 56 vesicular proteins from P. aurantium. Although completely unknown proteins were dominant, many of them could be categorised as hydrolytic enzymes targeting the fungal cell wall, indicating that fungal cell wall structures are under selection pressure by predatory phagocytes in natural environments. TAKE AWAY: The amoeba Protostelium aurantium feeds on fungi, such as Candida parapsilosis. Ingested yeast cells are exposed to reactive oxygen species. A copper exporter and a peroxiredoxin contribute to fungal defence. Yeast cells undergo intracellular lysis. Lysis occurs via a cocktail of hydrolytic enzymes from intracellular vesicles.

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.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Transient Mitochondria Dysfunction Confers Fungal Cross-Resistance against Phagocytic Killing and Fluconazole

Loss or inactivation of antivirulence genes is an adaptive strategy in pathogen evolution. Candida glabrata is an important opportunistic pathogen related to baker’s yeast, with the ability to both quickly increase its intrinsic high level of azole resistance and persist within phagocytes. During C. glabrata’s evolution as a pathogen, the mitochondrial DNA polymerase CgMip1 has been under positive selection. We show that CgMIP1 deletion not only triggers loss of mitochondrial function and a petite phenotype, but increases C. glabrata’s azole and endoplasmic reticulum (ER) stress resistance and, importantly, its survival in phagocytes. The same phenotype is induced by fluconazole and by exposure to macrophages, conferring a cross-resistance between antifungals and immune cells, and can be found in clinical isolates despite a slow growth of petite strains. This suggests that petite constitutes a bet-hedging strategy of C. glabrata and, potentially, a relevant cause of azole resistance. Mitochondrial function may therefore be considered a potential antivirulence factor.