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

Genome restructuring and lineage diversification of Cryptococcus neoformans during chronic infection of human hosts

Classified as a critical public health threat by the World Health Organization, Cryptococcus neo-formans infections with significant morbidity and mortality. Reports of cryptococcosis persistence, relapse, and reinfection date back to the 1950s, yet the factors driving chronic infections remain poorly understood. A major challenge is the scarcity of serial patient specimens and detailed medical records to study the simultaneous evolution of the pathogen and host health status. This study provides the first genomic and phenotypic analysis of in-host evolution of C. neoformans during chronic infections lasting over a year in six immunocompromised patients. We find fungal genome evolution during persistent infection is characterized by large-scale genome restructuring and increasing genomic heterogeneity. Phenotypic changes show diversification in virulence traits and antifungal susceptibility. Genotypically and phenotypically distinct sub-lineages arise and co-persist within the same tissues, consistent with a model of diversifying selection and niche partitioning in the complex environment of human hosts.

Quantifying Competitive Fitness in Yeast with High-Throughput Fluorescence Microscopy Imaging

Competitive fitness is a fundamental concept in evolutionary biology that captures the ability of organisms to survive, reproduce, and compete for resources in their environment. Competitive fitness is typically assessed in the lab by growing two or more competitors together and measuring the frequency of each at multiple time points. Traditional microbial competitive fitness assays are labor intensive and involve plating on solid medium and counting colonies. Here, we describe a method to quantitatively measure competitive fitness using fluorescence microscopic imaging and machine-learning-enabled image analysis to directly count the number of cells from each competitor in the mixed population. This high-throughput, primarily automated, and efficient process gives accurate and reproducible results for competitive fitness. Here, we describe the entire process, from sample preparation through microscopy to quantification, and provide instructions and scripts for the image analysis, fitness calculations, and sample data visualizations. © 2025 The Author(s). Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Sample preparation Basic Protocol 2: Photographing fluorescing and non-fluorescing cells using an EVOS microscope Basic Protocol 3: Counting fluorescing and non-fluorescing cells with Orbit Image Analysis Basic Protocol 4: Getting the average cell counts per well and changing the file names Basic Protocol 5: Calculating competitive fitness using R.

Massively parallel experimental interrogation of natural variants in ancient signaling pathways reveals both purifying selection and local adaptation

The nature of standing genetic variation remains a central debate in population genetics, with differing perspectives on whether common variants are mostly neutral or have functional effects. We address this question by directly mapping the fitness effects of over 9,000 natural variants in the Ras/PKA and TOR/Sch9 pathways-key regulators of cell proliferation in eukaryotes-across four conditions in Saccharomyces cerevisiae. While many variants are neutral in our assay, on the order of 3,500 exhibited significant fitness effects. These non-neutral variants tend to be missense and affect conserved, more densely packed, and less solvent-exposed protein regions. They are also typically younger, occur at lower frequencies, and more often found in heterozygous states, suggesting they are subject to purifying selection. A substantial fraction of non-neutral variants showing strong fitness effects in our experiments, however, is present at high frequencies in the population. These variants show signs of local adaptation as they tend to be found specifically in domesticated strains adapted to human-made environments. Our findings support the view that while common variants are often neutral, a significant proportion have adaptive functional consequences and are driven into the population by local positive selection. This study highlights the potential to explore the functional effects of natural genetic variation on a genome scale with quantitative fitness measurements in the laboratory, bridging the gap between population genetics and functional genomics to understand evolutionary dynamics in the wild.

Accumulating waves of random mutations before fixation

Mutations provide variation for evolution to emerge. A quantitative analysis of how mutations arising in single individuals expand and possibly fixate in a population is essential for studying evolutionary processes. While it is intuitive to expect that a continuous influx of mutations will lead to a continuous flow of mutations fixating in a stable constant population, joint fixation of multiple mutations occur frequently in stochastic simulations even under neutral selection. We quantitatively measure and analyze the distribution of joint fixation events of neutral mutations in constant populations and discussed the connection with previous results. We propose a new concept, the mutation "waves," where multiple mutations reach given frequencies simultaneously. We show that all but the lowest frequencies of the variant allele frequency distribution are dominated by single mutation "waves," which approximately follow an exponential distribution in terms of size. Consequently, large swaths of empty frequencies are observed in the variant allele frequency distributions, with a few frequencies having numbers of mutations far in excess of the expected average values over multiple realizations. We quantify the amount of time each frequency is empty of mutations and further show that the discrete mutation waves average out to a continuous distribution named as the wave frequency distribution, the shape of which is predictable based on few model parameters.

Effects of periodic bottlenecks on the dynamics of adaptive evolution in microbial populations

Population bottlenecks can impact the rate of adaptation in evolving populations. On the one hand, each bottleneck reduces the genetic variation that fuels adaptation. On the other hand, each founder that survives a bottleneck can undergo more generations and leave more descendants in a resource-limited environment, which allows surviving beneficial mutations to spread more quickly. A theoretical model predicted that the rate of fitness gains should be maximized using ~8-fold dilutions. Here we investigate the impact of repeated bottlenecks on the dynamics of adaptation using numerical simulations and experimental populations of Escherichia coli. Our simulations confirm the model's prediction when populations evolve in a regime where beneficial mutations are rare and waiting times between successful mutations are long. However, more extreme dilutions maximize fitness gains in simulations when beneficial mutations are common and clonal interference prevents most of them from fixing. To examine these predictions, we propagated 48 E. coli populations with 2-, 8-, 100-, and 1000-fold dilutions for 150 days. Adaptation began earlier and fitness gains were greater with 100- and 1000-fold dilutions than with 8-fold dilutions, consistent with the simulations when beneficial mutations are common. However, the selection pressures in the 2-fold treatment were qualitatively different from the other treatments, violating a critical assumption of the model and simulations. Thus, varying the dilution factor during periodic bottlenecks can have multiple effects on the dynamics of adaptation caused by differential losses of diversity, different numbers of generations, and altered selection.

Breeding and Selecting Corals Resilient to Global Warming

Selective breeding of resilient organisms is an emerging topic in marine conservation. It can help us predict how species will adapt in the future and how we can help restore struggling populations effectively in the present. Scleractinian corals represent a potential tractable model system given their widescale phenotypic plasticity across fitness-related traits and a reproductive life history based on mass synchronized spawning. Here, I explore the justification for breeding in corals, identify underutilized pathways of acclimation, and highlight avenues for quantitative targeted breeding from the coral host and symbiont perspective. Specifically, the facilitation of enhanced heat tolerance by targeted breeding of plasticity mechanisms is underutilized. Evidence from theoretical genetics identifies potential pitfalls, including inattention to physical and genetic characteristics of the receiving environment. Three criteria for breeding emerge from this synthesis: selection from warm, variable reefs that have survived disturbance. This information will be essential to protect what we have and restore what we can.

A dynamical limit to evolutionary adaptation

Natural selection makes evolutionary adaptation possible even if the overwhelming majority of new mutations are deleterious. However, in rapidly evolving populations where numerous linked mutations occur and segregate simultaneously, clonal interference and genetic hitchhiking can limit the efficiency of selection, allowing deleterious mutations to accumulate over time. This can in principle overwhelm the fitness increases provided by beneficial mutations, leading to an overall fitness decline. Here, we analyze the conditions under which evolution will tend to drive populations to higher versus lower fitness. Our analysis focuses on quantifying the boundary between these two regimes, as a function of parameters such as population size, mutation rates, and selection pressures. This boundary represents a state in which adaptation is precisely balanced by Muller's ratchet, and we show that it can be characterized by rapid molecular evolution without any net fitness change. Finally, we consider the implications of global fitness-mediated epistasis and find that under some circumstances, this can drive populations toward the boundary state, which can thus represent a long-term evolutionary attractor.

Species and condition shape the mutational spectrum in experimentally evolved biofilms

IMPORTANCE

Biofilm formation is a vital factor for the survival and adaptation of bacteria in diverse environmental niches. Experimental evolution combined with the advancement of whole-population genome sequencing provides us a powerful tool to understand the genomic dynamic of evolutionary adaptation to different environments, such as during biofilm development. Previous studies described the genetic and phenotypic changes of selected clones from experimentally evolved Bacillus thuringiensis and Bacillus subtilis that were adapted under abiotic and biotic biofilm conditions. However, the full understanding of the dynamic evolutionary landscapes was lacking. Furthermore, the differences and similarities of adaptive mechanisms in B. thuringiensis and B. subtilis were not identified. To overcome these limitations, we performed longitudinal whole-population genome sequencing to study the underlying genetic dynamics at high resolution. Our study provides the first comprehensive mutational landscape of two bacterial species' biofilms that is adapted to an abiotic and biotic surface.

Evolution of haploid and diploid populations reveals common, strong, and variable pleiotropic effects in non-home environments

Adaptation is driven by the selection for beneficial mutations that provide a fitness advantage in the specific environment in which a population is evolving. However, environments are rarely constant or predictable. When an organism well adapted to one environment finds itself in another, pleiotropic effects of mutations that made it well adapted to its former environment will affect its success. To better understand such pleiotropic effects, we evolved both haploid and diploid barcoded budding yeast populations in multiple environments, isolated adaptive clones, and then determined the fitness effects of adaptive mutations in 'non-home' environments in which they were not selected. We find that pleiotropy is common, with most adaptive evolved lineages showing fitness effects in non-home environments. Consistent with other studies, we find that these pleiotropic effects are unpredictable: they are beneficial in some environments and deleterious in others. However, we do find that lineages with adaptive mutations in the same genes tend to show similar pleiotropic effects. We also find that ploidy influences the observed adaptive mutational spectra in a condition-specific fashion. In some conditions, haploids and diploids are selected with adaptive mutations in identical genes, while in others they accumulate mutations in almost completely disjoint sets of genes.

A dynamical limit to evolutionary adaptation

Natural selection makes evolutionary adaptation possible even if the overwhelming majority of new mutations are deleterious. However, in rapidly evolving populations where numerous linked mutations occur and segregate simultaneously, clonal interference and genetic hitchhiking can limit the efficiency of selection, allowing deleterious mutations to accumulate over time. This can in principle overwhelm the fitness increases provided by beneficial mutations, leading to an overall fitness decline. Here, we analyze the conditions under which evolution will tend to drive populations to higher versus lower fitness. Our analysis focuses on quantifying the boundary between these two regimes, as a function of parameters such as population size, mutation rates, and selection pressures. This boundary represents a state in which adaptation is precisely balanced by Muller's ratchet, and we show that it can be characterized by rapid molecular evolution without any net fitness change. Finally, we consider the implications of global fitness-mediated epistasis, and find that under some circumstances this can drive populations towards the boundary state, which can thus represent a long-term evolutionary attractor.

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.

Fit-Seq2.0: An Improved Software for High-Throughput Fitness Measurements Using Pooled Competition Assays

The fitness of a genotype is defined as its lifetime reproductive success, with fitness itself being a composite trait likely dependent on many underlying phenotypes. Measuring fitness is important for understanding how alteration of different cellular components affects a cell's ability to reproduce. Here, we describe an improved approach, implemented in Python, for estimating fitness in high throughput via pooled competition assays.

Extreme Sensitivity of Fitness to Environmental Conditions: Lessons from #1BigBatch

The phrase "survival of the fittest" has become an iconic descriptor of how natural selection works. And yet, precisely measuring fitness, even for single-celled microbial populations growing in controlled laboratory conditions, remains a challenge. While numerous methods exist to perform these measurements, including recently developed methods utilizing DNA barcodes, all methods are limited in their precision to differentiate strains with small fitness differences. In this study, we rule out some major sources of imprecision, but still find that fitness measurements vary substantially from replicate to replicate. Our data suggest that very subtle and difficult to avoid environmental differences between replicates create systematic variation across fitness measurements. We conclude by discussing how fitness measurements should be interpreted given their extreme environment dependence. This work was inspired by the scientific community who followed us and gave us tips as we live tweeted a high-replicate fitness measurement experiment at #1BigBatch.

Best Practices in Microbial Experimental Evolution: Using Reporters and Long-Read Sequencing to Identify Copy Number Variation in Experimental Evolution

Copy number variants (CNVs), comprising gene amplifications and deletions, are a pervasive class of heritable variation. CNVs play a key role in rapid adaptation in both natural, and experimental, evolution. However, despite the advent of new DNA sequencing technologies, detection and quantification of CNVs in heterogeneous populations has remained challenging. Here, we summarize recent advances in the use of CNV reporters that provide a facile means of quantifying de novo CNVs at a specific locus in the genome, and nanopore sequencing, for resolving the often complex structures of CNVs. We provide guidance for the engineering and analysis of CNV reporters and practical guidelines for single-cell analysis of CNVs using flow cytometry. We summarize recent advances in nanopore sequencing, discuss the utility of this technology, and provide guidance for the bioinformatic analysis of these data to define the molecular structure of CNVs. The combination of reporter systems for tracking and isolating CNV lineages and long-read DNA sequencing for characterizing CNV structures enables unprecedented resolution of the mechanisms by which CNVs are generated and their evolutionary dynamics.

Hotspot of de novo telomere addition stabilizes linear amplicons in yeast grown in sulfate-limiting conditions

Evolution is driven by the accumulation of competing mutations that influence survival. A broad form of genetic variation is the amplification or deletion of DNA (≥50 bp) referred to as copy number variation (CNV). In humans, CNV may be inconsequential, contribute to minor phenotypic differences, or cause conditions such as birth defects, neurodevelopmental disorders, and cancers. To identify mechanisms that drive CNV, we monitored the experimental evolution of Saccharomyces cerevisiae populations grown under sulfate-limiting conditions. Cells with increased copy number of the gene SUL1, which encodes a primary sulfate transporter, exhibit a fitness advantage. Previously, we reported interstitial inverted triplications of SUL1 as the dominant rearrangement in a haploid population. Here, in a diploid population, we find instead that small linear fragments containing SUL1 form and are sustained over several generations. Many of the linear fragments are stabilized by de novo telomere addition within a telomere-like sequence near SUL1 (within the SNF5 gene). Using an assay that monitors telomerase action following an induced chromosome break, we show that this region acts as a hotspot of de novo telomere addition and that required sequences map to a region of <250 base pairs. Consistent with previous work showing that association of the telomere-binding protein Cdc13 with internal sequences stimulates telomerase recruitment, mutation of a four-nucleotide motif predicted to associate with Cdc13 abolishes de novo telomere addition. Our study suggests that internal telomere-like sequences that stimulate de novo telomere addition can contribute to adaptation by promoting genomic plasticity.

An improved algorithm for inferring mutational parameters from bar-seq evolution experiments

BACKGROUND

Genetic barcoding provides a high-throughput way to simultaneously track the frequencies of large numbers of competing and evolving microbial lineages. However making inferences about the nature of the evolution that is taking place remains a difficult task.

RESULTS

Here we describe an algorithm for the inference of fitness effects and establishment times of beneficial mutations from barcode sequencing data, which builds upon a Bayesian inference method by enforcing self-consistency between the population mean fitness and the individual effects of mutations within lineages. By testing our inference method on a simulation of 40,000 barcoded lineages evolving in serial batch culture, we find that this new method outperforms its predecessor, identifying more adaptive mutations and more accurately inferring their mutational parameters.

CONCLUSION

Our new algorithm is particularly suited to inference of mutational parameters when read depth is low. We have made Python code for our serial dilution evolution simulations, as well as both the old and new inference methods, available on GitHub ( https://github.com/FangfeiLi05/FitMut2 ), in the hope that it can find broader use by the microbial evolution community.

Paths to adaptation under fluctuating nitrogen starvation: The spectrum of adaptive mutations in Saccharomyces cerevisiae is shaped by retrotransposons and microhomology-mediated recombination

There are many mechanisms that give rise to genomic change: while point mutations are often emphasized in genomic analyses, evolution acts upon many other types of genetic changes that can result in less subtle perturbations. Changes in chromosome structure, DNA copy number, and novel transposon insertions all create large genomic changes, which can have correspondingly large impacts on phenotypes and fitness. In this study we investigate the spectrum of adaptive mutations that arise in a population under consistently fluctuating nitrogen conditions. We specifically contrast these adaptive alleles and the mutational mechanisms that create them, with mechanisms of adaptation under batch glucose limitation and constant selection in low, non-fluctuating nitrogen conditions to address if and how selection dynamics influence the molecular mechanisms of evolutionary adaptation. We observe that retrotransposon activity accounts for a substantial number of adaptive events, along with microhomology-mediated mechanisms of insertion, deletion, and gene conversion. In addition to loss of function alleles, which are often exploited in genetic screens, we identify putative gain of function alleles and alleles acting through as-of-yet unclear mechanisms. Taken together, our findings emphasize that how selection (fluctuating vs. non-fluctuating) is applied also shapes adaptation, just as the selective pressure (nitrogen vs. glucose) does itself. Fluctuating environments can activate different mutational mechanisms, shaping adaptive events accordingly. Experimental evolution, which allows a wider array of adaptive events to be assessed, is thus a complementary approach to both classical genetic screens and natural variation studies to characterize the genotype-to-phenotype-to-fitness map.

Widespread epistasis among beneficial genetic variants revealed by high-throughput genome editing

The phenotypic effect of any genetic variant can be altered by variation at other genomic loci. Known as epistasis, these genetic interactions shape the genotype-phenotype map of every species, yet their origins remain poorly understood. To investigate this, we employed high-throughput genome editing to measure the fitness effects of 1,826 naturally polymorphic variants in four strains of Saccharomyces cerevisiae. About 31% of variants affect fitness, of which 24% have strain-specific fitness effects indicative of epistasis. We found that beneficial variants are more likely to exhibit genetic interactions and that these interactions can be mediated by specific traits such as flocculation ability. This work suggests that adaptive evolution will often involve trade-offs where a variant is only beneficial in some genetic backgrounds, potentially explaining why many beneficial variants remain polymorphic. In sum, we provide a framework to understand the factors influencing epistasis with single-nucleotide resolution, revealing widespread epistasis among beneficial variants.

Long-Term Adaptation to Galactose as a Sole Carbon Source Selects for Mutations Outside the Canonical GAL Pathway

Galactose is a secondary fermentable sugar that requires specific regulatory and structural genes for its assimilation, which are under catabolite repression by glucose. When glucose is absent, the catabolic repression is attenuated, and the structural GAL genes are fully activated. In Saccharomyces cerevisiae, the GAL pathway is under selection in environments where galactose is present. However, it is unclear the adaptive strategies in response to long-term propagation in galactose as a sole carbon source in laboratory evolution experiments. Here, we performed a 4,000-generation evolution experiment using 48 diploid Saccharomyces cerevisiae populations to study adaptation in galactose. We show that fitness gains were greater in the galactose-evolved population than in identically evolved populations with glucose as a sole carbon source. Whole-genome sequencing of 96 evolved clones revealed recurrent de novo single nucleotide mutations in candidate targets of selection, copy number variations, and ploidy changes. We find that most mutations that improve fitness in galactose lie outside of the canonical GAL pathway. Reconstruction of specific evolved alleles in candidate target of selection, SEC23 and IRA1, showed a significant increase in fitness in galactose compared to glucose. In addition, most of our evolved populations (28/46; 61%) fixed aneuploidies on Chromosome VIII, suggesting a parallel adaptive amplification. Finally, we show greater loss of extrachromosomal elements in our glucose-evolved lineages compared with previous glucose evolution. Broadly, these data further our understanding of the evolutionary pressures that drive adaptation to less-preferred carbon sources.

Directed evolution as an approach to increase fructose utilization in synthetic grape juice by wine yeast AWRI 796

In winemaking, slow or stuck alcoholic fermentation can impact processing efficiency and wine quality. Residual fructose in the later stages of fermentation can leave the wine 'out of specification' unless removed, which requires reinoculation or use of a more fructophilic yeast. As such, robust, fermentation efficient strains are still highly desirable to reduce this risk. We report on a combined EMS mutagenesis and Directed Evolution (DE) approach as a 'proof of concept' to improve fructose utilization and decrease fermentation duration. One evolved isolate, Tee 9, was evaluated against the parent, AWRI 796 in defined medium (CDGJM) and Semillon juice. Interestingly, Tee 9 exhibited improved fermentation in CDGJM at several nitrogen contents, but not in juice. Genomic comparison between AWRI 796 and Tee 9 identified 371 mutations, but no chromosomal copy number variation. A total of 95 noncoding and 276 coding mutations were identified in 297 genes (180 of which encode proteins with one or more substitutions). Whilst introduction of two of these, Gid7 (E726K) or Fba1 (G135S), into AWRI 796 did not lead to the fermentation improvement seen in Tee 9, similar allelic swaps with the other mutations are needed to understand Tee 9's adaption to CDGJM. Furthermore, the 378 isolates, potentially mutagenized but with the same genetic background, are likely a useful resource for future phenotyping and genome-wide association studies.

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.

Population genetics of polymorphism and divergence in rapidly evolving populations

In rapidly evolving populations, numerous beneficial and deleterious mutations can arise and segregate within a population at the same time. In this regime, evolutionary dynamics cannot be analyzed using traditional population genetic approaches that assume that sites evolve independently. Instead, the dynamics of many loci must be analyzed simultaneously. Recent work has made progress by first analyzing the fitness variation within a population, and then studying how individual lineages interact with this traveling fitness wave. However, these "traveling wave" models have previously been restricted to extreme cases where selection on individual mutations is either much faster or much slower than the typical coalescent timescale Tc. In this work, we show how the traveling wave framework can be extended to intermediate regimes in which the scaled fitness effects of mutations (Tcs) are neither large nor small compared to one. This enables us to describe the dynamics of populations subject to a wide range of fitness effects, and in particular, in cases where it is not immediately clear which mutations are most important in shaping the dynamics and statistics of genetic diversity. We use this approach to derive new expressions for the fixation probabilities and site frequency spectra of mutations as a function of their scaled fitness effects, along with related results for the coalescent timescale Tc and the rate of adaptation or Muller's ratchet. We find that competition between linked mutations can have a dramatic impact on the proportions of neutral and selected polymorphisms, which is not simply summarized by the scaled selection coefficient Tcs. We conclude by discussing the implications of these results for population genetic inferences.

Intra- and inter-host evolution of H9N2 influenza A virus in Japanese quail

Influenza A viruses (IAVs) are constantly evolving. Crucial steps in the infection cycle, such as sialic acid (SA) receptor binding on the host cell surface, can either promote or hamper the emergence of new variants. We previously assessed the relative fitness in Japanese quail of H9N2 variant viruses differing at a single amino acid position, residue 216 in the hemagglutinin (HA) viral surface protein. This site is known to modulate SA recognition. Our prior study generated a valuable set of longitudinal samples from quail transmission groups where the inoculum comprised different mixed populations of HA 216 variant viruses. Here, we leveraged these samples to examine the evolutionary dynamics of viral populations within and between inoculated and naïve contact quails. We found that positive selection dominated HA gene evolution, but fixation of the fittest variant depended on the competition mixture. Analysis of the whole genome revealed further evidence of positive selection acting both within and between hosts. Positive selection drove fixation of variants in non-HA segments within inoculated and contact quails. Importantly, transmission bottlenecks were modulated by the molecular signature at HA 216, revealing viral receptor usage as a determinant of transmitted diversity. Overall, we show that selection strongly shaped the evolutionary dynamics within and between quails. These findings support the notion that selective processes act effectively on IAV populations in poultry hosts, facilitating rapid viral evolution in this ecological niche.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2^nd step of adaptation compared to the 1^st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2^nd adaptive step are contingent on the 1^st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Analyses of both natural and experimental evolution suggest that adaptation depends on the evolutionary past and adaptive potential decreases over time. Here, by tracking yeast adaptation with DNA barcoding, the authors show that such evolutionary phenomena can be observed even after a single adaptive step.

Adaptive response to wine selective pressures shapes the genome of a Saccharomyces interspecies hybrid

During industrial processes, yeasts are exposed to harsh conditions, which eventually lead to adaptation of the strains. In the laboratory, it is possible to use experimental evolution to link the evolutionary biology response to these adaptation pressures for the industrial improvement of a specific yeast strain. In this work, we aimed to study the adaptation of a wine industrial yeast in stress conditions of the high ethanol concentrations present in stopped fermentations and secondary fermentations in the processes of champagne production. We used a commercial Saccharomyces cerevisiae × S. uvarum hybrid and assessed its adaptation in a modified synthetic must (M-SM) containing high ethanol, which also contained metabisulfite, a preservative that is used during wine fermentation as it converts to sulfite. After the adaptation process under these selected stressful environmental conditions, the tolerance of the adapted strain (H14A7-etoh) to sulfite and ethanol was investigated, revealing that the adapted hybrid is more resistant to sulfite compared to the original H14A7 strain, whereas ethanol tolerance improvement was slight. However, a trade-off in the adapted hybrid was found, as it had a lower capacity to ferment glucose and fructose in comparison with H14A7. Hybrid genomes are almost always unstable, and different signals of adaptation on H14A7-etoh genome were detected. Each subgenome present in the adapted strain had adapted differently. Chromosome aneuploidies were present in S. cerevisiae chromosome III and in S. uvarum chromosome VII–XVI, which had been duplicated. Moreover, S. uvarum chromosome I was not present in H14A7-etoh and a loss of heterozygosity (LOH) event arose on S. cerevisiae chromosome I. RNA-sequencing analysis showed differential gene expression between H14A7-etoh and H14A7, which can be easily correlated with the signals of adaptation that were found on the H14A7-etoh genome. Finally, we report alterations in the lipid composition of the membrane, consistent with conserved tolerance mechanisms.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

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.

Quorum sensing provides a molecular mechanism for evolution to tune and maintain investment in cooperation

As selection frequently favors noncooperating defectors in mixed populations with cooperators, mechanisms that promote cooperation stability clearly exist. One potential mechanism is bacterial cell-to-cell communication, quorum sensing (QS), which can allow cooperators to prevent invasion by defectors. However, the impact of QS on widespread maintenance of cooperation in well-mixed conditions has not been experimentally demonstrated over extended evolutionary timescales. Here, we use wild-type (WT) Vibrio campbellii that regulates cooperation with QS and an unconditional cooperating (UC) mutant to examine the evolutionary origins and dynamics of novel defectors during a long-term evolution experiment. We found that UC lineages were completely outcompeted by defectors, whereas functioning QS enabled the maintenance of cooperative variants in most WT populations. Sequencing evolved populations revealed multiple luxR mutations that swept the UC lineages. However, the evolution of mutant lineages with reduced levels of bioluminescence (dims) occurred in many WT lineages. These dim variants also decreased other cooperative phenotypes regulated by QS, including protease production, indicating they result from changes to QS regulation. This diminished investment phenotype optimizes a tradeoff between cooperative input and growth output and suggests that decreasing the cost of QS could be a favorable strategy for maintaining the cooperative behaviors it regulates.

Evolutionary graph theory derived from eco-evolutionary dynamics

A biologically motivated individual-based framework for evolution in network-structured populations is developed that can accommodate eco-evolutionary dynamics. This framework is used to construct a network birth and death model. The evolutionary graph theory model, which considers evolutionary dynamics only, is derived as a special case, highlighting additional assumptions that diverge from real biological processes. This is achieved by introducing a negative ecological feedback loop that suppresses ecological dynamics by forcing births and deaths to be coupled. We also investigate how fitness, a measure of reproductive success used in evolutionary graph theory, is related to the life-history of individuals in terms of their birth and death rates. In simple networks, these ecologically motivated dynamics are used to provide new insight into the spread of adaptive mutations, both with and without clonal interference. For example, the star network, which is known to be an amplifier of selection in evolutionary graph theory, can inhibit the spread of adaptive mutations when individuals can die naturally.

Evolutionary dynamics and structural consequences of de novo beneficial mutations and mutant lineages arising in a constant environment

BACKGROUND

Microbial evolution experiments can be used to study the tempo and dynamics of evolutionary change in asexual populations, founded from single clones and growing into large populations with multiple clonal lineages. High-throughput sequencing can be used to catalog de novo mutations as potential targets of selection, determine in which lineages they arise, and track the fates of those lineages. Here, we describe a long-term experimental evolution study to identify targets of selection and to determine when, where, and how often those targets are hit.

RESULTS

We experimentally evolved replicate Escherichia coli populations that originated from a mutator/nonsense suppressor ancestor under glucose limitation for between 300 and 500 generations. Whole-genome, whole-population sequencing enabled us to catalog 3346 de novo mutations that reached > 1% frequency. We sequenced the genomes of 96 clones from each population when allelic diversity was greatest in order to establish whether mutations were in the same or different lineages and to depict lineage dynamics. Operon-specific mutations that enhance glucose uptake were the first to rise to high frequency, followed by global regulatory mutations. Mutations related to energy conservation, membrane biogenesis, and mitigating the impact of nonsense mutations, both ancestral and derived, arose later. New alleles were confined to relatively few loci, with many instances of identical mutations arising independently in multiple lineages, among and within replicate populations. However, most never exceeded 10% in frequency and were at a lower frequency at the end of the experiment than at their maxima, indicating clonal interference. Many alleles mapped to key structures within the proteins that they mutated, providing insight into their functional consequences.

CONCLUSIONS

Overall, we find that when mutational input is increased by an ancestral defect in DNA repair, the spectrum of high-frequency beneficial mutations in a simple, constant resource-limited environment is narrow, resulting in extreme parallelism where many adaptive mutations arise but few ever go to fixation.

Adaptation is influenced by the complexity of environmental change during evolution in a dynamic environment

The environmental conditions of microorganisms' habitats may fluctuate in unpredictable ways, such as changes in temperature, carbon source, pH, and salinity to name a few. Environmental heterogeneity presents a challenge to microorganisms, as they have to adapt not only to be fit under a specific condition, but they must also be robust across many conditions and be able to deal with the switch between conditions itself. While experimental evolution has been used to gain insight into the adaptive process, this has largely been in either unvarying or consistently varying conditions. In cases where changing environments have been investigated, relatively little is known about how such environments influence the dynamics of the adaptive process itself, as well as the genetic and phenotypic outcomes. We designed a systematic series of evolution experiments where we used two growth conditions that have differing timescales of adaptation and varied the rate of switching between them. We used lineage tracking to follow adaptation, and whole genome sequenced adaptive clones from each of the experiments. We find that both the switch rate and the order of the conditions influences adaptation. We also find different adaptive outcomes, at both the genetic and phenotypic levels, even when populations spent the same amount of total time in the two different conditions, but the order and/or switch rate differed. Thus, in a variable environment adaptation depends not only on the nature of the conditions and phenotypes under selection, but also on the complexity of the manner in which those conditions are combined to result in a given dynamic environment.

Whole-genome sequencing from the New Zealand Saccharomyces cerevisiae population reveals the genomic impacts of novel microbial range expansion

Saccharomyces cerevisiae is extensively utilized for commercial fermentation, and is also an important biological model; however, its ecology has only recently begun to be understood. Through the use of whole-genome sequencing, the species has been characterized into a number of distinct subpopulations, defined by geographical ranges and industrial uses. Here, the whole-genome sequences of 104 New Zealand (NZ) S. cerevisiae strains, including 52 novel genomes, are analyzed alongside 450 published sequences derived from various global locations. The impact of S. cerevisiae novel range expansion into NZ was investigated and these analyses reveal the positioning of NZ strains as a subgroup to the predominantly European/wine clade. A number of genomic differences with the European group correlate with range expansion into NZ, including 18 highly enriched single-nucleotide polymorphism (SNPs) and novel Ty1/2 insertions. While it is not possible to categorically determine if any genetic differences are due to stochastic process or the operations of natural selection, we suggest that the observation of NZ-specific copy number increases of four sugar transporter genes in the HXT family may reasonably represent an adaptation in the NZ S. cerevisiae subpopulation, and this correlates with the observations of copy number changes during adaptation in small-scale experimental evolution studies.

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.

Replicative aging is associated with loss of genetic heterogeneity from extrachromosomal circular DNA in Saccharomyces cerevisiae

Circular DNA can arise from all parts of eukaryotic chromosomes. In yeast, circular ribosomal DNA (rDNA) accumulates dramatically as cells age, however little is known about the accumulation of other chromosome-derived circles or the contribution of such circles to genetic variation in aged cells. We profiled circular DNA in Saccharomyces cerevisiae populations sampled when young and after extensive aging. Young cells possessed highly diverse circular DNA populations but 94% of the circular DNA were lost after ∼15 divisions, whereas rDNA circles underwent massive accumulation to >95% of circular DNA. Circles present in both young and old cells were characterized by replication origins including circles from unique regions of the genome and repetitive regions: rDNA and telomeric Y' regions. We further observed that circles can have flexible inheritance patterns: [HXT6/7circle] normally segregates to mother cells but in low glucose is present in up to 50% of cells, the majority of which must have inherited this circle from their mother. Interestingly, [HXT6/7circle] cells are eventually replaced by cells carrying stable chromosomal HXT6 HXT6/7 HXT7 amplifications, suggesting circular DNAs are intermediates in chromosomal amplifications. In conclusion, the heterogeneity of circular DNA offers flexibility in adaptation, but this heterogeneity is remarkably diminished with age.

Application of simultaneous selective pressures slows adaptation

Beneficial mutations that arise in an evolving asexual population may compete or interact in ways that alter the overall rate of adaptation through mechanisms such as clonal or functional interference. The application of multiple selective pressures simultaneously may allow for a greater number of adaptive mutations, increasing the opportunities for competition between selectively advantageous alterations, and thereby reducing the rate of adaptation. We evolved a strain of Saccharomyces cerevisiae that could not produce its own histidine or uracil for ~500 generations under one or three selective pressures: limitation of the concentration of glucose, histidine, and/or uracil in the media. The rate of adaptation was obtained by measuring evolved relative fitness using competition assays. Populations evolved under a single selective pressure showed a statistically significant increase in fitness on those pressures relative to the ancestral strain, but the populations evolved on all three pressures did not show a statistically significant increase in fitness over the ancestral strain on any single pressure. Simultaneously limiting three essential nutrients for a population of S. cerevisiae effectively slows the rate of evolution on any one of the three selective pressures applied, relative to the single selective pressure cases. We identify possible mechanisms for fitness changes seen between populations evolved on one or three limiting nutrient pressures by high-throughput sequencing. Adding multiple selective pressures to evolving disease like cancer and infectious diseases could reduce the rate of adaptation and thereby may slow disease progression, prolong drug efficacy and prevent deaths.

Current Status and Applications of Adaptive Laboratory Evolution in Industrial Microorganisms

Adaptive laboratory evolution (ALE) is an evolutionary engineering approach in artificial conditions that improves organisms through the imitation of natural evolution. Due to the development of multi-level omics technologies in recent decades, ALE can be performed for various purposes at the laboratory level. This review delineates the basics of the experimental design of ALE based on several ALE studies of industrial microbial strains and updates current strategies combined with progressed metabolic engineering, in silico modeling and automation to maximize the evolution efficiency. Moreover, the review sheds light on the applicability of ALE as a strain development approach that complies with non-recombinant preferences in various food industries. Overall, recent progress in the utilization of ALE for strain development leading to successful industrialization is discussed.

Exploiting strain diversity and rational engineering strategies to enhance recombinant cellulase secretion by Saccharomyces cerevisiae

Consolidated bioprocessing (CBP) of lignocellulosic material into bioethanol has progressed in the past decades; however, several challenges still exist which impede the industrial application of this technology. Identifying the challenges that exist in all unit operations is crucial and needs to be optimised, but only the barriers related to the secretion of recombinant cellulolytic enzymes in Saccharomyces cerevisiae will be addressed in this review. Fundamental principles surrounding CBP as a biomass conversion platform have been established through the successful expression of core cellulolytic enzymes, namely β-glucosidases, endoglucanases, and exoglucanases (cellobiohydrolases) in S. cerevisiae. This review will briefly address the challenges involved in the construction of an efficient cellulolytic yeast, with particular focus on the secretion efficiency of cellulases from this host. Additionally, strategies for studying enhanced cellulolytic enzyme secretion, which include both rational and reverse engineering approaches, will be discussed. One such technique includes bio-engineering within genetically diverse strains, combining the strengths of both natural strain diversity and rational strain development. Furthermore, with the advancement in next-generation sequencing, studies that utilise this method of exploiting intra-strain diversity for industrially relevant traits will be reviewed. Finally, future prospects are discussed for the creation of ideal CBP strains with high enzyme production levels.Key Points• Several challenges are involved in the construction of efficient cellulolytic yeast, in particular, the secretion efficiency of cellulases from the hosts.• Strategies for enhancing cellulolytic enzyme secretion, a core requirement for CBP host microorganism development, include both rational and reverse engineering approaches.• One such technique includes bio-engineering within genetically diverse strains, combining the strengths of both natural strain diversity and rational strain development.

A spontaneous complex structural variant in rcan-1 increases exploratory behavior and laboratory fitness of Caenorhabditis elegans

Over long evolutionary timescales, major changes to the copy number, function, and genomic organization of genes occur, however, our understanding of the individual mutational events responsible for these changes is lacking. In this report, we study the genetic basis of adaptation of two strains of C. elegans to laboratory food sources using competition experiments on a panel of 89 recombinant inbred lines (RIL). Unexpectedly, we identified a single RIL with higher relative fitness than either of the parental strains. This strain also displayed a novel behavioral phenotype, resulting in higher propensity to explore bacterial lawns. Using bulk-segregant analysis and short-read resequencing of this RIL, we mapped the change in exploration behavior to a spontaneous, complex rearrangement of the rcan-1 gene that occurred during construction of the RIL panel. We resolved this rearrangement into five unique tandem inversion/duplications using Oxford Nanopore long-read sequencing. rcan-1 encodes an ortholog to human RCAN1/DSCR1 calcipressin gene, which has been implicated as a causal gene for Down syndrome. The genomic rearrangement in rcan-1 creates two complete and two truncated versions of the rcan-1 coding region, with a variety of modified 5’ and 3’ non-coding regions. While most copy-number variations (CNVs) are thought to act by increasing expression of duplicated genes, these changes to rcan-1 ultimately result in the reduction of its whole-body expression due to changes in the upstream regions. By backcrossing this rearrangement into a common genetic background to create a near isogenic line (NIL), we demonstrate that both the competitive advantage and exploration behavioral changes are linked to this complex genetic variant. This NIL strain does not phenocopy a strain containing an rcan-1 loss-of-function allele, which suggests that the residual expression of rcan-1 is necessary for its fitness effects. Our results demonstrate how colonization of new environments, such as those encountered in the laboratory, can create evolutionary pressure to modify gene function. This evolutionary mismatch can be resolved by an unexpectedly complex genetic change that simultaneously duplicates and diversifies a gene into two uniquely regulated genes. Our work shows how complex rearrangements can act to modify gene expression in ways besides increased gene dosage.

Evolution acts on genetic variants that modify phenotypes that increase the likelihood of staying alive and passing on these genetic changes to subsequent generations (i.e. fitness). There is general interest in understanding the types of genetic variants that can increase fitness in specific environments. One route that fitness can be increased is through changes in behavior, such as finding new food sources. Here, we identify a spontaneous genetic change that increases exploration behavior and fitness of animals in laboratory environments. Interestingly, this genetic change is not a simple genetic change that deletes or changes the sequence of a protein product, but rather a complex structural variant that simultaneously duplicates the rcan-1 gene and also modifies its expression in a number of tissues. Our work demonstrates how a complex structural change can duplicate a gene, modify the DNA control regions that determine its cellular sites of action, and confer a fitness advantage that could lead to its spread in a population.

CgSTE11 mediates cross tolerance to multiple environmental stressors in Candida glabrata

Candida glabrata is a human commensal and an opportunistic human fungal pathogen. It is more closely related to the model yeast Saccharomyces cerevisiae than other Candida spp. Compared with S. cerevisiae, C. glabrata exhibits higher innate tolerance to various environmental stressors, including hyperthermal stress. Here we investigate the molecular mechanisms of C. glabrata adaptation to heat stress via adaptive laboratory evolution. We show that all parallel evolved populations readily adapt to hyperthermal challenge (from 47 °C to 50 °C) and exhibit convergence in evolved phenotypes with extensive cross-tolerance to various other environmental stressors such as oxidants, acids, and alcohols. Genome resequencing identified fixation of mutations in CgSTE11 in all parallel evolved populations. The CgSTE11 homolog in S. cerevisiae plays crucial roles in various mitogen-activated protein kinase (MAPK) signaling pathways, but its role is less understood in C. glabrata. Subsequent verification confirmed that CgSTE11 is important in hyperthermal tolerance and the observed extensive cross-tolerance to other environmental stressors. These results support the hypothesis that CgSTE11 mediates cross-talks between MAPK signaling pathways in C. glabrata in response to environmental challenges.

The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology

Harnessing the process of natural selection to obtain and understand new microbial phenotypes has become increasingly possible due to advances in culturing techniques, DNA sequencing, bioinformatics, and genetic engineering. Accordingly, Adaptive Laboratory Evolution (ALE) experiments represent a powerful approach both to investigate the evolutionary forces influencing strain phenotypes, performance, and stability, and to acquire production strains that contain beneficial mutations. In this review, we summarize and categorize the applications of ALE to various aspects of microbial physiology pertinent to industrial bioproduction by collecting case studies that highlight the multitude of ways in which evolution can facilitate the strain construction process. Further, we discuss principles that inform experimental design, complementary approaches such as computational modeling that help maximize utility, and the future of ALE as an efficient strain design and build tool driven by growing adoption and improvements in automation.

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.

Clonal interference can cause wavelet-like oscillations of multilocus linkage disequilibrium

Within-host adaptation of pathogens such as human immunodeficiency virus (HIV) often occurs at more than two loci. Multiple beneficial mutations may arise simultaneously on different genetic backgrounds and interfere, affecting each other's fixation trajectories. Here, we explore how these evolutionary dynamics are mirrored in multilocus linkage disequilibrium (MLD), a measure of multi-way associations between alleles. In the parameter regime corresponding to HIV, we show that deterministic early infection models induce MLD to oscillate over time in a wavelet-like fashion. We find that the frequency of these oscillations is proportional to the rate of adaptation. This signature is robust to drift, but can be eroded by high variation in fitness effects of beneficial mutations. Our findings suggest that MLD oscillations could be used as a signature of interference among multiple equally advantageous mutations and may aid the interpretation of MLD in data.

Designing Automated, High-throughput, Continuous Cell Growth Experiments Using eVOLVER

Continuous culture methods enable cells to be grown under quantitatively controlled environmental conditions, and are thus broadly useful for measuring fitness phenotypes and improving our understanding of how genotypes are shaped by selection. Extensive recent efforts to develop and apply niche continuous culture devices have revealed the benefits of conducting new forms of cell culture control. This includes defining custom selection pressures and increasing throughput for studies ranging from long-term experimental evolution to genome-wide library selections and synthetic gene circuit characterization. The eVOLVER platform was recently developed to meet this growing demand: a continuous culture platform with a high degree of scalability, flexibility, and automation. eVOLVER provides a single standardizing platform that can be (re)-configured and scaled with minimal effort to perform many different types of high-throughput or multi-dimensional growth selection experiments. Here, a protocol is presented to provide users of the eVOLVER framework a description for configuring the system to conduct a custom, large-scale continuous growth experiment. Specifically, the protocol guides users on how to program the system to multiplex two selection pressures - temperature and osmolarity - across many eVOLVER vials in order to quantify fitness landscapes of Saccharomyces cerevisiae mutants at fine resolution. We show how the device can be configured both programmatically, through its open-source web-based software, and physically, by arranging fluidic and hardware layouts. The process of physically setting up the device, programming the culture routine, monitoring and interacting with the experiment in real-time over the internet, sampling vials for subsequent offline analysis, and post experiment data analysis are detailed. This should serve as a starting point for researchers across diverse disciplines to apply eVOLVER in the design of their own complex and high-throughput cell growth experiments to study and manipulate biological systems.

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

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

Improved discovery of genetic interactions using CRISPRiSeq across multiple environments

Large-scale genetic interaction (GI) screens in yeast have been invaluable for our understanding of molecular systems biology and for characterizing novel gene function. Owing in part to the high costs and long experiment times required, a preponderance of GI data has been generated in a single environmental condition. However, an unknown fraction of GIs may be specific to other conditions. Here, we developed a pooled-growth CRISPRi-based sequencing assay for GIs, CRISPRiSeq, which increases throughput such that GIs can be easily assayed across multiple growth conditions. We assayed the fitness of approximately 17,000 strains encompassing approximately 7700 pairwise interactions in five conditions and found that the additional conditions increased the number of GIs detected nearly threefold over the number detected in rich media alone. In addition, we found that condition-specific GIs are prevalent and improved the power to functionally classify genes. Finally, we found new links during respiratory growth between members of the Ras nutrient-sensing pathway and both the COG complex and a gene of unknown function. Our results highlight the potential of conditional GI screens to improve our understanding of cellular genetic networks.

Microevolution of the pathogenic yeasts Candida albicans and Candida glabrata during antifungal therapy and host infection

Infections by the pathogenic yeasts Candida albicans and Candida glabrata are among the most common fungal diseases. The success of these species as human pathogens is contingent on their ability to resist antifungal therapy and thrive within the human host. C. glabrata is especially resilient to azole antifungal treatment, while C. albicans is best known for its wide array of virulence features. The core mechanisms that underlie antifungal resistance and virulence in these pathogens has been continuously addressed, but the investigation on how such mechanisms evolve according to each environment is scarcer. This review aims to explore current knowledge on micro-evolution experiments to several treatment and host-associated conditions in C. albicans and C. glabrata. The analysis of adaptation strategies that evolve over time will allow to better understand the mechanisms by which Candida species are able to achieve stable phenotypes in real-life scenarios, which are the ones that should constitute the most interesting drug targets.

Experimental evolution: its principles and applications in developing stress-tolerant yeasts

Stress tolerance and resistance in industrial yeast strains are important attributes for cost-effective bioprocessing. The source of stress-tolerant yeasts ranges from extremophilic environments to laboratory engineered strains. However, industrial stress-tolerant yeasts are very rare in nature as the natural environment forces them to evolve traits that optimize survival and reproduction and not the ability to withstand harsh habitat-irrelevant industrial conditions. Experimental evolution is a frequent method used to uncover the mechanisms of evolution and microbial adaption towards environmental stresses. It optimizes biological systems by means of adaptation to environmental stresses and thus has immense power of development of robust stress-tolerant yeasts. This mini-review briefly outlines the basics and implications of evolution experiments and their applications to industrial biotechnology. This work is meant to serve as an introduction to those new to the field of experimental evolution, and as a guide to biologists working in the field of yeast stress response. Future perspectives of experimental evolution for potential biotechnological applications have also been elucidated.

Effects of Seawater on Carotenoid Production and Lipid Content of Engineered Saccharomyces cerevisiae

The use of seawater in fermentation can potentially reduce the freshwater burden in the bio-based production of chemicals and fuels. We previously developed a Saccharomyces cerevisiae carotenoids hyperproducer SM14 capable of accumulating 18 mg g−1 DCW (DCW: dry cell weight) of β-carotene in rich media (YPD). In this work, the impacts of seawater on the carotenoid production of SM14 were investigated. When using nutrient-reduced media (0.1× YNB) in freshwater the β-carotene production of SM14 was 6.51 ± 0.37 mg g−1 DCW; however in synthetic seawater, the production was increased to 8.67 ± 0.62 mg g−1 DCW. We found that this improvement was partially due to the NaCl present in the synthetic seawater, since supplementation of 0.5 M NaCl in freshwater increased β-carotene production to 11.85 ± 0.77 mg g−1 DCW. The combination of synthetic seawater with higher carbon-to-nitrogen ratio (C:N = 50) further improved the β-carotene production to 10.44 ± 0.35 mg g−1 DCW. We further showed that the carotenoid production improvement in these conditions is related with lipid content and composition. These results demonstrated the benefit of using seawater to improve the production of carotenoids in S. cerevisiae, and have the potential to expand the utilization of seawater.