How to annotate a region in Heatmap?

[li1311139481]

Welcome to deepTools GitHub repository! Before opening the issue please check
that the following requirements are met :

• Search whether this issue (or a similar issue) has been solved before
  using the search tab above. Link the previous issue if appropriate below.

• Paste your deepTools version (deeptools --version) and your python
  version (python --version) below.
  deeptools 3.5.5
  Python 3.10.12

• Paste the full deepTools command that produces the issue below
  (ignore if you simply spotted the issue in the code/documentation).

/cluster/facility/hlhuang/miniconda3/envs/pip/bin/computeMatrix reference-point --referencePoint center -p ${THREADS} -b 3000 -a 3000 -R "/cluster/facility/hlhuang/zifeng001/1file/ATACregion/tss.bed" -S bigwig/*.bw --skipZeros -o tss/TSS.gz  
/cluster/facility/hlhuang/miniconda3/envs/pip/bin/plotHeatmap -m tss/TSS.gz  -out tss/TSS_Heatmap.pdf --plotFileFormat pdf  --dpi 720

• Paste the output printed on screen from the command that produces the issue
  below (ignore if you simply spotted the issue in the code/documentation).
  

Question:
I get the result. but i want anno the heatmap using some gene. As you can see, my input is differential peaks or Tss.bed, and i annotate the peaks to gene, Or I want to know which gene's TSS is different in different samples. I know featureAlignedHeatmap function in ChIPpeakAnno package in R. But I don't think it's clear enough, especially when it's presented as anno_block.

I want annotate the region in heatmap like this

Any help would be appreciated