tarSVM: Improving the accuracy of variant calls derived from microfluidic PCR-based targeted next generation sequencing using a support vector machine

BMC Bioinformatics

Methodology Article

Hyun Min Kang¹  Edgar A. Otto²  Simone Sanna-Cherchi³  Ali Gharavi³  Rajendra Bhimma⁴  Brendan Crawford⁵  Christopher E. Gillies⁵  Virginia Vega-Warner⁵  Catherine C. Robertson⁵  Matthew G. Sampson⁶  Cheryl Winkler⁷ 

[1] Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, MI, USA;Department of Internal Medicine-Nephrology, University of Michigan School of Medicine, Ann Arbor, MI, USA;Department of Medicine, Columbia University College of Physicians and Surgeons, New York, NY, USA;Department of Paediatrics and Child Health, University of KwaZulu Natal, Durban, South Africa;Department of Pediatrics-Nephrology, University of Michigan School of Medicine, Ann Arbor, MI, USA;Department of Pediatrics-Nephrology, University of Michigan School of Medicine, Ann Arbor, MI, USA;3560B MSRB2, 1150 West Medical Center Drive, 48109, Ann Arbor, MI, USA;NCI, Frederick National Lab for Cancer Research, Molecular Genetics Epidemiology Section, Frederick, MD, USA;

关键词: Microfluidic;    PCR;    Next-generation sequencing;    Variant calling;    Accuracy;    Support vector machine;    Nephrotic;   

DOI  :  10.1186/s12859-016-1108-4

 received in 2015-12-30, accepted in 2016-06-02,  发布年份 2016

来源: Springer

【 摘 要 】

BackgroundTargeted sequencing of discrete gene sets is a cost effective strategy to screen subjects for monogenic forms of disease. One method to achieve this pairs microfluidic PCR with next generation sequencing. The PCR step of this pipeline creates challenges in accurate variant calling. This includes that most reads targeting a specific exon are duplicates that have been amplified from the PCR step. To reduce false positive variant calls from these experiments, previous studies have used threshold-based filtering of alternative allele depth ratio and manual inspection of the alignments. However even after manual inspection and filtering, many variants fail to be validated via Sanger sequencing. To improve the accuracy of variant calling from these experiments, we are challenged to design a variant filtering strategy that sufficiently models microfluidic PCR-specific issues.ResultsWe developed an open source variant filtering pipeline, targeted sequencing support vector machine (“tarSVM”), that uses a Support Vector Machine (SVM) and a new score the normalized allele dosage test to identify high quality variants from microfluidic PCR data. tarSVM maximizes training knowledge by selecting variants that are likely true and likely false variants by incorporating knowledge from the 1000 Genomes and the Exome Aggregation Consortium projects. tarSVM improves on previous approaches by synthesizing variant features from the Genome Analysis Toolkit and allele dosage information. We compared the accuracy of tarSVM versus existing variant quality filtering strategies on two cohorts (n = 474 and n = 1152), and validated our method on a third cohort (n = 75). In the first cohort, our method achieved 84.5 % accuracy of predicting whether or not a variant would be validated with Sanger sequencing versus 78.8 % for the second most accurate method. In the second cohort, our method had an accuracy of 73.3 %, versus 61.5 % for the second best method. Finally, our method had a false discovery rate of 5 % for the validation cohort.ConclusionstarSVM increases the accuracy of variant calling when using microfluidic PCR based targeted sequencing approaches. This results in higher confidence downstream analyses, and ultimately reduces the costs Sanger validation. Our approach is less labor intensive than existing approaches, and is available as an open source pipeline for read trimming, aligning, variant calling, and variant quality filtering on GitHub at https://github.com/christopher-gillies/TargetSpecificGATKSequencingPipeline.

【 授权许可】

CC BY   
© Gillies et al. 2016