Question: Differentiating between stranded and unstranded RNA-seq data by mapping to genome
0
gravatar for amy.bashir
13 months ago by
amy.bashir80
amy.bashir80 wrote:

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 • 1.1k views
ADD COMMENTlink modified 13 months ago by Friederike2.3k • written 13 months ago by amy.bashir80
0
gravatar for Friederike
13 months ago by
Friederike2.3k
United States
Friederike2.3k wrote:

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.

ADD COMMENTlink written 13 months ago by Friederike2.3k

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

ADD REPLYlink written 13 months ago by amy.bashir80

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).

ADD REPLYlink written 13 months ago by Friederike2.3k

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?

ADD REPLYlink written 13 months ago by amy.bashir80
Please log in to add an answer.

Help
Access

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