joyplot for GSEA result

I am very glad to find that someone figure out how to use ggjoy with ggtree.

I really love ggjoy and believe it can be a good tool to visualize gene set enrichment (GSEA) result. DOSE/clusterProfiler support several visualization methods.

[Bioc 3.5] NEWS of my BioC packages

I have 8 packages published within the Bioconductor project.

A new package treeio was included in BioC 3.5 release.

dotplot for GSEA result

For GSEA analysis, we are familar with the above figure which shows the running enrichment score. But for most of the software, it lack of visualization method to summarize the whole enrichment result.

In DOSE (and related tools including clusterProfiler, ReactomePA and meshes), we provide enrichMap and cnetplot to summarize GSEA result.

[Bioc 34] NEWS of my BioC packages

I have 7 packages published within the Bioconductor project.

A new package meshes was included in BioC 3.4 release.

leading edge analysis

leading edge and core enrichment

Leading edge analysis reports Tags to indicate the percentage of genes contributing to the enrichment score, List to indicate where in the list the enrichment score is attained and Signal for enrichment signal strength.

It would also be very interesting to get the core enriched genes that contribute to the enrichment.

Now DOSE, clusterProfiler and ReactomePA all support leading edge analysis and report core enriched genes.

[Bioc 33] NEWS of my BioC packages

Today is my birthday and it happened to be the release day of Bioconductor 3.3. It’s again the time to reflect what I’ve done in the past year.

ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualization

My R/Bioconductor package, ReactomePA, published in Molecular BioSystems.

ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualization

G Yu, QY He. Molecular BioSystems 2016, 12:477-479

Received: 05 Oct 2015 Accepted: 20 Nov 2015 Online: 23 Nov 2015

[Bioc 32] NEWS of my BioC packages

In BioC 3.2 release, all my packages including GOSemSim, clusterProfiler, DOSE, ReactomePA, and ChIPseeker switch from Sweave to R Markdown for package vignettes.


To make it consistent between GOSemSim and clusterProfiler, ‘worm’ was deprecated and instead we should use ‘celegans’. As usual, information content data was updated.

functional enrichment for GTEx paper

The ENCODE consortium has recently published a great paper on Gene Expression from the GTEx dataset. A criticism raised on pubpeer is that the gene ontology enrichment analysis was done with DAVID which has not been updated in the last five years.

The result is shown below:

dotplot for enrichment result

This is a feature request from clusterProfiler user. It’s similar to what I implemented in clusterProfiler for comparing biological themes. For comparing different enrichment results, the x-axis represent different gene clusters while for a single enrichment result, the x-axis can be gene count or gene ratio. This is actually similar to traditional barplot, with dot position as bar height and dot color as bar color. But dotplot can represent one more feature nicely by dot size and it can be a good alternative to barplot.