feature count strand explained
2
1
Entering edit mode
7.0 years ago
popayekid55 ▴ 110

Hello, I am very new to rna-seq. I have got the RNA-seq data for bacteria (Hiseq 150*2). From lab it was told that its a directional library (NEB ultra directional kit). Now i assume the data is from Sense/Forward strand.

I have used featureCount with -s 1 and only 1% of the reads were assigned to feature. Then i tried with -s 2 (reverse strand) and >96% of reads were assigned. Could someone explain strand concept in feature count.

Thanks

RNA-Seq • 6.0k views
ADD COMMENT
2
Entering edit mode
7.0 years ago
lshepard ▴ 480

Hi, check out the following links for the description of strand specificity. From the output from featurecounts, looks like your first strand was the reverse (ISR in the first link). They really helped me:

http://salmon.readthedocs.io/en/latest/library_type.html

http://onetipperday.sterding.com/2012/07/how-to-tell-which-library-type-to-use.html

ADD COMMENT
1
Entering edit mode
7.0 years ago

You can look at the strandedness of your data in a browser such as igv. Are the reads in the same direction as your genes or in the reverse direction? For most recent stranded kits the reads are most often reverse stranded, which as you observed gave the best result.

ADD COMMENT
0
Entering edit mode

For most recent stranded kits the reads are most often reverse stranded,

Why is it so??

ADD REPLY
1
Entering edit mode

Because that's the way the library prep works, the second strand is not amplified (in the case of dUTP method).

ADD REPLY
0
Entering edit mode

thanks. If i count the data with unstranded option how does that impact the?

ADD REPLY
0
Entering edit mode

In all the excitement it seems you lost a word of your sentence.

ADD REPLY
0
Entering edit mode

If i count the data with unstranded option how does that impact the differential expression analysis/any down stream analysis?

ADD REPLY
0
Entering edit mode

Exons which are antisense will get counted incorrectly or not at all because the reads cannot unambiguously get attributed to one of both transcripts.

ADD REPLY

Login before adding your answer.

Traffic: 2707 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