Evolution in Saccharomyces cerevisiae: Identification of Mutations Increasing Fitness in Laboratory Populations

Template switching during DNA replication is a prevalent source of adaptive gene amplification

Copy number variants (CNVs) are an important source of genetic variation underlying rapid adaptation and genome evolution. Whereas point mutation rates vary with genomic location and local DNA features, the role of genome architecture in the formation and evolutionary dynamics of CNVs is poorly understood. Previously, we found the GAP1 gene in Saccharomyces cerevisiae undergoes frequent amplification and selection in glutamine-limitation. The gene is flanked by two long terminal repeats (LTRs) and proximate to an origin of DNA replication (autonomously replicating sequence, ARS), which likely promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we evolved engineered strains lacking either the adjacent LTRs, ARS, or all elements in glutamine-limited chemostats. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. Removal of local DNA elements significantly impacts the fitness effect of GAP1 CNVs and the rate of adaptation. In 177 CNV lineages, across all four strains, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, distal ones mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination can mediate gene amplification following de novo retrotransposon events. Our study reveals that template switching during DNA replication is a prevalent source of adaptive CNVs.

Template switching during DNA replication is a prevalent source of adaptive gene amplification

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.

Segregating Complete Tf2 Elements Are Largely Neutral in Fission Yeast

Transposable elements (TEs) comprise a large proportion of the eukaryote genomes. Yet it remains poorly understood how TEs influence the fitness of the hosts carrying them. Here, we empirically test the impact of TEs on the host fitness in the fission yeast Schizosaccharomyces pombe. We find that two families of TEs (Tf1 and Tf2 elements), both of which belong to long terminal repeat retrotransposons, are highly polymorphic among individual S. pombe strains. Only 13 complete Tf2 elements are identified in S. pombe laboratory strain 972. These 13 Tf2 elements integrated into host genomes in very recent time and are segregating within the S. pombe population. Through knocking out each of the 13 Tf2 elements in S. pombe strain 972, we find Tf2 knockout does not affect the host fitness, and Tf2 elements do not alter the expression of nearby genes. Challenged by diverse forms of stress, the Tf2 knockout strains do not exhibit different growth rates from wild-type strain. Together, we conclude that segregating complete Tf2 elements insertions are largely neutral to host fitness in the fission yeast. Our study provides genome-wide empirical support for the selfish nature of TEs in fission yeast.

The functional basis of adaptive evolution in chemostats

Two of the central problems in biology are determining the molecular basis of adaptive evolution and understanding how cells regulate their growth. The chemostat is a device for culturing cells that provides great utility in tackling both of these problems: it enables precise control of the selective pressure under which organisms evolve and it facilitates experimental control of cell growth rate. The aim of this review is to synthesize results from studies of the functional basis of adaptive evolution in long-term chemostat selections using Escherichia coli and Saccharomyces cerevisiae. We describe the principle of the chemostat, provide a summary of studies of experimental evolution in chemostats, and use these studies to assess our current understanding of selection in the chemostat. Functional studies of adaptive evolution in chemostats provide a unique means of interrogating the genetic networks that control cell growth, which complements functional genomic approaches and quantitative trait loci (QTL) mapping in natural populations. An integrated approach to the study of adaptive evolution that accounts for both molecular function and evolutionary processes is critical to advancing our understanding of evolution. By renewing efforts to integrate these two research programs, experimental evolution in chemostats is ideally suited to extending the functional synthesis to the study of genetic networks.

Gene dispensability

Genome-wide mutagenesis studies indicate that up to about 90% of genes in bacteria and 80% in eukaryotes can be inactivated individually leaving an organism viable, often seemingly unaffected. Several strategies are used to learn what these apparently dispensable genes contribute to fitness. Assays of growth under hundreds of physical and chemical stresses are among the most effective experimental approaches. Comparative studies of genomic DNA sequences continue to be valuable in discriminating between the core bacterial genome and the more variable niche-specific genes. The concept of the core genome appears currently unfeasible for eukaryotes but progress has been made in understanding why they contain numerous gene duplicates.

Adaptation of Saccharomyces cerevisiae to saline stress through laboratory evolution

Most laboratory evolution studies that characterize evolutionary adaptation genomically focus on genetically simple traits that can be altered by one or few mutations. Such traits are important, but they are few compared with complex, polygenic traits influenced by many genes. We know much less about complex traits, and about the changes that occur in the genome and in gene expression during their evolutionary adaptation. Salt stress tolerance is such a trait. It is especially attractive for evolutionary studies, because the physiological response to salt stress is well-characterized on the molecular and transcriptome level. This provides a unique opportunity to compare evolutionary adaptation and physiological adaptation to salt stress. The yeast Saccharomyces cerevisiae is a good model system to study salt stress tolerance, because it contains several highly conserved pathways that mediate the salt stress response. We evolved three replicate lines of yeast under continuous salt (NaCl) stress for 300 generations. All three lines evolved faster growth rate in high salt conditions than their ancestor. In these lines, we studied gene expression changes through microarray analysis and genetic changes through next generation population sequencing. We found two principal kinds of gene expression changes, changes in basal expression (82 genes) and changes in regulation (62 genes). The genes that change their expression involve several well-known physiological stress-response genes, including CTT1, MSN4 and HLR1. Next generation sequencing revealed only one high-frequency single-nucleotide change, in the gene MOT2, that caused increased fitness when introduced into the ancestral strain. Analysis of DNA content per cell revealed ploidy increases in all the three lines. Our observations suggest that evolutionary adaptation of yeast to salt stress is associated with genome size increase and modest expression changes in several genes.

Experimental evolution in yeast: a practical guide

Experimental evolution refers to a broad range of studies in which selection pressures are applied to populations. In some applications, particular traits are desired, while in others the subject of study is the mechanisms of evolution or the different modes of behavior between systems. This chapter will explore the range of studies falling under the experimental evolution umbrella, and their relative merits for different types of applications. Practical aspects of experimental evolution will also be discussed, including commercial suppliers, analysis methods, and best laboratory practices.

Low frequency of mutations with strongly deleterious but nonlethal fitness effects

Most experimentally detectable effects of mutations in cellular organisms are either lethal or mildly deleterious. A possible explanation for the paucity of strongly detrimental but nonlethal mutations is that the processes constituting cellular metabolism are either essential or largely redundant. Alternatively, the partition between lethal and inconspicuous mutations exists within important biological processes. To test this, we measured maximum growth rates of yeast strains each carrying the deletion of a single gene in one of 38 protein complexes. We also used relevant data from previous high-throughput phenotypic studies of the yeast gene-deletion collection. The complexes typified well-defined sets of genes engaged in a common process. Within virtually all essential complexes there were two clear modes of phenotypic effects, that is the cessation of growth or slowdown of growth by a few percent. This uniformity is striking given that complexes differ extensively in function, size, and proportion of essential proteins. The pattern of bimodality is observed both under optimal and suboptimal environmental conditions. The generic paucity of strong effects and abundance of small ones relates to the feasibility of analyses of quantitative traits and epidemiological surveys, irrespective of the particular element of metabolism under study.

Molecular characterization of clonal interference during adaptive evolution in asexual populations of Saccharomyces cerevisiae

The classical model of adaptive evolution in an asexual population postulates that each adaptive clone is derived from the one preceding it1. However, experimental evidence suggests more complex dynamics2-5 with theory predicting the fixation probability of a beneficial mutation as dependent on the mutation rate, population size, and the mutation's selection coefficient6. Clonal interference has been demonstrated in viruses7 and bacteria8, but has not been demonstrated in a eukaryote and a detailed molecular characterization is lacking. Here we use different fluorescent markers to visualize the dynamics of asexually evolving yeast populations. For each adaptive clone within one of our evolving populations, we have identified the underlying mutations, monitored their population frequencies and used microarrays to characterize changes in the transcriptome. These data provide the most detailed molecular characterization of an experimental evolution to date, and provide direct experimental evidence supporting both the clonal interference and the multiple mutation models.

The repertoire and dynamics of evolutionary adaptations to controlled nutrient-limited environments in yeast

The experimental evolution of laboratory populations of microbes provides an opportunity to observe the evolutionary dynamics of adaptation in real time. Until very recently, however, such studies have been limited by our inability to systematically find mutations in evolved organisms. We overcome this limitation by using a variety of DNA microarray-based techniques to characterize genetic changes—including point mutations, structural changes, and insertion variation—that resulted from the experimental adaptation of 24 haploid and diploid cultures of Saccharomyces cerevisiae to growth in either glucose, sulfate, or phosphate-limited chemostats for ∼200 generations. We identified frequent genomic amplifications and rearrangements as well as novel retrotransposition events associated with adaptation. Global nucleotide variation detection in ten clonal isolates identified 32 point mutations. On the basis of mutation frequencies, we infer that these mutations and the subsequent dynamics of adaptation are determined by the batch phase of growth prior to initiation of the continuous phase in the chemostat. We relate these genotypic changes to phenotypic outcomes, namely global patterns of gene expression, and to increases in fitness by 5–50%. We found that the spectrum of available mutations in glucose- or phosphate-limited environments combined with the batch phase population dynamics early in our experiments allowed several distinct genotypic and phenotypic evolutionary pathways in response to these nutrient limitations. By contrast, sulfate-limited populations were much more constrained in both genotypic and phenotypic outcomes. Thus, the reproducibility of evolution varies with specific selective pressures, reflecting the constraints inherent in the system-level organization of metabolic processes in the cell. We were able to relate some of the observed adaptive mutations (e.g., transporter gene amplifications) to known features of the relevant metabolic pathways, but many of the mutations pointed to genes not previously associated with the relevant physiology. Thus, in addition to answering basic mechanistic questions about evolutionary mechanisms, our work suggests that experimental evolution can also shed light on the function and regulation of individual metabolic pathways.

Adaptive evolution is a central biological process that underlies diverse phenomena from the acquisition of antibiotic resistance by microbes to the evolution of niche specialization. Two unresolved questions regarding adaptive evolution are what types of genomic variation are associated with adaptation and how repeatable is the process. We evolved yeast populations for more than 200 generations in nutrient-limited chemostats. We find that the phenotype of adapted individuals, as measured using global gene expression, is much less variable in clones adapted to sulfate limitation than either glucose or phosphate limitation. We comprehensively analyzed the genomes of adapted clones and found that those adapted to sulfate limitation almost invariably carry amplifications of the gene encoding a sulfur transporter, but the mutations in individuals adapted to glucose and phosphate limitation are much more diverse. This parallelism holds true at the level of single-nucleotide mutations. Although there may be other paths to adapt to sulfate limitation, one path confers a much greater advantage than all others so it dominates. By contrast, there are a number of ways to adapt to glucose and phosphate limitation that confer similar advantages. We conclude that the reproducibility of evolution depends on the specific selective pressure experienced by the organism.

Comparing whole genomes using DNA microarrays

The rapid accumulation of complete genomic sequences offers the opportunity to carry out an analysis of inter- and intra-individual genome variation within a species on a routine basis. Sequencing whole genomes requires resources that are currently beyond those of a single laboratory and therefore it is not a practical approach for resequencing hundreds of individual genomes. DNA microarrays present an alternative way to study differences between closely related genomes. Advances in microarray-based approaches have enabled the main forms of genomic variation (amplifications, deletions, insertions, rearrangements and base-pair changes) to be detected using techniques that are readily performed in individual laboratories using simple experimental approaches.

Genetic characterization of commercial Saccharomyces cerevisiae isolates recovered from vineyard environments

One hundred isolates of the commercial Saccharomyces cerevisiae strain Zymaflore VL1 were recovered from spontaneous fermentations carried out with grapes collected from vineyards located close to wineries in the Vinho Verde wine region of Portugal. Isolates were differentiated based on their mitochondrial DNA restriction patterns and the evaluation of genetic polymorphisms was carried out by microsatellite analysis, interdelta sequence typing and pulsed-field gel electrophoresis (PFGE). Genetic patterns were compared to those obtained for 30 isolates of the original commercialized Zymaflore VL1 strain. Among the 100 recovered isolates we found a high percentage of chromosomal size variations, most evident for the smaller chromosomes III and VI. Complete loss of heterozygosity was observed for two isolates that had also lost chromosomal heteromorphism; their growth and fermentative capacity in a synthetic must medium was also affected. A considerably higher number of variant patterns for interdelta sequence amplifications was obtained for grape-derived strains compared to the original VL1 isolates. Our data show that the long-term presence of strain VL1 in natural grapevine environments induced genetic changes that can be detected using different fingerprinting methods. The observed genetic changes may reflect adaptive mechanisms to changed environmental conditions that yeast cells encounter during their existence in nature.

High-resolution mutation mapping reveals parallel experimental evolution in yeast

Understanding the genetic basis of evolutionary adaptation is limited by our ability to efficiently identify the genomic locations of adaptive mutations. Here we describe a method that can quickly and precisely map the genetic basis of naturally and experimentally evolved complex traits using linkage analysis. A yeast strain that expresses the evolved trait is crossed to a distinct strain background and DNA from a large pool of progeny that express the trait of interest is hybridized to oligonucleotide microarrays that detect thousands of polymorphisms between the two strains. Adaptive mutations are detected by linkage to the polymorphisms from the evolved parent. We successfully tested our method by mapping five known genes to a precision of 0.2–24 kb (0.1–10 cM), and developed computer simulations to test the effect of different factors on mapping precision. We then applied this method to four yeast strains that had independently adapted to a fluctuating glucose–galactose environment. All four strains had acquired one or more missense mutations in GAL80, the repressor of the galactose utilization pathway. When transferred into the ancestral strain, the gal80 mutations conferred the fitness advantage that the evolved strains show in the transition from glucose to galactose. Our results show an example of parallel adaptation caused by mutations in the same gene.

An array hybridization method enables genetic mapping via linkage analysis; applied here this new method shows parallel adaptation to a fluctuating glucose-galactose environment evidenced by mutations in the GAL80 gene.

Global mapping of transposon location

Transposable genetic elements are ubiquitous, yet their presence or absence at any given position within a genome can vary between individual cells, tissues, or strains. Transposable elements have profound impacts on host genomes by altering gene expression, assisting in genomic rearrangements, causing insertional mutations, and serving as sources of phenotypic variation. Characterizing a genome's full complement of transposons requires whole genome sequencing, precluding simple studies of the impact of transposition on interindividual variation. Here, we describe a global mapping approach for identifying transposon locations in any genome, using a combination of transposon-specific DNA extraction and microarray-based comparative hybridization analysis. We use this approach to map the repertoire of endogenous transposons in different laboratory strains of Saccharomyces cerevisiae and demonstrate that transposons are a source of extensive genomic variation. We also apply this method to mapping bacterial transposon insertion sites in a yeast genomic library. This unique whole genome view of transposon location will facilitate our exploration of transposon dynamics, as well as defining bases for individual differences and adaptive potential.

Loss of dispensable genes is not adaptive in yeast

A substantial share of genes identified in yeast can be deleted without visible phenotypic effects. Current debate concentrates on the possible roles of seemingly dispensable genes. The costs of maintaining unnecessary functions has attracted little attention. The hypothesis of antagonistic pleiotropy postulates that adaptations to different constituents of the environment are likely to interfere with each other, and therefore loss of unnecessary functions is potentially advantageous. We tested an entire collection of nonessential yeast gene deletions in a benign and nutritionally rich environment in which the number of dispensable genes was particularly high. We applied a series of competition experiments that could detect differences in relative fitness of approximately 0.005. No beneficial deletions were found, except perhaps for the deletion of about a dozen genes that slightly improved competitive ability; however, a functional explanation of the fitness advantage is lacking. The paucity of beneficial gene deletions is striking because genetic adaptations to laboratory conditions are regularly observed in yeast. However, it accords with the finding that the gene contents of four species of Saccharomyces are nearly identical, despite up to 20 million years of independent evolution and extensive DNA sequence divergence. Such extreme conservation of functions would be improbable if there were periods of selection promoting the loss of temporarily dispensable genes. The evident cohesion of the yeast genomes may be their evolved feature or an intrinsic property of complex genetic systems.

Happy together: the life and times of Ty retrotransposons and their hosts

The aim of this review is to describe the level of intimacy between Ty retrotransposons (Ty1-Ty5) and their host the yeast Saccharomyces cerevisiae. The effects of Ty location in the genome and of host proteins on the expression and mobility of Ty elements are highlighted. After a brief overview of Ty diversity and evolution, we describe the factors that dictate Ty target-site preference and the impact of targeting on Ty and adjacent gene expression. Studies on Ty3 and Ty5 have been especially informative in unraveling the role of host factors (Pol III machinery and silencing proteins, respectively) and integrase in controlling the specificity of integration. In contrast, not much is known regarding Ty1, Ty2 and Ty4, except that their insertion depends on the transcriptional competence of the adjacent Pol III gene and might be influenced by some chromatin components. This review also brings together recent findings on the regulation of Ty1 retrotransposition. A large number of host proteins (over 30) involved in a wide range of cellular processes controls either directly or indirectly Ty1 mobility, primarily at post-transcriptional steps. We focus on several genes for which more detailed analyses have permitted the elaboration of regulatory models. In addition, this review describes new data revealing that repression of Ty1 mobility also involves two forms of copy number control that act at both the trancriptional and post-transcriptional levels. Since S. cerevisiae lacks the conserved pathways for copy number control via transcriptional and post-transcriptional gene silencing found in other eukaryotes, Ty1 copy number control must be via another mechanism whose features are outlined. Ty1 response to stress also implicates activation at both transcriptional and postranscriptional steps of Ty1. Finally, we provide several insights in the role of Ty elements in chromosome evolution and yeast adaptation and discuss the factors that might limit Ty ectopic recombination.

Microbial evolution in laboratory environments

The classical model of evolutionary change in microbial populations is evaluated in terms of recent evidence obtained from genetic and molecular analysis of evolving laboratory populations in simple, defined homogeneous environments. Patterns of change are observed which cannot be predicted from the classical model. These observations highlight the value of the study of evolving laboratory populations, which allow a degree of analysis typically not possible for natural populations.

The uses of genome-wide yeast mutant collections

We assess five years of usage of the major genome-wide collections of mutants from Saccharomyces cerevisiae: single deletion mutants, double mutants conferring 'synthetic' lethality and the 'TRIPLES' collection of mutants obtained by random transposon insertion. Over 100 experimental conditions have been tested and more than 5,000 novel phenotypic traits have been assigned to yeast genes using these collections.

Capturing the adaptive mutation in yeast

An accurate view of adaptive mutations is essential to evolutionary genetics, but their rarity makes them difficult to study. This can be partially overcome using the many tools of yeast genetics and the ability to study very large populations over many generations. Adaptation to laboratory environments has occurred primarily by chromosomal rearrangements, often involving retrotransposons and apparently selected for their effects on gene regulation. Estimated rates of adaptive mutation are on the order of 1 in 10(11) cell divisions. There remains great potential for the genomic study of variation within yeast species to contribute to our understanding of adaptive mutation.

Ty1 insertions in intergenic regions of the genome of Saccharomyces cerevisiae transcribed by RNA polymerase III have no detectable selective effect

The retrotransposon Ty1 of Saccharomyces cerevisiae inserts preferentially into intergenic regions in the vicinity of RNA polymerase III-transcribed genes. It has been suggested that this preference has evolved to minimize the deleterious effects of element transposition on the host genome, and thus to favor their evolutionary survival. This presupposes that such insertions have no selective effect. However, there has been no direct test of this hypothesis. Here we construct a series of strains containing single Ty1 insertions in the vicinity of tRNA genes, or in the rDNA cluster on chromosome XII, which are otherwise isogenic to strain 337, containing zero Ty1 elements. Competition experiments between 337 and the strains containing single Ty1 insertions revealed that in all cases, the Ty1 insertions have no selective effect in rich medium. These results are thus consistent with the hypothesis that the insertion site preference of Ty1 elements has evolved to minimize the deleterious effects of transposition on the host genome.