This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ensembl ID differences between assemblies?

Hi,

I have a list of transcripts, such as: ENST00000416191. I want to find the Ensembl IDs for the genes they are associated with. Would I get identical results, whether I query GRCh38 or GRCh37? Note that I am not interested in their coordinates, only the IDs.

Thank you

ensembl
• 8,370 views • link
updated by crisime crisime • written by jiab • jiab

In general, ensembl IDs do not change between assemblies. A new ID may come in (for eg. new finding) and an existing ID (for eg. deprecation) may disappear. But for an existing entity (gene, contig, transcript etc), IDs do not change between assemblies. This would be same for NCBI, UCSC, EBI, Uniprot etc.

• 0 views • link
written by cpad0112 cpad0112

The IDs may stay the same but the information (your entity) they are pointing to can change quite a lot!

Same ENST may point to a

in different builds. Depending on what you are using your ID and the collected data for, this can make a huge difference! Stay carefull and don't expect too much identity in the underlying data, when you search for the same ID in a new database version a few years later.

• 0 views • link
written by crisime crisime

It is unfortunate that in some case it does change even within the same genome version. Here is one example I have mentioned in my previous post.

• 0 views • link
written by EagleEye EagleEye

4 answers

Hi jiab, to give you an idea of the differences here is a little bash script, which gives you a list of ENSTs mapping to different ENSGs in grch37 and grch38:

#!/bin/bash
#download of all ENST mapped to ENSG from Ensembls grch37 biomart
wget -O mapping_37.txt 'http://grch37.ensembl.org/biomart/martservice?query=<Query virtualSchemaName="default" formatter="TSV" header="0" uniqueRows="1" count="" datasetConfigVersion="0.6"><Dataset name="hsapiens_gene_ensembl" interface="default"><Attribute name="ensembl_transcript_id"/><Attribute name="ensembl_gene_id"/></Dataset></Query>'

#download of all ENST mapped to ENSG from Ensembls grch38 biomart
wget -O mapping_38.txt 'http://www.ensembl.org/biomart/martservice?query=  <Query virtualSchemaName="default" formatter="TSV" header="0" uniqueRows="1" count="" datasetConfigVersion="0.6">  <Dataset name="hsapiens_gene_ensembl" interface="default"> <Attribute name="ensembl_transcript_id"/> <Attribute name="ensembl_gene_id"/> </Dataset> </Query>'

#join the two mappings on ENST
join <(sort mapping_37.txt) <(sort mapping_38.txt) >joined_on_ENST.txt

#get just the lines, where ENSGs differ for grch37 and grch38 for the same ENST
awk '$2!=$3 {print $1, $2, $3}' joined_on_ENST.txt >different_ENSG_37_38.txt

If you have a look at the number of lines in each file: 

  wc -l *
 1421 different_ENSG_37_38.txt
 190566 joined_on_ENST.txt
 215170 mapping_37.txt
 218207 mapping_38.txt

It shows you that the mappings in 37 and 38 share about 190k ENST having 25k (37) and 28k (38) unique ENSTs and 1421 ENSTs mapping to different ENSG in both builds. Out of about 243k different ENSTs (190566 +25k+28k) you could search your ENSG for, in the different versions, you would not get the same result for both builds (no ENSG in one build, or different ENSG) for about 54k queries (25k+28k+1421).

• 0 views • link
written by crisime crisime

Very on point answer, thanks!

• 0 views • link
written by jiab • jiab

Most likely the gene IDs would be different because annotations of GRCh37 are based on EnsEMBL v75 and are not updated whereas current annotations of GRCh38 are at version 89. Also, for the same reason, I wouldn't necessarily expect transcript IDs to be conserved. By the way, you can easily test this yourself using the API.

EDIT: Just to clarify in light of the other answers: A transcript or gene ID wouldn't change if the structure of the object it represents doesn't change but as the assemblies are different, genes/transcripts in regions that differ between the two assemblies may end up being different. In short, because the assemblies are different, I wouldn't expect their annotations to be identical.

• 0 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

Just check on the Ensembl browser... For this particular gene, hg38 and hg37 has the same gene-level ID as of Ensembl 89.

• 0 views • link
written by John Ma John Ma

Hello,

as a transcript is bind to a specific gene, you should get the same result in most cases. What can happen is, that a transcript that exists in hg37 does not have to exist in hg38 and vice versa. Also you cannot be sure that transcript with the same id has the same sequence and exons/introns in both reference genomes.

fin swimmer

• 0 views • link
written by finswimmer finswimmer

Log in to answer this question.

More posts like this