Abstract:Abstract Background Third-generation sequencing offers some advantages over next-generation sequencing predecessors, but with the caveat of harboring a much higher error rate. Clustering-related sequences is an essential task in modern biology. To accurately cluster sequences rich in errors, error type and frequency need to be accounted for. Levenshtein distance is a well-established mathematical algorithm for measuring the edit distance between words and can specifically weight insertions, deletions and substitutions. However, there are drawbacks to using Levenshtein distance in a biological context and hence has rarely been used for this purpose. We present novel modifications to the Levenshtein distance algorithm to optimize it for clustering error-rich biological sequencing data. Results We successfully introduced a bidirectional frameshift allowance with end-user determined accommodation caps combined with weighted error discrimination. Furthermore, our modifications dramatically improved the computational speed of Levenstein distance. For simulated ONT MinION and PacBio Sequel datasets, the average clustering sensitivity for 3GOLD was 41.45% (S.D. 10.39) higher than Sequence-Levenstein distance, 52.14% (S.D. 9.43) higher than Levenshtein distance, 55.93% (S.D. 8.67) higher than Starcode, 42.68% (S.D. 8.09) higher than CD-HIT-EST and 61.49% (S.D. 7.81) higher than DNACLUST. For biological ONT MinION data, 3GOLD clustering sensitivity was 27.99% higher than Sequence-Levenstein distance, 52.76% higher than Levenshtein distance, 56.39% higher than Starcode, 48% higher than CD-HIT-EST and 70.4% higher than DNACLUST. Conclusion Our modifications to Levenshtein distance have improved its speed and accuracy compared to the classic Levenshtein distance, Sequence-Levenshtein distance and other commonly used clustering approaches on simulated and biological third-generation sequenced datasets. Our clustering approach is appropriate for datasets of unknown cluster centroids, such as those generated with unique molecular identifiers as well as known centroids such as barcoded datasets. A strength of our approach is high accuracy in resolving small clusters and mitigating the number of singletons.

3GOLD: optimized Levenshtein distance for clustering third-generation sequencing data

Gene Sequence Alignment on a Public Computing Platform

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Interpreting Sequence-Levenshtein distance for determining error type and frequency between two embedded sequences of equal length

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

Efficient Frequency-Based De Novo Short-Read Clustering for Error Trimming in Next-Generation Sequencing

cascAGS: Comparative Analysis of SNP Calling Methods for Human Genome Data in the Absence of Gold Standard

De novo clustering of extensive long-read transcriptome datasets with isONclust3

SEED: efficient clustering of next-generation sequences.

Clustering pipeline for determining consensus sequences in targeted next-generation sequencing

Highly efficient clustering of long-read transcriptomic data with GeLuster

An Efficient Greedy Incremental Sequence Clustering Algorithm

Linear normalised hash function for clustering gene sequences and identifying reference sequences from multiple sequence alignments

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

A feasible roadmap to identifying significant intercellular genomic heterogeneity in deep sequencing data

Distangsd: Fast and Accurate Inference of Genetic Distances for Next-Generation Sequencing Data

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

A Crowdsourcing Method For Correcting Sequencing Errors For The Third-Generation Sequencing Data

GradHC: Highly Reliable Gradual Hash-based Clustering for DNA Storage Systems

Gclust:A Parallel Clustering Tool for Microbial Genomic Data

CLEVER: Clique-Enumerating Variant Finder