Related papers to read:
Model-based Analysis of ChIP-Seq (MACS)
Curr Protoc Bioinformatics. Author manuscript; available in PMC 2012 June 1.MACS empirically calculates FDR based on the number of peaks from control over ChIP that are called at the same p-value cutoff. Therefore if no control data is available, the FDR column does not exist in the output tabular file. Technically, MACS can also be applied to identify differential peaks between two conditions by treating one of the samples as the control. However, calculated FDR value should be ignored, as peaks from either sample are likely to be biologically meaningful in this case.
http://www.nature.com/nprot/journal/v7/n9/full/nprot.2012.101.htmlthe warning message ‘unbalanced reads between treatment and control’ means that the FDR of the resulting peaks will be overestimated when the control sample has more reads and will be underestimated when the ChIP-seq sample is sequenced more deeply.
MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets
|1||Supplementary Table 1. Enrichment of cell-type differentially expressed genes in genes near differential binding peaks defined by MAnorm, MACS and ChIPdiff. To compare the enrichment scores, we selected the same number of target genes associated with top differential binding regions identified by MAnorm with those identified by other two methods.|
|2||Supplementary Table 1A. Enrichment of genes more highly expressed in H1 ES cells (as compared to K562) in genes near H1 ES enriched peaks (as compared to K562) defined by MAnorm, MACS and ChIPdiff|
|3||H3K27ac H1 ES-enriched target genes||Number of Genes||Overlap with H1 ES up-regulated genes||Enrichment Score|
|6||MAnorm (top 1467 genes; same number of genes as identified by ChIPdiff with default settings)||1467||941||3.45|
|8||MAnorm (top 1589 genes; same number of genes as identified by MACS with P<1e-6)||1589||987||3.34|
|11||Supplementary Table 1B. Enrichment of K562 higher expressed genes (as compared to H1 ES) in genes near K562 enriched peaks (as compared to H1 ES) defined by MAnorm, MACS and ChIPdiff|
|12||H3K27ac K562-enriched target genes||Number of Genes||Overlap with K562 up-regulated genes||Enrichment Score|
|15||ChIPdiff (confidence threshold=0.9999999999)||2291||697||2.55|
|16||MAnorm (top 2291 genes; same number of genes as identified by ChIPdiff with confidence threshold=0.9999999999)||2291||820||3.00|
|19||MAnorm (top 1556 genes; same number of genes as identified by MACS with P<1e-6))||1556||644||3.47|
diffReps: Detecting Differential Chromatin Modification Sites from ChIP-seq Data with Biological Replicates
…shows that diffReps is the most sensitive method among all the methods compared, followed by edgeR, DESeq, ChIPDiff and, lastly, CCAT+DESeq. At each cutoff, diffReps (negative binomial test) typically detects a few thousands more differential sites than the secondly ranked method, edgeR.
The tool is in Perl (https://code.google.com/p/diffreps/)
Detecting differential binding of transcription factors with ChIP-seq, Bioinformatics, 28, 121-122, doi: 10.1093/bioinformatics/btr605
DiffBind : differential binding analysis of ChIP-Seq peak datahttp://bioconductor.org/packages/2.12/bioc/vignettes/DiffBind/inst/doc/DiffBind.pdf