Abstract
Scientists are assembling sequence data sets from increasing numbers of species and genes to build comprehensive timetrees. However, data are often unavailable for some species and gene combinations, and the proportion of missing data is often large for data sets containing many genes and species. Surprisingly, there has not been a systematic analysis of the effect of the degree of sparseness of the species-gene matrix on the accuracy of divergence time estimates. Here, we present results from computer simulations and empirical data analyses to quantify the impact of missing gene data on divergence time estimation in large phylogenies. We found that estimates of divergence times were robust even when sequences from a majority of genes for most of the species were absent. From the analysis of such extremely sparse data sets, we found that the most egregious errors occurred for nodes in the tree that had no common genes for any pair of species in the immediate descendant clades of the node in question. These problematic nodes can be easily detected prior to computational analyses based only on the input sequence alignment and the tree topology. We conclude that it is best to use larger alignments, because adding both genes and species to the alignment augments the number of genes available for estimating divergence events deep in the tree and improves their time estimates.
Keywords
divergence time
incomplete data
timetree
MeSH Terms
Computer Simulation
Evolution, Molecular
Genes
Models, Genetic
Phylogeny
Sequence Alignment/methods
Sequence Analysis, DNA
Authors & Affiliations
5 authors, click to expand affiliations / ORCID
Filipski Alan
Center for Evolutionary Medicine and Informatics, Biodesign Institute, Arizona State University.
Murillo Oscar
Center for Evolutionary Medicine and Informatics, Biodesign Institute, Arizona State UniversitySchool of Life Sciences, Arizona State University.
Freydenzon Anna
Center for Evolutionary Medicine and Informatics, Biodesign Institute, Arizona State University.
Tamura Koichiro
Department of Biological Sciences, Tokyo Metropolitan University, Tokyo, JapanResearch Center for Genomics and Bioinformatics, Tokyo Metropolitan University, Tokyo, Japan.
Kumar Sudhir
Center for Evolutionary Medicine and Informatics, Biodesign Institute, Arizona State UniversitySchool of Life Sciences, Arizona State UniversityCenter of Excellence in Genomic Medicine Research, King Abdulaziz University, Jeddah, Saudi ArabiaInstitute for Genomics and Evolutionary Medicine, Temple UniversityDepartment of Biology, Temple University s.kumar@temple.edu.
References (26)
26 references, click to expand
-
[Molecular dating in the genomic era].
Med Sci (Paris). 2006 Apr;22(4):374-80
PMID: 16597406
-
Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.
J Mol Evol. 1985;22(2):160-74
PMID: 3934395
-
Rate variation and estimation of divergence times using strict and relaxed clocks.
BMC Evol Biol. 2011 Sep 26;11:271
PMID: 21943087
-
Heterogeneity of nucleotide frequencies among evolutionary lineages and phylogenetic inference.
Mol Biol Evol. 2003 Apr;20(4):610-21
PMID: 12679548
-
Estimating divergence times in large molecular phylogenies.
Proc Natl Acad Sci U S A. 2012 Nov 20;109(47):19333-8
PMID: 23129628
-
Divergence time and evolutionary rate estimation with multilocus data.
Syst Biol. 2002 Oct;51(5):689-702
PMID: 12396584
-
MEGA6: Molecular Evolutionary Genetics Analysis version 6.0.
Mol Biol Evol. 2013 Dec;30(12):2725-9
PMID: 24132122
-
Divergence dates for Malagasy lemurs estimated from multiple gene loci: geological and evolutionary context.
Mol Ecol. 2004 Apr;13(4):757-73
PMID: 15012754
-
Performance of a divergence time estimation method under a probabilistic model of rate evolution.
Mol Biol Evol. 2001 Mar;18(3):352-61
PMID: 11230536
-
Missing data in phylogenetic analysis: reconciling results from simulations and empirical data.
Syst Biol. 2011 Oct;60(5):719-31
PMID: 21447483
-
Seq-Gen: an application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees.
Comput Appl Biosci. 1997 Jun;13(3):235-8
PMID: 9183526
-
Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach.
Mol Biol Evol. 2002 Jan;19(1):101-9
PMID: 11752195
-
Estimating divergence dates from molecular sequences.
Mol Biol Evol. 1998 Apr;15(4):442-8
PMID: 9549094
-
Performance of relaxed-clock methods in estimating evolutionary divergence times and their credibility intervals.
Mol Biol Evol. 2010 Jun;27(6):1289-300
PMID: 20093431
-
Phylogenomics with incomplete taxon coverage: the limits to inference.
BMC Evol Biol. 2010 May 25;10:155
PMID: 20500873
-
Impact of missing data on phylogenies inferred from empirical phylogenomic data sets.
Mol Biol Evol. 2013 Jan;30(1):197-214
PMID: 22930702
-
Efficiency of the neighbor-joining method in reconstructing deep and shallow evolutionary relationships in large phylogenies.
J Mol Evol. 2000 Dec;51(6):544-53
PMID: 11116328
-
BEAST: Bayesian evolutionary analysis by sampling trees.
BMC Evol Biol. 2007 Nov 08;7:214
PMID: 17996036
-
The timing of eukaryotic evolution: does a relaxed molecular clock reconcile proteins and fossils?
Proc Natl Acad Sci U S A. 2004 Oct 26;101(43):15386-91
PMID: 15494441
-
Phylogenomics of eukaryotes: impact of missing data on large alignments.
Mol Biol Evol. 2004 Sep;21(9):1740-52
PMID: 15175415
-
A stepwise algorithm for finding minimum evolution trees.
Mol Biol Evol. 1996 Apr;13(4):584-93
PMID: 8882501
-
Molecular clocks: four decades of evolution.
Nat Rev Genet. 2005 Aug;6(8):654-62
PMID: 16136655
-
Impacts of the Cretaceous Terrestrial Revolution and KPg extinction on mammal diversification.
Science. 2011 Oct 28;334(6055):521-4
PMID: 21940861
-
Highly incomplete taxa can rescue phylogenetic analyses from the negative impacts of limited taxon sampling.
PLoS One. 2012;7(8):e42925
PMID: 22900065
-
A molecular timescale for vertebrate evolution.
Nature. 1998 Apr 30;392(6679):917-20
PMID: 9582070
-
The effect of ambiguous data on phylogenetic estimates obtained by maximum likelihood and Bayesian inference.
Syst Biol. 2009 Feb;58(1):130-45
PMID: 20525573