Molecular clock

From Wikipedia, the free encyclopedia
Jump to: navigation, search

The molecular clock is a technique that uses the mutation rate of biomolecules to deduce the time in prehistory when two or more life forms diverged. The biomolecular data used for such calculations are usually nucleotide sequences for DNA or amino acid sequences for proteins. The benchmarks for determining the mutation rate are often fossil or archaeological dates. The molecular clock was first tested in 1962 on the haemoglobin protein variants of various animals, and is commonly used in molecular evolution to estimate times of speciation or radiation. It is sometimes called a gene clock or an evolutionary clock.

Early discovery and genetic equidistance[edit]

The notion of the existence of a so-called "molecular clock" was first attributed to Emile Zuckerkandl and Linus Pauling who, in 1962, noticed that the number of amino acid differences in hemoglobin between different lineages changes roughly linearly with time, as estimated from fossil evidence.[1] They generalized this observation to assert that the rate of evolutionary change of any specified protein was approximately constant over time and over different lineages (based on the molecular clock hypothesis (MCH)).

The genetic equidistance phenomenon was first noted in 1963 by Emanuel Margoliash, who wrote: "It appears that the number of residue differences between cytochrome c of any two species is mostly conditioned by the time elapsed since the lines of evolution leading to these two species originally diverged. If this is correct, the cytochrome c of all mammals should be equally different from the cytochrome c of all birds. Since fish diverges from the main stem of vertebrate evolution earlier than either birds or mammals, the cytochrome c of both mammals and birds should be equally different from the cytochrome c of fish. Similarly, all vertebrate cytochrome c should be equally different from the yeast protein."[2] For example, the difference between the cytochrome c of a carp and a frog, turtle, chicken, rabbit, and horse is a very constant 13% to 14%. Similarly, the difference between the cytochrome c of a bacterium and yeast, wheat, moth, tuna, pigeon, and horse ranges from 64% to 69%. Together with the work of Emile Zuckerkandl and Linus Pauling, the genetic equidistance result directly led to the formal postulation of the molecular clock hypothesis in the early 1960s.[3] Genetic equidistance has often been used to infer equal time of separation of different sister species from an outgroup.[4][5]

Later Allan Wilson and Vincent Sarich built upon this work.

Relationship with neutral theory[edit]

The observation of a clock-like rate of molecular change was originally purely phenomenological. Later, the work of Motoo Kimura[6] developed the neutral theory of molecular evolution, which predicted a molecular clock. Let there be N individuals, and to keep this calculation simple, let the individuals be haploid (i.e. have one copy of each gene). Let the rate of neutral mutations (i.e. mutations with no effect on fitness) in a new individual be . The probability that this new mutation will become fixed in the population is then 1/N, since each copy of the gene is as good as any other. Every generation, each individual can have new mutations, so there are N new neutral mutations in the population as a whole. That means that each generation, new neutral mutations will become fixed. If most changes seen during molecular evolution are neutral, then fixations in a population will accumulate at a clock-rate that is equal to the rate of neutral mutations in an individual.

Calibration[edit]

The molecular clock alone can only say that one time period is twice as long as another: it cannot assign concrete dates. To achieve this, the molecular clock must first be calibrated against independent evidence about dates, such as the fossil record.[7] For viral phylogenetics and ancient DNA studies – two areas of evolutionary biology where it is possible to sample sequences over an evolutionary timescale – the dates of the intermediate samples can be used to more precisely calibrate the molecular clock.

Non-constant rate of molecular clock[edit]

Sometimes only a single divergence date can be estimated from fossils, with all other dates inferred from that. Other sets of species have abundant fossils available, allowing the MCH of constant divergence rates to be tested. DNA sequences experiencing low levels of negative selection showed divergence rates of 0.7–0.8% per Myr in bacteria, mammals, invertebrates, and plants.[8] In the same study, genomic regions experiencing very high negative or purifying selection (encoding rRNA) were considerably slower (1% per 50 Myr).

In addition to such variation in rate with genomic position, since the early 1990s, variation among taxa has proven fertile ground for research too,[9] even over comparatively short periods of evolutionary time (for example mockingbirds[10]). Tube-nosed seabirds have molecular clocks that on average run at half speed of many other birds,[11] possibly due to long generation times, and many turtles have a molecular clock running at one-eighth the speed it does in small mammals or even slower.[12] Effects of small population size are also likely to confound molecular clock analyses. Researchers such as Francisco Ayala have more fundamentally challenged the molecular clock hypothesis.[13][14] According to Ayala's 1999 study, five factors combine to limit the application of molecular clock models:

  • Changing generation times (If the rate of new mutations depends at least partly on the number of generations rather than the number of years)
  • Population size (Genetic drift is stronger in small populations, and so more mutations are effectively neutral)
  • Species-specific differences (due to differing metabolism, ecology, evolutionary history,...)
  • Change in function of the protein studied (can be avoided in closely related species by utilizing non-coding DNA sequences or emphasizing silent mutations)
  • Changes in the intensity of natural selection.
Phylogram showing three groups, one of which has strikingly longer branches than the two others
Woody bamboos (tribes Arundinarieae and Bambuseae) have long generation times and lower mutation rates, as expressed by short branches in the phylogenetic tree, than the fast-evolving herbaceous bamboos (Olyreae).

Molecular clock users have developed workaround solutions using a number of statistical approaches including maximum likelihood techniques and later Bayesian modeling. In particular, models that take into account rate variation across lineages have been proposed in order to obtain better estimates of divergence times. These models are called relaxed molecular clocks[15] because they represent an intermediate position between the 'strict' molecular clock hypothesis and Joseph Felsenstein's many-rates model[16] and are made possible through MCMC techniques that explore a weighted range of tree topologies and simultaneously estimate parameters of the chosen substitution model. It must be remembered that divergence dates inferred using a molecular clock are based on statistical inference and not on direct evidence.

The molecular clock runs into particular challenges at very short and very long timescales. At long timescales, the problem is saturation. When enough time has passed, many sites have undergone more than one change, but it is impossible to detect more than one. This means that the observed number of changes is no longer linear with time, but instead flattens out.

At very short time scales, many differences between samples do not represent fixation of different sequences in the different populations. Instead, they represent alternative alleles that were both present as part of a polymorphism in the common ancestor. The inclusion of differences that have not yet become fixed leads to a potentially dramatic inflation of the apparent rate of the molecular clock at very short timescales.[17][18]

Methods[edit]

Bayesian methods can provide more appropriate estimates of divergence times, especially if large datasets—such as those yielded by phylogenomics—are employed.[19]

Uses[edit]

The molecular clock technique is an important tool in molecular systematics, the use of molecular genetics information to determine the correct scientific classification of organisms or to study variation in selective forces.

Knowledge of approximately constant rate of molecular evolution in particular sets of lineages also facilitates establishing the dates of phylogenetic events, including those not documented by fossils, such as the divergence of living taxa and the formation of the phylogenetic tree. In these cases—especially over long stretches of time—the limitations of MCH (above) must be considered; such estimates may be off by 50% or more.

See also[edit]

References[edit]

  1. ^ Zuckerkandl, E. and Pauling, L.B. (1962). "Molecular disease, evolution, and genic heterogeneity". In Kasha, M. and Pullman, B (editors). Horizons in Biochemistry. Academic Press, New York. pp. 189–225. 
  2. ^ Margoliash E (October 1963). "PRIMARY STRUCTURE AND EVOLUTION OF CYTOCHROME C". Proc. Natl. Acad. Sci. U.S.A. 50 (4): 672–9. doi:10.1073/pnas.50.4.672. PMC 221244free to read. PMID 14077496. 
  3. ^ Kumar S (August 2005). "Molecular clocks: four decades of evolution". Nat. Rev. Genet. 6 (8): 654–62. doi:10.1038/nrg1659. PMID 16136655. 
  4. ^ Pesole G, Gissi C, De Chirico A, Saccone C (April 1999). "Nucleotide substitution rate of mammalian mitochondrial genomes". J. Mol. Evol. 48 (4): 427–34. doi:10.1007/PL00006487. PMID 10079281. 
  5. ^ Huang, S. (2008) The genetic equidistance result of molecular evolution is independent of mutation rates. J. Comp. Sci. Syst. Biol., 1: 92-102. http://omicsonline.com/ArchiveJCSB/Ab01/JCSB1.092.html
  6. ^ Kimura, Motoo (1968). "Evolutionary rate at the molecular level". Nature. 217 (5129): 624–626. doi:10.1038/217624a0. PMID 5637732. 
  7. ^ Benton, M. J. & Donoghue, P. C. J. (2007). "Paleontological evidence to date the Tree of Life". Molecular Biology & Evolution. 24 (1): 26–53. doi:10.1093/molbev/msl150. PMID 17047029. 
  8. ^ Ochman H, Wilson AC (1987). "Evolution in bacteria: evidence for a universal substitution rate in cellular genomes". J Mol Evol. 26 (1–2): 74–86. doi:10.1007/BF02111283. PMID 3125340. 
  9. ^ Douzery, E.J.P., Delsuc, F., Stanhope, M.J. and Huchon, D. (2003). "Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals, and incompatibility among fossil calibrations". Journal of Molecular Evolution. 57: S201–S213. doi:10.1007/s00239-003-0028-x. PMID 15008417. 
  10. ^ Hunt, J.S., Bermingham, E., and Ricklefs, R.E. (2001). "Molecular systematics and biogeography of Antillean thrashers, tremblers, and mockingbirds (Aves: Mimidae)". Auk. 118 (1): 35–55. doi:10.1642/0004-8038(2001)118[0035:MSABOA]2.0.CO;2. ISSN 0004-8038. 
  11. ^ Rheindt, F. E. & Austin, J. (2005). "Major analytical and conceptual shortcomings in a recent taxonomic revision of the Procellariiformes - A reply to Penhallurick and Wink (2004)" (PDF). Emu. 105 (2): 181–186. doi:10.1071/MU04039. 
  12. ^ Avise, J.C., Bowen, W., Lamb, T., Meylan, A.B. and Bermingham, E. (1 May 1992). "Mitochondrial DNA Evolution at a Turtle's Pace: Evidence for Low Genetic Variability and Reduced Microevolutionary Rate in the Testudines". Molecular Biology and Evolution. 9 (3): 457–473. PMID 1584014. 
  13. ^ Ayala, F.J. (1999). "Molecular clock mirages". BioEssays. 21 (1): 71–75. doi:10.1002/(SICI)1521-1878(199901)21:1<71::AID-BIES9>3.0.CO;2-B. PMID 10070256. 
  14. ^ Schwartz, J. H. & Maresca, B. (2006). "Do Molecular Clocks Run at All? A Critique of Molecular Systematics". Biological Theory. 1 (4): 357–371. doi:10.1162/biot.2006.1.4.357. Lay summaryScience Daily. 
  15. ^ Drummond, A.J., Ho, S.Y.W., Phillips, M.J. and Rambaut A. (2006). "Relaxed Phylogenetics and Dating with Confidence". PLoS Biology. 4 (5): e88. doi:10.1371/journal.pbio.0040088. PMC 1395354free to read. PMID 16683862.  open access publication - free to read
  16. ^ "Taking variation of evolutionary rates between sites into account in inferring phylogenies.". J Mol Evol. 53: 447–55. doi:10.1007/s002390010234. PMID 11675604. 
  17. ^ Ho SY, Phillips MJ, Cooper A, Drummond AJ (2005). "Time dependency of molecular rate estimates and systematic overestimation of recent divergence times". Molecular Biology & Evolution. 22 (7): 1561–1568. doi:10.1093/molbev/msi145. PMID 15814826. 
  18. ^ Peterson GI, Masel J (2009). "Quantitative Prediction of Molecular Clock and Ka/Ks at Short Timescales". Molecular Biology & Evolution. 26 (11): 2595–2603. doi:10.1093/molbev/msp175. PMC 2912466free to read. PMID 19661199. 
  19. ^ Dos Reis, M.; Inoue, J.; Hasegawa, M.; Asher, R. J.; Donoghue, P. C. J.; Yang, Z. (2012). "Phylogenomic datasets provide both precision and accuracy in estimating the timescale of placental mammal phylogeny". Proceedings of the Royal Society B: Biological Sciences. 279 (1742): 3491–3500. doi:10.1098/rspb.2012.0683. 

Further reading[edit]

External links[edit]