Bayesian estimation of scaled mutation rate under the coalescent: a sequential Monte Carlo approach

BMC Bioinformatics

Research Article

Xiaodong Wang¹  Oyetunji E. Ogundijo¹ 

[1] Department of Electrical Engineering, Columbia University, 10027, New York, USA;

关键词: Coalescent;    Sequential Monte Carlo;    Genealogy;    Bayesian;   

DOI  :  10.1186/s12859-017-1948-6

 received in 2017-07-15, accepted in 2017-11-21,  发布年份 2017

来源: Springer

【 摘 要 】

BackgroundSamples of molecular sequence data of a locus obtained from random individuals in a population are often related by an unknown genealogy. More importantly, population genetics parameters, for instance, the scaled population mutation rate Θ=4Neμ for diploids or Θ=2Neμ for haploids (where Ne is the effective population size and μ is the mutation rate per site per generation), which explains some of the evolutionary history and past qualities of the population that the samples are obtained from, is of significant interest.ResultsIn this paper, we present the evolution of sequence data in a Bayesian framework and the approximation of the posterior distributions of the unknown parameters of the model, which include Θ via the sequential Monte Carlo (SMC) samplers for static models. Specifically, we approximate the posterior distributions of the unknown parameters with a set of weighted samples i.e., the set of highly probable genealogies out of the infinite set of possible genealogies that describe the sampled sequences. The proposed SMC algorithm is evaluated on simulated DNA sequence datasets under different mutational models and real biological sequences. In terms of the accuracy of the estimates, the proposed SMC method shows a comparable and sometimes, better performance than the state-of-the-art MCMC algorithms.ConclusionsWe showed that the SMC algorithm for static model is a promising alternative to the state-of-the-art approach for simulating from the posterior distributions of population genetics parameters.

【 授权许可】

CC BY   
© The Author(s) 2017

【 图 表 】

12864_2015_2297_Article_IEq4.gif

12864_2015_2129_Article_IEq18.gif

12864_2017_3920_Article_IEq2.gif

12864_2017_3733_Article_IEq36.gif

12864_2017_3645_Article_IEq1.gif

12864_2017_3733_Article_IEq35.gif

12864_2015_1944_Article_IEq4.gif

12864_2017_3605_Article_IEq2.gif

12864_2016_2793_Article_IEq36.gif

12864_2016_2793_Article_IEq35.gif

12864_2015_2055_Article_IEq32.gif

12864_2017_4132_Article_IEq28.gif

12864_2015_2296_Article_IEq52.gif

12864_2017_4179_Article_IEq35.gif

12864_2017_3990_Article_IEq8.gif

12864_2017_4179_Article_IEq34.gif

12864_2016_3425_Article_IEq2.gif

12880_2017_231_Article_IEq5.gif

12864_2017_3670_Article_IEq9.gif

12711_2016_271_Article_IEq59.gif

12864_2015_2192_Article_IEq19.gif

12864_2015_2174_Article_IEq1.gif

12864_2015_2273_Article_IEq14.gif

12864_2017_4132_Article_IEq14.gif

12864_2016_3098_Article_IEq70.gif

12864_2017_4030_Article_IEq29.gif

12864_2017_4226_Article_IEq2.gif

12864_2017_4030_Article_IEq28.gif

12864_2017_4004_Article_IEq15.gif

12864_2016_2392_Article_IEq2.gif

12864_2017_4004_Article_IEq14.gif

12864_2016_2392_Article_IEq1.gif

12864_2015_2300_Article_IEq26.gif

12864_2017_4004_Article_IEq12.gif

12864_2017_4116_Article_IEq2.gif

12864_2017_4166_Article_IEq3.gif

12864_2017_3547_Article_IEq1.gif

12864_2016_3105_Article_IEq15.gif

12864_2015_2192_Article_IEq11.gif

12888_2016_877_Article_IEq6.gif

12864_2016_3440_Article_IEq55.gif

12864_2015_2192_Article_IEq10.gif

12888_2016_877_Article_IEq4.gif

12864_2016_2789_Article_IEq16.gif

12888_2016_877_Article_IEq2.gif