Question: Adjusting Peak Calling For Broad Enrichment Of Histone Modifications.
8
gravatar for Dave Gerrard
8.4 years ago by
Dave Gerrard190
Manchester
Dave Gerrard190 wrote:

We are working on some ChIP-seq data targeting histone modifications which may be widely dispersed across the genome. We expect the enrichment to be much less peaky and more 'foothills vs. plains'. Does anyone have any tips or tricks for getting the most out of peak calling software in this situation or any alternative strategies? We are using Macs14 and QuEST but are open to trying other methods.

chip-seq peak-calling • 4.6k views
ADD COMMENTlink written 8.4 years ago by Dave Gerrard190
7
gravatar for Alastair Kerr
8.4 years ago by
Alastair Kerr5.2k
The University of Edinburgh, UK
Alastair Kerr5.2k wrote:

We have looked at this type of data in the following papers:

  1. DOI: 10.1038/nature08924 link
  2. DOI: 10.1016/j.molcel.2010.01.030 link
  3. DOI: 10.1371/journal.pgen.1001134 link

We ended up writing our own set of perl scripts and changed variables based on the type of data examined (full details and rationales in the methods section of these papers)

Our 'peak finder' defines a peak as an area where values are above the height threshold, less than the specified gap apart and the whole area falls within the length thresholds. It outputs the peaks as a Bed file with statistics. I'll put our peak finding script here but please be warned that is was not intended for release and it is still largely undocumented. If there is a need, we can tidy it up and put it on the galaxy toolshed.

ADD COMMENTlink written 8.4 years ago by Alastair Kerr5.2k
1

Can you please update the link to your perl script!

ADD REPLYlink written 2.3 years ago by Deepak Tanwar3.9k

I removed this script as there now exists a range of decent pack callers.. See this wikipedia page

ADD REPLYlink written 2.3 years ago by Alastair Kerr5.2k
2
gravatar for Brad Chapman
8.4 years ago by
Brad Chapman9.4k
Boston, MA
Brad Chapman9.4k wrote:

Nice discussion topic -- I've also struggled with large histone modification peaks. One approach which we had some success with is using MACS 1.4 with the call-subpeaks option to subdivide the larger peaks using PeakSplitter:

http://www.ebi.ac.uk/bertone/software.html

We then overlapped these with nucleosome position calls from NPS:

http://liulab.dfci.harvard.edu/NPS/

This was helpful to get a more refined set of reference nucleosome regions that could then be used for comparisons between experiments. Here's the python code used to combine the NPS and MACS calls:

https://github.com/chapmanb/mgh_projects/blob/master/cy_histone_chipseq/merge_nps_macs.py

ADD COMMENTlink modified 8.4 years ago • written 8.4 years ago by Brad Chapman9.4k
2
gravatar for Mikael Huss
8.4 years ago by
Mikael Huss4.6k
Stockholm
Mikael Huss4.6k wrote:

I've used SICER and CCAT with reasonable success for this kind of problem. The former is specifically designed for diffuse enrichment regions (but not TFs) and the latter has a "peak mode" for TFs and similar cases and a "region mode" for cases where you expect a more spread-out enrichment.

ADD COMMENTlink written 8.4 years ago by Mikael Huss4.6k
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: 1980 users visited in the last hour