How can I make RSEM to align ChIPseq signal to the nearest gene?
1
0
Entering edit mode
4.1 years ago
moxu ▴ 500

I guess I can make a GTF by myself, and let the GTF associate genomic coordinates with genes. But assigning coordinates to closest genes is not a simple task. Overlapping gene regions can be problematic. Is there such a GTF available or is there an option to turn on/off with rsem-prepare-reference/rsem-calculate-expression/etc.?

ChIP-Seq next-gen sequencing alignment gene • 887 views
0
Entering edit mode

What is your end goal with all of this?

0
Entering edit mode

Good question.

Assign ChIP signals to genes.

1
Entering edit mode
4.1 years ago

Exactly what kind of signal? I presume you either want to assign peaks to genes (bedtools closest) or look at signal distributions near genes (see computeMatrix and plotHeatmap from deepTools).

0
Entering edit mode

Could be Pol2, Histone modifications, ATAC, etc.

For Pol2, usually we would be interested in +-5k from TSS; for histone modifications, ATAC, double strand breaks, ... , often you want to check the whole genome regardless how far/close the signal is to a gene.

0
Entering edit mode

Then you want computeMatrix from deepTools followed by either plotProfile or plotHeatmap.