2
$\begingroup$

I have gene expression data from RNAseq, specifically: log2(x+1) transformed RSEM normalized counts.

How can I convert these data from multiple samples to determine fold change in gene expression?

$\endgroup$
3
$\begingroup$

You should use a proper statistical framework for RNA-seq dfferential analysis (which includes FC calculation). Standard tools for this are (among others) edgeR or DESeq2. You could use tximport to import RSEM outputs into R and then use its output for e.g. DESeq2. The linked manual provides example code for this. Please read the manuals of edgeR or DESeq2 towards the respective downstream analysis. Both are extensive and compreshensive.

$\endgroup$
3
  • $\begingroup$ Thanks a lot for this information! So, just to confirm, that means that there is no 'simple formula' that can convert RSEM to fold change? Thanks again. $\endgroup$ Mar 7 '20 at 23:28
  • 2
    $\begingroup$ No. By the way, fold changes without statistics are basically meaningless. Just run it through DESeq2 and then take the results for your analysis. $\endgroup$
    – ATpoint
    Mar 8 '20 at 10:44
  • $\begingroup$ Thank you very much. I will give this a try. $\endgroup$ Mar 8 '20 at 19:54

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.