Rapid expansion and functional divergence of subtelomeric gene families in yeasts

Computational analysis of polymorphic residues in maltose and maltotriose transporters of a wild Saccharomyces cerevisiae strain

The metabolism of maltose and maltotriose, the primary sugars in brewing wort, depends on an efficient transport system. However, most Saccharomyces cerevisiae strains transport maltotriose inefficiently, leaving residual α-glucosides in the final product. Proteins involved in maltotriose transport exhibit diverse polymorphic sequences linked to sugar transport efficiency. In this study, a wild S. cerevisiae strain was placed under adaptive selection, resulting in a strain with a 65 and 44% increase in maltose and maltotriose transport rates, respectively. Genes encoding maltose and maltotriose transporters, including MALx1, MPHx, and AGT1, were detected in both the native and adapted strains. One variant of Mal31p, carrying a polymorphism at position 371 in transmembrane helix 7, was identified. This helix has been reported to have a high likelihood of undergoing polymorphisms. Bioinformatics analysis revealed structural changes affecting substrate interactions and channel dynamics, with the polymorphism conferring greater protein flexibility and reducing electrostatic interactions. These results suggest that the residue at position 371 in maltose and maltotriose transporters is a key element distinct from those previously reported. Additionally, we propose a significant set of polymorphic residues within these transporters potentially resulting from the evolution of these proteins.

Genomic investigation of plant secondary metabolism: insights from synteny network analysis of oxidosqualene cyclase flanking genes

The clustered distribution of genes involved in metabolic pathways within the plant genome has garnered significant attention from researchers. By comparing and analyzing changes in the flanking regions of metabolic genes across a diverse array of species, we can enhance our understanding of the formation and distribution of biosynthetic gene clusters (BGCs). In this study, we have designed a workflow that uncovers and assesses conserved positional relationships between genes in various species by using synteny neighborhood networks (SNN). This workflow is then applied to the analysis of flanking genes associated with oxidosqualene cyclases (OSCs). The method allows for the recognition and comparison of homologous blocks with unique flanking genes accompanying different subfamilies of OSCs. The examination of the flanking genes of OSCs in 122 plant species revealed multiple genes with conserved positional relationships with OSCs in angiosperms. Specifically, the earliest adjacency of OSC genes and CYP716 genes first appeared in basal eudicots, and the nonrandom occurrence of CYP716 genes in the flanking region of OSC persists across different lineages of eudicots. Our study showed the substitution of genes in the flanking region of the OSC varies across different plant lineages, and our approach facilitates the investigation of flanking gene rearrangements in the formation of OSC-related BGCs.

SHIP identifies genomic safe harbors in eukaryotic organisms using genomic general feature annotation

Integrating foreign genes into loci, allowing their transcription without affecting endogenous gene expression, is the desirable strategy in genomic engineering. However, these loci, known as genomic safe harbors (GSHs), have been mainly identified by empirical methods. Furthermore, the most prominent available GSHs are localized within regions of high gene density, raising concerns about unstable expression. As synthetic biology is moving towards investigating polygenic modules rather than single genes, there is an increasing demand for tools to identify GSHs systematically. To expand the GSH repertoire, we present SHIP, an algorithm designed to detect potential GSHs in eukaryotes. Using the chassis organism Saccharomyces cerevisiae, five GSHs were experimentally curated based on data from DNA sequencing, stability, flow cytometry, qPCR, electron microscopy, RT-qPCR, and RNA-Seq assays. Our study places SHIP as a valuable tool for providing a list of promising candidates to assist in the experimental assessment of GSHs in eukaryotic organisms with available annotated genomes.

The type of carbon source not the growth rate it supports can determine diauxie in Saccharomyces cerevisiae

How cells choose between carbon sources is a classic example of cellular decision-making. Microbes often prioritise glucose, but there has been little investigation of whether other sugars are also preferred. Here we study budding yeast growing on mixtures of sugars with palatinose, a sucrose isomer that cells catabolise with the MAL regulon. We find that the decision-making involves more than carbon flux-sensing: yeast prioritise galactose over palatinose, but sucrose and fructose weakly if at all despite each allowing faster growth than palatinose. With genetic perturbations and transcriptomics, we show that the regulation is active with repression of the MAL genes via Gal4, the GAL regulon's master regulator. We argue, using mathematical modelling, that cells enforce their preference for galactose through weakening the MAL regulon's positive feedback. They do so through decreasing intracellular palatinose by repressing MAL11, the palatinose transporter, and expressing the isomaltases IMA1 and IMA5. Supporting these predictions, we show that deleting IMA1 abolishes diauxie. Our results demonstrate that budding yeast actively prioritises carbon sources other than glucose and that such priorities need not reflect differences in growth rates. They imply that carbon-sensing strategies even in model organisms are more complex than previously thought.

Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa

The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.

Chromosome-level subgenome-aware de novo assembly provides insight into Saccharomyces bayanus genome divergence after hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaptation, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitate adaptation, we have completed chromosome-level de novo assemblies of all chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380, using Oxford Nanopore Technologies' MinION long-read sequencing. We characterize the S. bayanus genome and compare it with its parent species, Saccharomyces uvarum and Saccharomyces eubayanus, and other S. bayanus genomes to better understand genome evolution after a relatively recent hybridization event. We observe multiple recombination events between the subgenomes in each chromosome, followed by loss of heterozygosity (LOH) in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of Saccharomyces cerevisiae, such as those involved in the maltose metabolism. Finally, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaptation to maltose brewing environments and mitigated the maladaptive effect of hybridization. This paper describes the first in-depth study using long-read sequencing technology of an S. bayanus hybrid genome, which may serve as an excellent reference for future studies of this important yeast and other yeast strains.

Characterization and organization of telomeric-linked helicase (tlh) gene families in Fusarium oxysporum

Telomere-linked RecQ helicase (tlh) genes have been reported in several fungi and a choanoflagellate in the regions adjacent to the terminal telomere repeats. In this study, we explored the Telomere-linked RecQ helicase (tlh) genes in four strains of Fusarium oxysporum, offering new insights into their genomic structure, functional motifs, and impact on chromosomal ends. We conducted a comprehensive analysis, comparing the tlh genes of F. oxysporum with those previously identified in other organisms and uncovering significant similarities. Through comparative genomics, we identified conserved protein motifs across these genes, including a TLH domain, C2H2, and RecQ helicase motifs. Our phylogenetic analysis positions the F. oxysporum tlh genes in a cluster with other known tlhs, suggesting a shared evolutionary origin. Mutation analysis revealed a relatively low level of deleterious mutations in tlh gene paralogs, with a notable proportion of full-size structures maintained across strains. Analysis of subtelomeric sequences indicates that a region with almost identical sequences flanks the majority of chromosome ends, termed tlh-containing region (TLHcr), across these strains. The presence of TLHcrs at chromosome ends, either as single entities or in arrays, underscores their potential role in telomere function and genome stability. Our findings provide a detailed examination of tlh genes in four strains of F. oxysporum, laying the groundwork for future studies on their biological significance and evolutionary history in fungal genomes.

Evaluating cellular roles and phenotypes associated with trehalose degradation genes in Saccharomyces cerevisiae

In the yeast Saccharomyces cerevisiae, 2 types of trehalase activities have been described. Neutral trehalases (Nth1 and Nth2) are considered to be the main proteins that catalyze intracellular trehalose mobilization. In addition to Nth1 and Nth2, studies have shown that acid trehalase Ath1 is required for extracellular trehalose degradation. Although both neutral and acid-type trehalases have been predominantly investigated in laboratory strains of S. cerevisiae, we sought to examine the phenotypic consequences of disrupting these genes in wild strains. In this study, we constructed mutants of the trehalose degradation pathway (NTH1, NTH2, and ATH1) in 5 diverse S. cerevisiae strains to examine whether published lab strain phenotypes are also exhibited by wild strains. For each mutant, we assessed a number of phenotypes for comparison to trehalose biosynthesis mutants, including trehalose production, glycogen production, cell size, acute thermotolerance, high-temperature growth, sporulation efficiency, and growth on a variety of carbon sources in rich and minimal medium. We found that all trehalase mutants including single deletion nth1Δ, nth2Δ, and ath1Δ, as well as double deletion nth1nth2Δ, accumulated higher intracellular trehalose levels compared to their isogenic wild-type cells. Also, nth1Δ and nth1Δnth2Δ mutants exhibited mild thermal sensitivity, suggesting a potential minor role for trehalose mobilization when cells recover from stress. In addition, we evaluated phenotypes more directly relevant to trehalose degradation, including both extracellular and intracellular trehalose utilization. We discovered that intracellular trehalose hydrolysis is critical for typical spore germination progression, highlighting a role for trehalose in cell cycle regulation, likely as a storage carbohydrate providing glycolytic fuel. Additionally, our work provides further evidence suggesting Ath1 is indispensable for extracellular trehalose utilization as a carbon source, even in the presence of AGT1.

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ∼400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

Selection for robust metabolism in domesticated yeasts is driven by adaptation to Hsp90 stress

Protein folding both promotes and constrains adaptive evolution. We uncover this surprising duality in the role of the protein-folding chaperone heat shock protein 90 (Hsp90) in maintaining the integrity of yeast metabolism amid proteotoxic stressors within industrial domestication niches. Ethanol disrupts critical Hsp90-dependent metabolic pathways and exerts strong selective pressure for redundant duplications of key genes within these pathways, yielding the classical genomic signatures of beer and bread domestication. This work demonstrates a mechanism of adaptive canalization in an ecology of major economic importance and highlights Hsp90-dependent variation as an important source of phantom heritability in complex traits.

Expression patterns of Mal genes and association with differential maltose and maltotriose transport rate of two Saccharomyces pastorianus yeasts

Beer brewing is a well-known process that still faces great challenges, such as the total consumption of sugars present in the fermentation media. Lager-style beer, a major worldwide beer type, is elaborated by Saccharomyces pastorianus (Sp) yeast, which must ferment high maltotriose content worts, but its consumption represents a notable problem, especially among Sp strains belonging to group I. Factors, such as fermentation conditions, presence of maltotriose transporters, transporter copy number variation, and genetic regulation variations contribute to this issue. We assess the factors affecting fermentation in two Sp yeast strains: SpIB1, with limited maltotriose uptake, and SpIB2, known for efficient maltotriose transport. Here, SpIB2 transported significantly more maltose (28%) and maltotriose (32%) compared with SpIB1. Furthermore, SpIB2 expressed all MAL transporters (ScMALx1, SeMALx1, ScAGT1, SeAGT1, MTT1, and MPHx) on the first day of fermentation, whereas SpIB1 only exhibited ScMalx1, ScAGT1, and MPH2/3 genes. Some SpIB2 transporters had polymorphic transmembrane domains (TMD) resembling MTT1, accompanied by higher expression of these transporters and its positive regulator genes, such as MAL63. These findings suggest that, in addition to the factors mentioned above, positive regulators of Mal transporters contribute significantly to phenotypic diversity in maltose and maltotriose consumption among the studied lager yeast strains.IMPORTANCEBeer, the third most popular beverage globally with a 90% market share in the alcoholic beverage industry, relies on Saccharomyces pastorianus (Sp) strains for lager beer production. These strains exhibit phenotypic diversity in maltotriose consumption, a crucial process for the acceptable organoleptic profile in lager beer. This diversity ranges from Sp group II strains with a notable maltotriose-consuming ability to Sp group I strains with limited capacity. Our study highlights that differential gene expression of maltose and maltotriose transporters and its upstream trans-elements, such as MAL gene-positive regulators, adds complexity to this variation. This insight can contribute to a more comprehensive analysis needed to the development of controlled and efficient biotechnological processes in the beer brewing industry.

Specialization restricts the evolutionary paths available to yeast sugar transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning approximately 400 million years of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

Mechanisms of adaptation and evolution in Toxoplasma gondii

Toxoplasma has high host flexibility, infecting all nucleated cells of mammals and birds. This implies that during its infective process the parasite must constantly adapt to different environmental situations, which in turn leads to modifications in its metabolism, regulation of gene transcription, translation of mRNAs and stage specific factors. There are conserved pathways that support these adaptations, which we aim to elucidate in this review. We begin by exploring the widespread epigenetic mechanisms and transcription regulators, continue with the supportive role of Heat Shock Proteins (Hsp), the translation regulation, stress granules, and finish with the emergence of contingency genes in highly variable genomic domains, such as subtelomeres. Within epigenetics, the discovery of a new histone variant of the H2B family (H2B.Z), contributing to T. gondii virulence and differentiation, but also gene expression regulation and its association with the metabolic state of the parasite, is highlighted. Associated with the regulation of gene expression are transcription factors (TFs). An overview of the main findings on TF and development is presented. We also emphasize the role of Hsp90 and Tgj1 in T. gondii metabolic fitness and the regulation of protein translation. Translation regulation is also highlighted as a mechanism for adaptation to conditions encountered by the parasite as well as stress granules containing mRNA and proteins generated in the extracellular tachyzoite. Another important aspect in evolution and adaptability are the subtelomeres because of their high variability and gene duplication rate. Toxoplasma possess multigene families of membrane proteins and contingency genes that are associated with different metabolic stresses. Among them parasite differentiation and environmental stresses stand out, including those that lead tachyzoite to bradyzoite conversion. Finally, we are interested in positioning protozoa as valuable evolution models, focusing on research related to the Extended Evolutionary Synthesis, based on models recently generated, such as extracellular adaptation and ex vivo cyst recrudescence.

Mutational robustness and the role of buffer genes in evolvability

Organisms rely on mutations to fuel adaptive evolution. However, many mutations impose a negative effect on fitness. Cells may have therefore evolved mechanisms that affect the phenotypic effects of mutations, thus conferring mutational robustness. Specifically, so-called buffer genes are hypothesized to interact directly or indirectly with genetic variation and reduce its effect on fitness. Environmental or genetic perturbations can change the interaction between buffer genes and genetic variation, thereby unmasking the genetic variation's phenotypic effects and thus providing a source of variation for natural selection to act on. This review provides an overview of our understanding of mutational robustness and buffer genes, with the chaperone gene HSP90 as a key example. It discusses whether buffer genes merely affect standing variation or also interact with de novo mutations, how mutational robustness could influence evolution, and whether mutational robustness might be an evolved trait or rather a mere side-effect of complex genetic interactions.

Structural characterization of a polymorphic repeat at the CACNA1C schizophrenia locus

Genetic variation within intron 3 of the CACNA1C calcium channel gene is associated with schizophrenia and bipolar disorder, but analysis of the causal variants and their effect is complicated by a nearby variable-number tandem repeat (VNTR). Here, we used 155 long-read genome assemblies from 78 diverse individuals to delineate the structure and population variability of the CACNA1C intron 3 VNTR. We categorized VNTR sequences into 7 Types of structural alleles using sequence differences among repeat units. Only 12 repeat units at the 5' end of the VNTR were shared across most Types, but several Types were related through a series of large and small duplications. The most diverged Types were rare and present only in individuals with African ancestry, but the multiallelic structural polymorphism Variable Region 2 was present across populations at different frequencies, consistent with expansion of the VNTR preceding the emergence of early hominins. VR2 was in complete linkage disequilibrium with fine-mapped schizophrenia variants (SNPs) from genome-wide association studies (GWAS). This risk haplotype was associated with decreased CACNA1C gene expression in brain tissues profiled by the GTEx project. Our work suggests that sequence variation within a human-specific VNTR affects gene expression, and provides a detailed characterization of new alleles at a flagship neuropsychiatric locus.

Machine learning enables identification of an alternative yeast galactose utilization pathway

How genomic differences contribute to phenotypic differences is a major question in biology. The recently characterized genomes, isolation environments, and qualitative patterns of growth on 122 sources and conditions of 1,154 strains from 1,049 fungal species (nearly all known) in the yeast subphylum Saccharomycotina provide a powerful, yet complex, dataset for addressing this question. We used a random forest algorithm trained on these genomic, metabolic, and environmental data to predict growth on several carbon sources with high accuracy. Known structural genes involved in assimilation of these sources and presence/absence patterns of growth in other sources were important features contributing to prediction accuracy. By further examining growth on galactose, we found that it can be predicted with high accuracy from either genomic (92.2%) or growth data (82.6%) but not from isolation environment data (65.6%). Prediction accuracy was even higher (93.3%) when we combined genomic and growth data. After the GALactose utilization genes, the most important feature for predicting growth on galactose was growth on galactitol, raising the hypothesis that several species in two orders, Serinales and Pichiales (containing the emerging pathogen Candida auris and the genus Ogataea, respectively), have an alternative galactose utilization pathway because they lack the GAL genes. Growth and biochemical assays confirmed that several of these species utilize galactose through an alternative oxidoreductive D-galactose pathway, rather than the canonical GAL pathway. Machine learning approaches are powerful for investigating the evolution of the yeast genotype-phenotype map, and their application will uncover novel biology, even in well-studied traits.

Identification of an active RNAi pathway in Candida albicans

RNA interference (RNAi) is a fundamental regulatory pathway with a wide range of functions, including regulation of gene expression and maintenance of genome stability. Although RNAi is widespread in the fungal kingdom, well-known species, such as the model yeast Saccharomyces cerevisiae, have lost the RNAi pathway. Until now evidence has been lacking for a fully functional RNAi pathway in Candida albicans, a human fungal pathogen considered critically important by the World Health Organization. Here, we demonstrated that the widely used C. albicans reference strain (SC5314) contains an inactivating missense mutation in the gene encoding for the central RNAi component Argonaute. In contrast, most other C. albicans isolates contain a canonical Argonaute protein predicted to be functional and RNAi-active. Indeed, using high-throughput small and long RNA sequencing combined with seamless CRISPR/Cas9-based gene editing, we demonstrate that an active C. albicans RNAi machinery represses expression of subtelomeric gene families. Thus, an intact and functional RNAi pathway exists in C. albicans, highlighting the importance of using multiple reference strains when studying this dangerous pathogen.

The Complete Genomic Sequence of Microbial Transglutaminase Producer, Streptomyces mobaraensis DSM40587

Actinomycetes are remarkable natural sources of active natural molecules and enzymes of considerable industrial value. Streptomyces mobaraensis is the first microorganism found to produce transglutaminase with broad industrial applications. Although transglutaminase in S. mobaraensis has been well studied over the past three decades, the genome of S. mobaraensis and its secondary metabolic potential were poorly reported. Here, we presented the complete genome of S. mobaraensis DSM40587 obtained from the German Collection of Microorganisms and Cell Cultures GmbH. It contains a linear chromosome of 7,633,041 bp and a circular plasmid of 23,857 bp. The chromosome with an average GC content of 73.49% was predicted to harbour 6683 protein-coding genes, seven rRNA and 69 tRNA genes. Comparative genomic analysis reveals its meaningful genomic characterisation. A comprehensive bioinformatics investigation identifies 35 putative BGCs (biosynthesis gene clusters) involved in synthesising various secondary metabolites. Of these, 13 clusters showed high similarity (> 55%) to known BGCs coding for polyketides, nonribosomal peptides, hopene, RiPP (Ribosomally synthesized and post-translationally modified peptides), and others. Furthermore, these BGCs with over 65% similarity to the known BGCs were analysed in detail. The complete genome of S. mobaraensis DSM40587 reveals its capacity to yield diverse bioactive natural products and provides additional insights into discovering novel secondary metabolites.

Chromosome-level Subgenome-aware de novo Assembly of Saccharomyces bayanus Provides Insight into Genome Divergence after Hybridization

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaption, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitated adaptation, we completed chromosome-level de novo assemblies of all 16 pairs chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380 (IFO11022), using Nanopore MinION long-read sequencing. Characterization of S. bayanus subgenomes and comparative analysis with the genomes of its parent species, S. uvarum and S. eubayanus, provide several new insights into understanding genome evolution after a relatively recent hybridization. For instance, multiple recombination events between the two subgenomes have been observed in each chromosome, followed by loss of heterozygosity (LOH) in most chromosomes in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of S. cerevisiae, such as those involved in the maltose metabolism. In addition, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus. The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaption to maltose brewing environments and mitigated the maladaptive effect of hybridization.

Physiological role of the EHL gene in sake yeast and its effects on quality of sake

The EHL1/2/3 genes were identified by whole-genome sequencing of Kyokai No. 7 (K7), which is a well-known representative Japanese sake yeast Saccharomyces cerevisiae. The genes are present in K7, but not in laboratory strain S288C. Although the genes were presumed to encode epoxide hydrolase based on homology analysis, their effect on cellular metabolism in sake yeast has not yet been clarified. We constructed ehl1/2/3 mutants harboring a stop codon in each gene using the haploid yeast strain H3 as the parental strain, which was derived from K701, and investigated the physiological role and effects of the EHL1/2/3 genes on sake quality. Metabolome analysis and vitamin requirement testing revealed that the EHL1/2/3 genes are partly responsible for the synthesis of pantothenate. For fermentation profiles, ethanol production by the ehl1/2/3 mutant was comparable with that of strain H3, but succinate production was decreased in the ehl1/2/3 mutant compared to strain H3 when cultured in yeast malt (YM) medium containing 10% glucose and during sake brewing. Ethyl hexanoate and isoamyl acetate levels in the ehl1/2/3 mutant strain were decreased compared to those of strain H3 during sake brewing. Thus, the EHL1/2/3 genes did not affect ethanol production but did affect the production of organic acids and aromatic components during sake brewing.

Ploidy evolution in a wild yeast is linked to an interaction between cell type and metabolism

Ploidy is an evolutionarily labile trait, and its variation across the tree of life has profound impacts on evolutionary trajectories and life histories. The immediate consequences and molecular causes of ploidy variation on organismal fitness are frequently less clear, although extreme mating type skews in some fungi hint at links between cell type and adaptive traits. Here, we report an unusual recurrent ploidy reduction in replicate populations of the budding yeast Saccharomyces eubayanus experimentally evolved for improvement of a key metabolic trait, the ability to use maltose as a carbon source. We find that haploids have a substantial, but conditional, fitness advantage in the absence of other genetic variation. Using engineered genotypes that decouple the effects of ploidy and cell type, we show that increased fitness is primarily due to the distinct transcriptional program deployed by haploid-like cell types, with a significant but smaller contribution from absolute ploidy. The link between cell-type specification and the carbon metabolism adaptation can be traced to the noncanonical regulation of a maltose transporter by a haploid-specific gene. This study provides novel mechanistic insight into the molecular basis of an environment-cell type fitness interaction and illustrates how selection on traits unexpectedly linked to ploidy states or cell types can drive karyotypic evolution in fungi.

Evolution and functional diversification of yeast sugar transporters

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

Improved Sugarcane-Based Fermentation Processes by an Industrial Fuel-Ethanol Yeast Strain

In Brazil, sucrose-rich broths (cane juice and/or molasses) are used to produce billions of liters of both fuel ethanol and cachaça per year using selected Saccharomyces cerevisiae industrial strains. Considering the important role of feedstock (sugar) prices in the overall process economics, to improve sucrose fermentation the genetic characteristics of a group of eight fuel-ethanol and five cachaça industrial yeasts that tend to dominate the fermentors during the production season were determined by array comparative genomic hybridization. The widespread presence of genes encoding invertase at multiple telomeres has been shown to be a common feature of both baker's and distillers' yeast strains, and is postulated to be an adaptation to sucrose-rich broths. Our results show that only two strains (one fuel-ethanol and one cachaça yeast) have amplification of genes encoding invertase, with high specific activity. The other industrial yeast strains had a single locus (SUC2) in their genome, with different patterns of invertase activity. These results indicate that invertase activity probably does not limit sucrose fermentation during fuel-ethanol and cachaça production by these industrial strains. Using this knowledge, we changed the mode of sucrose metabolism of an industrial strain by avoiding extracellular invertase activity, overexpressing the intracellular invertase, and increasing its transport through the AGT1 permease. This approach allowed the direct consumption of the disaccharide by the cells, without releasing glucose or fructose into the medium, and a 11% higher ethanol production from sucrose by the modified industrial yeast, when compared to its parental strain.

Ethanol Drives Evolution of Hsp90-Dependent Robustness by Redundancy in Yeast Domestication

Protein folding promotes and constrains adaptive evolution. We uncover this surprising duality in the role the protein-folding chaperone Hsp90 plays in mediating the interplay between proteome and the genome which acts to maintain the integrity of yeast metabolism in the face of proteotoxic stressors in anthropic niches. Of great industrial relevance, ethanol concentrations generated by fermentation in the making of beer and bread disrupt critical Hsp90-dependent nodes of metabolism and exert strong selective pressure for increased copy number of key genes encoding components of these nodes, yielding the classical genetic signatures of beer and bread domestication. This work establishes a mechanism of adaptive canalization in an ecology of major economic significance and highlights Hsp90-contingent variation as an important source of phantom heritability in complex traits.

Nuclear genome organization in fungi: from gene folding to Rabl chromosomes

Comparative genomics has recently provided unprecedented insights into the biology and evolution of the fungal lineage. In the postgenomics era, a major research interest focuses now on detailing the functions of fungal genomes, i.e. how genomic information manifests into complex phenotypes. Emerging evidence across diverse eukaryotes has revealed that the organization of DNA within the nucleus is critically important. Here, we discuss the current knowledge on the fungal genome organization, from the association of chromosomes within the nucleus to topological structures at individual genes and the genetic factors required for this hierarchical organization. Chromosome conformation capture followed by high-throughput sequencing (Hi-C) has elucidated how fungal genomes are globally organized in Rabl configuration, in which centromere or telomere bundles are associated with opposite faces of the nuclear envelope. Further, fungal genomes are regionally organized into topologically associated domain-like (TAD-like) chromatin structures. We discuss how chromatin organization impacts the proper function of DNA-templated processes across the fungal genome. Nevertheless, this view is limited to a few fungal taxa given the paucity of fungal Hi-C experiments. We advocate for exploring genome organization across diverse fungal lineages to ensure the future understanding of the impact of nuclear organization on fungal genome function.

Dynamics of synthetic yeast chromosome evolution shaped by hierarchical chromatin organization

Synthetic genome evolution provides a dynamic approach for systematically and straightforwardly exploring evolutionary processes. Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is an evolutionary system intrinsic to the synthetic yeast genome that can rapidly drive structural variations. Here, we detect over 260 000 rearrangement events after the SCRaMbLEing of a yeast strain harboring 5.5 synthetic yeast chromosomes (synII, synIII, synV, circular synVI, synIXR and synX). Remarkably, we find that the rearrangement events exhibit a specific landscape of frequency. We further reveal that the landscape is shaped by the combined effects of chromatin accessibility and spatial contact probability. The rearrangements tend to occur in 3D spatially proximal and chromatin-accessible regions. The enormous numbers of rearrangements mediated by SCRaMbLE provide a driving force to potentiate directed genome evolution, and the investigation of the rearrangement landscape offers mechanistic insights into the dynamics of genome evolution.

Highly complete long-read genomes reveal pangenomic variation underlying yeast phenotypic diversity

Understanding the genetic causes of trait variation is a primary goal of genetic research. One way that individuals can vary genetically is through variable pangenomic genes: genes that are only present in some individuals in a population. The presence or absence of entire genes could have large effects on trait variation. However, variable pangenomic genes can be missed in standard genotyping workflows, owing to reliance on aligning short-read sequencing to reference genomes. A popular method for studying the genetic basis of trait variation is linkage mapping, which identifies quantitative trait loci (QTLs), regions of the genome that harbor causative genetic variants. Large-scale linkage mapping in the budding yeast Saccharomyces cerevisiae has found thousands of QTLs affecting myriad yeast phenotypes. To enable the resolution of QTLs caused by variable pangenomic genes, we used long-read sequencing to generate highly complete de novo genome assemblies of 16 diverse yeast isolates. With these assemblies, we resolved QTLs for growth on maltose, sucrose, raffinose, and oxidative stress to specific genes that are absent from the reference genome but present in the broader yeast population at appreciable frequency. Copies of genes also duplicate onto chromosomes where they are absent in the reference genome, and we found that these copies generate additional QTLs whose resolution requires pangenome characterization. Our findings show the need for highly complete genome assemblies to identify the genetic basis of trait variation.

Structural Variations and Adaptations of Synthetic Chromosome Ends Driven by SCRaMbLE in Haploid and Diploid Yeasts

Variations and adaptations of chromosome ends play an important role in eukaryotic karyotype evolution. Traditional experimental studies of the adaptations of chromosome ends mainly rely on the strategy of introducing defects; thus, the adaptation methods of survivors may vary depending on the initial defects. Here, using the SCRaMbLE strategy, we obtained a library of haploid and diploid synthetic strains with variations in chromosome ends. Analysis of the SCRaMbLEd survivors revealed four routes of adaptation: homologous recombination between nonhomologous chromosome arms (haploids) or homologous chromosome arms (diploids), site-specific recombination between intra- or interchromosomal ends, circularization of chromosomes, and loss of whole chromosomes (diploids). We also found that circularization of synthetic chromosomes can be generated by SCRaMbLE. Our study of various adaptation routes of chromosome ends provides insight into eukaryotic karyotype evolution from the viewpoint of synthetic genomics.

TERRA and Telomere Maintenance in the Yeast Saccharomyces cerevisiae

Telomeres are structures made of DNA, proteins and RNA found at the ends of eukaryotic linear chromosomes. These dynamic nucleoprotein structures protect chromosomal tips from end-to-end fusions, degradation, activation of damage checkpoints and erroneous DNA repair events. Telomeres were thought to be transcriptionally silent regions because of their constitutive heterochromatin signature until telomeric long non-coding RNAs (LncRNAs) were discovered. One of them, TERRA (TElomeric Repeat-containing RNA), starts in the subtelomeric regions towards the chromosome ends from different telomeres and has been extensively studied in many evolutionarily distant eukaryotes. Changes in TERRA’s expression can lead to telomeric dysfunction, interfere with the replicative machinery and impact telomere length. TERRA also co-localizes in vivo with telomerase, and can form RNA:DNA hybrid structures called R-loops, which have been implicated in the onset of senescence and the alternative lengthening of telomere (ALT) pathway. Yet, the molecular mechanisms involving TERRA, as well as its function, remain elusive. Here, we review the current knowledge of TERRA transcription, structure, expression, regulation and its multiple telomeric and extra-telomeric functions in the budding yeast Saccharomyces cerevisiae.

Late Embryogenesis Abundant Proteins Contribute to the Resistance of Toxoplasma gondii Oocysts against Environmental Stresses

Toxoplasma gondii oocysts, which are shed in large quantities in the feces from infected felines, are very stable in the environment, resistant to most inactivation procedures, and highly infectious. The oocyst wall provides an important physical barrier for sporozoites contained inside oocysts, protecting them from many chemical and physical stressors, including most inactivation procedures. Furthermore, sporozoites can withstand large temperature changes, even freeze-thawing, as well as desiccation, high salinity, and other environmental insults; however, the genetic basis for this environmental resistance is unknown. Here, we show that a cluster of four genes encoding Late Embryogenesis Abundant (LEA)-related proteins are required to provide Toxoplasma sporozoites resistance to environmental stresses. Toxoplasma LEA-like genes (TgLEAs) exhibit the characteristic features of intrinsically disordered proteins, explaining some of their properties. Our in vitro biochemical experiments using recombinant TgLEA proteins show that they have cryoprotective effects on the oocyst-resident lactate dehydrogenase enzyme and that induced expression in E. coli of two of them leads to better survival after cold stress. Oocysts from a strain in which the four LEA genes were knocked out en bloc were significantly more susceptible to high salinity, freezing, and desiccation compared to wild-type oocysts. We discuss the evolutionary acquisition of LEA-like genes in Toxoplasma and other oocyst-producing apicomplexan parasites of the Sarcocystidae family and discuss how this has likely contributed to the ability of sporozoites within oocysts to survive outside the host for extended periods. Collectively, our data provide a first molecular detailed view on a mechanism that contributes to the remarkable resilience of oocysts against environmental stresses. IMPORTANCE Toxoplasma gondii oocysts are highly infectious and may survive in the environment for years. Their resistance against disinfectants and irradiation has been attributed to the oocyst and sporocyst walls by acting as physical and permeability barriers. However, the genetic basis for their resistance against stressors like changes in temperature, salinity, or humidity, is unknown. We show that a cluster of four genes encoding Toxoplasma Late Embryogenesis Abundant (TgLEA)-related proteins are important for this resistance to environmental stresses. TgLEAs have features of intrinsically disordered proteins, explaining some of their properties. Recombinant TgLEA proteins show cryoprotective effects on the parasite's lactate dehydrogenase, an abundant enzyme in oocysts, and expression in E. coli of two TgLEAs has a beneficial effect on growth after cold stress. Moreover, oocysts from a strain lacking all four TgLEA genes were more susceptible to high salinity, freezing, and desiccation compared to wild-type oocysts, highlighting the importance of the four TgLEAs for oocyst resilience.

Functional diversity and plasticity in the sugar preferences of Saccharomyces MALT transporters in domesticated yeasts

Maltose and maltotriose, together with glucose, are the major carbohydrates found in malts. Thus, brewing yeasts grown in malt-based brewing processes with serial re-pitching have likely increased their ability to uptake these sugars during domestication by modulating the expression and copy number of maltose transporter genes (MALT, also known as Malx1). However, the molecular basis for and structural insights into the sugar preferences of MALT proteins remain to be elucidated. Here we report the functional evaluation of two novel Saccharomyces cerevisiae MALT proteins, ScMalt#2p and ScMalt#5p, from industrial brewing yeasts, focusing on their maltose and maltotriose preferences. Structural models of the MALT proteins generated by AlphaFold2 and functional analyses of substitution mutants revealed that a very small number of amino acid residues in two spatially adjacent transmembrane helixes, TMH7 and TMH11, appear to be crucial for sugar preference. Thus, subtle conformational alterations conferred by a small number of amino acid polymorphisms within MALTs would contribute to the adaptation of domesticated brewing yeasts to the constrained carbohydrate environment of industrial wort during beer brewing.

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.

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

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

Untargeted metabolomics screening reveals unique secondary metabolite production from Alternaria section Alternaria

Alternaria section Alternaria is comprised of many species that infect a broad diversity of important crop plants and cause post-harvest spoilage. Alternaria section Alternaria species, such as A. alternata and A. arborescens, are prolific producers of secondary metabolites that act as virulence factors of disease and are mycotoxins that accumulate in infected tissues-metabolites that can vary in their spectrum of production between individuals from the same fungal species. Untargeted metabolomics profiling of secondary metabolite production using mass spectrometry is an effective means to detect phenotypic anomalies in secondary metabolism within a species. Secondary metabolite phenotypes from 36 Alternaria section Alternaria isolates were constructed to observe frequency of production patterns. A clear and unique mass feature pattern was observed for three of the strains that were linked with the production of the dehydrocurvularin family of toxins and associated detoxification products. Examination of corresponding genomes revealed the presence of the dehydrocurvularin biosynthesis gene cluster associated with a sub-telomeric accessory region. A comparison of sequence similarity and occurrences of the dehydrocurvularin biosynthetic gene cluster within Pleosporalean fungi is presented and discussed.

Adaptive sequence divergence forged new neurodevelopmental enhancers in humans

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

The pan-genome of Aspergillus fumigatus provides a high-resolution view of its population structure revealing high levels of lineage-specific diversity driven by recombination

Aspergillus fumigatus is a deadly agent of human fungal disease where virulence heterogeneity is thought to be at least partially structured by genetic variation between strains. While population genomic analyses based on reference genome alignments offer valuable insights into how gene variants are distributed across populations, these approaches fail to capture intraspecific variation in genes absent from the reference genome. Pan-genomic analyses based on de novo assemblies offer a promising alternative to reference-based genomics with the potential to address the full genetic repertoire of a species. Here, we evaluate 260 genome sequences of A. fumigatus including 62 newly sequenced strains, using a combination of population genomics, phylogenomics, and pan-genomics. Our results offer a high-resolution assessment of population structure and recombination frequency, phylogenetically structured gene presence-absence variation, evidence for metabolic specificity, and the distribution of putative antifungal resistance genes. Although A. fumigatus disperses primarily via asexual conidia, we identified extraordinarily high levels of recombination with the lowest linkage disequilibrium decay value reported for any fungal species to date. We provide evidence for 3 primary populations of A. fumigatus, with recombination occurring only rarely between populations and often within them. These 3 populations are structured by both gene variation and distinct patterns of gene presence-absence with unique suites of accessory genes present exclusively in each clade. Accessory genes displayed functional enrichment for nitrogen and carbohydrate metabolism suggesting that populations may be stratified by environmental niche specialization. Similarly, the distribution of antifungal resistance genes and resistance alleles were often structured by phylogeny. Altogether, the pan-genome of A. fumigatus represents one of the largest fungal pan-genomes reported to date including many genes unrepresented in the Af293 reference genome. These results highlight the inadequacy of relying on a single-reference genome-based approach for evaluating intraspecific variation and the power of combined genomic approaches to elucidate population structure, genetic diversity, and putative ecological drivers of clinically relevant fungi.

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

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

Functional Divergence in a Multi-gene Family Is a Key Evolutionary Innovation for Anaerobic Growth in Saccharomyces cerevisiae

The amplification and diversification of genes into large multi-gene families often mark key evolutionary innovations, but this process often creates genetic redundancy that hinders functional investigations. When the model budding yeast Saccharomyces cerevisiae transitions to anaerobic growth conditions, the cell massively induces the expression of seven serine/threonine-rich anaerobically-induced cell wall mannoproteins (anCWMPs): TIP1, TIR1, TIR2, TIR3, TIR4, DAN1, and DAN4. Here, we show that these genes likely derive evolutionarily from a single ancestral anCWMP locus, which was duplicated and translocated to new genomic contexts several times both prior to and following the budding yeast whole genome duplication (WGD) event. Based on synteny and their phylogeny, we separate the anCWMPs into four gene subfamilies. To resolve prior inconclusive genetic investigations of these genes, we constructed a set of combinatorial deletion mutants to determine their contributions toward anaerobic growth in S. cerevisiae. We found that two genes, TIR1 and TIR3, were together necessary and sufficient for the anCWMP contribution to anaerobic growth. Overexpressing either gene alone was insufficient for anaerobic growth, implying that they encode non-overlapping functional roles in the cell during anaerobic growth. We infer from the phylogeny of the anCWMP genes that these two important genes derive from an ancient duplication that predates the WGD event, whereas the TIR1 subfamily experienced gene family amplification after the WGD event. Taken together, the genetic and molecular evidence suggests that one key anCWMP gene duplication event, several auxiliary gene duplication events, and functional divergence underpin the evolution of anaerobic growth in budding yeasts.

Top-Down, Knowledge-Based Genetic Reduction of Yeast Central Carbon Metabolism

Saccharomyces cerevisiae, whose evolutionary past includes a whole-genome duplication event, is characterized by a mosaic genome configuration with substantial apparent genetic redundancy. This apparent redundancy raises questions about the evolutionary driving force for genomic fixation of “minor” paralogs and complicates modular and combinatorial metabolic engineering strategies. While isoenzymes might be important in specific environments, they could be dispensable in controlled laboratory or industrial contexts. The present study explores the extent to which the genetic complexity of the central carbon metabolism (CCM) in S. cerevisiae, here defined as the combination of glycolysis, the pentose phosphate pathway, the tricarboxylic acid cycle, and a limited number of related pathways and reactions, can be reduced by elimination of (iso)enzymes without major negative impacts on strain physiology. Cas9-mediated, groupwise deletion of 35 of the 111 genes yielded a “minimal CCM” strain which, despite the elimination of 32% of CCM-related proteins, showed only a minimal change in phenotype on glucose-containing synthetic medium in controlled bioreactor cultures relative to a congenic reference strain. Analysis under a wide range of other growth and stress conditions revealed remarkably few phenotypic changes from the reduction of genetic complexity. Still, a well-documented context-dependent role of GPD1 in osmotolerance was confirmed. The minimal CCM strain provides a model system for further research into genetic redundancy of yeast genes and a platform for strategies aimed at large-scale, combinatorial remodeling of yeast CCM.

Architecture, Chromatin and Gene Organization of Toxoplasma gondii Subtelomeres

Subtelomeres (ST) are chromosome regions that separate telomeres from euchromatin and play relevant roles in various biological processes of the cell. While their functions are conserved, ST structure and genetic compositions are unique to each species. This study aims to identify and characterize the subtelomeric regions of the 13 Toxoplasma gondii chromosomes of the Me49 strain. Here, STs were defined at chromosome ends based on poor gene density. The length of STs ranges from 8.1 to 232.4 kbp, with a gene density of 0.049 genes/kbp, lower than the Me49 genome (0.15 kbp). Chromatin organization showed that H3K9me3, H2A.X, and H3.3 are highly enriched near telomeres and the 5' end of silenced genes, decaying in intensity towards euchromatin. H3K4me3 and H2A.Z/H2B.Z are shown to be enriched in the 5' end of the ST genes. Satellite DNA was detected in almost all STs, mainly the sat350 family and a novel satellite named sat240. Beyond the STs, only short dispersed fragments of sat240 and sat350 were found. Within STs, there were 12 functional annotated genes, 59 with unknown functions (Hypothetical proteins), 15 from multigene FamB, and 13 from multigene family FamC. Some genes presented low interstrain synteny associated with the presence of satellite DNA. Orthologues of FamB and FamC were also detected in Neospora caninum and Hammondia hammondi. A re-analysis of previous transcriptomic data indicated that ST gene expression is strongly linked to the adaptation to different situations such as extracellular passage (evolve and resequencing study) and changes in metabolism (lack of acetyl-CoA cofactor). In conclusion, the ST region of the T. gondii chromosomes was defined, the STs genes were determined, and it was possible to associate them with high interstrain plasticity and a role in the adaptability of T. gondii to environmental changes.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Pseudo-Lager—Brewing with Lutra® Kveik Yeast

Brewers commonly produce ales since the ale yeast is more resilient, ferments quicker and requires higher temperatures, which are easier to ensure as opposed to lager and pilsner, which require lower temperatures and longer lagering time. However, Kveik yeasts are also resilient, ferment at fairly high temperatures (up to 35 °C), and can provide light, lager-like beers, but more quickly, in shorter lagering time, and with reduced off flavors. Diacetyl rest is not needed. The intention of this paper was to assess the possibility of producing pseudo-lager by using Lutra® Kveik. A batch (120 L) was divided into six fermenting vessels and inoculated with Lutra® yeast. To test its possibility to result in lager-like beer at higher temperature, we conducted fermentation at two temperatures (21 and 35 °C). Fermentation subjected to 21 °C lasted for 9 days, while at 35 °C, fermentation was finished in 2 days. After fermentation, both beers were stored in cold temperatures (4 °C) and then kegged, carbonized, and analyzed (pH, ethanol, polyphenols, color, bitterness, clarity). Alongside the sensory evaluation, a GC-MS analysis was also conducted in order to determine if there are any difference between the samples.

Intracellular trehalose accumulation via the Agt1 transporter promotes freeze-thaw tolerance in Saccharomyces cerevisiae

AIM

This study is to investigate the use of a constitutively expressed trehalose transport protein to directly control intracellular trehalose levels and protect baker's yeast (Saccharomyces cerevisiae) cells against freeze-thaw stress in vivo.

METHODS AND RESULTS

We used a constitutively overexpressed Agt1 transporter system to investigate the role of trehalose in the freeze-thaw tolerance of yeast cells by regulating intracellular trehalose concentrations independently of intracellular biosynthesis. Using this method, we found that increasing intracellular trehalose in yeast cells improved cell survival rate after 8 days of freezing at -80°C and -20°C. We also observed that freeze-thaw tolerance promoted by intracellular trehalose only occurs in highly concentrated cell pellets rather than cells in liquid suspension.

CONCLUSIONS

Trehalose is sufficient to provide freeze-thaw tolerance using our Agt1 overexpression system. Freeze-thaw tolerance can be further enhanced by deletion of genes encoding intracellular trehalose degradation enzymes.

SIGNIFICANCE AND IMPACT OF STUDY

These findings are relevant to improving the freeze-thaw tolerance of baker's yeast in the frozen baked goods industry through engineering strains that can accumulate intracellular trehalose via a constitutively expressed trehalose transporter and inclusion of trehalose into the growth medium.

The teenage years of yeast population genomics trace history, admixing and getting wilder

Population genomics studies the evolutionary processes that shape intraspecies genetic variations. In this review, I explore the insights into yeast-population genomics that have emerged from recent advances in sequencing. Genomes of the model Saccharomyces cerevisiae and many new yeast species from around the world are being used to address various aspects of population biology, including geographical origin, the level of introgression, domestication signatures, and outcrossing frequency. New long-read sequencing has enabled a greater capacity to quantify these variations at a finer resolution from complete de novo genomes at the population scale to phasing subgenomes of different origins. These resources provide a platform to dissect the relationship between phenotypes across environmental niches.

Maltose-Negative Yeast in Non-Alcoholic and Low-Alcoholic Beer Production

Although beer is a widely used beverage in many cultures, there is a need for a new drinking alternative in the face of rising issues such as health concerns or weight problems. However, non-alcoholic and low-alcoholic beers (NABLAB) still have some sensory problems that have not been fully remedied today, such as “wort-like”/”potato-like” flavours or a lack of aroma. These defects are due to the lack of alcohol (and the lack of the aldehyde-reducing effect of alcohol fermentation), as well as production techniques. The use of new yeast strains that cannot ferment maltose—the foremost sugar in the wort—is highly promising to produce a more palatable and sustainable NABLAB product because production with these yeast strains can be performed with standard brewery equipment. In the scientific literature, it is clear that interest in the production of NABLAB has increased recently, and experiments have been carried out with maltose-negative yeast strains isolated from many different environments. This study describes maltose-negative yeasts and their aromatic potential for the production of NABLAB by comprehensively examining recent academic studies.

An Integrated Approach Reveals DNA Damage and Proteotoxic Stress as Main Effects of Proton Radiation in S. cerevisiae

Proton radiotherapy (PRT) has the potential to reduce the normal tissue toxicity associated with conventional photon-based radiotherapy (X-ray therapy, XRT) because the active dose can be more directly targeted to a tumor. Although this dosimetric advantage of PRT is well known, the molecular mechanisms affected by PRT remain largely elusive. Here, we combined the molecular toolbox of the eukaryotic model Saccharomyces cerevisiae with a systems biology approach to investigate the physiological effects of PRT compared to XRT. Our data show that the DNA damage response and protein stress response are the major molecular mechanisms activated after both PRT and XRT. However, RNA-Seq revealed that PRT treatment evoked a stronger activation of genes involved in the response to proteotoxic stress, highlighting the molecular differences between PRT and XRT. Moreover, inhibition of the proteasome resulted in decreased survival in combination with PRT compared to XRT, not only further confirming that protons induced a stronger proteotoxic stress response, but also hinting at the potential of using proteasome inhibitors in combination with proton radiotherapy in clinical settings.

Genome-wide effect of non-optimal temperatures under anaerobic conditions on gene expression in Saccharomyces cerevisiae

Understanding of thermal adaptation mechanisms in yeast is crucial to develop better-adapted strains to industrial processes, providing more economical and sustainable products. We have analyzed the transcriptomic responses of three Saccharomyces cerevisiae strains, a commercial wine strain, ADY5, a laboratory strain, CEN.PK113-7D and a commercial bioethanol strain, Ethanol Red, grown at non-optimal temperatures under anaerobic chemostat conditions. Transcriptomic analysis of the three strains revealed a huge complexity of cellular mechanisms and responses. Overall, cold exerted a stronger transcriptional response in the three strains comparing with heat conditions, with a higher number of down-regulating genes than of up-regulating genes regardless the strain analyzed. The comparison of the transcriptome at both sub- and supra-optimal temperatures showed the presence of common genes up- or down-regulated in both conditions, but also the presence of common genes up- or down-regulated in the three studied strains. More specifically, we have identified and validated three up-regulated genes at sub-optimal temperature in the three strains, OPI3, EFM6 and YOL014W. Finally, the comparison of the transcriptomic data with a previous proteomic study with the same strains revealed a good correlation between gene activity and protein abundance, mainly at low temperature. Our work provides a global insight into the specific mechanisms involved in temperature adaptation regarding both transcriptome and proteome, which can be a step forward in the comprehension and improvement of yeast thermotolerance.

Full-Length Genome of an Ogataea polymorpha Strain CBS4732 ura3Δ Reveals Large Duplicated Segments in Subtelomeric Regions

BACKGROUND

Currently, methylotrophic yeasts (e.g., Pichia pastoris, Ogataea polymorpha, and Candida boindii) are subjects of intense genomics studies in basic research and industrial applications. In the genus Ogataea, most research is focused on three basic O. polymorpha strains-CBS4732, NCYC495, and DL-1. However, the relationship between CBS4732, NCYC495, and DL-1 remains unclear, as the genomic differences between them have not be exactly determined without their high-quality complete genomes. As a nutritionally deficient mutant derived from CBS4732, the O. polymorpha strain CBS4732 ura3Δ (named HU-11) is being used for high-yield production of several important proteins or peptides. HU-11 has the same reference genome as CBS4732 (noted as HU-11/CBS4732), because the only genomic difference between them is a 5-bp insertion.

RESULTS

In the present study, we have assembled the full-length genome of O. polymorpha HU-11/CBS4732 using high-depth PacBio and Illumina data. Long terminal repeat retrotransposons (LTR-rts), rDNA, 5' and 3' telomeric, subtelomeric, low complexity and other repeat regions were exactly determined to improve the genome quality. In brief, the main findings include complete rDNAs, complete LTR-rts, three large duplicated segments in subtelomeric regions and three structural variations between the HU-11/CBS4732 and NCYC495 genomes. These findings are very important for the assembly of full-length genomes of yeast and the correction of assembly errors in the published genomes of Ogataea spp. HU-11/CBS4732 is so phylogenetically close to NCYC495 that the syntenic regions cover nearly 100% of their genomes. Moreover, HU-11/CBS4732 and NCYC495 share a nucleotide identity of 99.5% through their whole genomes. CBS4732 and NCYC495 can be regarded as the same strain in basic research and industrial applications.

CONCLUSION

The present study preliminarily revealed the relationship between CBS4732, NCYC495, and DL-1. Our findings provide new opportunities for in-depth understanding of genome evolution in methylotrophic yeasts and lay the foundations for the industrial applications of O. polymorpha CBS4732, NCYC495, DL-1, and their derivative strains. The full-length genome of O. polymorpha HU-11/CBS4732 should be included into the NCBI RefSeq database for future studies of Ogataea spp.

Domestication reprogrammed the budding yeast life cycle

Domestication of plants and animals is the foundation for feeding the world human population but can profoundly alter the biology of the domesticated species. Here we investigated the effect of domestication on one of our prime model organisms, the yeast Saccharomyces cerevisiae, at a species-wide level. We tracked the capacity for sexual and asexual reproduction and the chronological life span across a global collection of 1,011 genome-sequenced yeast isolates and found a remarkable dichotomy between domesticated and wild strains. Domestication had systematically enhanced fermentative and reduced respiratory asexual growth, altered the tolerance to many stresses and abolished or impaired the sexual life cycle. The chronological life span remained largely unaffected by domestication and was instead dictated by clade-specific evolution. We traced the genetic origins of the yeast domestication syndrome using genome-wide association analysis and genetic engineering and disclosed causative effects of aneuploidy, gene presence/absence variations, copy number variations and single-nucleotide polymorphisms. Overall, we propose domestication to be the most dramatic event in budding yeast evolution, raising questions about how much domestication has distorted our understanding of the natural biology of this key model species.

Characteristics of New England India Pale Ale Beer Produced with the Use of Norwegian KVEIK Yeast

The aim of this research was to determine the potential of four unconventional Norwegian yeasts of the KVEIK type to produce NEIPA beer. The influence of yeast strains on fermentation process, physicochemical properties, antioxidant potential, volatile compounds, and sensory properties was investigated. The KVEIK-fermented beer did not differ in terms of physicochemical parameters from the beer produced with the commercial variants of US-05 yeast. The yeast strain influenced the sensory quality (taste and aroma) of the beers, with KVEIK-fermented beer rating significantly higher. The antioxidant activity of the tested beers also significantly depended on the yeast strain applied. The beers fermented with KVEIK had a significantly higher antioxidant potential (ABTS^•+) than those fermented with US-05. The strongest antioxidant activity was found in the beer brewed with the Lida KVEIK yeast. The use of KVEIK to produce NEIPA beer allowed enrichment of the finished products with volatile compounds isobutanol, 2-pentanol, 3-methylobutanol, ethyl octanoate, and ethyl decanoate.