how to select specific genes from heatmap of differentially accessible regions of ATAC-seq in DiffBind
1
0
Entering edit mode
Laura • 0
@5cec3508
Last seen 24 days ago
Austria

Hello everybody, I would like to do a heatmap of selected genes for my ATAC-seq data. When I use DiffBind and I do dba.PlotHeatmap of my differentially accessible regions I get a full heatmap of all the regions and all my samples. Is it possible to show only the heatmap for specific genes maybe by specifing the region around TSS? (for e.g. -1000bp,+500bp) Thank you

DiffBind • 131 views
0
Entering edit mode

I have the same question, have you found a solution yet? if you havent, let me know when you do. Thanks. Regards, Team Smokekitchen

0
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 19 days ago
CRUK, Cambridge, UK

I suspect the dba.plotProfile() function would be more appropriate for this. Here's a workbook showing how it works:

https://content.cruk.cam.ac.uk/bioinformatics/software/DiffBind/plotProfileDemo.html

The trick would be if the differential ATAC regions are not the same as the gene regions you want. In that case, you can provide a consensus peakset with the areas around the TSSs of interest.