This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GSEA analysis RNA-seq data

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!

rna-seq gsea pathway gene ontology
• 2,071 views • link
updated by GenoMax GenoMax • written by Marco Pannone • Marco Pannone

1 answer

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

• 0 views • link
written by MatthewP MatthewP

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?

• 0 views • link
written by Marco Pannone • Marco Pannone

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

• 0 views • link
written by MatthewP MatthewP

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.

• 0 views • link
written by Marco Pannone • Marco Pannone

Log in to answer this question.

More posts like this