The molecular basis of phenotypic variation in yeast

Trans-eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae. We discovered thousands of genetic correlations between gene expression and growth, suggesting potential functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans-acting regulatory loci. Trans-acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits, but with heterogeneous direction of effect. Our results reveal how trans-acting regulatory hotspots shape complex traits by altering cellular states.

Role of genetic architecture in phenotypic plasticity

Phenotypic plasticity, the ability of an organism to display different phenotypes across environments, is widespread in nature. Plasticity aids survival in novel environments. Herein, we review studies from yeast that allow us to start uncovering the genetic architecture of phenotypic plasticity. Genetic variants and their interactions impact the phenotype in different environments, and distinct environments modulate the impact of genetic variants and their interactions on the phenotype. Because of this, certain hidden genetic variation is expressed in specific genetic and environmental backgrounds. A better understanding of the genetic mechanisms of phenotypic plasticity will help to determine short- and long-term responses to selection and how wide variation in disease manifestation occurs in human populations.

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

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

High-throughput approaches to functional characterization of genetic variation in yeast

Expansion of sequencing efforts to include thousands of genomes is providing a fundamental resource for determining the genetic diversity that exists in a population. Now, high-throughput approaches are necessary to begin to understand the role these genotypic changes play in affecting phenotypic variation. Saccharomyces cerevisiae maintains its position as an excellent model system to determine the function of unknown variants with its exceptional genetic diversity, phenotypic diversity, and reliable genetic manipulation tools. Here, we review strategies and techniques developed in yeast that scale classic approaches of assessing variant function. These approaches improve our ability to better map quantitative trait loci at a higher resolution, even for rare variants, and are already providing greater insight into the role that different types of mutations play in phenotypic variation and evolution not just in yeast but across taxa.

Use of the Iron-Responsive RBT5 Promoter for Regulated Expression in Candida albicans

Engineered conditional gene expression is used in appraisal of gene function and pathway relationships. For pathogens like the fungus Candida albicans, conditional expression systems are most useful if they are active in the infection environment and if they can be utilized in multiple clinical isolates. Here, we describe such a system. It employs the RBT5 promoter and can be implemented with a few PCRs. We validated the system with RBT5 promoter fusions to two genes that promote filamentation and polarized growth, UME6 and HGC1, and with efg1Δ/Δ mutants, which are defective in an activator of filamentous growth. An RBT5 promoter fusion to either gene enabled filamentous growth of an efg1Δ/Δ mutant of strain SC5314 in iron-limited media, including RPMI with serum and yeast extract-peptone-dextrose with bathophenanthrolinedisulfonic acid. The RBT5-UME6 fusion promoted filamentation of efg1Δ/Δ mutants in RPMI with serum of four other clinical C. albicans isolates as well. In a mouse model of disseminated candidiasis, the RBT5-UME6 fusion promoted filamentation of the SC5314 efg1Δ/Δ mutant in kidney tissue, an indication that the RBT5 promoter is active in the iron-limited host environment. The RBT5 promoter expands the conditional expression toolkit for C. albicans genetics. IMPORTANCE Genetic strategies have been vital for mechanistic analysis of biological processes. Here, we describe a genetic tool for the fungal pathogen Candida albicans.

Going wild: ecology and genomics are crucial to understand yeast evolution

Improved and more accessible genome-sequencing approaches have allowed the analysis of large sets of natural yeast isolates. As a consequence, this unprecedented level of description of yeast-genome characteristics and variations in natural environments has provided crucial insights on yeast ecology and evolution. Here, we review some of the most relevant and intriguing aspects of yeast evolution pointed out, thanks to the combination of yeast ecology and genomics, and critically examine the resulting improvement of our knowledge on this field. Only integrated approaches, taking into consideration not only the characteristics of the microbe but also those of the hosting environment, will significantly move forward the exploration of yeast diversity, ecology, and evolution.

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.

Barcoded Bulk QTL mapping reveals highly polygenic and epistatic architecture of complex traits in yeast

Mapping the genetic basis of complex traits is critical to uncovering the biological mechanisms that underlie disease and other phenotypes. Genome-wide association studies (GWAS) in humans and quantitative trait locus (QTL) mapping in model organisms can now explain much of the observed heritability in many traits, allowing us to predict phenotype from genotype. However, constraints on power due to statistical confounders in large GWAS and smaller sample sizes in QTL studies still limit our ability to resolve numerous small-effect variants, map them to causal genes, identify pleiotropic effects across multiple traits, and infer non-additive interactions between loci (epistasis). Here, we introduce barcoded bulk quantitative trait locus (BB-QTL) mapping, which allows us to construct, genotype, and phenotype 100,000 offspring of a budding yeast cross, two orders of magnitude larger than the previous state of the art. We use this panel to map the genetic basis of eighteen complex traits, finding that the genetic architecture of these traits involves hundreds of small-effect loci densely spaced throughout the genome, many with widespread pleiotropic effects across multiple traits. Epistasis plays a central role, with thousands of interactions that provide insight into genetic networks. By dramatically increasing sample size, BB-QTL mapping demonstrates the potential of natural variants in high-powered QTL studies to reveal the highly polygenic, pleiotropic, and epistatic architecture of complex traits.

Predictive Potential of MALDI-TOF Analyses for Wine and Brewing Yeast

The potential of MALDI-TOF profiling for predicting potential applications of yeast strains in the beverage sector was assessed. A panel of 59 commercial yeasts (47 wine and 12 brewing yeasts) was used to validate the concept whereby 2 culture media (YPD agar and YPD broth), as well as two mass ranges m/z 500–4000 and m/z 2000–20,000, were evaluated for the best fit. Three machine learning-based algorithms, PCA, MDS, and UMAP, in addition to a hierarchical clustering method, were employed. Profiles derived from broth cultures yielded more peaks, but these were less well-defined compared with those from agar cultures. Hierarchical clustering more clearly resolved different species and gave a broad overview of potential strain utility, but more nuanced insights were provided by MDS and UMAP analyses. PCA-based displays were less informative. The potential of MALDI-TOF proteomics in predicting the utility of yeast strains of commercial benefit is supported in this study, provided appropriate approaches are used for data generation and analysis.

Multiple epistatic DNA variants in a single gene affect gene expression in trans

DNA variants that alter gene expression in trans are important sources of phenotypic variation. Nevertheless, the identity of trans-acting variants remains poorly understood. Single causal variants in several genes have been reported to affect the expression of numerous distant genes in trans. Whether these simple molecular architectures are representative of trans-acting variation is unknown. Here, we studied the large RAS signaling regulator gene IRA2, which contains variants with extensive trans-acting effects on gene expression in the yeast Saccharomyces cerevisiae. We used systematic CRISPR-based genome engineering and a sensitive phenotyping strategy to dissect causal variants to the nucleotide level. In contrast to the simple molecular architectures known so far, IRA2 contained at least seven causal nonsynonymous variants. The effects of these variants were modulated by nonadditive, epistatic interactions. Two variants at the 5'-end affected gene expression and growth only when combined with a third variant that also had no effect in isolation. Our findings indicate that the molecular basis of trans-acting genetic variation may be considerably more complex than previously appreciated.

Genomic characterization of a wild diploid isolate of Saccharomyces cerevisiae reveals an extensive and dynamic landscape of structural variation

The budding yeast Saccharomyces cerevisiae has been extensively characterized for many decades and is a critical resource for the study of numerous facets of eukaryotic biology. Recently, whole genome sequence analysis of over 1000 natural isolates of S. cerevisiae has provided critical insights into the evolutionary landscape of this species by revealing a population structure comprised of numerous genomically diverse lineages. These survey-level analyses have been largely devoid of structural genomic information, mainly because short read sequencing is not suitable for detailed characterization of genomic architecture. Consequently, we still lack a complete perspective of the genomic variation the exists within the species. Single molecule long read sequencing technologies, such as Oxford Nanopore and PacBio, provide sequencing-based approaches with which to rigorously define the structure of a genome, and have empowered yeast geneticists to explore this poorly described realm of eukaryotic genomics. Here, we present the comprehensive genomic structural analysis of a wild diploid isolate of S. cerevisiae, YJM311. We used long read sequence analysis to construct a haplotype-phased, telomere-to-telomere length assembly of the YJM311 genome and characterized the structural variations (SVs) therein. We discovered that the genome of YJM311 contains significant intragenomic structural variation, some of which imparts notable consequences to the genomic stability and developmental biology of the strain. Collectively, we outline a new methodology for creating accurate haplotype-phased genome assemblies and highlight how such genomic analyses can define the structural architectures of S. cerevisiae isolates. It is our hope that continued structural characterization of S. cerevisiae genomes, such as we have reported here for YJM311, will comprehensively advance our understanding of eukaryotic genome structure-function relationships, structural genomic diversity, and evolution.

Cis-regulatory variants affect gene expression dynamics in yeast

Evolution of cis-regulatory sequences depends on how they affect gene expression and motivates both the identification and prediction of cis-regulatory variants responsible for expression differences within and between species. While much progress has been made in relating cis-regulatory variants to expression levels, the timing of gene activation and repression may also be important to the evolution of cis-regulatory sequences. We investigated allele-specific expression (ASE) dynamics within and between Saccharomyces species during the diauxic shift and found appreciable cis-acting variation in gene expression dynamics. Within-species ASE is associated with intergenic variants, and ASE dynamics are more strongly associated with insertions and deletions than ASE levels. To refine these associations, we used a high-throughput reporter assay to test promoter regions and individual variants. Within the subset of regions that recapitulated endogenous expression, we identified and characterized cis-regulatory variants that affect expression dynamics. Between species, chimeric promoter regions generate novel patterns and indicate constraints on the evolution of gene expression dynamics. We conclude that changes in cis-regulatory sequences can tune gene expression dynamics and that the interplay between expression dynamics and other aspects of expression is relevant to the evolution of cis-regulatory sequences.

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.

Pervasive Phenotypic Impact of a Large Nonrecombining Introgressed Region in Yeast

To explore the origin of the diversity observed in natural populations, many studies have investigated the relationship between genotype and phenotype. In yeast species, especially in Saccharomyces cerevisiae, these studies are mainly conducted using recombinant offspring derived from two genetically diverse isolates, allowing to define the phenotypic effect of genetic variants. However, large genomic variants such as interspecies introgressions are usually overlooked even if they are known to modify the genotype–phenotype relationship. To have a better insight into the overall phenotypic impact of introgressions, we took advantage of the presence of a 1-Mb introgressed region, which lacks recombination and contains the mating-type determinant in the Lachancea kluyveri budding yeast. By performing linkage mapping analyses in this species, we identified a total of 89 loci affecting growth fitness in a large number of conditions and 2,187 loci affecting gene expression mostly grouped into two major hotspots, one being the introgressed region carrying the mating-type locus. Because of the absence of recombination, our results highlight the presence of a sexual dimorphism in a budding yeast for the first time. Overall, by describing the phenotype–genotype relationship in the Lachancea kluyveri species, we expanded our knowledge on how genetic characteristics of large introgression events can affect the phenotypic landscape.

Simultaneous quantification of mRNA and protein in single cells reveals post-transcriptional effects of genetic variation

Trans-acting DNA variants may specifically affect mRNA or protein levels of genes located throughout the genome. However, prior work compared trans-acting loci mapped in separate studies, many of which had limited statistical power. Here, we developed a CRISPR-based system for simultaneous quantification of mRNA and protein of a given gene via dual fluorescent reporters in single, live cells of the yeast Saccharomyces cerevisiae. In large populations of recombinant cells from a cross between two genetically divergent strains, we mapped 86 trans-acting loci affecting the expression of ten genes. Less than 20% of these loci had concordant effects on mRNA and protein of the same gene. Most loci influenced protein but not mRNA of a given gene. One locus harbored a premature stop variant in the YAK1 kinase gene that had specific effects on protein or mRNA of dozens of genes. These results demonstrate complex, post-transcriptional genetic effects on gene expression.

Pan-genomic open reading frames: A potential supplement of single nucleotide polymorphisms in estimation of heritability and genomic prediction

Pan-genomic open reading frames (ORFs) potentially carry protein-coding gene or coding variant information in a population. In this study, we suggest that pan-genomic ORFs are promising to be utilized in estimation of heritability and genomic prediction. A Saccharomyces cerevisiae dataset with whole-genome SNPs, pan-genomic ORFs, and the copy numbers of those ORFs is used to test the effectiveness of ORF data as a predictor in three prediction models for 35 traits. Our results show that the ORF-based heritability can capture more genetic effects than SNP-based heritability for all traits. Compared to SNP-based genomic prediction (GBLUP), pan-genomic ORF-based genomic prediction (OBLUP) is distinctly more accurate for all traits, and the predictive abilities on average are more than doubled across all traits. For four traits, the copy number of ORF-based prediction(CBLUP) is more accurate than OBLUP. When using different numbers of isolates in training sets in ORF-based prediction, the predictive abilities for all traits increased as more isolates are added in the training sets, suggesting that with very large training sets the prediction accuracy will be in the range of the square root of the heritability. We conclude that pan-genomic ORFs have the potential to be a supplement of single nucleotide polymorphisms in estimation of heritability and genomic prediction.

Mapping Ethanol Tolerance in Budding Yeast Reveals High Genetic Variation in a Wild Isolate

Ethanol tolerance, a polygenic trait of the yeast Saccharomyces cerevisiae, is the primary factor determining industrial bioethanol productivity. Until now, genomic elements affecting ethanol tolerance have been mapped only at low resolution, hindering their identification. Here, we explore the genetic architecture of ethanol tolerance, in the F6 generation of an Advanced Intercrossed Line (AIL) mapping population between two phylogenetically distinct, but phenotypically similar, S. cerevisiae strains (a common laboratory strain and a wild strain isolated from nature). Under ethanol stress, 51 quantitative trait loci (QTLs) affecting growth and 96 QTLs affecting survival, most of them novel, were identified, with high resolution, in some cases to single genes, using a High-Resolution Mapping Package of methodologies that provided high power and high resolution. We confirmed our results experimentally by showing the effects of the novel mapped genes: MOG1, MGS1, and YJR154W. The mapped QTLs explained 34% of phenotypic variation for growth and 72% for survival. High statistical power provided by our analysis allowed detection of many loci with small, but mappable effects, uncovering a novel “quasi-infinitesimal” genetic architecture. These results are striking demonstration of tremendous amounts of hidden genetic variation exposed in crosses between phylogenetically separated strains with similar phenotypes; as opposed to the more common design where strains with distinct phenotypes are crossed. Our findings suggest that ethanol tolerance is under natural evolutionary fitness-selection for an optimum phenotype that would tend to eliminate alleles of large effect. The study provides a platform for development of superior ethanol-tolerant strains using genome editing or selection.

DNA variants affecting the expression of numerous genes in trans have diverse mechanisms of action and evolutionary histories

DNA variants that alter gene expression contribute to variation in many phenotypic traits. In particular, trans-acting variants, which are often located on different chromosomes from the genes they affect, are an important source of heritable gene expression variation. However, our knowledge about the identity and mechanism of causal trans-acting variants remains limited. Here, we developed a fine-mapping strategy called CRISPR-Swap and dissected three expression quantitative trait locus (eQTL) hotspots known to alter the expression of numerous genes in trans in the yeast Saccharomyces cerevisiae. Causal variants were identified by engineering recombinant alleles and quantifying the effects of these alleles on the expression of a green fluorescent protein-tagged gene affected by the given locus in trans. We validated the effect of each variant on the expression of multiple genes by RNA-sequencing. The three variants differed in their molecular mechanism, the type of genes they reside in, and their distribution in natural populations. While a missense leucine-to-serine variant at position 63 in the transcription factor Oaf1 (L63S) was almost exclusively present in the reference laboratory strain, the two other variants were frequent among S. cerevisiae isolates. A causal missense variant in the glucose receptor Rgt2 (V539I) occurred at a poorly conserved amino acid residue and its effect was strongly dependent on the concentration of glucose in the culture medium. A noncoding variant in the conserved fatty acid regulated (FAR) element of the OLE1 promoter influenced the expression of the fatty acid desaturase Ole1 in cis and, by modulating the level of this essential enzyme, other genes in trans. The OAF1 and OLE1 variants showed a non-additive genetic interaction, and affected cellular lipid metabolism. These results demonstrate that the molecular basis of trans-regulatory variation is diverse, highlighting the challenges in predicting which natural genetic variants affect gene expression.

Differences in the DNA sequence of individual genomes contribute to differences in many traits, such as appearance, physiology, and the risk for common diseases. An important group of these DNA variants influences how individual genes across the genome are turned on or off. In this paper, we describe a strategy for identifying such “trans-acting” variants in different strains of baker’s yeast. We used this strategy to reveal three single DNA base changes that each influences the expression of dozens of genes. These three DNA variants were very different from each other. Two of them changed the protein sequence, one in a transcription factor and the other in a sugar sensor. The third changed the expression of an enzyme, a change that in turn caused other genes to alter their expression. One variant existed in only a few yeast isolates, while the other two existed in many isolates collected from around the world. This diversity of DNA variants that influence the expression of many other genes illustrates how difficult it is to predict which DNA variants in an individual’s genome will have effects on the organism.

Multiple Changes Underlie Allelic Divergence of CUP2 Between Saccharomyces Species

Under the model of micromutationism, phenotypic divergence between species is caused by accumulation of many small-effect changes. While mapping the causal changes to single nucleotide resolution could be difficult for diverged species, genetic dissection via chimeric constructs allows us to evaluate whether a large-effect gene is composed of many small-effect nucleotide changes. In a previously described non-complementation screen, we found an allele difference of CUP2, a copper-binding transcription factor, underlies divergence in copper resistance between Saccharomyces cerevisiae and S. uvarum. Here, we tested whether the allele effect of CUP2 was caused by multiple nucleotide changes. By analyzing chimeric constructs containing four separate regions in the CUP2 gene, including its distal promoter, proximal promoter, DNA binding domain and transcriptional activation domain, we found that all four regions of the S. cerevisiae allele conferred copper resistance, with the proximal promoter showing the largest effect, and that both additive and epistatic effects are likely involved. These findings support a model of multiple changes underlying evolution and suggest an important role of both protein coding and cis-regulatory changes in evolution.

Extensive impact of low-frequency variants on the phenotypic landscape at population-scale

Genome-wide association studies (GWAS) allow to dissect complex traits and map genetic variants, which often explain relatively little of the heritability. One potential reason is the preponderance of undetected low-frequency variants. To increase their allele frequency and assess their phenotypic impact in a population, we generated a diallel panel of 3025 yeast hybrids, derived from pairwise crosses between natural isolates and examined a large number of traits. Parental versus hybrid regression analysis showed that while most phenotypic variance is explained by additivity, a third is governed by non-additive effects, with complete dominance having a key role. By performing GWAS on the diallel panel, we found that associated variants with low frequency in the initial population are overrepresented and explain a fraction of the phenotypic variance as well as an effect size similar to common variants. Overall, we highlighted the relevance of low-frequency variants on the phenotypic variation.

Rare variants contribute disproportionately to quantitative trait variation in yeast

How variants with different frequencies contribute to trait variation is a central question in genetics. We use a unique model system to disentangle the contributions of common and rare variants to quantitative traits. We generated ~14,000 progeny from crosses among 16 diverse yeast strains and identified thousands of quantitative trait loci (QTLs) for 38 traits. We combined our results with sequencing data for 1011 yeast isolates to show that rare variants make a disproportionate contribution to trait variation. Evolutionary analyses revealed that this contribution is driven by rare variants that arose recently, and that negative selection has shaped the relationship between variant frequency and effect size. We leveraged the structure of the crosses to resolve hundreds of QTLs to single genes. These results refine our understanding of trait variation at the population level and suggest that studies of rare variants are a fertile ground for discovery of genetic effects.

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.

Quantitative Trait Nucleotides Impacting the Technological Performances of Industrial Saccharomyces cerevisiae Strains

The budding yeast Saccharomyces cerevisiae is certainly the prime industrial microorganism and is related to many biotechnological applications including food fermentations, biofuel production, green chemistry, and drug production. A noteworthy characteristic of this species is the existence of subgroups well adapted to specific processes with some individuals showing optimal technological traits. In the last 20 years, many studies have established a link between quantitative traits and single-nucleotide polymorphisms found in hundreds of genes. These natural variations constitute a pool of QTNs (quantitative trait nucleotides) that modulate yeast traits of economic interest for industry. By selecting a subset of genes functionally validated, a total of 284 QTNs were inventoried. Their distribution across pan and core genome and their frequency within the 1,011 Saccharomyces cerevisiae genomes were analyzed. We found that 150 of the 284 QTNs have a frequency lower than 5%, meaning that these variants would be undetectable by genome-wide association studies (GWAS). This analysis also suggests that most of the functional variants are private to a subpopulation, possibly due to their adaptive role to specific industrial environment. In this review, we provide a literature survey of their phenotypic impact and discuss the opportunities and the limits of their use for industrial strain selection.

Experimental Evolution of Yeast for High-Temperature Tolerance

Thermotolerance is a polygenic trait that contributes to cell survival and growth under unusually high temperatures. Although some genes associated with high-temperature growth (Htg+) have been identified, how cells accumulate mutations to achieve prolonged thermotolerance is still mysterious. Here, we conducted experimental evolution of a Saccharomyces cerevisiae laboratory strain with stepwise temperature increases for it to grow at 42 °C. Whole genome resequencing of 14 evolved strains and the parental strain revealed a total of 153 mutations in the evolved strains, including single nucleotide variants, small INDELs, and segmental duplication/deletion events. Some mutations persisted from an intermediate temperature to 42 °C, so they might be Htg+ mutations. Functional categorization of mutations revealed enrichment of exonic mutations in the SWI/SNF complex and F-type ATPase, pointing to their involvement in high-temperature tolerance. In addition, multiple mutations were found in a general stress-associated signal transduction network consisting of Hog1 mediated pathway, RAS-cAMP pathway, and Rho1-Pkc1 mediated cell wall integrity pathway, implying that cells can achieve Htg+ partly through modifying existing stress regulatory mechanisms. Using pooled segregant analysis of five Htg+ phenotype-orientated pools, we inferred causative mutations for growth at 42 °C and identified those mutations with stronger impacts on the phenotype. Finally, we experimentally validated a number of the candidate Htg+ mutations. This study increased our understanding of the genetic basis of yeast tolerance to high temperature.

Pervasive function and evidence for selection across standing genetic variation in S. cerevisiae

Quantitative genetics aims to map genotype to phenotype, often with the goal of understanding how organisms evolved. However, it remains unclear whether the genetic variants identified are exemplary of evolution. Here we analyzed progeny of two wild Saccharomyces cerevisiae isolates to identify 195 loci underlying complex metabolic traits, resolving 107 to single polymorphisms with diverse molecular mechanisms. More than 20% of causal variants exhibited patterns of emergence inconsistent with neutrality. Moreover, contrary to drift-centric expectation, variation in diverse wild yeast isolates broadly exhibited this property: over 30% of shared natural variants exhibited phylogenetic signatures suggesting that they are not neutral. This pattern is likely attributable to both homoplasy and balancing selection on ancestral polymorphism. Variants that emerged repeatedly were more likely to have done so in isolates from the same ecological niche. Our results underscore the power of super-resolution mapping of ecologically relevant traits in understanding adaptation and evolution.

Complex modifier landscape underlying genetic background effects

Genetic background impacts the phenotypic outcome of a mutation in different individuals; however, the underlying molecular mechanisms are often unclear. We characterized genes exhibiting conditional essentiality when mutated in two genetically distinct yeast strains. Hybrid crosses and whole-genome sequencing revealed that conditional essentiality can be associated with nonchromosomal elements or a single-modifier locus, but most involve a complex set of modifier loci. Detailed analysis of the cysteine biosynthesis pathway showed that independent, rare, single-gene modifiers, related to both up- and downstream pathway functions, can arise in multiple allelic forms from separate lineages. For several genes, we also resolved complex sets of modifying loci underlying conditional essentiality, revealing specific genetic interactions that drive an individual strain’s background effect.

The phenotypic consequence of a given mutation can be influenced by the genetic background. For example, conditional gene essentiality occurs when the loss of function of a gene causes lethality in one genetic background but not another. Between two individual Saccharomyces cerevisiae strains, S288c and Σ1278b, ∼1% of yeast genes were previously identified as “conditional essential.” Here, in addition to confirming that some conditional essential genes are modified by a nonchromosomal element, we show that most cases involve a complex set of genomic modifiers. From tetrad analysis of S288C/Σ1278b hybrid strains and whole-genome sequencing of viable hybrid spore progeny, we identified complex sets of multiple genomic regions underlying conditional essentiality. For a smaller subset of genes, including CYS3 and CYS4, each of which encodes components of the cysteine biosynthesis pathway, we observed a segregation pattern consistent with a single modifier associated with conditional essentiality. In natural yeast isolates, we found that the CYS3/CYS4 conditional essentiality can be caused by variation in two independent modifiers, MET1 and OPT1, each with roles associated with cellular cysteine physiology. Interestingly, the OPT1 allelic variation appears to have arisen independently from separate lineages, with rare allele frequencies below 0.5%. Thus, while conditional gene essentiality is usually driven by genetic interactions associated with complex modifier architectures, our analysis also highlights the role of functionally related, genetically independent, and rare variants.

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.

Mitochondrial Genome Variation Affects Multiple Respiration and Nonrespiration Phenotypes in Saccharomyces cerevisiae

Mitochondrial genome variation and its effects on phenotypes have been widely analyzed in higher eukaryotes but less so in the model eukaryote Saccharomyces cerevisiae Here, we describe mitochondrial genome variation in 96 diverse S. cerevisiae strains and assess associations between mitochondrial genotype and phenotypes as well as nuclear-mitochondrial epistasis. We associate sensitivity to the ATP synthase inhibitor oligomycin with SNPs in the mitochondrially encoded ATP6 gene. We describe the use of iso-nuclear F1 pairs, the mitochondrial genome equivalent of reciprocal hemizygosity analysis, to identify and analyze mitochondrial genotype-dependent phenotypes. Using iso-nuclear F1 pairs, we analyze the oligomycin phenotype-ATP6 association and find extensive nuclear-mitochondrial epistasis. Similarly, in iso-nuclear F1 pairs, we identify many additional mitochondrial genotype-dependent respiration phenotypes, for which there was no association in the 96 strains, and again find extensive nuclear-mitochondrial epistasis that likely contributes to the lack of association in the 96 strains. Finally, in iso-nuclear F1 pairs, we identify novel mitochondrial genotype-dependent nonrespiration phenotypes: resistance to cycloheximide, ketoconazole, and copper. We discuss potential mechanisms and the implications of mitochondrial genotype and of nuclear-mitochondrial epistasis effects on respiratory and nonrespiratory quantitative traits.

Comparative Genomics Approaches Accurately Predict Deleterious Variants in Plants

Recent advances in genome resequencing have led to increased interest in prediction of the functional consequences of genetic variants. Variants at phylogenetically conserved sites are of particular interest, because they are more likely than variants at phylogenetically variable sites to have deleterious effects on fitness and contribute to phenotypic variation. Numerous comparative genomic approaches have been developed to predict deleterious variants, but the approaches are nearly always assessed based on their ability to identify known disease-causing mutations in humans. Determining the accuracy of deleterious variant predictions in nonhuman species is important to understanding evolution, domestication, and potentially to improving crop quality and yield. To examine our ability to predict deleterious variants in plants we generated a curated database of 2,910 Arabidopsis thaliana mutants with known phenotypes. We evaluated seven approaches and found that while all performed well, their relative ranking differed from prior benchmarks in humans. We conclude that deleterious mutations can be reliably predicted in A. thaliana and likely other plant species, but that the relative performance of various approaches does not necessarily translate from one species to another.

Functional Genetic Variants Revealed by Massively Parallel Precise Genome Editing

A major challenge in genetics is to identify genetic variants driving natural phenotypic variation. However, current methods of genetic mapping have limited resolution. To address this challenge, we developed a CRISPR-Cas9-based high-throughput genome editing approach that can introduce thousands of specific genetic variants in a single experiment. This enabled us to study the fitness consequences of 16,006 natural genetic variants in yeast. We identified 572 variants with significant fitness differences in glucose media; these are highly enriched in promoters, particularly in transcription factor binding sites, while only 19.2% affect amino acid sequences. Strikingly, nearby variants nearly always favor the same parent's alleles, suggesting that lineage-specific selection is often driven by multiple clustered variants. In sum, our genome editing approach reveals the genetic architecture of fitness variation at single-base resolution and could be adapted to measure the effects of genome-wide genetic variation in any screen for cell survival or cell-sortable markers.

Genetics of trans-regulatory variation in gene expression

Heritable variation in gene expression forms a crucial bridge between genomic variation and the biology of many traits. However, most expression quantitative trait loci (eQTLs) remain unidentified. We mapped eQTLs by transcriptome sequencing in 1012 yeast segregants. The resulting eQTLs accounted for over 70% of the heritability of mRNA levels, allowing comprehensive dissection of regulatory variation. Most genes had multiple eQTLs. Most expression variation arose from trans-acting eQTLs distant from their target genes. Nearly all trans-eQTLs clustered at 102 hotspot locations, some of which influenced the expression of thousands of genes. Fine-mapped hotspot regions were enriched for transcription factor genes. While most genes had a local eQTL, most of these had no detectable effects on the expression of other genes in trans. Hundreds of non-additive genetic interactions accounted for small fractions of expression variation. These results reveal the complexity of genetic influences on transcriptome variation in unprecedented depth and detail.

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.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.

The natural diversity and ecology of fission yeast

While the fission yeast is a powerful model of eukaryote biology, there have been few studies of quantitative genetics, phenotypic or genetic diversity. Here I survey the small collection of fission yeast diversity research. I discuss what we can infer about the ecology and origins of Schizosaccharomyces pombe from microbiology field studies and the few strains that have been collected.

Mapping Causal Variants with Single-Nucleotide Resolution Reveals Biochemical Drivers of Phenotypic Change

Understanding the sequence determinants that give rise to diversity among individuals and species is the central challenge of genetics. However, despite ever greater numbers of sequenced genomes, most genome-wide association studies cannot distinguish causal variants from linked passenger mutations spanning many genes. We report that this inherent challenge can be overcome in model organisms. By pushing the advantages of inbred crossing to its practical limit in Saccharomyces cerevisiae, we improved the statistical resolution of linkage analysis to single nucleotides. This "super-resolution" approach allowed us to map 370 causal variants across 26 quantitative traits. Missense, synonymous, and cis-regulatory mutations collectively gave rise to phenotypic diversity, providing mechanistic insight into the basis of evolutionary divergence. Our data also systematically unmasked complex genetic architectures, revealing that multiple closely linked driver mutations frequently act on the same quantitative trait. Single-nucleotide mapping thus complements traditional deletion and overexpression screening paradigms and opens new frontiers in quantitative genetics.

Deciphering the Genic Basis of Yeast Fitness Variation by Simultaneous Forward and Reverse Genetics

The budding yeast Saccharomyces cerevisiae is the best studied eukaryote in molecular and cell biology, but its utility for understanding the genetic basis of phenotypic variation in natural populations is limited by inefficient association mapping due to strong and complex population structure. To overcome this challenge, we generated genome sequences for 85 strains and performed a comprehensive population genomic survey of a total of 190 diverse strains. We identified considerable variation in population structure among chromosomes and identified 181 genes that are absent from the reference genome. Many of these nonreference genes are expressed and we functionally confirmed that two of these genes confer increased resistance to antifungals. Next, we simultaneously measured the growth rates of over 4,500 laboratory strains, each of which lacks a nonessential gene, and 81 natural strains across multiple environments using unique DNA barcode present in each strain. By combining the genome-wide reverse genetic information gained from the gene deletion strains with a genome-wide association analysis from the natural strains, we identified genomic regions associated with fitness variation in natural populations. To experimentally validate a subset of these associations, we used reciprocal hemizygosity tests, finding that while the combined forward and reverse genetic approaches can identify a single causal gene, the phenotypic consequences of natural genetic variation often follow a complicated pattern. The resources and approach provided outline an efficient and reliable route to association mapping in yeast and significantly enhance its value as a model for understanding the genetic mechanisms underlying phenotypic variation and evolution in natural populations.

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.

Polymorphisms in the yeast galactose sensor underlie a natural continuum of nutrient-decision phenotypes

In nature, microbes often need to "decide" which of several available nutrients to utilize, a choice that depends on a cell's inherent preference and external nutrient levels. While natural environments can have mixtures of different nutrients, phenotypic variation in microbes' decisions of which nutrient to utilize is poorly studied. Here, we quantified differences in the concentration of glucose and galactose required to induce galactose-responsive (GAL) genes across 36 wild S. cerevisiae strains. Using bulk segregant analysis, we found that a locus containing the galactose sensor GAL3 was associated with differences in GAL signaling in eight different crosses. Using allele replacements, we confirmed that GAL3 is the major driver of GAL induction variation, and that GAL3 allelic variation alone can explain as much as 90% of the variation in GAL induction in a cross. The GAL3 variants we found modulate the diauxic lag, a selectable trait. These results suggest that ecological constraints on the galactose pathway may have led to variation in a single protein, allowing cells to quantitatively tune their response to nutrient changes in the environment.

Fitness Trade-Offs Lead to Suppressor Tolerance in Yeast

Genetic variation among individuals within a population provides the raw material for phenotypic diversity upon which natural selection operates. Some given variants can act on multiple standing genomic variations simultaneously and release previously inaccessible phenotypes, leading to increased adaptive potential upon challenging environments. Previously, we identified such a variant related to a tRNA nonsense suppressor in yeast. When introduced into other genetic backgrounds, the suppressor led to an increased population phenotypic variance on various culture conditions, conferring background and environment specific selective advantages. Nonetheless, most isolates are intolerant to the suppressor on rich media due to a severe fitness cost. Here, we found that the tolerance to suppressor is related to a surprising level of fitness outburst, showing a trade-off effect to accommodate the cost of carrying the suppressor. To dissect the genetic basis of such trade-offs, we crossed strains with contrasting tolerance levels on rich media, and analyzed the fitness distribution patterns in the offspring. Combining quantitative tetrad analysis and bulk segregant analysis, we identified two genes, namely MKT1 and RGA1, involved in suppressor tolerance. We showed that alleles from the tolerant parent for both genes conferred a significant gain of fitness, which increased the suppressor tolerance. Our results present a detailed dissection of suppressor tolerance in yeast and provide insights into the molecular basis of trade-offs between fitness and evolutionary potential.

The Hidden Complexity of Mendelian Traits across Natural Yeast Populations

Mendelian traits are considered to be at the lower end of the complexity spectrum of heritable phenotypes. However, more than a century after the rediscovery of Mendel's law, the global landscape of monogenic variants, as well as their effects and inheritance patterns within natural populations, is still not well understood. Using the yeast Saccharomyces cerevisiae, we performed a species-wide survey of Mendelian traits across a large population of isolates. We generated offspring from 41 unique parental pairs and analyzed 1,105 cross/trait combinations. We found that 8.9% of the cases were Mendelian. Further tracing of causal variants revealed background-specific expressivity and modified inheritances, gradually transitioning from Mendelian to complex traits in 30% of the cases. In fact, when taking into account the natural population diversity, the hidden complexity of traits could be substantial, confounding phenotypic predictability even for simple Mendelian traits.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.

Population genomics of yeasts: towards a comprehensive view across a broad evolutionary scale

With the advent of high-throughput technologies for sequencing, the complete description of the genetic variation that occurs in populations, also known as population genomics, is foreseeable but far from being reached. Explaining the forces that govern patterns of genetic variation is essential to elucidate the evolutionary history of species. Genetic variation results from a wide assortment of evolutionary forces, among which mutation, selection, recombination and drift play major roles in shaping genomes. In addition, exploring the genetic variation within a population also corresponds to the first step towards dissecting the genotype-phenotype relationship. In this context, yeast species are of particular interest because they represent a unique resource for studying the evolution of intraspecific genetic diversity in a phylum spanning a broad evolutionary scale. Here, we briefly review recent progress in yeast population genomics and provide some perspective on this rapidly evolving field. In fact, we truly believe that it is of interest to supplement comparative and early population genomic studies with the deep sequencing of more extensive sets of individuals from the same species. In parallel, it would be more than valuable to uncover the intraspecific variation of a large number of unexplored species, including those that are closely and more distantly related. Altogether, these data would enable substantially more powerful genomic scans for functional dissection.

Divergent MLS1 Promoters Lie on a Fitness Plateau for Gene Expression

Qualitative patterns of gene activation and repression are often conserved despite an abundance of quantitative variation in expression levels within and between species. A major challenge to interpreting patterns of expression divergence is knowing which changes in gene expression affect fitness. To characterize the fitness effects of gene expression divergence, we placed orthologous promoters from eight yeast species upstream of malate synthase (MLS1) in Saccharomyces cerevisiae. As expected, we found these promoters varied in their expression level under activated and repressed conditions as well as in their dynamic response following loss of glucose repression. Despite these differences, only a single promoter driving near basal levels of expression caused a detectable loss of fitness. We conclude that the MLS1 promoter lies on a fitness plateau whereby even large changes in gene expression can be tolerated without a substantial loss of fitness.

Population perspectives on functional genomic variation in yeast

Advances in high-throughput sequencing have facilitated large-scale surveys of genomic variation in the budding yeast,Saccharomyces cerevisiae These surveys have revealed extensive sequence variation between yeast strains. However, much less is known about how such variation influences the amount and nature of variation for functional genomic traits within and between yeast lineages. We review population-level studies of functional genomic variation, with a particular focus on how population functional genomic approaches can provide insights into both genome function and the evolutionary process. Although variation in functional genomics phenotypes is pervasive, our understanding of the consequences of this variation, either in physiological or evolutionary terms, is still rudimentary and thus motivates increased attention to appropriate null models. To date, much of the focus of population functional genomic studies has been on gene expression variation, but other functional genomic data types are just as likely to reveal important insights at the population level, suggesting a pressing need for more studies that go beyond transcription. Finally, we discuss how a population functional genomic perspective can be a powerful approach for developing a mechanistic understanding of the processes that link genomic variation to organismal phenotypes through gene networks.

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.

Temporal expression profiling identifies pathways mediating effect of causal variant on phenotype

Even with identification of multiple causal genetic variants for common human diseases, understanding the molecular processes mediating the causal variants’ effect on the disease remains a challenge. This understanding is crucial for the development of therapeutic strategies to prevent and treat disease. While static profiling of gene expression is primarily used to get insights into the biological bases of diseases, it makes differentiating the causative from the correlative effects difficult, as the dynamics of the underlying biological processes are not monitored. Using yeast as a model, we studied genome-wide gene expression dynamics in the presence of a causal variant as the sole genetic determinant, and performed allele-specific functional validation to delineate the causal effects of the genetic variant on the phenotype. Here, we characterized the precise genetic effects of a functional MKT1 allelic variant in sporulation efficiency variation. A mathematical model describing meiotic landmark events and conditional activation of MKT1 expression during sporulation specified an early meiotic role of this variant. By analyzing the early meiotic genome-wide transcriptional response, we demonstrate an MKT1-dependent role of novel modulators, namely, RTG1/3, regulators of mitochondrial retrograde signaling, and DAL82, regulator of nitrogen starvation, in additively effecting sporulation efficiency. In the presence of functional MKT1 allele, better respiration during early sporulation was observed, which was dependent on the mitochondrial retrograde regulator, RTG3. Furthermore, our approach showed that MKT1 contributes to sporulation independent of Puf3, an RNA-binding protein that steady-state transcription profiling studies have suggested to mediate MKT1-pleiotropic effects during mitotic growth. These results uncover interesting regulatory links between meiosis and mitochondrial retrograde signaling. In this study, we highlight the advantage of analyzing allele-specific transcriptional dynamics of mediating genes. Applications in higher eukaryotes can be valuable for inferring causal molecular pathways underlying complex dynamic processes, such as development, physiology and disease progression.

The causal path from a genetic variant to a complex phenotype such as disease progression is often not known. Studying gene expression variation is one approach to identify the mediating genes, however, it is difficult to distinguish causative from correlative genes. This becomes a challenge especially when studying developmental and physiological traits, since they involve dynamic processes contributing to the variation and only single static expression profiling is performed. As a proof of concept, we addressed this challenge here in yeast, by studying genome-wide gene expression in the presence of the causative polymorphism of MKT1 as the sole genetic variant, during the time phase when it contributes to sporulation efficiency variation. Our analysis during early sporulation identified mitochondrial retrograde signaling and nitrogen starvation as novel regulators, acting additively to regulate sporulation efficiency. Furthermore, we showed that PUF3, a known interactor of MKT1 had an independent role in sporulation. Our results highlight the role of differential mitochondrial signaling for efficient meiosis, providing insights into the factors regulating infertility. In addition, our study has implications for characterizing the molecular effects of causal genetic variants on dynamic biological processes during development and disease progression.

The fascinating and secret wild life of the budding yeast S. cerevisiae

The budding yeast Saccharomyces cerevisiae has been used in laboratory experiments for over a century and has been instrumental in understanding virtually every aspect of molecular biology and genetics. However, it wasn't until a decade ago that the scientific community started to realise how little was known about this yeast's ecology and natural history, and how this information was vitally important for interpreting its biology. Recent large-scale population genomics studies coupled with intensive field surveys have revealed a previously unappreciated wild lifestyle of S. cerevisiae outside the restrictions of human environments and laboratories. The recent discovery that Chinese isolates harbour almost twice as much genetic variation as isolates from the rest of the world combined suggests that Asia is the likely origin of the modern budding yeast.

DOI:http://dx.doi.org/10.7554/eLife.05835.001

Natural yeast promoter variants reveal epistasis in the generation of transcriptional-mediated noise and its potential benefit in stressful conditions

The increase in phenotypic variability through gene expression noise is proposed to be an evolutionary strategy in selective environments. Differences in promoter-mediated noise between Saccharomyces cerevisiae strains could have been selected for thanks to the benefit conferred by gene expression heterogeneity in the stressful conditions, for instance, those experienced by industrial strains. Here, we used a genome-wide approach to identify promoters conferring high noise levels in the industrial wine strain EC1118. Many promoters of genes related to environmental factors were identified, some of them containing genetic variations compared with their counterpart in the laboratory strain S288c. Each variant of eight promoters has been fused to yeast-Enhanced Green Fluorescent Protein and integrated in the genome of both strains. Some industrial variants conferred higher expression associated, as expected, with lower noise, but other variants either increased or decreased expression without modifying variability, so that they might exhibit different levels of transcriptional-mediated noise at equal mean. At different induction conditions giving similar expression for both variants of the CUP1 promoter, we indeed observed higher noise with the industrial variant. Nevertheless, this difference was only observed in the industrial strain, revealing epistasis in the generation of promoter-mediated noise. Moreover, the increased expression variability conferred by this natural yeast promoter variant provided a clear benefit in the face of an environmental stress. Thus, modulation of gene expression noise by a combination of promoter modifications and trans-influences might be a possible adaptation mechanism in yeast.

Differential regulation of antagonistic pleiotropy in synthetic and natural populations suggests its role in adaptation

Antagonistic pleiotropy (AP), the ability of a gene to show opposing effects in different phenotypes, has been identified in various life history traits and complex disorders, indicating its fundamental role in balancing fitness over the course of evolution. It is intuitive that natural selection might maintain AP to allow organisms phenotypic flexibility in different environments. However, despite several attempts, little evidence exists for its role in adaptation. We performed a meta-analysis in yeast to identify the genetic basis of AP in bi-parental segregants, natural isolates, and a laboratory strain genome-wide deletion collection, by comparing growth in favorable and stress conditions. We found that whereas AP was abundant in the synthetic populations, it was absent in the natural isolates. This finding indicated resolution of trade-offs, i.e., mitigation of trade-offs over evolutionary history, probably through accumulation of compensatory mutations. In the deletion collection, organizational genes showed AP, suggesting ancient resolutions of trade-offs in the basic cellular pathways. We find abundant AP in the segregants, greater than estimated in the deletion collection or observed in previous studies, with IRA2, a negative regulator of the Ras/PKA signaling pathway, showing trade-offs across diverse environments. Additionally, IRA2 and several other Ras/PKA pathway genes showed balancing selection in isolates of S. cerevisiae and S. paradoxus, indicating that multiple alleles maintain AP in this pathway in natural populations. We propose that during AP resolution, retaining the ability to vary signaling pathways such as Ras/PKA, may provide organisms with phenotypic flexibility. However, with increasing organismal complexity AP resolution may become difficult. A partial resolution of AP could manifest as complex human diseases, and the inability to resolve AP may play a role in speciation. Our findings suggest that testing a universal phenomenon like AP across multiple experimental systems may elucidate mechanisms underlying its regulation and evolution.