Fast comparison of DNA sequences by oligonucleotide profiling

Vicente Arnau¹ , Miguel Gallach² and Ignacio Marín³

¹Departamento de Informática. Universidad de Valencia, Spain

²Departamento de Genética. Universidad de Valencia, Spain

³Instituto de Biomedicina de Valencia. Consejo Superior de Investigaciones Científicas (IBV-CSIC), Spain

author email corresponding author email

BMC Research Notes 2008, 1:5doi:10.1186/1756-0500-1-5


Published:	28 February 2008

Abstract

Background

The comparison of DNA sequences is a traditional problem in genomics and bioinformatics. Many new opportunities emerge due to the improvement of personal computers, allowing the implementation of novel strategies of analysis.

Findings

We describe a new program, called UVWORD, which determines the number of times that each DNA word present in a sequence (target) is found in a second sequence (source), a procedure that we have called oligonucleotide profiling. On a standard computer, the user may search for words of a size ranging from k = 1 to k = 14 nucleotides. Average counts for groups of contiguous words may also be established. The rate of analysis on standard computers is from 3.4 (k = 14) to 16 millions of words per second (1 ≤ k ≤ 8). This makes feasible the fast screening of even the longest known DNA molecules.

Discussion

We show that the combination of the ability of analyzing words of relatively long size, which occur very rarely by chance, and the fast speed of the program allows to perform novel types of screenings, complementary to those provided by standard programs such as BLAST. This method can be used to determine oligonucleotide content, to characterize the distribution of repetitive sequences in chromosomes, to determine the evolutionary conservation of sequences in different species, to establish regions of similar DNA among chromosomes or genomes, etc.