Abstract:Targeted genomic sequencing (TS) greatly benefits precision oncology by rapidly detecting genetic variations with better accuracy and sensitivity owing to its high sequencing depth. Multiple sequencing platforms and variant calling tools are available for TS, making it excruciating for researchers to choose. Therefore, benchmarking study across different platforms and pipelines available for TS is imperative. In this study, we performed a TS of Reference OncoSpan FFPE (HD832) sample enriched by TSO500 panel using four commercially available sequencers, and analyzed the output 50 datasets using five commonly-used bioinformatics pipelines. We systematically investigated the sequencing quality and variant detection sensitivity, expecting to provide optimal recommendations for future research. Four sequencing platforms returned highly concordant results in terms of base quality (Q20 > 94%), sequencing coverage (>97%) and depth (>2000×). Benchmarking revealed good concordance of variant calling across different platforms and pipelines, among which, FASTASeq 300 platform showed the highest sensitivity (100%) and precision (100%) in high-confidence variants calling when analyzed by SNVer and VarScan 2 algorithms. Furthermore, this sequencer demonstrated the shortest sequencing time (∼21 h) at the sequencing mode PE150. Through the intersection of 50 datasets generated in this study, we recommended a novel set of variant genes outside the truth set published by HD832, expecting to replenish HD832 for future research on tumor variant diagnosis. Besides, we applied these five tools to another panel (TargetSeq One) for Twist cfDNA Pan-cancer Reference Standard, comprehensive consideration of SNP and InDel sensitivity, SNVer and VarScan 2 performed best among them. Furthermore, SNVer and VarScan 2 also performed best for six cancer cell lines samples regarding SNP and InDel sensitivity. Considering the dissimilarity of variant calls across different pipelines for datasets from the same platform, we recommended an integration of multiple tools to improve variant calling sensitivity and accuracy for the cancer genome. Illumina and GeneMind technologies can be used independently or together by public health laboratories performing tumor TS. SNVer and VarScan 2 perform better regarding variant detection sensitivity for three typical tumor samples. Our study provides a standardized target sequencing resource to benchmark new bioinformatics protocols and sequencing platforms.

Fast and accurate DNASeq variant calling workflow composed of LUSH toolkit

A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset

[A consensus on the standardization of the next generation sequencing process for the diagnosis of genetic diseases (3) - Data analysis].

From FastQ Data to High‐Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline

Large-scale Genotyping of Complex DNA

A study on fast calling variants from next-generation sequencing data using decision tree

Metapipeline-DNA: A Comprehensive Germline & Somatic Genomics Nextflow Pipeline

Diagnostic and Clinical Utility of Whole Genome Sequencing in a Cohort of Undiagnosed Chinese Families with Rare Diseases

wgatools: an ultrafast toolkit for manipulating whole genome alignments

DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly

Whole Animal Genome Sequencing: user-friendly, rapid, containerized pipelines for processing, variant discovery, and annotation of short-read whole genome sequencing data

DNAscan: a fast, computationally and memory efficient bioinformatics pipeline for the analysis of DNA next-generation-sequencing data

Systematic comparison of variant calling pipelines of target genome sequencing cross multiple next-generation sequencers

Fast-GBS: a new pipeline for the efficient and highly accurate calling of SNPs from genotyping-by-sequencing data

Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects

Kuura—An automated workflow for analyzing WES and WGS data

Variant calling pipelines for whole exome sequencing in clinical context

FLEXIBLE, PRODUCTION-SCALE, HUMAN WHOLE GENOME SEQUENCING ON A BENCHTOP SEQUENCER

OVarFlow: a resource optimized GATK 4 based Open source Variant calling workFlow

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

Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control