Question: GSEA analysis RNA-seq data
gravatar for marco
9 months ago by
Oslo, Norway
marco20 wrote:


I am performing the GO/pathway analysis of my RNA-seq data using GSEA. I have many DEGs, about 3000 in total, and I got many enriched pathways with both positive and negative NES. However, all or most of the enriched pathways have an FDR > 0.05 and therefore are not classified as statistically significant. Most of such enriched pathways are characterized by more than 30 DEGs, which is a reasonable number of differentially expressed genes associated with a certain pathway.

What could be the reason for that?

I would gladly appreciate any comment or suggestion and thanks in advance!

ADD COMMENTlink modified 9 months ago by GenoMax95k • written 9 months ago by marco20
gravatar for MatthewP
9 months ago by
MatthewP850 wrote:

How did you perform GSEA? What's your input data? You need to input ALL genes data to GSEA not only DEGs.

ADD COMMENTlink written 9 months ago by MatthewP850

I am running GSEA through the WebGestalt website, where it asks me to report a list of two columns: one should contain the gene identifier and another one the associated "scores", which I assume are the log2FC values of the retrieved DEGs. I do upload together both the UP and DOWN DEGs.

Is there something wrong I am doing?

ADD REPLYlink modified 9 months ago • written 9 months ago by marco20

Log2FC is fine. You need to upload all genes, including those not DEGs.

ADD REPLYlink written 9 months ago by MatthewP850

Thanks a lot! I will do so then. I was including only the DEGs with a log2FC threshold of > 0.5 for UP DEGs and < -0.5 for DOWN DEGs. Including all the genes means that I also have to remove the adjusted pvalue threshold? I have tried this way and some pathways now pass the FDR threshold and are classified as statistically significant.

ADD REPLYlink written 9 months ago by marco20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1124 users visited in the last hour