Development of a Comprehensive Genotype-to-Fitness Map of Adaptation-Driving Mutations in Yeast

The dynamics of loss of heterozygosity events in genomes

Genomic instability is a hallmark of tumorigenesis, yet it also plays an essential role in evolution. Large-scale population genomics studies have highlighted the importance of loss of heterozygosity (LOH) events, which have long been overlooked in the context of genetic diversity and instability. Among various types of genomic mutations, LOH events are the most common and affect a larger portion of the genome. They typically arise from recombination-mediated repair of double-strand breaks (DSBs) or from lesions that are processed into DSBs. LOH events are critical drivers of genetic diversity, enabling rapid phenotypic variation and contributing to tumorigenesis. Understanding the accumulation of LOH, along with its underlying mechanisms, distribution, and phenotypic consequences, is therefore crucial. In this review, we explore the spectrum of LOH events, their mechanisms, and their impact on fitness and phenotype, drawing insights from Saccharomyces cerevisiae to cancer. We also emphasize the role of LOH in genomic instability, disease, and genome evolution.

A high-resolution two-step evolution experiment in yeast reveals a shift from pleiotropic to modular adaptation

Evolution by natural selection is expected to be a slow and gradual process. In particular, the mutations that drive evolution are predicted to be small and modular, incrementally improving a small number of traits. However, adaptive mutations identified early in microbial evolution experiments, cancer, and other systems often provide substantial fitness gains and pleiotropically improve multiple traits at once. We asked whether such pleiotropically adaptive mutations are common throughout adaptation or are instead a rare feature of early steps in evolution that tend to target key signaling pathways. To do so, we conducted barcoded second-step evolution experiments initiated from 5 first-step mutations identified from a prior yeast evolution experiment. We then isolated hundreds of second-step mutations from these evolution experiments, measured their fitness and performance in several growth phases, and conducted whole genome sequencing of the second-step clones. Here, we found that while the vast majority of mutants isolated from the first-step of evolution in this condition show patterns of pleiotropic adaptation-improving both performance in fermentation and respiration growth phases-second-step mutations show a shift towards modular adaptation, mostly improving respiration performance and only rarely improving fermentation performance. We also identified a shift in the molecular basis of adaptation from genes in cellular signaling pathways towards genes involved in respiration and mitochondrial function. Our results suggest that the genes in cellular signaling pathways may be more likely to provide large, adaptively pleiotropic benefits to the organism due to their ability to coherently affect many phenotypes at once. As such, these genes may serve as the source of pleiotropic adaptation in the early stages of evolution, and once these become exhausted, organisms then adapt more gradually, acquiring smaller, more modular mutations.

Asynchronous abundance fluctuations can drive giant genotype frequency fluctuations

Large stochastic population abundance fluctuations are ubiquitous across the tree of life, impacting the predictability and outcomes of population dynamics. It is generally thought that abundance fluctuations with a Taylor's law exponent of two do not strongly impact evolution. However, we argue that such abundance fluctuations can lead to substantial genotype frequency fluctuations if different genotypes in a population experience these fluctuations asynchronously. By serially diluting mixtures of two closely related Escherichia coli strains, we show that such asynchrony can occur, leading to giant frequency fluctuations that far exceed expectations from genetic drift. We develop an effective model explaining that the abundance fluctuations arise from correlated offspring numbers between individuals, and the large frequency fluctuations result from (even slight) decoupling in offspring number correlations between genotypes. The model quantitatively predicts the observed abundance and frequency fluctuation scaling. Initially close trajectories diverge exponentially, suggesting that chaotic dynamics may underpin the excess frequency fluctuations. Our findings suggest that decoupling noise is also present in mixed-genotype Saccharomyces cerevisiae populations. Theoretical analyses demonstrate that decoupling noise can strongly influence evolutionary outcomes, in a manner distinct from genetic drift. Given the generic nature of these frequency fluctuations, we expect them to be widespread across biological populations.

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.

Gain- and loss-of-function alleles within signaling pathways lead to phenotypic diversity among individuals

Understanding how phenotypic diversity is generated is an important question in biology. We explored phenotypic diversity among wild yeast isolates (Saccharomyces cerevisiae) and found variation in the activity of MAPK signaling pathways as a contributing mechanism. To uncover the genetic basis of this mechanism, we identified 1957 SNPs in 62 candidate genes encoding signaling proteins from a MAPK signaling module within a large collection of yeast (>1500 individuals). Follow-up testing identified functionally relevant variants in key signaling proteins. Loss-of-function (LOF) alleles in a PAK kinase impacted protein stability and pathway specificity decreasing filamentous growth and mating phenotypes. In contrast, gain-of-function (GOF) alleles in G-proteins that were hyperactivating induced filamentous growth. Similar amino acid substitutions in G-proteins were identified in metazoans that in some cases were fixed in multicellular lineages including humans, suggesting hyperactivating GOF alleles may play roles in generating phenotypic diversity across eukaryotes. A mucin signaler that regulates MAPK activity was also found to contain a prevalance of presumed GOF alleles amoung individuals based on changes in mucin repeat numbers. Thus, genetic variation in signaling pathways may act as a reservoir for generating phenotypic diversity across eukaryotes.

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

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

Environment by environment interactions (ExE) differ across genetic backgrounds (ExExG)

While the terms "gene-by-gene interaction" (GxG) and "gene-by-environment interaction" (GxE) are widely recognized in the fields of quantitative and evolutionary genetics, "environment-byenvironment interaction" (ExE) is a term used less often. In this study, we find that environmentby-environment interactions are a meaningful driver of phenotypes, and moreover, that they differ across different genotypes (suggestive of ExExG). To support this conclusion, we analyzed a large dataset of roughly 1,000 mutant yeast strains with varying degrees of resistance to different antifungal drugs. Our findings reveal that the effectiveness of a drug combination, relative to single drugs, often differs across drug resistant mutants. Remarkably, even mutants that differ by only a single nucleotide change can have dramatically different drug × drug (ExE) interactions. We also introduce a new framework that more accurately predicts the direction and magnitude of ExE interactions for some mutants. Understanding how ExE interactions change across genotypes (ExExG) is crucial not only for modeling the evolution of pathogenic microbes, but also for enhancing our knowledge of the underlying cell biology and the sources of phenotypic variance within populations. While the significance of ExExG interactions has been overlooked in evolutionary and population genetics, these fields and others stand to benefit from understanding how these interactions shape the complex behavior of living systems.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Genetic conflicts in budding yeast: The 2μ plasmid as a model selfish element

Antagonistic coevolution, arising from genetic conflict, can drive rapid evolution and biological innovation. Conflict can arise both between organisms and within genomes. This review focuses on budding yeasts as a model system for exploring intra- and inter-genomic genetic conflict, highlighting in particular the 2-micron (2μ) plasmid as a model selfish element. The 2μ is found widely in laboratory strains and industrial isolates of Saccharomyces cerevisiae and has long been known to cause host fitness defects. Nevertheless, the plasmid is frequently ignored in the context of genetic, fitness, and evolution studies. Here, I make a case for further exploring the evolutionary impact of the 2μ plasmid as well as other selfish elements of budding yeasts, discuss recent advances, and, finally, future directions for the field.

Environmental memory alters the fitness effects of adaptive mutations in fluctuating environments

Evolution in a static laboratory environment often proceeds via large-effect beneficial mutations that may become maladaptive in other environments. Conversely, natural settings require populations to endure environmental fluctuations. A sensible assumption is that the fitness of a lineage in a fluctuating environment is the time average of its fitness over the sequence of static conditions it encounters. However, transitions between conditions may pose entirely new challenges, which could cause deviations from this time average. To test this, we tracked hundreds of thousands of barcoded yeast lineages evolving in static and fluctuating conditions and subsequently isolated 900 mutants for pooled fitness assays in 15 environments. Here we find that fitness in fluctuating environments indeed often deviates from the time average, leading to fitness non-additivity. Moreover, closer examination reveals that fitness in one component of a fluctuating environment is often strongly influenced by the previous component. We show that this environmental memory is especially common for mutants with high variance in fitness across tested environments. We use a simple mathematical model and whole-genome sequencing to propose mechanisms underlying this effect, including lag time evolution and sensing mutations. Our results show that environmental fluctuations impact fitness and suggest that variance in static environments can explain these impacts.

Asynchronous abundance fluctuations can drive giant genotype frequency fluctuations

Large stochastic population abundance fluctuations are ubiquitous across the tree of life1-7, impacting the predictability of population dynamics and influencing eco-evolutionary outcomes. It has generally been thought that these large abundance fluctuations do not strongly impact evolution, as the relative frequencies of alleles in the population will be unaffected if the abundance of all alleles fluctuate in unison. However, we argue that large abundance fluctuations can lead to significant genotype frequency fluctuations if different genotypes within a population experience these fluctuations asynchronously. By serially diluting mixtures of two closely related E. coli strains, we show that such asynchrony can occur, leading to giant frequency fluctuations that far exceed expectations from models of genetic drift. We develop a flexible, effective model that explains the abundance fluctuations as arising from correlated offspring numbers between individuals, and the large frequency fluctuations result from even slight decoupling in offspring numbers between genotypes. This model accurately describes the observed abundance and frequency fluctuation scaling behaviors. Our findings suggest chaotic dynamics underpin these giant fluctuations, causing initially similar trajectories to diverge exponentially; subtle environmental changes can be magnified, leading to batch correlations in identical growth conditions. Furthermore, we present evidence that such decoupling noise is also present in mixed-genotype S. cerevisiae populations. We demonstrate that such decoupling noise can strongly influence evolutionary outcomes, in a manner distinct from genetic drift. Given the generic nature of asynchronous fluctuations, we anticipate that they are widespread in biological populations, significantly affecting evolutionary and ecological dynamics.

Improving the Accuracy of Bulk Fitness Assays by Correcting Barcode Processing Biases

Measuring the fitnesses of genetic variants is a fundamental objective in evolutionary biology. A standard approach for measuring microbial fitnesses in bulk involves labeling a library of genetic variants with unique sequence barcodes, competing the labeled strains in batch culture, and using deep sequencing to track changes in the barcode abundances over time. However, idiosyncratic properties of barcodes can induce nonuniform amplification or uneven sequencing coverage that causes some barcodes to be over- or under-represented in samples. This systematic bias can result in erroneous read count trajectories and misestimates of fitness. Here, we develop a computational method, named REBAR (Removing the Effects of Bias through Analysis of Residuals), for inferring the effects of barcode processing bias by leveraging the structure of systematic deviations in the data. We illustrate this approach by applying it to two independent data sets, and demonstrate that this method estimates and corrects for bias more accurately than standard proxies, such as GC-based corrections. REBAR mitigates bias and improves fitness estimates in high-throughput assays without introducing additional complexity to the experimental protocols, with potential applications in a range of experimental evolution and mutation screening contexts.

Design principles of multi-map variation in biological systems

Complexity in biology is often described using a multi-map hierarchical architecture, where the genotype, representing the encoded information, is mapped to the functional level, known as the phenotype, which is then connected to a latent phenotype we refer to as fitness. This underlying architecture governs the processes driving evolution. Furthermore, natural selection, along with other neutral forces, can, in turn, modify these maps. At each level, variation is observed. Here, I propose the need to establish principles that can aid in understanding the transformation of variation within this multi-map architecture. Specifically, I will introduce three, related to the presence of modulators, constraints, and the modular channeling of variation. By comprehending these design principles in various biological systems, we can gain better insights into the mechanisms underlying these maps and how they ultimately contribute to evolutionary dynamics.

The loci of environmental adaptation in a model eukaryote

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into 6 classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Self-replicating artificial neural networks give rise to universal evolutionary dynamics

In evolutionary models, mutations are exogenously introduced by the modeler, rather than endogenously introduced by the replicator itself. We present a new deep-learning based computational model, the self-replicating artificial neural network (SeRANN). We train it to (i) copy its own genotype, like a biological organism, which introduces endogenous spontaneous mutations; and (ii) simultaneously perform a classification task that determines its fertility. Evolving 1,000 SeRANNs for 6,000 generations, we observed various evolutionary phenomena such as adaptation, clonal interference, epistasis, and evolution of both the mutation rate and the distribution of fitness effects of new mutations. Our results demonstrate that universal evolutionary phenomena can naturally emerge in a self-replicator model when both selection and mutation are implicit and endogenous. We therefore suggest that SeRANN can be applied to explore and test various evolutionary dynamics and hypotheses.

Quantifying the adaptive landscape of commensal gut bacteria using high-resolution lineage tracking

Gut microbiota can adapt to their host environment by rapidly acquiring new mutations. However, the dynamics of this process are difficult to characterize in dominant gut species in their complex in vivo environment. Here we show that the fine-scale dynamics of genome-wide transposon libraries can enable quantitative inferences of these in vivo evolutionary forces. By analyzing >400,000 lineages across four human Bacteroides strains in gnotobiotic mice, we observed positive selection on thousands of cryptic variants - most of which were unrelated to their original gene knockouts. The spectrum of fitness benefits varied between species, and displayed diverse tradeoffs over time and in different dietary conditions, enabling inferences of their underlying function. These results suggest that within-host adaptations arise from an intense competition between numerous contending variants, which can strongly influence their emergent evolutionary tradeoffs.

Horizontal gene transfer facilitates the molecular reverse-evolution of antibiotic sensitivity in experimental populations of H. pylori

Reversing the evolution of traits harmful to humans, such as antimicrobial resistance, is a key ambition of applied evolutionary biology. A major impediment to reverse evolution is the relatively low spontaneous mutation rates that revert evolved genotypes back to their ancestral state. However, the repeated re-introduction of ancestral alleles by horizontal gene transfer (HGT) could make reverse evolution likely. Here we evolve populations of an antibiotic-resistant strain of Helicobacter pylori in growth conditions without antibiotics while introducing an ancestral antibiotic-sensitive allele by HGT. We evaluate reverse evolution using DNA sequencing and find that HGT facilitates the molecular reverse evolution of the antibiotic resistance allele, and that selection for high rates of HGT drives the evolution of increased HGT rates in low-HGT treatment populations. Finally, we use a theoretical model and carry out simulations to infer how the fitness costs of antibiotic resistance, rates of HGT and effects of genetic drift interact to determine the probability and predictability of reverse evolution.

Disruption of Leishmania flagellum attachment zone architecture causes flagellum loss

Leishmania are flagellated eukaryotic parasites that cause leishmaniasis and are closely related to the other kinetoplastid parasites such as Trypanosoma brucei. In all these parasites there is a cell membrane invagination at the base of the flagellum called the flagellar pocket, which is tightly associated with and sculpted by cytoskeletal structures including the flagellum attachment zone (FAZ). The FAZ is a complex interconnected structure linking the flagellum to the cell body and has critical roles in cell morphogenesis, function and pathogenicity. However, this structure varies dramatically in size and organisation between these different parasites, suggesting changes in protein localisation and function. Here, we screened the localisation and function of the Leishmania orthologues of T. brucei FAZ proteins identified in the genome-wide protein tagging project TrypTag. We identified 27 FAZ proteins and our deletion analysis showed that deletion of two FAZ proteins in the flagellum, FAZ27 and FAZ34 resulted in a reduction in cell body size, and flagellum loss in some cells. Furthermore, after null mutant generation, we observed distinct and reproducible changes to cell shape, demonstrating the ability of the parasite to adapt to morphological perturbations resulting from gene deletion. This process of adaptation has important implications for the study of Leishmania mutants.

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.

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.

Strong environmental memory revealed by experimental evolution in static and fluctuating environments

Evolution in a static environment, such as a laboratory setting with constant and uniform conditions, often proceeds via large-effect beneficial mutations that may become maladaptive in other environments. Conversely, natural settings require populations to endure environmental fluctuations. A sensible assumption is that the fitness of a lineage in a fluctuating environment is the time-average of its fitness over the sequence of static conditions it encounters. However, transitions between conditions may pose entirely new challenges, which could cause deviations from this time-average. To test this, we tracked hundreds of thousands of barcoded yeast lineages evolving in static and fluctuating conditions and subsequently isolated 900 mutants for pooled fitness assays in 15 environments. We find that fitness in fluctuating environments indeed often deviates from the expectation based on static components, leading to fitness non-additivity. Moreover, closer examination reveals that fitness in one component of a fluctuating environment is often strongly influenced by the previous component. We show that this environmental memory is especially common for mutants with high variance in fitness across tested environments, even if the components of the focal fluctuating environment are excluded from this variance. We employ a simple mathematical model and whole-genome sequencing to propose mechanisms underlying this effect, including lag time evolution and sensing mutations. Our results demonstrate that environmental fluctuations have large impacts on fitness and suggest that variance in static environments can explain these impacts.

Intra-FCY1: a novel system to identify mutations that cause protein misfolding

Protein misfolding is a common intracellular occurrence. Most mutations to coding sequences increase the propensity of the encoded protein to misfold. These misfolded molecules can have devastating effects on cells. Despite the importance of protein misfolding in human disease and protein evolution, there are fundamental questions that remain unanswered, such as, which mutations cause the most misfolding? These questions are difficult to answer partially because we lack high-throughput methods to compare the destabilizing effects of different mutations. Commonly used systems to assess the stability of mutant proteins in vivo often rely upon essential proteins as sensors, but misfolded proteins can disrupt the function of the essential protein enough to kill the cell. This makes it difficult to identify and compare mutations that cause protein misfolding using these systems. Here, we present a novel in vivo system named Intra-FCY1 that we use to identify mutations that cause misfolding of a model protein [yellow fluorescent protein (YFP)] in Saccharomyces cerevisiae. The Intra-FCY1 system utilizes two complementary fragments of the yeast cytosine deaminase Fcy1, a toxic protein, into which YFP is inserted. When YFP folds, the Fcy1 fragments associate together to reconstitute their function, conferring toxicity in media containing 5-fluorocytosine and hindering growth. But mutations that make YFP misfold abrogate Fcy1 toxicity, thus strains possessing misfolded YFP variants rise to high frequency in growth competition experiments. This makes such strains easier to study. The Intra-FCY1 system cancels localization of the protein of interest, thus can be applied to study the relative stability of mutant versions of diverse cellular proteins. Here, we confirm this method can identify novel mutations that cause misfolding, highlighting the potential for Intra-FCY1 to illuminate the relationship between protein sequence and stability.

High-throughput library transgenesis in Caenorhabditis elegans via Transgenic Arrays Resulting in Diversity of Integrated Sequences (TARDIS)

High-throughput transgenesis using synthetic DNA libraries is a powerful method for systematically exploring genetic function. Diverse synthesized libraries have been used for protein engineering, identification of protein-protein interactions, characterization of promoter libraries, developmental and evolutionary lineage tracking, and various other exploratory assays. However, the need for library transgenesis has effectively restricted these approaches to single-cell models. Here, we present Transgenic Arrays Resulting in Diversity of Integrated Sequences (TARDIS), a simple yet powerful approach to large-scale transgenesis that overcomes typical limitations encountered in multicellular systems. TARDIS splits the transgenesis process into a two-step process: creation of individuals carrying experimentally introduced sequence libraries, followed by inducible extraction and integration of individual sequences/library components from the larger library cassette into engineered genomic sites. Thus, transformation of a single individual, followed by lineage expansion and functional transgenesis, gives rise to thousands of genetically unique transgenic individuals. We demonstrate the power of this system using engineered, split selectable TARDIS sites in Caenorhabditis elegans to generate (1) a large set of individually barcoded lineages and (2) transcriptional reporter lines from predefined promoter libraries. We find that this approach increases transformation yields up to approximately 1000-fold over current single-step methods. While we demonstrate the utility of TARDIS using C. elegans, in principle the process is adaptable to any system where experimentally generated genomic loci landing pads and diverse, heritable DNA elements can be generated.

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.

Resolving Deleterious and Near-Neutral Effects Requires Different Pooled Fitness Assay Designs

Pooled sequencing-based fitness assays are a powerful and widely used approach to quantifying fitness of thousands of genetic variants in parallel. Despite the throughput of such assays, they are prone to biases in fitness estimates, and errors in measurements are typically larger for deleterious fitness effects, relative to neutral effects. In practice, designing pooled fitness assays involves tradeoffs between the number of timepoints, the sequencing depth, and other parameters to gain as much information as possible within a feasible experiment. Here, we combined simulations and reanalysis of an existing experimental dataset to explore how assay parameters impact measurements of near-neutral and deleterious fitness effects using a standard fitness estimator. We found that sequencing multiple timepoints at relatively modest depth improved estimates of near-neutral fitness effects, but systematically biased measurements of deleterious effects. We showed that a fixed total number of reads, deeper sequencing at fewer timepoints improved resolution of deleterious fitness effects. Our results highlight a tradeoff between measurement of deleterious and near-neutral effect sizes for a fixed amount of data and suggest that fitness assay design should be tuned for fitness effects that are relevant to the specific biological question.

Barcoding Populations of Pseudomonas fluorescens SBW25

In recent years, evolutionary biologists have developed an increasing interest in the use of barcoding strategies to study eco-evolutionary dynamics of lineages within evolving populations and communities. Although barcoded populations can deliver unprecedented insight into evolutionary change, barcoding microbes presents specific technical challenges. Here, strategies are described for barcoding populations of the model bacterium Pseudomonas fluorescens SBW25, including the design and cloning of barcoded regions, preparation of libraries for amplicon sequencing, and quantification of resulting barcoded lineages. In so doing, we hope to aid the design and implementation of barcoding methodologies in a broad range of model and non-model organisms.

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.

Identifying Targets of Selection in Laboratory Evolution Experiments

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

Evolutionary repeatability of emergent properties of ecological communities

Most species belong to ecological communities where their interactions give rise to emergent community-level properties, such as diversity and productivity. Understanding and predicting how these properties change over time has been a major goal in ecology, with important practical implications for sustainability and human health. Less attention has been paid to the fact that community-level properties can also change because member species evolve. Yet, our ability to predict long-term eco-evolutionary dynamics hinges on how repeatably community-level properties change as a result of species evolution. Here, we review studies of evolution of both natural and experimental communities and make the case that community-level properties at least sometimes evolve repeatably. We discuss challenges faced in investigations of evolutionary repeatability. In particular, only a handful of studies enable us to quantify repeatability. We argue that quantifying repeatability at the community level is critical for approaching what we see as three major open questions in the field: (i) Is the observed degree of repeatability surprising? (ii) How is evolutionary repeatability at the community level related to repeatability at the level of traits of member species? (iii) What factors affect repeatability? We outline some theoretical and empirical approaches to addressing these questions. Advances in these directions will not only enrich our basic understanding of evolution and ecology but will also help us predict eco-evolutionary dynamics. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Mutational scans reveal differential evolvability of Drosophila promoters and enhancers

Rapid enhancer and slow promoter evolution have been demonstrated through comparative genomics. However, it is not clear how this information is encoded genetically and if this can be used to place evolution in a predictive context. Part of the challenge is that our understanding of the potential for regulatory evolution is biased primarily toward natural variation or limited experimental perturbations. Here, to explore the evolutionary capacity of promoter variation, we surveyed an unbiased mutation library for three promoters in Drosophila melanogaster. We found that mutations in promoters had limited to no effect on spatial patterns of gene expression. Compared to developmental enhancers, promoters are more robust to mutations and have more access to mutations that can increase gene expression, suggesting that their low activity might be a result of selection. Consistent with these observations, increasing the promoter activity at the endogenous locus of shavenbaby led to increased transcription yet limited phenotypic changes. Taken together, developmental promoters may encode robust transcriptional outputs allowing evolvability through the integration of diverse developmental enhancers.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

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.

Biomarkers of aging

Aging biomarkers are a combination of biological parameters to (i) assess age-related changes, (ii) track the physiological aging process, and (iii) predict the transition into a pathological status. Although a broad spectrum of aging biomarkers has been developed, their potential uses and limitations remain poorly characterized. An immediate goal of biomarkers is to help us answer the following three fundamental questions in aging research: How old are we? Why do we get old? And how can we age slower? This review aims to address this need. Here, we summarize our current knowledge of biomarkers developed for cellular, organ, and organismal levels of aging, comprising six pillars: physiological characteristics, medical imaging, histological features, cellular alterations, molecular changes, and secretory factors. To fulfill all these requisites, we propose that aging biomarkers should qualify for being specific, systemic, and clinically relevant.

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.

Novel evolved Yarrowia lipolytica strains for enhanced growth and lipid content under high concentrations of crude glycerol

BACKGROUND

Yarrowia lipolytica is a well-studied oleaginous yeast known for its ability to accumulate and store intracellular lipids, while growing on diverse, non-conventional substrates. Amongst them, crude glycerol, a low-cost by-product of the biodiesel industry, appears to be an interesting option for scaling up a sustainable single-cell oil production process. Adaptive laboratory evolution (ALE) is a powerful tool to force metabolic adaptations endowing tolerance to stressful environmental conditions, generating superior phenotypes with industrial relevance.

RESULTS

Y. lipolytica MUCL 28849 underwent ALE in a synthetic medium with increasing concentration of pure or crude glycerol as a stressing factor (9-20% v/v) for 520 generations. In one case of pure glycerol, chemical mutagenesis with ethyl methanesulfonate (EMS) was applied prior to ALE. Growth profile, biomass production and lipid content of 660 evolved strains (EVS), revealed 5 superior isolates; exhibiting from 1.9 to 3.6-fold increase of dry biomass and from 1.1 to 1.6-fold increase of lipid concentration compared to the parental strain, when grown in 15% v/v crude glycerol. NGS for differential gene expression analysis, showed induced expression in all EVS affecting nucleosomal structure and regulation of transcription. As strains differentiated, further changes accumulated in membrane transport and protein transport processes. Genes involved in glycerol catabolism and triacylglycerol biosynthesis were overexpressed in two EVS. Mismatches and gaps in the expressed sequences identified altered splicing and mutations in the EVS, with most of them, affecting different components of septin ring formation in the budding process. The selected YLE155 EVS, used for scale-up cultivation in a 3L benchtop bioreactor with 20% v/v crude glycerol, achieved extended exponential phase, twofold increase of dry biomass and lipid yields at 48 h, while citric acid secretion and glycerol consumption rates were 40% and 50% lower, respectively, compared to the parental strain, after 24 h of cultivation.

CONCLUSION

ALE and EMS-ALE under increasing concentrations of pure or crude glycerol generated novel Y. lipolytica strains with enhanced biomass and lipid content. Differential gene expression analysis and scale-up of YLE155, illustrated the potential of the evolved strains to serve as suitable "chassis" for rational engineering approaches towards both increased lipid accumulation, and production of high-added value compounds, through efficient utilization of crude glycerol.

Adaptation to Overflow Metabolism by Mutations That Impair tRNA Modification in Experimentally Evolved Bacteria

When microbes grow in foreign nutritional environments, selection may enrich mutations in unexpected pathways connecting growth and homeostasis. An evolution experiment designed to identify beneficial mutations in Burkholderia cenocepacia captured six independent nonsynonymous substitutions in the essential gene tilS, which modifies tRNA^Ile2 by adding a lysine to the anticodon for faithful AUA recognition. Further, five additional mutants acquired mutations in tRNA^Ile2, which strongly suggests that disrupting the TilS-tRNA^Ile2 interaction was subject to strong positive selection. Mutated TilS incurred greatly reduced enzymatic function but retained capacity for tRNA^Ile2 binding. However, both mutant sets outcompeted the wild type by decreasing the lag phase duration by ~3.5 h. We hypothesized that lysine demand could underlie fitness in the experimental conditions. As predicted, supplemental lysine complemented the ancestral fitness deficit, but so did the additions of several other amino acids. Mutant fitness advantages were also specific to rapid growth on galactose using oxidative overflow metabolism that generates redox imbalance, not resources favoring more balanced metabolism. Remarkably, 13 tilS mutations also evolved in the long-term evolution experiment with Escherichia coli, including four fixed mutations. These results suggest that TilS or unknown binding partners contribute to improved growth under conditions of rapid sugar oxidation at the predicted expense of translational accuracy. IMPORTANCE There is growing evidence that the fundamental components of protein translation can play multiple roles in maintaining cellular homeostasis. Enzymes that interact with transfer RNAs not only ensure faithful decoding of the genetic code but also help signal the metabolic state by reacting to imbalances in essential building blocks like free amino acids and cofactors. Here, we present evidence of a secondary function for the essential enzyme TilS, whose only prior known function is to modify tRNA^Ile(CAU) to ensure accurate translation. Multiple nonsynonymous substitutions in tilS, as well as its cognate tRNA, were selected in evolution experiments favoring rapid, redox-imbalanced growth. These mutations alone decreased lag phase and created a competitive advantage, but at the expense of most primary enzyme function. These results imply that TilS interacts with other factors related to the timing of exponential growth and that tRNA-modifying enzymes may serve multiple roles in monitoring metabolic health.

Optimizing the strain engineering process for industrial-scale production of bio-based molecules

Biomanufacturing could contribute as much as ${\$}$30 trillion to the global economy by 2030. However, the success of the growing bioeconomy depends on our ability to manufacture high-performing strains in a time- and cost-effective manner. The Design-Build-Test-Learn (DBTL) framework has proven to be an effective strain engineering approach. Significant improvements have been made in genome engineering, genotyping, and phenotyping throughput over the last couple of decades that have greatly accelerated the DBTL cycles. However, to achieve a radical reduction in strain development time and cost, we need to look at the strain engineering process through a lens of optimizing the whole cycle, as opposed to simply increasing throughput at each stage. We propose an approach that integrates all 4 stages of the DBTL cycle and takes advantage of the advances in computational design, high-throughput genome engineering, and phenotyping methods, as well as machine learning tools for making predictions about strain scale-up performance. In this perspective, we discuss the challenges of industrial strain engineering, outline the best approaches to overcoming these challenges, and showcase examples of successful strain engineering projects for production of heterologous proteins, amino acids, and small molecules, as well as improving tolerance, fitness, and de-risking the scale-up of industrial strains.

Quantifying the local adaptive landscape of a nascent bacterial community

The fitness effects of all possible mutations available to an organism largely shape the dynamics of evolutionary adaptation. Yet, whether and how this adaptive landscape changes over evolutionary times, especially upon ecological diversification and changes in community composition, remains poorly understood. We sought to fill this gap by analyzing a stable community of two closely related ecotypes ("L" and "S") shortly after they emerged within the E. coli Long-Term Evolution Experiment (LTEE). We engineered genome-wide barcoded transposon libraries to measure the invasion fitness effects of all possible gene knockouts in the coexisting strains as well as their ancestor, for many different, ecologically relevant conditions. We find consistent statistical patterns of fitness effect variation across both genetic background and community composition, despite the idiosyncratic behavior of individual knockouts. Additionally, fitness effects are correlated with evolutionary outcomes for a number of conditions, possibly revealing shifting patterns of adaptation. Together, our results reveal how ecological and epistatic effects combine to shape the adaptive landscape in a nascent ecological community.

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.

Mutualism-enhancing mutations dominate early adaptation in a two-species microbial community

Species interactions drive evolution while evolution shapes these interactions. The resulting eco-evolutionary dynamics and their repeatability depend on how adaptive mutations available to community members affect fitness and ecologically relevant traits. However, the diversity of adaptive mutations is not well characterized, and we do not know how this diversity is affected by the ecological milieu. Here we use barcode lineage tracking to address this question in a community of yeast Saccharomyces cerevisiae and alga Chlamydomonas reinhardtii that have a net commensal relationship that results from a balance between competitive and mutualistic interactions. We find that yeast has access to many adaptive mutations with diverse ecological consequences, in particular those that increase and reduce the yields of both species. The presence of the alga does not change which mutations are adaptive in yeast (that is, there is no fitness trade-off for yeast between growing alone or with alga), but rather shifts selection to favour yeast mutants that increase the yields of both species and make the mutualism stronger. Thus, in the presence of the alga, adaptative mutations contending for fixation in yeast are more likely to enhance the mutualism, even though cooperativity is not directly favoured by natural selection in our system. Our results demonstrate that ecological interactions not only alter the trajectory of evolution but also dictate its repeatability; in particular, weak mutualisms can repeatably evolve to become stronger.

On the hidden temporal dynamics of plant adaptation

Adaptation to a wide range of environments is a major driver of plant diversity. It is now possible to catalog millions of potential adaptive genomic differences segregating between environments within a plant species in a single experiment. Understanding which of these changes contributes to adaptive phenotypic divergence between plant populations is a major goal of evolutionary biologists and crop breeders. In this review, we briefly highlight the approaches frequently used to understand the genetic basis of adaptive phenotypes in plants, and we discuss some of the limitations of these methods. We propose that direct observation of the process of adaptation using multigenerational studies and whole genome sequencing is a crucial missing component of recent studies of plant adaptation because it complements several shortcomings of sampling-based techniques.

The Role of the Environment in Horizontal Gene Transfer

Gene-by-environment interactions play a crucial role in horizontal gene transfer by affecting how the transferred genes alter host fitness. However, how the environment modulates the fitness effect of transferred genes has not been tested systematically in an experimental study. We adapted a high-throughput technique for obtaining very precise estimates of bacterial fitness, in order to measure the fitness effects of 44 orthologs transferred from Salmonella Typhimurium to Escherichia coli in six physiologically relevant environments. We found that the fitness effects of individual genes were highly dependent on the environment, while the distributions of fitness effects across genes were not, with all tested environments resulting in distributions of same shape and spread. Furthermore, the extent to which the fitness effects of a gene varied between environments depended on the average fitness effect of that gene across all environments, with nearly neutral and nearly lethal genes having more consistent fitness effects across all environments compared to deleterious genes. Put together, our results reveal the unpredictable nature of how environmental conditions impact the fitness effects of each individual gene. At the same time, distributions of fitness effects across environments exhibit consistent features, pointing to the generalizability of factors that shape horizontal gene transfer of orthologous genes.

Using colony size to measure fitness in Saccharomyces cerevisiae

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

The landscape of transcriptional and translational changes over 22 years of bacterial adaptation

Organisms can adapt to an environment by taking multiple mutational paths. This redundancy at the genetic level, where many mutations have similar phenotypic and fitness effects, can make untangling the molecular mechanisms of complex adaptations difficult. Here, we use the Escherichia coli long-term evolution experiment (LTEE) as a model to address this challenge. To understand how different genomic changes could lead to parallel fitness gains, we characterize the landscape of transcriptional and translational changes across 12 replicate populations evolving in parallel for 50,000 generations. By quantifying absolute changes in mRNA abundances, we show that not only do all evolved lines have more mRNAs but that this increase in mRNA abundance scales with cell size. We also find that despite few shared mutations at the genetic level, clones from replicate populations in the LTEE are remarkably similar in their gene expression patterns at both the transcriptional and translational levels. Furthermore, we show that the majority of the expression changes are due to changes at the transcriptional level with very few translational changes. Finally, we show how mutations in transcriptional regulators lead to consistent and parallel changes in the expression levels of downstream genes. These results deepen our understanding of the molecular mechanisms underlying complex adaptations and provide insights into the repeatability of evolution.

Loss-of-function mutation survey revealed that genes with background-dependent fitness are rare and functionally related in yeast

In different individuals, the same mutation can lead to different phenotypes due to genetic background effects. This is commonly observed in various systems, including many human diseases. While isolated examples of such background effects have been observed, a systematic view across a large number of individuals is still lacking. Here, we surveyed genetic background effects associated with gene loss-of-function mutations across a population of natural isolates of the yeast Saccharomyces cerevisiae. We found that ∼15% of genes can display a background-dependent fitness change. Genes related to mitochondrial functions are significantly overrepresented, and showed reversed patterns of fitness gain or loss with genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport, suggesting a potential functional rewiring.

In natural populations, the same mutation can lead to different phenotypic outcomes due to the genetic variation that exists among individuals. Such genetic background effects are commonly observed, including in the context of many human diseases. However, systematic characterization of these effects at the species level is still lacking to date. Here, we sought to comprehensively survey background-dependent traits associated with gene loss-of-function (LoF) mutations in 39 natural isolates of Saccharomyces cerevisiae using a transposon saturation strategy. By analyzing the modeled fitness variability of a total of 4,469 genes, we found that 15% of them, when impacted by a LoF mutation, exhibited a significant gain- or loss-of-fitness phenotype in certain natural isolates compared with the reference strain S288C. Out of these 632 genes with predicted background-dependent fitness effects, around 2/3 impact multiple backgrounds with a gradient of predicted fitness change while 1/3 are specific to a single genetic background. Genes related to mitochondrial function are significantly overrepresented in the set of genes showing a continuous variation and display a potential functional rewiring with other genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport. Such rewiring effects are likely modulated by both the genetic background and the environment. While background-specific cases are rare and span diverse cellular processes, they can be functionally related at the individual level. All genes with background-dependent fitness effects tend to have an intermediate connectivity in the global genetic interaction network and have shown relaxed selection pressure at the population level, highlighting their potential evolutionary characteristics.

Inheritance of somatic mutations by animal offspring

Since 1892, it has been widely assumed that somatic mutations are evolutionarily irrelevant in animals because they cannot be inherited by offspring. However, some nonbilaterians segregate the soma and germline late in development or never, leaving the evolutionary fate of their somatic mutations unknown. By investigating uni- and biparental reproduction in the coral Acropora palmata (Cnidaria, Anthozoa), we found that uniparental, meiotic offspring harbored 50% of the 268 somatic mutations present in their parent. Thus, somatic mutations accumulated in adult coral animals, entered the germline, and were passed on to swimming larvae that grew into healthy juvenile corals. In this way, somatic mutations can increase allelic diversity and facilitate adaptation across habitats and generations in animals.

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

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

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

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