SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: A Strand-Specific Library Preparation Protocol for RNA Sequencing. Newsbot! Literature Watch 0 09-29-2011 07:00 AM
always reads on plus strand more than on minus strand tujchl Bioinformatics 4 04-29-2011 01:08 AM
different number of reads mapped to plus strand and minus strand gfmgfm Bioinformatics 2 02-03-2011 11:26 AM
cuffdiff not strand-specific? burkard Bioinformatics 0 11-15-2010 11:30 AM
FRT-seq: amplification-free, strand-specific transcriptome sequencing severin Literature Watch 0 03-22-2010 05:53 AM

Reply
 
Thread Tools
Old 11-15-2010, 09:08 AM   #1
nkwuji
Member
 
Location: Dublin

Join Date: Mar 2010
Posts: 19
Default cuffdiff for strand-unspecific sequencing

Hi,

I notice in my sequencing result that, the reads are not strand specific. There are lots of read matching the opposite strand of mRNA. There is a "--stranded=no" option in HTseq. But I didnot find any similar option in cuffdiff. May I ask whether there is such option in cuffdiff? And, how cuffdiff deals with reads matching opposite strand or genes with no strand information (e.g. "." for strand information)?

Cheers,
Jun
nkwuji is offline   Reply With Quote
Old 11-17-2010, 09:52 PM   #2
GKM
Member
 
Location: Pasadena, CA

Join Date: May 2009
Posts: 45
Default

Which version of cufflinks are you using?? 0.9 and above definitely include a number of options for library types

And, of course, what is your library type, you didn't mention it...
GKM is offline   Reply With Quote
Old 11-18-2010, 02:57 AM   #3
nkwuji
Member
 
Location: Dublin

Join Date: Mar 2010
Posts: 19
Default

Quote:
Originally Posted by GKM View Post
Which version of cufflinks are you using?? 0.9 and above definitely include a number of options for library types

And, of course, what is your library type, you didn't mention it...
Hi, GKM,

Thx for your reply. My library is Illumina RNA-seq. I just notice that there is an option in cufflinks of "--library-type", and the default setting is "fr-unstranded" for standard Illumina. I think I should use the default setting. But I still find a few assembled transcripts from cufflinks matched the opposite strand (or both strands) of RefSeq transcripts.

/Jun
nkwuji is offline   Reply With Quote
Old 11-18-2010, 03:03 AM   #4
GKM
Member
 
Location: Pasadena, CA

Join Date: May 2009
Posts: 45
Default

"Illumina RNA-Seq" doesn't help much. Is your library stranded or not and if it is, what kind of stranded protocol was used? If not, then your original worry about reads mapping to the opposite strand makes no sense since half of your reads will match the opposite strand anyway.

Cufflinks will assemble anti-sense transcripts even in unstranded data from time to time (usually those will have a slightly different splicing pattern from the annotated ones), those will be given very low abundance estimates though
GKM 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 01:37 PM.


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