Question: variant REF-ALT conversion to HGVS.c
0
gravatar for cocchi.e89
12 months ago by
cocchi.e8950
cocchi.e8950 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 • 470 views
ADD COMMENTlink modified 12 months ago by Pierre Lindenbaum127k • written 12 months ago by cocchi.e8950
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 12 months ago by manuel.belmadani1.2k

Thanks a lot, really helpful!

ADD REPLYlink written 12 months ago by cocchi.e8950
2
gravatar for Pierre Lindenbaum
12 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum127k 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 12 months ago by Pierre Lindenbaum127k
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: 1638 users visited in the last hour