Making sense of gene and proteins lists with functional enrichment analysis

9 Reporting results

Once we have done an enrichment analysis - how to communicate the results? Like everything - it depends what you need to say. Some examples are below.

9.1 In text

If you just need to emphasise that the differential expression relates to a condition of interest, you don’t need much:

Genes differentially expressed after SH-SY5Y cell differentiation were enriched for GO term “regulation of neurogenesis”, adjusted pval < 0.01).

And in the methods:

Enrichment calculated for differentially expressed genes with the g:GOSt enrichment tool (Raudvere et al, 2019) using a background of tested genes.

9.2 As a table

For a more complete view, a table of the significant or top n terms can be useful. E.g In a supplementary figure

The top 10 enriched GO terms for the differentially expressed genes;

GO:BP Term Term ID Adjusted p-value Term Size Num DE Genes
system process GO:0003008 2.90E-04 1243 45
nervous system development GO:0007399 4.19E-04 1985 60
regulation of cell development GO:0060284 1.98E-03 795 33
central nervous system development GO:0007417 2.67E-03 805 33
regulation of neurogenesis GO:0050767 5.59E-03 707 30
regulation of nervous system development GO:0051960 5.78E-03 790 32
regulation of cell differentiation GO:0045595 6.03E-03 1471 47
multicellular organismal process GO:0032501 2.04E-02 5300 111
system development GO:0048731 2.88E-02 3668 85
neurogenesis GO:0022008 4.36E-02 1375 43

9.3 As a figure

Plots of -log p-value are a popular option for a figure. By taking the exponent of the p-value bigger bars indicate more significance.

For instance, this figure is taken from the (Pezzini et al. 2017) paper. Note that in this figure shows 2 specifics categories from the IPA (ingenuity) database using the IPA tool (not covered here), so the terms themselves differ.

You can also build on such a plot with number of genes and the like. For an example, see the ‘bubble’ chart produced by a package called pathfindR here: https://www.biostars.org/p/322415/


There are also tools like cluego (a cytoscape plugin) that build enriched terms into a network : http://apps.cytoscape.org/apps/cluego