Impact of post-alignment processing in variant discovery from whole exome data

Abstract

BackgroundGATK Best Practices workflows are widely used in large-scale sequencing projects and recommend post-alignment processing before variant calling. Two key post-processing steps include the computationally intensive local realignment around known INDELs and base quality score recalibration (BQSR). Both have been shown to reduce erroneous calls; however, the findings are mainly supported by the analytical pipeline that incorporates BWA and GATK UnifiedGenotyper. It is not known whether there is any benefit of post-processing and to what extent the benefit might be for pipelines implementing other methods, especially given that both mappers and callers are typically updated. Moreover, because sequencing platforms are upgraded regularly and the new platforms provide better estimations of read quality scores, the need for post-processing is also unknown. Finally, some regions in the human genome show high sequence divergence from the reference genome; it is unclear whether there is benefit from post-processing in these regions.ResultsWe used both simulated and NA12878 exome data to comprehensively assess the impact of post-processing for five or six popular mappers together with five callers. Focusing on chromosome 6p21.3, which is a region of high sequence divergence harboring the human leukocyte antigen (HLA) system, we found that local realignment had little or no impact on SNP calling, but increased sensitivity was observed in INDEL calling for the Stampy + GATK UnifiedGenotyper pipeline. No or only a modest effect of local realignment was detected on the three haplotype-based callers and no evidence of effect on Novoalign. BQSR had virtually negligible effect on INDEL calling and generally reduced sensitivity for SNP calling that depended on caller, coverage and level of divergence. Specifically, for SAMtools and FreeBayes calling in the regions with low divergence, BQSR reduced the SNP calling sensitivity but improved the precision when the coverage is insufficient. However, in regions of high divergence (e.g., the HLA region), BQSR reduced the sensitivity of both callers with little gain in precision rate. For the other three callers, BQSR reduced the sensitivity without increasing the precision rate regardless of coverage and divergence level.ConclusionsWe demonstrated that the gain from post-processing is not universal; rather, it depends on mapper and caller combination, and the benefit is influenced further by sequencing depth and divergence level. Our analysis highlights the importance of considering these key factors in deciding to apply the computationally intensive post-processing to Illumina exome data.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1279-z) contains supplementary material, which is available to authorized users.

Genetic variants analysis of three dromedary camels using whole genome sequencing data

“…Picard tools ( http://broadinstitute.github.io/picard ) were then used to remove the duplicated reads. Duplicated reads were considered to be identical reads that arise during PCR and are mapped to same genomic position during mapping to the reference [ 19 ].…”

Machine learning approaches reveal genomic regions associated with sugarcane brown rust resistance

“…Indels, however, were not selected. These variants identified by in silico strategies do not provide reliable information, showing elevated divergence between the existent callers and a probability of producing spurious variants 114 .…”

Utility of pooled sequencing for association mapping in nonmodel organisms

“…Indels can cause misalignments that potentially lead to false‐positive proximate SNPs and can be removed (Kofler, Pandey, & Schlötterer, ). Other methods, such as local realignment around indels (e.g., GATK IndelRealigner , McKenna et al., ), can be used to retain SNPs around indel regions but are more computationally intensive and may not be necessary due to the small effect of indel regions on proper SNP calling (Tian, Yan, Kalmbach, & Slager, ). At the end of the alignment phase, a filtered sync file, VCF file, allele frequency table and mpileup file containing all pools of quality‐filtered aligned data are produced.…”