FVC As an Adaptive and Accurate Method for Filtering Variants from Popular NGS Analysis Pipelines

Overview

Journal Commun Biol

Specialty Biology

Date 2022 Sep 16

PMID 36114280

Authors

Yongyong Ren

Yan Kong

XiaoCheng Zhou

Georgi Z Genchev

Chao Zhou

Hongyu Zhao

Hui Lu

Affiliations

Soon will be listed here.

Abstract

The quality control of variants from whole-genome sequencing data is vital in clinical diagnosis and human genetics research. However, current filtering methods (Frequency, Hard-Filter, VQSR, GARFIELD, and VEF) were developed to be utilized on particular variant callers and have certain limitations. Especially, the number of eliminated true variants far exceeds the number of removed false variants using these methods. Here, we present an adaptive method for quality control on genetic variants from different analysis pipelines, and validate it on the variants generated from four popular variant callers (GATK HaplotypeCaller, Mutect2, Varscan2, and DeepVariant). FVC consistently exhibited the best performance. It removed far more false variants than the current state-of-the-art filtering methods and recalled ~51-99% true variants filtered out by the other methods. Once trained, FVC can be conveniently integrated into a user-specific variant calling pipeline.

References

Wei Q, Dunbrack Jr R . The role of balanced training and testing data sets for binary classifiers in bioinformatics. PLoS One. 2013; 8(7):e67863. PMC: 3706434. DOI: 10.1371/journal.pone.0067863. View

Ochoa D, Jarnuczak A, Vieitez C, Gehre M, Soucheray M, Mateus A . The functional landscape of the human phosphoproteome. Nat Biotechnol. 2019; 38(3):365-373. PMC: 7100915. DOI: 10.1038/s41587-019-0344-3. View

Poplin R, Chang P, Alexander D, Schwartz S, Colthurst T, Ku A . A universal SNP and small-indel variant caller using deep neural networks. Nat Biotechnol. 2018; 36(10):983-987. DOI: 10.1038/nbt.4235. View

Bentley D, Balasubramanian S, Swerdlow H, Smith G, Milton J, Brown C . Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008; 456(7218):53-9. PMC: 2581791. DOI: 10.1038/nature07517. View

Neums L, Suenaga S, Beyerlein P, Anders S, Koestler D, Mariani A . VaDiR: an integrated approach to Variant Detection in RNA. Gigascience. 2017; 7(2). PMC: 5827345. DOI: 10.1093/gigascience/gix122. View

Van der Auwera G, Carneiro M, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A . From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr Protoc Bioinformatics. 2014; 43:11.10.1-11.10.33. PMC: 4243306. DOI: 10.1002/0471250953.bi1110s43. View

Zhao E, Jones M, Jones S . Whole-Genome Sequencing in Cancer. Cold Spring Harb Perspect Med. 2018; 9(3). PMC: 6396343. DOI: 10.1101/cshperspect.a034579. View

Piskol R, Ramaswami G, Li J . Reliable identification of genomic variants from RNA-seq data. Am J Hum Genet. 2013; 93(4):641-51. PMC: 3791257. DOI: 10.1016/j.ajhg.2013.08.008. View

Jiang J, Huang J, Gu J, Cai X, Zhao H, Lu H . Genomic analysis of a spinal muscular atrophy (SMA) discordant family identifies a novel mutation in TLL2, an activator of growth differentiation factor 8 (myostatin): a case report. BMC Med Genet. 2020; 20(1):204. PMC: 6938020. DOI: 10.1186/s12881-019-0935-3. View

10.

Coudray A, Battenhouse A, Bucher P, Iyer V . Detection and benchmarking of somatic mutations in cancer genomes using RNA-seq data. PeerJ. 2018; 6:e5362. PMC: 6074801. DOI: 10.7717/peerj.5362. View

11.

Adelson R, Renton A, Li W, Barzilai N, Atzmon G, Goate A . Empirical design of a variant quality control pipeline for whole genome sequencing data using replicate discordance. Sci Rep. 2019; 9(1):16156. PMC: 6834861. DOI: 10.1038/s41598-019-52614-7. View

12.

Newell F, Kong Y, Wilmott J, Johansson P, Ferguson P, Cui C . Whole-genome landscape of mucosal melanoma reveals diverse drivers and therapeutic targets. Nat Commun. 2019; 10(1):3163. PMC: 6639323. DOI: 10.1038/s41467-019-11107-x. View

13.

Stranneheim H, Lagerstedt-Robinson K, Magnusson M, Kvarnung M, Nilsson D, Lesko N . Integration of whole genome sequencing into a healthcare setting: high diagnostic rates across multiple clinical entities in 3219 rare disease patients. Genome Med. 2021; 13(1):40. PMC: 7968334. DOI: 10.1186/s13073-021-00855-5. View

14.

McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A . The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20(9):1297-303. PMC: 2928508. DOI: 10.1101/gr.107524.110. View

15.

Ravasio V, Ritelli M, Legati A, Giacopuzzi E . GARFIELD-NGS: Genomic vARiants FIltering by dEep Learning moDels in NGS. Bioinformatics. 2018; 34(17):3038-3040. DOI: 10.1093/bioinformatics/bty303. View

16.

Werling D, Pochareddy S, Choi J, An J, Sheppard B, Peng M . Whole-Genome and RNA Sequencing Reveal Variation and Transcriptomic Coordination in the Developing Human Prefrontal Cortex. Cell Rep. 2020; 31(1):107489. PMC: 7295160. DOI: 10.1016/j.celrep.2020.03.053. View

17.

Mak A, White M, Eckalbar W, Szpiech Z, Oh S, Pino-Yanes M . Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma. Am J Respir Crit Care Med. 2018; 197(12):1552-1564. PMC: 6006403. DOI: 10.1164/rccm.201712-2529OC. View

18.

Zook J, Chapman B, Wang J, Mittelman D, Hofmann O, Hide W . Integrating human sequence data sets provides a resource of benchmark SNP and indel genotype calls. Nat Biotechnol. 2014; 32(3):246-51. DOI: 10.1038/nbt.2835. View

19.

Kendig K, Baheti S, Bockol M, Drucker T, Hart S, Heldenbrand J . Sentieon DNASeq Variant Calling Workflow Demonstrates Strong Computational Performance and Accuracy. Front Genet. 2019; 10:736. PMC: 6710408. DOI: 10.3389/fgene.2019.00736. View

20.

Chicco D, Jurman G . The advantages of the Matthews correlation coefficient (MCC) over F1 score and accuracy in binary classification evaluation. BMC Genomics. 2020; 21(1):6. PMC: 6941312. DOI: 10.1186/s12864-019-6413-7. View