Genetic variation in adaptability and pleiotropy in budding yeast

Environment-independent distribution of mutational effects emerges from microscopic epistasis

Predicting how new mutations alter phenotypes is difficult because mutational effects vary across genotypes and environments. Recently discovered global epistasis, where the fitness effects of mutations scale with the fitness of the background genotype, can improve predictions, but how the environment modulates this scaling is unknown. We measured the fitness effects of ~100 insertion mutations in 42 strains of Saccharomyces cerevisiae in six laboratory environments and found that the global-epistasis scaling is nearly invariant across environments. Instead, the environment tunes one global parameter, the background fitness at which most mutations switch sign. As a consequence, the distribution of mutational effects is predictable across genotypes and environments. Our results suggest that the effective dimensionality of genotype-to-phenotype maps across environments is surprisingly low.

Convergent genetic adaptation of Escherichia coli in minimal media leads to pleiotropic divergence

Adaptation in an environment can either be beneficial, neutral or disadvantageous in another. To test the genetic basis of pleiotropic behaviour, we evolved six lines of E. coli independently in environments where glucose and galactose were the sole carbon sources, for 300 generations. All six lines in each environment exhibit convergent adaptation in the environment in which they were evolved. However, pleiotropic behaviour was observed in several environmental contexts, including other carbon environments. Genome sequencing reveals that mutations in global regulators rpoB and rpoC cause this pleiotropy. We report three new alleles of the rpoB gene, and one new allele of the rpoC gene. The novel rpoB alleles confer resistance to Rifampicin, and alter motility. Our results show how single nucleotide changes in the process of adaptation in minimal media can lead to wide-scale pleiotropy, resulting in changes in traits that are not under direct selection.

Correlated Allele Frequency Changes Reveal Clonal Structure and Selection in Temporal Genetic Data

In evolving populations where the rate of beneficial mutations is large, subpopulations of individuals with competing beneficial mutations can be maintained over long times. Evolution with this kind of clonal structure is commonly observed in a wide range of microbial and viral populations. However, it can be difficult to completely resolve clonal dynamics in data. This is due to limited read lengths in high-throughput sequencing methods, which are often insufficient to directly measure linkage disequilibrium or determine clonal structure. Here, we develop a method to infer clonal structure using correlated allele frequency changes in time-series sequence data. Simulations show that our method recovers true, underlying clonal structures when they are known and accurately estimate linkage disequilibrium. This information can then be combined with other inference methods to improve estimates of the fitness effects of individual mutations. Applications to data suggest novel clonal structures in an E. coli long-term evolution experiment, and yield improved predictions of the effects of mutations on bacterial fitness and antibiotic resistance. Moreover, our method is computationally efficient, requiring orders of magnitude less run time for large data sets than existing methods. Overall, our method provides a powerful tool to infer clonal structures from data sets where only allele frequencies are available, which can also improve downstream analyses.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Experimental evolution for cell biology

Evolutionary cell biology explores the origins, principles, and core functions of cellular features and regulatory networks through the lens of evolution. This emerging field relies heavily on comparative experiments and genomic analyses that focus exclusively on extant diversity and historical events, providing limited opportunities for experimental validation. In this opinion article, we explore the potential for experimental laboratory evolution to augment the evolutionary cell biology toolbox, drawing inspiration from recent studies that combine laboratory evolution with cell biological assays. Primarily focusing on approaches for single cells, we provide a generalizable template for adapting experimental evolution protocols to provide fresh insight into long-standing questions in cell biology.

Detecting natural selection in trait-trait coevolution

No phenotypic trait evolves independently of all other traits, but the cause of trait-trait coevolution is poorly understood. While the coevolution could arise simply from pleiotropic mutations that simultaneously affect the traits concerned, it could also result from multivariate natural selection favoring certain trait relationships. To gain a general mechanistic understanding of trait-trait coevolution, we examine the evolution of 220 cell morphology traits across 16 natural strains of the yeast Saccharomyces cerevisiae and the evolution of 24 wing morphology traits across 110 fly species of the family Drosophilidae, along with the variations of these traits among gene deletion or mutation accumulation lines (a.k.a. mutants). For numerous trait pairs, the phenotypic correlation among evolutionary lineages differs significantly from that among mutants. Specifically, we find hundreds of cases where the evolutionary correlation between traits is strengthened or reversed relative to the mutational correlation, which, according to our population genetic simulation, is likely caused by multivariate selection. Furthermore, we detect selection for enhanced modularity of the yeast traits analyzed. Together, these results demonstrate that trait-trait coevolution is shaped by natural selection and suggest that the pleiotropic structure of mutation is not optimal. Because the morphological traits analyzed here are chosen largely because of their measurability and thereby are not expected to be biased with regard to natural selection, our conclusion is likely general.

Rapid evolutionary repair by secondary perturbation of a primary disrupted transcriptional network

The discrete steps of transcriptional rewiring have been proposed to occur neutrally to ensure steady gene expression under stabilizing selection. A conflict-free switch of a regulon between regulators may require an immediate compensatory evolution to minimize deleterious effects. Here, we perform an evolutionary repair experiment on the Lachancea kluyveri yeast sef1Δ mutant using a suppressor development strategy. Complete loss of SEF1 forces cells to initiate a compensatory process for the pleiotropic defects arising from misexpression of TCA cycle genes. Using different selective conditions, we identify two adaptive loss-of-function mutations of IRA1 and AZF1. Subsequent analyses show that Azf1 is a weak transcriptional activator regulated by the Ras1-PKA pathway. Azf1 loss-of-function triggers extensive gene expression changes responsible for compensatory, beneficial, and trade-off phenotypes. The trade-offs can be alleviated by higher cell density. Our results not only indicate that secondary transcriptional perturbation provides rapid and adaptive mechanisms potentially stabilizing the initial stage of transcriptional rewiring but also suggest how genetic polymorphisms of pleiotropic mutations could be maintained in the population.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.

Homo cerevisiae-Leveraging Yeast for Investigating Protein-Protein Interactions and Their Role in Human Disease

Understanding how genetic variation affects phenotypes represents a major challenge, particularly in the context of human disease. Although numerous disease-associated genes have been identified, the clinical significance of most human variants remains unknown. Despite unparalleled advances in genomics, functional assays often lack sufficient throughput, hindering efficient variant functionalization. There is a critical need for the development of more potent, high-throughput methods for characterizing human genetic variants. Here, we review how yeast helps tackle this challenge, both as a valuable model organism and as an experimental tool for investigating the molecular basis of phenotypic perturbation upon genetic variation. In systems biology, yeast has played a pivotal role as a highly scalable platform which has allowed us to gain extensive genetic and molecular knowledge, including the construction of comprehensive interactome maps at the proteome scale for various organisms. By leveraging interactome networks, one can view biology from a systems perspective, unravel the molecular mechanisms underlying genetic diseases, and identify therapeutic targets. The use of yeast to assess the molecular impacts of genetic variants, including those associated with viral interactions, cancer, and rare and complex diseases, has the potential to bridge the gap between genotype and phenotype, opening the door for precision medicine approaches and therapeutic development.

Epistasis and evolution: recent advances and an outlook for prediction

As organisms evolve, the effects of mutations change as a result of epistatic interactions with other mutations accumulated along the line of descent. This can lead to shifts in adaptability or robustness that ultimately shape subsequent evolution. Here, we review recent advances in measuring, modeling, and predicting epistasis along evolutionary trajectories, both in microbial cells and single proteins. We focus on simple patterns of global epistasis that emerge in this data, in which the effects of mutations can be predicted by a small number of variables. The emergence of these patterns offers promise for efforts to model epistasis and predict evolution.

Experimental approaches to study evolutionary cell biology using yeasts

The past century has witnessed tremendous advances in understanding how cells function. Nevertheless, how cellular processes have evolved is still poorly understood. Many studies have highlighted surprising molecular diversity in how cells from diverse species execute the same processes, and advances in comparative genomics are likely to reveal much more molecular diversity than was believed possible until recently. Extant cells remain therefore the product of an evolutionary history that we vastly ignore. Evolutionary cell biology has emerged as a discipline aiming to address this knowledge gap by combining evolutionary, molecular, and cellular biology thinking. Recent studies have shown how even essential molecular processes, such as DNA replication, can undergo fast adaptive evolution under certain laboratory conditions. These developments open new lines of research where the evolution of cellular processes can be investigated experimentally. Yeasts naturally find themselves at the forefront of this research line. Not only do they allow the observation of fast evolutionary adaptation, but they also provide numerous genomic, synthetic, and cellular biology tools already developed by a large community. Here we propose that yeasts can serve as an "evolutionary cell lab" to test hypotheses, principles, and ideas in evolutionary cell biology. We discuss various experimental approaches available for this purpose, and how biology at large can benefit from them.

Adaptation of the yeast gene knockout collection is near-perfectly predicted by fitness and diminishing return epistasis

Adaptive evolution of clonally dividing cells and microbes is the ultimate cause of cancer and infectious diseases. The possibility of constraining the adaptation of cell populations, by inhibiting proteins enhancing the evolvability, has therefore attracted interest. However, our current understanding of how genes influence adaptation kinetics is limited, partly because accurately measuring adaptation for many cell populations is challenging. We used a high-throughput adaptive laboratory evolution platform to track the adaptation of >18,000 cell populations corresponding to single-gene deletion strains in the haploid yeast deletion collection. We report that the preadaptation fitness of gene knockouts near-perfectly (R2= 0.91) predicts their adaptation to arsenic, leaving at the most a marginal role for dedicated evolvability gene functions. We tracked the adaptation of another >23,000 gene knockout populations to a diverse range of selection pressures and generalized the almost perfect (R2=0.72-0.98) capacity of preadaptation fitness to predict adaptation. We also reconstructed mutations in FPS1, ASK10, and ARR3, which together account for almost all arsenic adaptation in wild-type cells, in gene deletions covering a broad fitness range and show that the predictability of arsenic adaptation can be understood as a by global epistasis, where excluding arsenic is more beneficial to arsenic unfit cells. The paucity of genes with a meaningful evolvability effect on adaptation diminishes the prospects of developing adjuvant drugs aiming to slow antimicrobial and chemotherapy resistance.

Adaptive evolvability through direct selection instead of indirect, second-order selection

Can evolvability itself be the product of adaptive evolution? To answer this question is challenging, because any DNA mutation that alters only evolvability is subject to indirect, "second order" selection on the future effects of this mutation. Such indirect selection is weaker than "first-order" selection on mutations that alter fitness, in the sense that it can operate only under restrictive conditions. Here I discuss a route to adaptive evolvability that overcomes this challenge. Specifically, a recent evolution experiment showed that some mutations can enhance both fitness and evolvability through a combination of direct and indirect selection. Unrelated evidence from gene duplication and the evolution of gene regulation suggests that mutations with such dual effects may not be rare. Through such mutations, evolvability may increase at least in part because it provides an adaptive advantage. These observations suggest a research program on the adaptive evolution of evolvability, which aims to identify such mutations and to disentangle their direct fitness effects from their indirect effects on evolvability. If evolvability is itself adaptive, Darwinian evolution may have created more than life's diversity. It may also have helped create the very conditions that made the success of Darwinian evolution possible.

Challenges and potential solutions for studying the genetic and phenotypic architecture of adaptation in microbes

All organisms are defined by the makeup of their DNA. Over billions of years, the structure and information contained in that DNA, often referred to as genetic architecture, have been honed by a multitude of evolutionary processes. Mutations that cause genetic elements to change in a way that results in beneficial phenotypic change are more likely to survive and propagate through the population in a process known as adaptation. Recent work reveals that the genetic targets of adaptation are varied and can change with genetic background. Further, seemingly similar adaptive mutations, even within the same gene, can have diverse and unpredictable effects on phenotype. These challenges represent major obstacles in predicting adaptation and evolution. In this review, we cover these concepts in detail and identify three emerging synergistic solutions: higher-throughput evolution experiments combined with updated genotype-phenotype mapping strategies and physiological models. Our review largely focuses on recent literature in yeast, and the field seems to be on the cusp of a new era with regard to studying the predictability of evolution.

Changes in the growth rate of Chlamydomonas reinhardtii under long-term selection by temperature and salinity: Acclimation vs. evolution

We investigated the roles of acclimation and different components involved in evolution (adaptation, chance and history) on the changes in the growth rate of the model freshwater microalga Chlamydomonas reinhardtii P. A. Dang. exposed to selective temperature and salinity. Three C. reinhardtii strains previously grown during one year in freshwater medium and 20 °C were exposed to 5 °C temperature increase and a salinity of 5 g L^-1 NaCl. Cultures under each selective scenario and in combination (increase of salinity and temperature), were propagated until growth rate achieved an invariant mean value for 6 months (100-350 generations, varying as a function of scenario and strain). The changes of the growth rate under increased temperature were due to both adaptation and acclimation, as well as history. However, acclimation was the only mechanism detected under salinity increase as well as in the selective scenario of both temperature and salinity, suggesting that genetic variability would not allow survival at salinity higher than that to which experimental populations were exposed. Therefore, it could be hypothesized that under a global change scenario an increase in salinity would be a greater challenge than warming for some freshwater phytoplankton.

Overdominant and partially dominant mutations drive clonal adaptation in diploid Saccharomyces cerevisiae

Identification of adaptive targets in experimental evolution typically relies on extensive replication and genetic reconstruction. An alternative approach is to directly assay all mutations in an evolved clone by generating pools of segregants that contain random combinations of evolved mutations. Here, we apply this method to six Saccharomyces cerevisiae clones isolated from four diploid populations that were clonally evolved for 2,000 generations in rich glucose medium. Each clone contains 17-26 mutations relative to the ancestor. We derived intermediate genotypes between the founder and the evolved clones by bulk mating sporulated cultures of the evolved clones to a barcoded haploid version of the ancestor. We competed the resulting barcoded diploids en masse and quantified fitness in the experimental and alternative environments by barcode sequencing. We estimated average fitness effects of evolved mutations using barcode-based fitness assays and whole genome sequencing for a subset of segregants. In contrast to our previous work with haploid evolved clones, we find that diploids carry fewer beneficial mutations, with modest fitness effects (up to 5.4%) in the environment in which they arose. In agreement with theoretical expectations, reconstruction experiments show that all mutations with a detectable fitness effect manifest some degree of dominance over the ancestral allele, and most are overdominant. Genotypes with lower fitness effects in alternative environments allowed us to identify conditions that drive adaptation in our system.

Ploidy and recombination proficiency shape the evolutionary adaptation to constitutive DNA replication stress

In haploid budding yeast, evolutionary adaptation to constitutive DNA replication stress alters three genome maintenance modules: DNA replication, the DNA damage checkpoint, and sister chromatid cohesion. We asked how these trajectories depend on genomic features by comparing the adaptation in three strains: haploids, diploids, and recombination deficient haploids. In all three, adaptation happens within 1000 generations at rates that are correlated with the initial fitness defect of the ancestors. Mutations in individual genes are selected at different frequencies in populations with different genomic features, but the benefits these mutations confer are similar in the three strains, and combinations of these mutations reproduce the fitness gains of evolved populations. Despite the differences in the selected mutations, adaptation targets the same three functional modules in strains with different genomic features, revealing a common evolutionary response to constitutive DNA replication stress.

Temperature affects the repeatability of evolution in the microbial eukaryote Tetrahymena thermophila

Evolutionary biologists have long sought to understand what factors affect the repeatability of adaptive outcomes. To better understand the role of temperature in determining the repeatability of adaptive trajectories, we evolved populations of different genotypes of the ciliate Tetrahymena thermophila at low and high temperatures and followed changes in growth rate over 6,500 generations. As expected, growth rate increased with a decelerating rate for all populations; however, there were differences in the patterns of evolution at the two temperatures. The growth rates of the different genotypes tended to converge as evolution proceeded at both temperatures, but this convergence was quicker and more pronounced at the higher temperature. Additionally, over the first 4,000 generations we found greater repeatability of evolution, in terms of change in growth rate, among replicates of the same genotype at the higher temperature. Finally, we found limited evidence of trade-offs in fitness between temperatures, and an asymmetry in the correlated responses, whereby evolution in a high temperature increases growth rate at the lower temperature significantly more than the reverse. These results demonstrate the importance of temperature in determining the repeatability of evolutionary trajectories for the eukaryotic microbe Tetrahymena thermophila and may provide clues to how temperature affects evolution more generally.

Global epistasis emerges from a generic model of a complex trait

Epistasis between mutations can make adaptation contingent on evolutionary history. Yet despite widespread ‘microscopic’ epistasis between the mutations involved, microbial evolution experiments show consistent patterns of fitness increase between replicate lines. Recent work shows that this consistency is driven in part by global patterns of diminishing-returns and increasing-costs epistasis, which make mutations systematically less beneficial (or more deleterious) on fitter genetic backgrounds. However, the origin of this ‘global’ epistasis remains unknown. Here, we show that diminishing-returns and increasing-costs epistasis emerge generically as a consequence of pervasive microscopic epistasis. Our model predicts a specific quantitative relationship between the magnitude of global epistasis and the stochastic effects of microscopic epistasis, which we confirm by reanalyzing existing data. We further show that the distribution of fitness effects takes on a universal form when epistasis is widespread and introduce a novel fitness landscape model to show how phenotypic evolution can be repeatable despite sequence-level stochasticity.

Cancer cells employ an evolutionarily conserved polyploidization program to resist therapy

Unusually large cancer cells with abnormal nuclei have been documented in the cancer literature since 1858. For more than 100 years, they have been generally disregarded as irreversibly senescent or dying cells, too morphologically misshapen and chromatin too disorganized to be functional. Cell enlargement, accompanied by whole genome doubling or more, is observed across organisms, often associated with mitigation strategies against environmental change, severe stress, or the lack of nutrients. Our comparison of the mechanisms for polyploidization in other organisms and non-transformed tissues suggest that cancer cells draw from a conserved program for their survival, utilizing whole genome doubling and pausing proliferation to survive stress. These polyaneuploid cancer cells (PACCs) are the source of therapeutic resistance, responsible for cancer recurrence and, ultimately, cancer lethality.

High-throughput analysis of adaptation using barcoded strains of Saccharomyces cerevisiae

BACKGROUND

Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue. Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies. We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge.

RESULTS

We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays. Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations. In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics.

CONCLUSIONS

Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking. These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution. The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Effects of Genetic and Physiological Divergence on the Evolution of a Sulfate-Reducing Bacterium under Conditions of Elevated Temperature

Improving our understanding of how previous adaptation influences evolution has been a long-standing goal in evolutionary biology. Natural selection tends to drive populations to find similar adaptive solutions for the same selective conditions. However, variations in historical environments can lead to both physiological and genetic divergence that can make evolution unpredictable. Here, we assessed the influence of divergence on the evolution of a model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough, in response to elevated temperature and found a significant effect at the genetic but not the phenotypic level. Understanding how these influences drive evolution will allow us to better predict how bacteria will adapt to various ecological constraints.

Adaptation via natural selection is an important driver of evolution, and repeatable adaptations of replicate populations, under conditions of a constant environment, have been extensively reported. However, isolated groups of populations in nature tend to harbor both genetic and physiological divergence due to multiple selective pressures that they have encountered. How this divergence affects adaptation of these populations to a new common environment remains unclear. To determine the impact of prior genetic and physiological divergence in shaping adaptive evolution to accommodate a new common environment, an experimental evolution study with the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough (DvH) was conducted. Two groups of replicate populations with genetic and physiological divergence, derived from a previous evolution study, were propagated in an elevated-temperature environment for 1,000 generations. Ancestor populations without prior experimental evolution were also propagated in the same environment as a control. After 1,000 generations, all the populations had increased growth rates and all but one had greater fitness in the new environment than the ancestor population. Moreover, improvements in growth rate were moderately affected by the divergence in the starting populations, while changes in fitness were not significantly affected. The mutations acquired at the gene level in each group of populations were quite different, indicating that the observed phenotypic changes were achieved by evolutionary responses that differed between the groups. Overall, our work demonstrated that the initial differences in fitness between the starting populations were eliminated by adaptation and that phenotypic convergence was achieved by acquisition of mutations in different genes.

Evolutionary stalling and a limit on the power of natural selection to improve a cellular module

Cells consist of molecular modules which perform vital biological functions. Cellular modules are key units of adaptive evolution because organismal fitness depends on their performance. Theory shows that in rapidly evolving populations, such as those of many microbes, adaptation is driven primarily by common beneficial mutations with large effects, while other mutations behave as if they are effectively neutral. As a consequence, if a module can be improved only by rare and/or weak beneficial mutations, its adaptive evolution would stall. However, such evolutionary stalling has not been empirically demonstrated, and it is unclear to what extent stalling may limit the power of natural selection to improve modules. Here we empirically characterize how natural selection improves the translation machinery (TM), an essential cellular module. We experimentally evolved populations of Escherichia coli with genetically perturbed TMs for 1,000 generations. Populations with severe TM defects initially adapted via mutations in the TM, but TM adaptation stalled within about 300 generations. We estimate that the genetic load in our populations incurred by residual TM defects ranges from 0.5 to 19%. Finally, we found evidence that both epistasis and the depletion of the pool of beneficial mutations contributed to evolutionary stalling. Our results suggest that cellular modules may not be fully optimized by natural selection despite the availability of adaptive mutations.

The Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses

Pleiotropy, the concept that a single gene controls multiple distinct traits, is prevalent in most organisms and has broad implications for medicine and agriculture. The identification of the molecular mechanisms underlying pleiotropy has the power to reveal previously unknown biological connections between seemingly unrelated traits. Additionally, the discovery of pleiotropic genes increases our understanding of both genetic and phenotypic complexity by characterizing novel gene functions. Quantitative trait locus (QTL) mapping has been used to identify several pleiotropic regions in many organisms. However, gene knockout studies are needed to eliminate the possibility of tightly linked, non-pleiotropic loci. Here, we use a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to identify a single large-effect QTL on the center of chromosome V associated with variation in responses to eight chemotherapeutics. We validate this QTL with near-isogenic lines and pair genome-wide gene expression data with drug response traits to perform mediation analysis, leading to the identification of a pleiotropic candidate gene, scb-1, for some of the eight chemotherapeutics. Using deletion strains created by genome editing, we show that scb-1, which was previously implicated in response to bleomycin, also underlies responses to other double-strand DNA break-inducing chemotherapeutics. This finding provides new evidence for the role of scb-1 in the nematode drug response and highlights the power of mediation analysis to identify causal genes.

Candida albicans Genetic Background Influences Mean and Heterogeneity of Drug Responses and Genome Stability during Evolution in Fluconazole

Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

The importance of within-species diversity in determining the evolutionary potential of a population to evolve drug resistance or tolerance is not well understood, including in eukaryotic pathogens. To examine the influence of genetic background, we evolved replicates of 20 different clinical isolates of Candida albicans, a human fungal pathogen, in fluconazole, the commonly used antifungal drug. The isolates hailed from the major C. albicans clades and had different initial levels of drug resistance and tolerance to the drug. The majority of replicates rapidly increased in fitness in the evolutionary environment, with the degree of improvement inversely correlated with parental strain fitness in the drug. Improvement was largely restricted to up to the evolutionary level of drug: only 4% of the evolved replicates increased resistance (MIC) above the evolutionary level of drug. Prevalent changes were altered levels of drug tolerance (slow growth of a subpopulation of cells at drug concentrations above the MIC) and increased diversity of genome size. The prevalence and predominant direction of these changes differed in a strain-specific manner, but neither correlated directly with parental fitness or improvement in fitness. Rather, low parental strain fitness was correlated with high levels of heterogeneity in fitness, tolerance, and genome size among evolved replicates. Thus, parental strain background is an important determinant in mean improvement to the evolutionary environment as well as the diversity of evolved phenotypes, and the range of possible responses of a pathogen to an antimicrobial drug cannot be captured by in-depth study of a single strain background.

IMPORTANCE Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

Predicting Evolution Using Regulatory Architecture

The limits of evolution have long fascinated biologists. However, the causes of evolutionary constraint have remained elusive due to a poor mechanistic understanding of studied phenotypes. Recently, a range of innovative approaches have leveraged mechanistic information on regulatory networks and cellular biology. These methods combine systems biology models with population and single-cell quantification and with new genetic tools, and they have been applied to a range of complex cellular functions and engineered networks. In this article, we review these developments, which are revealing the mechanistic causes of epistasis at different levels of biological organization-in molecular recognition, within a single regulatory network, and between different networks-providing first indications of predictable features of evolutionary constraint.

Chance and necessity in the pleiotropic consequences of adaptation for budding yeast

Mutations that a population accumulates during evolution in one 'home' environment may cause fitness gains or losses in other environments. Such pleiotropic fitness effects determine the evolutionary fate of the population in variable environments and can lead to ecological specialization. It is unclear how the pleiotropic outcomes of evolution are shaped by the intrinsic randomness of the evolutionary process and by the deterministic variation in selection pressures across environments. Here, to address this question, we evolved 20 replicate populations of the yeast Saccharomyces cerevisiae in 11 laboratory environments and measured their fitness across multiple conditions. We found that evolution led to diverse pleiotropic fitness gains and losses, driven by multiple types of mutations. Approximately 60% of this variation is explained by the home environment of a clone and the most common parallel genetic changes, whereas about 40% is attributed to the stochastic accumulation of mutations whose pleiotropic effects are unpredictable. Although populations are typically specialized to their home environment, generalists also evolved in almost all of the conditions. Our results suggest that the mutations that accumulate during evolution incur a variety of pleiotropic costs and benefits with different probabilities. Thus, whether a population evolves towards a specialist or a generalist phenotype is heavily influenced by chance.

Evolutionary repair: Changes in multiple functional modules allow meiotic cohesin to support mitosis

The role of proteins often changes during evolution, but we do not know how cells adapt when a protein is asked to participate in a different biological function. We forced the budding yeast, Saccharomyces cerevisiae, to use the meiosis-specific kleisin, recombination 8 (Rec8), during the mitotic cell cycle, instead of its paralog, Scc1. This perturbation impairs sister chromosome linkage, advances the timing of genome replication, and reduces reproductive fitness by 45%. We evolved 15 parallel populations for 1,750 generations, substantially increasing their fitness, and analyzed the genotypes and phenotypes of the evolved cells. Only one population contained a mutation in Rec8, but many populations had mutations in the transcriptional mediator complex, cohesin-related genes, and cell cycle regulators that induce S phase. These mutations improve sister chromosome cohesion and delay genome replication in Rec8-expressing cells. We conclude that changes in known and novel partners allow cells to use an existing protein to participate in new biological functions.

Higher-fitness yeast genotypes are less robust to deleterious mutations

Natural selection drives populations toward higher fitness, but second-order selection for adaptability and mutational robustness can also influence evolution. In many microbial systems, diminishing-returns epistasis contributes to a tendency for more-fit genotypes to be less adaptable, but no analogous patterns for robustness are known. To understand how robustness varies across genotypes, we measure the fitness effects of hundreds of individual insertion mutations in a panel of yeast strains. We find that more-fit strains are less robust: They have distributions of fitness effects with lower mean and higher variance. These differences arise because many mutations have more strongly deleterious effects in faster-growing strains. This negative correlation between fitness and robustness implies that second-order selection for robustness will tend to conflict with first-order selection for fitness.

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.

Adaptive walks on high-dimensional fitness landscapes and seascapes with distance-dependent statistics

The dynamics of evolution is intimately shaped by epistasis - interactions between genetic elements which cause the fitness-effect of combinations of mutations to be non-additive. Analyzing evolutionary dynamics that involves large numbers of epistatic mutations is intrinsically difficult. A crucial feature is that the fitness landscape in the vicinity of the current genome depends on the evolutionary history. A key step is thus developing models that enable study of the effects of past evolution on future evolution. In this work, we introduce a broad class of high-dimensional random fitness landscapes for which the correlations between fitnesses of genomes are a general function of genetic distance. Their Gaussian character allows for tractable computational as well as analytic understanding. We study the properties of these landscapes focusing on the simplest evolutionary process: random adaptive (uphill) walks. Conventional measures of "ruggedness" are shown to not much affect such adaptive walks. Instead, the long-distance statistics of epistasis cause all properties to be highly conditional on past evolution, determining the statistics of the local landscape (the distribution of fitness-effects of available mutations and combinations of these), as well as the global geometry of evolutionary trajectories. In order to further explore the effects of conditioning on past evolution, we model the effects of slowly changing environments. At long times, such fitness "seascapes" cause a statistical steady state with highly intermittent evolutionary dynamics: populations undergo bursts of rapid adaptation, interspersed with periods in which adaptive mutations are rare and the population waits for more new directions to be opened up by changes in the environment. Finally, we discuss prospects for studying more complex evolutionary dynamics and on broader classes of high-dimensional landscapes and seascapes.

Cryptic genetic variation accelerates evolution by opening access to diverse adaptive peaks

Cryptic genetic variation can facilitate adaptation in evolving populations. To elucidate the underlying genetic mechanisms, we used directed evolution in Escherichia coli to accumulate variation in populations of yellow fluorescent proteins and then evolved these proteins toward the new phenotype of green fluorescence. Populations with cryptic variation evolved adaptive genotypes with greater diversity and higher fitness than populations without cryptic variation, which converged on similar genotypes. Populations with cryptic variation accumulated neutral or deleterious mutations that break the constraints on the order in which adaptive mutations arise. In doing so, cryptic variation opens paths to adaptive genotypes, creates historical contingency, and reduces the predictability of evolution by allowing different replicate populations to climb different adaptive peaks and explore otherwise-inaccessible regions of an adaptive landscape.

Microbial Experimental Evolution - a proving ground for evolutionary theory and a tool for discovery

Microbial experimental evolution uses controlled laboratory populations to study the mechanisms of evolution. The molecular analysis of evolved populations enables empirical tests that can confirm the predictions of evolutionary theory, but can also lead to surprising discoveries. As with other fields in the life sciences, microbial experimental evolution has become a tool, deployed as part of the suite of techniques available to the molecular biologist. Here, I provide a review of the general findings of microbial experimental evolution, especially those relevant to molecular microbiologists that are new to the field. I also relate these results to design considerations for an evolution experiment and suggest future directions for those working at the intersection of experimental evolution and molecular biology.

The causes of evolvability and their evolution

Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.

IAMBEE: a web-service for the identification of adaptive pathways from parallel evolved clonal populations

IAMBEE is a web server designed for the Identification of Adaptive Mutations in Bacterial Evolution Experiments (IAMBEE). Input data consist of genotype information obtained from independently evolved clonal populations or strains that show the same adapted behavior (phenotype). To distinguish adaptive from passenger mutations, IAMBEE searches for neighborhoods in an organism-specific interaction network that are recurrently mutated in the adapted populations. This search for recurrently mutated network neighborhoods, as proxies for pathways is driven by additional information on the functional impact of the observed genetic changes and their dynamics during adaptive evolution. In addition, the search explicitly accounts for the differences in mutation rate between the independently evolved populations. Using this approach, IAMBEE allows exploiting parallel evolution to identify adaptive pathways. The web-server is freely available at http://bioinformatics.intec.ugent.be/iambee/ with no login requirement.

Patterns and Mechanisms of Diminishing Returns from Beneficial Mutations

Diminishing returns epistasis causes the benefit of the same advantageous mutation smaller in fitter genotypes and is frequently observed in experimental evolution. However, its occurrence in other contexts, environment dependence, and mechanistic basis are unclear. Here, we address these questions using 1,005 sequenced segregants generated from a yeast cross. Under each of 47 examined environments, 66-92% of tested polymorphisms exhibit diminishing returns epistasis. Surprisingly, improving environment quality also reduces the benefits of advantageous mutations even when fitness is controlled for, indicating the necessity to revise the global epistasis hypothesis. We propose that diminishing returns originates from the modular organization of life where the contribution of each functional module to fitness is determined jointly by the genotype and environment and has an upper limit, and demonstrate that our model predictions match empirical observations. These findings broaden the concept of diminishing returns epistasis, reveal its generality and potential cause, and have important evolutionary implications.

Modular epistasis and the compensatory evolution of gene deletion mutants

Screens for epistatic interactions have long been used to characterize functional relationships corresponding to protein complexes, metabolic pathways, and other functional modules. Although epistasis between adaptive mutations is also common in laboratory evolution experiments, the functional basis for these interactions is less well characterized. Here, we quantify the extent to which gene function (as determined by a genome-wide screen for epistasis among deletion mutants) influences the rate and genetic basis of compensatory adaptation in a set of 37 gene deletion mutants nested within 16 functional modules. We find that functional module has predictive power: mutants with deletions in the same module tend to adapt more similarly, on average, than those with deletions in different modules. At the same time, initial fitness also plays a role: independent of the specific functional modules involved, adaptive mutations tend to be systematically more beneficial in less-fit genetic backgrounds, consistent with a general pattern of diminishing returns epistasis. We measured epistatic interactions between initial gene deletion mutations and the mutations that accumulate during compensatory adaptation and found a general trend towards positive epistasis (i.e. mutations tend to be more beneficial in the background in which they arose). In two functional modules, epistatic interactions between the initial gene deletions and the mutations in their descendant lines caused evolutionary entrenchment, indicating an intimate functional relationship. Our results suggest that genotypes with similar epistatic interactions with gene deletion mutations will also have similar epistatic interactions with adaptive mutations, meaning that genome scale maps of epistasis between gene deletion mutations can be predictive of evolutionary dynamics.

The effects of mutations often depend on the presence or absence of other mutations. This phenomenon, known as epistasis, has been used extensively to infer functional associations between genes. For example, genes that participate in the same functional module will often show a characteristic pattern of positive epistasis where the knockout of one gene will mask the deleterious effects of knockouts in the other genes. In the context of adaptation, epistasis can cause the outcomes of evolution to depend strongly on the initial genotype. Although studies have found that epistasis is common in laboratory populations, we do not know the extent to which the patterns of epistasis that reveal functional associations overlap with the patterns of epistasis that are important in evolution. Here, by comparing evolution in strains with gene deletions in different functional modules, we quantify the effect of functional epistasis on evolutionary outcomes. We find that mutants with deletions in the same module have more similar evolutionary outcomes, on average, than mutants with deletions in different modules. This suggests that screens for epistasis between gene deletion mutations will not only reveal functional interactions between those genes but may also predict evolutionary dynamics.

Existing Host Range Mutations Constrain Further Emergence of RNA Viruses

RNA viruses are capable of rapid host shifting, typically due to a point mutation that confers expanded host range. As additional point mutations are necessary for further expansions, epistasis among host range mutations can potentially affect the mutational neighborhood and frequency of niche expansion. We mapped the mutational neighborhood of host range expansion using three genotypes of the double-stranded RNA (dsRNA) bacteriophage φ6 (wild type and two isogenic host range mutants) on the novel host Pseudomonas syringae pv. atrofaciens. Both Sanger sequencing of 50 P. syringae pv. atrofaciens mutant clones for each genotype and population Illumina sequencing revealed the same high-frequency mutations allowing infection of P. syringae pv. atrofaciens. Wild-type φ6 had at least nine different ways of mutating to enter the novel host, eight of which are in p3 (host attachment protein gene), and 13/50 clones had unchanged p3 genes. However, the two isogenic mutants had dramatically restricted neighborhoods: only one or two mutations, all in p3. Deep sequencing revealed that wild-type clones without mutations in p3 likely had changes in p12 (morphogenic protein), a region that was not polymorphic for the two isogenic host range mutants. Sanger sequencing confirmed that 10/13 of the wild-type φ6 clones had nonsynonymous mutations in p12, and 2 others had point mutations in p9 and p5. None of these genes had previously been associated with host range expansion in φ6. We demonstrate, for the first time, epistatic constraint in an RNA virus due to host range mutations themselves, which has implications for models of serial host range expansion.IMPORTANCE RNA viruses mutate rapidly and frequently expand their host ranges to infect novel hosts, leading to serial host shifts. Using an RNA bacteriophage model system (Pseudomonas phage φ6), we studied the impact of preexisting host range mutations on another host range expansion. Results from both clonal Sanger and Illumina sequencing show that extant host range mutations dramatically narrow the neighborhood of potential host range mutations compared to that of wild-type φ6. This research suggests that serial host-shifting viruses may follow a small number of molecular paths to enter additional novel hosts. We also identified new genes involved in φ6 host range expansion, expanding our knowledge of this important model system in experimental evolution.

The Genetic Basis of Mutation Rate Variation in Yeast

Mutations are the root source of genetic variation and underlie the process of evolution. Although the rates at which mutations occur vary considerably between species, little is known about differences within species, or the genetic and molecular basis of these differences. Here, we leveraged the power of the yeast Saccharomyces cerevisiae as a model system to uncover natural genetic variants that underlie variation in mutation rate. We developed a high-throughput fluctuation assay and used it to quantify mutation rates in seven natural yeast isolates and in 1040 segregant progeny from a cross between BY, a laboratory strain, and RM, a wine strain. We observed that mutation rate varies among yeast strains and is heritable (H² = 0.49). We performed linkage mapping in the segregants and identified four quantitative trait loci underlying mutation rate variation in the cross. We fine-mapped two quantitative trait loci to the underlying causal genes, RAD5 and MKT1, that contribute to mutation rate variation. These genes also underlie sensitivity to the DNA-damaging agents 4NQO and MMS, suggesting a connection between spontaneous mutation rate and mutagen sensitivity.

Contingency and determinism in evolution: Replaying life’s tape

Historical processes display some degree of “contingency,” meaning their outcomes are sensitive to seemingly inconsequential events that can fundamentally change the future. Contingency is what makes historical outcomes unpredictable. Unlike many other natural phenomena, evolution is a historical process. Evolutionary change is often driven by the deterministic force of natural selection, but natural selection works upon variation that arises unpredictably through time by random mutation, and even beneficial mutations can be lost by chance through genetic drift. Moreover, evolution has taken place within a planetary environment with a particular history of its own. This tension between determinism and contingency makes evolutionary biology a kind of hybrid between science and history. While philosophers of science examine the nuances of contingency, biologists have performed many empirical studies of evolutionary repeatability and contingency. Here, we review the experimental and comparative evidence from these studies. Replicate populations in evolutionary “replay” experiments often show parallel changes, especially in overall performance, although idiosyncratic outcomes show that the particulars of a lineage’s history can affect which of several evolutionary paths is taken. Comparative biologists have found many notable examples of convergent adaptation to similar conditions, but quantification of how frequently such convergence occurs is difficult. On balance, the evidence indicates that evolution tends to be surprisingly repeatable among closely related lineages, but disparate outcomes become more likely as the footprint of history grows deeper. Ongoing research on the structure of adaptive landscapes is providing additional insight into the interplay of fate and chance in the evolutionary process.

Effective models and the search for quantitative principles in microbial evolution

Microbes evolve rapidly. Yet they do so in idiosyncratic ways, which depend on the specific mutations that are beneficial or deleterious in a given situation. At the same time, some population-level patterns of adaptation are strikingly similar across different microbial systems, suggesting that there may also be simple, quantitative principles that unite these diverse scenarios. We review the search for simple principles in microbial evolution, ranging from the biophysical level to emergent evolutionary dynamics. A key theme has been the use of effective models, which coarse-grain over molecular and cellular details to obtain a simpler description in terms of a few effective parameters. Collectively, these theoretical approaches provide a set of quantitative principles that facilitate understanding, prediction, and potentially control of evolutionary phenomena, though formidable challenges remain due to the ecological complexity of natural populations.

The complex underpinnings of genetic background effects

Genetic interactions between mutations and standing polymorphisms can cause mutations to show distinct phenotypic effects in different individuals. To characterize the genetic architecture of these so-called background effects, we genotype 1411 wild-type and mutant yeast cross progeny and measure their growth in 10 environments. Using these data, we map 1086 interactions between segregating loci and 7 different gene knockouts. Each knockout exhibits between 73 and 543 interactions, with 89% of all interactions involving higher-order epistasis between a knockout and multiple loci. Identified loci interact with as few as one knockout and as many as all seven knockouts. In mutants, loci interacting with fewer and more knockouts tend to show enhanced and reduced phenotypic effects, respectively. Cross–environment analysis reveals that most interactions between the knockouts and segregating loci also involve the environment. These results illustrate the complicated interactions between mutations, standing polymorphisms, and the environment that cause background effects.

Mutations often show distinct phenotypic effects across different genetic backgrounds. Here the authors describe the genetic basis of these so-called background effects using data on genotype and growth in 10 environments from 1411 segregants from a cross of two strains of budding yeast.

Power law fitness landscapes and their ability to predict fitness

Whether or not evolution by natural selection is predictable depends on the existence of general patterns shaping the way mutations interact with the genetic background. This interaction, also known as epistasis, has been observed during adaptation (macroscopic epistasis) and in individual mutations (microscopic epistasis). Interestingly, a consistent negative correlation between the fitness effect of beneficial mutations and background fitness (known as diminishing returns epistasis) has been observed across different species and conditions. We tested whether the adaptation pattern of an additional species, Schizosaccharomyces pombe, followed the same trend. We used strains that differed by the presence of large karyotype differences and observed the same pattern of fitness convergence. Using these data along with published datasets, we measured the ability of different models to describe adaptation rates. We found that a phenotype-fitness landscape shaped like a power law is able to correctly predict adaptation dynamics in a variety of species and conditions. Furthermore we show that this model can provide a link between the observed macroscopic and microscopic epistasis. It may be very useful in the development of algorithms able to predict the adaptation of microorganisms from measures of the current phenotypes. Overall, our results suggest that even though adaptation quickly slows down, populations adapting to lab conditions may be quite far from a fitness peak.

Experimental Studies of Evolutionary Dynamics in Microbes

Evolutionary dynamics in laboratory microbial evolution experiments can be surprisingly complex. In the past two decades, observations of these dynamics have challenged simple models of adaptation and have shown that clonal interference, hitchhiking, ecological diversification, and contingency are widespread. In recent years, advances in high-throughput strain maintenance and phenotypic assays, the dramatically reduced cost of genome sequencing, and emerging methods for lineage barcoding have made it possible to observe evolutionary dynamics at unprecedented resolution. These new methods can now begin to provide detailed measurements of key aspects of fitness landscapes and of evolutionary outcomes across a range of systems. These measurements can highlight challenges to existing theoretical models and guide new theoretical work towards the complications that are most widely important.

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.