gene_exp.diff file from different versions of cuffdiff

narges · 09-16-2012, 04:20 AM

Hi all,
Regarding my question I have found another threads but it could not help and I really need your help because my work deadline is approaching soon.
I have an rna-seq dataset (6 samples for each labels) and I should announce differentially expressed genes.
Since my read sequence length is very short (read length = 36), we decided to use tophat version 1 instead of version 2 or the newest one (now I am not sure whether my decision was true or not?)
Then I applied cuffdiff over the bam output file from tophat (I did not run cufflinks since in the manual it is mentioned that there is no need to run it if I only need the differentially expressed genes without any intentions to find new genes). First I used cuffdiff version 2 but I got very few numbers of DEgenes so my college offered to apply cuffdiff version 2 which he believes is less restrict. Below are the results in gene_exp.diff file:

with cufflinks/2.0.2 :
30 DE genes (those which have pvalue <0.01)
status:
fail 116
lowdata 122
notest 10789
ok 12228

with cuflinks cufflinks/1.3.0 and additive parameters: cuffdiff -c 0 -u (I used this version plus these parameters to make cuffdiff less restrict):
74 DE genes
fail 4050
lowdata 104
notest 2618
ok 16513

Now first I want to know if you think it was a correct decision to use tophat version 1 due to the short reads and second do you have any suggestion to which cuffdiff result should I announce or if you have any advice to get more reliable result.
Thanks a lot in advance.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
multiple lines per single gene in gene_exp.diff output of cuffdiff	reut	Bioinformatics	2	03-09-2015 05:16 AM
Cuffdiff splicing.diff	Greg	Bioinformatics	1	11-04-2013 05:04 PM
Cuffdiff NOTEST for splicing.diff	Greg	Bioinformatics	5	11-14-2011 06:22 AM
multiple FPKM problem for single gene in gene_exp.diff after running cuffdiff	ngs	RNA Sequencing	4	03-30-2011 02:55 PM
Strange cuffdiff result (disagreement between tracking and diff files)	jebe	Bioinformatics	0	07-01-2010 11:42 AM

09-16-2012, 04:20 AM	#1
narges Member Location: Finland Join Date: Aug 2012 Posts: 29	gene_exp.diff file from different versions of cuffdiff Hi all, Regarding my question I have found another threads but it could not help and I really need your help because my work deadline is approaching soon. I have an rna-seq dataset (6 samples for each labels) and I should announce differentially expressed genes. Since my read sequence length is very short (read length = 36), we decided to use tophat version 1 instead of version 2 or the newest one (now I am not sure whether my decision was true or not?) Then I applied cuffdiff over the bam output file from tophat (I did not run cufflinks since in the manual it is mentioned that there is no need to run it if I only need the differentially expressed genes without any intentions to find new genes). First I used cuffdiff version 2 but I got very few numbers of DEgenes so my college offered to apply cuffdiff version 2 which he believes is less restrict. Below are the results in gene_exp.diff file: with cufflinks/2.0.2 : 30 DE genes (those which have pvalue <0.01) status: fail 116 lowdata 122 notest 10789 ok 12228 with cuflinks cufflinks/1.3.0 and additive parameters: cuffdiff -c 0 -u (I used this version plus these parameters to make cuffdiff less restrict): 74 DE genes fail 4050 lowdata 104 notest 2618 ok 16513 Now first I want to know if you think it was a correct decision to use tophat version 1 due to the short reads and second do you have any suggestion to which cuffdiff result should I announce or if you have any advice to get more reliable result. Thanks a lot in advance.