Top-down and Middle-down Protein Analysis Reveals that Intact and Clipped Human Histones Differ in Post-translational Modification Patterns*

Post-translational modifications (PTMs) of histone proteins play a fundamental role in regulation of DNA-templated processes. There is also growing evidence that proteolytic cleavage of histone N-terminal tails, known as histone clipping, influences nucleosome dynamics and functional properties. Using top-down and middle-down protein analysis by mass spectrometry, we report histone H2B and H3 N-terminal tail clipping in human hepatocytes and demonstrate a relationship between clipping and co-existing PTMs of histone H3. Histones H2B and H3 undergo proteolytic processing in primary human hepatocytes and the hepatocellular carcinoma cell line HepG2/C3A when grown in spheroid (3D) culture, but not in a flat (2D) culture. Using tandem mass spectrometry we localized four different clipping sites in H3 and one clipping site in H2B. We show that in spheroid culture clipped H3 proteoforms are mainly represented by canonical histone H3, whereas in primary hepatocytes over 90% of clipped H3 correspond to the histone variant H3.3. Comprehensive analysis of histone H3 modifications revealed a series of PTMs, including K14me1, K27me2/K27me3, and K36me1/me2, which are differentially abundant in clipped and intact H3. Analysis of co-existing PTMs revealed negative crosstalk between H3K36 methylation and H3K23 acetylation in clipped H3. Our data provide the first evidence of histone clipping in human hepatocytes and demonstrate that clipped H3 carry distinct co-existing PTMs different from those in intact H3.

Chromatin is a highly dynamic structure that must respond to different stimuli in order to orchestrate all DNA-dependent processes. Post translational modifications (PTMs) 1 of histones play a major role in regulation of chromatin functionality. Evidence is emerging that not only "classical" histone PTMs, such as methylation, acetylation, and phosphorylation at distinct residues, but also proteolytic processing of nucleosome proteins, known as "histone clipping," can be involved in regulation of key cellular processes, such as transcriptional regulation, cell differentiation, and senescence (1)(2)(3)(4)(5)(6)(7).
Clipping of the histone H3 N-terminal tail was reported to be associated with gene activation in yeast. Santos-Rosa et al. demonstrated a serine protease activity in S. cerevisiae that cleaves histone H3 after residue Ala21 (A21) during sporulation and stationary phase (1). H3 clipping took place specifically within the promoters of sporulation-induced genes following the induction of transcription and prior to histone eviction from these DNA regions. Prevention of H3 N-tail cleavage by amino acid substitution at the endoproteinase recognition site (H3 Q19A, L20A) abolished expression of these genes, indicating that H3 clipping is essential for productive transcription.
The biological significance of histone clipping in higher eukaryotes is not yet understood but also appears to be related to functional commitment by the cell. Duncan et al. demonstrated that histone H3 is proteolytically cleaved by the enzyme Cathepsin L1 (CTSL1) at several sites between residues A21 and S28 during mouse ESC differentiation (5). The "in vitro" proteolytic activity of CTSL1 was found to be dependent on the H3 N-tail PTM status. H3K27me2 increased H3 cleavage by CTSL1, whereas H3K23ac reduced H3 cleavage. H3 clipping was also demonstrated in human ESCs (6).
A recent study by Duarte et al. proposed a role for H3 clipping in cellular senescence (7). Histone variant H3.3 was found to be proteolytically processed by CTSL1 upon oncogene-induced and replicative senescence in human fibroblasts and melanocytes. Ectopic expression of H3.3 and particularly its clipped proteoform was sufficient to induce senescence in fibroblasts presumably via transcriptional silencing of cell cycle regulatory genes.
Although the mechanism of regulation of histone clipping remains unclear, several studies suggested that this process might be affected by canonical histone PTMs (1,3,5,8). However, because of technical challenges in the characterization of co-existing histone modifications the relation between histone clipping and covalent histone PTMs has remained poorly defined. In the present study, we address this question by using a middle-down proteomic workflow optimized in our laboratory for efficient characterization of combinatorial histone modifications (9).
First, we demonstrate that the N-terminal tails of two core histones, H2B and H3, undergo proteolytic processing in human hepatocytes both in vitro in hepatocarcinoma cell line HepG2/C3A and in vivo in primary hepatocytes and liver tissue. We find that cell culture conditions have profound effect on this process. Histone clipping takes place in HepG2/3CA cell line cultivated as a spheroid 3D culture (i.e. when cells are at their metabolic equilibrium (10)) but not when grown in a flat 2D culture using conventional cell culture techniques (when cells are in exponential growth). By using middle-and topdown proteomic approaches optimized for histone analysis we localize four different H3 cleavage sites and identify the position of H2B clipping. Finally we provide a comprehensive analysis of the PTM status of clipped H3 proteoforms and show that clipped H3 contain distinct PTM patterns enriched in K3K36 mono-and dimethylation.

MATERIALS AND METHODS
Cell Culturing-Standard culture conditions (2D culture): the immortal human hepatocellular carcinoma cell line, HepG2/C3A, was grown in monolayer cultures in Dulbecco's modified Eagle medium (DMEM) as previously described (10).
Spheroid culture conditions (3D culture): the HepG2/C3A cell spheroids were prepared using AggreWell™ 400 plates (Stemcel Technologies cat. no. 27845) and grown in ProtoTissue™ bioreactors as previously described (10). Spheroids were collected at specific time points, washed 5 times with PBS, snap-frozen in liquid nitrogen and stored at Ϫ80°C until analysis.
Primary Human Hepatocytes-Primary adult human hepatocytes were obtained from BioreclamationIVT (Brussels, Belgium) Human Liver Tissue-Protein lysate of normal liver tissue was obtained from Acris Antibodies GmbH (Herford, Germany) Euchromatin and Heterochromatin Extraction Through Partial MNase Digestion-Nuclei were isolated from the 21 day old C3A 3D culture followed by chromatin fractionation and immunoblot analysis as described (11). Briefly, nuclei were isolated by cell lysis with Nuclear Isolation Buffer (NIB, 15 mM Tris pH ϭ 7.5, 15 mM NaCl, 60 mM KCl, 5 mM MgCl 2 , 1 mM CaCl 2 , 250 mM sucrose and "complete" protease inhibitor mixture (Roche)) supplemented with 0.3% Nonidet P-40. The resulting nuclei pellet was separated by centrifugation (600 rpm for 5 min at 4°C) and washed twice with NIB. Nuclei were resuspended in NIB to a concentration of ϳ10 7 nuclei/ml and chromatin was partially digested with MNase (Sigma 1.6 U/ml) at 37°C for 10 min. Reaction was quenched by 1 mM EGTA on ice for 10 min. The sample was centrifuged at 1000 rpm for 5 min at 4°C to generate the first supernatant (S1). The pellet was resuspended in 2 mM EDTA at pH ϭ 7.2 in the same volume as S 1 and incubated on ice for 10 min. The sample was then centrifuged at 12,000 rpm at 4°C to yield a second supernatant (S2) and a pellet (P).
Histone Purification-Histones were acid extracted as described previously (13). One hundred micrograms of bulk histone mixture was separated by reverse-phase high performance liquid chromatography (RP-HPLC) as previously described (14). Twelve fractions were collected between 28 -72 min of elution, lyophilized and stored at Ϫ80°C until further analysis.
Chromatography and Mass Spectrometry-Top-down Protein Analysis-RP-HPLC fractions containing H2B were combined and fractions containing H3 were combined. Then three g of protein sample (H2B or H3) from each pool were analyzed by our hybrid weak cation exchange/hydrophilic interaction liquid chromatography-tandem mass spectrometry (WCX/HILIC-MS/MS) method (9,14). Protein samples were separated by using an Ulti-Mate3000 (Dionex) HPLC system, equipped with a two-column setup, consisting of a reversed-phase trap column (3 cm, 100 m i.d., 360 m o.d., packed with ReproSil, Pur C18AQ 3 m; Dr. Maisch) and a WCX/HILIC analytical column (18 cm, 75 m i.d. 360 m o.d., packed with PolycatA 1500 Å, 3 m, PolyLC Inc). Solvent A and solvent B were composed as follows: solvent A 75% MeCN (Fisher Scientific), 20 mM propionic acid (Sigma) at pH 6; solvent B 25% MeCN and 0.3% formic acid (Sigma) at pH 2.5. Solvent C used for sample loading was composed of 0.1% formic acid. Proteins were eluted directly into an Orbitrap Fusion mass spectrometer (Thermo Fisher Scientific) with a 70-min linear gradient of 50 -90% buffer B at flow rate of 300 nL/min. MS was controlled by Xcalibur software (Thermo Fisher Scientific). The nanoelectrospray ion source (Thermo Fisher Scientific) was used with a spray voltage of 1.9 kV. The ion transfer tube temperature was 275°C. MS survey scans were obtained for the m/z range of 700 -1200 in the orbitrap with maximum ion injection time of 50 ms, auto gain control target 5 ϫ 10 5 , resolution 60,000 and 10 microscans per each full MS event. Four most abundant precursor ions with charge state in the range 9 to 21 from each MS scan were selected for ETD with 20 ms reaction time. MS/MS spectra were measured in the orbitrap with maximum ion injection time 500 ms, auto gain control target 2 ϫ 10 5 , mass resolution 30,000 and 8 microscans per MS/MS. Dynamic exclusion was set to a repeat count of one with a repeat duration of 30 s, exclusion duration 45 s, and exclusion m/z width 10.
Middle-down Protein Analysis-Acid-precipitated bulk histones were resuspended in 100 mM NH 4 HCO 3 . Then a 50-g aliquot of histones was subjected to endoproteinase GluC digestion (Calbiochem, 1:20 enzyme to substrate ratio 1:20) performed for 8 h at room temperature. The resulting peptide mixture was analyzed in triplicate by the WCX/HILIC-MS/MS method adapted from our previous studies (9,14) and described above in a previous section. Peptides were eluted directly into the ESI source of a LTQ-Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) using a 180-min linear gradient of 40 -80% buffer B at flow rate of 300 nL/min. The mass spectrometer was controlled by Xcalibur software (Thermo Fisher Scientific). A nanoelectrospray ion source (Proxeon, Odense, Denmark) was used with a ESI voltage of 2.2 kV. Capillary temperature was 270°C. Dynamic exclusion (DE) was disabled during data acquisition. Data acquisition was performed in the Orbitrap for both precursor ions and product ions, with a mass resolution of 60,000 (fullwidth at half-height) for MS and 30,000 for MS/MS. Precursor charge states 1ϩ, 2ϩ, and 3ϩ were excluded. Isolation width was set at 2 m/z. The six most intense ions with MS signal higher than 5000 counts were isolated for fragmentation using ETD with an activation Q value of 0.25, activation time of 90 ms with charge state dependent ETD time enabled, and Supplemental activation. Three microscans were used for each MS/MS spectrum, and the AGC target was set to 2 ϫ 10 5 . Acquisition window was set at m/z 450 -750, including charge states 5ϩ to 11ϩ.
Data Analysis-Middle-down protein analysis. Raw files were processed and searched with Mascot (version 2.3.2, Matrix Science, London, UK) using the Proteome Discoverer (version 1.4.0.288, Thermo Fischer Scientific) data analysis package. Spectra were deconvoluted by using Xtract (Thermo Fisher Scientific) with the following parameters: S/N threshold was set to 0, resolution at 400 m/z was 30,000 and monoisotopic mass only was true. The following parameters were used in the Mascot search: MS mass tolerance: 2.1 Da, to include possible errors in isotopic recognition; MS/MS mass tolerance: 0.01 Da; semi GluC specificity with no missed cleavage allowed; mono-and dimethylation (KR), trimethylation (K), acetylation (K), and phosphorylation (ST) were set as variable modifications. MS/MS spectra were searched against histone human database, consisted of histone H3.1/H3.2 and H3.3 N-terminal sequences, covering residues 1-50, downloaded from Uniprot (updated September 2011). XML results file from Mascot were exported by ticking all the query level information and the start and end from the peptide match information (plus all the default settings). XML files were imported and processed using an in-house developed software Histone Coder (9) (see supplemental methods, supplemental Table S3). Histone Coder was used with a tolerance of 30 ppm and only c/z fragment ions allowed. Only PTMs with at least one site-determining ion before and after the assigned PTM site were accepted. Peptide quantification was performed using in-house developed software IsoScale (9) (see supplemental methods, supplemental Table S4). By using IsoScale we extracted the total ion intensity from MS/MS spectra of identified peptides. We used the output from Histone Coder after eliminating the matches that contain 0 (zeros) in the site determining ion cell, excluding all PTMs with ambiguous assignments. In addition, we loaded the .csv result file(s) from Mascot with all the query information included. The total ion intensity obtained from the Mascot result files was summed for all peptide-spectrum matches (PSMs) and divided by the sum of all peptides in all modified forms for each histone to obtain the relative abundance of each PTM form. Different clipped histone proteoform were considered as different proteins. IsoScale also calculated the fragment ion relative ratio (FIRR) of isobaric cofragmented peptides and divided the total ion intensity of the given MS/MS spectrum between the two species (see supplemental Methods). The MS/MS tolerance to identify fragment ion for FIRR calculation was set to 0.05 Da. The software tools Histone Coder and isoScale are freely available (9).
The relative abundance of each individual modification was calculated as sum of relative abundances of all PTM forms containing this modification. Different histone proteoforms were considered as a different proteins and relative abundance of each individual modification was calculated for each H3 proteoform separately.
The "interplay score" was used to evaluate the likelihood of two histone marks to coexist on the same polypeptide as previously described (15). Interplay score was calculated by using the following equation: Interplay xy ϭ log2 (F xy /(F x * F y )), where xy is the coexistence frequency of the binary mark, Fx is the frequency or relative abundance of the mark x and Fy is the frequency of the mark y.
For a detailed description of the top-down protein data analysis, please see supplemental Methods.  Tables (supplemental Table S1 and Table  S2) and supplemental Methods. Supplemental tables.xml contains supplemental Table S3 and supplemental Table S4.

Histones H2B and H3 Undergo Proteolytic Processing in a Hepatocarcinoma Cell Line Grown in 3D
Culture-We first demonstrated the occurrence of histone clipping in the human hepatocarcinoma cell line HepG2/C3A (hereafter C3A). We used immunoblotting to probe whole-cell extracts (WCEs) of the C3A cell line grown in flat culture (2D) and as spheroids (3D). By using antibodies (Abs) specific for the C-terminal part of histone H3 we reproducibly observed two faster migrating H3 bands in WCE isolated from 3D culture, but not in C3A or HeLa grown in conventional 2D culture (Fig. 1A). To confirm the obtained results we analyzed histone fraction isolated from 3D culture using SDS-PAGE and liquid chromatographytandem mass spectrometry (LC-MS/MS) after limited in-gel trypsin digestion. As expected the faster migrating H3 band contained peptides corresponding to histone H3. Interest- ingly, in addition to H3 this band also contained peptides derived from the histone H2B, suggesting that this protein might also be proteolytically processed (supplemental Fig.  S1A). We confirmed the presence of clipped H2B in WCE from C3A cells by Western blot analysis using anti-H2B Ab. The additional low molecular weight H2B band was detected only in 3D culture, but not in C3A or Hela 2D cultures (Fig. 1A).
To rule out the possibility that proteolysis of histone H3 occurred during sample preparation, we carried out a control experiment. Purified intact H3 from C3A cells grown in 2D culture was labeled with cysteine-reactive Tandem Mass Tag (TMT) reagent and added to 3D culture cell pellet with lysis buffer. Subsequently the cells were lysed and the resulting WCE was analyzed by Western blot using anti-TMT Ab (supplemental Fig. S2A). Uncontrolled degradation of TMT-tagged H3 by endogenous proteases during sample preparation would lead to multiple bands in the Western blot analysis. Thus the analysis revealed that histone H3 clipping did not occur during sample preparation (supplemental Fig. S2B). This procedure was not applicable to H2B because of the absence of cysteine residues in its sequence.
Clipped H2B and H3 Histone Forms are Associated with Euand Heterochromatin-To test whether clipped H2B and H3 are associated with distinct chromatin states, we evaluated the relative abundance of clipped histones in crude euchromatin and heterochromatin subfractions (11,16,17). Briefly, chromatin was isolated from 3D culture, partially digested with MNase and subsequently fractionated as described previously (11). The presence of clipped histones was examined by immunoblotting (Fig. 1B, supplemental Fig. S1B). We detected clipped H2B and H3 forms in both soluble (S1) and insoluble (S2) MNase-treated chromatin fractions, indicating that clipped histones are localized within both eu-and heterochromatin. Faster-migrating H2B and H3 bands were also observed in the final EDTA-insoluble chromatin pellet (P) containing matrix-associated chromatin (11).
Abundance of Clipped Histones Increases During C3A Cultivation as a 3D Culture-H3 clipping is involved in transcriptional silencing of cell cycle-promoting genes upon senescence in fibroblasts (7). Previously we demonstrated that C3A cells in 3D culture grow significantly slower than in classical 2D culture (18). The doubling time for the C3A cell line grown in 2D cultures is about 3 days, whereas in 3D it increases drastically to 17 days in 21-day old spheroids and to 60 days in 42-day old spheroids (18). To test whether histone clipping can be correlated with proliferation activity of the C3A cell line, we assessed the presence of clipped histones in 3D culture at different stages of cultivation ( Fig. 2A). A signal from a faster migrating H3 band was first detected after 10 days of cultivation and increased up to day 15-24 (Fig. 2B). The similar trend was observed for histone H2B; the faster migrating H2B band was first observed after 10 days of C3A cultivation and increased in amount up to day 47 (Fig. 2B). To check whether histone clipping is associated with 3D culture or it occurs in cells growing by both cultivation techniques but gets diluted in 2D culture because of the faster cell doubling time we tested 2D culture cells exhibiting growth arrest caused by contact inhibition. C3A cells were grown in conventional 2D culture to confluence and then cultivated for 10 more days, prolonging the growth in arrested duplication state. Cells were collected at 1, 5, and 10 days post confluence and corresponding WCEs were probed with anti-H2B and anti-H3 Abs. No faster migrating histone band was detected in any of the post-confluence cultures (supplemental Fig. S3), demonstrating that histone clipping is not evenly accumulated in 2D culture, but is likely related to the functional differences be- tween cells in 2D and 3D cultures. The possible contribution of cell death on proteolytic processing of histone proteins was eliminated by using TUNEL staining (data not shown), which confirmed that there are very few apoptotic cells (Ͻ3%) in the spheroids C3A culture, indicating that histone clipping is not caused by cell death-associated protein degradation.
Histone Clipping Occurs in Human Primary Hepatocytes and Human Liver Tissue-Tumor-derived cell lines differ from both normal and tumor tissues (of the same origin) in many ways, including protein expression profiles (19). Much less is known about the epigenetic differences between those cell types, including differences in histone PTM abundance. To determine whether histone clipping is restricted to C3A cells and whether it also occurs in normal liver cells we analyzed human primary hepatocytes and liver tissue. Western blotting analysis confirmed the presence of clipped forms of both histones H3 and H2B in primary human hepatocytes (Fig. 2C). The immunoblot pattern for histone H2B was similar to those seen in hepatocarcinoma cell line, whereas for histone H3 more bands were detected in primary cells as compared with 3D culture. The similar pattern for H3 was observed in liver tissue lysate, however no faster migrating species for H2B were detected (Fig. 2C).
Identification of Histone Clipping Sites Using MS-To define the nature of the lower molecular weight histone proteoforms and determine the sequence localization of histone clipping sites we sequenced the faster migrating histone species using mass spectrometry (20). Histones were isolated from 3D culture and separated by reversed-phase liquid chromatography (RPLC) (13). The presence of clipped histones was determined by immunoblotting (supplemental Fig. S4C). Fractions containing faster migrating H2B and H3 species were analyzed using a top-down proteomics strategy, i.e. characterization of intact histones by using nanoliter flow LC interfaced to electrospray ionization tandem mass spectrometry (LC-MS/MS) (9,21). We identified both the intact and clipped histone H2B proteoforms (supplemental Fig. S5, supplemental Table S1). The clipped H2B lacked 17 N-terminal residues demonstrating that proteolytic processing of H2B occurs between residues K17 and A18 ( supplemental Fig.  S5C).
We also observed truncated H3 N-tails starting at residue L20 (H3⌬Q19). The presence of such peptides can be explained by glutamine deamidation which results in the conversion of Gln19 to Glu19 (22). This, consequently, may lead to generation of artificial truncated H3 tails during histone digestion by endoproteinase GluC. We confirmed this hypothesis by analysis of Arg-C digested H3 using selected reaction monitoring (SRM) (supplemental Fig. S7).
PTMs Interplay Within Intact and Clipped Histone H3-Next, using the data obtained by the middle-down proteomics analysis of residues 1-50 of the histone H3 tail we investigated the interplay between the different PTMs within intact and clipped histone H3. We calculated the "interplay score" of binary PTMs by comparing observed PTM pair frequencies with theoretical PTM pair frequencies (15). The "observed" co-occurrence frequencies of each co-existing PTM pairs were calculated based on the frequencies of peptides with combinatorial PTMs containing these PTM pairs. Theoretical PTM pair co-occurrence frequencies were calculated based on the single PTM relative abundances on the assumption that they are independent (14). Briefly, such interplay score considers how often two PTMs should co-exist on the same protein based on their relative abundance and compares this value with the experimental value (15). We examined which PTMs are unlikely to co-exist (negative interplay) and which are most likely co-dependent (positive interplay) for intact and clipped H3 N-terminal tails. The most abundant PTMs including K14ac, K23ac, K27me2/me3, and K36me1/me2 were examined (Fig. 5A). We found that K27me3 rarely co-exist with K36me1/me2 (except for the H3⌬G12 tail). This result conforms to previous studies suggesting antagonism between high methylation states of H3K27 and H3K36 (23,24). K14ac and K23ac were found to be mutually exclusive by interplay analysis, in agreement with a previous report (15). Interestingly, the calculated interplay scores for some of the H3 PTM pairs were different for intact and clipped H3 tails. For instance, K36me1/me2 was found with a negative interplay score with K23ac on clipped H3 tails, whereas this was not so pronounced in intact H3 tails (Fig. 5A, Fig. 6).
Collectively, our results demonstrate that clipped H3 forms contain distinct PTM patterns that are different from those in intact H3.

DISCUSSION
Histone Clipping Activity in Human Hepatocytes-A number of studies demonstrated that clipping of the N-terminal tail of histone H3 takes place in various organisms, including Tetrahymena thermophila, Saccharomyces cerevisiae, chicken, mouse, and human (1, 2, 5-7, 25, 26). Although this phenomenon is likely associated with transcriptional activation in yeast, the significance of H3 clipping in higher eukaryotes remains poorly understood (1)(2)(3)(4). Here, we demonstrate that site-specific clipping of core histones H2B and H3 takes place in human hepatocytes both in vitro and in vivo. Histone clip- ping was induced during the cultivation of human hepatocarcinoma cell line, HepG2/C3A, in 3D culture but not in 2D culture. There is growing evidence that cell culture techniques can have a profound effect on cell metabolism (27,28). In particular, cells cultured in 3D spheroids more closely mimic original tissues compared with those grown using traditional 2D cultures (29). Recently Duarte et al. proposed a role for histone H3.3 clipping in down-regulation of cell cycle-promoting genes during senescence in human fibroblasts (30). Previously we demonstrated that proliferation rate of hepatocarcinoma cells constitutively decreases during the cultivation in 3D culture and after 21 days of growth the majority of cells are either arrested in G 1 or have entered G o phase (31). We now show that the abundance of clipped H2B and H3 consistently increases in C3A cells during their cultivation in 3D culture. By using middle-down MS analysis of histones isolated from human primary hepatocytes we identified four different clipped H3 proteoforms-namely H3⌬K9, H3⌬S10, H3⌬G12, and H3⌬K23-and demonstrate that clipping preferentially occurs in histone variant H3.3 at the same regions within N-terminal tail (between lysine 9 and glycine 12 and after lysine 23) as was observed in senescent fibroblasts (30). Although the role of histone clipping in hepatocytes remains undefined, we hypothesize that it might be involved in regulation of cell-cycle progression as was previously proposed for human fibroblasts.
Histone modifications are known to play an important role in the regulation of many chromatin-associated processes. Since Allis and colleagues have proposed the "histone code hypothesis" in 2001, many findings provided strong evidence that histone modifications can act in a cooperative manner (32)(33)(34). Different combinations of histone PTMs lead to different functional outcomes, providing a strong and flexible platform for the regulation of key cellular processes. Previously, Duncan et al. found that different H3 PTMs can modulate proteolytic activity of the Cathepsin L enzyme, which was found to be responsible for histone H3 clipping during mouse ESCs differentiation, in vitro (5). The combinatorial PTM status of cellular clipped histones, however, remains poorly defined. By using an optimized middle-down MS platform that is highly efficient for the characterization of histone PTM isoforms we carried out for the first time a comprehensive comparison of histone PTMs between intact and clipped H3 N-terminal tails isolated from the C3A cell line grown in 3D culture. First, by using MS and MS/MS we identified two different clipped proteoforms of H3, namely H3⌬S10 and H3⌬G12. Furthermore we show that the relative abundance of distinct H3 PTMs, such as K14me1, K27me2, K27me3, K36me1, and K36me2, is significantly different between clipped and intact H3 N-terminal tails. H3⌬S10 and H3⌬G12 forms were found to contain extremely low levels of K14me1, suggesting a strong negative correlation between the H3K14 methylation and clipping of H3. This finding supports the previously reported idea that H3 clipping can be regulated by other H3 PTMs (5). However which of these events occurs first is not known. Thus, there are several possible explanations of the relationship between H3 clipping and H3 lysine 14 methylation. First, H3K14me1, as a mark near the H3 clipping site, can decrease the efficiency of the cleavage site recognition by the "clipping protease," thus, preventing H3 processing. Alternatively, H3K14me1 can be specifically demethylated in clipped H3. It is also possible that these events do not have a direct link and H3 clipping just occurs in chromatin regions which have low level of H3K14 methylation. These questions need to be addressed in future studies.
We found that clipped H3 proteoforms H3⌬S10 and H3⌬G12 are highly enriched in both H3K36me1 and H3K36me2 as compared with intact H3. Histone H3 K36 methylation is known as an active mark which has several well established roles in different nuclear processes, including transcription, DNA repair and regulation of Polycomb repressive complex 2 (PRC2) (23,(35)(36)(37).
Numerous studies demonstrated that H3K36 methylation is associated with transcribed chromatin in various organisms (38 -40). In human TKO2 myeloma cell line H3K36me2 was enriched near the transcription start sites (TSSs) and the level of H3K36me2 correlated with gene activity (41). In yeast cotranscriptional dimethylation of H3K36, catalyzed by RNA polymerase II (pol II)-associated methyltransferase Set2, recruits Rpd3S deacetylase complex and consequently prevents aberrant transcriptional initiation within coding sequences by maintaining the hypo-acetylated state of transcribed chromatin (42,43). Our data reveal that clipped H3 proteoforms are over 4 times enriched in K36me2 as compared with intact H3. Furthermore, H3 tails methylated at K36 exhibit a low acetylation level primarily because of the negative crosstalk between K36me1/me2 and H3K23ac. Taken together, these findings demonstrate that clipped H3 proteoforms have a PTMs pattern similar to those known to occur in H3 located in transcribed chromatin regions. Although the role of histone clipping in transcriptional regulation still needs to be explored previous studies suggested that this irreversible modification can act both as a transcriptional activator and repressor (1,30,44,45). Here we provide firs time comprehensive characterization of H3 clipping sites in human cells. We mapped four different cleavage sites within the H3 N-terminal tail: after lysine 9 (H3⌬K9), serine 10 (H3⌬S10), glycine 12 (H3⌬G12), and lysine 23 (H3⌬23). This result can be used for generation of antibodies specific to distinct clipped H3 proteoforms and further characterization of proteolytic processing of histone H3 and its genomic distribution. Although proteolytic processing of histone H3 has been demonstrated in several species there was no prior evidence of clipping of H2B, except for nickel-induced H2B truncation (46). We here provide the first evidence of H2B clipping as a part of normal cellular functionality. Using top-down MS we identified two distinct intact H2B variants, H2B type 1-C/E/F/ G/I and H2B type 2-E. The clipped proteoform lacking 17 N-terminal residues was identified only for H2B type 1-C/E/ F/G/I, suggesting that H2B clipping might be variant-specific. Positive and negative interplay is represented in red and blue, respectively. PTM pairs with negative (positive) interplay values co-occur less (more) frequently than if they were independent. A negative interplay score corresponds therefore to mutual exclusion indicating competition between the two marks whereas positive interplay reveals enhanced, correlated crosstalk. PTMs which appear to be mutually exclusive are indicated in bold. B, The total relative abundance of acetylation marks in peptides carrying K36me1 or K36me2 and the total relative abundance of acetylation marks calculated based on all identified peptides. Histone H3 containing K36me1 or K36me2 demonstrate hypoacetylated status. In summary, our study provides the first observation of the proteolytic processing of histones H2B and H3 in human hepatocytes and the first detailed PTM characterization of clipped H3, including co-existing PTMs. Our data highlight the relationship between H3 clipping and canonical H3 PTMs and provide a basis for further investigation of this phenomenon.