Efficient estimation of pairwise distances between genomes.
BACKGROUND
Genome comparison is central to contemporary genomics and typically relies on sequence alignment. However, genome-wide alignments are difficult to compute. We have, therefore, recently developed an accurate alignment-free estimator of the number of substitutions per site based on the lengths of exact matches between pairs of sequences. The previous implementation of this measure requires n(n-1) suffix tree constructions and traversals, where n is the number of sequences analyzed. This does not scale well for large n.


RESULTS
We present an algorithm to extract pairwise distances in a single traversal of a single suffix tree containing n sequences. As a result, the run time of the suffix tree construction phase of our algorithm is reduced from O(n(2)L) to O(nL), where L is the length of each sequence. We implement this algorithm in the program kr version 2 and apply it to 825 HIV genomes, 13 genomes of enterobacteria and the complete genomes of 12 Drosophila species. We show that, depending on the input dataset, the new program is at least 10 times faster than its predecessor.


BACKGROUND
Version 2 of kr can be tested via a web interface at http://guanine.evolbio.mpg.de/kr2/. It is written in standard C and its source code is available under the GNU General Public License from the same web site.


BACKGROUND
haubold@evolbio.mpg.de Supplementary informations: Supplementary data are available at Bioinformatics online.
Domazet-Loso
Mirjana
M
Department of Evolutionary Genetics, Max-Planck-Institute for Evolutionary Biology, 24306 Plön, Germany.
Haubold
Bernhard
B
D000465:Algorithms; D000818:Animals; D001483:Base Sequence; D030541:Databases, Genetic; D016678:Genome; D023281:Genomics / Q000379:methods*; D006801:Humans; D016415:Sequence Alignment; D017422:Sequence Analysis, DNA