title: Inference of demographic history from genealogical trees using reversible jump Markov chain Monte Carlo.
abstract: BACKGROUND
  Coalescent theory is a general framework to model genetic variation in a population. Specifically, it allows inference about population parameters from sampled DNA sequences. However, most currently employed variants of coalescent theory only consider very simple demographic scenarios of population size changes, such as exponential growth.
  
  
  RESULTS
  Here we develop a coalescent approach that allows Bayesian non-parametric estimation of the demographic history using genealogies reconstructed from sampled DNA sequences. In this framework inference and model selection is done using reversible jump Markov chain Monte Carlo (MCMC). This method is computationally efficient and overcomes the limitations of related non-parametric approaches such as the skyline plot. We validate the approach using simulated data. Subsequently, we reanalyze HIV-1 sequence data from Central Africa and Hepatitis C virus (HCV) data from Egypt.
  
  
  CONCLUSIONS
  The new method provides a Bayesian procedure for non-parametric estimation of the demographic history. By construction it additionally provides confidence limits and may be used jointly with other MCMC-based coalescent approaches.
authors.affiliation: Department of Statistics, University of Munich, Ludwigstr, 33, D-80539 Munich, Germany. opgen@stat.uni-muenchen.de
authors.forename: Rainer
authors.initials: R
authors.lastname: Opgen-Rhein
authors.forename: Ludwig
authors.initials: L
authors.lastname: Fahrmeir
authors.forename: Korbinian
authors.initials: K
authors.lastname: Strimmer
mesh_terms: D000465:Algorithms; D001499:Bayes Theorem; D004247:DNA / Q000235:genetics; D003627:Data Interpretation, Statistical; D019143:Evolution, Molecular; D005814:Genes, Viral; D014644:Genetic Variation; D005828:Genetics, Population; D015497:HIV-1 / Q000235:genetics; D016174:Hepacivirus / Q000235:genetics; D008390:Markov Chains; D008957:Models, Genetic; D015233:Models, Statistical; D009010:Monte Carlo Method; D012984:Software
chemical_list: D004247:DNA