Abstract:Abstract Background Accurate variant detection in the coding regions of the human genome is a key requirement for molecular diagnostics of Mendelian disorders. Efficiency of variant discovery from next-generation sequencing (NGS) data depends on multiple factors, including reproducible coverage biases of NGS methods and the performance of read alignment and variant calling software. Although variant caller benchmarks are published constantly, no previous publications have leveraged the full extent of available gold standard whole-genome (WGS) and whole-exome (WES) sequencing datasets. Results In this work, we systematically evaluated the performance of 4 popular short read aligners (Bowtie2, BWA, Isaac, and Novoalign) and 9 novel and well-established variant calling and filtering methods (Clair3, DeepVariant, Octopus, GATK, FreeBayes, and Strelka2) using a set of 14 “gold standard” WES and WGS datasets available from Genome In A Bottle (GIAB) consortium. Additionally, we have indirectly evaluated each pipeline’s performance using a set of 6 non-GIAB samples of African and Russian ethnicity. In our benchmark, Bowtie2 performed significantly worse than other aligners, suggesting it should not be used for medical variant calling. When other aligners were considered, the accuracy of variant discovery mostly depended on the variant caller and not the read aligner. Among the tested variant callers, DeepVariant consistently showed the best performance and the highest robustness. Other actively developed tools, such as Clair3, Octopus, and Strelka2, also performed well, although their efficiency had greater dependence on the quality and type of the input data. We have also compared the consistency of variant calls in GIAB and non-GIAB samples. With few important caveats, best-performing tools have shown little evidence of overfitting. Conclusions The results show surprisingly large differences in the performance of cutting-edge tools even in high confidence regions of the coding genome. This highlights the importance of regular benchmarking of quickly evolving tools and pipelines. We also discuss the need for a more diverse set of gold standard genomes that would include samples of African, Hispanic, or mixed ancestry. Additionally, there is also a need for better variant caller assessment in the repetitive regions of the coding genome.

Variant calling in polyploids for population and quantitative genetics

Genotyping Polyploids from Messy Sequencing Data

Single Nucleotide Polymorphism Identification in Polyploids: A Review, Example, and Recommendations

Benchmarking small-variant genotyping in polyploids

Tools for Genetic Studies in Experimental Populations of Polyploids

Impact of genotype‐calling methodologies on genome‐wide association and genomic prediction in polyploids

Multiple Variant Calling Pipelines in Wheat Whole Exome Sequencing

Sequence coverage required for accurate genotyping by sequencing in polyploid species

Genotyping single nucleotide polymorphisms and inferring ploidy by amplicon sequencing for polyploid, ploidy‐variable organisms

Towards Better Understanding of Artifacts in Variant Calling from High-Coverage Samples

The Analysis of Polyploid Genetic Data

Systematic benchmark of state-of-the-art variant calling pipelines identifies major factors affecting accuracy of coding sequence variant discovery

Benchmarking variant callers in next-generation and third-generation sequencing analysis

Genotype Calling from Population-Genomic Sequencing Data

Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array

A Benchmark of Genetic Variant Calling Pipelines Using Metagenomic Short-Read Sequencing

Large-scale Genotyping of Complex DNA

Scalable approaches for functional analyses of whole-genome sequencing non-coding variants

Detection of Genomic Variations and DNA Polymorphisms and Impact on Analysis of Meiotic Recombination and Genetic Mapping.

polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids