Strong Selection Significantly Increases Epistatic Interactions in the Long-Term Evolution of a Protein

Epistatic interactions between residues determine a protein's adaptability and shape its evolutionary trajectory. When a protein experiences a changed environment, it is under strong selection to find a peak in the new fitness landscape. It has been shown that strong selection increases epistatic interactions as well as the ruggedness of the fitness landscape, but little is known about how the epistatic interactions change under selection in the long-term evolution of a protein. Here we analyze the evolution of epistasis in the protease of the human immunodeficiency virus type 1 (HIV-1) using protease sequences collected for almost a decade from both treated and untreated patients, to understand how epistasis changes and how those changes impact the long-term evolvability of a protein. We use an information-theoretic proxy for epistasis that quantifies the co-variation between sites, and show that positive information is a necessary (but not sufficient) condition that detects epistasis in most cases. We analyze the "fossils" of the evolutionary trajectories of the protein contained in the sequence data, and show that epistasis continues to enrich under strong selection, but not for proteins whose environment is unchanged. The increase in epistasis compensates for the information loss due to sequence variability brought about by treatment, and facilitates adaptation in the increasingly rugged fitness landscape of treatment. While epistasis is thought to enhance evolvability via valley-crossing early-on in adaptation, it can hinder adaptation later when the landscape has turned rugged. However, we find no evidence that the HIV-1 protease has reached its potential for evolution after 9 years of adapting to a drug environment that itself is constantly changing.Comment: 25 pages, 9 figures, plus Supplementary Material including Supplementary Text S1-S7, Supplementary Tables S1-S2, and Supplementary Figures S1-2. Version that appears in PLoS Genetic

Compensatory mutations cause excess of antagonistic epistasis in RNA secondary structure folding

BACKGROUND: The rate at which fitness declines as an organism's genome accumulates random mutations is an important variable in several evolutionary theories. At an intuitive level, it might seem natural that random mutations should tend to interact synergistically, such that the rate of mean fitness decline accelerates as the number of random mutations is increased. However, in a number of recent studies, a prevalence of antagonistic epistasis (the tendency of multiple mutations to have a mitigating rather than reinforcing effect) has been observed. RESULTS: We studied in silico the net amount and form of epistatic interactions in RNA secondary structure folding by measuring the fraction of neutral mutants as a function of mutational distance d. We found a clear prevalence of antagonistic epistasis in RNA secondary structure folding. By relating the fraction of neutral mutants at distance d to the average neutrality at distance d, we showed that this prevalence derives from the existence of many compensatory mutations at larger mutational distances. CONCLUSIONS: Our findings imply that the average direction of epistasis in simple fitness landscapes is directly related to the density with which fitness peaks are distributed in these landscapes

Environmental change exposes beneficial epistatic interactions in a catalytic RNA

Natural selection drives populations of individuals towards local peaks in a fitness landscape. These peaks are created by the interactions between individual mutations. Fitness landscapes may change as an environment changes. In a previous contribution, we discovered a variant of the Azoarcus group I ribozyme that represents a local peak in the RNA fitness landscape. The genotype at this peak is distinguished from the wild-type by four point mutations. We here report ribozyme fitness data derived from constructing all possible combinations of these point mutations. We find that these mutations interact epistatically. Importantly, we show that these epistatic interactions change qualitatively in the three different environments that we studied. We find examples where the relative fitness of a ribozyme can change from neutral or negative in one environment, to positive in another. We also show that the fitness effect of a specific GC–AU base pair switch is dependent on both the environment and the genetic context. Moreover, the mutations that we study improve activity at the cost of decreased structural stability. Environmental change is ubiquitous in nature. Our results suggest that such change can facilitate adaptive evolution by exposing new peaks of a fitness landscape. They highlight a prominent role for genotype–environment interactions in doing so