Rapid parallel genome indexing with MapReduce

Menon, R. K., Bhat, G. P., Schatz, M. C. (2011) Rapid parallel genome indexing with MapReduce. MapReduce '11 Proceedings of the second international workshop on MapReduce and its applications . pp. 51-58.

Abstract

Sequence alignment is one of the most important applications in computational biology, and is used for such diverse tasks as identifying homologous proteins, analyzing gene expression, mapping variations between individuals, or assembling de novo the genome of organism. Except for trivial cases involving just a small number of short sequences, virtually all other sequence alignment tasks rely on a precomputed index of the sequence to accelerate the alignment. Two of the most important index structures are the suffix array, which consists of the lexicographically sorted list of suffixes of a genome, and the closely related Burrows-Wheeler Transform (BWT), which is a permutation of the genome based on the suffix array. Constructing these structures on large sequences, such as the human genome, requires several hours of serial computation and must be performed for each genome, or genome assembly, to be analyzed. Here we present a novel parallel algorithm for constructing the suffix array and the BWT of a sequence leveraging the unique features of the MapReduce parallel programming model. We demonstrate the performance of the algorithm by greatly accelerating suffix array and BWT construction on five significant genomes using as many as 120 cores leased from the Amazon Elastic Compute Cloud (EC2), reducing the end-to-end runtime from hours to mere minutes. The source code is available under an open source GPL License at: http://code.google.com/p/genome-indexing/ © 2011 ACM.

Item Type: Paper
Additional Information: Meeting Abstract
Uncontrolled Keywords: burrows-wheeler transform cloud computing DNA sequence analysis MapReduce parallel suffix array construction
Subjects: bioinformatics > genomics and proteomics > alignment
bioinformatics > genomics and proteomics > alignment > sequence alignment
bioinformatics > genomics and proteomics > analysis and processing > Sequence Data Processing
bioinformatics > genomics and proteomics > Mapping and Rendering > Sequence Rendering
bioinformatics > computational biology
Publication Type > Meeting Abstract
CSHL Authors:
Communities: CSHL labs > Schatz lab
Depositing User: CSHL Librarian
Date: 2011
Date Deposited: 16 Mar 2012 14:18
Last Modified: 28 Feb 2018 15:43
URI: https://repository.cshl.edu/id/eprint/25377

Actions (login required)

Administrator's edit/view item Administrator's edit/view item