Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
cuffdiff for differential gene testing Paul Walker RNA Sequencing 4 01-26-2015 08:04 AM
Cufflinks (Cuffdiff) v0.9.3 multiple samples keebs42 Bioinformatics 4 11-06-2013 02:20 AM
New differential testing of cuffdiff/cufflinks since 1.3.0 DerSeb Bioinformatics 46 09-12-2013 05:46 AM
multiple testing correction question gene_x Bioinformatics 4 11-08-2012 11:51 PM
edgeR multiple testing colindaven RNA Sequencing 2 04-26-2012 05:04 AM

Thread Tools
Old 02-26-2013, 08:37 AM   #1
Location: Freiburg

Join Date: Oct 2012
Posts: 56
Default Cuffdiff - 1vs1 or multiple testing?

Dear all,

I am analyzing my RNA-seq data from mice that are isogenic, but yet exhibit different phenotypes. I am interesting in getting DEG and have therefore, in parallel with DESeq, analyzed my data with Cuffdiff (no gene/transcript discovery mode).

In essence, I have a control wildtype sample (WT) and then I have samples from two mice that are isogenic but different phenotypes (Pheno01 and Pheno02).

Firstly, I am interested in knowing what genes are DE between:
WT vs Pheno01
WT vs Pheno02

but more interestingly, what genes are DE between the two isogenic mice:
Pheno01 vs Pheno02

My question then goes on if it would be most correct to run Cuffdiff three times, or just run Cuffdiff one time but then include Wt, Pheno01 and Pheno02 and then do the pairwise comparison using cummeRbund.
The reason for asking is that it gives completely different results regarding the amount of DEG.

The same holds true for a colleagues experiment where he has progenitor cells which he has differentiated in to five different lineages, where he wants to compare them with each other to see how different/similar they are.

Thanks a lot! I appreciate your help.
DonDolowy is offline   Reply With Quote
Old 02-26-2013, 05:15 PM   #2
Location: china

Join Date: Dec 2011
Posts: 48

I encountered the same proplem, for multicondition's comparision, I first use cuffdiff to do reads counting and normalizition (with fpkm output) by --nondiff option, then extracted those expression information and test them manually(test method : anova or GLM). hope it helps!
zinky is offline   Reply With Quote
Old 02-26-2013, 10:15 PM   #3
Senior Member
Location: Vienna

Join Date: Mar 2010
Posts: 107
Default you need biological replicates!

my opinion is that you can't (by no ways) get a meaningful DE list without biological replicates.

you never know if your DE genes are (really) differentially expressed (from your conditions) or are only fluctuating...

if you would have, let's say, at least three replicates from each type you could use programs which can deal with your design, like edgeR (w/ limma) or limma, using contrasts.

cuffdiff is always the worst choice...
dietmar13 is offline   Reply With Quote
Old 02-26-2013, 11:11 PM   #4
Location: china

Join Date: Dec 2011
Posts: 48

well , I mean , I just employed cuffdiff to do reads counting as well as normalizition(fpkm , the most classic method), and of course there are replicates on my hand. DonDolowy's question is not about wether the DElist is meanful but on how to avoid three times cufdiff running , actually, you can put all your data into one command (with " " delimed for each condition ,ps: you must have replicates), cuffdiff can do comparision of each two condition for you.
zinky is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:44 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO