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 25 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 • 133 views
ADD COMMENT
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

ADD REPLY
0
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 20 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.

ADD COMMENT

Login before adding your answer.

Traffic: 323 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