enrichKEGG results for non-model organism
0
0
Entering edit mode
AP • 0
@ap-14194
Last seen 2.4 years ago
United States

Hey everyone,

In order to do gene enrichment analysis on non-model organism (fungi), I obtained the KEGG ids using KOfam analysis. Using clusterProfiler enrichKEGG function I did the enrichment analysis. But the problem now is the enriched pathways I am getting is human diseases like Parkinsons disease, Huntington diseases etc. or human related pathways. Is there a way to filter these out or I am doing something wrong. My codes are given below:


file <- read.delim("T161_core.txt", header = FALSE, sep="\t")
colnames(file) <- c("ProteinID", "KO")
file1<- (file$KO)

kegg <- enrichKEGG(  file1,  organism = "ko",  keyType = "kegg",  pvalueCutoff = 0.05,
                     pAdjustMethod = "BH", qvalueCutoff = 0.05)

p<- dotplot(kegg, showCategory=20) + ggtitle("T161 Core ")
p + theme(
  plot.title = element_text(color="black", size=14, face="bold", hjust=0.5))



sessionInfo( )

Top 20 pathways

clusterProfiler enrichment KEGG • 2.4k views
ADD COMMENT
1
Entering edit mode

If anyone still wondering about this:

You can filter your 'kegg' object based on a list of target pathways.

Look at these fields:

kegg@result[["Description"]], which has pathway names

and

kegg@result[["ID"]], which has respective pathway IDs

ADD REPLY
0
Entering edit mode

can you tell me the formate of output file from KOfam please? I don't understand the formate of input file for ClusterProfiler.

ADD REPLY
0
Entering edit mode

Hello, I'm also analyzing a non-model organism. Could you tell me which input files you used? Also, why did you set `organism = 'ko' ?

ADD REPLY

Login before adding your answer.

Traffic: 768 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6