The degeneration of asexual haploid populations and the speed of Muller's ratchet

The accumulation of deleterious mutations due to the process known as Muller's ratchet can lead to the degeneration of nonrecombining populations. We present an analytical approximation for the rate at which this process is expected to occur in a haploid population. The approximation is based on a diffusion equation and is valid when N exp(-u/s) > 1, where N is the population size, u is the rate at which deleterious mutations occur, and s is the effect of each mutation on fitness. Simulation results are presented to show that the approximation estimates the rate of the process better than previous approximations for values of mutation rates and selection coefficients that are compatible with the biological data. Under certain conditions, the ratchet can turn at a biologically significant rate when the deterministic equilibrium number of individuals free of mutations is substantially >100. The relevance of this process for the degeneration of Y or neo-Y chromosomes is discussed.     

Hemiclonal reproduction slows down the speed of Muller's ratchet in the hybridogenetic frog Rana esculenta

Rare recombination in otherwise asexually reproducing organisms is known to beneficially influence the fitness in small populations. In most of the investigated organisms, asexual and rare sexual generations with recombination follow each other sequentially. Here we present a case where clonal reproduction and rare recombination occur simultaneously in the same population. The hybridogenetic water frog Rana esculenta (E), a hybrid between R. lessonae (L) and R. ridibunda (R) produces gametes that only contain the unaltered maternal R part of their genome. New generations of R. esculenta usually arise from E x L matings. Intraspecific E x E matings produce mostly inviable offspring, but in rare cases, female R. ridibunda arise from such matings which are capable of recombination. In the absence of conspecific males, these R females have to mate with E males, which results in further R females, or with L males, which produces new E lineages. This indirect mechanism reintroduces recombination into the otherwise clonally transmitted R genomes in R. esculenta populations. In this study, we show through Monte Carlo simulations that, in most cases, it is sufficient that only between 1 % and 10 % of mixed water frog populations consist of R females to prevent or significantly reduce the fixation and accumulation of deleterious mutations.     

Muller's ratchet and the accumulation of favourable mutations

Under the influence of recurrent deleterious mutation and selection, asexual and sexual populations reach a deterministic equilibrium with individuals carrying 0,1,2,. . . harmful mutations. When a favourable mutation (aA) occurs in an asexual population it will usually occur in an individual who has one or more (k) deleterious mutations. Muller's ratchet then applies as A will thereafter never occur in an individual with less than k mutations. If the selective advantage of A is less than the selective disadvantage of k harmful mutations then A will not spread. If it is greater it may spread carrying k deleterious mutations to fixation. Sexual populations are not affected in this way. A will spread through the population experiencing genomes with 0,1,2,. . . deleterious mutations in accordance with the deterministic equilibrium.     

The speed of Muller's ratchet with background selection, and the degeneration of Y chromosomes.

The rate of accumulation of deleterious mutations by Muller's ratchet is investigated in large asexual haploid populations, for a range of parameters with potential biological relevance. The rate of this process is studied by considering a very simple model in which mutations can have two types of effect: either strongly deleterious or mildly deleterious. It is shown that the rate of accumulation of mildly deleterious mutations can be greatly increased by the presence of strongly deleterious mutations, and that this can be predicted from the associated reduction in effective population size (the background selection effect). We also examine the rate of the ratchet when there are two classes of mutation of similar but unequal effects on fitness. The accuracy of analytical approximations for the rate of this process is analysed. Its possible role in causing the degeneration of Y and neo-Y chromosomes is discussed in the light of our present knowledge of deleterious mutation rates and selection coefficients.     

Muller's ratchet and the pattern of variation at a neutral locus

The levels and patterns of variation at a neutral locus are analyzed in a haploid asexual population undergoing accumulation of deleterious mutations due to Muller's ratchet. We find that the movement of Muller's ratchet can be associated with a considerable reduction in genetic diversity below classical neutral expectation. The extent to which variability is reduced is a function of the deleterious mutation rate, the fitness effects of the mutations, and the population size. Approximate analytical expressions for the expected genetic diversity are compared with simulation results under two different models of deleterious mutations: a model where all deleterious mutations have equal effects and a model where there are two classes of deleterious mutations. We also find that Muller's ratchet can produce a considerable distortion in the neutral frequency spectrum toward an excess of rare variants.     

A stochastic model for a single click of Muller's ratchet

This work presents a new approach to Muller's ratchet, where Haigh's model is approximately mapped into a simpler model that describes the behaviour of a population after a click of the ratchet, i.e., after loss of what was the fittest class. This new model predicts the distribution of times to the next click of the ratchet and is equivalent to a Wright-Fisher model for a population of haploid asexual individuals with one locus and two alleles. Within this model, the fittest members of a population correspond to carriers of one allele, while all other individuals have suboptimal fitness and are represented as carriers of the other allele. In this way, all suboptimal fitness individuals are amalgamated into a single “mutant” class.The approach presented here has some limitations and the potential for improvement. However, it does lead to results for the rate of the ratchet that, over a wide range of parameters, are accurate within one order of magnitude of simulation results. This contrasts with existing approaches, which are designed for only one or other of the two different parameter regimes known for the ratchet and are more accurate only in the parameter regime they were designed for.Numerical results are presented for the mean time between clicks of the ratchet for (i) the Wright-Fisher model, (ii) a diffusion approximation of this model and (iii) individually based simulations of a full model. The diffusion approximation is validated over a wide range of parameters by its close agreement with the Wright-Fisher model.The present work predicts that: (a) the time between clicks of the ratchet is insensitive to the value of the selection coefficient when the genomic mutation rate is large compared with the selection coefficient against a deleterious mutation, (b) the time interval between clicks of the ratchet has, approximately, an exponential distribution (or its discrete analogue). It is thus possible to determine the variance in times between clicks, given the expected time between clicks. Evidence for both (a) and (b) is seen in simulations.     

Temporal variation in selection accelerates mutational decay by Muller's ratchet

Asexual species accumulate deleterious mutations through an irreversible process known as Muller's ratchet. Attempts to quantify the rate of the ratchet have ignored the role of temporal environmental heterogeneity even though it is common in nature and has the potential to affect overall ratchet rate. Here we examine Muller's ratchet in the context of conditional neutrality (i.e., mutations that are deleterious in some environmental conditions but neutral in others) as well as more subtle changes in the strength (but not sign) of selection. We find that temporal variation increases the rate of the ratchet (mutation accumulation) and the rate of fitness decline over that of populations experiencing constant selection of equivalent average strength. Temporal autocorrelation magnifies the effects of temporal heterogeneity and can allow the ratchet to operate at large population sizes in which it would be halted under constant selection. Classic studies of Muller's ratchet show that the rate of fitness decline is maximized when selection is of a low but intermediate strength. This relationship changes quantitatively with all forms of temporal heterogeneity studied and changes qualitatively when there is temporal autocorrelation in selection. In particular, the rate of fitness decline can increase indefinitely with the strength of selection with some forms of temporal heterogeneity. Our finding that temporal autocorrelation in selection dramatically increases ratchet rate and rate of fitness decline may help to explain the paucity of asexual taxa.     

Muller's ratchet and the degeneration of Y chromosomes: a simulation study

A typical pattern in sex chromosome evolution is that Y chromosomes are small and have lost many of their genes. One mechanism that might explain the degeneration of Y chromosomes is Muller's ratchet, the perpetual stochastic loss of linkage groups carrying the fewest number of deleterious mutations. This process has been investigated theoretically mainly for asexual, haploid populations. Here, I construct a model of a sexual population where deleterious mutations arise on both X and Y chromosomes. Simulation results of this model demonstrate that mutations on the X chromosome can considerably slow down the ratchet. On the other hand, a lower mutation rate in females than in males, background selection, and the emergence of dosage compensation are expected to accelerate the process.     

Absorbing phenomena and escaping time for Muller's ratchet in adaptive landscape

Background

The accumulation of deleterious mutations of a population directly contributes to the fate as to how long the population would exist, a process often described as Muller's ratchet with the absorbing phenomenon. The key to understand this absorbing phenomenon is to characterize the decaying time of the fittest class of the population. Adaptive landscape introduced by Wright, a re-emerging powerful concept in systems biology, is used as a tool to describe biological processes. To our knowledge, the dynamical behaviors for Muller's ratchet over the full parameter regimes are not studied from the point of the adaptive landscape. And the characterization of the absorbing phenomenon is not yet quantitatively obtained without extraneous assumptions as well.

Methods

We describe how Muller's ratchet can be mapped to the classical Wright-Fisher process in both discrete and continuous manners. Furthermore, we construct the adaptive landscape for the system analytically from the general diffusion equation. The constructed adaptive landscape is independent of the existence and normalization of the stationary distribution. We derive the formula of the single click time in finite and infinite potential barrier for all parameters regimes by mean first passage time.

Results

We describe the dynamical behavior of the population exposed to Muller's ratchet in all parameters regimes by adaptive landscape. The adaptive landscape has rich structures such as finite and infinite potential, real and imaginary fixed points. We give the formula about the single click time with finite and infinite potential. And we find the single click time increases with selection rates and population size increasing, decreases with mutation rates increasing. These results provide a new understanding of infinite potential. We analytically demonstrate the adaptive and unadaptive states for the whole parameters regimes. Interesting issues about the parameters regions with the imaginary fixed points is demonstrated. Most importantly, we find that the absorbing phenomenon is characterized by the adaptive landscape and the single click time without any extraneous assumptions. These results suggest a graphical and quantitative framework to study the absorbing phenomenon.

     

Quantifying the genomic decay paradox due to Muller's ratchet in human mitochondrial DNA

The observation of high mitochondrial mutation rates in human pedigrees has led to the question of how such an asexual genetic system can survive the accumulation of slightly deleterious mutations caused by Muller's ratchet. I define a null model to quantify in unprecedented detail the threat from extinction caused by Muller's ratchet. This model is general enough to explore the biological significance of Muller's ratchet in various species where its operation has been suspected. For increased precision over a wide range of parameter space I employ individual-based simulations run by evolution@home, the first global computing system for evolutionary biology. After compiling realistic values for the key parameters in human mitochondrial DNA (mtDNA) I find that a surprisingly large range of biologically realistic parameter combinations would lead to the extinction of the human line over a period of 20 million years - if accepted wisdom about mtDNA and Muller's ratchet is correct. The resulting genomic decay paradox complements a similar threat from extinction due to mutation accumulation in nuclear DNA and suggests evaluation of unconventional explanations for long-term persistence. A substantial list of potential solutions is given, including compensatory back mutations, mutation rate heterogeneity and occasional recombination in mtDNA. Future work will have to explore which of these actually solves the paradox. Nonetheless, the results presented here provide yet another reason to minimize anthropogenic increase of mutation rates.     

Many-trillionfold amplification of single RNA virus particles fails to overcome the Muller's ratchet effect.

We showed earlier that transfers of large populations of RNA viruses lead to fitness gains and that repeated genetic bottleneck transfers result in fitness losses due to Muller's ratchet. In the present study, we examined the effects of genetic bottleneck passages intervening between population passages, a process akin to some natural viral transmissions, using vesicular stomatitis virus as a model. Our findings show that the pronounced fitness increases that occur during two successive population passages cannot overcome the fitness decreases caused by a single intervening genetic bottleneck passage. The implications for natural transmissions of RNA viruses are discussed.     

Mutation accumulation in transfer RNAs: molecular evidence for Muller's ratchet in mitochondrial genomes

The accumulation of deleterious mutations is thought to be a major factor preventing the long-term persistence of obligately asexual lineages relative to their sexual ancestors. This phenomenon is also of potential relevance to sexual species that harbor asexually propagating organelle genomes. A comparative study of the transfer RNA genes in animal mitochondrial and nuclear genomes demonstrates that the former accumulate nucleotide substitutions much more rapidly than do the latter, and several lines of evidence are consistent with the idea that the excess substitutions are mildly deleterious. First, the average binding stability between complementary strands in the stems of mitochondrial tRNAs is less than half that in nuclear tRNAs. Second, most loop sizes in the mitochondrial tRNAs have experienced a net reduction in size over evolutionary time, and they are nearly 50 times more variable in the mitochondrial than in the nuclear genome. Third, although nearly 20% of the nucleotides in nuclear tRNA genes (particularly those involved in tertiary interactions) are invariant across all animal taxa and all tRNA species, there are no invariant sites in the mitochondrial tRNAs. These observations, as well as results from recent laboratory experiments, are consistent with the hypothesis that nonrecombining organelle genomes are subject to gradual loss of fitness due to the cumulative chance fixation of mildly deleterious mutations.      

Quantifying the threat of extinction from Muller's ratchet in the diploid Amazon molly (<Emphasis Type="Italic">Poecilia formosa</Emphasis>)

Background  

The Amazon molly (Poecilia formosa) is a small unisexual fish that has been suspected of being threatened by extinction from the stochastic accumulation of slightly deleterious mutations that is caused by Muller's ratchet in non-recombining populations. However, no detailed quantification of the extent of this threat is available.     

Fluctuations of Fitness Distributions and the Rate of Muller's Ratchet

The accumulation of deleterious mutations is driven by rare fluctuations that lead to the loss of all mutation free individuals, a process known as Muller's ratchet. Even though Muller's ratchet is a paradigmatic process in population genetics, a quantitative understanding of its rate is still lacking. The difficulty lies in the nontrivial nature of fluctuations in the fitness distribution, which control the rate of extinction of the fittest genotype. We address this problem using the simple but classic model of mutation selection balance with deleterious mutations all having the same effect on fitness. We show analytically how fluctuations among the fittest individuals propagate to individuals of lower fitness and have dramatically amplified effects on the bulk of the population at a later time. If a reduction in the size of the fittest class reduces the mean fitness only after a delay, selection opposing this reduction is also delayed. This delayed restoring force speeds up Muller's ratchet. We show how the delayed response can be accounted for using a path-integral formulation of the stochastic dynamics and provide an expression for the rate of the ratchet that is accurate across a broad range of parameters.