A statistical method for alignment-free comparison of regulatory sequences
1Department of Computer Science, 2Department of Entomology, 3Institute for Genomic Biology, University of Illinois, Urbana-Champaign, Illinois, USA
*To whom correspondence should be addressed.
| Abstract |
|---|
Motivation: The similarity of two biological sequences has traditionally been assessed within the well-established framework of alignment. Here we focus on the task of identifying functional relationships between cis-regulatory sequences that are non-orthologous or greatly diverged. Alignment-free measures of sequence similarity are required in this regime.
Results: We investigate the use of a new score for alignment-free sequence comparison, called the
score. It is based on comparing the frequencies of all fixed-length words in the two sequences. An important, novel feature of the score is that it is comparable across sequence pairs drawn from arbitrary background distributions. We present a method that gives quadratic improvement in the time complexity of calculating the
score, over the naïve method. We then evaluate the score on several tissue-specific families of cis-regulatory modules (in Drosophila and human). The new score is highly successful in discriminating functionally related regulatory sequences from unrelated sequence pairs. The performance of the
score is compared to five other alignment-free similarity measures, and shown to be consistently superior to all of these measures.
Availability: Our implementation of the
score will be made freely available as source code, upon publication of this article, at: http://veda.cs.uiuc.edu/d2z/
Contact: sinhas{at}cs.uiuc.edu
Supplementary information: Supplementary data are available at Bioinformatics online.
This article has been cited by other articles:
![]() |
Q. Dai, Y. Yang, and T. Wang Markov model plus k-word distributions: a synergy that produces novel statistical measures for sequence comparison Bioinformatics, October 15, 2008; 24(20): 2296 - 2302. [Abstract] [Full Text] [PDF] |
||||
