Question

DeepVariant and Gene Panels

1

Entering edit mode

4.8 years ago

vaktson ▴ 10

Dear *,

I currently try out DeepVariant, which was designed for Whole Exome Sequencing and Whole Genome Sequencing.

Now I wonder if there is any problem, if I would use it for gene panels, so a subset of whole exome sequencing (more or less).

Best regards, Daniel

deepvariant genepanel panel ngs • 1.9k views

ADD COMMENT • link updated 4.8 years ago by Charles Warden 8.2k • written 4.8 years ago by vaktson ▴ 10

score 2 · Answer 1 · 2019-07-03

In general, I would probably try GATK before DeepVariant:

C: How Much Does It Cost To Align A Flowcell In The Cloud?

GATK performance also seems similar (or slightly better) than DeepVariant in this paper (with BWA alignments):

https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2928-9

enter image description here

In the specific case of targeted gene panels, I have previously found indications that there can be a relatively greater need for modification for targeted gene panels than Exome sequencing (where I was using VarScan with a set of non-default parameters for targeted gene panels):

https://peerj.com/articles/600/

targeted gene versus Exome variant calling

[please note that the y-axis is different in the panels from the figure in the Warden et al. 2014 paper]

So, without directly testing DeepVariant, I would definitely recommend checking for your variants (and I think it is possible you may want to use an alternative model, from a large panel of samples with your own target design).

score 1 · Answer 2 · 2019-07-03

I don't think there would be a problem using DeepVariant for targeted gene panels. make_examples accepts a bed file defining your regions with the --regions flag; you would have to decide whether to use the WGS model or the WES model for call_variants. I'd recommend using the WES model, assuming your targeted data has the kind of non-uniform coverage patterns seen in WES data.