Related Experiment Video
Updated: Jul 4, 2025
Author Spotlight: Investigating the Role of Repetitive DNA Misregulation in Cancer Initiation and Immunotherapy Resistance
Published on: December 13, 2024
A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset
Yong Zhou1,2, Nagarajan Kathiresan3, Zhichao Yu1,4
1Center for Desert Agriculture (CDA), Biological and Environmental Sciences & Engineering Division (BESE), King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia.
A new open-source workflow (HPC-GVCW) enables faster, affordable SNP calling using GATK on various platforms. This accelerates genetic variation studies across multiple crop species, discovering millions of novel SNPs for functional research.
Area of Science:
- Genomics
- Bioinformatics
- Computational Biology
Background:
- Single-nucleotide polymorphisms (SNPs) are crucial for molecular genetic variation studies.
- Exponential growth in reference genomes and resequencing data necessitates faster SNP calling tools.
- Existing Genome Analysis Toolkit (GATK) versions lack affordable, high-performance computing (HPC) accessibility.
Purpose of the Study:
- To develop an open-source, high-performance computing genome variant calling workflow (HPC-GVCW) for GATK.
- To enable efficient SNP calling across diverse computing platforms, from supercomputers to desktops.
- To accelerate SNP discovery and facilitate functional and breeding studies.
Main Methods:
- Developed and implemented the HPC-GVCW pipeline for GATK.
- Benchmarked HPC-GVCW performance and accuracy on multiple crop species.
- Applied HPC-GVCW in production mode using a subpopulation-aware 16-genome rice reference panel with ~3000 resequenced accessions.
Main Results:
- HPC-GVCW demonstrated comparable performance and accuracy to GATK alone across crop species.
- SNP calling on the rice panel took ~16 weeks, identifying an average of 27.3 million SNPs/genome.
- Discovered ~2.3 million novel SNPs in rice, not present in the IRGSP RefSeq reference genome.
Conclusions:
- The open-source HPC-GVCW pipeline significantly enhances SNP calling speed on HPC platforms.
- The workflow is broadly applicable, successfully tested on four major crop species with varying genome sizes.
- Combining HPC-GVCW with subpopulation-aware reference panels enables rapid SNP discovery and public release, revealing functionally relevant novel SNPs.
