Question: Problem With Repeated Genes In Cuffdiff'S Output
gravatar for Hoang, Thanh
5.4 years ago by
Hoang, Thanh200
Hoang, Thanh200 wrote:
Hi all, I am working on RNA-seq using TopHat/Cufflink/Cuffdiff for differential gene expression and new gene discovery ( this is what I am interested in). However, I found many genes that are repeated in the Cuffdiff's ouput. These are the same genes and at the exact the same locus. There should be only one gene for 1 line. Something like this: *Genes* *Locus* *Status* *q1* *q2* *Log2 Folg change* *Significance* * Lnp* *2:74517521-74584544* *OK* *8.91501* *85.2735* *3.25779* * yes* *Lnp * *2:74517521-74584544* *OK* *12.0044* *171.352* *3.83533* * yes* If I re-run the Cuffdiff for differential gene expression only ( No gene discovery), the problem is fixed. Anyone knows how o explain and fix this? Thank you so much
rna-seq cuffdiff • 1.5k views
ADD COMMENTlink modified 5.4 years ago by Jennifer Hillman Jackson25k • written 5.4 years ago by Hoang, Thanh200
gravatar for Jennifer Hillman Jackson
5.4 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hello, The whole output line is not included, so it is difficult to tell you exactly what is going on in the case you are looking at, but most likely there is a difference below the gene bound level. Look for an alternate transcript start site or coding region - this should be indicated on the line in one of the other columns. (When not doing gene discovery, these alternate TSS/CDS in your data were not utilized, this is the difference between restricting to only known vs including/discovering novel in the input samples). Also, just in case you want to confirm some of these (or visualize them) the first column of each file generally has a test id/identifier that can be used to look up what is being testing in the tracking files. For the actual transcripts, going back to through the data in these files to map into the Cuffmerge/Cufflinks transcripts files, where you can review the details. Pulling up the input BAM with the Cufflinks transcripts files, along with the reference annotation, in trackster may help you confirm/decide if you agree with the results or not (you know the coordinates, so can zoom into the region(s)). The Cuffdiff manual describes each file in detail - this is a great reference. The Galaxy wrapper doesn't alter the output. Hopefully this helps, Jen Galaxy team -- Jennifer Hillman-Jackson Galaxy Support and Training
ADD COMMENTlink written 5.4 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: 168 users visited in the last hour