MUTATIONAL MELTDOWNS IN SEXUAL POPULATIONS

Although it is widely acknowledged that the gradual accumulation of mildly deleterious mutations is an important source of extinction for asexual populations, it is generally assumed that this process is of little relevance to sexual species. Here we present results, based on computer simulations and supported by analytical approximations, that indicate that mutation accumulation in small, random-mating monoecious populations can lead to mean extinction times less than a few hundred to a few thousand generations. Unlike the situation in obligate asexuals in which the mean time to extinction (t?_e) increases more slowly than linearly with the population carrying capacity (K), t?_e increases approximately exponentially with K in outcrossing sexual populations. The mean time to extinction for obligately selfing populations is shown to be equivalent to that for asexual populations of the same size, but with half the mutation rate and twice the mutational effect; this suggests that obligate selfing, like obligate asexuality, is inviable as a long-term reproductive strategy. Under all mating systems, the mean time to extinction increases relatively slowly with the logarithm of fecundity, and mutations with intermediate effects (similar to those observed empirically) cause the greatest risk of extinction. Because our analyses ignore sources of demographic and environmental stochasticity, which have synergistic effects that exacerbate the accumulation of deleterious mutations, our results should yield liberal upper bounds to the mean time to extinction caused by mutational degradation. Thus, deleterious mutation accumulation cannot be ruled out generally as a significant source of extinction vulnerability in small sexual populations or as a selective force influencing mating-system evolution.     

Mutation accumulation in a selfing population: consequences of different mutation rates between selfers and outcrossers

Currently existing theories predict that because deleterious mutations accumulate at a higher rate, selfing populations suffer from more intense genetic degradation relative to outcrossing populations. This prediction may not always be true when we consider a potential difference in deleterious mutation rate between selfers and outcrossers. By analyzing the evolutionary stability of selfing and outcrossing in an infinite population, we found that the genome-wide deleterious mutation rate would be lower in selfing than in outcrossing organisms. When this difference in mutation rate was included in simulations, we found that in a small population, mutations accumulated more slowly under selfing rather than outcrossing. This result suggests that under frequent and intense bottlenecks, a selfing population may have a lower risk of genetic extinction than an outcrossing population.     

The double edged sword: The demographic consequences of the evolution of self‐fertilization

Phylogenies indicate that the transition from outcrossing to selfing is frequent, with selfing populations being more prone to extinction. The rates of transition to selfing and extinction, acting on different timescales, could explain the observed distributions of extant selfing species among taxa. However, phylogenetic and theoretical studies consider these mechanisms independently, that is transitions do not cause extinction. Here, we theoretically explore the demographic consequences of the evolution of self‐fertilization. Deleterious mutations and mutations modifying the selfing rate are recurrently introduced and the number of offspring depends on individual fitness, allowing for a demographic feedback. We show that mutational meltdowns can be triggered in populations evolving near strict selfing. Populations having survived a demographic crash are more stable than ancestral outcrossing populations once deleterious mutations are purged. The relatively rapid time‐scales at which extinctions occur indicate that during evolutionary transitions the accumulation of deleterious mutations may not be the cause of extinctions observed on longer time scales, but could lead to the underestimation of transition rates from outcrossing to selfing.     

MUTATION LOAD AND THE SURVIVAL OF SMALL POPULATIONS

Previous attempts to model the joint action of selection and mutation in finite populations have treated population size as being independent of the mutation load. However, the accumulation of deleterious mutations is expected to cause a gradual reduction in population size. Consequently, in small populations random genetic drift will progressively overpower selection making it easier to fix future mutations. This synergistic interaction, which we refer to as a mutational melt-down, ultimately leads to population extinction. For many conditions, the coefficient of variation of extinction time is less than 0.1, and for species that reproduce by binary fission, the expected extinction time is quite insensitive to population carrying capacity. These results are consistent with observations that many cultures of ciliated protozoans and vertebrate fibroblasts have characteristic extinction times. The model also predicts that clonal lineages are unlikely to survive more than 10⁴ to 10⁵ generations, which is consistent with existing data on parthenogenetic animals. Contrary to the usual view that Muller's ratchet does more damage when selection is weak, we show that the mean extinction time declines as mutations become more deleterious. Although very small sexual populations, such as self-fertilized lines, are subject to mutational meltdowns, recombination effectively eliminates the process when the effective population size exceeds a dozen or so. The concept of the effective mutation load is developed, and several procedures for estimating it are described. It is shown that this load can be reduced substantially when mutational effects are highly variable.     

MUTATIONAL MELTDOWN IN SELFING ARABIDOPSIS LYRATA

The majority of plant species and many animals are hermaphrodites, with individuals expressing both female and male function. Although hermaphrodites can potentially reproduce by self‐fertilization, they have a high prevalence of outcrossing. The genetic advantages of outcrossing are described by two hypotheses: avoidance of inbreeding depression because selfing leads to immediate expression of recessive deleterious mutations, and release from drift load because self‐fertilization leads to long‐term accumulation of deleterious mutations due to genetic drift and, eventually, to extinction. I tested both hypotheses by experimentally crossing Arabidopsis lyrata plants (self‐pollinated, cross‐pollinated within the population, or cross‐pollinated between populations) and measuring offspring performance over 3 years. There were 18 source populations, each of which was either predominantly outcrossing, mixed mating, or predominantly selfing. Contrary to predictions, outcrossing populations had low inbreeding depression, which equaled that of selfing populations, challenging the central role of inbreeding depression in mating system shifts. However, plants from selfing populations showed the greatest increase in fitness when crossed with plants from other populations, reflecting higher drift load. The results support the hypothesis that extinction by mutational meltdown is why selfing hermaphroditic taxa are rare, despite their frequent appearance over evolutionary time.     

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.     

An explicit model for the inbreeding load in the evolutionary analysis of selfing

I present analytical predictions for the equilibrium inbreeding load expected in a population under mutation, selection, and a regular mating system for any population size and for any magnitude and recessivity of the deleterious effects. Using this prediction, I deduce the relative fitness of mutant alleles with small effect on selfing to explore the situations where selfing or outcrossing are expected to evolve. The results obtained are in agreement with previous literature, showing that natural selection is expected to lead to stable equilibria where populations show either complete outcrossing or complete selfing, and that selfing is promoted by large deleterious mutation rates. I find that the evolution of selfing is favored by a large recessivity of deleterious effects, while the magnitude of homozygous deleterious effects only becomes relevant in relatively small populations. This result contradicts the standard assumption that purging in large populations will only promote selfing when homozygous deleterious effects are large, and implies that previously published results obtained assuming lethal mutations in large populations can be extrapolated to nonlethal alleles of similar recessivity. This conclusion and the general approach used in this analysis can be useful in the study of the evolution of mating systems.     

Testing for epistasis between deleterious mutations in a parasitoid wasp

Abstract.— Determining the way in which deleterious mutations interact to effect fitness is crucial to numerous areas in evolutionary biology. For example, if each additional mutation leads to a greater decrease in log fitness than the last, termed synergistic epistasis, then sex and recombination provide an advantage because they enable deleterious mutations to be eliminated more efficiently. However, there is a severe shortage of relevant empirical data, especially of the form that can help test mutational explanations for the widespread occurrence of sex. Here, we test for epistasis in the parasitic wasp Nasonia vitripennis , examining the fitness consequences of chemically induced deleterious mutations. We examine two components of fitness, both of which are thought to be important in natural populations of parasitic wasps: longevity and egg production. Our results show synergistic epistasis for longevity, but not for egg production.     

Cyclical dynamics under constant selection against mutations in haploid and diploid populations with facultative selfing

We have found that constant selection against mutations can cause cyclical dynamics in a population with facultative selfing. When this happens, the distribution of the number of deleterious mutations per genotype fluctuates with the period approximately 1/sHe generations, where sHe is the coefficient of selection against a heterozygous mutation. The amplitude of oscillations of the mean population fitness often exceeds an order of magnitude. Cyclical dynamics can occur under intermediate selfing rates if selection against heterozygous mutations is weak and selection against homozygous mutations is much stronger. Cycling is possible without epistasis or with diminishing-returns epistasis, but not with synergistic epistasis. Under multiplicative selection, cycling might happen if the haploid mutation rate exceeds 1.9 in the case of selfing of haploids, and if this diploid mutation rate exceeds 4.5 in the case of selfing of diploids. We propose a heuristic explanation for cycling under facultative selfing and discuss its possible relevance.     

Mitosis, stature and evolution of plant mating systems: low-Phi and high-Phi plants

There is a long-recognized association in plants between small stature and selfing, and large stature and outcrossing. Inbreeding depression is central to several hypotheses for this association, but differences in the evolutionary dynamics of inbreeding depression associated with differences in stature are rarely considered. Here, we propose and test the Phi model of plant mating system evolution, which assumes that the per-generation mutation rate of a plant is a function of the number of mitoses (Phi) that occur from zygote to gamete, and predicts fundamental differences between low-Phi (small-statured) and high-Phi (large-statured) plants in the outcomes of the joint evolution of outcrossing rate and inbreeding depression. Using a large dataset of published population genetic studies of angiosperms and conifers, we compute fitted values of inbreeding depression and deleterious mutation rates for small- and large-statured plants. Consistent with our Phi model, we find that populations of small-statured plants exhibit a range of mating systems, significantly lower mutation rates, and intermediate inbreeding depression, while large-statured plants exhibit very high mutation rates and the maximum inbreeding depression of unity. These results indicate that (i) inbred progeny typically observed in large-statured plant populations are completely lost prior to maturity in nearly all populations; (ii) evolutionary shifts from outcrossing to selfing are generally not possible in large-statured species, rather, large-statured species are more likely to evolve mating systems that avoid selfing such as self-incompatibility and dioecy; (iii) destabilization of the mating system-high selfing rate with high-inbreeding depression-might be a common occurrence in large-statured species; and (iv) large-statured species in fragmented populations might be at higher risk of extinction than previously thought. Our results help to unify and simplify a large and diverse field of research, and serve to emphasize the importance that developmental and genetic constraints play in the evolution of plant mating systems.     

Can females benefit from selfing avoidance? Genetic associations and the evolution of plant gender

Evolutionary stability of dioecy and nuclear gynodioecy in higher plants requires that females produce over twice as many successful seeds as hermaphrodites. This fitness differential is widely thought to derive primarily from the advantages of outcrossing caused by high selfing rates and inbreeding depression in the hermaphrodite. This study hypothesized that (i) extraordinarily high deleterious mutation rates are necessary to double female seed success due to outcrossing, and (ii) the large difference in outcrossing rates between sex morphs causes differential purging of these mutations, resulting in additional genetic selection on male sterility. Using genetically explicit models, I showed that the phenotypic outcrossing advantage requires at least one new highly recessive deleterious mutation per genome per generation, regardless of selection coefficient. However, under this mutational regime, differential purging created strong genetic selection against recessive male sterility that overwhelmed the phenotypic selection in favour of outcrossing. In very small populations and for dominant male sterility, this genetic selection was weaker or absent. This first genetically explicit study of the outcrossing advantage of unisexual females may shed new light on both the genetic and selective conditions for the evolution of gynodioecy and dioecy.     

Muller''s Ratchet, Epistasis and Mutation Effects

In this study, computer simulation is used to show that despite synergistic epistasis for fitness, Muller's ratchet can lead to lethal fitness loss in a population of asexuals through the accumulation of deleterious mutations. This result contradicts previous work that indicated that epistasis will halt the ratchet. The present results show that epistasis will not halt the ratchet provided that rather than a single deleterious mutation effect, there is a distribution of deleterious mutation effects with sufficient density near zero. In addition to epistasis and mutation distribution, the ability of Muller's ratchet to lead to the extinction of an asexual population under epistasis for fitness depends strongly on the expected number of offspring that survive to reproductive age. This strong dependence is not present in the nonepistatic model and suggests that interpreting the population growth parameter as fecundity is inadequate. Because a continuous distribution of mutation effects is used in this model, an emphasis is placed on the dynamics of the mutation effect distribution rather than on the dynamics of the number of least mutation loaded individuals. This perspective suggests that current models of gene interaction are too simple to apply directly to long-term prediction for populations undergoing the ratchet.     

THE EVOLUTION OF SELF-FERTILIZATION AND INBREEDING DEPRESSION IN PLANTS. I. GENETIC MODELS

The amounts of inbreeding depression upon selfing and of heterosis upon outcrossing determine the strength of selection on the selfing rate in a population when this evolves polygenically by small steps. Genetic models are constructed which allow inbreeding depression to change with the mean selfing rate in a population by incorporating both mutation to recessive and partially dominant lethal and sublethal alleles at many loci and mutation in quantitative characters under stabilizing selection. The models help to explain observations of high inbreeding depression (> 50%) upon selfing in primarily outcrossing populations, as well as considerable heterosis upon outcrossing in primarily selfing populations. Predominant selfing and predominant outcrossing are found to be alternative stable states of the mating system in most plant populations. Which of these stable states a species approaches depends on the history of its population structure and the magnitude of effect of genes influencing the selfing rate.     

Fixation of new alleles and the extinction of small populations: drift load, beneficial alleles, and sexual selection

With a small effective population size, random genetic drift is more important than selection in determining the fate of new alleles. Small populations therefore accumulate deleterious mutations. Left unchecked, the effect of these fixed alleles is to reduce the reproductive capacity of a species, eventually to the point of extinction. New beneficial mutations, if fixed by selection, can restore some of this lost fitness. This paper derives the overall change in fitness due to fixation of new deleterious and beneficial alleles, as a function of the distribution of effects of new mutations and the effective population size. There is a critical effective size below which a population will on average decline in fitness, but above which beneficial mutations allow the population to persist. With reasonable estimates of the relevant parameters, this critical effective size is likely to be a few hundred. Furthermore, sexual selection can act to reduce the fixation probability of deleterious new mutations and increase the probability of fixing new beneficial mutations. Sexual selection can therefore reduce the risk of extinction of small populations.     

Epistasis,inbreeding depression,and the evolution of self-fertilization

Inbreeding depression resulting from partially recessive deleterious alleles is thought to be the main genetic factor preventing self-fertilizing mutants from spreading in outcrossing hermaphroditic populations. However, deleterious alleles may also generate an advantage to selfers in terms of more efficient purging, while the effects of epistasis among those alleles on inbreeding depression and mating system evolution remain little explored. In this article, we use a general model of selection to disentangle the effects of different forms of epistasis (additive-by-additive, additive-by-dominance, and dominance-by-dominance) on inbreeding depression and on the strength of selection for selfing. Models with fixed epistasis across loci, and models of stabilizing selection acting on quantitative traits (generating distributions of epistasis) are considered as special cases. Besides its effects on inbreeding depression, epistasis may increase the purging advantage associated with selfing (when it is negative on average), while the variance in epistasis favors selfing through the generation of linkage disequilibria that increase mean fitness. Approximations for the strengths of these effects are derived, and compared with individual-based simulation results.     

Estimation of parameters of deleterious mutations in partial selfing or partial outcrossing populations and in nonequilibrium populations

The Deng-Lynch method was developed to estimate the rate and effects of deleterious genomic mutations (DGM) in natural populations under the assumption that populations are either completely outcrossing or completely selfing and that populations are at mutation-selection (M-S) balance. However, in many plant and animal populations, selfing or outcrossing is often incomplete in that a proportion of populations undergo inbreeding while the rest are outcrossing. In addition, the degrees of deviation of populations from M-S balance are often not known. Through computer simulations, we investigated the robustness and the applicability of the Deng-Lynch method under different degrees of partial selfing or partial outcrossing and for nonequilibrium populations approaching M-S balance at different stages. The investigation was implemented under constant, variable, and epistatic mutation effects. We found that, generally, the estimation by the Deng-Lynch method is fairly robust if the selfing rate (S) is <0.10 in outcrossing populations and if S > 0.8 in selfing populations. The estimation may be unbiased under partial selfing with variable and epistatic mutation effects in predominantly outcrossing populations. The estimation is fairly robust in nonequilibrium populations at different stages approaching M-S balance. The dynamics of populations approaching M-S balance under various parameters are also studied. Under mutation and selection, populations approach balance at a rapid pace. Generally, it takes 400-2000 generations to reach M-S balance even when starting from homogeneous individuals free of DGM. Our investigation here provides a basis for characterizing DGM in partial selfing or outcrossing populations and for nonequilibrium populations.     

FIVE GENERATIONS OF ENFORCED SELFING AND OUTCROSSING IN MIMULUS GUTTATUS: INBREEDING DEPRESSION VARIATION AT THE POPULATION AND FAMILY LEVEL

The focus of this study was to examine the consequences of five sequential generations of enforced selfing and outcrossing in two annual populations of the mixed-mating Mimulus guttatus. Our primary goal was to determine whether purging of deleterious recessive alleles occurs uniformly between populations and among families, and thus gain insights into the mode of gene action (dominance, overdominance, and/or epistasis) governing the expression of inbreeding depression at both the population and family levels across the life cycle.     

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.     

Hitchhiking of Deleterious Alleles and the Cost of Adaptation in Partially Selfing Species

Self-fertilization is generally seen to be disadvantageous in the long term. It increases genetic drift, which subsequently reduces polymorphism and the efficiency of selection, which also challenges adaptation. However, high selfing rates can increase the fixation probability of recessive beneficial mutations, but existing theory has generally not accounted for the effect of linked sites. Here, we analyze a model for the fixation probability of deleterious mutants that hitchhike with selective sweeps in diploid, partially selfing populations. Approximate analytical solutions show that, conditional on the sweep not being lost by drift, higher inbreeding rates increase the fixation probability of the deleterious allele, due to the resulting reduction in polymorphism and effective recombination. When extending the analysis to consider a distribution of deleterious alleles, as well as the average fitness increase after a sweep, we find that beneficial alleles generally need to be more recessive than the previously assumed dominance threshold (h < 1/2) for selfing to be beneficial from one-locus theory. Our results highlight that recombination aiding the efficiency of selection on multiple loci amplifies the fitness benefits of outcrossing over selfing, compared to results obtained from one-locus theory. This effect additionally increases the parameter range under which obligate outcrossing is beneficial over partial selfing.     

EFFECTS OF DIFFERENT LEVELS OF INBREEDING ON PROGENY FITNESS IN PLANTAGO CORONOPUS

Inbreeding depression (δ) is a major selective force favoring outcrossing in flowering plants. Many phenotypic and genetic models of the evolution of selfing conclude that complete outcrossing should evolve whenever inbreeding depression is greater than one-half, otherwise selfing should evolve. Recent theoretical work, however, has challenged this view and emphasized (1) the importance of variation in inbreeding depression among individuals within a population; and (2) the nature of gene action between deleterious mutations at different loci (epistasis) as important determinants for the evolution of plant mating systems. The focus of this study was to examine the maintenance of inbreeding depression and the relationship between inbreeding level and inbreeding depression at both the population and the individual level in one population of the partially self-fertilizing plant Plantago coronopus (L.). Maternal plants, randomly selected from an area of about 50 m² in a natural population, were used to establish lines with expected inbreeding coefficients (f) of 0, 0.25, 0.50, 0.75, and 0.875. Inbreeding depression was estimated both in the greenhouse and at the site of origin of the maternal plants by comparing growth, survival, flowering, and seed production of the progeny with different inbreeding coefficients. No significant inbreeding depression for these fitness traits was detected in the greenhouse after 16 weeks. This was in strong contrast to the field, where the traits all displayed significant inbreeding depression and declined with increased inbreeding. The results were consistent with the view that mutation to mildly deleterious alleles is the primary cause of inbreeding depression. At the family level, significantly different maternal line responses (maternal parent × inbreeding level interaction) provide a mechanism for the invasion of a selfing variant into the population through any maternal line exhibiting purging of its genetic load. At the population level, evidence for synergistic epistasis was detected for the probability of flowering, but not for total seed production. At the family level, however, a significant interaction between inbreeding level and maternal families for both traits was observed, indicating that epistasis could play a role in the expression of inbreeding depression among maternal lines.