Question: Cuffcompare or Cuffmerge for tissue specificity?
gravatar for colin.kern
3.5 years ago by
United States
colin.kern220 wrote:

I have RNA-seq libraries from 8 tissues (2 biological replicates for each). I want to determine which transcripts have tissue-specific expression. After I build cufflinks annotations for each of the individual tissues, would it be better to combine them into a single annotation using cuffmerge or cuffcompare? I'm concerned that if only one of the tissues has an assembled transcript at a specific loci, cuffmerge may consider it to not have enough evidence to include. But that's exactly the type of thing I'm interested in finding. Am I misunderstanding how cuffmerge works? It seems like cuffcompare might be better in this situation as it just does a straightforward union of all fragments.

rna-seq cufflinks • 885 views
ADD COMMENTlink modified 3.5 years ago by vivekbhr510 • written 3.5 years ago by colin.kern220
gravatar for vivekbhr
3.5 years ago by
vivekbhr510 wrote:

You can safely use cuffmerge for merging. It causes no loss of detected transcripts. 

Also you will find this post useful..

ADD COMMENTlink written 3.5 years ago by vivekbhr510

Thanks for the answer. I tried both and I get 33886 XLOCs with cuffmerge and 51144 with cuffcmp. Also almost twice the number of TCONS. It says in the cuffmerge manual that it "automatically filters a number of transfrags that are probably artfifacts." Could this be what's causing the difference?

ADD REPLYlink written 3.5 years ago by colin.kern220

Did you use the reference annotation (ensembl gtf) with cuffcompare?

ADD REPLYlink written 3.5 years ago by vivekbhr510
Please log in to add an answer.


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