An Eulerian path approach to local multiple alignment for DNA sequences

Yu Zhang, Michael S. Waterman

Research output: Contribution to journalArticle

22 Scopus citations

Abstract

Expensive computation in handling a large number of sequences limits the application of local multiple sequence alignment. We present an Eulerian path approach to local multiple alignment for DNA sequences. The computational time and memory usage of this approach is approximately linear to the total size of sequences analyzed; hence, it can handle thousands of sequences or millions of letters simultaneously. By constructing a De Bruijn graph, most of the conserved segments are amplified as heavy Eulerian paths in the graph, and the original patterns distributed in sequences are recovered even if they do not exist in any single sequence. This approach can accurately detect unknown conserved regions, for both short and long, conserved and degenerate patterns. We further present a Poisson heuristic to estimate the significance of a local multiple alignment. The performance of our method is demonstrated by finding Alu repeats in the human genome. We compare the results with Alus marked by REPEATMASKER, where the two programs are in good agreement. Our method is robust under various conditions and superior to other methods in terms of efficiency and accuracy.

Original languageEnglish (US)
Pages (from-to)1285-1290
Number of pages6
JournalProceedings of the National Academy of Sciences of the United States of America
Volume102
Issue number5
DOIs
StatePublished - Feb 1 2005

    Fingerprint

All Science Journal Classification (ASJC) codes

  • General

Cite this