Comparative genome analysis and gene finding in Candida species using CGOB

The Candida Genome Database: annotation and visualization updates

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

De novo SSR markers development and genetic diversity assessment in Ailanthus excelsa in India

Ailanthus excelsa is a fast-growing, multipurpose agroforestry tree species of Indian Arid Regions (IAR). It is widely cultivated as tree outside forests (TOFs) on farm lands, roadside, canal banks, etc., where the genetic stocks were randomly planted. To ensure the availability of quality planting materials (QPM) for industrial profitability, the germplasm must undergo a systematic genetic improvement program. Genetic variability in the base population is crucial for effective selection, but the lack of genomic resources and marker impedes this process. This study aimed to generate genome sequence information and de novo development of simple sequence repeats (SSRs) in A. excelsa. About 96 million raw reads were generated using Illumina platform, assembled into ~ 183,000 contigs with 33% GC content and an N50 value of 641 bp. A total of 7,667 microsatellite repeats were identified, with di-nucleotides being the most abundant. AT rich repeats were more prevalent than GC rich motifs. A total of 3,696 primer pairs were designed, and 150 of these were selected for validation. In PCR, 145 SSRs were positively amplified and 15 showed polymorphic banding pattern. These polymorphic SSRs were used to characterize 213 individuals from northern and central India. SSR analysis revealed high gene diversity (He = 0.71; Ar = 9.12) with negligible genetic differentiation in populations. The study presents a comprehensive set of de novo SSR markers and provides baseline knowledge of genetic structure of A. excelsa, essential for conservation and long-term genetic improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-025-01566-6.

When HSFs bring the heat-mapping the transcriptional circuitries of HSF-type regulators in Candida albicans

Heat shock factor (HSF)-type regulators are stress-responsive transcription factors widely distributed among eukaryotes, including fungi. They carry a four-stranded winged helix-turn-helix DNA-binding domain considered as the signature domain for HSFs. The genome of the opportunistic yeast Candida albicans encodes four HSF members, namely, Sfl1, Sfl2, Skn7, and the essential regulator, Hsf1. C. albicans HSFs do not only respond to heat shock and/or temperature variation but also to CO2 levels, oxidative stress, and quorum sensing, acting this way as central decision makers. In this minireview, I follow on the heels of my mSphere of Influence commentary (2020) to provide an overview of the repertoire of HSF regulators in Saccharomyces cerevisiae and C. albicans and describe how their genetic perturbation in C. albicans, coupled with genome-wide expression and location analyses, allow to map their transcriptional circuitry. I highlight how they can regulate, in common, a crucial developmental program: filamentous growth.

Regulation of Catalase Expression and Activity by DhHog1 in the Halotolerant Yeast Debaryomyces hansenii Under Saline and Oxidative Conditions

Efficient transcriptional regulation of the stress response is critical for microorganism survival. In yeast, stress-related gene expression, particularly for antioxidant enzymes like catalases, mitigates reactive oxygen species such as hydrogen peroxide (H2O2), preventing cell damage. The halotolerant yeast Debaryomyces hansenii shows oxidative stress tolerance, largely due to high catalase activity from DhCTA and DhCTT genes. This study evaluates D. hansenii's response to oxidative stress caused by H2O2 under saline conditions, focusing on cell viability, gene expression, and catalase activity. Chromatin organization in the promoter of DhCTA and DhCTT was analyzed, revealing low nucleosome occupancy in promoter regions, correlating with active gene expression. Stress-related motifs for transcription factors like Msn2/4 and Sko1 were found, suggesting regulation by the DhHog1 MAP kinase. Analysis of a Dhhog1Δ mutant showed DhHog1's role in DhCTA expression under H2O2 or NaCl conditions. These findings highlight DhHog1's critical role in regulating the stress response in D. hansenii, offering insights for enhancing stress tolerance in halotolerant yeasts, particularly for industrial applications in saline wastewater management.

Alternative sulphur metabolism in the fungal pathogen Candida parapsilosis

Candida parapsilosis is an opportunistic fungal pathogen commonly isolated from the environment and associated with nosocomial infection outbreaks worldwide. We describe here the construction of a large collection of gene disruptions, greatly increasing the molecular tools available for probing gene function in C. parapsilosis. We use these to identify transcription factors associated with multiple metabolic pathways, and in particular to dissect the network regulating the assimilation of sulphur. We find that, unlike in other yeasts and filamentous fungi, the transcription factor Met4 is not the main regulator of methionine synthesis. In C. parapsilosis, assimilation of inorganic sulphur (sulphate) and synthesis of cysteine and methionine is regulated by Met28, a paralog of Met4, whereas Met4 regulates expression of a wide array of transporters and enzymes involved in the assimilation of organosulfur compounds. Analysis of transcription factor binding sites suggests that Met4 is recruited by the DNA-binding protein Met32, and Met28 is recruited by Cbf1. Despite having different target genes, Met4 and Met28 have partial functional overlap, possibly because Met4 can contribute to assimilation of inorganic sulphur in the absence of Met28.

RFX transcription factor in the human-associated yeast Candida albicans regulates adhesion to oral epithelium

Adhesion to mucosal surfaces is a critical step in many bacterial and fungal infections. Here, using a mouse model of oral infection by the human fungal pathobiont Candida albicans, we report the identification of a novel regulator of C. albicans adhesion to the oral mucosa. The regulator is a member of the regulatory factor X (RFX) family of transcription factors, which control cellular processes ranging from genome integrity in model yeasts to tissue differentiation in vertebrates. Mice infected with the C. albicans rfx1 deletion mutant displayed increased fungal burden in tongues compared to animals infected with the reference strain. High-resolution imaging revealed RFX1 transcripts being expressed by C. albicans cells during infection. Concomitant with the increase in fungal burden, the rfx1 mutant elicited an enhanced innate immune response. Transcriptome analyses uncovered HWP1, a gene encoding an adhesion protein that mediates covalent attachment to buccal cells, as a major RFX1-regulated locus. Consistent with this result, we establish that C. albicans adhesion to oral cells is modulated by RFX1 in an HWP1-dependent manner. Our findings expand the repertoire of biological processes controlled by the RFX family and illustrate a mechanism whereby C. albicans can adjust adhesion to the oral epithelium.

A gain-of-function mutation in zinc cluster transcription factor Rob1 drives Candida albicans adaptive growth in the cystic fibrosis lung environment

Candida albicans chronically colonizes the respiratory tract of patients with Cystic Fibrosis (CF). It competes with CF-associated pathogens (e.g. Pseudomonas aeruginosa) and contributes to disease severity. We hypothesize that C. albicans undergoes specific adaptation mechanisms that explain its persistence in the CF lung environment. To identify the underlying genetic and phenotypic determinants, we serially recovered 146 C. albicans clinical isolates over a period of 30 months from the sputum of 25 antifungal-naive CF patients. Multilocus sequence typing analyses revealed that most patients were individually colonized with genetically close strains, facilitating comparative analyses between serial isolates. We strikingly observed differential ability to filament and form monospecies and dual-species biofilms with P. aeruginosa among 18 serial isolates sharing the same diploid sequence type, recovered within one year from a pediatric patient. Whole genome sequencing revealed that their genomes were highly heterozygous and similar to each other, displaying a highly clonal subpopulation structure. Data mining identified 34 non-synonymous heterozygous SNPs in 19 open reading frames differentiating the hyperfilamentous and strong biofilm-former strains from the remaining isolates. Among these, we detected a glycine-to-glutamate substitution at position 299 (G299E) in the deduced amino acid sequence of the zinc cluster transcription factor ROB1 (ROB1G299E), encoding a major regulator of filamentous growth and biofilm formation. Introduction of the G299E heterozygous mutation in a co-isolated weak biofilm-former CF strain was sufficient to confer hyperfilamentous growth, increased expression of hyphal-specific genes, increased monospecies biofilm formation and increased survival in dual-species biofilms formed with P. aeruginosa, indicating that ROB1G299E is a gain-of-function mutation. Disruption of ROB1 in a hyperfilamentous isolate carrying the ROB1G299E allele abolished hyperfilamentation and biofilm formation. Our study links a single heterozygous mutation to the ability of C. albicans to better survive during the interaction with other CF-associated microbes and illuminates how adaptive traits emerge in microbial pathogens to persistently colonize and/or infect the CF-patient airways.

Candida auris-macrophage cellular interactions and transcriptional response

The pathogenic yeast Candida auris represents a global threat of the utmost clinical relevance. This emerging fungal species is remarkable in its resistance to commonly used antifungal agents and its persistence in the nosocomial settings. The innate immune system is one the first lines of defense preventing the dissemination of pathogens in the host. C. auris is susceptible to circulating phagocytes, and understanding the molecular details of these interactions may suggest routes to improved therapies. In this work, we examined the interactions of this yeast with macrophages. We found that macrophages avidly phagocytose C. auris; however, intracellular replication is not inhibited, indicating that C. auris resists the killing mechanisms imposed by the phagocyte. Unlike Candida albicans, phagocytosis of C. auris does not induce macrophage lysis. The transcriptional response of C. auris to macrophage phagocytosis is very similar to other members of the CUG clade (C. albicans, C. tropicalis, C. parapsilosis, C. lusitaniae), i.e., downregulation of transcription/translation and upregulation of alternative carbon metabolism pathways, transporters, and induction of oxidative stress response and proteolysis. Gene family expansions are common in this yeast, and we found that many of these genes are induced in response to macrophage co-incubation. Among these, amino acid and oligopeptide transporters, as well as lipases and proteases, are upregulated. Thus, C. auris shares key transcriptional signatures shared with other fungal pathogens and capitalizes on the expansion of gene families coding for potential virulence attributes that allow its survival, persistence, and evasion of the innate immune system.

A Candida auris-specific adhesin, Scf1, governs surface association, colonization, and virulence

Candida auris is an emerging fungal pathogen responsible for health care-associated outbreaks that arise from persistent surface and skin colonization. We characterized the arsenal of adhesins used by C. auris and discovered an uncharacterized adhesin, Surface Colonization Factor (Scf1), and a conserved adhesin, Iff4109, that are essential for the colonization of inert surfaces and mammalian hosts. SCF1 is apparently specific to C. auris, and its expression mediates adhesion to inert and biological surfaces across isolates from all five clades. Unlike canonical fungal adhesins, which function through hydrophobic interactions, Scf1 relies on exposed cationic residues for surface association. SCF1 is required for C. auris biofilm formation, skin colonization, virulence in systemic infection, and colonization of inserted medical devices.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

Genome survey sequencing-based SSR marker development and their validation in Dendrocalamus longispathus

Bamboo is an important genetic resource in India, supporting rural livelihood and industries. Unfortunately, most Indian bamboo taxa are devoid of basic genomic or marker information required to comprehend the genetic processes for further conservation and management. In this study, we perform genome survey sequencing for development of de novo genomic SSRs in Dendrocalamus longispathus, a socioeconomically important bamboo species of northeast India. Using Illumina platform, 69.49 million raw reads were generated and assembled into 1,145,321 contig with GC content 43% and N50 1228 bp. In total, 46,984 microsatellite repeats were mined-out wherein di-nucleotide repeats were most abundant (54.71%) followed by mono- (31.91%) and tri-repeats (9.85%). Overall, AT-rich repeats were predominant in the genome, but GC-rich motifs were more frequent in tri-repeats. Afterwards, 21,596 SSR loci were successfully tagged with the primer pairs, and a subset of 50 were validated through polymerase chain reaction amplification. Of these, 36 SSR loci were successfully amplified, and 16 demonstrated polymorphism. Using 13 polymorphic SSRs, a moderate level of gene diversity (He = 0.480; Ar = 3.52) was recorded in the analysed populations of D. longispathus. Despite the high gene flow (Nm = 4.928) and low genetic differentiation (F_ST = 0.119), severe inbreeding (F_IS = 0.407) was detected. Further, genetic clustering and STRUCTURE analysis revealed that the entire genetic variability is captured under two major gene pools. Conclusively, we present a comprehensive set of novel SSR markers in D. longispathus as well as other taxa of tropical woody bamboos.

An Approximate Bayesian Computation Approach for Modeling Genome Rearrangements

The inference of genome rearrangement events has been extensively studied, as they play a major role in molecular evolution. However, probabilistic evolutionary models that explicitly imitate the evolutionary dynamics of such events, as well as methods to infer model parameters, are yet to be fully utilized. Here, we developed a probabilistic approach to infer genome rearrangement rate parameters using an Approximate Bayesian Computation (ABC) framework. We developed two genome rearrangement models, a basic model, which accounts for genomic changes in gene order, and a more sophisticated one which also accounts for changes in chromosome number. We characterized the ABC inference accuracy using simulations and applied our methodology to both prokaryotic and eukaryotic empirical datasets. Knowledge of genome-rearrangement rates can help elucidate their role in evolution as well as help simulate genomes with evolutionary dynamics that reflect empirical genomes.

The DNA replication protein Orc1 from the yeast Torulaspora delbrueckii is required for heterochromatin formation but not as a silencer-binding protein

To understand the process by which new protein functions emerge, we examined how the yeast heterochromatin protein Sir3 arose through gene duplication from the conserved DNA replication protein Orc1. Orc1 is a subunit of the origin recognition complex (ORC), which marks origins of DNA replication. In Saccharomyces cerevisiae, Orc1 also promotes heterochromatin assembly by recruiting the structural proteins Sir1-4 to silencer DNA. In contrast, the paralog of Orc1, Sir3, is a nucleosome-binding protein that spreads across heterochromatic loci in conjunction with other Sir proteins. We previously found that a non-duplicated Orc1 from the yeast Kluyveromyces lactis behaved like ScSir3 but did not have a silencer-binding function like ScOrc1. Moreover, K. lactis lacks Sir1, the protein that interacts directly with ScOrc1 at the silencer. Here, we examined whether the emergence of Sir1 coincided with Orc1 acting as a silencer-binding protein. In the non-duplicated species Torulaspora delbrueckii, which has an ortholog of Sir1 (TdKos3), we found that TdOrc1 spreads across heterochromatic loci independently of ORC, as ScSir3 and KlOrc1 do. This spreading is dependent on the nucleosome binding BAH domain of Orc1 and on Sir2 and Kos3. However, TdOrc1 does not have a silencer-binding function: T. delbrueckii silencers do not require ORC binding sites to function, and Orc1 and Kos3 do not appear to interact. Instead, Orc1 and Kos3 both spread across heterochromatic loci with other Sir proteins. Thus, Orc1 and Sir1/Kos3 originally had different roles in heterochromatin formation than they do now in S. cerevisiae.

Development and validation of monoclonal antibodies specific for Candida albicans Als2, Als9-1, and Als9-2

Fungal agglutinin-like sequence (Als) cell-surface glycoproteins, best characterized in Candida albicans, mediate adhesive and aggregative interactions with host cells, other microbes, and abiotic surfaces. Monoclonal antibodies (MAbs) specific for each C. albicans Als protein are valuable reagents for gaining insight into Als protein localization and function. This manuscript describes development and validation of MAbs specific for C. albicans Als2, as well as for C. albicans Als9-1 and Als9-2, two protein variants produced from the ALS9 locus. Native C. albicans ALS9 expression levels were not sufficiently high to produce detectable Als9 protein on the wild-type cell surface so MAb validation required production of overexpression strains, each featuring one of the two ALS9 alleles. An anti-Als2 MAb was raised against an N-glycosylated form of the protein immunogen, as well as an Endoglycosidase H-treated immunogen. The MAb raised against the N-glycosylated immunogen proved superior and immunolabeled C. albicans yeast cells and germ tubes, and the surface of Candida dubliniensis and Candida tropicalis yeasts. Als2 was visible on C. albicans yeast cells recovered from a murine model of oral candidiasis, demonstrating Als2 production both in vivo and in vitro. These new MAbs add to the collection of anti-Als MAbs that are powerful tools to better understand the role of Als proteins in C. albicans biology and pathogenesis.

Development and Use of a Monoclonal Antibody Specific for the Candida albicans Cell-Surface Protein Hwp1

The Candida albicans cell-surface protein Hwp1 functions in adhesion to the host and in biofilm formation. A peptide from the Gln-Pro-rich adhesive domain of Hwp1 was used to raise monoclonal antibody (MAb) 2-E8. MAb 2-E8 specificity for Hwp1 was demonstrated using a hwp1/hwp1 C. albicans isolate and strains that expressed at least one HWP1 allele. Immunofluorescence and atomic force microscopy experiments using MAb 2-E8 confirmed C. albicans germ-tube-specific detection of the Hwp1 protein. MAb 2-E8 also immunolabeled the tips of some Candida dubliniensis germ tubes grown under conditions that maximized HWP1 expression. The phylogeny of HWP1 and closely related genes suggested that the Gln-Pro-rich adhesive domain was unique to C. albicans and C. dubliniensis focusing the utility of MAb 2-E8 on these species. This new reagent can be used to address unanswered questions about Hwp1 and its interactions with other proteins in the context of C. albicans biology and pathogenesis.

Peering Into Candida albicans Pir Protein Function and Comparative Genomics of the Pir Family

The fungal cell wall, comprised primarily of protein and polymeric carbohydrate, maintains cell structure, provides protection from the environment, and is an important antifungal drug target. Pir proteins (proteins with internal repeats) are linked to cell wall β-1,3-glucan and are best studied in Saccharomyces cerevisiae. Sequential deletion of S. cerevisiae PIR genes produces strains with increasingly notable cell wall damage. However, a true null mutant lacking all five S. cerevisiae PIR genes was never constructed. Because only two PIR genes (PIR1, PIR32) were annotated in the Candida albicans genome, the initial goal of this work was to construct a true Δpir/Δpir null strain in this species. Unexpectedly, the phenotype of the null strain was almost indistinguishable from its parent, leading to the search for other proteins with Pir function. Bioinformatic approaches revealed nine additional C. albicans proteins that share a conserved Pir functional motif (minimally DGQ). Examination of the protein sequences revealed another conserved motif (QFQFD) toward the C-terminal end of each protein. Sequence similarities and presence of the conserved motif(s) were used to identify a set of 75 proteins across 16 fungal species that are proposed here as Pir proteins. The Pir family is greatly expanded in C. albicans and C. dubliniensis compared to other species and the orthologs are known to have specialized function during chlamydospore formation. Predicted Pir structures showed a conserved core of antiparallel beta-sheets and sometimes-extensive loops that contain amino acids with the potential to form linkages to cell wall components. Pir phylogeny demonstrated emergence of specific ortholog groups among the fungal species. Variation in gene expression patterns was noted among the ortholog groups during growth in rich medium. PIR allelic variation was quite limited despite the presence of a repeated sequence in many loci. Results presented here demonstrate that the Pir family is larger than previously recognized and lead to new hypotheses to test to better understand Pir proteins and their role in the fungal cell wall.

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.

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.

How to Use the Candida Genome Database

The Candida Genome Database provides access to biological information about genes and proteins of several medically important Candida species. The website is organized into easily navigable pages that enable data retrieval and analysis. This chapter shows how to explore the CGD Home page and Locus Summary pages, which are the main access points to the database. It also provides a description of how to use the GO analysis tools, GO Term Finder, and GO Slim Mapper and how to browse large-scale datasets using the JBrowse genome browser. Finally, it shows how to search and retrieve data for user-defined sets of genes using the Advanced Search and Batch Download tools.

Plasmid-Based CRISPR-Cas9 Editing in Multiple Candida Species

CRISPR-Cas9 technology radically changed the approach to genetic manipulation of both medically and industrially relevant Candida species, as attested by the ever-increasing number of applications to the study of pathogenesis, drug resistance, gene expression, and host pathogen interaction and drug discovery. Here, we describe the use of plasmid-based systems for high efficiency CRISPR-Cas9 gene editing into any strain of four non-albicans Candida species, namely, Candida parapsilosis, Candida orthopsilosis, Candida metapsilosis, and Candida tropicalis. The plasmids pCP-tRNA and pCT-tRNA contain all the elements necessary for expressing the CRISPR-Cas9 machinery, and they can be used in combination with a repair template for disrupting gene function by insertion of a premature stop codon or by gene deletion. The plasmids are easily lost in the absence of selection, allowing scarless gene editing and minimizing detrimental effects of prolonged Cas9 expression.

Interactions of Both Pathogenic and Nonpathogenic CUG Clade Candida Species with Macrophages Share a Conserved Transcriptional Landscape

Candida species are a leading cause of opportunistic, hospital-associated bloodstream infections with high mortality rates, typically in immunocompromised patients. Several species, including Candida albicans, the most prevalent cause of infection, belong to the monophyletic CUG clade of yeasts. Innate immune cells such as macrophages are crucial for controlling infection, and C. albicans responds to phagocytosis by a coordinated induction of pathways involved in catabolism of nonglucose carbon sources, termed alternative carbon metabolism, which together are essential for virulence. However, the interactions of other CUG clade species with macrophages have not been characterized. Here, we analyzed transcriptional responses to macrophage phagocytosis by six Candida species across a range of virulence and clinical importance. We define a core induced response common to pathogenic and nonpathogenic species alike, heavily weighted to alternative carbon metabolism. One prominent pathogen, Candida parapsilosis, showed species-specific expansion of phagocytosis-responsive genes, particularly metabolite transporters. C. albicans and Candida tropicalis, the other prominent pathogens, also had species-specific responses, but these were largely comprised of functionally uncharacterized genes. Transcriptional analysis of macrophages also demonstrated highly correlated proinflammatory transcriptional responses to different Candida species that were largely independent of fungal viability, suggesting that this response is driven by recognition of conserved cell wall components. This study significantly broadens our understanding of host interactions in CUG clade species, demonstrating that although metabolic plasticity is crucial for virulence in Candida, it alone is not sufficient to confer pathogenicity. Instead, we identify sets of mostly uncharacterized genes that may explain the evolution of pathogenicity.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Using Genomics to Shape the Definition of the Agglutinin-Like Sequence (ALS) Family in the Saccharomycetales

The Candida albicans agglutinin-like sequence (ALS) family is studied because of its contribution to cell adhesion, fungal colonization, and polymicrobial biofilm formation. The goal of this work was to derive an accurate census and sequence for ALS genes in pathogenic yeasts and other closely related species, while probing the boundaries of the ALS family within the Order Saccharomycetales. Bioinformatic methods were combined with laboratory experimentation to characterize 47 novel ALS loci from 8 fungal species. AlphaFold predictions suggested the presence of a conserved N-terminal adhesive domain (NT-Als) structure in all Als proteins reported to date, as well as in S. cerevisiae alpha-agglutinin (Sag1). Lodderomyces elongisporus, Meyerozyma guilliermondii, and Scheffersomyces stipitis were notable because each species had genes with C. albicans ALS features, as well as at least one that encoded a Sag1-like protein. Detection of recombination events between the ALS family and gene families encoding other cell-surface proteins such as Iff/Hyr and Flo suggest widespread domain swapping with the potential to create cell-surface diversity among yeast species. Results from the analysis also revealed subtelomeric ALS genes, ALS pseudogenes, and the potential for yeast species to secrete their own soluble adhesion inhibitors. Information presented here supports the inclusion of SAG1 in the ALS family and yields many experimental hypotheses to pursue to further reveal the nature of the ALS family.

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 H₂O₂, 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.

Opportunistic yeast pathogen Candida spp.: Secreted and membrane-bound virulence factors

Candidiasis is a fungal infection caused by Candida spp. especially Candida albicans, C. glabrata, C. parapsilosis and C. tropicalis. Although the medicinal therapeutic strategies have rapidly improved, the mortality rate due to candidiasis has continuously increased. The secreted and membrane-bound virulence factors (VFs) are responsible for fungal invasion, damage and translocation through the host enterocytes besides the evasion from host immune system. VFs such as agglutinin-like sequences (Als), heat shock protein 70, phospholipases, secreted aspartyl proteinases (Sap), lipases, enolases and phytases are mostly hydrolases which degrade the enterocyte membrane components except for candidalysin, the VF acts as a peptide toxin to induce necrotic cell lysis. To date, structural studies of the VFs remain underexplored, hindering their functional analyses. Among the VFs, only secreted aspartyl proteinases and agglutinin-like sequences have their structures deposited in Protein Data Bank (PDB). Therefore, this review scrutinizes the mechanisms of these VFs by discussing the VF-deficient studies of several Candida spp. and their abilities to produce these VFs. Nonetheless, their latest reported sequential and structural analyses are discussed to impart a wider perception of the host-pathogen interactions and potential vaccine or antifungal drug targets. This review signifies that more VFs structural investigations and mining in the emerging Candida spp. are required to decipher their pathogenicity and virulence mechanisms compared to the prominent C. albicans.

LAY ABSTRACT

Candida virulence factors (VFs) including mainly enzymes and proteins play vital roles in breaching the human intestinal barrier and causing deadly candidiasis. Limited VFs' structural studies hinder deeper comprehension of their mechanisms and thus the design of vaccines and antifungal drugs against fungal infections.

Conservation of a DNA Replication Motif among Phylogenetically Distant Budding Yeast Species

Eukaryotic DNA replication begins at genomic loci termed origins, which are bound by the origin recognition complex (ORC). Although ORC is conserved across species, the sequence composition of origins is more varied. In the budding yeast Saccharomyces cerevisiae, the ORC-binding motif consists of an A/T-rich 17 bp “extended ACS” sequence adjacent to a B1 element composed of two 3-bp motifs. Similar sequences occur at origins in closely related species, but it is not clear when this type of replication origin arose and whether it predated a whole-genome duplication that occurred around 100 Ma in the budding yeast lineage. To address these questions, we identified the ORC-binding sequences in the nonduplicated species Torulaspora delbrueckii. We used chromatin immunoprecipitation followed by sequencing and identified 190 ORC-binding sites distributed across the eight T. delbrueckii chromosomes. Using these sites, we identified an ORC-binding motif that is nearly identical to the known motif in S. cerevisiae. We also found that the T. delbrueckii ORC-binding sites function as origins in T. delbrueckii when cloned onto a plasmid and that the motif is required for plasmid replication. Finally, we compared an S. cerevisiae origin with two T. delbrueckii ORC-binding sites and found that they conferred similar stabilities to a plasmid. These results reveal that the ORC-binding motif arose prior to the whole-genome duplication and has been maintained for over 100 Myr.

Genetic Manipulation as a Tool to Unravel Candida parapsilosis Species Complex Virulence and Drug Resistance: State of the Art

An increase in the rate of isolation of Candida parapsilosis in the past decade, as well as increased identification of azole-resistant strains are concerning, and require better understanding of virulence-like factors and drug-resistant traits of these species. In this regard, the present review “draws a line” on the information acquired, thus far, on virulence determinants and molecular mechanisms of antifungal resistance in these opportunistic pathogens, mainly derived from genetic manipulation studies. This will provide better focus on where we stand in our understanding of the C. parapsilosis species complex–host interaction, and how far we are from defining potential novel targets or therapeutic strategies—key factors to pave the way for a more tailored management of fungal infections caused by these fungal pathogens.

Candida pathogens induce protective mitochondria-associated type I interferon signalling and a damage-driven response in vaginal epithelial cells

Vaginal candidiasis is an extremely common disease predominantly caused by four phylogenetically diverse species: Candida albicans; Candida glabrata; Candida parapsilosis; and Candida tropicalis. Using a time course infection model of vaginal epithelial cells and dual RNA sequencing, we show that these species exhibit distinct pathogenicity patterns, which are defined by highly species-specific transcriptional profiles during infection of vaginal epithelial cells. In contrast, host cells exhibit a homogeneous response to all species at the early stages of infection, which is characterized by sublethal mitochondrial signalling inducing a protective type I interferon response. At the later stages, the transcriptional response of the host diverges in a species-dependent manner. This divergence is primarily driven by the extent of epithelial damage elicited by species-specific mechanisms, such as secretion of the toxin candidalysin by C. albicans. Our results uncover a dynamic, biphasic response of vaginal epithelial cells to Candida species, which is characterized by protective mitochondria-associated type I interferon signalling and a species-specific damage-driven response.

Evolution of the complex transcription network controlling biofilm formation in Candida species

We examine how a complex transcription network composed of seven 'master' regulators and hundreds of target genes evolved over a span of approximately 70 million years. The network controls biofilm formation in several Candida species, a group of fungi that are present in humans both as constituents of the microbiota and as opportunistic pathogens. Using a variety of approaches, we observed two major types of changes that have occurred in the biofilm network since the four extant species we examined last shared a common ancestor. Master regulator 'substitutions' occurred over relatively long evolutionary times, resulting in different species having overlapping but different sets of master regulators of biofilm formation. Second, massive changes in the connections between the master regulators and their target genes occurred over much shorter timescales. We believe this analysis is the first detailed, empirical description of how a complex transcription network has evolved.

Transcription Profiles Associated with Inducible Adhesion in Candida parapsilosis

Candida parapsilosis has emerged as a frequent cause of invasive candidiasis with increasing evidence of unique biological features relative to C. albicans As it adapts to conditions within a mammalian host, rapid changes in gene expression are necessary to facilitate colonization and persistence in this environment. Adhesion of the organism to biological surfaces is a key first step in this process and is the focus of this study. Building on previous observations showing the importance of a member of the ALS gene family in C. parapsilosis adhesion, three clinical isolates were cultured under two conditions that mimic the mammalian host and promote adhesion, incubation at 37°C in tissue culture medium 199 or in human plasma. Transcriptional profiles using RNA-seq were obtained in these adhesion-inducing conditions and compared to profiles following growth in yeast media that suppress adhesion to identify gene expression profiles associated with adhesion. Overall gene expression profiles among the three strains were similar in both adhesion-inducing conditions and distinct from adhesion-suppressing conditions. Pairwise analysis among the three growth conditions identified 133 genes that were differentially expressed at a cutoff of ±4-fold, with the most upregulated genes significantly enriched in iron acquisition and transmembrane transport, while the most downregulated genes were enriched in oxidation-reduction processes. Gene family enrichment analysis identified gene families with diverse functions that may have an important role in this important step for colonization and disease.IMPORTANCE Invasive Candida infections are frequent complications of the immunocompromised and are associated with substantive morbidity and mortality. Although C. albicans is the best-studied species, emerging infections by non-albicans Candida species have led to increased efforts to understand aspects of their pathogenesis that are unique from C. albicans C. parapsilosis is a frequent cause of invasive infections, particularly among premature infants. Recent efforts have identified important virulence mechanisms that have features distinct from C. albicans C. parapsilosis can exist outside a host environment and therefore requires rapid modifications when it encounters a mammalian host to prevent its clearance. An important first step in the process is adhesion to host surfaces. This work takes a global, nonbiased approach to investigate broad changes in gene expression that accompany efficient adhesion. As such, biological pathways and individual protein targets are identified that may be amenable to manipulation to reduce colonization and disease from this organism.

Systematic Analysis of Functionally Related Gene Clusters in the Opportunistic Pathogen, Candida albicans

The proper balance of gene expression is essential for cellular health, organismal development, and maintaining homeostasis. In response to complex internal and external signals, the cell needs to modulate gene expression to maintain proteostasis and establish cellular identity within its niche. On a genome level, single-celled prokaryotic microbes display clustering of co-expressed genes that are regulated as a polycistronic RNA. This phenomenon is largely absent from eukaryotic microbes, although there is extensive clustering of co-expressed genes as functional pairs spread throughout the genome in Saccharomyces cerevisiae. While initial analysis demonstrated conservation of clustering in divergent fungal lineages, a comprehensive analysis has yet to be performed. Here we report on the prevalence, conservation, and significance of the functional clustering of co-regulated genes within the opportunistic human pathogen, Candida albicans. Our analysis reveals that there is extensive clustering within this organism-although the identity of the gene pairs is unique compared with those found in S. cerevisiae-indicating that this genomic arrangement evolved after these microbes diverged evolutionarily, rather than being the result of an ancestral arrangement. We report a clustered arrangement in gene families that participate in diverse molecular functions and are not the result of a divergent orientation with a shared promoter. This arrangement coordinates the transcription of the clustered genes to their neighboring genes, with the clusters congregating to genomic loci that are conducive to transcriptional regulation at a distance.

Pursuing Advances in DNA Sequencing Technology to Solve a Complex Genomic Jigsaw Puzzle: The Agglutinin-Like Sequence (ALS) Genes of Candida tropicalis

The agglutinin-like sequence (ALS) gene family encodes cell-surface adhesins that interact with host and abiotic surfaces, promoting colonization by opportunistic fungal pathogens such as Candida tropicalis. Studies of Als protein contribution to C. tropicalis adhesion would benefit from an accurate catalog of ALS gene sequences as well as insight into relative gene expression levels. Even in the genomics era, this information has been elusive: genome assemblies are often broken within ALS genes because of their extensive regions of highly conserved, repeated DNA sequences and because there are many similar ALS genes at different chromosomal locations. Here, we describe the benefit of long-read DNA sequencing technology to facilitate characterization of C. tropicalis ALS loci. Thirteen ALS loci in C. tropicalis strain MYA-3404 were deduced from a genome assembly constructed from Illumina MiSeq and Oxford Nanopore MinION data. Although the MinION data were valuable, PCR amplification and Sanger sequencing of ALS loci were still required to complete and verify the gene sequences. Each predicted Als protein featured an N-terminal binding domain, a central domain of tandemly repeated sequences, and a C-terminal domain rich in Ser and Thr. The presence of a secretory signal peptide and consensus sequence for addition of a glycosylphosphatidylinositol (GPI) anchor was consistent with predicted protein localization to the cell surface. TaqMan assays were designed to recognize each ALS gene, as well as both alleles at the divergent CtrALS3882 locus. C. tropicalis cells grown in five different in vitro conditions showed differential expression of various ALS genes. To place the C. tropicalis data into a larger context, TaqMan assays were also designed and validated for analysis of ALS gene expression in Candida albicans and Candida dubliniensis. These comparisons identified the subset of highly expressed C. tropicalis ALS genes that were predicted to encode proteins with the most abundant cell-surface presence, prioritizing them for subsequent functional analysis. Data presented here provide a solid foundation for future experimentation to deduce ALS family contributions to C. tropicalis adhesion and pathogenesis.

Pentose metabolism and conversion to biofuels and high-value chemicals in yeasts

Pentose sugars are widespread in nature and two of them, D-xylose and L-arabinose belong to the most abundant sugars being the second and third by abundance sugars in dry plant biomass (lignocellulose) and in general on planet. Therefore, it is not surprising that metabolism and bioconversion of these pentoses attract much attention. Several different pathways of D-xylose and L-arabinose catabolism in bacteria and yeasts are known. There are even more common and really ubiquitous though not so abundant pentoses, D-ribose and 2-deoxy-D-ribose, the constituents of all living cells. Thus, ribose metabolism is example of endogenous metabolism whereas metabolism of other pentoses, including xylose and L-arabinose, represents examples of the metabolism of foreign exogenous compounds which normally are not constituents of yeast cells. As a rule, pentose degradation by the wild-type strains of microorganisms does not lead to accumulation of high amounts of valuable substances; however, productive strains have been obtained by random selection and metabolic engineering. There are numerous reviews on xylose and (less) L-arabinose metabolism and conversion to high value substances; however, they mostly are devoted to bacteria or the yeast Saccharomyces cerevisiae. This review is devoted to reviewing pentose metabolism and bioconversion mostly in non-conventional yeasts, which naturally metabolize xylose. Pentose metabolism in the recombinant strains of S. cerevisiae is also considered for comparison. The available data on ribose, xylose, L-arabinose transport, metabolism, regulation of these processes, interaction with glucose catabolism and construction of the productive strains of high-value chemicals or pentose (ribose) itself are described. In addition, genome studies of the natural xylose metabolizing yeasts and available tools for their molecular research are reviewed. Metabolism of other pentoses (2-deoxyribose, D-arabinose, lyxose) is briefly reviewed.

Parallel Concerted Evolution of Ribosomal Protein Genes in Fungi and Its Adaptive Significance

Ribosomal protein (RP) genes encode structural components of ribosomes, the cellular machinery for protein synthesis. A single functional copy has been maintained in most of 78–80 RP families in animals due to evolutionary constraints imposed by gene dosage balance. Some fungal species have maintained duplicate copies in most RP families. The mechanisms by which the RP genes were duplicated and maintained and their functional significance are poorly understood. To address these questions, we identified all RP genes from 295 fungi and inferred the timing and nature of gene duplication events for all RP families. We found that massive duplications of RP genes have independently occurred by different mechanisms in three distantly related lineages: budding yeasts, fission yeasts, and Mucoromycota. The RP gene duplicates in budding yeasts and Mucoromycota were mainly created by whole genome duplication events. However, duplicate RP genes in fission yeasts were likely generated by retroposition, which is unexpected considering their dosage sensitivity. The sequences of most RP paralogs have been homogenized by repeated gene conversion in each species, demonstrating parallel concerted evolution, which might have facilitated the retention of their duplicates. Transcriptomic data suggest that the duplication and retention of RP genes increased their transcript abundance. Physiological data indicate that increased ribosome biogenesis allowed these organisms to rapidly consume sugars through fermentation while maintaining high growth rates, providing selective advantages to these species in sugar-rich environments.

The Methylotroph Gene Order Browser (MGOB) reveals conserved synteny and ancestral centromere locations in the yeast family Pichiaceae

The yeast family Pichiaceae, also known as the 'methylotrophs clade', is a relatively little studied group of yeasts despite its economic and clinical relevance. To explore the genome evolution and synteny relationships within this family, we developed the Methylotroph Gene Order Browser (MGOB, http://mgob.ucd.ie) similar to our previous gene order browsers for other yeast families. The dataset contains genome sequences from nine Pichiaceae species, including our recent reference sequence of Pichia kudriavzevii. As an example, we demonstrate the conservation of synteny around the MOX1 locus among species both containing and lacking the MOX1 gene for methanol assimilation. We found ancient clusters of genes that are conserved as adjacent between Pichiaceae and Saccharomycetaceae. Surprisingly, we found evidence that the locations of some centromeres have been conserved among Pichiaceae species, and between Pichiaceae and Saccharomycetaceae, even though the centromeres fall into different structural categories-point centromeres, inverted repeats and retrotransposon cluster centromeres.

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.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

Biological and genomic analyses of a clinical isolate of Yarrowia galli from China

Infections caused by emerging fungal pathogens represent a new threat to human health. The yeast Yarrowia (Candida) galli was first described from chicken breast and liver in 2004 and has occasionally been isolated in clinical settings. In this study, we present the first report of a Y. galli isolate from a face granuloma of a woman. Y. galli is unable to grow at human physiological temperature (37 °C). Phenotypic analysis demonstrates that Y. galli can exist as several morphological types, namely fluffy, sticky, tight, and yeast forms, based on their cellular and colony appearances. Interestingly, Y. galli is able to undergo switching among different morphologies. These morphological changes are similar to the switching systems in pathogenic Candida species such as Candida albicans and Candida tropicalis. We further sequenced the genome of the Y. galli isolate. A comparative analysis with pathogenic yeast species indicated that a set of lipid metabolism genes were enriched in Y. galli. Domain enrichment analysis demonstrated that, similar to Candida clade species, the genome of Y. galli maintained several gene families required for virulence. Our biological and genomic analyses provide new insights into the understanding of the biology of Y. galli as either an environmental isolate or a potential human pathogen.

Correlating Genotype and Phenotype in the Asexual Yeast Candida orthopsilosis Implicates ZCF29 in Sensitivity to Caffeine

Candida orthopsilosis is diploid asexual yeast that causes human disease. Most C. orthopsilosis isolates arose from at least four separate hybridizations between related, but not identical, parents. Here, we used population genomics data to correlate genotypic and phenotypic variation in 28 C. orthopsilosis isolates. We used cosine similarity scores to identify 65 variants with potential high-impact (deleterious effects) that correlated with specific phenotypes. Of these, 19 were Single Nucleotide Polymorphisms (SNPs) that changed stop or start codons, or splice sites. One variant resulted in a premature stop codon in both alleles of the gene ZCF29 in C. orthopsilosis isolate 185, which correlated with sensitivity to nystatin and caffeine. We used CRISPR-Cas9 editing to introduce this polymorphism into two resistant C. orthopsilosis isolates. Introducing the stop codon resulted in sensitivity to caffeine and to ketoconazole, but not to nystatin. Our analysis shows that it is possible to associate genomic variants with phenotype in asexual Candida species, but that only a small amount of genomic variation can be easily explored.

Recent trends in molecular diagnostics of yeast infections: from PCR to NGS

The incidence of opportunistic yeast infections in humans has been increasing over recent years. These infections are difficult to treat and diagnose, in part due to the large number and broad diversity of species that can underlie the infection. In addition, resistance to one or several antifungal drugs in infecting strains is increasingly being reported, severely limiting therapeutic options and showcasing the need for rapid detection of the infecting agent and its drug susceptibility profile. Current methods for species and resistance identification lack satisfactory sensitivity and specificity, and often require prior culturing of the infecting agent, which delays diagnosis. Recently developed high-throughput technologies such as next generation sequencing or proteomics are opening completely new avenues for more sensitive, accurate and fast diagnosis of yeast pathogens. These approaches are the focus of intensive research, but translation into the clinics requires overcoming important challenges. In this review, we provide an overview of existing and recently emerged approaches that can be used in the identification of yeast pathogens and their drug resistance profiles. Throughout the text we highlight the advantages and disadvantages of each methodology and discuss the most promising developments in their path from bench to bedside.

Conservation of location of several specific inhibitory codon pairs in the Saccharomyces sensu stricto yeasts reveals translational selection

Synonymous codons provide redundancy in the genetic code that influences translation rates in many organisms, in which overall codon use is driven by selection for optimal codons. It is unresolved if or to what extent translational selection drives use of suboptimal codons or codon pairs. In Saccharomyces cerevisiae, 17 specific inhibitory codon pairs, each comprised of adjacent suboptimal codons, inhibit translation efficiency in a manner distinct from their constituent codons, and many are translated slowly in native genes. We show here that selection operates within Saccharomyces sensu stricto yeasts to conserve nine of these codon pairs at defined positions in genes. Conservation of these inhibitory codon pairs is significantly greater than expected, relative to conservation of their constituent codons, with seven pairs more highly conserved than any other synonymous pair. Conservation is strongly correlated with slow translation of the pairs. Conservation of suboptimal codon pairs extends to two related Candida species, fungi that diverged from Saccharomyces ∼270 million years ago, with an enrichment for codons decoded by I•A and U•G wobble in both Candida and Saccharomyces. Thus, conservation of inhibitory codon pairs strongly implies selection for slow translation at particular gene locations, executed by suboptimal codon pairs.

AYbRAH: a curated ortholog database for yeasts and fungi spanning 600 million years of evolution

Budding yeasts inhabit a range of environments by exploiting various metabolic traits. The genetic bases for these traits are mostly unknown, preventing their addition or removal in a chassis organism for metabolic engineering. Insight into the evolution of orthologs, paralogs and xenologs in the yeast pan-genome can help bridge these genotypes; however, existing phylogenomic databases do not span diverse yeasts, and sometimes cannot distinguish between these homologs. To help understand the molecular evolution of these traits in yeasts, we created Analyzing Yeasts by Reconstructing Ancestry of Homologs (AYbRAH), an open-source database of predicted and manually curated ortholog groups for 33 diverse fungi and yeasts in Dikarya, spanning 600 million years of evolution. OrthoMCL and OrthoDB were used to cluster protein sequence into ortholog and homolog groups, respectively; MAFFT and PhyML reconstructed the phylogeny of all homolog groups. Ortholog assignments for enzymes and small metabolite transporters were compared to their phylogenetic reconstruction, and curated to resolve any discrepancies. Information on homolog and ortholog groups can be viewed in the AYbRAH web portal (https://lmse.github.io/aybrah/), including functional annotations, predictions for mitochondrial localization and transmembrane domains, literature references and phylogenetic reconstructions. Ortholog assignments in AYbRAH were compared to HOGENOM, KEGG Orthology, OMA, eggNOG and PANTHER. PANTHER and OMA had the most congruent ortholog groups with AYbRAH, while the other phylogenomic databases had greater amounts of under-clustering, over-clustering or no ortholog annotations for proteins. Future plans are discussed for AYbRAH, and recommendations are made for other research communities seeking to create curated ortholog databases.

Unique, Diverged, and Conserved Mitochondrial Functions Influencing Candida albicans Respiration

Candida albicans is an opportunistic fungal pathogen of major clinical concern. The virulence of this pathogen is intimately intertwined with its metabolism. Mitochondria, which have a central metabolic role, have undergone many lineage-specific adaptations in association with their eukaryotic host. A screen for lineage-specific genes identified seven such genes specific to the CTG clade of fungi, of which C. albicans is a member. Each is required for respiratory growth and is integral to expression of complex I, III, or IV of the electron transport chain. Two genes, NUO3 and NUO4, encode supernumerary subunits of complex I, whereas NUE1 and NUE2 have nonstructural roles in expression of complex I. Similarly, the other three genes have nonstructural roles in expression of complex III (QCE1) or complex IV (COE1 and COE2). In addition to these novel additions, an alternative functional assignment was found for the mitochondrial protein encoded by MNE1 MNE1 was required for complex I expression in C. albicans, whereas the distantly related Saccharomyces cerevisiae ortholog participates in expression of complex III. Phenotypic analysis of deletion mutants showed that fermentative metabolism is unable to support optimal growth rates or yields of C. albicans However, yeast-hypha morphogenesis, an important virulence attribute, did not require respiratory metabolism under hypoxic conditions. The inability to respire also resulted in hypersensitivity to the antifungal fluconazole and in attenuated virulence in a Galleria mellonella infection model. The results show that lineage-specific adaptations have occurred in C. albicans mitochondria and highlight the significance of respiratory metabolism in the pathobiology of C. albicans IMPORTANCE Candida albicans is an opportunistic fungal pathogen of major clinical concern. The virulence of this pathogen is intimately intertwined with its metabolic behavior, and mitochondria have a central role in that metabolism. Mitochondria have undergone many evolutionary changes, which include lineage-specific adaptations in association with their eukaryotic host. Seven lineage-specific genes required for electron transport chain function were identified in the CTG clade of fungi, of which C. albicans is a member. Additionally, examination of several highly diverged orthologs encoding mitochondrial proteins demonstrated functional reassignment for one of these. Deficits imparted by deletion of these genes revealed the critical role of respiration in virulence attributes of the fungus and highlight important evolutionary adaptations in C. albicans metabolism.

Snf2 controls pulcherriminic acid biosynthesis and antifungal activity of the biocontrol yeast Metschnikowia pulcherrima

Metschnikowia pulcherrima synthesises the pigment pulcherrimin, from cyclodileucine (cyclo(Leu-Leu)) as a precursor, and exhibits strong antifungal activity against notorious plant pathogenic fungi. This yeast therefore has great potential for biocontrol applications against fungal diseases; particularly in the phyllosphere where this species is frequently found. To elucidate the molecular basis of the antifungal activity of M. pulcherrima, we compared a wild-type strain with a spontaneously occurring, pigmentless, weakly antagonistic mutant derivative. Whole genome sequencing of the wild-type and mutant strains identified a point mutation that creates a premature stop codon in the transcriptional regulator gene SNF2 in the mutant. Complementation of the mutant strain with the wild-type SNF2 gene restored pigmentation and recovered the strong antifungal activity. Mass spectrometry (UPLC HR HESI-MS) proved the presence of the pulcherrimin precursors cyclo(Leu-Leu) and pulcherriminic acid and identified new precursor and degradation products of pulcherriminic acid and/or pulcherrimin. All of these compounds were identified in the wild-type and complemented strain, but were undetectable in the pigmentless snf2 mutant strain. These results thus identify Snf2 as a regulator of antifungal activity and pulcherriminic acid biosynthesis in M. pulcherrima and provide a starting point for deciphering the molecular functions underlying the antagonistic activity of this yeast.

Precise genome editing using a CRISPR-Cas9 method highlights the role of CoERG11 amino acid substitutions in azole resistance in Candida orthopsilosis

Background

Azoles are one of the main antifungal classes for the treatment of candidiasis. In the current context of emerging drug resistance, most studies have focused on Candida albicans, Candida glabrata or Candida auris but, so far, less is known about the underlying mechanisms of resistance in other species, including Candida orthopsilosis.

Objectives

We investigated azole resistance in a C. orthopsilosis clinical isolate recovered from a patient with haematological malignancy receiving fluconazole prophylaxis.

Methods

Antifungal susceptibility to fluconazole was determined in vitro (CLSI M27-A3) and in vivo (in a Galleria mellonella model of invasive candidiasis). The CoERG11 gene was then sequenced and amino acid substitutions identified were mapped on the predicted 3D structure of CoErg11p. A clustered regularly interspaced short palindromic repeat-Cas9 (CRISPR-Cas9) genome-editing strategy was used to introduce relevant mutations into a fluconazole-susceptible C. orthopsilosis isolate.

Results

Compared with unrelated C. orthopsilosis isolates, the clinical isolate exhibited both in vitro and in vivo fluconazole resistance. Sequencing of the CoERG11 gene identified several amino acid substitutions, including two possibly involved in fluconazole resistance (L376I and G458S). Both mutations mapped close to the active site of CoErg11p. Engineering these mutations in a different genetic background using CRISPR-Cas9 demonstrated that G458S, but not L376I, confers resistance to fluconazole and voriconazole.

Conclusions

Our data show that the G458S amino acid substitution in CoERG11p, but not L376I, contributes to azole resistance in C. orthopsilosis. In addition to highlighting the potential of CRISPR-Cas9 technology for precise genome editing in the field of antifungal resistance, we discuss some points that are critical to improving its efficiency.

Agglutinin-Like Sequence (ALS) Genes in the Candida parapsilosis Species Complex: Blurring the Boundaries Between Gene Families That Encode Cell-Wall Proteins

The agglutinin-like sequence (Als) proteins are best-characterized in Candida albicans and known for their role in adhesion of the fungal cell to host and abiotic surfaces. ALS sequences are often misassembled in whole-genome sequence data because each species has multiple ALS loci that contain similar sequences, most notably tandem copies of highly conserved repeated sequences. The Candida parapsilosis species complex includes Candida parapsilosis, Candida orthopsilosis, and Candida metapsilosis, three distinct but closely related species. Using publicly available genome resources, de novo genome assemblies, and laboratory experimentation including Sanger sequencing, five ALS genes were characterized in C. parapsilosis strain CDC317, three in C. orthopsilosis strain 90-125, and four in C. metapsilosis strain ATCC 96143. The newly characterized ALS genes shared similar features with the well-known C. albicans ALS family, but also displayed unique attributes such as novel short, imperfect repeat sequences that were found in other genes encoding fungal cell-wall proteins. Evidence of recombination between ALS sequences and other genes was most obvious in CmALS2265, which had the 5' end of an ALS gene and the repeated sequences and 3' end from the IFF/HYR family. Together, these results blur the boundaries between the fungal cell-wall families that were defined in C. albicans. TaqMan assays were used to quantify relative expression for each ALS gene. Some measurements were complicated by the assay location within the ALS gene. Considerable variation was noted in relative gene expression for isolates of the same species. Overall, however, there was a trend toward higher relative gene expression in saturated cultures rather than younger cultures. This work provides a complete description of the ALS genes in the C. parapsilosis species complex and a toolkit that promotes further investigations into the role of the Als proteins in host-fungal interactions.

Characterization of the Candida orthopsilosis agglutinin-like sequence (ALS) genes

Agglutinin like sequence (Als) cell-wall proteins play a key role in adhesion and virulence of Candida species. Compared to the well-characterized Candida albicans ALS genes, little is known about ALS genes in the Candida parapsilosis species complex. Three incomplete ALS genes were identified in the genome sequence for Candida orthopsilosis strain 90–125 (GenBank assembly ASM31587v1): CORT0C04210 (named CoALS4210), CORT0C04220 (CoALS4220) and CORT0B00800 (CoALS800). To complete the gene sequences, new data were derived from strain 90–125 using Illumina (short-read) and Oxford Nanopore (long-read) methods. Long-read sequencing analysis confirmed the presence of 3 ALS genes in C. orthopsilosis 90–125 and resolved the gaps located in repetitive regions of CoALS800 and CoALS4220. In the new genome assembly (GenBank PQBP00000000), the CoALS4210 sequence was slightly longer than in the original assembly. C. orthopsilosis Als proteins encoded features well-known in C. albicans Als proteins such as a secretory signal peptide, N-terminal domain with a peptide-binding cavity, amyloid-forming region, repeated sequences, and a C-terminal site for glycosylphosphatidylinositol anchor addition that, in yeast, suggest localization of the proteins in the cell wall. CoAls4210 and CoAls800 lacked the classic C. albicans Als tandem repeats, instead featuring short, imperfect repeats with consensus motifs such as SSSEPP and GSGN. Quantitative RT-PCR showed differential regulation of CoALS genes by growth stage in six genetically diverse C. orthopsilosis clinical isolates, which also exhibited length variation in the ALS alleles, and strain-specific gene expression patterns. Overall, long-read DNA sequencing methodology was instrumental in generating an accurate assembly of CoALS genes, thus revealing their unconventional features and first insights into their allelic variability within C. orthopsilosis clinical isolates.

Genetic Modification of Closely Related Candida Species

Species from the genus Candida are among the most important human fungal pathogens. Several of them are frequent commensals of the human microbiota but are also able to cause a variety of opportunistic infections, especially when the human host becomes immunocompromised. By far, most of the research to understand the molecular underpinnings of the pathogenesis of these species has focused on Candida albicans, the most virulent member of the genus. However, epidemiological data indicates that related Candida species are also clinically important. Here, we describe the generation of a set of strains and plasmids to genetically modify C. dubliniensis and C. tropicalis, the two pathogenic species most closely related to C. albicans. C. dubliniensis is an ideal model to understand C. albicans pathogenesis since it is the closest species to C. albicans but considerably less virulent. On the other hand, C. tropicalis is ranked among the four most common causes of infections by Candida species. Given that C. dubliniensis and C. tropicalis are obligate diploids with no known conventional sexual cycle, we generated strains that are auxotrophic for at least two amino acids which allows the tandem deletion of both alleles of a gene by complementing the two auxotrophies. The strains were generated in two different genetic backgrounds for each species — one for which the genomic sequence is available and a second clinically important one. In addition, we have adapted plasmids developed to delete genes and epitope/fluorophore tag proteins in C. albicans so that they can be employed in C. tropicalis. The tools generated here allow for efficient genetic modification of C. dubliniensis and C. tropicalis, and thus facilitate the study of the molecular basis of pathogenesis in these medically relevant fungi.

Candida parapsilosis: from Genes to the Bedside

Patients with suppressed immunity are at the highest risk for hospital-acquired infections. Among these, invasive candidiasis is the most prevalent systemic fungal nosocomial infection. Over recent decades, the combined prevalence of non-albicans Candida species outranked Candida albicans infections in several geographical regions worldwide, highlighting the need to understand their pathobiology in order to develop effective treatment and to prevent future outbreaks. Candida parapsilosis is the second or third most frequently isolated Candida species from patients. Besides being highly prevalent, its biology differs markedly from that of C. albicans, which may be associated with C. parapsilosis' increased incidence. Differences in virulence, regulatory and antifungal drug resistance mechanisms, and the patient groups at risk indicate that conclusions drawn from C. albicans pathobiology cannot be simply extrapolated to C. parapsilosis Such species-specific characteristics may also influence their recognition and elimination by the host and the efficacy of antifungal drugs. Due to the availability of high-throughput, state-of-the-art experimental tools and molecular genetic methods adapted to C. parapsilosis, genome and transcriptome studies are now available that greatly contribute to our understanding of what makes this species a threat. In this review, we summarize 10 years of findings on C. parapsilosis pathogenesis, including the species' genetic properties, transcriptome studies, host responses, and molecular mechanisms of virulence. Antifungal susceptibility studies and clinician perspectives are discussed. We also present regional incidence reports in order to provide an updated worldwide epidemiology summary.

Development and genomic elucidation of hybrid yeast with improved glucose-xylose co-fermentation at high temperature

Enhanced capability of co-fermenting glucose and xylose at high temperature is highly desirable for yeast application in second-generation bioethanol production. Here, we obtained hybrid strains with improved glucose-xylose co-fermentation properties at high temperature by combining genome shuffling and adaptive evolution. Genome resequencing of these strains suggested predominantly inherited genetic information from one parental strain Spathaspora passalidarum SP rather than the other parental strain Saccharomyces cerevisiae ScY01, possibly due to that the CUG codon system of S. passalidarum might have systematically eliminated most of the functional proteins from S. cerevisiae through misfolding. Compared to SP, one-copy loss of a 146-kb fragment was found in the hybrid strain and regained after being evolved for a while, whereas one-copy loss of an 11-kb fragment was only found after being evolved for a longer time. Besides, the genes affected by nonsynonymous variants were also identified, especially the mutation S540F in the endoplasmic reticulum chaperon Kar2. Structural prediction indicated that S540F might change the substrate binding activity of Kar2, and thus play a role in preventing protein aggregation in yeast at high temperature. Our results illustrated genomic alterations during this process and revealed some genomic factors that might be involved to determine yeast thermotolerance.