GO and KEGG enrichment analysis based on cuffdiff

Roger Z · 10-27-2014, 06:25 PM

Hi,everyone!I used cufflinks and cuffdiff to detect differentially expressed genes, and I am going to do GO and KEGG pathway enrichment analysis based on the result. The problem is that some genes of cuffdiff output have no corresponding GeneID(they are assembled by cufflinks), which makes it impossible to do GO enrichment analysis. Could you please give me some advice? How can we solve this problem?

Best regards,

Parharn · 10-28-2014, 02:08 AM

I do DESeq2, have you tried that?
And what is your organism?

kiwistew · 11-13-2014, 03:02 PM

Option one - do this part of the analysis without the novel transcripts. Just do the ontology analysis with expression data for known transcripts.

Option two - try and assign gene ontologies to the novel transcripts. Then included them in the analysis.

Option two would clearly be a lot more work but might be neccesary if you're working on a poorly annotated organism.

Roger Z · 11-17-2014, 10:34 PM

Thank you, kiwistew. I have chosen option 2 and I used blast2GO.

Roger Z · 11-17-2014, 10:35 PM

Thank you,Parharn! Organism is zebrafish and the problem has been solved using cuffdiff.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Gene enrichment on cuffdiff output	Parharn	Bioinformatics	3	10-28-2014 02:05 AM
Mapping transcripts to Kegg gene IDs / Using Paintomics on a non-Kegg species	anth	Bioinformatics	0	06-27-2014 03:26 PM
Enrichment analysis based on genetic location	LouiseT	Bioinformatics	1	03-04-2013 07:37 PM
Editing a Kegg Ortholgy file for metagenomic analysis in Qiime	Giorgio C	Bioinformatics	0	11-19-2012 12:49 PM
Pubmed: Microarray-based multicycle-enrichment	read2009	Literature Watch	0	09-27-2009 01:20 PM

10-27-2014, 06:25 PM	#1
Roger Z Junior Member Location: China Join Date: Oct 2014 Posts: 8	GO and KEGG enrichment analysis based on cuffdiff Hi,everyone!I used cufflinks and cuffdiff to detect differentially expressed genes, and I am going to do GO and KEGG pathway enrichment analysis based on the result. The problem is that some genes of cuffdiff output have no corresponding GeneID(they are assembled by cufflinks), which makes it impossible to do GO enrichment analysis. Could you please give me some advice? How can we solve this problem? Best regards,

10-28-2014, 02:08 AM	#2
Parharn Member Location: Sweden Join Date: Jul 2013 Posts: 84	I do DESeq2, have you tried that? And what is your organism?

11-13-2014, 03:02 PM	#3
kiwistew Junior Member Location: New Zealand Join Date: Nov 2014 Posts: 2	Option one - do this part of the analysis without the novel transcripts. Just do the ontology analysis with expression data for known transcripts. Option two - try and assign gene ontologies to the novel transcripts. Then included them in the analysis. Option two would clearly be a lot more work but might be neccesary if you're working on a poorly annotated organism.

11-17-2014, 10:34 PM	#4
Roger Z Junior Member Location: China Join Date: Oct 2014 Posts: 8	Thank you, kiwistew. I have chosen option 2 and I used blast2GO.

11-17-2014, 10:35 PM	#5
Roger Z Junior Member Location: China Join Date: Oct 2014 Posts: 8	Thank you,Parharn! Organism is zebrafish and the problem has been solved using cuffdiff.