1
$\begingroup$

3 solutions that I had to unfortunately rule out.

Potential solution #1

As far as I know, the conventional way is to do this is by using CCDS ids (e.g. using biomart). However, the issue is that, as I see on the CCDS project page, the last release of CCDS ids was made almost 3 years ago (!) and it seems to be compatible with Ensembl release 92. Then how should one convert ids to the latest Ensembl release e.g. 100?

CCDS Release Species NCBI Annotation Release Ensembl Annotation Release Assembly Name Assembly ID Made Public
22 Homo sapiens 109 92 GRCh38.p12 GCF_000001405.38 06/14/2018

Potential solution #2

I am also aware of Ensembl's ID History Converter. However, I am looking for a programmatic approach. Also, I noticed that if I try to convert ~>20K ids with ID History Converter, in my experience, there is always some issue. e.g. an error I got when I tried to convert ~70K ensembl transcript ids. enter image description here

As a side note, I wonder why this ID conversion functionality is not available via Ensembl's REST API. To me, id conversion seems to be a basic enough task that almost every bioinformatician has to do. So it would be immensely helpful to be someday able to do this programmatically using Ensembl's mighty a REST API.

Potential solution #3

Finally, as another potential solution, over at Ensembl's FTP website, located within the directory of the release to which I am interested in converting my ids to (release 100, target release), I saw that according to a README file, "ensembl_archive_" directory should contain "Data on historical Ensembl releases". So, I wondered if this directory contains a map between ids from different releases. However, without any extra info about column names, contents, file formats etc, I couldn't really figure out how to use those files. So, I had to rule this one out too.

Index of /pub/release-100/mysql/ensembl_archive_100/
../
CHECKSUMS                                          26-Mar-2020 14:31                 132
ens_release.txt.gz                                 26-Mar-2020 14:31                1098
ensembl_archive_100.sql.gz                         26-Mar-2020 14:31                 943
release_species.txt.gz                             26-Mar-2020 14:31               39620
species.txt.gz                                     26-Mar-2020 14:31                6669

Of note, I saw a similar "mysql" directory over at the directory of the release whose ids I want to convert (source release): http://ftp.ensembl.org/pub/grch37/release-100/mysql/homo_sapiens_core_100_37/ . However, again, I do not understand its contents without any extra information.

So, finally, I am not sure what is the best way to programmatically convert Ensembl ids from between different releases. As a test case scenario, I want to convert all the protein_coding transcript ids from GRCh37 (release 75) to GRCh38 (release 100). If anybody has any suggestions about this, please kindly let me know.

$\endgroup$
1
+50
$\begingroup$

The ID history converter is available as a Perl script.

This script accesses these MySQL tables, so you could query these tables directly.

Note that the mapping session is only one release at a time, so you would need to run 25 queries to go from release 75 to 100.

$\endgroup$
2
  • $\begingroup$ Thank you for your answer. To use the perl script, I wonder if there is any tutorial to refer to. I am particularly interested in knowing how to "tell" the script to convert from a given specific release to other specific release e.g. from 75 to 76. $\endgroup$
    – Ramirez
    Mar 22 at 14:38
  • 1
    $\begingroup$ It just gives you the whole history, you would need to parse out the releases you want. $\endgroup$ Mar 22 at 14:57
0
$\begingroup$

Here's one way to do it programmatically using the biomaRt package (your potential solution #1):

library(biomaRt)
library(dplyr)

mart37 <- useMart(biomart = "ENSEMBL_MART_ENSEMBL",
                  host    = "grch37.ensembl.org",
                  path    = "/biomart/martservice",
                  dataset = "hsapiens_gene_ensembl")

mart38 <- useMart(biomart = "ENSEMBL_MART_ENSEMBL",
                  path    = "/biomart/martservice",
                  dataset = "hsapiens_gene_ensembl")

#pick all genes on chr22 from build 37 as a test dataset::
out37 <- getBM(attributes = c('hgnc_symbol', 'chromosome_name',
                              'ensembl_gene_id', 'ensembl_transcript_id', 'ccds'),
                  filters = c('chromosome_name'),
                  values = list("22"),
                  mart = mart37) %>% #empty string CCDS ID's
  filter(ccds != "")

#NOTE:: to check which filters can be used to query mart37, run this command:
#View(listFilters(mart37))

out38 <- getBM(attributes = c('hgnc_symbol', 'chromosome_name',
                              'ensembl_gene_id', 'ensembl_transcript_id', 'ccds'),
                        filters = c('ccds'),
                        values = list(out37$ccds),
                        mart = mart38)

What this does is query the build 38 DB via the CCDS ID's of build 37. You can include other variables in this out38 dataframe by simply adding to the attributes argument above. And in order to see all of the potential variables that can be included, you just have to run View(listFilters(mart37)) or View(listFilters(mart38)), depending on which build you are interested in.

Does this solve your problem?

$\endgroup$
5
  • $\begingroup$ your solution could indeed be useful. I am just curious which ensembl release numbers are used when mart37 and mart38 were built. Is there a way to specify that? $\endgroup$
    – Ramirez
    Mar 25 at 22:04
  • $\begingroup$ mart37 corresponds to the GRCh37 release, while mart38 corresponds to GRCh38. This is specified in the host= argument in the code above. mart38 is the default, and so is not specified. Does this answer your question? $\endgroup$
    – h3ab74
    Mar 26 at 13:10
  • $\begingroup$ P.S. if you're referring to the version number of the ensembl release, you can obtain this by adding 'ensembl_gene_id_version' to the attributes array in the code above. $\endgroup$
    – h3ab74
    Mar 26 at 13:56
  • $\begingroup$ By Ensembl release number I meant version numbers of the assemblies maintained by ensembl, e.g. release 75 corresponds to GRCh37 for human genome. $\endgroup$
    – Ramirez
    Mar 26 at 22:53
  • $\begingroup$ ensembl_gene_id_version could be useful if it is possible to know to which release that belongs to. e.g. how would I know if x.3 id belongs to genome release 90 or 91? To get this information, I would have to run 'Ensembl IDmapper' which I would like to avoid, because of the issue I mentioned in my question. $\endgroup$
    – Ramirez
    Mar 26 at 22:57

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.