Redesigned and improved ChIP-seq broad peak calling tool SICER
Chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-seq) can be used to map binding sites of a protein of interest in the genome. Histone modifications usually occupy broad chromatin domains and result in diffuse patterns in ChIP-seq data that make it difficult to identify signal enrichment. SICER, a spatial clustering approach for the identification of ChIP-enriched regions, was developed for calling broad peaks from ChIP-seq data.
Usability of the original SICER software has been affected by increased throughputs of ChIP-seq experiments over the years. We now present SICER2 a more user-friendly version of SICER that has been redisgned and streamlined to handle large ChIP-seq data sets. This new Python package supports multiple job submissions on cluster systems and parallel processing on multicore architectures.
For more information about the original SICER algorithm, please see
“A clustering approach for identification of enriched domains from histone modification ChIP-Seq data" Chongzhi Zang, Dustin E. Schones, Chen Zeng, Kairong Cui, Keji Zhao, and Weiqun Peng, Bioinformatics 25, 1952 - 1958 (2009)
In addition, we present an alternative algorithm for identification of broad domains from ChIP-seq data called RECOGNICER.
(More detail about RECOGNICER TBD)
Easiest way to install SICER2 is through
pip. Simply open the terminal and type
pip install SICER2.
Unlike the original version of SICER, SICER2 runs in Python 3. Please use Python 3 to install and run SICER 2.0.
Numpy and Scipy are required to run SICER2. Please have these installed before installing SICER2.
This can be done by simply typing
pip install numpy scipy under command line (if python2.7 is your default python version, use
C compiler is required to compile C codes that are part of the SICER2 package. This also means that python header files (e.g. Python.h) are needed. For Linux users, make sure to have python-dev installed. For Mac OS X users, it is recommended that you install Xcode.
Lastly, if you would like to directly pass BAM files as input files for SICER2, you need to have bedtools installed. Please refer to this link for more details on installing bedtools. This is not required if you will intend to only pass BED files as input files.
For local installation, the source distribution file is available at Zang Lab website (link)
The terminal command to run SICER is
sicer. The command to run RECOGNICER is
sicer command and
recognicer command take several flag-based arguments as inputs. Two arguments are required: (1) name/path of treatment file, (2) species of reference genome. Arguments for each command are explained below.
The file must either be in BED or BAM format (note that you need bedtools installed to directly enter BAM files). The file name can either the relative path or the absolute path of the file.
Like the treatment file, control file must be in BED or BAM format and can be the relative path or the absolute path of the file. However, control library input is optional.
The number of copies of indentical reads allowed in a library. Default value is 1
Resolution of SICER. Default value is 200 (bp)
The amount of shift from the beginning of a read to the center of the DNA fragment represented by the read. Default value is 150 (bp).
Effective genome as fraction of the genome size. Default value is 0.74.
Remove all islands with an false-discovery-rate below cutoff. Default value is 0.01.
The minimum length of a "gap" such that neighboring window is an "island." Please note that this value must be a multiple of the window size. Default value is 600 (bp).
E-value. Requires user input when no control library is provided. Default value is 1000
Path of the directory in which results will be stored. Default output directory is the current working directory.
The number of CPU cores SICER program will use when executing multi-processing tasks. Optimal number of cores is the species' number of chromosomes. Default value is the maximum number of cores avaiable in the system.
Significant Reads: Type "--significant_reads" flag to have SICER produce a BED file of treatment reads filtered by significant islands and WIG file of filtered reads binned into windows.
All of the arguments for RECOGNICER are identical to those of SICER except for
Instead of these two arguments, RECOGNICER has two arguments called
The number of windows in one graining unit. Default value is 3.
The minimum number of positive elements in the graining unit to call the unit positive. Default value is 2.
Using SICER2 for differential peak calling
The commands for differential peak calling are
The arguments for both SICER and RECOGNICER differential peak calling are identical to those of the regular peak callings except for the following arguments specified below.
Also, differential peak calling has one additional argument called
Two files must be given as input. The first file must be the knockout (KO) file and the second file must be the wild-type (WT) file. Both files must either be in BED or BAM format.
While optional, two files must be given as input if you decide to provide the input. The first file must be the control library corresponding to the knockout (KO) treatment file and the second file must be the control library corresponding to the wild-type (WT) treatment file. Both files must either be in BED or BAM format.
Cutoff for identification of significant changes been wild-type library and knockout library. Default value is 0.01.
- Calling SICER with a control library. Default parameters are explicitly entered for the sake of demonstration.
sicer -t treatment.bed -c control.bed -s hg38 -w 200 -rt 1 -f 150 -egf 0.74 -fdr 0.01 -g 600 -e 1000
- Calling SICER without a control library
sicer -t treatment.bed -s hg38
- Calling SICER with control libraries for differential peak calling.
sicer_df -t treatment1.bed treatment2.bed -c control1.bed control2.bed -s hg38
- Calling SICER without control libraries for differential peak calling.
sicer_df -t treatment1.bed treatment2.bed -s hg38