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,
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:
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.
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?
Ok, got it.
Could you recommend me any alternative tool?
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....
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.
Aha, AFAIK in that case you can use the argument species.KEGG = "ko" when using the function kegga().
species.KEGG = "ko"
Assuming DEGs is a character vector with the differentially expressed KO IDs, something like:
k<- kegga(DEGs, species.KEGG = "ko", FDR = 0.05)
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy