Question: Computing the significance of overlapping ChIP-seq peaks
0
gravatar for Marcel
13 months ago by
Marcel0
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?

Thanks.

chip-seq • 722 views
ADD COMMENTlink modified 13 months ago by colin.kern930 • written 13 months ago by Marcel0

Have you checked these previous threads?

ADD REPLYlink modified 13 months ago • written 13 months ago by igor11k

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 13 months ago by Marcel0
1
gravatar for colin.kern
13 months ago by
colin.kern930
United States
colin.kern930 wrote:

There's a tool called IntervalStats that does this. Here is their paper: https://academic.oup.com/bioinformatics/article/28/5/607/247792

ADD COMMENTlink written 13 months ago by colin.kern930

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 13 months ago by Marcel0
Please log in to add an answer.

Help
Access

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