How to create a customized OrgDb in order to perform kegga analysis?
1
0
Entering edit mode
baus87 • 0
@baus87-13757
Last seen 6.5 years ago

Hi,

i'm working on a set of DEGs and i'd like to analyze them using the kegga function. Since i'm working on a species whose OrgDb is not implemented in, how can i create an in-house dedicated OrgDb?
thanks in advance,

Antonio

limma kegga orgdb • 2.1k views
ADD COMMENT
2
Entering edit mode
@gordon-smyth
Last seen 2 hours ago
WEHI, Melbourne, Australia

You don't need an OrgDb package to run kegga(). kegga simply reads annotation from the KEGG website.

However KEGG does not provide annotation for the Trichoderma gamsii species, as you can see from their website:

http://www.kegg.jp/kegg/catalog/org_list.html

The closest species is Trichoderma reesei.

So you can't do a KEGG enrichment analysis for Trichoderma gamsii because KEGG pathways don't exist for that species.

ADD COMMENT
0
Entering edit mode

Dear Gordon,
i'm working with Trichoderma gamsii (filamentous fungus). I'd like to perform a kegg enrichment analysis upon DEGs. Is it kegga suitable for this purpose?

ADD REPLY
0
Entering edit mode

Ok, got it.
Could you recommend me any alternative tool?

ADD REPLY
0
Entering edit mode

Let me ask you in turn: is there any pathway / annotation info publicly available for this organism? At KEGG apparently not... Being a person working with T gamsii I suppose you should know this.... If there isn't, you cannot expect members of this (Bioconductor) support site to give you a useful answer to your last question....

ADD REPLY
0
Entering edit mode

Dear Guido,
nowadays i'm annotating my reference genomes with ko codes from KEGG. Given a set o differentially expressed genes, I was just looking for the proper enrichment analysis to be applied. Maybe my first post has been misleading, sorry for that.
Antonio

ADD REPLY
1
Entering edit mode

Aha, AFAIK in that case you can use the argument species.KEGG = "ko" when using the function kegga().

Assuming DEGs is a character vector with the differentially expressed KO IDs, something like:

k<- kegga(DEGs, species.KEGG = "ko", FDR = 0.05)
ADD REPLY

Login before adding your answer.

Traffic: 669 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