This paper describes a semiautomated, microscaled ChIP-Seq protocol of DNA regions associated with a specific histone modification (H3K27ac) using a ChIP liquid-handler platform, followed by library preparation using tagmentation. The procedure includes control assessment for quality and quantity purposes and can be adopted to other histone modifications or transcription factors.
Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) is a powerful and widely used approach to profile chromatin DNA associated with specific histone modifications, such as H3K27ac, to help identify cis-regulatory DNA elements. The manual process to complete a ChIP-Seq is labor intensive, technically challenging, and often requires large-cell numbers (>100,000 cells). The method described here helps to overcome those challenges. A complete semiautomated, microscaled H3K27ac ChIP-Seq procedure including cell fixation, chromatin shearing, immunoprecipitation, and sequencing library preparation, for batch of 48 samples for cell number inputs less than 100,000 cells is described in detail. The semiautonomous platform reduces technical variability, improves signal-to-noise ratios, and drastically reduces labor. The system can thereby reduce costs by allowing for reduced reaction volumes, limiting the number of expensive reagents such as enzymes, magnetic beads, antibodies, and hands-on time required. These improvements to the ChIP-Seq method suit perfectly for large-scale epigenetic studies of clinical samples with limited cell numbers in a highly reproducible manner.
The wide use of ChIP-Seq assays for determining fragments of DNA associated with specific histone modifications is in part due to its ability to identify cis-regulatory DNA elements, including active enhancers, promoters, silencers, heterochromatin, and others1,2,3,4. Identification of non-coding regulatory regions across the genome has shown valuable insight to better understand gene regulation in health and diseases4. Previous work from the lab has used ChIP-Seq to show that cis-regulatory elements can play important roles in different cell types5. Transcription factor (TF) ChIP assays has been utilized to show disease associated risk single-nucleotide polymorphisms6.
The use of ChIP-Seq with human clinical samples is challenging, mainly due to the limitation of cell numbers or the desired tissue sample. As a result, there has been a concerted effort in the field to improve and microscale these techniques and as a result, several assays have emerged, such as CUT&TAG5,7,8,9,10,11,12. This assay utilizes a transposase to tagment and isolate genomic regions bound by a specific antibody9. This technique has been able to reduce the cell numbers down to 1,000s and in some cases to a single-cell, however, the use of this technique in translational research and clinical set-up has shown limitations due to the requirements of using live cells for this method9,12. The live cell requirement makes clinical samples logistically difficult to handle and can introduce batch effects if the samples are not processed at the same time. Others have optimized microscaled techniques for formaldehyde-fixed cells, including the development of ChIPmentation11, which is adapted here in a high-throughput manner. The use of fixed cells allows samples to be stored until collection and subsequent processing of all samples together to minimize batch effects.
Here, a semiautomated microscaled ChIP-Seq assay is described which reduces experimental hands-on time to profile histone modifications10. The semiautomated method allows for high-throughput ChIP-Seq assays, allowing for up to 48 samples to be fully processed and ready for sequencing in as little as 5 days, for as few as 10,000 cells per sample using a ChIP liquid-handler. The handler completes the immunoprecipitation (IP) and subsequent washes in an autonomous manner, which helps to reduce variability between samples. The semiautomated method lowers both the hands-on time by over 15 h for 48 samples and the technical variability, enabling large-scale epigenetic studies to be conducted in a reproducible and rapid manner for either primary or cultured cells. The protocol explains the process from start to finish for high quality ChIP-Seq. If the specific machines are not available, the protocol will still be a useful resource to set up and trouble-shoot ChIP-Seq experiments manually.
The assay was performed with three different primary human immune cell types and one cultured cell line (HUT78 – ATCC: TIB-161). For clarity, the protocol has been divided into seven sections: cell fixation, chromatin shearing via sonication, automated chromatin immunoprecipitation, library preparation by DNA fragment tagmentation, library amplification, library purification, followed by DNA quantification. For buffer recipes please refer to Supplementary Table 1.
The Institutional Review Board (IRB) of the La Jolla Institute for Allergy and Immunology (LJI; IRB protocol no. SGE-121-0714) approved the study. Healthy volunteers were recruited and provided leukapheresis samples after written informed consent.
1. Cell fixation
2. Chromatin shearing
NOTE: This protocol is optimized for the chromatin shearing of pellets with 0.3 to 3 x 106 cells in 0.65 mL low binding tubes.
3. Automated ChIP-Seq for histone modification
NOTE: This protocol is designed to run on a ChIP liquid-handler. Although the system can use customized buffers, all buffers are provided with the ChIP kit. The ChIP strips with 8 tubes used in this section are specific to the ChIP liquid-handler.
4. Transposase integration of library adaptors for library preparation
5. Tagmented DNA fragments purification
6. Amplification and size-selection of the purified samples
7. Quantify the final libraries and QC samples using a fluorescence assay
As proof of concept, ChIP-Seq was completed for six human donors with three sets of immune cell types: naive CD4 T cells (CD4), classical monocytes (MO) and natural killer cells (NK), enriched by FACS sorting as described before13. The underlined procedure consists of nine distinct procedures as represented in Figure 1.
Figure 1: General Flowchart for the procedure. (A) A cartoon of the overall procedure (generated in BioRender). (B) Flow-chart for all the major steps for the protocol and the estimated hands-on and total time associated with each day. The sequencing could happen at the end of Day 5 or later with multiple rounds. The timeline can also be staggered throughout the week, where sequential Day 3-4 can be completed multiple times in a week to generate 48 ChIP samples. Please click here to view a larger version of this figure.
After cell isolation by flow cytometry13, sorted cells were centrifuged and cells fixed and stored as described above. Once all the samples were collected, the samples were lysed and prepared for chromatic shearing in batches of 12 as described above. For each sample, the number of cycles to reach optimal sonication was completed10. Quantitative measurement, as well as sheared chromatin fragment size measurements showed great reproducibility of our method on the three sets of immune cells (Figure 2A). The different human immune cells were sonicated in separate batches and yielded very consistently with > 70% of the sample between 100 - 500 bp for 14 cycles (16 s ON, 32 s OFF per cycle). At this point, samples with large fragments after sonication (< 70% of the sample between 100 - 500 bp) were considered as failed. These samples could either sonicated for 1-2 additional cycles or were discarded and replaced later with cells from another pellet. Our method showed none of the samples required more sonication or were eliminated, suggesting absolute robustness of the procedure.
Figure 2: Pre-sequencing QC examples. (A) 1.2% agarose gels show reproducibility of sonication. Sonication samples for 6 donors in three cell types: naive CD4 T cells (CD4), Classical monocytes (MO), and Natural killer cells (NK). The samples were sonicated for 14 cycles (16 s ON, 32 s OFF per cycle). For each sample about 200 ng of decrosslinked chromatin were loaded on a 1% agarose gel. Samples were considered good if more than 70 % of the fragments are within 100-500 bp. (B) Top - Analysis qPCR amplification curves to determine the optimal number of cycles for amplification (Ct where there is ½ the max intensity). The ideal samples have a Ct of about 15 and amplification can be completed up to 2 cycle more of the measured Ct. The arrow is an example of a bad example where the Ct is greater than 18. Bottom - An example of a poor set of samples is shown which have a Ct greater than 18. These samples also showed lower fluorescence intensity. (C) Left - Fragment analyzer electrophoresis traces showed the distribution of final tagmented libraries after amplification and size-selection. Samples with more than 85 % of fragment library lies within 200-1,000 bp were considered as good samples. Measurement of peak intensity of fluorescence is also considered as an important QC parameter, indeed if signal is low, the sample is unlikely to sequence well. Right - Examples for positive samples in CD4, MO, and NK are shown. Please click here to view a larger version of this figure.
After quantification, the samples were run on a ChIP liquid-handler with H3K27ac antibodies, followed by tagmentation with Tn5 transposase enzyme. To determine the appropriate number of amplification cycles by qPCR, 10% of tagmented samples were used. For the determination of the number of cycles for the amplification of the samples, we find the cycle at which the intensity of the sample is half the average maximum for cycle determination (Figure 2B). Samples with Ct values of more than 18 did not perform well post sequencing and their Ct value was thus indicative of a failed ChIP sample. These samples generally also yielded a lower amount of DNA after amplification. Samples (100,000 cells input) with a Ct value equal or lesser than 15 were ideal and samples between 15 and 18 were acceptable but less consistent post sequencing. For samples with less than 100,000 input cells, the Ct values were usually found between 15 and 18 but did not need more than 18 cycles to yield enough product for sequencing.
After DNA-tagmented amplification, libraries were purified and size-selected to obtain an ideal size distribution, ranging from 200 to 1,000 bp, for the NextGen sequencing. Size distribution assessment on each of the libraries was completed because best sequencing data were obtained when more than 85% of the DNA fragments ranged between 200 to 1,000 bp (Figure 2C). Notably, as the same quantity of DNA (measured by fluorescence quantification) was loaded, it was noticed the samples with lower fluorescence intensity generally sequenced poorly (Figure 2C).
Post sequencing, standard quality controls based on the ENCODE ChIP-Seq guidelines were applied5,14,15.
Figure 3: Reproducibility of the immune-cell samples. (A) H3K27ac tracks (UCSC Genome Browser, maximum intensity, smoothing function of 4, all with equally scaled Y-axis) for 6 donors (100,000 cells per replicate) in each cell type (CD4, MO, and NK). Four exemplary loci are shown, two with (IL2RA locus and PTPRC) and two without enrichment for H3K27ac (CCR4 and MS4A1). (B) Pearson correlation between the donors and corresponding correlation plots generated using a 300 bp extension and 500 bp window within the MEDIPS package for each of the cell type replicates16. (C) Merged donor files for each cell type showing H3K27ac tracks (UCSC Genome Browser maximum intensity, smoothing function of 4) in cell type-specific regions (IL17R for CD4, CCR2 for MO, and KLRC1 for NK) and the house-keeping gene B2M, present in all cell types. Please click here to view a larger version of this figure.
For visual quality control, H3K27ac enrichment tracks for display in the UCSC genome browser were prepared. For four gene loci, individual tracks for each sample showed high mapping quality and signal-to-noise ratio reflecting the high consistency and robustness of our assay (Figure 3A). The two loci to the left harbor well-expressed genes in these cell types, while the genes in the two loci to the right are not expressed and served as background controls13 (Figure 3A). Further, the MEDIPS analysis package was used as post-sequencing variable to assess the correlation index between technical replicates (Figure 3B)5,16,17, establishing the degree of correlation for reads enrichment level for 500 bp bins16. For the majority of the pairwise comparisons, Pearson correlations indexes showed more than 90% correlation suggesting high level of consistence between the biological replicates (Figure 3B). Replicates with acceptable correlation were merged to increase signal-to-noise ratio. While cell type-specific loci showed high enrichment in the appropriate cells, a house-keeping gene (B2M) showed very consistent histone modification (Figure 3C). For the analysis, merging tracks from replicates will increase the enrichment, reinforce the specific signal, including for important cell type-specific enhancers, and reduces the inter-individual variability inherent to human samples5.
Although 100,000 cells were used for this study, there was high reproducibility for as few as 10,000 cells in a human cultured T-cell line (HUT78). Correlation analysis between ChIP-Seq dataset performed from samples with less than 100,000 cells showed high reproducibility and correlation down to 10,000 cells (Figure 4A).
Figure 4: Reproducibility of low input samples. (A) Examples of the consistency of H3K27ac ChIP-Seq for cells 100,000 down to 10,000 in HUT-78 cells (a T-cell lymphoma cell line). The tracks (UCSC Genome Browser, maximum intensity, smoothing function of 4, all with equally scaled Y-axis) show the IL4 locus. (B) Pearson correlations of the replicates using a 300 bp extension and 500 bp window within the MEDIPS package16. (C) Pearson correlations between the different cell number groups (100,000, 50,000, and 10,000 cells) using the same MEDIPS parameters as in (B)16. Please click here to view a larger version of this figure.
Pearson correlation analysis showed high correlation index (83% to 92%), suggesting maintenance of signal in low cell number samples. However, there was increased background as the cell numbers were reduced as well as a dropping of the correlation coefficients (Figure 4B). To maintain low background signals, technical duplicates were merged, and the correlation was tested between groups (Figure 4C).
10X Cell Fixation Buffer | |
Compound | Final Concentration |
Formaldehyde solution | 11% |
NaCl | 100 mM |
EDTA, pH 8.0 | 1 mM |
EGTA, pH 8.0 | 0.5 mM |
HEPES, pH 7.5 | 50 mM |
Complete Lysis Buffer | |
Compound | Final Concentration |
Tris-HCI, pH 8.0 | 50 mM |
EDTA, pH 8.0 | 10 mM |
SDS | 0.25% |
Sodium Butyrate | 20 mM |
Protease Inhibitor Cocktail | 1X |
Short-term Lysis Buffer | |
Compound | Final Concentration |
Tris-HCI, pH 8.0 | 50 mM |
EDTA, pH 8.0 | 10 mM |
SDS | 0.25% |
Tagmentation Mix | |
Compound | Final Concentration |
Tris-HCI, pH 8.0 | 10 mM |
MgCl2 | 5 mM |
N,N-dimethylformamide | 10% |
Illumina tagmentation enzyme | 1:24 vol:vol |
CtD Mix | |
Compound | Per sample (µL) |
NextEra Index Primer A (25 µM) | 0.275 |
NextEra Index Primer B (25 µM) | 0.275 |
2X KAPA HiFi HotStart Ready Mix | 2.75 |
1:1000 SYBR Green dye | 0.11 |
ROX passive dye | 0.11 |
Water | Fill to 4 µL |
AMP Mix | |
Compound | Per sample (µL) |
2X KAPA HiFi HotStart Ready Mix | 27.5 |
Water | Fill to 31 µL |
Supplementary Table 1: Buffer recipes.
Supplementary Table 2: Spearman and Pearson sample correlations for the 6 donors and each cell type. Please click here to download this table.
The method described here expands on the ChIPmentation procedure11, which implements a tagmentation library preparation protocol prior to DNA purification, by automating and microscaling the protocol. Since the onset of ChIP-Seq, the required cell numbers have been reduced drastically, from about 20 million cells for histones down to hundreds and even single-cells1,7,10,12,18,19,20,21. These newly developed methods have allowed for a deeper understanding of how cis-regulatory mechanisms are working in cells by increasing the sensitivity and allowing for rare clinical cell populations to be tested5,6,12,17. For instance, one of the more recent and popular procedure, called CUT&TAG, as robust and sensitive ChIP-Seq alternative9. It produces an excellent signal-to-noise ratio as the Tn5 enzyme is covalently bound to protein A and recognizes the Fc chain of the ChIP antibody with high specificity9. Background activity of Tn5-enzyme is reduced as the enzyme is not functional before binding to the target antibody9. However, the implementation of this method in a clinical context is limited since it requires non-fixed, live cells. Also, the removal of DNA fragments from the hypotonic nucleus could have negative effects on the chromatin as it is removed from during the assay. The necessary requirement to work with fresh and living cells is a source of issues for rare clinical samples and for large cohorts of samples, since large cohorts can take numerous years to collect5. Another type of method, drop-ChIP, elegantly uses a microfluidics device to generate droplet based tagmention prior to processing the ChIP19. However, it uses a highly specialized microfluidic device and, while it is possible to complete single-cell ChIP-Seq, it is also limited to the use of live cells7,8,9,18,19. Newer methods relying on ChIP-Seq such as PLAC-Seq or HiChIP, attempt to understanding 3-dimension (3D) interactions between the ChIP-Seq peaks22,23. These 3D methods are exciting as they are identifying cis-regulatory or TF mediated interactions across the genome and better the understanding of the regulation of gene expression in cell types of interest, in healthy tissues and in the context of disease.
There are a few critical steps to consider for the protocol to be successful such as quality of the sonicated chromatin and quality of the antibody. Shearing efficiency is critical, if the chromatin is not sonicated well, the efficiency of the assay decreases drastically24. Sonication is a challenging aspect of ChIP-Seq due to the cell numbers required. On the sonicator used in the protocol, efficiency was drastically reduced under 300,000 cells. This is a challenging aspect in ChIP-Seq as to sonicate under that level would often require enzymatic fragmentation, which is less impartial. As a result, sonication is a major limited factor for true microscaled ChIP-Seq. Other sonication platforms and commercially available kits were tested for sonicating chromatin, but the sonicator used here had the most robust and reproducible results. Another advantage of the sonicator is not having to purchase specialized tubes to run the sonication, which reduces costs when dealing with large number of samples. For optimal sonication, firstly, it is important to pre-warm the sonicator as described above. Second, to lyse the pellet, it is recommended to have the pipette tip touching the bottom of the tube while lysing to break up the cells with more physical constrains. Third, any bubble formation prior to sonication hinders the ability of the sample to be sonicated evenly. If there are any bubbles formed during the lysis, it is important to remove them with a pipette. This can be challenging without removing a lot of sample, but if the tip is lightly pressed against the bubble it can be slowly drawn up without loss of much sample. Lastly, when determining the number of cycles, complete a time-course where every three cycles, sample is removed, purified, and ran on an agarose gel. Avoid over/under sonication of samples as this decreases the ChIP efficiency. If the sample is under sonicated, the large fragments can have a negative effect on the ChIP-Seq quality24. On the other hand, if the sample is over sonicated, there is a risk of the target epitope getting lost in the process.
Another essential part of ChIP-Seq is the quality of the antibody. Prior to running any large-scale study, it is necessary to optimize the antibody which will be used. The goal is to obtain a significantly high signal to noise ratio of known regions of the genome and another is the reproducibility. If the antibody is pulling a lot of background signal, it might be recommended to use a larger input or try a different lot/supplier. This will add time before starting a large-scale experiment, but it is an essential step. To test for the signal-to-noise it is recommended to use qPCR with regions known to be a target of your antibody and another region known to be absent. It has been noticed histone modifications are more robust and easier to optimize than TFs.
The protocol described above provides a robust method for high-throughput histone-modification ChIP-Seq in a semiautonomous, microscaled manner. The method limits the amount of hands-on time and increases the reproducibility over manual ChIP-Seq. Previous studies completed in the lab used manual ChIP on technical replicates and obtained a Spearman correlation average of 0.505, however, with the semiautomated system, the Spearman correlation between different donors with an NK cells average of 0.66 (Supplementary Table 2). This was also completed with about 40% less hands-on time. The method described here has been optimized for histone-modifications (H3K27ac shown here, but the protocol should not need any modification for others) and would only require minor modifications for to be implemented for TF ChIP-Seq. Despite the quality of the antibody, the main modification would be for the sonication time and potentially the buffers used during the IP. Usually, for TF ChIP assays, the method may work better with slightly longer fragments of chromatin (with a range of around 350-800 bp) as TF:DNA complexes are likely less able to be maintained through rigorous sonication6. The buffers might also need to change to a custom mix or other industry available kits, as TFs can behave differently than histone modifications.
Although the automated ChIP liquid-handler was been tested for as few as 10,000 cells, there was a noticeable decrease in reproducibility at lower chromatin concentrations. Due to this, the protocol was not recommended to less than 10,000 cells, with 100,000 cells being the optimal conditions. The protocol was also completed using industry ChIP buffers, which was an added expense but provided higher quality data. The protocol could be modified with regard to the sonication conditions (as long as the sheared chromatin is kept within the same range), buffers could be customized for the immunoprecipitation (IP; optimization may be required), or the ChIP liquid-handler may not be used. A limitation of the protocol is the use of the ChIP liquid-handler, which can be an expensive investment and can only run 16 samples at once. The ChIP liquid-handler is limited to small-scale reactions and cell numbers greater than one million are not recommended. However, the protocol could be completed without it, by completing the IP and wash steps manually. If the IP and washes were completed by hand, the time to complete the assay will increase and the reproducibility may decrease, but this guide will still be useful running a high-quality ChIP-Seq experiment. Of note, other liquid handlers could be adapted to run semiautomated ChIP reactions.
To summarize, the major benefits of this system is the high-throughput nature, since the IP and washing steps are completed autonomously. As such, sequential rounds of ChIP experiments can be completed, allowing up to 48 samples to be fully processed and ready for sequencing in 5 days, with limited hands-on time compared to manual ChIP-Seq experiments. Another benefit is the increased reproducibility since ChIP-Seq can be difficult to obtain highly reproducible results. Other methods either require live cells, complex micro-pipetting systems, or the work to be completed all by hand. This system will have to be optimized for low-input samples (<10,000 cells), ultimately allowing single-cell ChIP reactions. The system is also capable of being adapted for the newer ChIP methods, such as PLAC-Seq and HiChIP22,23.
We thank the Vijayanand lab members for technical help and constructive discussions and Dr. Sharron Squazzo and Mr. Geoffrey Berguet from Diagenode for technical assistance with the sonicator and ChIP liquid-handler machine and protocols. This work was supported by NIH grants AI108564, R01HL114093, S10RR027366 (BD FACSAria II), and S10OD016262 (Illumina HiSeq 2500).
Name | Company | Catalog Number | Comments |
200 µl tube strips (8 tubes/strip) + cap strips | Diagenode | C30020002 | Strip tubes for use on the IP Star; ChIP 8-tube strip |
AMPure XP for PCR Purification | Beckman Coulter | A63880 | SPRI beads |
Axygen 0.6 mL MaxyClear Snaplock Microcentrifuge Tube | Corning | MCT-060-C | 0.65 mL low binding tube |
Bioruptor Pico Sonicator | Diagenode | B01060010 | Sonicator used in the lab but others can be used |
ChIP DNA Clean & Concentrator (Capped Columns) | Zymo Research | D5205 | DNA clean-up kit |
Dynabeads Protein A for Immunoprecipitation | ThermoFisher | 10001D | |
EDTA (0.5 M), pH 8.0, RNase-free | ThermoFisher | AM9260G | |
EGTA pH 8.0 | Millipore Sigma | E3889-25G | |
Eppendorf ThermoMixer C | Eppendorf | 2231000667 | |
Formaldehyde solution | Millipore Sigma | 252549-1L | |
Glycine | Millipore Sigma | 50046-250G | |
H3K27ac polyclonal antibody - Premium | Diagenode | C15410196 | |
HEPES (1 M) pH 7.5 | ThermoFisher | 15630080 | |
IDT for Illumina Nextera DNA Unique Dual Indexes | Illumina | 20027213 | |
Illumina Tagment DNA Enzyme and Buffer Small Kit | Illumina | 20034197 | |
IP-Star Compact Automated System | Diagenode | B03000002 | Automated system for ChIP-Seq studies; ChIP liquid handler |
KAPA HiFi HotStart ReadyMix | Roche | KK2601 | PCR mix |
Medium reagent container for SX-8G IP-Star Compact | Diagenode | C30020003 | |
MgCl2 (magnesium chloride) (25 mM) | ThermoFisher | R0971 | |
N,N-Dimethylformamide | Millipore Sigma | D4551-250ML | CAUTION - low flash point |
NaCl (5 M), RNase-free | ThermoFisher | AM9760G | |
PBS (10X), pH 7.4 | ThermoFisher | 70011044 | |
PCR Flex-free 8-tube stripes, attached individual optical caps | USA Scientific | 1402-4700 | 8 strip tubes, 0.2 mL 8-tube strip |
Proteinase Inhibitor Cocktail | Millipore Sigma | P8340 | |
Proteinase K Solution (20 mg/mL), RNA grade | ThermoFisher | 25530049 | |
PureLink RNase A (20 mg/mL) | ThermoFisher | 12091021 | |
Quant-iT PicoGreen dsDNA Reagent | ThermoFisher | P7581 | Used in the flourescence quantification |
QuantStudio 6 Flex Real-Time PCR System | ThermoFisher | 4485699 | qPCR |
ROX Reference Dye | ThermoFisher | 12223012 | |
Sodium butyrate | Millipore Sigma | 303410-100G | |
SYBR Gold Nucleic Acid Gel Stain (10,000X Concentrate in DMSO) | ThermoFisher | S11494 | nucleic acid dye |
SYBR Green I Nucleic Acid Gel Stain - 10,000X concentrate in DMSO | ThermoFisher | S7563 | |
TE Buffer | ThermoFisher | 12090015 | |
Tips (bulk) | Diagenode | C30040020 | Tips for the IP Star |
True MicroChIP Kit | Diagenode | C01010130 | Contains all the buffers for the IP; ChIP kit |
UltraPure 1M Tris-HCI, pH 8.0 | ThermoFisher | 15568025 | |
UltraPure SDS Solution, 10% | ThermoFisher | 24730020 |
This article has been published
Video Coming Soon
ABOUT JoVE
Copyright © 2024 MyJoVE Corporation. All rights reserved