I am trying to interpret the following annotation from SnpEff:
chr10 96633300 rs4919045 G A 15277.3 PASS ...ANN=A|sequence_feature|MODERATE|PIK3AP1|ENSG00000155629|modified-residue:Phosphotyrosine|ENST00000339364|protein_coding|8/16|c.1376-4807C>T||||||,A|upstream_gene_variant|MODIFIER|PIK3AP1|ENSG00000155629|transcript|ENST00000371109.3|protein_coding||c.-241C>T|||||59|,A|intron_variant|MODIFIER|PIK3AP1|ENSG00000155629|transcript|ENST00000339364.9|protein_coding|8/16|c.1376-4807C>T||||||,A|intron_variant|MODIFIER|PIK3AP1|ENSG00000155629|transcript|ENST00000371110.6|protein_coding|7/15|c.842-4807C>T||||||,A|intron_variant|MODIFIER|PIK3AP1|ENSG00000155629|transcript|ENST00000468783.1|processed_transcript|7/7|n.1022-4807C>T|||||| ...
In particular, I am interested in understanding the MODERATE
effect annotation with the Sequence Ontology feature type of modified-residue:Phosphotyrosine
as the location seems to be intronic. I assume there must be an alternatively-spliced transcript where this is not the case but cannot find it anywhere. How do I interpret/track down the origin of this annotation? Thank you!
what's the reference build you have used for annotation? For GRCh 37 and 38. VEP classifies rs4919045 as modifier for all the transcripts. I got his annotation with SNPeff.
Here is the relevant header from the VCF:
Which does reproduce the above
MODERATE
effect but it doesn't with GRCh38.99. What would be the/a likely rationale for this change? I am guessing the more recent GRCh should be taken as the correct annotation.with SNPeff version SnpEff 5.0e (build 2021-03-09 06:01), for 38.13.Refseq
with GRCh38.99,
Seems so far all the GRCh37, and 38 (refseq and ENSEMBL), result is consistent with your understanding, however it differs with your result. Can you try using some other database of GRCh38 within SNPeff in annotating this variant? Extract just this variant with headers and re-annotate with some other GRCh38 version.