【单细胞转录组】KEGG通路反推基因列表

#寻找通路和基因的关系
library(clusterProfiler)
mmu_kegg <- clusterProfiler::download_KEGG("mmu")
names(mmu_kegg)

head(mmu_kegg$KEGGPATHID2NAME)

head(mmu_kegg$KEGGPATHID2EXTID)

PATH2ID <- mmu_kegg$KEGGPATHID2EXTID
PATH2NAME <- mmu_kegg$KEGGPATHID2NAME
PATH_ID_NAME <- merge(PATH2ID, PATH2NAME, by="from")
colnames(PATH_ID_NAME) <- c("KEGGID", "ENTREZID", "DESCRPTION")

head(PATH_ID_NAME)

#转换EntrezID为Gene Symbol
PATH_ID_NAME$geneid <- PATH_ID_NAME$ENTREZID
head(PATH_ID_NAME)
View(map_dt)
# transform id  
map_dt <- bitr(PATH_ID_NAME$ENTREZID, fromType = "ENTREZID",
               toType = c( "SYMBOL"),
               OrgDb = org.Mm.eg.db)

dt_merge <- merge(map_dt, PATH_ID_NAME, by.y = "geneid", by.x = "ENTREZID")

View(dt_merge)