The Community for Technology Leaders
Green Image
Issue No. 02 - April-June (2009 vol. 6)
ISSN: 1545-5963
pp: 281-295
Ana Arribas-Gil , Universidad Carlos III de Madrid, Spain
Dirk Metzler , Johann Wolfgang Goethe-Universität, Germany
Jean-Louis Plouhinec , Institut de Transgénose, CNRS-IEM, France
We present a stochastic sequence evolution model to obtain alignments and estimate mutation rates between two homologous sequences. The model allows two possible evolutionary behaviors along a DNA sequence in order to determine conserved regions and take its heterogeneity into account. In our model, the sequence is divided into slow and fast evolution regions. The boundaries between these sections are not known. It is our aim to detect them. The evolution model is based on a fragment insertion and deletion process working on fast regions only and on a substitution process working on fast and slow regions with different rates. This model induces a pair hidden Markov structure at the level of alignments, thus making efficient statistical alignment algorithms possible. We propose two complementary estimation methods, namely, a Gibbs sampler for Bayesian estimation and a stochastic version of the EM algorithm for maximum likelihood estimation. Both algorithms involve the sampling of alignments. We propose a partial alignment sampler, which is computationally less expensive than the typical whole alignment sampler. We show the convergence of the two estimation algorithms when used with this partial sampler. Our algorithms provide consistent estimates for the mutation rates and plausible alignments and sequence segmentations on both simulated and real data.
Markov processes, probabilistic algorithms, mathematics and statistics, sequence evolution, biology and genetics.

A. Arribas-Gil, J. Plouhinec and D. Metzler, "Statistical Alignment with a Sequence Evolution Model Allowing Rate Heterogeneity along the Sequence," in IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 6, no. , pp. 281-295, 2007.
95 ms
(Ver 3.3 (11022016))