SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: Comparative Analysis of RNA-Seq Alignment Algorithms and the RNA-Seq Unified Newsbot! Literature Watch 3 07-31-2011 07:08 PM

Reply
 
Thread Tools
Old 08-29-2016, 10:56 PM   #1
mike.t
Member
 
Location: Spain

Join Date: Mar 2010
Posts: 36
Default PE vs SE for RNA-Seq

I have used paired end sequences for RNA-Seq experiments in the past with the understanding that PE reads are 'more accurate' somehow. One paper that I found says that the false positive rate is lower with PE. Single reads are a lot less expensive and I may be able to afford to sequence the libraries deeper. What do others on this forum use? PE or SE for rna-seq?
mike.t is offline   Reply With Quote
Old 08-30-2016, 12:52 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

For normal differential expression single-end 50 base reads are sufficient. We tend to use paired-end reads whenever we need to look at splicing, though, since inevitably that works better.
dpryan is offline   Reply With Quote
Old 08-30-2016, 04:38 AM   #3
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 179
Default

Actually, if the organism you work with is danio rerio (zebrafish), it may help to use paired-end even for single DE (a lot of duplications in this genome). For mouse or Humans, single-end are fine...
SylvainL is offline   Reply With Quote
Old 09-16-2016, 05:20 AM   #4
prathima
Junior Member
 
Location: Bangalore

Join Date: Aug 2016
Posts: 6
Default

Paired-end sequencing facilitates detection of genomic rearrangements and repetitive sequence elements, as well as gene fusions and novel transcripts.
prathima is offline   Reply With Quote
Old 09-17-2016, 10:15 AM   #5
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

I would always recommend PE for anything other than quantification studies. For quantification, it depends. But otherwise - per base, PE should always be cheaper than SE, and give massive advantages for accuracy.

This:
Quote:
Single reads are a lot less expensive and I may be able to afford to sequence the libraries deeper.
Should, in the general case, be false. I guess it depends on how you define "deeper", though. Normally I think of it in terms of fold coverage rather than the total quantity of unique fragments.

So, it depends on the goal of your RNA-seq experiment! If you want to identify isoforms, or want maximal specificity, absolutely you should do PE. If you are just interested in differential expression of genes in general and don't care which isoforms are expressed, then SE will be more cost-effective.
Brian Bushnell is offline   Reply With Quote
Old 09-19-2016, 10:12 PM   #6
MikhailFokin
Member
 
Location: NZ

Join Date: Mar 2014
Posts: 15
Default

Sorry, may be a bit off-topic. I personally use PE libraries, but just recently found the issue with read-through (when forward and reverse reads overlap completely). The question is - to drop reverse ones or leave them? Could it bias coverage and downstream DE analysis? If drop, then how to deal with less overlapped pairs? Ideally I would like to preserve non-overlapped fragments only, but not sure if there is a software capable for...
MikhailFokin is offline   Reply With Quote
Old 09-20-2016, 11:23 AM   #7
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

If you preserve non-overlapping fragments only, you will get a bias against shorter transcripts. So I don't suggest it. Though if you really want to, you can do so by adapter-trimming and setting the minimum length to the initial read length (filtering out anything with adapter sequence).
Brian Bushnell 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 07:50 AM.


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