This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Mutation Rate In Human And Calculating Time From Diverged Sequences.

I have two sequences which used to be identical. During time they experienced several mutations, insertions or deletions. I would like to count time that passed.

I read about:

T=k/2r, where k=divergence, r=substitution rate, T=time

However, I read very ambiguous information about human mutation rate (I am interested mostly in mutation rate in retrotransposons) and about the calculation formula (there are many and many).

I would like to start with something really simple at the beginning.

Could you please advice me how to do it, supposing I have two sequences and substitution rate? Or some literature considering this topic? In many articles there stands just "we calculated" instead of exact formula.

Thanks a lot.

mutation

1 answer

Under the neutral theory, the basic divergence formula is d = 2uT (divergence is equal to twice the mutation rate times time)

The factor of 2 comes from the fact that when comparing 2 sequences that diverged time T in the past, there are 2 branches on which independent evolutionary changes can occur.

    X       |
   X X      |
  X   X     | T
 X     X    |
X       X   |

Rearranging this expression gives an estimate of the divergence time:

T = d/2u

This assumes a molecular clock, which effectively implies neutral sequence divergence and a constant mutation rate u. The reason that no explicit method is given is because this is (for a point estimates) a simple calculation.

Mutation rates can be estimated directly via mutation accumulation experiments or by measuring new mutations in parent-offspring comparisons. A good review on the topic can be found here: http://rstb.royalsocietypublishing.org/content/365/1544/1169.full

If you are interested dating the age of insertions using intra-element LTR-LTR comparisons, I suggest you look into using REannotate, which can perform the calculation of d (its standard deviation) automatically: http://www.bioinformatics.org/reannotate/

Log in to answer this question.