This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Most variable gene from the differentially expressed genes

I m getting deseq2 output of differentially expressed genes with p ,0.05 .So from that list of differentially expressed i want to take out the most variable genes getting expressed .Is there any way to do that .

Any suggestion or help would be highly appreciated.

r rna-seq
• 3,858 views • link
updated by Satyajeet Khare Satyajeet Khare • written by 1769mkc 1769mkc

1 answer

I think, by "most variable" you mean in terms of "differences in expression between the two samples". You need to have fold change values for that. Generally, genes with p<0.05 and log2FC > 1 (more than two fold) are selected.

• 0 views • link
written by Satyajeet Khare Satyajeet Khare

well yes i did that too, but when i plot heatmap i dont see much differences,and lets say i have differentially expressed gene list of 1000 ,i dont think that would be wise to take into downstream analysis .I want to narrow it down..make it more stringent

• 0 views • link
written by 1769mkc 1769mkc

The p-value and log2FC thresholds are entirely arbitrary. Just rank your genes in decreasing order of differential expression and take the top N, N being whatever is reasonable in your context. Sometimes, it may help to plot the values in decreasing order. You may notice an elbow which may suggest a convenient cut-off.

• 0 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

okay.I will try it, I have another query would you suggest to use Z score normalisation of the rna-seq data which are already normalised to plot heatmap ?

• 0 views • link
written by 1769mkc 1769mkc

It depends on what your RNA seq data is. Search for related posts here on Biostars.

• 0 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

Log in to answer this question.

More posts like this