This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get pathway information from KEGG using Entry id

Hi, I have a list of ids from KEGG, is there any way to get the corresponding pathway information in a single search?

Kegg_ids are - stm:STM2395, sau:SA0132, efa:EF2297

Expected output stm:STM2395 : Pathway : Cationic antimicrobial peptide (CAMP) resistance

Thanks in advance

r sequencing gene next-gen snp
• 4,678 views • link
updated by EagleEye EagleEye • written by sinumolgeorge • sinumolgeorge 
> library(KEGGREST)
> kegg_ids=read.csv("kegg_ids", header = F, stringsAsFactors = F)
> kegg_ids
            V1
1 stm:STM2395 
2 sau:SA0132
3 efa:EF2297

> data.frame(Pathway=unlist(sapply(keggGet(kegg_ids[,1]), "[[", "PATHWAY")),stringsAsFactors = F)
                                                      Pathway
    stm01503 Cationic antimicrobial peptide (CAMP) resistance
    stm02020                             Two-component system
    efa00550                       Peptidoglycan biosynthesis
    efa01100                               Metabolic pathways
    efa01502                            Vancomycin resistance
    efa02020                             Two-component system
• 0 views • link
written by cpad0112 cpad0112

Thank you, but its giving pathway id, how we know the corresponding kegg id ?

• 0 views • link
written by sinumolgeorge • sinumolgeorge

I want to get the kegg_id along with the result. Are there any options?

• 0 views • link
written by sinumolgeorge • sinumolgeorge 
library(KEGGREST)
library(purrr)
library(magrittr)


kegg_ids=read.csv("kegg_ids", header = F, stringsAsFactors = F)
kegg_ids

kegg_pathways=data.frame(Pathway=unlist(sapply(keggGet(kegg_ids[,1]), "[[", c("PATHWAY"))),stringsAsFactors = F)
kegg_list=map(keggGet(kegg_ids[,1]), extract, c("ENTRY", "PATHWAY"))
library(dplyr)
kegg_df=bind_rows(lapply(kegg_list, function (x) data.frame(t(unlist(x)),stringsAsFactors = F)))
library(tidyr)
kegg_df1=na.omit(gather(kegg_df,key = "", value = Pathway, -ENTRY.CDS))[,c(1,3)]

Note: beware of library loading. Some of the libraries mask the function of others and results in execution issues. output:

> na.omit(gather(kegg_df,key = "", value = Pathway, -ENTRY.CDS))[,c(1,3)]
   ENTRY.CDS                                          Pathway
1    STM2395 Cationic antimicrobial peptide (CAMP) resistance
4    STM2395                             Two-component system
9     EF2297                       Peptidoglycan biosynthesis
12    EF2297                               Metabolic pathways
15    EF2297                            Vancomycin resistance
18    EF2297                             Two-component system
• 0 views • link
written by cpad0112 cpad0112

Error in UseMethod("extract_") : no applicable method for 'extract_' applied to an object of class "list"

• 0 views • link
written by sinumolgeorge • sinumolgeorge

This is because of one library masking the function of another library. A note was added in between regarding the same. Following is the updated code with correct order of libraries:

library(KEGGREST)
kegg_ids=read.csv("test.txt", header = F, stringsAsFactors = F)
library(purrr)
library(magrittr)
kegg_list=map(keggGet(kegg_ids[,1]), extract, c("ENTRY", "PATHWAY"))
library(dplyr)
kegg_df=bind_rows(lapply(kegg_list, function (x) data.frame(t(unlist(x)),stringsAsFactors = F)))
library(tidyr)
kegg_df1=na.omit(gather(kegg_df,key = "", value = Pathway, -ENTRY.CDS))[,c(1,3)]
• 0 views • link
written by cpad0112 cpad0112

I have used GENE NAME instead of pathway information for getting corresponding gene name. But it's not fetching output. Is there any other keyword for fetching corresponding gene name of KEGG id?

• 0 views • link
written by sinumolgeorge • sinumolgeorge

NAME for gene name

library(KEGGREST)
library(purrr)
library(magrittr)

kegg_ids=read.csv("test.txt", header = F, stringsAsFactors = F)
library(biomaRt)

kegg_pathways=data.frame(Pathway=unlist(sapply(keggGet(kegg_ids[,1]), "[[", c("PATHWAY"))),stringsAsFactors = F)
kegg_list=map(keggGet(kegg_ids[,1]), extract, c("ENTRY", "PATHWAY", "NAME"))
library(dplyr)
kegg_df=bind_rows(lapply(kegg_list, function (x) data.frame(t(unlist(x)),stringsAsFactors = F)))
library(tidyr)
na.omit(gather(kegg_df,key = "", value = Pathway, -c(ENTRY.CDS,NAME)))[,c(1,2,4)]

output:

> na.omit(gather(kegg_df,key = "", value = Pathway, -c(ENTRY.CDS,NAME)))[,c(1,2,4)]
   ENTRY.CDS  NAME                                          Pathway
1    STM2395  pgtE Cationic antimicrobial peptide (CAMP) resistance
4    STM2395  pgtE                             Two-component system
9     EF2297 vanYB                       Peptidoglycan biosynthesis
12    EF2297 vanYB                               Metabolic pathways
15    EF2297 vanYB                            Vancomycin resistance
18    EF2297 vanYB                             Two-component system
• 0 views • link
written by cpad0112 cpad0112

Thanks! Very useful. How if I want to list out all the results even with the ones that do not have any pathways but still have name or definition?

• 1 views • link
written by hadiaziz • hadiaziz

1 answer

You can download all KEGG pathways with ids, description and corresponding genes involved as a simple table (plain text file) using GeneSCF.

• 0 views • link
written by EagleEye EagleEye

Log in to answer this question.

More posts like this