Generating heatmaps of pre-selected candidate genes after deseq2
2
1
Entering edit mode
3.6 years ago

Hi biostars, Is it possible to make a heatmap for a preselected list of genes after deseq2 DGE analysis? I wanted to represent the log2FC of 20 differentially expressed candidate genes as a heatmap.

RNA-Seq DESEq2 heatmap • 3.4k views
0
Entering edit mode

Hi, did you solve this? I want to do a similar thing with a number of candidate genes but using rlog data to include in a pheatmap

3
Entering edit mode
3.6 years ago
e.rempel ★ 1.0k

Hi,

it is absolutely possible. Have a look at this tutorial. You could slice the matrix of normalized counts to contain only interesting genes:

pheatmap(assay(ntd)[intersting_genes,], cluster_rows=FALSE, show_rownames=FALSE,
cluster_cols=FALSE, annotation_col=df)

0
Entering edit mode

How do you "slice the matrix of normalized counts to contain only interesting genes", considering a list/data frame named "GENES_TO_LABEL" (containing the names of the genes you wish to label) which you've imported from a txt file?

0
Entering edit mode

Hi, did you solve this? I want to do a similar thing with a number of candidate genes but using rlog data to include in a pheatmap

3
Entering edit mode
6 months ago
svp ▴ 440

You can use as follows: use a vector that contains list of genes you wish.

genes <- c("gene_1", "gene_2", ...)
mat <- assay(rld)[genes, ]
mat <- mat - rowMeans(mat)
anno <- as.data.frame(colData(rld)[c("Condition")])
pheatmap(mat, annotation_col = anno)