Bayesian Inference of Sampled Ancestor Trees for Epidemiology and Fossil Calibration

Abstract

Phylogenetic analyses which include fossils or molecular sequences that are sampled through time require models that allow one sample to be a direct ancestor of another sample. As previously available phylogenetic inference tools assume that all samples are tips, they do not allow for this possibility. We have developed and implemented a Bayesian Markov Chain Monte Carlo (MCMC) algorithm to infer what we call sampled ancestor trees, that is, trees in which sampled individuals can be direct ancestors of other sampled individuals. We use a family of birth-death models where individuals may remain in the tree process after sampling, in particular we extend the birth-death skyline model [Stadler et al., 2013] to sampled ancestor trees. This method allows the detection of sampled ancestors as well as estimation of the probability that an individual will be removed from the process when it is sampled. We show that even if sampled ancestors are not of specific interest in an analysis, failing to account for them leads to significant bias in parameter estimates. We also show that sampled ancestor birth-death models where every sample comes from a different time point are non-identifiable and thus require one parameter to be known in order to infer other parameters. We apply our phylogenetic inference accounting for sampled ancestors to epidemiological data, where the possibility of sampled ancestors enables us to identify individuals that infected other individuals after being sampled and to infer fundamental epidemiological parameters. We also apply the method to infer divergence times and diversification rates when fossils are included along with extant species samples, so that fossilisation events are modelled as a part of the tree branching process. Such modelling has many advantages as argued in the literature. The sampler is available as an open-source BEAST2 package (https://github.com/CompEvol/sampled-ancestors).

Full text

Its licence does not allow us to show the full text. Read it at PubMed Central:

References 48

1. 1.

   Yang Z, Rannala B (1997 ) Bayesian phylogenetic inference using DNA sequences: a Markov chain Monte Carlo method. Mol Biol Evol 14: 717–24.

   Not in our library10.1093/oxfordjournals.molbev.a025811

2. 2.

   Mau B, Newton MA, Larget B (1999 ) Bayesian phylogenetic inference via Markov chain Monte Carlo methods. Biometrics 55: 1–12.

   Not in our library10.1111/j.0006-341x.1999.00001.x

3. 3.

   Huelsenbeck JP, Ronquist F (2001 ) MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754–755.

   Not in our library10.1093/bioinformatics/17.8.754

4. 4.

   Drummond AJ, Suchard MA, Xie D, Rambaut A (2012 ) Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 29: 1969–73.

   Not in our library10.1093/molbev/mss075

5. 5.

   Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, et al (2012 ) MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61: 539–42.

   Not in our library10.1093/sysbio/sys029

6. 6.
   BEAST 2: A Software Platform for Bayesian Evolutionary Analysis · 2014

   No notice
7. 7.

   Lartillot N, Lepage T, Blanquart S (2009 ) PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating. Bioinformatics 25: 2286–8.

   Not in our library10.1093/bioinformatics/btp368

8. 8.

   Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, Teller E (1953 ) Equations of state calculations by fast computing machines. The Journal of Chemical Physics 21: 1087–1092.

   Not in our library

9. 9.

   Hastings W (1970 ) Monte Carlo sampling methods using Markov chains and their applications. Biometrika 57: 97–109.

   Not in our library

10. 10.

    Lewis PO, Holder MT, Holsinger KE (2005 ) Polytomies and Bayesian phylogenetic inference. Syst Biol 54: 241–253.

    Not in our library10.1080/10635150590924208

11. 11.

    Gavryushkina A, Welch D, Drummond AJ (2013 ) Recursive algorithms for phylogenetic tree counting. Algorithms for Molecular Biology 8: 26.

    Not in our library10.1186/1748-7188-8-26

12. 12.

    Drummond AJ, Pybus OG, Rambaut A, Forsberg R, Rodrigo AG (2003 ) Measurably evolving populations. Trends in Ecology & Evolution 18: 481–488.

    Not in our library

13. 13.

    Pybus OG, Charleston MA, Gupta S, Rambaut A, Holmes EC, et al (2001 ) The epidemic behavior of the hepatitis C virus. Science 292: 2323–2325.

    Not in our library10.1126/science.1058321

14. 14.

    Grenfell BT, Pybus OG, Gog JR, Wood JL, Daly JM, et al (2004 ) Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303: 327–32.

    Not in our library10.1126/science.1090727

15. 15.

    Stadler T, Kouyos R, von Wyl V, Yerly S, Böni J, et al (2011 ) Estimating the basic reproductive number from viral sequence data. Mol Biol Evol 29: 347–357.

    Not in our library10.1093/molbev/msr217

16. 16.

    Ypma RJ, van Ballegooijen WM, Wallinga J (2013 ) Relating phylogenetic trees to transmission trees of infectious disease outbreaks. Genetics 195: 1055–1062.

    Not in our library10.1534/genetics.113.154856

17. 17.

    Kühnert D, Stadler T, Vaughan TG, Drummond AJ (2014 ) Simultaneous reconstruction of evolutionary history and epidemiological dynamics from viral sequences with the birth-death SIR model. J R Soc Interface 11: 20131106.

    Not in our library10.1098/rsif.2013.1106

18. 18.

    Volz EM, Frost SD (2013 ) Inferring the source of transmission with phylogenetic data. PLoS Comput Biol 9: e1003397.

    Not in our library10.1371/journal.pcbi.1003397

19. 19.

    Teunis P, Heijne JCM, Sukhrie F, van Eijkeren J, Koopmans M, et al (2013 ) Infectious disease transmission as a forensic problem: who infected whom? J R Soc Interface 10: 20120955.

    Not in our library10.1098/rsif.2012.0955

20. 20.

    Vrancken B, Rambaut A, Suchard MA, Drummond A, Baele G, et al (2014 ) The genealogical population dynamics of HIV-1 in a large transmission chain: Bridging within and among host evolutionary rates. PLoS Comput Biol 10: e1003505.

    Not in our library10.1371/journal.pcbi.1003505

21. 21.

    Stadler T (2010 ) Sampling-through-time in birth-death trees. Journal of Theoretical Biology 267: 396–404.

    Not in our library10.1016/j.jtbi.2010.09.010

22. 22.

    Didier G, Royer-Carenzi M, Laurin M (2012 ) The reconstructed evolutionary process with the fossil record. J Theor Biol 315: 26–37.

    Not in our library10.1016/j.jtbi.2012.08.046

23. 23.

    Stadler T, Kühnert D, Bonhoeffer S, Drummond AJ (2013 ) Birth-death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis C virus (HCV). Proc Natl Acad Sci USA 110: 228–33.

    Not in our library10.1073/pnas.1207965110

24. 24.

    Sanderson MJ (1997 ) A nonparametric approach to estimating divergence times in the absence of rate consistency. Mol Biol Evol 14: 1218–1231.

    Not in our library

25. 25.

    Thorne JL, Kishino H, Painter IS (1998 ) Estimating the rate of evolution of the rate of molecular evolution. Mol Biol Evol 15: 1647–1657.

    Not in our library10.1093/oxfordjournals.molbev.a025892

26. 26.
    Relaxed Phylogenetics and Dating with Confidence · 2006

    No notice
27. 27.

    Rannala B, Yang Z (2007 ) Inferring speciation times under an episodic molecular clock. Syst Biol 56: 453–466.

    Not in our library10.1080/10635150701420643

28. 28.

    Ho SY, Phillips MJ (2009 ) Accounting for calibration uncertainty in phylogenetic estimation of evolutionary divergence times. Syst Biol 58: 367–380.

    Not in our library10.1093/sysbio/syp035

29. 29.

    Ronquist F, Klopfstein S, Vilhelmsen L, Schulmeister S, Murray DL, et al (2012 ) A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera. Syst Biol 61: 973–999.

    Not in our library10.1093/sysbio/sys058

30. 30.

    Heath TA, Huelsenbeck JP, Stadler T (2014 ) The fossilized birth–death process for coherent calibration of divergence-time estimates. Proc Natl Acad Sci U S A 111: E2957–E2966.

    Not in our library10.1073/pnas.1319091111

31. 31.

    Heled J, Drummond AJ (2012 ) Calibrated tree priors for relaxed phylogenetics and divergence time estimation. Syst Biol 61: 138–49.

    Not in our library10.1093/sysbio/syr087

32. 32.
    Recent progress in paleontological methods for dating the Tree of Life · 2012

    Not checked
33. 33.

    Pyron RA (2011 ) Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia. Syst Biol 60: 466–81.

    Not in our library10.1093/sysbio/syr047

34. 34.

    Wood HM, Matzke NJ, Gillespie RG, Griswold CE (2013 ) Treating fossils as terminal taxa in divergence time estimation reveals ancient vicariance patterns in the palpimanoid spiders. Syst Biol 62: 264–284.

    Not in our library10.1093/sysbio/sys092

35. 35.

    Schrago CG, Mello B, Soares AER (2013 ) Combining fossil and molecular data to date the diversification of New World Primates. J Evol Biol 26: 2438–2446.

    Not in our library10.1111/jeb.12237

36. 36.

    Silvestro D, Schnitzler J, Liow LH, Antonelli A, Salamin N (2014 ) Bayesian estimation of speciation and extinction from incomplete fossil occurrence data. Syst Biol 63: 349–367.

    Not in our library10.1093/sysbio/syu006

37. 37.

    Foote M (1996 ) On the probability of ancestors in the fossil record. Paleobiology 22: 141–151.

    Not in our library

38. 38.

    Wilkinson RD, Tavaré S (2009 ) Estimating primate divergence times by using conditioned birth-and-death processes. Theor Popul Biol 75: 278–285.

    Not in our library10.1016/j.tpb.2009.02.003

39. 39.

    Green PJ (1995 ) Reversible jump Markov chain Monte Carlo computation and Bayesian model determination. Biometrica 82: 771–732.

    Not in our library

40. 40.

    Wilson IJ, Balding DJ (1998 ) Genealogical inference from microsatellite data. Genetics 150: 499–510.

    Not in our library10.1093/genetics/150.1.499

41. 41.

    Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W (2002 ) Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics 161: 1307–20.

    Not in our library10.1093/genetics/161.3.1307

42. 42.

    Hué S, Pillay D, Clewley JP, Pybus OG (2005 ) Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups. Proc Natl Acad Sci USA 102: 4425–4429.

    Not in our library10.1073/pnas.0407534102

43. 43.

    Dawid AP (1982 ) The well-calibrated Bayesian. J Am Stat Assoc 77: 605–610.

    Not in our library

44. 44.

    Swets JA (1996) Signal detection theory and ROC analysis in psychology and diagnostics: collected papers. Lawrence Erlbaum Associates, Mahwah, NJ.

    Not in our library

45. 45.

    Lewis PO (2001 ) A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol 50: 913–25.

    Not in our library10.1080/106351501753462876

46. 46.

    Tavaré S, Marshall CR, Will O, Soligo C, Martin RD (2002 ) Using the fossil record to estimate the age of the last common ancestor of extant primates. Nature 416: 726–9.

    Not in our library10.1038/416726a

47. 47.

    Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, et al (1999 ) Consistent viral evolutionary changes associated with the disease progression of human immunodeficiency virus type 1 infection. Journal of virology 73: 10489–10502.

    Not in our library10.1128/jvi.73.12.10489-10502.1999

48. 48.

    Wu CH (2014) Bayesian approaches to model uncertainty in phylogenetics. Ph.D. thesis, University of Auckland, Auckland, New Zealand.

    Not in our library

Cited by papers we hold 6

• BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis · 2019
  No notice
• Macroecology and macroevolution of the latitudinal diversity gradient in ants · 2018
  No notice
• Bayesian Divergence-Time Estimation with Genome-Wide Single-Nucleotide Polymorphism Data of Sea Catfishes (Ariidae) Supports Miocene Closure of the Panamanian Isthmus · 2018
  No notice
• bModelTest: Bayesian phylogenetic site model averaging and model comparison · 2017
  No notice
• Inferences from tip‐calibrated phylogenies: a review and a practical guide · 2016
  No notice
• Total-Evidence Dating under the Fossilized Birth–Death Process · 2015
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.