Optimized pipeline of MuTect and GATK tools to improve the detection of somatic single nucleotide polymorphisms in whole-exome sequencing data

Abstract

BackgroundDetecting somatic mutations in whole exome sequencing data of cancer samples has become a popular approach for profiling cancer development, progression and chemotherapy resistance. Several studies have proposed software packages, filters and parametrizations. However, many research groups reported low concordance among different methods. We aimed to develop a pipeline which detects a wide range of single nucleotide mutations with high validation rates. We combined two standard tools – Genome Analysis Toolkit (GATK) and MuTect – to create the GATK-LODN method. As proof of principle, we applied our pipeline to exome sequencing data of hematological (Acute Myeloid and Acute Lymphoblastic Leukemias) and solid (Gastrointestinal Stromal Tumor and Lung Adenocarcinoma) tumors. We performed experiments on simulated data to test the sensitivity and specificity of our pipeline.ResultsThe software MuTect presented the highest validation rate (90 %) for mutation detection, but limited number of somatic mutations detected. The GATK detected a high number of mutations but with low specificity. The GATK-LODN increased the performance of the GATK variant detection (from 5 of 14 to 3 of 4 confirmed variants), while preserving mutations not detected by MuTect. However, GATK-LODN filtered more variants in the hematological samples than in the solid tumors. Experiments in simulated data demonstrated that GATK-LODN increased both specificity and sensitivity of GATK results.ConclusionWe presented a pipeline that detects a wide range of somatic single nucleotide variants, with good validation rates, from exome sequencing data of cancer samples. We also showed the advantage of combining standard algorithms to create the GATK-LODN method, that increased specificity and sensitivity of GATK results. This pipeline can be helpful in discovery studies aimed to profile the somatic mutational landscape of cancer genomes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1190-7) contains supplementary material, which is available to authorized users.

PDX models recapitulate the genetic and epigenetic landscape of pediatric T‐cell leukemia

“…Briefly, alignment was performed with bwa (Burrows‐Wheeler Aligner; Li & Durbin, ) mem algorithm and followed by removing duplicate reads (Picard) and local realignment by Genome Analysis Toolkit (GATK; McKenna et al , ). Mutation callers Mutect (do Valle et al , ) and Strelka (Saunders et al , ) were next applied on the recalibrated data for detection of somatic mutations. The presence of single nucleotide variants (SNV) and small insertions and deletions (InDels) in the germ line was excluded by analysis of remission samples.…”

Mutations in the coat complex II component SEC23B promote colorectal cancer metastasis

“…Somatic mutations were called with Genome Analysis Toolkit (GATK) 18 (https://software.broadinstitute. org/gatk/) best practice pipeline and Mutect 19 . Variants were annotated using ANNOVAR 20 (http://www.…”

Informatics for Cancer Immunotherapy

“…Most labs still manually examine somatic mutation calls with tools like IGV [64] or pileup.js [65] to confirm their accuracy, as in [32] . Whole pipelines have been published that can serve as a model for investigators setting up their own computational programs [66][67][68] . As of 2017, there is no consensus on the most sensitive and specific pipeline.…”