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

Hello!

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
2
gravatar for MatthewP
9 months ago by
MatthewP850
China
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
2

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.

Help
Access

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
_