SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
how does stranded and un-stranded quantification work? weasteam Bioinformatics 2 05-10-2013 10:15 AM
tool for checking the stranded vs unstranded library MichalO Illumina/Solexa 5 04-29-2013 08:49 AM
Mapping paired-end stranded RNA-seq data hubin.keio Bioinformatics 3 03-28-2013 02:10 PM
Bias adjustment in cufflinks: stranded and unstranded protocol question... adeonari Bioinformatics 0 03-21-2013 03:53 AM
stranded or non-stranded protocol? enelkinsan RNA Sequencing 7 12-18-2012 12:42 PM

Reply
 
Thread Tools
Old 02-27-2014, 06:36 PM   #1
jake13
Member
 
Location: San Francisco

Join Date: Aug 2013
Posts: 23
Default Cuffdiff stranded and unstranded data

I have generated stranded RNA-seq data and have been comparing it using cuffdiff. However, I now want to compare some of my samples to a previously published data set that does not have stranded data. Can I compare these samples in cuffdiff? If so, what do I specify for the library type?
Thanks
jake13 is offline   Reply With Quote
Old 02-27-2014, 07:22 PM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

No, you can't really do a valid comparison unless the organism does not have overlapping opposite-strand genes. Fortunately... outside of prokaryotes, organisms generally don't have overlapping opposite-strand genes. So if you're working with multicellular eukaryotes, just declare both libraries as unstranded and compare the results directly. Try evaluating the same library against itself, processing it stranded versus unstranded. If the results indicate no difference, you might be OK. The different protocols may have different biases, though. I would consider the data as informative, but I would certainly not publish anything based on comparing datasets with different methodologies.

Last edited by Brian Bushnell; 02-27-2014 at 08:08 PM.
Brian Bushnell is offline   Reply With Quote
Old 03-03-2014, 12:05 AM   #3
wenhua
Junior Member
 
Location: NanJing

Join Date: Oct 2013
Posts: 1
Default cuffdiff for DEGs

I use tophat ,cufflinks and cuffdiff analysis of RNA-seq 2 samples of three biological replicates .
I want to know why the value of log2(fold_change) between two samples is relatively large(such as 6) and q>0.05 resulting in significant value is "no".
How to use "log2(fold_change)"and "q value" in screening differentially expressed genes between the two samples.

Another case is that one gene FPKM value in sample1 is 0 and,in sample2 is greater than 0 (such as 14), the result from cuffdiff also no difference between two samples.Why?
Anyone else notice this?

Thanks
wenhua is offline   Reply With Quote
Old 03-03-2014, 06:32 AM   #4
jwfoley
Senior Member
 
Location: Stanford

Join Date: Jun 2009
Posts: 179
Default

You could just accumulate per-transcript read counts and use a more flexible program like DESeq2.

Of course there is the risk of ambiguity for overlapping genes, as Brian Bushnell points out (and partial overlap does happen frequently in multicellular organisms), but if you're worried about that then you have to discard your unstranded data altogether.
jwfoley 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 06:14 AM.


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