This is a test version of Biostars. For the public version, visit https://www.biostars.org.
TruSeq strand-specificity in rsem-calculate-expression

Hello,

I'm trying to figure out the settings for analysing strand-specific RNA-Seq data from Illumina TruSeq chemistry with rsem-calculate-expression

Specifically, regarding the --forward-prob parameter, the manual has this to say:

....Set to 1 for a strand-specific protocol where all (upstream) reads are derived from the forward strand, 0 for a strand-specific protocol where all (upstream) read are derived from the reverse strand, or 0.5 for a non-strand-specific protocol. (Default: 0.5)

Question: based upon my understanding of the molecular biology of the TruSeq strand-specific kit, the strand actually getting sequenced in read 1 of a pair is the cDNA sequence, (or in RSEM manual terminology, the reverse strand).

As such, my hunch is that I should be using --forward-prob=0, but I'd like to check with the community to see if I'm getting things right here.

Any help much appreciated

Thanks,
Shurjo

rna-seq rsem strand-specific
• 7,943 views • link
updated by Ram Ram • written by shurjo.sen • shurjo.sen

3 answers

The Tru-seq strand-specific protocol is specific to the REVERSE STRAND.

Instead of RSEM's --strand-specific option, Shurjo needs to use --forward-prob 0.

RSEM will automatically translate this into the --nofw option for bowtie.

• 1 views • link
updated by Ram Ram • written by scott.swanson • scott.swanson

After reading the docs, I would guess that you want to use --strand-specific or, equivalently, --forward-prob=1.0.

http://deweylab.biostat.wisc.edu/rsem/rsem-calculate-expression.html

--strand-specific

The RNA-Seq protocol used to generate the reads is strand specific, i.e., all (upstream) reads are derived from the forward strand. This option is equivalent to --forward-prob=1.0. With this option set, if RSEM runs the Bowtie/Bowtie 2 aligner, the '--norc' Bowtie/Bowtie 2 option will be used, which disables alignment to the reverse strand of transcripts. (Default: off)

• 0 views • link
updated by Ram Ram • written by Kamil Kamil

I just stumbled upon this while struggling with rsem strand directions. Do not simply use --strand-specific for stranded reads, as you would need to specify direction of the read. If you've used TopHat you'll be familiar with parameters such as --library-type fr-firststrand, do not assume fr-firststrand is equivalent to --strand-specific or --forward-prob=1.0, it turns out that as @scott.swanson suggested you will need to use --forward-prob 0 instead!

Check out this blog for more iformation:

• 0 views • link
updated by Ram Ram • written by Yuka Takemon • Yuka Takemon

Log in to answer this question.

More posts like this