Heatmap for ChIP-seq data
1
0
Entering edit mode
6 months ago
Marco Pannone ▴ 240

Hello everybody!

I have some data from ChIP-seq analysis, I have performed peak calling and I have annotated the peaks to the nearby genes. So basically I have a .bed file with unique peaks, genomic coordinates of the peaks, and a final column with the gene name of the nearby annotated gene for each peak.

I look for a way to produce a heatmap (and a profile plot on top of it) where I can plot all the peaks accordingly to their distance to the TSS, while having the gene names on each row of the heatmap.

I have checked some functions such as plotHeatmap from deepTools, but it does not seem to provide me an option for keeping the gene names as row names of the heatmap when plotting.

Any help would be very appreciated!

Thanks in advance :)

heatmap chip-seq • 301 views
ADD COMMENT
2
Entering edit mode
6 months ago
Papyrus ★ 2.1k

I think you can use EnrichedHeatmap (R package) with the show_row_names = T argument, like in the example shown here by the author: example.

ADD COMMENT
0
Entering edit mode

Thanks for the suggestion! I did not know the existence of this package and I will for sure give it a try! :)

ADD REPLY

Login before adding your answer.

Traffic: 2186 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6