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.     

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.     

Environment Determines Epistatic Patterns for a ssDNA Virus

Despite the accumulation of substantial quantities of information about epistatic interactions among both deleterious and beneficial mutations in a wide array of experimental systems, neither consistent patterns nor causal explanations for these interactions have yet emerged. Furthermore, the effects of mutations depend on the environment in which they are characterized, implying that the environment may also influence epistatic interactions. Recent work with beneficial mutations for the single-stranded DNA bacteriophage ID11 demonstrated that interactions between pairs of mutations could be understood by means of a simple model that assumes that mutations have additive phenotypic effects and that epistasis arises through a nonlinear phenotype–fitness map with a single intermediate optimum. To determine whether such a model could also explain changes in epistatic patterns associated with changes in environment, we measured epistatic interactions for these same mutations under conditions for which we expected to find the wild-type ID11 at different distances from its phenotypic optimum by assaying fitnesses at three different temperatures: 33°, 37°, and 41°. Epistasis was present and negative under all conditions, but became more pronounced as temperature increased. We found that the additive-phenotypes model explained these patterns as changes in the parameters of the phenotype–fitness map, but that a model that additionally allows the phenotypes to vary across temperatures performed significantly better. Our results show that ostensibly complex patterns of fitness effects and epistasis across environments can be explained by assuming a simple structure for the genotype–phenotype relationship.     

Strong Epistatic Selection on the RNA Secondary Structure of HIV

A key question in evolutionary genomics is how populations navigate the adaptive landscape in the presence of epistasis, or interactions among loci. This problem can be directly addressed by studying the evolution of RNA secondary structures, for which there is constraint to maintain pairing between Watson-Crick (WC) sites. Replacement of a nucleotide at one site of a WC pair reduces fitness by disrupting binding, which can be restored via a compensatory replacement at the interacting site. Here, I present the first genome-scale analysis of epistasis on the RNA secondary structure of human immunodeficiency virus type 1 (HIV-1). Comparison of polymorphism frequencies at ancestrally conserved sites reveals that selection against replacements is ∼2.7 times stronger at WC than at non-WC sites, such that nearly 50% of constraint can be attributed to epistasis. However, almost all epistatic constraint is due to selection against conversions of WC pairs to unpaired (UP) nucleotides, whereas conversions to GU wobbles are only slightly deleterious. This disparity is also evident in pairs with second-site compensatory replacements; conversions from UP nucleotides to WC pairs increase median fitness by ∼4.2%, whereas conversions from GU wobbles to WC pairs only increase median fitness by ∼0.3%. Moreover, second-site replacements that convert UP nucleotides to GU wobbles also increase median fitness by ∼4%, indicating that such replacements are nearly as compensatory as those that restore WC pairing. Thus, WC peaks of the HIV-1 epistatic adaptive landscape are connected by high GU ridges, enabling the viral population to rapidly explore distant peaks without traversing deep UP valleys.     

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.     

HIV-1 protease catalytic efficiency effects caused by random single amino acid substitutions

Protein evolution has occurred by successive fixation of individual mutations. The probability of fixation depends on the fitness of the mutation, and the arising variant can be deleterious, neutral, or beneficial. Despite its relevance, only few studies have estimated the distribution of fitness effects caused by random single mutations on protein function. The human immunodeficiency virus type 1 (HIV-1) protease was chosen as a model protein to quantify protein's tolerability to random single mutations. After determining the enzymatic activity of 107 single random mutants, we found that 86% of single mutations were deleterious for the enzyme catalytic efficiency and 54% lethal. Only 2% of the mutations significantly increased the catalytic efficiency of the enzyme. These data demonstrate the vulnerability of HIV-1 protease to single random mutations. When a second random mutagenesis library was constructed from an HIV-1 protease carrying a highly deleterious single mutation (D30N), a higher proportion of mutations with neutral or beneficial effect were found, 26% and 9%, respectively. Importantly, antagonist epistasis was observed between deleterious mutations. In particular, the mutation N88D, lethal for the wild-type protease, restored the wild-type catalytic efficiency when combined with the highly deleterious mutation D30N. The low tolerability to single random substitutions shown here for the wild-type HIV-1 protease contrasts with its in vivo ability to generate an adaptive variation. Thus, the antagonist epistasis between deleterious or lethal mutations may be responsible for increasing the protein mutational robustness and evolvability.     

Fitness ranking of individual mutants drives patterns of epistatic interactions in HIV-1

Fitness interactions between mutations, referred to as epistasis, can strongly impact evolution. For RNA viruses and retroviruses with their high mutation rates, epistasis may be particularly important to overcome fitness losses due to the accumulation of deleterious mutations and thus could influence the frequency of mutants in a viral population. As human immunodeficiency virus type 1 (HIV-1) resistance to azidothymidine (AZT) requires selection of sequential mutations, it is a good system to study the impact of epistasis. Here we present a thorough analysis of a classical AZT-resistance pathway (the 41-215 cluster) of HIV-1 variants by fitness measurements in single round infection assays covering physiological drug concentrations ex vivo. The sign and value of epistasis varied and did not predict the epistatic effect on the mutant frequency. This complex behavior is explained by the fitness ranking of the variants that strongly depends on environmental factors, i.e., the presence and absence of drugs and the host cells used. Although some interactions compensate fitness losses, the observed small effect on the relative mutant frequencies suggests that epistasis might be inefficient as a buffering mechanism for fitness losses in vivo. While the use of epistasis-based hypotheses to make general assumptions on the evolutionary dynamics of viral populations is appealing, our data caution their interpretation without further knowledge on the characteristics of the viral mutant spectrum under different environmental conditions.     

Epistasis as a Determinant of the HIV-1 Protease's Robustness to Mutation

The robustness of phenotypes to mutation is critical to protein evolution; robustness may be an adaptive trait if it promotes evolution. We hypothesised that native proteins subjected to natural selection in vivo should be more robust than proteins generated in vitro in the absence of natural selection. We compared the mutational robustness of two human immunodeficiency virus type 1 (HIV-1) proteases with comparable catalytic efficiencies, one isolated from an infected individual and the second generated in vitro via random mutagenesis. Single mutations in the protease (82 and 60 in the wild-type and mutant backgrounds, respectively) were randomly generated in vitro and the catalytic efficiency of each mutant was determined. No differences were observed between these two protease variants when lethal, neutral, and deleterious mutations were compared (P = 0.8025, chi-squared test). Similarly, average catalytic efficiency (−72.6% and −64.5%, respectively) did not significantly differ between protease mutant libraries (P = 0.3414, Mann Whitney test). Overall, the two parental proteins displayed similar mutational robustness. Importantly, strong and widespread epistatic interactions were observed when the effect of the same mutation was compared in both proteases, suggesting that epistasis can be a key determinant of the robustness displayed by the in vitro generated protease.     

Long-term adaptation of epistatic genetic networks

Gene networks are likely to govern most traits in nature. Mutations at these genes often show functional epistatic interactions that lead to complex genetic architectures and variable fitness effects in different genetic backgrounds. Understanding how epistatic genetic systems evolve in nature remains one of the great challenges in evolutionary biology. Here we combine an analytical framework with individual-based simulations to generate novel predictions about long-term adaptation of epistatic networks. We find that relative to traits governed by independently evolving genes, adaptation with epistatic gene networks is often characterized by longer waiting times to selective sweeps, lower standing genetic variation, and larger fitness effects of adaptive mutations. This may cause epistatic networks to either adapt more slowly or more quickly relative to a nonepistatic system. Interestingly, epistatic networks may adapt faster even when epistatic effects of mutations are on average deleterious. Further, we study the evolution of epistatic properties of adaptive mutations in gene networks. Our results show that adaptive mutations with small fitness effects typically evolve positive synergistic interactions, whereas adaptive mutations with large fitness effects evolve positive synergistic and negative antagonistic interactions at approximately equal frequencies. These results provide testable predictions for adaptation of traits governed by epistatic networks and the evolution of epistasis within networks.     

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.     

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.     

Interaction between directional epistasis and average mutational effects.

We investigate the relationship between the average fitness decay due to single mutations and the strength of epistatic interactions in genetic sequences. We observe that epistatic interactions between mutations are correlated to the average fitness decay, both in RNA secondary structure prediction as well as in digital organisms replicating in silico. This correlation implies that, during adaptation, epistasis and average mutational effect cannot be optimized independently. In experiments with RNA sequences evolving on a neutral network, the selective pressure to decrease the mutational load then leads to a reduction in the amount of sequences with strong antagonistic interactions between deleterious mutations in the population.     

Recombination Accelerates Adaptation on a Large-Scale Empirical Fitness Landscape in HIV-1

Recombination has the potential to facilitate adaptation. In spite of the substantial body of theory on the impact of recombination on the evolutionary dynamics of adapting populations, empirical evidence to test these theories is still scarce. We examined the effect of recombination on adaptation on a large-scale empirical fitness landscape in HIV-1 based on in vitro fitness measurements. Our results indicate that recombination substantially increases the rate of adaptation under a wide range of parameter values for population size, mutation rate and recombination rate. The accelerating effect of recombination is stronger for intermediate mutation rates but increases in a monotonic way with the recombination rates and population sizes that we examined. We also found that both fitness effects of individual mutations and epistatic fitness interactions cause recombination to accelerate adaptation. The estimated epistasis in the adapting populations is significantly negative. Our results highlight the importance of recombination in the evolution of HIV-I.     

Epistatic interactions promote persistence of NS3-Q80K in HCV infection by compensating for protein folding instability

The Q80K polymorphism in the NS3-4A protease of the hepatitis C virus is associated with treatment failure of direct-acting antiviral agents. This polymorphism is highly prevalent in genotype 1a infections and stably transmitted between hosts. Here, we investigated the underlying molecular mechanisms of evolutionarily conserved coevolving amino acids in NS3-Q80K and revealed potential implications of epistatic interactions in immune escape and variants persistence. Using purified protein, we characterized the impact of epistatic amino acid substitutions on the physicochemical properties and peptide cleavage kinetics of the NS3-Q80K protease. We found that Q80K destabilized the protease protein fold (p < 0.0001). Although NS3-Q80K showed reduced peptide substrate turnover (p < 0.0002), replicative fitness in an H77S.3 cell culture model of infection was not significantly inferior to the WT virus. Epistatic substitutions at residues 91 and 174 in NS3-Q80K stabilized the protein fold (p < 0.0001) and leveraged the WT protease stability. However, changes in protease stability inversely correlated with enzymatic activity. In infectious cell culture, these secondary substitutions were not associated with a gain of replicative fitness in NS3-Q80K variants. Using molecular dynamics, we observed that the total number of residue contacts in NS3-Q80K mutants correlated with protein folding stability. Changes in the number of contacts reflected the compensatory effect on protein folding instability by epistatic substitutions. In summary, epistatic substitutions in NS3-Q80K contribute to viral fitness by mechanisms not directly related to RNA replication. By compensating for protein-folding instability, epistatic interactions likely protect NS3-Q80K variants from immune cell recognition.     

Pleiotropic Effects of a Mitochondrial–Nuclear Incompatibility Depend upon the Accelerating Effect of Temperature in Drosophila

Interactions between mitochondrial and nuclear gene products that underlie eukaryotic energy metabolism can cause the fitness effects of mutations in one genome to be conditional on variation in the other genome. In ectotherms, the effects of these interactions are likely to depend upon the thermal environment, because increasing temperature accelerates molecular rates. We find that temperature strongly modifies the pleiotropic phenotypic effects of an incompatible interaction between a Drosophila melanogaster polymorphism in the nuclear-encoded, mitochondrial tyrosyl-transfer (t)RNA synthetase and a D. simulans polymorphism in the mitochondrially encoded tRNA^Tyr. The incompatible mitochondrial–nuclear genotype extends development time, decreases larval survivorship, and reduces pupation height, indicative of decreased energetic performance. These deleterious effects are ameliorated when larvae develop at 16° and exacerbated at warmer temperatures, leading to complete sterility in both sexes at 28°. The incompatible genotype has a normal metabolic rate at 16° but a significantly elevated rate at 25°, consistent with the hypothesis that inefficient energy metabolism extends development in this genotype at warmer temperatures. Furthermore, the incompatibility decreases metabolic plasticity of larvae developed at 16°, indicating that cooler development temperatures do not completely mitigate the deleterious effects of this genetic interaction. Our results suggest that the epistatic fitness effects of metabolic mutations may generally be conditional on the thermal environment. The expression of epistatic interactions in some environments, but not others, weakens the efficacy of selection in removing deleterious epistatic variants from populations and may promote the accumulation of incompatibilities whose fitness effects will depend upon the environment in which hybrids occur.     

The role of deleterious mutations in allopatric speciation

Allopatric speciation is often assumed to occur as a consequence of adaptive divergence between two isolated populations. However, there are some scenarios in which reproductive isolation can be favored due to accumulated unconditionally deleterious mutations. If deleterious mutations have synergistic epistatic effects, it is shown here that the average fitness of recombinants between two parental lines with a given number of fixed mutations is lower than that of the parents in both the F1 and F2 generations. If individual mutations are only slightly deleterious, then they will tend to fixation at a high enough rate to cause lower hybrid fitness. If the fitness effects of mutation give rise to antagonistic epistasis, the hybrids tend to have a higher average fitness than the parental lines, suggesting a possible scenario for the origin of hybrid vigor. The other model of deleterious mutations investigated is the accumulation of knockout mutants in a duplicated gene family. While neutral in the parental lines, upon contact the F1 and later generations have a significant probability of carrying double knockouts. Under this scenario, selection may also favor reproductive isolation between the two lines. Even when the selection coefficients generated are too low to drive speciation, epistatic interactions between deleterious mutations offer a possible explanation for both outbreeding depression and hybrid vigor.     

Epistasis and the adaptability of an RNA virus

We have explored the patterns of fitness recovery in the vesicular stomatitis RNA virus. We show that, in our experimental setting, reversions to the wild-type genotype were rare and fitness recovery was at least partially driven by compensatory mutations. We compared compensatory adaptation for genotypes carrying (1) mutations with varying deleterious fitness effects, (2) one or two deleterious mutations, and (3) pairs of mutations showing differences in the strength and sign of epistasis. In all cases, we found that the rate of fitness recovery and the proportion of reversions were positively affected by population size. Additionally, we observed that mutations with large fitness effect were always compensated faster than mutations with small fitness effect. Similarly, compensatory evolution was faster for genotypes carrying a single deleterious mutation than for those carrying pairs of mutations. Finally, for genotypes carrying two deleterious mutations, we found evidence of a negative correlation between the epistastic effect and the rate of compensatory evolution.     

Prevalence of epistasis in the evolution of influenza A surface proteins

The surface proteins of human influenza A viruses experience positive selection to escape both human immunity and, more recently, antiviral drug treatments. In bacteria and viruses, immune-escape and drug-resistant phenotypes often appear through a combination of several mutations that have epistatic effects on pathogen fitness. However, the extent and structure of epistasis in influenza viral proteins have not been systematically investigated. Here, we develop a novel statistical method to detect positive epistasis between pairs of sites in a protein, based on the observed temporal patterns of sequence evolution. The method rests on the simple idea that a substitution at one site should rapidly follow a substitution at another site if the sites are positively epistatic. We apply this method to the surface proteins hemagglutinin and neuraminidase of influenza A virus subtypes H3N2 and H1N1. Compared to a non-epistatic null distribution, we detect substantial amounts of epistasis and determine the identities of putatively epistatic pairs of sites. In particular, using sequence data alone, our method identifies epistatic interactions between specific sites in neuraminidase that have recently been demonstrated, in vitro, to confer resistance to the drug oseltamivir; these epistatic interactions are responsible for widespread drug resistance among H1N1 viruses circulating today. This experimental validation demonstrates the predictive power of our method to identify epistatic sites of importance for viral adaptation and public health. We conclude that epistasis plays a large role in shaping the molecular evolution of influenza viruses. In particular, sites with , which would normally not be identified as positively selected, can facilitate viral adaptation through epistatic interactions with their partner sites. The knowledge of specific interactions among sites in influenza proteins may help us to predict the course of antigenic evolution and, consequently, to select more appropriate vaccines and drugs.     

Resistance mechanism revealed by crystal structures of unliganded nelfinavir-resistant HIV-1 protease non-active site mutants N88D and N88S

Nelfinavir is an inhibitor of HIV-1 protease, and is used for treatment of patients suffering from HIV/AIDS. However, treatment results in drug resistant mutations in HIV-1 protease. N88D and N88S are two such mutations which occur in the non-active site region of the enzyme. We have determined crystal structures of unliganded N88D and N88S mutants of HIV-1 protease to resolution of 1.65 Å and 1.8 Å, respectively. These structures refined against synchrotron data lead to R-factors of 0.1859 and 0.1780, respectively. While structural effects of N88D are very subtle, the mutation N88S has caused a significant conformational change in D30, an active site residue crucial for substrate and inhibitor binding.     

Asymptotic distribution for epistatic tests in case-control studies

We propose a statistical model for dissecting a multilocus genotypic value into its main (additive and dominant) effects and epistatic effects between different loci in a case-control association study. The model can discern four different kinds of epistasis, additive × additive, additive × dominant, dominant × additive, and dominant × dominant interactions. To test each kind of epistasis, a χ² test statistic was computed for a two by two contingency table derived from combined genotypes in both case and control groups. We derived an analytical approach for estimating the asymptotic distribution of the χ² test statistic for epistatic tests under the null hypothesis, with the result being consistent with that from Monte Carlo simulations. The new model was used to analyze a case-control data set for candidate gene studies of stroke, leading to the identification of several significant interactions between causal SNPs on this disease.