Pervasive,yet idiosyncratic,epistatic pleiotropy during adaptation in a behaviourally complex microbe

Understanding how multiple mutations interact to jointly impact multiple ecologically important traits is critical for creating a robust picture of organismal fitness and the process of adaptation. However, this is complicated by both environmental heterogeneity and the complexity of genotype‐to‐phenotype relationships generated by pleiotropy and epistasis. Moreover, little is known about how pleiotropic and epistatic relationships themselves change over evolutionary time. The soil bacterium Myxococcus xanthus employs several distinct social traits across a range of environments. Here, we use an experimental lineage of M. xanthus that evolved a novel form of social motility to address how interactions between epistasis and pleiotropy evolve. Specifically, we test how mutations accumulated during selection on soft agar pleiotropically affect several other social traits (hard agar motility, predation and spore production). Relationships between changes in swarming rate in the selective environment and the four other traits varied greatly over time in both direction and magnitude, both across timescales of the entire evolutionary lineage and individual evolutionary time steps. We also tested how a previously defined epistatic interaction is pleiotropically expressed across these traits. We found that phenotypic effects of this epistatic interaction were highly correlated between soft and hard agar motility, but were uncorrelated between soft agar motility and predation, and inversely correlated between soft agar motility and spore production. Our results show that ‘epistatic pleiotropy’ varied greatly in magnitude, and often even in sign, across traits and over time, highlighting the necessity of simultaneously considering the interacting complexities of pleiotropy and epistasis when studying the process of adaptation.     

EPISTATIC INTERACTIONS DETERMINE THE MUTATIONAL PATHWAYS AND COEXISTENCE OF LINEAGES IN CLONAL ESCHERICHIA COLI POPULATIONS

Understanding how diversity emerges in a single niche is not fully understood. Rugged fitness landscapes and epistasis between beneficial mutations could explain coexistence among emerging lineages. To provide an experimental test of this notion, we investigated epistasis among four pleiotropic mutations in rpoS, mglD, malT, and hfq present in two coexisting lineages that repeatedly fixed in experimental populations of Escherichia coli. The mutations were transferred into the ancestral background individually or in combination of double or triple alleles. The combined competitive fitness of two or three beneficial mutations from the same lineage was consistently lower than the sum of the competitive fitness of single mutants—a clear indication of negative epistasis within lineages. We also found sign epistasis (i.e., the combined fitness of two beneficial mutations lower than the ancestor), not only from two different lineages (i.e., hfq and rpoS) but also from the same lineage (i.e., mglD and malT). The sign epistasis between loci of different lineages indeed indicated a rugged fitness landscape, providing an epistatic explanation for the coexistence of distinct rpoS and hfq lineages in evolving populations. The negative and sign epistasis between beneficial mutations within the same lineage can further explain the order of mutation acquisition.     

A SHIFT FROM MAGNITUDE TO SIGN EPISTASIS DURING ADAPTIVE EVOLUTION OF A BACTERIAL SOCIAL TRAIT

Although the importance of epistasis in evolution has long been recognized, remarkably little is known about the processes by which epistatic interactions evolve in real time in specific biological systems. Here, we have characterized how the epistatic fitness relationship between a social gene and an adapting genome changes radically over a short evolutionary time frame in the social bacterium Myxococcus xanthus. We show that a highly beneficial effect of this social gene in the ancestral genome is gradually reduced—and ultimately reversed into a deleterious effect—over the course of an experimental adaptive trajectory in which a primitive form of novel cooperation evolved. This reduction and reversal of a positive social allelic effect is driven solely by changes in the genetic context in which the gene is expressed as new mutations are sequentially fixed during adaptive evolution, and explicitly demonstrates a significant evolutionary change in the genetic architecture of an ecologically important social trait.     

Expression profiles reveal parallel evolution of epistatic interactions involving the CRP regulon in Escherichia coli

The extent and nature of epistatic interactions between mutations are issues of fundamental importance in evolutionary biology. However, they are difficult to study and their influence on adaptation remains poorly understood. Here, we use a systems-level approach to examine epistatic interactions that arose during the evolution of Escherichia coli in a defined environment. We used expression arrays to compare the effect on global patterns of gene expression of deleting a central regulatory gene, crp. Effects were measured in two lineages that had independently evolved for 20,000 generations and in their common ancestor. We found that deleting crp had a much more dramatic effect on the expression profile of the two evolved lines than on the ancestor. Because the sequence of the crp gene was unchanged during evolution, these differences indicate epistatic interactions between crp and mutations at other loci that accumulated during evolution. Moreover, a striking degree of parallelism was observed between the two independently evolved lines; 115 genes that were not crp-dependent in the ancestor became dependent on crp in both evolved lines. An analysis of changes in crp dependence of well-characterized regulons identified a number of regulatory genes as candidates for harboring beneficial mutations that could account for these parallel expression changes. Mutations within three of these genes have previously been found and shown to contribute to fitness. Overall, these findings indicate that epistasis has been important in the adaptive evolution of these lines, and they provide new insight into the types of genetic changes through which epistasis can evolve. More generally, we demonstrate that expression profiles can be profitably used to investigate epistatic interactions.     

Epistasis among Deleterious Mutations in the HIV-1 Protease

A central goal in molecular evolution is to understand how genetic interactions between protein mutations shape protein function and fitness. While intergenic epistasis has been extensively explored in eukaryotes, bacteria, and viruses, intragenic epistatic interactions have been insufficiently studied. Here, we employ a model system in which lambda phage fitness correlates with the enzymatic activity of human immunodeficiency virus type 1 (HIV-1) protease to systematically determine the epistatic interactions between intragenic pairs of deleterious protein substitutions. We generated 114 genotypes of the HIV-1 protease, each carrying pairs of nucleotide substitution mutations whose separated and combined deleterious effects on fitness were then determined. A high proportion (39%) of pairs displayed lethality. Several pairs exhibited significant interactions for fitness, including positive and negative epistasis. Significant negative epistatic interactions predominated (15%) over positive interactions (2%). However, the average ± SD epistatic effect, ē = 0.0025 ± 0.1334, was not significantly different from zero (p = 0.8368). Notably, epistatic interactions, regardless of epistatic direction, tend to be more frequent in the context of less deleterious mutations. In the present study, the high frequencies of lethality and negative epistasis indicate that the HIV-1 protease is highly sensitive to the effects of deleterious mutations. Therefore, proteins may not be as robust to mutational change as is usually expected.     

The Impact of Macroscopic Epistasis on Long-Term Evolutionary Dynamics

Genetic interactions can strongly influence the fitness effects of individual mutations, yet the impact of these epistatic interactions on evolutionary dynamics remains poorly understood. Here we investigate the evolutionary role of epistasis over 50,000 generations in a well-studied laboratory evolution experiment in Escherichia coli. The extensive duration of this experiment provides a unique window into the effects of epistasis during long-term adaptation to a constant environment. Guided by analytical results in the weak-mutation limit, we develop a computational framework to assess the compatibility of a given epistatic model with the observed patterns of fitness gain and mutation accumulation through time. We find that a decelerating fitness trajectory alone provides little power to distinguish between competing models, including those that lack any direct epistatic interactions between mutations. However, when combined with the mutation trajectory, these observables place strong constraints on the set of possible models of epistasis, ruling out many existing explanations of the data. Instead, we find that the data are consistent with a “two-epoch” model of adaptation, in which an initial burst of diminishing-returns epistasis is followed by a steady accumulation of mutations under a constant distribution of fitness effects. Our results highlight the need for additional DNA sequencing of these populations, as well as for more sophisticated models of epistasis that are compatible with all of the experimental data.     

Plasticity and epistasis strongly affect bacterial fitness after losing multiple metabolic genes

Many bacterial lineages lack seemingly essential metabolic genes. Previous work suggested selective benefits could drive the loss of biosynthetic functions from bacterial genomes when the corresponding metabolites are sufficiently available in the environment. However, the factors that govern this “genome streamlining” remain poorly understood. Here we determine the effect of plasticity and epistasis on the fitness of Escherichia coli genotypes from whose genome biosynthetic genes for one, two, or three different amino acids have been deleted. Competitive fitness experiments between auxotrophic mutants and prototrophic wild‐type cells in one of two carbon environments revealed that plasticity and epistasis strongly affected the mutants’ fitness individually and interactively. Positive and negative epistatic interactions were prevalent, yet on average cancelled each other out. Moreover, epistasis correlated negatively with the expected effects of combined auxotrophy‐causing mutations, thus producing a pattern of diminishing returns. Moreover, computationally analyzing 1,432 eubacterial metabolic networks revealed that most pairs of auxotrophies co‐occurred significantly more often than expected by chance, suggesting epistatic interactions and/or environmental factors favored these combinations. Our results demonstrate that both the genetic background and environmental conditions determine the adaptive value of a loss‐of‐biochemical‐function mutation and that fitness gains decelerate, as more biochemical functions are lost.     

SIGN EPISTASIS LIMITS EVOLUTIONARY TRADE‐OFFS AT THE CONFLUENCE OF SINGLE‐ AND MULTI‐CARBON METABOLISM IN METHYLOBACTERIUM EXTORQUENS AM1

Adaptation of one set of traits is often accompanied by attenuation of traits important in other selective environments, leading to fitness trade‐offs. The mechanisms that either promote or prevent the emergence of trade‐offs remain largely unknown, and are difficult to discern in most systems. Here, we investigate the basis of trade‐offs that emerged during experimental evolution of Methylobacterium extorquens AM1 to distinct growth substrates. After 1500 generations of adaptation to a multi‐carbon substrate, succinate (S), many lineages had lost the ability to use one‐carbon compounds such as methanol (M), generating a mixture of M⁺ and M^? evolved phenotypes. We show that trade‐offs in M^? strains consistently arise via antagonistic pleiotropy through recurrent selection for loss‐of‐function mutations to ftfL (formate‐tetrahydrofolate ligase), which improved growth on S while simultaneously eliminating growth on M. But if loss of FtfL was beneficial, why were M trade‐offs not found in all populations? We discovered that eliminating FtfL was not universally beneficial on S, as it was neutral or even deleterious in certain evolved lineages that remained M⁺. This suggests that sign epistasis with earlier arising mutations prevented the emergence of mutations that drove trade‐offs through antagonistic pleiotropy, limiting the evolution of metabolic specialists in some populations.     

Resource Level Affects Relative Performance of the Two Motility Systems of Myxococcus xanthus

The adventurous (A) and social (S) motility systems of the microbial predator Myxococcus xanthus show differential swarming performance on distinct surface types. Under standard laboratory conditions, A-motility performs well on hard agar but poorly on soft agar, whereas the inverse pattern is shown by S-motility. These properties may allow M. xanthus to swarm effectively across a greater diversity of natural surfaces than would be possible with one motility system alone. Nonetheless, the range of ecological conditions under which dual motility enhances effective swarming across distinct surfaces and how ecological parameters affect the complementarity of A-motility and S-motility remain unclear. Here we have examined the role of nutrient concentration in determining swarming patterns driven by dual motility on distinct agar surfaces, as well as the relative contributions of A-motility and S-motility to these patterns. Swarm expansion rates of dually motile (A⁺S⁺), solely A-motile (A⁺S⁻), and solely S-motile (A⁻S⁺) strains were compared on hard and soft agar across a wide range of casitone concentrations. At low casitone concentrations (0–0.1%), swarming on soft agar driven by S-motility is very poor, and is significantly slower than swarming on hard agar driven by A-motility. This reverses at high casitone concentration (1–3.2%) such that swarming on soft agar is much faster than swarming on hard agar. This pattern greatly constrained the ability of M. xanthus to encounter patches of prey bacteria on a soft agar surface when nutrient levels between the patches were low. The swarming patterns of a strain that is unable to produce extracellular fibrils indicate that these appendages are responsible for the elevated swarming of S-motility at high resource levels. Together, these data suggest that large contributions by S-motility to predatory swarming in natural soils may be limited to soft, wet, high-nutrient conditions that may be uncommon. Several likely benefits of S-motility to the M. xanthus life cycle are discussed, including synergistic interactions with A-motility across a wide variety of conditions.     

The Evolution of Epistasis and Its Links With Genetic Robustness, Complexity and Drift in a Phenotypic Model of Adaptation

The epistatic interactions among mutations have a large effect on the evolution of populations. In this article we provide a formalism under which epistatic interactions among pairs of mutations have a distribution whose mean can be modulated. We find that the mean epistasis is correlated to the effect of mutations or genetic robustness, which suggests that such formalism is in good agreement with most in silico models of evolution where the same pattern is observed. We further show that the evolution of epistasis is highly dependant on the intensity of drift and of how complex the organisms are, and that either positive or negative epistasis could be selected for, depending on the balance between the efficiency of selection and the intensity of drift.     

Magnitude and sign epistasis among deleterious mutations in a positive-sense plant RNA virus

How epistatic interactions between mutations determine the genetic architecture of fitness is of central importance in evolution. The study of epistasis is particularly interesting for RNA viruses because of their genomic compactness, lack of genetic redundancy, and apparent low complexity. Moreover, interactions between mutations in viral genomes determine traits such as resistance to antiviral drugs, virulence and host range. In this study we generated 53 Tobacco etch potyvirus genotypes carrying pairs of single-nucleotide substitutions and measured their separated and combined deleterious fitness effects. We found that up to 38% of pairs had significant epistasis for fitness, including both positive and negative deviations from the null hypothesis of multiplicative effects. Interestingly, the sign of epistasis was correlated with viral protein-protein interactions in a model network, being predominantly positive between linked pairs of proteins and negative between unlinked ones. Furthermore, 55% of significant interactions were cases of reciprocal sign epistasis (RSE), indicating that adaptive landscapes for RNA viruses maybe highly rugged. Finally, we found that the magnitude of epistasis correlated negatively with the average effect of mutations. Overall, our results are in good agreement to those previously reported for other viruses and further consolidate the view that positive epistasis is the norm for small and compact genomes that lack genetic robustness.     

Pervasive Pairwise Intragenic Epistasis among Sequential Mutations in TEM-1 β-Lactamase

Interactions between mutations play a central role in shaping the fitness landscape, but a clear picture of intragenic epistasis has yet to emerge. To further reveal the prevalence and patterns of intragenic epistasis, we present a survey of epistatic interactions between sequential mutations in TEM-1 β-lactamase. We measured the fitness effect of ~ 12,000 pairs of consecutive amino acid substitutions and used our previous study of the fitness effects of single amino acid substitutions to calculate epistasis for over 8000 mutation pairs. Since sequential mutations are prone to physically interact, we postulated that our study would be surveying specific epistasis instead of nonspecific epistasis. We found widespread negative epistasis, especially in beta-strands, and a high frequency of negative sign epistasis among individually beneficial mutations. Negative epistasis (52%) occurred 7.6 times as frequently as positive epistasis (6.8%). Buried residues experienced more negative epistasis that surface-exposed residues. However, TEM-1 exhibited a couple of hotspots for positive epistasis, most notably L221/ R222 at which many combinations of mutations positively interacted. This study is the first to systematically examine pairwise epistasis throughout an entire protein performing its native function in its native host.     

Complex genotype interactions influence social fitness during the developmental phase of the social amoeba Dictyostelium discoideum

When individuals interact, phenotypic variation can be partitioned into direct genetic effects (DGEs) of the individuals’ own genotypes, indirect genetic effects (IGEs) of their social partners’ genotypes and epistatic interactions between the genotypes of interacting individuals (‘genotype‐by‐genotype (G×G) epistasis’). These components can all play important roles in evolutionary processes, but few empirical studies have examined their importance. The social amoeba Dictyostelium discoideum provides an ideal system to measure these effects during social interactions and development. When starved, free‐living amoebae aggregate and differentiate into a multicellular fruiting body with a dead stalk that holds aloft viable spores. By measuring interactions among a set of natural strains, we quantify DGEs, IGEs and G×G epistasis affecting spore formation. We find that DGEs explain most of the phenotypic variance (57.6%) whereas IGEs explain a smaller (13.3%) but highly significant component. Interestingly, G×G epistasis explains nearly a quarter of the variance (23.0%), highlighting the complex nature of genotype interactions. These results demonstrate the large impact that social interactions can have on development and suggest that social effects should play an important role in developmental evolution in this system.     

Initial mutations direct alternative pathways of protein evolution

Whether evolution is erratic due to random historical details, or is repeatedly directed along similar paths by certain constraints, remains unclear. Epistasis (i.e. non-additive interaction between mutations that affect fitness) is a mechanism that can contribute to both scenarios. Epistasis can constrain the type and order of selected mutations, but it can also make adaptive trajectories contingent upon the first random substitution. This effect is particularly strong under sign epistasis, when the sign of the fitness effects of a mutation depends on its genetic background. In the current study, we examine how epistatic interactions between mutations determine alternative evolutionary pathways, using in vitro evolution of the antibiotic resistance enzyme TEM-1 β-lactamase. First, we describe the diversity of adaptive pathways among replicate lines during evolution for resistance to a novel antibiotic (cefotaxime). Consistent with the prediction of epistatic constraints, most lines increased resistance by acquiring three mutations in a fixed order. However, a few lines deviated from this pattern. Next, to test whether negative interactions between alternative initial substitutions drive this divergence, alleles containing initial substitutions from the deviating lines were evolved under identical conditions. Indeed, these alternative initial substitutions consistently led to lower adaptive peaks, involving more and other substitutions than those observed in the common pathway. We found that a combination of decreased enzymatic activity and lower folding cooperativity underlies negative sign epistasis in the clash between key mutations in the common and deviating lines (Gly238Ser and Arg164Ser, respectively). Our results demonstrate that epistasis contributes to contingency in protein evolution by amplifying the selective consequences of random mutations.     

Epistasis buffers the fitness effects of rifampicin- resistance mutations in Pseudomonas aeruginosa

Epistatic interactions between resistance mutations in antibiotic-free environments potentially play a crucial role in the spread of resistance in pathogen populations by determining the fitness cost associated with resistance. We used an experimental evolution approach to test for epistatic interactions between 14 different pairs of rifampicin mutations in the pathogenic bacterium Pseudomonas aeruginosa in 42 different rifampicin-free environments. First, we show that epistasis between rifampicin-resistance mutations tends to be antagonistic: the fitness effect of having two mutations is generally smaller than that predicted from the effects of individual mutations on the wild-type. Second, we show that sign epistasis between resistance mutations is both common and strong; most notably, pairs of deleterious resistance mutations often partially or completely compensate for each others' costs, revealing a novel mechanism for compensatory adaptation. These results suggest that antagonistic epistasis between intragenic resistance mutations may be a key determinant of the cost of antibiotic resistance and compensatory adaptation in pathogen populations.     

SELECTION BIASES THE PREVALENCE AND TYPE OF EPISTASIS ALONG ADAPTIVE TRAJECTORIES

The contribution to an organism's phenotype from one genetic locus may depend upon the status of other loci. Such epistatic interactions among loci are now recognized as fundamental to shaping the process of adaptation in evolving populations. Although little is known about the structure of epistasis in most organisms, recent experiments with bacterial populations have concluded that antagonistic interactions abound and tend to deaccelerate the pace of adaptation over time. Here, we use the NK model of fitness landscapes to examine how natural selection biases the mutations that substitute during evolution based on their epistatic interactions. We find that, even when beneficial mutations are rare, these biases are strong and change substantially throughout the course of adaptation. In particular, epistasis is less prevalent than the neutral expectation early in adaptation and much more prevalent later, with a concomitant shift from predominantly antagonistic interactions early in adaptation to synergistic and sign epistasis later in adaptation. We observe the same patterns when reanalyzing data from a recent microbial evolution experiment. These results show that when the order of substitutions is not known, standard methods of analysis may suggest that epistasis retards adaptation when in fact it accelerates it.     

Mistranslation Reduces Mutation Load in Evolving Proteins through Negative Epistasis with DNA Mutations

Translational errors during protein synthesis cause phenotypic mutations that are several orders of magnitude more frequent than DNA mutations. Such phenotypic mutations may affect adaptive evolution through their interactions with DNA mutations. To study how mistranslation may affect the adaptive evolution of evolving proteins, we evolved populations of green fluorescent protein (GFP) in either high-mistranslation or low-mistranslation Escherichia coli hosts. In both hosts, we first evolved GFP under purifying selection for the ancestral phenotype green fluorescence, and then under directional selection toward the new phenotype yellow fluorescence. High-mistranslation populations evolved modestly higher yellow fluorescence during each generation of evolution than low-mistranslation populations. We demonstrate by high-throughput sequencing that elevated mistranslation reduced the accumulation of deleterious DNA mutations under both purifying and directional selection. It did so by amplifying the fitness effects of deleterious DNA mutations through negative epistasis with phenotypic mutations. In contrast, mistranslation did not affect the incidence of beneficial mutations. Our findings show that phenotypic mutations interact epistatically with DNA mutations. By reducing a population’s mutation load, mistranslation can affect an important determinant of evolvability.     

EPISTASIS BETWEEN NEW MUTATIONS AND GENETIC BACKGROUND AND A TEST OF GENETIC CANALIZATION

Abstract The importance for fitness of epistatic interactions among mutations is poorly known, yet epistasis can exert important effects on the dynamics of evolving populations. We showed previously that epistatic interactions are common between pairs of random insertion mutations in the bacterium Escherichia coli . In this paper, we examine interactions between these mutations and other mutations by transducing each of twelve insertion mutations into two genetic backgrounds, one ancestral and the other having evolved in, and adapted to, a defined laboratory environment for 10,000 generations. To assess the effect of the mutation on fitness, we allowed each mutant to compete against its unmutated counterpart in that same environment. Overall, there was a strong positive correlation between the mutational effects on the two genetic backgrounds. Nonetheless, three of the twelve mutations had significantly different effects on the two backgrounds, indicating epistasis. There was no significant tendency for the mutations to be less harmful on the derived background. Thus, there is no evidence supporting the hypothesis that the derived bacteria had adapted, in part, by becoming buffered against the harmful effects of mutations.     

Dynamic Epistasis under Varying Environmental Perturbations

Epistasis describes the phenomenon that mutations at different loci do not have independent effects with regard to certain phenotypes. Understanding the global epistatic landscape is vital for many genetic and evolutionary theories. Current knowledge for epistatic dynamics under multiple conditions is limited by the technological difficulties in experimentally screening epistatic relations among genes. We explored this issue by applying flux balance analysis to simulate epistatic landscapes under various environmental perturbations. Specifically, we looked at gene-gene epistatic interactions, where the mutations were assumed to occur in different genes. We predicted that epistasis tends to become more positive from glucose-abundant to nutrient-limiting conditions, indicating that selection might be less effective in removing deleterious mutations in the latter. We also observed a stable core of epistatic interactions in all tested conditions, as well as many epistatic interactions unique to each condition. Interestingly, genes in the stable epistatic interaction network are directly linked to most other genes whereas genes with condition-specific epistasis form a scale-free network. Furthermore, genes with stable epistasis tend to have similar evolutionary rates, whereas this co-evolving relationship does not hold for genes with condition-specific epistasis. Our findings provide a novel genome-wide picture about epistatic dynamics under environmental perturbations.