Abstract:The cost of DNA sequencing has dropped exponentially over the past decade, making genomic data accessible to a growing number of scientists. In bioinformatics, localization of short DNA sequences (reads) within large genomic sequences is commonly facilitated by constructing index data structures which allow for efficient querying of substrings. Recent metagenomic classification pipelines annotate reads with taxonomic labels by analyzing their $k$-mer histograms with respect to a reference genome database. CPU-based index construction is often performed in a preprocessing phase due to the relatively high cost of building irregular data structures such as hash maps. However, the rapidly growing amount of available reference genomes establishes the need for index construction and querying at interactive speeds. In this paper, we introduce MetaCache-GPU -- an ultra-fast metagenomic short read classifier specifically tailored to fit the characteristics of CUDA-enabled accelerators. Our approach employs a novel hash table variant featuring efficient minhash fingerprinting of reads for locality-sensitive hashing and their rapid insertion using warp-aggregated operations. Our performance evaluation shows that MetaCache-GPU is able to build large reference databases in a matter of seconds, enabling instantaneous operability, while popular CPU-based tools such as Kraken2 require over an hour for index construction on the same data. In the context of an ever-growing number of reference genomes, MetaCache-GPU is the first metagenomic classifier that makes analysis pipelines with on-demand composition of large-scale reference genome sets practical. The source code is publicly available at <a class="link-external link-https" href="https://github.com/muellan/metacache" rel="external noopener nofollow">this https URL</a> .

Genetically Improved BarraCUDA

Gene Sequence Alignment on a Public Computing Platform

Implementation of Short Read Alignment Algorithm in OpenCL on Xeon Phi Coprocessor

GPU Accelerated Biological Sequence Alignment

Accelerating Genome-Wide Association Studies Using CUDA Compatible Graphics Processing Units

gpuPairHMM: High-speed Pair-HMM Forward Algorithm for DNA Variant Calling on GPUs

Fast Variation-Aware Read Alignment with Debga-Vara

CGAP-align: a high performance DNA short read alignment tool.

Arioc: High-concurrency short-read alignment on multiple GPUs

Fast Read Alignment with Incorporation of Known Genomic Variants.

GBC: a parallel toolkit based on highly addressable byte-encoding blocks for extremely large-scale genotypes of species

Accelerating Biological Sequence Alignment Algorithm on GPU with CUDA

High-Performance Genomic Analysis Heterogeneous System Using OpenCL

Parallel Accelerated Custom Correlation Coefficient Calculations for Genomics Applications

High-Throughput Sequence Translation Using Cuda

MetaCache-GPU: Ultra-Fast Metagenomic Classification

Pacybara: accurate long-read sequencing for barcoded mutagenized allelic libraries

CUDAMPF++: A Proactive Resource Exhaustion Scheme for Accelerating Homologous Sequence Search on CUDA-enabled GPU

CURC: a CUDA-based reference-free read compressor

CUDA ClustalW: An efficient parallel algorithm for progressive multiple sequence alignment on Multi-GPUs

GSGP-CUDA -- a CUDA framework for Geometric Semantic Genetic Programming