NGS ChIP-seq / Find peaks using MACS, treatment only

Description

This tool will identify statistically significantly enriched genomic regions in sequencing data from a ChIP-seq experiment. The analysis is performed on one treatment samples alone, without taking into account samples from a control experiment.

Parameters


Details

This tool leverages the functionality of the Model-Based Analysis of ChIP-Seq (MACS) algorithm, which can be controlled by a number of parameters. While file format, genome and read length are defined by the experiment parameters chosen at the time of the experiment the band with, p-value cutoff, whether to use a peak model and the M-fold cutoffs should be titrated to obtain optimal results.

A good starting point for optimizing the band with, which determines the window size used by MACS to scan the genomic regions, is to input a value roughly half the size of the average fragment length of the DNA.

Lowering the p/value cutoff will increase the specificity of the found peak regions but may result in a restrictively low number of peaks. Conversely, applying a less stringent cutoff will improve sensitivity but may yield too high proportion of false positives.

By default MACS tries to build a peek model, based on the band with and the M-fold cutoff limits, that is later used to determine statistically significantly enriched genomic regions. A small value for the lower cutoff limit will provide a greater number of peak regions for model building but could potentially include spurious data into the model that will adversely affect the peak finding results. Setting it too high improves the quality of the model but could have the downside of excluding so many candidate model regions that model building fails. If a good compromise setting is difficult to obtain there is the option to turn off the model building altogether, resulting in MACS simply using a shift-size of 100 bp when searching for peaks.

Output

The analysis output consists of the following:


References

This tool uses the MACS package for peak detection. Please read the following article for more detailed information.

Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. (2008) Model-based Analysis of ChIP-Seq (MACS), Genome Biology, 2008;9(9):R137.