This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Deeptools for ChIP-Seq

Hello,

I am new to deeptools and have referred to this post Y axis PlotHeatMap Deeptools. I am still not clear on what the Y axis represents in my case enter image description here

The figure shows 2 ChIP-Seq replicates (1_normalized and 2_normalized); these are generated using bamCompare. Also, is there a way to combine my ChIP-seq replicates in deeptools?

chip-seq deeptools heatmap

Your image is broken and should be embedded directly.

Well, what files did you feed into bamCompare? Both of your replicates or each replicate individually against a control (IGG, input, etc)? What's in your bed files being fed to plotHeatMap?

You can use samtools to merge bam files, though it's usually best to show replicates individually when possible.

Okay. Thank you for your replies! I fed in the bam files - each replicate against its control - into bamCompare. The bigwig files were then used in computeMatrix and finally the matrix file was used to generate the heatmap. How can I obtain the top genes with the peaks in/ near the TSS from this?

were you able to figure out what the Y-axis meant?

1 answer

The y-axis is whatever value was used to generate the bigwig files. Probably a kind of normalized read count, but it technically can be any numeric value. bamCompare in deeptools can combine replicates.

Log in to answer this question.