This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GSEA top genes vs all platform

Hello,

I have done an analysis in LIMMA to obtain the differential expressed genes. However, now I have 2 questions when I perform GSEA:

  • Must I order my samples by fold change (first negative and after positive forld change, por example) in a preranked analysis? Is it better to order by p.val? Can I give the list obtained by limma directly in GSEA?

  • GSEA is usually done with the top genes (let's say, the top 50-100 genes) or must I use THE WHOLE platform (all the 30K probes)?

Thanks in advance

rna-seq functional analysis gsea
• 2,920 views • link
updated by Charles Warden Charles Warden • written by ra2967 • ra2967

Have a loot at this thread, may be helpful

• 0 views • link
written by Prakash Prakash

You can also use gene set tests included in the limma package, i.e. CAMERA and ROAST. CAMERA, in particular, is way, way faster than GSEA. Prof. Smyth has made the MSigDB gene sets available in easy to ingest .rda files at his website: http://bioinf.wehi.edu.au/software/MSigDB/index.html

• 0 views • link
written by glocke01 glocke01

2 answers

GSEA is typically performed on the whole platform, without the differential expression step.

However, you are correct that a pre-ranked list is another option.

• 0 views • link
written by Charles Warden Charles Warden

Instead of using pre-ranked with 50-100 genes, you can use entire set. Then check how many genes have core GSEA enrichment and are differentially expressed.

• 0 views • link
written by Rashedul Islam Rashedul Islam

Log in to answer this question.

More posts like this