Abstract:Background: Clustering is a fundamental operation in the analysis of biological sequence data. New DNA sequencing technologies have dramatically increased the rate at which we can generate data, resulting in datasets that cannot be efficiently analyzed by traditional clustering methods.This is particularly true in the context of taxonomic profiling of microbial communities through direct sequencing of phylogenetic markers (e.g. 16S rRNA) - the domain that motivated the work described in this paper. Many analysis approaches rely on an initial clustering step aimed at identifying sequences that belong to the same operational taxonomic unit (OTU). When defining OTUs (which have no universally accepted definition), scientists must balance a trade-off between computational efficiency and biological accuracy, as accurately estimating an environment's phylogenetic composition requires computationally-intensive analyses. We propose that efficient and mathematically well defined clustering methods can benefit existing taxonomic profiling approaches in two ways: (i) the resulting clusters can be substituted for OTUs in certain applications; and (ii) the clustering effectively reduces the size of the data-sets that need to be analyzed by complex phylogenetic pipelines (e.g., only one sequence per cluster needs to be provided to downstream analyses). Results: To address the challenges outlined above, we developed DNACLUST, a fast clustering tool specifically designed for clustering highly-similar DNA sequences.Given a set of sequences and a sequence similarity threshold, DNACLUST creates clusters whose radius is guaranteed not to exceed the specified threshold. Underlying DNACLUST is a greedy clustering strategy that owes its performance to novel sequence alignment and k-mer based filtering algorithms.DNACLUST can also produce multiple sequence alignments for every cluster, allowing users to manually inspect clustering results, and enabling more detailed analyses of the clustered data. Conclusions: We compare DNACLUST to two popular clustering tools: CD-HIT and UCLUST. We show that DNACLUST is about an order of magnitude faster than CD-HIT and UCLUST (exact mode) and comparable in speed to UCLUST (approximate mode). The performance of DNACLUST improves as the similarity threshold is increased (tight clusters) making it well suited for rapidly removing duplicates and near-duplicates from a dataset, thereby reducing the size of the data being analyzed through more elaborate approaches.

CD-HIT: accelerated for clustering the next-generation sequencing data

Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences

Gene Sequence Alignment on a Public Computing Platform

Gclust:A Parallel Clustering Tool for Microbial Genomic Data

An Efficient Greedy Incremental Sequence Clustering Algorithm

Dace: A Scalable Dp-Means Algorithm for Clustering Extremely Large Sequence Data

An efficient parallel algorithm for multiple sequence similarities calculation using a low complexity method.

DGCF: A Distributed Greedy Clustering Framework for Large-scale Genomic Sequences.

TreQ-CG: Clustering Accelerates High-Throughput Sequencing Read Mapping

SEED: efficient clustering of next-generation sequences.

EdClust: A Heuristic Sequence Clustering Method with Higher Sensitivity

Accelerating Gene Clustering on Heterogeneous Clusters

Sequence Clustering in Bioinformatics: an Empirical Study.

DNACLUST: accurate and efficient clustering of phylogenetic marker genes

Comparison of Methods for Biological Sequence Clustering.

SaAlign: Multiple DNA/RNA Sequence Alignment and Phylogenetic Tree Construction Tool for Ultra-Large Datasets and Ultra-Long Sequences Based on Suffix Array

Ngia: A Novel Greedy Incremental Alignment Based Algorithm for Gene Sequence Clustering

A data parallel strategy for aligning multiple biological sequences on multi-core computers

Cluster-Distribute-Align-Merge: A General Algorithm to Speed Up Multiple Sequence Alignment on Multi-Core Computers

Interpolative Multidimensional Scaling Techniques for the Identification of Clusters in Very Large Sequence Sets

ClustalXeed: a GUI-based grid computation version for high performance and terabyte size multiple sequence alignment