Question: kegg enrichment using clusterProfiler could not read KEGG Orthology IDs
0
gravatar for jfo
10 days ago by
jfo0
jfo0 wrote:

I wanted to perform the kegg enrichment method. Here's what I did: ran my fasta to kaas to get the KEGG Orthology IDs (K#####) then used these ID to perform enrichKEGG.

  1. An example of the command was performed to test the function.

keggID <- c("K15179", "K03377") x <- enrichKEGG(keggID, organism='ko', keyType='kegg')

--> No gene can be mapped.... --> Expected input gene ID: K10780,K16305,K08093,K14275,K18819,K22554 --> return NULL...

When I use other KO ID like K05164, the program ran without any errors. This could mean that the KEGG IDs (in my example above) could not be mapped because they are not enriched? Also for the IDs that were accepted (e.g. K05164), it didn't have qvalue. I wonder if what I am doing wrong.

  1. My other question is that, are the Kegg orthologous IDs invalid for bitr_kegg() function? I could not convert my K##### IDs. This shows up everytime I try different IDs (except for the example, ofc):

Warning message: In bitr_kegg("K15179", fromType = "kegg", toType = "uniprot", organism = "hsa") : 100% of input gene IDs are fail to map...

This is also true for ncbi protein and geneID. Am I using the wrong organism?

  1. Also, what is the proper way to perform KEGG enrichment for non-model organism? I am really confused as to whether I would just use the closest model organism with available db to perform KEGG or use all available databases?

Thank you for the help!

clusterprofiler • 88 views
ADD COMMENTlink modified 10 days ago by Guangchuang Yu1.1k • written 10 days ago by jfo0
Answer: kegg enrichment using clusterProfiler could not read KEGG Orthology IDs
0
gravatar for Guangchuang Yu
10 days ago by
Guangchuang Yu1.1k
China/Guangzhou/Southern Medical University
Guangchuang Yu1.1k wrote:

keggID <- c("K15179", "K03377") x <- enrichKEGG(keggID, organism='ko', keyType='kegg')

--> No gene can be mapped.... --> Expected input gene ID: K10780,K16305,K08093,K14275,K18819,K22554 --> return NULL...

Usually, this means that you are using the wrong ID.

In your case, you are indeed using the correct ID type. Then such error means that all your input ID have no pathway annotation.

You can verify by visiting the KEGG website.

For example, no pathway annotated in https://www.genome.jp/dbget-bin/www_bget?ko:K15179.

And have pathway annotated for the example ID, https://www.genome.jp/dbget-bin/www_bget?ko:K10780.


Warning message: In bitr_kegg("K15179", fromType = "kegg", toType = "uniprot", organism = "hsa") : 100% of input gene IDs are fail to map...

bitr_kegg is for converting ID within the same species, and K number is not a valid ID type of hsa, albeit that we can find ortholog ID.

Look at the gene session on https://www.genome.jp/dbget-bin/www_bget?ko:K15179, and try:

> bitr_kegg("7469", fromType = "kegg", toType = "uniprot", organism = "hsa")
  kegg    uniprot
1 7469     Q9H3P2
2 7469 A0A0C4DFX9
ADD COMMENTlink modified 10 days ago • written 10 days ago by Guangchuang Yu1.1k

Is it safe to use hsa database for a non-model crustacean? I am still confused how everybody use KEGG analysis for non-model organism. Thank you for the answer!

ADD REPLYlink modified 10 days ago • written 10 days ago by jfo0
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 16.09
Traffic: 260 users visited in the last hour