Question: variant REF-ALT conversion to HGVS.c
0
gravatar for cocchi.e89
16 months ago by
cocchi.e8980
cocchi.e8980 wrote:

Dear all,

I am quite new so sorry if this is a silly question. I am working on a VCF files and comparing my variant with GnomAD DB. I noticed that the VEP column of Gnomad HGVS.c bases are different from the REF-ALT bases.

As example, GnomAD DB row:

21 10944679 rs770820681 C A 3696.48 RF AC=1 ... vep=A|missense_variant|MODERATE|TPTE|ENSG00000166157|Transcript|ENST00000298232|protein_coding|10/23||ENST00000298232.7:c.501G>T|ENSP00000298232.7:p.Arg167Ser...

The variant is a C>A and the VEP column is referring to A allele (first VEP field), but it points to a transcript with a change G>T.

How is that possible? What the "conversion" between them?

Thanks a lot in advance for any help!

variants nomeclature hgvs gene • 591 views
ADD COMMENTlink modified 16 months ago by Pierre Lindenbaum129k • written 16 months ago by cocchi.e8980
1

You might find a tool like "TransVar" to be useful: https://bioinformatics.mdanderson.org/transvar/

Say you start with the HGVS.c notation (ENST00000298232.7:c.501G>T), select "Reverse annotation: cDNA", some reference genome (e.g. GRCh37) and in this case use Ensembl as the annotation database. Enter "ENST00000298232.7:c.501G>T" as the query and it'll give you the different representations between genomic, cDNA, protein variant formats:

input   transcript  gene    strand  coordinates(gDNA/cDNA/protein)  region  info
ENST00000298232.7:c.501G>T  ENST00000298232 (protein_coding)    TPTE    -   chr21:g.10944679C>A/c.501G>T/p.R167S    inside_[cds_in_exon_10] CSQN=Missense;dbsnp=rs770820681(chr21:10944679C>A);reference_codon=AGG;alternative_codon=AGT;aliases=ENSP00000298232;source=Ensembl
ADD REPLYlink written 16 months ago by manuel.belmadani1.2k

Thanks a lot, really helpful!

ADD REPLYlink written 16 months ago by cocchi.e8980
2
gravatar for Pierre Lindenbaum
16 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

How is that possible? What the "conversion" between them?

ENSG00000166157 is mapped on the reverse strand in grch37 http://grch37.ensembl.org/Homo_sapiens/Gene/Summary?g=ENSG00000166157;r=21:10906201-11029719

ADD COMMENTlink written 16 months ago by Pierre Lindenbaum129k
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: 1784 users visited in the last hour