Genome evolution in yeasts

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

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

Summary

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

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Cost-effective in-house COVID-19 reverse transcription-polymerase chain reaction testing with yeast-derived Taq polymerase

BACKGROUND

Despite the decline of the COVID-19 pandemic, there continues to be a persistent requirement for reliable testing methods that can be adapted to future outbreaks and areas with limited resources. While the standard approach of using reverse transcription-polymerase chain reaction (RT-PCR) with Taq polymerase is effective, it faces challenges such as limited access to high-quality enzymes and the presence of bacterial DNA contamination in commercial kits, which can impact the accuracy of test results.

METHODS

This study investigates the production of recombinant Taq polymerase in yeast cells and assesses its crude lysate in a multiplex RT-PCR assay for detecting the SARS-CoV-2 RNA-dependent RNA polymerase (RdRP) and N genes, with human Ribonuclease P serving as an internal control.

RESULTS

The unpurified yeast Taq polymerase demonstrates sensitivity comparable to commercially purified bacterial Taq polymerase and unpurified bacterial counterparts in detecting the RdRP and N genes. It exhibits the highest specificity, with 100% accuracy, for the N gene. The specificity for the RdRP gene closely aligns with that of commercially purified bacterial Taq polymerase and unpurified bacterial Taq polymerase.

CONCLUSIONS

The use of unpurified recombinant yeast Taq polymerase shows promise as a cost-effective approach for conducting in-house COVID-19 RT-PCR testing. By eliminating the need for chromatography purification steps, the production of RT-PCR kits can be streamlined, potentially improving accessibility and scalability, especially in resource-limited settings and future pandemics.

Sister chromatid cohesion halts DNA loop expansion

Eukaryotic genomes are folded into DNA loops mediated by structural maintenance of chromosomes (SMC) complexes such as cohesin, condensin, and Smc5/6. This organization regulates different DNA-related processes along the cell cycle, such as transcription, recombination, segregation, and DNA repair. During the G2 stage, SMC-mediated DNA loops coexist with cohesin complexes involved in sister chromatid cohesion (SCC). However, the articulation between the establishment of SCC and the formation of SMC-mediated DNA loops along the chromatin remains unknown. Here, we show that SCC is indeed a barrier to cohesin-mediated DNA loop expansion along G2/M Saccharomyces cerevisiae chromosomes.

Genome-scale transcriptional activation by non-homologous end joining-mediated integration in Yarrowia lipolytica

BACKGROUND

Genome-scale screening can be applied to efficiently mine for unknown genes with phenotypes of interest or special functions. It is also useful to identify new targets for engineering desirable properties of cell factories.

RESULTS

Here, we designed a new approach for genome-scale transcription activation using non-homologous end joining (NHEJ)-mediated integration in Yarrowia lipolytica. We utilized this approach to screen for genes that, upon activation, confer phenotypes including improved acetic acid tolerance and xylose metabolism. The candidates were validated using gene overexpression, and functional changes including improved growth performance under multiple stressors and activated pentose metabolism were identified.

CONCLUSIONS

This study provides a simple and effective approach to randomly activate endogenous genes and mine for key targets associated with phenotypes of interest. The specific gene targets identified here will be useful for cell factory construction and biorefining lignocellulose.

Application of New Yarrowia lipolytica Transformants in Production of Citrates and Erythritol from Glycerol

Citric acid and erythritol are obtained on an industrial scale using biotechnological methods. Due to the growing market demand for these products, research is underway to improve the process economics by introducing new microorganisms, in particular of the species Yarrowia lipolytica. The aim of this study was to evaluate transformants of Y. lipolytica for growth and ability to overproduce citric acids and erythritol from glycerol. The transformants were constructed by overexpressing glycerol kinase, methylcitrate synthase and mitochondrial succinate-fumarate transporter in the mutant Wratislavia 1.31. Next, strains were assessed for biosynthesis of citrate (pH 5.5; nitrogen limitation) and erythritol (pH 3.0; high osmotic pressure) from glycerol. Regardless of culture conditions strains, 1.31.GUT1/6 and 1.31.GUT1/6.CIT1/3 exhibited high rates of substrate utilization. Under conditions favoring citrate biosynthesis, both strains produced several percent more citrates, accompanied by higher erythritol production compared to the parental strain. During erythritol biosynthesis, the strain 1.31.GUT1/6.CIT1/3.E34672g obtained as a result of co-expression of all three genes stood out, producing 84.0 g/L of erythritol with yield and productivity of 0.54 g/g and 0.72 g/Lh, respectively, which places it in the group of the highest-ranked producers of erythritol among Y. lipolytica species.

'Mother(Nature) knows best' - hijacking nature-designed transcriptional programs for enhancing stress resistance and protein production in Yarrowia lipolytica; presentation of YaliFunTome database

BACKGROUND

In the era of rationally designed synthetic biology, heterologous metabolites production, and other counter-nature engineering of cellular metabolism, we took a step back and recalled that 'Mother(-Nature) knows best'. While still aiming at synthetic, non-natural outcomes of generating an 'over-production phenotype' we dug into the pre-designed transcriptional programs evolved in our host organism-Yarrowia lipolytica, hoping that some of these fine-tuned orchestrated programs could be hijacked and used. Having an interest in the practical outcomes of the research, we targeted industrially-relevant functionalities-stress resistance and enhanced synthesis of proteins, and gauged them over extensive experimental design's completion.

RESULTS

Technically, the problem was addressed by screening a broad library of over 120 Y. lipolytica strains under 72 combinations of variables through a carefully pre-optimized high-throughput cultivation protocol, which enabled actual phenotype development. The abundance of the transcription program elicitors-transcription factors (TFs), was secured by their overexpression, while challenging the strains with the multitude of conditions was inflicted to impact their activation stratus. The data were subjected to mathematical modeling to increase their informativeness. The amount of the gathered data prompted us to present them in the form of a searchable catalog - the YaliFunTome database ( https://sparrow.up.poznan.pl/tsdatabase/ )-to facilitate the withdrawal of biological sense from numerical data. We succeeded in the identification of TFs that act as omni-boosters of protein synthesis, enhance resistance to limited oxygen availability, and improve protein synthesis capacity under inorganic nitrogen provision.

CONCLUSIONS

All potential users are invited to browse YaliFunTome in the search for homologous TFs and the TF-driven phenotypes of interest.

Dissecting an ancient stress resistance trait syndrome in the compost yeast Kluyveromyces marxianus

In the search to understand how evolution builds new traits, ancient events are often the hardest to dissect. Species-unique traits pose a particular challenge for geneticists-cases in which a character arose long ago and, in the modern day, is conserved within a species, distinguishing it from reproductively isolated relatives. In this work, we have developed the budding yeast genus Kluyveromyces as a model for mechanistic dissection of trait variation across species boundaries. Phenotypic profiling revealed robust heat and chemical-stress tolerance phenotypes that distinguished the compost yeast K. marxianus from the rest of the clade. We used culture-based, transcriptomic, and genetic approaches to characterize the metabolic requirements of the K. marxianus trait syndrome. We then generated a population-genomic resource for K. marxianus and harnessed it in molecular-evolution analyses, which found hundreds of housekeeping genes with evidence for adaptive protein variation unique to this species. Our data support a model in which, in the distant past, K. marxianus underwent a vastly complex remodeling of its proteome to achieve stress resistance. Such a polygenic architecture, involving nucleotide-level allelic variation on a massive scale, is consistent with theoretical models of the mechanisms of long-term adaptation, and suggests principles of broad relevance for interspecies trait genetics.

Metabolic cross-feeding enhances branched-chain aldehydes production in a synthetic community of fermented sausages

Microbial interactions play an important role in regulating the metabolic function of fermented food communities, especially the production of key flavor compounds. However, little is known about specific molecular mechanisms that regulate the production of key flavor compounds through microbial interactions. Here, we designed a synthetic consortium containing Debaryomyces hansenii D1, Staphylococcus xylosus S1, and Pediococcus pentosaceus PP1 to explore the mechanism of the microbial interactions underlying the branched-chain aldehydes production. In this consortium, firstly, D. hansenii secreted amino acids that promoted the growth of P. pentosaceus and S. xylosus. Specifically, D. hansenii D1 secreted alanine, aspartate, glutamate, glutamine, glycine, phenylalanine, serine, and threonine, which were the primary nutrients for bacterial growth. P. pentosaceus PP1 utilized all these eight amino acids through cross-feeding, whereas S. xylosus S1 did not utilize aspartate and serine. Furthermore, D. hansenii D1 promoted the production of branched-chain aldehydes from S. xylosus and P. pentosaceus through cross-feeding of α-keto acids (intermediate metabolites). Thus, the accumulation of 2-methyl-butanal was promoted in all co-culture. Overall, this work revealed the mechanism by which D. hansenii and bacteria cross-feed to produce branched-chain aldehydes in fermented sausages.

Zinc prevents vaginal candidiasis by inhibiting expression of an inflammatory fungal protein

Candida causes an estimated half-billion cases of vulvovaginal candidiasis (VVC) every year. VVC is most commonly caused by Candida albicans, which, in this setting, triggers nonprotective neutrophil infiltration, aggressive local inflammation, and symptomatic disease. Despite its prevalence, little is known about the molecular mechanisms underpinning the immunopathology of this fungal infection. In this study, we describe the molecular determinant of VVC immunopathology and a potentially straightforward way to prevent disease. In response to zinc limitation, C. albicans releases a trace mineral binding molecule called Pra1 (pH-regulated antigen). Here, we show that the PRA1 gene is strongly up-regulated during vaginal infections and that its expression positively correlated with proinflammatory cytokine concentrations in women. Genetic deletion of PRA1 prevented vaginal inflammation in mice, and application of a zinc solution down-regulated expression of the gene and also blocked immunopathology. We also show that treatment of women suffering from recurrent VVC with a zinc gel prevented reinfections. We have therefore identified a key mediator of symptomatic VVC, giving us an opportunity to develop a range of preventative measures for combatting this disease.

Taxogenomic placement of Rasporella oleae and Rasporella dianae gen. and spp. nov., two insect associated yeast species

During the course of independent studies in Europe, North America, and Africa, seven yeast strains were isolated from insect frass, decaying wood, tree flux, and olive oil sediment. Phylogenetic analysis of two barcoding DNA regions (internal transcribed spacer and the D1/D2 domain of the LSU rRNA gene) revealed that they belong to two closely related undescribed species distinct from all genera in the family Debaryomycetaceae. For reliable taxonomic placement the genomes of four strains of the two novel species and six type strains of closely related species were sequenced. Orthologous genes from 54 genomes of representatives of the Pichiomycetes and 23 outgroup taxa were concatenated to construct a fully supported phylogenetic tree. Consistent with the assumptions, we found that the two new species belong to a novel genus. In addition, the delimitation of the novel species was supported by genetic distance calculations from average nucleotide identity (ANI) and digital DNA:DNA hybridization (dDDH) values. The physiological characterization of the novel species was generally consistent with their genomic content. All strains had two alleles encoding secretory lipase in either two or three copies depending on the species. However, lipolytic activity was detected only in strains with three copies of the secretory lipase gene. Nevertheless, lipolytic activity might be related to their association with the insect gut. Based on these results, formal descriptions of the new genus Rasporella gen. nov. and of two new species Rasporella dianae sp. nov. (holotype UCDFST 68-643T , MycoBank no.: 850238) and Rasporella oleae sp. nov. (holotype ZIM 2471T , MycoBank no.: 850126) are provided.

Multi-omics view of recombinant Yarrowia lipolytica: Enhanced ketogenic amino acid catabolism increases polyketide-synthase-driven docosahexaenoic production to high selectivity at the gram scale

DHA is a marine PUFA of commercial value, given its multiple health benefits. The worldwide emerging shortage in DHA supply has increased interest in microbial cell factories that can provide the compound de novo. In this regard, the present work aimed to improve DHA production in the oleaginous yeast strain Y. lipolytica Af4, which synthetized the PUFA via a heterologous myxobacterial polyketide synthase (PKS)-like gene cluster. As starting point, we used transcriptomics, metabolomics, and 13C-based metabolic pathway profiling to study the cellular dynamics of Y. lipolytica Af4. The shift from the growth to the stationary DHA-production phase was associated with fundamental changes in carbon core metabolism, including a strong upregulation of the PUFA gene cluster, as well as an increase in citrate and fatty acid degradation. At the same time, the intracellular levels of the two DHA precursors acetyl-CoA and malonyl-CoA dropped by up to 98% into the picomolar range. Interestingly, the degradation pathways for the ketogenic amino acids l-lysine, l-leucine, and l-isoleucine were transcriptionally activated, presumably to provide extra acetyl-CoA. Supplementation with small amounts of these amino acids at the beginning of the DHA production phase beneficially increased the intracellular CoA-ester pools and boosted the DHA titer by almost 40%. Isotopic 13C-tracer studies revealed that the supplements were efficiently directed toward intracellular CoA-esters and DHA. Hereby, l-lysine was found to be most efficient, as it enabled long-term activation, due to storage within the vacuole and continuous breakdown. The novel strategy enabled DHA production in Y. lipolytica at the gram scale for the first time. DHA was produced at a high selectivity (27% of total fatty acids) and free of the structurally similar PUFA DPA, which facilitates purification for high-value medical applications that require API-grade DHA. The assembled multi-omics picture of the central metabolism of Y. lipolytica provides valuable insights into this important yeast. Beyond our work, the enhanced catabolism of ketogenic amino acids seems promising for the overproduction of other compounds in Y. lipolytica, whose synthesis is limited by the availability of CoA ester precursors.

Extending the reach of homology by using successive computational filters to find yeast pheromone genes

The mating of fungi depends on pheromones that mediate communication between two mating types. Most species use short peptides as pheromones, which are either unmodified (e.g., α-factor in Saccharomyces cerevisiae) or C-terminally farnesylated (e.g., a-factor in S. cerevisiae). Peptide pheromones have been found by genetics or biochemistry in a small number of fungi, but their short sequences and modest conservation make it impossible to detect homologous sequences in most species. To overcome this problem, we used a four-step computational pipeline to identify candidate a-factor genes in sequenced genomes of the Saccharomycotina, the fungal clade that contains most of the yeasts: we require that candidate genes have a C-terminal prenylation motif, are shorter than 100 amino acids long, and contain a proteolytic-processing motif upstream of the potential mature pheromone sequence and that closely related species contain highly conserved homologs of the potential mature pheromone sequence. Additional manual curation exploits the observation that many species carry more than one a-factor gene, encoding identical or nearly identical pheromones. From 332 Saccharomycotina genomes, we identified strong candidate pheromone genes in 241 genomes, covering 13 clades that are each separated from each other by at least 100 million years, the time required for evolution to remove detectable sequence homology among small pheromone genes. For one small clade, the Yarrowia, we demonstrated that our algorithm found the a-factor genes: deleting all four related genes in the a-mating type of Yarrowia lipolytica prevents mating.

Deubiquitination module is critical for oxidative stress response and biofilm formation in Candida glabrata

Candidiasis is one of the most important fungal diseases and generally refers to diseases of the skin or mucosal tissues caused by Candida species. Candida glabrata is an opportunistic human fungal pathogen. Infection with C. glabrata has significantly increased due to innate antifungal drug tolerance and the ability to adhere to mucocutaneous surfaces. Spt-Ada-Gcn5 acetyltransferase complex contains two different post-translational modifications, histone acetylation (HAT) module and deubiquitination (DUB) module, which are decisive in gene regulation and highly conserved in eukaryotes. Previous research in our laboratory found that the HAT module ADA2 could regulate C. glabrata oxidative stress tolerance, drug tolerance, cell wall integrity, and virulence. However, the roles of the DUB module that is comprised of UBP8, SGF11, SGF73, and SUS1 genes in those phenotypes are not yet understood. In this study, we found that DUB module genes UBP8, SGF11, and SUS1, but not SGF73 positively regulate histone H2B DUB. Furthermore, ubp8, sgf11, and sus1 mutants exhibited decreased biofilm formation and sensitivity to cell wall-perturbing agent sodium dodecyl sulfate and antifungal drug amphotericin B. In addition, the sgf73 mutant showed increased biofilm formation but was susceptible to oxidative stresses, antifungal drugs, and cell wall perturbing agents. The ubp8, sgf11, and sus1 mutants showed marginal hypovirulence, whereas the sgf73 mutant exhibited virulence similar to the wild type in a murine systemic infection model. In conclusion, the C. glabrata DUB module plays distinct roles in H2B ubiquitination, oxidative stress response, biofilm formation, cell wall integrity, and drug tolerance, but exhibits minor roles in virulence.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Divergence of TORC1-mediated stress response leads to novel acquired stress resistance in a pathogenic yeast

Acquired stress resistance (ASR) enables organisms to prepare for environmental changes that occur after an initial stressor. However, the genetic basis for ASR and how the underlying network evolved remain poorly understood. In this study, we discovered that a short phosphate starvation induces oxidative stress response (OSR) genes in the pathogenic yeast C. glabrata and protects it against a severe H2O2 stress; the same treatment, however, provides little benefit in the low pathogenic-potential relative, S. cerevisiae. This ASR involves the same transcription factors (TFs) as the OSR, but with different combinatorial logics. We show that Target-of-Rapamycin Complex 1 (TORC1) is differentially inhibited by phosphate starvation in the two species and contributes to the ASR via its proximal effector, Sch9. Therefore, evolution of the phosphate starvation-induced ASR involves the rewiring of TORC1's response to phosphate limitation and the repurposing of TF-target gene networks for the OSR using new regulatory logics.

Genotypic diversity and unrecognized antifungal resistance among populations of Candida glabrata from positive blood cultures

The longstanding model is that most bloodstream infections (BSIs) are caused by a single organism. We perform whole genome sequencing of five-to-ten strains from blood culture (BC) bottles in each of ten patients with Candida glabrata BSI. We demonstrate that BCs contain mixed populations of clonal but genetically diverse strains. Genetically distinct strains from two patients exhibit phenotypes that are potentially important during BSIs, including differences in susceptibility to antifungal agents and phagocytosis. In both patients, the clinical microbiology lab recovered a fluconazole-susceptible index strain, but we identify mixed fluconazole-susceptible and -resistant populations. Diversity in drug susceptibility is likely clinically relevant, as fluconazole-resistant strains were subsequently recovered by the clinical laboratory during persistent or relapsing infections. In one patient, unrecognized respiration-deficient small colony variants are fluconazole-resistant and significantly attenuated for virulence during murine candidiasis. Our data suggest a population-based model of C. glabrata genotypic and phenotypic diversity during BSIs.

Divergence of TORC1-mediated Stress Response Leads to Novel Acquired Stress Resistance in a Pathogenic Yeast

Acquired stress resistance (ASR) enables organisms to prepare for environmental changes that occur after an initial stressor. However, the genetic basis for ASR and how the underlying network evolved remain poorly understood. In this study, we discovered that a short phosphate starvation induces oxidative stress response (OSR) genes in the pathogenic yeast C. glabrata and protects it against a severe H2O2 stress; the same treatment, however, provides little benefit in the low pathogenic-potential relative, S. cerevisiae. This ASR involves the same transcription factors (TFs) as the OSR, but with different combinatorial logics. We show that Target-of-Rapamycin Complex 1 (TORC1) is differentially inhibited by phosphate starvation in the two species and contributes to the ASR via its proximal effector, Sch9. Therefore, evolution of the phosphate starvation-induced ASR involves the rewiring of TORC1's response to phosphate limitation and the repurposing of TF-target gene networks for the OSR using new regulatory logics.

Fungal Alcohol Dehydrogenases: Physiological Function, Molecular Properties, Regulation of Their Production, and Biotechnological Potential

Fungal alcohol dehydrogenases (ADHs) participate in growth under aerobic or anaerobic conditions, morphogenetic processes, and pathogenesis of diverse fungal genera. These processes are associated with metabolic operation routes related to alcohol, aldehyde, and acid production. The number of ADH enzymes, their metabolic roles, and their functions vary within fungal species. The most studied ADHs are associated with ethanol metabolism, either as fermentative enzymes involved in the production of this alcohol or as oxidative enzymes necessary for the use of ethanol as a carbon source; other enzymes participate in survival under microaerobic conditions. The fast generation of data using genome sequencing provides an excellent opportunity to determine a correlation between the number of ADHs and fungal lifestyle. Therefore, this review aims to summarize the latest knowledge about the importance of ADH enzymes in the physiology and metabolism of fungal cells, as well as their structure, regulation, evolutionary relationships, and biotechnological potential.

Pyruvate decarboxylase and thiamine biosynthetic genes are regulated differently by Pdc2 in S. cerevisiae and C. glabrata

Understanding metabolism in the pathogen Candida glabrata is key to identifying new targets for antifungals. The thiamine biosynthetic (THI) pathway is partially defective in C. glabrata, but the transcription factor CgPdc2 upregulates some thiamine biosynthetic and transport genes. One of these genes encodes a recently evolved thiamine pyrophosphatase (CgPMU3) that is critical for accessing external thiamine. Here, we demonstrate that CgPdc2 primarily regulates THI genes. In Saccharomyces cerevisiae, Pdc2 regulates both THI and pyruvate decarboxylase (PDC) genes, with PDC proteins being a major thiamine sink. Deletion of PDC2 is lethal in S. cerevisiae in standard growth conditions, but not in C. glabrata. We uncover cryptic cis elements in C. glabrata PDC promoters that still allow for regulation by ScPdc2, even when that regulation is not apparent in C. glabrata. C. glabrata lacks Thi2, and it is likely that inclusion of Thi2 into transcriptional regulation in S. cerevisiae allows for a more complex regulation pattern and regulation of THI and PDC genes. We present evidence that Pdc2 functions independent of Thi2 and Thi3 in both species. The C-terminal activation domain of Pdc2 is intrinsically disordered and critical for species differences. Truncation of the disordered domains leads to a gradual loss of activity. Through a series of cross species complementation assays of transcription, we suggest that there are multiple Pdc2-containing complexes, and C. glabrata appears to have the simplest requirement set for THI genes, except for CgPMU3. CgPMU3 has different cis requirements, but still requires Pdc2 and Thi3 to be upregulated by thiamine starvation. We identify the minimal region sufficient for thiamine regulation in CgTHI20, CgPMU3, and ScPDC5 promoters. Defining the cis and trans requirements for THI promoters should lead to an understanding of how to interrupt their upregulation and provide targets in metabolism for antifungals.

Transcriptome analysis reveals a common adaptive transcriptional response of Candida glabrata to diverse environmental stresses

Candida glabrata, an opportunistic fungal pathogen, causes superficial and life-threatening infections in humans. In the host microenvironment, C. glabrata encounters a variety of stresses, and its ability to cope with these stresses is crucial for its pathogenesis. To gain insights into how C. glabrata adapts to adverse environmental conditions, we examined its transcriptional landscape under heat, osmotic, cell wall, oxidative, and genotoxic stresses using RNA sequencing and reveal that C. glabrata displays a diverse transcriptional response involving ∼75% of its genome for adaptation to different environmental stresses. C. glabrata mounts a central common adaptation response wherein ∼25% of all genes (n = 1370) are regulated in a similar fashion at different environmental stresses. Elevated cellular translation and diminished mitochondrial activity-associated transcriptional signature characterize the common adaptation response. Transcriptional regulatory association networks of common adaptation response genes revealed a set of 29 transcription factors acting as potential activators and repressors of associated adaptive response genes. Overall, the current work delineates the adaptive responses of C. glabrata to diverse environmental stresses and reports the existence of a common adaptive transcriptional response upon prolonged exposure to environmental stresses.

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.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.

Engineering Strategies for Efficient Bioconversion of Glycerol to Value-Added Products by Yarrowia lipolytica

Yarrowia lipolytica has been a valuable biotechnological workhorse for the production of commercially important biochemicals for over 70 years. The knowledge gained so far on the native biosynthetic pathways, as well as the availability of numerous systems and synthetic biology tools, enabled not only the regulation and the redesign of the existing metabolic pathways, but also the introduction of novel synthetic ones; further consolidating the position of the yeast in industrial biotechnology. However, for the development of competitive and sustainable biotechnological production processes, bioengineering should be reinforced by bioprocess optimization strategies. Although there are many published reviews on the bioconversion of various carbon sources to value-added products by Yarrowia lipolytica, fewer works have focused on reviewing up-to-date strain, medium, and process engineering strategies with an aim to emphasize the significance of integrated engineering approaches. The ultimate goal of this work is to summarize the necessary knowledge and inspire novel routes to manipulate at a systems level the yeast biosynthetic machineries by combining strain and bioprocess engineering. Due to the increasing surplus of biodiesel-derived waste glycerol and the favored glycerol-utilization metabolic pathways of Y. lipolytica over other carbon sources, the present review focuses on pure and crude glycerol-based biomanufacturing.

Metabolic Engineering of Yarrowia lipolytica for Terpenoid Production: Tools and Strategies

Terpenoids are a diverse group of compounds with isoprene units as basic building blocks. They are widely used in the food, feed, pharmaceutical, and cosmetic industries due to their diverse biological functions such as antioxidant, anticancer, and immune enhancement. With an increase in understanding the biosynthetic pathways of terpenoids and advances in synthetic biology techniques, microbial cell factories have been built for the heterologous production of terpenoids, with the oleaginous yeast Yarrowia lipolytica emerging as an outstanding chassis. In this paper, recent progress in the development of Y. lipolytica cell factories for terpenoid production with a focus on the advances in novel synbio tools and metabolic engineering strategies toward enhanced terpenoid biosynthesis is reviewed.

Draft Genome Sequence of the Fluconazole-Resistant Yarrowia lipolytica Clinical Isolate CBS 18115

Yarrowia lipolytica is a nonconventional yeast of industrial interest that can sometimes act as an opportunistic pathogen and is responsible for invasive fungal infections. We report the draft genome sequence of the fluconazole-resistant strain CBS 18115, which was isolated from a blood culture. The Y132F substitution in ERG11, which was previously described in fluconazole-resistant Candida isolates, was identified.

High-throughput functional profiling of the human fungal pathogen Candida albicans genome

Candida albicans is a major fungal pathogen of humans. Although its genome has been sequenced more than two decades ago, there are still over 4300 uncharacterized C. albicans genes. We previously generated an ORFeome as well as a collection of destination vectors to facilitate overexpression of C. albicans ORFs. Here, we report the construction of ∼2500 overexpression mutants and their evaluation by in vitro spotting on rich medium and in a liquid pool experiment in rich medium, allowing the identification of genes whose overexpression has a fitness cost. The candidates were further validated at the individual strain level. This new resource allows large-scale screens in different growth conditions to be performed routinely. Altogether, based on the concept of identifying functionally related genes by cluster analysis, the availability of this overexpression mutant collection will facilitate the characterization of gene functions in C. albicans.

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

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

The Gcn5-Ada2-Ada3 histone acetyltransferase module has divergent roles in pathogenesis of Candida glabrata

Candida glabrata is an opportunistic fungal pathogen and the second most prevalent species isolated from candidiasis patients. C. glabrata has intrinsic tolerance to antifungal drugs and oxidative stresses and the ability to adhere to mucocutaneous surfaces. However, knowledge about the regulation of its virulence traits is limited. The Spt-Ada-Gcn5 acetyltransferase (SAGA) complex modulates gene transcription by histone acetylation through the histone acetyltransferase (HAT) module comprised of Gcn5-Ada2-Ada3. Previously, we showed that the ada2 mutant was hypervirulent but displayed decreased tolerance to antifungal drugs and cell wall perturbing agents. In this study, we further characterized the functions of Ada3 and Gcn5 in C. glabrata. We found that single, double, or triple deletions of the HAT module, as expected, resulted in a decreased level of acetylation on histone H3 lysine 9 (H3K9) and defective growth. These mutants were more susceptible to antifungal drugs, oxidative stresses, and cell wall perturbing agents compared with the wild-type. In addition, HAT module mutants exhibited enhanced agar invasion and upregulation of adhesin and proteases encoding genes, whereas the biofilm formation of those mutants was impaired. Interestingly, HAT module mutants exhibited enhanced induction of catalases (CTA1) expression upon treatment with H2O2 compared with the wild-type. Lastly, although ada3 and gcn5 exhibited marginal hypervirulence, the HAT double and triple mutants were hypervirulent in a murine model of candidiasis. In conclusion, the HAT module of the SAGA complex plays unique roles in H3K9 acetylation, drug tolerance, oxidative stress response, adherence, and virulence in C. glabrata.

Heterologous expression reveals unique properties of Trk K+ importers from nonconventional biotechnologically relevant yeast species together with their potential to support Saccharomyces cerevisiae growth

In the model yeast Saccharomyces cerevisiae, Trk1 is the main K⁺ importer. It is involved in many important physiological processes, such as the maintenance of ion homeostasis, cell volume, intracellular pH, and plasma-membrane potential. The ScTrk1 protein can be of great interest to industry, as it was shown that changes in its activity influence ethanol production and tolerance in S. cerevisiae and also cell performance in the presence of organic acids or high ammonium under low K⁺ conditions. Nonconventional yeast species are attracting attention due to their unique properties and as a potential source of genes that encode proteins with unusual characteristics. In this work, we aimed to study and compare Trk proteins from Debaryomyces hansenii, Hortaea werneckii, Kluyveromyces marxianus, and Yarrowia lipolytica, four biotechnologically relevant yeasts that tolerate various extreme environments. Heterologous expression in S. cerevisiae cells lacking the endogenous Trk importers revealed differences in the studied Trk proteins' abilities to support the growth of cells under various cultivation conditions such as low K⁺ or the presence of toxic cations, to reduce plasma-membrane potential or to take up Rb⁺ . Examination of the potential of Trks to support the stress resistance of S. cerevisiae wild-type strains showed that Y. lipolytica Trk1 is a promising tool for improving cell tolerance to both low K⁺ and high salt and that the overproduction of S. cerevisiae's own Trk1 was the most efficient at improving the growth of cells in the presence of highly toxic Li⁺ ions.

Contrasting Genomic Evolution Between Domesticated and Wild Kluyveromyces lactis Yeast Populations

The process of domestication has variable consequences on genome evolution leading to different phenotypic signatures. Access to the complete genome sequences of a large number of individuals makes it possible to explore the different facets of this domestication process. Here, we sought to explore the genome evolution of Kluyveromyces lactis, a yeast species well known for its involvement in dairy processes and also present in natural environments. Using a combination of short- and long-read sequencing strategies, we investigated the genomic variability of 41 K. lactis isolates and found that the overall genetic diversity of this species is very high (θw = 3.3 × 10-2) compared with other species such as Saccharomyces cerevisiae (θw = 1.6 × 10-2). However, the domesticated dairy population shows a reduced level of diversity (θw = 1 × 10-3), probably due to a domestication bottleneck. In addition, this entire population is characterized by the introgression of the LAC4 and LAC12 genes, responsible for lactose fermentation and coming from the closely related species, Kluyveromyces marxianus, as previously described. Our results highlighted that the LAC4/LAC12 gene cluster was acquired through multiple and independent introgression events. Finally, we also identified several genes that could play a role in adaptation to dairy environments through copy number variation. These genes are involved in sugar consumption, flocculation, and drug resistance, and may play a role in dairy processes. Overall, our study illustrates contrasting genomic evolution and sheds new light on the impact of domestication processes on it.

Bidirectional hybrid erythritol-inducible promoter for synthetic biology in Yarrowia lipolytica

BACKGROUND

The oleaginous yeast Yarrowia lipolytica is increasingly used as a chassis strain for generating bioproducts. Several hybrid promoters with different strengths have been developed by combining multiple copies of an upstream activating sequence (UAS) associated with a TATA box and a core promoter. These promoters display either constitutive, phase-dependent, or inducible strong expression. However, there remains a lack of bidirectional inducible promoters for co-expressing genes in Y. lipolytica.

RESULTS

This study built on our previous work isolating and characterizing the UAS of the erythritol-induced genes EYK1 and EYD1 (UAS-eyk1). We found an erythritol-inducible bidirectional promoter (BDP) located in the EYK1-EYL1 intergenic region. We used the BDP to co-produce YFP and RedStarII fluorescent proteins and demonstrated that the promoter's strength was 2.7 to 3.5-fold stronger in the EYL1 orientation compared to the EYK1 orientation. We developed a hybrid erythritol-inducible bidirectional promoter (HBDP) containing five copies of UAS-eyk1 in both orientations. It led to expression levels 8.6 to 19.2-fold higher than the native bidirectional promoter. While the BDP had a twofold-lower expression level than the strong constitutive TEF promoter, the HBDP had a 5.0-fold higher expression level when oriented toward EYL1 and a 2.4-fold higher expression level when oriented toward EYK1. We identified the optimal media for BDP usage by exploring yeast growth under microbioreactor conditions. Additionally, we constructed novel Golden Gate biobricks and a destination vector for general use.

CONCLUSIONS

In this research, we developed novel bidirectional and hybrid bidirectional promoters of which expression can be fine-tuned, responding to the need for versatile promoters in the yeast Y. lipolytica. This study provides effective tools that can be employed to smoothly adjust the erythritol-inducible co-expression of two target genes in biotechnology applications. BDPs developed in this study have potential applications in the fields of heterologous protein production, metabolic engineering, and synthetic biology.

Insights into the Genomic and Phenotypic Landscape of the Oleaginous Yeast Yarrowia lipolytica

Although Yarrowia lipolytica is a model yeast for the study of lipid metabolism, its diversity is poorly known, as studies generally consider only a few standard laboratory strains. To extend our knowledge of this biotechnological workhorse, we investigated the genomic and phenotypic diversity of 56 natural isolates. Y. lipolytica is classified into five clades with no correlation between clade membership and geographic or ecological origin. A low genetic diversity (π = 0.0017) and a pan-genome (6528 genes) barely different from the core genome (6315 genes) suggest Y. lipolytica is a recently evolving species. Large segmental duplications were detected, totaling 892 genes. With three new LTR-retrotransposons of the Gypsy family (Tyl4, Tyl9, and Tyl10), the transposable element content of genomes appeared diversified but still low (from 0.36% to 3.62%). We quantified 34 traits with substantial phenotypic diversity, but genome-wide association studies failed to evidence any associations. Instead, we investigated known genes and found four mutational events leading to XPR2 protease inactivation. Regarding lipid metabolism, most high-impact mutations were found in family-belonging genes, such as ALK or LIP, and therefore had a low phenotypic impact, suggesting that the huge diversity of lipid synthesis and accumulation is multifactorial or due to complex regulations.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Engineering of Yarrowia lipolytica for terpenoid production

Terpenoids are a group of chemicals of great importance for human health and prosperity. Terpenoids can be used for human and animal nutrition, treating diseases, enhancing agricultural output, biofuels, fragrances, cosmetics, and flavouring. However, due to the rapid depletion of global natural resources and manufacturing practices relying on unsustainable petrochemical synthesis, there is a need for economic alternatives to supply the world's demand for these essential chemicals. Microbial biosynthesis offers the means to develop scalable and sustainable bioprocesses for terpenoid production. In particular, the non-conventional yeast Yarrowia lipolytica demonstrates excellent potential as a chassis for terpenoid production due to its amenability to industrial production scale-up, genetic engineering, and high accumulation of terpenoid precursors. This review aims to illustrate the scientific progress in developing Y. lipolytica terpenoid cell factories, focusing on metabolic engineering approaches for strain improvement and cultivation optimization.

Engineering Yarrowia lipolytica to produce nutritional fatty acids: Current status and future perspectives

Due to their vital physiological functions, nutritional fatty acids have great potential as nutraceutical food supplements for preventing an array of diseases such as inflammation, depression, arthritis, osteoporosis, diabetes and cancer. Microbial biosynthesis of fatty acids follows the trend of sustainable development, as it enables green, environmentally friendly and efficient production. As a natural oleaginous yeast, Yarrowia lipolytica is especially well-suited for the production of fatty acids. Moreover, it has a variety of genetic engineering tools and novel metabolic engineering strategies that make it a robust workhorse for the production of an array of value-added products. In this review, we summarize recent advances in metabolic engineering strategies for accumulating nutritional fatty acids in Y. lipolytica, including conjugated fatty acids and polyunsaturated fatty acids. In addition, the future prospects of nutritional fatty acid production using the Y. lipolytica platform are discussed in light of the current progress, challenges, and trends in this field. Finally, guidelines for future studies are also emphasized.

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

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

Rapid Gene Target Tracking for Enhancing β-Carotene Production Using Flow Cytometry-Based High-Throughput Screening in Yarrowia lipolytica

β-Carotene is a provitamin A precursor and an important antioxidant that is used widely in the aquaculture, food, cosmetic, and pharmaceutical industries. Oleaginous Yarrowia lipolytica has been demonstrated as a competitive producer microorganism for the production of hydrophobic β-carotene through rational engineering strategies. However, the limited understanding of the complexity of the metabolic network between carotenoid biosynthesis and other cellular processes has hampered further advancement. Genome-scale mutagenesis and high-throughput screening of mutagenesis libraries have been extensively employed in gene mining or in the identification of key targets associated with particular phenotypes. In this study, we developed a fluorescence-activated cell-sorting approach for the effective high-throughput screening of high-β-carotene-producing strains. Using this approach, millions of mutants were screened rapidly, and new gene targets involved in lipid metabolism, sterol metabolism, signal transduction, and stress response were identified. The disruption of the genes affecting fatty acid oxidation, lipid composition, and sterol transcriptional regulation (4CL-8, GCS, and YIsterTF) increased β-carotene significantly. By engineering these targets in a high-β-carotene production, a strain that produced 9.4 g/L β-carotene was constructed. Here, we used a flow cytometry approach to improve screening efficiency and eliminate the interference of intermediate metabolites. The targets obtained in this study can be used in studies focusing on metabolic engineering in the future for improving carotenoid production. IMPORTANCE β-Carotene is a high-value-added product that is widely used in the aquaculture, food, cosmetic, and pharmaceutical industries. In our previous study, Yarrowia lipolytica has been engineered extensively to produce β-carotene. To further improve its production, high-throughput screening and the identification of new beneficial gene targets are required. Herein, we developed a fluorescence-activated cell-sorting approach for the effective high-throughput screening of high-β-carotene-producing strains. Using this approach, millions of mutants were screened rapidly, and new gene targets involved in lipid metabolism, sterol metabolism, signal transduction, and stress response were identified. The disruption of the genes affecting fatty acid oxidation, lipid composition, and sterol transcriptional regulation (4CL-8, GCS, and YIsterTF) increased β-carotene significantly. By engineering these targets in a high-β-carotene production, a strain that produced 9.4 g/L β-carotene was constructed.

Distinct chromosomal “niches” in the genome of Saccharomyces cerevisiae provide the background for genomic innovation and shape the fate of gene duplicates

Nearly one third of Saccharomyces cerevisiae protein coding sequences correspond to duplicate genes, equally split between small-scale duplicates (SSD) and whole-genome duplicates (WGD). While duplicate genes have distinct properties compared to singletons, to date, there has been no systematic analysis of their positional preferences. In this work, we show that SSD and WGD genes are organized in distinct gene clusters that occupy different genomic regions, with SSD being more peripheral and WGD more centrally positioned close to centromeric chromatin. Duplicate gene clusters differ from the rest of the genome in terms of gene size and spacing, gene expression variability and regulatory complexity, properties that are also shared by singleton genes residing within them. Singletons within duplicate gene clusters have longer promoters, more complex structure and a higher number of protein–protein interactions. Particular chromatin architectures appear to be important for gene evolution, as we find SSD gene-pair co-expression to be strongly associated with the similarity of nucleosome positioning patterns. We propose that specific regions of the yeast genome provide a favourable environment for the generation and maintenance of small-scale gene duplicates, segregating them from WGD-enriched genomic domains. Our findings provide a valuable framework linking genomic innovation with positional genomic preferences.

Papiliotrema laurentii: general features and biotechnological applications

Papiliotrema laurentii, previously classified as Cryptococcus laurentii, is an oleaginous yeast that has been isolated from soil, plants, and agricultural and industrial residues. This variety of habitats reflects the diversity of carbon sources that it can metabolize, including monosaccharides, oligosaccharides, glycerol, organic acids, and oils. Compared to other oleaginous yeasts, such as Yarrowia lipolytica and Rhodotorula toruloides, there is little information regarding its genetic and physiological characteristics. From a biotechnological point of view, P. laurentii can produce surfactants, enzymes, and high concentrations of lipids, which can be used as feedstock for fatty acid-derived products. Moreover, it can be applied for the biocontrol of phytopathogenic fungi, contributing to quality maintenance in post- and pre-harvest fruits. It can also improve mycorrhizal colonization, nitrogen nutrition, and plant growth. P. laurentii is also capable of degrading polyester and diesel derivatives and acting in the bioremediation of heavy metals. In this review, we present the current knowledge about the basic and applied aspects of P. laurentii, underscoring its biotechnological potential and future perspectives. KEY POINTS: • The physiological characteristics of P. laurentii confer a wide range of biotechnological applications. • The regulation of the acetyl-CoA carboxylase in P. laurentii is different from most other oleaginous yeasts. • The GEM is a valuable tool to guide the construction of engineered P. laurentii strains with improved features for bio-based products.

Characterization of the Secretome of Pathogenic Candida glabrata and Their Effectiveness against Systemic Candidiasis in BALB/c Mice for Vaccine Development

Infections by non-albicans Candida species have increased drastically in the past few decades. Candida glabrata is one of the most common opportunistic fungal pathogens in immunocompromised individuals, owing to its capability to attach to various human cell types and medical devices and being intrinsically weakly susceptible to azoles. Immunotherapy, including the development of antifungal vaccines, has been recognized as an alternative approach for preventing and treating fungal infections. Secretory proteins play a crucial role in establishing host–pathogen interactions and are also responsible for eliciting an immune response in the host during candidiasis. Therefore, fungal secretomes can provide promising protein candidates for antifungal vaccine development. This study attempts to uncover the presence of immunodominant antigenic proteins in the C. glabrata secretome and delineate their role in various biological processes and their potency in the development of antifungal vaccines. LC–MS/MS results uncovered that C. glabrata secretome consisted of 583 proteins, among which 33 were identified as antigenic proteins. The protection ability of secretory proteins against hematogenously disseminated infection caused by C. glabrata was evaluated in BALB/c mice. After immunization and booster doses, all the animals were challenged with a lethal dose of C. glabrata. All the mice showing signs of distress were sacrificed post-infection, and target organs were collected, followed by histopathology and C. glabrata (CFU/mg) estimation. Our results showed a lower fungal burden in target organs and increased survival in immunized mice compared to the infection control group, thus revealing the immunogenic property of secreted proteins. Thus, identified secretome proteins of C. glabrata have the potential to act as antigenic proteins, which can serve as potential candidates for the development of antifungal vaccines. This study also emphasizes the importance of a mass-spectrometry approach to identifying the antigenic proteins in C. glabrata secretome.

Current Ethanol Production Requirements for the Yeast Saccharomyces cerevisiae

An increase in global energy demand has caused oil prices to reach record levels in recent times. High oil prices together with concerns over CO₂ emissions have resulted in renewed interest in renewable energy. Nowadays, ethanol is the principal renewable biofuel. However, the industrial need for increased productivity, wider substrate range utilization, and the production of novel compounds leads to renewed interest in further extending the use of current industrial strains by exploiting the immense, and still unknown, potential of natural yeast strains. This review seeks to answer the following questions: (a) which characteristics should S. cerevisiae have for the current production of first- and second-generation ethanol? (b) Why are alcohol-tolerance and thermo-tolerance characteristics required? (c) Which genes are related to these characteristics? (d) What are the advances that can be achieved with the isolation of new organisms from the environment?

Genomic features, aroma profiles, and probiotic potential of the Debaryomyces hansenii species complex strains isolated from Korean soybean fermented food

Fermented soybean products are gaining attention in the food industry owing to their nutritive value and health benefits. In this study, we performed genomic analysis and physiological characterization of two Debaryomyces spp. yeast isolates obtained from a Korean traditional fermented soy sauce "ganjang". Both Debaryomyces hansenii ganjang isolates KD2 and C11 showed halotolerance to concentrations of up to 15% NaCl and improved growth in the presence of salt. Ploidy and whole-genome sequencing analyses indicated that the KD2 genome is haploid, whereas the C11 genome is heterozygous diploid with two distinctive subgenomes. Interestingly, phylogenetic analysis using intron sequences indicated that the C11 strain was generated via hybridization between D. hansenii and D. tyrocola ancestor strains. The D. hansenii KD2 and D. hansenii-hybrid C11 produced various volatile flavor compounds associated with butter, caramel, cheese, and fruits, and showed high bioconversion activity from ferulic acid to 4-vinylguaiacol, a characteristic flavor compound of soybean products. Both KD2 and C11 exhibited viability in the presence of bile salts and at low pH and showed immunomodulatory activity to induce high levels of the anti-inflammatory cytokine IL-10. The safety of the yeast isolates was confirmed by analyzing virulence and acute oral toxicity. Together, the D. hansenii ganjang isolates possess physiological properties beneficial for improving the flavor and nutritional value of fermented products.

Gene Coexpression Connectivity Predicts Gene Targets Underlying High Ionic-Liquid Tolerance in Yarrowia lipolytica

Microbial tolerance to organic solvents such as ionic liquids (ILs) is a robust phenotype beneficial for novel biotransformation. While most microbes become inhibited in 1% to 5% (vol/vol) IL (e.g., 1-ethyl-3-methylimidazolium acetate), we engineered a robust Yarrowia lipolytica strain (YlCW001) that tolerates a record high of 18% (vol/vol) IL via adaptive laboratory evolution. Yet, genotypes conferring high IL tolerance in YlCW001 remain to be discovered. In this study, we shed light on the underlying cellular processes that enable robust Y. lipolytica to thrive in inhibitory ILs. By using dynamic transcriptome sequencing (RNA-Seq) data, we introduced Gene Coexpression Connectivity (GeCCo) as a metric to discover genotypes conferring desirable phenotypes that might not be found by the conventional differential expression (DE) approaches. GeCCo selects genes based on their number of coexpressed genes in a subnetwork of upregulated genes by the target phenotype. We experimentally validated GeCCo by reverse engineering a high-IL-tolerance phenotype in wild-type Y. lipolytica. We found that gene targets selected by both DE and GeCCo exhibited the best statistical chance at increasing IL tolerance when individually overexpressed. Remarkably, the best combination of dual-overexpression genes was genes selected by GeCCo alone. This nonintuitive combination of genes, BRN1 and OYE2, is involved in guiding/regulating mitotic cell division, chromatin segregation/condensation, microtubule and cytoskeletal organization, and Golgi vesicle transport.

IMPORTANCE Cellular robustness to cope with stressors is an important phenotype. Y. lipolytica is an industrial robust oleaginous yeast that has recently been discovered to tolerate record high concentrations of ILs, beneficial for novel biotransformation in organic solvents. However, genotypes that link to IL tolerance in Y. lipolytica are largely unknown. Due to the complex IL-tolerant phenotype, conventional gene discovery and validation based on differential gene expression approaches are time-consuming due to a large search space and might encounter a high false-discovery rate. Here, using the developed Gene Coexpression Connectivity (GeCCo) method, we identified and validated a subset of most promising gene targets conferring the IL-tolerant phenotypes and shed light on their potential mechanisms. We anticipate GeCCo being a useful method to discover the genotype-to-phenotype link.

Evolution of holocentric chromosomes: Drivers, diversity, and deterrents

Centromeres are specialized chromosomal regions that recruit kinetochore proteins and mediate spindle microtubule attachment to ensure faithful chromosome segregation during mitosis and meiosis. Centromeres can be restricted to one region of the chromosome. Named "monocentromere", this type represents the most commonly found centromere organization across eukaryotes. Alternatively, centromeres can also be assembled at sites chromosome-wide. This second type is called "holocentromere". Despite their early description over 100 years ago, research on holocentromeres has lagged behind that of monocentromeres. Nevertheless, the application of next generation sequencing approaches and advanced microscopic technologies enabled recent advances understanding the molecular organization and regulation of holocentromeres in different organisms. Here we review the current state of research on holocentromeres focusing on evolutionary considerations. First, we provide a brief historical perspective on the discovery of holocentric chromosomes. We then discuss models/drivers that have been proposed over the years to explain the evolutionary transition from mono- to holocentric chromosomes. We continue to review the description of holocentric chromosomes in diverse eukaryotic groups and then focus our discussion on a specific and recently characterized type of holocentromere organization in insects that functions independently of the otherwise essential centromeric marker protein CenH3, thus providing novel insights into holocentromere evolution in insects. Finally, we propose reasons to explain why the holocentric trait is not more frequent across eukaryotes despite putative selective advantages.

A roadmap to understanding diversity and function of coral reef-associated fungi

Tropical coral reefs are hotspots of marine productivity, owing to the association of reef-building corals with endosymbiotic algae and metabolically diverse bacterial communities. However, the functional importance of fungi, well-known for their contribution to shaping terrestrial ecosystems and global nutrient cycles, remains underexplored on coral reefs. We here conceptualize how fungal functional traits may have facilitated the spread, diversification, and ecological adaptation of marine fungi on coral reefs. We propose that functions of reef-associated fungi may be diverse and go beyond their hitherto described roles of pathogens and bioeroders, including but not limited to reef-scale biogeochemical cycles and the structuring of coral-associated and environmental microbiomes via chemical mediation. Recent technological and conceptual advances will allow the elucidation of the physiological, ecological, and chemical contributions of understudied marine fungi to coral holobiont and reef ecosystem functioning and health and may help provide an outlook for reef management actions.