Friday, November 17, 2017

Calculating time of divergence using genome sequences and mutation rates (humans vs other apes)

There are several ways to report a mutation rate. You can state it as the number of mutations per base pair per year in which case a typical mutation rate for humans is about 5 × 10-10. Or you can express it as the number of mutations per base pair per generation (~1.5 × 10-8).

You can use the number of mutations per generation or per year if you are only discussing one species. In humans, for example, you can describe the mutation rate as 100 mutations per generation and just assume that everyone knows the number of base pairs (6.4 × 109).

The intrinsic mutation rate depends on the error rate of DNA replication. We don't know the exact value of this error rate but it's pretty close to 10-10 per base pair when you take repair into account [Estimating the Human Mutation Rate: Biochemical Method]. For single-cell species you simply multiply this number by the number of base pairs in the genome to get a good estimate of the mutation rate per generation.

The calculation for multicellular species is much more complicated because you have to know the number of cell divisions between zygote and mature germ cells. In some cases it's impossible to know this number (e.g.flowering plants, yeast). In other cases we have a pretty good estimate: for example, in humans there are about 400 cell divisions between zygote and mature sperm and about 30 cell divisions between zygote and mature egg cells. The number of cell divisions depends on the age of the parent, especially in males [Parental age and the human mutation rate]. This effect is significant—older parents pass on twice as many mutations as the youngest parents.

The parental age effect is comparable to the extremes in estimations of the human mutation rate based on different ways of measuring it [Human mutation rates - what's the right number?] [Human mutation rates ]. Those values range from about 60 mutations per generation to about 160 mutations per generation.

Thus, in the case of humans, we're dealing with estimates that differ by a factor of two depending on method and parental age.

Theme

Mutation

-definition
-mutation types
-mutation rates
-phylogeny
-controversies
Let's assume that each child is born with 100 new mutations. This seems like a reasonable number. It's on the high end of direct counts by sequencing parents and siblings but there are reasons to believe these counts are underestimated (Scally, 2016). On the other hand, this value (100 mutations) is on the low end of the estimates using the biochemical method and the phylogenetic method.

Most of these mutations occur in the father but some were contributed by the mother. Since the child is diploid, we calculate the mutation rate per bp as: 100 ÷ 6.4 × 109 = 1.56 × 10-8 per base pair per generation. Assuming an average generation time of 30 years, this gives 1.56 × 10-8 ÷ 30 = 5.2 × 10-10 mutations per bp per year. That's the value given above (rounded to 5 × 10-10). Scally (2016) uses this same value except he assumes a generation time of 29 years.

There are many who think this value is considerably lower than previous estimates and this casts doubt on the traditional times of divergence chimps and human and the other great apes. For example, Scally (2016) says that prior to the availability of direct sequencing date the "consensus value" was 10 × 10-10 per bp per year.1 That's twice the value he prefers today. It works out to 186 mutations per generation!

I think it's been a long time since workers in the field assumed such a high mutation rate but let's assume he is correct and current estimates are considerably lower than those from twenty years ago.

You can calculate a time of divergence (t) between any two species if you know the genetic distance (d) between them measured in base pairs and the mutation rate (μ) in mutations per year.2 The genetic distance can be estimated by comparing genome sequences and counting the differences. It represents the number of mutations that have become fixed in the two lineages since they shared a common ancestor. Haploid reference genome sequences are sufficient for this estimate.

The mutation rate (μ) is 100 mutations per generation divided by 30 years = 3.3 mutations per year.

The time of divergence is then calculated by dividing half that distance (in nucleotides) by the mutation rate (t = d/2 ÷ μ). (There are all kinds of "corrections" that can be applied to these values but let's ignore them for now and see what the crude data says.)

Human and chimp genomes differ by about 1.4%, which corresponds to 44.8 million nucleotide differences and d/2 = 22.4 million. Using 100 mutations per generation as the mutation rate means 5 × 10-10 per bp per year. From t = d/2 ÷ μ we get t = 6.8 million years. (Or, 22.4 ÷ 3.3 = 6.8.)

This is a reasonable number. It's consistent with the known fossil record and it's in line with the current views of a divergence time for chimps and humans.

However, there are reasons to believe that some of the assumptions in this calculation are wrong. For example, the average generation time is probably not 30 years in both lineages over the last few million years. It's probably shorter, at least in the chimp lineage where the current generation time is 25 years. Using a generation time of 25 years gives a divergence time of 5.6 million years.

In addition, the overall differences between the human and chimp genomes may be only 1.2% instead of 1.4% (see Moorjani et al., 2016). If you combine this value with the shorter generation time, you get 4.25 million years for the time of divergence.

Given the imprecision of the mutation rate, the question of real generation time, and problems in estimating the overall difference between humans and chimps, we can't know for certain what time of divergence is predicted by a molecular clock. On the other hand, the range of values (e.g. 4.25 - 6.8 million years) isn't cause for great concern.

So, what's the problem? The problem is that applying the human mutation rate (100 mutations per generation) to more distantly related species gives strange results. For example. Scally (2016) uses this mutation rate and a difference of 2.6% to estimate the time of divergence of humans and orangutans. The calculation yields a value of 26 million years. This is far too old according to the fossil record.

Several recent papers have addressed this issue (Scally, 2016; Moorjani et al., 2016a; Moorjani et al., 2016b). Most of the problem is solved by assuming a much higher mutation rate in the past. The biggest effect is the generation time in years. It may have been as low as 15 years for much of the past ten million years. Many of the problems go away when you adjust for this effect.

What puzzles me is the approach taken by Moorjani et al. in their two recent papers. They say that the "new" mutation rate is 5 × 10-10 per bp per year. That's exactly the value I use above. It's roughly 100 new mutations per child (per generation). Moorjani et al. (2016a) think this value is surprisingly low because it leads to a surprising result. They explain it in a section titled "The Puzzle."

They assume that the human and chimp genomes differ by 1.2%. That works out to 38 million mutations over the entire genome. This is 19 million fixed mutated alleles in each lineage if the mutation rate in both lineages is equal and constant.

If the mutation rate is 5 × 10-10 per bp per year then for a haploid genome this is 1.6 mutations per year. Dividing 19 million by 1.6 gives 11.9 million years (rounded to 12 million) for the time of divergence. This is the value quoted by the authors.
Taken at face value, this mutation rate suggests that African and non-African populations split over 100,000 years and a human-chimpanzee divergence time of 12 million years ago (Mya) (for a human–chimpanzee average nucleotide divergence of 1.2% at putatively neutral sites). These estimates are older than previously believed, but not necessarily at odds with the existing—and very limited—paleontological evidence for Homininae. More clearly problematic are the divergence times that are obtained for humans and orangutans or humans and OWMs [Old World Monkeys]. As an illustration, using whole genome divergence estimates for putatively neutral sites suggests a human–orangutan divergence time of 31 Mya and human–OWM divergence time of 62 Mya. These estimates are implausibly old, implying a human-oraguntan divergence well into the Oligocene and OWM-hominoid divergence well into or beyond the Eocene. Thus, the yearly mutation rates obtained from pedigrees seem to suggest dates that are too ancient to be readily reconciled with the current understanding of the fossil record.
Here's the problem. If the mutation rate is 100 mutations per generation then this applies to DIPLOID genomes. Some of the mutations are contribute by the mother and some (more) by the father. If you apply this rate to a DIPLOID genome then the number of mutations per year is 3.1 (100/30 years). Or,

         5 × 10-10 per bp per year × 6.4 × 109 bp (diploid) = 3.2 mutations per year

Dividing 19 million mutations by 3.2 give a time of divergence of 5.9 million years. This is a reasonable number but it's half the value calculated by Moorjani et al. (2016a).

They also calculate a value of 12.1 million years for the human-chimp divergence in their second paper (and 15.1 million years for the divergence of humans and gorillas) (Moorjani et al., 2016b).

I think their calculations are wrong because they used the haploid genome size rather than the diploid genome where the mutations are accumulating. Both these papers appear in good journals and both were peer-reviewed. Furthermore, the senior author, Molly Przeworski, is a Professor at Columbia University (New York, NY, USA) and she's an expert in this field.

What am I doing wrong? Is it true that a mutation rate of ~100 mutations per generation means that human and chimpanzees must have been separated for 12 million years as Moorjani et al. say? Or is the real value 5.9 million years as I've calculated above?


Image Credit: The chromosome image is from Wikipedia: Creative Commons Attribution 2.0 Generic license. The chimp photo is also from Wikipedia.

1. Scally takes this value from Nachman and Crowell (2000) who claim that the mutation rate is ~2.5 × 1008 mutations per bp in humans. This works out to 160 mutations per generation and an overall mutation rate of 8 × 10-10 based on a generation time of 30 years, not 10 × 10-10 as Scally states.

2. This assumes that all mutations are neutral. The rate of fixation of neutral alleles over time is equal to the mutation rate. Since 8% of the genome is under selection, it's not true that all mutations are neutral but to a first approximation it's not far off.

Moorjani, P., Gao, Z., and Przeworski, M. (2016) Human germline mutation and the erratic evolutionary clock. PLoS Biology, 14:e2000744. [doi: 10.1371/journal.pbio.2000744]

Moorjani, P., Amorim, C.E.G., Arndt, P.F., and Przeworski, M. (2016b) Variation in the molecular clock of primates. Proc. Nat. Acad. Sci. (USA) 113:10607-10612. [doi: 10.1073/pnas.1600374113 ]

Nachman, M.W., and Crowell, S.L. (2000) Estimate of the mutation rate per nucleotide in humans. Genetics, 156:297-304. [PDF]

Scally, A. (2016) The mutation rate in human evolution and demographic inference. Current opinion in genetics & development, 41:36-43. [doi: 10.1016/j.gde.2016.07.008]

52 comments:

  1. If you assume a mutation rate of 5e-10 mutations per base pair per year, and a human-chimp distance of 1.2% or 0.012, the divergence time is about 0.5 * 0.012 / 5e-10 = 12My. That is the puzzle, and it remains whatever convention you use (haploid vs diploid) for stating a number of mutations per generation.

    ReplyDelete
  2. The first term in your equation (0.5 * 5e-10) is half the percentage of the HAPLOID reference genomes that differ between humans and chimps.

    The second term in your equation is the mutation rate expressed as the number of mutations per DIPLOID genome per year. We know that the mutation rate is about 100 mutations per generation in a DIPLOID genome. The mutation rate in your equation is based on the number of mutations per year per 6.4e9 base pairs.

    5e-10 X 6.4e9 X 30 years = 96 mutations/year

    ReplyDelete
    Replies
    1. We agree that the first term in your equation is 19 million ...

      0.5 X 0.012 X haploid genome size = 19 million

      The second term in your equation is either 1.6 mutations per year or 3.2 mutations per year depending on who does the calculation.

      The data says about 100 mutations per generation so the second value is the correct one.

      Delete
    2. i think here is another interesting case: according to evolution fly and mosquito split off about 250 my ago. fly generation is about less then one month. so even if one generation mean only 1 new mutation (i think in reality they get much more) we will need only 10^8 month to change their entire genome. or about less than 10^7 years. so fly and mosquito should be different in about their entire genomes from each other (or at least 75% difference). far from reality as far as im aware about.

      Delete
    3. https://en.wikipedia.org/wiki/Conserved_sequence

      Delete
    4. are you saying that the entire genome is conserve?

      Delete
    5. You're not accounting for saturation effects, either. If there's a 10% difference between genomes, then there's a 0.1 chance of a mutation occuring where the genomes are already divergent. I'm also not sure where you get the 205ma figure from. Misof et al. 2014 has the split at ~160Ma and the refined value we got (preprint: https://www.biorxiv.org/content/early/2017/04/18/128314) is at ~120Ma. That figure is half of the divergence time you assume, which is closer to the Diptera - (Mecoptera+Siphoneraptera) split in both Misof et al. and the paper in prep.

      Delete
    6. And then there's the fact that mutations don't happen uniformly across the genome either. The chance of mutation is not equiprobable across all nucleotides.

      Delete
    7. "2014 has the split at ~160Ma and the refined value we got"-

      but it doesnt matter at all. even 100my is too much. i just assume about 100 new mutations per generation. and this is the number you will get. so something is very weird here.


      "And then there's the fact that mutations don't happen uniformly across the genome either."-

      but we are talking about the entire genome here. are you saying that the majority of the genome doesnt get any mutations?

      Delete
    8. All the genome gets mutations. However, mutations happen more often at certain "hotspots" than elsewhere. (Why? I don't know.)

      Only part of the genome is conserved. DNA sequences for genes, regulatory regions, telomeres, and other useful functions are conserved. Mutations happen there, but most of the mutations are selected out because they do harm. In the "junk" regions of the DNA (about 90% of the DNA in humans, sequence is not conserved.

      Delete
    9. What regulates the "hot spots"?
      What regulates mutations in overlapping genes?

      Delete
  3. The amount of difference between human and chimp DNA should be the product of the mutation rate per year with the sum of the divergence time plus the coalescence time of two gene copies, times 2. One must be careful not to leave out the coalescence time, which is a number of generations equal to twice the effective population size.

    ReplyDelete
  4. This is from the PNAS paper by Moorjani et al. I've lost track of what the puzzle is.

    "we estimate that humans diverged from chim-
    panzees ~12.1 Mya and from gorillas ~15.1 Mya (Fig. 4). Assuming
    further that the effective population size of the human–ape an-
    cestor was five times the current population size (as estimated by
    refs. 43, 44), the human–chimpanzee split time is ~7.9 Mya, and
    the human–gorilla split time is 10.8 Mya. We note that there is
    substantial uncertainty in estimates of ancestral population size of
    apes, with previous estimates ranging between 50,000 and 100,000
    (43–45). Accounting for this uncertainty provides estimates of
    human–chimpanzee split time in the range of 6.5–9.3 Mya, and
    human–gorilla split time in the range of 9.4–12.2 Mya."

    ReplyDelete
  5. It sounds as if they have taken coalescence effects properly into account.

    ReplyDelete
  6. Fossil record should not be a factor in figuring genetic operations.! this must stop about studying rocks and using it to confirm biology hypothesis.
    Anyways.
    This divergence concept is entirely based on observing modern mutation rates which do not produce much mutation in people.
    this idea allows for no other options.
    If a creator created biology in finished form and hen mutations were going on in trivial ways IT WOULD be a math error, or a concept error, to extrap[olate back from this mutation rate.
    mutation rates do not demonstrate man/ape divergence but only show a divergence score, AFTER A PRESUMPTION, and so they show nothing of genetic evidence for common descent.
    its still all reinforcing presumptions without scientific methodology.
    you can't falsify any of these rates conclusions unless already agreeing that these mutations have been going on.

    ReplyDelete
    Replies
    1. Robert, you ever hear the phrase from carpentry "Measure twice, cut once"? The biology and the rocks are measuring twice, that's all.

      Delete
    2. But you presume the existence of a creator god, despite the utter nonexistence of any testable evidence to that effect.

      Mote, beam, Robert Byers.

      Delete
  7. Going back to the beginning of the post, the difference between the mutation rates calculated by biochemical measurements and the direct counts from parent/sibling sequences seems to me likely to be a product of differential mortality. Parents and siblings are survivor, no, therefore a biased sample. (Yes, direct counts from non-survivors are impossible.) For the purposes of estimating time of evolutionary divergence, though, it appears that the direct counts are the best estimate.

    The thing is, if there is a period in the evolutionary history when the constraints of natural selection are slackened, so that variations in phenotype that previously suffered differential mortality instead prosper, the mutation rate as measured in direct counts would be higher. And it seems to me that a population moving into a new area could indeed tolerate a higher effective mutation rate (as I suppose you could call it.) Obviously this is just a version of radiative adaptation.

    Doesn't it seem as though these divergences in estimates are really suggesting a period of radiative adaptation in those lineages?

    ReplyDelete
    Replies
    1. There are plenty of ways to modify the mutation rate in order to account for various effect that the crude values ignore. There are several ways of adjusting the formula for calculating times of divergence. One of them involves taking coalescence and population size into account as Joe pointed out above.

      These are the main points of the papers I referenced but it's not the main point of my post. I'm asking whether Moorjani et al. have done the arithmetic correctly.

      Delete
  8. Why do you have negative signs in front of you exponents? e.g.

    "In humans, for example, you can describe the mutation rate as 100 mutations per generation and just assume that everyone knows the number of base pairs (6.4 × 10-9). "

    I'm pretty sure you mean 10 to the 9th power or billions here, right?

    Also this page says, "The human genome contains approximately 3 billion of these base pairs..."
    https://www.genome.gov/11006943/human-genome-project-completion-frequently-asked-questions/

    Your quote above says the number of base pairs is 6.4 billion. What are you referring to?

    ReplyDelete
    Replies
    1. Thank-you for spotting the typo. I fixed it.

      The current version of the human reference genome (GRCh38.p11) contains 3,092,480,053 base pairs. There are still quite a few gaps but the approximate size of most gaps is known. The total genome size is estimated to be 3,253,848,404 bp.

      Diploid cells (e.g. human zygote) contain twice this amount of DNA.

      Delete
  9. Larry Moran-
    From the paper it appears Moorjani et al. does use the haploid rather than diploid genome size.
    Your arithmetic is correct.
    It seems you have solved ‘the puzzle’.

    ReplyDelete
  10. Funnily enough, we use the human mutation rate and haploid genome size in an exercise for first year undergraduates in our Evolution 101 course. The aim is to calculate the number of new mutations per mitosis. None of the students ever accounts for the genome being diploid.

    ReplyDelete
  11. There's currently a discussion going on between Dr Richard Buggs (Reader in Evolutionary Genomics at Queen Mary, University of London), and Dennis Venema (Professor of biology at Trinity Western University), about whether or not there has been a human bottleneck down to only two people in the last 500,000 years. Dr Buggs is arguing that there is no evidence against such a bottleneck, and that none of the extant population studies have ruled it out.

    http://www.richardbuggs.com/email-to-dennis-venema-adam-eve-genomics.html

    ReplyDelete
  12. Important issue there are:
    1. Really? "No evidence against?" MHC (H:A) diversity and shared alleles with chimpanzees?
    2. "Ruled it out". How "ruled out" do you need it to be to be "ruled out"? Science rarely, absolutely certainly, rules anything in or out.
    Dennis Venema has properly qualified his statements on these issues. Furthermore you, Furtigern, have put words in Buggs's mouth -- the letter by Buggs that you link to does not say "no evidence against".

    ReplyDelete
    Replies
    1. "1. Really? "No evidence against?" MHC (H:A) diversity and shared alleles with chimpanzees?"

      Dr Buggs claims neither of those provide evidence against a bottleneck of just two people.

      "2. "Ruled it out". How "ruled out" do you need it to be to be "ruled out"? Science rarely, absolutely certainly, rules anything in or out."

      Agreed. Instead we have "persistently resistant to falsfiability". There does come a point at which that resistance seems futile to deny, which is why mainstream science doesn't say the sun goes around the earth. Heliocentrism has proved fairly consistently resistant to falsification. At this point geocentrism is unlikely to make a comeback, even if we argue it hasn't yet been "ruled out".

      "Dennis Venema has properly qualified his statements on these issues."

      Yes. And Dr Buggs has continued to resist his conclusions.

      "Furthermore you, Furtigern, have put words in Buggs's mouth -- the letter by Buggs that you link to does not say "no evidence against"."

      I've read every post he had made in the discussion on the subject (which is over 30 so far), and no I don't think I am putting words in his mouth. You're only looking at a single blog post.

      Dennis says that none of the evidence even remotely suggests a bottleneck down to two people. Buggs replies thus.

      * "I think we are getting ahead of ourselves if we say that none of the evidence even remotely suggests it, given that the hypothesis has not been directly tested"

      He clearly disagrees, and claims the hypothesis has not even been tested. Later in the same post he says "they never set out to test the bottleneck hypothesis" and "the scientific literature does not specifically address the question of whether or not humans have passed through a bottleneck of two".

      Read these comments of his and tell me if you think he's arguing that there's evidence against a bottleneck of two.

      * "This is fine, so long as mainstream science really is showing that a single-pair bottleneck somewhere in the history of humankind is an absolute impossibility. But having looked at the evidence that Venema describes and cites, I am not yet convinced."

      * "I realise I am acting as a defense attorney in a case that seems to you to be a bit hopeless,"

      * "I would also note that thus far, no mainstream scientist has disputed my blog which has been read by many scientists."

      * "I think that the Zhao et al (2000) data is explainable from a bottleneck of two within the last few hundred thousand years, taking their own coalescence analysis as a starting point"

      * "Does genomic evidence make it scientifically impossible that the human lineage could have ever passed through a population bottleneck of just two individuals? This is a question I am asked semi-frequently by religious friends. With my current understanding of the genetic evidence, I can’t state categorically that it’s impossible."

      * "this hypothesis still needs testing and has not been explicitly and reliably tested in an already-published analysis"

      * " I therefore cannot conclude from the models that you have run that a bottleneck of two in the history of the human lineage is not possible"

      * "I have pointed out to you that neither of the parts you highlighted are good evidence against the bottleneck of two hypothesis"

      * "I struggle to see this paper as evidence that there was never a short sharp bottleneck in human history."

      * "Please could you explain to me why you think that data excludes a bottleneck of two at some point in the human lineage?"

      * "This would appear to leave open the possibility of a bottleneck in the previous 5.8 million years."

      Delete
    2. Dr Buggs claims neither of those provide evidence against a bottleneck of just two people.

      What is the basis of Dr. Buggs's interesting contention?

      Delete
    3. Ostensibly science, but ultimately theology. He is a proponent of Intelligent Design, which he describes as a science. Predictably, the Discovery Institute is already appealing to this discussion as supposed evidence that every human who has ever lived, descended from Adam and Eve.

      https://uncommondescent.com/genetics/adam-eve-richard-buggs-and-dennis-venema-could-adam-and-eve-have-existed/

      Delete
    4. It seems to me that it is irrelevant whether an observation was made specifically to test the A&E hypothesis. If it ends up having relevant information we cannot casually dismiss that information.

      And there is no evidence, as far as I know, that specifically points to a bottleneck of two people. (Leaving aside whether we should be asking about a bottleneck that is effectively of one person, Eve being a clone of Adam).

      Delete
    5. Both of those points have been put to Dr Buggs several times. He insists that since no test has been made of a specific bottleneck of two people, the hypothesis has not yet been tested. On the matter of Eve being taken from Adam, he has not commented (though this has also been raised).

      Delete
    6. That is like a flat-earther saying that the roundness of the shadow of the earth on the moon is not evidence of the earth being round, because when people observed it, they were not specifically testing that hypothesis.

      Delete
    7. Yes, but what about the existence of more than four HLA alleles shared between humans and chimps? Doesn't that preclude a bottleneck of 2? What is Buggs's argument to the contrary?

      Delete
    8. He hasn't addressed that. It's one of a number of points he hasn't addressed.

      Delete
    9. And even though he hasn't addressed it, he claims that it doesn't preclude a bottleneck of two?

      Delete
    10. (Based on the discussion in this thread) it sounds as if he dismisses the evidence from HLA alleles because it was not collected with the specific intent of testing the existence of A&E. That is facepalm territory -- please, someone tell me he isn't saying that.

      Delete
    11. Hi guys, I know this is not really pertinent to this particular post, but I can't think of a better place to ask this.

      The Discovery Institute is bragging about an article that has been published in a legit journal by John C Sanford (of "Genetic Entropy" infamy) and someone named William F. Basener. They claim to have falsified Fisher's Theorem. The math is way beyond me, but as best as I can tell they calculated the results assuming that deleterious mutations outnumber beneficial mutations by 1000:1 and determine that only extinction results. However, as far as I can tell, they neglected to include neutral/nearly neutral mutations in the model, which would obviously make the results meaningless. Am I right in suspecting that is what they have done?

      Here is the link:

      https://link.springer.com/article/10.1007/s00285-017-1190-x

      Delete
    12. lutesuite-
      The paper you linked to gives an updated version of Fisher’s model that includes mutations. (The new model will give Fisher’s if one sets the effects of mutations to zero).
      The new model is then used to test what will happen given what are called ‘realistic’ assumptions about mutations- including deleterious mutations, nearly- neutral mutations and beneficial mutations.
      They use values from Kimura for the deleterious and nearly- neutral mutations. That is covered in section 5 of the paper.
      The model runs do not show the ‘mean fitness’ increasing, as Fisher suggested, but rather that the ‘mean fitness’ likely declines overtime.

      I think it observably true that species are more likely to go extinct than to continue- this result would be in keeping with that observation, so I’m not sure if there is any problem here.

      Delete
    13. Lutesuite,

      I had a quick look at the paper. I do not think there's a serious mistake in their maths or simulations, and they certainly have not ignored neutral/nearly neutral mutations. Their main point is about very slightly deleterious mutations (VSDMs), too small in effect to be removed by selection (so usually called neutral/nearly neutral), but possibly important because there could be vast numbers of them.

      When a new mutation arises, it has some selection coefficient x. The controversy is about what x is likely to be. What is the distribution of x? We can only measure x very approximately. Every species (and the same species in different situations) will have a different answer to this fundamental question: What is the distribution of x? Search for 'Distribution of fitness effects' for more information. The basic problem is that if x is close to zero, we can't tell if it is very slightly deleterious or very slightly beneficial. We do know that big positive values of x (say bigger than 0.01) are way more likely (1000x or 1000000x more likely) than small values (less than -0.01), but that doesn't tell us about the distribution of x very close to zero.

      You can choose many distributions for x which don't contradict the available evidence. Sanford and Basener chose one. Their distribution is assymetric by a factor of 1000 regardless of how far x is from zero (see their Fig 1). I think (but of course do not know) that it should be almost symmetric near zero, and increasingly assymetric further away.


      Delete
    14. Graham,
      The authors say they overestimate the probability of beneficial mutations (.0001).
      They also state-
      “By assuming the same parameterization for beneficial and deleterious mutations, we also overestimate the range of the beneficial mutations because beneficial mutations tend to consistently have more modest fitness effects.”

      Are you thinking this statement is in error?
      Would this be due to beneficial mutations having an exponential distribution?

      Delete
    15. It seems to me impossible that the fraction of beneficial to deleterious mutations could ever be 1:1000 unless organisms were already extremely well adapted. Because every deleterious mutation would imply a beneficial one if reversed.

      Delete
    16. @Mikkel Rumraket Rasmussen: You noted that "every deleterious mutation would imply a beneficial one if reversed".

      @Mikkel Rumraket Rasmussed: You noted that "every deleterious mutation would imply a beneficial one if reversed".

      An important point. Whatever it implies about the ratio of beneficial to deleterious mutations, it also completely demolishes another common creationist argument. They claim that there is some mysterious informational / theological rule that there is no such thing as a beneficial mutation. Which puts their whole worldview in a perilous state if a beneficial mutation happens to be observed. Your point, which I've also made myself also here at Sandwalk, neatly applies a needle to that balloon.

      Delete
    17. Readers, kindly ignore the duplicate opening sentence of the preceding comment -- it was an editing mistake which I cannot reverse now.

      Delete
    18. Well, they did say the presumed ratio of 1000:1 was "generous." Maybe they actually think it's infinity.

      Delete
    19. Jack Jackson,

      I think they overestimate the probability of beneficial mutations among those of large effect (and know this and say so), but underestimate the probability of beneficial mutations among those of small effect.

      In the sentence you quoted, I don't know what by they mean by 'range'. I think it is likely that they were trying to say something I'd agree with, but I can't really tell, it is too imprecise.

      Delete
    20. Graham-
      Thank-you for that. Deciding the issue observationally seems to be very difficult when mutations have little effect.
      I think the same mutation can be beneficial or deleterious depending on the situation- making the adjudication even more perplexing.
      It would be interesting to see how much a change like you mention would effect the overall result mathematically.

      Delete
    21. So, I finally had time to look at the paper and the first thing I checked was if Sanford was repeating his old mistake (conflating relative fitness and absolute fitness). And (of course) I find precisely that. It's a fundamental error to assume that if an allele with s<0 is fixed, this implies that the mean absolute fitness of the population has decreased. s<0 means that the individuals carrying the allele have a lower mean fitness than those that do not, *in the presence of those that do not*. Most of what s measures is competition within a population for some resource. And that means that if a detrimental allele gets fixed, it also means the better competitors are gone. In the medium to long term the means absolute fitness of every population is 1.

      Delete
    22. Michael Lynch and I have recently posted at The Skeptical Zone a more comprehensive reply to the paper by Basener and Sanford. We argue that the omission of recombination from their model, together with the assumption of a particular distribution of net whole genome effects of mutation, create a situation where mutation will typically overwhelm selection. If there is recombination, and if mutation effects at different loci are relatively independent, then the conventional result of mutation-selection balance will be expected.

      Delete
  13. This comment has been removed by the author.

    ReplyDelete
  14. Dear Professor Moran, I hope you are in good health and your absence from your blog is because you are enjoying yourself elsewhere. Your regular posts are missed. Looking forward to your next posts. Wishing you a happy new year.

    ReplyDelete
  15. This comment has been removed by the author.

    ReplyDelete