Gene Symbol to KEGG
3
0
Entering edit mode
8.0 years ago
Folder40g ▴ 190

Hi!

I have a long list (R) which contains genes, compounds and other things (gex; glycan).

Example:

"CDK4" "C01245" "PIK3R5" "C05981" "C01245" "MAPK1" "E2F1"

I'm working with some pathways retrieved from KEGG and I want to know which elements from this list are genes and which are not.

I need to convert gene symbols to KEGG id, so far compounds and other that are not genes can be identified, because those are already in KEGG IDs.

names((keggGet("C00044"))[[1]]$ENTRY)

[1] "Compound"

names((keggGet("G13064"))[[1]]$ENTRY)

[1] "Glycan"

But this does not work with gene symbols.

names((keggGet("KRAS"))[[1]]$ENTRY)

Error in .getUrl(url, .flatFileParser) : Bad Request (HTTP 400).

I have tried using library(org.Hs.eg.db), but so far I have not been able to get the results that I expected.

Any hint?

Thanks for your time.

KEGG Symbol Gene • 11k views
ADD COMMENT
0
Entering edit mode
8.0 years ago
ivivek_ngs ★ 5.2k

You should convert the gene symbols to ENTREZ ids as they are not unique and then try to map them to the corresponding KEGG IDs Take a look at this answer here, this should serve the purpose for you.

ADD COMMENT
0
Entering edit mode
8.0 years ago
Girolamo ▴ 140

If the genes you are interested in are few you can use: http://rest.kegg.jp/find/genes/KRAS otherwise, you will save time by converting the gene symbol in ensembl/uniprot/entrez/ncbi id and then retrieve the pathway.

ADD COMMENT
0
Entering edit mode
8.0 years ago
natasha.sernova ★ 4.0k

See this post:

A: Ensembl ID to Gene Symbol

and this link inside:

https://biodbnet-abcc.ncifcrf.gov/db/db2db.php

It may help you.

ADD COMMENT

Login before adding your answer.

Traffic: 2018 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6