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

    • 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
Published December 7, 2023. https://doi.org/10.1101/gr.278070.123
Download PDF Cite Article Permissions Share
cover of Genome Research Vol 36 Issue 5
Current Issue:

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.

Loading
Loading
Loading
Back to top