Single end RNAseq data strandedness infer-experiment-py
2
0
Entering edit mode
5.3 years ago
nilus1432 ▴ 30

Hi all,

I have a single end RNAseq data and would like to understand the strandedness of the data.

  1. From wet-lab input I know "Stranded cDNA library was generated by reverse transcribing the RNA molecules".

  2. I used infer-experiment-py The output is:

This is SingleEnd Data
Fraction of reads failed to determine: 0.0795
Fraction of reads explained by "++,--": 0.0703
Fraction of reads explained by "+-,-+": 0.8502
  

So it's stranded but is it forward or reverse? I do not understand the help given here:

Does it means that its reverse stranded and I have to use s -2 option in featureCounts, "reverse" strandedness in htseq, --rf in StringTie?

Thank you

rna-seq • 2.4k views
ADD COMMENT
2
Entering edit mode
9 months ago
DareDevil ★ 4.3k

Based on the output provided by infer_experiment.py, the strandedness can be inferred as follows:

  • Fraction of reads failed to determine: 0.0795 This indicates the fraction of reads for which the strandedness could not be determined. It could be due to various reasons such as low-quality reads or issues with the library preparation. This value is relatively small, indicating that a majority of the reads could be properly assigned a strandedness.

  • Fraction of reads explained by "++,--": 0.0703 This fraction represents the reads that align in the forward orientation ("++") or the reverse orientation ("--"). In other words, the reads are mapping to the same strand as their corresponding reference transcripts. In this case, it suggests that 7.03% of the reads are aligned in the forward orientation and 7.03% are aligned in the reverse orientation.

  • Fraction of reads explained by "+-,-+": 0.8502 This fraction represents the reads that align in the forward-reverse ("+-") or reverse-forward ("-+") orientation. In other words, the reads are mapping to the opposite strand compared to their corresponding reference transcripts. In this case, it suggests that 85.02% of the reads are aligned in the forward-reverse orientation ("-+") and 85.02% are aligned in the reverse-forward orientation ("+-").

Based on this information, we can conclude that the library appears to be stranded in a reverse manner. The majority of the reads align to the reverse-complementary strand compared to their corresponding reference transcripts, while a smaller fraction aligns to the same strand.

ADD COMMENT
1
Entering edit mode
5.3 years ago
h.mon 35k

The majority (85%) of your reads falls in the case described by:

+-,-+
read mapped to ‘+’ strand indicates parental gene on ‘-‘ strand
read mapped to ‘-‘ strand indicates parental gene on ‘+’ strand
  

So yes, your library is reverse stranded. As featureCounts is very fast, you can quickly confirm this by running it with both settings (-s 1 and -s 2), the correct setting will have a much lower count for NoFeature than the incorrect one.

ADD COMMENT

Login before adding your answer.

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