Somatic Structural Variant Discovery from Tumor Specific Strings Computed Against a Draft Composite Germline Reference
- 1University of Milano-Bicocca
Abstract
Somatic structural variants (SVs) are among the most clinically relevant genomic alterations in cancer. However, their accurate detection remains challenging due to complex genomic rearrangements, repetitive genomic regions, difficulties in distinguishing somatic variants from germline variation and sequencing artifacts, and limited sequencing support for low-frequency events. Existing mapping-based somatic SV callers are susceptible to alignment ambiguities in repetitive regions and near complex SV breakpoints, while approaches relying on comparisons of independently generated tumor and normal variant callsets may propagate errors from both callsets. Building on the concept of specific strings, previously applied to germline SV discovery, we introduce VESPA-SS, a novel approach for detecting somatic SVs from paired tumor-normal long-read whole-genome sequencing data. VESPA-SS directly extracts somatic SV signals using specific strings computed from tumor reads against a Draft Composite Germline Reference, without requiring initial read mapping or a high-quality germline assembly. These signals are subsequently clustered and refined through consensus sequence generation and alignment-based procedures to identify somatic SVs. Evaluation on real PacBio HiFi datasets demonstrates that VESPA-SS outperforms state-of-the-art somatic SV callers, achieving F1 scores exceeding 0.92 and 0.90 on the well-characterized HG008 and COLO829 cancer cell lines, respectively, compared with the highest F1 scores of 0.86 and 0.85 achieved by the evaluated methods. These results demonstrate the effectiveness of specific string based signal extraction for accurate somatic SV discovery and establish a new computational paradigm for identifying complex genomic alterations in cancer. VESPA-SS is freely available at https://github.com/dcmonti/VESPA.
Related
- Vritra: gene-specific reference construction for species-resolved functional profiling of metagenomic and metatranscriptomic data
- The Genetic Discovery (GeDi) resource: Deep characterization of a multi-ancestry panel of human reference iPSCs for functional genomic studies
- Reference-free isoform discovery links unannotated RNA processing to disease risk
- A Plasmodium falciparum Pangenome Resource to Drive Structural Variant Discovery and to assist Malaria Control
- A scalable Zebrafish Platform for genetic validation and functional discovery of stage-specific breast cancer metastasis regulators
BibTeX
@misc{monti2026somatic,
title = {{Somatic Structural Variant Discovery from Tumor Specific Strings Computed Against a Draft Composite Germline Reference}},
author = {D. C. Monti and S. Ciccolella and D. Cozzi and G. Della Vedova and P. Bonizzoni and L. Denti and F. Hormozdiari},
howpublished = {bioRxiv},
year = {2026},
month = nov,
doi = {10.64898/2026.11.06.764802},
url = {https://www.biorxiv.org/content/10.64898/2026.11.06.764802v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Somatic Structural Variant Discovery from Tumor Specific Strings Computed Against a Draft Composite Germline Reference
AU - D. C. Monti
AU - S. Ciccolella
AU - D. Cozzi
AU - G. Della Vedova
AU - P. Bonizzoni
AU - L. Denti
AU - F. Hormozdiari
PY - 2026
DA - 2026/11/06
PB - bioRxiv
DO - 10.64898/2026.11.06.764802
UR - https://www.biorxiv.org/content/10.64898/2026.11.06.764802v1
AB - Somatic structural variants (SVs) are among the most clinically relevant genomic alterations in cancer. However, their accurate detection remains challenging due to complex genomic rearrangements, repetitive genomic regions, difficulties in distinguishing somatic variants from germline variation and sequencing artifacts, and limited sequencing support for low-frequency events. Existing mapping-based somatic SV callers are susceptible to alignment ambiguities in repetitive regions and near complex SV breakpoints, while approaches relying on comparisons of independently generated tumor and normal variant callsets may propagate errors from both callsets. Building on the concept of specific strings, previously applied to germline SV discovery, we introduce VESPA-SS, a novel approach for detecting somatic SVs from paired tumor-normal long-read whole-genome sequencing data. VESPA-SS directly extracts somatic SV signals using specific strings computed from tumor reads against a Draft Composite Germline Reference, without requiring initial read mapping or a high-quality germline assembly. These signals are subsequently clustered and refined through consensus sequence generation and alignment-based procedures to identify somatic SVs. Evaluation on real PacBio HiFi datasets demonstrates that VESPA-SS outperforms state-of-the-art somatic SV callers, achieving F1 scores exceeding 0.92 and 0.90 on the well-characterized HG008 and COLO829 cancer cell lines, respectively, compared with the highest F1 scores of 0.86 and 0.85 achieved by the evaluated methods. These results demonstrate the effectiveness of specific string based signal extraction for accurate somatic SV discovery and establish a new computational paradigm for identifying complex genomic alterations in cancer. VESPA-SS is freely available at https://github.com/dcmonti/VESPA.
ER -