DEGreport to plot nice RNA-seq figures
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
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”