Reducing INDEL calling errors in whole genome and exome sequencing data

Fang, Han, Wu, Yiyang, Narzisi, G., O'Rawe, Jason, Jimenez Barron, Laura, Rosenbaum, J., Ronemus, M., Iossifov, I., Schatz, M. C., Lyon, Gholson J. (2014) Reducing INDEL calling errors in whole genome and exome sequencing data. Genome Med, 6 (10). p. 89.

[thumbnail of Paper]
PDF (Paper)
Lyon And Schatz and Iossifov Genome Medicine 2014.pdf - Published Version

Download (2MB) | Preview
DOI: 10.1186/s13073-014-0089-z


BACKGROUND: INDELs, especially those disrupting protein-coding regions of the genome, have been strongly associated with human diseases. However, there are still many errors with INDEL variant calling, driven by library preparation, sequencing biases, and algorithm artifacts. METHODS: We characterized whole genome sequencing (WGS), whole exome sequencing (WES), and PCR-free sequencing data from the same samples to investigate the sources of INDEL errors. We also developed a classification scheme based on the coverage and composition to rank high and low quality INDEL calls. We performed a large-scale validation experiment on 600 loci, and find high-quality INDELs to have a substantially lower error rate than low-quality INDELs (7% vs. 51%). RESULTS: Simulation and experimental data show that assembly based callers are significantly more sensitive and robust for detecting large INDELs (>5 bp) than alignment based callers, consistent with published data. The concordance of INDEL detection between WGS and WES is low (53%), and WGS data uniquely identifies 10.8-fold more high-quality INDELs. The validation rate for WGS-specific INDELs is also much higher than that for WES-specific INDELs (84% vs. 57%), and WES misses many large INDELs. In addition, the concordance for INDEL detection between standard WGS and PCR-free sequencing is 71%, and standard WGS data uniquely identifies 6.3-fold more low-quality INDELs. Furthermore, accurate detection with Scalpel of heterozygous INDELs requires 1.2-fold higher coverage than that for homozygous INDELs. Lastly, homopolymer A/T INDELs are a major source of low-quality INDEL calls, and they are highly enriched in the WES data. CONCLUSIONS: Overall, we show that accuracy of INDEL detection with WGS is much greater than WES even in the targeted region. We calculated that 60X WGS depth of coverage from the HiSeq platform is needed to recover 95% of INDELs detected by Scalpel. While this is higher than current sequencing practice, the deeper coverage may save total project costs because of the greater accuracy and sensitivity. Finally, we investigate sources of INDEL errors (for example, capture deficiency, PCR amplification, homopolymers) with various data that will serve as a guideline to effectively reduce INDEL errors in genome sequencing.

Item Type: Paper
Subjects: bioinformatics
bioinformatics > genomics and proteomics
bioinformatics > computational biology
Investigative techniques and equipment > whole exome sequencing
Investigative techniques and equipment > assays > whole exome sequencing
Investigative techniques and equipment > assays > whole genome sequencing
CSHL Authors:
Communities: CSHL Cancer Center Program > Cancer Genetics
CSHL labs > Iossifov lab
CSHL labs > Lyon lab
CSHL labs > Schatz lab
CSHL labs > Wigler lab
Stanley Institute for Cognitive Genomics
Depositing User: Andrea Potopsingh
Date: 2014
Date Deposited: 12 Jan 2015 20:01
Last Modified: 06 Nov 2015 20:43
PMCID: PMC4240813
Related URLs:

Actions (login required)

Administrator's edit/view item Administrator's edit/view item
CSHL HomeAbout CSHLResearchEducationNews & FeaturesCampus & Public EventsCareersGiving