Question: Meaningful DE using rpkm values in a RNA-Seq analysis
0
gravatar for Seq225
4 weeks ago by
Seq22590
Seq22590 wrote:

Is it conventional/required to use a cut-off value to make rpkm values meaningful? For example, in a transcriptome analysis if a contig/gene (mRNA) has rpkm values of 0.025 and 0.0025 in two different conditions, the DE is 10. The ratio is high but the actual values (0.025 vs 0.0025) are very low. Should I consider this DE as significant (with multiple reps and p-value <.05)?

Is there a guideline to follow/decide what is meaningful and what not?

Thanks!

ADD COMMENTlink modified 4 weeks ago by c.chakraborty140 • written 4 weeks ago by Seq22590
3
gravatar for swbarnes2
4 weeks ago by
swbarnes26.5k
United States
swbarnes26.5k wrote:

The ratio is high but the actual values (0.025 vs 0.0025) are very low.

No one can really understand what those values mean, because they are so corrected. If you had fanatically deep coverage, that difference might be perfectly real and valid.

Should I consider this DE as significant (with multiple reps and p-value <.05)?

Don't try to assign a p-value to differences in RPKM. Put the raw counts through DESeq or edgeR.

ADD COMMENTlink written 4 weeks ago by swbarnes26.5k
1
gravatar for c.chakraborty
4 weeks ago by
c.chakraborty140
c.chakraborty140 wrote:

Use raw read counts for genes and put them in DESeq2 or EdgeR for analysis. With RPKM, your expression is normalized to the length of the transcript. Thus if you have a very long transcript, the normalized value will be lower. As it is a practice in RPKM based analysis, to use a cutoff for expression, for example 0.1 RPKM, hence there are chances that you will lose a lot many interesting transcripts. Therefore use the read count values generated by htseq counts or feature counts.

ADD COMMENTlink written 4 weeks ago by c.chakraborty140
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: 774 users visited in the last hour