open-cravat icon indicating copy to clipboard operation
open-cravat copied to clipboard

Annotation confusion

Open cheanney opened this issue 1 year ago • 1 comments

Hi,

First of all, we want to thank you for this wonderful application. It has streamlined our variant annotation workflow tremendously.

We submitted the HG19 chr16-11349332-CC-TT variant to OpenCravat and received the following results:

cDNA change: c.3_4delinsAA Protein change: p.Val2_Met119del As you can see, the cDNA change does not match the protein change p.Val2_Met119del. Not sure why this variant which is in the initiator methionine is being called a large deletion Val2_Met119del. It is only a 2-nucleotide deletion.

Screenshot 2024-07-30 at 10 55 52 AM

Then, we submitted this variant to VEP and VarSome, and the results were:

cDNA change: c.3_4delGGinsAA Protein change: p.(Met1_Val2delinsIleIle) VarSome Result Screenshot 2024-07-30 at 12 35 43 PM

Here are the cDNA and protein locations from vep:

Screenshot 2024-07-30 at 12 38 20 PM

Let me know if we can provide additional information.

Best,

Anney

cheanney avatar Jul 30 '24 19:07 cheanney

Hi, we will take a look at this issue! Thank you for the information.

jasminebro avatar Jul 31 '24 21:07 jasminebro

Hi Open-Cravat team,

First, we want to express our gratitude for this fantastic tool. Open-Cravat is a vital component of our annotation workflow. That said, we would like to ask if you have had the chance to look into this issue.

Thanks,

Anney

cheanney avatar Dec 09 '24 14:12 cheanney

Hi cheanney,

Apologies for the delay on our end for looking into this.

Our opinion is that the call of Val2_Met119del is correct. chr16-11349332-CC-TT results in a start loss by changing the (reverse strand) sequence from ATGGTA to ATAATA. The next start codon is at position 119.

HGVS recommendations are that, "variants affecting the translation initiation site (Met1), activating an upstream (N-terminal) translation initiation site, are described as deletion-insertion; those activating a downstream (C-terminal) initiation site as a deletion." This specific case is addressed more explicitly here.

Here is the site of the variant, and here is the site of the predicted new start codon.

We think that the VarSome call is incomplete because, while the new sequence ATAATA would translate to IleIle, it's important to annotate that the translation start is lost.

We're very interested in what you think about our observations. Please let us know if you see something wrong with our reasoning in this case. We really appreciate you bringing this to our attention (this was fun to chase down).

kmoad avatar Dec 13 '24 22:12 kmoad