Question: How to perform KEGG enrichment analysis for DE genes from a de-novo transcriptome?
0
gravatar for Light92
15 months ago by
Light9250
Italy
Light9250 wrote:

Hello, I'm running a pipeline to analyze differential expression for transcripts predicted and BLAST-annotated from a de-novo assembled transcriptome... At the end of the pipeline, I'm interested in a KEGG Enrichment analysis, but I have no idea about how to perform it. I used Trinity till now but no functions for this kind of analysis are supported, so I guess I need an external script. But I was unable to find a reliable strategy to perform it with external functions.

More specifically, most functions are designed to work for a specific organism, when my transcripts are annotated against a lot of species...

Any suggestions?Thanks in advance!

kegg kegg enrichment • 647 views
ADD COMMENTlink modified 15 months ago • written 15 months ago by Light9250
2
gravatar for GenoMax
15 months ago by
GenoMax94k
United States
GenoMax94k wrote:

when my transcripts are annotated against a lot of species...

Transcripts should be annotated for function not other species. Sure the "hits" from blast searches or whatever you may have done are to other species but based on the homology you need to decide if the sequence codes for the particular gene it seems to be pointing to.

I believe GeneSCF v.1.0 is able to use custom databases. Take a look to see if it will work in your case.

ADD COMMENTlink modified 15 months ago • written 15 months ago by GenoMax94k

It worked like a bless! Thank you!

ADD REPLYlink written 15 months ago by Light9250
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1012 users visited in the last hour