Question: SVtyper on lumpyexpress VCF: somatic-germline
gravatar for ATpoint
3.2 years ago by
ATpoint42k wrote:

Does someone has experience in filtering somatic variants from SVtyper-genotyped VCF files produced with lumpy (WGS, 30-50x)? Lumpyexpress was run as instructed in the manual, and the VCF was genotyped with SVtyper, everything on default settings. There are some issues on Git on how to discriminate somatic from germline variants, but still I did not really make progress in how to filter out the somatics. Anyone experienced in this?

EDIT: Ryan Layer was so kind to give this response on Github:

Select the variants that are non reference in your tumor and have no evidence in the normal. You can use SnpSift to do this with something like: GEN[0].GT != 0/0 && GEN[1].AO == 0 The syntax is not exactly right. Check here for the exact details.

somatic variants lumpy svtyper sv • 2.0k views
ADD COMMENTlink modified 3.1 years ago • written 3.2 years ago by ATpoint42k

Which with proper SnpSift syntax would be the following, given that the tumor column comes before the normal column:

SnpSift.jar filter --file in.vcf "( isVariant( GEN[0] ) && ( GEN[1].AO == 0 ) )" > out.vcf
ADD REPLYlink written 3.1 years ago by ATpoint42k

In the absence of matched 'normal' DNA, such as that from leukocytes in the plasma buffy coat or buccal swab DNA, Why not just build your own 'in house' database of normal DNA by downloading all 1000 Genomes Phase III FASTQ files, processing them, and then creating an easy lookup in order to filter out all likely germline variants?

From what I've seen so far, the major cancer centers (in the USA) each has their own 'panel of normals', which they use for filtering.

ADD REPLYlink written 3.1 years ago by Kevin Blighe68k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1178 users visited in the last hour