MBE Advance Access originally published online on September 21, 2007
Molecular Biology and Evolution 2007 24(12):2669-2680; doi:10.1093/molbev/msm193
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
A General Comparison of Relaxed Molecular Clock Models
* Department of Mathematics and Statistics, McGill University, Montréal, Québec, Canada
Department of Mathematics, University of Auckland, New Zealand
Département de biochimie, Université de Montréal, Québec, Canada
Laboratoire d'Informatique, de Robotique et de Microélectronique de Montpellier, UMR 5506, CNRS-Université de Montpellier 2, 161, rue Ada, 34392 Montpellier Cedex 5, France
E-mail: nicolas.lartillot{at}lirmm.fr.
Accepted for publication September 7, 2007.
Several models have been proposed to relax the molecular clock in order to estimate divergence times. However, it is unclear which model has the best fit to real data and should therefore be used to perform molecular dating. In particular, we do not know whether rate autocorrelation should be considered or which prior on divergence times should be used. In this work, we propose a general bench mark of alternative relaxed clock models. We have reimplemented most of the already existing models, including the popular lognormal model, as well as various prior choices for divergence times (birth–death, Dirichlet, uniform), in a common Bayesian statistical framework. We also propose a new autocorrelated model, called the "CIR" process, with well-defined stationary properties. We assess the relative fitness of these models and priors, when applied to 3 different protein data sets from eukaryotes, vertebrates, and mammals, by computing Bayes factors using a numerical method called thermodynamic integration. We find that the 2 autocorrelated models, CIR and lognormal, have a similar fit and clearly outperform uncorrelated models on all 3 data sets. In contrast, the optimal choice for the divergence time prior is more dependent on the data investigated. Altogether, our results provide useful guidelines for model choice in the field of molecular dating while opening the way to more extensive model comparisons.
Key Words: relaxed clock • Bayes factor • molecular dating • CIR process • phylogeny • Markov chain Monte Carlo
Spencer Muse, Associate Editor
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  J. Inoue, P. C. J. Donoghue, and Z. Yang The Impact of the Representation of Fossil Calibrations on Bayesian Estimation of Species Divergence Times Syst Biol, January 1, 2010; 59(1): 74 - 89. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. P. Brown and Z. Yang Bayesian Dating of Shallow Phylogenies with a Relaxed Clock Syst Biol, December 10, 2009; (2009) syp082v1. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 L. Himmelmann and D. Metzler TreeTime: an extensible C++ software package for Bayesian phylogeny reconstruction with time-calibration Bioinformatics, September 15, 2009; 25(18): 2440 - 2441. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
N. Lartillot, T. Lepage, and S. Blanquart PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating Bioinformatics, September 1, 2009; 25(17): 2286 - 2288. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. Y. W. Ho and M. J. Phillips Accounting for Calibration Uncertainty in Phylogenetic Estimation of Evolutionary Divergence Times Syst Biol, July 3, 2009; (2009) syp035v1. [Full Text] [PDF]
|
|
|
|
 |
|
 S. Y.W. Ho An examination of phylogenetic models of substitution rate variation among lineages Biol Lett, June 23, 2009; 5(3): 421 - 424. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 O. Akerborg, B. Sennblad, L. Arvestad, and J. Lagergren Simultaneous Bayesian gene tree reconstruction and reconciliation analysis PNAS, April 7, 2009; 106(14): 5714 - 5719. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 N. Rodrigue, N. Lartillot, and H. Philippe Bayesian Comparisons of Codon Substitution Models Genetics, November 1, 2008; 180(3): 1579 - 1591. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. S. Renner, G. W. Grimm, G. M. Schneeweiss, T. F. Stuessy, and R. E. Ricklefs Rooting and Dating Maples (Acer) with an Uncorrelated-Rates Molecular Clock: Implications for North American/Asian Disjunctions Syst Biol, October 1, 2008; 57(5): 795 - 808. [Abstract] [Full Text] [PDF]
|
|