A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset

Zhou, Yong, Kathiresan, Nagarajan, Yu, Zhichao, Rivera, Luis F, Yang, Yujian, Thimma, Manjula, Manickam, Keerthana, Chebotarov, Dmytro, Mauleon, Ramil, Chougule, Kapeel, Wei, Sharon, Gao, Tingting, Green, Carl D, Zuccolo, Andrea, Xie, Weibo, Ware, Doreen, Zhang, Jianwei, McNally, Kenneth L, Wing, Rod A (January 2024) A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset. BMC Biology, 22 (1). p. 13. ISSN 1741-7007

	PDF s12915-024-01820-5.pdf - Published Version Download (5MB)
	PDF s12915-024-01820-5.pdf Available under License Creative Commons Attribution. Download (5MB)

Abstract

BACKGROUND: Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable. RESULTS: Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq). CONCLUSIONS: This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

Item Type:	Paper
Subjects:	bioinformatics bioinformatics > genomics and proteomics > genetics & nucleic acid processing > DNA, RNA structure, function, modification bioinformatics > genomics and proteomics > genetics & nucleic acid processing bioinformatics > genomics and proteomics bioinformatics > genomics and proteomics > genetics & nucleic acid processing > genomes organism description > plant bioinformatics > genomics and proteomics > genetics & nucleic acid processing > DNA, RNA structure, function, modification > single nucleotide polymorphism
CSHL Authors:	Chougule, Kapeel Wei, Sharon Ware, Doreen H.
Communities:	CSHL labs > Ware lab
SWORD Depositor:	CSHL Elements
Depositing User:	CSHL Elements
Date:	25 January 2024
Date Deposited:	29 Jan 2024 13:53
Last Modified:	29 Jan 2024 13:53
PMCID:	PMC10809545
Related URLs:	Publisher
URI:	https://repository.cshl.edu/id/eprint/41423

Actions (login required)

Administrator's edit/view item