SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
50 bp paired end reads vs. 100 bp single end reads efoss Bioinformatics 12 01-15-2014 08:05 PM
cuffdiff single end reads madsaan Bioinformatics 7 08-06-2013 08:32 PM
TopHat -paired end vs single end reads adarshjose RNA Sequencing 10 06-12-2012 06:15 PM
Can paired-end mapping produce more reads than single-end ? warrenemmett Bioinformatics 13 03-20-2012 11:10 PM
CuffLinks/CuffDiff frag length for single-end reads mbblack Bioinformatics 2 12-17-2011 01:47 PM

Reply
 
Thread Tools
Old 11-14-2011, 05:53 PM   #1
zun
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 26
Default Can Cuffdiff treat paired-end and single-end reads at the same time?

Hi there,

I run cuffdiff with two paired-end reads samples and 2 single-end reads samples at the same time the other day(after running tophat toward each samples), and got results which seemed normal and correct with no error messages.

As I read cufflinks, cuffcompare, cuffdiff manuals and didn't find any statements that say they can't treat single-end and paired-end samples together, so I think my results are reliable(I don't know how accurate they are, though).

But I just want to know other RNA-seq users' opinions over cuffdiff's treatment with single/paired-end reads together.

Have anyone done the same analysis before?
I would be very happy if anyone would reply to me and give some advises!

zun
zun is offline   Reply With Quote
Old 01-22-2012, 11:38 PM   #2
zun
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 26
Default

Hi there,
I am now wondering how I can treat pair-end reads as two single-end read, and integrate them together when doing cufflinks.

I got a low mapping rate(40%) after running tophat to my pair-end reads, but when I checked log files(ex. logs/bowtie.left_kept_reads.fixmap.log), the rate was more than 75%!! I think it is because this high rate was just a result from Bowtie, after screening the left and right reads which are not mapped on the same chromosome and disposing these.

Then, I come to think whether I can separate pair-end reads to single-end ones , run tophat for each of them and finally run cufflinks for them together....

What do you think of this idea?
Does anyone have other idea on this ?
zun is offline   Reply With Quote
Old 02-20-2012, 01:14 AM   #3
zorph
Member
 
Location: FL

Join Date: May 2010
Posts: 40
Default

hello out there,
I too was also wondering whether the cufflinks calculations would be correct if I merged a single read data set with a paired-end read data set and ran cufflinks on the merge file?

Would everything be ok? or would the numbers be off?

Thanks!
zorph is offline   Reply With Quote
Old 06-12-2012, 05:37 PM   #4
zun
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 26
Default

Hello again,

Now Tophat2 was released and has a new option called "--report-discordant-pair-alignments", which allows mate pairs to map to different chromosomes.
So my problems was solved described above, in fact my mapping rate improved up to 70% from 40%! thanks!

BUT, on the other hand, when it comes to FPKM calculation, how are the mate pairs which are mapped on different chromosomes treated in Cufflinks??
Is there any new calibration to it? or these are treated as single read?
I don't know whether I can trust the FPKM resulted from this option....

I will be appreciated any help!

zun
zun is offline   Reply With Quote
Reply

Tags
cuffdiff, cufflinks, differential, paired, single

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:37 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO