Best Peak Caller For Histone Modifications (Mainly H3K4Me3) Chip-Seq Datasets?
5
7
Entering edit mode
13.6 years ago

I've heard of MACS for peak/region calling of histone modification chip-seq datasets (mainly H3K4me3), but I am wondering if there is any other tool I should consider as well.

chip-seq peak-calling • 12k views
ADD COMMENT
0
Entering edit mode

This BioStar thread is similar.

ADD REPLY
4
Entering edit mode
13.6 years ago

SICER is especially designed for histone marks. The recently published PeakRanger also does diffuse profiles (like histone marks) and seems interesting.

ADD COMMENT
0
Entering edit mode

New version of SICER: https://github.com/endrebak/epic

ADD REPLY
3
Entering edit mode
13.6 years ago
Dataminer ★ 2.8k

MACS is best although you can try PICS as well.

Reason for using PICS: Two binding sites are separated into two disjoint regions by PICS, whereas MACS combines these two sites into a single region. This makes PICS particularly attractive for subsequent motif based analysis.

:)

ADD COMMENT
2
Entering edit mode

"Best" in terms of "user experience": having first-hand experience with your own data and liking one as the best in comparison to others.

ADD REPLY
1
Entering edit mode

Personal experience shows that MACS is not good for histones. But excellent for TFs, although not terribly sensitive.

ADD REPLY
0
Entering edit mode

MACS is best, based on what criterion?

ADD REPLY
0
Entering edit mode

Personally no algorithm is based, it was my personal biased opinion. Although I also gave a second option of PICS and guess you are big fan of SICER .. :)

ADD REPLY
0
Entering edit mode

Personally no algorithm is best, it was my personal biased opinion. I also gave a second option of PICS and I guess you are a big fan of SICER .. :)

ADD REPLY
0
Entering edit mode

Haven't used SICER much really, but it does have an underlying model of histone mark enrichment (based on nucleosome spacing). CCAT is the package I've used the most for histone modifications

ADD REPLY
0
Entering edit mode

@avilella: Well said

ADD REPLY
0
Entering edit mode

Sorry for being grumpy ... I guess what I was trying to say was that the answer did not address the issue of how MACS and PICS are addressing histone modification profiles specifically. This is a scenario where motif based analysis (as you suggested for PICS) is not very relevant.

ADD REPLY
0
Entering edit mode

Its ok Mikael, once in a while we all are grumpy :) To be honest no algorithm is perfect (as far as I know) they all have drawbacks. I agree that SICER could be a better choice depending upon what kind of biological question you want to address after marking the histone modification profiles. And, I never said you were wrong and I was right.

Friends?

ADD REPLY
2
Entering edit mode
13.6 years ago
Ian 6.1k

I have not had chance to look at this yet, but RSEG looks very promising! I expect this to be good as Andrew Smith has historically produced very useful software tools, e.g. CREAD and DME.

We have also had a good experience using SICER.

Although i routinely use MACS for ChIP-seq i would not recommend its use with histone data.

ADD COMMENT
1
Entering edit mode
13.6 years ago
Benm ▴ 710

Maybe you can find the answer from this article: Shirley Pepke, Barbara Wold & Ali Mortazavi, Computation for ChIP-seq and RNA-seq studies. Nature Methods 6, S22 - S32 (2009).

And this question had been discussed in SEQanswers before, just try to google "ChIP-Seq Challenge".

ADD COMMENT
0
Entering edit mode

It would be great to be able to pull SEQanswers content in the "Related questions" area in Biostar.

ADD REPLY
1
Entering edit mode
13.6 years ago

I have not tried this myself, but you can have a look at the ChipSeqR Bioconductor package. There is a fairly extensive section on testing this method in the associated paper by Humburg et al. that may get to the issue of "best" peak caller for nucleosomes/histon modifications. This paper also cites two other papers by Albert et al and Kharchenko et al that are probably worth having a look at as well.

ADD COMMENT

Login before adding your answer.

Traffic: 1506 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6