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.     

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.     

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.     

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.     

Purging deleterious mutations under self fertilization: paradoxical recovery in fitness with increasing mutation rate in Caenorhabditis elegans

Background

The accumulation of deleterious mutations can drastically reduce population mean fitness. Self-fertilization is thought to be an effective means of purging deleterious mutations. However, widespread linkage disequilibrium generated and maintained by self-fertilization is predicted to reduce the efficacy of purging when mutations are present at multiple loci.

Methodology/Principal Findings

We tested the ability of self-fertilizing populations to purge deleterious mutations at multiple loci by exposing obligately self-fertilizing populations of Caenorhabditis elegans to a range of elevated mutation rates and found that mutations accumulated, as evidenced by a reduction in mean fitness, in each population. Therefore, purging in obligate selfing populations is overwhelmed by an increase in mutation rate. Surprisingly, we also found that obligate and predominantly self-fertilizing populations exposed to very high mutation rates exhibited consistently greater fitness than those subject to lesser increases in mutation rate, which contradicts the assumption that increases in mutation rate are negatively correlated with fitness. The high levels of genetic linkage inherent in self-fertilization could drive this fitness increase.

Conclusions

Compensatory mutations can be more frequent under high mutation rates and may alleviate a portion of the fitness lost due to the accumulation of deleterious mutations through epistatic interactions with deleterious mutations. The prolonged maintenance of tightly linked compensatory and deleterious mutations facilitated by self-fertilization may be responsible for the fitness increase as linkage disequilibrium between the compensatory and deleterious mutations preserves their epistatic interaction.     

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.     

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.     

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.     

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.     

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 adaptation to the deleterious effect of antibiotic resistance in Salmonella typhimurium

Most chromosomal mutations that cause antibiotic resistance impose fitness costs on the bacteria. This biological cost can often be reduced by compensatory mutations. In Salmonella typhimurium, the nucleotide substitution AAA42 --> AAC in the rpsL gene confers resistance to streptomycin. The resulting amino acid substitution (K42N) in ribosomal protein S12 causes an increased rate of ribosomal proofreading and, as a result, the rate of protein synthesis, bacterial growth and virulence are decreased. Eighty-one independent lineages of the low-fitness, K42N mutant were evolved in the absence of antibiotic to ameliorate the costs. From the rate of fixation of compensated mutants and their fitness, the rate of compensatory mutations was estimated to be > or = 10-7 per cell per generation. The size of the population bottleneck during evolution affected fitness of the adapted mutants: a larger bottleneck resulted in higher average fitness. Only four of the evolved lineages contained streptomycin-sensitive revertants. The remaining 77 lineages contained mutants that were still fully streptomycin resistant, had retained the original resistance mutation and also acquired compensatory mutations. Most of the compensatory mutations, resulting in at least 35 different amino acid substitutions, were novel single-nucleotide substitutions in the rpsD, rpsE, rpsL or rplS genes encoding the ribosomal proteins S4, S5, S12 and L19 respectively. Our results show that the deleterious effects of a resistance mutation can be compensated by an unexpected variety of mutations.     

AMELIORATION OF THE DELETERIOUS PLEIOTROPIC EFFECTS OF AN ADAPTIVE MUTATION IN BACILLUS SUBTILIS

The deleterious pleiotropic effects of an adaptive mutation may be ameliorated by one of two modes of evolution: (1) by replacement, in which an adaptive mutation with harmful pleiotropic effects is replaced by one that confers an equal benefit but at less cost; or (2) by compensatory evolution, in which natural selection favors modifiers at other loci that compensate for the deleterious effects of the mutant allele. In this study, we have measured the potential of these two modes of evolution to ameliorate the deleterious pleiotropic effects of resistance to the antibiotic rifampicin in the soil bacterium Bacillus subtilis. One approach was to measure the fitness cost of a series of spontaneous rifampicin-resistance mutations from each of several strains. The potential for amelioration by the replacement mode was estimated by the variation in fitness cost among the mutants of a single strain. Another approach was to introduce a series of different rifampicin-resistance alleles into a diversity of strains, and to measure the fitness cost of rifampicin resistance for each allele-by-strain combination. The potential for amelioration by the replacement mode was estimated by the variation in fitness costs among rifampicin-resistance alleles; the potential for compensatory evolution was estimated by variation in the fitness cost of rifampicin resistance among strains. This study has shown that the cost of rifampicin resistance may be ameliorated by both the compensatory and replacement modes.     

Note on the computation of critical effective population sizes.

This work extends the work of Whitlock in examining the critical effective population sizes from the fixation of both deleterious and beneficial mutations under drift and selection to prevent mutation breakdown of the population. The validity of approximations for the probability of fixation depends on the nature of the assumed distribution for the fitness effect of both types of mutations. Using no approximation for the probability of fixation and assuming a heavy tailed fitness effect distribution, the current model indicates that the coefficients of variation for the fitness effect distributions of both types of mutations and the fitness effect distribution mean for the beneficial mutations are important predictors of the critical effective population size. The current model further predicts that very small populations can be sustained if the fitness effect variances for both types of mutations and the mean for beneficial mutations are large.     

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.     

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.     

Mutation, multilevel selection, and the evolution of propagule size during the origin of multicellularity

Evolutionary transitions require the organization of genetic variation at two (or more) levels of selection so that fitness heritability may emerge at the new level. In this article, we consider the consequences for fitness variation and heritability of two of the main modes of reproduction used in multicellular organisms: vegetative reproduction and single-cell reproduction. We study a model where simple cell colonies reproduce by fragments or propagules of differing size, with mutations occurring during colony growth. Mutations are deleterious at the colony level but can be advantageous or deleterious at the cell level ("selfish" or "uniformly deleterious" mutants). Fragment size affects fitness in two ways: through a direct effect on adult group size (which in turn affects fitness) and by affecting the within- and between-group variances and opportunity for selection on mutations at the two levels. We show that the evolution of fragment size is determined primarily by its direct effects on group size except when mutations are selfish. When mutations are selfish, smaller propagule size may be selected, including single-cell reproduction, even though smaller propagule size has a direct fitness cost by virtue of producing smaller organisms, that is, smaller adult cell groups.     

Characterization of deleterious mutations in outcrossing populations.

Deng and Lynch recently proposed estimating the rate and effects of deleterious genomic mutations from changes in the mean and genetic variance of fitness upon selfing/outcrossing in outcrossing/highly selfing populations. The utility of our original estimation approach is limited in outcrossing populations, since selfing may not always be feasible. Here we extend the approach to any form of inbreeding in outcrossing populations. By simulations, the statistical properties of the estimation under a common form of inbreeding (sib mating) are investigated under a range of biologically plausible situations. The efficiencies of different degrees of inbreeding and two different experimental designs of estimation are also investigated. We found that estimation using the total genetic variation in the inbred generation is generally more efficient than employing the genetic variation among the mean of inbred families, and that higher degree of inbreeding employed in experiments yields higher power for estimation. The simulation results of the magnitude and direction of estimation bias under variable or epistatic mutation effects may provide a basis for accurate inferences of deleterious mutations. Simulations accounting for environmental variance of fitness suggest that, under full-sib mating, our extension can achieve reasonably well an estimation with sample sizes of only approximately 2000-3000.     

Mutation accumulation in populations of varying size: large effect mutations cause most mutational decline in the rotifer Brachionus calyciflorus under UV‐C radiation

Theory predicts that fitness decline via mutation accumulation will depend on population size, but there are only a few direct tests of this key idea. To gain a qualitative understanding of the fitness effect of new mutations, we performed a mutation accumulation experiment with the facultative sexual rotifer Brachionus calyciflorus at six different population sizes under UV‐C radiation. Lifetime reproduction assays conducted after ten and sixteen UV‐C radiations showed that while small populations lost fitness, fitness losses diminished rapidly with increasing population size. Populations kept as low as 10 individuals were able to maintain fitness close to the nonmutagenized populations throughout the experiment indicating that selection was able to remove the majority of large effect mutations in small populations. Although our results also seem to imply that small populations are effectively immune to mutational decay, we caution against this interpretation. Given sufficient time, populations of moderate to large size can experience declines in fitness from accumulating weakly deleterious mutations as demonstrated by fitness estimates from simulations and, tentatively, from a long‐term experiment with populations of moderate size. There is mounting evidence to suggest that mutational distributions contain a heavier tail of large effects. Our results suggest that this is also true when the mutational spectrum is altered by UV radiation.     

Sexual selection on spontaneous mutations strengthens the between‐sex genetic correlation for fitness

A proposed benefit to sexual selection is that it promotes purging of deleterious mutations from populations. For this benefit to be realized, sexual selection, which is usually stronger on males, must purge mutations deleterious to both sexes. Here, we experimentally test the hypothesis that sexual selection on males purges deleterious mutations that affect both male and female fitness. We measured male and female fitness in two panels of spontaneous mutation‐accumulation lines of the fly, Drosophila serrata, each established from a common ancestor. One panel of mutation accumulation lines limited both natural and sexual selection (LS lines), whereas the other panel limited natural selection, but allowed sexual selection to operate (SS lines). Although mutation accumulation caused a significant reduction in male and female fitness in both the LS and SS lines, sexual selection had no detectable effect on the extent of the fitness reduction. Similarly, despite evidence of mutational variance for fitness in males and females of both treatments, sexual selection had no significant impact on the amount of mutational genetic variance for fitness. However, sexual selection did reshape the between‐sex correlation for fitness: significantly strengthening it in the SS lines. After 25 generations, the between‐sex correlation for fitness was positive but considerably less than one in the LS lines, suggesting that, although most mutations had sexually concordant fitness effects, sex‐limited, and/or sex‐biased mutations contributed substantially to the mutational variance. In the SS lines this correlation was strong and could not be distinguished from unity. Individual‐based simulations that mimick the experimental setup reveal two conditions that may drive our results: (1) a modest‐to‐large fraction of mutations have sex‐limited (or highly sex‐biased) fitness effects, and (2) the average fitness effect of sex‐limited mutations is larger than the average fitness effect of mutations that affect both sexes similarly.     

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.