Open growland2 opened 1 month ago
Hi @growland2,
Sorry to hear about this inconvenience.
While running VEP with ClinVar as a custom file (similar to your command), I do get the same results consistently for the expected variants, regardless of their Ensembl/RefSeq transcript.
I am going to try and see if any of the options you are using could be affecting the results.
Best regards, Nuno
Thanks, FYI running these variants through VEP GRCh37 online failed to annotate with any clinical significance for any transcript for both example variants:
VEP online results for Example variant 1 and Example variant 2.
Kind regards,
Greg
Describe the issue
When running VEP using a local implementation via Docker, we see that variants are annotated with different ClinVar clinical significance values depending on the corresponding transcript or feature. Examples of such variants are shown below:
Example variant 1:
When annotated via VEP locally using VEP version=107, VEP Cache=RefSeq and ClinVar GRCh37 version=20240317 (as a custom annotation source), this variant was not annotated with a ClinVar_CLNSIG for the given transcripts (see excerpt from VEP output below):
VEP v107 and RefSeq cache output:
Command:
We then re-ran locally using VEP version=v112, VEP Cache=RefSeq and ClinVar GRCh37 version=20240317, to ensure the prior output was not specific to v107. This returned the same results as above.
Next, we ran specifying the merged cache containing RefSeq and Ensembl transcripts. Here we saw that we only get the expected ClinVar_CLNSIG annotation
"Pathogenic/Likely_pathogenic"
for some of the Ensembl (ENS) transcripts (ENST00000469403, ENST00000477585 and ENST00000479707) but not for the remaining transcripts (ENST00000379370, NM_001305275.2, NM_001364727.2 and NM_198576.4) (see output from VEP below).VEP v112 and merged cache output:
Command:
Example variant 2:
As was seen for "Example variant 1", when annotated via VEP locally using VEP version=107, VEP Cache=RefSeq and ClinVar GRCh37 version=20240317, this variant received different ClinVar_CLNSIG values depending on the trascript/feature with NM_022834.5 and NM_199121.3 receiving no ClinVar_CLNSIG annotation, and NR_125994.1, NR_125995.1 and NR_125996.1 being annotated with the expected
"Pathogenic/Likely_pathogenic"
(see excerpt from VEP output below):VEP v107 and RefSeq cache output:
What decides which transcript receives the expected ClinVar_CLNSIG value?
Additional information
System