Differential Expression Analysis Pipeline with Reference

pcg · 01-27-2011, 07:34 AM

Hi All,

I am wondering what are people's thoughts on the following way to look for DE genes between samples because cuffdiff does not seem to be summing up the FPKM transcripts scores for each gene. SO because I know the reference genome (human) I am using tophat with the -G option and then running cuffdiff as follows with the human genome annotation file downloaded from ensembl.

cuffdiff hg19.gtf accepted_hits_time1 accepted_hits_time 2

Do I need to go through the entire cufflinks/cuffcompare pipeline or will the above suffice

Many thanks in advancce everyone

honey · 01-27-2011, 10:04 AM

I read somewhere in the post if you run TopHat with anotation file it will give one gene one row. I have not tried as yet and has same problem.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
RNA-Seq: Differential expression analysis for sequence count data.	Newsbot!	Literature Watch	3	09-25-2012 10:19 AM
Differential Expression analysis without replicates	polsum	Bioinformatics	1	08-05-2011 04:40 AM
RNA-Seq Quantification and Differential Expression Analysis	days369	RNA Sequencing	2	04-06-2011 01:24 AM
Differential expression analysis workflow in Cufflinks	anna_vt	Bioinformatics	4	12-19-2010 03:04 AM
GobyWeb web application for alignment and differential expression analysis	Fabien Campagne	Core Facilities	0	11-04-2010 03:43 PM

01-27-2011, 07:34 AM	#1
pcg Junior Member Location: USA Join Date: Jan 2010 Posts: 8	Differential Expression Analysis Pipeline with Reference Hi All, I am wondering what are people's thoughts on the following way to look for DE genes between samples because cuffdiff does not seem to be summing up the FPKM transcripts scores for each gene. SO because I know the reference genome (human) I am using tophat with the -G option and then running cuffdiff as follows with the human genome annotation file downloaded from ensembl. cuffdiff hg19.gtf accepted_hits_time1 accepted_hits_time 2 Do I need to go through the entire cufflinks/cuffcompare pipeline or will the above suffice Many thanks in advancce everyone

01-27-2011, 10:04 AM	#2
honey Senior Member Location: Pittsburgh Join Date: Feb 2010 Posts: 151	TopHat with annotation I read somewhere in the post if you run TopHat with anotation file it will give one gene one row. I have not tried as yet and has same problem.