The coupon collector and the suppressor mutation: estimating the number of compensatory mutations by maximum likelihood

Compensatory mutation occurs when a loss of fitness caused by a deleterious mutation is restored by its epistatic interaction with a second mutation at a different site in the genome. How many different compensatory mutations can act on a given deleterious mutation? Although this quantity is fundamentally important to understanding the evolutionary consequence of mutation and the genetic complexity of adaptation, it remains poorly understood. To determine the shape of the statistical distribution for the number of compensatory mutations per deleterious mutation, we have performed a maximum-likelihood analysis of experimental data collected from the suppressor mutation literature. Suppressor mutations are used widely to assess protein interactions and are under certain conditions equivalent to compensatory mutations. By comparing the maximum likelihood of a variety of candidate distribution functions, we established that an L-shaped gamma distribution (alpha=0.564, theta=21.01) is the most successful at explaining the collected data. This distribution predicts an average of 11.8 compensatory mutations per deleterious mutation. Furthermore, the success of the L-shaped gamma distribution is robust to variation in mutation rates among sites. We have detected significant differences among viral, prokaryotic, and eukaryotic data subsets in the number of compensatory mutations and also in the proportion of compensatory mutations that are intragenic. This is the first attempt to characterize the overall diversity of compensatory mutations, identifying a consistent and accurate prior distribution of compensatory mutation diversity for theoretical evolutionary models.     

The rate of compensatory mutation in the DNA bacteriophage phiX174

A compensatory mutation occurs when the fitness loss caused by one mutation is remedied by its epistatic interaction with a second mutation at a different site in the genome. This poorly understood biological phenomenon has important implications, not only for the evolutionary consequences of mutation, but also for the genetic complexity of adaptation. We have carried out the first direct experimental measurement of the average rate of compensatory mutation. An arbitrary selection of 21 missense substitutions with deleterious effects on fitness was introduced by site-directed mutagenesis into the bacteriophage phiX174. For each deleterious mutation, we evolved 8-16 replicate populations to determine the frequency at which a compensatory mutation, instead of the back mutation, was acquired to recover fitness. The overall frequency of compensatory mutation was approximately 70%. Deleterious mutations that were more severe were significantly more likely to be compensated for. Furthermore, experimental reversion of deleterious mutations revealed that compensatory mutations have deleterious effects in a wild-type background. A large diversity of intragenic compensatory mutations was identified from sequencing fitness-recovering genotypes. Subsequent analyses of intragenic mutation diversity revealed a significant degree of clustering around the deleterious mutation in the linear sequence and also within folded protein structures. Moreover, a likelihood analysis of mutation diversity predicts that, on average, a deleterious mutation can be compensated by about nine different intragenic compensatory mutations. We estimate that about half of all compensatory mutations are located extragenically in this organism.     

The nearly neutral and selection theories of molecular evolution under the fisher geometrical framework: substitution rate, population size, and complexity

The general theories of molecular evolution depend on relatively arbitrary assumptions about the relative distribution and rate of advantageous, deleterious, neutral, and nearly neutral mutations. The Fisher geometrical model (FGM) has been used to make distributions of mutations biologically interpretable. We explored an FGM-based molecular model to represent molecular evolutionary processes typically studied by nearly neutral and selection models, but in which distributions and relative rates of mutations with different selection coefficients are a consequence of biologically interpretable parameters, such as the average size of the phenotypic effect of mutations and the number of traits (complexity) of organisms. A variant of the FGM-based model that we called the static regime (SR) represents evolution as a nearly neutral process in which substitution rates are determined by a dynamic substitution process in which the population's phenotype remains around a suboptimum equilibrium fitness produced by a balance between slightly deleterious and slightly advantageous compensatory substitutions. As in previous nearly neutral models, the SR predicts a negative relationship between molecular evolutionary rate and population size; however, SR does not have the unrealistic properties of previous nearly neutral models such as the narrow window of selection strengths in which they work. In addition, the SR suggests that compensatory mutations cannot explain the high rate of fixations driven by positive selection currently found in DNA sequences, contrary to what has been previously suggested. We also developed a generalization of SR in which the optimum phenotype can change stochastically due to environmental or physiological shifts, which we called the variable regime (VR). VR models evolution as an interplay between adaptive processes and nearly neutral steady-state processes. When strong environmental fluctuations are incorporated, the process becomes a selection model in which evolutionary rate does not depend on population size, but is critically dependent on the complexity of organisms and mutation size. For SR as well as VR we found that key parameters of molecular evolution are linked by biological factors, and we showed that they cannot be fixed independently by arbitrary criteria, as has usually been assumed in previous molecular evolutionary models.     

The Rate of Compensatory Evolution

A two-locus model is presented to analyze the evolution of compensatory mutations occurring in stems of RNA secondary structures. Single mutations are assumed to be deleterious but harmless (neutral) in appropriate combinations. In proceeding under mutation pressure, natural selection and genetic drift from one fitness peak to another one, a population must therefore pass through a valley of intermediate deleterious states of individual fitness. The expected time for this transition is calculated using diffusion theory. The rate of compensatory evolution, k(c), is then defined as the inverse of the expected transition time. When selection against deleterious single mutations is strong, k(c) depends on the recombination fraction r between the two loci. Recombination generally reduces the rate of compensatory evolution because it breaks up favorable combinations of double mutants. For complete linkage, k(c) is given by the rate at which favorable combinations of double mutants are produced by compensatory mutation. For r>0, k(c) decreases exponentially with r. In contrast, k(c) becomes independent of r for weak selection. We discuss the dynamics of evolutionary substitutions of compensatory mutants in relation to WRIGHT's shifting balance theory of evolution and use our results to analyze the substitution process in helices of mRNA secondary structures.     

Viral RNA and evolved mutational robustness.

Many properties of organisms show great robustness against mutations. Whether this robustness is an evolved property or intrinsic to genetic systems is by and large unknown. An evolutionary origin of robustness would require a rethinking of key concepts in the field of molecular evolution, such as gene-specific neutral mutation rates, or the context-independence of deleterious mutations. We provide evidence that mutational robustness of the genome of RNA viruses to mutational changes in secondary structure has evolved. J. Exp. Zool. ( Mol. Dev. Evol.) 285:119-127, 1999.     

Parasites and deleterious mutations: interactions influencing the evolutionary maintenance of sex

The restrictive assumptions associated with purely genetic and purely ecological mechanisms suggest that neither of the two forces, in isolation, can offer a general explanation for the evolutionary maintenance of sex. Consequently, attention has turned to pluralistic models (i.e. models that apply both ecological and genetic mechanisms). Existing research has shown that combining mutation accumulation and parasitism allows restrictive assumptions about genetic and parasite parameter values to be relaxed while still predicting the maintenance of sex. However, several empirical studies have shown that deleterious mutations and parasitism can reduce fitness to a greater extent than would be expected if the two acted independently. We show how interactions between these genetic and ecological forces can completely reverse predictions about the evolution of reproductive modes. Moreover, we demonstrate that synergistic interactions between infection and deleterious mutations can render sex evolutionarily stable even when there is antagonistic epistasis among deleterious mutations, thereby widening the conditions for the evolutionary maintenance of sex.     

Evolution by Gene Duplication and Compensatory Advantageous Mutations

Relaxation of selective constraint is thought to play an important role for evolution by gene duplication, in connection with compensatory advantageous mutant substitutions. Models were investigated by incorporating gene duplication by unequal crossing over, selection, mutation and random genetic drift into Monte Carlo simulations. Compensatory advantageous mutations were introduced, and simulations were carried out with and without relaxation, when genes are redundant on chromosomes. Relaxation was introduced by assuming that deleterious mutants have no effect on fitness, so long as one or more genes free of such mutations remain in the array. Compensatory mutations are characterized by the intermediate deleterious step of their substitutions, and therefore relaxation by gene redundancy is important. Through extensive Monte Carlo simulations, it was found that compensatory mutant substitutions require relaxation in addition to gene duplication, when mutant effects are large. However when mutant effects are small, such that the product of selection coefficient and population size is around unity, evolution by compensatory mutation is enhanced by gene duplication even without relaxation.     

Compensatory mutations are repeatable and clustered within proteins

Compensatory mutations improve fitness in genotypes that contain deleterious mutations but have no beneficial effects otherwise. As such, compensatory mutations represent a very specific form of epistasis. We show that intragenic compensatory mutations occur non-randomly over gene sequence. Compensatory mutations are more likely to appear at some sites than others. Moreover, the sites of compensatory mutations are more likely than expected by chance to be near the site of the original deleterious mutation. Furthermore, compensatory mutations tend to occur more commonly in certain regions of the protein even when controlling for clustering around the site of the deleterious mutation. These results suggest that compensatory evolution at the protein level is partially predictable and may be convergent.     

Moving forward in determining the causes of mutations: the features of plants that make them suitable for assessing the impact of environmental factors and cell age

Currently, the types of factors that impact the mutation rate is a controversial issue. The marked attention towards identifying the factors that impact the genomic mutation rate is justified because mutations are the source of genetic variation underlying evolution and because many mutations have deleterious effects and can cause diseases. Although data showing correlations between germ cell division number and mutation rates (from epidemiological studies and molecular evolutionary rate analyses) have suggested that most mutations in animals are replication errors, this notion is highly debated and inconsistencies in the correlations suggest that other, replication-independent factors, could play an important role. Likely candidates include environmental parameters and cell age, but these issues have proved to be difficult to study using animals and in vitro systems, and consequently, very few or no data currently exist. The specific features of plants that make them powerful model systems for revealing the influence of the environment (natural environmental factors) and cell age on the spontaneous genomic mutation rate are discussed here. Overall, the evidence suggests that plants could be key biological systems for advancing our knowledge about how and why heritable mutations arise.     

Characterization of compensated mutations in terms of structural and physico-chemical properties

The study of the evolution of compensatory mechanisms among amino acids is paramount to our understanding of intramolecular epistatic interactions. It has been addressed from different points of view, for example much effort has been devoted to establish the number of compensatory mutations required per deleterious mutation. However, we still do not know how the nature of the compensated mutation determines the existence of compensatory mutations. Within this context, recent studies have produced several instances of an interesting phenomenon: human disease-associated residues may sometimes appear as wild-type residues in non-human proteins. This can be explained in terms of compensatory mutations, present in the non-human protein, which would neutralize the damage caused by the disease-associated residue. Therefore, comparison between these compensated mutations and non-compensated pathological mutations provides a simple approach to understand how the nature of the compensated deleterious mutation determines the existence of compensatory mutations. To address this issue, we have obtained a large set of compensated mutations and characterised them with a series of different properties. When comparing the resulting distributions with those from pathological mutations we find that in general compensated mutations are milder than pathological mutations. More precisely, we find that the probability that a compensatory mutation will evolve is directly related (i) to the location in the protein structure and (ii) to changes in physico-chemical properties (e.g. amino acid volume or hydrophobicity) of the compensated mutation.     

Rate of adaptive peak shifts with partial genetic robustness

How adaptive evolution occurs with individually deleterious but jointly beneficial mutations has been one of the major problems in population genetics theory. Adaptation in this case is commonly described as a population's escape from a local peak to a higher peak on Sewall Wright's fitness landscape. Recent molecular genetic and computational studies have suggested that genetic robustness can facilitate such peak shifts. If phenotypic expressions of new mutations are suppressed under genetic robustness, mutations that are otherwise deleterious can accumulate in the population as neutral variants. When the robustness is perturbed by an environmental change or a major mutation, these variants become exposed to natural selection. It is argued that this process promotes adaptation because allelic combinations enriched under genetic robustness can then be positively selected. Here, I propose simple two- and three-locus models of adaptation with partial genetic robustness as suggested by recent studies. The waiting time until the fixation of an adaptive haplotype was observed in stochastic simulations and compared to the expectation without robustness. It is shown that peak shifts can be delayed or accelerated depending on the conditions of genetic robustness. The evolutionary significance of these processes is discussed.     

Compensatory vs. pseudocompensatory evolution in molecular and developmental interactions

The evolution of molecules, developmental circuits, and new species are all characterized by the accumulation of incompatibilities between ancestors and descendants. When specific interactions between components are necessary at any of these levels, this requires compensatory coevolution. Theoretical treatments of compensatory evolution that only consider the endpoints predict that it should be rare because intermediate states are deleterious. However, empirical data suggest that compensatory evolution is common at all levels of molecular interaction. A general solution to this paradox is provided by plausible neutral or nearly neutral intermediates that possess informational redundancy. These intermediates provide an evolutionary path between coadapted allelic combinations. Although they allow incompatible end points to evolve, at no point was a deleterious mutation ever in need of compensation. As a result, what appears to be compensatory evolution may often actually be “pseudocompensatory.” Both theoretical and empirical studies indicate that pseudocompensation can speed the evolution of intergenic incompatibility, especially when driven by adaptation. However, under strong stabilizing selection the rate of pseudocompensatory evolution is still significant. Important examples of this process at work discussed here include the evolution of rRNA secondary structures, intra- and inter-protein interactions, and developmental genetic pathways. Future empirical work in this area should focus on comparing the details of intra- and intergenic interactions in closely related organisms.     

Functional origins of fitness effect-sizes of compensatory mutations in the DNA bacteriophage phiX174

Abstract Epistasis is an important and poorly understood aspect of mutations and strongly influences the evolutionary impact of genetic variation on adaptation and fitness. Although recent studies have begun to characterize the distribution of epistatic effects between mutations affecting fitness, there is currently a lack of empirical information on the underlying biological causes of these epistatic interactions. What are the functional constraints that determine the effectiveness of a compensatory mutation at restoring fitness? We have measured the effect‐sizes of 52 compensatory mutations affecting nine different deleterious mutations in the major capsid and spike proteins of the DNA bacteriophage X174. On average, an experimentally detectable compensatory mutation recovers about two‐thirds of the fitness cost of the preceding deleterious mutation. Variation in fitness effect‐sizes is only weakly associated with measures of the distance separating the deleterious and compensatory mutations in the amino acid sequence or the folded protein structure. However, there is a strong association of fitness effect‐size with the correlation in the effects of the mutations on the biochemical properties of amino acids. A compensatory mutation has the largest effect‐size, on average, when both the compensatory and deleterious mutations have radical effects on the overall biochemical make‐up of the amino acids. By examining the relative contributions of specific biochemical properties to variation in fitness effect‐size, we find that the area and charge of amino acids have a major influence, which suggests that the complexity of the amino acid phenotype is simplified by selection into a reduced number of phenotypic components.     

Evolution by small steps and rugged landscapes in the RNA virus phi6

Fisher's geometric model of adaptive evolution argues that adaptive evolution should generally result from the substitution of many mutations of small effect because advantageous mutations of small effect should be more common than those of large effect. However, evidence for both evolution by small steps and for Fisher's model has been mixed. Here we report supporting results from a new experimental test of the model. We subjected the bacteriophage phi6 to intensified genetic drift in small populations and caused viral fitness to decline through the accumulation of a deleterious mutation. We then propagated the mutated virus at a range of larger population sizes and allowed fitness to recover by natural selection. Although fitness declined in one large step, it was usually recovered in smaller steps. More importantly, step size during recovery was smaller with decreasing size of the recovery population. These results confirm Fisher's main prediction that advantageous mutations of small effect should be more common. We also show that the advantageous mutations of small effect are compensatory mutations whose advantage is conditional (epistatic) on the presence of the deleterious mutation, in which case the adaptive landscape of phi6 is likely to be very rugged.     

Genetic constraints on protein evolution

Evolution requires the generation and optimization of new traits ("adaptation") and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

Compensatory evolution in diploid populations

Compensatory mutations are individually deleterious but harmless in appropriate combinations either at more than two sites within a gene or on separate genes. Considering that dominance effects of selection and heterodimer formation of gene products may affect the rate of compensatory evolution, we investigate compensatory neutral mutation models for diploid populations. Our theoretical analysis on the average time until fixation of compensatory mutations shows that these factors play an important role in reducing the fixation time of compensatory mutations if mutation rates are not low. Compensatory evolution of heterodimers is shown to occur more easily if the deleterious effects of single mutants are recessive.     

Fixation probability of rare nonmutator and evolution of mutation rates

Although mutations drive the evolutionary process, the rates at which the mutations occur are themselves subject to evolutionary forces. Our purpose here is to understand the role of selection and random genetic drift in the evolution of mutation rates, and we address this question in asexual populations at mutation‐selection equilibrium neglecting selective sweeps. Using a multitype branching process, we calculate the fixation probability of a rare nonmutator in a large asexual population of mutators and find that a nonmutator is more likely to fix when the deleterious mutation rate of the mutator population is high. Compensatory mutations in the mutator population are found to decrease the fixation probability of a nonmutator when the selection coefficient is large. But, surprisingly, the fixation probability changes nonmonotonically with increasing compensatory mutation rate when the selection is mild. Using these results for the fixation probability and a drift‐barrier argument, we find a novel relationship between the mutation rates and the population size. We also discuss the time to fix the nonmutator in an adapted population of asexual mutators, and compare our results with experiments.     

PERSPECTIVE: SPONTANEOUS DELETERIOUS MUTATION

Mildly deleterious mutation has been invoked as a leading explanation for a diverse array of observations in evolutionary genetics and molecular evolution and is thought to be a significant risk of extinction for small populations. However, much of the empirical evidence for the deleterious-mutation process derives from studies of Drosophila melanogaster, some of which have been called into question. We review a broad array of data that collectively support the hypothesis that deleterious mutations arise in flies at rate of about one per individual per generation, with the average mutation decreasing fitness by about only 2% in the heterozygous state. Empirical evidence from microbes, plants, and several other animal species provide further support for the idea that most mutations have only mildly deleterious effects on fitness, and several other species appear to have genomic mutation rates that are of the order of magnitude observed in Drosophila. However, there is mounting evidence that some organisms have genomic deleterious mutation rates that are substantially lower than one per individual per generation. These lower rates may be at least partially reconciled with the Drosophila data by taking into consideration the number of germline cell divisions per generation. To fully resolve the existing controversy over the properties of spontaneous mutations, a number of issues need to be clarified. These include the form of the distribution of mutational effects and the extent to which this is modified by the environmental and genetic background and the contribution of basic biological features such as generation length and genome size to interspecific differences in the genomic mutation rate. Once such information is available, it should be possible to make a refined statement about the long-term impact of mutation on the genetic integrity of human populations subject to relaxed selection resulting from modern medical procedures.