Environment changes epistasis to alter trade‐offs along alternative evolutionary paths

Data-driven discovery of the interplay between genetic and environmental factors in bacterial growth

A complex interplay of genetic and environmental factors influences bacterial growth. Understanding these interactions is crucial for insights into complex living systems. This study employs a data-driven approach to uncover the principles governing bacterial growth changes due to genetic and environmental variation. A pilot survey is conducted across 115 Escherichia coli strains and 135 synthetic media comprising 45 chemicals, generating 13,944 growth profiles. Machine learning analyzes this dataset to predict the chemicals' priorities for bacterial growth. The primary gene-chemical networks are structured hierarchically, with glucose playing a pivotal role. Offset in bacterial growth changes is frequently observed across 1,445,840 combinations of strains and media, with its magnitude correlating to individual alterations in strains or media. This counterbalance in the gene-chemical interplay is supposed to be a general feature beneficial for bacterial population growth.

Competition-driven eco-evolutionary feedback reshapes bacteriophage lambda's fitness landscape and enables speciation

A major challenge in evolutionary biology is explaining how populations navigate rugged fitness landscapes without getting trapped on local optima. One idea illustrated by adaptive dynamics theory is that as populations adapt, their newly enhanced capacities to exploit resources alter fitness payoffs and restructure the landscape in ways that promote speciation by opening new adaptive pathways. While there have been indirect tests of this theory, to our knowledge none have measured how fitness landscapes deform during adaptation, or test whether these shifts promote diversification. Here, we achieve this by studying bacteriophage [Formula: see text], a virus that readily speciates into co-existing receptor specialists under controlled laboratory conditions. We use a high-throughput gene editing-phenotyping technology to measure [Formula: see text]'s fitness landscape in the presence of different evolved-[Formula: see text] competitors and find that the fitness effects of individual mutations, and their epistatic interactions, depend on the competitor. Using these empirical data, we simulate [Formula: see text]'s evolution on an unchanging landscape and one that recapitulates how the landscape deforms during evolution. [Formula: see text] heterogeneity only evolves in the shifting landscape regime. This study provides a test of adaptive dynamics, and, more broadly, shows how fitness landscapes dynamically change during adaptation, potentiating phenomena like speciation by opening new adaptive pathways.

Environmental modulation of global epistasis in a drug resistance fitness landscape

Interactions between mutations (epistasis) can add substantial complexity to genotype-phenotype maps, hampering our ability to predict evolution. Yet, recent studies have shown that the fitness effect of a mutation can often be predicted from the fitness of its genetic background using simple, linear relationships. This phenomenon, termed global epistasis, has been leveraged to reconstruct fitness landscapes and infer adaptive trajectories in a wide variety of contexts. However, little attention has been paid to how patterns of global epistasis may be affected by environmental variation, despite this variation frequently being a major driver of evolution. This is particularly relevant for the evolution of drug resistance, where antimicrobial drugs may change the environment faced by pathogens and shape their adaptive trajectories in ways that can be difficult to predict. By analyzing a fitness landscape of four mutations in a gene encoding an essential enzyme of P. falciparum (a parasite cause of malaria), here we show that patterns of global epistasis can be strongly modulated by the concentration of a drug in the environment. Expanding on previous theoretical results, we demonstrate that this modulation can be quantitatively explained by how specific gene-by-gene interactions are modified by drug dose. Importantly, our results highlight the need to incorporate potential environmental variation into the global epistasis framework in order to predict adaptation in dynamic environments.

Have you tried turning it off and on again? Oscillating selection to enhance fitness-landscape traversal in adaptive laboratory evolution experiments

Adaptive Laboratory Evolution (ALE) is a powerful tool for engineering and understanding microbial physiology. ALE relies on the selection and enrichment of mutations that enable survival or faster growth under a selective condition imposed by the experimental setup. Phenotypic fitness landscapes are often underpinned by complex genotypes involving multiple genes, with combinatorial positive and negative effects on fitness. Such genotype relationships result in mutational fitness landscapes with multiple local fitness maxima and valleys. Traversing local maxima to find a global maximum often requires an individual or sub-population of cells to traverse fitness valleys. Traversing involves gaining mutations that are not adaptive for a given local maximum but are necessary to 'peak shift' to another local maximum, or eventually a global maximum. Despite these relatively well understood evolutionary principles, and the combinatorial genotypes that underlie most metabolic phenotypes, the majority of applied ALE experiments are conducted using constant selection pressures. The use of constant pressure can result in populations becoming trapped within local maxima, and often precludes the attainment of optimum phenotypes associated with global maxima. Here, we argue that oscillating selection pressures is an easily accessible mechanism for traversing fitness landscapes in ALE experiments, and provide theoretical and practical frameworks for implementation.

Competition-driven eco-evolutionary feedback reshapes bacteriophage lambda's fitness landscape and enables speciation

A major challenge in evolutionary biology is explaining how populations navigate rugged fitness landscapes without getting trapped on local optima. One idea illustrated by adaptive dynamics theory is that as populations adapt, their newly enhanced capacities to exploit resources alter fitness payoffs and restructure the landscape in ways that promote speciation by opening new adaptive pathways. While there have been indirect tests of this theory, none have measured how fitness landscapes deform during adaptation, or test whether these shifts promote diversification. Here, we achieve this by studying bacteriophage λ, a virus that readily speciates into co-existing receptor specialists under controlled laboratory conditions. We used a high-throughput gene editing-phenotyping technology to measure λ's fitness landscape in the presence of different evolved-λ competitors and found that the fitness effects of individual mutations, and their epistatic interactions, depend on the competitor. Using these empirical data, we simulated λ's evolution on an unchanging landscape and one that recapitulates how the landscape deforms during evolution. λ heterogeneity only evolved in the shifting landscape regime. This study provides a test of adaptive dynamics, and, more broadly, shows how fitness landscapes dynamically change during adaptation, potentiating phenomena like speciation by opening new adaptive pathways.

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.

Epistasis decreases with increasing antibiotic pressure but not temperature

Predicting mutational effects is essential for the control of antibiotic resistance (ABR). Predictions are difficult when there are strong genotype-by-environment (G × E), gene-by-gene (G × G or epistatic) or gene-by-gene-by-environment (G × G × E) interactions. We quantified G × G × E effects in Escherichia coli across environmental gradients. We created intergenic fitness landscapes using gene knock-outs and single-nucleotide ABR mutations previously identified to vary in the extent of G × E effects in our environments of interest. Then, we measured competitive fitness across a complete combinatorial set of temperature and antibiotic dosage gradients. In this way, we assessed the predictability of 15 fitness landscapes across 12 different but related environments. We found G × G interactions and rugged fitness landscapes in the absence of antibiotic, but as antibiotic concentration increased, the fitness effects of ABR genotypes quickly overshadowed those of gene knock-outs, and the landscapes became smoother. Our work reiterates that some single mutants, like those conferring resistance or susceptibility to antibiotics, have consistent effects across genetic backgrounds in stressful environments. Thus, although epistasis may reduce the predictability of evolution in benign environments, evolution may be more predictable in adverse environments. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Impact of population size on early adaptation in rugged fitness landscapes

Owing to stochastic fluctuations arising from finite population size, known as genetic drift, the ability of a population to explore a rugged fitness landscape depends on its size. In the weak mutation regime, while the mean steady-state fitness increases with population size, we find that the height of the first fitness peak encountered when starting from a random genotype displays various behaviours versus population size, even among small and simple rugged landscapes. We show that the accessibility of the different fitness peaks is key to determining whether this height overall increases or decreases with population size. Furthermore, there is often a finite population size that maximizes the height of the first fitness peak encountered when starting from a random genotype. This holds across various classes of model rugged landscapes with sparse peaks, and in some experimental and experimentally inspired ones. Thus, early adaptation in rugged fitness landscapes can be more efficient and predictable for relatively small population sizes than in the large-size limit. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Phage-resistant mutations impact bacteria susceptibility to future phage infections and antibiotic response

Bacteriophage (phage) therapy in combination with antibiotic treatment serves as a potential strategy to overcome the continued rise in antibiotic resistance across bacterial pathogens. Understanding the impacts of evolutionary and ecological processes to the phage-antibiotic-resistance dynamic could advance the development of such combinatorial therapy. We tested whether the acquisition of mutations conferring phage resistance may have antagonistically pleiotropic consequences for antibiotic resistance. First, to determine the robustness of phage resistance across different phage strains, we infected resistant Escherichia coli cultures with phage that were not previously encountered. We found that phage-resistant E. coli mutants that gained resistance to a single phage strain maintain resistance to other phages with overlapping adsorption methods. Mutations underlying the phage-resistant phenotype affects lipopolysaccharide (LPS) structure and/or synthesis. Because LPS is implicated in both phage infection and antibiotic response, we then determined whether phage-resistant trade-offs exist when challenged with different classes of antibiotics. We found that only 1 out of the 4 phage-resistant E. coli mutants yielded trade-offs between phage and antibiotic resistance. Surprisingly, when challenged with novobiocin, we uncovered evidence of synergistic pleiotropy for some mutants allowing for greater antibiotic resistance, even though antibiotic resistance was never selected for. Our results highlight the importance of understanding the role of selective pressures and pleiotropic interactions in the bacterial response to phage-antibiotic combinatorial therapy.

Fitness of evolving bacterial populations is contingent on deep and shallow history but only shallow history creates predictable patterns

Long-term evolution experiments have tested the importance of genetic and environmental factors in influencing evolutionary outcomes. Differences in phylogenetic history, recent adaptation to distinct environments and chance events, all influence the fitness of a population. However, the interplay of these factors on a population's evolutionary potential remains relatively unexplored. We tracked the outcome of 2000 generations of evolution of four natural isolates of Escherichia coli bacteria that were engineered to also create differences in shallow history by adding previously identified mutations selected in a separate long-term experiment. Replicate populations started from each progenitor evolved in four environments. We found that deep and shallow phylogenetic histories both contributed significantly to differences in evolved fitness, though by different amounts in different selection environments. With one exception, chance effects were not significant. Whereas the effect of deep history did not follow any detectable pattern, effects of shallow history followed a pattern of diminishing returns whereby fitter ancestors had smaller fitness increases. These results are consistent with adaptive evolution being contingent on the interaction of several evolutionary forces but demonstrate that the nature of these interactions is not fixed and may not be predictable even when the role of chance is small.

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.

The genetic architecture underlying prey-dependent performance in a microbial predator

Natural selection should favour generalist predators that outperform specialists across all prey types. Two genetic solutions could explain why intraspecific variation in predatory performance is, nonetheless, widespread: mutations beneficial on one prey type are costly on another (antagonistic pleiotropy), or mutational effects are prey-specific, which weakens selection, allowing variation to persist (relaxed selection). To understand the relative importance of these alternatives, we characterised natural variation in predatory performance in the microbial predator Dictyostelium discoideum. We found widespread nontransitive differences among strains in predatory success across different bacterial prey, which can facilitate stain coexistence in multi-prey environments. To understand the genetic basis, we developed methods for high throughput experimental evolution on different prey (REMI-seq). Most mutations (~77%) had prey-specific effects, with very few (~4%) showing antagonistic pleiotropy. This highlights the potential for prey-specific effects to dilute selection, which would inhibit the purging of variation and prevent the emergence of an optimal generalist predator.

What prevents a generalist predator from evolving and outperforming specialist predators? By combing analyses of natural variation with experimental evolution, Stewart et al. suggest that predator variation persists because most mutations have prey-specific effects, which results in relaxed selection

The mutation effect reaction norm (mu-rn) highlights environmentally dependent mutation effects and epistatic interactions

Since the modern synthesis, the fitness effects of mutations and epistasis have been central yet provocative concepts in evolutionary and population genetics. Studies of how the interactions between parcels of genetic information can change as a function of environmental context have added a layer of complexity to these discussions. Here I introduce the "mutation effect reaction norm" (Mu-RN), a new instrument through which one can analyze the phenotypic consequences of mutations and interactions across environmental contexts. It embodies the fusion of measurements of genetic interactions with the reaction norm, a classic depiction of the performance of genotypes across environments. I demonstrate the utility of the Mu-RN through the signature of a "compensatory ratchet" mutation that undermines reverse evolution of antimicrobial resistance. More broadly, I argue that the mutation effect reaction norm may help us resolve the dynamism and unpredictability of evolution, with implications for theoretical biology, genetic modification technology, and public health. This article is protected by copyright. All rights reserved.

Role of bacterial motility in differential resistance mechanisms of silver nanoparticles and silver ions

Unlike conventional antimicrobials, the study of bacterial resistance to silver nanoparticles (AgNPs) remains in its infancy and the mechanism(s) through which it evolves are limited and inconclusive. The central question remains whether bacterial resistance is driven by the AgNPs, released Ag(I) ions or a combination of these and other factors. Here, we show a specific resistance in an Escherichia coli K-12 MG1655 strain to subinhibitory concentrations of AgNPs, and not Ag(I) ions, as indicated by a statistically significant greater-than-twofold increase in the minimum inhibitory concentration occurring after eight repeated passages that was maintained after the AgNPs were removed and reintroduced. Whole-population genome sequencing identified a cusS mutation associated with the heritable resistance that possibly increased silver ion efflux. Finally, we rule out the effect of particle aggregation on resistance and suggest that the mechanism of resistance may be enhanced or mediated by flagellum-based motility.

Evolution of resistance under alternative models of selective interference

The use of multiple pesticides or drugs can lead to a simultaneous selection pressure for resistance alleles at different loci. Models of resistance evolution focus on how this can delay the spread of resistance through a population, but often neglect how this can also reduce the probability that a resistance allele spreads. This neglected factor has been studied in a parallel literature as selective interference. Models of interference use alternative constructions of fitness, where selection coefficients from different loci either add or multiply. Although these are equivalent under weak selection, the two constructions make alternative predictions under the strong selection that characterizes resistance evolution. Here, simulations are used to examine the effects of interference on the probability of fixation and time to fixation of a new and strongly beneficial mutation in the presence of another strongly beneficial allele with variable starting frequency. The results from simulations show a complicated pattern of effects. The key result is that, under multiplicativity, the presence of the strongly beneficial allele leads to a small reduction in the probability of fixation for the new beneficial mutation up to ~10%, and a negligible increase in the average time to fixation up to ~2%, whereas under additivity, the effect is more substantial at up to ~50% for the probability of fixation and ~100% for the average time to fixation. Consequently, the effect of interference is only an important feature of resistance evolution under additivity. Current evidence from studies of experimental evolution provides widespread support for the basic features of additivity, which suggests that interference may afford resistance a different pattern of evolution than other adaptations: rather than the gradual and simultaneous selection of many alleles with small effects, the rapid evolution of resistance may involve the sequential selection of alleles with large effects.

The somatic molecular evolution of cancer: Mutation, selection, and epistasis

Cancer progression has been attributed to somatic changes in single-nucleotide variants, copy-number aberrations, loss of heterozygosity, chromosomal instability, epistatic interactions, and the tumor microenvironment. It is not entirely clear which of these changes are essential and which are ancillary to cancer. The dynamic nature of cancer evolution in a patient can be illuminated using several concepts and tools from classical evolutionary biology. Neutral mutation rates in cancer cells are calculable from genomic data such as synonymous mutations, and selective pressures are calculable from rates of fixation occurring beyond the expectation by neutral mutation and drift. However, these cancer effect sizes of mutations are complicated by epistatic interactions that can determine the likely sequence of gene mutations. In turn, longitudinal phylogenetic analyses of somatic cancer progression offer an opportunity to identify key moments in cancer evolution, relating the timing of driver mutations to corresponding landmarks in the clinical timeline. These analyses reveal temporal aspects of genetic and phenotypic change during tumorigenesis and across clinical timescales. Using a related framework, clonal deconvolution, physical locations of clones, and their phylogenetic relations can be used to infer tumor migration histories. Additionally, genetic interactions with the tumor microenvironment can be analyzed with longstanding approaches applied to organismal genotype-by-environment interactions. Fitness landscapes for cancer evolution relating to genotype, phenotype, and environment could enable more accurate, personalized therapeutic strategies. An understanding of the trajectories underlying the evolution of neoplasms, primary, and metastatic tumors promises fundamental advances toward accurate and personalized predictions of therapeutic response.

Extreme Antagonism Arising from Gene-Environment Interactions

Antagonistic interactions in biological systems, which occur when one perturbation blunts the effect of another, are typically interpreted as evidence that the two perturbations impact the same cellular pathway or function. Yet, this interpretation ignores extreme antagonistic interactions wherein an otherwise deleterious perturbation compensates for the function lost because of a prior perturbation. Here, we report on gene-environment interactions involving genetic mutations that are deleterious in a permissive environment but beneficial in a specific environment that restricts growth. These extreme antagonistic interactions constitute gene-environment analogs of synthetic rescues previously observed for gene-gene interactions. Our approach uses two independent adaptive evolution steps to address the lack of experimental methods to systematically identify such extreme interactions. We apply the approach to Escherichia coli by successively adapting it to defined glucose media without and with the antibiotic rifampicin. The approach identified multiple mutations that are beneficial in the presence of rifampicin and deleterious in its absence. The analysis of transcription shows that the antagonistic adaptive mutations repress a stringent response-like transcriptional program, whereas nonantagonistic mutations have an opposite transcriptional profile. Our approach represents a step toward the systematic characterization of extreme antagonistic gene-drug interactions, which can be used to identify targets to select against antibiotic resistance.