A comparison of methods for estimating substitution rates from ancient DNA sequence data.
Abstract: Phylogenetic analysis of DNA from modern and ancient samples allows the reconstruction of important demographic and evolutionary processes. A critical component of these analyses is the estimation of evolutionary rates, which can be calibrated using information about the ages of the samples. However, the reliability of these rate estimates can be negatively affected by among-lineage rate variation and non-random sampling. Using a simulation study, we compared the performance of three phylogenetic methods for inferring evolutionary rates from time-structured data sets: regression of root-to-tip distances, least-squares dating, and Bayesian inference. We also applied these three methods to time-structured mitogenomic data sets from six vertebrate species. Our results from 12 simulation scenarios show that the three methods produce reliable estimates when the substitution rate is high, rate variation is low, and samples of similar ages are not all grouped together in the tree (i.e., low phylo-temporal clustering). The interaction of these factors is particularly important for least-squares dating and Bayesian estimation of evolutionary rates. The three estimation methods produced consistent estimates of rates across most of the six mitogenomic data sets, with sequence data from horses being an exception. We recommend that phylogenetic studies of ancient DNA sequences should use multiple methods of inference and test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data.
Publication Date: 2018-05-16 PubMed ID: 29769015PubMed Central: PMC5956955DOI: 10.1186/s12862-018-1192-3Google Scholar: Lookup
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
- Comparative Study
- Journal Article
- Research Support
- Non-U.S. Gov't
Summary
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
The research in this article focuses on comparing three phylogenetic methods for calculating evolutionary rates using DNA from modern and ancient samples. The aim was to establish which method provided the most reliable results in different circumstances. The authors suggested muti-method inference and testing for certain data characteristics for optimal evaluation of ancient DNA sequences.
Phylogenetic analysis of DNA samples
- The researchers’ main task was to analyze DNA from modern and ancient samples to reconstruct demographic and evolutionary processes. They were particularly interested in estimating evolutionary rates, which can provide insights into the genetic changes in species over time.
- These rates can be calibrated using information about the ages of the samples. However, two issues can negatively affect the rate estimates: variation in rates among different lineages and non-random sampling.
Comparison of different methods
- The study compared three methods for determining evolutionary rates from time-structured data sets. These methods include: regression of root-to-tip distances, least-squares dating, and Bayesian inference.
- These were applied to time-structured mitogenomic data sets from six different vertebrate species to evaluate their effectiveness for different types and amounts of data.
Findings from the simulation scenarios
- The results from 12 simulation scenarios demonstrated that all three methods tend to make reliable estimates when the following conditions are met: the substitution rate is high, rate variation is low, and samples of similar ages are scattered throughout the tree, rather than being grouped together (also known as low phylo-temporal clustering).
- The interaction between these factors was noted to be particularly significant for the least-squares dating and Bayesian estimation methods. This means that if these conditions are not met, these latter two methods may produce less accurate rate estimates.
- The three methods provided largely consistent estimates of rates across most of the six mitogenomic data sets. However, the sequence data from horses was a notable exception, indicating potential limitations of these methods for certain types of data.
Recommendations for future research
- The authors recommend using a mixed-methods approach in future phylogenetic studies of ancient DNA sequence data to provide a more robust analysis.
- Additionally, researchers should also test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data as these factors can significantly affect the accuracy of the evolutionary rate estimates.
Cite This Article
APA
Tong KJ, Duchêne DA, Duchêne S, Geoghegan JL, Ho SYW.
(2018).
A comparison of methods for estimating substitution rates from ancient DNA sequence data.
BMC Evol Biol, 18(1), 70.
https://doi.org/10.1186/s12862-018-1192-3 Publication
Researcher Affiliations
- School of Life and Environmental Sciences, University of Sydney, Sydney, Australia.
- School of Life and Environmental Sciences, University of Sydney, Sydney, Australia.
- Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Melbourne, Australia.
- Department of Biological Sciences, Macquarie University, Sydney, Australia.
- School of Life and Environmental Sciences, University of Sydney, Sydney, Australia. simon.ho@sydney.edu.au.
MeSH Terms
- Animals
- Base Sequence
- Bayes Theorem
- Computer Simulation
- DNA, Ancient
- Evolution, Molecular
- Genome, Mitochondrial
- Genomics / methods
- Horses
- Mutation / genetics
- Phylogeny
- Time Factors
- Uncertainty
- Vertebrates / genetics
Grant Funding
- FT160100167 / Australian Research Council
Conflict of Interest Statement
ETHICS APPROVAL AND CONSENT TO PARTICIPATE: Not applicable. COMPETING INTERESTS: The authors declare that they have no competing interests. PUBLISHER’S NOTE: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
This article includes 59 references
- de Bruyn M, Hoelzel AR, Carvalho GR, Hofreiter M. Faunal histories from Holocene ancient DNA.. Trends Ecol Evol 2011 Aug;26(8):405-13.
- Rambaut A. Estimating the rate of molecular evolution: incorporating non-contemporaneous sequences into maximum likelihood phylogenies.. Bioinformatics 2000 Apr;16(4):395-9.
- Li WH, Tanimura M, Sharp PM. Rates and dates of divergence between AIDS virus nucleotide sequences.. Mol Biol Evol 1988 Jul;5(4):313-30.
- Biek R, Pybus OG, Lloyd-Smith JO, Didelot X. Measurably evolving pathogens in the genomic era.. Trends Ecol Evol 2015 Jun;30(6):306-13.
- Hagelberg E, Hofreiter M, Keyser C. Introduction. Ancient DNA: the first three decades.. Philos Trans R Soc Lond B Biol Sci 2015 Jan 19;370(1660):20130371.
- Drummond AJ, Pybus OG, Rambaut A, Forsberg R, Rodrigo AG. Measurably evolving populations.. Trends Ecol Evol 2003;18(9):481–488.
- Bromham L. Why do species vary in their rate of molecular evolution?. Biol Lett 2009 Jun 23;5(3):401-4.
- Baer CF, Miyamoto MM, Denver DR. Mutation rate variation in multicellular eukaryotes: causes and consequences.. Nat Rev Genet 2007 Aug;8(8):619-31.
- Ho SY, Lanfear R, Phillips MJ, Barnes I, Thomas JA, Kolokotronis SO, Shapiro B. Bayesian estimation of substitution rates from ancient DNA sequences with low information content.. Syst Biol 2011 May;60(3):366-75.
- Rieux A, Balloux F. Inferences from tip-calibrated phylogenies: a review and a practical guide.. Mol Ecol 2016 May;25(9):1911-24.
- Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen).. Virus Evol 2016 Jan;2(1):vew007.
- Buonagurio DA, Nakada S, Parvin JD, Krystal M, Palese P, Fitch WM. Evolution of human influenza A viruses over 50 years: rapid, uniform rate of change in NS gene.. Science 1986 May 23;232(4753):980-2.
- To TH, Jung M, Lycett S, Gascuel O. Fast Dating Using Least-Squares Criteria and Algorithms.. Syst Biol 2016 Jan;65(1):82-97.
- Langley CH, Fitch WM. An examination of the constancy of the rate of molecular evolution.. J Mol Evol 1974;3(3):161-77.
- Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W. Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data.. Genetics 2002 Jul;161(3):1307-20.
- Shapiro B, Ho SY, Drummond AJ, Suchard MA, Pybus OG, Rambaut A. A Bayesian phylogenetic method to estimate unknown sequence ages.. Mol Biol Evol 2011 Feb;28(2):879-87.
- Molak M, Suchard MA, Ho SY, Beilman DW, Shapiro B. Empirical calibrated radiocarbon sampler: a tool for incorporating radiocarbon-date and calibration error into Bayesian phylogenetic analyses of ancient DNA.. Mol Ecol Resour 2015 Jan;15(1):81-6.
- Bromham L, Duchêne S, Hua X, Ritchie AM, Duchêne DA, Ho SYW. Bayesian molecular dating: opening up the black box.. Biol Rev Camb Philos Soc 2018 May;93(2):1165-1191.
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence.. PLoS Biol 2006 May;4(5):e88.
- Rambaut A, Ho SY, Drummond AJ, Shapiro B. Accommodating the effect of ancient DNA damage on inferences of demographic histories.. Mol Biol Evol 2009 Feb;26(2):245-8.
- Ho SY, Heupink TH, Rambaut A, Shapiro B. Bayesian estimation of sequence damage in ancient DNA.. Mol Biol Evol 2007 Jun;24(6):1416-22.
- Ho SY, Kolokotronis SO, Allaby RG. Elevated substitution rates estimated from ancient DNA sequences.. Biol Lett 2007 Dec 22;3(6):702-5.
- Subramanian S, Denver DR, Millar CD, Heupink T, Aschrafi A, Emslie SD, Baroni C, Lambert DM. High mitogenomic evolutionary rates and time dependency.. Trends Genet 2009 Nov;25(11):482-6.
- Molak M, Ho SY. Prolonged decay of molecular rate estimates for metazoan mitochondrial DNA.. PeerJ 2015;3:e821.
- Ho SY, Lanfear R, Bromham L, Phillips MJ, Soubrier J, Rodrigo AG, Cooper A. Time-dependent rates of molecular evolution.. Mol Ecol 2011 Aug;20(15):3087-101.
- Duchêne D, Duchêne S, Ho SY. Tree imbalance causes a bias in phylogenetic estimation of evolutionary timescales using heterochronous sequences.. Mol Ecol Resour 2015 Jul;15(4):785-94.
- Duchêne S, Duchêne D, Holmes EC, Ho SY. The Performance of the Date-Randomization Test in Phylogenetic Analyses of Time-Structured Virus Data.. Mol Biol Evol 2015 Jul;32(7):1895-906.
- Murray GG, Wang F, Harrison EM, Paterson GK, Mather AE, Harris SR, Holmes MA, Rambaut A, Welch JJ. The effect of genetic structure on molecular dating and tests for temporal signal.. Methods Ecol Evol 2016 Jan;7(1):80-89.
- Navascués M, Emerson BC. Elevated substitution rate estimates from ancient DNA: model violation and bias of Bayesian methods.. Mol Ecol 2009 Nov;18(21):4390-7.
- Wertheim JO, Fourment M, Kosakovsky Pond SL. Inconsistencies in estimating the age of HIV-1 subtypes due to heterotachy.. Mol Biol Evol 2012 Feb;29(2):451-6.
- Duchêne S, Geoghegan JL, Holmes EC, Ho SY. Estimating evolutionary rates using time-structured data: a general comparison of phylogenetic methods.. Bioinformatics 2016 Nov 15;32(22):3375-3379.
- Molak M, Lorenzen ED, Shapiro B, Ho SY. Phylogenetic estimation of timescales using ancient DNA: the effects of temporal sampling scheme and uncertainty in sample ages.. Mol Biol Evol 2013 Feb;30(2):253-62.
- Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ. BEAST 2: a software platform for Bayesian evolutionary analysis.. PLoS Comput Biol 2014 Apr;10(4):e1003537.
- Paijmans JL, Gilbert MT, Hofreiter M. Mitogenomic analyses from ancient DNA.. Mol Phylogenet Evol 2013 Nov;69(2):404-16.
- Ho SY, Gilbert MT. Ancient mitogenomics.. Mitochondrion 2010 Jan;10(1):1-11.
- Ho SY, Duchêne S, Duchêne D. Simulating and detecting autocorrelation of molecular evolutionary rates among lineages.. Mol Ecol Resour 2015 Jul;15(4):688-96.
- Rieux A, Eriksson A, Li M, Sobkowiak B, Weinert LA, Warmuth V, Ruiz-Linares A, Manica A, Balloux F. Improved calibration of the human mitochondrial clock using ancient genomes.. Mol Biol Evol 2014 Oct;31(10):2780-92.
- Lepage T, Bryant D, Philippe H, Lartillot N. A general comparison of relaxed molecular clock models.. Mol Biol Evol 2007 Dec;24(12):2669-80.
- Schliep KP. phangorn: phylogenetic analysis in R.. Bioinformatics 2011 Feb 15;27(4):592-3.
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.. Bioinformatics 2014 May 1;30(9):1312-3.
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7.. Mol Biol Evol 2012 Aug;29(8):1969-73.
- Ferreira MAR, Suchard MA. Bayesian analysis of elapsed times in continuous-time Markov chains.. Can J Stat 2008;36:355–368.
- Miller W, Schuster SC, Welch AJ, Ratan A, Bedoya-Reina OC, Zhao F, Kim HL, Burhans RC, Drautz DI, Wittekindt NE, Tomsho LP, Ibarra-Laclette E, Herrera-Estrella L, Peacock E, Farley S, Sage GK, Rode K, Obbard M, Montiel R, Bachmann L, Ingólfsson O, Aars J, Mailund T, Wiig O, Talbot SL, Lindqvist C. Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change.. Proc Natl Acad Sci U S A 2012 Sep 4;109(36):E2382-90.
- Thalmann O, Shapiro B, Cui P, Schuenemann VJ, Sawyer SK, Greenfield DL, Germonpré MB, Sablin MV, López-Giráldez F, Domingo-Roura X, Napierala H, Uerpmann HP, Loponte DM, Acosta AA, Giemsch L, Schmitz RW, Worthington B, Buikstra JE, Druzhkova A, Graphodatsky AS, Ovodov ND, Wahlberg N, Freedman AH, Schweizer RM, Koepfli KP, Leonard JA, Meyer M, Krause J, Pääbo S, Green RE, Wayne RK. Complete mitochondrial genomes of ancient canids suggest a European origin of domestic dogs.. Science 2013 Nov 15;342(6160):871-4.
- Achilli A, Olivieri A, Soares P, Lancioni H, Hooshiar Kashani B, Perego UA, Nergadze SG, Carossa V, Santagostino M, Capomaccio S, Felicetti M, Al-Achkar W, Penedo MC, Verini-Supplizi A, Houshmand M, Woodward SR, Semino O, Silvestrelli M, Giulotto E, Pereira L, Bandelt HJ, Torroni A. Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication.. Proc Natl Acad Sci U S A 2012 Feb 14;109(7):2449-54.
- Lippold S, Matzke NJ, Reissmann M, Hofreiter M. Whole mitochondrial genome sequencing of domestic horses reveals incorporation of extensive wild horse diversity during domestication.. BMC Evol Biol 2011 Nov 14;11:328.
- Orlando L, Ginolhac A, Zhang G, Froese D, Albrechtsen A, Stiller M, Schubert M, Cappellini E, Petersen B, Moltke I, Johnson PL, Fumagalli M, Vilstrup JT, Raghavan M, Korneliussen T, Malaspinas AS, Vogt J, Szklarczyk D, Kelstrup CD, Vinther J, Dolocan A, Stenderup J, Velazquez AM, Cahill J, Rasmussen M, Wang X, Min J, Zazula GD, Seguin-Orlando A, Mortensen C, Magnussen K, Thompson JF, Weinstock J, Gregersen K, Røed KH, Eisenmann V, Rubin CJ, Miller DC, Antczak DF, Bertelsen MF, Brunak S, Al-Rasheid KA, Ryder O, Andersson L, Mundy J, Krogh A, Gilbert MT, Kjær K, Sicheritz-Ponten T, Jensen LJ, Olsen JV, Hofreiter M, Nielsen R, Shapiro B, Wang J, Willerslev E. Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse.. Nature 2013 Jul 4;499(7456):74-8.
- Brotherton P, Haak W, Templeton J, Brandt G, Soubrier J, Jane Adler C, Richards SM, Der Sarkissian C, Ganslmeier R, Friederich S, Dresely V, van Oven M, Kenyon R, Van der Hoek MB, Korlach J, Luong K, Ho SYW, Quintana-Murci L, Behar DM, Meller H, Alt KW, Cooper A. Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans.. Nat Commun 2013;4:1764.
- Gilbert MT, Drautz DI, Lesk AM, Ho SY, Qi J, Ratan A, Hsu CH, Sher A, Dalén L, Götherström A, Tomsho LP, Rendulic S, Packard M, Campos PF, Kuznetsova TV, Shidlovskiy F, Tikhonov A, Willerslev E, Iacumin P, Buigues B, Ericson PG, Germonpré M, Kosintsev P, Nikolaev V, Nowak-Kemp M, Knight JR, Irzyk GP, Perbost CS, Fredrikson KM, Harkins TT, Sheridan S, Miller W, Schuster SC. Intraspecific phylogenetic analysis of Siberian woolly mammoths using complete mitochondrial genomes.. Proc Natl Acad Sci U S A 2008 Jun 17;105(24):8327-32.
- Minin VN, Bloomquist EW, Suchard MA. Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics.. Mol Biol Evol 2008 Jul;25(7):1459-71.
- Xie W, Lewis PO, Fan Y, Kuo L, Chen MH. Improving marginal likelihood estimation for Bayesian phylogenetic model selection.. Syst Biol 2011 Mar;60(2):150-60.
- Kass RE, Raftery AE. Bayes factors.. J Am Stat Assoc 1995;90:773–795.
- Ramsden C, Holmes EC, Charleston MA. Hantavirus evolution in relation to its rodent and insectivore hosts: no evidence for codivergence.. Mol Biol Evol 2009 Jan;26(1):143-53.
- Fiala KL, Sokal RR. FACTORS DETERMINING THE ACCURACY OF CLADOGRAM ESTIMATION: EVALUATION USING COMPUTER SIMULATION.. Evolution 1985 May;39(3):609-622.
- Keis M, Remm J, Ho SYW, Davison J, Tammeleht E, Tumanov IL, Saveljev AP, Männil P, Kojola I, Abramov AV. Complete mitochondrial genomes and a novel spatial genetic method reveal cryptic phylogeographical structure and migration patterns among brown bears in north-western Eurasia.. J Biogeogr 2013;40:915–927.
- Millar CD, Dodd A, Anderson J, Gibb GC, Ritchie PA, Baroni C, Woodhams MD, Hendy MD, Lambert DM. Mutation and evolutionary rates in adélie penguins from the antarctic.. PLoS Genet 2008 Oct 3;4(10):e1000209.
- Williamson S, Orive ME. The genealogy of a sequence subject to purifying selection at multiple sites.. Mol Biol Evol 2002 Aug;19(8):1376-84.
- Mourad R, Chevennet F, Dunn DT, Fearnhill E, Delpech V, Asboe D, Gascuel O, Hue S. A phylotype-based analysis highlights the role of drug-naive HIV-positive individuals in the transmission of antiretroviral resistance in the UK.. AIDS 2015 Sep 24;29(15):1917-25.
- Volz EM, Frost SDW. Scalable relaxed clock phylogenetic dating.. Virus Evol 2017;3(2):vex025.
Citations
This article has been cited 9 times.- Nishimura L, Fujito N, Sugimoto R, Inoue I. Detection of Ancient Viruses and Long-Term Viral Evolution.. Viruses 2022 Jun 18;14(6).
- Översti S, Palo JU. Variation in the Substitution Rates among the Human Mitochondrial Haplogroup U Sublineages.. Genome Biol Evol 2022 Jul 2;14(7).
- Connell JR, Benton MC, Lea RA, Sutherland HG, Chaseling J, Haupt LM, Wright KM, Griffiths LR. Pedigree derived mutation rate across the entire mitochondrial genome of the Norfolk Island population.. Sci Rep 2022 Apr 26;12(1):6827.
- Franceschi VB, Ferrareze PAG, Zimerman RA, Cybis GB, Thompson CE. Mutation hotspots and spatiotemporal distribution of SARS-CoV-2 lineages in Brazil, February 2020-2021.. Virus Res 2021 Oct 15;304:198532.
- Wilson PJ, Rutledge LY. Considering Pleistocene North American wolves and coyotes in the eastern Canis origin story.. Ecol Evol 2021 Jul;11(13):9137-9147.
- Franceschi VB, Caldana GD, de Menezes Mayer A, Cybis GB, Neves CAM, Ferrareze PAG, Demoliner M, de Almeida PR, Gularte JS, Hansen AW, Weber MN, Fleck JD, Zimerman RA, Kmetzsch L, Spilki FR, Thompson CE. Genomic epidemiology of SARS-CoV-2 in Esteio, Rio Grande do Sul, Brazil.. BMC Genomics 2021 May 20;22(1):371.
- Rocha-Méndez A, Sánchez-González LA, González C, Navarro-Sigüenza AG. The geography of evolutionary divergence in the highly endemic avifauna from the Sierra Madre del Sur, Mexico.. BMC Evol Biol 2019 Dec 30;19(1):237.
- Menardo F, Duchêne S, Brites D, Gagneux S. The molecular clock of Mycobacterium tuberculosis.. PLoS Pathog 2019 Sep;15(9):e1008067.
- Garcia-Porta J, Irisarri I, Kirchner M, Rodríguez A, Kirchhof S, Brown JL, MacLeod A, Turner AP, Ahmadzadeh F, Albaladejo G, Crnobrnja-Isailovic J, De la Riva I, Fawzi A, Galán P, Göçmen B, Harris DJ, Jiménez-Robles O, Joger U, Jovanović Glavaš O, Karış M, Koziel G, Künzel S, Lyra M, Miles D, Nogales M, Oğuz MA, Pafilis P, Rancilhac L, Rodríguez N, Rodríguez Concepción B, Sanchez E, Salvi D, Slimani T, S'khifa A, Qashqaei AT, Žagar A, Lemmon A, Moriarty Lemmon E, Carretero MA, Carranza S, Philippe H, Sinervo B, Müller J, Vences M, Wollenberg Valero KC. Environmental temperatures shape thermal physiology as well as diversification and genome-wide substitution rates in lizards.. Nat Commun 2019 Sep 9;10(1):4077.
Use Nutrition Calculator
Check if your horse's diet meets their nutrition requirements with our easy-to-use tool Check your horse's diet with our easy-to-use tool
Talk to a Nutritionist
Discuss your horse's feeding plan with our experts over a free phone consultation Discuss your horse's diet over a phone consultation
Submit Diet Evaluation
Get a customized feeding plan for your horse formulated by our equine nutritionists Get a custom feeding plan formulated by our nutritionists