Towards standards for analysing structural variants from short-read whole genome sequencing data
- 1Menzies Institute for Medical Research, University of Tasmania
Abstract
No widely accepted standard exists for analysing structural variants from short-read whole genome sequencing (WGS) data. We use comprehensive SV truth sets (HG002 T2T-Q100, Platinum Pedigree) and long-read sequencing databases (1KGP-ONT, CoLoRS) to benchmark SV breakpoints and sizes at base pair resolution. Using 30X Illumina WGS data, only five of sixteen individual tools tested reliably output deletions, insertions and duplications that perfectly match the truth sets. Ensemble pipelines (e.g. GATK-SV, MOPline, Parliament2), frequently used in population-scale datasets, reduced perfect matching accuracy by up to 98% compared to their highest-performing constituent tools. This inaccuracy propagates into a short-read database (gnomAD-SV) constructed from an ensemble pipeline (GATK-SV) and likely affects other databases (TOPMed, All of Us, Biobank Japan). As a proof-of-concept, we demonstrate an intra-sample merging strategy that generates the highest F1-scores for deletions and insertions on the HG002 sample, recovering more perfect matches than total matches from the best performing ensemble pipeline tested. Common HG002 truth SVs (allele frequency > 1%) showed extremely high concordance for exact matches between our ensemble callset and two orthogonal long-read sequencing databases (CoLoRS, ONT-1KGP). This study provides a vital step towards robust standards for SV analysis, guiding short-read WGS ensemble pipeline development and improving cross-study comparisons.
Related
- The Genetic Landscape of Structural Variants in Aspergillus flavus using ONT and PacBio Sequencing
- Whole genome sequencing identifies genetic variants associated with senescence-accelerated phenotype of OXYS rats.
- Long-read sequencing reveals genomic landscape of structural variants in a major disease vector
- 3D-genome impacts of TAD-boundary-disrupting ultra-rare structural variants identified in schizophrenia
- Genome-wide associations of structural variants with human traits through imputation from long-read assemblies
BibTeX
@misc{slimmer2026towards,
title = {{Towards standards for analysing structural variants from short-read whole genome sequencing data}},
author = {J. C. Slimmer and J. L. Jones and W. R. Reay and K. P. Burdon and N. B. Blackburn},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.13.758836},
url = {https://www.biorxiv.org/content/10.64898/2026.10.13.758836v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Towards standards for analysing structural variants from short-read whole genome sequencing data
AU - J. C. Slimmer
AU - J. L. Jones
AU - W. R. Reay
AU - K. P. Burdon
AU - N. B. Blackburn
PY - 2026
DA - 2026/10/13
PB - bioRxiv
DO - 10.64898/2026.10.13.758836
UR - https://www.biorxiv.org/content/10.64898/2026.10.13.758836v1
AB - No widely accepted standard exists for analysing structural variants from short-read whole genome sequencing (WGS) data. We use comprehensive SV truth sets (HG002 T2T-Q100, Platinum Pedigree) and long-read sequencing databases (1KGP-ONT, CoLoRS) to benchmark SV breakpoints and sizes at base pair resolution. Using 30X Illumina WGS data, only five of sixteen individual tools tested reliably output deletions, insertions and duplications that perfectly match the truth sets. Ensemble pipelines (e.g. GATK-SV, MOPline, Parliament2), frequently used in population-scale datasets, reduced perfect matching accuracy by up to 98% compared to their highest-performing constituent tools. This inaccuracy propagates into a short-read database (gnomAD-SV) constructed from an ensemble pipeline (GATK-SV) and likely affects other databases (TOPMed, All of Us, Biobank Japan). As a proof-of-concept, we demonstrate an intra-sample merging strategy that generates the highest F1-scores for deletions and insertions on the HG002 sample, recovering more perfect matches than total matches from the best performing ensemble pipeline tested. Common HG002 truth SVs (allele frequency > 1%) showed extremely high concordance for exact matches between our ensemble callset and two orthogonal long-read sequencing databases (CoLoRS, ONT-1KGP). This study provides a vital step towards robust standards for SV analysis, guiding short-read WGS ensemble pipeline development and improving cross-study comparisons.
ER -