» Articles » PMID: 39345270

NCBench: Providing an Open, Reproducible, Transparent, Adaptable, and Continuous Benchmark Approach for DNA-sequencing-based Variant Calling

Abstract

We present the results of the human genomic small variant calling benchmarking initiative of the German Research Foundation (DFG) funded Next Generation Sequencing Competence Network (NGS-CN) and the German Human Genome-Phenome Archive (GHGA). In this effort, we developed NCBench, a continuous benchmarking platform for the evaluation of small genomic variant callsets in terms of recall, precision, and false positive/negative error patterns. NCBench is implemented as a continuously re-evaluated open-source repository. We show that it is possible to entirely rely on public free infrastructure (Github, Github Actions, Zenodo) in combination with established open-source tools. NCBench is agnostic of the used dataset and can evaluate an arbitrary number of given callsets, while reporting the results in a visual and interactive way. We used NCBench to evaluate over 40 callsets generated by various variant calling pipelines available in the participating groups that were run on three exome datasets from different enrichment kits and at different coverages. While all pipelines achieve high overall quality, subtle systematic differences between callers and datasets exist and are made apparent by NCBench.These insights are useful to improve existing pipelines and develop new workflows. NCBench is meant to be open for the contribution of any given callset. Most importantly, for authors, it will enable the omission of repeated re-implementation of paper-specific variant calling benchmarks for the publication of new tools or pipelines, while readers will benefit from being able to (continuously) observe the performance of tools and pipelines at the time of reading instead of at the time of writing.

References
1.
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

2.
Barbitoff Y, Abasov R, Tvorogova V, Glotov A, Predeus A . Systematic benchmark of state-of-the-art variant calling pipelines identifies major factors affecting accuracy of coding sequence variant discovery. BMC Genomics. 2022; 23(1):155. PMC: 8862519. DOI: 10.1186/s12864-022-08365-3. View

3.
Pedersen B, Quinlan A . Mosdepth: quick coverage calculation for genomes and exomes. Bioinformatics. 2017; 34(5):867-868. PMC: 6030888. DOI: 10.1093/bioinformatics/btx699. View

4.
Krusche P, Trigg L, Boutros P, Mason C, De La Vega F, Moore B . Best practices for benchmarking germline small-variant calls in human genomes. Nat Biotechnol. 2019; 37(5):555-560. PMC: 6699627. DOI: 10.1038/s41587-019-0054-x. View

5.
Zhao S, Agafonov O, Azab A, Stokowy T, Hovig E . Accuracy and efficiency of germline variant calling pipelines for human genome data. Sci Rep. 2020; 10(1):20222. PMC: 7678823. DOI: 10.1038/s41598-020-77218-4. View