Natural proteome diversity links aneuploidy tolerance to protein turnover

Accessing the natural genetic diversity of species unveils hidden genetic traits, clarifies gene functions and allows the generalizability of laboratory findings to be assessed. One notable discovery made in natural isolates of Saccharomyces cerevisiae is that aneuploidy-an imbalance in chromosome copy numbers-is frequent1,2 (around 20%), which seems to contradict the substantial fitness costs and transient nature of aneuploidy when it is engineered in the laboratory3-5. Here we generate a proteomic resource and merge it with genomic1 and transcriptomic6 data for 796 euploid and aneuploid natural isolates. We find that natural and lab-generated aneuploids differ specifically at the proteome. In lab-generated aneuploids, some proteins-especially subunits of protein complexes-show reduced expression, but the overall protein levels correspond to the aneuploid gene dosage. By contrast, in natural isolates, more than 70% of proteins encoded on aneuploid chromosomes are dosage compensated, and average protein levels are shifted towards the euploid state chromosome-wide. At the molecular level, we detect an induction of structural components of the proteasome, increased levels of ubiquitination, and reveal an interdependency of protein turnover rates and attenuation. Our study thus highlights the role of protein turnover in mediating aneuploidy tolerance, and shows the utility of exploiting the natural diversity of species to attain generalizable molecular insights into complex biological processes.

Ancient and recent origins of shared polymorphisms in yeast

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Engineering heterothallic strains in fission yeast

In poor nitrogen conditions, fission yeast cells mate, undergo meiosis and form spores that are resistant to deleterious environments. Natural isolates of Schizosaccharomyces pombe are homothallic. This allows them to naturally switch between the two h- and h+ mating types with a high frequency, thereby ensuring the presence of both mating partners in a population of cells. However, alteration of the mating type locus can abolish mating type switching or reduce it to a very low frequency. Such heterothallic strains have been isolated and are common in research laboratories due to the simplicity of their use for Mendelian genetics. In addition to the standard laboratory strains, a large collection of natural S. pombe isolates is now available, representing a powerful resource for investigating the genetic diversity and biology of fission yeast. However, most of these strains are homothallic, and only tedious or mutagenic strategies have been described to obtain heterothallic cells from a homothallic parent. Here, we describe a simple approach to generate heterothallic strains. It takes advantage of an alteration of the mating type locus that was previously identified in a mating type switching-deficient strain and the CRISPR-Cas9 editing tool, allowing for a one-step engineering of heterothallic cells with high efficiency.

Genome instability footprint under rapamycin and hydroxyurea treatments

The mutational processes dictating the accumulation of mutations in genomes are shaped by genetic background, environment and their interactions. Accurate quantification of mutation rates and spectra under drugs has important implications in disease treatment. Here, we used whole-genome sequencing and time-resolved growth phenotyping of yeast mutation accumulation lines to give a detailed view of the mutagenic effects of rapamycin and hydroxyurea on the genome and cell growth. Mutation rates depended on the genetic backgrounds but were only marginally affected by rapamycin. As a remarkable exception, rapamycin treatment was associated with frequent chromosome XII amplifications, which compensated for rapamycin induced rDNA repeat contraction on this chromosome and served to maintain rDNA content homeostasis and fitness. In hydroxyurea, a wide range of mutation rates were elevated regardless of the genetic backgrounds, with a particularly high occurrence of aneuploidy that associated with dramatic fitness loss. Hydroxyurea also induced a high T-to-G and low C-to-A transversion rate that reversed the common G/C-to-A/T bias in yeast and gave rise to a broad range of structural variants, including mtDNA deletions. The hydroxyurea mutation footprint was consistent with the activation of error-prone DNA polymerase activities and non-homologues end joining repair pathways. Taken together, our study provides an in-depth view of mutation rates and signatures in rapamycin and hydroxyurea and their impact on cell fitness, which brings insights for assessing their chronic effects on genome integrity.

phyBWT2: phylogeny reconstruction via eBWT positional clustering

BACKGROUND

Molecular phylogenetics studies the evolutionary relationships among the individuals of a population through their biological sequences. It may provide insights about the origin and the evolution of viral diseases, or highlight complex evolutionary trajectories. A key task is inferring phylogenetic trees from any type of sequencing data, including raw short reads. Yet, several tools require pre-processed input data e.g. from complex computational pipelines based on de novo assembly or from mappings against a reference genome. As sequencing technologies keep becoming cheaper, this puts increasing pressure on designing methods that perform analysis directly on their outputs. From this viewpoint, there is a growing interest in alignment-, assembly-, and reference-free methods that could work on several data including raw reads data.

RESULTS

We present phyBWT2, a newly improved version of phyBWT (Guerrini et al. in 22nd International Workshop on Algorithms in Bioinformatics (WABI) 242:23-12319, 2022). Both of them directly reconstruct phylogenetic trees bypassing both the alignment against a reference genome and de novo assembly. They exploit the combinatorial properties of the extended Burrows-Wheeler Transform (eBWT) and the corresponding eBWT positional clustering framework to detect relevant blocks of the longest shared substrings of varying length (unlike the k-mer-based approaches that need to fix the length k a priori). As a result, they provide novel alignment-, assembly-, and reference-free methods that build partition trees without relying on the pairwise comparison of sequences, thus avoiding to use a distance matrix to infer phylogeny. In addition, phyBWT2 outperforms phyBWT in terms of running time, as the former reconstructs phylogenetic trees step-by-step by considering multiple partitions, instead of just one partition at a time, as previously done by the latter.

CONCLUSIONS

Based on the results of the experiments on sequencing data, we conclude that our method can produce trees of quality comparable to the benchmark phylogeny by handling datasets of different types (short reads, contigs, or entire genomes). Overall, the experiments confirm the effectiveness of phyBWT2 that improves the performance of its previous version phyBWT, while preserving the accuracy of the results.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Single nucleotide polymorphisms associated with wine fermentation and adaptation to nitrogen limitation in wild and domesticated yeast strains

For more than 20 years, Saccharomyces cerevisiae has served as a model organism for genetic studies and molecular biology, as well as a platform for biotechnology (e.g., wine production). One of the important ecological niches of this yeast that has been extensively studied is wine fermentation, a complex microbiological process in which S. cerevisiae faces various stresses such as limited availability of nitrogen. Nitrogen deficiencies in grape juice impair fermentation rate and yeast biomass production, leading to sluggish or stuck fermentations, resulting in considerable economic losses for the wine industry. In the present work, we took advantage of the "1002 Yeast Genomes Project" population, the most complete catalogue of the genetic variation in the species and a powerful resource for genotype-phenotype correlations, to study the adaptation to nitrogen limitation in wild and domesticated yeast strains in the context of wine fermentation. We found that wild and domesticated yeast strains have different adaptations to nitrogen limitation, corroborating their different evolutionary trajectories. Using a combination of state-of-the-art bioinformatic (GWAS) and molecular biology (CRISPR-Cas9) methodologies, we validated that PNP1, RRT5 and PDR12 are implicated in wine fermentation, where RRT5 and PDR12 are also involved in yeast adaptation to nitrogen limitation. In addition, we validated SNPs in these genes leading to differences in fermentative capacities and adaptation to nitrogen limitation. Altogether, the mapped genetic variants have potential applications for the genetic improvement of industrial yeast strains.

Building pangenome graphs

Pangenome graphs can represent all variation between multiple genomes, but existing methods for constructing them are biased due to reference-guided approaches. In response, we have developed PanGenome Graph Builder (PGGB), a reference-free pipeline for constructing unbi-ased pangenome graphs. PGGB uses all-to-all whole-genome alignments and learned graph embeddings to build and iteratively refine a model in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Spontaneous Mutation Rates and Spectra of Respiratory-Deficient Yeast

The yeast petite mutant was first discovered in the yeast Saccharomyces cerevisiae, which shows growth stress due to defects in genes encoding the respiratory chain. In a previous study, we described that deletion of the nuclear-encoded gene MRPL25 leads to mitochondrial genome (mtDNA) loss and the petite phenotype, which can be rescued by acquiring ATP3 mutations. The mrpl25Δ strain showed an elevated SNV (single nucleotide variant) rate, suggesting genome instability occurred during the crisis of mtDNA loss. However, the genome-wide mutation landscape and mutational signatures of mitochondrial dysfunction are unknown. In this study we profiled the mutation spectra in yeast strains with the genotype combination of MRPL25 and ATP3 in their wildtype and mutated status, along with the wildtype and cytoplasmic petite rho0 strains as controls. In addition to the previously described elevated SNV rate, we found the INDEL (insertion/deletion) rate also increased in the mrpl25Δ strain, reinforcing the occurrence of genome instability. Notably, although both are petites, the mrpl25Δ and rho0 strains exhibited different INDEL rates and transition/transversion ratios, suggesting differences in the mutational signatures underlying these two types of petites. Interestingly, the petite-related mutagenesis effect disappeared when ATP3 suppressor mutations were acquired, suggesting a cost-effective mechanism for restoring both fitness and genome stability. Taken together, we present an unbiased genome-wide characterization of the mutation rates and spectra of yeast strains with respiratory deficiency, which provides valuable insights into the impact of respiratory deficiency on genome instability.

Telomeres are shorter in wild Saccharomyces cerevisiae isolates than in domesticated ones

Telomeres are ribonucleoproteins that cap chromosome-ends and their DNA length is controlled by counteracting elongation and shortening processes. The budding yeast Saccharomyces cerevisiae has been a leading model to study telomere DNA length control and dynamics. Its telomeric DNA is maintained at a length that slightly varies between laboratory strains, but little is known about its variation at the species level. The recent publication of the genomes of over 1,000 S. cerevisiae strains enabled us to explore telomere DNA length variation at an unprecedented scale. Here, we developed a bioinformatic pipeline (YeaISTY) to estimate telomere DNA length from whole-genome sequences and applied it to the sequenced S. cerevisiae collection. Our results revealed broad natural telomere DNA length variation among the isolates. Notably, telomere DNA length is shorter in those derived from wild rather than domesticated environments. Moreover, telomere DNA length variation is associated with mitochondrial metabolism, and this association is driven by wild strains. Overall, these findings reveal broad variation in budding yeast's telomere DNA length regulation, which might be shaped by its different ecological life-styles.

Genetically controlled mtDNA deletions prevent ROS damage by arresting oxidative phosphorylation

Deletion of mitochondrial DNA in eukaryotes is currently attributed to rare accidental events associated with mitochondrial replication or repair of double-strand breaks. We report the discovery that yeast cells arrest harmful intramitochondrial superoxide production by shutting down respiration through genetically controlled deletion of mitochondrial oxidative phosphorylation genes. We show that this process critically involves the antioxidant enzyme superoxide dismutase 2 and two-way mitochondrial-nuclear communication through Rtg2 and Rtg3. While mitochondrial DNA homeostasis is rapidly restored after cessation of a short-term superoxide stress, long-term stress causes maladaptive persistence of the deletion process, leading to complete annihilation of the cellular pool of intact mitochondrial genomes and irrevocable loss of respiratory ability. This shows that oxidative stress-induced mitochondrial impairment may be under strict regulatory control. If the results extend to human cells, the results may prove to be of etiological as well as therapeutic importance with regard to age-related mitochondrial impairment and disease.

Resistance to Arsenite and Arsenate in Saccharomyces cerevisiae Arises through the Subtelomeric Expansion of a Cluster of Yeast Genes

Arsenic is one of the most prevalent toxic elements in the environment, and its toxicity affects every organism. Arsenic resistance has mainly been observed in microorganisms, and, in bacteria, it has been associated with the presence of the Ars operon. In Saccharomyces cerevisiae, three genes confer arsenic resistance: ARR1, ARR2, and ARR3. Unlike bacteria, in which the presence of the Ars genes confers per se resistance to arsenic, most of the S. cerevisiae isolates present the three ARR genes, regardless of whether the strain is resistant or sensitive to arsenic. To assess the genetic features that make natural S. cerevisiae strains resistant to arsenic, we used a combination of comparative genomic hybridization, whole-genome sequencing, and transcriptomics profiling with microarray analyses. We observed that both the presence and the genomic location of multiple copies of the whole cluster of ARR genes were central to the escape from subtelomeric silencing and the acquisition of resistance to arsenic. As a result of the repositioning, the ARR genes were expressed even in the absence of arsenic. In addition to their relevance in improving our understanding of the mechanism of arsenic resistance in yeast, these results provide evidence for a new cluster of functionally related genes that are independently duplicated and translocated.

Extensive simulations assess the performance of genome-wide association mapping in various Saccharomyces cerevisiae subpopulations

With the advent of high throughput sequencing technologies, genome-wide association studies (GWAS) have become a powerful paradigm for dissecting the genetic origins of the observed phenotypic variation. We recently completely sequenced the genome of 1011 Saccharomyces cerevisiae isolates, laying a strong foundation for GWAS. To assess the feasibility and the limits of this approach, we performed extensive simulations using five selected subpopulations as well as the total set of 1011 genomes. We measured the ability to detect the causal genetic variants involved in Mendelian and more complex traits using a linear mixed model approach. The results showed that population structure is well accounted for and is not the main problem when the sample size is high enough. While the genetic determinant of a Mendelian trait is easily mapped in all studied subpopulations, discrepancies are seen between datasets when performing GWAS on a complex trait in terms of detection, false positive and false negative rate. Finally, we performed GWAS on the different defined subpopulations using a real quantitative trait (resistance to copper sulfate) and showed the feasibility of this approach. The performance of each dataset depends simultaneously on several factors such as sample size, relatedness and population evolutionary history.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

RecombineX: A generalized computational framework for automatic high-throughput gamete genotyping and tetrad-based recombination analysis

Meiotic recombination is an essential biological process that ensures faithful chromosome segregation and promotes parental allele shuffling. Tetrad analysis is a powerful approach to quantify the genetic makeups and recombination landscapes of meiotic products. Here we present RecombineX (https://github.com/yjx1217/RecombineX), a generalized computational framework that automates the full workflow of marker identification, gamete genotyping, and tetrad-based recombination profiling based on any organism or genetic background with batch processing capability. Aside from conventional reference-based analysis, RecombineX can also perform analysis based on parental genome assemblies, which facilitates analyzing meiotic recombination landscapes in their native genomic contexts. Additional features such as copy number variation profiling and missing genotype inference further enhance downstream analysis. RecombineX also includes a dedicate module for simulating the genomes and reads of recombinant tetrads, which enables fine-tuned simulation-based hypothesis testing. This simulation module revealed the power and accuracy of RecombineX even when analyzing tetrads with very low sequencing depths (e.g., 1-2X). Tetrad sequencing data from the budding yeast Saccharomyces cerevisiae and green alga Chlamydomonas reinhardtii were further used to demonstrate the accuracy and robustness of RecombineX for organisms with both small and large genomes, manifesting RecombineX as an all-around one stop solution for future tetrad analysis. Interestingly, our re-analysis of the budding yeast tetrad sequencing data with RecombineX and Oxford Nanopore sequencing revealed two unusual structural rearrangement events that were not noticed before, which exemplify the occasional genome instability triggered by meiosis.

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.

Extensive sampling of Saccharomyces cerevisiae in Taiwan reveals ecology and evolution of predomesticated lineages

The ecology and genetic diversity of the model yeast Saccharomyces cerevisiae before human domestication remain poorly understood. Taiwan is regarded as part of this yeast's geographic birthplace, where the most divergent natural lineage was discovered. Here, we extensively sampled the broadleaf forests across this continental island to probe the ancestral species’ diversity. We found that S. cerevisiae is distributed ubiquitously at low abundance in the forests. Whole-genome sequencing of 121 isolates revealed nine distinct lineages that diverged from Asian lineages during the Pleistocene, when a transient continental shelf land bridge connected Taiwan to other major landmasses. Three lineages are endemic to Taiwan and six are widespread in Asia, making this region a focal biodiversity hotspot. Both ancient and recent admixture events were detected between the natural lineages, and a genetic ancestry component associated with isolates from fruits was detected in most admixed isolates. Collectively, Taiwanese isolates harbor genetic diversity comparable to that of the whole Asia continent, and different lineages have coexisted at a fine spatial scale even on the same tree. Patterns of variations within each lineage revealed that S. cerevisiae is highly clonal and predominantly reproduces asexually in nature. We identified different selection patterns shaping the coding sequences of natural lineages and found fewer gene family expansion and contractions that contrast with domesticated lineages. This study establishes that S. cerevisiae has rich natural diversity sheltered from human influences, making it a powerful model system in microbial ecology.

Yeasts from temperate forests

Yeasts are ubiquitous in temperate forests. While this broad habitat is well-defined, the yeasts inhabiting it and their life cycles, niches, and contributions to ecosystem functioning are less understood. Yeasts are present on nearly all sampled substrates in temperate forests worldwide. They associate with soils, macroorganisms, and other habitats, and no doubt contribute to broader ecosystem-wide processes. Researchers have gathered information leading to hypotheses about yeasts' niches and their life cycles based on physiological observations in the laboratory as well as genomic analyses, but the challenge remains to test these hypotheses in the forests themselves. Here we summarize the habitat and global patterns of yeast diversity, give some information on a handful of well-studied temperate forest yeast genera, discuss the various strategies to isolate forest yeasts, and explain temperate forest yeasts' contributions to biotechnology. We close with a summary of the many future directions and outstanding questions facing researchers in temperate forest yeast ecology. Yeasts present an exciting opportunity to better understand the hidden world of microbial ecology in this threatened and global habitat.

Aborting meiosis allows recombination in sterile diploid yeast hybrids

Hybrids between diverged lineages contain novel genetic combinations but an impaired meiosis often makes them evolutionary dead ends. Here, we explore to what extent an aborted meiosis followed by a return-to-growth (RTG) promotes recombination across a panel of 20 Saccharomyces cerevisiae and S. paradoxus diploid hybrids with different genomic structures and levels of sterility. Genome analyses of 275 clones reveal that RTG promotes recombination and generates extensive regions of loss-of-heterozygosity in sterile hybrids with either a defective meiosis or a heavily rearranged karyotype, whereas RTG recombination is reduced by high sequence divergence between parental subgenomes. The RTG recombination preferentially arises in regions with low local heterozygosity and near meiotic recombination hotspots. The loss-of-heterozygosity has a profound impact on sexual and asexual fitness, and enables genetic mapping of phenotypic differences in sterile lineages where linkage analysis would fail. We propose that RTG gives sterile yeast hybrids access to a natural route for genome recombination and adaptation.

Lager Yeast Design Through Meiotic Segregation of a Saccharomyces cerevisiae × Saccharomyces eubayanus Hybrid

Yeasts in the lager brewing group are closely related and consequently do not exhibit significant genetic variability. Here, an artificial Saccharomyces cerevisiae × Saccharomyces eubayanus tetraploid interspecies hybrid was created by rare mating, and its ability to sporulate and produce viable gametes was exploited to generate phenotypic diversity. Four spore clones obtained from a single ascus were isolated, and their brewing-relevant phenotypes were assessed. These F1 spore clones were found to differ with respect to fermentation performance under lager brewing conditions (15°C, 15 °Plato), production of volatile aroma compounds, flocculation potential and temperature tolerance. One spore clone, selected for its rapid fermentation and acetate ester production was sporulated to produce an F2 generation, again comprised of four spore clones from a single ascus. Again, phenotypic diversity was introduced. In two of these F2 clones, the fermentation performance was maintained and acetate ester production was improved relative to the F1 parent and the original hybrid strain. Strains also performed well in comparison to a commercial lager yeast strain. Spore clones varied in ploidy and chromosome copy numbers, and faster wort fermentation was observed in strains with a higher ploidy. An F2 spore clone was also subjected to 10 consecutive wort fermentations, and single cells were isolated from the resulting yeast slurry. These isolates also exhibited variable fermentation performance and chromosome copy numbers, highlighting the instability of polyploid interspecific hybrids. These results demonstrate the value of this natural approach to increase the phenotypic diversity of lager brewing yeast strains.

Population Size, Sex and Purifying Selection: Comparative Genomics of Two Sister Taxa of the Wild Yeast Saccharomyces paradoxus

This study uses population genomic data to estimate demographic and selection parameters in two sister lineages of the wild yeast Saccharomyces paradoxus and compare their evolution. We first estimate nucleotide and recombinational diversities in each of the two lineages to infer their population size and frequency of sex and then analyze the rate of mutation accumulation since divergence from their inferred common ancestor to estimate the generation time and efficacy of selection. We find that one of the lineages has significantly higher silent nucleotide diversity and lower linkage disequilibrium, indicating a larger population with more frequent sexual generations. The same lineage also shows shorter generation time and higher efficacy of purifying selection, the latter consistent with the finding of larger population size and more frequent sex. Similar analyses are also performed on the ancestries of individual strains within lineages and we find significant differences between strains implying variation in rates of mitotic cell divisions. Our sample includes some strains originating in the Chernobyl nuclear-accident exclusion zone, which has been subjected to high levels of radiation for nearly 30 years now. We find no evidence, however, for increased rates of mutation. Finally, there is a positive correlation between rates of mutation accumulation and length of growing period, as measured by latitude of the place of origin of strains. Our study illustrates the power of genomic analyses in estimating population and life history parameters and testing predictions based on population genetic theory.

Natural variants suppress mutations in hundreds of essential genes

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature‐sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain‐of‐function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature‐sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

The histone H3K9M mutation synergizes with H3K14 ubiquitylation to selectively sequester histone H3K9 methyltransferase Clr4 at heterochromatin

Oncogenic histone lysine-to-methionine mutations block the methylation of their corresponding lysine residues on wild-type histones. One attractive model is that these mutations sequester histone methyltransferases, but genome-wide studies show that mutant histones and histone methyltransferases often do not colocalize. Using chromatin immunoprecipitation sequencing (ChIP-seq), here, we show that, in fission yeast, even though H3K9M-containing nucleosomes are broadly distributed across the genome, the histone H3K9 methyltransferase Clr4 is mainly sequestered at pericentric repeats. This selective sequestration of Clr4 depends not only on H3K9M but also on H3K14 ubiquitylation (H3K14ub), a modification deposited by a Clr4-associated E3 ubiquitin ligase complex. In vitro, H3K14ub synergizes with H3K9M to interact with Clr4 and potentiates the inhibitory effects of H3K9M on Clr4 enzymatic activity. Moreover, binding kinetics show that H3K14ub overcomes the Clr4 aversion to H3K9M and reduces its dissociation. The selective sequestration model reconciles previous discrepancies and demonstrates the importance of protein-interaction kinetics in regulating biological processes.

The budding yeast life cycle: More complex than anticipated?

The budding yeast, Saccharomyces cerevisiae, has served as a model for nearly a century to understand the principles of the eukaryotic life cycle. The canonical life cycle of S. cerevisiae comprises a regular alternation between haploid and diploid phases. Haploid gametes generated by sporulation are expected to quickly restore the diploid phase mainly through inbreeding via intratetrad mating or haploselfing, thereby promoting genome homozygotization. However, recent large population genomics data unveiled that heterozygosity and polyploidy are unexpectedly common. This raises the interesting paradox of a haplo-diplobiontic species being well-adapted to inbreeding and able to maintain high levels of heterozygosity and polyploidy, thereby suggesting an unanticipated complexity of the yeast life cycle. Here, we propose that unprogrammed mating type switching, heterothallism, reduced spore formation and viability, cell-cell fusion and dioecy could play key and uncharted contributions to generate and maintain heterozygosity through polyploidization.

Meiotic Cells Counteract Programmed Retrotransposon Activation via RNA-Binding Translational Repressor Assemblies

Retrotransposon proliferation poses a threat to germline integrity. While retrotransposons must be activated in developing germ cells in order to survive and propagate, how they are selectively activated in the context of meiosis is unclear. We demonstrate that the transcriptional activation of Ty3/Gypsy retrotransposons and host defense are controlled by master meiotic regulators. We show that budding yeast Ty3/Gypsy co-opts binding sites of the essential meiotic transcription factor Ndt80 upstream of the integration site, thereby tightly linking its transcriptional activation to meiotic progression. We also elucidate how yeast cells thwart Ty3/Gypsy proliferation by blocking translation of the retrotransposon mRNA using amyloid-like assemblies of the RNA-binding protein Rim4. In mammals, several inactive Ty3/Gypsy elements are undergoing domestication. We show that mammals utilize equivalent master meiotic regulators (Stra8, Mybl1, Dazl) to regulate Ty3/Gypsy-derived genes in developing gametes. Our findings inform how genes that are evolving from retrotransposons can build upon existing regulatory networks during domestication.

Slow Growth and Increased Spontaneous Mutation Frequency in Respiratory Deficient afo1- Yeast Suppressed by a Dominant Mutation in ATP3

A yeast deletion mutation in the nuclear-encoded gene, AFO1, which codes for a mitochondrial ribosomal protein, led to slow growth on glucose, the inability to grow on glycerol or ethanol, and loss of mitochondrial DNA and respiration. We noticed that afo1- yeast readily obtains secondary mutations that suppress aspects of this phenotype, including its growth defect. We characterized and identified a dominant missense suppressor mutation in the ATP3 gene. Comparing isogenic slowly growing rho-zero and rapidly growing suppressed afo1- strains under carefully controlled fermentation conditions showed that energy charge was not significantly different between strains and was not causal for the observed growth properties. Surprisingly, in a wild-type background, the dominant suppressor allele of ATP3 still allowed respiratory growth but increased the petite frequency. Similarly, a slow-growing respiratory deficient afo1- strain displayed an about twofold increase in spontaneous frequency of point mutations (comparable to the rho-zero strain) while the suppressed strain showed mutation frequency comparable to the respiratory-competent WT strain. We conclude, that phenotypes that result from afo1- are mostly explained by rapidly emerging mutations that compensate for the slow growth that typically follows respiratory deficiency.

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of guide RNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Finally, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

simuG: a general-purpose genome simulator

SUMMARY

Simulated genomes with pre-defined and random genomic variants can be very useful for benchmarking genomic and bioinformatics analyses. Here we introduce simuG, a lightweight tool for simulating the full-spectrum of genomic variants (single nucleotide polymorphisms, Insertions/Deletions, copy number variants, inversions and translocations) for any organisms (including human). The simplicity and versatility of simuG make it a unique general-purpose genome simulator for a wide-range of simulation-based applications.

AVAILABILITY AND IMPLEMENTATION

Code in Perl along with user manual and testing data is available at https://github.com/yjx1217/simuG. This software is free for use under the MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential Gene Expression and Allele Frequency Changes Favour Adaptation of a Heterogeneous Yeast Population to Nitrogen-Limited Fermentations

Alcoholic fermentation is fundamentally an adaptation process, in which the yeast Saccharomyces cerevisiae outperforms its competitors and takes over the fermentation process itself. Although wine yeast strains appear to be adapted to the stressful conditions of alcoholic fermentation, nitrogen limitations in grape must cause stuck or slow fermentations, generating significant economic losses for the wine industry. One way to discover the genetic bases that promote yeast adaptation to nitrogen-deficient environments are selection experiments, where a yeast population undergoes selection under conditions of nitrogen restriction for a number of generations, to then identify by sequencing the molecular characteristics that promote this adaptation. In this work, we carried out selection experiments in bioreactors imitating wine fermentation under nitrogen-limited fermentation conditions (SM60), using the heterogeneous SGRP-4X yeast population, to then sequence the transcriptome and the genome of the population at different time points of the selection process. The transcriptomic results showed an overexpression of genes from the NA strain (North American/YPS128), a wild, non-domesticated isolate. In addition, genome sequencing and allele frequency results allowed several QTLs to be mapped for adaptation to nitrogen-limited fermentation. Finally, we validated the ECM38 allele of NA strain as responsible for higher growth efficiency under nitrogen-limited conditions. Taken together, our results revealed a complex pattern of molecular signatures favouring adaptation of the yeast population to nitrogen-limited fermentations, including differential gene expression, allele frequency changes and loss of the mitochondrial genome. Finally, the results suggest that wild alleles from a non-domesticated isolate (NA) may have a relevant role in the adaptation to the assayed fermentation conditions, with the consequent potential of these alleles for the genetic improvement of wine yeast strains.

An Out-of-Patagonia migration explains the worldwide diversity and distribution of Saccharomyces eubayanus lineages

Population‐level sampling and whole‐genome sequences of different individuals allow one to identify signatures of hybridization, gene flow and potential molecular mechanisms of environmental responses. Here, we report the isolation of 160 Saccharomyces eubayanus strains, the cryotolerant ancestor of lager yeast, from ten sampling sites in Patagonia along 2,000 km of Nothofagus forests. Frequency of S. eubayanus isolates was higher towards southern and colder regions, demonstrating the cryotolerant nature of the species. We sequenced the genome of 82 strains and, together with 23 available genomes, performed a comprehensive phylogenetic analysis. Our results revealed the presence of five different lineages together with dozens of admixed strains. Various analytical methods reveal evidence of gene flow and historical admixture between lineages from Patagonia and Holarctic regions, suggesting the co-occurrence of these ancestral populations. Analysis of the genetic contribution to the admixed genomes revealed a Patagonian genetic origin of the admixed strains, even for those located in the North Hemisphere. Overall, the Patagonian lineages, particularly the southern populations, showed a greater global genetic diversity compared to Holarctic and Chinese lineages, in agreement with a higher abundance in Patagonia. Thus, our results are consistent with a likely colonization of the species from peripheral glacial refugia from South Patagonia. Furthermore, fermentative capacity and maltose consumption resulted negatively correlated with latitude, indicating better fermentative performance in northern populations. Our genome analysis, together with previous reports in the sister species S. uvarum suggests that a S. eubayanus ancestor was adapted to the harsh environmental conditions of Patagonia, a region that provides the ecological conditions for the diversification of these ancestral lineages.

Lager yeast history has intrigued scientists for decades. The recent isolation of S. eubayanus, the lager yeast ancestor, represents an unprecedented opportunity to extend our knowledge on yeast phylogeography and the origins of the S. pastorianus lager hybrid. However, the genetic, phenotypic and evolutionary history of this species remains poorly known. Our work demonstrates that S. eubayanus isolates from Patagonia have the greatest genetic diversity, comprising the largest number of lineages within a single geographic region and experienced ancestral and recent admixture between lineages, likely suggesting co-occurrence in Patagonia. Importantly, some isolates exhibited significant phenotypic differences for traits such as high temperature and ethanol tolerance, together with fermentation performance, demonstrating their potential in the brewing industry for the generation of new styles of lager beers. Furthermore, our results support the idea of colonization from peripheral glacial refugia from the South, as responsible for the high genetic diversity observed in southern Chilean Patagonia. Our results allow hypothesizing a successful physiological adjustment of the species to the local conditions in Patagonia, explaining its wide distribution in the southern hemisphere.

Intragenic repeat expansion in the cell wall protein gene HPF1 controls yeast chronological aging

Aging varies among individuals due to both genetics and environment, but the underlying molecular mechanisms remain largely unknown. Using a highly recombined Saccharomyces cerevisiae population, we found 30 distinct quantitative trait loci (QTLs) that control chronological life span (CLS) in calorie-rich and calorie-restricted environments and under rapamycin exposure. Calorie restriction and rapamycin extended life span in virtually all genotypes but through different genetic variants. We tracked the two major QTLs to the cell wall glycoprotein genes FLO11 and HPF1 We found that massive expansion of intragenic tandem repeats within the N-terminal domain of HPF1 was sufficient to cause pronounced life span shortening. Life span impairment by HPF1 was buffered by rapamycin but not by calorie restriction. The HPF1 repeat expansion shifted yeast cells from a sedentary to a buoyant state, thereby increasing their exposure to surrounding oxygen. The higher oxygenation altered methionine, lipid, and purine metabolism, and inhibited quiescence, which explains the life span shortening. We conclude that fast-evolving intragenic repeat expansions can fundamentally change the relationship between cells and their environment with profound effects on cellular lifestyle and longevity.

Transcriptional Activity and Protein Levels of Horizontally Acquired Genes in Yeast Reveal Hallmarks of Adaptation to Fermentative Environments

In the past decade, the sequencing of large cohorts of Saccharomyces cerevisiae strains has revealed a landscape of genomic regions acquired by Horizontal Gene Transfer (HGT). The genes acquired by HGT play important roles in yeast adaptation to the fermentation process, improving nitrogen and carbon source utilization. However, the functional characterization of these genes at the molecular level has been poorly attended. In this work, we carried out a systematic analysis of the promoter activity and protein level of 30 genes contained in three horizontally acquired regions commonly known as regions A, B, and C. In three strains (one for each region), we used the luciferase reporter gene and the mCherry fluorescent protein to quantify the transcriptional and translational activity of these genes, respectively. We assayed the strains generated in four different culture conditions; all showed low levels of transcriptional and translational activity across these environments. However, we observed an increase in protein levels under low nitrogen culture conditions, suggesting a possible role of the horizontally acquired genes in the adaptation to nitrogen-limited environments. Furthermore, since the strains carrying the luciferase reporter gene are null mutants for the horizontally acquired genes, we assayed growth parameters (latency time, growth rate, and efficiency) and the fermentation kinetics in this set of deletion strains. The results showed that single deletion of 20 horizontally acquired genes modified the growth parameters, whereas the deletion of five of them altered the maximal CO₂ production rate (Vmax). Interestingly, we observed a correlation between growth parameters and Vmax for an ORF within region A, encoding an ortholog to a thiamine (vitamin B1) transporter whose deletion decreased the growth rate, growth efficiency, and CO₂ production. Altogether, our results provided molecular and phenotypic evidence highlighting the importance of horizontally acquired genes in yeast adaptation to fermentative environments.

Human RAP1 specifically protects telomeres of senescent cells from DNA damage

Repressor/activator protein 1 (RAP1) is a highly evolutionarily conserved protein found at telomeres. Although yeast Rap1 is a key telomere capping protein preventing non-homologous end joining (NHEJ) and consequently telomere fusions, its role at mammalian telomeres in vivo is still controversial. Here, we demonstrate that RAP1 is required to protect telomeres in replicative senescent human cells. Downregulation of RAP1 in these cells, but not in young or dividing pre-senescent cells, leads to telomere uncapping and fusions. The anti-fusion effect of RAP1 was further explored in a HeLa cell line where RAP1 expression was depleted through an inducible CRISPR/Cas9 strategy. Depletion of RAP1 in these cells gives rise to telomere fusions only when telomerase is inhibited. We further show that the fusions triggered by RAP1 loss are dependent upon DNA ligase IV. We conclude that human RAP1 is specifically involved in protecting critically short telomeres. This has important implications for the functions of telomeres in senescent cells.

Accurate Tracking of the Mutational Landscape of Diploid Hybrid Genomes

Mutations, recombinations, and genome duplications may promote genetic diversity and trigger evolutionary processes. However, quantifying these events in diploid hybrid genomes is challenging. Here, we present an integrated experimental and computational workflow to accurately track the mutational landscape of yeast diploid hybrids (MuLoYDH) in terms of single-nucleotide variants, small insertions/deletions, copy-number variants, aneuploidies, and loss-of-heterozygosity. Pairs of haploid Saccharomyces parents were combined to generate ancestor hybrids with phased genomes and varying levels of heterozygosity. These diploids were evolved under different laboratory protocols, in particular mutation accumulation experiments. Variant simulations enabled the efficient integration of competitive and standard mapping of short reads, depending on local levels of heterozygosity. Experimental validations proved the high accuracy and resolution of our computational approach. Finally, applying MuLoYDH to four different diploids revealed striking genetic background effects. Homozygous Saccharomyces cerevisiae showed a ∼4-fold higher mutation rate compared with its closely related species S. paradoxus. Intraspecies hybrids unveiled that a substantial fraction of the genome (∼250 bp per generation) was shaped by loss-of-heterozygosity, a process strongly inhibited in interspecies hybrids by high levels of sequence divergence between homologous chromosomes. In contrast, interspecies hybrids exhibited higher single-nucleotide mutation rates compared with intraspecies hybrids. MuLoYDH provided an unprecedented quantitative insight into the evolutionary processes that mold diploid yeast genomes and can be generalized to other genetic systems.

Comparison of Phylogenetic Tree Topologies for Nitrogen Associated Genes Partially Reconstruct the Evolutionary History of Saccharomyces cerevisiae

Massive sequencing projects executed in Saccharomyces cerevisiae have revealed in detail its population structure. The recent "1002 yeast genomes project" has become the most complete catalogue of yeast genetic diversity and a powerful resource to analyse the evolutionary history of genes affecting specific phenotypes. In this work, we selected 22 nitrogen associated genes and analysed the sequence information from the 1011 strains of the "1002 yeast genomes project". We constructed a total evidence (TE) phylogenetic tree using concatenated information, which showed a 27% topology similarity with the reference (REF) tree of the "1002 yeast genomes project". We also generated individual phylogenetic trees for each gene and compared their topologies, identifying genes with similar topologies (suggesting a shared evolutionary history). Furthermore, we pruned the constructed phylogenetic trees to compare the REF tree topology versus the TE tree and the individual genes trees, considering each phylogenetic cluster/subcluster within the population, observing genes with cluster/subcluster topologies of high similarity to the REF tree. Finally, we used the pruned versions of the phylogenetic trees to compare four strains considered as representatives of S. cerevisiae clean lineages, observing for 15 genes that its cluster topologies match 100% the REF tree, supporting that these strains represent main lineages of yeast population. Altogether, our results showed the potential of tree topologies comparison for exploring the evolutionary history of a specific group of genes.

Genetic variants of TORC1 signaling pathway affect nitrogen consumption in Saccharomyces cerevisiae during alcoholic fermentation

In the alcoholic fermentation process, Saccharomyces cerevisiae strains present differences in their nitrogen consumption profiles, these phenotypic outcomes have complex genetic and molecular architectures. In this sense, variations in nitrogen signaling pathways regulated by TORC1 represent one of the main sources of phenotypic diversity in nitrogen consumption. This emphasizes the possible roles that allelic variants from the TORC1 pathway have in the nitrogen consumption differences observed in yeast during the alcoholic fermentation. Here, we studied the allelic diversity in the TORC1 pathway across four yeast strains and determined how these polymorphisms directly impact nitrogen consumption during alcoholic fermentation. Using a reciprocal hemizygosity approach combined with phenotyping under fermentative conditions, we found that allelic variants of GTR1, TOR2, SIT4, SAP185, EAP1, NPR1 and SCH9 underlie differences in the ammonium and amino acids consumption phenotypes. Among these, GTR1 alleles from the Wine/European and West African genetic backgrounds showed the greatest effects on ammonium and amino acid consumption, respectively. Furthermore, we identified allelic variants of SAP185, TOR2, SCH9 and NPR1 from an oak isolate that increased the amino acid consumption preference over ammonium; representing putative candidates coming from a non-domesticated strain that could be used for genetic improvement programs. In conclusion, our results demonstrated that a large number of allelic variants within the TORC1 pathway significantly impacts on regulatory mechanisms of nitrogen assimilation during alcoholic fermentation.

Shared Molecular Targets Confer Resistance over Short and Long Evolutionary Timescales

Pre-existing and de novo genetic variants can both drive adaptation to environmental changes, but their relative contributions and interplay remain poorly understood. Here we investigated the evolutionary dynamics in drug-treated yeast populations with different levels of pre-existing variation by experimental evolution coupled with time-resolved sequencing and phenotyping. We found a doubling of pre-existing variation alone boosts the adaptation by 64.1% and 51.5% in hydroxyurea and rapamycin, respectively. The causative pre-existing and de novo variants were selected on shared targets: RNR4 in hydroxyurea and TOR1, TOR2 in rapamycin. Interestingly, the pre-existing and de novo TOR variants map to different functional domains and act via distinct mechanisms. The pre-existing TOR variants from two domesticated strains exhibited opposite rapamycin resistance effects, reflecting lineage-specific functional divergence. This study provides a dynamic view on how pre-existing and de novo variants interactively drive adaptation and deepens our understanding of clonally evolving populations.

André Goffeau's imprinting on second generation yeast "genomologists"

All authors of the present paper have worked in labs that participated to the sequencing effort of the Saccharomyces cerevisiae reference genome, and we owe to this the fact that we have all chosen to work on genomics of yeasts. S. cerevisiae has been a popular model species for genetics since the 20th century as well as being a model for general eukaryotic cellular processes. Although it has also been used empirically in fermentation for millennia, there was until recently, a lack of knowledge about the natural and evolutionary history of this yeast. The achievement of the international effort to sequence its genome was the foundation for understanding many eukaryotic biological processes but also represented the first step towards the study of the genome and ecological diversity of yeast populations worldwide. We will describe recent advances in yeast comparative and population genomics that find their origins in the S. cerevisiae genome project initiated and pursued by André Goffeau.

Clonal Heterogeneity Influences the Fate of New Adaptive Mutations

The joint contribution of pre-existing and de novo genetic variation to clonal adaptation is poorly understood but essential to designing successful antimicrobial or cancer therapies. To address this, we evolve genetically diverse populations of budding yeast, S. cerevisiae, consisting of diploid cells with unique haplotype combinations. We study the asexual evolution of these populations under selective inhibition with chemotherapeutic drugs by time-resolved whole-genome sequencing and phenotyping. All populations undergo clonal expansions driven by de novo mutations but remain genetically and phenotypically diverse. The clones exhibit widespread genomic instability, rendering recessive de novo mutations homozygous and refining pre-existing variation. Finally, we decompose the fitness contributions of pre-existing and de novo mutations by creating a large recombinant library of adaptive mutations in an ensemble of genetic backgrounds. Both pre-existing and de novo mutations substantially contribute to fitness, and the relative fitness of pre-existing variants sets a selective threshold for new adaptive mutations.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

Isolation and Laboratory Domestication of Natural Yeast Strains

The process from yeast isolation to their use in laboratory experiments is lengthy. Historically, Saccharomyces strains were easily obtained by sampling alcoholic fermentation processes or other substrates associated with human activity in which Saccharomyces was heavily enriched. In contrast, wild Saccharomyces yeasts are found in complex microbial communities and small population sizes, making isolation challenging. We have overcome this problem by enriching yeast on media favoring the growth of Saccharomyces over other microorganisms. The isolation process is usually followed by molecular characterization that allows the strain identification. Finally, yeast isolated from domestic or wild environments need to be genetically manipulated before they can be used in laboratory experiments.

Yeast Reciprocal Hemizygosity to Confirm the Causality of a Quantitative Trait Loci-Associated Gene

Pinpointing causal alleles within a quantitative trait loci region is a key challenge when dissecting the genetic basis of natural variation. In yeast, homing in on culprit genes is often achieved using engineered reciprocal hemizygotes as outlined here. Based on prior information on gene-trait associations, candidate genes are identified. In haploid versions of both founder strains, a candidate gene is then deleted. Gene knockouts are independently mated to a wild-type version of the other strain, such that two diploid hybrid strains are obtained. These strains are identical with regard to the nuclear genome, except for that they are hemizygotic at the locus of interest and contain different alleles of the candidate gene. If correctly measured, a trait difference between these reciprocal hemizygotes can confidently be ascribed to allelic variation at the target locus.

Mapping Quantitative Trait Loci in Yeast

Natural Saccharomyces strains isolated from the wild differ quantitatively in molecular and organismal phenotypes. Quantitative trait loci (QTL) mapping is a powerful approach for identifying sequence variants that alter gene function. In yeast, QTL mapping has been used in designed crosses to map functional polymorphisms. This approach, outlined here, is often the first step in understanding the molecular basis of quantitative traits. New large-scale sequencing surveys have the potential to directly associate genotypes with organismal phenotypes, providing a broader catalog of causative genetic variants. Additional analysis of intermediate phenotypes (e.g., RNA, protein, or metabolite levels) can produce a multilayered and integrated view of individual variation, producing a high-resolution view of the genotype-phenotype map.

Budding Yeast Strains and Genotype-Phenotype Mapping

A small number of well-studied laboratory strains of Saccharomyces cerevisiae, mostly derived from S288C, are used in yeast research. Although powerful, studies for understanding S288C do not always capture the phenotypic essence or the genetic complexity of S. cerevisiae biology. This is particularly problematic for multilocus phenotypes identified in laboratory strains because these loci have never been jointly exposed to natural selection and the corresponding phenotypes do not represent optimization for any particular purpose or environment. Isolation and sequencing of new natural yeast strains also reveal that the total sequence diversity of the S. cerevisiae global population is poorly sampled in common laboratory strains. Here we discuss methodologies required for using the natural genetic variation in yeast to complete a genotype-phenotype map.

de novo assembly and population genomic survey of natural yeast isolates with the Oxford Nanopore MinION sequencer

Background: Oxford Nanopore Technologies Ltd (Oxford, UK) have recently commercialized MinION, a small single-molecule nanopore sequencer, that offers the possibility of sequencing long DNA fragments from small genomes in a matter of seconds. The Oxford Nanopore technology is truly disruptive; it has the potential to revolutionize genomic applications due to its portability, low cost, and ease of use compared with existing long reads sequencing technologies. The MinION sequencer enables the rapid sequencing of small eukaryotic genomes, such as the yeast genome. Combined with existing assembler algorithms, near complete genome assemblies can be generated and comprehensive population genomic analyses can be performed. Results: Here, we resequenced the genome of the Saccharomyces cerevisiae S288C strain to evaluate the performance of nanopore-only assemblers. Then we de novo sequenced and assembled the genomes of 21 isolates representative of the S. cerevisiae genetic diversity using the MinION platform. The contiguity of our assemblies was 14 times higher than the Illumina-only assemblies and we obtained one or two long contigs for 65 % of the chromosomes. This high contiguity allowed us to accurately detect large structural variations across the 21 studied genomes. Conclusion: Because of the high completeness of the nanopore assemblies, we were able to produce a complete cartography of transposable elements insertions and inspect structural variants that are generally missed using a short-read sequencing strategy. Our analyses show that the Oxford Nanopore technology is already usable for de novo sequencing and assembly; however, non-random errors in homopolymers require polishing the consensus using an alternate sequencing technology.

Contrasting evolutionary genome dynamics between domesticated and wild yeasts

Structural rearrangements have long been recognized as an important source of genetic variation, with implications in phenotypic diversity and disease, yet their detailed evolutionary dynamics remain elusive. Here we use long-read sequencing to generate end-to-end genome assemblies for 12 strains representing major subpopulations of the partially domesticated yeast Saccharomyces cerevisiae and its wild relative Saccharomyces paradoxus. These population-level high-quality genomes with comprehensive annotation enable precise definition of chromosomal boundaries between cores and subtelomeres and a high-resolution view of evolutionary genome dynamics. In chromosomal cores, S. paradoxus shows faster accumulation of balanced rearrangements (inversions, reciprocal translocations and transpositions), whereas S. cerevisiae accumulates unbalanced rearrangements (novel insertions, deletions and duplications) more rapidly. In subtelomeres, both species show extensive interchromosomal reshuffling, with a higher tempo in S. cerevisiae. Such striking contrasts between wild and domesticated yeasts are likely to reflect the influence of human activities on structural genome evolution.

Mismatch Repair Incompatibilities in Diverse Yeast Populations

An elevated mutation rate can provide cells with a source of mutations to adapt to changing environments. We identified a negative epistatic interaction involving naturally occurring variants in the MLH1 and PMS1 mismatch repair (MMR) genes of Saccharomyces cerevisiae We hypothesized that this MMR incompatibility, created through mating between divergent S. cerevisiae, yields mutator progeny that can rapidly but transiently adapt to an environmental stress. Here we analyzed the MLH1 and PMS1 genes across 1010 S. cerevisiae natural isolates spanning a wide range of ecological sources (tree exudates, Drosophila, fruits, and various fermentation and clinical isolates) and geographical sources (Europe, America, Africa, and Asia). We identified one homozygous clinical isolate and 18 heterozygous isolates containing the incompatible MMR genotype. The MLH1-PMS1 gene combination isolated from the homozygous clinical isolate conferred a mutator phenotype when expressed in the S288c laboratory background. Using a novel reporter to measure mutation rates, we showed that the overall mutation rate in the homozygous incompatible background was similar to that seen in compatible strains, indicating the presence of suppressor mutations in the clinical isolate that lowered its mutation rate. This observation and the identification of 18 heterozygous isolates, which can lead to MMR incompatible genotypes in the offspring, are consistent with an elevated mutation rate rapidly but transiently facilitating adaptation. To avoid long-term fitness costs, the incompatibility is apparently buffered by mating or by acquiring suppressors. These observations highlight effective strategies in eukaryotes to avoid long-term fitness costs associated with elevated mutation rates.

The genetic architecture of low-temperature adaptation in the wine yeast Saccharomyces cerevisiae

Background

Low-temperature growth and fermentation of wine yeast can enhance wine aroma and make them highly desirable traits for the industry. Elucidating response to cold in Saccharomyces cerevisiae is, therefore, of paramount importance to select or genetically improve new wine strains. As most enological traits of industrial importance in yeasts, adaptation to low temperature is a polygenic trait regulated by many interacting loci.

Results

In order to unravel the genetic determinants of low-temperature fermentation, we mapped quantitative trait loci (QTLs) by bulk segregant analyses in the F13 offspring of two Saccharomyces cerevisiae industrial strains with divergent performance at low temperature. We detected four genomic regions involved in the adaptation at low temperature, three of them located in the subtelomeric regions (chromosomes XIII, XV and XVI) and one in the chromosome XIV. The QTL analysis revealed that subtelomeric regions play a key role in defining individual variation, which emphasizes the importance of these regions’ adaptive nature.

Conclusions

The reciprocal hemizygosity analysis (RHA), run to validate the genes involved in low-temperature fermentation, showed that genetic variation in mitochondrial proteins, maintenance of correct asymmetry and distribution of phospholipid in the plasma membrane are key determinants of low-temperature adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3572-2) contains supplementary material, which is available to authorized users.