Obtaining gene sets for GSEA
 

Hi All,

I work with Bacillus subtilis and I want to perform some GSEA on an RNA-seq data set. Looking at the literature, I've seen that it is very much possible to do this. I know that gene sets can be constructed as well, but I'm not super great at automation and it seems like it would require a massive amount of time to do something like that manually.

In some of the literature I've read, I've seen that gene sets could be obtained from KEGG. However, browsing the website, I'm not certain how gene sets could be extracted from KEGG.

Can anyone help me in figuring out how to obtain gene sets for GSEA?

Thanks!

Hi,

I'm not sure what you have now. But if you have a list of differentially expressed genes with their fold change, such as the results from edgeR or DESeq2 package, I think the package `clusterProfiler` is good for you to do GO/GESA/KEGG directly. You could find their detailed manual here -> http://yulab-smu.top/clusterProfiler-book/chapter1.html.

best