This quick tutorial will guide you through the generation of an Enrichment Map for an analysis performed using g:Profiler (Functional Profiling of Gene List from large-scale experiments).
Download the test data:
Description of the tutorial files contained in the gProfilerTutorial folder:
12hr_topgenes.txtList of top genes expressed in Estrogen dataset at 12hr - Official Gene Symbol.
24hr_topgenes.txtList of top genes expressed in Estrogen dataset at 24hr - Official Gene Symbol.
Step 1: Generate g:Profiler output files¶
Go to g:Profiler website
Select and copy all genes in the tutorial file 12hr_topgenes.txt in the Query box. Make sure that your list contains only official gene symbol (HUGO).
In Options, check Significant only, No electronic GO annotations
Set the Output type to Generic EnrichmentMap
Show advanced options
Set Min and Max size of functional category to 3 and 500 respectively.
Select 2 for Size of Q&T
On the right panel, choose the Gene Ontology Biological process and Reactome
Set Significance threshold to Benjamini-Hochberg FDR
Click on g:Profile! to run the analysis
[Note] - if some of your identifiers in your query have multiple mappings in g:Profiler! by default they get excluded. If this happens you will see the following above the g:Profiler! results:
Click on above link to manually map each gene to its correct annotation
Click on Resubmit query to update your results with the specified mappings.
If the identifier discrepancy warning is ignored there might be differences between the number of genes g:Profiler attributes to a particular gene set and those associated with it in the Enrichment Map.
Download g:Profiler data as gmt name Note, you will have to unzip the folder
Download the result file: Download data in Generic Enrichment Map (GEM) format
Repeat these steps for the 24hrs time-point and the file 24hr_topgenes.txt
Link to a step by step tutorial:
Step 2: Generate Enrichment Map with g:Profiler Output¶
g:Profiler output files from Step 1:
- Open Cytoscape
- In the menu bar, locate the App tab and then select > EnrichmentMap > Create Enrichment Map
- Make sure the Analysis Type is set to generic(ex:gProfiler)
- Please select the following files by clicking on the respective (…) button and selecting the file in the Dialog:
- GMT / hsapiens.pathways.NAME.gmt
- Dataset 1 / Enrichments: gProfiler_results_12hr.txt
- Tune Parameters
- P-value cut-off: 1
- Q-value cut-off: 1
- Overlap Coefficient cut-off: 0.5
- Click on the Build radio button at the bottom of the panel to create the Enrichment Map
- In the menu bar, Go to View, and activate Show Graphics Details
- In the control panel, go to Style, click on Label and select EM1_GS_DESCR in the Column dropdown. This will label nodes with names rather thsn GO IDs. The selected value may be EM2_GS_DESCR or other if you have more than one Enrichment Map open in Cytoscape.
Step 3: Examining Results¶
- Node size corresponds to the number of genes in dataset 1 within the geneset
- Colour of the node corresponds to the significance of the geneset for dataset 1.
- Edge size corresponds to the number of genes that overlap between two connected genesets.