Quantifying single nucleotide variant detection sensitivity in exome sequencing

Abstract

BackgroundThe targeted capture and sequencing of genomic regions has rapidly demonstrated its utility in genetic studies. Inherent in this technology is considerable heterogeneity of target coverage and this is expected to systematically impact our sensitivity to detect genuine polymorphisms. To fully interpret the polymorphisms identified in a genetic study it is often essential to both detect polymorphisms and to understand where and with what probability real polymorphisms may have been missed.ResultsUsing down-sampling of 30 deeply sequenced exomes and a set of gold-standard single nucleotide variant (SNV) genotype calls for each sample, we developed an empirical model relating the read depth at a polymorphic site to the probability of calling the correct genotype at that site. We find that measured sensitivity in SNV detection is substantially worse than that predicted from the naive expectation of sampling from a binomial. This calibrated model allows us to produce single nucleotide resolution SNV sensitivity estimates which can be merged to give summary sensitivity measures for any arbitrary partition of the target sequences (nucleotide, exon, gene, pathway, exome). These metrics are directly comparable between platforms and can be combined between samples to give “power estimates” for an entire study. We estimate a local read depth of 13X is required to detect the alleles and genotype of a heterozygous SNV 95% of the time, but only 3X for a homozygous SNV. At a mean on-target read depth of 20X, commonly used for rare disease exome sequencing studies, we predict 5–15% of heterozygous and 1–4% of homozygous SNVs in the targeted regions will be missed.ConclusionsNon-reference alleles in the heterozygote state have a high chance of being missed when commonly applied read coverage thresholds are used despite the widely held assumption that there is good polymorphism detection at these coverage levels. Such alleles are likely to be of functional importance in population based studies of rare diseases, somatic mutations in cancer and explaining the “missing heritability” of quantitative traits.

Highly Sensitive Diagnosis of 43 Monogenic Forms of Diabetes or Obesity Through One-Step PCR-Based Enrichment in Combination With Next-Generation Sequencing

“…A limitation of our assay may be the use of a coverage threshold of 83 for variant detection. A recent study has reported that raising the coverage threshold to 133 markedly limited the number of missing variant calls (25). However, in the current study, a concordance rate of 100% was found between WES and PCR-based technology in microfluidic droplets in combination with NGS (in target genes that are present in the panel) using a coverage threshold of 83.…”

Validation of Next Generation Sequencing Cancer Panels for Clinical Somatic Mutation Profiling-- Identification of Source of Variations and Artifacts using FFPE Tissues

“…Whether or not another round of new sequencing technologies comes soon to replace current "next generation" sequencing technologies, their revolutionary impact is here to stay. Recently a huge impulse of NGS publications flooded into the scientific journals [1][2][3][4][5][6][7][8][9][10][11][12]. This is understandable as scientists try to make their mark in the NGS field.…”

Exome Sequencing to Detect Rare Variants Associated With General Cognitive Ability: A Pilot Study

“…The cumulative coverage counts (reported at depths 0-500), together with the genomic target size, were used to obtain the proportion of sites at each depth. The results of the product between estimated values for heterozygous sites from Meynert et al, (2013) and the proportion of sites, at each depth, were summed to get an overall SNP detection sensitivity per sample. The average of detection sensitivities per cohort was obtained by averaging the individual values across samples in each cohort.…”