- Open Access
ChIP-less analysis of chromatin states
Epigenetics & Chromatinvolume 7, Article number: 7 (2014)
Histone post-translational modifications (PTMs) are key epigenetic regulators in chromatin-based processes. Increasing evidence suggests that vast combinations of PTMs exist within chromatin histones. These complex patterns, rather than individual PTMs, are thought to define functional chromatin states. However, the ability to interrogate combinatorial histone PTM patterns at the nucleosome level has been limited by the lack of direct molecular tools.
Here we demonstrate an efficient, quantitative, antibody-free, chromatin immunoprecipitation-less (ChIP-less) method for interrogating diverse epigenetic states. At the heart of the workflow are recombinant chromatin reader domains, which target distinct chromatin states with combinatorial PTM patterns. Utilizing a newly designed combinatorial histone peptide microarray, we showed that three reader domains (ATRX-ADD, ING2-PHD and AIRE-PHD) displayed greater specificity towards combinatorial PTM patterns than corresponding commercial histone antibodies. Such specific recognitions were employed to develop a chromatin reader-based affinity enrichment platform (matrix-assisted reader chromatin capture, or MARCC). We successfully applied the reader-based platform to capture unique chromatin states, which were quantitatively profiled by mass spectrometry to reveal interconnections between nucleosomal histone PTMs. Specifically, a highly enriched signature that harbored H3K4me0, H3K9me2/3, H3K79me0 and H4K20me2/3 within the same nucleosome was identified from chromatin enriched by ATRX-ADD. This newly reported PTM combination was enriched in heterochromatin, as revealed by the associated DNA.
Our results suggest the broad utility of recombinant reader domains as an enrichment tool specific to combinatorial PTM patterns, which are difficult to probe directly by antibody-based approaches. The reader affinity platform is compatible with several downstream analyses to investigate the physical coexistence of nucleosomal PTM states associated with specific genomic loci. Collectively, the reader-based workflow will greatly facilitate our understanding of how distinct chromatin states and reader domains function in gene regulatory mechanisms.
Post-translational modifications (PTMs) on histones constitute an intricate ‘language’ that instructs chromatin-based processes . Different histone PTMs deposited across multiple residues on histone proteins generate a combinatorial pattern. The number of combinations can be further amplified by the presence of two copies of each histone (H3, H4, H2A and H2B) within a nucleosome. One intriguing hypothesis is that PTM combinations better categorize unique chromatin states than an assessment of an individual PTM. Evidence of the importance of combinatorial PTMs comes from studies on embryonic stem cells, where the coexistence of ‘active mark’ H3K4me3 and ‘repressive mark’ H3K27me3 at differentiation-related gene promoters poises genes for subsequent regulation . Also, the coexistence of H3S10ph and H3K9me3 is a signal to evict heterochromatin protein 1 during the M phase . How combinatorial PTM patterns reflect specific chromatin states or biological processes remains incompletely understood. Decoding the enormous combinatorial capacity of PTM patterns [4, 5] has been hindered by a lack of proper tools to distinguish different combinatorial PTM patterns.
Current methodologies to define chromatin states rely almost exclusively on histone antibodies, which have been developed to recognize individual PTMs. As summarized in Table 1, each of these approaches has its own limitations for analysis of the combinatorial histone PTM patterns. Proper interpretation of such analyses relies on a thorough understanding of the specificity of these antibodies. However, recent studies reveal cross-reactivity to unexpected antigens (‘off-targets’) of many commercial histone antibodies [6–11]. Lot-to-lot variability, high costs and engineering difficulties are critical concerns for the development and commercialization of histone antibodies [12, 13]. All these have challenged the use of histone antibodies in accurately probing chromatin states. Therefore an antibody-free approach to directly enrich combinatorial PTM patterns that maintains native nucleosome structure and combinatorial PTM complexity for downstream quantitation and genome mapping of physical intranucleosomal PTM connections will greatly simplify attempts to characterize chromatin states.
To address this bottleneck, we developed an epigenetic platform utilizing ‘reader’ domains to enrich, identify and quantify chromatin states. Reader domains are a diverse collection of histone-interacting protein modules that can interpret histone PTM language in cells . Previously, we and other groups described the ability of reader domains to discriminate combinatorial PTM patterns [24, 25], suggesting their potential for probing native chromatin states. Using three readers as proof of concept, we demonstrated a powerful workflow involving a combinatorial PTM histone peptide microarray and a reader-based affinity enrichment platform (matrix-assisted reader chromatin capture, or MARCC) to capture and quantify chromatin states maintained on native nucleosomes.
Results and discussion
The chromatin affinity enrichment requires highly specific ‘baits’ that can distinguish different PTM states, making readers such as malignant brain tumor domains unsuitable for such application due to their insensitivity to residues surrounding the targeted methyl lysine group . Three reader domains were selected (Additional file 1: Figure S1) based on their reported chromatin-binding specificities: plant homeodomain (PHD) zinc fingers from ING2 and AIRE, and the ADD (ATRX-DNMT3A-DNMT3L) domain from ATRX [27–32]. To evaluate the specificity of reader domains systematically, we designed a histone peptide library that featured a more comprehensive collection of histone peptides than previous platforms [6, 7] (compared in Table 2). The new library contained 746 tiled peptide species, ≈60% of which comprised new PTMs or novel combinations (Figure 1A and Additional file 2). Such extensive PTM combinations were critical in evaluating the effects of combinatorial PTM patterns on reader binding specificity.
Reader domains were recombinantly expressed and purified as HaloTag fusions, which permitted sensitive detection by covalent labeling with fluorescence  when screened on microarrays containing ≈ 4,500 peptides (Figure 1B). An additional fluorescence dye was introduced during peptide printing to enable dual-channel fluorescence detection (532 nm/635 nm) and allowing easy identification of misprinting (Figure 1C). Binding events were quantified from fluorescence signal intensities at 635 nm of spatially addressed peptides. Each peptide was printed six times on each microarray for statistical analysis and proper control peptides were included in each subarray (Figure 1C). For comparison, two ChIP-grade commercial antibodies for H3K9me3 and H3K4me3 were analyzed using the microarrays.
Microarray analysis revealed high selectivity for amino-acid sequence and combinatorial PTM states with reader domains (for complete dataset, see Additional file 3). Sequence specificity was assessed by comparing binding signals for a panel of methylated or unmodified peptides within the library (Figure 2A, Additional file 4: Figure S2 and Additional file 5: Figure S3). ATRX-ADD, AIRE-PHD and ING2-PHD bound three different PTM states, H3K9me3, H3K4unmod and H3K4me3, respectively, with no evident off-targets. By contrast, the two antibodies exhibited unsatisfactory sequence specificities. Most notably, there was dramatic off-target binding to H3K79me3 peptide using the H3K9me3 antibody (ab8898) and unexpected binding to H3K18me3 using the H3K4me3 antibody (ab8580) (Figure 2A and Additional file 5: Figure S3). Such cross-reactivity of histone antibodies at multiple sites is particularly concerning because the biological functions of histone PTMs are site-dependent and therefore confound the interpretation of antibody-dependent data. Collectively, our binding results further underscore specificity issues with histone antibodies and highlight the high specificity of the tested reader domains.
Within the same amino-acid sequence, the combinatorial PTM pattern is an important component of the chromatin landscape, and an appropriate enrichment strategy should select among diverse PTM states. Binding signals resulting from the same sequence were extracted from the microarray and the effects on coexisting (combinatorial) modifications in the H3 tail were assessed (Figure 2B). Modifications at R2 and T3 diminished binding of AIRE-PHD and ING2-PHD to the H3 tail (Figure 2B). For ATRX-ADD, K9me3 binding was abrogated by modifications at T3, T6 and R8 (Figure 2B). Most importantly, ATRX-ADD but not the antibody ab8898 discriminated against the bivalent peptide H3K4me3K9me3 (Figure 2B), suggesting that ATRX-ADD can uniquely enrich K4me0K9me3-containing chromatin. The strict specificity for K9me3 within the proper surrounding PTM state and a correct N-terminal context (Figure 2B and Additional file 5: Figure S3) can be explained by the extensive interactions of ATRX-ADD with H3 1-10aa [27, 28]. By comparison, the antibodies were less affected by surrounding residues and PTMs, failing to distinguish different PTM combinations (Figure 2B). We reasoned that the recognition specificity is largely dependent on the contact surface between the binding module (readers or antibodies) and its histone target. Reader modules with more surface interaction with the histone molecules, as we have shown here with ATRX-ADD, ING2-PHD, AIRE-PHD, will have longer sequence discrimination and therefore greater sensitivity to combinatorial PTMs within those contacts. For example, the linked domains of ATRX-ADD revealed increased binding specificity that results from the larger interaction surface with histone substrates (Figure 2).
Previous studies have demonstrated that reader domains interact with preferred histone substrates at low micromolar affinity [27–32], a range equal to or weaker than antibody recognition, which varies from micromolar to nanomolar . To test whether such micromolar affinity is adequate for enriching specific nucleosomal PTMs, we next assessed the ability of ATRX-ADD to discriminate PTMs on a folded nucleosome structure, consisting of reconstituted mononucleosomes bearing either H3Kc4me3 or H3Kc9me3 (chemical analogs of H3K4me3 and H3K9me3 that permit stoichiometric modification). Consistent with the peptide microarray analysis, ATRX-ADD only exhibited binding with H3Kc9me3-containing nucleosomes (Additional file 6: Figure S4). Furthermore, using standard Western blotting protocols, we demonstrate that ATRX-ADD can function as a highly specific probe to detect endogenous H3K9me3 from cell lysate, showing sufficient affinity and specificity that rivals the corresponding antibody (Additional file 7: Figure S5).
We then evaluated whether the immobilized reader domains (MARCC approach) (Figure 3) captured unique chromatin states in native mononucleosomes from MCF-7 cells (Additional file 8: Figure S6). In this application, HaloTag provided an effective means of covalently linking the reader domains to a solid support, allowing extensive washing . By Western blot analysis, MARCC revealed specific enrichment with each reader domain (Figure 4A,B): ING2-PHD enriched for H3K4me3-containing nucleosomes, AIRE-PHD discriminated against H3K4me3, and ATRX-ADD enriched H3K9me3-containing chromatin. To further define the chromatin states captured by MARCC, sufficient quantities of histone protein are required for quantitative mass spectrometry (qMS) while maintaining the native state of nucleosomes during enrichment. For this purpose, resin-captured nucleosomes from MARCC were released with tobacco etch virus (TEV) protease (Additional file 9: Figure S7) and analyzed by qMS. Owing to the low abundance of H3K4me3 (<0.1%) in the input chromatin, ING2-PHD MARCC enriched from ≈ 107 cells was below the detection limit for qMS, though the high enrichment was evident by Western quantification and qPCR (Figure 4C).
Quantitative mass spectrometry (qMS) revealed unique nucleosomal PTM signatures of MARCC-enriched chromatin with excellent reproducibility (Additional file 10: Figure S8 B-D). Consistent with their direct binding sites, MARCC with AIRE-PHD enriched H3K4unmod-containing chromatin to 90% and MARCC with ATRX-ADD enriched both unmodified H3K4 and hypermethylated H3K9 to more than 90% (Figure 4). Importantly, we were able to identify and quantify other PTM states that were co-enriched with the targeted PTM states (Figure 4 and Additional file 11: Table S1). Specifically, ATRX-ADD, but not AIRE-PHD, enriched H3K79unmod to over 90% and H4K20me2/3 to ≈ 80%, implying that these two modifications are more strongly associated with the combination of H3K4unmod and H3K9me2/3, but not H3K4unmod alone. In other words, ATRX-ADD captured specific heterochromatin state with four coexisting PTM signatures (H3K4unmod, H3K9me2/3, H4K20me2/3 and H3K79unmod), all of which were present at >80% abundance, suggesting that most of these states exist on both histone copies within the mononucleosome. Previous ChIP-sequencing data have correlated H4K20 and H3K9 hypermethylation . Here, our data suggest that these two modifications exist within the same nucleosome. One possibility is that the interaction between the H4K20 methyltransferase Suv4-20 and H3K9me3-binding protein HP1  bridges the methyltransferase for intranucleosomal H4K20 methylation. H4K20 methylation is located near the basic batch of H4  and might interfere with Dot1 (H3K79 methyltransferase), explaining the co-enrichment of H4K20 hypermethylation with unmodified H3K79. Still, the functional importance of these four PTMs coexisting at the mononucleosome level needs further investigation. We speculate that this combination might be enriched for a particular subtype of heterochromatin.
To support the hypothesis that ATRX-ADD enriched unique silent chromatin, relative enrichment of additional PTM signatures was quantified by qMS (Figure 5A). Similar analyses were performed with AIRE-PHD from the MARCC enrichment (Figure 5B). We observed positive enrichment of repressive marks by ATRX-ADD, including H3K9 hypermethylation, H3K27 hypermethylation and H4K20 hypermethylation (Figure 5A, above x-axis). Consistent with silent chromatin, active marks such as H4 hyperacetylation and H3K36me3 were depleted by ATRX-ADD (Figure 5A, below x-axis). MARCC also provided information on less characterized PTMs, for example, methylations at H3K79 and H3K18. H3K79 methylation was dramatically depleted using MARCC with ATRX-ADD, indicating a negative correlation with this specific heterochromatin state. Strong enrichment of H3K18me1 by ATRX-ADD and AIRE-PHD (Figure 5A,B), which both strongly prefer H3K4unmod, suggests a role of this PTM in connection with H3K4unmod.
Finally, qPCR analysis of MARCC-eluted chromatin revealed highly enriched heterochromatic DNA by ATRX-ADD (Figure 5C). Collectively, these results suggest that the combination of four coexisting nucleosomal PTMs enriched by ATRX-ADD represents a heterochromatic nucleosome species. The specificity of the enrichment and the combinatorial PTM signatures captured using MARCC with ATRX-ADD corroborates the application of MARCC as a tool to enrich, identify and quantify chromatin states with combinatorial PTM patterns and associated DNA.
In summary, we demonstrate an efficient, quantitative, reader-based platform for interrogating epigenetic states of chromatin (Table 1). The workflow involves (i) using a combinatorial PTM histone peptide microarray to define reader specificity and suitability for MARCC and (ii) utilization of MARCC to enrich, identify and quantify unique chromatin states. Utilizing this workflow, we describe three different reader domains that can recognize combinatorial PTM patterns and enrich unique chromatin states. Histone PTMs and their corresponding DNA associated with specific chromatin states are comprehensively profiled for as few as 107 cells. Future systematic analyses of uncharacterized or genetically engineered readers will expand the capacity of MARCC to reveal undiscovered chromatin states. We envision that reader-based epigenetic platforms can function as a new paradigm for customizable chromatin enrichment reagents, avoiding issues of antibody variability and specificity.
Plasmids encoding mouse ING2-PHD (201 to 281) and human AIRE-PHD1 (293 to 354) were kindly provided by T Kutateladze (University of Colorado-Denver) and G Musco (Dulbecco Telethon Institute), respectively. cDNA encoding human ATRX-ADD (163 to 292) was synthesized (Integrated DNA Technologies, Inc.). A bacterial expression vector encoding an N-terminal (HQ)5-HaloTag followed by a TEV protease cleavage site was kindly supplied by Dr. M Slater and Dr. J Hartnett (Promega). A similar vector (pFN29) is now commercially available (Promega). ING2-PHD (200 to 280), AIRE-PHD (293 to 354) and ATRX-ADD (163 to 292) were amplified from the parent vectors and cloned into the (HQ)5-HaloTag N-terminal vector with SgfI and PmeI sites.
Protein expression and purification
All recombinant (HQ)5-HaloTag fusion proteins and (HQ)5-HaloTag alone were expressed in BL21 (DE3) pLysS cells. When OD600nm reached 0.6 to 0.8, expression was induced by addition of 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 20 to 25°C for 4 hours or at 18°C overnight. 150 μM ZnCl2 was supplemented into the media after induction. The cell pellets were resuspended in 30 mM HEPES, 500 mM NaCl, pH 7.4, with 1 mM dithiotheritol (DTT), 1 mM phenylmethylsulfonyl fluoride (PMSF), 10 μg/ml leupeptin, 10 μg/ml aprotinin and 1 mg/ml lysozyme. After agitation at 4°C for 30 min, the cell resuspension was further lysed by three mild sonication cycles (20% amplitude; 5 s on; 5 s off; 1 min total per sonication cycle) with 1/4″ probe (Thermo Fisher). Clarified lysate supernatant was incubated with Ni-NTA resin (GE Life Sciences) in batches for 2 h, 4°C. Resin-bound proteins were washed three times with 30 mM HEPES, 150 mM NaCl, 20 mM imidazole, pH 7.4, followed by batch elution with 30 mM HEPES, 150 mM NaCl, 300 mM imidazole, pH 7.4. Eluents were dialyzed into 30 mM HEPES, 150 mM NaCl, pH 7.4, supplemented with 3 mM DTT and up to 10% (v/v) glycerol. After being concentrated in 10 kDa molecular weight cut-off (MWCO) Amicon (Millipore), protein was quantified by Bradford assay and analyzed on 12% SDS-PAGE followed by Coomassie staining to check protein purity and integrity (Additional file 1: Figure S1A). Proteins were aliquoted and stored at -80°C.
Combinatorial PTM histone peptide microarray synthesis
Peptides were synthesized on modified cellulose with a Respep SL automated synthesizer (Intavis AG, Köln, Germany). A spotting control cellulose-Cy3 dye conjugate was created by simply removing Fmoc from the celluspot and coupling Cy3 (Lumiprobe, corp.). A blank celluspot control was created by deprotecting the Fmoc of the spot and subjecting the paper to the dissolution conditions. For peptides, the initial coupling consisted of equimolar amounts of Fmoc-Ala-OH and Boc-Ala-OH to decrease resin loading and improve synthesis quality. Each synthesis included peptides with an acid cleavable rink linker for quality assessment by HPLC and mass spectrometry (Additional file 12: Table S2). A 20-atom polyethylene linker (PE) was also included (Novabiochem), which improved synthesis and protein binding in the array. A streptavidin control was created by coupling Fmoc-Glu-(biotinyl-PEG)-OH (Novabiochem). Standard Fmoc/tBu chemistry was used to create the 13-amino-acid peptides. Peptide PTMs were introduced as suitably protected monomers, as described previously . On completion of the synthesis, the side chain protecting groups were removed with 82.5% trifluoroacetic acid (TFA): 5% thioanisole: 5% water: 5% saturated phenol in dichloromethane: 2.5% ethane dithiol for 90 minutes (150 μl per spot). The solution was removed and replaced with 250 μl of 88.5% TFA, 4% trifluoromethanesulfonic acid, 2.5% triisopropylsilane, 5% water and agitated gently overnight to dissolve the cellulose. Peptide cellulose conjugates were precipitated twice with cold diethyl ether, allowed to dry for several minutes, and redissolved into 100 μl dimethyl sulfoxide (DMSO).
The arrays were generated by combining 75% cellulose peptide DMSO solution with 10% cy3 cellulose and 15% 7× saline-sodium citrate (SSC) buffer (1.05 M NaCl, 105 mM sodium citrate). The biotin positive control was created by using 7.5% of the DMSO peptide stock solution with 67.5% DMSO in addition to the Cy3 cellulose conjugate and SSC buffer. Each solution was mixed in separate wells of 384-well polypropylene plates and spotted in triplicate onto 75 mm × 25 mm nitrocellulose sides (Intuitive Biosciences) with a Gene Machines OmniGrid Arrayer (Genomic Solutions). The dip time was 0 ms and the print time varied between 30 and 50 ms. Humidity and room temperature were controlled between 40% to 60% relative humidity and 20°C to 25°C. A test microarray was printed and tested with HaloTag only or streptavidin only as a control to identify non-specific binding peptides, which were removed for subsequent microarrays (summarized in Additional file 2 and Additional file 12, Table S2). The final library contained 746 distinct histone peptide species. On each nitrocellulose slide, two identical libraries were printed, with each library containing 16 blocks and an overall number of 2,337 spots. In each block, blank cellulose spots and biotin peptide spots were included as negative and positive controls. The physical map of the peptide array is described in Additional file 2.
Peptide array binding assay
Peptide arrays were handled with care and protected from light. All solutions were filtered before putting on the array. Assays were performed with a modified two-chamber simplex gasket (Intuitive Bioscience). After tightening the gasket onto the slide surface, the peptide array was blocked with blocking solution (1 × PBS, 0.05% Tween-20, pH 7.4, 1% BSA) at 4°C overnight to reduce non-specific binding.
For reader binding assay
Purified recombinant (HQ)5-HaloTag proteins (10 nM to 1 μM) were labeled with excessive HaloTag biotin ligand (Promega) in blocking solution and incubated with the peptide array under mild rocking at 4°C for 1 hour. The slide was washed with PBS with 0.05% Tween (PBST) three times and incubated with 1:2,000 streptavidin-conjugated Alexa-Fluor647 (Invitrogen) in blocking solution at room temperature for 1 hour.
For antibody binding assay
Different dilutions (1:10,000 to 1:2,000) of rabbit primary polyclonal antibodies (H3K9me3 antibody ab8898 lot#GR102573-1, H3K4me3 antibody ab8580 lot#GR56122-1) in blocking solution were incubated with the array under mild rocking at 4°C for 1 hour. The slide was washed with PBST three times and incubated with 1:1,000 anti-rabbit IgG Alexa 647 (Cell Signaling) in blocking solution at room temperature for 1 hour.
For both assays
After three washes in PBST and a final wash in distilled water, the slide was dried by centrifugation and imaged at dual wavelengths of 532 nm and 635 nm on Axon GenePix 4000B (Molecular Devices). The laser power was set to 100%, with automatic gain adjustment (0.05% saturation tolerance) for dual photomultipliers. Image was obtained at 5 μm pixel resolution. Features in each block were defined by manual adjustment of 13 × 13 grid (feature diameter, 280 μm; column spacing and row spacing, 320 μm) to cover every spot. Signal intensities were quantified by GenePix Pro 6.1 software (Molecular Devices). For each spot (feature), the mean intensities for 635 nm wavelength were used for subsequent analysis. For each peptide species, an average was calculated from three replicate spots. The averaged intensities were normalized to the range between 0 and 1 by the lowest and highest values of selected peptides on each library. For comparison between peptides with combinatorial modifications, the intensities were normalized to the peptide with single modification as log2 ratio. Heat maps were generated by Java Treeview version 1.1. The signal at 532 nm wavelength was used to identify misprinting events. For a complete dataset, see Additional file 3.
Nucleosome reconstitution and MLA (methyl lysine analog) production
Methyl lysine analog (MLA), H3Kc4me3 (Kc = aminoethyl cysteine), was prepared according to the literature . Briefly, Histone H3 (Xenopus laevis) with K4C and C110A mutations was expressed in Escherichia coli and purified by Superdex200 column followed by ion-exchange columns. The purified histone was reduced, and the cysteine residue was alkylated with excess (2-bromoethyl) trimethylammonium bromide (Sigma Aldrich) at elevated temperature. After the reaction was quenched with 2-mercaptoethanol, excess reagents and salts were removed using a PD-10 desalting device. Along with other core histones (H2A, H2B, and H4) that are recombinantly prepared, H3Kc4me3 MLA was refolded into histone octamers. The purified octamer and 146-bp DNA fragments (shown to exhibit strong positioning to the histone octamers) were reconstituted into nucleosome core particle by a salt-gradient dialysis method described previously .
Reader-probe Western blot
HEK293 cells were lysed in radioimmunoprecipitation assay (RIPA) buffer with protease inhibitors, nicotinamide and trichostatin A. The protein level was quantified by Bradford assay. 10 μg total cell lysates were separated using 12% SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) membrane. After blocking with 5% BSA for 1 hour at room temperature, the membrane was incubated with 100 nM HaloTag ATRX-ADD at 4° for 3 hours. Before incubation, ATRX-ADD was reacted with HaloTag ligand (biotin conjugate or AlexaFluor 660 conjugate) (Promega). For biotin conjugate, the membrane was further incubated with streptavidin-AlexaFluor647 conjugate (Invitrogen) at room temperature for 1 hour. After three washes, the membrane was directly exposed at Cy5 setting (GE ImageQuant LAS 4000).
Native mononucleosome isolation from MCF-7 cells
MCF-7 cells were cultured in DMEM supplemented with 10% FBS. For each nucleosome isolation, MCF-7 cells at ≈ 90% confluency from two 10-cm plates (≈2 × 107 cells in total) were collected and washed for three times in ice-cold Buffer M (10 mM HEPES, 10 mM KCl, 1.5 mM MgCl2, 340 mM sucrose, pH 7.9, 10% glycerol, v/v), supplemented with 1 μg/ml trichostatin A, 1 mM DTT, 0.5 mM PMSF, 10 mM β-glycerophosphate, 1 mM leupeptin, and 1 mM aprotinin. At the last wash, the cell resuspension was lysed with 0.1% Triton X-100 on ice for 10 min. After lysis, nuclei pellets were resuspended in Buffer M and centrifuged at 1,300 g for 12 min through chilled sucrose cushion buffer (10 mM HEPES, pH 7.9, 30% sucrose, w/v, 1.5 mM MgCl2) to further purify the nuclei pellets. After three washes with Buffer M, nuclei was diluted to 1.2 to 1.6 mg/ml DNA concentration and digested with 2,000 gel units of micrococcal nuclease (New England Biolabs) at 37°C for 12 min with constant mixing in the presence of 1 mM final concentration of CaCl2. Prior to assay, the amount of enzyme and digestion time was optimized to obtain above 90% purity of mononucleosomes. MNase activity was stopped with 10 mM ethylenediaminetetraacetic acid (EDTA) and spun down. Soluble chromatin from the supernatant (S1) was collected. Less soluble chromatin (S2) was recovered from the nuclei pellets resuspended in 5 mM HEPES, 0.2 mM EDTA at 4°C overnight. Pooled chromatin extract (S1 plus S2) was concentrated to ≈ 10 μM and dialyzed into 30 mM HEPES, 150 mM NaCl, pH 7.4, 10% (v/v) glycerol. Nucleosomes were analyzed on 1.2% agarose gel (with 0.1% sodium dodecyl sulfate (SDS) in the sample) and 18% SDS-PAGE to check the size and quantity of DNA and histones separately (Additional file 1: Figure S1). A yield of ≈ 2 nmole mononucleosomes (or 200 μg DNA) was usually obtained from a single isolation.
MARCC (matrix-assisted reader chromatin capture)
The MARCC resins were prepared by incubation of saturating amounts (more than 10 nmole) of purified recombinant (HQ)5-HaloTag proteins with 200 μl HaloLink resin slurry (Promega) in MARCC buffer (30 mM HEPES, 150 mM NaCl, pH 7.4, 0.01% NP-40, 10% glycerol) at 4°C for 1 hour. (HQ)5-HaloTag protein alone was included as a negative control. Excessive proteins were removed by three brief washes with MARCC buffer. Chromatin capture was achieved by incubating 1 nmole native mononucleosomes or reconstituted nucleosomes with MARCC resins under constant rotation at 4°C overnight. Bound chromatin was further washed with MARCC buffer and eluted by Halo-TEV protease (Promega) cleavage in the presence of 150 μl 1 mM Tris-HCl, pH 7.4 at room temperature for 2 hours. Eluted chromatin was combined with another 150 μl resin-recovered chromatin and could be used for downstream analysis. For small-scale enrichment, 100 pmole nucleosomes and 20 μl MARCC resin were used.
Native chromatin immunoprecipitation (NChIP)
100 pmole native mononucleosome extract was mixed with 2 μg antibody (H3K4me3 - Ab8580, lot# GR56122-1; H3K9me3 - Ab8898, lot# GR102573-1; normal rabbit IgG - sc2027) and incubated with 20 μL Magna ChIP Protein G magnetic beads (Millipore) in buffer N (30 mM HEPES, 150 mM NaCl, pH 7.4, 0.01% NP-40, 5 mM EDTA) under constant rotation at 4°C overnight. Bound chromatin was further washed with buffer N three times and eluted twice with 150 μl elution buffer containing 1% SDS.
Quantitative mass spectrometry for histone PTMs
Procedures for quantitative mass spectrometry sample preparation and data analysis were as previous described [42, 43]. Briefly, histones were subjected to chemical derivatization using propionic anhydride (Sigma-Aldrich) and digested with sequencing grade trypsin (Promega) at a 10:1 substrate to enzyme ratio for 6 hours at 37°C. The digested peptides were treated with an additional round of propionylation for the purpose of adding propionyl group to the newly generated N-terminus. Peptides were desalted using C18 extracted mini disk (Empore 3 M, MN, USA) and dissolved in 0.1% formic acid. Approximately 1 μg of each sample was loaded via an autosampler (EASY-nLC, Thermo Fisher Scientific Inc.) onto a Thermo Scientific Acclaim PepMap 100 pre-column (75 μm × 2 cm, C18 resin, 3 μm particle size, 100 Å pore sizes). Peptides were chromatographically resolved via an Acclaim PepMap rapid separation liquid chromatography (RSLC) analytical column (50 μm × 15 cm, C18 resin, 2 μm particle sizes, 100 Å pore sizes), using a 67-min 1-98% solvent B gradient (solvent A = 0.1% formic acid, solvent B = 100% acetonitrile) at a flow rate of 300 nl/min. The eluted peptides were electrosprayed into and detected by a Q Exactive mass spectrometer (Thermo Fisher Scientific Inc.) with a resolution of 70,000 for full MS spectrum followed by MS/MS spectra obtained in a higher-energy collisional dissociation cell. The relative abundance of each modification, expressed as a percentage on a histone peptide sequence, was quantified by analyzing its MS and MS/MS spectra via a program developed in-house. All results were also manually verified.
DNA isolation and real-time PCR
DNA was isolated from chromatin enriched by MARCC or NChIP along with 10% input chromatin using GeneJet PCR purification kit (Fermentas) in the presence of isopropanol. Real-time PCR was performed in duplicate with SsoFast EvaGreen supermix (Bio-Rad) on CFX96 (Bio-rad). For each specific PCR primer set, a standard curve was created from a serial dilution covering 0.2 to 5% of input DNA. For linear regression, Ct (cycle threshold) values were plotted with log starting quantity. Primer sets with R2 over 0.99 and amplification efficiency over 90% were chosen
Nuclease-free water was used as a no-template control to monitor contaminants interfering with reactions. To quantify a specific DNA fragment in enriched chromatin samples, individual Ct values were converted into starting DNA quantity using corresponding standard curve and the average was calculated from two individual MARCC or ChIP experiments.
Primer sets for qPCR
Chr1 SAT2-forward: 5′-CATCGATGGAAATGAAAGGAGTC-3′;
Chr1 SAT2-reverse: 5′- ACCATTGGATGATTGCAGTCAA-3′.
alpha-thalassemia/mental retardation, X-linked
bovine serum albumin
Dulbecco’s modified Eagle’s medium
fetal bovine serum
histone H3 trimethylated at lysine 4
histone H3 trimethylated at lysine 9
high-performance liquid chromatography
inhibitor of growth
matrix-assisted reader chromatin capture
Michigan Cancer Foundation-7
methyl lysine analog
molecular weight cut-off
native chromatin immunoprecipitation
phosphate-buffered saline with Tween
polymerase chain reaction
quantitative mass spectrometry
quantitative polymerase chain reaction
rapid separation liquid chromatography
sodium dodecyl sulfate
tobacco etch virus
Strahl BD, Allis CD: The language of covalent histone modifications. Nature. 2000, 403 (6765): 41-45. 10.1038/47412.
Voigt P, LeRoy G, Drury WJ, Zee BM, Son J, Beck DB, Young NL, Garcia BA, Reinberg D: Asymmetrically modified nucleosomes. Cell. 2012, 151 (1): 181-193. 10.1016/j.cell.2012.09.002.
Fischle W, Tseng BS, Dormann HL, Ueberheide BM, Garcia BA, Shabanowitz J, Hunt DF, Funabiki H, Allis CD: Regulation of HP1-chromatin binding by histone H3 methylation and phosphorylation. Nature. 2005, 438 (7071): 1116-1122. 10.1038/nature04219.
Young NL, DiMaggio PA, Plazas-Mayorca MD, Baliban RC, Floudas CA, Garcia BA: High throughput characterization of combinatorial histone codes. Mol Cell Proteomics. 2009, 8 (10): 2266-2284. 10.1074/mcp.M900238-MCP200.
Phanstiel D, Brumbaugh J, Berggren WT, Conard K, Feng X, Levenstein ME, McAlister GC, Thomson JA, Coon JJ: Mass spectrometry identifies and quantifies 74 unique histone H4 isoforms in differentiating human embryonic stem cells. Proc Natl Acad Sci U S A. 2008, 105 (11): 4093-4098. 10.1073/pnas.0710515105.
Fuchs SM, Krajewski K, Baker RW, Miller VL, Strahl BD: Influence of combinatorial histone modifications on antibody and effector protein recognition. Curr Biol. 2011, 21 (1): 53-58. 10.1016/j.cub.2010.11.058.
Bock I, Dhayalan A, Kudithipudi S, Brandt O, Rathert P, Jeltsch A: Detailed specificity analysis of antibodies binding to modified histone tails with peptide arrays. Epigenetics. 2011, 6 (2): 256-263. 10.4161/epi.6.2.13837.
Bua DJ, Kuo AJ, Cheung P, Liu CL, Migliori V, Espejo A, Casadio F, Bassi C, Amati B, Bedford MT, Guccione E, Gozani O: Epigenome microarray platform for proteome-wide dissection of chromatin-signaling networks. PLoS One. 2009, 4 (8): e6789-10.1371/journal.pone.0006789.
Egelhofer TA, Minoda A, Klugman S, Lee K, Kolasinska-Zwierz P, Alekseyenko AA, Cheung MS, Day DS, Gadel S, Gorchakov AA, Gu T, Kharchenko PV, Kuan S, Latorre I, Linder-Basso D, Luu Y, Ngo Q, Perry M, Rechtsteiner A, Riddle NC, Schwartz YB, Shanower GA, Vielle A, Ahringer J, Elgin SC, Kuroda MI, Pirrotta V, Ren B, Strome S, Park PJ: An assessment of histone-modification antibody quality. Nat Struct Mol Biol. 2011, 18 (1): 91-93. 10.1038/nsmb.1972.
Heubach Y, Planatscher H, Sommersdorf C, Maisch D, Maier J, Joos TO, Templin MF, Poetz O: From spots to beads-PTM-peptide bead arrays for the characterization of anti-histone antibodies. Proteomics. 2013, 13 (6): 1010-1015. 10.1002/pmic.201200383.
Nishikori S, Hattori T, Fuchs SM, Yasui N, Wojcik J, Koide A, Strahl BD, Koide S: Broad ranges of affinity and specificity of anti-histone antibodies revealed by a quantitative peptide immunoprecipitation assay. J Mol Biol. 2012, 424 (5): 391-399. 10.1016/j.jmb.2012.09.022.
Hattori T, Taft JM, Swist KM, Luo H, Witt H, Slattery M, Koide A, Ruthenburg AJ, Krajewski K, Strahl BD, White KP, Farnham PJ, Zhao Y, Koide S: Recombinant antibodies to histone post-translational modifications. Nat Methods. 2013, 10 (10): 992-995. 10.1038/nmeth.2605.
Fuchs SM, Strahl BD: Antibody recognition of histone post-translational modifications: emerging issues and future prospects. Epigenomics. 2011, 3 (3): 247-249. 10.2217/epi.11.23.
Ernst J, Kellis M: Discovery and characterization of chromatin states for systematic annotation of the human genome. Nat Biotechnol. 2010, 28 (8): 817-825. 10.1038/nbt.1662.
Kharchenko PV, Alekseyenko AA, Schwartz YB, Minoda A, Riddle NC, Ernst J, Sabo PJ, Larschan E, Gorchakov AA, Gu T, Linder-Basso D, Plachetka A, Shanower G, Tolstorukov MY, Luquette LJ, Xi R, Jung YL, Park RW, Bishop EP, Canfield TK, Sandstrom R, Thurman RE, MacAlpine DM, Stamatoyannopoulos JA, Kellis M, Elgin SC, Kuroda MI, Pirrotta V, Karpen GH, Park PJ: Comprehensive analysis of the chromatin landscape in Drosophila melanogaster. Nature. 2011, 471 (7339): 480-485. 10.1038/nature09725.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, Ku M, Durham T, Kellis M, Bernstein BE: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473 (7345): 43-49. 10.1038/nature09906.
Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, Lee W, Mendenhall E, O’Donovan A, Presser A, Russ C, Xie X, Meissner A, Wernig M, Jaenisch R, Nusbaum C, Lander ES, Bernstein BE: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007, 448 (7153): 553-560. 10.1038/nature06008.
Zhu J, Adli M, Zou JY, Verstappen G, Coyne M, Zhang X, Durham T, Miri M, Deshpande V, De Jager PL, Bennett DA, Houmard JA, Muoio DM, Onder TT, Camahort R, Cowan CA, Meissner A, Epstein CB, Shoresh N, Bernstein BE: Genome-wide chromatin state transitions associated with developmental and environmental cues. Cell. 2013, 152 (3): 642-654. 10.1016/j.cell.2012.12.033.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, Jaenisch R, Wagschal A, Feil R, Schreiber SL, Lander ES: A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006, 125 (2): 315-326. 10.1016/j.cell.2006.02.041.
van Rossum B, Fischle W, Selenko P: Asymmetrically modified nucleosomes expand the histone code. Nat Struct Mol Biol. 2012, 19 (11): 1064-1066. 10.1038/nsmb.2433.
Peach SE, Rudomin EL, Udeshi ND, Carr SA, Jaffe JD: Quantitative assessment of chromatin immunoprecipitation grade antibodies directed against histone modifications reveals patterns of co-occurring marks on histone protein molecules. Mol Cell Proteomics. 2012, 11 (5): 128-137. 10.1074/mcp.M111.015941.
Soldi M, Bonaldi T: The proteomic investigation of chromatin functional domains reveals novel synergisms among distinct heterochromatin components. Mol Cell Proteomics. 2013, 12 (3): 764-780. 10.1074/mcp.M112.024307.
Musselman CA, Lalonde ME, Cote J, Kutateladze TG: Perceiving the epigenetic landscape through histone readers. Nat Struct Mol Biol. 2012, 19 (12): 1218-1227. 10.1038/nsmb.2436.
Garske AL, Oliver SS, Wagner EK, Musselman CA, LeRoy G, Garcia BA, Kutateladze TG, Denu JM: Combinatorial profiling of chromatin binding modules reveals multisite discrimination. Nat Chem Biol. 2010, 6 (4): 283-290. 10.1038/nchembio.319.
Vermeulen M, Mulder KW, Denissov S, Pijnappel WW, van Schaik FM, Varier RA, Baltissen MP, Stunnenberg HG, Mann M, Timmers HT: Selective anchoring of TFIID to nucleosomes by trimethylation of histone H3 lysine 4. Cell. 2007, 131 (1): 58-69. 10.1016/j.cell.2007.08.016.
Moore KE, Carlson SM, Camp ND, Cheung P, James RG, Chua KF, Wolf-Yadlin A, Gozani O: A general molecular affinity strategy for global detection and proteomic analysis of lysine methylation. Mol Cell. 2013, 50 (3): 444-456. 10.1016/j.molcel.2013.03.005.
Eustermann S, Yang JC, Law MJ, Amos R, Chapman LM, Jelinska C, Garrick D, Clynes D, Gibbons RJ, Rhodes D, Higgs DR, Neuhaus D: Combinatorial readout of histone H3 modifications specifies localization of ATRX to heterochromatin. Nat Struct Mol Biol. 2011, 18 (7): 777-782. 10.1038/nsmb.2070.
Iwase S, Xiang B, Ghosh S, Ren T, Lewis PW, Cochrane JC, Allis CD, Picketts DJ, Patel DJ, Li H, Shi Y: ATRX ADD domain links an atypical histone methylation recognition mechanism to human mental-retardation syndrome. Nat Struct Mol Biol. 2011, 18 (7): 769-776. 10.1038/nsmb.2062.
Pena PV, Davrazou F, Shi X, Walter KL, Verkhusha VV, Gozani O, Zhao R, Kutateladze TG: Molecular mechanism of histone H3K4me3 recognition by plant homeodomain of ING2. Nature. 2006, 442 (7098): 100-103.
Shi X, Hong T, Walter KL, Ewalt M, Michishita E, Hung T, Carney D, Pena P, Lan F, Kaadige MR, Lacoste N, Cayrou C, Davrazou F, Saha A, Cairns BR, Ayer DE, Kutateladze TG, Shi Y, Côté J, Chua KF, Gozani O: ING2 PHD domain links histone H3 lysine 4 methylation to active gene repression. Nature. 2006, 442 (7098): 96-99.
Chakravarty S, Zeng L, Zhou MM: Structure and site-specific recognition of histone H3 by the PHD finger of human autoimmune regulator. Structure. 2009, 17 (5): 670-679. 10.1016/j.str.2009.02.017.
Chignola F, Gaetani M, Rebane A, Org T, Mollica L, Zucchelli C, Spitaleri A, Mannella V, Peterson P, Musco G: The solution structure of the first PHD finger of autoimmune regulator in complex with non-modified histone H3 tail reveals the antagonistic role of H3R2 methylation. Nucleic Acids Res. 2009, 37 (9): 2951-2961. 10.1093/nar/gkp166.
Bock I, Kudithipudi S, Tamas R, Kungulovski G, Dhayalan A, Jeltsch A: Application of Celluspots peptide arrays for the analysis of the binding specificity of epigenetic reading domains to modified histone tails. BMC Biochem. 2011, 12: 48-10.1186/1471-2091-12-48.
Rothbart SB, Dickson BM, Ong MS, Krajewski K, Houliston S, Kireev DB, Arrowsmith CH, Strahl BD: Multivalent histone engagement by the linked tandem Tudor and PHD domains of UHRF1 is required for the epigenetic inheritance of DNA methylation. Genes Dev. 2013, 27 (11): 1288-1298. 10.1101/gad.220467.113.
Cai L, Rothbart SB, Lu R, Xu B, Chen WY, Tripathy A, Rockowitz S, Zheng D, Patel DJ, Allis CD, Strahl BD, Song J, Wang GG: An H3K36 methylation-engaging Tudor motif of polycomb-like proteins mediates PRC2 complex targeting. Mol Cell. 2013, 49 (3): 571-582. 10.1016/j.molcel.2012.11.026.
Los GV, Encell LP, McDougall MG, Hartzell DD, Karassina N, Zimprich C, Wood MG, Learish R, Ohana RF, Urh M, Simpson D, Mendez J, Zimmerman K, Otto P, Vidugiris G, Zhu J, Darzins A, Klaubert DH, Bulleit RF, Wood KV: HaloTag: a novel protein labeling technology for cell imaging and protein analysis. ACS Chem Biol. 2008, 3 (6): 373-382. 10.1021/cb800025k.
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009.
Schotta G, Lachner M, Sarma K, Ebert A, Sengupta R, Reuter G, Reinberg D, Jenuwein T: A silencing pathway to induce H3-K9 and H4-K20 trimethylation at constitutive heterochromatin. Genes Dev. 2004, 18 (11): 1251-1262. 10.1101/gad.300704.
Fingerman IM, Li HC, Briggs SD: A charge-based interaction between histone H4 and Dot1 is required for H3K79 methylation and telomere silencing: identification of a new trans-histone pathway. Genes Dev. 2007, 21 (16): 2018-2029. 10.1101/gad.1560607.
Simon MD, Chu F, Racki LR, de la Cruz CC, Burlingame AL, Panning B, Narlikar GJ, Shokat KM: The site-specific installation of methyl-lysine analogs into recombinant histones. Cell. 2007, 128 (5): 1003-1012. 10.1016/j.cell.2006.12.041.
Luger K, Rechsteiner TJ, Richmond TJ: Preparation of nucleosome core particle from recombinant histones. Methods Enzymol. 1999, 304: 3-19.
Plazas-Mayorca MD, Zee BM, Young NL, Fingerman IM, LeRoy G, Briggs SD, Garcia BA: One-pot shotgun quantitative mass spectrometry characterization of histones. J Proteome Res. 2009, 8 (11): 5367-5374. 10.1021/pr900777e.
Garcia BA, Mollah S, Ueberheide BM, Busby SA, Muratore TL, Shabanowitz J, Hunt DF: Chemical derivatization of histones for facilitated analysis by mass spectrometry. Nat Protoc. 2007, 2 (4): 933-938. 10.1038/nprot.2007.106.
Leroy G, Dimaggio PA, Chan EY, Zee BM, Blanco MA, Bryant B, Flaniken IZ, Liu S, Kang Y, Trojer P, Garcia BA: A quantitative atlas of histone modification signatures from human cancer cells. Epigenetics Chromatin. 2013, 6 (1): 20-10.1186/1756-8935-6-20.
We are grateful to former laboratory members E Wagner and B Albaugh for cloning the corresponding constructs. We thank R Sridharan, X Zhong, P Lewis and E Wagner (University of Wisconsin-Madison) for helpful discussion and critical reading of the manuscript. We also thank Gene Expression Center (University of Wisconsin Biotechnology Center) for technical assistance with peptide array printing and scanning and W Xu (University of Wisconsin-Madison) for sharing MCF-7 cells. This research is supported by NIH R01 grant #2R37GM059785-15/P250VA (JMD), National Science Foundation Early Faculty CAREER award (BAG) and an NIH Innovator grant (DP2OD007447) from the Office of the Director (BAG).
The authors declare that they have no competing interests.
ZS and JMD conceived and designed the overall study. SSO, ZS, MDB and JHL designed the histone peptide array layout. MDB synthesized the peptides and printed the peptide arrays. ZS carried out the experiments and analyzed the results. JHL synthesized the MLA nucleosomes. SL and BAG performed the qMS analysis. ZS and JMD wrote the manuscript. All authors read and approved the final manuscript.