Comprehensive and accurate genome analysis at scale using DRAGEN accelerated algorithms

Research and medical genomics require comprehensive and scalable solutions to drive the discovery of novel disease targets, evolutionary drivers, and genetic markers with clinical significance. This necessitates a framework to identify all types of variants independent of their size (e.g., SNV/SV) or location (e.g., repeats). Here we present DRAGEN that utilizes novel methods based on multigenomes, hardware acceleration, and machine learning based variant detection to provide novel insights into individual genomes with ~30min computation time (from raw reads to variant detection). DRAGEN outperforms all other state-of-the-art methods in speed and accuracy across all variant types (SNV, indel, STR, SV, CNV) and further incorporates specialized methods to obtain key insights in medically relevant genes (e.g., HLA, SMN, GBA). We showcase DRAGEN across 3,202 genomes and demonstrate its scalability, accuracy, and innovations to further advance the integration of comprehensive genomics for research and medical applications.

PrecisionFDA Truth Challenge V2: Calling variants from short and long reads in difficult-to-map regions

The precisionFDA Truth Challenge V2 aimed to assess the state of the art of variant calling in challenging genomic regions. Starting with FASTQs, 20 challenge participants applied their variant-calling pipelines and submitted 64 variant call sets for one or more sequencing technologies (Illumina, PacBio HiFi, and Oxford Nanopore Technologies). Submissions were evaluated following best practices for benchmarking small variants with updated Genome in a Bottle benchmark sets and genome stratifications. Challenge submissions included numerous innovative methods, with graph-based and machine learning methods scoring best for short-read and long-read datasets, respectively. With machine learning approaches, combining multiple sequencing technologies performed particularly well. Recent developments in sequencing and variant calling have enabled benchmarking variants in challenging genomic regions, paving the way for the identification of previously unknown clinically relevant variants.

Abstract 5463: Accuracy improvements in somatic whole-genome small-variant calling with the DRAGEN platform

Introduction: Next-generation whole-genome sequencing promises to enable dramatic expansion of precision oncology and personalized cancer care. To support this effort, it is critical to develop computational tools that can analyze sequence data accurately and with rapid turn-around-time. One particularly challenging problem is the calling of somatic variants in matched tumor and normal samples. We present the DRAGENTM 3.5 somatic pipeline that performs 110x/40x whole-genome end-to-end analysis in under two hours with accuracy superior to that of all other tools we have benchmarked, including Mutect2/GATK4 and Strelka2. In addition, DRAGEN 3.5 is robust against variations in coverage, sequencing platform, sample preparation chemistry, and tumor purity. It also tolerates tumor-in-normal contamination, thereby making the pipeline applicable to late-stage solid tumors or hematological cancers.

Methods: DRAGEN 3.5 replaces the legacy genotyping model (originally developed as part of MuTect2) with that of Strelka2. Whereas the MuTect2 model performs separate analyses on the tumor and normal samples, the Strelka2 model performs a joint analysis, allowing (1) detection of systematic errors that affect both samples simultaneously and (2) modeling of tumor-in-normal contamination, where the amount of contamination in the normal sample depends on the allele frequency observed in the tumor sample at the locus in question. In addition, we improved the probabilistic model of systematic error by incorporating models of strand bias and mismapping.

Results: We benchmarked DRAGEN against Mutect2 (included in GATK 4.1.2) and Strelka2 (version 2.9.9) on five public synthetic and real datasets with known truth sets. DRAGEN greatly outperformed the other methods on all five datasets, producing 14-67% and 22-91% fewer false SNV calls, and 35-86% and 48-89% fewer false indel calls than Strelka2 and Mutect2 respectively. DRAGEN also exhibits higher tolerance to tumor-in-normal (TiN) contamination than Strelka2 which is already equipped with a model tolerating TiN contamination. The average end-to-end workflow runtime of the DRAGEN somatic pipeline was 77 minutes, 75% and 830% faster than Strelka2 and Mutect2 taking DRAGEN alignments as input.

Conclusion: The DRAGEN pipeline enables reliable whole-genome analysis that can be scaled to large numbers of samples, leading to better tumor characterization and improved interpretation. We anticipate that it will ultimately fuel progress in oncology, cancer research and precision medicine. The DRAGEN 3.5 somatic pipeline can be run either locally on a DRAGEN server or remotely in the cloud via https://basespace.illumina.com.

Citation Format: Konrad Scheffler, Sangtae Kim, Varun Jain, Jeffrey Yuan, Westley Sherman, Taylor O'Connell, Eric Ojard, Lisa Murray, Rami Mehio, Severine Catreux. Accuracy improvements in somatic whole-genome small-variant calling with the DRAGEN platform [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5463.

A 26-hour system of highly sensitive whole genome sequencing for emergency management of genetic diseases

While the cost of whole genome sequencing (WGS) is approaching the realm of routine medical tests, it remains too tardy to help guide the management of many acute medical conditions. Rapid WGS is imperative in light of growing evidence of its utility in acute care, such as in diagnosis of genetic diseases in very ill infants, and genotype-guided choice of chemotherapy at cancer relapse. In such situations, delayed, empiric, or phenotype-based clinical decisions may meet with substantial morbidity or mortality. We previously described a rapid WGS method, STATseq, with a sensitivity of >96 % for nucleotide variants that allowed a provisional diagnosis of a genetic disease in 50 h. Here improvements in sequencing run time, read alignment, and variant calling are described that enable 26-h time to provisional molecular diagnosis with >99.5 % sensitivity and specificity of genotypes. STATseq appears to be an appropriate strategy for acutely ill patients with potentially actionable genetic diseases.