Evolution of sex in RNA viruses

The distribution of deleterious mutations in a population of organisms is determined by the opposing effects of two forces, mutation pressure and selection. If mutation rates are high, the resulting mutation-selection balance can generate a substantial mutational load in the population. Sex can be advantageous to organisms experiencing high mutation rates because it can either buffer the mutation-selection balance from genetic drift, thus preventing any increases in the mutational load (Muller, 1964: Mut. Res. 1, 2), or decrease the mutational load by increasing the efficiency of selection (Crow, 1970: Biomathematics 1, 128). Muller's hypothesis assumes that deleterious mutations act independently, whereas Crow's hypothesis assumes that deleterious mutations interact synergistically, i.e., the acquisition of a deleterious mutation is proportionately more harmful to a genome with many mutations than it is to a genome with a few mutations. RNA viruses provide a test for these two hypotheses because they have extremely high mutation rates and appear to have evolved specific adaptations to reproduce sexually. Population genetic models for RNA viruses show that Muller's and Crow's hypotheses are also possible explanations for why sex is advantageous to these viruses. A re-analysis of published data on RNA viruses that are cultured by undiluted passage suggests that deleterious mutations in such viruses interact synergistically and that sex evolved there as a mechanism to reduce the mutational load.     

Deleterious mutations and the genetic variance of male fitness components in Mimulus guttatus

Deleterious mutations are relevant to a broad range of questions in genetics and evolutionary biology. I present an application of the "biometric method" for estimating mutational parameters for male fitness characters of the yellow monkeyflower, Mimulus guttatus. The biometric method rests on two critical assumptions. The first is that experimental inbreeding changes genotype frequencies without changing allele frequencies; i.e., there is no genetic purging during the experiment. I satisfy this condition by employing a breeding design in which the parents are randomly extracted, fully homozygous inbred lines. The second is that all genetic variation is attributable to deleterious mutations maintained in mutation-selection balance. I explicitly test this hypothesis using likelihood ratios. Of the three deleterious mutation models tested, the first two are rejected for all characters. The failure of these models is due to an excess of additive genetic variation relative to the expectation under mutation-selection balance. The third model is not rejected for either of two log-transformed male fitness traits. However, this model imposes only "weak conditions" and is not sufficiently detailed to provide estimates for mutational parameters. The implication is that, if biometric methods are going to yield useful parameter estimates, they will need to consider mutational models more complicated than those typically employed in experimental studies.     

Dynamic mutation-selection balance as an evolutionary attractor

The vast majority of mutations are deleterious and are eliminated by purifying selection. Yet in finite asexual populations, purifying selection cannot completely prevent the accumulation of deleterious mutations due to Muller's ratchet: once lost by stochastic drift, the most-fit class of genotypes is lost forever. If deleterious mutations are weakly selected, Muller's ratchet can lead to a rapid degradation of population fitness. Evidently, the long-term stability of an asexual population requires an influx of beneficial mutations that continuously compensate for the accumulation of the weakly deleterious ones. Hence any stable evolutionary state of a population in a static environment must involve a dynamic mutation-selection balance, where accumulation of deleterious mutations is on average offset by the influx of beneficial mutations. We argue that such a state can exist for any population size N and mutation rate U and calculate the fraction of beneficial mutations, ε, that maintains the balanced state. We find that a surprisingly low ε suffices to achieve stability, even in small populations in the face of high mutation rates and weak selection, maintaining a well-adapted population in spite of Muller's ratchet. This may explain the maintenance of mitochondria and other asexual genomes.     

Dynamics of inbreeding depression due to deleterious mutations in small populations: mutation parameters and inbreeding rate

A multilocus stochastic model is developed to simulate the dynamics of mutational load in small populations of various sizes. Old mutations sampled from a large ancestral population at mutation-selection balance and new mutations arising each generation are considered jointly, using biologically plausible lethal and deleterious mutation parameters. The results show that inbreeding depression and the number of lethal equivalents due to partially recessive mutations can be partly purged from the population by inbreeding, and that this purging mainly involves lethals or detrimentals of large effect. However, fitness decreases continuously with inbreeding, due to increased fixation and homozygosity of mildly deleterious mutants, resulting in extinctions of very small populations with low reproductive rates. No optimum inbreeding rate or population size exists for purging with respect to fitness (viability) changes, but there is an optimum inbreeding rate at a given final level of inbreeding for reducing inbreeding depression or the number of lethal equivalents. The interaction between selection against partially recessive mutations and genetic drift in small populations also influences the rate of decay of neutral variation. Weak selection against mutants relative to genetic drift results in apparent overdominance and thus an increase in effective size (Ne) at neutral loci, and strong selection relative to drift leads to a decrease in Ne due to the increased variance in family size. The simulation results and their implications are discussed in the context of biological conservation and tests for purging.     

Antagonistic coevolution with parasites increases the cost of host deleterious mutations

The fitness consequences of deleterious mutations are sometimes greater when individuals are parasitized, hence parasites may result in the more rapid purging of deleterious mutations from host populations. The significance of host deleterious mutations when hosts and parasites antagonistically coevolve (reciprocal evolution of host resistance and parasite infectivity) has not previously been experimentally investigated. We addressed this by coevolving the bacterium Pseudomonas fluorescens and a parasitic bacteriophage in laboratory microcosms, using bacteria with high and low mutation loads. Directional coevolution between bacterial resistance and phage infectivity occurred in all populations. Bacterial population fitness, as measured by competition experiments with ancestral genotypes in the absence of phage, declined with time spent coevolving. However, this decline was significantly more rapid in bacteria with high mutation loads, suggesting the cost of bacterial resistance to phage was greater in the presence of deleterious mutations (synergistic epistasis). As such, resistance to phage was more costly to evolve in the presence of a high mutation load. Consistent with these data, bacteria with high mutation loads underwent less rapid directional coevolution with their phage populations, and showed lower levels of resistance to their coevolving phage populations. These data suggest that coevolution with parasites increases the rate at which deleterious mutations are purged from host populations.     

Deleterious Passengers in Adapting Populations

Most new mutations are deleterious and are eventually eliminated by natural selection. But in an adapting population, the rapid amplification of beneficial mutations can hinder the removal of deleterious variants in nearby regions of the genome, altering the patterns of sequence evolution. Here, we analyze the interactions between beneficial “driver” mutations and linked deleterious “passengers” during the course of adaptation. We derive analytical expressions for the substitution rate of a deleterious mutation as a function of its fitness cost, as well as the reduction in the beneficial substitution rate due to the genetic load of the passengers. We find that the fate of each deleterious mutation varies dramatically with the rate and spectrum of beneficial mutations and the deleterious substitution rate depends nonmonotonically on the population size and the rate of adaptation. By quantifying this dependence, our results allow us to estimate which deleterious mutations will be likely to fix and how many of these mutations must arise before the progress of adaptation is significantly reduced.     

Mutation-selection balance: ancestry,load, and maximum principle

We analyze the equilibrium behavior of deterministic haploid mutation-selection models. To this end, both the forward and the time-reversed evolution processes are considered. The stationary state of the latter is called the ancestral distribution, which turns out as a key for the study of mutation-selection balance. We find that the ancestral genotype frequencies determine the sensitivity of the equilibrium mean fitness to changes in the corresponding fitness values and discuss implications for the evolution of mutational robustness. We further show that the difference between the ancestral and the population mean fitness, termed mutational loss, provides a measure for the sensitivity of the equilibrium mean fitness to changes in the mutation rate. The interrelation of the loss and the mutation load is discussed. For a class of models in which the number of mutations in an individual is taken as the trait value, and fitness is a function of the trait, we use the ancestor formulation to derive a simple maximum principle, from which the mean and variance of fitness and the trait may be derived; the results are exact for a number of limiting cases, and otherwise yield approximations which are accurate for a wide range of parameters. These results are applied to threshold phenomena caused by the interplay of selection and mutation (known as error thresholds). They lead to a clarification of concepts, as well as criteria for the existence of error thresholds.     

Mate choice, monogamy and mutational load

The classical Haldane mutational load for dominant mutations is 1 - e-alpha (where alpha = the mutation rate of the genome). This is independent of the fitness of the mutant heterozygotes. As alpha may be between 0.9 and 1.8 this may be too high for low fecundity species to bear. The load may be greatly reduced for mildly deleterious mutations if female choice (for males with few or no mutations) is present in a polygamous species. It is shown here that if females and males prefer partners with few or no mutations then the load may also be reduced in monogamous species. As in the polygamous model the effect is greatest if mutations are mildly deleterious.     

The Interaction between Selection,Demography and Selfing and How It Affects Population Viability

Population extinction due to the accumulation of deleterious mutations has only been considered to occur at small population sizes, large sexual populations being expected to efficiently purge these mutations. However, little is known about how the mutation load generated by segregating mutations affects population size and, eventually, population extinction. We propose a simple analytical model that takes into account both the demographic and genetic evolution of populations, linking population size, density dependence, the mutation load, and self-fertilisation. Analytical predictions were found to be relatively good predictors of population size and probability of population viability when verified using an explicit individual based stochastic model. We show that initially large populations do not always reach mutation-selection balance and can go extinct due to the accumulation of segregating deleterious mutations. Population survival depends not only on the relative fitness and demographic stochasticity, but also on the interaction between the two. When deleterious mutations are recessive, self-fertilisation affects viability non-monotonically and genomic cold-spots could favour the viability of outcrossing populations.     

Most rare missense alleles are deleterious in humans: implications for complex disease and association studies

The accumulation of mildly deleterious missense mutations in individual human genomes has been proposed to be a genetic basis for complex diseases. The plausibility of this hypothesis depends on quantitative estimates of the prevalence of mildly deleterious de novo mutations and polymorphic variants in humans and on the intensity of selective pressure against them. We combined analysis of mutations causing human Mendelian diseases, of human-chimpanzee divergence, and of systematic data on human genetic variation and found that ~20% of new missense mutations in humans result in a loss of function, whereas ~27% are effectively neutral. Thus, the remaining 53% of new missense mutations have mildly deleterious effects. These mutations give rise to many low-frequency deleterious allelic variants in the human population, as is evident from a new data set of 37 genes sequenced in >1,500 individual human chromosomes. Surprisingly, up to 70% of low-frequency missense alleles are mildly deleterious and are associated with a heterozygous fitness loss in the range 0.001-0.003. Thus, the low allele frequency of an amino acid variant can, by itself, serve as a predictor of its functional significance. Several recent studies have reported a significant excess of rare missense variants in candidate genes or pathways in individuals with extreme values of quantitative phenotypes. These studies would be unlikely to yield results if most rare variants were neutral or if rare variants were not a significant contributor to the genetic component of phenotypic inheritance. Our results provide a justification for these types of candidate-gene (pathway) association studies and imply that mutation-selection balance may be a feasible evolutionary mechanism underlying some common diseases.     

Estimating selection on nonsynonymous mutations

The distribution of mutational effects on fitness is of fundamental importance for many aspects of evolution. We develop two methods for characterizing the fitness effects of deleterious, nonsynonymous mutations, using polymorphism data from two related species. These methods also provide estimates of the proportion of amino acid substitutions that are selectively favorable, when combined with data on between-species sequence divergence. The methods are applicable to species with different effective population sizes, but that share the same distribution of mutational effects. The first, simpler, method assumes that diversity for all nonneutral mutations is given by the value under mutation-selection balance, while the second method allows for stronger effects of genetic drift and yields estimates of the parameters of the probability distribution of mutational effects. We apply these methods to data on populations of Drosophila miranda and D. pseudoobscura and find evidence for the presence of deleterious nonsynonymous mutations, mostly with small heterozygous selection coefficients (a mean of the order of 10(-5) for segregating variants). A leptokurtic gamma distribution of mutational effects with a shape parameter between 0.1 and 1 can explain observed diversities, in the absence of a separate class of completely neutral nonsynonymous mutations. We also describe a simple approximate method for estimating the harmonic mean selection coefficient from diversity data on a single species.     

Haldane's sieve and adaptation from the standing genetic variation

We consider populations that adapt to a sudden environmental change by fixing alleles found at mutation-selection balance. In particular, we calculate probabilities of fixation for previously deleterious alleles, ignoring the input of new mutations. We find that "Haldane's sieve"--the bias against the establishment of recessive beneficial mutations--does not hold under these conditions. Instead probabilities of fixation are generally independent of dominance. We show that this result is robust to patterns of sex expression for both X-linked and autosomal loci. We further show that adaptive evolution is invariably slower at X-linked than autosomal loci when evolution begins from mutation-selection balance. This result differs from that obtained when adaptation uses new mutations, a finding that may have some bearing on recent attempts to distinguish between hitchhiking and background selection by contrasting the molecular population genetics of X-linked vs. autosomal loci. Last, we suggest a test to determine whether adaptation used new mutations or previously deleterious alleles from the standing genetic variation.     

Beneficial mutation selection balance and the effect of linkage on positive selection

When beneficial mutations are rare, they accumulate by a series of selective sweeps. But when they are common, many beneficial mutations will occur before any can fix, so there will be many different mutant lineages in the population concurrently. In an asexual population, these different mutant lineages interfere and not all can fix simultaneously. In addition, further beneficial mutations can accumulate in mutant lineages while these are still a minority of the population. In this article, we analyze the dynamics of such multiple mutations and the interplay between multiple mutations and interference between clones. These result in substantial variation in fitness accumulating within a single asexual population. The amount of variation is determined by a balance between selection, which destroys variation, and beneficial mutations, which create more. The behavior depends in a subtle way on the population parameters: the population size, the beneficial mutation rate, and the distribution of the fitness increments of the potential beneficial mutations. The mutation-selection balance leads to a continually evolving population with a steady-state fitness variation. This variation increases logarithmically with both population size and mutation rate and sets the rate at which the population accumulates beneficial mutations, which thus also grows only logarithmically with population size and mutation rate. These results imply that mutator phenotypes are less effective in larger asexual populations. They also have consequences for the advantages (or disadvantages) of sex via the Fisher-Muller effect; these are discussed briefly.     

Beneficial mutations, hitchhiking and the evolution of mutation rates in sexual populations

Natural selection acts in three ways on heritable variation for mutation rates. A modifier allele that increases the mutation rate is (i) disfavored due to association with deleterious mutations, but is also favored due to (ii) association with beneficial mutations and (iii) the reduced costs of lower fidelity replication. When a unique beneficial mutation arises and sweeps to fixation, genetic hitchhiking may cause a substantial change in the frequency of a modifier of mutation rate. In previous studies of the evolution of mutation rates in sexual populations, this effect has been underestimated. This article models the long-term effect of a series of such hitchhiking events and determines the resulting strength of indirect selection on the modifier. This is compared to the indirect selection due to deleterious mutations, when both types of mutations are randomly scattered over a given genetic map. Relative to an asexual population, increased levels of recombination reduce the effects of beneficial mutations more rapidly than those of deleterious mutations. However, the role of beneficial mutations in determining the evolutionarily stable mutation rate may still be significant if the function describing the cost of high-fidelity replication has a shallow gradient.     

MUTATION AND EXTINCTION: THE ROLE OF VARIABLE MUTATIONAL EFFECTS,SYNERGISTIC EPISTASIS,BENEFICIAL MUTATIONS,AND DEGREE OF OUTCROSSING

Recent theoretical studies have illustrated the potential role of spontaneous deleterious mutation as a cause of extinction in small populations. However, these studies have not addressed several genetic issues, which can in principle have a substantial influence on the risk of extinction. These include the presence of synergistic epistasis, which can reduce the rate of mutation accumulation by progressively magnifying the selective effects of mutations, and the occurrence of beneficial mutations, which can offset the effects of previous deleterious mutations. In stochastic simulations of small populations (effective sizes on the order of 100 or less), we show that both synergistic epistasis and the rate of beneficial mutation must be unrealistically high to substantially reduce the risk of extinction due to random fixation of deleterious mutations. However, in analytical calculations based on diffusion theory, we show that in large, outcrossing populations (effective sizes greater than a few hundred), very low levels of beneficial mutation are sufficient to prevent mutational decay. Further simulation results indicate that in populations small enough to be highly vulnerable to mutational decay, variance in deleterious mutational effects reduces the risk of extinction, assuming that the mean deleterious mutational effect is on the order of a few percent or less. We also examine the magnitude of outcrossing that is necessary to liberate a predominantly selfing population from the threat of long-term mutational deterioration. The critical amount of outcrossing appears to be greater than is common in near-obligately selfing plant species, supporting the contention that such species are generally doomed to extinction via random drift of new mutations. Our results support the hypothesis that a long-term effective population size in the neighborhood of a few hundred individuals defines an approximate threshold, below which outcrossing populations are vulnerable to extinction via fixation of deleterious mutations, and above which immunity is acquired.     

Mutation-selection balance accounting for genetic variation for viability in Drosophila melanogaster as deduced from an inbreeding and artificial selection experiment

We carried out an experiment of inbreeding and upward artificial selection for egg-to-adult viability in a recently captured population of Drosophila melanogaster, as well as computer simulations of the experimental design, in order to obtain information on the nature of genetic variation for this important fitness component. The inbreeding depression was linear with a rate of 0.70 +/- 0.11% of the initial mean per 1% increase in inbreeding coefficient, and the realized heritability was 0.06 +/- 0.07. We compared the empirical observations of inbreeding depression and selection response with computer simulations assuming a balance between the occurrence of partially recessive deleterious mutations and their elimination by selection. Our results suggest that a model assuming mutation-selection balance with realistic mutational parameters can explain the genetic variation for viability in the natural population studied. Several mutational models are incompatible with some observations and can be discarded. Mutational models assuming a low rate of mutations of large average effect and highly recessive gene action, and others assuming a high rate of mutations of small average effect and close to additive gene action, are compatible with all the observations.     

Effects of spontaneous mutation accumulation on sex ratio traits in a parasitoid wasp

Sex allocation theory has proved extremely successful at predicting when individuals should adjust the sex of their offspring in response to environmental conditions. However, we know rather little about the underlying genetics of sex ratio or how genetic architecture might constrain adaptive sex-ratio behavior. We examined how mutation influenced genetic variation in the sex ratios produced by the parasitoid wasp Nasonia vitripennis. In a mutation accumulation experiment, we determined the mutability of sex ratio, and compared this with the amount of genetic variation observed in natural populations. We found that the mutability (h(2)(m)) ranges from 0.001 to 0.002, similar to estimates for life-history traits in other organisms. These estimates suggest one mutation every 5-60 generations, which shift the sex ratio by approximately 0.01 (proportion males). In this and other studies, the genetic variation in N. vitripennis sex ratio ranged from 0.02 to 0.17 (broad-sense heritability, H(2)). If sex ratio is maintained by mutation-selection balance, a higher genetic variance would be expected given our mutational parameters. Instead, the observed genetic variance perhaps suggests additional selection against sex-ratio mutations with deleterious effects on other fitness traits as well as sex ratio (i.e., pleiotropy), as has been argued to be the case more generally.     

Nonequivalent Loci and the distribution of mutant effects

It has been observed repeatedly that the distribution of new mutations of a quantitative trait has a kurtosis (a statistical measure of the distribution's shape) that is systematically larger than that of a normal distribution. Here we suggest that rather than being a property of individual loci that control the trait, the enhanced kurtosis is highly likely to be an emergent property that arises directly from the loci being mutationally nonequivalent. We present a method of incorporating nonequivalent loci into quantitative genetic modeling and give an approximate relation between the kurtosis of the mutant distribution and the degree of mutational nonequivalence of loci. We go on to ask whether incorporating the experimentally observed kurtosis through nonequivalent loci, rather than at locus level, affects any biologically important conclusions of quantitative genetic modeling. Concentrating on the maintenance of quantitative genetic variation by mutation-selection balance, we conclude that typically nonequivalent loci yield a genetic variance that is of order 10% smaller than that obtained from the previous approaches. For large populations, when the kurtosis is large, the genetic variance may be <50% of the result of equivalent loci, with Gaussian distributions of mutant effects.     

Epistasis in polygenic traits and the evolution of genetic architecture under stabilizing selection

We consider the effects of epistasis in a polygenic trait in the balance of mutation and stabilizing selection. The main issues are the genetic variation maintained in equilibrium and the evolution of the mutational effect distribution. The model assumes symmetric mutation and a continuum of alleles at all loci. Epistasis is modeled proportional to pairwise products of the single-locus effects. A general analytical formalism is developed. Assuming linkage equilibrium, we derive results for the equilibrium mutation load and the genetic and mutational variance in the house of cards and the Gaussian approximation. The additive genetic variation maintained in mutation-selection balance is reduced by any pattern of the epistatic interactions. The mutational variance, in contrast, is often increased. Large differences in mutational effects among loci emerge, and a negative correlation among (standard mean) locus mutation effects and mutation rates is predicted. Contrary to the common view since Waddington, we find that stabilizing selection in general does not lead to canalization of the trait. We propose that canalization as a target of selection instead occurs at the genic level. Here, primarily genes with a high mutation rate are buffered, often at the cost of decanalization of other genes. An intuitive interpretation of this view is given in the discussion.     

Probing the evolution of senescence inDrosophila melanogaster withP-element tagging

Natural populations host a wealth of genetic variation in longevity and age-specific schedules of reproduction. This variation provides critical information for inferring the evolutionary origin of senescence. Patterns of mutational effects on age-specific fecundity and survival provide additional insight to distinguish alternative models of senescence. In this study,P-elements bearing thewhite minigene were inserted at random into a common genetic background, generating lines ofD. melanogaster with single, stable transposon inserts. A series of 48 single-P-element lines revealed statistically significant heterogeneity in both longevity and fecundity. Longevity and early fecundity were only weakly positively correlated (r=0.286,P=0.0398). Both the pooled sample and 30 of the individual lines exhibited a leveling of age-specific mortality at advanced ages, in opposition to the classical demographic models. To the extent that these mutational effects are representative of naturally-occurring mutations in heterogeneous populations, this result presents a problem for the evolutionary theory of senescence. Natural selection is inefficient at removing deleterious mutations that are expressed only at late ages, and selection may not differentiate between mutations whose effects on longevity are post-reproductive. A leveling of the mortality rate would also be seen if mutations whose expression is delayed until very late simply do not occur. A simulation of mutation-selection balance among the 48P-element tagged lines shows that the mean longevity declines monotonically with increasing mutation rate, consistent with the mutation-accumulation model.