Evaluating the use of GPUs in liver image segmentation and HMMER database searches

In this paper we present the results of parallelizing two life sciences applications, Markov random fields-based (MRF) liver segmentation and HMMER's Viterbi algorithm, using GPUs. We relate our experiences in porting both applications to the GPU as well as the techniques and optimizations that are most beneficial. The unique characteristics of both algorithms are demonstrated by implementations on an NVIDIA 8800 GTX Ultra using the CUDA programming environment. We test multiple enhancements in our GPU kernels in order to demonstrate the effectiveness of each strategy. Our optimized MRF kernel achieves over 130× speedup, and our hmmsearch implementation achieves up to 38× speedup. We show that the differences in speedup between MRF and hmmsearch is due primarily to the frequency at which the hmmsearch must read from the GPU's DRAM.

[1]  Pat Hanrahan,et al.  ClawHMMER: A Streaming HMMer-Search Implementation , 2005, SC.

[2]  Sean R. Eddy,et al.  Biological Sequence Analysis: Probabilistic Models of Proteins and Nucleic Acids , 1998 .

[3]  Pat Hanrahan,et al.  Brook for GPUs: stream computing on graphics hardware , 2004, SIGGRAPH 2004.

[4]  Weiguo Liu,et al.  Bio-sequence database scanning on a GPU , 2006, Proceedings 20th IEEE International Parallel & Distributed Processing Symposium.

[5]  Kurt Keutzer,et al.  Data-Parallel Large Vocabulary Continuous Speech Recognition on Graphics Processors , 2008 .

[6]  Binsheng Zhao,et al.  Liver segmentation for CT images using GVF snake. , 2005, Medical physics.

[7]  Vipin Chaudhary,et al.  Design of a benchmark dataset, similarity metrics, and tools for liver segmentation , 2008, SPIE Medical Imaging.

[8]  J.J. Rodriguez,et al.  Snake-based liver lesion segmentation , 2004, 6th IEEE Southwest Symposium on Image Analysis and Interpretation, 2004..

[9]  Kuni Ohtomo,et al.  Liver CT image processing: a short introduction of the technical elements. , 2006, European journal of radiology.

[10]  Klaus D. Tönnies,et al.  A New Approach for Model-Based Adaptive Region Growing in Medical Image Analysis , 2001, CAIP.

[11]  Pat Hanrahan,et al.  ClawHMMER: A Streaming HMMer-Search Implementatio , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[12]  Bashar Qudah,et al.  Accelerating the HMMER sequence analysis suite using conventional processors , 2006, 20th International Conference on Advanced Information Networking and Applications - Volume 1 (AINA'06).

[13]  DePaul University Cphilips.  Segmentation of Soft Tissue Using Texture Features and Gradient Snakes , 2008 .

[14]  Shaohui Huang,et al.  Using GVF Snake to Segment Liver from CT Images , 2006, 2006 3rd IEEE/EMBS International Summer School on Medical Devices and Biosensors.

[15]  Kunio Doi,et al.  Automated hepatic volumetry for living related liver transplantation at multisection CT. , 2006, Radiology.

[16]  Jacob D. Furst,et al.  A comparison of texture models for automatic liver segmentation , 2007, SPIE Medical Imaging.

[17]  Demetri Terzopoulos,et al.  Deformable models in medical image analysis: a survey , 1996, Medical Image Anal..

[18]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[19]  Bertil Schmidt,et al.  Accelerating the Viterbi Algorithm for Profile Hidden Markov Models Using Reconfigurable Hardware , 2006, International Conference on Computational Science.

[20]  Sean R. Eddy,et al.  Profile hidden Markov models , 1998, Bioinform..

[21]  Brandon Harris,et al.  Accelerator design for protein sequence HMM search , 2006, ICS '06.

[22]  Brian E. Smith,et al.  An Efficient Parallel Implementation of the Hidden Markov Methods for Genomic Sequence-Search on a Massively Parallel System , 2008, IEEE Transactions on Parallel and Distributed Systems.

[23]  E-Liang Chen,et al.  An automatic diagnostic system for CT liver image classification , 1998, IEEE Transactions on Biomedical Engineering.