The diverse gene expression programs that allow for differentiation and response to environmental stimuli result from the regulated binding of transcription factors to DNA. The prevalent technique used in chromatin biology for mapping protein–DNA interactions is chromatin immunoprecipitation (ChIP), but little has changed since it was first described 27 years ago (Solomon et al., 1988). Despite recent advances in read-out technologies for ChIP, such as high-throughput sequencing (ChIP-seq), the basic chromatin preparation protocol remains the same and has a number of limitations. For example, sonication is typically used to fragment the chromatin. This however, has been shown to be non-random, with heterochromatic regions showing increased resistance to fragmentation leading to bias in the experiment (Teytelman et al., 2009). In addition, sonication typically produces chromatin fragments between 200 and 500 bp, whereas the footprint of a typical chromatin-associated protein is ∼10-fold smaller, indicating a lack of resolution currently obtained by ChIP-seq. Even extensive sonication only yields fragments with an average length of 200 bp, suggesting that sonication is of limited use in generating high-resolution maps of genome-wide protein binding (Fan et al., 2008). In a previous study, we were interested in how RNA Polymerase II transcribes through nucleosomes at the promoter (Skene et al., 2014). Answering this question required the precise mapping of PolII with respect to the position of nucleosomes, but conventional ChIP-seq that uses sonication to fragment the chromatin, yields fragments approximately twice the size of a nucleosome. Additionally, it has been shown that PolII can crosslink to nearby nucleosomes and therefore mapping the immunoprecipitated DNA fragments from these composite PolII:nucleosome:DNA complexes fails to precisely map the position of PolII on the DNA (Koerber et al., 2009; Skene et al., 2014). By using micrococcal nuclease (MNase) to digest unprotected DNA, we were able to achieve high resolution in a ChIP experiment, mapping the precise location of PolII and chromatin remodelers on the DNA (Figure 1A; a detailed protocol is provided as a Supplementary file 1) (Skene et al., 2014). Optimization of this simple protocol for the high-resolution mapping of protein–DNA interactions has the potential to revolutionize our understanding of genome-wide protein binding. High resolution is especially a requirement at closely spaced transcription factor binding sites, such as locus control regions and super-enhancers, which have been shown to be vital to cell fate decisions and human diseases (Hnisz et al., 2013; Pott and Lieb, 2014).

Figure 1 with 1 supplement see all

Download asset Open asset

To evaluate high-resolution crosslinking ChIP-seq (X-ChIP-seq), we compared it to existing methodologies. Initially, we focused on mapping PolII at the transcriptional start site (TSS) in Drosophila S2 cells, where there are existing data sets at both low and high resolution. We performed high-resolution X-ChIP-seq using the same antibody against total PolII (Rpb3 subunit) as previously used by a conventional sonication ChIP experiment (Figure 1B) (Core et al., 2012). Using this cell line also allowed a comparison with the single base-pair resolution technique that maps the last ribonucleotide incorporated into the nascent RNA chain (3′NT), thereby mapping the exact position of the PolII active site (Weber et al., 2014). By using paired-end sequencing, we can selectively study specific lengths of immunoprecipitated fragments. Analyzing sequenced fragments with lengths 20–70 bp, which more closely represent the footprint of PolII (Samkurashvili and Luse, 1996), avoids the complication of mapping fragments consisting of PolII crosslinked to adjacent nucleosomes. Using this technique, we find that the maximal peak of PolII signal coincides with the position of the polymerase's active site at ∼+35 bp, as measured by 3′NT. This is consistent with evidence suggesting that the vast majority of Drosophila genes have a productively engaged PolII enzyme stalled just downstream of the promoter rather than PolII stably bound at the pre-initiation complex (Core et al., 2012). In contrast, PolII distribution as measured by conventional ChIP with the chromatin fragmented by sonication, shows a distinct distribution at the promoter with a broader peak centered at the TSS with maximal density at −5 bp. This discrepancy likely comes from biases in the probability of sonication breaking the DNA at the nucleosome-depleted region of the promoter, as accessible regions such as DNase I sites and promoters of active genes have been shown to be sonicated at higher probability than inactive genomic regions (Teytelman et al., 2009). Analysis of a published sonicated input chromatin sample indicates a strong sonication bias at the promoter region (Figure 1—figure supplement 1). In contrast, by predominantly fragmenting the chromatin with MNase, it is possible to generate footprints corresponding to nucleosomes and other DNA-bound factors (Henikoff et al., 2011; Skene et al., 2014). Overall, this shows that using a high-resolution ChIP technique to map the protected footprint of PolII achieves comparable resolution to the single base-pair resolution achieved by mapping the position of the active site of PolII via nascent chain mapping. In comparison to conventional ChIP-seq, using high-resolution X-ChIP-seq achieves both higher resolution, as indicated by the width of the ChIP peak and higher accuracy by avoiding sonication bias, as shown by high similarity to 3′NT. Moreover, the depth of sequencing indicates the cost-effectiveness of this high-resolution ChIP approach, with the 3′NT profile based on 150 million reads (Weber et al., 2014), whereas our method required only 7 million paired-end reads with a fragment length of 20–70 bp. For comparison, the PolII profile generated by conventional ChIP was based on 13 million mapped reads (Core et al., 2012).

A limitation of high-resolution X-ChIP-seq is that a minority of the immunoprecipitated fragments represent the footprint of PolII on DNA, likely as a consequence of formaldehyde readily forming protein–protein crosslinks generating complexes such as PolII crosslinked to nucleosomes (Koerber et al., 2009; Skene et al., 2014). In our previous study, mapping murine PolII, only 10% of the fragments were 20–70 bp in length and less than 3% were under 50 bp (Figure 1C) (Skene et al., 2014). Therefore, to improve the cost-effectiveness of this technique and make it more applicable to transcription factors, which typically have a <50-bp footprint, we have further optimized the method to enrich for short fragments prior to sequencing. Previously, Agencourt AMpure beads have been used to select for short fragments prior to linker ligation (Orsi et al., 2015). In agreement, initial attempts indicated that Agencourt AMpure beads could enrich for DNA fragments below 100 bp from a complex mixture, but were unable to selectively purify fragments of ∼50 bp. However, by adjusting the volumetric ratio of beads to DNA, we could reproducibly control the selection within the 100–200 bp range with a ratio of 1.1×, leaving fragments of ∼170 bp in the unbound fraction (Figure 2A). Given that the ligation of the Illumina adapters to the immunoprecipitated DNA adds ∼125 bp, by using this ratio of AMpure beads, we could selectively enrich for ligated products containing short inserts (Figure 2B). Using this approach on input DNA from a MNase ChIP experiment, where the vast majority of the DNA fragments are from mono-nucleosomes, we find a 25-fold enrichment of fragments below 50 bp (Figure 2C). Therefore, combining this modification to the existing library preparation protocol with the MNase X-ChIP approach yields cost-effective high-resolution data.

Figure 2

Download asset Open asset

To assess the resolution of this method, we chose the well-characterized transcription factor CCCTC-binding factor (CTCF). We performed high-resolution X-ChIP-seq in K562 cells and analyzed 20–50 bp fragments, and compared this to conventional ChIP as used in the ENCODE project (Figure 3A). To avoid the complexities of peak-calling algorithms that might be biased for differing data types, we used an unbiased approach of centering the data on CTCF motifs that were found within DNase I sites and therefore most likely bound by CTCF. High-resolution X-ChIP-seq of CTCF yielded a more focused distribution of reads centered over the CTCF motif. To quantify this, we measured the width of the ChIP peak at its half-height for each individual CTCF site (Figure 3B). A conventional ChIP approach using sonication gave a half-height width of 200 bp. In contrast, analysis of 20–50 bp fragments from MNase ChIP gave much higher resolution, with a half-height width of only 50 bp, suggesting that genome-wide MNase is chewing back to a minimal footprint of CTCF bound to the DNA. By analyzing different ranges of fragment lengths, it was possible to see that shorter fragments gave the highest resolution and smallest range in peak widths (Figure 3B and Figure 3—figure supplement 1).

Figure 3 with 2 supplements see all

Download asset Open asset

We also compared our CTCF high-resolution X-ChIP-seq results to CTCF profiles obtained using ChIP-exo, which is based upon the sonication of crosslinked chromatin, followed by exonuclease digestion of the immunoprecipitated complexes (Rhee and Pugh, 2011). In ChIP-exo, sequential ligation reactions allow the demarcation of 5′ and 3′ ends and bioinformatic analysis is used to identify ‘peak pairs’ that flank the transcription factor binding site. Figure 3C shows profiles based on ENCODE X-ChIP-seq data as processed by the ENCODE uniform processing pipelines and downloaded as ‘raw signal’, our high-resolution X-ChIP-seq stacked read data, and raw ChIP-exo data around a representative CTCF motif. Due to the low amounts of noise, high-resolution X-ChIP-seq is amenable to a very simple thresholding algorithm to identify peaks (Kasinathan et al., 2014), requiring only 13 million paired-end reads to obtain a crisp peak feature. This is in contrast to ChIP-exo, where more complex analysis is required, including dedicated software (Rhee and Pugh, 2011; Starick et al., 2015). Based on 82 million reads, the ChIP-exo raw data show significant signal at a distance from the CTCF motif. This might be a consequence of immunoprecipitating sonicated 200–300 bp chromatin fragments containing more than one protein, which would block the subsequent exonuclease cleavage. However, by using MNase to fragment the chromatin, which has both endo- and exo-nuclease activity, high-resolution X-ChIP-seq should be able to discriminate between nearby proteins. An additional limitation of ChIP-exo is that the input chromatin is not subjected to the same exonuclease treatment and therefore subsequent analyses cannot be normalized to input. With high-resolution X-ChIP-seq, however, all the steps in chromatin fragmentation are prepared prior to immunoprecipitation, thereby allowing input normalization. Moreover, high-resolution X-ChIP-seq requires only minimal alteration to the existing conventional ChIP workflow and library preparation, whereas other techniques require more extensive changes (Starick et al., 2015). To more directly compare to ChIP-exo, we plotted the end positions of each of our 20–50 bp paired-end reads (Figure 3D). We find two predominant sharp peaks on either side of the 19-bp CTCF motif that are separated by 19 bp, indicating that on average for each of our immunoprecipitated fragments, MNase has chewed back to one side of the minimal sequence motif. In contrast, the signal for ChIP-exo is relatively diffuse when centered around the CTCF motif, with an average distance of 52 bp between peak pairs for the peak-called sites (Rhee and Pugh, 2011). DNase I footprinting is often used to generate maps of global transcription factor binding at nucleotide resolution, with the drawback that the technique is not targeted to a specific transcription factor (Hesselberth et al., 2009; Neph et al., 2012). By comparing the ends of the DNA fragments released by DNase I footprinting and that of high-resolution X-ChIP-seq, we find that both techniques identify protected fragments with ends separated by 19 bp at the 19 bp consensus CTCF motifs (Figure 3—figure supplement 2). This therefore suggests that high-resolution X-ChIP-seq can achieve single nucleotide resolution, and by using immunoprecipitation, has the advantage that it can be used to interrogate individual transcription factors.

By harnessing the endo- and exo-nuclease activity of MNase to fragment chromatin, high-resolution X-ChIP-seq has key advantages over conventional ChIP-seq and ChIP-exo in terms of the resolution obtained (Figure 4). Overall, the combination of the improvements to enrich for short immunoprecipitated fragments and the unparalleled ChIP resolution for PolII and transcription factor binding indicate that high-resolution X-ChIP-seq is a cost-effective and simple approach that easily fits within existing ChIP-seq pipelines for determining precise genome-wide maps of protein–DNA binding.

Figure 4

Download asset Open asset

1. 1. Skene PJ
   2. Henikoff S

   (2015) A simple method for generating high-resolution maps of genome wide protein binding

   Publicly available at the NCBI Gene Expression Omnibus (Accession no: GSE67454).

   http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67454
2. 1. Skene PJ
   2. Hernandez AE
   3. Groudine M
   4. Henikoff S

   (2014) The nucleosomal barrier to promoter escape by RNA Polymerase II is overcome by the chromatin remodeler Chd1

   Publicly available at the NCBI Gene Expression Omnibus (Accession no: GSE52349).

   http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE52349

1. 1. Core LJ
   2. Waterfall JJ
   3. Gilchrist DA
   4. Fargo DC
   5. Kwak H
   6. Adelman K
   7. Lis JT

   (2012) Defining the status of RNA polymerase at promoters

   Cell Reports 2:1025–1035.

   https://doi.org/10.1016/j.celrep.2012.08.034

   • Google Scholar
2. 1. Fan X
   2. Lamarre-Vincent N
   3. Wang Q
   4. Struhl K

   (2008) Extensive chromatin fragmentation improves enrichment of protein binding sites in chromatin immunoprecipitation experiments

   Nucleic Acids Research 36:e125.

   https://doi.org/10.1093/nar/gkn535

   • Google Scholar
3. 1. Grant CE
   2. Bailey TL
   3. Noble WS

   (2011) FIMO: Scanning for occurrences of a given motif

   Bioinformatics 27:1017–1018.

   https://doi.org/10.1093/bioinformatics/btr064

   • Google Scholar
4. 1. Henikoff JG
   2. Belsky JA
   3. Krassovsky K
   4. MacAlpine DM
   5. Henikoff S

   (2011) Epigenome characterization at single base-pair resolution

   Proceedings of the National Academy of Sciences of USA 108:18318–18323.

   https://doi.org/10.1073/pnas.1110731108

   • Google Scholar
5. 1. Hesselberth JR
   2. Chen X
   3. Zhang Z
   4. Sabo PJ
   5. Sandstrom R
   6. Reynolds AP
   7. Thurman RE
   8. Neph S
   9. Kuehn MS
   10. Noble WS
   11. Fields S
   12. Stamatoyannopoulos JA

   (2009) Global mapping of protein-DNA interactions in vivo by digital genomic footprinting

   Nature Methods 6:283–289.

   https://doi.org/10.1038/nmeth.1313

   • Google Scholar
6. 1. Hnisz D
   2. Abraham BJ
   3. Lee TI
   4. Lau A
   5. Saint-André V
   6. Sigova AA
   7. Hoke HA
   8. Young RA

   (2013) Super-enhancers in the control of cell identity and disease

   Cell 155:934–947.

   https://doi.org/10.1016/j.cell.2013.09.053

   • Google Scholar
7. 1. Kasinathan S
   2. Orsi GA
   3. Zentner GE
   4. Ahmad K
   5. Henikoff S

   (2014) High-resolution mapping of transcription factor binding sites on native chromatin

   Nature Methods 11:203–209.

   https://doi.org/10.1038/nmeth.2766

   • Google Scholar
8. 1. Koerber RT
   2. Rhee HS
   3. Jiang C
   4. Pugh BF

   (2009) Interaction of transcriptional regulators with specific nucleosomes across the Saccharomyces genome

   Molecular Cell 35:889–902.

   https://doi.org/10.1016/j.molcel.2009.09.011

   • Google Scholar
9. 1. Neph S
   2. Vierstra J
   3. Stergachis AB
   4. Reynolds AP
   5. Haugen E
   6. Vernot B
   7. Thurman RE
   8. John S
   9. Sandstrom R
   10. Johnson AK
   11. Maurano MT
   12. Humbert R
   13. Rynes E
   14. Wang H
   15. Vong S
   16. Lee K
   17. Bates D
   18. Diegel M
   19. Roach V
   20. Dunn D
   21. Neri J
   22. Schafer A
   23. Hansen RS
   24. Kutyavin T
   25. Giste E
   26. Weaver M
   27. Canfield T
   28. Sabo P
   29. Zhang M
   30. Balasundaram G
   31. Byron R
   32. MacCoss MJ
   33. Akey JM
   34. Bender MA
   35. Groudine M
   36. Kaul R
   37. Stamatoyannopoulos JA

   (2012) An expansive human regulatory lexicon encoded in transcription factor footprints

   Nature 489:83–90.

   https://doi.org/10.1038/nature11212

   • Google Scholar
10. 1. Orsi GA
    2. Kasinathan S
    3. Zentner GE
    4. Henikoff S
    5. Ahmad K

    (2015) Mapping regulatory factors by immunoprecipitation from native chromatin

    Current Protocols in Molecular Biology 110:21.31.1–21.31.25.

    https://doi.org/10.1002/0471142727.mb2131s110

    • Google Scholar
11. 1. Pott S
    2. Lieb JD

    (2014) What are super-enhancers?

    Nature Genetics 47:8–12.

    https://doi.org/10.1038/ng.3167

    • Google Scholar
12. 1. Rhee HS
    2. Pugh BF

    (2011) Comprehensive genome-wide protein-DNA interactions detected at single-nucleotide resolution

    Cell 147:1408–1419.

    https://doi.org/10.1016/j.cell.2011.11.013

    • Google Scholar
13. 1. Samkurashvili I
    2. Luse DS

    (1996) Translocation and transcriptional arrest during transcript elongation by RNA polymerase II

    The Journal of Biological Chemistry 271:23495–23505.

    https://doi.org/10.1074/jbc.271.38.23495

    • Google Scholar
14. 1. Skene PJ
    2. Hernandez AE
    3. Groudine M
    4. Henikoff S

    (2014) The nucleosomal barrier to promoter escape by RNA polymerase II is overcome by the chromatin remodeler Chd1

    eLife 3:e02042.

    https://doi.org/10.7554/eLife.02042

    • Google Scholar
15. 1. Solomon MJ
    2. Larsen PL
    3. Varshavsky A

    (1988) Mapping protein-DNA interactions in vivo with formaldehyde: evidence that histone H4 is retained on a highly transcribed gene

    Cell 53:937–947.

    https://doi.org/10.1016/S0092-8674(88)90469-2

    • Google Scholar
16. 1. Starick SR
    2. Ibn-Salem J
    3. Jurk M
    4. Hernandez C
    5. Love MI
    6. Chung HR
    7. Vingron M
    8. Thomas-Chollier M
    9. Meijsing SH

    (2015) ChIP-exo signal associated with DNA-binding motifs provides insight into the genomic binding of the glucocorticoid receptor and cooperating transcription factors

    Genome Research 25:825–835.

    https://doi.org/10.1101/gr.185157.114

    • Google Scholar
17. 1. Teytelman L
    2. Ozaydin B
    3. Zill O
    4. Lefrançois P
    5. Snyder M
    6. Rine J
    7. Eisen MB

    (2009) Impact of chromatin structures on DNA processing for genomic analyses

    PLOS ONE 4:e6700.

    https://doi.org/10.1371/journal.pone.0006700

    • Google Scholar
18. 1. Weber CM
    2. Ramachandran S
    3. Henikoff S

    (2014) Nucleosomes are context-specific, H2A.Z-modulated barriers to RNA polymerase

    Molecular Cell 53:819–830.

    https://doi.org/10.1016/j.molcel.2014.02.014

    • Google Scholar

Author details

1. Peter J Skene

   Fred Hutchinson Cancer Research Center, Seattle, United States

   Contribution

   PJS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

2. Steven Henikoff

   1. Fred Hutchinson Cancer Research Center, Seattle, United States
   2. Howard Hughes Medical Institute, Fred Hutchinson Cancer Research Center, Seattle, United States

   Contribution

   SH, Conception and design, Drafting or revising the article

   For correspondence

   steveh@fhcrc.org

   Competing interests

   The authors declare that no competing interests exist.

• 13,277

  views

• 2,578

  downloads

• 108

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.