SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Collaborator for 3C-related methodology analysis Sarah_Tulin General 0 08-21-2013 12:09 PM
Methodology used by Fastx clipper babi2305 Bioinformatics 0 02-13-2013 11:50 AM
In Sequence: ABI, Illumina, Former Lawyer Fail to Settle IP-Ownership Suit; ABI Seeks Newsbot! Illumina/Solexa 0 04-01-2008 01:14 PM
In Sequence: Court ‘Encourages’ ABI to Settle Suit Against Former Counsel, Illumina Newsbot! Illumina/Solexa 0 03-11-2008 03:04 PM

Reply
 
Thread Tools
Old 03-10-2015, 05:26 AM   #1
lre1234
Senior Member
 
Location: US

Join Date: Aug 2011
Posts: 105
Default help settle a methodology discussion.

Hi all,

My boss and I are having a debate on a protocol that I need a little clarity on, and hopefully someone can give the best answer.

The problem: we are interested in using RNA-seq data and looking at strand specific expression data. But, all of the sequence data that we have is from paired-end RNA sequence.

My boss thinks that this isn't a problem and we can do one of two things. either a) only use the forward sequence reads and look at the data from that, or b) simply take the reverse reads, and use the reverse-compliment of them and treat them as forward reads.

My answer is that he is wrong in both cases. By doing paired-end sequencing, you inherently loose the strand information and cannot do any strand specific information. Simply reversing the two reads will not work, and you do not know which strand the forward (or reverse read came from). It is not always the case that the forward read is the plus strand and the reverse read is the minus strand, and these can be swapped. Essentially, we do not know where (strand and orientation) that the read came from.

Can someone help confirm that I am correct .
lre1234 is offline   Reply With Quote
Old 03-10-2015, 05:45 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

You're both wrong actually. Whether strand information is conveyed is determined by the nature of the library prep, not whether you sequence single-end or paired-end. You can certainly do strand-specific paired-end sequencing, that's not even all that uncommon. So, the question is purely whether you have a strand-specific dataset. If so, then just determine which of the reads is conferring the strand information (normally read #2 does, since dUTP methods are the most popular) and proceed accordingly.
dpryan is offline   Reply With Quote
Old 03-10-2015, 08:46 AM   #3
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 311
Default

Quote:
Originally Posted by dpryan View Post
You're both wrong actually. Whether strand information is conveyed is determined by the nature of the library prep, not whether you sequence single-end or paired-end. You can certainly do strand-specific paired-end sequencing, that's not even all that uncommon. So, the question is purely whether you have a strand-specific dataset. If so, then just determine which of the reads is conferring the strand information (normally read #2 does, since dUTP methods are the most popular) and proceed accordingly.
As a side comment... Paired end sequencing is such a source of confusion. Sometimes I wish it was never invented!
dariober is offline   Reply With Quote
Old 03-10-2015, 10:23 AM   #4
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,315
Default

Quote:
Originally Posted by dariober View Post
As a side comment... Paired end sequencing is such a source of confusion. Sometimes I wish it was never invented!
Do you also wish you didn't have a nose every time you catch a cold?

--
Phillip
pmiguel is offline   Reply With Quote
Old 03-10-2015, 10:27 AM   #5
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,315
Default

Quote:
Originally Posted by dpryan View Post
You're both wrong actually. Whether strand information is conveyed is determined by the nature of the library prep, not whether you sequence single-end or paired-end. You can certainly do strand-specific paired-end sequencing, that's not even all that uncommon. So, the question is purely whether you have a strand-specific dataset. If so, then just determine which of the reads is conferring the strand information (normally read #2 does, since dUTP methods are the most popular) and proceed accordingly.
Completely agree.

In addition I don't think informaticians should trust that data sets are of a particular strandedness. They should test data sets to determine that. Just align 1000 reads to the large rRNA subunit sequence for the species in question. If you have paired reads this will also answer the question "what insert size was sequenced".

--
Phillip
pmiguel is offline   Reply With Quote
Old 03-10-2015, 10:56 AM   #6
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 311
Default

Quote:
Originally Posted by pmiguel View Post
Do you also wish you didn't have a nose every time you catch a cold?
In case it wasn't obvious... My comment was ironic. I just wanted to point out that PE can make things confusing. (And if I really hated PE I would just use only the 1st-in-pair reads).
dariober is offline   Reply With Quote
Old 03-10-2015, 11:05 AM   #7
pmiguel
Senior Member
 
Location: Purdue University, West Lafayette, Indiana

Join Date: Aug 2008
Posts: 2,315
Default

Quote:
Originally Posted by dariober View Post
In case it wasn't obvious... My comment was ironic. I just wanted to point out that PE can make things confusing. (And if I really hated PE I would just use only the 1st-in-pair reads).
Hey, I'm not here to judge...

--
Phillip
pmiguel is offline   Reply With Quote
Old 03-10-2015, 11:11 AM   #8
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 311
Default

Quote:
Originally Posted by pmiguel View Post
Hey, I'm not here to judge...

--
Phillip
Sure, no worries!
dariober is offline   Reply With Quote
Reply

Tags
paired-end reads, rna-seq advice, strand direction

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:27 PM.


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