Normalization and analysis of microarray data was as described previously (28 (link)). For CpG analysis, CpG islands (CGIs) were identified by finding probes with a minimum of 25 bp overlap with CGI found ±1 kb within genes using Galaxy software. CGI positions were taken from the University of California Santa Cruz (UCSC) table browser.
Enrichment analysis for ChIP and run-on data was performed for probes ±1 kb from transcription start sites (TSS) or transcription end site (TES). The smoothed conditional mean plots were generated using the R package ggplot2 and the geom_smooth function.
The following mm9 coordinates were used for quantification of ChIP enrichment; non-expressed 3′ Hoax(Hoxa1 to Hoxa7) genes chr6:52,101,011–52,172,728, expressed 5′ Hoxa genes (Hoxa9-Hoxa13) chr6:52,171,296–52,211,033, 3′ non-expressed Hoxd (Hoxd1 to Hoxd9 genes) chr2:74,534,258–74,606,421, expressed 5′ Hoxd genes (Hoxd9-Hoxd13) chr2:74,484,916–74,537,448. To test the significance of differential ChIP enrichment at genomic regions a Wilcoxon rank-sum test was performed with a correction for multiple testing (Holm method) using the R statistical program.
For sequential ChIP (SeqChIP), antibodies were covalently coupled to Dynabeads with antibody coupling kit, (Invitrogen Cat. 14311D), using the manufacture's protocol. The first ChIP was eluted with 10 mM DTT and the elute was diluted 30 times with Radio-Immunoprecipitation Assay (RIPA) (50 mMTris, pH 7.5, 150 mMNaCl, 1% IGEPAL CA-630, 0.5% deoxycholate) buffer before continuing with the second ChIP. Primers used for ChIP qPCR are given in Supplementary Table S1.