Cited by 546Genomics and Phylogenetic StudiesEvolution and Paleontology StudiesEvolution and Genetic Dynamics

The DOI matches its Crossref record; 46 of the 46 references it cites were found in a registry or an index; none of the 27 sources its claims cite is checked for support yet.

A General Comparison of Relaxed Molecular Clock Models

Open access
Lepage, T., Bryant, D., Philippe, H., Lartillot, N.·2007·Molecular Biology and Evolution

Abstract

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.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper benchmarks relaxed molecular clock models and divergence-time priors in a common Bayesian framework, adding a new autocorrelated CIR process. Across three eukaryotic protein data sets, autocorrelated rate models fit better than uncorrelated ones, while the best divergence-time prior depends on the data set.

Key findings

  • The strict clock was strongly rejected relative to the unconstrained model on all data sets, with log Bayes factors ranging from about −157 to −3,072.
  • Autocorrelated models, LOGN and CIR, had similar fit and clearly outperformed uncorrelated models (UEXP, UGAM, WN) on nearly all data sets; LOGN was slightly favored except on MITO100, where they were equivalent, and MITO150, where CIR was significantly preferred.
  • The Dirichlet prior was generally close to the uniform prior, whereas the birth-death prior was rejected in most cases (ln BF < −6) but strongly favored on MITO100 and MITO150 (ln BF > 45).
  • Dense taxonomic sampling was important for distinguishing autocorrelated from uncorrelated models; with only 10 taxa, the uncorrelated models were comparable to the autocorrelated ones.

Methods

  • The authors reimplemented seven relaxed-clock models—STRICT, UEXP, UGAM, WN, LOGN, CIR, and an unconstrained reference—and combined them with uniform, Dirichlet, and birth-death priors on divergence times.
  • Data sets were MOC (129 genes, 30,399 positions, 36 eukaryotes), MAM (3 nuclear genes, 3,700 positions, 60 mammals), and MITO (11 mitochondrial protein-coding genes, 3,220 positions, 155 osteichthyans plus 5 chondrichthyans, with ingroup subsamples of 10, 50, 100, and 150 taxa).
  • Posterior sampling used MCMC in an extension of PhyloBayes under the normal approximation, and model fit was compared by Bayes factors computed with thermodynamic integration.

Limitations

  • The paper states that no fossil calibrations were used, so only relative divergence times were examined and no conclusions about absolute dates or root age are provided.
  • The analyses fixed the tree topology, used the normal approximation exclusively, and did not account for heterogeneity across sites or data partitions.
  • The authors note that the reasons for the birth-death prior’s data-dependent performance are unclear, that sampling error in Bayes factor estimates was on the order of a few log units, and that substitution-model effects and other relaxed-clock models were not investigated.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1093/molbev/msm193

Cites

46 of 46 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Verification

Claims the paper makes

0 of 27 cited sources checked · 0 supported · 0 flagged

Show all 48 claims

  • It may state two separate claims as one. 28
  • Probably not a substantial claim. 14
38 of 38
  • The authors state this as an assumption.

    The rate process is Markovian.

    “Assuming that the rate process is Markovian, the joint probability density of the values of the process at the nodes is given by:”
  • The paper relies on this without support.

    The strict molecular clock is not biologically realistic and substitution rates change over time, requiring explicit accounting in molecular dating.

    “However, it has become clear that the strict molecular clock hypothesis is not biologically realistic and that the rate at which substitutions accumulate with time is subject to changes that need to be explicitly accounted for in molecular dating techniques.”
  • The paper relies on this without support.

    Few attempts have been made to compare the performance of available clock models.

    “Remarkably, there have been few attempts thus far to compare the performances of available clock models.”
  • The paper relies on this without support.

    The random subsampling of osteichthyans in the MITO data sets produces taxon samplings that differ only in number of taxa, so differences in model performance can be attributed to taxonomic sampling density.

    “considered 4 different taxon samplings for the ingroup, obtained by randomly subsampling 10, 50, 100, and 150 osteichthyans.”
  • The paper relies on this without support.

    The data sets used are representative of real sequence evolution, so conclusions about the strict molecular clock and autocorrelation generalize to real sequences.

    “indicating that the strict molecular clock is not adapted to describing the evolution of real sequences.”
5 of 5

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

  • Autocorrelated models do not outperform UGAM on MITO10

    Minor

    The abstract says both autocorrelated models clearly outperform uncorrelated models on all three data sets. For MITO10, table 1 gives UGAM a higher log Bayes factor (3.9-4.0) than either LOGN (3.2-3.7) or CIR (2.5-3.0).

  • MITO100 birth-death Bayes factor is not above 45

    Minor

    The stated threshold holds for MITO150, but not MITO100. Table 2 reports MITO100 log Bayes factors of 13.3-17.4 for CIR and 13.3-17.1 for LOGN under the birth-death prior.

  • Figure 1b uses a different prior from the cited model

    Minor

    The text points to figure 1b when discussing the CIR plus birth-death model, but the figure caption identifies its prior as uniform.

3 of 3

Citations

Cites

1 retracted

Flagged first, then by citations

Research works printed but not linked: 1

  1. Kimura M. 1983. The neutral allele theory of molecular evolution. Cambridge: Cambridge University Press.

Cited by

Sorted by citations

See all 546

Metadata

Authors
Lepage, T., Bryant, D., Philippe, H., & Lartillot, N.
Affiliations
  • McGill University (CA)
  • University of Auckland (NZ)
  • Université de Montréal (CA)
  • Centre National de la Recherche Scientifique (FR)
Venue
Molecular Biology and Evolution
Published
2007-06-29
Volume
24
Issue
12
Pages
2669-2680
ISSN
  • 0737-4038 (ISSN-L)
  • 1537-1719
OpenAlex
W2116474217
Cited by
546
References
46
Citations by year
2026: 27 · 2025: 28 · 2024: 24 · 2023: 33 · 2022: 39 · 2021: 34 · 2020: 38 · 2019: 29 · 2018: 43 · 2017: 33 · 2016: 32 · 2015: 38 · 2014: 50 · 2013: 18 · 2012: 19
FWCI
4.98
Type
journal article
Language
en
Open access
Locations
  • Molecular Biology and Evolution (publishedVersion)
  • PubMed (publishedVersion)
  • HAL (submittedVersion, cc-by)
Funders
  • Génome Québec
  • Centre National de la Recherche Scientifique
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
  • Physical Sciences › Earth and Planetary Sciences › Paleontology › Evolution and Paleontology Studies
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Evolution and Genetic Dynamics
Keywords
relaxed molecular clock, divergence time estimation, CIR process, Dirichlet prior, model comparison, Bayes factor, thermodynamic integration, protein sequence data
MeSH terms
Bayes Theorem, Models, Molecular, Phylogeny, Time Factors, Genetic Variation

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-05, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-05 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review