Question: Cufflink Workflow
gravatar for Ianiri, Giuseppe
5.6 years ago by
Ianiri, Giuseppe90 wrote:
Hi, I want to compare gene deferentially expressed in two conditions (1 and 2). I don't have a GTF reference annotation. In my first attempt, I mapped my Illumina reads with TopHat. Then I run Cufflink, and Cuffmerge on my Cufflink outputs for the conditions 1 and 2. At the end, I run Cuffdiff using the Cuffmerge output (1 + 2) and the TopHat output for the two samples of interest. I got genes and transcripts differentially expressed, but I was expecting a quite different output (based on the literature). Do you guys think I did anything wrong? Also, I am not sure the question makes sense, but is the Cuffmerge step essential for the comparison of my two samples? Or could I use my Cufflink outputs for samples 1 and 2 for running Cuffdiff? Regards, Giuseppe Ianiri Email:
rna-seq cuffmerge • 1.8k views
ADD COMMENTlink modified 5.6 years ago by Jennifer Hillman Jackson25k • written 5.6 years ago by Ianiri, Giuseppe90
gravatar for Jennifer Hillman Jackson
5.6 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hi Giuseppe, Your workflow sounds correct. Cuffmerge or Cuffcompare is needed and no, you do not want to use Cufflinks as sample input - these should be mapped reads as you state that you have used them so far. Are you using the public Main Galaxy server at If so and you would like to share a history, I can take a quick look to see if any usage problems are present. For the overall results - there are many factors that can contribute to different results vs what is published. It may be that parameters for some of these steps need to be tuned, or that obtaining the same results without a reference annotation to guide assembly is not possible. We can offer help torule out anything technically wrong with the workflow, so you can concentrate on these scientific questions (test/try alternate parameter settings and such). The paper from the tool authors (in the link below) would be a good place to start once you get to this point, if you are not already using it. To share, use the gear icon at the top of the history panel (upper right corner), select "Share or Publish", click on the first button to generate the share link, copy and paste that into a reply email to just me (not the entire list, to keep your data private). If instead you want to just explore this more on your own, links to help for RNA-seq tools are in this wiki: See Tools on the Main server: Example ? RNA-seq analysis tools. Best, Jen Galaxy team -- Jennifer Hillman-Jackson Galaxy Support and Training
ADD COMMENTlink written 5.6 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.


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