<div dir="ltr"><div>Hi Andrew,</div><div><br></div><div>I have been already in contact with someone from the VEP Ensembl team, and she indeed super kindly addressed some of my issues and keep me updated of your roadmap about HGVS conversions. I however thought that my question here was not part of that 'todo' list (and that the whole forum may benefit from the potential answer), apologies for the double check, it was not on purpose.</div><div><br></div><div>In any case, and just FYI, I m using VEP for that conversions because i need the result in 'vcf  representation', meaning that i need chr-pos-ref-alt specifically stated (eg I would need additional steps to retrieve the ref from <span style="font-family:Menlo;font-size:11px;background-color:rgb(255,255,255)"><font size="2"><span style="font-family:arial,helvetica,sans-serif">the 'g.</span></font><font size="2"><span style="font-family:arial,helvetica,sans-serif">32336340_32336341delinsAA' result). I have been not able to find a tool supporting that among your tools (hope I m wrong?).</span></font></span></div><div><span style="font-family:Menlo;font-size:11px;background-color:rgb(255,255,255)"><font size="2"><span style="font-family:arial,helvetica,sans-serif"><br></span></font></span></div><div><span style="font-family:Menlo;font-size:11px;background-color:rgb(255,255,255)"><font size="2"><span style="font-family:arial,helvetica,sans-serif">Thanks for your answer anyway!</span></font></span></div><div><span style="font-family:Menlo;font-size:11px;background-color:rgb(255,255,255)"><font size="2"><span style="font-family:arial,helvetica,sans-serif">and have a good weekend<br></span></font></span></div><div><span style="font-family:Menlo;font-size:11px;background-color:rgb(255,255,255)"><font size="2"><span style="font-family:arial,helvetica,sans-serif">d</span></font><br></span></div></div><br><div class="gmail_quote"><div dir="ltr">El vie., 7 dic. 2018 a las 18:16, Andrew Parton (<<a href="mailto:aparton@ebi.ac.uk">aparton@ebi.ac.uk</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">Hi David,<div><br></div><div>Thanks for your query - while VEP supports HGVS input, due to the complexity of HGVS and the variety of ways in which people use it, then we require that input HGVS is relative to genomic or transcript coordinates. We have some documentation (that we are in the process of improving) on the matter here: <a href="http://www.ensembl.org/info/docs/tools/vep/vep_formats.html#hgvs" target="_blank">http://www.ensembl.org/info/docs/tools/vep/vep_formats.html#hgvs</a></div><div><br></div><div>It may be more appropriate for you to use our variant recoder tool, also contained within the ensembl-vep repository, rather than VEP itself, as this will give you HGVSg output. However, without a particular transcript within the input, then it is possible that the variant could map to multiple locations. Variant Recoder will often guess at a solution (in the BRCA2 example you gave then it suggests <span style="font-family:Menlo;font-size:11px;background-color:rgb(255,255,255)">NC_000013.11:g.32336340_32336341delinsAA)</span>, however as you have noticed, unless the correct genomic/transcript nomenclature is used then a result cannot be guaranteed.</div><div><br></div><div>We are currently in the process of improving how we handle HGVS inputs, however this work is still in development. Sorry that I couldn’t give you a more beneficial response.</div><div><br></div><div>Kind Regards,</div><div>Andrew</div><div><br></div><div><div><br><blockquote type="cite"><div>On 5 Dec 2018, at 09:27, David Tamborero <<a href="mailto:david.tamborero@gmail.com" target="_blank">david.tamborero@gmail.com</a>> wrote:</div><br class="gmail-m_-3339877667212711762Apple-interchange-newline"><div><div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi there,</div><div><br></div><div>I m badly interested in converting protein changes to genomic changes via the VEP  standalone perl script. Of note, I have the gene (but not the specific transcript) in which such a change is annotated (eg NRAS:Q61L).</div><div><br></div><div>I m experiencing some problems that I cannot figure it out how to solve; I ve tried to find any documentation to address those in the dev list archives etc but I failed. So I hope this is a good way to look for help (thanks in advance!). <br></div><div><br></div><div>I m using VEP 93.3 with hg19 (I do not think that the vep parameters are relevant, but note that I m using vcf as output format since I need the coordinates in chr-pos-ref-alt   ---and not HGVS---  format). <br></div><div><br></div><div>-- insertion of stop codons--</div><div><br></div><div>so for cases as CDKN2A:p.Q50* or PTEN:p.R233*, VEP works smooth; however, for other cases (eg BRCA2:p.S662*) VEP does not seem happy and does not give an output. <br></div><div><br></div><div>My guess is that this has to be with the fact that the ones that work are stop codons caused by a specific nucleotide change (eg. CDKN2A:p.Q50* -> chr10 89717672 C>T) whereas those that do not work are caused by insertions, which I guess it gives a larger universe of possible nucleotide changes. <br></div><div><br></div><div> However, why not to have a 'more possible' guess? e.g. for BRCA2:p.S662*, TransVar gives a fair chr13:g.32910477_32910478delCTinsAA try.</div><div><br></div><div>--indels--</div><div><br></div><div>here I m failing badly. I m not even sure that I m using the correct nomenclature. For instance, for frameshifts, DTX1:p.P11fs*2 is giving a 'Unable to parse HGVS notation' error. I went to the variant recoder API REST, and this specific variant in cdna nomenclature (which i retrieved also by TransVar), is represented as <span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">"ENSP00000257600.3:p.Pro11LeufsTer2" in protein coordinates<br></span></span></div><div><br></div><div><a href="https://rest.ensembl.org/variant_recoder/human/DTX1:c.32delC?content-type=application/json" target="_blank">https://rest.ensembl.org/variant_recoder/human/DTX1:c.32delC?content-type=application/json</a></div><div><br></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">but when I input this representation to the VEP command, I m having the format error anyway, even after changing <span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">ENSP00000257600 for the gene  (DTX1) or transcript (ENST00000257600) <span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">name</span></span></span></span> plus some alterntives to the <span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">p.Pro11LeufsTer2 representation.</span></span></span></span></span></span></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><br></span></span></span></span></span></span></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">Am I missing some point?</span></span></span></span></span></span></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><br></span></span></span></span></span></span></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">many thanks in advance (and congratulations for such a great tool!)</span></span></span></span></span></span></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">br</span></span></span></span></span></span></div><div><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox"><span><span class="gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox-string gmail-m_-3339877667212711762m_-5679131016372609052gmail-objectBox">d<br></span></span></span></span></span></span></div></div></div></div></div></div></div></div></div>
</div></div>
_______________________________________________<br>Dev mailing list    <a href="mailto:Dev@ensembl.org" target="_blank">Dev@ensembl.org</a><br>Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br></div></blockquote></div><br></div></div>_______________________________________________<br>
Dev mailing list    <a href="mailto:Dev@ensembl.org" target="_blank">Dev@ensembl.org</a><br>
Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" rel="noreferrer" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" rel="noreferrer" target="_blank">http://www.ensembl.info/</a><br>
</blockquote></div>