A Markov Chain Monte Carlo expectation maximization algorithm for statistical analysis of DNA sequence evolution with neighbor-dependent substitution rates

Journal of Computational and Graphical Statistics

Volumn 17, Issue 1, 2008, Pages 138-162

  Hobolth, Asger ^a  

^a NORTH CAROLINA STATE UNIVERSITY   (United States)

Author keywords

EM algorithm; Gibbs sampling; Likelihood inference; Molecular evolution; Neighbor dependence; Path sampling

Indexed keywords

EID: 37549017437     PISSN: 10618600     EISSN: None     Source Type: Journal    
DOI: 10.1198/106186008X289010     Document Type: Article

References (36)

1. Albert, B., Johnson, A., Lewis, J., Raff, M., Roberts, K. and Walter, P. (2002), Molecular Biology of the Cell, New York: Garland Science.
2. Arndt, P.F., and Hwa, T. (2005), "Identification and Measurement of Neighbour-Dependent Nucleotide Substitution Processes," Bioinformatics, 21, 2322-2328.
3. Bladt, M., and Sørensen, M. (2005), "Statistical Inference for Discretely Observed Markov Jump Processes," Journal of the Royal Statistical Society, Ser. B, 67, 395-410.
4. Blake, R.D., Hess, S.T., and Nicholson-Tuell, J. (1992), "The Influence of Nearest Neighbors on the Rate and Pattern of Spontaneous Point Mutations," Journal of Molecular Evolution, 34, 189-200.
5. Caffo, B.S., Jank, W., and Jones, G.L. (2005), "Ascent-Based Monte Carlo Expectation-Maximization," Journal of the Royal Statistical Society, Ser. B, 67, 235-251.
6. Christensen, O.F., Hobolth, A., and Jensen, J.L. (2005), "Pseudo-Likelihood Analysis of Context-Dependent Codon Substitution Models," Journal of Computational Biology, 12, 1166-1182.
7. Clark, A.G., Glanowski, S., Nielsen, R., Thomas, P.D., Kejariwal, A., Todd, M.A., Tanenbaum, D.M., Civello, D., Lu, F., Murphy, B., Ferriera, S., Wang, G., Zheng, X., White, T.J., Sninsky, J.J., Adams, M.D., Cargill, M. (2003), "Inferring Nonneutral Evolution from Human-Chimp-Mouse Orthologous Gene Trios," Science, 302, 1960-1963.
8. Dempster, A.P., Laird, N.M., and Rubin, D.B. (1977), "Maximum Likelihood from Incomplete Data via the EM Algorithm," Journal of the Royal Statistical Society, Ser. B, 39, 1-22.
9. Diaconis, P., and Rolles, S.W.W. (2006), "Bayesian Analysis for Reversible Markov Chains," The Annals of Statistics, 34, 1270-1292.
10. Drton, M. (2004), "Maximum Likelihood Estimation in Gaussian AMP Chain Graph Models and Gaussian Ancestral Graph Models," unpublished Ph.D. thesis, Department of Statistics, University of Washington.
11. Ewens, W. J., and Grant, G.R. (2001), Statistical Methods in Bioinformatics, New York: Springer.
12. Fort, G., and Moulines, E. (2003), "Convergence of the Monte Carlo Expectation Maximization for Curved Exponential Families," The Annals of Statistics, 31, 1220-1259.
13. Guindon, S., and Gascuel, O. (2003), "A Simple, Fast, and Accurate Algorithm to Estimate Large Phylogenies by Maximum Likelihood," Systematic Biology, 52, 696-704.
14. Hess, S.T., Blake, J.D., and Blake, R.D. (1994), "Wide Variations in Neighbour-Dependent Substitution Rates," Journal of Molecular Biology, 236, 1022-1033.
15. Holmes, I., and Rubin, G.M. (2002), "An Expectation Maximization Algorithm for Training Hidden Substitution Models," Journal of Molecular Biology, 317, 757-768.
16. Hobolth, A., and Jensen, J.L. (2005a), "Statistical Inference in Evolutionary Models of DNA Sequences via the EM Algorithm," Statistical applications in Genetics and Molecular Biology, 4, 18.
17. (2005b), "Applications of Hidden Markov Models for Characterization of Homologous DNA Sequences with a Common Gene," Journal of Computational Biology, 12, 186-203.
18. Huelsenbeck, J.P., Nielsen, R., and Bollback, J.P. (2003), "Stochastic Mapping of Morphological Characters," Systematic Biology, 52, 131-158.
19. Hwang, D.G., and Green, P. (2004), "Bayesian Markov Chain Monte Carlo Sequence Analysis Reveals Varying Neutral Substitution Patterns in Mammalian Evolution," PNAS, 101, 13994-14001.
20. Jensen, J.L. (2005), "Context Dependent DNA Evolutionary Models," Research Report 458, Department of Mathematical Sciences, Aarhus University.
21. Jensen, J.L., and Pedersen, A.K. (2000), "Probabilistic Models of DNA Sequence Evolution with Context Dependent Rates of Substitution," Advances in Applied Probability, 32, 499-517.
22. Kimura, M. (1980), "A Simple Method for Estimating Evolutionary Rate in a Finite Population due to Mutational Production of Neutral and Nearly Neutral Base Substitution through Comparative Studies of Nucleotide Sequences," Journal of Molecular Evolution, 16, 111-120.
23. Lauritzen, S.L. (1996), Graphical Models, Oxford, UK: Clarendon Press.
24. Louis, A.T. (1982), "Finding the Observed Information Matrix When using the EM Algorithm," Journal of the Royal Statistical Society, Ser. B, 44, 226-233.
25. Lunter, G.A., and Hein, J. (2004), "A Nucleotide Substitution Model with Nearest-Neighbour Interactions," Bioinformatics, special issue for ISMB2004, 20, i216-i223.
26. Lunter, G. A., Miklos, I., Drummond, A.J., Jensen, J.L., and Hein, J. (2005), "Bayesian Coestimation of Phylogeny and Sequence Alignment," BMC Bioinformatics, 6, 83.
27. Nielsen, R. (2002), "Mapping Mutations on Phylogenies," Systematic Biology, 51, 729-739.
28. Oehlert, G.W. (1992), "A Note on the Delta Method," The American Statistician, 46, 27-29.
29. Pedersen, J.S., and Hein, J. (2003), "Gene Finding with a Hidden Markov Model of Genome Structure and Evolution," Bioinformatics, 19, 219-227.
30. Redelings, B. D., and Suchard, M.A. (2005), "Joint Bayesian Estimation of Alignment and Phylogeny," Systematic Biology, 54, 401-418.
31. Ren, F., and Yang, Z. (2005), "An Empirical Examination of the Utility of Codon-Substitution Models in Phylogeny Reconstruction," Systematic Biology, 54, 808-818.
32. Robinson, D.M., Jones, D.T., Kishino, H., Goldman, N., and Thorne, J.L. (2003), "Protein Evolution with Dependence Among Codons due to Tertiary Structure," Molecular Biology and Evolution, 20, 1692-1704.
33. Siepel, A., and Haussler, D. (2004), "Phylogenetic Estimation of Context-Dependent Substitution Rates by Maximum Likelihood," Molecular Biology and Evolution, 21, 468-488.
34. Wei, G.C.G., and Tanner, M.A. (1990), "A Monte Carlo Implementation of the EM Algorithm and the Poor Man's Data Augmentation Algorithms," Journal of the American Statistical Association, 85, 699-704.
35. Yap, V.B., and Speed, T.P. (2004), "Modeling DNA Base Substitution in Large Genomic Regions from Two Organisms," Journal of Molecular Evolution, 58, 12-18.
36. (2005), "Estimating Substitution Matrices," in Statistical Methods in Molecular Evolution, ed. R. Nielsen, New York: Springer.