This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Meaningful DE using rpkm values in a RNA-Seq analysis

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!

rna-seq r chip-seq sequencing assembly
• 2,139 views • link
updated by c.chakraborty c.chakraborty • written by Seq225 Seq225

2 answers

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.

• 0 views • link
written by swbarnes2 swbarnes2

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.

• 0 views • link
written by c.chakraborty c.chakraborty

Log in to answer this question.

More posts like this