Bioinformatics Advance Access published online on September 5, 2007
Bioinformatics, doi:10.1093/bioinformatics/btm422
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Multiple spaced seeds for homology search
aDepartment of Computer Science, University of Western Ontario, N6A 5B7, London, Ontario, CANADA, bNumerical Analysis, Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 Lund, SWEDEN
*To whom correspondence should be addressed. Dr. Silvana Ilie, E-mail: ilie{at}csd.uwo.ca
| Abstract |
|---|
Motivation: Homology search finds similar segments between two biological sequences, such as DNA or protein sequences. The introduction of optimal spaced seeds in PatternHunter, (Ma et al., 2002), has increased both the sensitivity and the speed of homology search and it has been adopted by many alignment programs such as BLAST. With the further improvement provided by multiple spaced seeds in PatternHunterII, (Li et al., 2004), Smith-Waterman sensitivity is approached at BLASTn speed. However, computing optimal multiple spaced seeds was proved to be NP-hard and current heuristic algorithms are all very slow (exponential).
Results: We give a simple algorithm which computes good multiple seeds in polynomial time. Due to a completely different approach, the difference with respect to the previous methods is dramatic. The multiple spaced seed of PatternHunterII, with 16 weight 11 seeds, (Li et al., 2004), was computed in 12 days. It takes us 17 seconds to find a better one. Our approach changes the way of looking at multiple spaced seeds.
Contact: ilie{at}csd.uwo.ca
Associate Editor: Prof. John Quackenbush
Received on May 14, 2007; revised on July 19, 2007; accepted on August 11, 2007