- Methodology article
- Open Access
High-throughput ChIPmentation: freely scalable, single day ChIPseq data generation from very low cell-numbers
BMC Genomics volume 20, Article number: 59 (2019)
Chromatin immunoprecipitation coupled to sequencing (ChIP-seq) is widely used to map histone modifications and transcription factor binding on a genome-wide level.
We present high-throughput ChIPmentation (HT-ChIPmentation) that eliminates the need for DNA purification prior to library amplification and reduces reverse-crosslinking time from hours to minutes.
The resulting workflow is easily established, extremely rapid, and compatible with requirements for very low numbers of FACS sorted cells, high-throughput applications and single day data generation.
The combination of chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq) has become the method of choice for mapping chromatin-associated proteins and histone-modifications on a genome-wide level.
The ChIP-seq methodology has rapidly developed [1,2,3,4]. Despite this, performing ChIP-seq on limited cell-numbers and in a high-throughput manner remains technically challenging. This is largely due to decreasing input material leading to progressively increasing losses of material during DNA preparation and inefficiencies of enzymatic reactions used for library preparation. While elegant strategies have been developed to resolve these issues, they remain laborious and have not seen wider use [5,6,7,8,9,10,11,12].
ChIPmentation  effectively alleviates the issues associated with traditional library preparation methodologies by introducing sequencing-compatible adapters to bead-bound chromatin using Tn5 transposase (tagmentation). While fast and convenient, the methodology still relies on the use of traditional reverse crosslinking and DNA purification procedures prior to library amplification, hampering processing time, DNA recovery, and limiting scalability for high-throughput applications.
Here, we present freely scalable high-throughput ChIPmentation (HT-ChIPmentation) that by eliminating the need for DNA purification and traditional reverse-crosslinking prior to library amplification, dramatically reduces required time and input cell numbers. In comparison with current ChIP-seq variants [3, 5,6,7,8,9,10,11,12], HT-ChIPmentation is technically simple, extremely rapid and widely applicable, being compatible with both very low cell number requirements and high-throughput applications.
The adapters introduced by Tn5 are covalently linked only to one strand of the tagmented DNA. The complete adapters, compatible with PCR amplification, are created through a subsequent extension reaction. With this in mind, we reasoned that performing adapter extension of tagmented bead-bound chromatin and high-temperature reverse crosslinking , would allow us to bypass the DNA purification step.
To validate this approach and benchmark it against standard ChIPmentation (Fig. 1a and Additional file 1: Figure S1), we FACS sorted defined numbers of formaldehyde fixed cells and performed ChIP with subsequent library preparation on cell numbers ranging from 0.1 to 150 k cells. HT-ChIPmentation indeed produced excellent sequencing profiles (Fig. 1b), and a consistent library size over > 100-fold difference in input cell numbers (Additional file 1: Figure S2A).
Looking specifically at H3K27Ac (a histone modification demarcating active promoters and enhancers ) HT-ChIPmentation and ChIPmentation samples generated in parallel from high cell-numbers (50–150 k cells), both methods generated high-quality data that is comparable in regard to: concordance of library profiles (Fig. 1b); mappability of sequencing reads (Additional file 1: Table S1); correlation between samples (Fig. 1c); number, quality scores and signal range of identified peaks (Additional file 1: Figure S2B–D); and peak overlap (Fig. 1d).
To perform accurate peak calling, input controls were generated by direct tagmentation of 500 cell equivalents of sonicated chromatin (5% of 10 k sonicated cells), subsequently processed in parallel with corresponding 10 k HT-ChIPmentation samples (Fig. 1a). The HT-ChIPmentation compatible input controls produced similar results as input controls prepared using traditional library preparation methodology, in terms of library profiles and even genomic coverage (Fig. 1b and e).
We next compared H3K27Ac HT-ChIPmentation and ChIPmentation samples from progressively lower input cell-numbers. As expected, eliminating losses associated with DNA purification allowed HT-ChIPmentation samples to maintain much higher library complexity (> 75% unique reads down to 2.5 k cells) than ChIPmentation samples generated from the same number of cells (Fig. 1f). This difference in library quality was directly reflected in HT-ChIPmentation samples generated from a few thousand cells maintaining: consistent high quality library profiles (Fig. 1b); mappability (Additional file 1: Table S1); number, quality scores and signal range of identified peaks (Additional file 1: Figure S2B–D); high correlation between samples (Fig. 1g); and high peak overlap (Fig. 1h). Similar results were obtained for H3K27Ac HT-ChIPmentation data generated in a single day (Fig. 1b, g and Additional file 1: Figure S2B–D). Based on the same metrics, CTCF (a chromatin organizing protein ) HT-ChIPmentation experiments further verified the robustness of the method with cell numbers in the range of a few thousands cells (Fig. 1b, g, h; Additional file 1: Figures S2B–D and S3A–B).
Here we present HT-ChIPmentation, an improved and simplified tagmentation based approach to produce ChIP-seq libraries. We demonstrate that the adapters introduced by Tn5 can be extended directly on the bead-bound chromatin. Through this, we can combine ChIPmentation  with high-temperature reverse crosslinking and direct library amplification without prior DNA purification . Even compared to the already technically simple and fast ChIPmentation method, HT-ChIPmentation is easier to perform and greatly reduces the time needed to produce sequencing ready libraries (Fig. 1i). In fact, HT-ChIPmentation together with sequencing can be performed in a single day (Fig. 1b, g and i; Additional file 1: Figure S2B–D). This makes the protocol ideal for rapid data generation and compatible with the development of clinical diagnostic/prognostic applications relying on chromatin associated features to distinguish, for example, tumor subtypes [15, 16].
The removal of the DNA purification step, allows for fully taking advantage of that tagmentation of chromatin – as opposed to traditional adapter ligation [6, 8] – remains highly effective even with very limited input material ( and Additional file 1: Figure S2A). Together, the reduced losses of material and effective addition of adapters, allows HT-ChIPmentation to be performed on just a few thousand FACS sorted cells with maintained quality and library complexity. Hence, HT-ChIPmentation provides a robust and technically simple workflow for characterizing epigenetic changes and transcription factor binding in rare subsets of cells.
Input controls are commonly used to exclude biases in the input material and as a negative control for identification of peak regions. Here we show that input controls can be prepared in parallel with HT-ChIPmentation samples, through direct tagmentation and library amplification of sonicated chromatin. The protocol requires very limited material (500 cell equivalents of sonicated chromatin), making it both feasible and convenient to directly prepare adequate controls for peak finding, also from rare subsets of cells.
The simplicity of the HT-ChIPmentation protocol – allowing for performing all steps from cells to amplified sequencing ready library without DNA purification – makes it perfectly suited for epigenetic characterization at any scale. While HT-ChIPmentation is directly compatible with full automation, experiments presented here were simply performed in 96-well plates using a multi-channel pipette, demonstrating that HT-ChIPmentation makes it highly feasible to perform epigenome scale projects in a matter of days using standard laboratory equipment.
Here we introduce HT-ChIPmentation, an improved tagmentation based ChIP-seq protocol that through the extension of the Tn5-inserted adapters on bead-bound chromatin, allows for direct library amplification without prior DNA purification. In comparison to current state-of-the-art ChIP-seq protocols [3, 5,6,7,8,9,10,11,12], HT-ChIPmentation is technically simple, extremely rapid and widely applicable, being compatible with very low cell number requirements, high-throughput applications and single day data generation. Taken together, HT-ChIPmentation provides a versatile and simplistic workflow attractive as the mainstay protocol for epigenome projects of any scale.
Cultured MEC1 cells were stained with LIVE/DEAD fixable Aqua stain (Invitrogen) to allow for excluding cells dead already prior to fixation (during subsequent FACS sorting) and fixed using 1% PFA (Pierce). Aliquots of 10 k cells were FACS sorted directly into 100 μl SDS lysis buffer (50 mM Tris/HCl, 0.5% SDS, and 10 mM EDTA) supplemented with 1X cOmplete EDTA-free protease inhibitor (Roche) and stored at − 80 °C until use. For aliquots of cells (50 and 150 k), where the sheath fluid volume is non-negligible, cells were sorted into PBS, spun down (2000 g 5 min) and resuspended in 100 μl SDS lysis buffer prior to freezing. Sorting was performed using a BD FACSAriaIIu cell sorter (BD Biosciences) with an 85 μm nozzle.
Chromatin immunoprecipitation and tagmentation
For ChIP, polyclonal anti-H3K27Ac (Diagenode, cat# C15410196, lot# A1723-0041D) antibody or anti-CTCF (Diagenode, cat# C15410210, lot# A2359-00234P) antibody was added to Protein G-coupled Dynabeads (ThermoFisher) in PBS with 0.5% BSA and incubated with rotation for 4 h at 4 °C (0.5 h at RT for HT-ChIPmentation samples processed in a single day). For 50–150 k cells, 10 μl beads incubated with 3 μg H3K27Ac or 1.5 μg CTCF antibody were used per ChIP. For 0.1–10 k cells, 2 μl beads incubated with 0.6 μg H3K27Ac or 0.3 μg CTCF antibody were used per ChIP. Fixed cells (FACS sorted) frozen in SDS lysis buffer were thawed at room temperature. To perform ChIP on < 10 k cells, aliquots were diluted with SDS lysis buffer and 100 μl containing the appropriate number of cells were processed. Cells were sonicated for 12 cycles of 30 s on/30 s off on high power using a Bioruptor Plus (Diagenode). To neutralize the SDS, Triton X100 was added to a final concentration of 1% along with 2 μl 50x cOmplete protease inhibitor (final 1x). Samples were incubated at room temperature for 10 min and when applicable 5% aliquots were saved for preparation of input controls. Antibody-coated Dynabeads were washed with PBS with 0.5% FCS and mixed with cell lysate in PCR tubes. Tubes were incubated rotating overnight (or 4 h for HT-ChIPmentation samples processed in a single day) at 4 °C.
Immunoprecipitated chromatin was washed with 150 μl of low-salt buffer (50 mM Tris/HCl, 150 mM NaCl, 0.1% SDS, 0.1% NaDOC, 1% Triton X-100, and 1 mM EDTA), high-salt buffer (50 mM Tris/HCl, 500 mM NaCl, 0.1% SDS, 0.1% NaDoc, 1% Triton X-100, and 1 mM EDTA) and LiCl buffer (10 mM Tris/HCl, 250 mM LiCl, 0.5% IGEPAL CA-630, 0.5% NaDOC, and 1 mM EDTA), followed by two washes with TE buffer (10 mM Tris/HCl and 1 mM EDTA) and two washes with ice cold Tris/HCl pH 8. For tagmentation, bead bound chromatin was resuspended in 30 μl of tagmentation buffer, 1 μl of transposase (Nextera, Illumina) was added and samples were incubated at 37 °C for 10 min followed by two washes with low-salt buffer.
High-throughput ChIPmentation library preparation
For High-throughput ChIPmentation (HT-CM) samples, bead bound tagmented chromatin was diluted in 20 μl of water. PCR master mix (Nextera, Illumina) and indexed amplification primers  (0.125uM final concentration) was added and libraries prepared using the following program: 72 °C 5 min (adapter extension); 95 °C 5 min (reverse cross-linking); followed by 11 cycles of 98 °C 10s, 63 °C 30s and 72 °C 3 min.
For preparation of HT-CM compatible input controls, 1 μl of 50 mM MgCl2 was added to 5 μl sonicated lysate (5% aliquot of 10 k samples) to neutralize the EDTA in the SDS lysis buffer. Thirty microliters of tagmentation buffer and 1 μl transposase (Nextera, Illumina) was added, and samples were incubated at 37 °C for 10 min. 22.5 μl of the transposition reaction were combined with 15 μl of PCR master mix and 2.5 μl of primer mix (Nextera, Illumina). Libraries were subsequently amplified as described for HT-ChIPmentation samples.
ChIPmentation library preparation
For standard reverse crosslinking, chromatin complexes were diluted with 200 μl ChIP elution buffer (10 mM Tris/HCl, 0.5% SDS, 300 mM NaCl, and 5 mM EDTA) and 2 μl of 20 μg/ml proteinase K (Thermo Scientific). Samples were vortexed and incubated with shaking overnight at 65 °C. After reverse crosslinking, 1 μl 20 μg/ml RNase (Sigma) was added and incubated at 37 °C for 30 min. After another 2 h of incubation with 2 μl of proteinase K (20 mg/ml) at 55 °C, samples were placed in a magnet to trap magnetic beads and supernatants were collected. DNA purification was carried out using Qiagen MinElute PCR Purification Kit. Fifteen microliters of PCR master mix and 5 μl of primer mix (Nextera, Illumina) was added to 20 μl of eluted DNA, and libraries were amplified as described for HT-ChIPmentation libraries.
Preparation of conventional input control
Sonicated material from 50 k cells was reverse crosslinked as described for ChIPmentation. Two nanograms of DNA was used for library preparation using the ThruPLEX DNA-seq kit (Rubicon Genomics) with 11 cycles of PCR amplification.
Post-PCR library cleanup and sequencing
After PCR amplification, library cleanup was done using Agencourt AmPureXP beads (Beckman Coulter) at a ratio of 1:1. DNA concentrations in purified samples were measured using the Qubit dsDNA HS Kit (Invitrogen). Libraries were pooled and single-end sequenced (50 cycles) using the Nextseq500 platform (Illumina).
Basic processing of ChIP-seq and input control sequencing data
Quality of the sequenced samples was assessed using FastQC v0.11.5 . Samples were mapped to the human reference genome (hg19) using Bowtie2 v2.2.3  with default settings. Further basic processing was performed using HOMER v4.8.3 . Specifically, mapped reads were converted into tagdirectories by the makeTagDirectory command using settings for the human genome (-genome hg19) and removing duplicate reads by allowing only one tag to start per base pair (-tbp 1).
Genome browser visualizations
Bedgraphs were created for each sample using HOMER’s makeUCSCfile. Tracks were uploaded and visualized using the UCSC genome browser .
Peak finding and plotting peak metrics
Peak finding was performed using the findPeaks command in HOMER. Peaks were called using default settings for histone modifications (-style histone) and transcription factors (-style factor) for H3K27Ac and CTCF respectively with input (-i) as a control. Visualization was done in R v3.1.0 , using the built in barplot and boxplot R-functions to plot peak numbers and peak quality scores, respectively.
Making and annotating peak catalogs
Peak catalogs were created by merging all peak files of samples analyzed using HOMER’s mergePeaks command. Setting used (-size given) ensured that peaks with literal overlap were merged to one peak while peaks unique to one sample were directly added to the peak catalog. Subsequently, peak catalogs were annotated with unnormalized (-raw) read counts within peaks in the catalog for each individual sample using HOMER’s annotatePeaks.pl script.
Plotting peak read distributions and correlation between samples
Raw counts were log normalized in R as follows: log(df[,countsCols]+1,2). Log2 counts were subsequently plotted using the build in boxplot R-function. These same Log2 counts were used to calculate sample correlations, using the build-in cor R-function with spearman correlation. Correlation matrices were visualized with the pheatmap function from the pheatmap R-package using color scales generated with the build-in colorRampPalette R-function.
Plotting reads within 1 kb bins for input control samples
A file containing 1 kb bins covering the whole genome was created using the makewindows command from bedtools v2.26.0  using a window size of 1 kb (-w 1000). Chromosome sizes were retrieved as follows: mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -e "select chrom, size from hg19.chromInfo" > hg19.genome. Raw reads in each 1 kb bin for each input control were counted using HOMER’s annotatePeaks.pl script, as described above. Raw read distributions were converted to RPKM in R based on the standard RPKM formula. Resulting RPKM distributions were plotted with the build-in boxplot R-function.
Determining top peak overlap
Peaks identified in individual samples were overlapped with in-house code using the IRanges  R-package. Top peaks overlap was considered to be the percentage of high quality peaks (50% of peaks with highest quality scores) in the reference sample that overlap (≥1 bp) with a peak in the second sample. For purposes of determining peak overlap, CTCF peaks were extended with 50 bp up and downstream, considering findPeaks with -style factor only calls a small region around the peak maximum. Peak overlaps were visualized using the pheatmap function from the pheatmap R-package using color scales generated with the build-in colorRampPalette R-function.
Comparing library complexity
To compare duplication rates between HT-ChIPmentation and ChIPmentation samples, fastq files were randomly down-sampled to the total number of reads in the smallest file for each cell number. Down sampling was performed using the fastq-sample script from fastq-tools v0.8 . Fraction of unique reads was subsequently determined for each file using FastQC v0.11.5.
Motif enrichment analysis
Enrichments of known transcription factor binding motifs in peaks were identified using HOMER’s findMotifsGenome.pl script with default settings.
Chromatin immunoprecipitation with high-throughput sequencing
Furey TS. ChIP–seq and beyond: new and improved methodologies to detect and characterize protein–DNA interactions. Nat Rev Genet. 2012;13:1–13.
Sundaram AYM, Hughes T, Biondi S, Bolduc N, Bowman SK, Camilli A, Chew YC, Couture C, Farmer A, Jerome JP, et al. A comparative study of ChIP-seq sequencing library preparation methods. BMC Genomics. 2016;17:1–12.
Schmidl C, Rendeiro AF, Sheffield NC, Bock C. ChIPmentation: fast, robust, low-input ChIP-seq for histones and transcription factors. Nat Methods. 2015;12:963–5.
Rotem A, Ram O, Shoresh N, Sperling RA, Goren A, Weitz DA, Bernstein BE. Single-cell ChIP-seq reveals cell subpopulations defined by chromatin state. Nat Biotechnol. 2015;33:1165–72.
Lara-Astiaso D, Weiner A, Lorenzo-Vivas E, Zaretsky I, Jaitin DA, David E, Keren-Shaul H, Mildner A, Winter D, Jung S, et al. Immunogenetics. Chromatin state dynamics during blood formation. Science. 2014;345:943–9.
Wallerman O, Nord H, Bysani M, Borghini L, Wadelius C. lobChIP: from cells to sequencing ready ChIP libraries in a single day. Epigenetics Chromatin. 2015;8:25.
Adli M, Zhu J, Bernstein BE. Genome-wide chromatin maps derived from limited numbers of hematopoietic progenitors. Nat Publ Group. 2010;7:615–8.
Jakobsen JS, Bagger FO, Hasemann MS, Schuster MB, Frank A-K, Waage J, Vitting-Seerup K, Porse BT. Amplification of pico-scale DNA mediated by bacterial carrier DNA for small-cell-number transcription factor ChIP-seq. BMC Genomics. 2015;16:46.
Blecher-Gonen R, Barnett-Itzhaki Z, Jaitin D, Amann-Zalcenstein D, Lara-Astiaso D, Amit I. High-throughput chromatin immunoprecipitation for genome-wide mapping of in vivo protein-DNA interactions and epigenomic states. Nat Protoc. 2013;8:539–54.
Sachs M, Onodera C, Blaschke K, Ebata KT, Song JS, Ramalho-Santos M. Bivalent chromatin Marks developmental regulatory genes in the mouse embryonic germline in vivo. Cell Rep. 2013;3:1777–84.
Shankaranarayanan P, Mendoza-Parra M-A, Walia M, Wang L, Li N, Trindade LM, Gronemeyer H. Single-tube linear DNA amplification (LinDA) for robust ChIP-seq. Nat Methods. 2011;8:565–7.
Zwart W, Koornstra R, Wesseling J, Rutgers E, Linn S, Carroll JS. A carrier-assisted ChIP-seq method for estrogen receptor-chromatin interactions from breast cancer core needle biopsy samples. BMC Genomics. 2013;14:232.
Creyghton MP, Cheng AW, Welstead GG, Kooistra T, Carey BW, Steine EJ, Hanna J, Lodato MA, Frampton GM, Sharp PA, et al. Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc Natl Acad Sci. 2010;107:21931–6.
Phillips JE, Corces VG. CTCF: master weaver of the genome. Cell. 2009;137:1194–211.
Muller-Tidow C, Klein HU, Hascher A, Isken F, Tickenbrock L, Thoennissen N, Agrawal-Singh S, Tschanter P, Disselhoff C, Wang Y, et al. Profiling of histone H3 lysine 9 trimethylation levels predicts transcription factor activity and survival in acute myeloid leukemia. Blood. 2010;116:3564–71.
Lin CY, Erkek S, Tong Y, Yin L, Federation AJ, Zapatka M, Haldipur P, Kawauchi D, Risch T, Warnatz HJ, et al. Active medulloblastoma enhancers reveal subgroup-specific cellular origins. Nature. 2016;530:57–62.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10:1213–8.
FastQC. www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 1 Feb 2018.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Publ Group. 2012;9:357–9.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
James Kent W, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D. The human genome browser at UCSC. J Genome Res. 2002;12:996–1006.
The R Project for Statistical Computing. www.r-project.org/. Accessed 1 Feb 2018.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
Lawrence M, Huber W, Pages H, Aboyoun P, Carlson M, Gentleman R, Morgan M, Carey V. Software for computing and annotating genomic ranges. PLoS Comput Biol. 9(8):e1003118.
fastq-tools. https://homes.cs.washington.edu/~dcjones/fastq-tools/. Accessed 1 Feb 2018.
European Nucleotide Archive (ENA). www.ebi.ac.uk/ena. Accessed 1 Feb 2018.
We thank: Prof. Joakim Dillner and colleagues for access to the NextSeq 500 system; Prof. Anders Rosén for providing the MEC1 cell-line; Dr. Yin C. Lin for critical reading of the manuscript; and UPPMAX Next Generation Sequencing Cluster & Storage (UPPNEX) for computational resources.
This work was funded by the Swedish Cancer Foundation (Cancerfonden), the Swedish Research Council (VR), the Knut and Alice Wallenberg Foundation (KAW) and the Swedish Foundation for Strategic Research (SSF). The funding bodies had no role in the design of the study, collection, analysis and interpretation of data, or in writing the manuscript.
Availability of data and materials
The generated data sets are available from the European Nucleotide Archive  under the study accession number: PRJEB23059.
Ethics approval and consent to participate
Consent for publication
CS has a pending patent application for ChIPmentation. CG, ADP, and RM have no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Schematic comparison of the ChIPmentation and high-throughput ChIPmentation protocols. Figure S2. High-throughput ChIPmentation (HT-CM) samples maintain library quality over progressively lower input cell numbers. Figure S3. High-throughput ChIPmentation (HT-CM) maintains high library complexity in CTCF samples. Table S1. Sample details. (PDF 323 kb)
About this article
Cite this article
Gustafsson, C., De Paepe, A., Schmidl, C. et al. High-throughput ChIPmentation: freely scalable, single day ChIPseq data generation from very low cell-numbers. BMC Genomics 20, 59 (2019). https://doi.org/10.1186/s12864-018-5299-0
- Chromatin immunoprecipitation
- High-throughput genomics