Question: High number of reads with the incorrect strand designation.
0
gravatar for Ati
9 months ago by
Ati30
Ati30 wrote:

I have used CollectRnaSeqMetrics to check the numbers of reads with the correct/incorrect strand designation. The results showed that more than 80% of reads in almost all of the samples are mapped to the incorrect strand!

What would be a good explanation for these results?

Thank you in advance!

ADD COMMENTlink modified 9 months ago • written 9 months ago by Ati30

Every RNAseq kit out there does not capture the reverse strand. There are some that capture coding strand as well. Have you checked into what kind of kit was used?

ADD REPLYlink written 9 months ago by genomax91k

I'm not sure about the exact kit but it was a stranded RNA-seq library prep kit.

ADD REPLYlink modified 9 months ago • written 9 months ago by Ati30

Make sure to set the right flags and use the tool correctly, depending on the library prep either the second in pair or the first in pair will indicate the correct strand. So you can easily get the opposite result.

Do visualize your data as well.

ADD REPLYlink written 9 months ago by Istvan Albert ♦♦ 85k

Thanks! As my data is paired I have specified STRAND=SECOND_READ_TRANSCRIPTION_STRAND. That's why I'm a bit surprise with the results!

ADD REPLYlink written 9 months ago by Ati30
2

That setting has nothing do with the data being paired-end. You need to identify which strand was captured by library prep and got sequenced.

See this image: A: Tophat Library-Type : Illumina Truseq Stranded Total Rna Sample Prep Kit

ADD REPLYlink modified 9 months ago • written 9 months ago by genomax91k
1

Plus the nomenclature is non standardized and quite confusing. Different tools may designate the same protocol with different names.

For example, the TrueSeq Illumina protocol when analyzed with TopHat you need to specify fr-firststrand designation, but as a matter of fact, in the read pair, the second in pair will match in the original orientation.

ADD REPLYlink modified 9 months ago • written 9 months ago by Istvan Albert ♦♦ 85k

That's now clear to me! Thank you!

ADD REPLYlink written 9 months ago by Ati30
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: 793 users visited in the last hour