Question: Chip-Seq - Peak caller for broad enrichment of histone modifications?
gravatar for ej
5.1 years ago by
European Union
ej60 wrote:

 I'm analyzing Chip-Seq data and am interested in histone modifications that have very broad peaks such as H3K27me3 and H3K36me3. I have tried peak calling using Homer but it does not seem to work well. Is there a better peak caller for these modifications that have broad enrichment across the genome? Thanks

chip-seq • 6.5k views
ADD COMMENTlink modified 4.0 years ago by Ming Tang2.6k • written 5.1 years ago by ej60
gravatar for Manvendra Singh
5.1 years ago by
Manvendra Singh2.1k
Berlin, Germany
Manvendra Singh2.1k wrote:

Why aren't you using macs for calling peaks. Its very good for Histone peaks , including broad ones

ADD COMMENTlink written 5.1 years ago by Manvendra Singh2.1k
gravatar for Ming Tang
4.0 years ago by
Ming Tang2.6k
Houston/MD Anderson Cancer Center
Ming Tang2.6k wrote:

I have some collections here

ADD COMMENTlink written 4.0 years ago by Ming Tang2.6k
gravatar for Gary
5.1 years ago by
Taiwan/Taichung/China Medical University Hospital
Gary480 wrote:

SICER is a specialized peak caller for broad enriched histone modification (Zang, et al 2009). A comment in OMICtools also said that SICER is better than MACS for broad enrichment antibodies. However, I use MACS 1.4, the most stable version, to call H3K27me3 peaks with parameters suggested by Dr. Liu as below (Feng, et al 2012). For our ChIP-Seq data, MACS 1.4 is good enough. I used MACS, but not SICER, because we would like to use the same peak calling tool to call all ChIP-Seq data, including broad and sharp enriched peaks. In addition, MACS 2 has added broad enriched parameters (--broad) that could be helpful for your purpose.


Zang, et al. 2009. A clustering approach for identification of enriched domains from histone modification ChIP-Seq data. Bioinformatics

A comment for SICER in OMICtools

Ideal for broad peaks like H3K27me3 and H3K36me3. Where MACS will truncate peaks in smaller fragments, generating dozens of small peaks, the strenght of SICER is that you can take gaps into account. You just have to perform various tests to choose a good combination between gap size and window size. Typically, I used a window of 500nt and a gap of 3500nt for a good compromise between sensitivity and specificity.

Posted by Fabien Pichon on Oct 17, 2014

Feng JX, Liu T, Qin B, Zhang Y, Liu XS. 2012. Identifying ChIP-seq enrichment using MACS. Nature Protocols 7:1728-1740.

> macs14 -t BROAD_GM12878_H3K4me3.bam –c BROAD_GM12878_H3K4me3_Control.bam -g hs -n BROAD_GM12878_H3K4me3 --nomodel --shiftsize 73 -B -S --call-subpeaks

ADD COMMENTlink modified 5.1 years ago • written 5.1 years ago by Gary480

I have tested both, and even though I like MACS a lot, easy to use and had good results in the past, I found that for broad regions SICER gave me some better results after setting the Window and gap size appropriately. By good results I mean more enriched regions matching expectations and what could be observed visually from the normalized coverage.    

ADD REPLYlink written 5.1 years ago by A. Domingues2.2k

New version of SICER:

ADD REPLYlink written 4.0 years ago by endrebak810
gravatar for endrebak
4.0 years ago by
endrebak810 wrote:

I am biased, but is good.

Also, if you prefer R and know some linear modelling csaw is an option:

ADD COMMENTlink written 4.0 years ago by endrebak810
gravatar for Ido Tamir
5.1 years ago by
Ido Tamir5.0k
Ido Tamir5.0k wrote:

You could try ZINBA "ZINBA integrates local covariates with DNA-seq data to identify broad and narrow regions of enrichment, even within amplified genomic regions" It also uses local mappability information. I never used it because it has no CLI (only R).

ADD COMMENTlink written 5.1 years ago by Ido Tamir5.0k
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: 2129 users visited in the last hour