Question: (No Subject)
gravatar for Davide Degli Esposti
5.7 years ago by
Davide Degli Esposti80 wrote:
Dear Galaxy team, I have a question about RNA analysis with the cufflinks package. I have some bam files to analyze from a SOLiD platform. Some previous tests show that these bam/sam files are different from those coming from Tophat and cufflinks cannot assemble them using a reference annotation (XS attribute lacking in spliced alignments). (see datasetscufflinks). An apparent solution is to include the reference annotation in the cuffmerge (see degliesposti/h/rna-seqtest-datasetsapril-20132) or cuffcompare (see datasetsjan-2013-1) steps. Doing like this allowed me to run cuffdiff on my datasets without apparent technical errors. However, when I compare the list of differentially expressed transcripts (DETs), these results extremely different: using cuffcompare, I got 390 DETs and using cuffmerge I got 770 DETs, but just 60 genes are shared between the two lists. The parameters used in cuffdiff (FDR, Min Alignement counts, etc.) are the same for the two analyses. Do you have any explanation about that? I expected that cuffcompare and cuffmerge did not lead to outputs quantitatively different. Where may the source of this difference be? I thank you for your cooperation Davide Davide Degli Esposti, PhD Epigenetic (EGE) Group International Agency for Research on Cancer Tel. +33 4 72738036 Fax. +33 4 72738322 150, cours Albert Thomas 69372 Lyon Cedex 08 France This message and its attachments are strictly confidential. If you are not the intended recipient of this message, please immediately notify the sender and delete it. Since its integrity cannot be guaranteed, its content cannot involve the sender's responsibility. Any misuse, any disclosure or publication of its content, either whole or partial, is prohibited, exception made of formally approved use
rna-seq cuffmerge cufflinks • 1.0k views
ADD COMMENTlink modified 5.7 years ago by Jeremy Goecks2.2k • written 5.7 years ago by Davide Degli Esposti80
gravatar for Jeremy Goecks
5.7 years ago by
Jeremy Goecks2.2k
Jeremy Goecks2.2k wrote:
Cuffmerge does some additional steps that Cuffcompare does not; specifically, Cuffmerge attempts to remove assembly artifacts: It's likely that the (presumed) artifacts removed by Cuffmerge account for the differences that you're seeing. Best, J.
ADD COMMENTlink written 5.7 years ago by Jeremy Goecks2.2k
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: 172 users visited in the last hour