Abstract
Next generation sequencing is extensively applied to catalogue somatic mutations in cancer, in research settings and increasingly in clinical settings for molecular diagnostics, guiding therapy decisions. Somatic variant callers perform paired comparisons of sequencing data from cancer tissue and matched normal tissue in order to detect somatic mutations. The advent of many new somatic variant callers creates a need for comparison and validation of the tools, as no de facto standard for detection of somatic mutations exists and only limited comparisons have been reported. We have performed a comprehensive evaluation using exome sequencing and targeted deep sequencing data of paired tumor-normal samples from five breast cancer patients to evaluate the performance of nine publicly available somatic variant callers: EBCall, Mutect, Seurat, Shimmer, Indelocator, Somatic Sniper, Strelka, VarScan 2 and Virmid for the detection of single nucleotide mutations and small deletions and insertions. We report a large variation in the number of calls from the nine somatic variant callers on the same sequencing data and highly variable agreement. Sequencing depth had markedly diverse impact on individual callers, as for some callers, increased sequencing depth highly improved sensitivity. For SNV calling, we report EBCall, Mutect, Virmid and Strelka to be the most reliable somatic variant callers for both exome sequencing and targeted deep sequencing. For indel calling, EBCall is superior due to high sensitivity and robustness to changes in sequencing depths.
MeSH Terms
Algorithms
Base Sequence
Carcinoma, Ductal, Breast/genetics
Computational Biology/methods
Exome/genetics
Female
High-Throughput Nucleotide Sequencing/methods
Humans
Mutation/genetics
Sequence Analysis, DNA/methods
Authors & Affiliations
5 authors, click to expand affiliations / ORCID
Krøigård Anne Bruun
Department of Clinical Genetics, Odense University Hospital, Sdr. Boulevard 29, 5000, Odense, Denmark. | Human Genetics, Institute of Clinical Research, University of Southern Denmark, Winsløvparken 19, 5000, Odense, Denmark.
Thomassen Mads
Department of Clinical Genetics, Odense University Hospital, Sdr. Boulevard 29, 5000, Odense, Denmark. | Human Genetics, Institute of Clinical Research, University of Southern Denmark, Winsløvparken 19, 5000, Odense, Denmark.
Lænkholm Anne-Vibeke
Department of Pathology, Slagelse Hospital, Ingemannsvej 18, 4200, Slagelse, Denmark.
Kruse Torben A
Department of Clinical Genetics, Odense University Hospital, Sdr. Boulevard 29, 5000, Odense, Denmark. | Human Genetics, Institute of Clinical Research, University of Southern Denmark, Winsløvparken 19, 5000, Odense, Denmark.
Larsen Martin Jakob
Department of Clinical Genetics, Odense University Hospital, Sdr. Boulevard 29, 5000, Odense, Denmark. | Human Genetics, Institute of Clinical Research, University of Southern Denmark, Winsløvparken 19, 5000, Odense, Denmark.
References (20)
20 references, click to expand
-
A comparative analysis of algorithms for somatic SNV detection in cancer.
Bioinformatics. 2013 Sep 15;29(18):2223-30
PMID: 23842810
-
Combining calls from multiple somatic mutation-callers.
BMC Bioinformatics. 2014;15:154
PMID: 24885750
-
An empirical Bayesian framework for somatic mutation detection from cancer genome sequencing data.
Nucleic Acids Res. 2013 Apr;41(7):e89
PMID: 23471004
-
Cake: a bioinformatics pipeline for the integrated analysis of somatic variants in cancer genomes.
Bioinformatics. 2013 Sep 1;29(17):2208-10
PMID: 23803469
-
VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing.
Genome Res. 2012 Mar;22(3):568-76
PMID: 22300766
-
Detecting somatic point mutations in cancer genome sequencing data: a comparison of mutation callers.
Genome Med. 2013 Oct 11;5(10):91
PMID: 24112718
-
Genotype and SNP calling from next-generation sequencing data.
Nat Rev Genet. 2011 Jun;12(6):443-51
PMID: 21587300
-
Sequence-specific error profile of Illumina sequencers.
Nucleic Acids Res. 2011 Jul;39(13):e90
PMID: 21576222
-
Accurately identifying low-allelic fraction variants in single samples with next-generation sequencing: applications in tumor subclone resolution.
Hum Mutat. 2013 Oct;34(10):1432-8
PMID: 23766071
-
Identification of somatic mutations in cancer through Bayesian-based analysis of sequenced genome pairs.
BMC Genomics. 2013;14:302
PMID: 23642077
-
Optimizing cancer genome sequencing and analysis.
Cell Syst. 2015 Sep 23;1(3):210-223
PMID: 26645048
-
A framework for variation discovery and genotyping using next-generation DNA sequencing data.
Nat Genet. 2011 May;43(5):491-8
PMID: 21478889
-
SomaticSniper: identification of somatic point mutations in whole genome sequencing data.
Bioinformatics. 2012 Feb 1;28(3):311-7
PMID: 22155872
-
Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs.
Bioinformatics. 2012 Jul 15;28(14):1811-7
PMID: 22581179
-
Clonal expansion and linear genome evolution through breast cancer progression from pre-invasive stages to asynchronous metastasis.
Oncotarget. 2015 Mar 20;6(8):5634-49
PMID: 25730902
-
A survey of tools for variant analysis of next-generation genome sequencing data.
Brief Bioinform. 2014 Mar;15(2):256-78
PMID: 23341494
-
Comparison of somatic mutation calling methods in amplicon and whole exome sequence data.
BMC Genomics. 2014;15:244
PMID: 24678773
-
Shimmer: detection of genetic alterations in tumors using next-generation sequence data.
Bioinformatics. 2013 Jun 15;29(12):1498-503
PMID: 23620360
-
Virmid: accurate detection of somatic mutations with sample impurity inference.
Genome Biol. 2013;14(8):R90
PMID: 23987214
-
Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples.
Nat Biotechnol. 2013 Mar;31(3):213-9
PMID: 23396013