Group-based variant calling leveraging next-generation supercomputing for large-scale whole-genome sequencing studies
Authors:
- Kristopher A. StandishKSKristopher A. Standish
- Tristan M. CarlandTCTristan M. Carland
- Glenn K. LockwoodGLGlenn K. Lockwood
- Wayne PfeifferWPWayne Pfeiffer
- Mahidhar TatineniMTMahidhar Tatineni
- C Chris HuangCHC Chris Huang
- Sarah LamberthSLSarah Lamberth
- Yauheniya CherkasYCYauheniya Cherkas
- Carrie BrodmerkelCBCarrie Brodmerkel
- Ed JaegerEJEd Jaeger
- Lance SmithLSLance Smith
- Gunaretnam RajagopalGRGunaretnam Rajagopal
Open Access
- 22 September 2015
- journal article
- research article
- Published by Springer Nature in BMC Bioinformatics
- Vol. 16 (1) , 1-14
Abstract
Next-generation sequencing (NGS) technologies have become much more efficient, allowing whole human genomes to be sequenced faster and cheaper than ever before. However, processing the raw sequence reads associated with NGS technologies requires care and sophistication in order to draw compelling inferences about phenotypic consequences of variation in human genomes. It has been shown that different approaches to variant calling from NGS data can lead to different conclusions. Ensuring appropriate accuracy and quality in variant calling can come at a computational cost. We describe our experience implementing and evaluating a group-based approach to calling variants on large numbers of whole human genomes. We explore the influence of many factors that may impact the accuracy and efficiency of group-based variant calling, including group size, the biogeographical backgrounds of the individuals who have been sequenced, and the computing environment used. We make efficient use of the Gordon supercomputer cluster at the San Diego Supercomputer Center by incorporating job-packing and parallelization considerations into our workflow while calling variants on 437 whole human genomes generated as part of large association study. We ultimately find that our workflow resulted in high-quality variant calls in a computationally efficient manner. We argue that studies like ours should motivate further investigations combining hardware-oriented advances in computing systems with algorithmic developments to tackle emerging ‘big data’ problems in biomedical research brought on by the expansion of NGS technologies.Keywords
This publication has 17 references indexed in Scilit:
- Radiographic benefit and maintenance of clinical benefit with intravenous golimumab therapy in patients with active rheumatoid arthritis despite methotrexate therapy: results up to 1 year of the phase 3, randomised, multicentre, double blind, placebo controlled GO-FURTHER trialAnnals of the Rheumatic Diseases, 2014
- A Method for Inferring an Individual’s Genetic Ancestry and Degree of Admixture Associated with Six Major Continental PopulationsFrontiers in Genetics, 2013
- Clinical Implications of Human Population Differences in Genome-Wide Rates of Functional GenotypesFrontiers in Genetics, 2012
- Five Years of GWAS DiscoveryThe American Journal of Human Genetics, 2012
- A framework for variation discovery and genotyping using next-generation DNA sequencing dataNature Genetics, 2011
- The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing dataGenome Research, 2010
- The Sequence Alignment/Map format and SAMtoolsBioinformatics, 2009
- Fast and accurate short read alignment with Burrows–Wheeler transformBioinformatics, 2009
- The Sequence of the Human GenomeScience, 2001
- Initial sequencing and analysis of the human genomeNature, 2001