DEGreport to plot nice RNA-seq figures

Lorena Pantano bio photo By Lorena Pantano

Differentially gene expression analysis with RNA-seq data is quite common nowadays, and there are pretty good Bioconductor packages for that: limma::voom, DESeq2

The code for that part is quite simple, being super quick to get a list of de-regulated genes. However, downstream analyses vary a lot depending on the project itself. But I found myself doing the same plots and analyses many times for different project, so I put together a bunch of plots and analyses using code from my colleagues at work (@HSPH bioinformatics core) and myself.

There are many possibilities, including QC figures, but the features that seem more interesting are:

  • plot of individual genes using multiple variables from the metadata. This would be an evolution of plotCounts function from DESeq2. Fig1
  • plot volcano plots labeling some genes. Fig2
  • clustering set of genes to find common patterns expression. Most useful for time seria data. You get as well the list of genes that follows that pattern, so you can use it for functional analysis, for instance. Fig3

note: maybe some features only available with development version

fig1 fig2 fig3

Package available at BioC page

The code is available here

Hope you enjoy!

Any suggestion and/or issue is welcome here

“sincerity in bioinformatics: try hard to explain the limitation of your analysis to non-computational researchers”