Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles

Daniel G. Mulcahy, Brice P. Noonan, Travis Moss, Ted M. Townsend, Tod W. Reeder, Jack W. Sites, John J Wiens

Research output: Contribution to journalArticle

60 Citations (Scopus)

Abstract

Recently, phylogenetics has expanded to routinely include estimation of clade ages in addition to their relationships. Various dating methods have been used, but their relative performance remains understudied. Here, we generate and assemble an extensive phylogenomic data set for squamate reptiles (lizards and snakes) and evaluate two widely used dating methods, penalized likelihood in r8s (r8s-PL) and Bayesian estimation with uncorrelated relaxed rates among lineages (BEAST). We obtained sequence data from 25 nuclear loci (∼500-1000. bp per gene; 19,020. bp total) for 64 squamate species and nine outgroup taxa, estimated the phylogeny, and estimated divergence dates using 14 fossil calibrations. We then evaluated how well each method approximated these dates using random subsets of the nuclear loci (2, 5, 10, 15, and 20; replicated 10 times each), and using ∼1. kb of the mitochondrial ND2 gene. We find that estimates from r8s-PL based on 2, 5, or 10 loci can differ considerably from those based on 25 loci (mean absolute value of differences between 2-locus and 25-locus estimates were 9.0. Myr). Estimates from BEAST are somewhat more consistent given limited sampling of loci (mean absolute value of differences between 2 and 25-locus estimates were 5.0. Myr). Most strikingly, age estimates using r8s-PL for ND2 were ∼68-82. Myr older (mean = 73.1) than those using 25 nuclear loci with r8s-PL. These results show that dates from r8s-PL with a limited number of loci (and especially mitochondrial data) can differ considerably from estimates derived from a large number of nuclear loci, whereas estimates from BEAST derived from fewer nuclear loci or mitochondrial data alone can be surprisingly similar to those from many nuclear loci. However, estimates from BEAST using relatively few loci and mitochondrial data could still show substantial deviations from the full data set (>50. Myr), suggesting the benefits of sampling many nuclear loci. Finally, we found that confidence intervals on ages from BEAST were not significantly different when sampling 2 vs. 25 loci, suggesting that adding loci decreased errors but did not increase confidence in those estimates.

Original languageEnglish (US)
Pages (from-to)974-991
Number of pages18
JournalMolecular Phylogenetics and Evolution
Volume65
Issue number3
DOIs
StatePublished - Dec 2012
Externally publishedYes

Fingerprint

Reptiles
dating method
Squamata
reptile
reptiles
divergence
loci
sampling
Lizards
Mitochondrial Genes
Snakes
Phylogeny
Calibration
gene
snake
methodology
lizard
Confidence Intervals
confidence interval
phylogeny

Keywords

  • BEAST
  • Lizards
  • Penalized likelihood
  • Phylogeny
  • R8s
  • Snakes

ASJC Scopus subject areas

  • Ecology, Evolution, Behavior and Systematics
  • Molecular Biology
  • Genetics

Cite this

Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles. / Mulcahy, Daniel G.; Noonan, Brice P.; Moss, Travis; Townsend, Ted M.; Reeder, Tod W.; Sites, Jack W.; Wiens, John J.

In: Molecular Phylogenetics and Evolution, Vol. 65, No. 3, 12.2012, p. 974-991.

Research output: Contribution to journalArticle

Mulcahy, Daniel G. ; Noonan, Brice P. ; Moss, Travis ; Townsend, Ted M. ; Reeder, Tod W. ; Sites, Jack W. ; Wiens, John J. / Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles. In: Molecular Phylogenetics and Evolution. 2012 ; Vol. 65, No. 3. pp. 974-991.
@article{9a097ec2ca40418f902c61532e8c5309,
title = "Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles",
abstract = "Recently, phylogenetics has expanded to routinely include estimation of clade ages in addition to their relationships. Various dating methods have been used, but their relative performance remains understudied. Here, we generate and assemble an extensive phylogenomic data set for squamate reptiles (lizards and snakes) and evaluate two widely used dating methods, penalized likelihood in r8s (r8s-PL) and Bayesian estimation with uncorrelated relaxed rates among lineages (BEAST). We obtained sequence data from 25 nuclear loci (∼500-1000. bp per gene; 19,020. bp total) for 64 squamate species and nine outgroup taxa, estimated the phylogeny, and estimated divergence dates using 14 fossil calibrations. We then evaluated how well each method approximated these dates using random subsets of the nuclear loci (2, 5, 10, 15, and 20; replicated 10 times each), and using ∼1. kb of the mitochondrial ND2 gene. We find that estimates from r8s-PL based on 2, 5, or 10 loci can differ considerably from those based on 25 loci (mean absolute value of differences between 2-locus and 25-locus estimates were 9.0. Myr). Estimates from BEAST are somewhat more consistent given limited sampling of loci (mean absolute value of differences between 2 and 25-locus estimates were 5.0. Myr). Most strikingly, age estimates using r8s-PL for ND2 were ∼68-82. Myr older (mean = 73.1) than those using 25 nuclear loci with r8s-PL. These results show that dates from r8s-PL with a limited number of loci (and especially mitochondrial data) can differ considerably from estimates derived from a large number of nuclear loci, whereas estimates from BEAST derived from fewer nuclear loci or mitochondrial data alone can be surprisingly similar to those from many nuclear loci. However, estimates from BEAST using relatively few loci and mitochondrial data could still show substantial deviations from the full data set (>50. Myr), suggesting the benefits of sampling many nuclear loci. Finally, we found that confidence intervals on ages from BEAST were not significantly different when sampling 2 vs. 25 loci, suggesting that adding loci decreased errors but did not increase confidence in those estimates.",
keywords = "BEAST, Lizards, Penalized likelihood, Phylogeny, R8s, Snakes",
author = "Mulcahy, {Daniel G.} and Noonan, {Brice P.} and Travis Moss and Townsend, {Ted M.} and Reeder, {Tod W.} and Sites, {Jack W.} and Wiens, {John J}",
year = "2012",
month = "12",
doi = "10.1016/j.ympev.2012.08.018",
language = "English (US)",
volume = "65",
pages = "974--991",
journal = "Molecular Phylogenetics and Evolution",
issn = "1055-7903",
publisher = "Academic Press Inc.",
number = "3",

}

TY - JOUR

T1 - Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles

AU - Mulcahy, Daniel G.

AU - Noonan, Brice P.

AU - Moss, Travis

AU - Townsend, Ted M.

AU - Reeder, Tod W.

AU - Sites, Jack W.

AU - Wiens, John J

PY - 2012/12

Y1 - 2012/12

N2 - Recently, phylogenetics has expanded to routinely include estimation of clade ages in addition to their relationships. Various dating methods have been used, but their relative performance remains understudied. Here, we generate and assemble an extensive phylogenomic data set for squamate reptiles (lizards and snakes) and evaluate two widely used dating methods, penalized likelihood in r8s (r8s-PL) and Bayesian estimation with uncorrelated relaxed rates among lineages (BEAST). We obtained sequence data from 25 nuclear loci (∼500-1000. bp per gene; 19,020. bp total) for 64 squamate species and nine outgroup taxa, estimated the phylogeny, and estimated divergence dates using 14 fossil calibrations. We then evaluated how well each method approximated these dates using random subsets of the nuclear loci (2, 5, 10, 15, and 20; replicated 10 times each), and using ∼1. kb of the mitochondrial ND2 gene. We find that estimates from r8s-PL based on 2, 5, or 10 loci can differ considerably from those based on 25 loci (mean absolute value of differences between 2-locus and 25-locus estimates were 9.0. Myr). Estimates from BEAST are somewhat more consistent given limited sampling of loci (mean absolute value of differences between 2 and 25-locus estimates were 5.0. Myr). Most strikingly, age estimates using r8s-PL for ND2 were ∼68-82. Myr older (mean = 73.1) than those using 25 nuclear loci with r8s-PL. These results show that dates from r8s-PL with a limited number of loci (and especially mitochondrial data) can differ considerably from estimates derived from a large number of nuclear loci, whereas estimates from BEAST derived from fewer nuclear loci or mitochondrial data alone can be surprisingly similar to those from many nuclear loci. However, estimates from BEAST using relatively few loci and mitochondrial data could still show substantial deviations from the full data set (>50. Myr), suggesting the benefits of sampling many nuclear loci. Finally, we found that confidence intervals on ages from BEAST were not significantly different when sampling 2 vs. 25 loci, suggesting that adding loci decreased errors but did not increase confidence in those estimates.

AB - Recently, phylogenetics has expanded to routinely include estimation of clade ages in addition to their relationships. Various dating methods have been used, but their relative performance remains understudied. Here, we generate and assemble an extensive phylogenomic data set for squamate reptiles (lizards and snakes) and evaluate two widely used dating methods, penalized likelihood in r8s (r8s-PL) and Bayesian estimation with uncorrelated relaxed rates among lineages (BEAST). We obtained sequence data from 25 nuclear loci (∼500-1000. bp per gene; 19,020. bp total) for 64 squamate species and nine outgroup taxa, estimated the phylogeny, and estimated divergence dates using 14 fossil calibrations. We then evaluated how well each method approximated these dates using random subsets of the nuclear loci (2, 5, 10, 15, and 20; replicated 10 times each), and using ∼1. kb of the mitochondrial ND2 gene. We find that estimates from r8s-PL based on 2, 5, or 10 loci can differ considerably from those based on 25 loci (mean absolute value of differences between 2-locus and 25-locus estimates were 9.0. Myr). Estimates from BEAST are somewhat more consistent given limited sampling of loci (mean absolute value of differences between 2 and 25-locus estimates were 5.0. Myr). Most strikingly, age estimates using r8s-PL for ND2 were ∼68-82. Myr older (mean = 73.1) than those using 25 nuclear loci with r8s-PL. These results show that dates from r8s-PL with a limited number of loci (and especially mitochondrial data) can differ considerably from estimates derived from a large number of nuclear loci, whereas estimates from BEAST derived from fewer nuclear loci or mitochondrial data alone can be surprisingly similar to those from many nuclear loci. However, estimates from BEAST using relatively few loci and mitochondrial data could still show substantial deviations from the full data set (>50. Myr), suggesting the benefits of sampling many nuclear loci. Finally, we found that confidence intervals on ages from BEAST were not significantly different when sampling 2 vs. 25 loci, suggesting that adding loci decreased errors but did not increase confidence in those estimates.

KW - BEAST

KW - Lizards

KW - Penalized likelihood

KW - Phylogeny

KW - R8s

KW - Snakes

UR - http://www.scopus.com/inward/record.url?scp=84867715558&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84867715558&partnerID=8YFLogxK

U2 - 10.1016/j.ympev.2012.08.018

DO - 10.1016/j.ympev.2012.08.018

M3 - Article

C2 - 22982760

AN - SCOPUS:84867715558

VL - 65

SP - 974

EP - 991

JO - Molecular Phylogenetics and Evolution

JF - Molecular Phylogenetics and Evolution

SN - 1055-7903

IS - 3

ER -