I'm currently doing RNAseq differential expression analysis, on which I've got no significative genes with FDR < 0.05. (I'm working on human tumor biopsies data, with 111 samples.)
However, when I perform GO enrichment analysis on the top hits (p-value < 0.05, logFC > 1 or < -1), it gives significative enriched pathways, which seems consistent from a biology point a view.
This bring me two questions :
- Could those results be relevant? Would the biological signal detected with GO enrichment in the top hits counteract the fact that there is no significative genes detected?
- If so, how could I illustrate those findings, I've tried to do heatmaps on a subset of genes belonging to a specific pathway, but, as the gene expression differential is rather low between the two studied conditions, there is no relevant clustering of the samples in the figure between the two conditions. (see below for example the type of expression patterns that I get for one gene between the two studied conditions)
Thank you in advance for any input