Question: RNA-seq Infer experiment wrong strandeness information
0
gravatar for abedkurdi10
14 days ago by
abedkurdi1030
abedkurdi1030 wrote:

Hello everyone,

I have received recently eight RNA-seq samples. We know that the reads are reversely stranded. After aligning those samples to their reference genome, I run infer_experiment on the alignment data just to make sure everything is fine.

Five samples out of eight showed that they are really reversely stranded, however, the remaining three samples showed that the reads are forward stranded. In addition, those three samples show that they contain high percentage of rRNA sequences (using bbduk ~30 percent of the reads).

Removing those reads did not change the result of infer_experiment.

So what could be the reason(s) behind such confusion? Thank you in advance.

ADD COMMENTlink modified 14 days ago • written 14 days ago by abedkurdi1030

This is weird. Can you extract the reads mapped to forward strand and see what are they? Maybe you have a lot of anti-sense RNA? (long shot, I know). Also, what organism?

ADD REPLYlink written 14 days ago by Asaf8.2k

I will have a look on them. Actually, It's mouse data.

ADD REPLYlink written 14 days ago by abedkurdi1030
0
gravatar for abedkurdi10
14 days ago by
abedkurdi1030
abedkurdi1030 wrote:

It seems that there is more rRNA sequences. I downloaded rmsk bed file from UCSC and removed all reads covering those Repeating Elements (including and not limited to rRNA, tRNA...), then I run again infer_experiment and everything went well.

ADD COMMENTlink written 14 days ago by abedkurdi1030

Why did you delete and close the post? It is relevant to the forum, please revert it

ADD REPLYlink written 13 days ago by Asaf8.2k
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: 1891 users visited in the last hour