Question: PE or SE reads for alternative splicing / isoform quantification?
gravatar for nevev
2.7 years ago by
nevev90 wrote:

Hi All,

I've had a long break, and last time I learnt about it, I thought PE reads are required for detection and quantification of differential isoform expression. A few days ago someone was trying to convince me that long SE reads (100-150bp) can do as well. Any comments/ideas/opinions?

I am planning an RNA-seq experiment now, and need to decide if I want to pay 2x for PE, or can settle for SE. My main interest will be in DE of genes, but if I already do it, I'd rather have a dataset that'll allow me to look at isoform expression, TSS usage etc. too later on (if I recall correctly, CuffDiff could do all that).

Thanks in advance for good advice,


ADD COMMENTlink modified 2.7 years ago by Brian Bushnell16k • written 2.7 years ago by nevev90
gravatar for Brian Bushnell
2.7 years ago by
Walnut Creek, USA
Brian Bushnell16k wrote:

SE reads are absolutely capable of differential isoform quantification. The best example, of course, would be PacBio reads-of-insert, which give you whole-isoform reads with fairly low error rates (the error rate depends on the isoform length). For Illumina reads, you have two choices:

  1. Merge overlapping paired reads, then map them as single-ended, using a splice-aware tool such as BBMap; then, use any tool to identify isoforms.


  2. Map the paired reads (again using a splice-aware aligner) and rely on post-processing to identify isoforms via a pair-aware protocol.

PE is strictly better than SE for isoform identification, whether or not you decide to aim for overlapping reads that you merge prior to quantification. But it PE is not required, particularly with longer reads.

It's important to decide whether you care more about isoforms or genes. For gene expression quantification, more fragments are better; for isoform expression quantification, longer and paired reads are better.

ADD COMMENTlink modified 5 weeks ago by RamRS18k • written 2.7 years ago by Brian Bushnell16k

Thank you for the answer!

ADD REPLYlink written 2.7 years ago by nevev90
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2204 users visited in the last hour