More specifically, after I get the KO numbers for each gene, how do I make the actual assignment of each gene to a specific pathway? I suppose I will have to download a database and parse it somehow, but since there are quite a few of them I'm a bit confused!
Do you have any other method to suggest (except KEGG)?
2 answers
You can do this mapping using custom scripts by downloading the flat file of KO list here. You can use that list to get pathway associated with each KO numbers. If you want the pathway mapping you can check the mapping file here.
If you want to perform enrichment analysis, you may take a look at SubPathwayMiner
AFAIK, these two are the most recent data from KEGG. I am not sure about the reason behind missing KO number. I would recommend you to check with KEGG curators regarding this.
Just for your info... I asked the KEGG curators as you told me Khader and they told me that this is normal; the "ko.list" file contains only KOs that are assigned to a certain pathway. Hence, if a KO doesn't appear in "ko.list" then this is because it wasn't possible for it to be assigned to a pathway (poor characterization).
Thanks for the feedback Panos. This will be useful for the future reference.
You should check the KEGG API through SOAP, also because the access to the KEGG ftp is now restricted without a licence.
This API contains various commands that may help you in mapping a KO ID to a pathway, like the command "get_pathways_by_kos"
Log in to answer this question.
-
KEGG pathway annotation and enrichment analysis on an unannotated organism
written by Raito92
-
How to get KEGG pathway names in KEGG enrichment analysis using clusterprofiler
written by tianshenbio
-
How to find KO number for each protein represented by STRING ID?
written by larcwind •
-
How to generate KEGG pathway classification chart
written by dongchenchen •
-
syntax questions about getting sequences from KO numbers using KEGGREST
written by jon.sy.tarn •
-
JGI files - CSP
written by lmach •
-
Pathway annotation using KEGG
written by sharmatina189059 •
-
How To Batch Grab All Ko (Kegg Ontology) Terms Associated With Each Kegg Pathway?
written by JacobS
-
How To Blast Kegg-Only Data In The 'Nr' Database
written by Shuixia100 •
-
How To See Which Proteins / Enzymes Are Present In A Pathway (Using Resources In Kegg)?
written by Panos
You know that you can also download species-specific pathways?
No, I didn't know that! Nevertheless, I want to find pathways in a metagenome that doesn't exist in KEGG so I have to do it from the beginning.