Does It Makes Sense To Assemble All My Rna-Seq Samples As Though They Were One?
2
2
Entering edit mode
12.3 years ago
Ryan Thompson ★ 3.6k

I am doing reference-guided transcript assembly and quantitifaction using Tophat and Cufflinks. Currently, I am running Cufflinks on each sample individually and then merging the assemblies with the cuffmerge tool. Would I get a better assembly by combining all the BAM files from all my samples into one and running Cufflinks on that? Or should it theoretically produce an identical assembly?

(Obviously the second strategy will yield nonsense FPKM values, but I can always requantify against each BAM file individually, so that's not a problem.)

cufflinks transcriptome transcript assembly • 3.3k views
ADD COMMENT
1
Entering edit mode
12.3 years ago
Bioinfosm ▴ 620

I guess using all samples to make a topHat assembly would not make sense. I have seen a lot of groups combine the assembly from different samples, and then use that as the guide for tophat to get a better handle on transcripts (similar to what you are suggesting). Using data from all samples could actually confuse the assembly and yield inefficient transcript models to use with topHat!

ADD COMMENT
0
Entering edit mode

Even if the samples are from the same biological group ?

ADD REPLY
0
Entering edit mode

Unless it is the same sample, you might be better off to run independent tophat analysis to find all the transcripts separately. Then there is always the option of combining all those to form the inclusive 'transcriptome'

ADD REPLY
0
Entering edit mode
11.4 years ago

Hi,

Little additional question : Is it a good idea if the samples are from the same biological group ?

Thanks

N.

ADD COMMENT

Login before adding your answer.

Traffic: 2855 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6