Question: TruSeq strand-specificity in rsem-calculate-expression
1
gravatar for shurjo.sen
3.4 years ago by
shurjo.sen10
United States
shurjo.sen10 wrote:

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 • 2.7k views
ADD COMMENTlink modified 2.2 years ago by Yuka Takemon20 • written 3.4 years ago by shurjo.sen10
2
gravatar for scott.swanson
3.3 years ago by
United States
scott.swanson20 wrote:

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.

ADD COMMENTlink written 3.3 years ago by scott.swanson20
0
gravatar for Kamil
3.4 years ago by
Kamil1.9k
Boston
Kamil1.9k wrote:

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

ADD COMMENTlink written 3.4 years ago by Kamil1.9k
0
gravatar for Yuka Takemon
2.2 years ago by
Yuka Takemon20
The Jackson Laboratory
Yuka Takemon20 wrote:

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: http://likit.github.io/category/next-gen-sequencing.html

ADD COMMENTlink written 2.2 years ago by Yuka Takemon20
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: 1020 users visited in the last hour