Finding similar regions in many sequences

  • Authors:
  • Ming Li;Bin Ma;Lusheng Wang

  • Affiliations:
  • Department of Computer Science, University of California, Santa Barbara, California;Department of Computer Science, University of Western Ontario, London, Ontario N6A5B7, Canada;City University of Hong Kong, Kowloon, Hong Kong

  • Venue:
  • Journal of Computer and System Sciences - STOC 1999
  • Year:
  • 2002

Quantified Score

Hi-index 0.00

Visualization

Abstract

Algorithms for finding similar, or highly conserved, regions in a group of sequences are at the core of many molecular biology problems. Assume that we are given n DNA sequences s1, ...., sn. The Consensus Patterns problem, which has been widely studied in bioinformatics research, in its simplest form, asks for a region of length L in each si, and a median string s of length L so that the total Hamming distance from s to these regions is minimized. We show that the problem is NP-hard and give a polynomial time approximation scheme (PTAS) for it. We then present an efficient approximation algorithm for the consensus pattern problem under the original relative entropy measure. As an interesting application of our analysis, we further obtain a PTAS for a restricted (but still NP-hard) version of the important consensus alignment problem allowing at most constant number of gaps, each of arbitrary length, in each sequence.