ChIP-seq libraries were prepared using standard protocols provided by Illumina (www.illumina.com). Libraries were sequenced (Solexa, Illumina) in two successive runs and output from each sample was used for analysis. The resulting sequence outputs (bases 2-42) were aligned to the reference mouse genome (version mm9) using bowtie 0.12.7 (19 (link)).
The peaks were filtered for P < 10−20, next, peak calling and annotation were done by using Model-based Analysis of ChIP-Seq (MACS, version 1.4) (20 (link)) and peakAnnotator (version 1.4) (21 (link)), respectively. The published ChIP-seq data from Gene Expression Omnibus (GEO) database were acquired to analyze the binding targets of H3K9me2 and H3K9me3 (accession number GSE82144) (22 (link)), RELA (GSE36099) (23 (link)), and KDM4A and KDM4C in mouse acute myeloid leukemia (AML) cells (GSE81300) (24 (link)).