Question

How to get a mapping between KEGG module and KEGG orthologs?

0

Entering edit mode

3.7 years ago

O.rka ▴ 740

I'm looking for the most simple of tables but it's difficult to find. A table of KEGG orthologs in the following format: [MODULE]\t[KO_1, KO_2, ..., KO_N]

I downloaded a weird formatted flat file from KEGG but for some of the KEGG modules there were other KEGG modules in the hierarchy (yes, I know KEGG is hierarchical) such as https://www.genome.jp/kegg-bin/show_module?M00615

Does anyone know where I can find this? I just need a very simple table for set comprehension.

KEGG Annotation Ortholog Module • 2.5k views

ADD COMMENT • link updated 3.7 years ago by Elucidata ▴ 270 • written 3.7 years ago by O.rka ▴ 740

score 4 · Accepted Answer · 2021-03-18

You can use R to construct the table. You would need to install the package and load it. You can use the code below to connect to the KEGG database, retrieve module information, map to get corresponding ortholog information, and construct the table.

#Install package to get relevant information from KEGG database

if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("KEGGREST")

#Note: Uncomment the above code once the installation is successful

#Load package
library(KEGGREST)

#Get list of modules in KEGG
mod <- keggList("module")

#Loop through each module and get corresponding orthologs
#Return module ID, corresponding list of orthologs
#obj is a list of list(moduleID,orthologs)
#[[moduleID1,orthologs list 1],[moduleID2,orthologs list 2].etc]

obj<-lapply(names(mod),function(x)
{
  module<-strsplit(x,"md:")[[1]][2]
  #Search for corresponding ortholog
  ko<-keggGet(x)
  #Save list of orthologs as a string separated by ","
  orthologs<-paste(names(ko[[1]]$ORTHOLOGY),collapse = ",")
  list(module,orthologs)
})

#Convert list to dataframe
df<-do.call(rbind,obj)

#Name columns
colnames(df)<-c("Module","KO")

#Display first few entries in the table
head(df)

#Save table to csv file
write.csv(df,path to file/filename.csv)