An alignment-free heuristic for fast sequence comparisons with applications to phylogeny reconstruction.

BACKGROUND

Alignment-free methods for sequence comparisons have become popular in many bioinformatics applications, specifically in the estimation of sequence similarity measures to construct phylogenetic trees. Recently, the average common substring measure, ACS, and its k-mismatch counterpart, ACS, have been shown to produce results as effective as multiple-sequence alignment based methods for reconstruction of phylogeny trees. Since computing ACS takes O(n logkn) time and hence impractical for large datasets, multiple heuristics that can approximate ACS have been introduced.

RESULTS

In this paper, we present a novel linear-time heuristic to approximate ACS, which is faster than computing the exact ACS while being closer to the exact ACS values compared to previously published linear-time greedy heuristics. Using four real datasets, containing both DNA and protein sequences, we evaluate our algorithm in terms of accuracy, runtime and demonstrate its applicability for phylogeny reconstruction. Our algorithm provides better accuracy than previously published heuristic methods, while being comparable in its applications to phylogeny reconstruction.

CONCLUSIONS

Our method produces a better approximation for ACS and is applicable for the alignment-free comparison of biological sequences at highly competitive speed. The algorithm is implemented in Rust programming language and the source code is available at https://github.com/srirampc/adyar-rs .