This is a test version of Biostars. For the public version, visit https://www.biostars.org.
somatic variant calling without matched normal in long-reads

When calling somatic variants in ONT cancer data (from tumours), is there grounds for using matched short-read (illumina) WGS from normal (blood) to remove germline variants if matched normal sequencing from long-reads isn't available?

Leading on from this, would it be sensible, given a wider cohort of short-read germline variant calls from many individuals, to generate a larger panel of normals to subtract from long-read tumour sequencing variant calls?

The primary reason being the expense of sequencing normal tissue if matched SR-WGS is already available...

Most studies have matched tumour-normal pairs, but I have found an example where a panel of normals is constructed from a mixture of long- and short-read sequencing data.

And another where SVs from multiple sequencing technologies were filtered against 15 healthy genomes sequenced with pacbio.

Some tools such as nanomonsv have a panel of normal function included, which makes use of 30 healthy ONT normal samples from the human pangenome reference consortium.

Then there is the question of population databases such as dbSNP and gnomAD, which are also not based on long-read data (to my knowledge)

variants cancer illumina ont nanopore

I like this question. Thank you for digging into this!

See also https://github.com/KolmogorovLab/Severus for another tool using tumor-normal pairs.

I guess you could use a PoN based on short reads, but that will inherently be incomplete so you will always miss things that are only detected with long reads.

What about using a healthy genome sequenced with long-reads e.g. HG002 ?

1 answer

We are working on such a panel-of-normal SV set for 1000 Genomes samples: ONT SV analysis. We hope to release delly's, sniffles's and SVarp's long-read SV calls for this in April and then it should be straight forward to filter somatic SVs against this.

Great! I know the latest gnomad release (v4.0) has been benchmarked with long-reads but actual long-read calls is even better!

Severus have now included this panel in their somatic variant caller for tumour-only samples

Log in to answer this question.