SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
comparing results by cuffdiff, edgeR, DESeq PFS Bioinformatics 5 03-12-2014 03:01 AM
Comparing output from Bowtie and BWA maasha Bioinformatics 15 10-25-2012 05:56 AM
Questions about Top-Hat and Paired End Reads Cirno Bioinformatics 4 12-03-2011 11:07 AM
DESeq results give extremely small p-values? chris Bioinformatics 11 08-29-2011 06:33 AM
comparing results from two different reference genomes BAJ Bioinformatics 2 02-24-2009 07:38 AM

Reply
 
Thread Tools
Old 03-12-2012, 08:04 AM   #1
maryb
Junior Member
 
Location: CA

Join Date: Jan 2012
Posts: 6
Default comparing Bowtie/DESeq and Top-Hat/Cufflinks results

Dear all, I have 2 RNA-seq libraries (40 bp single end) and a genome annotation, I am interested in differential expression.

I run:
1- Bowtie/DESeq at the gene level
2-TopHat/Cufflinks at the transcript level.

I got very different results (not only in terms of quantification, but also "direction" of changes)-see example below. I was expecting differences, but not this much.

Which method do you think best suites the type of data I have?
Is is appropriate to try to run TopHat with 40bp-single end reads?
The mean N. of reads given by DESeq does not account for transcript length, would this prevent comparison of transcript quantification levels within a library?

thanks in advance for any reply,


Bowtie
Raw N. reads
gene Transcript Transcript length conditionA conditionB
1 1 1590 297 242
2 2 198 0 0
3 3 2048 383 500
4 4 2034 283 109
5 5-a 788 86 137
5 5-b 1268
6 6 2087 303 640
7 7 1656 0 0
8 8 1809 316 335
9 9 761 0 0
10 10-a 735 658 386
10 10-b 524
TopHat-Cufflinks
FPKM-A FPKM-B
gene Transcript Transcript length conditionA conditionB ln(fold_ch) AvB
1 1 1590 20.526 11.7229 0.560149
2 2 198 45.8285 0 1.79769e+308
3 3 2048 17.5533 9.35482 0.62935
4 4 2034 28.2751 9.71151 1.06867
5 5-a 788 6.67631 1.6504 1.39755
5 5-b 1268 32.5 4.01143 2.09209
6 6 2087 53.4758 3.36856 2.76474
7 7 1656 0.110199 0 1.79769e+308
8 8 1809 16.365 15.5165 0.0532368
9 9 761 2.85777 0 1.79769e+308
10 10-a 735 6.11169 3.07078 0.688272
10 10-b 524 818.778 1315.66 -0.474281
Bowtie-DESeq
MeanReadsA MeanReadsB
gene Transcript Transcript length conditionA conditionB Log2FCAvB
1 1 1590 284.9435568 252.2394288 -0.175228351
2 2 198 0 0 0
3 3 2048 367.4524655 521.1558445 0.503001955
4 4 2034 271.511874 113.6119741 -1.249561315
5 5-a 788 82.50890871 142.7967014 0.784028564
5 5-b 1268
6 6 2087 290.6999923 667.079481 1.195534401
7 7 1656 0 0 0
8 8 1809 303.1722692 349.1744158 0.203185051
9 9 761 0 0 0
10 10-a 735 631.2890922 402.332312 -0.648615343
10 10-b 524
maryb is offline   Reply With Quote
Old 03-12-2012, 09:23 AM   #2
kopi-o
Senior Member
 
Location: Stockholm, Sweden

Join Date: Feb 2008
Posts: 319
Default

I don't get why you are using Bowtie here for the gene level analysis. Why not calculate read counts for the TopHat output and feed into DESeq? With Bowtie, you will miss splice junction spanning alignments. This is one reason why the two cases may not be comparable.
kopi-o is offline   Reply With Quote
Old 03-12-2012, 09:25 AM   #3
maryb
Junior Member
 
Location: CA

Join Date: Jan 2012
Posts: 6
Default

thank you Kopi-o for your reply. I have 40 bp single end reads, i used Bowtie because i ma not sure that this type of read will be ok with TopHat. will they? thanks
maryb is offline   Reply With Quote
Old 03-12-2012, 09:30 AM   #4
kopi-o
Senior Member
 
Location: Stockholm, Sweden

Join Date: Feb 2008
Posts: 319
Default

I think that should be OK. To be honest I haven't run TopHat on shorter reads than 1x50 bp, which worked fine.
kopi-o is offline   Reply With Quote
Old 03-13-2012, 04:01 AM   #5
pbluescript
Senior Member
 
Location: Boston

Join Date: Nov 2009
Posts: 224
Default

Quote:
Originally Posted by maryb View Post
thank you Kopi-o for your reply. I have 40 bp single end reads, i used Bowtie because i ma not sure that this type of read will be ok with TopHat. will they? thanks
You can use Tophat with these reads, but be sure to change the --segment-length option to 20. With a default of 25, you'll have issues with reads that short.

If you want an opinion on which method for DE analysis is better, you should settle in for a long read on this board. People debate it frequently. Personally, I see no clear winner that will work for every situation and generally end up using Cufflinks, DESeq, and RSEM. Whatever pops up consistently is what I take on for further analysis.
pbluescript is offline   Reply With Quote
Old 03-13-2012, 07:46 AM   #6
maryb
Junior Member
 
Location: CA

Join Date: Jan 2012
Posts: 6
Default

thank you very much
maryb is offline   Reply With Quote
Reply

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 12:11 PM.


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