Question: Strand specific in RNA-seq
0
gravatar for swarajgaikwad15
7 months ago by
Tezpur University, Assam, India
swarajgaikwad150 wrote:

One small confusion is there, I am using paired-end RNA-seq data and Platform for my sequencing was 'Illumina NextSeq 500 (Homo sapiens)' I have used the following command for quantification,

 featureCounts -p -t exon -g gene_id -a hg38ucsc.gtf -o SRRXXX.bam.txt SRRXXX.bam

The output which I do have looks like

Geneid   Chr      start              End                             Strand   Length    
NMXYZ    Chr1;Chr1;...  67092176;67096252;...   67093604;67096321;...   +;+;  2203,..  0

Now here I can see 'Strand' column which contains + and - both options in different rows, Is it correct? and if not then am I supposed to use -s 0 or -s 1 or -s 2 option in the above command?

rna-seq subread featurecounts • 289 views
ADD COMMENTlink modified 7 months ago by caggtaagtat650 • written 7 months ago by swarajgaikwad150

You did not tell us which kit was used for generating the library, and that determines if you are sequencing strand-specific.

ADD REPLYlink written 7 months ago by WouterDeCoster40k

SMARTer Ultra Low Input RNA Kit for Sequencing - v4 kit (Clontech Laboratories).

ADD REPLYlink written 7 months ago by swarajgaikwad150
1

Based on the information available from the Clontech website, this is not a stranded kit and therefore you have a normal unstranded library, and any inference of the strand-of-origin is impossible.

ADD REPLYlink modified 7 months ago • written 7 months ago by ATpoint19k
1
gravatar for caggtaagtat
7 months ago by
caggtaagtat650
caggtaagtat650 wrote:

Hi, every row stands for a certain feature which is either encoded on the plus or minus strand, this is therefor correct. The option -s 0 on the other hand would state, that the data originates from the sequencing of an "unstranded" library. This means, that you don't know if the synthesized read comes from the forward or reverse strand of your cDNA. These days, its a standart to use stranded library protocolls, which enables to generate reads which for example are reverse complementary to the original fragment. Again, this does not give you information about the strand your feature was located on, but which strand of your cDNA was sequenced by synthesis during read generation.

ADD COMMENTlink modified 7 months ago • written 7 months ago by caggtaagtat650

Thank you very much. That's a really good explanation. So finally I don't have to use either -s 1 or -s 2 which is typically meant for specifically for forward or reverse strand library? And other small question what is meant by '0' column present after length column?

Geneid   Chr      start              End                             Strand   Length    
NMXYZ    Chr1;Chr1;...  67092176;67096252;...   67093604;67096321;...   +;+;  2203,..  0
ADD REPLYlink written 7 months ago by swarajgaikwad150

Like ATpoint already said, you dont have to state -s 1 or -s 2, since your library seems to be unstranded. Im not familier with feature counts and the heading seems to be missing, but if there would be a heading like "Name_of_your_imput.sam" the 0 would be the read count.

ADD REPLYlink written 7 months ago by caggtaagtat650

Yes, It's correct Got it Thank you.

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