Question: How to draw gene expression heatmap ?
0
gravatar for scheme4193
18 months ago by
scheme419380
scheme419380 wrote:

I want to draw heatmap using differential expression genes.

First, I mapped using hisat2, and then counted genes using htseq. From this expression count matrix, I used DESeq2 to detect DEGs.

Now I have DEG list. I want to draw heatmap about these DEGs.

Which expression value do many people use to draw heatmap ? In my pipeline (hisat2 → htseq-count → DESeq2), which data can be used to draw genes expression heatmap ?

rna-seq next-gen R gene • 1.1k views
ADD COMMENTlink modified 18 months ago by Kevin Blighe69k • written 18 months ago by scheme419380

Depends of what you want to show. Usually heatmap are made on the log of the counts

http://bioconductor.org/packages/devel/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#heatmap-of-the-count-matrix

Something like :

pheatmap(log2(counts(dds, normalized=TRUE))...)
ADD REPLYlink written 18 months ago by Bastien Hervé4.9k

Thank you for your reply ! I appreciate you sharing URL. I will use it as a reference.

ADD REPLYlink written 18 months ago by scheme419380
1
gravatar for Kevin Blighe
18 months ago by
Kevin Blighe69k
Republic of Ireland
Kevin Blighe69k wrote:

You just need to subset your normalised expression data (in DESeq2, it is returned by counts(dds, normalized=TRUE), as pointed out by Bastien), but you can also use the normalised and then transformed expression data - these are produced via rlog() or vst(), then accessed via assay(). For example:

rlogcounts <- rlog(dds, blind = FALSE)
dataHeatmap <- assay(rlogcounts)

There is no standard for what to actually use in a heatmap. You can just use whatever you want. I have seen a wide range of things:

  • binary data
  • raw counts
  • normalised counts
  • log-normalised counts
  • normalised, transformed counts
  • log [base 2] ratios
  • Z-scores

You just need to specifically state what it is that you are using.

There is a wealth of information at your fingertips via a search engine, too.

Kevin

ADD COMMENTlink written 18 months ago by Kevin Blighe69k
1

Thank you for your polite reply. So, there are no uniform indicator. I'm going to try some scores. Thank you !

ADD REPLYlink written 18 months ago by scheme419380
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1983 users visited in the last hour