Almost all the tools below allow to define a custom background. Use a custom background if you can.
For single cell RNASeq analysis, the background should be the list of genes present in the Seurat object. For more info on the background I refer to the Function analysis course on the e-learning system.
If you’re a VIB scientist who works on human, mouse and rat, check Ingenuity Pathway Analysis. VIB offers this software for free to VIB scientists. You can find more info in the Functional analysis course on the e-learning system. Contact email@example.com if you want us to organize a training on this software.
Other useful websites for ORA and GSEA:
Easy to use software to install on your computer for GSEA:
More information on these tools can be found on:
R/Bioconductor packages for ORA and GSEA:
Load required packages
I’m not sure it works on Gene Symbols you might need to map them first to EntrezIDs (see above)
For the GO enrichment analysis you need:
> GO2 <- enrichGO(gene=MarkersList,universe=Background,keyType=”ENTREZID” (or “SYMBOL”), OrgDb=org.Hs.eg.db,ont=”BP”,pAdjustMethod=”BH”,qvalueCutoff=0.05,readable=TRUE)
Output results from GO analysis to a table
> go_summary <- data.frame(GO2)
The package also has a function enrichKEGG() for pathways. The package has additional functions gseGO(), gseKEGG() for GSEA but then your input also needs to contain values (log fold changes, p-values…), and not only gene names.