Cryptic fitness advantage: diploids invade haploid populations despite lacking any apparent advantage as measured by standard fitness assays

Ploidy evolution in a wild yeast is linked to an interaction between cell type and metabolism

Ploidy is an evolutionarily labile trait, and its variation across the tree of life has profound impacts on evolutionary trajectories and life histories. The immediate consequences and molecular causes of ploidy variation on organismal fitness are frequently less clear, although extreme mating type skews in some fungi hint at links between cell type and adaptive traits. Here, we report an unusual recurrent ploidy reduction in replicate populations of the budding yeast Saccharomyces eubayanus experimentally evolved for improvement of a key metabolic trait, the ability to use maltose as a carbon source. We find that haploids have a substantial, but conditional, fitness advantage in the absence of other genetic variation. Using engineered genotypes that decouple the effects of ploidy and cell type, we show that increased fitness is primarily due to the distinct transcriptional program deployed by haploid-like cell types, with a significant but smaller contribution from absolute ploidy. The link between cell-type specification and the carbon metabolism adaptation can be traced to the noncanonical regulation of a maltose transporter by a haploid-specific gene. This study provides novel mechanistic insight into the molecular basis of an environment-cell type fitness interaction and illustrates how selection on traits unexpectedly linked to ploidy states or cell types can drive karyotypic evolution in fungi.

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.

A Saccharomyces eubayanus haploid resource for research studies

Since its identification, Saccharomyces eubayanus has been recognized as the missing parent of the lager hybrid, S. pastorianus. This wild yeast has never been isolated from fermentation environments, thus representing an interesting candidate for evolutionary, ecological and genetic studies. However, it is imperative to develop additional molecular genetics tools to ease manipulation and thus facilitate future studies. With this in mind, we generated a collection of stable haploid strains representative of three main lineages described in S. eubayanus (PB-1, PB-2 and PB-3), by deleting the HO gene using CRISPR-Cas9 and tetrad micromanipulation. Phenotypic characterization under different conditions demonstrated that the haploid derivates were extremely similar to their parental strains. Genomic analysis in three strains highlighted a likely low frequency of off-targets, and sequencing of a single tetrad evidenced no structural variants in any of the haploid spores. Finally, we demonstrate the utilization of the haploid set by challenging the strains under mass-mating conditions. In this way, we found that S. eubayanus under liquid conditions has a preference to remain in a haploid state, unlike S. cerevisiae that mates rapidly. This haploid resource is a novel set of strains for future yeast molecular genetics studies.

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.

The genetic basis of differential autodiploidization in evolving yeast populations

Spontaneous whole-genome duplication, or autodiploidization, is a common route to adaptation in experimental evolution of haploid budding yeast populations. The rate at which autodiploids fix in these populations appears to vary across strain backgrounds, but the genetic basis of these differences remains poorly characterized. Here, we show that the frequency of autodiploidization differs dramatically between two closely related laboratory strains of Saccharomyces cerevisiae, BY4741 and W303. To investigate the genetic basis of this difference, we crossed these strains to generate hundreds of unique F1 segregants and tested the tendency of each segregant to autodiplodize across hundreds of generations of laboratory evolution. We find that variants in the SSD1 gene are the primary genetic determinant of differences in autodiploidization. We then used multiple laboratory and wild strains of S. cerevisiae to show that clonal populations of strains with a functional copy of SSD1 autodiploidize more frequently in evolution experiments, while knocking out this gene or replacing it with the W303 allele reduces autodiploidization propensity across all genetic backgrounds tested. These results suggest a potential strategy for modifying rates of spontaneous whole-genome duplications in laboratory evolution experiments in haploid budding yeast. They may also have relevance to other settings in which eukaryotic genome stability plays an important role, such as biomanufacturing and the treatment of pathogenic fungal diseases and cancers.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

Frequency-dependent interactions determine outcome of competition between two breast cancer cell lines

Tumors are highly dynamic ecosystems in which diverse cancer cell subpopulations compete for space and resources. These complex, often non-linear interactions govern continuous spatial and temporal changes in the size and phenotypic properties of these subpopulations. Because intra-tumoral blood flow is often chaotic, competition for resources may be a critical selection factor in progression and prognosis. Here, we quantify resource competition using 3D spheroid cultures with MDA-MB-231 and MCF-7 breast cancer cells. We hypothesized that MCF-7 cells, which primarily rely on efficient aerobic glucose metabolism, would dominate the population under normal pH and low glucose conditions; and MDA-MB-231 cells, which exhibit high levels of glycolytic metabolism, would dominate under low pH and high glucose conditions. In spheroids with single populations, MCF-7 cells exhibited equal or superior intrinsic growth rates (density-independent measure of success) and carrying capacities (density-dependent measure of success) when compared to MDA-MB-231 cells under all pH and nutrient conditions. Despite these advantages, when grown together, MCF-7 cells do not always outcompete MDA-MB-231 cells. MDA-MB-231 cells outcompete MCF-7 cells in low glucose conditions and coexistence is achieved in low pH conditions. Under all conditions, MDA-MB-231 has a stronger competitive effect (frequency-dependent interaction) on MCF-7 cells than vice-versa. This, and the inability of growth rate or carrying capacity when grown individually to predict the outcome of competition, suggests a reliance on frequency-dependent interactions and the need for competition assays. We frame these results in a game-theoretic (frequency-dependent) model of cancer cell interactions and conclude that competition assays can demonstrate critical density-independent, density-dependent and frequency-dependent interactions that likely contribute to in vivo outcomes.

The population genetics of ploidy change in unicellular fungi

Changes in ploidy are a significant type of genetic variation, describing the number of chromosome sets per cell. Ploidy evolves in natural populations, clinical populations, and lab experiments, particularly in fungi. Despite a long history of theoretical work on this topic, predicting how ploidy will evolve has proven difficult, as it is often unclear why one ploidy state outperforms another. Here, we review what is known about contemporary ploidy evolution in diverse fungal species through the lens of population genetics. As with typical genetic variants, ploidy evolution depends on the rate that new ploidy states arise by mutation, natural selection on alternative ploidy states, and random genetic drift. However, ploidy variation also has unique impacts on evolution, with the potential to alter chromosomal stability, the rate and patterns of point mutation, and the nature of selection on all loci in the genome. We discuss how ploidy evolution depends on these general and unique factors and highlight areas where additional experimental evidence is required to comprehensively explain the ploidy transitions observed in the field and the lab.

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.

Host-Induced Genome Instability Rapidly Generates Phenotypic Variation across Candida albicans Strains and Ploidy States

Candida albicans is an opportunistic fungal pathogen of humans. The ability to generate genetic variation is essential for adaptation and is a strategy that C. albicans and other fungal pathogens use to change their genome size. Stressful environments, including the host, induce C. albicans genome instability. Here, we investigated how C. albicans genetic background and ploidy state impact genome instability, both in vitro and in a host environment. We show that the host environment induces genome instability, but the magnitude depends on C. albicans genetic background. Furthermore, we show that tetraploid C. albicans is highly unstable in host environments and rapidly reduces in genome size. These reductions in genome size often resulted in reduced virulence. In contrast, diploid C. albicans displayed modest host-induced genome size changes, yet these frequently resulted in increased virulence. Such studies are essential for understanding how opportunistic pathogens respond and potentially adapt to the host environment.

Candida albicans is an opportunistic fungal pathogen of humans that is typically diploid yet has a highly labile genome tolerant of large-scale perturbations including chromosomal aneuploidy and loss-of-heterozygosity events. The ability to rapidly generate genetic variation is crucial for C. albicans to adapt to changing or stressful environments, like those encountered in the host. Genetic variation occurs via stress-induced mutagenesis or can be generated through its parasexual cycle, in which tetraploids arise via diploid mating or stress-induced mitotic defects and undergo nonmeiotic ploidy reduction. However, it remains largely unknown how genetic background contributes to C. albicans genome instability in vitro or in the host environment. Here, we tested how genetic background, ploidy, and the host environment impacts C. albicans genome stability. We found that host association induced both loss-of-heterozygosity events and genome size changes, regardless of genetic background or ploidy. However, the magnitude and types of genome changes varied across C. albicans strain background and ploidy state. We then assessed if host-induced genomic changes resulted in fitness consequences on growth rate and nonlethal virulence phenotypes and found that many host-derived isolates significantly changed relative to their parental strain. Interestingly, diploid host-associated C. albicans predominantly decreased host reproductive fitness, whereas tetraploid host-associated C. albicans increased host reproductive fitness. Together, these results are important for understanding how host-induced genomic changes in C. albicans alter its relationship with the host.

IMPORTANCECandida albicans is an opportunistic fungal pathogen of humans. The ability to generate genetic variation is essential for adaptation and is a strategy that C. albicans and other fungal pathogens use to change their genome size. Stressful environments, including the host, induce C. albicans genome instability. Here, we investigated how C. albicans genetic background and ploidy state impact genome instability, both in vitro and in a host environment. We show that the host environment induces genome instability, but the magnitude depends on C. albicans genetic background. Furthermore, we show that tetraploid C. albicans is highly unstable in host environments and rapidly reduces in genome size. These reductions in genome size often resulted in reduced virulence. In contrast, diploid C. albicans displayed modest host-induced genome size changes, yet these frequently resulted in increased virulence. Such studies are essential for understanding how opportunistic pathogens respond and potentially adapt to the host environment.

The genome-wide rate and spectrum of spontaneous mutations differ between haploid and diploid yeast

By altering the dynamics of DNA replication and repair, alternative ploidy states may experience different rates and types of new mutations, leading to divergent evolutionary outcomes. We report a direct comparison of the genome-wide spectrum of spontaneous mutations arising in haploids and diploids following a mutation-accumulation experiment in the budding yeast Saccharomyces cerevisiae Characterizing the number, types, locations, and effects of thousands of mutations revealed that haploids were more prone to single-nucleotide mutations (SNMs) and mitochondrial mutations, while larger structural changes were more common in diploids. Mutations were more likely to be detrimental in diploids, even after accounting for the large impact of structural changes, contrary to the prediction that mutations would have weaker effects, due to masking, in diploids. Haploidy is expected to reduce the opportunity for conservative DNA repair involving homologous chromosomes, increasing the insertion-deletion rate, but we found little support for this idea. Instead, haploids were more susceptible to SNMs in late-replicating genomic regions, resulting in a ploidy difference in the spectrum of substitutions. In diploids, we detect mutation rate variation among chromosomes in association with centromere location, a finding that is supported by published polymorphism data. Diploids are not simply doubled haploids; instead, our results predict that the spectrum of spontaneous mutations will substantially shape the dynamics of genome evolution in haploid and diploid populations.

Adaptive genome duplication affects patterns of molecular evolution in Saccharomyces cerevisiae

Genome duplications are important evolutionary events that impact the rate and spectrum of beneficial mutations and thus the rate of adaptation. Laboratory evolution experiments initiated with haploid Saccharomyces cerevisiae cultures repeatedly experience whole-genome duplication (WGD). We report recurrent genome duplication in 46 haploid yeast populations evolved for 4,000 generations. We find that WGD confers a fitness advantage, and this immediate fitness gain is accompanied by a shift in genomic and phenotypic evolution. The presence of ploidy-enriched targets of selection and structural variants reveals that autodiploids utilize adaptive paths inaccessible to haploids. We find that autodiploids accumulate recessive deleterious mutations, indicating an increased susceptibility for nonadaptive evolution. Finally, we report that WGD results in a reduced adaptation rate, indicating a trade-off between immediate fitness gains and long-term adaptability.

Whole genome duplications—the simultaneous doubling of each chromosome—can have a profound influence on evolution. Evidence of ancient whole genome duplications can be seen in most modern genomes. Experimental evolution, the long-term propagation of organisms under well-controlled laboratory conditions, yields valuable insight into the processes of adaptation and genome evolution. One interesting, and common, outcome of laboratory evolution experiments that start with haploid yeast populations is the emergence of diploid lineages via whole genome duplication. We show that, under our laboratory conditions, whole genome duplication provides a direct fitness benefit, and we identify several consequences of whole genome duplication on adaptation. Following whole-genome duplication, the rate of adaptation slows, the biological targets of selection change, and aneuploidies, copy-number variants and recessive lethal mutations accumulate. By studying the effect of whole genome duplication on adaptation, we can better understand how selection acts on ploidy, a fundamental biological parameter that varies considerably across life.

Increase in viability due to the accumulation of X chromosome mutations in Drosophila melanogaster males

To increase our understanding of the role of new X-chromosome mutations in adaptive evolution, single-X Drosophila melanogaster males were mated with attached-X chromosome females, allowing the male X chromosome to accumulate mutations over 28 generations. Contrary to our hypothesis that male viability would decrease over time, due to the accumulation and expression of X-linked recessive deleterious mutations in hemizygous males, viability significantly increased. This increase may be attributed to germinal selection and to new X-linked beneficial or compensatory mutations, possibly supporting the faster-X hypothesis.

Feedback between environment and traits under selection in a seasonal environment: consequences for experimental evolution

Batch cultures are frequently used in experimental evolution to study the dynamics of adaptation. Although they are generally considered to simply drive a growth rate increase, other fitness components can also be selected for. Indeed, recurrent batches form a seasonal environment where different phases repeat periodically and different traits can be under selection in the different seasons. Moreover, the system being closed, organisms may have a strong impact on the environment. Thus, the study of adaptation should take into account the environment and eco-evolutionary feedbacks. Using data from an experimental evolution on yeast Saccharomyces cerevisiae, we developed a mathematical model to understand which traits are under selection, and what is the impact of the environment for selection in a batch culture. We showed that two kinds of traits are under selection in seasonal environments: life-history traits, related to growth and mortality, but also transition traits, related to the ability to react to environmental changes. The impact of environmental conditions can be summarized by the length of the different seasons which weight selection on each trait: the longer a season is, the higher the selection on associated traits. Since phenotypes drive season length, eco-evolutionary feedbacks emerge. Our results show how evolution in successive batches can affect season lengths and strength of selection on different traits.

Evolution: Zeroing In on the Rate of Genome Doubling

Changes in genome copy number have occurred numerous times throughout the history of life, with profound evolutionary consequences. New experiments with budding yeast shed light on how frequently spontaneous genome doubling occurs within populations and the environmental conditions that favour cells with doubled genomes.

Frequent ploidy changes in growing yeast cultures

Ploidy is considered a very stable cellular characteristic. Although rare, changes in ploidy play important roles in the acquisition of long-term adaptations. Since these duplications allow the subsequent loss of individual chromosomes and accumulation of mutations, changes in ploidy can also cause genomic instability, and have been found to promote cancer. Despite the importance of the subject, measuring the rate of whole-genome duplications has proven extremely challenging. We have recently measured the rate of diploidization in yeast using long-term, in-lab experiments. We found that spontaneous diploidization occurs frequently, by two different mechanisms: endoreduplication and mating type switching. Despite its common occurrence, spontaneous diploidization is usually selected against, although it can be advantageous under some stressful conditions. Our results have implications for the understanding of evolutionary processes, as well as for the use of yeast cells in biotechnological applications.

Spontaneous Changes in Ploidy Are Common in Yeast

Changes in ploidy are relatively rare, but play important roles in the development of cancer and the acquisition of long-term adaptations. Genome duplications occur across the tree of life, and can alter the rate of adaptive evolution. Moreover, by allowing the subsequent loss of individual chromosomes and the accumulation of mutations, changes in ploidy can promote genomic instability and/or adaptation. Although many studies have been published in the last years about changes in chromosome number and their evolutionary consequences, tracking and measuring the rate of whole-genome duplications have been extremely challenging. We have systematically studied the appearance of diploid cells among haploid yeast cultures evolving for over 100 generations in different media. We find that spontaneous diploidization is a relatively common event, which is usually selected against, but under certain stressful conditions may become advantageous. Furthermore, we were able to detect and distinguish between two different mechanisms of diploidization, one that requires whole-genome duplication (endoreduplication) and a second that involves mating-type switching despite the use of heterothallic strains. Our results have important implications for our understanding of evolution and adaptation in fungal pathogens and the development of cancer, and for the use of yeast cells in biotechnological applications.

Hidden Complexity of Yeast Adaptation under Simple Evolutionary Conditions

Few studies have "quantitatively" probed how adaptive mutations result in increased fitness. Even in simple microbial evolution experiments, with full knowledge of the underlying mutations and specific growth conditions, it is challenging to determine where within a growth-saturation cycle those fitness gains occur. A common implicit assumption is that most benefits derive from an increased exponential growth rate. Here, we instead show that, in batch serial transfer experiments, adaptive mutants' fitness gains can be dominated by benefits that are accrued in one growth cycle, but not realized until the next growth cycle. For thousands of evolved clones (most with only a single mutation), we systematically varied the lengths of fermentation, respiration, and stationary phases to assess how their fitness, as measured by barcode sequencing, depends on these phases of the growth-saturation-dilution cycles. These data revealed that, whereas all adaptive lineages gained similar and modest benefits from fermentation, most of the benefits for the highest fitness mutants came instead from the time spent in respiration. From monoculture and high-resolution pairwise fitness competition experiments for a dozen of these clones, we determined that the benefits "accrued" during respiration are only largely "realized" later as a shorter duration of lag phase in the following growth cycle. These results reveal hidden complexities of the adaptive process even under ostensibly simple evolutionary conditions, in which fitness gains can accrue during time spent in a growth phase with little cell division, and reveal that the memory of those gains can be realized in the subsequent growth cycle.

Ploidy tug-of-war: Evolutionary and genetic environments influence the rate of ploidy drive in a human fungal pathogen

Variation in baseline ploidy is seen throughout the tree of life, yet the factors that determine why one ploidy level is maintained over another remain poorly understood. Experimental evolution studies using asexual fungal microbes with manipulated ploidy levels intriguingly reveals a propensity to return to the historical baseline ploidy, a phenomenon that we term "ploidy drive." We evolved haploid, diploid, and polyploid strains of the human fungal pathogen Candida albicans under three different nutrient limitation environments to test whether these conditions, hypothesized to select for low ploidy levels, could counteract ploidy drive. Strains generally maintained or acquired smaller genome sizes (measured as total nuclear DNA through flow cytometry) in minimal medium and under phosphorus depletion compared to in a complete medium, while mostly maintained or acquired increased genome sizes under nitrogen depletion. Improvements in fitness often ran counter to changes in genome size; in a number of scenarios lines that maintained their original genome size often increased in fitness more than lines that converged toward diploidy (the baseline ploidy of C. albicans). Combined, this work demonstrates a role for both the environment and genotype in determination of the rate of ploidy drive, and highlights questions that remain about the force(s) that cause genome size variation.

Adaptation to High Ethanol Reveals Complex Evolutionary Pathways

Tolerance to high levels of ethanol is an ecologically and industrially relevant phenotype of microbes, but the molecular mechanisms underlying this complex trait remain largely unknown. Here, we use long-term experimental evolution of isogenic yeast populations of different initial ploidy to study adaptation to increasing levels of ethanol. Whole-genome sequencing of more than 30 evolved populations and over 100 adapted clones isolated throughout this two-year evolution experiment revealed how a complex interplay of de novo single nucleotide mutations, copy number variation, ploidy changes, mutator phenotypes, and clonal interference led to a significant increase in ethanol tolerance. Although the specific mutations differ between different evolved lineages, application of a novel computational pipeline, PheNetic, revealed that many mutations target functional modules involved in stress response, cell cycle regulation, DNA repair and respiration. Measuring the fitness effects of selected mutations introduced in non-evolved ethanol-sensitive cells revealed several adaptive mutations that had previously not been implicated in ethanol tolerance, including mutations in PRT1, VPS70 and MEX67. Interestingly, variation in VPS70 was recently identified as a QTL for ethanol tolerance in an industrial bio-ethanol strain. Taken together, our results show how, in contrast to adaptation to some other stresses, adaptation to a continuous complex and severe stress involves interplay of different evolutionary mechanisms. In addition, our study reveals functional modules involved in ethanol resistance and identifies several mutations that could help to improve the ethanol tolerance of industrial yeasts.

Organisms can evolve resistance to specific stress factors, which allows them to thrive in environments where non-adapted organisms fail to grow. However, the molecular mechanisms that underlie adaptation to complex stress factors that interfere with basic cellular processes are poorly understood. In this study, we reveal how yeast populations adapt to high ethanol concentrations, an ecologically and industrially relevant stress that is still poorly understood. We exposed six independent populations of genetically identical yeast cells to gradually increasing ethanol levels, and we monitored the changes in their DNA sequence over a two-year period. Together with novel computational analyses, we could identify the mutational dynamics and molecular mechanisms underlying increased ethanol resistance. Our results show how adaptation to high ethanol is complex and can be reached through different mutational pathways. Together, our study offers a detailed picture of how populations adapt to a complex continuous stress and identifies several mutations that increase ethanol resistance, which opens new routes to obtain superior biofuel yeast strains.

Genomic investigations of evolutionary dynamics and epistasis in microbial evolution experiments

Microbial evolution experiments enable us to watch adaptation in real time, and to quantify the repeatability and predictability of evolution by comparing identical replicate populations. Further, we can resurrect ancestral types to examine changes over evolutionary time. Until recently, experimental evolution has been limited to measuring phenotypic changes, or to tracking a few genetic markers over time. However, recent advances in sequencing technology now make it possible to extensively sequence clones or whole-population samples from microbial evolution experiments. Here, we review recent work exploiting these techniques to understand the genomic basis of evolutionary change in experimental systems. We first focus on studies that analyze the dynamics of genome evolution in microbial systems. We then survey work that uses observations of sequence evolution to infer aspects of the underlying fitness landscape, concentrating on the epistatic interactions between mutations and the constraints these interactions impose on adaptation.

Shift and adapt: the costs and benefits of karyotype variations

Variation is the spice of life or, in the case of evolution, variation is the necessary material on which selection can act to enable adaptation. Karyotypic variation in ploidy (the number of homologous chromosome sets) and aneuploidy (imbalance in the number of chromosomes) are fundamentally different than other types of genomic variants. Karyotypic variation emerges through different molecular mechanisms than other mutational events, and unlike mutations that alter the genome at the base pair level, rapid reversion to the wild type chromosome number is often possible. Although karyotypic variation has long been noted and discussed by biologists, interest in the importance of karyotypic variants in evolutionary processes has spiked in recent years, and much remains to be discovered about how karyotypic variants are produced and subsequently selected.

Parasexual Ploidy Reduction Drives Population Heterogeneity Through Random and Transient Aneuploidy in Candida albicans

The opportunistic pathogen Candida albicans has a large repertoire of mechanisms to generate genetic and phenotypic diversity despite the lack of meiosis in its life cycle. Its parasexual cycle enables shifts in ploidy, which in turn facilitate recombination, aneuploidy, and homozygosis of whole chromosomes to fuel rapid adaptation. Here we show that the tetraploid state potentiates ploidy variation and drives population heterogeneity. In tetraploids, the rate of losing a single heterozygous marker [loss of heterozygosity (LOH)] is elevated ∼30-fold higher than the rate in diploid cells. Furthermore, isolates recovered after selection for LOH of one, two, or three markers were highly aneuploid, with a broad range of karyotypes including strains with a combination of di-, tri-, and tetrasomic chromosomes. We followed the ploidy trajectories for these tetraploid- and aneuploid-derived isolates, using a combination of flow cytometry and double-digestion restriction-site-associated DNA analyzed with next-generation sequencing. Isolates derived from either tetraploid or aneuploid isolates predominately resolved to a stable euploid state. The majority of isolates reduced to the conventional diploid state; however, stable triploid and tetraploid states were observed in ∼30% of the isolates. Notably, aneuploid isolates were more transient than tetraploid isolates, resolving to a euploid state within a few passages. Furthermore, the likelihood that a particular isolate will resolve to the same ploidy state in replicate evolution experiments is only ∼50%, supporting the idea that the chromosome loss process of the parasexual cycle is random and does not follow trajectories involving specific combinations of chromosomes. Together, our results indicate that tetraploid progenitors can produce populations of progeny cells with a high degree of genomic diversity, from altered ploidy to homozygosis, providing an excellent source of genetic variation upon which selection can act.

Efficient purging of deleterious mutations in plants with haploid selfing

In diploid organisms, selfing reduces the efficiency of selection in removing deleterious mutations from a population. This need not be the case for all organisms. Some plants, for example, undergo an extreme form of selfing known as intragametophytic selfing, which immediately exposes all recessive deleterious mutations in a parental genome to selective purging. Here, we ask how effectively deleterious mutations are removed from such plants. Specifically, we study the extent to which deleterious mutations accumulate in a predominantly selfing and a predominantly outcrossing pair of moss species, using genome-wide transcriptome data. We find that the selfing species purge significantly more nonsynonymous mutations, as well as a greater proportion of radical amino acid changes which alter physicochemical properties of amino acids. Moreover, their purging of deleterious mutation is especially strong in conserved regions of protein-coding genes. Our observations show that selfing need not impede but can even accelerate the removal of deleterious mutations, and do so on a genome-wide scale.

Loss-of-heterozygosity facilitates passage through Haldane's sieve for Saccharomyces cerevisiae undergoing adaptation

Haldane's sieve posits that the majority of beneficial mutations that contribute to adaptation should be dominant, as these are the mutations most likely to establish and spread when rare. It has been argued, however, that if the dominance of mutations in their current and previous environments are correlated, Haldane's sieve could be eliminated. We constructed heterozygous lines of Saccharomyces cerevisiae containing single adaptive mutations obtained during exposure to the fungicide nystatin. Here we show that no clear dominance relationship exists across environments: mutations exhibited a range of dominance levels in a rich medium, yet were exclusively recessive under nystatin stress. Surprisingly, heterozygous replicates exhibited variable-onset rapid growth when exposed to nystatin. Targeted Sanger sequencing demonstrated that loss-of-heterozygosity (LOH) accounted for these growth patterns. Our experiments demonstrate that recessive beneficial mutations can avoid Haldane's sieve in clonal organisms through rapid LOH and thus contribute to rapid evolutionary adaptation.

Mapping small effect mutations in Saccharomyces cerevisiae: impacts of experimental design and mutational properties

Genetic variants identified by mapping are biased toward large phenotypic effects because of methodologic challenges for detecting genetic variants with small phenotypic effects. Recently, bulk segregant analysis combined with next-generation sequencing (BSA-seq) was shown to be a powerful and cost-effective way to map small effect variants in natural populations. Here, we examine the power of BSA-seq for efficiently mapping small effect mutations isolated from a mutagenesis screen. Specifically, we determined the impact of segregant population size, intensity of phenotypic selection to collect segregants, number of mitotic generations between meiosis and sequencing, and average sequencing depth on power for mapping mutations with a range of effects on the phenotypic mean and standard deviation as well as relative fitness. We then used BSA-seq to map the mutations responsible for three ethyl methanesulfonate−induced mutant phenotypes in Saccharomyces cerevisiae. These mutants display small quantitative variation in the mean expression of a fluorescent reporter gene (−3%, +7%, and +10%). Using a genetic background with increased meiosis rate, a reliable mating type marker, and fluorescence-activated cell sorting to efficiently score large segregating populations and isolate cells with extreme phenotypes, we successfully mapped and functionally confirmed a single point mutation responsible for the mutant phenotype in all three cases. Our simulations and experimental data show that the effects of a causative site not only on the mean phenotype, but also on its standard deviation and relative fitness should be considered when mapping genetic variants in microorganisms such as yeast that require population growth steps for BSA-seq.

An experimental test on the probability of extinction of new genetic variants

In 1927, J.B.S. Haldane reasoned that the probability of fixation of new beneficial alleles is twice their fitness effect. This result, later generalized by M. Kimura, has since become the cornerstone of modern population genetics. There is no experimental test of Haldane's insight that new beneficial alleles are lost with high probability. Here we demonstrate that extinction rates decrease with increasing initial numbers of beneficial alleles, as expected, by performing invasion experiments with inbred lines of the nematode Caenorhabditis elegans. We further show that the extinction rates of deleterious alleles are higher than those of beneficial alleles, also as expected. Interestingly, we also find that for these inbred lines, when at intermediate frequencies, the fate of invaders might not result in their ultimate fixation or loss but on their maintenance. Our study confirms the key results from classical population genetics and highlights that the nature of adaptation can be complex.

Molecular specificity, convergence and constraint shape adaptive evolution in nutrient-poor environments

One of the central goals of evolutionary biology is to explain and predict the molecular basis of adaptive evolution. We studied the evolution of genetic networks in Saccharomyces cerevisiae (budding yeast) populations propagated for more than 200 generations in different nitrogen-limiting conditions. We find that rapid adaptive evolution in nitrogen-poor environments is dominated by the de novo generation and selection of copy number variants (CNVs), a large fraction of which contain genes encoding specific nitrogen transporters including PUT4, DUR3 and DAL4. The large fitness increases associated with these alleles limits the genetic heterogeneity of adapting populations even in environments with multiple nitrogen sources. Complete identification of acquired point mutations, in individual lineages and entire populations, identified heterogeneity at the level of genetic loci but common themes at the level of functional modules, including genes controlling phosphatidylinositol-3-phosphate metabolism and vacuole biogenesis. Adaptive strategies shared with other nutrient-limited environments point to selection of genetic variation in the TORC1 and Ras/PKA signaling pathways as a general mechanism underlying improved growth in nutrient-limited environments. Within a single population we observed the repeated independent selection of a multi-locus genotype, comprised of the functionally related genes GAT1, MEP2 and LST4. By studying the fitness of individual alleles, and their combination, as well as the evolutionary history of the evolving population, we find that the order in which these mutations are acquired is constrained by epistasis. The identification of repeatedly selected variation at functionally related loci that interact epistatically suggests that gene network polymorphisms (GNPs) may be a frequent outcome of adaptive evolution. Our results provide insight into the mechanistic basis by which cells adapt to nutrient-limited environments and suggest that knowledge of the selective environment and the regulatory mechanisms important for growth and survival in that environment greatly increase the predictability of adaptive evolution.

We studied adaptive evolution in different nitrogen-limited environments using long-term selection of asexually reproducing Saccharomyces cerevisiae populations in chemostats. Using next generation sequencing and DNA microarrays, we identified all acquired genetic variation associated with increased fitness, in both individual lineages and entire populations. We find that amplification alleles that include nutrient transporter genes specific to the molecular form of the nitrogen present in the environment are a common mechanism underlying increased fitness. In addition, we identified a general strategy for adaptation to nitrogen-limited environments that entails remodeling of phospholipid biogenesis required for producing important cellular components including vacuoles and autophagosomes. More general strategies for adaptation to nutrient-limited environments point to a role for re-wiring of signaling pathways that coordinate cell growth with nutrient availability. We reconstructed the evolutionary dynamics of a population evolving in ammonium-limited conditions and find that a multi-locus genotype is repeatedly selected within the population and constrained by epistasis. We propose that this genotype constitutes a “gene network polymorphism (GNP),” which may be a common outcome of adaptive evolution. Our study suggests that when the selective pressure is understood the molecular basis of adaptive evolution in large microbial populations may be predicted with reasonable precision.

Rapid evolution of quantitative traits: theoretical perspectives

An increasing number of studies demonstrate phenotypic and genetic changes in natural populations that are subject to climate change, and there is hope that some of these changes will contribute to avoiding species extinctions ('evolutionary rescue'). Here, we review theoretical models of rapid evolution in quantitative traits that can shed light on the potential for adaptation to a changing climate. Our focus is on quantitative-genetic models with selection for a moving phenotypic optimum. We point out that there is no one-to-one relationship between the rate of adaptation and population survival, because the former depends on relative fitness and the latter on absolute fitness. Nevertheless, previous estimates that sustainable rates of genetically based change usually do not exceed 0.1 haldanes (i.e., phenotypic standard deviations per generation) are probably correct. Survival can be greatly facilitated by phenotypic plasticity, and heritable variation in plasticity can further speed up genetic evolution. Multivariate selection and genetic correlations are frequently assumed to constrain adaptation, but this is not necessarily the case and depends on the geometric relationship between the fitness landscape and the structure of genetic variation. Similar conclusions hold for adaptation to shifting spatial gradients. Recent models of adaptation in multispecies communities indicate that the potential for rapid evolution is strongly influenced by interspecific competition.

An evaluation of high-throughput approaches to QTL mapping in Saccharomyces cerevisiae

Dissecting the molecular basis of quantitative traits is a significant challenge and is essential for understanding complex diseases. Even in model organisms, precisely determining causative genes and their interactions has remained elusive, due in part to difficulty in narrowing intervals to single genes and in detecting epistasis or linked quantitative trait loci. These difficulties are exacerbated by limitations in experimental design, such as low numbers of analyzed individuals or of polymorphisms between parental genomes. We address these challenges by applying three independent high-throughput approaches for QTL mapping to map the genetic variants underlying 11 phenotypes in two genetically distant Saccharomyces cerevisiae strains, namely (1) individual analysis of >700 meiotic segregants, (2) bulk segregant analysis, and (3) reciprocal hemizygosity scanning, a new genome-wide method that we developed. We reveal differences in the performance of each approach and, by combining them, identify eight polymorphic genes that affect eight different phenotypes: colony shape, flocculation, growth on two nonfermentable carbon sources, and resistance to two drugs, salt, and high temperature. Our results demonstrate the power of individual segregant analysis to dissect QTL and address the underestimated contribution of interactions between variants. We also reveal confounding factors like mutations and aneuploidy in pooled approaches, providing valuable lessons for future designs of complex trait mapping studies.

Phenotypic and genotypic convergences are influenced by historical contingency and environment in yeast

Different organisms have independently and recurrently evolved similar phenotypic traits at different points throughout history. This phenotypic convergence may be caused by genotypic convergence and in addition, constrained by historical contingency. To investigate how convergence may be driven by selection in a particular environment and constrained by history, we analyzed nine life-history traits and four metabolic traits during an experimental evolution of six yeast strains in four different environments. In each of the environments, the population converged toward a different multivariate phenotype. However, the evolution of most traits, including fitness components, was constrained by history. Phenotypic convergence was partly associated with the selection of mutations in genes involved in the same pathway. By further investigating the convergence in one gene, BMH1, mutated in 20% of the evolved populations, we show that both the history and the environment influenced the types of mutations (missense/nonsense), their location within the gene itself, as well as their effects on multiple traits. However, these effects could not be easily predicted from ancestors' phylogeny or past selection. Combined, our data highlight the role of pleiotropy and epistasis in shaping a rugged fitness landscape.

Parasexuality and ploidy change in Candida tropicalis

Candida species exhibit a variety of ploidy states and modes of sexual reproduction. Most species possess the requisite genes for sexual reproduction, recombination, and meiosis, yet only a few have been reported to undergo a complete sexual cycle including mating and sporulation. Candida albicans, the most studied Candida species and a prevalent human fungal pathogen, completes its sexual cycle via a parasexual process of concerted chromosome loss rather than a conventional meiosis. In this study, we examine ploidy changes in Candida tropicalis, a closely related species to C. albicans that was recently revealed to undergo sexual mating. C. tropicalis diploid cells mate to form tetraploid cells, and we show that these can be induced to undergo chromosome loss to regenerate diploid forms by growth on sorbose medium. The diploid products are themselves mating competent, thereby establishing a parasexual cycle in this species for the first time. Extended incubation (>120 generations) of C. tropicalis tetraploid cells under rich culture conditions also resulted in instability of the tetraploid form and a gradual reduction in ploidy back to the diploid state. The fitness levels of C. tropicalis diploid and tetraploid cells were compared, and diploid cells exhibited increased fitness relative to tetraploid cells in vitro, despite diploid and tetraploid cells having similar doubling times. Collectively, these experiments demonstrate distinct pathways by which a parasexual cycle can occur in C. tropicalis and indicate that nonmeiotic mechanisms drive ploidy changes in this prevalent human pathogen.

Ancient evolutionary trade-offs between yeast ploidy states

The number of chromosome sets contained within the nucleus of eukaryotic organisms is a fundamental yet evolutionarily poorly characterized genetic variable of life. Here, we mapped the impact of ploidy on the mitotic fitness of baker's yeast and its never domesticated relative Saccharomyces paradoxus across wide swaths of their natural genotypic and phenotypic space. Surprisingly, environment-specific influences of ploidy on reproduction were found to be the rule rather than the exception. These ploidy–environment interactions were well conserved across the 2 billion generations separating the two species, suggesting that they are the products of strong selection. Previous hypotheses of generalizable advantages of haploidy or diploidy in ecological contexts imposing nutrient restriction, toxin exposure, and elevated mutational loads were rejected in favor of more fine-grained models of the interplay between ecology and ploidy. On a molecular level, cell size and mating type locus composition had equal, but limited, explanatory power, each explaining 12.5%–17% of ploidy–environment interactions. The mechanism of the cell size–based superior reproductive efficiency of haploids during Li⁺ exposure was traced to the Li⁺ exporter ENA. Removal of the Ena transporters, forcing dependence on the Nha1 extrusion system, completely altered the effects of ploidy on Li⁺ tolerance and evoked a strong diploid superiority, demonstrating how genetic variation at a single locus can completely reverse the relative merits of haploidy and diploidy. Taken together, our findings unmasked a dynamic interplay between ploidy and ecology that was of unpredicted evolutionary importance and had multiple molecular roots.

Organisms vary in the number of chromosome sets contained within the nucleus of each cell, but neither the reasons nor the consequences of this variation are well understood. We designed yeasts that differed in the number of chromosome sets but were otherwise identical and mapped the consequences of such ploidy variations during exposure to a large palette of environments. Contrary to commonly held assumptions, we found ploidy effects on the mitotic reproductive capacity of yeast to be the rule rather than the exception and to be highly evolutionarily conserved. Furthermore, our data rejected previously contemplated hypotheses of generalizable advantages of haploidy or diploidy when cells face nutrient starvation or are exposed to toxins or increased mutation rates. We also mapped the molecular processes mediating ploidy–environment interactions, showing that cell size and mating type locus composition had equal explanatory power. Finally we show that ploidy effects can be mechanistically very subtle, as a designed shift from one plasma membrane Li⁺ transporter to another completely altered the relative merits of having one or two chromosome sets when exposed to high Li⁺ concentrations. This complex and dynamic interplay between the number of chromosomes sets and the fluctuating environment must be taken into account when considering organismal form and behavior.

Mutational effects depend on ploidy level: all else is not equal

Ploidy is predicted to influence adaptation directly, yet whether single mutations behave the same in different ploidy backgrounds has not been well studied. It has often been assumed theoretically that aside from dominance, selective parameters do not differ between cells of varying ploidy. Using the budding yeast Saccharomyces cerevisiae, I compared the effect size of 20 adaptive mutations in haploids and homozygous diploids and found, surprisingly, that the same mutations often had a much larger effect in haploids than homozygous diploids. This empirical result demonstrates that it cannot be assumed that mutations will have the same effect in haploids and homozygous diploids.

Parallel genetic changes and nonparallel gene-environment interactions characterize the evolution of drug resistance in yeast

Beneficial mutations are required for adaptation to novel environments, yet the range of mutational pathways that are available to a population has been poorly characterized, particularly in eukaryotes. We assessed the genetic changes of the first mutations acquired during adaptation to a novel environment (exposure to the fungicide, nystatin) in 35 haploid lines of Saccharomyces cerevisiae. Through whole-genome resequencing we found that the genomic scope for adaptation was narrow; all adapted lines acquired a mutation in one of four late-acting genes in the ergosterol biosynthesis pathway, with very few other mutations found. Lines that acquired different ergosterol mutations in the same gene exhibited very similar tolerance to nystatin. All lines were found to have a cost relative to wild type in an unstressful environment; the level of this cost was also strongly correlated with the ergosterol gene bearing the mutation. Interestingly, we uncovered both positive and negative effects on tolerance to other harsh environments for mutations in the different ergosterol genes, indicating that these beneficial mutations have effects that differ in sign among environmental challenges. These results demonstrate that although the genomic target was narrow, different adaptive mutations can lead populations down different evolutionary pathways, with respect to their ability to tolerate (or succumb to) other environmental challenges.