- Open Access
Quantitative proteomic analysis of histone modifications in decitabine sensitive and resistant leukemia cell lines
Clinical Proteomicsvolume 13, Article number: 14 (2016)
The refractory nature of many cancers remains the main health challenge over the past century. The epigenetic drug, decitabine (DAC), represents one of the most promising therapeutic agents in cancers particularly in acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS). However, its ambiguous anti-tumor mechanism and the unpredictable drug-resistant nature in some population compromise its application in cancer therapy. In crosstalk with DNA methylation, histone post-translational modifications (PTMs) are the key players in modulating the downstream epigenetic status of tumor suppressor genes. This study targets the role of decitabine in epigenetic regulation in leukemia therapy and searches responsive predictors and therapeutic targets for pretreatment evaluation and drug development.
A simple, fast, and robust proteomic strategy identified 15 novel PTMs and 60 PTM combinations in two leukemia cell lines (MDS-L and TF-1). Histone modification profiles have been generated and compared between DAC sensitive and resistant groups (n = 3) in response to DAC treatment. Among these histone PTMs, five of which were found differentially upon DAC treatment in drug sensitive and resistant cells: H3.3K36me3, H4K8acK12acK16ac in MDS-L cells; and H3.1K27me1, H3.1K36me1, H3.1K27me1K36me1 in TF-1 cells. They may serve as biomarkers in predicting leukemia and drug responsiveness. In addition, we also explored PTM differences in two cell lines which were developed from early and advanced stages of AML. Three PTMs (H3.1K27me3, H3.1K27me2K36me2 and H3.3K27me2K36me2) are highly abundant in TF-1 cells (advanced AML cell line), suggesting their relevance to leukemogenesis. Our method allowed deep analysis of histone proteins and elucidation of a large number of histone PTMs with high precision and sensitivity.
DAC-induced DNA hypomethylation has wide impact on chromatin modifications. This study represents first effort to investigate the undefined epigenetic mechanism of decitabine in leukemia therapy. The identification of 15 novel PTMs and the discovery of several marks have relevance to epigenetic directed therapies.
Epigenetics is an active field of cancer research. DNA methylation and histone modifications are two major contributors to epigenetics. As a crucial epigenetic mark on the genome, DNA methylation modulates many important cellular processes including embryonic development, transcription, mammalian X chromosome inactivation, genomic imprinting and chromosome stability [1, 2]. About 60–80 % of CpGs in the human genome undergo DNA methylation . There are two classes of DNA methyltransferases (DNMTs) that catalyze the transfer of the methyl group onto DNA: (1) DNMT1 maintains the global methylation pattern; and (2) DNMT3A and DNMT3B perform de novo methylation during embryonic development [3–7]. DNA methylation is generally considered a repressive mark and is often associated with gene silencing [1, 8]. The balance between DNA methylation and demethylation is required to be precisely maintained and dysregulation of the balance may lead to human diseases notably cancer [3, 9, 10]. Global DNA hypomethylation as well as DNA hypermethylation in the promoter regions of tumor suppressor genes are common hallmarks of cancer cells [10–12].
Histone modifications, on the other hand, provide another layer in modulating DNA replication and gene transcription. Histones are relatively small proteins with high content of basic residues, lysine and arginine, making these proteins strongly positively charged which contributes to their tight interactions with the negatively charged DNA backbone . There are two types of histones: (1) Linker histone H1; and (2) Core histones (H2A, H2B, H3, and H4) [14, 15]. A protein octamer formed by two of each copy of core histones constitutes the basic unit of eukaryotic chromatin [13, 14]. Both linker histones and core histones undergo a large number of chemical modifications including methylation, acetylation, phosphorylation, biotinylation, citrullination, ADP-ribosylation, and ubiquitination . Histone modifying enzymes are responsible for addition or removal of these different types of chemical modifications [13, 17, 18]. A wide range of dynamic histone marks (histone code), their modifying enzymes (code writers or erasers), and their downstream effector molecules (code readers) are key players in regulating eukaryotic chromatin structure and functions [19–22]. Multiple regulatory layers can be achieved by changing the levels, types, and positions of PTMs on these proteins. As with DNA methylation, a wrong histone code (i.e. an aberrant histone modification pattern) written or erased by histone modifying enzymes may result in disease [23–26].
As DNA and histones are physically intertwined, changes in chromatin states often require synergistic actions affecting both DNA methylation and histone modifications involving proteins including DNMTs, histone modifying enzymes, and chromatin binding proteins [27–29]. Little is known, however, about the crosstalk between these two types of modifications in many diseases. Some preliminary studies have pointed to a connection between these two types of modifications and their relevance to human cancers [8, 11, 28]. DNA methylation may affect histone modification patterns and altered histone modifications also impact the DNA methylome and vice versa [30–37]. Therefore, all of these marks and proteins mentioned above provide a rich source to mine for cancer biomarkers and therapeutic targets.
The recent advancements in cancer epigenetics has led to the development of new epigenetic drugs, most of which are DNA methyltransferase inhibitors (DMNTi) or histone modifying enzyme inhibitors (e.g. HATi, HDACi, HMTi, and HDMi) [38–42]. Decitabine (DAC), a DNA methylation inhibitor, has potential as a therapy for myelodysplastic syndromes (MDS) and certain types of leukemia [43–48]. DAC depletes DNMT1 and reverses aberrant epigenetic repression of tumor suppressor genes via an unknown mechanism. Low-dose administration of DAC is less toxic and improves efficacy in MDS and other cancers [49–52]. On the other hand, ~50 % of patients are non-responders and most of patients eventually develop resistance to the drug . Patients who are resistant to DAC have limited alternative options and have high mortality. The ambiguous anti-tumor mechanism and unpredictable basis for drug-resistance represent a current challenge. Therefore, knowledge of DAC anti-tumor effects and mechanisms of resistance has essential relevance to cancer therapy.
Liquid chromatography tandem mass spectrometry (LC–MS/MS) technique has become the most popular and powerful tool for large-scale protein identification, PTM analysis, protein–protein interactions, and etc. [53, 54]. The quantitative feature gives it extra credit in comparing thousands of components across multiple samples in different cell states [55–57]. However, analysis of histones and their posttranslational modifications are challenging owing the high complexity and combinatory manners of the numerous chemical modifications, which often involves the multi-step histone purifications such as HPLC and gel electrophoresis [58–60]. As a result a large number of cells or tissues are required to obtain enough histones for mass spectrometry analysis, complicating analysis of clinical samples available in limited quantities. Since histone proteins have excess basic amino acid residues (lysine and arginine), trypsin digestion will generate undersized, hydrophilic peptides which are unfavorable when conventional proteomics strategy is adopted. Chemical derivation of histones, such as propionylation , is a popular protocol which greatly facilitates LC–MS/MS analysis by introducing an artificial modification (propionyl group) onto lysine residues generating longer, more hydrophobic tryptic peptides. However, operational variations during sample preparation are introduced since propionylation of histones is usually incomplete and varies from batch to batch . For comparison of histone PTMs in different biological states, chemical or metabolic labeling techniques often provide more accurate quantification than label-free methods which require more measurements and repeats to achieve the power of statistical tests [63–65]. Although middle-down analysis of large-size peptides (3–9 kDa) or top-down analysis of intact histones gives more comprehensive assessment of the histone code, the use of these approaches is largely limited since they require more separation steps and large amounts of histones . Targeted proteomics allows absolute quantification of histone PTMs by spiking samples with synthetic, isotope labeled peptides . Nonetheless, this strategy is designed upon biomarker validation not for discovery purpose. To dissect the role of DAC in epigenetic regulation in leukemia therapy, we have implemented a simple, fast, highly sensitive proteomic method for mapping and quantifying histone modifications and their combination patterns in leukemia cells based on high-resolution mass spectrometry. Prior to enzyme digestion and LC–MS/MS analysis, this method relies on a simple histone enrichment step without HPLC or gel separation, thus greatly simplifying the overall workflow. Propionylation in combination with stable isotope labeled histones as internal standard further improves the PTM identification and quantification. Using this strategy, we have identified 61 individual histone marks and quantified the relative levels of 60 PTM combinations in two leukemia cell lines. We also investigated effects of treatment with thus linking DNA methylation and chromatin modifications.
Cell culture and experimental design
Two cell lines, TF-1 (derived from human erythroleukemia, purchased from ATCC) and MDS-L (derived from patient MDS) were chosen for this study. TF-1 cells were maintained in RPMI-1640 containing 10 % fetal bovine serum and 10 ng/mL interleukin 3. MDS-L cells were maintained in RPMI-1640 containing 20 % fetal bovine serum and 10 ng/mL interleukin 3. An additional cell line U937 was maintained in RPMI-1640 medium containing 10 % fetal bovine serum. To develop resistance to DAC, parental (DAC-sensitive) TF-1 and MDS-L cell lines were challenged with decitabine over a year. Both DAC-sensitive and DAC-resistant cell lines were then cultured in DAC-free medium and analyzed with or without treatment with DAC (1-3 µM) for 72 h in the same medium. To obtain the global internal standard, the three parental cell lines (MDS-L, TF-1, U937) were cultured in SILAC medium (R1780 SIGMA, RPMI-1640 Medium) supplemented with heavy amino acids-l-ARGININE: HCL (13C6, 99 %, CLM-2265-H-0.5) and l-LYSINE: 2HCL (U-13C6, 99 %; CLM-2247-H-0.5) purchased from Cambridge Isotope Laboratories. After SILAC labeling, they were mixed in equal numbers prior to nuclear isolation and histone acid extraction. Meanwhile, regularly cultured (unlabeled) cell lines also underwent nuclear isolation and acid extraction separately (Fig. 1).
Nuclear isolation and bulk histone acid extraction
About 1 × 106 cells were collected by centrifugation and washed twice with 10 mL ice-cold PBS (Phosphate-buffered saline) supplemented with 5 mM Sodium Butyrate (Sigma-Aldrich). Cell pellets were suspended in 1 mL ice-cold TEB Buffer (Triton Extraction Buffer: PBS containing 0.5 % Triton X 100 (v/v), 2 mM PMSF) supplemented with 5 mM Sodium Butyrate, 1× Protease inhibitor (Roche), 1× PhoStop (Roche), and incubated in a rotator at 4 °C for 10 min. Intact nuclei released from cells were collected by spinning at 3000g for 10 min at 4 °C. Nuclei were washed with 1 mL TEB buffer and spun down again. For histone acid extraction, nuclear pellets were resuspended in 1 mL of 0.4 N H2SO4 with gentle spin in a rotator overnight at 4 °C. Nuclear debris was removed by centrifugation at 16,000g for 10 min at 4 °C in a cooled table-top centrifuge. The supernatant containing bulk histones was transferred into a new Eppendorf tube for precipitation by addition of TCA (Trichloroacetic acid, Sigma-Aldrich) drop by drop to the histone supernatant (final concentration of TCA = 33 %). Histones were collected by spinning at 16,000g, for 10 min at 4 °C. Pellets containing histones were washed with ice-cold acetone twice. Once acetone was removed, pellets were air dried for 20 min at room temperature, dissolved in 400 µL H2O, and stored at −20 °C. Both linker histone and core histones could be effectively enriched after this simple step. Typically, more than 100 µg proteins can be extracted from 1 × 106 cells and up to 80 % of the proteins are histones. Purity of crude histones after acid extraction was checked by means of SDS-PAGE (Sodium dodecyl sulfate polyacrylamide gel electrophoresis) and gel staining with Coomassie Brilliant Blue.
Chemical derivatization, trypsinization, and LC–MS/MS
Unlabeled bulk histones extracted from TF-1 or MDS-L (sensitive or resistant cells treated with or without DAC) were equally mixed with SILAC labeled crude histones extracted from a three cell line mixture which serves as internal standard for comparison across samples collected under different conditions (Fig. 1). The light and heavy protein mixture was chemically propionylated twice before and after trypsin digestion as previously described [61, 68]. Propionylation significantly increases the hydrophobicity and retention times of tryptic peptides on a reversed-phase column, which improves the number of PTMs identified and quantified by LC–MS/MS. Samples were then mixed in 100 mM ammonium bicarbonate buffer and digested with sequencing-grade trypsin (Promega, Madison, WI) at a ratio of 1:20 (enzyme:substrate) at 37 °C for 6 h. The digested mixture was vacuum-dried, reconstituted in 0.1 % formic acid and filtered by 10 kDa cut-off centrifugal filter unit (Millipore Ultracel YM-10). About ~1 μg of protein digest was loaded onto a C18 trap column (Waters, 180 μm ID × 20 mL) and resolved on a 25 cm long, capillary column (75 μm ID × 360 μm OD, New Objective) packed with 5 μm, 200 Å C18 silica-bonded material (Magic C18 AQ, New Objective). Peptides were eluted using a linear gradient as follows: run 100 % solvent A (0.1 % formic acid) over 5 min; run a 0–40 % gradient against solvent B (0.1 % formic acid in acetonitrile) over 90 min; and finally run 10 min, 95 % solvent B at a flow rate of 300 nL/min. The eluted peptides were analyzed using a Thermo Fisher Q Exactive Hybrid Quadrupole-Orbitrap Mass Spectrometer under data dependent acquisition mode. After a full scan (m/z 300–1800), tandem spectra were collected by selecting the 10 most intensive peaks with normalized collision energy of 35 %. Resolution on MS1 and MS2 was 70,000 and 17,500 respectively. A dynamic exclusion of 35 s and internal mass lock (445.12002 m/z) were also executed to obtain best performance.
A quantitative proteomics software package MaxQuant (version 126.96.36.199) was used for raw data processing, peptide identification, and quantification . Collected spectra were searched against Uniprot Human proteome FASTA database (released in April 2015) containing 90,411 sequences. Mass tolerance for precursor ion was 10 and 20 ppm for the fragment ions. N-terminal propionylation was considered as fixed modification while lysine acetylation, lysine and arginine methylation, and lysine propionylation were searched as variable modifications. Up to five missed cleavages were allowed during digestion as trypsin doesn’t cleave propionylated lysine. Data were normalized so that the median of the logarithms of all peptide ratios in each LC–MS run is zero . We accepted 1 % false discovery rate at both peptide and protein level estimated from decoyed sequences. All spectra identified with PTMs were manually checked and validated. The peptide ratios of light/heavy were evaluated in biological triplicates. The true ratios of peptides and proteins in any two different conditions can be directly converted from their light/heavy ratios. Unpaired, two-tailed Student’s t test was analyzed in Microsoft Excel.
Results and discussion
Research strategy and method development
Owing to the challenges and difficulties in analyzing histone PTMs (discussed in “Background” section), we have developed a simple, fast, and robust quantitative LC–MS/MS proteomic method. Compared to the existing methods, our entire experimental platform has several advantageous features for confident PTM identification and quantification in human cells: (1) The simple enrichment procedure with acid extraction effectively isolated bulk histones containing both linker histone and core histones with some minor non-histone contaminants, which greatly reduces the workload, and improves sensitivity, throughput and reproducibility. The whole enrichment procedure took less than 1 day before trypsin digestion and LC–MS/MS analysis. (2) In comparison to label-free quantification, the use of internal standard derived from three SILAC-labeled cell lines reduces sample sizes and gives a better PTM coverage owing to PTM pattern varies from cell line to cell line. In comparison to multiplex SILAC labeling, our protocol is more cost-effective and also allows longitudinal study of histone modifications in clinical samples. (3) The chemical derivatization took advantage of the fact that propionylation blocks trypsin cleavage at lysine and neutralizes positive charges on lysine residues , thus generates adequately sized, more hydrophobic peptides suitable for reverse-phase mass spectrometry analysis. The variations caused by incomplete propionylation during sample preparation are greatly minimized by using the internal standard which enables reliable and consistent comparison of protein and PTM abundance in different samples . (4) The use of Q Exactive mass spectrometers that offer high-resolution option on both MS1 and MS2, allows for high confidence identification of peptides and PTMs. All together, these efforts greatly improve the number of PTMs identified and quantified by LC–MS/MS [63, 68, 70, 71]. Our optimized method allowed profiling of histones and their PTMs from as few as 105 leukemia cells.
In total, 22 histones or histone variants were identified with at least one unique peptide (Additional file 1: Table S1) and 108 modified peptide species (Additional file 2: Table S3) consisting of 61 distinct histone marks at 39 sites in both linker and core histones (Additional file 3: Figure S1, Additional file 4: Figure S2). Importantly, conserved histone PTMs were detected and quantified in this study. Among these marks, 15 are novel PTMs (Additional file 3: Figure S1) according to the UniProt Knowledgebase. They are: H1.0 Lys-40 mono-methylation, H1.2 Lys-52 (or H1.3 Lys-53, H1.4 Lys-52) mono-methylation, H1.5 Lys-55 mono-methylation, H1.5 Lys-78 mono-methylation, H1.4 Lys-119 (or H1.5 Lys-122) mono-methylation, H1.4 Lys-121 (or H1.5 Lys-124) mono-methylation, H1.2/H1.5 Lys-172 mono-methylation, H2A type 1-D Lys-43 mono-methylation, H2A.1/H2A type 1-C/H2A type 1-D Lys-126 mono-methylation, H2B type 2-F Lys-21 mono-methylation, H2B type F-S Lys-25 mono-methylation, H2B type F-S Lys-29 mono-methylation, H4 Lys-9 mono-methylation, H4 Arg-56 mono-methylation, and H4 Lys-78 mono-methylation. In addition, 40 individual PTMs or 60 PTM combinations were quantified in leukemia sensitive and resistant cell lines (Additional file 5: Table S2). After normalization, most proteins and PTMs quantified had Light/Heavy ratios close to one with small statistical variations (Fig. 3; Additional file 1: Table S1, Additional file 5: Table S2) suggesting reliability of the method.
PTMs related to decitabine resistance
As our primary effort, we first targeted the DAC resistance in MDS-L and TF-1 cells which represent early and late stage of the disease respectively. To identify PTMs that are related to DAC resistance, we investigated histone modification profiles in DAC-sensitive and DAC-resistant cells upon drug treatment.
In MDS-L cells, two peptide species (Additional file 5: Table S2, Figs. 2, 3b) showed significantly difference among sensitive and resistant groups in response to DAC treatment (p < 0.05): 27-KSAPSTGGVKme3KPHR-40 (H3.3K36me3, Additional file 6: Figure S3. Note: N-terminal methionine is not counted in the following text) and 4-GKGGKacGLGKacGGAKacR-17 (H4K8acK12acK16ac, Additional file 6: Figure S3). H3.3K36me3 was moderately induced in MDS-L resistant cells after DAC treatment (final ratio = 1.39, p < 0.01) as compared to sensitive cells (final ratio = 1.09). Whereas tri-acetylation on the H4 N-terminal tail (H4K8acK12acK16ac) was reduced after DAC treatment in MDS-L sensitive cells (final ratio: 0.68, p < 0.05) but this modification remained mostly unchanged in resistant cells upon drug exposure (final ratio: 0.93). Alteration of histone lysine methylation and their enzymes are often associated with leukemogenesis or tumor progression and drugs targeting disordered patterns of histone methylation are being developed . On the other hand, deregulation of histone acetylation and their modifying enzymes (HATs and HDACs) are common observations in AML and MDS . H4K8acK12acK16ac may be relevant to the development of inhibitors of HAT and HDAC as in the case of Vorinostat and Romidepsin .
In TF-1 cells, three mono-methylation states on H3 Lys-27, Lys-36, and Lys-37 changed differently among sensitive and resistant groups in response to DAC treatment (Figs. 2, 3b). 27-Kme1SAPATGGVKKPHR-40 (H3.1K27me1, Additional file 6: Figure S3) was reduced by about 50 % (final ratio = 0.51, p < 0.05) in resistant cells while this modification maintained stable in TF-1 sensitive cells (final ratio = 1.08). The other two modification states: 27-KSAPATGGVKme1KPHR-40 (H3.1K36me1, Additional file 6: Figure S3) and 27-Kme1SAPATGGVKme1KPHR-40 (H3.1K27me1K36me1, Additional file 6: Figure S3) were moderately decreased in resistant cells after drug treatment. These methylation states were also stable in sensitive cells in response to the drug treatment. EZH2, one of the H3K27 methyltransferases and the principal component of Polycomb Repressive Complex 2 (PRC2), is a key player (code writer) in regulating the methylation states of Lys-27 in histone H3 . EZH2 also controls DNA methylation by interacting with DNMTs . Given its key function in maintaining gene transcription, mutation or deregulation of EZH2 is associated with cancer development . Development of inhibitors of EZH2 represents a therapeutic strategy for cancer . Our study has provided additional evidence that H3K27 methylation may play a role in DAC treatment effect. However, no common individual marks related to drug resistance were found in both cell lines. This may be explained by cell line difference or disease stages.
Highly abundant PTMs in TF-1 cells
As our secondary effort, we investigated and compared the PTM patterns in these two cell lines since alterations of epigenetic marks may be related to disease stages. The TF-1 cell line was first established from a patient with erythroleukemia, a rare form of acute myeloid leukemia (AML) . The MDS-L cell line was derived from myelodysplastic cells . MDS could transform into AML thus is generally regarded as a pre-cancer stage of AML . Therefore, studying these clinically closely related cell lines could provide information about leukemia development and progression.
By comparing steady-state level of histone PTMs in the sensitive cell lines, we found that di- and tri-methylation states of H3K27 were significantly more abundant in TF-1 cells compared to MDS-L cells (Fig. 2). Three species were identified in total consisting of 27-Kme3SAPATGGVK-36 (H3.1K27me3, Additional file 6: Figure S3), 27-Kme2SAPATGGVKme2KPHR-40 (H3.1K27me2K36me2, Additional file 6: Figure S3), and 27-Kme2SAPSTGGVKme2KPHR-40 (H3.3K27me2K36me2, Additional file 6: Figure S3) that were substantially more abundant in TF-1 or missing in MDS-L (Additional file 5: Table S2). H3.1K27me3 and H3.3K27me2K36me2 were identified in all samples from the TF-1 cell line but were not detectable in all MDS-L cell line samples. H3.1K27me2K36me2 were found in all TF-1 samples but were only detected in two MDS-L samples with very low Light/Heavy ratios (0.068 on average). Of note the dimethylation state of H3K27 (27-Kme2SAPATGGVK-36) was also significantly reduced in MDS-L cells compared to TF1 cells (Light/Heavy ratio: 0.23 vs 1.61, p < 0.01).
The role of H3K27 methylation has been extensively studied and a precise balance of this modification is a key to maintaining normal cell growth . There are several histone methyltransferases that impact H3K27 methylation including EZH1, EZH2, and WHSC1LI . Moreover, there are also three histone demethylases that reverse the methylation status of H3K27 . Mutations of these chromatin modifiers are often linked to cancer . H3K27 methylation is generally considered as a repressive mark and alteration of H3K27 methylation has been found in several cancer types including leukemia . Histone modifications are currently being explored as potential biomarkers for disease progression and prognosis . The induction of H3K27me2/me3 and H3K36me2 in TF-1 cells suggests a potential for these epigenetic marks to serve as biomarkers that differentiate the early (MDS) and late stage (AML) of leukemia. However, if the alteration of H3K27 methylation contributes to myeloid leukemogenesis remains to be fully elucidated.
Interplay between DNA methylation and histone modifications
Prior studies have pointed to the relatedness of DNA methylation and histone modifications in impacting tumor pathogenesis [8, 11, 28, 30–35]. Three major DNMTs (DNMT1, DNMT3A, and DNMT3B) are involved in DNA methylation in cells. All three DNMTs are post-translationally modified with a variety of PTMs including acetylation, methylation, phosphorylation, SUMOylation, and ubiquitination . Reversible covalent modifications of DNMTs may affect their stability, enzyme activity, DNA binding, and interactions with other partners [85–87]. For example, SUMOylation of DNMT3A affects protein binding to HDACs and changes DNA methylation profiles [88, 89]. DNMT1 is the primary target trapped and depleted by DAC. Among its listed known interacting proteins are chromatin modifiers including HDAC1/2 (HDACs), KDM1A (HDM), SUV39H1 (HMT), EHMT2 (HMT), EZH2 (Polycomb group protein), EED (Polycomb group protein), SETD7 (HMT), KAT5 (HAT), and etc. .
The role of DNMTs in linking these two types of modifications is of particular interest. In this study, we have demonstrated that DAC-induced DNA hypomethylation has significant impact on chromatin modifications. We have revealed five particular histone marks (H3K27me1, H3K36me1, H3K27me1K36me1, H3.3K36me3, and H4K8acK12acK16ac) that were differentially expressed in MDS-L and TF-1 sensitive and resistant cells in response to DAC treatment, pointing to a potential role for these modifications in drug resistance mechanisms. These findings suggest a DNMT-dependent pathway through which DAC inhibits DNMTs and re-activates downstream tumor suppressor genes via histone modifying enzymes or other unknown factors.
We have developed a fast proteomic method for robust quantitative analysis of histone PTMs. (1) The proteomic strategy described here enhances data quality and acquisition sensitivity resulting in comprehensive analysis of histone modifications in relatively small number of cells; (2) The systematic analysis of epigenetic profiles in drug sensitive and resistant cell lines have identified 61 PTMs, of which 15 are novel histone modifications; (3) The identification of H3.1K27me3 and H3.1/H3.3K27me2K36me2 as signature in TF-1 cells suggests their potential role in leukemogenesis; (4) The discovery of five additional histone marks that were differentially impacted by DAC in the sensitive and resistant cells suggests their potential relevance to the development of drug resistance. Our results also suggest a conserved, DNMT-dependent pathway in DAC-mediated leukemia treatment. More investigations, however, are needed to elucidate how DAC exerts its anti-tumor effect and how tumor cells develop a resistant strategy to escape the DAC-mediated anti-tumor therapy.
propionylation and monomethylation
DNA (cytosine-5)-methyltransferase 1
DNA (cytosine-5)-methyltransferase 3A
DNA (cytosine-5)-methyltransferase 3 beta
histone acetyltransferase inhibitor
histone deacetylase inhibitor
histone methyltransferase inhibitor
histone demethylase inhibitor
a first stage of mass spectrometry
tandem mass spectrometry or MS/MS
stable isotope labeling by amino acids in cell culture
Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005;6:597–610.
Jones PA, Takai D. The role of DNA methylation in mammalian epigenetics. Science. 2001;293:1068–70.
Wu H, Zhang Y. Reversing DNA methylation: mechanisms, genomics, and biological functions. Cell. 2014;156:45–68.
Robertson KD. DNA methylation, methyltransferases, and cancer. Oncogene. 2001;20:3139–55.
Bestor TH. The DNA methyltransferases of mammals. Hum Mol Genet. 2000;9:2395–402.
Siedlecki P, Zielenkiewicz P. Mammalian DNA methyltransferases. Acta Biochim Pol. 2006;53:245–56.
Cheng X, Blumenthal RM. Mammalian DNA methyltransferases: a structural perspective. Structure. 2008;16:341–50.
Baylin SB. DNA methylation and gene silencing in cancer. Nat Clin Pract Oncol. 2005;2(Suppl 1):S4–11.
Song CX, He C. Balance of DNA methylation and demethylation in cancer development. Genome Biol. 2012;13:173.
Ehrlich M. DNA methylation in cancer: too much, but also too little. Oncogene. 2002;21:5400–13.
Tsai HC, Baylin SB. Cancer epigenetics: linking basic biology to clinical medicine. Cell Res. 2011;21:502–17.
Herman JG, Baylin SB. Gene silencing in cancer in association with promoter hypermethylation. N Engl J Med. 2003;349:2042–54.
Kouzarides T. Chromatin modifications and their function. Cell. 2007;128:693–705.
Campos EI, Reinberg D. Histones: annotating chromatin. Annu Rev Genet. 2009;43:559–99.
Peterson CL, Laniel MA. Histones and histone modifications. Curr Biol. 2004;14:R546–51.
Sidoli S, Cheng L, Jensen ON. Proteomics in chromatin biology and epigenetics: elucidation of post-translational modifications of histone proteins by mass spectrometry. J Proteomics. 2012;75:3419–33.
Allis CD, Berger SL, Cote J, Dent S, Jenuwien T, Kouzarides T, Pillus L, Reinberg D, Shi Y, Shiekhattar R, et al. New nomenclature for chromatin-modifying enzymes. Cell. 2007;131:633–6.
Zhang K, Dent SY. Histone modifying enzymes and cancer: going beyond histones. J Cell Biochem. 2005;96:1137–48.
Strahl BD, Allis CD. The language of covalent histone modifications. Nature. 2000;403:41–5.
Jenuwein T, Allis CD. Translating the histone code. Science. 2001;293:1074–80.
Berger SL. The complex language of chromatin regulation during transcription. Nature. 2007;447:407–12.
Krebs JE. Moving marks: dynamic histone modifications in yeast. Mol BioSyst. 2007;3:590–7.
Chi P, Allis CD, Wang GG. Covalent histone modifications–miswritten, misinterpreted and mis-erased in human cancers. Nat Rev Cancer. 2010;10:457–69.
Portela A, Esteller M. Epigenetic modifications and human disease. Nat Biotechnol. 2010;28:1057–68.
Greer EL, Shi Y. Histone methylation: a dynamic mark in health, disease and inheritance. Nat Rev Genet. 2012;13:343–57.
Muntean AG, Hess JL. Epigenetic dysregulation in cancer. Am J Pathol. 2009;175:1353–61.
Cameron EE, Bachman KE, Myohanen S, Herman JG, Baylin SB. Synergy of demethylation and histone deacetylase inhibition in the re-expression of genes silenced in cancer. Nat Genet. 1999;21:103–7.
Esteller M. Cancer epigenomics: DNA methylomes and histone-modification maps. Nat Rev Genet. 2007;8:286–98.
Cedar H, Bergman Y. Linking DNA methylation and histone modification: patterns and paradigms. Nat Rev Genet. 2009;10:295–304.
Weisenberger DJ. Characterizing DNA methylation alterations from The Cancer Genome Atlas. J Clin Invest. 2014;124:17–23.
Clements EG, Mohammad HP, Leadem BR, Easwaran H, Cai Y, Van Neste L, Baylin SB. DNMT1 modulates gene expression without its catalytic activity partially through its interactions with histone-modifying enzymes. Nucleic Acids Res. 2012;40:4334–46.
Fuks F, Burgers WA, Brehm A, Hughes-Davies L, Kouzarides T. DNA methyltransferase Dnmt1 associates with histone deacetylase activity. Nat Genet. 2000;24:88–91.
Robertson KD, Ait-Si-Ali S, Yokochi T, Wade PA, Jones PL, Wolffe AP. DNMT1 forms a complex with Rb, E2F1 and HDAC1 and represses transcription from E2F-responsive promoters. Nat Genet. 2000;25:338–42.
Rountree MR, Bachman KE, Baylin SB. DNMT1 binds HDAC2 and a new co-repressor, DMAP1, to form a complex at replication foci. Nat Genet. 2000;25:269–77.
Peng L, Yuan Z, Ling H, Fukasawa K, Robertson K, Olashaw N, Koomen J, Chen J, Lane WS, Seto E. SIRT1 deacetylates the DNA methyltransferase 1 (DNMT1) protein and alters its activities. Mol Cell Biol. 2011;31:4720–34.
Kondo Y. Epigenetic cross-talk between DNA methylation and histone modifications in human cancers. Yonsei Med J. 2009;50:455–63.
Vaissiere T, Sawan C, Herceg Z. Epigenetic interplay between histone modifications and DNA methylation in gene silencing. Mutat Res. 2008;659:40–8.
Rodriguez-Paredes M, Esteller M. Cancer epigenetics reaches mainstream oncology. Nat Med. 2011;17:330–9.
Szyf M. Epigenetics, DNA methylation, and chromatin modifying drugs. Annu Rev Pharmacol Toxicol. 2009;49:243–63.
Di Costanzo A, Del Gaudio N, Migliaccio A, Altucci L. Epigenetic drugs against cancer: an evolving landscape. Arch Toxicol. 2014;88:1651–68.
Simo-Riudalbas L, Esteller M. Targeting the histone orthography of cancer: drugs for writers, erasers and readers. Br J Pharmacol. 2015;172:2716–32.
Azad N, Zahnow CA, Rudin CM, Baylin SB. The future of epigenetic therapy in solid tumours—lessons from the past. Nat Rev Clin Oncol. 2013;10:256–66.
Montalban-Bravo G, Garcia-Manero G. Novel drugs for older patients with acute myeloid leukemia. Leukemia. 2015;29:760–9.
Saunthararajah Y, Sekeres M, Advani A, Mahfouz R, Durkin L, Radivoyevitch T, Englehaupt R, Juersivich J, Cooper K, Husseinzadeh H, et al. Evaluation of noncytotoxic DNMT1-depleting therapy in patients with myelodysplastic syndromes. J Clin Invest. 2015;125:1043–55.
Santos FP, Kantarjian H, Garcia-Manero G, Issa JP, Ravandi F. Decitabine in the treatment of myelodysplastic syndromes. Expert Rev Anticancer Ther. 2010;10:9–22.
Santini V. Novel therapeutic strategies: hypomethylating agents and beyond. Hematology Am Soc Hematol Educ Program. 2012;2012:65–73.
Li X, Mei Q, Nie J, Fu X, Han W. Decitabine: a promising epi-immunotherapeutic agent in solid tumors. Expert Rev Clin Immunol. 2015;11:363–75.
Kadia TM, Thomas XG, Dmoszynska A, Wierzbowska A, Minden M, Arthur C, Delaunay J, Ravandi F, Kantarjian H. Decitabine improves outcomes in older patients with acute myeloid leukemia (AML) and higher blast counts. Am J Hematol. 2015;90:E139–41.
Tsai HC, Li H, Van Neste L, Cai Y, Robert C, Rassool FV, Shin JJ, Harbom KM, Beaty R, Pappou E, et al. Transient low doses of DNA-demethylating agents exert durable antitumor effects on hematological and epithelial tumor cells. Cancer Cell. 2012;21:430–46.
Shen H, Laird PW. In epigenetic therapy, less is more. Cell Stem Cell. 2012;10:353–4.
Nie J, Liu L, Li X, Han W. Decitabine, a new star in epigenetic therapy: the clinical application and biological mechanism in solid tumors. Cancer Lett. 2014;354:12–20.
Lubbert M, Suciu S, Baila L, Ruter BH, Platzbecker U, Giagounidis A, Selleslag D, Labar B, Germing U, Salih HR, et al. Low-dose decitabine versus best supportive care in elderly patients with intermediate- or high-risk myelodysplastic syndrome (MDS) ineligible for intensive chemotherapy: final results of the randomized phase III study of the European Organisation for Research and Treatment of Cancer Leukemia Group and the German MDS Study Group. J Clin Oncol. 2011;29:1987–96.
Cravatt BF, Simon GM, Yates JR 3rd. The biological impact of mass-spectrometry-based proteomics. Nature. 2007;450:991–1000.
Aebersold R, Mann M. Mass spectrometry-based proteomics. Nature. 2003;422:198–207.
Cox J, Mann M. Quantitative, high-resolution proteomics for data-driven systems biology. Annu Rev Biochem. 2011;80:273–99.
Schulze WX, Usadel B. Quantitation in mass-spectrometry-based proteomics. Annu Rev Plant Biol. 2010;61:491–516.
Ong SE, Mann M. Mass spectrometry-based proteomics turns quantitative. Nat Chem Biol. 2005;1:252–62.
Britton LM, Gonzales-Cope M, Zee BM, Garcia BA. Breaking the histone code with quantitative mass spectrometry. Expert Rev Proteomics. 2011;8:631–43.
Ueberheide BM, Mollah S. Deciphering the histone code using mass spectrometry. Int J Mass Spectrom. 2007;259:46–56.
Su X, Ren C, Freitas MA. Mass spectrometry-based strategies for characterization of histones and their post-translational modifications. Expert Rev Proteomics. 2007;4:211–25.
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:933–8.
Meert P, Govaert E, Scheerlinck E, Dhaenens M, Deforce D. Pitfalls in histone propionylation during bottom-up mass spectrometry analysis. Proteomics. 2015;15:2966–71.
Zhang C, Liu Y, Andrews PC. Quantification of histone modifications using 15N metabolic labeling. Methods. 2013;61:236–43.
Lin S, Garcia BA. Examining histone posttranslational modification patterns by high-resolution mass spectrometry. Methods Enzymol. 2012;512:3–28.
Huang H, Lin S, Garcia BA, Zhao Y. Quantitative proteomic analysis of histone modifications. Chem Rev. 2015;115:2376–418.
Moradian A, Kalli A, Sweredoski MJ, Hess S. The top-down, middle-down, and bottom-up mass spectrometry approaches for characterization of histone variants and their post-translational modifications. Proteomics. 2014;14:489–97.
Gao J, Liao R, Yu Y, Zhai H, Wang Y, Sack R, Peters AH, Chen J, Wu H, Huang Z, et al. Absolute quantification of histone PTM marks by MRM-based LC–MS/MS. Anal Chem. 2014;86:9679–86.
Zhang C, Molascon AJ, Gao S, Liu Y, Andrews PC. Quantitative proteomics reveals that the specific methyltransferases Txr1p and Ezl2p differentially affect the mono-, di- and trimethylation states of histone H3 lysine 27 (H3K27). Mol Cell Proteomics. 2013;12:1678–88.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26:1367–72.
Zhang C, Gao S, Molascon AJ, Liu Y, Andrews PC. Quantitative proteomics reveals histone modifications in crosstalk with H3 lysine 27 methylation. Mol Cell Proteomics. 2014;13:749–59.
Zhang C, Gao S, Molascon AJ, Wang Z, Gorovsky MA, Liu Y, Andrews PC. Bioinformatic and proteomic analysis of bulk histones reveals PTM crosstalk and chromatin features. J Proteome Res. 2014;13:3330–7.
Wouters BJ, Delwel R. Epigenetics and approaches to targeted epigenetic therapy in acute myeloid leukemia. Blood. 2016;127:42–52.
Bug G, Ottmann OG. The DAC system and associations with acute leukemias and myelodysplastic syndromes. Invest New Drugs. 2010;28(Suppl 1):S36–49.
Cao R, Zhang Y. The functions of E(Z)/EZH2-mediated methylation of lysine 27 in histone H3. Curr Opin Genet Dev. 2004;14:155–64.
Vire E, Brenner C, Deplus R, Blanchon L, Fraga M, Didelot C, Morey L, Van Eynde A, Bernard D, Vanderwinden JM, et al. The Polycomb group protein EZH2 directly controls DNA methylation. Nature. 2006;439:871–4.
McCabe MT, Ott HM, Ganji G, Korenchuk S, Thompson C, Van Aller GS, Liu Y, Graves AP, Della Pietra A 3rd, Diaz E, et al. EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations. Nature. 2012;492:108–12.
Kitamura T, Tange T, Terasawa T, Chiba S, Kuwaki T, Miyagawa K, Piao YF, Miyazono K, Urabe A, Takaku F. Establishment and characterization of a unique human cell line that proliferates dependently on GM-CSF, IL-3, or erythropoietin. J Cell Physiol. 1989;140:323–34.
Cogle CR, Saki N, Khodadi E, Li J, Shahjahani M, Azizidoost S. Bone marrow niche in the myelodysplastic syndromes. Leuk Res. 2015;39:1020–7.
Martinez-Garcia E, Licht JD. Deregulation of H3K27 methylation in cancer. Nat Genet. 2010;42:100–1.
Schneider R, Bannister AJ, Kouzarides T. Unsafe SETs: histone lysine methyltransferases and cancer. Trends Biochem Sci. 2002;27:396–402.
Klose RJ, Zhang Y. Regulation of histone methylation by demethylimination and demethylation. Nat Rev Mol Cell Biol. 2007;8:307–18.
Vu LP, Luciani L, Nimer SD. Histone-modifying enzymes: their role in the pathogenesis of acute leukemia and their therapeutic potential. Int J Hematol. 2013;97:198–209.
Hock H. A complex Polycomb issue: the two faces of EZH2 in cancer. Genes Dev. 2012;26:751–5.
Mulero-Navarro S, Esteller M. Epigenetic biomarkers for human cancer: the time is now. Crit Rev Oncol Hematol. 2008;68:1–11.
Denis H, Ndlovu MN, Fuks F. Regulation of mammalian DNA methyltransferases: a route to new mechanisms. EMBO Rep. 2011;12:647–56.
Kar S, Deb M, Sengupta D, Shilpi A, Parbin S, Torrisani J, Pradhan S, Patra S. An insight into the various regulatory mechanisms modulating human DNA methyltransferase 1 stability and function. Epigenetics. 2012;7:994–1007.
Qin W, Leonhardt H, Pichler G. Regulation of DNA methyltransferase 1 by interactions and modifications. Nucleus. 2011;2:392–402.
Geiman TM, Sankpal UT, Robertson AK, Zhao Y, Zhao Y, Robertson KD. DNMT3B interacts with hSNF2H chromatin remodeling enzyme, HDACs 1 and 2, and components of the histone methylation system. Biochem Biophys Res Commun. 2004;318:544–55.
Ling Y, Sankpal UT, Robertson AK, McNally JG, Karpova T, Robertson KD. Modification of de novo DNA methyltransferase 3a (Dnmt3a) by SUMO-1 modulates its interaction with histone deacetylases (HDACs) and its capacity to repress transcription. Nucleic Acids Res. 2004;32:598–610.
CZ, GGM, and SH conceived and designed the project; CZ and JS processed samples; YW developed DAC resistant cell lines; HK carried out LC–MS/MS run; CZ performed data analysis; CZ and SH wrote the manuscript. All authors read and approved the manuscript.
This work was supported by The University of Texas MD Anderson Cancer Center’s Moon Shots Program.
The authors declare that they have no competing interests