This is a test version of Biostars. For the public version, visit https://www.biostars.org.
From cuffmerge to cuffdiff: What input for cuffdiff?

I want to compare 2 conditions: controls (4 biological replicates) v treated (3 biological replicates). I have run cufflinks successfully on all 7 BAM files independently,

I then ran cuffmerge twice to produce 2 transctipt.gtf files:

1) merged_control_transcrpts.gtf

2) merged_treated_transcripts.gtf

I was expecting that the two transctipts.gtf files would then be used as input (with BAM files) for cuffdiff. But, cuffdiff only allows 1 transctipts.gtf. Where have I gone wrong? Am I supposed to merge everything? Can cuffdiff tell the treated group from the control group?

Thanks in advance. Kenneth

rna-seq cuffmerge cuffdiff treated controls
• 3,075 views • link
updated by Biostar Biostar • written by BioinfGuru BioinfGuru

1 answer

You need to merge all your cufflinks assemblies into a single gtf file using cuffmerge. Both conditions will be contained in a single file. You can then use this within cuffdiff and your original BAM files to test for your conditions.

• 1 views • link
written by onspotproductions onspotproductions

You sir....are a legend.

• 1 views • link
written by BioinfGuru BioinfGuru

Log in to answer this question.

More posts like this