Take a sneak peek at the new NIST.gov and let us know what you think!
(Please note: some content may not be complete on the beta site.).

View the beta site
NIST logo

Publication Citation: Synthetic spike-in standards improve run-specific systematic error analysis for DNA and RNA sequencing

NIST Authors in Bold

Author(s): Justin M. Zook; Daniel V. Samarov; Jennifer H. McDaniel; Shurjo Sen; Marc L. Salit;
Title: Synthetic spike-in standards improve run-specific systematic error analysis for DNA and RNA sequencing
Published: July 31, 2012
Abstract: While the importance of random sequencing errors decreases at higher DNA or RNA sequencing depths, systematic sequencing errors (SSEs) dominate at high sequencing depths and can be difficult to distinguish from biological variants. These SSEs can cause base quality scores to underestimate the probability of error at certain genomic positions, resulting in false positive variant calls, particularly in mixtures such as samples with RNA editing, tumors, circulating tumor cells, bacteria, mitochondrial heteroplasmy, or pooled DNA. Most algorithms proposed for correction of SSEs require a training data set, which is typically either from a part of the data set being recalibrated (Genome Analysis ToolKit, or GATK) or from a separate data set with special characteristics (SysCall). Here, we combine the advantages of these approaches by adding synthetic RNA spike-in standards to human RNA, and use GATK to recalibrate base quality scores with reads mapped to the spike-in standards. Compared to conventional GATK recalibration that uses reads mapped to the genome, spike-ins improve the accuracy of Illumina base quality scores by a mean of 5 units, and by as much as 13 units at CpG sites. In addition, since reads mapping to the genome are not used for recalibration, our method allows run-specific recalibration even for the many species without a comprehensive and accurate SNP database. We also use GATK with the spike-in standards to demonstrate that the Illumina RNA sequencing runs overestimate quality scores for AC, CC, GC, GG, and TC dinucleotides, while SOLiD has less dinucleotide bias but more bias for certain cycles. We conclude that using these DNA and RNA spike-in standards with GATK improves base quality score recalibration.
Citation: PLoS One
Volume: 7
Issue: 7
Keywords: bioscience and health, DNA and RNA sequencing, systematic sequencing errors
Research Areas: Molecular Pathology, Clinical Diagnostics, Standard Reference Materials, Bioscience & Health
DOI: http://dx.doi.org/10.1371/journal.pone.0041356  (Note: May link to a non-U.S. Government webpage)
PDF version: PDF Document Click here to retrieve PDF version of paper (920KB)