Abstract
Next-generation sequencing technologies have led to the high-throughput production of sequence data (reads) at low cost. However, these reads are significantly shorter and more error-prone than conventional Sanger shotgun reads. This poses a challenge for the de novo assembly in terms of assembly quality and scalability for large-scale short read datasets. We present DecGPU, the first parallel and distributed error correction algorithm for high-throughput short reads (HTSRs) using a hybrid combination of CUDA and MPI parallel programming models. DecGPU provides CPU-based and GPU-based versions, where the CPU-based version employs coarse-grained and fine-grained parallelism using the MPI and OpenMP parallel programming models, and the GPU-based version takes advantage of the CUDA and MPI parallel programming models and employs a hybrid CPU+GPU computing model to maximize the performance by overlapping the CPU and GPU computation. The distributed feature of our algorithm makes it feasible and flexible for the error correction of large-scale HTSR datasets. Using simulated and real datasets, our algorithm demonstrates superior performance, in terms of error correction quality and execution speed, to the existing error correction algorithms. Furthermore, when combined with Velvet and ABySS, the resulting DecGPU-Velvet and DecGPU-ABySS assemblers demonstrate the potential of our algorithm to improve de novo assembly quality for de-Bruijn-graph-based assemblers. DecGPU is publicly available open-source software, written in CUDA C++ and MPI. The experimental results suggest that DecGPU is an effective and feasible error correction algorithm to tackle the flood of short reads produced by next-generation sequencing technologies.
MeSH Terms
Algorithms
Computer Simulation
Reproducibility of Results
Sequence Analysis, DNA/methods
Software
Authors & Affiliations
3 authors, click to expand affiliations / ORCID
Liu Yongchao
School of Computer Engineering, Nanyang Technological University, 639798, Singapore. liuy0039@ntu.edu.sg
Schmidt Bertil
Maskell Douglas L
References (20)
20 references, click to expand
-
Ultrafast and memory-efficient alignment of short DNA sequences to the human genome.
Genome Biol. 2009;10(3):R25
PMID: 19261174
-
Correction of sequencing errors in a mixed set of reads.
Bioinformatics. 2010 May 15;26(10):1284-90
PMID: 20378555
-
Velvet: algorithms for de novo short read assembly using de Bruijn graphs.
Genome Res. 2008 May;18(5):821-9
PMID: 18349386
-
A fast hybrid short read fragment assembly algorithm.
Bioinformatics. 2009 Sep 1;25(17):2279-80
PMID: 19535537
-
An Eulerian path approach to DNA fragment assembly.
Proc Natl Acad Sci U S A. 2001 Aug 14;98(17):9748-53
PMID: 11504945
-
ARACHNE: a whole-genome shotgun assembler.
Genome Res. 2002 Jan;12(1):177-89
PMID: 11779843
-
A parallel algorithm for error correction in high-throughput short-read data on CUDA-enabled graphics hardware.
J Comput Biol. 2010 Apr;17(4):603-15
PMID: 20426693
-
SHREC: a short-read error correction method.
Bioinformatics. 2009 Sep 1;25(17):2157-63
PMID: 19542152
-
ABySS: a parallel assembler for short read sequence data.
Genome Res. 2009 Jun;19(6):1117-23
PMID: 19251739
-
A whole-genome assembly of Drosophila.
Science. 2000 Mar 24;287(5461):2196-204
PMID: 10731133
-
Assembling millions of short DNA sequences using SSAKE.
Bioinformatics. 2007 Feb 15;23(4):500-1
PMID: 17158514
-
SHARCGS, a fast and highly accurate short-read assembly algorithm for de novo genomic sequencing.
Genome Res. 2007 Nov;17(11):1697-706
PMID: 17908823
-
Extending assembly of short DNA sequences to handle error.
Bioinformatics. 2007 Nov 1;23(21):2942-4
PMID: 17893086
-
Short read fragment assembly of bacterial genomes.
Genome Res. 2008 Feb;18(2):324-30
PMID: 18083777
-
CUDASW++2.0: enhanced Smith-Waterman protein database search on CUDA-enabled GPUs based on SIMT and virtualized SIMD abstractions.
BMC Res Notes. 2010 Apr 06;3:93
PMID: 20370891
-
CUDASW++: optimizing Smith-Waterman sequence database searches for CUDA-enabled graphics processing units.
BMC Res Notes. 2009 May 06;2:73
PMID: 19416548
-
De novo assembly of human genomes with massively parallel short read sequencing.
Genome Res. 2010 Feb;20(2):265-72
PMID: 20019144
-
PCAP: a whole-genome assembly program.
Genome Res. 2003 Sep;13(9):2164-70
PMID: 12952883
-
ALLPATHS: de novo assembly of whole-genome shotgun microreads.
Genome Res. 2008 May;18(5):810-20
PMID: 18340039
-
The Atlas genome assembly system.
Genome Res. 2004 Apr;14(4):721-32
PMID: 15060016