Issue No. 02 - April-June (2009 vol. 6)
DOI Bookmark: http://doi.ieeecomputersociety.org/10.1109/TCBB.2009.9
Todd J. Treangen , Institut Pasteur, CNRS and UPMC University, France
Aaron E. Darling , University of California-Davis, Davis
Guillaume Achaz , UPMC University, CNRS, MNHN, IRD, UMR, Paris
Mark A. Ragan , University of Queensland and ARC Centre of Excellence, Brisbane
Xavier Messeguer , Technical University of Catalonia, Barcelona
Eduardo P.C. Rocha , Institut Pasteur, CNRS and UPMC University, France
Pairwise local sequence alignment methods have been the prevailing technique to identify homologous nucleotides between related species. However, existing methods that identify and align all homologous nucleotides in one or more genomes have suffered from poor scalability and limited accuracy. We propose a novel method that couples a gapped extension heuristic with an efficient filtration method for identifying interspersed repeats in genome sequences. During gapped extension, we use the MUSCLE implementation of progressive global multiple alignment with iterative refinement. The resulting gapped extensions potentially contain alignments of unrelated sequence. We detect and remove such undesirable alignments using a hidden Markov model (HMM) to predict the posterior probability of homology. The HMM emission frequencies for nucleotide substitutions can be derived from any time-reversible nucleotide substitution matrix. We evaluate the performance of our method and previous approaches on a hybrid data set of real genomic DNA with simulated interspersed repeats. Our method outperforms a related method in terms of sensitivity, positive predictive value, and localizing boundaries of homology. The described methods have been implemented in freely available software, Repeatoire, available from: http://wwwabi.snv.jussieu.fr/public/Repeatoire.
Sequence alignment, genome comparison, DNA repeats, local multiple alignment, hidden Markov model, gapped extension.
T. J. Treangen, G. Achaz, M. A. Ragan, X. Messeguer, E. P. Rocha and A. E. Darling, "A Novel Heuristic for Local Multiple Alignment of Interspersed DNA Repeats," in IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 6, no. , pp. 180-189, 2009.