Abstract:Genomic sequence alignment is an important research topic in bioinformatics and continues to attract significant efforts. As genomic data grow exponentially, however, most of alignment methods face challenges due to their huge computational costs. HMMER, a suite of bioinformatics tools, is widely used for the analysis of homologous protein and nucleotide sequences with high sensitivity, based on profile hidden Markov models (HMMs). Its latest version, HMMER3, introdues a heuristic pipeline to accelerate the alignment process, which is carried out on central processing units (CPUs) with the support of streaming SIMD extensions (SSE) instructions. Few acceleration results have since been reported based on HMMER3. In this paper, we propose a five-tiered parallel framework, CUDAMPF++, to accelerate the most computationally intensive stages of HMMER3's pipeline, multiple/single segment Viterbi (MSV/SSV), on a single graphics processing unit (GPU). As an architecture-aware design, the proposed framework aims to fully utilize hardware resources via exploiting finer-grained parallelism (multi-sequence alignment) compared with its predecessor (CUDAMPF). In addition, we propose a novel method that proactively sacrifices L1 Cache Hit Ratio (CHR) to get improved performance and scalability in return. A comprehensive evaluation shows that the proposed framework outperfroms all existig work and exhibits good consistency in performance regardless of the variation of query models or protein sequence datasets. For MSV (SSV) kernels, the peak performance of the CUDAMPF++ is 283.9 (471.7) GCUPS on a single K40 GPU, and impressive speedups ranging from 1.x (1.7x) to 168.3x (160.7x) are achieved over the CPU-based implementation (16 cores, 32 threads).

The accelerating implementation of BLAST with stream processor

Research on Parallel BLAST Algorithm Based on Multi-core Stream Processors

Gene Sequence Alignment on a Public Computing Platform

Large Scale DNA sequence alignment and Kernel method implemented with GPUs

GPU Accelerated Biological Sequence Alignment

Exploiting Parallelization of BLAST on Dawning 4000A

Fast Alignment of Biological Sequences Based on General Purpose Graphic Processor Unit

Mega-base Biological Sequence Alignment Targeting OpenCL Architecture

Accelerating Biological Sequence Alignment Algorithm on GPU with CUDA

Exploiting the Task-Pipelined Parallelism of Stream Programs on Many-Core GPUs.

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

AGAThA: Fast and Efficient GPU Acceleration of Guided Sequence Alignment for Long Read Mapping

Distributed Sequence Alignment Applications for the Public Computing Architecture

BLVector: Fast BLAST-Like Algorithm for Manycore CPU With Vectorization

BLASX: A High Performance Level-3 BLAS Library for Heterogeneous Multi-GPU Computing

Massively Parallel Implementation of Sequence Alignment with Basic Local Alignment Search Tool Using Parallel Computing in Java Library

Stream-K++: Adaptive GPU GEMM Kernel Scheduling and Selection using Bloom Filters

ScalaBLAST: A Scalable Implementation of BLAST for High-Performance Data-Intensive Bioinformatics Analysis

High Performance Computing Via a GPU

Parallelization of BLAST with MapReduce for Long Sequence Alignment.

Accelerating BWA-MEM Read Mapping on GPUs