Question: InterVar and Varsome discrepancies in ACMG classification
0
gravatar for cocchi.e89
12 months ago by
cocchi.e8950
cocchi.e8950 wrote:

I am working on ACMG classification of variants in a cohort of patients, I found that InterVar classification differs (generally more strict for what does it seems to me at the moment) form Varsome, e.g. var X-107939570-C-T is classified as "Pathogenic" by Varsome (https://varsome.com/variant/hg19/X-107939570-C-T) and Uncertain Significance by InterVar (command line tool: "InterVar: Uncertain significance PVS1=0 PS=[0, 0, 0, 0, 0] PM=[0, 1, 0, 0, 0, 0, 0] PP=[0, 0, 1, 0, 1, 0] BA1=0 BS=[0, 0, 0, 0, 0] BP=[0, 0, 0, 0, 0, 0, 0, 0]")

At the same manner var: X-107898624-G-T. Varsome: Likely Pathogenic (https://varsome.com/variant/hg19/X-107898624-G-T) and Intervar: "InterVar: Uncertain significance PVS1=0 PS=[0, 0, 0, 0, 0] PM=[0, 1, 0, 0, 0, 0, 0] PP=[0, 0, 1, 0, 1, 0] BA1=0 BS=[0, 0, 0, 0, 0] BP=[1, 0, 0, 0, 0, 0, 0, 0]"

Does somebody has got a clue for this difference? There are other command line tools to make this annotation automated (make more similar to Varsome "relaxed" interpretation?

Thanks a lot in advance for any help!

ADD COMMENTlink modified 10 months ago by manuel.belmadani1.2k • written 12 months ago by cocchi.e8950

Hey I am also working on ACMG classification and was wondering what command you used with InterVar and also how your vcf is structured since I am struggling to get an output in general.

Thanks in advance for your help!

ADD REPLYlink written 12 months ago by zankov0

Just clone the InterVar repo on GitHub:

git clone https://github.com/WGLab/InterVar.git

and then access the InterVar directory in which you find the python script InterVar.py Call it on a normal VCF file with, as example (I use python3.7 but any distribution >3.3 should work properly):

python3.7 InterVar.py  \
    -b hg19 \
    -i <your_input.vcf> \
     --input_type=VCF \
     -o <your_output>

And then you will find your output files where you specified. Please note, in the same InterVar directory you must download Annovar executables (annotate_variation.pl, convert2annovar.pl, retrieve_seq_from_fasta.pl, coding_change.pl, table_annovar.pl, variants_reduction.pl) and the first time you run it will take long in order for Annovar to download all its reference databases in the "humandb" directory that will appear in the same "InterVar" directory where you launch the script.

ADD REPLYlink modified 12 months ago • written 12 months ago by cocchi.e8950
0
gravatar for manuel.belmadani
10 months ago by
Canada
manuel.belmadani1.2k wrote:

You can actually replicate what InterVar tells you with the VarSome UI.

Example:

PVS1=0 PS=[0, 0, 0, 0, 0] PM=[0, 1, 0, 0, 0, 0, 0] PP=[0, 0, 1, 0, 1, 0] BA1=0 BS=[0, 0, 0, 0, 0] BP=[1, 0, 0, 0, 0, 0, 0, 0]

enter image description here

Why do they get different information? I'm not too familiar with InterVar, but from the looks of it the data is static from: https://github.com/WGLab/InterVar/tree/master/intervardb

Some files have been recently updated (~3 months ago) while others are a couple years old. Whereas VarSome uses ClinVar from this month's release. My guess is that one of difference is just that VarSome uses more recent information, though presumably you would get the same thing if you updated the InterVar files to equivalent versions.

There's seems to be a difference with BP1:

"Missense variant in a gene for which primarily truncating variants are known to cause disease"

Intevar considers this to be true while Varsome does not, because:

Reason it failed
Missense variant in gene COL4A5 that has 488 pathogenic variants of which 392 pathogenic missense variants = 80.3% which is more than maximum of 10.0%.

ADD COMMENTlink modified 10 months ago • written 10 months ago by manuel.belmadani1.2k

Please check out our automated tool for clinical characterization of germline variants following ACMG guidelines: https://github.com/ding-lab/CharGer

Fernanda Martins Rodrigues | Washington University in St. Louis | Ding Lab

ADD REPLYlink written 6 months ago by mrodrigues.fernanda40
Please log in to add an answer.

Help
Access

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