Search
Question: org.Hs.eg.db - Error in mapping gene symbol to gene id?
0
gravatar for af2202
16 months ago by
af22020
af22020 wrote:

Please look at the following code:

> library("org.Hs.eg.db")

> get("MAK", org.Hs.egALIAS2EG)
$MAK
[1] "4117"  "57538"

> get("ALPK3", org.Hs.egALIAS2EG)
$ALPK3
[1] "57538"

As far as I can tell, "MAK" is correctly mapped to the gene ID 4177 and incorrectly mapped to the gene id 57538 (the latter is the id corresponding to gene symbol "ALPK3").

Am I misunderstanding/misusing something?

Aris Floratos
Columbia University

ADD COMMENTlink modified 16 months ago by James W. MacDonald45k • written 16 months ago by af22020
0
gravatar for James W. MacDonald
16 months ago by
United States
James W. MacDonald45k wrote:

You are confusing ALIAS and SYMBOL. The gene that currently has the HUGO symbol of ALPK3 has been known in the past as MAK. So if you ask for the ALIAS, you get all the genes that have ever been called MAK, and that includes the one that is currently called MAK, as well as ALPK3, which used to be called MAK.

If you just want the current SYMBOL, then you should ask for that. And you should also stop using get, and start using select.

> select(org.Hs.eg.db, c("MAK","ALPK3"), "ENTREZID","SYMBOL")
'select()' returned 1:1 mapping between keys and columns
  SYMBOL ENTREZID
1    MAK     4117
2  ALPK3    57538
> select(org.Hs.eg.db, c("MAK","ALPK3"), "ENTREZID","ALIAS")
'select()' returned 1:many mapping between keys and columns
  ALIAS ENTREZID
1   MAK     4117
2   MAK    57538
3 ALPK3    57538

 

ADD COMMENTlink written 16 months ago by James W. MacDonald45k

I got it, thank you...

ADD REPLYlink written 16 months ago by af22020
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 120 users visited in the last hour