This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VEP vs webserver Polyphen scores

Hello, I have recently noted a discrepancy the scores given by running PolyPhen at its website vs getting the PolyPhen score from ensembl's variant effect predictor (VEP, hg19). Here is one example, with the polyphen scores bolded:

POLYPHEN

VEP

  • Pasted data: AKT2:c.136G>A
  • 1000 Genomes continental allele frequencies: Enabled
  • 1000 Genomes global minor allele frequency: Enabled
  • APPRIS: Enabled
  • BLOSUM62(p): Enabled
  • Condel(p): Enabled
  • Condel Score/prediction(p): Prediction and score
  • CSN(p): Enabled
  • ExAC allele frequencies: Enabled
  • Exon and intron numbers: Enabled
  • Filter by frequency: Exclude variants with MAF greater than 0.05 in 1000 genomes (1KG) combined population
  • Find co-located known variants: Enabled
  • Gene symbol: Enabled
  • HGVS: Enabled LoFtool(p): Enabled
  • MaxEntScan(p): Enabled
  • PolyPhen: Prediction and score
  • Protein: Enabled
  • PubMed IDs for citations of co-located variants: Enabled
  • Get regulatory region consequences: Yes
  • Restrict results: Disabled
  • SIFT: Prediction and score
  • Transcript biotype: Enabled
  • Transcript database to use: Ensembl transcripts
  • Transcript support level: Enabled
  • Polyphen Score for transcript ENST00000392038: 0.805

Which score is more trustworthy?

-edited, as my initial question (why they are different) has been answered (differing protein databases)

polyphen vep annotation missense
• 5,613 views • link
updated by Emily Emily • written by regrant • regrant

1 answer

VEP PolyPhen scores are pre-computed rather than run on the fly. This means that the scores may have used an older version of PolyPhen, and that the protein database that the proteins were compared to will be older. If there has been a massive expansion of a particular protein family, this would result in a change in the PolyPhen score. You can see which version of PolyPhen and when the protein database snapshot was taken in the Ensembl documentation.

• 1 views • link
written by Emily Emily

Sorry,

I uploaded one of my pheno.RNA.passed.somatic.indel.vcf in Vep, I tried to filtrate by Pathogenicity predictions based on SIFT (Prediction and score) and PolyPhen (Prediction and score). But, after analysis done, results shows only empty columns related to SIFT or Polyphen with no score. I am just wondering why?

Thank you

• 1 views • link
written by zizigolu zizigolu

Were any of the variants missense?

• 1 views • link
written by Emily Emily

Sorry I am not sure but I am seeing most is frameshift_variation and frameshift_deletion ; This is link of VEP for one of my .vcf files. Could you please have a look?

I was trying Vep in Linux but I am getting error so that would be nice if I could get results from web tool. However, by web tool I would need to filter by SIFT or polyphen

http://grch37.ensembl.org/Homo_sapiens/Tools/VEP/Ticket?tl=2H9vfkjI8DGQReZk

Thanks a lot

• 1 views • link
written by zizigolu zizigolu

No, there weren't any. I used the filter at the top of the table and filtered by consequence is missense_variant. You only get SIFT and PolyPhen scores for missense variants, and the columns only display if they have data in them.

• 1 views • link
written by Emily Emily

Thanks a lot

Mostly tools visualizing mutation results like lollipop plots needs 

Chromosome  Start_Position  End_Position    Reference_Allele    Variant_Allele

In VEP output I am not able to find protein change column. Is there any protein change column in Vep results please?

• 1 views • link
written by zizigolu zizigolu

Yes, there is a column for amino acid and codon change, which you should be able to see data for for your inframe and frameshift deletions. The columns are only shown when there's data in them, so while you're just viewing the first few variants, which don't have these data associated, you can't see them. Filter by consequences to see the relevant variants and the columns should appear.

You could also get protein domains by ticking the box on the input form. Coming in the next release will also be the option to see the variants on the protein structure, if there is one available.

• 1 views • link
written by Emily Emily

Yeah thanks I saw that for few genes. Is the Vep working on GRCh37 (hg19)? Because when I am loading some columns of my Vep results in cBioPortal for visualization in mutation plot says 

Critical annotation error. Annotation services might be temporarily down, or your input format might be invalid.

These are some of my inputs manipulated from Vep

Chromosome  Start_Position  End_Position    Reference_Allele    Variant_Allele
1   56987   56988   A   -
1   56987   56988   A   -
1   56987   56988   -   -
1   4133415 4133416 -   -
1   4611030 4611030 -   -
1   4872820 4872820 -   -
1   10442174    10442175    A   -
1   10442174    10442175    A   -
1   10442174    10442175    A   T
1   10442174    10442175    A   T
1   10442174    10442175    A   T
1   10442174    10442175    A   T
1   10442174    10442175    A   T
1   10442174    10442175    A   T
1   10442174    10442175    A   T
1   10442174    10442175    A   T
• 1 views • link
written by zizigolu zizigolu

VEP is working just fine. I have no idea about cBioPortal

• 1 views • link
written by Emily Emily

Really thank you so much

Sorry, because I don't have any missense variants, so I don't need to any filtration? I mean can I carry on with all of my variants as I don't have any SIFT or PolyPhen score?

• 1 views • link
written by zizigolu zizigolu

It depends what you're trying to do. Some people filter by the consequences themselves, or the known variant frequencies, or by lists of genes known to be linked to their phenotype of interest.

• 1 views • link
written by Emily Emily

If you export as VCF you'll get the reference in that.

• 0 views • link
written by Emily Emily

Sorry,

I have uploaded my separate SNV and indel as vcf in Vep. I want to score the functional impact of each missense non-synonomous mutation (from 0, non-impactful to 1 highly impactful), synonymous mutation a score of 0 impact and truncating mutations (Non-sense and frameshift mutations) a score of 1 by SIFT and PolyPhen. Finally removing each gene with less than 7 mutations. Is it possible to do these filtrations all at once?

• 0 views • link
written by zizigolu zizigolu

Your job says that you've filtered by frequency:

Filter by frequency: Exclude variants with MAF greater than 0.01 in 1000 genomes (1KG) combined population

Also, please do not hijack old posts to ask continued questions. regrant opened this post over 2.5 years ago: I'm sure they do not want to be receiving notification after notification on your particular problems that are mostly unrelated to their original post. You should have started a new post with your first question, and I should have told you about this back then too.

• 0 views • link
written by Emily Emily

I am sorry

You righ

• 0 views • link
written by zizigolu zizigolu

Log in to answer this question.

More posts like this