Cell state-dependent chromatin targeting in NUT carcinoma.

Aberrant transcriptional programming and chromatin dysregulation are common to most cancers. Whether by deranged cell signaling or environmental insult, the resulting oncogenic phenotype is typically manifested in transcriptional changes characteristic of undifferentiated cell growth. Here we analyze targeting of an oncogenic fusion protein, BRD4-NUT, composed of two normally independent chromatin regulators. The fusion causes the formation of large hyperacetylated genomic regions or megadomains, mis-regulation of c-MYC , and an aggressive carcinoma of squamous cell origin. Our previous work revealed largely distinct megadomain locations in different NUT carcinoma patient cell lines. To assess whether this was due to variations in individual genome sequences or epigenetic cell state, we expressed BRD4-NUT in a human stem cell model and found that megadomains formed in dissimilar patterns when comparing cells in the pluripotent state with the same cell line following induction along a mesodermal lineage. Thus, our work implicates initial cell state as the critical factor in the locations of BRD4-NUT megadomains. These results, together with our analysis of c-MYC protein-protein interactions in a patient cell line, are consistent with a cascade of chromatin misregulation underlying NUT carcinoma.

Cell state-dependent chromatin targeting in NUT carcinoma.

Aberrant transcriptional programming and chromatin dysregulation are common to most cancers. Whether by deranged cell signaling or environmental insult, the resulting oncogenic phenotype is typically manifested in transcriptional changes characteristic of undifferentiated cell growth. Here we analyze targeting of an oncogenic fusion protein, BRD4-NUT, composed of 2 normally independent chromatin regulators. The fusion causes the formation of large hyperacetylated genomic regions or megadomains, mis-regulation of c-MYC, and an aggressive carcinoma of squamous cell origin. Our previous work revealed largely distinct megadomain locations in different NUT carcinoma patient cell lines. To assess whether this was due to variations in individual genome sequences or epigenetic cell state, we expressed BRD4-NUT in a human stem cell model and found that megadomains formed in dissimilar patterns when comparing cells in the pluripotent state with the same cell line following induction along a mesodermal lineage. Thus, our work implicates initial cell state as the critical factor in the locations of BRD4-NUT megadomains. These results, together with our analysis of c-MYC protein-protein interactions in a patient cell line, are consistent with a cascade of chromatin misregulation underlying NUT carcinoma.

Variant Polycomb complexes in Drosophila consistent with ancient functional diversity.

Polycomb group (PcG) mutants were first identified in Drosophila on the basis of their failure to maintain proper Hox gene repression during development. The proteins encoded by the corresponding fly genes mainly assemble into one of two discrete Polycomb repressive complexes: PRC1 or PRC2. However, biochemical analyses in mammals have revealed alternative forms of PRC2 and multiple distinct types of noncanonical or variant PRC1. Through a series of proteomic analyses, we identify analogous PRC2 and variant PRC1 complexes in Drosophila, as well as a broader repertoire of interactions implicated in early development. Our data provide strong support for the ancient diversity of PcG complexes and a framework for future analysis in a longstanding and versatile genetic system.

Combined epigenetic and metabolic treatments overcome differentiation blockade in acute myeloid leukemia.

A hallmark of acute myeloid leukemia (AML) is the inability of self-renewing malignant cells to mature into a non-dividing terminally differentiated state. This differentiation block has been linked to dysregulation of multiple cellular processes, including transcriptional, chromatin, and metabolic regulation. The transcription factor HOXA9 and the histone demethylase LSD1 are examples of such regulators that promote differentiation blockade in AML. To identify metabolic targets that interact with LSD1 inhibition to promote myeloid maturation, we screened a small molecule library to identify druggable substrates. We found that differentiation caused by LSD1 inhibition is enhanced by combined perturbation of purine nucleotide salvage and de novo lipogenesis pathways, and identified multiple lines of evidence to support the specificity of these pathways and suggest a potential basis of how perturbation of these pathways may interact synergistically to promote myeloid differentiation. In sum, these findings suggest potential drug combination strategies in the treatment of AML.

Re-programing Chromatin with a Bifunctional LSD1/HDAC Inhibitor Induces Therapeutic Differentiation in DIPG.

H3K27M mutations resulting in epigenetic dysfunction are frequently observed in diffuse intrinsic pontine glioma (DIPGs), an incurable pediatric cancer. We conduct a CRISPR screen revealing that knockout of KDM1A encoding lysine-specific demethylase 1 (LSD1) sensitizes DIPG cells to histone deacetylase (HDAC) inhibitors. Consistently, Corin, a bifunctional inhibitor of HDACs and LSD1, potently inhibits DIPG growth in vitro and in xenografts. Mechanistically, Corin increases H3K27me3 levels suppressed by H3K27M histones, and simultaneously increases HDAC-targeted H3K27ac and LSD1-targeted H3K4me1 at differentiation-associated genes. Corin treatment induces cell death, cell-cycle arrest, and a cellular differentiation phenotype and drives transcriptional changes correlating with increased survival time in DIPG patients. These data suggest a strategy for treating DIPG by simultaneously inhibiting LSD1 and HDACs.

Intersection of Epigenetic and Metabolic Regulation of Histone Modifications in Acute Myeloid Leukemia.

Acute myeloid leukemia (AML) is one of the most lethal blood cancers, accounting for close to a quarter of a million annual deaths worldwide. Even though genetically heterogeneous, all AMLs are characterized by two interrelated features-blocked differentiation and high proliferative capacity. Despite significant progress in our understanding of the molecular and genetic basis of AML, the treatment of AMLs with chemotherapeutic regimens has remained largely unchanged in the past 30 years. In this review, we will consider the role of two cellular processes, metabolism and epigenetics, in the development and progression of AML and highlight the studies that suggest an interconnection of therapeutic importance between the two. Large-scale whole-exome sequencing of AML patients has revealed the presence of mutations, translocations or duplications in several epigenetic effectors such as DNMT3, MLL, ASXL1, and TET2, often times co-occuring with mutations in metabolic enzymes such as IDH1 and IDH2. These mutations often result in impaired enzymatic activity which leads to an altered epigenetic landscape through dysregulation of chromatin modifications such as DNA methylation, histone acetylation and methylation. We will discuss the role of enzymes that are responsible for establishing these modifications, namely histone acetyl transferases (HAT), histone methyl transferases (HMT), demethylases (KDMs), and deacetylases (HDAC), and also highlight the merits and demerits of using inhibitors that target these enzymes. Furthermore, we will tie in the metabolic regulation of co-factors such as acetyl-CoA, SAM, and α-ketoglutarate that are utilized by these enzymes and examine the role of metabolic inhibitors as a treatment option for AML. In doing so, we hope to stimulate interest in this topic and help generate a rationale for the consideration of the combinatorial use of metabolic and epigenetic inhibitors for the treatment of AML.

Bivalent complexes of PRC1 with orthologs of BRD4 and MOZ/MORF target developmental genes in Drosophila.

Regulatory decisions in Drosophila require Polycomb group (PcG) proteins to maintain the silent state and Trithorax group (TrxG) proteins to oppose silencing. Since PcG and TrxG are ubiquitous and lack apparent sequence specificity, a long-standing model is that targeting occurs via protein interactions; for instance, between repressors and PcG proteins. Instead, we found that Pc-repressive complex 1 (PRC1) purifies with coactivators Fs(1)h [female sterile (1) homeotic] and Enok/Br140 during embryogenesis. Fs(1)h is a TrxG member and the ortholog of BRD4, a bromodomain protein that binds to acetylated histones and is a key transcriptional coactivator in mammals. Enok and Br140, another bromodomain protein, are orthologous to subunits of a mammalian MOZ/MORF acetyltransferase complex. Here we confirm PRC1-Br140 and PRC1-Fs(1)h interactions and identify their genomic binding sites. PRC1-Br140 bind developmental genes in fly embryos, with analogous co-occupancy of PRC1 and a Br140 ortholog, BRD1, at bivalent loci in human embryonic stem (ES) cells. We propose that identification of PRC1-Br140 "bivalent complexes" in fly embryos supports and extends the bivalency model posited in mammalian cells, in which the coexistence of H3K4me3 and H3K27me3 at developmental promoters represents a poised transcriptional state. We further speculate that local competition between acetylation and deacetylation may play a critical role in the resolution of bivalent protein complexes during development.

Ectopic protein interactions within BRD4-chromatin complexes drive oncogenic megadomain formation in NUT midline carcinoma.

To investigate the mechanism that drives dramatic mistargeting of active chromatin in NUT midline carcinoma (NMC), we have identified protein interactions unique to the BRD4-NUT fusion oncoprotein compared with wild-type BRD4. Using cross-linking, affinity purification, and mass spectrometry, we identified the EP300 acetyltransferase as uniquely associated with BRD4 through the NUT fusion in both NMC and non-NMC cell types. We also discovered ZNF532 associated with BRD4-NUT in NMC patient cells but not detectable in 293T cells. EP300 and ZNF532 are both implicated in feed-forward regulatory loops leading to propagation of the oncogenic chromatin complex in BRD4-NUT patient cells. Adding key functional significance to our biochemical findings, we independently discovered a ZNF532-NUT translocation fusion in a newly diagnosed NMC patient. ChIP sequencing of the major players NUT, ZNF532, BRD4, EP300, and H3K27ac revealed the formation of ZNF532-NUT-associated hyperacetylated megadomains, distinctly localized but otherwise analogous to those found in BRD4-NUT patient cells. Our results support a model in which NMC is dependent on ectopic NUT-mediated interactions between EP300 and components of BRD4 regulatory complexes, leading to a cascade of misregulation.

upSET, the Drosophila homologue of SET3, Is Required for Viability and the Proper Balance of Active and Repressive Chromatin Marks.

Chromatin plays a critical role in faithful implementation of gene expression programs. Different post-translational modifications (PTMs) of histone proteins reflect the underlying state of gene activity, and many chromatin proteins write, erase, bind, or are repelled by, these histone marks. One such protein is UpSET, the Drosophila homolog of yeast Set3 and mammalian KMT2E (MLL5). Here, we show that UpSET is necessary for the proper balance between active and repressed states. Using CRISPR/Cas-9 editing, we generated S2 cells that are mutant for upSET We found that loss of UpSET is tolerated in S2 cells, but that heterochromatin is misregulated, as evidenced by a strong decrease in H3K9me2 levels assessed by bulk histone PTM quantification. To test whether this finding was consistent in the whole organism, we deleted the upSET coding sequence using CRISPR/Cas-9, which we found to be lethal in both sexes in flies. We were able to rescue this lethality using a tagged upSET transgene, and found that UpSET protein localizes to transcriptional start sites (TSS) of active genes throughout the genome. Misregulated heterochromatin is apparent by suppressed position effect variegation of the wm4 allele in heterozygous upSET-deleted flies. Using nascent-RNA sequencing in the upSET-mutant S2 lines, we show that this result applies to heterochromatin genes generally. Our findings support a critical role for UpSET in maintaining heterochromatin, perhaps by delimiting the active chromatin environment.

The Oncoprotein BRD4-NUT Generates Aberrant Histone Modification Patterns.

Defects in chromatin proteins frequently manifest in diseases. A striking case of a chromatin-centric disease is NUT-midline carcinoma (NMC), which is characterized by expression of NUT as a fusion partner most frequently with BRD4. ChIP-sequencing studies from NMC patients revealed that BRD4-NUT (B4N) covers large genomic regions and elevates transcription within these domains. To investigate how B4N modulates chromatin, we performed affinity purification of B4N when ectopically expressed in 293-TREx cells and quantified the associated histone posttranslational modifications (PTM) using proteomics. We observed significant enrichment of acetylation particularly on H3 K18 and of combinatorial patterns such as H3 K27 acetylation paired with K36 methylation. We postulate that B4N complexes override the preexisting histone code with new PTM patterns that reflect aberrant transcription and that epigenetically modulate the nucleosome environment toward the NMC state.

Streamlined discovery of cross-linked chromatin complexes and associated histone modifications by mass spectrometry.

Posttranslational modifications (PTMs) are key contributors to chromatin function. The ability to comprehensively link specific histone PTMs with specific chromatin factors would be an important advance in understanding the functions and genomic targeting mechanisms of those factors. We recently introduced a cross-linked affinity technique, BioTAP-XL, to identify chromatin-bound protein interactions that can be difficult to capture with native affinity techniques. However, BioTAP-XL was not strictly compatible with similarly comprehensive analyses of associated histone PTMs. Here we advance BioTAP-XL by demonstrating the ability to quantify histone PTMs linked to specific chromatin factors in parallel with the ability to identify nonhistone binding partners. Furthermore we demonstrate that the initially published quantity of starting material can be scaled down orders of magnitude without loss in proteomic sensitivity. We also integrate hydrophilic interaction chromatography to mitigate detergent carryover and improve liquid chromatography-mass spectrometric performance. In summary, we greatly extend the practicality of BioTAP-XL to enable comprehensive identification of protein complexes and their local chromatin environment.

Sex comb on midleg (Scm) is a functional link between PcG-repressive complexes in Drosophila.

The Polycomb group (PcG) proteins are key regulators of development in Drosophila and are strongly implicated in human health and disease. How PcG complexes form repressive chromatin domains remains unclear. Using cross-linked affinity purifications of BioTAP-Polycomb (Pc) or BioTAP-Enhancer of zeste [E(z)], we captured all PcG-repressive complex 1 (PRC1) or PRC2 core components and Sex comb on midleg (Scm) as the only protein strongly enriched with both complexes. Although previously not linked to PRC2, we confirmed direct binding of Scm and PRC2 using recombinant protein expression and colocalization of Scm with PRC1, PRC2, and H3K27me3 in embryos and cultured cells using ChIP-seq (chromatin immunoprecipitation [ChIP] combined with deep sequencing). Furthermore, we found that RNAi knockdown of Scm and overexpression of the dominant-negative Scm-SAM (sterile α motif) domain both affected the binding pattern of E(z) on polytene chromosomes. Aberrant localization of the Scm-SAM domain in long contiguous regions on polytene chromosomes revealed its independent ability to spread on chromatin, consistent with its previously described ability to oligomerize in vitro. Pull-downs of BioTAP-Scm captured PRC1 and PRC2 and additional repressive complexes, including PhoRC, LINT, and CtBP. We propose that Scm is a key mediator connecting PRC1, PRC2, and transcriptional silencing. Combined with previous structural and genetic analyses, our results strongly suggest that Scm coordinates PcG complexes and polymerizes to produce broad domains of PcG silencing.

Proteogenomics analysis reveals specific genomic orientations of distal regulatory regions composed by non-canonical histone variants.

BACKGROUND: Histone variants play further important roles in DNA packaging and controlling gene expression. However, our understanding about their composition and their functions is limited. RESULTS: Integrating proteomic and genomic approaches, we performed a comprehensive analysis of the epigenetic landscapes containing the four histone variants H3.1, H3.3, H2A.Z, and macroH2A. These histones were FLAG-tagged in HeLa cells and purified using chromatin immunoprecipitation (ChIP). By adopting ChIP followed by mass spectrometry (ChIP-MS), we quantified histone post-translational modifications (PTMs) and histone variant nucleosomal ratios in highly purified mononucleosomes. Subsequent ChIP followed by next-generation sequencing (ChIP-seq) was used to map the genome-wide localization of the analyzed histone variants and define their chromatin domains. Finally, we included in our study large datasets contained in the ENCODE database. We newly identified a group of regulatory regions enriched in H3.1 and the histone variant associated with repressive marks macroH2A. Systematic analysis identified both symmetric and asymmetric patterns of histone variant occupancies at intergenic regulatory regions. Strikingly, these directional patterns were associated with RNA polymerase II (PolII). These asymmetric patterns correlated with the enhancer activities measured using global run-on sequencing (GRO-seq) data. CONCLUSIONS: Our studies show that H2A.Z and H3.3 delineate the orientation of transcription at enhancers as observed at promoters. We also showed that enhancers with skewed histone variant patterns well facilitate enhancer activity. Collectively, our study indicates that histone variants are deposited at regulatory regions to assist gene regulation.

BioTAP-XL: Cross-linking/Tandem Affinity Purification to Study DNA Targets, RNA, and Protein Components of Chromatin-Associated Complexes.

In order to understand how chromatin complexes function in the nucleus, it is important to obtain a comprehensive picture of their protein, DNA, and RNA components, as well as their mutual interactions. This unit presents a chromatin cross-linking approach (BioTAP-XL) that utilizes a special BioTAP-tagged transgenic protein bait along with mass spectrometry to identify protein complex components, and high-throughput sequencing to identify RNA components and DNA binding sites. Full protocols are provided for Drosophila cells and for human cells in culture, along with an additional protocol for Drosophila embryos as the source material. A key element of the approach in all cases is the generation of control data from input chromatin samples.

High resolution is not a strict requirement for characterization and quantification of histone post-translational modifications.

Mass spectrometry (MS) is a powerful tool to accurately identify and quantify histone post-translational modifications (PTMs). High-resolution mass analyzers have been regarded as essential for these PTM analyses because the mass accuracy afforded is sufficient to differentiate trimethylation versus acetylation (42.0470 and 42.0106 Da, respectively), whereas lower-resolution mass analyzers cannot. Noting this limitation, we sought to determine whether lower-resolution detectors are nonetheless adequate for histone PTM analysis by comparing the low-resolution LTQ Velos Pro with the high-resolution LTQ-Orbitrap Velos Pro. We first determined that the optimal scan mode on the LTQ Velos Pro is the Enhanced scan mode with respect to apparent resolution, number of MS and MS/MS scans per run, and reproducibility of label-free quantifications. We next compared the performance of the LTQ Velos Pro to the LTQ-Orbitrap Velos Pro using the same criteria for comparison, and we found that the main difference is that the LTQ-Orbitrap Velos Pro is able to resolve the difference between acetylation and trimethylation while the LTQ Velos Pro cannot. However, using heavy isotope labeled synthetic peptide standards and retention time information enables confident assignment of these modifications and comparable quantification between the instruments. Therefore, lower-resolution instruments can confidently be utilized for histone PTM analysis.

Histone methylation has dynamics distinct from those of histone acetylation in cell cycle reentry from quiescence.

Cell growth is attuned to nutrient availability to sustain homeostatic biosynthetic processes. In unfavorable environments, cells enter a nonproliferative state termed quiescence but rapidly return to the cell cycle once conditions support energetic needs. Changing cellular metabolite pools are proposed to directly alter the epigenome via histone acetylation. Here we studied the relationship between histone modification dynamics and the dramatic transcriptional changes that occur during nutrient-induced cell cycle reentry from quiescence in the yeast Saccharomyces cerevisiae. SILAC (stable isotope labeling by amino acids in cell culture)-based mass spectrometry showed that histone methylation-in contrast to histone acetylation-is surprisingly static during quiescence exit. Chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq) revealed genome-wide shifts in histone acetylation at growth and stress genes as cells exit quiescence and transcription dramatically changes. Strikingly, however, the patterns of histone methylation remain intact. We conclude that the functions of histone methylation and acetylation are remarkably distinct during quiescence exit: acetylation rapidly responds to metabolic state, while methylation is independent. Thus, the initial burst of growth gene reactivation emerging from quiescence involves dramatic increases of histone acetylation but not of histone methylation.

Heterochromatin-associated interactions of Drosophila HP1a with dADD1, HIPP1, and repetitive RNAs.

Heterochromatin protein 1 (HP1a) has conserved roles in gene silencing and heterochromatin and is also implicated in transcription, DNA replication, and repair. Here we identify chromatin-associated protein and RNA interactions of HP1a by BioTAP-XL mass spectrometry and sequencing from Drosophila S2 cells, embryos, larvae, and adults. Our results reveal an extensive list of known and novel HP1a-interacting proteins, of which we selected three for validation. A strong novel interactor, dADD1 (Drosophila ADD1) (CG8290), is highly enriched in heterochromatin, harbors an ADD domain similar to human ATRX, displays selective binding to H3K9me2 and H3K9me3, and is a classic genetic suppressor of position-effect variegation. Unexpectedly, a second hit, HIPP1 (HP1 and insulator partner protein-1) (CG3680), is strongly connected to CP190-related complexes localized at putative insulator sequences throughout the genome in addition to its colocalization with HP1a in heterochromatin. A third interactor, the histone methyltransferase MES-4, is also enriched in heterochromatin. In addition to these protein-protein interactions, we found that HP1a selectively associated with a broad set of RNAs transcribed from repetitive regions. We propose that this rich network of previously undiscovered interactions will define how HP1a complexes perform their diverse functions in cells and developing organisms.

De novo peptide design and experimental validation of histone methyltransferase inhibitors.

Histones are small proteins critical to the efficient packaging of DNA in the nucleus. DNA­protein complexes, known as nucleosomes, are formed when the DNA winds itself around the surface of the histones. The methylation of histone residues by enhancer of zeste homolog 2 (EZH2) maintains gene repression over successive cell generations. Overexpression of EZH2 can silence important tumor suppressor genes leading to increased invasiveness of many types of cancers. This makes the inhibition of EZH2 an important target in the development of cancer therapeutics. We employed a three-stage computational de novo peptide design method to design inhibitory peptides of EZH2. The method consists of a sequence selection stage and two validation stages for fold specificity and approximate binding affinity. The sequence selection stage consists of an integer linear optimization model that was solved to produce a rank-ordered list of amino acid sequences with increased stability in the bound peptide-EZH2 structure. These sequences were validated through the calculation of the fold specificity and approximate binding affinity of the designed peptides. Here we report the discovery of novel EZH2 inhibitory peptides using the de novo peptide design method. The computationally discovered peptides were experimentally validated in vitro using dose titrations and mechanism of action enzymatic assays. The peptide with the highest in vitro response, SQ037, was validated in nucleo using quantitative mass spectrometry-based proteomics. This peptide had an IC50 of 13.5 mM, demonstrated greater potency as an inhibitor when compared to the native and K27A mutant control peptides, and demonstrated competitive inhibition versus the peptide substrate. Additionally, this peptide demonstrated high specificity to the EZH2 target in comparison to other histone methyltransferases. The validated peptides are the first computationally designed peptides that directly inhibit EZH2. These inhibitors should prove useful for further chromatin biology investigations.

De novo peptide design and experimental validation of histone methyltransferase inhibitors.

Histones are small proteins critical to the efficient packaging of DNA in the nucleus. DNA-protein complexes, known as nucleosomes, are formed when the DNA winds itself around the surface of the histones. The methylation of histone residues by enhancer of zeste homolog 2 (EZH2) maintains gene repression over successive cell generations. Overexpression of EZH2 can silence important tumor suppressor genes leading to increased invasiveness of many types of cancers. This makes the inhibition of EZH2 an important target in the development of cancer therapeutics. We employed a three-stage computational de novo peptide design method to design inhibitory peptides of EZH2. The method consists of a sequence selection stage and two validation stages for fold specificity and approximate binding affinity. The sequence selection stage consists of an integer linear optimization model that was solved to produce a rank-ordered list of amino acid sequences with increased stability in the bound peptide-EZH2 structure. These sequences were validated through the calculation of the fold specificity and approximate binding affinity of the designed peptides. Here we report the discovery of novel EZH2 inhibitory peptides using the de novo peptide design method. The computationally discovered peptides were experimentally validated in vitro using dose titrations and mechanism of action enzymatic assays. The peptide with the highest in vitro response, SQ037, was validated in nucleo using quantitative mass spectrometry-based proteomics. This peptide had an IC50 of 13.5 [Formula: see text]M, demonstrated greater potency as an inhibitor when compared to the native and K27A mutant control peptides, and demonstrated competitive inhibition versus the peptide substrate. Additionally, this peptide demonstrated high specificity to the EZH2 target in comparison to other histone methyltransferases. The validated peptides are the first computationally designed peptides that directly inhibit EZH2. These inhibitors should prove useful for further chromatin biology investigations.

A quantitative atlas of histone modification signatures from human cancer cells.

BACKGROUND: An integral component of cancer biology is the understanding of molecular properties uniquely distinguishing one cancer type from another. One class of such properties is histone post-translational modifications (PTMs). Many histone PTMs are linked to the same diverse nuclear functions implicated in cancer development, including transcriptional activation and epigenetic regulation, which are often indirectly assayed with standard genomic technologies. Thus, there is a need for a comprehensive and quantitative profiling of cancer lines focused on their chromatin modification states. RESULTS: To complement genomic expression profiles of cancer lines, we report the proteomic classification of 24 different lines, the majority of which are cancer cells, by quantifying the abundances of a large panel of single and combinatorial histone H3 and H4 PTMs, and histone variants. Concurrent to the proteomic analysis, we performed transcriptomic analysis on histone modifying enzyme abundances as a proxy for quantifying their activity levels. While the transcriptomic and proteomic results were generally consistent in terms of predicting histone PTM abundance from enzyme abundances, several PTMs were regulated independently of the modifying enzyme expression. In addition, combinatorial PTMs containing H3K27 methylation were especially enriched in breast cell lines. Knockdown of the predominant H3K27 methyltransferase, enhancer of zeste 2 (EZH2), in a mouse mammary xenograft model significantly reduced tumor burden in these animals and demonstrated the predictive utility of proteomic techniques. CONCLUSIONS: Our proteomic and genomic characterizations of the histone modification states provide a resource for future investigations of the epigenetic and non-epigenetic determinants for classifying and analyzing cancer cells.