Perspective: Sign epistasis and genetic constraint on evolutionary trajectories

Epistasis for fitness means that the selective effect of a mutation is conditional on the genetic background in which it appears. Although epistasis is widely observed in nature, our understanding of its consequences for evolution by natural selection remains incomplete. In particular, much attention focuses only on its influence on the instantaneous rate of changes in frequency of selected alleles via epistatic contribution to the additive genetic variance for fitness. Thus, in this framework epistasis only has evolutionary importance if the interacting loci are simultaneously segregating in the population. However, the selective accessibility of mutational trajectories to high fitness genotypes may depend on the genetic background in which novel mutations appear, and this effect is independent of population polymorphism at other loci. Here we explore this second influence of epistasis on evolution by natural selection. We show that it is the consequence of a particular form of epistasis, which we designate sign epistasis. Sign epistasis means that the sign of the fitness effect of a mutation is under epistatic control; thus, such a mutation is beneficial on some genetic backgrounds and deleterious on others. Recent experimental innovations in microbial systems now permit assessment of the fitness effects of individual mutations on multiple genetic backgrounds. We review this literature and identify many examples of sign epistasis, and we suggest that the implications of these results may generalize to other organisms. These theoretical and empirical considerations imply that strong genetic constraint on the selective accessibility of trajectories to high fitness genotypes may exist and suggest specific areas of investigation for future research.     

The distribution of fitness effects among beneficial mutations in Fisher's geometric model of adaptation

Recent models of adaptation at the DNA sequence level assume that the fitness effects of new mutations show certain statistical properties. In particular, these models assume that the distribution of fitness effects among new mutations is in the domain of attraction of the so-called Gumbel-type extreme value distribution. This assumption has not, however, been justified on any biological or theoretical grounds. In this note, I study random mutation in one of the simplest models of mutation and adaptation-Fisher's geometric model. I show that random mutation in this model yields a distribution of mutational effects that belongs to the Gumbel type. I also show that the distribution of fitness effects among rare beneficial mutations in Fisher's model is asymptotically exponential. I confirm these analytic findings with exact computer simulations. These results provide some support for the use of Gumbel-type extreme value theory in studies of adaptation and point to a surprising connection between recent phenotypic- and sequence-based models of adaptation: in both, the distribution of fitness effects among rare beneficial mutations is approximately exponential.     

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

The fitness effect of a mutation can depend on both its genetic background, known as epistasis, and the prevailing external environment. Many examples of these dependencies are known, but few studies consider both aspects in combination, especially as they affect mutations that have been selected together. We examine interactions between five coevolved mutations in eight diverse environments. We find that mutations are, on average, beneficial across environments, but that there is high variation in their fitness effects, including many examples of mutations conferring a cost in some, but not other, genetic background‐environment combinations. Indeed, even when global interaction trends are accounted for, specific local mutation interactions are common and differed across environments. One consequence of this dependence is that the range of trade‐offs in genotype fitness across selected and alternative environments are contingent on the particular evolutionary path followed over the mutation landscape. Finally, although specific interactions were common, there was a consistent pattern of diminishing returns epistasis whereby mutation effects were less beneficial when added to genotypes of higher fitness. Our results underline that specific mutation effects are highly dependent on the combination of genetic and external environments, and support a general relationship between a genotype's current fitness and its potential to increase in fitness.     

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life‐history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation‐accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.     

A resolution of the mutation load paradox in humans

Current information on the rate of mutation and the fraction of sites in the genome that are subject to selection suggests that each human has received, on average, at least two new harmful mutations from its parents. These mutations were subsequently removed by natural selection through reduced survival or fertility. It has been argued that the mutation load, the proportional reduction in population mean fitness relative to the fitness of an idealized mutation-free individual, allows a theoretical prediction of the proportion of individuals in the population that fail to reproduce as a consequence of these harmful mutations. Application of this theory to humans implies that at least 88% of individuals should fail to reproduce and that each female would need to have more than 16 offspring to maintain population size. This prediction is clearly at odds with the low reproductive excess of human populations. Here, we derive expressions for the fraction of individuals that fail to reproduce as a consequence of recurrent deleterious mutation () for a model in which selection occurs via differences in relative fitness, such as would occur through competition between individuals. We show that is much smaller than the value predicted by comparing fitness to that of a mutation-free genotype. Under the relative fitness model, we show that depends jointly on U and the selective effects of new deleterious mutations and that a species could tolerate 10's or even 100's of new deleterious mutations per genome each generation.     

Complexity, pleiotropy, and the fitness effect of mutations

One of the assumptions underlying many theoretical predictions in evolutionary biology concerns the distribution of the fitness effect of mutations. Approximations to this distribution have been derived using various theoretical approaches, of which Fisher's geometrical model is among the most popular ones. Two key concepts in this model are complexity and pleiotropy. Recent studies have proposed different methods for estimating how complexity varies across species, but their results have been contradictory. Here, we show that contradictory results are to be expected when the assumption of universal pleiotropy is violated. We develop a model in which the two key parameters are the total number of traits and the mean number of traits affected by a single mutation. We derive approximations for the distribution of the fitness effect of mutations when populations are either well-adapted or away from the optimum. We also consider drift load in a well-adapted population and show that it is independent of the distribution of the fitness effect of mutations. We show that mutation accumulation experiments can only measure the effect of the mean number of traits affected by mutations, whereas drift load only provides information about the total number of traits. We discuss the plausibility of the model.     

Consequences of life history for inbreeding depression and mating system evolution in plants

Many plants are perennial, but most studies of inbreeding depression and mating system evolution focus on annuals. This paper extends a population genetic model of inbreeding depression due to recessive deleterious mutations to perennials. The model incorporates life history and mating system variation, and multiplicative selection across many genetic loci. In the absence of substantial mitotic mutation, perennials have higher mean fitness and lower, or even negative, inbreeding depression than annuals with the same mating system. As in annuals, self fertilization exposes deleterious recessive mutations to selection, increasing mean fitness and decreasing inbreeding depression. Including mitotic mutation decreases mean fitness while increasing inbreeding depression. Perenniality introduces a kind of selective sieve, such that strongly recessive mutations contribute disproportionately to mean fitness and inbreeding depression. In the presence of high mitotic mutation, this selective sieve may provide a mechanistic basis for high inbreeding depression observed in some long lived perennials. Without substantial mitotic mutation, it is difficult to reconcile genetically based models of inbreeding depression with the empirical generalization that perennials outcross while related annuals self fertilize.     

Imperfect Genes, Fisherian Mutation and the Evolution of Sex

In this paper we present a mathematical model of mutation and selection that allows for the coexistence of multiple alleles at a locus with very small selective differences between alleles. The model also allows for the determination of fitness by multiple loci. Models of this sort are biologically plausible. However, some previous attempts to construct similar models have assumed that all mutations produce a decrease in fitness, and this has led to a tendency for the average fitness of population members to decline when population numbers are finite. In our model we incorporate some of the ideas of R. A. FISHER, so that both deleterious and beneficial mutations are possible. As a result, average fitness tends to approach a stationary distribution. We have used computer simulation methods to apply the Fisherian mutation model to the problem of the evolution of sex and recombination. The results suggest that sex and recombination can provide very large benefits in terms of average fitness. The results also suggest that obligately sexual species will win ecological competitions with species that produce a substantial fraction of their offspring asexually, so long as the number of sites under selection within the genomes of the competing species is not too small and the population sizes are not too large. Our model focuses on fertility selection in an hermaphroditic plant. However, the results are likely to generalize to a wide variety of other situations as well.     

Estimation of Deleterious-Mutation Parameters in Natural Populations

The rate and average effects of spontaneous deleterious mutations are important determinants of the evolution of breeding systems and of the vulnerability of small populations to extinction. Nevertheless, few attempts have been made to estimate the properties of such mutations, and those studies that have been performed have been extremely labor intensive, relying on long-term, laboratory mutation-accumulation experiments. We present an alternative to the latter approach. For populations in which the genetic variance for fitness is a consequence of selection-mutation balance, the mean fitness and genetic variance of fitness in outbred and inbred generations can be expressed as simple functions of the genomic mutation rate, average homozygous effect and average dominance coefficient of new mutations. Using empirical estimates for the mean and genetic variance of fitness, these expressions can then be solved to obtain joint estimates of the deleterious-mutation parameters. We employ computer simulations to evaluate the degree of bias of the estimators and present some general recommendations on the application of the technique. Our procedures provide some hope for obtaining estimates of the properties of deleterious mutations from a wide phylogenetic range of species as well as a mechanism for testing the validity of alternative models for the maintenance of genetic variance for fitness.     

Models of Frequency-Dependent Selection with Mutation from Parental Alleles

Frequency-dependent selection (FDS) remains a common heuristic explanation for the maintenance of genetic variation in natural populations. The pairwise-interaction model (PIM) is a well-studied general model of frequency-dependent selection, which assumes that a genotype’s fitness is a function of within-population intergenotypic interactions. Previous theoretical work indicated that this type of model is able to sustain large numbers of alleles at a single locus when it incorporates recurrent mutation. These studies, however, have ignored the impact of the distribution of fitness effects of new mutations on the dynamics and end results of polymorphism construction. We suggest that a natural way to model mutation would be to assume mutant fitness is related to the fitness of the parental allele, i.e., the existing allele from which the mutant arose. Here we examine the numbers and distributions of fitnesses and alleles produced by construction under the PIM with mutation from parental alleles and the impacts on such measures due to different methods of generating mutant fitnesses. We find that, in comparison with previous results, generating mutants from existing alleles lowers the average number of alleles likely to be observed in a system subject to FDS, but produces polymorphisms that are highly stable and have realistic allele-frequency distributions.     

Computational alanine scanning with linear scaling semiempirical quantum mechanical methods

Alanine scanning is a powerful experimental tool for understanding the key interactions in protein–protein interfaces. Linear scaling semiempirical quantum mechanical calculations are now sufficiently fast and robust to allow meaningful calculations on large systems such as proteins, RNA and DNA. In particular, they have proven useful in understanding protein–ligand interactions. Here we ask the question: can these linear scaling quantum mechanical methods developed for protein–ligand scoring be useful for computational alanine scanning? To answer this question, we assembled 15 protein–protein complexes with available crystal structures and sufficient alanine scanning data. In all, the data set contains ΔΔGs for 400 single point alanine mutations of these 15 complexes. We show that with only one adjusted parameter the quantum mechanics‐based methods outperform both buried accessible surface area and a potential of mean force and compare favorably to a variety of published empirical methods. Finally, we closely examined the outliers in the data set and discuss some of the challenges that arise from this examination. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Assessing in vivo mutation frequencies and creating a high-resolution genome-wide map of fitness costs of Hepatitis C virus

Like many viruses, Hepatitis C Virus (HCV) has a high mutation rate, which helps the virus adapt quickly, but mutations come with fitness costs. Fitness costs can be studied by different approaches, such as experimental or frequency-based approaches. The frequency-based approach is particularly useful to estimate in vivo fitness costs, but this approach works best with deep sequencing data from many hosts are. In this study, we applied the frequency-based approach to a large dataset of 195 patients and estimated the fitness costs of mutations at 7957 sites along the HCV genome. We used beta regression and random forest models to better understand how different factors influenced fitness costs. Our results revealed that costs of nonsynonymous mutations were three times higher than those of synonymous mutations, and mutations at nucleotides A or T had higher costs than those at C or G. Genome location had a modest effect, with lower costs for mutations in HVR1 and higher costs for mutations in Core and NS5B. Resistance mutations were, on average, costlier than other mutations. Our results show that in vivo fitness costs of mutations can be site and virus specific, reinforcing the utility of constructing in vivo fitness cost maps of viral genomes.     

Fixation probabilities depend on life history: fecundity, generation time and survival in a burst-death model

The burst-death model has been developed to describe the life history of organisms with variable generation times and a burst of a fixed number of offspring. The model also includes an optional constant clearance rate, such as washout from a chemostat, and the possibility of sustained periods of population growth followed by severe bottlenecks, as in serial passaging. In this model, a beneficial mutation can either increase the burst rate or the burst size, or reduce the clearance rate, thus increasing survival. In this article we examine the effects of these three possible mechanisms on both the Malthusian fitness and the fixation probability of the lineage. We find that equivalent relative increases in the burst rate or burst size confer equivalent increases in the Malthusian fitness of a lineage, whereas increasing survival typically has a more moderate effect on Malthusian fitness. In contrast, for beneficial mutations that confer the same increase in fitness, mutations that increase survival are the most likely to fix, followed by mutations that increase the burst rate. Mutations that increase the burst size are the least likely to fix. These results imply that mutant lineages with the highest Malthusian fitness are not, in many cases, the most likely to escape extinction.     

Frequent beneficial mutations during single-colony serial transfer of Streptococcus pneumoniae

The appearance of new mutations within a population provides the raw material for evolution. The consistent decline in fitness observed in classical mutation accumulation studies has provided support for the long-held view that deleterious mutations are more common than beneficial mutations. Here we present results of a study using a mutation accumulation design with the bacterium Streptococcus pneumoniae in which the fitness of the derived populations increased. This rise in fitness was associated specifically with adaptation to survival during brief stationary phase periods between single-colony population bottlenecks. To understand better the population dynamics behind this unanticipated adaptation, we developed a maximum likelihood model describing the processes of mutation and stationary-phase selection in the context of frequent population bottlenecks. Using this model, we estimate that the rate of beneficial mutations may be as high as 4.8×10(-4) events per genome for each time interval corresponding to the pneumococcal generation time. This rate is several orders of magnitude higher than earlier estimates of beneficial mutation rates in bacteria but supports recent results obtained through the propagation of small populations of Escherichia coli. Our findings indicate that beneficial mutations may be relatively frequent in bacteria and suggest that in S. pneumoniae, which develops natural competence for transformation, a steady supply of such mutations may be available for sampling by recombination.     

Statistical Mechanics and Thermodynamics of Viral Evolution

This paper uses methods drawn from physics to study the life cycle of viruses. The paper analyzes a model of viral infection and evolution using the "grand canonical ensemble" and formalisms from statistical mechanics and thermodynamics. Using this approach we enumerate all possible genetic states of a model virus and host as a function of two independent pressures–immune response and system temperature. We prove the system has a real thermodynamic temperature, and discover a new phase transition between a positive temperature regime of normal replication and a negative temperature “disordered” phase of the virus. We distinguish this from previous observations of a phase transition that arises as a function of mutation rate. From an evolutionary biology point of view, at steady state the viruses naturally evolve to distinct quasispecies. This paper also reveals a universal relationship that relates the order parameter (as a measure of mutational robustness) to evolvability in agreement with recent experimental and theoretical work. Given that real viruses have finite length RNA segments that encode proteins which determine virus fitness, the approach used here could be refined to apply to real biological systems, perhaps providing insight into immune escape, the emergence of novel pathogens and other results of viral evolution.     

The pace of evolution across fitness valleys

How fast does a population evolve from one fitness peak to another? We study the dynamics of evolving, asexually reproducing populations in which a certain number of mutations jointly confer a fitness advantage. We consider the time until a population has evolved from one fitness peak to another one with a higher fitness. The order of mutations can either be fixed or random. If the order of mutations is fixed, then the population follows a metaphorical ridge, a single path. If the order of mutations is arbitrary, then there are many ways to evolve to the higher fitness state. We address the time required for fixation in such scenarios and study how it is affected by the order of mutations, the population size, the fitness values and the mutation rate.     

RATES OF FITNESS DECLINE AND REBOUND SUGGEST PERVASIVE EPISTASIS

Unraveling the factors that determine the rate of adaptation is a major question in evolutionary biology. One key parameter is the effect of a new mutation on fitness, which invariably depends on the environment and genetic background. The fate of a mutation also depends on population size, which determines the amount of drift it will experience. Here, we manipulate both population size and genotype composition and follow adaptation of 23 distinct Escherichia coli genotypes. These have previously accumulated mutations under intense genetic drift and encompass a substantial fitness variation. A simple rule is uncovered: the net fitness change is negatively correlated with the fitness of the genotype in which new mutations appear—a signature of epistasis. We find that Fisher's geometrical model can account for the observed patterns of fitness change and infer the parameters of this model that best fit the data, using Approximate Bayesian Computation. We estimate a genomic mutation rate of 0.01 per generation for fitness altering mutations, albeit with a large confidence interval, a mean fitness effect of mutations of ?0.01, and an effective number of traits nine in mutS^? E. coli. This framework can be extended to confront a broader range of models with data and test different classes of fitness landscape models.