Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Quantifying the threat of extinction from Muller's ratchet in the diploid Amazon molly (Poecilia formosa)

Laurence Loewe12* and Dunja K Lamatsch3456

Author affiliations

1 Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Ashworth Laboratories, King's Buildings, Edinburgh EH9 3JT, UK

2 Centre for Systems Biology Edinburgh, School of Biological Sciences, University of Edinburgh, Darwin Building, King's Buildings, Edinburgh EH9 3JU, UK

3 Universität Würzburg, Institute of Physiological Chemistry I, Biocenter, Würzburg, 97074 Würzburg, Germany

4 Freshwater Biology, Royal Belgian Institute of Natural Sciences, Vautierstraat 29, B – 1000 Brussels, Belgium

5 University of Sheffield, Department of Animal and Plant Sciences, Alfred Denny Building, Western Bank, Sheffield, S10 2TN, UK

6 Austrian Academy of Sciences, Institute for Limnology, Mondseestrasse 9, 5310 Mondsee, Austria

For all author emails, please log on.

Citation and License

BMC Evolutionary Biology 2008, 8:88  doi:10.1186/1471-2148-8-88


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/8/88


Received:28 June 2007
Accepted:19 March 2008
Published:19 March 2008

© 2008 Loewe and Lamatsch; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

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.

Results

Here we quantify genomic decay in this fish by using a simple model of Muller's ratchet with the most realistic parameter combinations available employing the evolution@home global computing system. We also describe simple extensions of the standard model of Muller's ratchet that allow us to deal with selfing diploids, triploids and mitotic recombination. We show that Muller's ratchet creates a threat of extinction for the Amazon molly for many biologically realistic parameter combinations. In most cases, extinction is expected to occur within a time frame that is less than previous estimates of the age of the species, leading to a genomic decay paradox.

Conclusion

How then does the Amazon molly survive? Several biological processes could individually or in combination solve this genomic decay paradox, including paternal leakage of undamaged DNA from sexual sister species, compensatory mutations and many others. More research is needed to quantify the contribution of these potential solutions towards the survival of the Amazon molly and other (ancient) asexual species.

Background

It is a general observation that asexual lineages do not last over very long periods of time, but the precise reasons for this are less clear [1]. Ancient asexuals are the rare exceptions to this rule and it is of considerable biological interest to know what mechanisms allow them to survive for so long. They must have found a way to overcome the long-term fitness-degrading consequences of genetic processes like Muller's ratchet and/or ecological processes like Red-Queen dynamics [1-6]. Muller's ratchet describes the long-term accumulation of slightly deleterious mutations in asexual populations and has been suggested as a key mechanism for the extinction of asexual species on the long term [3,7,8]. Unfortunately most of these predictions remain at a stage of verbal argument [1], making it very difficult to rule out that Muller's ratchet may not have had enough time to cause extinction. As the behavior of Muller's ratchet can be very sensitive to model parameters [7], realistic values need to be used to predict the consequences of mutation accumulation for a given ancient asexual species. We employ a simple null hypothesis [see [7]] for testing the threat of extinction from Muller's ratchet in the unisexual fish Poecilia formosa, the Amazon molly.

Overcoming the lack of quantification

Quantifications of the threat of extinction from Muller's ratchet are often not trivial theoretical work that requires either challenging mathematics [9] or complex computer simulations [7] or both [7]. Therefore some adopt the pragmatic approach that any system with no recombination and a potential for appreciable slightly deleterious mutation rates could be driven to extinction by Muller's ratchet [3]. Such arguments frequently overlook the fact that the particular combination of parameters in that species might not be expected to lead to extinction within the known time of its existence, even if Muller's ratchet is clearly operating [1,3,7]. Thus statements about the evolutionary short lives of asexuals are often less quantitative than would be desirable [1]. In other words there may not be a genomic decay paradox that calls for any special solutions [7]. We advocate the use of a simple model for predicting extinction times caused by Muller's ratchet in order to make current discussions about ancient asexuals more quantitative [1,7], even if that model cannot capture the full complexity of our study species and therefore only leads to tentative predictions. We believe that small steps in model development will allow future models to benefit from the experiences with simpler models. Hence our use of a simple model of Muller's ratchet that ignores all complications like potential Red Queen dynamics that might accelerate the rate of Muller's ratchet (see discussion below and [6]). We focus on testing the null hypothesis that Muller's ratchet could not have led to extinction in a given time frame, as described by Loewe [7]. While this is an important advance over the purely verbal stage, we want to encourage future work to model the various processes that increase or decrease the predicted speed of genomic decay. We also want to encourage more empirical work to establish the precise values of parameters in these models. Such work is needed for other asexual species as well.

The Amazon molly and Muller's ratchet

The fish Poecilia formosa was the first unisexual vertebrate that was discovered [10]. This all-female species resulted from hybridization between relatives of Poecilia mexicana and Poecilia latipinna [11,12] that probably happened between 40,000 and 100,000 years ago (see section 'Age ...' below and [11,12]). Reproduction normally occurs by sperm-dependent parthenogenesis, i.e. diploid eggs are produced, which need to be activated for embryonic development by sperm of closely related species. It has been argued that paternal leakage, leading to the expression of paternal genes, plays a pivotal role to stop Muller's ratchet [13] that otherwise would have driven the species to extinction in less time than its current estimates of existence [11]. Paternal leakage and other processes that may slow down genomic decay are discussed below. We want to determine if these processes are necessary to explain the survival of this fish into our times.

Habitat and population structure

The Amazon molly is a small fish (3 – 7 cm) that lives in a rather limited range from the Nueces River in south-east Texas southward to the mouth of the Rio Tuxpan, north of the Sierra del Abra in Mexico. All these river systems flow from west to east and have no connection other than the sea. The population on such a large scale may have some structure, as populations from south Texas, for example, have no reasonable connection with those in the Río Purificación. However, as the Amazon molly tolerates marine conditions [14], migration cannot be entirely excluded. A study of FST in subpopulations that span a distance of about 100 km in the same river system did not find significant population subdivision [15]. Some simple models of population subdivision do not affect the effective population size Ne and probabilities of the fixation [such as some island models [16,17]]. However, more realistic models of population structure that allow for extinctions and recolonizations can have a substantial impact on deleterious mutation accumulation [18,19]. To simplify our theoretical treatment, we will assume that the whole species has no substructure that is not already accounted for by our assumed Ne.

Genetics

The Amazon molly reproduces gynogenetically, i.e. its eggs contain an unreduced set of chromosomes, that need the sperm of one of the sister species Poecilia mexicana or Poecilia latipinna as a mechanical trigger to start development [20-22]. Usually, only the diploid set of maternal genes is expressed and the paternal genome is expelled. However, occasionally, the paternal genome remains, giving rise to a triploid clone, that reproduces as a triploid gynogen [15,23]. In other cases, only traces of the paternal genome (so-called micro-chromosomes) escape the enzymatic machinery that clears the egg from the nucleus that arrived with the sperm [24]. In all cases, the full diploid or triploid set of chromosomes (with or without micro-chromosomes [24,25]) is clonally passed on to offspring, without an obvious opportunity for recombination [26]. Based on this lack of recombination it was hypothesized that Muller's ratchet should have driven the Amazon molly to extinction within the presumed time of its existence, unless processes like paternal leakage would stop genomic decay [13].

Aims

Here we aim to quantify the verbally predicted effects of Muller's ratchet in order to see, whether there really is a genomic decay paradox as defined by Loewe [7] that calls for an explanation. Results show that indeed a range of realistic parameter combinations should have led to the extinction of Amazon molly within the time of its presumed existence.

Results

We quantified the rate of mutation accumulation due to Muller's ratchet using the best available analytical approximations [9,27,28] and globally distributed individual-based simulations run by Simulator005r6 of evolution@home [7,29-31] assuming our estimates of the most realistic parameter combinations. To this end we used the standard null model of Muller's ratchet described elsewhere [7] and extended it to accommodate the slowdown in fitness decay that can be caused by polyploidy and mitotic recombination. We employ the U-shaped plot of extinction times against selection coefficients to allow easy visualization of situations that lead to a genomic decay paradox. The frequency of these situations can be measured by specifying the range of critical selection coefficients sc that are defined by the prediction of corresponding extinction times that are below Tage, the presumed age of asexuality in an evolutionary line. For more detailed explanations of this plot see Loewe [7].

The results show that values for Usdm, the slightly deleterious mutation rate, that are above Usdm ≈ 0.1, lead to the extinction of the Amazon molly within the estimated Tage = 81,000 years of its existence, even if lower and upper limits for Tage are considered (see Figure 1 and Tage estimates below). This is also true for our best estimate of Usdm in the unlikely case of extremely high levels of mitotic recombination (then Usdm is scaled to 0.2 deleterious mutations with critical effects/diploid genome/generation). These findings are rather independent of the effective population size Ne, as even Ne = 107 (certainly larger than the true Ne of the Amazon molly) will not help against Muller's ratchet if mutation rates are too high.

thumbnailFigure 1. Predicted extinction times of the Amazon molly. Muller's ratchet might cause extinctions of the unisexual Amazon molly due to the accumulation of slightly deleterious mutations. (A) analytical results only, (B) analytical results and simulations combined. The upper bar denotes the assumed age of the line (70 Ky, min 40 Ky, max 100 Ky). The lower bar marks the border to neutrality for the effective population sizes (Ne) used by spanning the selection coefficients from Nes = 1 for the largest (107 females) to the smallest conceivable Ne (104 females). The lines represent the analytic predictions of the extinction time for different deleterious genomic mutation rates (Usdm) with Ne = 316,000, generation time Tgen = 1 year and maximal reproductive capacity Rmax = 500 offspring/generation. Our current best overall estimate for Usdm is 0.4 critically deleterious mutations/diploid genome/generation. Our upper limit is close to Usdm = 1. The dashed lines indicate the variability of the extinction time estimates for a value that is close to our lowest credible mutation rate estimate (Usdm = 0.05, green diamantes) using the corresponding upper and lower limits of Ne, Tgen and Rmax; variability in extinction time is similar for other Usdm. Large symbols denote valid extinction time estimates from simulations with at least 2 observed clicks of Muller's ratchet (usually many more, up to 500). Small symbols denote lower limits for extinction times from simulations without observed clicks, based on the (usually wrong) assumption that the ratchet would have clicked just after stopping the simulation. Each symbol denotes an independent simulation with a different random seed and assumes the same mean Tgen and Rmax as analytic predictions (different Ne have been plotted on top of each other to avoid a series of similar plots). This plot contains 24,251 simulations with a total of 14.78 years of computing time. See main text for a discussion of parameter combinations and Loewe [7] for an in-depth explanation of the U-shaped plot shown here. The location of the wall of background selection for a particular parameter combination is approximately given by the vertical part of the corresponding line: all mutations with effects larger than the location of this wall are removed deterministically.

A deleterious mutation rate of about Usdm ≈ 0.05 per generation results in extinction times that border with upper limits of estimates of the age of the Amazon molly of 105 years. At this mutation rate some values of Ne, Rmax and Tgen permit persistence and others cause extinction within this timeframe, where Rmax is the maximal reproductive capacity of the non-degenerated ancestors and Tgen is the generation time. The variability due to Ne, Rmax and Tgen is indicated for Usdm = 0.05 by the dashed black lines in Figure 1. These factors lead to similar variability of extinction times for other mutation rates. If the deleterious genomic mutation rate were only 0.01, then extinction due to Muller's ratchet in the known age of the P. formosa lineage could be excluded, even with the current uncertainty in other parameters.

If deleterious mutation rates are greater than our best estimate of Usdm ≈ 0.4, then Muller's ratchet will cause extinction of the Amazon molly in a few thousand years from the origin of asexuality under the standard model in the absence of mitigating forces. Such high mutation rates seem to be supported by our best approximation of the genetic architecture in the Amazon molly (many approximately codominant mutations of small effects in a genome with very low levels of mitotic recombination lead to an effective doubling of Usdm as shown in the Equal-Contribution-Model in Table 1 and discussed in the Models section below).

Table 1. Useful transformations for computing expectations and bounds for the rate of Muller's ratchet in diploids.

Following Loewe [7], the threat of extinction from Muller's ratchet in Amazon molly can be quantified in detail as follows:

Usdm = 0.5 (and a mean of all other values) yields a minimal extinction time, Tex, of Tex ≈ 5 Kyr due to most critical selection coefficients, scm, in the range from scm ≈ 0.005 to 0.05. A genomic decay paradox for this mutation rate exists at Tage = 100 Kyr in the range of critical selection coefficients sc ≈ 0.0002 to 0.06.

Usdm = 0.1 (and a mean of all other values) yields a minimal extinction time of Tex ≈ 50 Kyr due to most critical selection coefficients in the range from scm ≈ 0.005 to 0.01. A genomic decay paradox for this mutation rate exists at Tage = 100 Kyr in the range of critical selection coefficients sc ≈ 0.001 to 0.01.

Lower mutation rates under the same conditions lead to minimal extinction times that are longer than the assumed age of the line (Usdm = 0.05 leads to extinction in Tex ≈ 200 Kyr at scm ≈ 0.005 and Usdm = 0.01 to Tex ≈ 3 Myr at scm ≈ 0.001). See the coloured lines in Figure 1 for a visual overview.

If mutation rates are increased beyond the natural levels that we estimated (e.g. by mutagenic pollution), then the most damaging mutational effects are all in the range of several percent or more and resulting minimal extinction times can be surprisingly short. For example, Usdm = 1 can lead to extinction times of about 300 years by accumulating deleterious mutational effects of s ≈ 10%. If mutagenic pollution leads to a further 10-fold increase of Usdm, then extinction times are expected to be less than 30 years, as increasingly harmful mutations start to accumulate as well (extrapolation from Figure 1).

Discussion

This is the first detailed assessment of the threat of extinction from deleterious mutation accumulation through Muller's ratchet in the Amazon molly. Combining our best estimates of the haploid deleterious mutation rate (Usdm ≈ 0.2) with our best approximation of the genetic architecture in the Amazon molly (many approximately codominant mutations of small effects in a genome with very low levels of mitotic recombination) leads to an effective deleterious mutation rate of Usdm ≈ 0.4 mutations with critical effects/diploid genome/generation. In the absence of significant mitigating forces this would cause extinction of the Amazon molly in a few thousand years. More precise results can be taken from Figure 1 if needed. This represents a big step forward for understanding the asexuality in the Amazon molly, since never before have the times to extinction been quantified in such detail for this fish.

Genomic decay paradox

Combining predicted extinction times with current estimates of the age of asexuality in the Amazon molly leads to a genomic decay paradox, as defined by Loewe [7]. Thus mechanisms that extend our standard model of Muller's ratchet are needed to explain why this fish has thus far escaped extinction. Such mechanisms are not needed if either our most plausible mutation rate estimates are too high or the Amazon molly is younger than current age estimates suggest (effective Usdm = 0.04 – 0.9; Tage = 40,000 – 100,000 years, see section 'Age...' below and [11,12]). As current evidence seems hardly compatible with such low mutation rates or such a young age of P. formosa, the quest for mechanisms that help this fish escape genomic decay seems to be valid. This is corroborated by our observation of a fitness correlate, the number of embryos found in females, which does not show significant differences between P. formosa and P. mexicana (one of its parental species, see section on life history below). While we have no such information for other correlates such as longevity, number of broods per life, offspring survival, etc., our current limited evidence suggests that P. formosa experiences only little or no genomic decay despite our predictions of substantial deleterious mutation accumulation. Below we will discuss processes that might be of special importance for maintaining fitness in the Amazon molly.

Mitotic recombination

Some reports of apomixis in the Amazon molly suggest a little debated potential solution to the mystery of its long-term survival. Rasch et al. [26] reported low, but consistent, levels of tissue graft rejections after prolonged periods (up to one year) within certain sibships, suggesting that not all inheritance is strictly isogenic in the Amazon molly. In the absence of a meiotic prophase these observations have been interpreted as the result of either a mutation rate that exceeds expectations or as the result of somatic cell crossing-over [26,32]. This process is also known as mitotic recombination and is most likely an inevitable result of the way that cells organize mitosis. It was first discovered in Drosophila [33,34] and has been intensively studied in yeast [35-37], mice [38], humans [39] and Daphnia [40]. The frequency of mitotic recombination in some fungi was found to be between 102 and 104 times less than that of meiotic recombination [41] and recent experiments in yeast reported a factor of ≈ 25,000 [37]. Estimates suggest that about 10 sister chromatid exchanges per cellcycle can occur in mammalian cells [42]. Such evidence suggests that the Amazon molly experiences mitotic recombination, even if the effective rate of segregation of different mutations is probably much lower than in selfers with meiotic recombination. The fact that asexual Daphnia have been shown to experience mitotic recombination [40] means that the Amazon molly would not be the only asexual to experience mitotic recombination.

If this is true, the resulting segregation might slow down Muller's ratchet for some selection coefficients [43-45], compared to expectations from the Equal-Contributions-Model described below. However, the fact that the distribution of mutational effects is expected to be very wide on a log scale [46,47] means that corresponding smaller selection coefficients will continue to drive Muller's ratchet. Thus it is difficult to see how mitotic recombination could stop the ratchet on its own without the contribution of other processes. We can use a simple model to put an upper limit on the maximal potential of mitotic recombination to stop Muller's ratchet. Mitotic recombination cannot possibly do more to stop fitness decay than in the case of completely free mitotic recombination. This allows for mutations to segregate at a maximal rate and has been used to model Muller's ratchet in selfers [45]. As shown below (see section on Muller's ratchet with selfing), a simple scaling of mutation rate and selection coefficient is enough to extend the standard model of Muller's ratchet to this case. Applying such a scaling to our results does not remove the genomic decay paradox that we find (see Table 1; it is questionable whether mitotic recombination will have a large effect, since the rates of mitotic recombination are probably far from free recombination). Similarly an analysis of the various levels of dominance that are possible for diploids and triploids as shown in Table 1 does not remove the paradox. Thus we will have to search for other solutions.

Rare recombination events during oogenesis in the Amazon molly are not conceptually different from mitotic recombination in the germ line or selfing. To limit the largest possible contribution of these processes towards stopping genomic decay we assume the most generous form of selfing, which is free recombination between both diploid copies. As our analysis above shows, the largest possible amount of recombination during oogenesis cannot stop Muller's ratchet – if no fresh genes are introduced by outcrossing of some sort.

Paternal leakage

Occasional paternal leakage of fresh genetic material from sexual sister species could serve as a template for DNA repair [48] or restore genes that had been destroyed by Muller's ratchet [13,49]. Paternal leakage has been demonstrated to lead to the expression of paternal genes [13,50], suggesting a compensatory role in the Amazon molly [13]. There are two mechanistic scenarios that may facilitate this.

(i) Micro-chromosomes with a size of about 1% of the genome have been observed to leak from a sexual sister species to the Amazon molly [13,24]. If they carry an occasional random sample of genes from the non-degraded sexual genome into the degrading Amazon molly population, then the possibility exists that they might restore an ancient, non-degraded state of fitness. This could slow down Muller's ratchet enough to solve the genomic decay paradox [13,25]. Micro-chromosomes have been demonstrated to be stably inherited over many generations [25,51].

(ii) The finding of triploid clones might encourage the speculation that the third genome copy that is contributed by the sexual sister species might help restore fitness [49]. One might speculate that triploids should produce more offspring or survive better and ultimately substitute the diploid individuals. However, observations suggest the contrary, as triploids are rather limited in their range, young, and much less frequently produced than in other asexuals [15,23]. An alternative possibility is that triploids might occasionally lose one of their three genome copies and give rise to secondary diploids that can carry presumably fitter genes from their sexual sister species. If such a gene flow exists from the sexual sister species over paternal leakage and an intermediate triploid stage to a final diploid stage, then such a flow might contribute towards solving a genomic decay paradox [49]. However, there is no evidence for such a directed gene flow and triploids are produced much rarer in the Amazon molly than in other asexuals or than micro-chromosomes in P. formosa [23].

In any case, paternal leakage should not be confused with true recombination, which has stimulated discussions about the ratchet stopping potential of paternal leakage [52,53]. It has also been speculated that the paternal genome might be used as a template for DNA repair, but its precise role remains unclear [48].

Other processes

There is a long list of other potential solutions for the genomic decay paradox that has been given elsewhere [7] and thus shall not be discussed in detail here. This list includes the unlikely possibility that the true deleterious mutation rate might be much lower, either because mutation rates in the Amazon molly happen to be generally much lower for yet unknown reasons, or because the distribution of mutational effects happens to be strongly bimodal with almost no mutations in the critical intermediate range. These hypotheses are not well supported by comparative analyses of mutation rates and effects in different species (see Parameter Estimates below). Synergistic epistatic effects have also been argued to have the capacity to stop genomic decay [54,55], but this is only true in combination with very specific distributions of mutational effects. If these distributions are reasonably wide, the potential for epistatic effects to decelerate Muller's ratchet is virtually non-existent [56].

We have also ignored advantageous and compensatory mutations, which have a substantial potential to stop genomic decay completely, if they are frequent enough [57-59]. Recent work has suggested that a substantial fraction of all mutations is advantageous [47,60]. If the underlying patterns are not caused by other processes and the selection coefficients of the corresponding mutations are large enough, then advantageous mutations can stop Muller's ratchet. Recent work in viruses has suggested that the ratio of beneficial to deleterious mutations can increase as the mean fitness decreases [61]. If similar dynamics hold for fish then Muller's ratchet may operate much slower, if at all. A more detailed discussion of back mutations and compensatory mutations can be found elsewhere ([7] and references therein).

To increase the precision of extinction time estimates it would be desirable to have more direct empirical estimates of mutation rates and effects in the Amazon molly. Such more precise estimates are needed when specific potential solutions for this genomic decay paradox are to be tested. The evidence presented here makes it seem unlikely that such added information will change our main conclusion that a genomic decay paradox exists for the current age estimates. For more details and additional potential solutions, please see [7].

Red Queen

It is possible that Muller's ratchet is not the only process that leads to genomic decay and that the speed of Muller's ratchet may be significantly increased by other processes. The Amazon molly and its closely related sister species are known to harbor parasites [62,63]. These parasites probably decrease the fitness of their host substantially in the wild and co-evolve with it in an evolutionary arms race for survival. This scenario is described by the Red Queen hypothesis [3,5,64,65] and may in itself lead to the extinction of a species. A Red Queen scenario can also be caused by antagonistic co-evolution in general, which may occur in many circumstances, including evolving predator-prey or plant-herbivore relationships or intra-specific co-evolution. It is also known that the speed of mutation accumulation caused by Muller's ratchet is enhanced in a population that experiences Red Queen dynamics [6,66]. Empirical support for a prerequisite of the Red Queen hypothesis could be found in another member of the family Poeciliidae to which the Amazon molly belongs [67]. Thus the Amazon molly might participate in such an arms race, however we currently do not know the range of biologically realistic parameter combinations that are needed to quantify Red Queen dynamics here. The Red Queen Hypothesis predicts an increased load from parasites in asexuals, because they cannot adapt as fast as sexuals to newly evolving parasites. Based on this prediction one might not expect substantial Red Queen dynamics in the Amazon molly, since it seems to have about as many parasites as its sexual sister species [68]. Another attempt to discover Red Queen dynamics in the Amazon molly was also negative [69,70]. If the Amazon Molly is forced to constantly evolve as under Red Queen dynamics, then the genomic decay paradox might be more extreme.

Implications for the origin of the Amazon molly

Many attempts to produce fertile asexual hybrids in the laboratory have been unsuccessful [71-73], but under natural conditions such hybridization attempts between the sympatric parental species of the Amazon molly, P. mexicana and P. latipinna, might happen occasionally. If such hybridization events occurred regularly, they could have led to a stable existence of the Amazon molly form, even though all individual hybrids are on their way to extinction as they will soon be replaced by fresh hybrids. In this case we expect multiple different young hybrid lines in random samples from the overall population of the Amazon molly, as each lineage of hybrids will be closer to its parental species than to other independent hybrids. The corresponding phylogenetic tree is expected to be polyphyletic.

Such a scenario is not supported by existing mtDNA data [12,72]. Amazon molly individuals sampled from a wide range of locations show a paraphyletic tree that is compatible with a single hybridization event in the distant past. This is based on the observation that Amazon molly mtDNA sequences cluster either with each other or with exactly the same ancestral sequence from the parental species [12]. If all sequences sampled in the future follow the same pattern, this rules out polyphyly and casts serious doubts over the hypothesis of multiple hybridizations in the past [12]. Combining this with the difficulties to produce asexual hybrids in the lab [71-73] suggests that the Amazon molly probably comes from a singularly rare hybridization success. This is in marked difference to some other species, where asexual hybrids are much easier to obtain [73]. Thus our quantification of Muller's ratchet only applies to the descendents of this one clone that appears to have an age of about 80 Ky.

We cannot distinguish whether the Amazon molly originated by a single ancient hybridization event or whether there was a small series of such events involving very similar parental individuals in the distant past. This is of no importance for our conclusions, as the subsequent course of mutation accumulation is not expected to be different from that in a population that goes back to a single hybridization event.

In order for other hybridization events to affect our conclusions, there would have to be a repeated production of fresh clones that could then constantly replace the decaying genomes in the population. This is expected to show in phylogenetic analyses. Current analyses show the absence of polyphyly and indicate that there is no ongoing hybridization in the wild [12,72]. Thus our analysis is appropriate for the time that came after the population was cut off from geneflow. This is the time Tagethat we estimate as the time that Muller's ratchet had for degrading the Amazon molly. In the case that future work uncovers multiple hybridization events, our analysis is valid for the lifespan of each clone, assuming no additional factors like interclonal competition that is independent of Muller's ratchet.

As current genetic data cannot infer what happened before the origin of our current lineage of the Amazon molly, this lineage might well be only the last unisexual Amazon molly that did not yet go extinct, with many others preceding it. As we know nothing about potential previous clones our ability to test the so-called "Frozen-Niche-Variation model" [73,74] is rather limited in this case.

Conservation genetics

Since extinction time is very sensitive to changes in the mutation rate, it is conceivable that the anthropogenic release of mutagenic substances could lead to such a strong increase in mutation rate that extinction times are predicted to be in time frames that are frequently considered by conservation biologists. One of us (DKL) could actually observe a typical Poeciliid habitat that had lost almost all vertebrate life due to apparent water pollution (Altamira, Tamaulipas, Mexico, 2002). Although we have no evidence to decide whether this incident was mutagenic or not, pollution in general often has mutagenic side effects. This suggests the possibility of considerable pollution at least of parts of the habitat of the Amazon molly and would not be the first instance where pollution in rivers leads to a several-fold increase of mutation rates (see observations in ferns [75-78]).

Quantifying Muller's ratchet in other ancient asexuals is easier now

It is well known that asexual lineages are typically short-lived. While this observation is central to theories about the origin of sex, it still needs to be properly quantified [1]. The evolution@home results used in this work make it manageable to quantify the threat from Muller's ratchet for the long list of other putatively ancient asexuals that have been suspected of being threatened by genomic decay [1,4,79,80], including the Amazon molly's close relative Poeciliopsis [74,81]. We suspect that it might not be possible to explain the existence of all these asexuals by the surprising discovery of a recent ancestor with a very young age as in the case of the clonal, hybrid, gynogenetic mole salamander Ambystoma [82-84]. Even if that were the case, then it would still be interesting to quantify Muller's ratchet, as this would shed more light on how fast it actually clicks in clonal lineages, where some evidence is consistent with its operation [85-89]. Examples of ancient asexuals that could benefit from a more rigorous quantification of the effects of Muller's ratchet include:

(i) Darwinula stevensoni, a small non-marine ostracod. Darwinulidae are believed to have lived for about 200 Myr without sex as the fossil record shows only females [80] – but see Smith et al. [90]. The species Darwinula stevensoni is a member of this group and is thought to exist for more than 20 Myr now [80,91,92].

(ii) Bdelloid rotifers. The Class Bdelloidea of the Phylum Rotifera is the largest taxonomic group that has apparently lived completely without sex for at least 40 Myr [93-95]. The ancient asexuality of these 0.1 to 1 mm long animals appears to be as well established as it can possibly be. Various special features of this group have been discussed as the reason behind its long-term survival [96-99].

(iii) Oribatid mites. Parthenogenetic automicts are frequent among oribatid mites [100,101]. Current estimates suggest that some lines have lived parthenogenetically for perhaps 100 Myr and that their extant distribution was strongly affected by continental drift [102].

Practical aspects

To quantify the possibilities of extinction from Muller's ratchet in a given asexual system using the evolution@home results database, please contact one of us (LL). You will receive help in completing a survey of various details about your study system (see online questionnaire [103]) and if enough data is available, a preliminary report will be produced (including plots similar to Figure 1). If necessary, the existing evolution@home infrastructure can be used to compute new parameter combinations. Our experience shows that the prediction of extinction times in ancient asexuals can often be simplified by using rough upper and lower limits for important parameters. For some parameters the lack of precision will not be critical, where as for others it will help focus further empirical work towards parameter estimation.

Why quantify the ratchet as often as possible?

There are many open questions that can be asked about the general accuracy of the simple standard model of Muller's ratchet used for the quantifications presented here [7]. We believe that the following reasons justify a series of analyses of extinction times in various asexuals based on our null-model. We expect such work to contribute towards a mature, quantitative discussion about the evolutionary biology of asexuals.

(i) Critical predictions of extinction time help biologists to look for the key data that is also needed for more realistic quantifications of the effects of mutation accumulation. Therefore such predictions help in the design of empirical work.

(ii) It can be expected that at least some of the putative ancient asexuals are not examples of the genomic decay paradox, because the known age of their asexuality is smaller than their predicted extinction time. This does not deny the age of their asexuality, it just takes the 'scandal' out of the observations [4].

(iii) Those species with an apparent genomic decay paradox can be subject to a more detailed search for mechanisms that solve the paradox and help them to avoid extinction [7].

(iv) Experiences with the present simple system for predicting extinction times can be expected to lead to the development of more realistic systems for the quantification of genomic decay paradoxes. Such improved systems might include the processes discussed here and might measure their mitigating effects on mutation accumulation.

Conclusion

A genomic decay paradox is predicted by a large number of biologically realistic parameter combinations for the unisexual Amazon molly. This is based on a simple model of Muller's ratchet that accounts for the distribution of mutational effects on fitness, the availability of multiple copies of the genome and mitotic recombination. Our prediction of a genomic decay paradox strengthens the conclusions of earlier work that suggested the existence of additional biological processes that slow down or halt the mutational decay of fitness in this fish. This conclusion is consistent with our observation that the Amazon molly carries approximately the same number of embryos per adult as its sexual sister species. Future work will have to establish whether paternal leakage of micro-chromosomes or still other processes have helped the Amazon molly to survive until today. If these mitigating processes are weak enough, an increase of mutagenic substances in the environment could easily lead to a rate of mutation accumulation that might allow extinction within time frames that are frequently of interest to conservation biologists.

Methods

Standard null model

To test the hypothesis that Muller's ratchet does not threaten the Amazon molly with extinction during the known time of its existence, we used the null model for quantifying the threat of extinction from Muller's ratchet, as described elsewhere in detail [7]. This model is a simple extension of the standard model of Haigh [104] combined with mutational meltdowns [105]. In short, we combined a multiplicative fitness model with an upper limit for Rmax, the maximal effective number of offspring that can be produced by an individual. This allows computing Cmm, the number of clicks of Muller's ratchet that are needed to start mutational meltdown

Cmm = log(1/Rmax)/log(1 - s),(1)

where each click is the stochastic extinction of the genotype carrying the fewest deleterious mutations that have the constant positive selection coefficient s [see [7,105]]. While this extinction does not require the fixation of a new mutation in the population, such a fixation typically follows shortly after the loss of the best class [106,107], unless there are special circumstances. Combining the effective population size Ne with Usdm, the genomic rate of the origin of

    s
lightly
    d
eleterious
    m
utations with effect s, allows us to compute Tcl, the average effective time between two clicks of Muller's ratchet in generations that are assumed to be discrete. Since computation of this click time is rather difficult, we combine two of the best analytical approximations [9,27,28] with extensive individual-based computer simulations that were distributed over the Internet using evolution@home, the first global computing system for evolutionary biology [7,29-31]. We do not quantify mutational meltdown, because this demographic process at the end of genomic decay is so fast that it can be neglected. Therefore, we approximate Tex, the time to extinction by

Tex = Cmm * Tcl * Tgen(2)

where Tgen is the time between generations that are assumed to be discrete [7]. Extinction times are computed for a large number of parameter combinations that span the whole realistic range of parameters, to assess how many parameter combinations could lead to an extinction by Tage, the known age of existence of the asexual lines of descent. Our simplifying assumption here is that no significant mutation accumulation had occurred in the sexual ancestors, since regular recombination would have facilitated the selective removal of all slightly deleterious mutations of critical effects (i.e. effects that could endanger the long-term survival of the Amazon molly in the presence of Muller's ratchet). To account for the distribution of mutational effects we scale the total genomic mutation rate appropriately to obtain Usdm [7] and we estimate Ne from diversity data [108]. This approach is justified in the corresponding sections below.

Mutation rate estimates

Before we can estimate Usdm, we need to estimate Utot, the mutation rate at all potentially deleterious sites in a haploid genome per generation. To do this we focus on synonymous point mutations for estimating the rate per base pair and later extrapolate to potentially harmful sites that change amino acids or affect regulatory functions. We focus on synonymous point mutations for several reasons.

(i) Synonymous mutations are expected to be mostly effectively neutral (or only under weak selection) and therefore their rate of fixation can be predicted from the neutral theory [109] (the rate is slightly lower if there is very weak selection). Thus synonymous substitution rates are probably very close to the true mutation rates of neighbouring non-synonymous sites that accumulate the actual mutational damage.

(ii) We cannot apply observations from microsatellites, since their mutational model is very different from that of normal point mutations and it is not clear, how an estimate of one rate could be converted to a direct estimate of the other rate.

(iii) One might want to ignore chromosome rearrangements and frame-shift mutations like indels or TE insertions, if one wants to obtain a conservative estimate of extinction time. Many of these mutations are so strongly selected against that they have no chance of accumulating in a large population (their selection coefficients are firmly behind the 'wall of background selection', see Figure 1 and [7]). Thus, such drastic mutations are probably only of importance in the context of corrections for polyploidy that were discussed above and if mutation rates are already very high.

In addition to mutation rates in the nuclear genome, one might want to consider mutation rates in the mitochondrial genome, since mitochondria are probably as essential to fish as they are to humans [7].

Distribution of mutational effects

To quantify Muller's ratchet in the presence of a distribution of mutational effects on fitness we partition this distribution in three as suggested elsewhere [7,108]:

(i) Very deleterious mutations are implicitly dealt with by our method of determining Ne, see below.

(ii) Slightly deleterious mutations with effects in the critical range or close to the critical range are the main focus of our attention here, as these contribute most to extinction, see below.

(iii) Effectively neutral mutations that accumulate like neutral mutations are ignored, as their effects are too small for impacting fitness.

In order to quantify Muller's ratchet we need a good estimate of Usdm, the slightly deleterious mutation rate, which is given by

Usdm = Utot * fsdm,(8)

where Utot is the genomic mutation rate at sites that are in a functional category with potentially deleterious effects like non-synonymous mutations, and fsdm is the fraction of sites with slightly deleterious, critical selection coefficients sc among all potentially deleterious mutations (for more details on this approach, see [7]). The corresponding range of critical selection coefficients approximates all s that lead to extinction within a minimal time or within a given time Tage; this can be determined from Figure 1. To compute fsdm we need to combine the range of all sc with estimates of the distribution of mutational effects on fitness, which traditionally have been difficult to obtain. However recent progress has been substantial and we now know that the distribution of non-synonymous mutational effects is very leptokurtic and spans many orders of magnitude in many different species [46,47,110-112].

Effective population size

There are three population sizes that are potentially relevant here:

(i) total census population size Nt,

(ii) effective population size in the absence of any deleterious mutations Ne0,

(iii) effective population size in the presence of background selection Neb.

Here we argue that Ne = Neb for the purposes of quantifying the rate of Muller's ratchet; this is how we use Ne outside of this section. In order to see this, we need to consider how to approximate the rate of Muller's ratchet in the presence of a wide distribution of mutational effects. New work by Söderberg & Berg [108] has shown that the rate of the ratchet in the presence of such a distribution can be approximated reasonably by dividing mutational effects into different categories (see above). Here we consider the category that contains mutations with very strong effects (see background selection theory [113,114]). Simulations [108] show that the effect of these mutations on the rate of the ratchet is well approximated by using a scaled effective population size Ne0 that is appropriately reduced by the factor Neb/Ne0 so that it only contains individuals that are free from strongly deleterious mutations (see equation 7 in [108]). This is consistent with classical Hill-Robertson effect theory, which states that a locus that is linked to another locus under selection will experience a reduction in effective population size [109].

In the absence of deleterious mutations, estimates of Ne = Ne0 can be made from diversity data and mutation rates, but we cannot estimate Ne0 in our system, since deleterious mutations are present. Surprisingly, we do not need such an estimate, since background selection theory predicts that any estimate of Ne based on diversity and mutation rates will in reality be an estimate of Neb, the effective size of the class that is free from strongly deleterious mutatons [113,114]. This means that two difficulties in our analysis cancel each other out: we do not need to estimate Ne0 and we do not need to simulate background selection for a population of size Ne0 to quantify the rate of the ratchet for mutations with critical effects. We only need to use standard approaches for measuring Ne from DNA sequence diversity; this will approximately give us Neb for our simulations, which then have to ignore all mutations with effects in the background selection range [108]. Thus our analysis uses information from two categories of sites:

(i) Neutral sites are important for estimating Ne = Neb under background selection, but not for the ratchet itself (ignore when computing Usdm).

(ii) Selected sites contribute towards operating the ratchet (use for computing Usdm), but cannot be used to estimate Ne.

Since mtDNA and all nuclear chromosomes are completely linked with the same selective unit in the absence of outcrossing, Ne is the same for both systems like in pure selfers [115,116].

Simplifications overview

Here we approximate the effects of Muller's ratchet in a real diploid genome by proposing to choose (i) effective deleterious mutation rates that are adjusted for the corresponding ploidy level and exclude mutations with effects that are either too large to accumulate or too small to cause effective harm and (ii) effective selection coefficients that are adjusted for arbitrary levels of ploidy, dominance and mitotic recombination. The corresponding adjustments are compiled in Table 1 and explained below. Given the large degree of uncertainty about mutation rates and selection coefficients in our study organism, the errors from these approximations seem insignificant. Also, it is not possible to convincingly decide between the "Core-Genome-Model" and the alternative "Equal-Contribution-Model", although the latter seems to be much more realistic (both models are explained below). Therefore all possible options are explored by including a wide range of possible mutation rates and effects.

Extensions for diploids and polyploids

To apply the model above to non-haploid organisms requires some adjustments to account for the fact that duplicate copies of genes can buffer deleterious effects. To avoid the complexity associated with recent models of the evolution of gene duplicates [117-121], we propose two simple models that allow the reduction of polyploids to an effectively haploid genome model: (i) The Core-Genome-Model assumes that one functional copy of each essential gene is sufficient and all other alleles can be discarded without problems, so that fitness degradation from Muller's ratchet can only start when the 'backup alleles' have been deactivated. (ii) The Equal-Contribution-Model assumes that each allelic copy contributes equally to fitness, so that selection coefficients are reduced and mutation rates are increased proportionally to the ploidy level. Both models are discussed below for diploid and triploid Amazon molly along with the impact of a variable dominance coefficient h on quantifications of Muller's ratchet.

Core-Genome-Model

One might assume that only one copy of each important gene is really needed and that all additional copies from higher ploidy levels of the genome could be regarded as 'neutral backups' that can be deleted without any negative effect, as long as the last functional copy is still working. An extreme interpretation of this assumption implies that all deleterious mutations are completely recessive (dominance coefficient h = 0), which is rather unrealistic in many settings [122-130]. However, in less extreme cases, recessivity may be strong enough to reduce the effects of a complete gene knockout to effective selective neutrality (Nesh < 1). Under the Core-Genome-Model otherwise strongly deleterious mutations like indels or transposable element insertions can disrupt gene function in all additional copies and thus accumulate effectively like neutral mutations. Eventually the core genome will be distributed across all the different haploid copies of the original genome, as it is really only a combination of all the least damaged parts of the genome. Since one frame-shifting mutation is enough to transform a copy of gene i into a pseudogene, we can approximate its time to inactivation, Tko, by

Tko ≈ 1/(li * μbp * fko)(3)

where fko is the factor that specifies the frequency of frame-shift or other knock-out mutations relative to μbp, the synonymous point mutation rate/basepair/generation, and li is the total length of gene i. Lets assume a typical gene has 2,000 base pairs and frame-shifts occur at a rate of μbp * fko ≈ 1*10-9 (probably a lower limit if compared to ≈ 1*10-8 mutations/bp/generation in mice and fko ≈ 1/3 as observed in bacteria [see the 'C' parameter in [131,132]]). In this case one may have to wait on average for about 500,000 generations, before inactivation of such a gene can be expected. Since these events most probably follow a Poisson distribution, it can be expected that a substantial fraction of essential genes would be inactivated very quickly, leaving only one copy that is then maintained by selection. It is these last copies that are then slowly degraded by Muller's ratchet, since they can most probably mutate to such slightly deleterious states that allow the operation of the ratchet. Early inactivation of enough 'backup copies' will give Muller's ratchet extensive periods of time for degrading the core genome (late inactivation may allow for slightly deleterious point mutations in the non-deactivated allele). For example, the rates above suggest that frame-shift accumulation in 20,000 diploid genes would inactivate on average about 1,000 genes (sd = 32) over 25,000 generations. A precise calculation of extinction time in this model would require the integration of the effects of increasing mutation accumulation over time, where the increase is caused by rising deleterious mutation rates due to the progressive inactivation of genes. To avoid these complex calculations, the following two extreme simplifications are proposed.

A lower limit of Tex can be obtained by ignoring all additional copies of the essential core genes. In this case the resulting Usdm to use for the computation of click time is that of the haploid core copy of the genome. An upper limit can be obtained by partitioning time into (i) the accumulation of knockouts while ignoring Muller's ratchet and (ii) the operation of Muller's ratchet while neglecting any further increase of mutation rates due to additional knockouts. Testing several possible ratios of these two times allows minimization of extinction time between the extreme of no time for frameshifts (no extinction due to Muller's ratchet, since mutation rates are too low) and all time for frame-shifts (no extinction, because Muller's ratchet will always need some time to cause an extinction, even with high mutation rates).

In the case of triploids or higher ploidy levels, the Core-Genome-Model needs correspondingly longer for genomic decay, but in no case can additional copies stop genomic decay. Higher ploidy levels are not discussed here, since we assume that the first Amazon molly was diploid. Triploids lose their advantage in the Core-Genome-Model, if they occasionally lose one of their three sets of chromosomes, as they might do [49]. The fact that this is possible indicates that the lost set of chromosomes does not carry any last functional copy of a vital gene. For more details on triploids see the Equal-Contribution-Model.

Equal-Contribution-Model

Diploids

Many mutations of small effects do not seem to be completely recessive, but rather close to codominance, especially if their effects are very small [122-130]. If that is approximately true for most genes, then both copies in a diploid genome are actually needed to produce the required dose of proteins to maintain fitness. In this case, all copies of the genome of the Amazon molly will contribute approximately the same amount of functionality, implying that the size of the mutational target increases two-fold, while effective selection coefficients decrease two-fold, relative to the haploid case (there are twice as many sites to hit and selection only operates against heterozygotes with strength sh, where h = 0.5).

Triploids

It is not clear whether such reasoning can be extended to triploids, as ancestral genomes may have been fine-tuned for diploidy. The fact that triploids can have significant disadvantages if compared directly to diploids [133,134] cautions against a positive functional role of the additional copy of the genome. On the other hand, one could argue that the surprising flexibility of fish in tolerating various ploidy levels [133] allows for an easy incorporation of an additional dosage of proteins provided by a third set of chromosomes. This is not contradicted by the fact that Amazon molly triploids can be observed in the wild [23], as strong purifying selection would predict vanishing frequencies, since the rates of origin of triploids are rather low in this fish [23]. Existing data on expression patterns of muscle proteins in Amazon molly triploids are inconclusive [50]. These complications suggest that triploids may face purifying selection and therefore they are not likely to reach 100% in a population. Depending on the specific mechanism of how triploids incur a possible disadvantage, the inactivation of one copy might occasionally even be beneficial for triploids. In addition, triploids may revert back to diploidy by occasionally losing one of their three sets of chromosomes [49]. The fraction of triploids in the wild seems to fluctuate substantially, but no population with 100% triploids has been observed. Frequencies of triploids for two different samples were in the range of 4%–15% (Lamatsch et al., unpublished) and 3% – 46% [135], indicating that selection against triploids is probably weak. Therefore, equal-contribution-corrections to haploid estimates of Usdm (and s if available) have to be treated with caution in the case of triploids, but appear to be reasonable approximations for diploids.

Varying dominance

The two models above assume either complete recessivity or codominance. Observed dominance levels in other species are frequently somewhere in-between, where mutations with smaller effects tend to be closer to codominance [122-130]. While we know very little about the precise values of the dominance coefficient, h, and the selection coefficient, s, in the Amazon molly, it appears to be highly probable that general findings from other species apply here too. Of special interest here is the fact that the distribution of mutational effects on fitness appears to be very wide on a log scale [46,47]. Thus we can propose a simplification that allows for arbitrary dominance coefficients between h = 0 and h = 1 by making slight adjustments to the effective s used in simulations. The simplification is based on the observation that

(i) deleterious mutations with smaller effects have a higher probability of accumulating in the population than those with larger effects and

(ii) mutations with large enough effects do not accumulate (see the 'wall of background selection' in Figure 1 and in [7]).

To simplify the treatment we speak of the ratchet as if it fixes mutations, although strictly speaking the ratchet will only cause the mutation-free class to go extinct and fixing happens shortly after that by genetic drift (see [106,107]). The following qualitative analysis is supported by published simulation results [106].

We can treat dominant mutations as follows. If and only if the ratchet can fix the first mutation at a site (i.e. individuals that are homozygous for that size go extinct and all individuals in the population become heterozygous) then it will eventually also fix the second mutation (i.e. all individuals become homozygous for the deleterious allele). If a heterozygote cannot be fixed, there is no way the homozygote state can be fixed and the effective mutation rate used to compute the speed of the ratchet should be reduced accordingly to exclude such sites.

In the case of recessive deleterious mutations, there are the two following extremes. If s is small (most mutations), then click times are hardly affected and almost no correction is necessary. If s is large enough to prevent fixation of heterozygotes, then again the effective mutation rate may be reduced, as fixation of the homozygous state is impossible. Only few intermediate s can become heterozygous but not homozygous, as they have to be close to the switch-like transition between 'can accumulate' and 'cannot accumulate' (see the 'wall of background selection' in Figure 1 and [7]) and their homozygous effects have to be too deleterious to accumulate. This simplification relies on the almost switch-like transition between mutation accumulation by Muller's ratchet and complete mutation removal by background selection.

We ignore overdominant sites, assuming that these stay in their optimal heterozygous state. This will lead to a conservative estimate of extinction time, as these sites do not contribute to fitness decay under this assumption. Future models will have to investigate to what extent such an approximation is justified in the presence of a small fraction of overdominant sites [129,130].

Extensions for mitotic recombination

Non-meiotic recombination can play a substantial role in reducing diversity in asexual lineages [40]. To assess the extent to which the rate of Muller's ratchet could be reduced by this process [43,44,136], the following approach was used. We assume that the maximal possible rate of mitotic recombination is equivalent to selfing with free recombination between loci. The true reduction of the rate of Muller's ratchet is probably much smaller, since mitotic recombination is certainly not as effective as free recombination in a selfing population. Thus the strongest possible reduction of the rate of Muller's ratchet from mitotic recombination can be estimated from the analytical work of Heller and Maynard Smith [45]. As explained below, a simple scaling of U and s is enough to reduce predictions to the haploid case.

Muller's ratchet with selfing

Heller and Maynard Smith [45] derived an equation that describes the distribution of deleterious mutations in a selfing population under Muller's ratchet. We assume here that the expectation of the size of the 'best class' in mutation-selection balance mainly determines the rate of Muller's ratchet. This appears to be supported by some simulation results [43,44] although Gordo & Charlesworth [137] found that the ratchet could click at different rates for the same size of the best class if other parameters were different. Our assumption allows us to use a simple scaling of U and s to reduce predictions for diploid selfers to the haploid case.

In detail, Heller and Maynard Smith [45] show that in a selfing population xij, the fraction of individuals carrying i homozygous and j heterozygous mutations is given by

<a onClick="popup('http://www.biomedcentral.com/1471-2148/8/88/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/8/88/mathml/M11">View MathML</a>

(4)

where Usdm, d is the slightly deleterious mutation rate/diploid genome/generation that changes wildtype homozygotes into heterozygotes, h is the dominance coefficient and s is the homozygous selection coefficient (positive for harmful mutations). This equation assumes that (i) recombination is free between sites, so that the ratchet only operates at homozygous sites, since selfing can restore the best class on a heterozygous site, (ii) Usdm, d is constant and independent from the number of mutations per individual, (iii) mutations at heterozygous sites are negligible, since selfing produces many more homozygotes, (iv) the standard model of Muller's ratchet without back mutations is valid [7,104]. Then the number of individuals in the 'best class', N0, can be computed by adding up all heterozygotes that are free from homozygous deleterious mutations and then scaling by the effective population size

<a onClick="popup('http://www.biomedcentral.com/1471-2148/8/88/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/8/88/mathml/M12">View MathML</a>

(5)

Using the well known result that

<a onClick="popup('http://www.biomedcentral.com/1471-2148/8/88/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/8/88/mathml/M13">View MathML</a>

(6)

equations (4) – (6) yield

<a onClick="popup('http://www.biomedcentral.com/1471-2148/8/88/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/8/88/mathml/M14">View MathML</a>

(7)

For complete recessivity (h = 0) this reduces to the haploid asexual case (<a onClick="popup('http://www.biomedcentral.com/1471-2148/8/88/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/8/88/mathml/M10">View MathML</a>), because the diploid mutation rate Usdm, d equals twice the haploid rate Usdm. It is easy to see from (7) that the degree of dominance of mutations with small selection coefficients does not significantly affect numerical values and hence the rate of the ratchet. Since the ratchet does not click for large s, one might as well neglect dominance altogether.

Parameter estimates

Effective population size

For our simulations we need an estimate of the effective population size in the presence of background selection (Ne = Neb, see above). To arrive at such an estimate we obtain corresponding diversity data from the control region of mtDNA, since it is linked to the rest of the genome by unisexual inheritance (see the Models section above for a justification of this approach). We use the DNA sequence diversity data from Möller (see pp. 46, 47 and 55 in [12]), who observed 14 haplotypes with S = 13 sites segregating for single base pair polymorphisms in a total of n = 63 sequences of length L = 886 bp from 17 geographic locations. Since the control region of mtDNA is likely to evolve neutrally at most sites, we can use Watterson's [138] equation to compute Ne = S/(L * μbp * a), where Ne is the effective population size of mitochondrial DNA in the presence of background selection, μbp is the mutation rate per site/generation and a = 1 + 1/2 + 1/3 + ... + 1/(n-1) for n sequences in the sample. Here we assume an mtDNA divergence rate of 3.6% ± 0.46%/site/Myr as estimated in the species groups of the fish Centropomus [139] from divergence since the Panama seaway closed 3.0–3.5 Myr ago [140]. This suggests a value for Ne between 153,000 and 198,000 that will have to be scaled downward, if actual mutation rates are higher, like in human pedigrees [141]. The nuclear value is identical to this mtDNA estimate, if there is no outcrossing [115,116].

Comparative analyses suggest that the resulting values are reasonable. The total population size Nt may be used to place upper bounds on Ne in most situations. Some mark and recapture experiments could only estimate that Nt > 10,000 in that local study area (M. Doebler, personal communication). One may speculate that about 10 million individuals of P. formosa may exist in total. If this is combined with the typical finding that vertebrates frequently have Ne-values that are about 10% of their census population sizes [142], then Ne would be about a million. Given that a species as abundant as the fruitfly Drosophila melanogaster has effective population sizes of little more than a million [143], suggests that the Amazon molly probably has Ne < 106. On the other extreme it is highly unlikely that the Amazon molly would have an effective population size smaller than that of primates, suggesting Ne > 104.

We use deliberately wide upper and lower bounds of Ne = 104 to 107 to demonstrate that population size has a rather small influence on the operation of Muller's ratchet in the presence of a distribution of mutational effects.

Life history

Based on observations we estimate that generation time Tgen is between 4 month (min) and 2.5 years (max), so we assume 1 year as our middle value. The maximal reproductive capacity in the wild with a non-degraded genome Rmax is assumed to be 50 (min) 500 (mid) and 2000 (max) offspring per lifetime (see below). While Tgen scales extinction times linearly, Rmax has a much smaller effect. It is possible that our values for Rmax are too large since they are more based on the highest possible values conceivable under laboratory conditions than on any real life survival situation. Once such data are available, results from Figure 1 can be easily scaled appropriately (see equation 2).

Unfortunately there is very little information about life history from natural habitats. Most of the knowledge about P. formosa comes from extensive aquarium cultures. In captivity the Amazon molly reaches maturity between 4 and 6 months, produces more offspring at larger body sizes and may live for up to 3 years under optimal laboratory conditions. Here we report that females caught in the wild contain a mean of 28 embryos (sd ± 16; n = 20; observed maximum in the wild = 70 embryos, slightly less than the maximum of 90 observed in the lab [144]). Females will probably produce this number of newborn fish about every 28 days, depending on water temperature. They stop reproduction about 3 to 6 months before death in the lab. The number of surviving newborn fish can be reduced considerably by predation or disease. Thus under good conditions Rmax might be 18 periods * 28 offspring = 504, assuming a life span of 2 years. As a lower limit one might use 50 offspring per lifetime, assuming that many offspring die due to predation or disease and many parents are small and frequently trapped in little ponds without sister species that provide the sperm that they need as a trigger for development. An upper limit of Rmax = 2,030 might be derived from combining the maximal number of offspring-producing periods observed in the lab (29 = 36-4-3) with the maximal number of embryos per period observed in the wild (70 in n = 20, the upper expectation for a fish of this size [10]). We did not find that the sexual sister species, P. mexicana limantouri differs significantly from these values (observed mean = 27 embryos ± sd = 13; n = 11).

Age of the evolutionary line

The ancestor of the Amazon molly clone known to us was formed by a hybridisation event between a Poecilia mexicana female and a Poecilia latipinna male [11,12] and is believed to have reproduced asexually since then. A date Tage for this event may be derived from comparing alleles in the Amazon molly with alleles in their corresponding parent species using a molecular clock. At the moment there are two datasets that are large enough to infer an evolutionary age that is different from zero. Comparisons of 1,377 bp in various nuclear genes (mostly introns) lead to a point estimate of Tage = 100,000 years assuming a divergence rate of 2%/Myr from standard mitochondrial clocks [11]. This age may be an underestimate, since nuclear mutation rates are usually much lower than mitochondrial rates. However, the confidence intervals for this estimate range from Tage = 8,930 to Tage = 209,000 years due to the small number of substitutions observed [69,70]. Thus it is desirable to arrive at a more precise estimate. This can be done using the same set of mtDNA diversity data that was used to estimate Ne above [12]. Analyses of the median joining network show that there is only one mtDNA haplotype that is shared between P. formosa and P. mexicana and that this haplotype is at the center of P. formosa diversity (p.49 in [12]). Thus it is reasonable to assume that this haplotype was the most recent common ancestor (MRCA) of all Amazon mollies and one can use the number of mutational steps back to the MRCA as an indicator of the time since divergence of Amazon mollies from their MRCA. Using the method described by Saillard et al. [145] and the data of Möller [see p.69 [12]] leads to an estimate of 1.294 (± sd = 0.159) mutational steps back to the MRCA in the 886 bp of the control region of mtDNA that were sequenced for this purpose. If this is combined with the divergence rate of 3.6% ± 0.46%/site/Myr as estimated in the species groups of the fish Centropomus [139] from divergence since the Panama seaway closed 3.0–3.5 Myr ago [140], then we arrive at an estimate of Tage = 81,000 years for the time to the MRCA of P. formosa (1.294/(886 * 3.6 *10-8/2); lower and upper limits range from 63,000 to 104,000 years). This time will have to be scaled downward, if actual mutation rates are higher, as has been found in human pedigrees [141]. Thus we have marked Tage = 40 Kyr, 70 Kyr and 100 Kyr as estimates of the age of the Amazon molly in Figure 1. It is easy to check our conclusions, if more precise information becomes available. A more rigorous analysis involving multiple genes from multiple individuals from each species would be desirable [see [146]].

Mutation rate estimates

The last decades have seen a large increase of our understanding of mutation rates across a wide range of organisms [132,147,148] and we have no reason to assume hat the Amazon molly might have extraordinarily low mutation rates. Thus it is possible in the absence of more specific data to extrapolate from other species to arrive at a credible estimate for the Amazon molly specific mutation rate at all potentially deleterious sites in a haploid genome per generation, Utot.

The size of a haploid genome is about 950 Mbp, as derived from DNA content of cells [149]. If this is combined with the mutation rate counterpart of the divergence rate that was used to date the age of Amazon molly [11], then about 9.5 new mutations per haploid genome can be expected in each new generation (950 Mbp * 0.01 subst/bp/Myr assuming Tgen = 1 year). Since many of these mutations will be synonymous or affect only non-functional DNA, we have to scale this rate by the effective genome size [132]. Fish are presumably more complex than flies. Therefore, we may use the effective genome size of Drosophila melanogaster as a lower limit (13,379 genes with 54,934 exons and a total length of 27.8 Mbp in exons [150] suggest about 19 Mbp as the target size for non-synonymous, non-frameshifting mutations). This is just 2% of the actual genome size and has to be increased to about 8% to account for non-coding functional sequences that seem to be about three times as abundant as functional sequences in exons in organisms as diverse as fruitflies and mice [151,152]. This would result in Utot ≈ 0.74 potentially deleterious mutations/haploid genome/generation and may be an underestimate, as the amount of coding DNA in fish might be higher than in Drosophila. If a similar calculation is based on the 33 Mbp exons that we can estimate from the Fugu genome project, then Utot ≈ 0.88 seems plausible (all genes occupy 108 Mbp and we estimate that this includes roughly 75 Mbp introns from Figure 2 in [153]).

These estimates do not change much, even if we use a completely different way of obtaining the mutation rate. Observations in mice have lead to estimates of 1.8 × 10-10/site/cell replication, respectively [132]. If an estimated 88 Mbp functional sites in P. formosa experience about 25 germ line cell divisions per adult generation (as female mice [132] that have a similar body weight to the Amazon molly) then we arrive at Utot ≈ 0.4 potentially deleterious mutations/haploid genome/generation. Obviously, P. formosa is not a mouse and there is much room for improving the precision of such mutation rate estimates, even in mice. Recent high-precision measurements of the deleterious mutation rate in mutation accumulation experiments in Drosophila found Utot ≈ 0.6 deleterious mutations/haploid genome/generation [154]. Scaling this estimate by the relative effective genome size of Fugu suggests Utot ≈ 0.7.

It is encouraging to find such general agreement between these estimates that are based on very different approaches. This suggests that we can have some confidence in our estimates, even in the absence of direct observations in the Amazon molly. Since mutations are the inevitable consequence of DNA replication errors and the number of mitotic cell divisions per generation is a key determinant of mutation rates per generation [132,155], we conclude that the absence of recombination in the Amazon molly is unlikely to cause a major change in mutation rate when compared to its sexual sister species.

Distribution of mutational effects

Recent progress has shown that the distribution of non-synonymous mutational effects is very leptokurtic and spans many orders of magnitude in many different species [46,47,110-112]. While we have no specific data for the Amazon molly, we have no reason to believe that this result does not hold here as well. One can use data on the ratio of non-synonymous to synonymous substitution rates, KA/KS, to infer that the fraction 1-KA/KS of all mutations must be more deleterious than about s = 1/Ne, since mutations with effects below that threshold accumulate as if they were neutral [47,156]. If the data in Drosophila can be used as a landmark in this new territory, then about 90% of all mutations are more deleterious than the limit of 1/Ne ≈ 10-6 [112]. Estimates of the distribution of mutational effects in Drosophila have shown that this distribution may well follow a lognormal law and there may be good theoretical reasons for this [46]. Robust features are its large width on a logscale and the fact that most probability mass is between effective neutrality and lethality [46,47]. From this and from Figure 1 we may infer that the fraction of slightly deleterious mutations, fsdm, could be perhaps around 30%, probably larger than 10%, but probably not much larger than 50%. We use these values as a point estimate with lower and upper bound. Thus Usdm probably shares the same order of magnitude with Utot. As this calculation has large errors, an array of other mutation rates is included as well when quantifying the ratchet. We assume that selection coefficients did not change much for most genes in the time since asexuality arose in the Amazon molly. This is compatible with the assumption that most genes in the Amazon molly are generic to similar fish species and thus well adapted, while only a small fraction is actually responsible for the specific adaptations of the Amazon molly.

Deleterious mutation rates summary

As discussed above, our best estimate for the total mutation rate at potentially deleterious sites per haploid genome per generation is expected to be between Utot ≈ 0.4 and Utot ≈ 0.9, while our best estimate of the fraction of mutations with critical slightly deleterious effects is expected to be between fsdm ≈ 10% and fsdm ≈ 50%. Combining these values makes us expect Usdm ≈ 0.2 slightly deleterious mutations with critical effects/haploid genome/generation with a lower limit of Usdm ≈ 0.04 and an upper limit of Usdm ≈ 0.45. These values may have to be doubled, depending on the genome model that is used for accounting for the effects of diploidy. To err on the side of caution and provide a better feeling for the effects of mutation rates we plot values from Usdm = 0.01 to Usdm = 1 in Figure 1.

Authors' contributions

LL conceived the study, wrote the software, managed the computation of simulation results, developed the theoretical aspects of the paper, estimated many parameters, analyzed the results and wrote the manuscript. DKL provided most details on the Amazon molly biology, contributed to parameter estimation and helped writing the manuscript. Both authors read and approved the final manuscript.

Acknowledgements

Special thanks go to the more than 300 participants of evolution@home (especially to David Corbett, Rechenkraft.net and Seth A. Keel) for more than 80 years of CPU-time that are the basis of this analysis. We thank Kathrin Lampert and Ute Hornung for help in the field and Manfred Schartl, Brian Charlesworth and Asher Cutter for fruitful discussions and helpful comments on the manuscript. We also thank four anonymous reviewers for their feedback. Financial support came from the Leverhulme Trust, the Deutsche Forschungsgemeinschaft (SFB 567 "Mechanismen der interspezifischen Interaktion von Organismen") and EU Marie Curie Research Training Network 'SexAsex' (MCRTN-CT-2004-512492).

References

  1. Butlin R: Opinion – evolution of sex: The costs and benefits of sex: new insights from old asexual lineages.

    Nat Rev Genet 2002, 3:311-317. PubMed Abstract | Publisher Full Text OpenURL

  2. Martens K, Schon I: Parasites, predators and the Red Queen.

    Trends Ecol Evol 2000, 15:392-393. PubMed Abstract | Publisher Full Text OpenURL

  3. Maynard Smith J: The evolution of sex. New York: Cambridge University Press; 1978.

  4. Judson OP, Normark BB: Ancient asexual scandals.

    Trends Ecol Evol 1996, 11:41-46. OpenURL

  5. West SA, Lively CM, Read AF: A pluralist approach to sex and recombination.

    J Evol Biol 1999, 12:1003-1012. OpenURL

  6. Howard RS, Lively CM: Parasitism, mutation accumulation and the maintenance of sex.

    Nature 1994, 367:554-557. PubMed Abstract | Publisher Full Text OpenURL

  7. Loewe L: Quantifying the genomic decay paradox due to Muller's ratchet in human mitochondrial DNA.

    Genet Res 2006, 87:133-159. PubMed Abstract | Publisher Full Text OpenURL

  8. Kondrashov AS: Classification of hypotheses on the advantage of amphimixis.

    J Hered 1993, 84:372-387. PubMed Abstract | Publisher Full Text OpenURL

  9. Stephan W, Kim Y: Recent applications of diffusion theory to population genetics. In Modern developments in theoretical population genetics. Edited by Slatkin M, Veuille M. Oxford: Oxford University Press; 2002:72-93. OpenURL

  10. Hubbs CL, Hubbs LC: Apparent parthogenesis in nature, in a form of fish of hybrid origin.

    Science 1932, 76:728-630. OpenURL

  11. Schartl M, Wilde B, Schlupp I, Parzefall J: Evolutionary origin of a parthenoform, the Amazon molly Poecilia formosa, on the basis of a molecular genealogy.

    Evolution 1995, 49:827-835. OpenURL

  12. Möller D: Aspekte zur Populationsgenetik des eingeschlechtlichen Amazonenkärpflings Poecilia formosa (GIRARD 1859) unter Berücksichtigung der genetischen parentalen Arten, dem Breitflossenkärpfling Poecilia latipinna (LESUEUR 1821) und dem Atlantikkärpfling Poecilia mexicana (STEINDACHNER 1863). Hamburg, Germany: University of Hamburg; 2001.

  13. Schartl M, Nanda I, Schlupp I, Wilde B, Epplen JT, Schmid M, Parzefall J: Incorporation of subgenomic amounts of DNA as compensation for mutational load in a gynogenetic fish.

    Nature 1995, 373:68-71. OpenURL

  14. Schlupp I, Parzefall J, Schartl M: Biogeography of the Amazon molly, Poecilia formosa.

    J Biogeogr 2002, 29:1-6. OpenURL

  15. Lampert KP, Lamatsch DK, Epplen JT, Schartl M: Evidence for a monophyletic origin of triploid clones of the Amazon molly, Poecilia formosa.

    Evolution 2005, 59:881-889. PubMed Abstract OpenURL

  16. Maruyama T: On fixation probability of mutant genes in a subdivided population.

    Genet Res 1970, 15:221-225. PubMed Abstract OpenURL

  17. Maruyama T: Simple proof that certain quantities are independent of geographical structure of population.

    Theor Popul Biol 1974, 5:148-154. PubMed Abstract OpenURL

  18. Whitlock MC, Barton NH: The effective size of a subdivided population.

    Genetics 1997, 146:427-441. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Whitlock MC: Fixation probability and time in subdivided populations.

    Genetics 2003, 164:767-779. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Hubbs CL: Interactions between bisexual fish species and its gynogenetic sexual parasite.

    Bull Tex Mem Mus 1964, 8:1-72. OpenURL

  21. Schlupp I, et al.: Dispensable and indispensable genes in an ameiotic fish, the Amazon molly Poecilia formosa.

    Cytogenet Cell Genet 1998, 80:193-198. PubMed Abstract | Publisher Full Text OpenURL

  22. Stöck M, Lamatsch DK: Triploide Wirbeltiere: Wege aus der Unfruchtbarkeit oder Eingeschlechtigkeit.

    Naturw Rdsch 2002, 55:349-358. OpenURL

  23. Schories S, Lampert KP, Lamatsch DK, de Leon FJ, Schartl M: Analysis of a possible independent origin of triploid P. formosa outside of the Río Purificación river system.

    Front Zool 2007, 4:13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Lamatsch DK, Nanda I, Schlupp I, Epplen JT, Schmid M, Schartl M: Distribution and stability of supernumerary microchromosomes in natural populations of the Amazon molly, Poecilia formosa.

    Cytogenet Genome Res 2004, 106:189-194. PubMed Abstract | Publisher Full Text OpenURL

  25. Nanda I, Schlupp I, Lamatsch D, Lampert K, Schmid M, Schartl M: Stable inheritance of host species-derived microchromosomes in the gynogenetic fish Poecilia formosa.

    Genetics 2007, 177:917-926. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Rasch EM, Monaco PJ, Balsano JS: Cytophotometric and autoradiographic evidence for functional apomixis in a gynogenetic fish, Poecilia formosa and its related, triploid unisexuals.

    Histochemistry 1982, 73:515-533. PubMed Abstract OpenURL

  27. Gessler DDG: The constraints of finite size in asexual populations and the rate of the ratchet.

    Genet Res 1995, 66:241-253. PubMed Abstract OpenURL

  28. Stephan W, Chao L, Smale JG: The advance of Muller's ratchet in a haploid asexual population: Approximate solutions based on diffusion theory.

    Genet Res 1993, 61:225-231. PubMed Abstract OpenURL

  29. Loewe L: evolution@home: Experiences with work units that span more than 7 orders of magnitude in computational complexity. [http://archive.evolution.ws/pdf/loewe-l/Loewe2002-EaHworkunits.pdf] webcite

    In Proceedings of the 2nd IEEE/ACM International Symposium on Cluster Computing and the Grid (CCGrid2002) Edited by Bal HE, Löhr K-P, Reinefeld A. Berlin, Germany: IEEE Computer Society; 2002, 425-431. OpenURL

  30. Loewe L: Evolution@home: observations on participant choice, work unit variation and low-effort global computing.

    Softw Pract Exper 2007, 37:1289-1318. OpenURL

  31. Loewe L: The evolution@home website. [http://evolutionary-research.net] webcite

    2008.

  32. Sola L, Galetti PM, Monaco PJ, Rasch EM: Cytogenetics of bisexual/unisexual species of Poecilia. 6. Additional nucleolus organizer region chromosomal clones of Poecilia formosa (Amazon molly) from Texas, with a survey of chromosomal clones detected in the Amazon molly.

    Heredity 1997, 78:612-619. OpenURL

  33. Stern C: Somatic crossing over and segregation in Drosophila melanogaster.

    Genetics 1936, 21:625-730. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Ashburner M: Drosophila: A laboratory handbook. Cold Spring Harbor, N.Y.: Cold Spring Harbor Laboratory; 1989.

  35. Prado F, Cortes-Ledesma F, Huertas P, Aguilera A: Mitotic recombination in Saccharomyces cerevisiae.

    Curr Genet 2003, 42:185-198. PubMed Abstract | Publisher Full Text OpenURL

  36. Aguilera A, Chavez S, Malagon F: Mitotic recombination in yeast: elements controlling its incidence.

    Yeast 2000, 16:731-754. PubMed Abstract | Publisher Full Text OpenURL

  37. Barbera MA, Petes TD: Selection and analysis of spontaneous reciprocal mitotic cross-overs in Saccharomyces cerevisiae.

    Proc Natl Acad Sci USA 2006, 103:12819-12824. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Shao C, Stambrook PJ, Tischfield JA: Mitotic recombination is suppressed by chromosomal divergence in hybrids of distantly related mouse strains.

    Nat Genet 2001, 28:169-172. PubMed Abstract | Publisher Full Text OpenURL

  39. Holt D, Dreimanis M, Pfeiffer M, Firgaira F, Morley A, Turner D: Interindividual variation in mitotic recombination.

    Am J Hum Genet 1999, 65:1423-1427. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Omilian AR, Cristescu ME, Dudycha JL, Lynch M: Ameiotic recombination in asexual lineages of Daphnia.

    Proc Natl Acad Sci USA 2006, 103:18638-18643. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Holliday R: Recombination and meiosis.

    Philos Trans R Soc Lond B Biol Sci 1977, 277:359-370. PubMed Abstract OpenURL

  42. Helleday T: Pathways for mitotic homologous recombination in mammalian cells.

    Mutat Res 2003, 532:103-115. PubMed Abstract | Publisher Full Text OpenURL

  43. Antezana MA, Hudson RR: Before crossing over: The advantages of eukaryotic sex in genomes lacking chiasmatic recombination.

    Genet Res 1997, 70:7-25. PubMed Abstract OpenURL

  44. Charlesworth D, Morgan MT, Charlesworth B: Mutation accumulation in finite outbreeding and inbreeding populations.

    Genet Res 1993, 61:39-56. OpenURL

  45. Heller R, Maynard Smith J: Does Muller's ratchet work with selfing?

    Genet Res 1978, 32:289-293. OpenURL

  46. Loewe L, Charlesworth B: Inferring the distribution of mutational effects on fitness in Drosophila.

    Biol Lett 2006, 2:426-430. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Eyre-Walker A, Keightley PD: The distribution of fitness effects of new mutations.

    Nat Rev Genet 2007, 8:610-618. PubMed Abstract | Publisher Full Text OpenURL

  48. Beukeboom LW, Vrijenhoek RC: Evolutionary genetics and ecology of sperm-dependent parthenogenesis.

    J Evol Biol 1998, 11:755-782. OpenURL

  49. Lamatsch DK, Schmid M, Schartl M: A somatic mosaic of the gynogenetic Amazon molly.

    J Fish Biol 2002, 60:1417-1422. OpenURL

  50. Monaco PJ, Rasch EM, Balsano JS, Turner BJ: Muscle protein phenotypes and the probable evolutionary origin of a unisexual fish, Poecilia formosa, and its triploid derivatives.

    J Exp Zool 1982, 221:265-274. OpenURL

  51. Monaco P, Rasch E, Balsano J: Apomictic reproduction in the Amazon molly, Poecilia formosa and its triploid hybrids. In Evolutionary Genetics of Fishes. Edited by Turner BJ. New York: Plenum Press; 1984:311-318. OpenURL

  52. Beukeboom LW, van Batenburg FHD: The effect of paternal leakage on the rate of Muller's ratchet.

    7th Conference of the European Society for Evolutionary Biology: 1999; Barcelona, Spain 1999.

    handout for poster.

    OpenURL

  53. Beukeboom LW, Weinzierl RP, Michiels NK: Amazon molly and Muller's ratchet.

    Nature 1995, 375:111-112. OpenURL

  54. Crow JF: The high spontaneous mutation rate: Is it a health risk?

    Proc Natl Acad Sci USA 1997, 94:8380-8386. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Kondrashov AS: Muller's Ratchet under epistatic selection.

    Genetics 1994, 136:1469-1473. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Butcher D: Muller's ratchet, epistasis and mutation effects.

    Genetics 1995, 141:431-437. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Poon A, Otto SP: Compensating for our load of mutations: Freezing the meltdown of small populations.

    Evolution 2000, 54:1467-1479. PubMed Abstract OpenURL

  58. Bachtrog D, Gordo I: Adaptive evolution of asexual populations under Muller's ratchet.

    Evolution 2004, 58:1403-1413. PubMed Abstract OpenURL

  59. Schultz ST, Lynch M: Mutation and extinction: The role of variable mutational effects, synergistic epistasis, beneficial mutations, and the degree of outcrossing.

    Evolution 1997, 51:1363-1371. OpenURL

  60. Eyre-Walker A: The genomic rate of adaptive evolution.

    Trends Ecol Evol 2006, 21:569-575. PubMed Abstract | Publisher Full Text OpenURL

  61. Silander OK, Tenaillon O, Chao L: Understanding the evolutionary fate of finite populations: the dynamics of mutational effects.

    PLoS Biol 2007, 5:e94. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Plath M: Cave molly females (Poecilia mexicana) avoid parasitised males.

    Acta Ethol 2004, 6:47-51. OpenURL

  63. Tobler M, Plath M, Burmeister H, Schlupp I: Black spots and female association preferences in a sexual/asexual mating complex (Poecilia, Poeciliidae, Teleostei).

    Behav Ecol Sociobiol 2006, 60:159-165. OpenURL

  64. Jaenike J: An hypothesis to account for the maintenance of sex within populations.

    Evolutionary Theory 1978, 3:191-194. OpenURL

  65. Sasaki A, Hamilton WD, Ubeda F: Clone mixtures and a pacemaker: new facets of Red-Queen theory and ecology.

    Proceedings of the Royal Society of London Series B-Biological Sciences 2002, 269:761-772. OpenURL

  66. Meirmans S, Neiman M: Methodologies for testing a pluralist idea for the maintenance of sex.

    Biol J Linn Soc 2006, 89:605-613. OpenURL

  67. Lively CM, Craddock C, Vrijenhoek RC: Red Queen hypothesis supported by parasitism in sexual and clonal fish.

    Nature 1990, 344:864-866. OpenURL

  68. Tobler M, Schlupp I: Parasites in sexual and asexual mollies (Poecilia, Poeciliidae, Teleostei): a case for the Red Queen?

    Biol Lett 2005, 1:166-168. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  69. Dries LA: The evolutionary persistence of the gynogenetic Amazon molly, Poecilia formosa. Austin, Texas: University of Texas; 2000.

  70. Dries LA: Peering through the looking glass at a sexual parasite: Are Amazon mollies red queens?

    Evolution 2003, 57:1387-1396. PubMed Abstract OpenURL

  71. Lampert KP, Lamatsch DK, Fischer P, Epplen JT, Nanda I, Schmid M, Schartl M: Automictic reproduction in interspecific hybrids of poeciliid fish.

    Curr Biol 2007, 17:1948-1953. PubMed Abstract | Publisher Full Text OpenURL

  72. Schlupp I: The evolutionary ecology of gynogenesis.

    Annual Review of Ecology Evolution and Systematics 2005, 36:399-417. OpenURL

  73. Vrijenhoek RC: Unisexual Fish – Model Systems for Studying Ecology and Evolution.

    Annu Rev Ecol Syst 1994, 25:71-96. OpenURL

  74. Vrijenhhoek RC: The origin and evolution of clones versus the maintenance of sex in Poeciliopsis.

    J Hered 1993, 84:388-395. OpenURL

  75. Schoen DJ, Brownstein JS, White PA: Mutation load in natural populations of the sensitive fern Onoclea sensibilis exposed to soil mutagens.

    Ecol Applic 2002, 12:124-137. OpenURL

  76. Klekowski EJ: Mutational load in a fern population growing in a polluted environment.

    Am J Bot 1976, 63:1024-1030. OpenURL

  77. Klekowski EJ, Klekowski E: Mutation in ferns growing in an environment contaminated with polychlorinated-biphenyls.

    Am J Bot 1982, 69:721-727. OpenURL

  78. Majer BJ, Grummt T, Uhl M, Knasmuller S: Use of plant bioassays for the detection of genotoxins in the aquatic environment.

    Acta Hydrochim Hydrobiol 2005, 33:45-55. OpenURL

  79. Bell G: The Masterpiece of Nature: The Evolution and Genetics of Sexuality. Berkeley: University of California Press; 1982.

  80. Martens K, Rossetti G, Horne DJ: How ancient are ancient asexuals?

    Proc R Soc Lond B Biol Sci 2003, 270:723-729. OpenURL

  81. Leslie JF, Vrijenhoek RC: Consideration of Muller's ratchet mechanism through studies of genetic-linkage and genomic compatibilities in clonally reproducing Poeciliopsis.

    Evolution 1980, 34:1105-1115. OpenURL

  82. Spolsky CM, Phillips CA, Uzzell T: Antiquity of clonal salamander lineages revealed by mitochondrial DNA.

    Nature 1992, 356:706-708. PubMed Abstract | Publisher Full Text OpenURL

  83. Hedges SB, Bogart JP, Maxson LR: Ancestry of unisexual salamanders.

    Nature 1992, 356:708-710. PubMed Abstract | Publisher Full Text OpenURL

  84. Robertson AV, Ramsden C, Niedzwiecki J, Fu J, Bogart JP: An unexpected recent ancestor of unisexual Ambystoma.

    Mol Ecol 2006, 15:3339-3351. PubMed Abstract | Publisher Full Text OpenURL

  85. Guex GD, Hotz H, Semlitsch RD: Deleterious alleles and differential viability in progeny of natural hemiclonal frogs.

    Evolution 2002, 56:1036-1044. PubMed Abstract OpenURL

  86. Som C, Bagheri HC, Reyer HU: Mutation accumulation and fitness effects in hybridogenetic populations: a comparison to sexual and asexual systems.

    BMC Evol Biol 2007, 7:80. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  87. Som C, Reyer HU: Hemiclonal reproduction slows down the speed of Muller's ratchet in the hybridogenetic frog Rana esculenta.

    J Evol Biol 2007, 20:650-660. PubMed Abstract | Publisher Full Text OpenURL

  88. Vorburger C: Fixation of deleterious mutations in clonal lineages: evidence from hybridogenetic frogs.

    Evolution 2001, 55:2319-2332. PubMed Abstract OpenURL

  89. Klekowski EJ: Plant clonality, mutation, diplontic selection and mutational meltdown.

    Biol J Linn Soc 2003, 79:61-67. OpenURL

  90. Smith RJ, Kamiya T, Horne DJ: Living males of the 'ancient asexual' Darwinulidae (Ostracoda: Crustacea).

    Proc R Soc Biol Sci Ser B 2006, 273:1569-1578. OpenURL

  91. Martens K: Sex and ostracods: a new synthesis. In Sex and parthenogenesis: Evolutionary ecology of reproductive modes in non-marine ostracods. Edited by Martens K. Leiden, Netherlands: Backhuys Publishers; 1998:295-321. OpenURL

  92. Butlin R, Schön I, Martens K: Asexual reproduction in non-marine ostracods.

    Heredity 1998, 81:473-480. OpenURL

  93. Welch DM, Meselson M: Evidence for the evolution of Bdelloid Rotifers without sexual reproduction or genetic exchange.

    Science 2000, 288:1211-1215. PubMed Abstract | Publisher Full Text OpenURL

  94. Judson OP, Normark BB: Evolutionary genetics. Sinless originals.

    Science 2000, 288:1185-1186. PubMed Abstract | Publisher Full Text OpenURL

  95. Butlin RK: Virgin rotifers.

    Trends Ecol Evol 2000, 15:389-390. PubMed Abstract | Publisher Full Text OpenURL

  96. Dolgin ES, Charlesworth B: The fate of transposable elements in asexual populations.

    Genetics 2006, 174:817-827. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  97. Arkhipova I, Meselson M: Deleterious transposable elements and the extinction of asexuals.

    Bioessays 2005, 27:76-85. PubMed Abstract | Publisher Full Text OpenURL

  98. Pouchkina-Stantcheva NN, McGee BM, Boschetti C, Tolleter D, Chakrabortee S, Popova AV, Meersman F, Macherel D, Hincha DK, Tunnacliffe A: Functional divergence of former alleles in an ancient asexual invertebrate.

    Science 2007, 318:268-271. PubMed Abstract | Publisher Full Text OpenURL

  99. Meselson M, Welch DM: Evolution. Stable heterozygosity?

    Science 2007, 318:202-203. PubMed Abstract | Publisher Full Text OpenURL

  100. Maraun M, Heethoff M, Scheu S, Norton RA, Weigmann G, Thomas RH: Radiation in sexual and parthenogenetic oribatid mites (Oribatida, Acari) as indicated by genetic divergence of closely related species.

    Experimental and Applied Acarology 2003, 29:265-277. OpenURL

  101. Maraun M, Heethoff M, Schneider K, Scheu S, Weigmann G, Cianciolo J, Thomas RH, Norton RA: Molecular phylogeny of oribatid mites (Oribatida, Acari): evidence for multiple radiations of parthenogenetic lineages.

    Experimental and Applied Acarology 2004, 33:183-201. OpenURL

  102. Heethoff M, Domes K, Laumann M, Maraun M, Norton RA, Scheu S: High genetic divergences indicate ancient separation of parthenogenetic lineages of the oribatid mite Platynothrus peltifer (Acari, Oribatida).

    J Evol Biol 2007, 20:392-402. PubMed Abstract | Publisher Full Text OpenURL

  103. Loewe L: How to quantify effects of Muller's ratchet. [http://evolutionary-research.net/science/mullers-ratchet/] webcite

    2008.

  104. Haigh J: The accumulation of deleterious genes in a population – Muller's ratchet.

    Theor Popul Biol 1978, 14:251-267. PubMed Abstract | Publisher Full Text OpenURL

  105. Lynch M, Butcher RBD, Gabriel W: The mutational meltdown in asexual populations.

    J Hered 1993, 84:339-344. PubMed Abstract | Publisher Full Text OpenURL

  106. Charlesworth B, Charlesworth D: Rapid fixation of deleterious alleles can be caused by Muller's ratchet.

    Genet Res 1997, 70:63-73. PubMed Abstract OpenURL

  107. Higgs PG, Woodcock G: The accumulation of mutations in asexual populations and the structure of genealogical trees in the presence of selection.

    J Math Biol 1995, 33:677-702. OpenURL

  108. Söderberg RJ, Berg OG: Mutational interference and the progression of Muller's ratchet when mutations have a broad range of effects.

    Genetics 2007, 177:971-986. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  109. Birky CW Jr, Walsh JB: Effects of linkage on rates of molecular evolution.

    Proc Natl Acad Sci USA 1988, 85:6414-6418. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  110. Eyre-Walker A, Woolfit M, Phelps T: The distribution of fitness effects of new deleterious amino acid mutations in humans.

    Genetics 2006, 173:891-900. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  111. Piganeau G, Eyre-Walker A: Estimating the distribution of fitness effects from DNA sequence data: Implications for the molecular clock.

    Proc Natl Acad Sci USA 2003, 100:10335-10340. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  112. Loewe L, Charlesworth B, Bartolomé C, Nöel V: Estimating selection on non-synonymous mutations.

    Genetics 2006, 172:1079-1092. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  113. Charlesworth D, Charlesworth B, Morgan MT: The pattern of neutral molecular variation under the background selection model.

    Genetics 1995, 141:1619-1632. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  114. Charlesworth B, Morgan MT, Charlesworth D: The effect of deleterious mutations on neutral molecular variation.

    Genetics 1993, 134:1289-1303. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  115. Pollak E: On the theory of partially inbreeding finite populations. I. Partial selfing.

    Genetics 1987, 117:353-360. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  116. Nordborg M, Donnelly P: The coalescent process with selfing.

    Genetics 1997, 146:1185-1195. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  117. Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations.

    Genetics 1999, 151:1531-1545. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  118. Lynch M, O'Hely M, Walsh B, Force A: The probability of preservation of a newly arisen gene duplicate.

    Genetics 2001, 159:1789-1804. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  119. Lynch M, Katju V: The altered evolutionary trajectories of gene duplicates.

    Trends Genet 2004, 20:544-549. PubMed Abstract | Publisher Full Text OpenURL

  120. Lynch M, Force A: The probability of duplicate gene preservation by subfunctionalization.

    Genetics 2000, 154:459-473. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  121. Lynch M: The origins of genome architecture. Sunderland, Mass.: Sinauer Associates; 2007.

  122. Greenberg R, Crow JF: A Comparison of the Effect of Lethal and Detrimental Chromosomes from Drosophila Populations.

    Genetics 1960, 45:1153-1168. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  123. Gillespie JH: The heterozygous effects of alleles. In Population genetics: A concise guide. Edited by Gillespie JH. Baltimore, Maryland: John Hopkins University Press; 1998:62-71. OpenURL

  124. Simmons MJ, Crow JF: Mutations affecting fitness in Drosophila populations.

    Annu Rev Genet 1977, 11:49-78. PubMed Abstract | Publisher Full Text OpenURL

  125. Higgins K, Lynch M: Metapopulation extinction caused by mutation accumulation.

    Proc Natl Acad Sci USA 2001, 98:2928-2933. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  126. Johnston MO, Schoen DJ: Mutation rates and dominance levels of genes affecting total fitness in two angiosperm species.

    Science 1995, 267:226-229. PubMed Abstract | Publisher Full Text OpenURL

  127. Willis JH: Inbreeding load, average dominance and the mutation rate for mildly deleterious alleles in Mimulus guttatus.

    Genetics 1999, 153:1885-1898. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  128. Charlesworth B, Hughes KA: The maintenance of genetic variation in life-history traits. In Evolutionary genetics: from molecules to morphology. Edited by Singh RS, Krimbas CB. Cambridge, U.K.: Cambridge University Press; 2000:369-392. OpenURL

  129. Fry JD: How common are overdominant mutations?

    Genetics 2004, 167:1031-1032. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  130. Peters AD, Halligan DL, Whitlock MC, Keightley PD: Dominance and overdominance of mildly deleterious induced mutations for fitness traits in Caenorhabditis elegans.

    Genetics 2003, 165:589-599. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  131. Drake JW: A constant rate of spontaneous mutation in DNA-based microbes.

    Proc Natl Acad Sci USA 1991, 88:7160-7164. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  132. Drake JW, Charlesworth B, Charlesworth D, Crow JF: Rates of spontaneous mutation.

    Genetics 1998, 148:1667-1686. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  133. Pandian TJ, Koteeswaran R: Ploidy induction and sex control in fish.

    Hydrobiologia 1998, 384:167-243. OpenURL

  134. Comai L: The advantages and disadvantages of being polyploid.

    Nat Rev Genet 2005, 6:836-846. PubMed Abstract | Publisher Full Text OpenURL

  135. Balsano JS, Rasch EM, Monaco PJ: The evolutionary ecology of Poecilia formosa and its triploid associate. In Ecology and evolution of livebearing fishes (Poeciliidae). Edited by Meffe GK, Snelson FF. Englewood Cliffs, N.J.: Prentice Hall; 1989. OpenURL

  136. Pamilo P, Nei M, Li WH: Accumulation of mutations in sexual and asexual populations.

    Genet Res 1987, 49:135-146. PubMed Abstract OpenURL

  137. Gordo I, Charlesworth B: The degeneration of asexual haploid populations and the speed of Muller's ratchet.

    Genetics 2000, 154:1379-1387. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  138. Watterson GA: On the number of segregating sites in genetical models without recombination.

    Theor Popul Biol 1975, 7:256-276. PubMed Abstract | Publisher Full Text OpenURL

  139. Donaldson KA, Wilson RR Jr: Amphi-panamic geminates of snook (Percoidei: Centropomidae) provide a calibration of the divergence rate in the mitochondrial DNA control region of fishes.

    Mol Phylogenet Evol 1999, 13:208-213. PubMed Abstract | Publisher Full Text OpenURL

  140. Knowlton N, Weigt LA, Solorzano LA, Mills DK, Bermingham E: Divergence in proteins, mitochondrial-DNA, and reproductive compatibility across the Isthmus of Panama.

    Science 1993, 260:1629-1632. PubMed Abstract | Publisher Full Text OpenURL

  141. Howell N, Smejkal CB, Mackey DA, Chinnery PF, Turnbull DM, Herrnstadt C: The pedigree rate of sequence divergence in the human mitochondrial genome: there is a difference between phylogenetic and pedigree rates.

    Am J Hum Genet 2003, 72:659-670. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  142. Frankham R: Effective population size/adult population size ratios in wildlife – a review.

    Genet Res 1995, 66:95-107. OpenURL

  143. McVean GAT, Vieira J: Inferring parameters of mutation, selection and demography from patterns of synonymous site evolution in Drosophila.

    Genetics 2001, 157:245-257. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  144. Woodhead AD, Setlow RB, Pond V: The Amazon molly, Poecilia formosa, as a test animal in carcinogenicity studies: chronic exposures to physical agents.

    Natl Cancer Inst Monogr 1984, 65:45-52. PubMed Abstract OpenURL

  145. Saillard J, Forster P, Lynnerup N, Bandelt HJ, Norby S: mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion.

    Am J Hum Genet 2000, 67:718-726. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  146. Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision.

    Trends Genet 2004, 20:80-86. PubMed Abstract | Publisher Full Text OpenURL

  147. Lynch M, Blanchard J, Houle D, Kibota T, Schultz S, Vassilieva L, Willis J: Perspective: Spontaneous deleterious mutation.

    Evolution 1999, 53:645-663. OpenURL

  148. Baer CF, Miyamoto MM, Denver DR: Mutation rate variation in multicellular eukaryotes: causes and consequences.

    Nat Rev Genet 2007, 8:619-631. PubMed Abstract | Publisher Full Text OpenURL

  149. Lamatsch DK, Steinlein C, Schmid M, Schartl M: Non-invasive determination of genome size and ploidy level in fishes by flow cytometry: Detection of triploid Poecilia formosa.

    Cytometry 2000, 39:91-95. PubMed Abstract | Publisher Full Text OpenURL

  150. Misra S, Crosby MA, Mungall CJ, Matthews BB, Campbell KS, Hradecky P, Huang Y, Kaminker JS, Millburn GH, Prochnik SE, Smith CD, Tupy JL, Whitfied EJ, Bayraktaroglu L, Berman BP, Bettencourt BR, Celniker SE, de Grey AD, Drysdale RA, Harris NL, Richter J, Russo S, Schroeder AJ, Shu SQ, Stapleton M, Yamada C, Ashburner M, Gelbart WM, Rubin GM, Lewis SE: Annotation of the Drosophila melanogaster euchromatic genome: a systematic review.

    Genome Biol 2002, 3:research0083.0081-0022. OpenURL

  151. Gaffney DJ, Keightley PD: Genomic selective constraints in murid noncoding DNA.

    PLoS Genetics 2006, 2:e204. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  152. Halligan DL, Keightley PD: Ubiquitous selective constraints in the Drosophila genome revealed by a genome-wide interspecies comparison.

    Genome Res 2006, 16:875-884. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  153. Aparicio S, Chapman J, Stupka E, Putnam N, Chia JM, Dehal P, Christoffels A, Rash S, Hoon S, Smit A, Gelpke MD, Roach J, Oh T, Ho IY, Wong M, Detter C, Verhoef F, Predki P, Tay A, Lucas S, Richardson P, Smith SF, Clark MS, Edwards YJ, Doggett N, Zharkikh A, Tavtigian SV, Pruss D, Barnstead M, Evans C, Baden H, Powell J, Glusman G, Rowen L, Hood L, Tan YH, Elgar G, Hawkins T, Venkatesh B, Rokhsar D, Brenner S: Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes.

    Science 2002, 297:1301-1310. PubMed Abstract | Publisher Full Text OpenURL

  154. Haag-Liautard C, Dorris M, Maside X, Macaskill S, Halligan DL, Charlesworth B, Keightley PD: Direct estimation of per nucleotide and genomic deleterious mutation rates in Drosophila.

    Nature 2007, 445:82-85. PubMed Abstract | Publisher Full Text OpenURL

  155. Kibota TT, Lynch M: Estimate of the genomic mutation rate deleterious to overall fitness in E. coli.

    Nature 1996, 381:694-696. PubMed Abstract | Publisher Full Text OpenURL

  156. Kimura M: The neutral theory of molecular evolution. Cambridge: Cambridge University Press; 1983.