Question: Wrong predict consequence in variants using VEP Ensembl
1
gravatar for theodorobiotec
3.3 years ago by
theodorobiotec10 wrote:

Hi,

I've been using Ensembl Variant Effect Predictor (VEP) to determine the effect of variants in non-model organism. Despite of being a collection of variants identified in non-model organisms it has a reference genome and a VEP database hosted in ensembl fungi portal.

So, I ran vep using a multi-sample VCF as input and got a final VCF with annotations about the consequence of each variant. When I was looking at the results I realised a variant that supposedly result in a stop codon insertion in some samples as you can check in first line bellow.

Allele      Consequence   IMPACT     SYMBOL             Gene Feature_type          Feature        BIOTYPE EXON
2055      T      stop_gained     HIGH ATEG_09961 CADATEAG00000507   Transcript CADATEAT00000507 protein_coding  8/8
2059      C missense_variant MODERATE ATEG_09961 CADATEAG00000507   Transcript CADATEAT00000507 protein_coding  8/8
     INTRON HGVSc HGVSp cDNA_position CDS_position Protein_position Amino_acids  Codons Existing_variation DISTANCE
2055                             8260         8260             2754         Q/* Cag/Tag                            
2059                             8262         8262             2754         Q/H caG/caC                            
     STRAND SYMBOL_SOURCE
2055      1         CADRE
2059      1         CADRE

 

But if you pay attention there's also a second variant into same codon (same Protein_position) but two bases apart from the first variant. Don't VEP should take into account the variants into same codon before output the consequence results?

Because if we consider the second's variant role and supposing that a sample is homozygote for variant allele in both variants the consequence will be the insertion of a TaC (translate to tyrosine).

I don't know if I did something wrong or miss any parameter in my analysis but it sounds a bit strange for me that VEP's algorithm doesnt handle with issues related to variants located at same codon.

snp vep ensembl software error • 1.8k views
ADD COMMENTlink modified 3.3 years ago by Christian2.8k • written 3.3 years ago by theodorobiotec10
3
gravatar for Emily_Ensembl
3.3 years ago by
Emily_Ensembl18k
EMBL-EBI
Emily_Ensembl18k wrote:

The VEP does not consider variant phasing and only considers each variant individually.

ADD COMMENTlink written 3.3 years ago by Emily_Ensembl18k
3
gravatar for Christian
3.3 years ago by
Christian2.8k
Cambridge, US
Christian2.8k wrote:

We built CooVar a while back for exactly that purpose/reason (and uncovered a few discrepancies with complex variants in VEP of which I don't know if they have been resolved).

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by Christian2.8k

Very cool - well done for creating this! :)

ADD REPLYlink written 3.3 years ago by John12k
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: 996 users visited in the last hour