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.
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.