Entering edit mode
John St. John
▴
30
@john-st-john-5628
Last seen 9.6 years ago
I have a gene for example, CDK4
http://www.ncbi.nlm.nih.gov/gene?term=1019, and I was testing its
mapping between entrez ID and ENSEMBL id in R. It's entrez id is 1019.
library(org.Hs.eg.db)
revmap(org.Hs.egENSEMBL2EG)[['1019']]
[1] "ENSG00000135446"
#this appears to be the expected output (see
http://www.ncbi.nlm.nih.gov/gene?term=1019)
org.Hs.egENSEMBL[['11019']]
[1] "ENSG00000121897"
#wrong output, this is a different unrelated gene, LIAS
org.Hs.egSYMBOL[['11019']]
[1] "LIAS"
#expected "CDK4"
Help description for org.Hs.egENSEMBL:
org.Hs.egENSEMBL is an R object that contains mappings between Entrez
Gene identifiers and Ensembl gene accession numbers.
Help description for org.Hs.egENSEMBL2EG (it is the same help
message):
org.Hs.egENSEMBL is an R object that contains mappings between Entrez
Gene identifiers and Ensembl gene accession numbers.
Help description for org.Hs.egSYMBOL:
org.Hs.egSYMBOL is an R object that provides mappings between entrez
gene identifiers and gene abbreviations.
There is something wrong here, right? From the helps it looks like
these different annotation databases do the mappings I expect, and
should all work with entrez ids as input. The fact that one of these
works, and the rest do not is concerning. Is there a more widely used,
less buggy mapping between ENSEMBL and Entrez that I should be using
instead of this one?
Thanks,
John
[[alternative HTML version deleted]]