AGAPE (Automated Genome Analysis PipelinE) for pan-genome analysis of Saccharomyces cerevisiae

Footprints of Human Migration in the Population Structure of Wild Baker's Yeast

Humans have a long history of fermenting food and beverages that led to domestication of the baker's yeast, Saccharomyces cerevisiae. Despite their tight companionship with humans, yeast species that are domesticated or pathogenic can also live on trees. Here we used over 300 genomes of S. cerevisiae from oaks and other trees to determine whether tree-associated populations are genetically distinct from domesticated lineages and estimate the timing of forest lineage divergence. We found populations on trees are highly structured within Europe, Japan, and North America. Approximate estimates of when forest lineages diverged out of Asia and into North America and Europe coincide with the end of the last ice age, the spread of agriculture, and the onset of fermentation by humans. It appears that migration from human-associated environments to trees is ongoing. Indeed, patterns of ancestry in the genomes of three recent migrants from the trees of North America to Europe could be explained by the human response to the Great French Wine Blight. Our results suggest that human-assisted migration affects forest populations, albeit rarely. Such migration events may even have shaped the global distribution of S. cerevisiae. Given the potential for lasting impacts due to yeast migration between human and natural environments, it seems important to understand the evolution of human commensals and pathogens in wild niches.

Isolation and characterization of haploid heterothallic beer yeasts

Improving ale or lager yeasts by conventional breeding is a non-trivial task. Domestication of lager yeasts, which are hybrids between Saccharomyces cerevisiae and Saccharomyces eubayanus, has led to evolved strains with severely reduced or abolished sexual reproduction capabilities, due to, e.g. postzygotic barriers. On the other hand, S. cerevisiae ale yeasts, particularly Kveik ale yeast strains, were shown to produce abundant viable spores (~ 60%; Dippel et al. Microorganisms 10(10):1922, 2022). This led us to investigate the usefulness of Kveik yeasts for conventional yeast breeding. Surprisingly, we could isolate heterothallic colonies from germinated spores of different Kveik strains. These strains presented stable mating types in confrontation assays with pheromone-sensitive tester strains. Heterothallism was due to inactivating mutations in their HO genes. These led to amino acid exchanges in the Ho protein, revealing a known G223D mutation and also a novel G217R mutation, both of which abolished mating type switching. We generated stable MATa or MATα lines of four different Kveik yeasts, named Odin, Thor, Freya and Vör. Analyses of bud scar positions in these strains revealed both axial and bipolar budding patterns. However, the ability of Freya and Vör to form viable meiotic offspring with haploid tester strains demonstrated that these strains are haploid. Fermentation analyses indicated that all four yeast strains were able to ferment maltose and maltotriose. Odin was found to share not only mutations in the HO gene, but also inactivating mutations in the PAD1 and FDC1 genes with lager yeasts, which makes this strain POF-, i.e. not able to generate phenolic off-flavours, a key feature of lager yeasts. These haploid ale yeast-derived strains may open novel avenues also for generating novel lager yeast strains by breeding or mutation and selection utilizing the power of yeast genetics, thus lifting a block that domestication of lager yeasts has brought about. KEY POINTS: • Haploid Kveik ale yeasts with stable MATa and MATα mating types were isolated. • Heterothallic strains bear mutant HO alleles leading to a novel inactivating G217R amino acid change. • One strain was found to be POF- due to inactivating mutations in the PAD1 and FDC1 gene rendering it negative for phenolic off-flavor production. • These strains are highly accessible for beer yeast improvements by conventional breeding, employing yeast genetics and mutation and selection regimes.

Horizontal Transfer and Recombination Fuel Ty4 Retrotransposon Evolution in Saccharomyces

Horizontal transposon transfer (HTT) plays an important role in the evolution of eukaryotic genomes; however, the detailed evolutionary history and impact of most HTT events remain to be elucidated. To better understand the process of HTT in closely related microbial eukaryotes, we studied Ty4 retrotransposon subfamily content and sequence evolution across the genus Saccharomyces using short- and long-read whole genome sequence data, including new PacBio genome assemblies for two Saccharomyces mikatae strains. We find evidence for multiple independent HTT events introducing the Tsu4 subfamily into specific lineages of Saccharomyces paradoxus, Saccharomyces cerevisiae, Saccharomyces eubayanus, Saccharomyces kudriavzevii and the ancestor of the S. mikatae/Saccharomyces jurei species pair. In both S. mikatae and S. kudriavzevii, we identified novel Ty4 clades that were independently generated through recombination between resident and horizontally transferred subfamilies. Our results reveal that recurrent HTT and lineage-specific extinction events lead to a complex pattern of Ty4 subfamily content across the genus Saccharomyces. Moreover, our results demonstrate how HTT can lead to coexistence of related retrotransposon subfamilies in the same genome that can fuel evolution of new retrotransposon clades via recombination.

Overview of the Saccharomyces cerevisiae population structure through the lens of 3,034 genomes

With the rise of high-throughput sequencing technologies, a holistic view of genetic variation within populations-through population genomics studies-appears feasible, although it remains an ongoing effort. Genetic variation arises from a diverse range of evolutionary forces, with mutation and recombination being key drivers in shaping genomes. Studying genetic variation within a population represents a crucial first step in understanding the relationship between genotype and phenotype and the evolutionary history of species. In this context, the budding yeast Saccharomyces cerevisiae has been at the forefront of population genomic studies. In addition, it has a complex history that involves adaptation to a wide range of wild and human-related ecological niches. Although to date more than 3,000 diverse isolates have been sequenced, there is currently a lack of a resource bringing together sequencing data and associated metadata for all sequenced isolates. To perform a comprehensive analysis of the population structure of S. cerevisiae, we collected genome sequencing data from 3,034 natural isolates and processed the data uniformly. We determined ploidy levels, identified single nucleotide polymorphisms (SNPs), small insertion-deletions (InDels), copy number variations (CNVs), and aneuploidies across the population, creating a publicly accessible resource for the yeast research community. Interestingly, we showed that this population captures ∼93% of the species diversity. Using neighbor-joining and Bayesian methods, we redefined the populations, revealing clustering patterns primarily based on ecological origin. This work represents a valuable resource for the community and efforts have been made to make it evolvable and integrable to future yeast population studies.

The Saccharomyces cerevisiae cell lysis mutant cly8 is a temperature-sensitive allele of ESP1

The Saccharomyces cerevisiae mutants designated cly (cell lysis) cause cell lysis at elevated temperatures. The cly8 mutation, previously localized to an 80 kilobase region between GCD2 and SPT6 on the right arm of chromosome VII, has been used for mutation mapping and in recombination assays, but its genetic identity has remained unknown. Whole genome sequencing of cly8 mutant and CLY8 wild-type strains revealed four missense mutations specific to the cly8 mutant strain in the GCD2-SPT6 interval, three of these missense mutations affected essential genes. The esp1-G543E mutation, but not the two other missense mutations, recapitulated the temperature-sensitive growth, the cell lysis phenotype, and recessive nature of the cly8 mutation. The ESP1 gene encodes S. cerevisiae separase and is required for normal chromosome segregation during mitosis. Shifting cells containing the esp1-G543E mutation to the non-permissive temperature caused chromosome missegregation based on the analysis of DNA content by flow cytometry analysis. Taken together, these results indicate that the temperature-sensitive cly8 mutation is an allele of the essential ESP1 gene.

Histone deacetylase Hos2 regulates protein expression noise by potentially modulating the protein translation machinery

Non-genetic variations derived from expression noise at transcript or protein levels can result in cell-to-cell heterogeneity within an isogenic population. Although cells have developed strategies to reduce noise in some cellular functions, this heterogeneity can also facilitate varying levels of regulation and provide evolutionary benefits in specific environments. Despite several general characteristics of cellular noise having been revealed, the detailed molecular pathways underlying noise regulation remain elusive. Here, we established a dual-fluorescent reporter system in Saccharomyces cerevisiae and performed experimental evolution to search for mutations that increase expression noise. By analyzing evolved cells using bulk segregant analysis coupled with whole-genome sequencing, we identified the histone deacetylase Hos2 as a negative noise regulator. A hos2 mutant down-regulated multiple ribosomal protein genes and exhibited partially compromised protein translation, indicating that Hos2 may regulate protein expression noise by modulating the translation machinery. Treating cells with translation inhibitors or introducing mutations into several Hos2-regulated ribosomal protein genes-RPS9A, RPS28B and RPL42A-enhanced protein expression noise. Our study provides an effective strategy for identifying noise regulators and also sheds light on how cells regulate non-genetic variation through protein translation.

Horizontal transfer and recombination fuel Ty4 retrotransposon evolution in Saccharomyces

Horizontal transposon transfer (HTT) plays an important role in the evolution of eukaryotic genomes, however the detailed evolutionary history and impact of most HTT events remain to be elucidated. To better understand the process of HTT in closely-related microbial eukaryotes, we studied Ty4 retrotransposon subfamily content and sequence evolution across the genus Saccharomyces using short- and long-read whole genome sequence data, including new PacBio genome assemblies for two S. mikatae strains. We find evidence for multiple independent HTT events introducing the Tsu4 subfamily into specific lineages of S. paradoxus, S. cerevisiae, S. eubayanus, S. kudriavzevii and the ancestor of the S. mikatae/S. jurei species pair. In both S. mikatae and S. kudriavzevii, we identified novel Ty4 clades that were independently generated through recombination between resident and horizontally-transferred subfamilies. Our results reveal that recurrent HTT and lineage-specific extinction events lead to a complex pattern of Ty4 subfamily content across the genus Saccharomyces. Moreover, our results demonstrate how HTT can lead to coexistence of related retrotransposon subfamilies in the same genome that can fuel evolution of new retrotransposon clades via recombination.

Differential Hsp90-dependent gene expression is strain-specific and common among yeast strains

Enhanced phenotypic diversity increases a population's likelihood of surviving catastrophic conditions. Hsp90, an essential molecular chaperone and a central network hub in eukaryotes, has been observed to suppress or enhance the effects of genetic variation on phenotypic diversity in response to environmental cues. Because many Hsp90-interacting genes are involved in signaling transduction pathways and transcriptional regulation, we tested how common Hsp90-dependent differential gene expression is in natural populations. Many genes exhibited Hsp90-dependent strain-specific differential expression in five diverse yeast strains. We further identified transcription factors (TFs) potentially contributing to variable expression. We found that on Hsp90 inhibition or environmental stress, activities or abundances of Hsp90-dependent TFs varied among strains, resulting in differential strain-specific expression of their target genes, which consequently led to phenotypic diversity. We provide evidence that individual strains can readily display specific Hsp90-dependent gene expression, suggesting that the evolutionary impacts of Hsp90 are widespread in nature.

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.

Overexpression profiling reveals cellular requirements in the context of genetic backgrounds and environments

Overexpression can help life adapt to stressful environments, making an examination of overexpressed genes valuable for understanding stress tolerance mechanisms. However, a systematic study of genes whose overexpression is functionally adaptive (GOFAs) under stress has yet to be conducted. We developed a new overexpression profiling method and systematically identified GOFAs in Saccharomyces cerevisiae under stress (heat, salt, and oxidative). Our results show that adaptive overexpression compensates for deficiencies and increases fitness under stress, like calcium under salt stress. We also investigated the impact of different genetic backgrounds on GOFAs, which varied among three S. cerevisiae strains reflecting differing calcium and potassium requirements for salt stress tolerance. Our study of a knockout collection also suggested that calcium prevents mitochondrial outbursts under salt stress. Mitochondria-enhancing GOFAs were only adaptive when adequate calcium was available and non-adaptive when calcium was deficient, supporting this idea. Our findings indicate that adaptive overexpression meets the cell's needs for maximizing the organism's adaptive capacity in the given environment and genetic context.

Strain-dependent differences in coordination of yeast signalling networks

The yeast mitogen-activated protein kinase pathways serve as a model system for understanding how network interactions affect the way in which cells coordinate the response to multiple signals. We have quantitatively compared two yeast strain backgrounds YPH499 and ∑1278b (both of which have previously been used to study these pathways) and found several important differences in how they coordinate the interaction between the high osmolarity glycerol (HOG) and mating pathways. In the ∑1278b background, in response to simultaneous stimulus, mating pathway activation is dampened and delayed in a dose-dependent manner. In the YPH499 background, only dampening is dose-dependent. Furthermore, leakage from the HOG pathway into the mating pathway (crosstalk) occurs during osmostress alone in the ∑1278b background only. The mitogen-activated protein kinase Hog1p suppresses crosstalk late in an induction time course in both strains but does not affect the early crosstalk seen in the ∑1278b background. Finally, the kinase Rck2p plays a greater role suppressing late crosstalk in the ∑1278b background than in the YPH499 background. Our results demonstrate that comparisons between laboratory yeast strains provide an important resource for understanding how signalling network interactions are tuned by genetic variation without significant alteration to network structure.

Foci-forming regions of pyruvate kinase and enolase at the molecular surface incorporate proteins into yeast cytoplasmic metabolic enzymes transiently assembling (META) bodies

Spatial reorganization of metabolic enzymes to form the "metabolic enzymes transiently assembling (META) body" is increasingly recognized as a mechanism contributing to regulation of cellular metabolism in response to environmental changes. A number of META body-forming enzymes, including enolase (Eno2p) and phosphofructokinase, have been shown to contain condensate-forming regions. However, whether all META body-forming enzymes have condensate-forming regions or whether enzymes have multiple condensate-forming regions remains unknown. The condensate-forming regions of META body-forming enzymes have potential utility in the creation of artificial intracellular enzyme assemblies. In the present study, the whole sequence of yeast pyruvate kinase (Cdc19p) was searched for condensate-forming regions. Four peptide fragments comprising 27-42 amino acids were found to form condensates. Together with the fragment previously identified from Eno2p, these peptide regions were collectively termed "META body-forming sequences (METAfos)." METAfos-tagged yeast alcohol dehydrogenase (Adh1p) was found to co-localize with META bodies formed by endogenous Cdc19p under hypoxic conditions. The effect of Adh1p co-localization with META bodies on cell metabolism was further evaluated. Expression of Adh1p fused with a METAfos-tag increased production of ethanol compared to acetic acid, indicating that spatial reorganization of metabolic enzymes affects cell metabolism. These results contribute to understanding of the mechanisms and biological roles of META body formation.

Understanding the Role of Yeast Yme1 in Mitochondrial Function Using Biochemical and Proteomics Analyses

Mitochondrial i-AAA proteinase Yme1 is a multifunctional protein that plays important roles in maintaining mitochondrial protein homeostasis and regulating biogenesis and function of mitochondrial proteins. However, due to the complex interplay of mitochondria and the multifunctional nature of Yme1, how Yme1 affects mitochondrial function and protein homeostasis is still poorly understood. In this study, we investigated how YME1 deletion affects yeast Saccharomyces cerevisiae growth, chronological life span, mitochondrial protein homeostasis and function, with a focus on the mitochondrial oxidative phosphorylation (OXPHOS) complexes. Our results show that whilst the YME1 deleted cells grow poorly under respiratory conditions, they grow similar to wild-type yeast under fermentative conditions. However, the chronological life span is impaired, indicating that Yme1 plays a key role in longevity. Using highly enriched mitochondrial extract and proteomic analysis, we show that the abundances of many mitochondrial proteins are altered by YME1 deletion. Several components of the respiratory chain complexes II, III, IV and V were significantly decreased, suggesting that Yme1 plays an important role in maintaining the level and function of complexes II-V. This result was confirmed using blue native-PAGE and in-solution-based enzyme activity assays. Taken together, this study shows that Yme1 plays an important role in the chronological life span and mitochondrial protein homeostasis and has deciphered its function in maintaining the activity of mitochondrial OXPHOS complexes.

New data and collaborations at the Saccharomyces Genome Database: updated reference genome, alleles, and the Alliance of Genome Resources

Saccharomyces cerevisiae is used to provide fundamental understanding of eukaryotic genetics, gene product function, and cellular biological processes. Saccharomyces Genome Database (SGD) has been supporting the yeast research community since 1993, serving as its de facto hub. Over the years, SGD has maintained the genetic nomenclature, chromosome maps, and functional annotation, and developed various tools and methods for analysis and curation of a variety of emerging data types. More recently, SGD and six other model organism focused knowledgebases have come together to create the Alliance of Genome Resources to develop sustainable genome information resources that promote and support the use of various model organisms to understand the genetic and genomic bases of human biology and disease. Here we describe recent activities at SGD, including the latest reference genome annotation update, the development of a curation system for mutant alleles, and new pages addressing homology across model organisms as well as the use of yeast to study human disease.

Exceptional origin activation revealed by comparative analysis in two laboratory yeast strains

We performed a comparative analysis of replication origin activation by genome-wide single-stranded DNA mapping in two yeast strains challenged by hydroxyurea, an inhibitor of the ribonucleotide reductase. We gained understanding of the impact on origin activation by three factors: S-phase checkpoint control, DNA sequence polymorphisms, and relative positioning of origin and transcription unit. Wild type W303 showed a significant reduction of fork progression accompanied by an elevated level of Rad53 phosphorylation as well as physical presence at origins compared to A364a. Moreover, a rad53K227A mutant in W303 activated more origins, accompanied by global reduction of ssDNA across all origins, compared to A364a. Sequence polymorphism in the consensus motifs of origins plays a minor role in determining strain-specific activity. Finally, we identified a new class of origins only active in checkpoint-proficient cells, which we named “Rad53-dependent origins”. Our study presents a comprehensive list of differentially used origins and provide new insights into the mechanisms of origin activation.

Integrated multi-omics reveals common properties underlying stress granule and P-body formation

Non-membrane-bound compartments such as P-bodies (PBs) and stress granules (SGs) play important roles in the regulation of gene expression following environmental stresses. We have systematically and quantitatively determined the protein and mRNA composition of PBs and SGs formed before and after nutrient stress. We find that high molecular weight (HMW) complexes exist prior to glucose depletion that we propose may act as seeds for further condensation of proteins forming mature PBs and SGs. We identify an enrichment of proteins with low complexity and RNA binding domains, as well as long, structured mRNAs that are poorly translated following nutrient stress. Many proteins and mRNAs are shared between PBs and SGs including several multivalent RNA binding proteins that promote condensate interactions during liquid-liquid phase separation. We uncover numerous common protein and RNA components across PBs and SGs that support a complex interaction profile during the maturation of these biological condensates. These interaction networks represent a tuneable response to stress, highlighting previously unrecognized condensate heterogeneity. These studies therefore provide an integrated and quantitative understanding of the dynamic nature of key biological condensates.

Draft genome sequence and potential identification of a biosurfactant from Brevibacterium casei strain LS14 an isolate from fresh water Loktak Lake

This study reports the whole-genome sequencing and sequence analysis of a bacterial isolate Brevibacterium casei strain LS14, isolated from Loktak Lake, Imphal, India. The de novo assembled genome reported in this paper featured a size of 3,809,532 bp, has GC content of 68% and contains 3602 genomic features, including 3551 protein-coding genes, 46 tRNA and 5rRNA. A biosurfactant biosynthesis gene cluster in the genome of the isolated strain was identified using AntiSMASH online tool V3.0.5 and KAAS (KEGG Automatic Annotation Server). The presence of biosurfactant was demonstrated by drop collapse, oil displacement and emulsification index. Subsequent chemical characterization using FTIR and LC-MS analyses revealed surfactin and terpene containing biosurfactant moieties. Also, the presence of genes involved in terpenoid synthesis pathway in the genome sequence may account for biosurfactant terpenoid backbone, but genes for later-stage conversion of terpenoid to biosurfactant were not ascertained.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02867-9.

Design of an improved universal signal peptide based on the α-factor mating secretion signal for enzyme production in yeast

Saccharomyces cerevisiae plays an important role in the heterologous expression of an array of proteins due to its easy manipulation, low requirements and ability for protein post-translational modifications. The implementation of the preproleader secretion signal of the α-factor mating pheromone from this yeast contributes to increase the production yields by targeting the foreign protein to the extracellular environment. The use of this signal peptide combined with enzyme-directed evolution allowed us to achieve the otherwise difficult functional expression of fungal laccases in S. cerevisiae, obtaining different evolved α-factor preproleader sequences that enhance laccase secretion. However, the design of a universal signal peptide to enhance the production of heterologous proteins in S. cerevisiae is a pending challenge. We describe here the optimisation of the α-factor preproleader to improve recombinant enzyme production in S. cerevisiae through two parallel engineering strategies: a bottom-up design over the native α-factor preproleader (αnat) and a top-down design over the fittest evolved signal peptide obtained in our lab (α9H2 leader). The goal was to analyse the effect of mutations accumulated in the signal sequence throughout iterations of directed evolution, or of other reported mutations, and their possible epistatic interactions. Both approaches agreed in the positive synergism of four mutations (Aα9D, Aα20T, Lα42S, Dα83E) contained in the final optimised leader (αOPT), which notably enhanced the secretion of several fungal oxidoreductases and hydrolases. Additionally, we suggest a guideline to further drive the heterologous production of a particular enzyme based on combinatorial saturation mutagenesis of positions 86th and 87th of the αOPT leader fused to the target protein.

The pan-genome of Saccharomyces cerevisiae

Understanding genotype-phenotype relationship is fundamental in biology. With the benefit from next-generation sequencing and high-throughput phenotyping methodologies, there have been generated much genome and phenome data for Saccharomyces cerevisiae. This makes it an excellent model system to understand the genotype-phenotype relationship. In this paper, we presented the reconstruction and application of the yeast pan-genome in resolving genotype-phenotype relationship by a machine learning-assisted approach.

Resistance Mechanisms of Saccharomyces cerevisiae to Commercial Formulations of Glyphosate Involve DNA Damage Repair, the Cell Cycle, and the Cell Wall Structure

The use of glyphosate-based herbicides is widespread and despite their extensive use, their effects are yet to be deciphered completely. The additives in commercial formulations of glyphosate, though labeled inert when used individually, have adverse effects when used in combination with other additives along with the active ingredient. As a species, Saccharomyces cerevisiae has a wide range of resistance to glyphosate-based herbicides. To investigate the underlying genetic differences between sensitive and resistant strains, global changes in gene expression were measured, when yeast were exposed to a glyphosate-based herbicide (GBH). Expression of genes involved in numerous pathways crucial to the cell's functioning, such as DNA replication, MAPK signaling, meiosis, and cell wall synthesis changed. Because so many diverse pathways were affected, these strains were then subjected to in-lab-evolutions (ILE) to select mutations that confer increased resistance. Common fragile sites were found to play a role in adaptation to resistance to long-term exposure of GBHs. Copy number increased in approximately 100 genes associated with cell wall proteins, mitochondria, and sterol transport. Taking ILE and transcriptomic data into account it is evident that GBHs affect multiple biological processes in the cell. One such component is the cell wall structure which acts as a protective barrier in alleviating the stress caused by exposure to inert additives in GBHs. Sed1, a GPI-cell wall protein, plays an important role in tolerance of a GBH. Hence, a detailed study of the changes occurring at the genome and transcriptome levels is essential to better understand the effects of an environmental stressor such as a GBH, on the cell as a whole.

MCHM Acts as a Hydrotrope, Altering the Balance of Metals in Yeast

While drugs and other industrial chemicals are routinely studied to assess risks, many widely used chemicals have not been thoroughly evaluated. One such chemical, 4-methylcyclohexane methanol (MCHM), is an industrial coal-cleaning chemical that contaminated the drinking water supply in Charleston, WV, USA in 2014. While a wide range of ailments was reported following the spill, little is known about the molecular effects of MCHM exposure. We used the yeast model to explore the impacts of MCHM on cellular function. Exposure to MCHM dramatically altered the yeast transcriptome and the balance of metals in yeast. Underlying genetic variation in the response to MCHM, transcriptomics and, mutant analysis uncovered the role of the metal transporters, Arn2 and Yke4, to MCHM response. Expression of Arn2, which is involved in iron uptake, was lower in MCHM-tolerant yeast and loss of Arn2 further increased MCHM tolerance. Genetic variation within Yke4, an ER zinc transporter, also mediated response to MCHM, and loss of Yke4 decreased MCHM tolerance. The addition of zinc to MCHM-sensitive yeast rescued growth inhibition. In vitro assays demonstrated that MCHM acted as a hydrotrope and prevented protein interactions, while zinc induced the aggregation of proteins. We hypothesized that MCHM altered the structures of extracellular domains of proteins, and the addition of zinc stabilized the structure to maintain metal homeostasis in yeast exposed to MCHM.

Challenges and opportunities for strain verification by whole-genome sequencing

Laboratory strains, cell lines, and other genetic materials change hands frequently in the life sciences. Despite evidence that such materials are subject to mix-ups, contamination, and accumulation of secondary mutations, verification of strains and samples is not an established part of many experimental workflows. With the plummeting cost of next generation technologies, it is conceivable that whole genome sequencing (WGS) could be applied to routine strain and sample verification in the future. To demonstrate the need for strain validation by WGS, we sequenced haploid yeast segregants derived from a popular commercial mutant collection and identified several unexpected mutations. We determined that available bioinformatics tools may be ill-suited for verification and highlight the importance of finishing reference genomes for commonly used laboratory strains.

The Yeast Atlas of Appalachia: Species and Phenotypic Diversity of Herbicide Resistance in Wild Yeast

Glyphosate and copper-based herbicides/fungicides affect non-target organisms, and these incidental exposures can impact microbial populations. In this study, glyphosate resistance was found in the historical collection of S. cerevisiae, which was collected over the last century, but only in yeast isolated after the introduction of glyphosate. Although herbicide application was not recorded, the highest glyphosate-resistant S. cerevisiae were isolated from agricultural sites. In an effort to assess glyphosate resistance and impact on non-target microorganisms, different yeast species were harvested from 15 areas with known herbicidal histories, including an organic farm, conventional farm, remediated coal mine, suburban locations, state park, and a national forest. Yeast representing 23 genera were isolated from 237 samples of plant, soil, spontaneous fermentation, nut, flower, fruit, feces, and tree material samples. Saccharomyces, Candida, Metschnikowia, Kluyveromyces, Hanseniaspora, and Pichia were other genera commonly found across our sampled environments. Managed areas had less species diversity, and at the brewery only Saccharomyces and Pichia were isolated. A conventional farm growing RoundUp Ready™ corn had the lowest phylogenetic diversity and the highest glyphosate resistance. The mine was sprayed with multiple herbicides including a commercial formulation of glyphosate; however, the S. cerevisiae did not have elevated glyphosate resistance. In contrast to the conventional farm, the mine was exposed to glyphosate only one year prior to sample isolation. Glyphosate resistance is an example of the anthropogenic selection of nontarget organisms.

The Polymorphic PolyQ Tail Protein of the Mediator Complex, Med15, Regulates the Variable Response to Diverse Stresses

The Mediator is composed of multiple subunits conserved from yeast to humans and plays a central role in transcription. The tail components are not required for basal transcription but are required for responses to different stresses. While some stresses are familiar, such as heat, desiccation, and starvation, others are exotic, yet yeast can elicit a successful stress response. 4-Methylcyclohexane methanol (MCHM) is a hydrotrope that induces growth arrest in yeast. We found that a naturally occurring variation in the Med15 allele, a component of the Mediator tail, altered the stress response to many chemicals in addition to MCHM. Med15 contains two polyglutamine repeats (polyQ) of variable lengths that change the gene expression of diverse pathways. The Med15 protein existed in multiple isoforms and its stability was dependent on Ydj1, a protein chaperone. The protein level of Med15 with longer polyQ tracts was lower and turned over faster than the allele with shorter polyQ repeats. MCHM sensitivity via variation of Med15 was regulated by Snf1 in a Myc-tag-dependent manner. Tagging Med15 with Myc altered its function in response to stress. Genetic variation in transcriptional regulators magnified genetic differences in response to environmental changes. These polymorphic control genes were master variators.

Evolution of Ty1 copy number control in yeast by horizontal transfer and recombination

Transposable elements constitute a large fraction of most eukaryotic genomes. Insertion of mobile DNA sequences typically has deleterious effects on host fitness, and thus diverse mechanisms have evolved to control mobile element proliferation. Mobility of the Ty1 retrotransposon in Saccharomyces yeasts is regulated by copy number control (CNC) mediated by a self-encoded restriction factor derived from the Ty1 gag capsid gene that inhibits virus-like particle function. Here, we survey a panel of wild and human-associated strains of S. cerevisiae and S. paradoxus to investigate how genomic Ty1 content influences variation in Ty1 mobility. We observe high levels of mobility for a tester element with a gag sequence from the canonical Ty1 subfamily in permissive strains that either lack full-length Ty1 elements or only contain full-length copies of the Ty1' subfamily that have a divergent gag sequence. In contrast, low levels of canonical Ty1 mobility are observed in restrictive strains carrying full-length Ty1 elements containing a canonical gag sequence. Phylogenomic analysis of full-length Ty1 elements revealed that Ty1' is the ancestral subfamily present in wild strains of S. cerevisiae, and that canonical Ty1 in S. cerevisiae is a derived subfamily that acquired gag from S. paradoxus by horizontal transfer and recombination. Our results provide evidence that variation in the ability of S. cerevisiae and S. paradoxus strains to repress canonical Ty1 transposition via CNC is regulated by the genomic content of different Ty1 subfamilies, and that self-encoded forms of transposon control can spread across species boundaries by horizontal transfer.

Saccharomyces arboricola and Its Hybrids’ Propensity for Sake Production: Interspecific Hybrids Reveal Increased Fermentation Abilities and a Mosaic Metabolic Profile

The use of interspecific hybrids during the industrial fermentation process has been well established, positioning the frontier of advancement in brewing to capitalize on the potential of Saccharomyces hybridization. Interspecific yeast hybrids used in modern monoculture inoculations benefit from a wide range of volatile metabolites that broaden the organoleptic complexity. This is the first report of sake brewing by Saccharomyces arboricola and its hybrids. S. arboricola x S. cerevisiae direct-mating generated cryotolerant interspecific hybrids which increased yields of ethanol and ethyl hexanoate compared to parental strains, important flavor attributes of fine Japanese ginjo sake rice wine. Hierarchical clustering heatmapping with principal component analysis for metabolic profiling was used in finding low levels of endogenous amino/organic acids clustered S. arboricola apart from the S. cerevisiae industrial strains. In sake fermentations, hybrid strains showed a mosaic profile of parental strains, while metabolic analysis suggested S. arboricola had a lower amino acid net uptake than S. cerevisiae. Additionally, this research found an increase in ethanolic fermentation from pyruvate and increased sulfur metabolism. Together, these results suggest S. arboricola is poised for in-depth metabolomic exploration in sake fermentation.

Pan-genome analyses of model fungal species

The concept of the species ‘pan-genome’, the union of ‘core’ conserved genes and all ‘accessory’ non-conserved genes across all strains of a species, was first proposed in prokaryotes to account for intraspecific variability. Species pan-genomes have been extensively studied in prokaryotes, but evidence of species pan-genomes has also been demonstrated in eukaryotes such as plants and fungi. Using a previously published methodology based on sequence homology and conserved microsynteny, in addition to bespoke pipelines, we have investigated the pan-genomes of four model fungal species: Saccharomyces cerevisiae, Candida albicans, Cryptococcus neoformans var. grubii and Aspergillus fumigatus. Between 80 and 90 % of gene models per strain in each of these species are core genes that are highly conserved across all strains of that species, many of which are involved in housekeeping and conserved survival processes. In many of these species, the remaining ‘accessory’ gene models are clustered within subterminal regions and may be involved in pathogenesis and antimicrobial resistance. Analysis of the ancestry of species core and accessory genomes suggests that fungal pan-genomes evolve by strain-level innovations such as gene duplication as opposed to wide-scale horizontal gene transfer. Our findings lend further supporting evidence to the existence of species pan-genomes in eukaryote taxa.

Diversity analysis of Saccharomyces cerevisiae isolated from natural sources by multilocus sequence typing (MLST)

We used multilocus sequence typing (MLST) to analyze the diversity of natural isolates of Saccharomyces cerevisiae, the most important microorganism in alcoholic fermentation. Six loci, ADP1, RPN2, GLN4, ACC1, MET4, and NUP116, in S. cerevisiae genome were selected as MLST markers. To investigate genetic diversity within S. cerevisiae, 42 S. cerevisiae isolated from natural sources in Korea as well as six S. cerevisiae obtained from Genbank and four industrial S. cerevisiae were examined using MLST. Twenty-six polymorphic sites were found in the six loci. Among them, ACC1 had the most genetic variation with eight polymorphic sites. MLST differentiated the 52 strains into three clades. Alcohol fermentation results revealed that S. cerevisiae in Clade III produced less alcohol than those in Clades I and II. These results suggested that MLST is a powerful tool to differentiate S. cerevisiae and can potentially be used to select S. cerevisiae suitable for industrial use.

Whole-Genome Analysis of Three Yeast Strains Used for Production of Sherry-Like Wines Revealed Genetic Traits Specific to Flor Yeasts

Flor yeast strains represent a specialized group of Saccharomyces cerevisiae yeasts used for biological wine aging. We have sequenced the genomes of three flor strains originated from different geographic regions and used for production of sherry-like wines in Russia. According to the obtained phylogeny of 118 yeast strains, flor strains form very tight cluster adjacent to the main wine clade. SNP analysis versus available genomes of wine and flor strains revealed 2,270 genetic variants in 1,337 loci specific to flor strains. Gene ontology analysis in combination with gene content evaluation revealed a complex landscape of possibly adaptive genetic changes in flor yeast, related to genes associated with cell morphology, mitotic cell cycle, ion homeostasis, DNA repair, carbohydrate metabolism, lipid metabolism, and cell wall biogenesis. Pangenomic analysis discovered the presence of several well-known "non-reference" loci of potential industrial importance. Events of gene loss included deletions of asparaginase genes, maltose utilization locus, and FRE-FIT locus involved in iron transport. The latter in combination with a flor-yeast-specific mutation in the Aft1 transcription factor gene is likely to be responsible for the discovered phenotype of increased iron sensitivity and improved iron uptake of analyzed strains. Expansion of the coding region of the FLO11 flocullin gene and alteration of the balance between members of the FLO gene family are likely to positively affect the well-known propensity of flor strains for velum formation. Our study provides new insights in the nature of genetic variation in flor yeast strains and demonstrates that different adaptive properties of flor yeast strains could have evolved through different mechanisms of genetic variation.

KinFin: Software for Taxon-Aware Analysis of Clustered Protein Sequences

The field of comparative genomics is concerned with the study of similarities and differences between the information encoded in the genomes of organisms. A common approach is to define gene families by clustering protein sequences based on sequence similarity, and analyze protein cluster presence and absence in different species groups as a guide to biology. Due to the high dimensionality of these data, downstream analysis of protein clusters inferred from large numbers of species, or species with many genes, is nontrivial, and few solutions exist for transparent, reproducible, and customizable analyses. We present KinFin, a streamlined software solution capable of integrating data from common file formats and delivering aggregative annotation of protein clusters. KinFin delivers analyses based on systematic taxonomy of the species analyzed, or on user-defined, groupings of taxa, for example, sets based on attributes such as life history traits, organismal phenotypes, or competing phylogenetic hypotheses. Results are reported through graphical and detailed text output files. We illustrate the utility of the KinFin pipeline by addressing questions regarding the biology of filarial nematodes, which include parasites of veterinary and medical importance. We resolve the phylogenetic relationships between the species and explore functional annotation of proteins in clusters in key lineages and between custom taxon sets, identifying gene families of interest. KinFin can easily be integrated into existing comparative genomic workflows, and promotes transparent and reproducible analysis of clustered protein data.

Genome Diversity and Evolution in the Budding Yeasts (Saccharomycotina)

Considerable progress in our understanding of yeast genomes and their evolution has been made over the last decade with the sequencing, analysis, and comparisons of numerous species, strains, or isolates of diverse origins. The role played by yeasts in natural environments as well as in artificial manufactures, combined with the importance of some species as model experimental systems sustained this effort. At the same time, their enormous evolutionary diversity (there are yeast species in every subphylum of Dikarya) sparked curiosity but necessitated further efforts to obtain appropriate reference genomes. Today, yeast genomes have been very informative about basic mechanisms of evolution, speciation, hybridization, domestication, as well as about the molecular machineries underlying them. They are also irreplaceable to investigate in detail the complex relationship between genotypes and phenotypes with both theoretical and practical implications. This review examines these questions at two distinct levels offered by the broad evolutionary range of yeasts: inside the best-studied Saccharomyces species complex, and across the entire and diversified subphylum of Saccharomycotina. While obviously revealing evolutionary histories at different scales, data converge to a remarkably coherent picture in which one can estimate the relative importance of intrinsic genome dynamics, including gene birth and loss, vs. horizontal genetic accidents in the making of populations. The facility with which novel yeast genomes can now be studied, combined with the already numerous available reference genomes, offer privileged perspectives to further examine these fundamental biological questions using yeasts both as eukaryotic models and as fungi of practical importance.

Non-invasive real time monitoring of yeast volatilome by PTR-ToF-MS

INTRODUCTION

Producing a wide range of volatile secondary metabolites Saccharomyces cerevisiae influences wine, beer, and bread sensory quality and hence selection of strains based on their volatilome becomes pivotal. A rapid on-line method for volatilome assessing of strains growing on standard solid media is still missing.

OBJECTIVES

Methodologically, the aim of this study was to demonstrate the automatic, real-time, direct, and non-invasive monitoring of yeast volatilome in order to rapidly produce a robust large data set encompassing measurements relative to many strains, replicates and time points. The fundamental scope was to differentiate volatilomes of genetically similar strains of oenological relevance during the whole growing process.

METHOD

Six different S. cerevisiae strains (four meiotic segregants of a natural strain and two laboratory strains) inoculated onto a solid medium have been monitored on-line by Proton Transfer Reaction-Time-of-Flight-Mass Spectrometry for 11 days every 4 h (3540 time points). FastGC PTR-ToF-MS was performed during the stationary phase on the 5th day.

RESULTS

More than 300 peaks have been extracted from the average spectra associated to each time point, 70 have been tentatively identified. Univariate and multivariate analyses have been performed on the data matrix (3640 measurements × 70 peaks) highlighting the volatilome evolution and strain-specific features. Laboratory strains with opposite mating type, and meiotic segregants of the same natural strain showed significantly different profiles.

CONCLUSIONS

The described set-up allows the on-line high-throughput screening of yeast volatilome of S. cerevisiae strains and the identification of strain specific features and new metabolic pathways, discriminating also genetically similar strains, thus revealing a novel method for strain phenotyping, identification, and quality control.

Moss Phylogeny Reconstruction Using Nucleotide Pangenome of Complete Mitogenome Sequences

Stability of composition and sequence of genes was shown earlier in 13 mitochondrial genomes of mosses (Rensing, S. A., et al. (2008) Science, 319, 64-69). It is of interest to study the evolution of mitochondrial genomes not only at the gene level, but also on the level of nucleotide sequences. To do this, we have constructed a "nucleotide pangenome" for mitochondrial genomes of 24 moss species. The nucleotide pangenome is a set of aligned nucleotide sequences of orthologous genome fragments covering the totality of all genomes. The nucleotide pangenome was constructed using specially developed new software, NPG-explorer (NPGe). The stable part of the mitochondrial genome (232 stable blocks) is shown to be, on average, 45% of its length. In the joint alignment of stable blocks, 82% of positions are conserved. The phylogenetic tree constructed with the NPGe program is in good correlation with other phylogenetic reconstructions. With the NPGe program, 30 blocks have been identified with repeats no shorter than 50 bp. The maximal length of a block with repeats is 140 bp. Duplications in the mitochondrial genomes of mosses are rare. On average, the genome contains about 500 bp in large duplications. The total length of insertions and deletions was determined in each genome. The losses and gains of DNA regions are rather active in mitochondrial genomes of mosses, and such rearrangements presumably can be used as additional markers in the reconstruction of phylogeny.

Integration of new alternative reference strain genome sequences into the Saccharomyces genome database

The Saccharomyces Genome Database (SGD; http://www.yeastgenome.org/) is the authoritative community resource for the Saccharomyces cerevisiae reference genome sequence and its annotation. To provide a wider scope of genetic and phenotypic variation in yeast, the genome sequences and their corresponding annotations from 11 alternative S. cerevisiae reference strains have been integrated into SGD. Genomic and protein sequence information for genes from these strains are now available on the Sequence and Protein tab of the corresponding Locus Summary pages. We illustrate how these genome sequences can be utilized to aid our understanding of strain-specific functional and phenotypic differences.

Database URL:www.yeastgenome.org

Genome Sequence and Analysis of a Stress-Tolerant, Wild-Derived Strain of Saccharomyces cerevisiae Used in Biofuels Research

The genome sequences of more than 100 strains of the yeast Saccharomyces cerevisiae have been published. Unfortunately, most of these genome assemblies contain dozens to hundreds of gaps at repetitive sequences, including transposable elements, tRNAs, and subtelomeric regions, which is where novel genes generally reside. Relatively few strains have been chosen for genome sequencing based on their biofuel production potential, leaving an additional knowledge gap. Here, we describe the nearly complete genome sequence of GLBRCY22-3 (Y22-3), a strain of S. cerevisiae derived from the stress-tolerant wild strain NRRL YB-210 and subsequently engineered for xylose metabolism. After benchmarking several genome assembly approaches, we developed a pipeline to integrate Pacific Biosciences (PacBio) and Illumina sequencing data and achieved one of the highest quality genome assemblies for any S. cerevisiae strain. Specifically, the contig N50 is 693 kbp, and the sequences of most chromosomes, the mitochondrial genome, and the 2-micron plasmid are complete. Our annotation predicts 92 genes that are not present in the reference genome of the laboratory strain S288c, over 70% of which were expressed. We predicted functions for 43 of these genes, 28 of which were previously uncharacterized and unnamed. Remarkably, many of these genes are predicted to be involved in stress tolerance and carbon metabolism and are shared with a Brazilian bioethanol production strain, even though the strains differ dramatically at most genetic loci. The Y22-3 genome sequence provides an exceptionally high-quality resource for basic and applied research in bioenergy and genetics.

CRISPR-directed mitotic recombination enables genetic mapping without crosses

Linkage and association studies have mapped thousands of genomic regions that contribute to phenotypic variation, but narrowing these regions to the underlying causal genes and variants has proven much more challenging. Resolution of genetic mapping is limited by the recombination rate. We developed a method that uses CRISPR (clustered, regularly interspaced, short palindromic repeats) to build mapping panels with targeted recombination events. We tested the method by generating a panel with recombination events spaced along a yeast chromosome arm, mapping trait variation, and then targeting a high density of recombination events to the region of interest. Using this approach, we fine-mapped manganese sensitivity to a single polymorphism in the transporter Pmr1. Targeting recombination events to regions of interest allows us to rapidly and systematically identify causal variants underlying trait differences.

Genome Sequencing and Comparative Analysis of Saccharomyces cerevisiae Strains of the Peterhof Genetic Collection

The Peterhof genetic collection of Saccharomyces cerevisiae strains (PGC) is a large laboratory stock that has accumulated several thousands of strains for over than half a century. It originated independently of other common laboratory stocks from a distillery lineage (race XII). Several PGC strains have been extensively used in certain fields of yeast research but their genomes have not been thoroughly explored yet. Here we employed whole genome sequencing to characterize five selected PGC strains including one of the closest to the progenitor, 15V-P4, and several strains that have been used to study translation termination and prions in yeast (25-25-2V-P3982, 1B-D1606, 74-D694, and 6P-33G-D373). The genetic distance between the PGC progenitor and S288C is comparable to that between two geographically isolated populations. The PGC seems to be closer to two bakery strains than to S288C-related laboratory stocks or European wine strains. In genomes of the PGC strains, we found several loci which are absent from the S288C genome; 15V-P4 harbors a rare combination of the gene cluster characteristic for wine strains and the RTM1 cluster. We closely examined known and previously uncharacterized gene variants of particular strains and were able to establish the molecular basis for known phenotypes including phenylalanine auxotrophy, clumping behavior and galactose utilization. Finally, we made sequencing data and results of the analysis available for the yeast community. Our data widen the knowledge about genetic variation between Saccharomyces cerevisiae strains and can form the basis for planning future work in PGC-related strains and with PGC-derived alleles.

From one to many: expanding the Saccharomyces cerevisiae reference genome panel

In recent years, thousands of Saccharomyces cerevisiae genomes have been sequenced to varying degrees of completion. The Saccharomyces Genome Database (SGD) has long been the keeper of the original eukaryotic reference genome sequence, which was derived primarily from S. cerevisiae strain S288C. Because new technologies are pushing S. cerevisiae annotation past the limits of any system based exclusively on a single reference sequence, SGD is actively working to expand the original S. cerevisiae systematic reference sequence from a single genome to a multi-genome reference panel. We first commissioned the sequencing of additional genomes and their automated analysis using the AGAPE pipeline. Here we describe our curation strategy to produce manually reviewed high-quality genome annotations in order to elevate 11 of these additional genomes to Reference status.

Database URL: http://www.yeastgenome.org/

The Saccharomyces Genome Database Variant Viewer

The Saccharomyces Genome Database (SGD; http://www.yeastgenome.org) is the authoritative community resource for the Saccharomyces cerevisiae reference genome sequence and its annotation. In recent years, we have moved toward increased representation of sequence variation and allelic differences within S. cerevisiae. The publication of numerous additional genomes has motivated the creation of new tools for their annotation and analysis. Here we present the Variant Viewer: a dynamic open-source web application for the visualization of genomic and proteomic differences. Multiple sequence alignments have been constructed across high quality genome sequences from 11 different S. cerevisiae strains and stored in the SGD. The alignments and summaries are encoded in JSON and used to create a two-tiered dynamic view of the budding yeast pan-genome, available at http://www.yeastgenome.org/variant-viewer.

Review of current methods for characterizing virulence and pathogenicity potential of industrial Saccharomyces cerevisiae strains towards humans

Most industrial Saccharomyces cerevisiae strains used in food or biotechnology processes are benign. However, reports of S. cerevisiae infections have emerged and novel strains continue to be developed. In order to develop recommendations for the human health risk assessment of S. cerevisiae strains, we conducted a literature review of current methods used to characterize their pathogenic potential and evaluated their relevance towards risk assessment. These studies revealed that expression of virulence traits in S. cerevisiae is complex and depends on many factors. Given the opportunistic nature of this organism, an approach using multiple lines of evidence is likely necessary for the reasonable prediction of the pathogenic potential of a particular strain. Risk assessment of S. cerevisiae strains would benefit from more research towards the comparison of virulent and non-virulent strains in order to better understand those genotypic and phenotypic traits most likely to be associated with pathogenicity.

Clade- and species-specific features of genome evolution in the Saccharomycetaceae

Many aspects of the genomes of yeast species in the family Saccharomycetaceae have been well conserved during evolution. They have similar genome sizes, genome contents, and extensive collinearity of gene order along chromosomes. Gene functions can often be inferred reliably by using information from Saccharomyces cerevisiae. Beyond this conservative picture however, there are many instances where a species or a clade diverges substantially from the S. cerevisiae paradigm—for example, by the amplification of a gene family, or by the absence of a biochemical pathway or a protein complex. Here, we review clade-specific features, focusing on genomes sequenced in our laboratory from the post-WGD genera Naumovozyma, Kazachstania and Tetrapisispora, and from the non-WGD species Torulaspora delbrueckii. Examples include the loss of the pathway for histidine synthesis in the cockroach-associated species Tetrapisispora blattae; the presence of a large telomeric GAL gene cluster in To. delbrueckii; losses of the dynein and dynactin complexes in several independent yeast lineages; fragmentation of the MAT locus and loss of the HO gene in Kazachstania africana; and the patchy phylogenetic distribution of RNAi pathway components.

The authors review species-specific evolutionary attributes of yeast genomes.