SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
stranded or non-stranded protocol? enelkinsan RNA Sequencing 7 12-18-2012 12:42 PM
Scriptseq v2 vs TruSeq Stranded new2seq Sample Prep / Library Generation 0 11-28-2012 11:20 AM
How does new Truseq Stranded Protocol actually work? Is it dUTP based?? seqseq123 Sample Prep / Library Generation 3 10-02-2012 03:52 AM
question about stranded RNA seq StephaniePi83 RNA Sequencing 11 09-02-2012 04:25 PM
Stranded TruSeq TonyBrooks Sample Prep / Library Generation 1 06-23-2012 08:09 AM

Reply
 
Thread Tools
Old 05-08-2013, 07:16 AM   #1
weasteam
Member
 
Location: FL

Join Date: Dec 2008
Posts: 26
Default how does stranded and un-stranded quantification work?

I have one stranded library and it was aligned to genome in both un-stranded (--library-type fr-unstranded) and stranded (--library-type fr-firststrand) methods (with tophat). The gene level counts were generated with easyRNASeq separately for the two different bam files.

The raw counts are different between un-stranded and stranded aligned bam files. For some genes, the reads from the un-stranded are a little bit higher than the stranded counting. I am wondering how does this happen? When easyRNASeq quantify from the stranded bam file, does it only count the reads from the correct strand and discard the reads from the other strand? I am trying to verify this by counting the reading in IGV, but failed.How does stranded and un-stranded quantification work?
weasteam is offline   Reply With Quote
Old 05-10-2013, 12:59 AM   #2
kokonech
Curious Character
 
Location: Berlin, Germany

Join Date: Sep 2010
Posts: 13
Default

The protocol strand-specificity is well-explained in this post:
http://onetipperday.blogspot.de/2012...pe-to-use.html

According TopHat FAQ the library information is used when calculating jucntions:
http://tophat.cbcb.umd.edu/faq.shtml#library_type

I am not aware how the quantification with easyRNASeq works, but in general the for strand specific library the read is counted only if the alignment is concordant with transcript strand according to the protocol rule. For unstranded library the strand of the alignment in transcript does not have any restrictions. Usually gene coutning tools have an option to specify strand-specificity of the protocol. See for example, HTSeq-count.
kokonech is offline   Reply With Quote
Old 05-10-2013, 10:15 AM   #3
weasteam
Member
 
Location: FL

Join Date: Dec 2008
Posts: 26
Default

Quote:
Originally Posted by kokonech View Post
The protocol strand-specificity is well-explained in this post:
http://onetipperday.blogspot.de/2012...pe-to-use.html

According TopHat FAQ the library information is used when calculating jucntions:
http://tophat.cbcb.umd.edu/faq.shtml#library_type

I am not aware how the quantification with easyRNASeq works, but in general the for strand specific library the read is counted only if the alignment is concordant with transcript strand according to the protocol rule. For unstranded library the strand of the alignment in transcript does not have any restrictions. Usually gene coutning tools have an option to specify strand-specificity of the protocol. See for example, HTSeq-count.
Thank you kokonech, I will compare the result outputed from htseq-count and easyrnaseq.
weasteam 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:52 PM.


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