This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Differentiating between stranded and unstranded RNA-seq data by mapping to genome

Hello everyone!

I have some paired-end RNA-seq data from a collaborator. I was told that they are "stranded", but nothing about whether it's RF or FR.

I pulled out 5 pairs of reads (read 1 and read 2) from each of the four RNA-seq libraries and aligned them to the reference genome of the organism. I see that 50% of the read 1s map to the "+" strand, the other 50% to the "-" strand. Same as the read 2s.

Does this mean that the data is actually unstranded?

Any help would be greatly appreciated!

rna-seq
• 4,256 views • link
updated by Friederike Friederike • written by amy.bashir amy.bashir

1 answer

I recommend you run QoRTs on your data. This will tell you a lot about the quality of the data including which protocol for the strandedness was used or if the data stems from unstranded library prep.

• 0 views • link
written by Friederike Friederike

Is it not possible to tell from the information that I provided?

• 0 views • link
written by amy.bashir amy.bashir

I cannot tell from the information you provide above, no. And honestly, I also wouldn't want to do it, not based on 5 read pairs. There are established tools for doing that basing their verdict on thousands of reads. If you don't want to use QoRTs, maybe RSeQC tickles your fancy (see this post).

• 0 views • link
written by Friederike Friederike

Great, thank you very much for the clarification. I did try RSeQC, and got the output:

Fraction of reads explained by "1+-,1-+,2++,2--": 0.9803

From what I see around online, this corresponds to a RF library?

• 0 views • link
written by amy.bashir amy.bashir

Log in to answer this question.

More posts like this