LRSomatic: a highly scalable and robust pipeline for somatic variant calling in long-read sequencing data

Abstract

MotivationLong-read sequencing is increasingly used in cancer research and clinical genomics due to its ability to resolve complex genomic variation and previously inaccessible regions of the genome. However, dedidated workflows for comprehensive somatic variant analysis from long-read whole-genome data remain scarce, limiting uptake in cancer genomics. ResultsWe present LRSomatic, a Nextflow-based, nf-core-compliant pipeline supporting somatic SNV, indel, structural variant, and copy number calling from PacBio HiFi and ONT data. LRSomatic supports paired tumor-normal and tumor-only designs, as well as integration of epigenetic integration via Fiber-seq. Benchmarked on COLO829 and HG008 reference cell lines, LRSomatic achieves state-of-the-art performance across both platforms and variant types. Applied to a case of clear cell sarcoma, it recovers all identified driver alterations, including the pathognomonic EWSR1::ATF1 fusion, and resolves haplotype-specific chromatin accessibility via Fiber-seq. Availability and ImplementationFreely available at https://github.com/intgenomicslab/lrsomatic, implemented in Nextflow DSL2, supported via Docker and Singularity.