Whole-genome long-read sequencing downsampling and its effect on variant-calling precision and recall [RESEARCH]

William T. Harvey1, Peter Ebert2,3,4, Jana Ebler2,4, Peter A. Audano5, Katherine M. Munson1, Kendra Hoekzema1, David Porubsky1, Christine R. Beck5,6, Tobias Marschall2,4, Kiran Garimella7 and Evan E. Eichler1,8 1Department of Genome Sciences, University of Washington School of Medicine, Seattle, Washington 98195-5065, USA; 2Institute for Medical Biometry and Bioinformatics, Medical Faculty, Heinrich Heine University, 40225 Düsseldorf, Germany; 3Core Unit Bioinformatics, Medical Faculty, Heinrich Heine University, 40225 Düsseldorf, Germany; 4Center for Digital Medicine, Heinrich Heine University, 40225 Düsseldorf, Germany; 5The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut 06032, USA; 6Department of Genetics and Genome Sciences, University of Connecticut Health Center, Farmington, Connecticut 06030-6403, USA; 7Data Sciences Platform, Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, USA; 8Howard Hughes Medical Institute, University of Washington, Seattle, Washington 98195, USA Corresponding author: eeegs.washington.edu Abstract

Advances in long-read sequencing (LRS) technologies continue to make whole-genome sequencing more complete, affordable, and accurate. LRS provides significant advantages over short-read sequencing approaches, including phased de novo genome assembly, access to previously excluded genomic regions, and discovery of more complex structural variants (SVs) associated with disease. Limitations remain with respect to cost, scalability, and platform-dependent read accuracy and the tradeoffs between sequence coverage and sensitivity of variant discovery are important experimental considerations for the application of LRS. We compare the genetic variant-calling precision and recall of Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio) HiFi platforms over a range of sequence coverages. For read-based applications, LRS sensitivity begins to plateau around 12-fold coverage with a majority of variants called with reasonable accuracy (F1 score above 0.5), and both platforms perform well for SV detection. Genome assembly increases variant-calling precision and recall of SVs and indels in HiFi data sets with HiFi outperforming ONT in quality as measured by the F1 score of assembly-based variant call sets. While both technologies continue to evolve, our work offers guidance to design cost-effective experimental strategies that do not compromise on discovering novel biology.

Footnotes

[Supplemental material is available for this article.]

Article published online before print. Article, supplemental material, and publication date are at https://www.genome.org/cgi/doi/10.1101/gr.278070.123.

Freely available online through the Genome Research Open Access option.

Received May 4, 2023. Accepted November 3, 2023.

留言 (0)

沒有登入
gif