Biostar Beta. Not for public use.
Question: Help with HOMER motif profile plots
0
Entering edit mode

Hi,

I'm trying to understand profile plots for motifs. Homer has a function annotatePeaks.pl) to generate motif density profiles that looks like following (lifted from here)

enter image description here

While I understand what these profile plots convey for ChIP-seq data (e.g. those generated by deepTools plotProfile), but I'm not sure what exactly is the y-axis in above plot.

So the question is, is it just the probability of finding given motif at each position (as MEME does here) in a set of given sequences or is it the actual count normalized to something?

I didn't find any explanation on homer website. It would be great if anyone can explain what exactly is going on in above plot.

Thank you.

ADD COMMENTlink 12 months ago venu 6.2k
Entering edit mode
1

Convert a BED file of motif positions to a bigWig with a 1 where there's a motif and a 0 otherwise and it's the same thing you'd then get with plotProfile I think.

ADD REPLYlink 12 months ago
Devon Ryan
90k
Entering edit mode
0

I created a bigiwg ( using bedtools genomecov + kent utils) and plotted with plotProfile. The motif I'm interested in is has very sparse hits, applying loess fit on plotProfile output data is more/less looking closer to the above plot. I think I understood what's happening inside.

ADD REPLYlink 12 months ago
venu
6.2k
3
Entering edit mode

For future googlers, here is the response from lead author of HOMER, Chris Benner

The motif plots are normalized by taking the actual motif totals, then normalizing by the bin size used in the histogram (so that the motif densities are comparable across different bin sizes) and the total number of peaks/regions used to create the histogram (so that the motif densities are independent of the number of regions used to search).

ADD COMMENTlink 12 months ago venu 6.2k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0