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

Horizontal transfer and recombination fuel Ty4 retrotransposon evolution in Saccharomyces

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

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Contributions of mutation and selection to regulatory variation: lessons from the Saccharomyces cerevisiae TDH3 gene

Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting either cis- or trans-regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of the TDH3 gene in Saccharomyces cerevisiae and comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties of cis- and trans-regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Network Topology Can Explain Differences in Pleiotropy Between Cis- and Trans-regulatory Mutations

A mutation's degree of pleiotropy (i.e., the number of traits it alters) is predicted to impact the probability of the mutation being detrimental to fitness. For mutations that impact gene expression, mutations acting in cis have been hypothesized to generally be less pleiotropic than mutations affecting the same gene's expression in trans, suggesting that cis-regulatory mutations should be less deleterious and more likely to fix over evolutionary time. Here, we use expression and fitness data from Saccharomyces cerevisiae gene deletion strains to test these hypotheses. By treating deletion of each gene as a cis-regulatory mutation affecting its own expression and deletions of other genes affecting expression of this focal gene as trans-regulatory mutations, we find that cis-acting mutations do indeed tend to be less pleiotropic than trans-acting mutations affecting expression of the same gene. This pattern was observed for the vast majority of genes in the data set and could be explained by the topology of the regulatory network controlling gene expression. Comparing the fitness of cis- and trans-acting mutations affecting expression of the same gene also confirmed that trans-acting deletions tend to be more deleterious. These findings provide strong support for pleiotropy playing a role in the preferential fixation of cis-regulatory alleles over evolutionary time.

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.

Fluorescence lifetime imaging of NAD(P)H upon oxidative stress in Kluyveromyces marxianus

K. marxianus is a promising cell factory for producing heterologous proteins. Oxidative stresses were raised during overexpression of heterologous proteins, leading to the shift of the redox state. How to measure the redox state of live K. marxianus cells without perturbing their growth remains a big challenge. Here, a fluorescence lifetime imaging (FLIM)-based method was developed in live K. marxianus cells. During the early exponential growth, K. marxianus cells exhibited an increased mean fluorescence lifetime (τ-mean) of NAD(P)H compared with Saccharomyces cerevisiae cells, which was consistent with the preference for respiration in K. marxianus cells and that for fermentation in S. cerevisiae cells. Upon oxidative stresses induced by high temperature or H2O2, K. marxianus cells exhibited an increased τ-mean in company with decreased intracellular NAD(P)H/NAD(P)+, suggesting a correlation between an increased τ-mean and a more oxidized redox state. The relationship between τ-mean and the expression level of a heterologous protein was investigated. There was no difference between the τ-means of K. marxianus strains which were not producing a heterologous protein. The τ-mean of a strain yielding a high level of a heterologous protein was higher than that of a low-yielding strain. The results suggested the potential application of FLIM in the non-invasive screen of high-yielding cells.

Pleiotropic effects of trans-regulatory mutations on fitness and gene expression

Variation in gene expression arises from cis- and trans-regulatory mutations, which contribute differentially to expression divergence. We compare the impacts on gene expression and fitness resulting from cis- and trans-regulatory mutations in Saccharomyces cerevisiae, with a focus on the TDH3 gene. We use the effects of cis-regulatory mutations to infer effects of trans-regulatory mutations attributable to impacts beyond the focal gene, revealing a distribution of pleiotropic effects. Cis- and trans-regulatory mutations had different effects on gene expression with pleiotropic effects of trans-regulatory mutants affecting expression of genes both in parallel to and downstream of the focal gene. The more widespread and deleterious effects of trans-regulatory mutations we observed are consistent with their decreasing relative contribution to expression differences over evolutionary time.

Epistatic drift causes gradual decay of predictability in protein evolution

Epistatic interactions can make the outcomes of evolution unpredictable, but no comprehensive data are available on the extent and temporal dynamics of changes in the effects of mutations as protein sequences evolve. Here, we use phylogenetic deep mutational scanning to measure the functional effect of every possible amino acid mutation in a series of ancestral and extant steroid receptor DNA binding domains. Across 700 million years of evolution, epistatic interactions caused the effects of most mutations to become decorrelated from their initial effects and their windows of evolutionary accessibility to open and close transiently. Most effects changed gradually and without bias at rates that were largely constant across time, indicating a neutral process caused by many weak epistatic interactions. Our findings show that protein sequences drift inexorably into contingency and unpredictability, but that the process is statistically predictable, given sufficient phylogenetic and experimental data.

Transposon insertional mutagenesis of diverse yeast strains suggests coordinated gene essentiality polymorphisms

Due to epistasis, the same mutation can have drastically different phenotypic consequences in different individuals. This phenomenon is pertinent to precision medicine as well as antimicrobial drug development, but its general characteristics are largely unknown. We approach this question by genome-wide assessment of gene essentiality polymorphism in 16 Saccharomyces cerevisiae strains using transposon insertional mutagenesis. Essentiality polymorphism is observed for 9.8% of genes, most of which have had repeated essentiality switches in evolution. Genes exhibiting essentiality polymorphism lean toward having intermediate numbers of genetic and protein interactions. Gene essentiality changes tend to occur concordantly among components of the same protein complex or metabolic pathway and among a group of over 100 mitochondrial proteins, revealing molecular machines or functional modules as units of gene essentiality variation. Most essential genes tolerate transposon insertions consistently among strains in one or more coding segments, delineating nonessential regions within essential genes.

The rate and molecular spectrum of mutation are selectively maintained in yeast

What determines the rate (μ) and molecular spectrum of mutation is a fundamental question. The prevailing hypothesis asserts that natural selection against deleterious mutations has pushed μ to the minimum achievable in the presence of genetic drift, or the drift barrier. Here we show that, contrasting this hypothesis, μ substantially exceeds the drift barrier in diverse organisms. Random mutation accumulation (MA) in yeast frequently reduces μ, and deleting the newly discovered mutator gene PSP2 nearly halves μ. These results, along with a comparison between the MA and natural yeast strains, demonstrate that μ is maintained above the drift barrier by stabilizing selection. Similar comparisons show that the mutation spectrum such as the universal AT mutational bias is not intrinsic but has been selectively preserved. These findings blur the separation of mutation from selection as distinct evolutionary forces but open the door to alleviating mutagenesis in various organisms by genome editing.

How natural selection shapes the rate and molecular spectrum of mutations is debated. Yeast mutation accumulation experiments identify a gene promoting mutagenesis and show stabilizing selection maintaining the mutation rate above the drift barrier. Selection also preserves the mutation spectrum.

Asexual Experimental Evolution of Yeast Does Not Curtail Transposable Elements

Compared with asexual reproduction, sex facilitates the transmission of transposable elements (TEs) from one genome to another, but boosts the efficacy of selection against deleterious TEs. Thus, theoretically, it is unclear whether sex has a positive net effect on TE’s proliferation. An empirical study concluded that sex is at the root of TE’s evolutionary success because the yeast TE load was found to decrease rapidly in approximately 1,000 generations of asexual but not sexual experimental evolution. However, this finding contradicts the maintenance of TEs in natural yeast populations where sexual reproduction occurs extremely infrequently. Here, we show that the purported TE load reduction during asexual experimental evolution is likely an artifact of low genomic sequencing coverages. We observe stable TE loads in both sexual and asexual experimental evolution from multiple yeast data sets with sufficient coverages. To understand the evolutionary dynamics of yeast TEs, we turn to asexual mutation accumulation lines that have been under virtually no selection. We find that both TE transposition and excision rates per generation, but not their difference, tend to be higher in environments where yeast grows more slowly. However, the transposition rate is not significantly higher than the excision rate and the variance of the TE number among natural strains is close to its neutral expectation, suggesting that selection against TEs is at best weak in yeast. We conclude that the yeast TE load is maintained largely by a transposition–excision balance and that the influence of sex remains unclear.

Competition experiments in a soil microcosm reveal the impact of genetic and biotic factors on natural yeast populations

The ability to measure microbial fitness directly in natural conditions and in interaction with other microbes is a challenge that needs to be overcome if we want to gain a better understanding of microbial fitness determinants in nature. Here we investigate the influence of the natural microbial community on the relative fitness of the North American populations SpB, SpC and SpC* of the wild yeast Saccharomyces paradoxus using DNA barcodes and a soil microcosm derived from soil associated with oak trees. We find that variation in fitness among these genetically distinct groups is influenced by the microbial community. Altering the microbial community load and diversity with an irradiation treatment significantly diminishes the magnitude of fitness differences among populations. Our findings suggest that microbial interactions could affect the evolution of yeast lineages in nature by modulating variation in fitness.

Antagonistic pleiotropy conceals molecular adaptations in changing environments

The importance of positive selection in molecular evolution is debated. Evolution experiments under invariant laboratory conditions typically show a higher rate of nonsynonymous nucleotide changes than the rate of synonymous changes, demonstrating prevalent molecular adaptations. Natural evolution inferred from genomic comparisons, however, almost always exhibits the opposite pattern even among closely related conspecifics, which is indicative of a paucity of positive selection. Here we hypothesize that this apparent contradiction is at least in part attributable to ubiquitous and frequent environmental changes in nature, causing nonsynonymous mutations that are beneficial at one time to become deleterious soon after because of antagonistic pleiotropy and hindering their fixations relative to synonymous mutations despite continued population adaptations. To test this hypothesis, we performed yeast evolution experiments in changing and corresponding constant environments, followed by genome sequencing of the evolving populations. We observed a lower nonsynonymous to synonymous rate ratio in antagonistic changing environments than in the corresponding constant environments, and the population dynamics of mutations supports our hypothesis. These findings and the accompanying population genetic simulations suggest that molecular adaptation is consistently underestimated in nature due to the antagonistic fitness effects of mutations in changing environments.

Evaluation of Saccharomyces cerevisiae Wine Yeast Competitive Fitness in Enologically Relevant Environments by Barcode Sequencing

When a wine yeast is inoculated into grape juice the potential variation in juice composition that confronts it is huge. Assessing the performance characteristics of the many commercially available wine yeasts in the many possible grape juice compositions is a daunting task. To this end we have developed a barcoded Saccharomyces cerevisiae wine yeast collection to facilitate the task of performance assessment that will contribute to a broader understanding of genotype-phenotype relations. Barcode sequencing of mixed populations is used to monitor strain abundance in different grape juices and grape juice-like environments. Choice of DNA extraction method is shown to affect strain-specific barcode count in this highly related set of S. cerevisiae strains; however, the analytical approach is shown to be robust toward strain dependent variation in DNA extraction efficiency. Of the 38 unique compositional variables assessed, resistance to copper and SO₂ are found to be dominant discriminatory factors in wine yeast performance. Finally, a comparison of competitive fitness profile with performance in single inoculum fermentations reveal strain dependent correspondence of yeast performance using these two different approaches.

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.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

Empirical measures of mutational effects define neutral models of regulatory evolution in Saccharomyces cerevisiae

Understanding how phenotypes evolve requires disentangling the effects of mutation generating new variation from the effects of selection filtering it. Tests for selection frequently assume that mutation introduces phenotypic variation symmetrically around the population mean, yet few studies have tested this assumption by deeply sampling the distributions of mutational effects for particular traits. Here, we examine distributions of mutational effects for gene expression in the budding yeast Saccharomyces cerevisiae by measuring the effects of thousands of point mutations introduced randomly throughout the genome. We find that the distributions of mutational effects differ for the 10 genes surveyed and are inconsistent with normality. For example, all 10 distributions of mutational effects included more mutations with large effects than expected for normally distributed phenotypes. In addition, some genes also showed asymmetries in their distribution of mutational effects, with new mutations more likely to increase than decrease the gene's expression or vice versa. Neutral models of regulatory evolution that take these empirically determined distributions into account suggest that neutral processes may explain more expression variation within natural populations than currently appreciated.

Compensatory trans-regulatory alleles minimizing variation in TDH3 expression are common within Saccharomyces cerevisiae

Heritable variation in gene expression is common within species. Much of this variation is due to genetic differences outside of the gene with altered expression and is trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture and evolution of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including the variability of trans-regulatory variation within a species, how selection affects trans-regulatory variation, and how trans-regulatory variants are distributed throughout the genome and within a species. To address these questions, we isolated and measured trans-regulatory differences affecting TDH3 promoter activity among 56 strains of Saccharomyces cerevisiae, finding that trans-regulatory backgrounds varied approximately twofold in their effects on TDH3 promoter activity. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that despite this variability in the effects of trans-regulatory backgrounds, stabilizing selection has constrained trans-regulatory differences within this species. Using a powerful quantitative trait locus mapping method, we identified ∼100 trans-acting expression quantitative trait locus in each of three crosses to a common reference strain, indicating that regulatory variation is more polygenic than previous studies have suggested. Loci altering expression were located throughout the genome, and many loci were strain specific. This distribution and prevalence of alleles is consistent with recent theories about the genetic architecture of complex traits. In all mapping experiments, the nonreference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. This variation may provide the raw material for compensatory evolution and larger scale regulatory rewiring observed in developmental systems drift among species.

Overdosage of Balanced Protein Complexes Reduces Proliferation Rate in Aneuploid Cells

Cells with complex aneuploidies display a wide range of phenotypic abnormalities. However, the molecular basis for this has been mainly studied in trisomic (2n + 1) and disomic (n + 1) cells. To determine how karyotype affects proliferation in cells with complex aneuploidies, we generated 92 2n + x yeast strains in which each diploid cell has between 3 and 12 extra chromosomes. Genome-wide and, for individual protein complexes, proliferation defects are caused by the presence of protein complexes in which all subunits are balanced at the 3-copy level. Proteomics revealed that over 50% of 3-copy members of imbalanced complexes were expressed at only 2n protein levels, whereas members of complexes in which all subunits are stoichiometrically balanced at 3 copies per cell had 3n protein levels. We validated this finding using orthogonal datasets from yeast and from human cancers. Taken together, our study provides an explanation of how aneuploidy affects phenotype.

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.

A collection of barcoded natural isolates of Saccharomyces paradoxus to study microbial evolutionary ecology

While the use of barcoded collections of laboratory microorganisms and the development of barcode-based cell tracking are rapidly developing in genetics and genomics research, tools to track natural populations are still lacking. The yeast Saccharomyces paradoxus is an emergent microbial model in ecology and evolution. More than five allopatric and sympatric lineages have been identified and hundreds of strains have been isolated for this species, allowing to assess the impact of natural diversity on complex traits. We constructed a collection of 550 barcoded and traceable strains of S. paradoxus, including all three North American lineages SpB, SpC, and SpC*. These strains are diploid, many have their genome fully sequenced and are barcoded with a unique 20 bp sequence that allows their identification and quantification. This yeast collection is functional for competitive experiments in pools as the barcodes allow to measure each lineage's and individual strains' fitness in common conditions. We used this tool to demonstrate that in the tested conditions, there are extensive genotype-by-environment interactions for fitness among S. paradoxus strains, which reveals complex evolutionary potential in variable environments. This barcoded collection provides a valuable resource for ecological genomics studies that will allow gaining a better understanding of S. paradoxus evolution and fitness-related traits.

The optimal mating distance resulting from heterosis and genetic incompatibility

Theory predicts that the fitness of an individual is maximized when the genetic distance between its parents (i.e., mating distance) is neither too small nor too large. However, decades of research have generally failed to validate this prediction or identify the optimal mating distance (OMD). Respectively analyzing large numbers of crosses of fungal, plant, and animal model organisms, we indeed find the hybrid phenotypic value a humped quadratic polynomial function of the mating distance for the vast majority of fitness-related traits examined, with different traits of the same species exhibiting similar OMDs. OMDs are generally slightly greater than the nucleotide diversities of the species concerned but smaller than the observed maximal intraspecific genetic distances. Hence, the benefit of heterosis is at least partially offset by the harm of genetic incompatibility even within species. These results have multiple theoretical and practical implications for speciation, conservation, and agriculture.

Genome-Wide Screen for Saccharomyces cerevisiae Genes Contributing to Opportunistic Pathogenicity in an Invertebrate Model Host

Environmental opportunistic pathogens can exploit vulnerable hosts through expression of traits selected for in their natural environments. Pathogenicity is itself a complicated trait underpinned by multiple complex traits, such as thermotolerance, morphology, and stress response. The baker's yeast, Saccharomyces cerevisiae, is a species with broad environmental tolerance that has been increasingly reported as an opportunistic pathogen of humans. Here we leveraged the genetic resources available in yeast and a model insect species, the greater waxmoth Galleria mellonella, to provide a genome-wide analysis of pathogenicity factors. Using serial passaging experiments of genetically marked wild-type strains, a hybrid strain was identified as the most fit genotype across all replicates. To dissect the genetic basis for pathogenicity in the hybrid isolate, bulk segregant analysis was performed which revealed eight quantitative trait loci significantly differing between the two bulks with alleles from both parents contributing to pathogenicity. A second passaging experiment with a library of deletion mutants for most yeast genes identified a large number of mutations whose relative fitness differed in vivovs.in vitro, including mutations in genes controlling cell wall integrity, mitochondrial function, and tyrosine metabolism. Yeast is presumably subjected to a massive assault by the innate insect immune system that leads to melanization of the host and to a large bottleneck in yeast population size. Our data support that resistance to the innate immune response of the insect is key to survival in the host and identifies shared genetic mechanisms between S. cerevisiae and other opportunistic fungal pathogens.

Effects of mutation and selection on plasticity of a promoter activity in Saccharomyces cerevisiae

Phenotypic plasticity is an evolvable property of biological systems that can arise from environment-specific regulation of gene expression. To better understand the evolutionary and molecular mechanisms that give rise to plasticity in gene expression, we quantified the effects of 235 single-nucleotide mutations in the Saccharomyces cerevisiae TDH3 promoter (P_TDH3 ) on the activity of this promoter in media containing glucose, galactose, or glycerol as a carbon source. We found that the distributions of mutational effects differed among environments because many mutations altered the plastic response exhibited by the wild-type allele. Comparing the effects of these mutations with the effects of 30 P_TDH3 polymorphisms on expression plasticity in the same environments provided evidence of natural selection acting to prevent the plastic response in P_TDH3 activity between glucose and galactose from becoming larger. The largest changes in expression plasticity were observed between fermentable (glucose or galactose) and nonfermentable (glycerol) carbon sources and were caused by mutations located in the RAP1 and GCR1 transcription factor binding sites. Mutations altered expression plasticity most frequently between the two fermentable environments, with mutations causing significant changes in plasticity between glucose and galactose distributed throughout the promoter, suggesting they might affect chromatin structure. Taken together, these results provide insight into the molecular mechanisms underlying gene-by-environment interactions affecting gene expression as well as the evolutionary dynamics affecting natural variation in plasticity of gene expression.

Most m6A RNA Modifications in Protein-Coding Regions Are Evolutionarily Unconserved and Likely Nonfunctional

Methylation of the adenosine base at the nitrogen-6 position (m6A) is the most prevalent internal posttranscriptional modification of mRNAs in many eukaryotes. Despite the rapid progress in the transcriptome-wide mapping of m6As, identification of proteins responsible for writing, reading, and erasing m6As, and elucidation of m6A functions in splicing, RNA stability, translation, and other processes, it is unknown whether most observed m6A modifications are functional. To address this question, we respectively analyze the evolutionary conservation of yeast and human m6As in protein-coding regions. Relative to comparable unmethylated As, m6As are overall no more conserved in yeasts and only slightly more conserved in mammals. Furthermore, yeast m6As and comparable unmethylated As have no significant difference in single nucleotide polymorphism (SNP) density or SNP site frequency spectrum. The same is true in human. The methylation status of a gene, not necessarily the specific sites methylated in the gene, is subject to purifying selection for no more than ∼20% of m6A-modified genes. These observations suggest that most m6A modifications in protein-coding regions are nonfunctional and nonadaptive, probably resulting from off-target activities of m6A methyltransferases. In addition, our reanalysis invalidates the recent claim of positive selection for newly acquired m6A modifications in human evolution. Regarding the small number of evolutionarily conserved m6As, evidence suggests that a large proportion of them are likely functional; they should be prioritized in future functional characterizations of m6As. Together, these findings have important implications for understanding the biological significance of m6A and other posttranscriptional modifications.

Testing the neutral hypothesis of phenotypic evolution

Although evolution by natural selection is widely regarded as the most important principle of biology, it is unknown whether phenotypic variations within and between species are mostly adaptive or neutral due to the lack of relevant studies of large, unbiased samples of phenotypic traits. Here, we examine 210 yeast morphological traits chosen because of experimental feasibility irrespective of their potential adaptive values. Our analysis is based on the premise that, under neutrality, the rate of phenotypic evolution measured in the unit of mutational size declines as the trait becomes more important to fitness, analogous to the neutral paradigm that functional genes evolve more slowly than functionless pseudogenes. However, we find faster evolution of more important morphological traits within and between species, rejecting the neutral hypothesis. By contrast, an analysis of 3,466 gene expression traits fails to refute neutrality. Thus, at least in yeast, morphological evolution appears largely adaptive, but the same may not apply to other classes of phenotypes. Our neutrality test is applicable to other species, especially genetic model organisms, for which estimations of mutational size and trait importance are relatively straightforward.

Intra and Interspecific Variations of Gene Expression Levels in Yeast Are Largely Neutral: (Nei Lecture, SMBE 2016, Gold Coast)

It is commonly, although not universally, accepted that most intra and interspecific genome sequence variations are more or less neutral, whereas a large fraction of organism-level phenotypic variations are adaptive. Gene expression levels are molecular phenotypes that bridge the gap between genotypes and corresponding organism-level phenotypes. Yet, it is unknown whether natural variations in gene expression levels are mostly neutral or adaptive. Here we address this fundamental question by genome-wide profiling and comparison of gene expression levels in nine yeast strains belonging to three closely related Saccharomyces species and originating from five different ecological environments. We find that the transcriptome-based clustering of the nine strains approximates the genome sequence-based phylogeny irrespective of their ecological environments. Remarkably, only ∼0.5% of genes exhibit similar expression levels among strains from a common ecological environment, no greater than that among strains with comparable phylogenetic relationships but different environments. These and other observations strongly suggest that most intra and interspecific variations in yeast gene expression levels result from the accumulation of random mutations rather than environmental adaptations. This finding has profound implications for understanding the driving force of gene expression evolution, genetic basis of phenotypic adaptation, and general role of stochasticity in evolution.

Evolutionary Dynamics of Regulatory Changes Underlying Gene Expression Divergence among Saccharomyces Species

Heritable changes in gene expression are important contributors to phenotypic differences within and between species and are caused by mutations in cis-regulatory elements and trans-regulatory factors. Although previous work has suggested that cis-regulatory differences preferentially accumulate with time, technical restrictions to closely related species and limited comparisons have made this observation difficult to test. To address this problem, we used allele-specific RNA-seq data from Saccharomyces species and hybrids to expand both the evolutionary timescale and number of species in which the evolution of regulatory divergence has been investigated. We find that as sequence divergence increases, cis-regulatory differences do indeed become the dominant type of regulatory difference between species, ultimately becoming a better predictor of expression divergence than trans-regulatory divergence. When both cis- and trans-regulatory differences accumulate for the same gene, they more often have effects in opposite directions than in the same direction, indicating widespread compensatory changes underlying the evolution of gene expression. The frequency of compensatory changes within and between species and the magnitude of effect for the underlying cis- and trans-regulatory differences suggests that compensatory changes accumulate primarily due to selection against divergence in gene expression as a result of weak stabilizing selection on gene expression levels. These results show that cis-regulatory differences and compensatory changes in regulation play increasingly important roles in the evolution of gene expression as time increases.