Question: Computing the significance of overlapping ChIP-seq peaks
gravatar for Marcel
16 months ago by
Marcel0 wrote:

Hello everyone,

I have ChIP-seq data of stalled Pol2 (13.000 peaks) and a protein that mostly binds in promoter regions (250 peaks). Now I want to know whether there is a significant overlap between the two.

What I did try so far, is using "bedtools fisher" with the two datasets. However, as far as I understand, this will not take care of the fact, that both datasets mostly bind in promoter regions. Of course, it's telling me that the overlap is extremely significant.

In other words: I want to know whether my protein of interest binds significantly more often than expected to a region where Pol2 is sitting, considering that Pol2 already occupies a large fraction of promoters with its 13.000 peaks.

Any ideas on how to do this properly?


chip-seq • 864 views
ADD COMMENTlink modified 16 months ago by colin.kern940 • written 16 months ago by Marcel0

Have you checked these previous threads?

ADD REPLYlink modified 16 months ago • written 16 months ago by igor12k

Sorry for the late reply. I've had a look at these threats and they mention IntervalStats (which was also suggested in the other reply to my thread) which appears to fit my purpose.. however, I can't get it running. Other than that, I'll have a look at Cooccur, which was also mentioned. Thanks for your reply!

ADD REPLYlink written 16 months ago by Marcel0
gravatar for colin.kern
16 months ago by
United States
colin.kern940 wrote:

There's a tool called IntervalStats that does this. Here is their paper:

ADD COMMENTlink written 16 months ago by colin.kern940

Sorry for the late reply and thanks for the suggestion. This tool indeed seems to be doing exactly what I want. However, it's throwing me errors saying "Malformed BED line".. but my .bed files are allright. I've been using them for other analyses already.

ADD REPLYlink written 16 months ago by Marcel0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2612 users visited in the last hour