This is a test version of Biostars. For the public version, visit https://www.biostars.org.
heatmap for coding and noncoding genes

Hi friends, I did differential expression for rna seq data: I did for coding genes also I did for all genes(coding+noncoding). I want to show significant genes in heatmap. Do you think heatmap is good for significant coding genes or all genes?

rna-seq
• 2,017 views • link
updated by cristopher.biotec • cristopher.biotec • written by Rob • Rob

2 answers

It is your choice and it depends on the message that you want to convey in your work. You do not really need anybody else's verification for this task.

  • All genes == unsupervised clustering
  • Statistically significant genes == supervised clustering

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

I may have missed your point slightly in relation to non-coding versus protein coding genes. If you have conducted your entire analysis separately for protein coding genes and non-coding genes (that is, normalisation and differential expression analysis separately), then it is my view that they also should be shown separately in different heatmaps.

There is no problem conducting your analysis by just focusing on protein coding genes alone, and there is neither a problem focusing on protein coding + non-coding combined.

As you can see, there are no rules here... you can present your results in the way that helps to get the information best to the reader.

• 0 views • link
written by Kevin Blighe Kevin Blighe

Thanks Kevin I have differentially expressed genes for coding genes alone and also all genes together(Coding+noncoding). I have volcano plot and now I want to plot heatmap for them.

• 0 views • link
written by Rob • Rob

For the heatmap, you could try to follow this: A simple tutorial for a complex ComplexHeatmap

, or, at the most basic level, you just need:

heatmap(assay(vst)[sigGenes,])

, assuming here that you have used DESeq2 and that your normalsied+transformed data is contained in assay(vst) and that your vector of statistically significant genes, whose ID type should match that of the rownames of assay(vst), are held in sigGenes.

• 0 views • link
written by Kevin Blighe Kevin Blighe

You can show those results in the enhanced volcano plot using a different shape for coding and non-coding genes

• 0 views • link
written by cristopher.biotec • cristopher.biotec

Good suggestion - I am also a developer of EnhancedVolcano.

• 0 views • link
written by Kevin Blighe Kevin Blighe

Hi Cristopher Thanks good suggestion. I did volcano plot before just for significant genes all same shape. But how can I do this with different shape for coding and non-coding genes?

• 0 views • link
written by Rob • Rob

Thanks Kevin I use your EnhancedVolcano a lot. This is really amazing. Appreciate that.

• 0 views • link
written by Rob • Rob

Thank you Sir. You got it working?

• 0 views • link
written by Kevin Blighe Kevin Blighe

Yes dear Kevin Your Enhancedvolcano is awesome and very useful. I wish you would do something like this for heatmap.

thanks

• 0 views • link
written by Rob • Rob

Log in to answer this question.

More posts like this