HTSeq-count TruSeq RNA Exome Lib Prep
1
0
Entering edit mode
15 months ago
tdfyoder ▴ 30

Hello,

I observed a high percentage of "no features" while running HTseq w/ the --stranded yes option enabled (>80%). The library prep kit I am using is Illumina TruSeq RNA Exome which generates stranded data. If I run HTseq-count w/ strand == "no" or "reverse" I observe expected levels of "no features" <2%.

I have checked for compatibility between reference genome and gtf coordinates (chr1 vs 1) and mapped to a second reference genome to assure gtf ranges are concordant.

To my knowledge this library prep kit does not use dUTP.

CFLAR (+ strand) IGV export seen below.

I came across this paper describing the best option to call based on strand synthesis. The Illumina guide describes the molecular reactions as,

The RNA is fragmented using divalent cations under elevated temperature. cDNA is generated from the
cleaved RNA fragments using random priming during first and second strand synthesis. Then, sequencing
adapters are ligated to the resulting double-stranded cDNA fragments. The coding regions of the
transcriptome are captured from this library using sequence-specific probes to create the final library.

Can you please advise on the best --stranded option for this chemistry? Thank you.

-T

HTseq strandness TruSeqRNAExome nofeatures • 723 views
ADD COMMENT
0
Entering edit mode

CFLAR (+ strand)

ADD REPLY
1
Entering edit mode
15 months ago
tdfyoder ▴ 30

from the developer:

https://github.com/simon-anders/htseq/issues/106

If I recall correctly, "reverse" is correct for TruSeq, because TrueSeq reads are 3' to 5', i.e. reverse to the "biological" reading direction.

That the non-reversed direction is called "yes" rather than "forward" hast historical reasons: When I write that, TruSeq did not exist yet, and the first stranded protocols where 5' to 3'.

To be more precise: Read 1 is 3' to 5'. Read 2 is the opposite direction, of course. For paired-end reads, the "stranded" parameter should describe the orientation of read 1.
ADD COMMENT

Login before adding your answer.

Traffic: 1174 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6