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

Evolution of evolvability in rapidly adapting populations

Mutations can alter the short-term fitness of an organism, as well as the rates and benefits of future mutations. While numerous examples of these evolvability modifiers have been observed in rapidly adapting microbial populations, existing theory struggles to predict when they will be favoured by natural selection. Here we develop a mathematical framework for predicting the fates of genetic variants that modify the rates and benefits of future mutations in linked genomic regions. We derive analytical expressions showing how the fixation probabilities of these variants depend on the size of the population and the diversity of competing mutations. We find that competition between linked mutations can dramatically enhance selection for modifiers that increase the benefits of future mutations, even when they impose a strong direct cost on fitness. However, we also find that modest direct benefits can be sufficient to drive evolutionary dead ends to fixation. Our results suggest that subtle differences in evolvability could play an important role in shaping the long-term success of genetic variants in rapidly evolving microbial populations.

Environment-independent distribution of mutational effects emerges from microscopic epistasis

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

Changing fitness effects of mutations through long-term bacterial evolution

The distribution of fitness effects of new mutations shapes evolution, but it is challenging to observe how it changes as organisms adapt. Using Escherichia coli lineages spanning 50,000 generations of evolution, we quantify the fitness effects of insertion mutations in every gene. Macroscopically, the fraction of deleterious mutations changed little over time whereas the beneficial tail declined sharply, approaching an exponential distribution. Microscopically, changes in individual gene essentiality and deleterious effects often occurred in parallel; altered essentiality is only partly explained by structural variation. The identity and effect sizes of beneficial mutations changed rapidly over time, but many targets of selection remained predictable because of the importance of loss-of-function mutations. Taken together, these results reveal the dynamic-but statistically predictable-nature of mutational fitness effects.

Gene birth in a model of non-genic adaptation

BACKGROUND

Over evolutionary timescales, genomic loci can switch between functional and non-functional states through processes such as pseudogenization and de novo gene birth. Particularly, de novo gene birth is a widespread process, and many examples continue to be discovered across diverse evolutionary lineages. However, the general mechanisms that lead to functionalization are poorly understood, and estimated rates of de novo gene birth remain contentious. Here, we address this problem within a model that takes into account mutations and structural variation, allowing us to estimate the likelihood of emergence of new functions at non-functional loci.

RESULTS

Assuming biologically reasonable mutation rates and mutational effects, we find that functionalization of non-genic loci requires the realization of strict conditions. This is in line with the observation that most de novo genes are localized to the vicinity of established genes. Our model also provides an explanation for the empirical observation that emerging proto-genes are often lost despite showing signs of adaptation.

CONCLUSIONS

Our work elucidates the properties of non-genic loci that make them fertile for adaptation, and our results offer mechanistic insights into the process of de novo gene birth.

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.

A null model for the distribution of fitness effects of mutations

The distribution of fitness effects (DFE) of new mutations is key to our understanding of many evolutionary processes. Theoreticians have developed several models to help understand the patterns seen in empirical DFEs. Many such models reproduce the broad patterns seen in empirical DFEs but these models often rely on structural assumptions that cannot be tested empirically. Here, we investigate how much of the underlying "microscopic" biological processes involved in the mapping of new mutations to fitness can be inferred from "macroscopic" observations of the DFE. We develop a null model by generating random genotype-to-fitness maps and show that the null DFE is that with the largest possible information entropy. We further show that, subject to one simple constraint, this null DFE is a Gompertz distribution. Finally, we illustrate how the predictions of this null DFE match empirically measured DFEs from several datasets, as well as DFEs simulated from Fisher's geometric model. This suggests that a match between models and empirical data is often not a very strong indication of the mechanisms underlying the mapping of mutation to fitness.

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.

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.

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.

Epistasis and evolution: recent advances and an outlook for prediction

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

Transient mutation bias increases the predictability of evolution on an empirical genotype-phenotype landscape

Predicting how a population will likely navigate a genotype-phenotype landscape requires consideration of selection in combination with mutation bias, which can skew the likelihood of following a particular trajectory. Strong and persistent directional selection can drive populations to ascend toward a peak. However, with a greater number of peaks and more routes to reach them, adaptation inevitably becomes less predictable. Transient mutation bias, which operates only on one mutational step, can influence landscape navigability by biasing the mutational trajectory early in the adaptive walk. This sets an evolving population upon a particular path, constraining the number of accessible routes and making certain peaks and routes more likely to be realized than others. In this work, we employ a model system to investigate whether such transient mutation bias can reliably and predictably place populations on a mutational trajectory to the strongest selective phenotype or usher populations to realize inferior phenotypic outcomes. For this we use motile mutants evolved from ancestrally non-motile variants of the microbe Pseudomonas fluorescens SBW25, of which one trajectory exhibits significant mutation bias. Using this system, we elucidate an empirical genotype-phenotype landscape, where the hill-climbing process represents increasing strength of the motility phenotype, to reveal that transient mutation bias can facilitate rapid and predictable ascension to the strongest observed phenotype in place of equivalent and inferior trajectories. 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.

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.

Mutational robustness changes during long-term adaptation in laboratory budding yeast populations

As an adapting population traverses the fitness landscape, its local neighborhood (i.e., the collection of fitness effects of single-step mutations) can change shape because of interactions with mutations acquired during evolution. These changes to the distribution of fitness effects can affect both the rate of adaptation and the accumulation of deleterious mutations. However, while numerous models of fitness landscapes have been proposed in the literature, empirical data on how this distribution changes during evolution remains limited. In this study, we directly measure how the fitness landscape neighborhood changes during laboratory adaptation. Using a barcode-based mutagenesis system, we measure the fitness effects of 91 specific gene disruption mutations in genetic backgrounds spanning 8000–10,000 generations of evolution in two constant environments. We find that the mean of the distribution of fitness effects decreases in one environment, indicating a reduction in mutational robustness, but does not change in the other. We show that these distribution-level patterns result from differences in the relative frequency of certain patterns of epistasis at the level of individual mutations, including fitness-correlated and idiosyncratic epistasis.

Influence of genomic structural variations and nutritional conditions on the emergence of quorum sensing-dependent gene regulation defects in Burkholderia glumae

Bacteria often change their genetic and physiological traits to survive in harsh environments. To determine whether, in various strains of Burkholderia glumae, genomic diversity is associated with the ability to adapt to ever-changing environments, whole genomes of 44 isolates from different hosts and regions were analyzed. Whole-genome phylogenetic analysis of the 44 isolates revealed six clusters and two divisions. While all isolates possessed chromosomes 1 and 2, strains BGR80S and BGR81S had one chromosome resulting from the merging of the two chromosomes. Upon comparison of genomic structures to the prototype BGR1, inversions, deletions, and rearrangements were found within or between chromosomes 1 and/or 2 in the other isolates. When three isolates—BGR80S, BGR15S, and BGR21S, representing clusters III, IV, and VI, respectively—were grown in Luria-Bertani medium, spontaneous null mutations were identified in qsmR encoding a quorum-sensing master regulator. Six days after subculture, qsmR mutants were found at detectable frequencies in BGR15S and BGR21S, and reached approximately 40% at 8 days after subculture. However, the qsmR mutants appeared 2 days after subculture in BGR80S and dominated the population, reaching almost 80%. No qsmR mutant was detected at detectable frequency in BGR1 or BGR13S. The spontaneous qsmR mutants outcompeted their parental strains in the co-culture. Daily addition of glucose or casamino acids to the batch cultures of BGR80S delayed emergence of qsmR mutants and significantly reduced their incidence. These results indicate that spontaneous qsmR mutations are correlated with genomic structures and nutritional conditions.

The teenage years of yeast population genomics trace history, admixing and getting wilder

Population genomics studies the evolutionary processes that shape intraspecies genetic variations. In this review, I explore the insights into yeast-population genomics that have emerged from recent advances in sequencing. Genomes of the model Saccharomyces cerevisiae and many new yeast species from around the world are being used to address various aspects of population biology, including geographical origin, the level of introgression, domestication signatures, and outcrossing frequency. New long-read sequencing has enabled a greater capacity to quantify these variations at a finer resolution from complete de novo genomes at the population scale to phasing subgenomes of different origins. These resources provide a platform to dissect the relationship between phenotypes across environmental niches.

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

The population genetics of collateral resistance and sensitivity

Resistance mutations against one drug can elicit collateral sensitivity against other drugs. Multi-drug treatments exploiting such trade-offs can help slow down the evolution of resistance. However, if mutations with diverse collateral effects are available, a treated population may evolve either collateral sensitivity or collateral resistance. How to design treatments robust to such uncertainty is unclear. We show that many resistance mutations in Escherichia coli against various antibiotics indeed have diverse collateral effects. We propose to characterize such diversity with a joint distribution of fitness effects (JDFE) and develop a theory for describing and predicting collateral evolution based on simple statistics of the JDFE. We show how to robustly rank drug pairs to minimize the risk of collateral resistance and how to estimate JDFEs. In addition to practical applications, these results have implications for our understanding of evolution in variable environments.

Drugs known as antibiotics are the main treatment for most serious infections caused by bacteria. However, many bacteria are acquiring genetic mutations that make them resistant to the effects of one or more types of antibiotics, making them harder to eliminate.

One way to tackle drug-resistant bacteria is to develop new types of antibiotics; however, in recent years, the rate at which new antibiotics have become available has been dwindling. Using two or more existing drugs, one after another, can also be an effective way to eliminate resistant bacteria. The success of any such ‘multi-drug’ treatment lies in being able to predict whether mutations that make the bacteria resistant to one drug simultaneously make it sensitive to another, a phenomenon known as collateral sensitivity.

Different resistance mutations may have different collateral effects: some may increase the bacteria’s sensitivity to the second drug, while others might make the bacteria more resistant. However, it is currently unclear how to design robust multi-drug treatments that take this diversity of collateral effects into account.

Here, Ardell and Kryazhimskiy used a concept called JDFE (short for the joint distribution of fitness effects) to describe the diversity of collateral effects in a population of bacteria exposed to a single drug. This information was then used to mathematically model how collateral effects evolved in the population over time.

Ardell and Kryazhimskiy showed that this approach can predict how likely a population is to become collaterally sensitive or collaterally resistant to a second antibiotic. Drug pairs can then be ranked according to the risk of collateral resistance emerging, so long as information on the variety of resistance mutations available to the bacteria are included in the model.

Each year, more than 700,000 people die from infections caused by bacteria that are resistant to one or more antibiotics. The findings of Ardell and Kryazhimskiy may eventually help clinicians design multi-drug treatments that effectively eliminate bacterial infections and help to prevent more bacteria from evolving resistance to antibiotics. However, to achieve this goal, more research is needed to fully understand the range collateral effects caused by resistance mutations.

Joint effects of genes underlying a temperature specialization tradeoff in yeast

A central goal of evolutionary genetics is to understand, at the molecular level, how organisms adapt to their environments. For a given trait, the answer often involves the acquisition of variants at unlinked sites across the genome. Genomic methods have achieved landmark successes in pinpointing these adaptive loci. To figure out how a suite of adaptive alleles work together, and to what extent they can reconstitute the phenotype of interest, requires their transfer into an exogenous background. We studied the joint effect of adaptive, gain-of-function thermotolerance alleles at eight unlinked genes from Saccharomyces cerevisiae, when introduced into a thermosensitive sister species, S. paradoxus. Although the loci damped each other’s beneficial impact (that is, they were subject to negative epistasis), most boosted high-temperature growth alone and in combination, and none was deleterious. The complete set of eight genes was sufficient to confer ~15% of the S. cerevisiae thermotolerance phenotype in the S. paradoxus background. The same loci also contributed to a heretofore unknown advantage in cold growth by S. paradoxus. Together, our data establish temperature resistance in yeasts as a model case of a genetically complex evolutionary tradeoff, which can be partly reconstituted from the sequential assembly of unlinked underlying loci.

Organisms adapt to threats in the environment by acquiring DNA sequence variants that tweak traits to improve fitness. Experimental studies of this process have proven to be a particular challenge when they involve manipulation of a suite of genes, all on different chromosomes. We set out to understand how so many loci could work together to confer a trait. We used as a model system eight genes that govern the ability of the unicellular yeast Saccharomyces cerevisiae to grow at high temperature. We introduced these variant loci stepwise into a non-thermotolerant sister species, and found that the more S. cerevisiae alleles we added, the better the phenotype. We saw no evidence for toxic interactions between the genes as they were combined. We also used the eight-fold transgenic to dissect the biological mechanism of thermotolerance. And we discovered a tradeoff: the same alleles that boosted growth at high temperature eroded the organism’s ability to deal with cold conditions. These results serve as a case study of modular construction of a trait from nature, by assembling the genes together in one genome.