Dual EZH2 and EHMT2 histone methyltransferase inhibition increases biological efficacy in breast cancer cells
- Edward Curry1,
- Ian Green2,
- Nadine Chapman-Rothe1,
- Elham Shamsaei1,
- Sarah Kandil1,
- Fanny L Cherblanc3,
- Luke Payne2,
- Emma Bell2,
- Thota Ganesh4,
- Nitipol Srimongkolpithak3,
- Joachim Caron3,
- Fengling Li5,
- Anthony G. Uren6,
- James P. Snyder7,
- Masoud Vedadi5,
- Matthew J. Fuchter3Email author and
- Robert Brown1, 8Email author
© Curry et al. 2015
Received: 12 December 2014
Accepted: 28 July 2015
Published: 21 August 2015
Many cancers show aberrant silencing of gene expression and overexpression of histone methyltransferases. The histone methyltransferases (HKMT) EZH2 and EHMT2 maintain the repressive chromatin histone methylation marks H3K27me and H3K9me, respectively, which are associated with transcriptional silencing. Although selective HKMT inhibitors reduce levels of individual repressive marks, removal of H3K27me3 by specific EZH2 inhibitors, for instance, may not be sufficient for inducing the expression of genes with multiple repressive marks.
We report that gene expression and inhibition of triple negative breast cancer cell growth (MDA-MB-231) are markedly increased when targeting both EZH2 and EHMT2, either by siRNA knockdown or pharmacological inhibition, rather than either enzyme independently. Indeed, expression of certain genes is only induced upon dual inhibition. We sought to identify compounds which showed evidence of dual EZH2 and EHMT2 inhibition. Using a cell-based assay, based on the substrate competitive EHMT2 inhibitor BIX01294, we have identified proof-of-concept compounds that induce re-expression of a subset of genes consistent with dual HKMT inhibition. Chromatin immunoprecipitation verified a decrease in silencing marks and an increase in permissive marks at the promoter and transcription start site of re-expressed genes, while Western analysis showed reduction in global levels of H3K27me3 and H3K9me3. The compounds inhibit growth in a panel of breast cancer and lymphoma cell lines with low to sub-micromolar IC50s. Biochemically, the compounds are substrate competitive inhibitors against both EZH2 and EHMT1/2.
We have demonstrated that dual inhibition of EZH2 and EHMT2 is more effective at eliciting biological responses of gene transcription and cancer cell growth inhibition compared to inhibition of single HKMTs, and we report the first dual EZH2-EHMT1/2 substrate competitive inhibitors that are functional in cells.
EZH2 along with EED and SUZ12 are the indispensible core components of the Polycomb Repressive Complex (PRC2) responsible for maintenance of the repressive epigenetic mark H3K27me3: trimethylation of lysine 27 of histone 3 . High expression of the histone methyltransferase (HKMT) EZH2, in some cases associated with gene amplification, has been well documented in a variety of cancers , . EZH2 overexpression has been linked to poor prognosis [4, 5] and shown to be a marker of aggressive breast cancer , associated with difficult-to-treat basal or triple negative breast cancer . Gene knockdown of EZH2 reduces growth of a variety of tumour cell types [5, 8, 9]. Several groups have reported specific co-factor competitive EZH2 inhibitors [10–16], which have shown a strong capacity to reduce growth of cells expressing mutated forms of EZH2 (such as certain non-Hodgkin’s lymphoma ). However, removal of the repressive mark H3K27me3 alone may not always be sufficient for reversal of gene silencing. Indeed, it has been shown that highly specific EZH2 inhibitors require a mutant EZH2 status to inhibit cell growth, being less effective in cells solely expressing wild type EZH2 [5, 8, 9]. Elimination of further repressive methylation marks by inhibition of additional HKMTs may be required to fully realise the epigenetic potential of HKMT inhibitors.
EHMT2 (also known as G9a) and the highly homologous EHMT1 (also known as GLP) are HKMTs partly responsible for mono- and di-methylation of lysine nine of histone 3 (H3K9me1 and H3K9me2, respectively); repressive chromatin marks found on the promoter regions of genes that are often aberrantly silenced in cancer . EHMT2 is overexpressed and amplified in various cancers including leukaemia, prostate carcinoma, and lung cancer, with gene knockdown of EHMT2 inhibiting cancer cell growth in these tumour types [18, 19]. BIX-01294 (see Fig. 2) was previously identified as an inhibitor of the HKMTs EHMT2 and EHMT1, and subsequent medicinal chemistry studies around the 2, 4-diamino-6, 7-dimethoxyquinazoline template of BIX-01294 have yielded a number of follow-up EHMT2 inhibitors [20–25].
In addition to its role in methylating H3K9, EHMT2 has been shown to be able to methylate H3K27 [26, 27]. It has been suggested that this could provide cells with a mechanism to compensate in part for a loss of EZH2 . The picture is further complicated by recent evidence that EHMT2 and EZH2 (via the PRC2 complex) interact physically and share targets for epigenetic silencing . Combining this evidence, it would again suggest that specifically targeting either EZH2 or EHMT2 alone may not be sufficient to reverse epigenetic silencing of genes, but rather combined inhibition may be required. To this end, we have examined the effect of dual EZH2 and EHMT2 gene knockdown or enzyme inhibition in breast cancer cells. Consistent with the requirement for removal of both repressive H3K9 and H3K27 methylation marks, we show that dual inhibition of EHMT2 and EZH2 pharmacologically or by SiRNA is necessary for reactivation of certain genes and induces greater inhibition of cell growth than targeting either HKMT alone in triple negative breast cancer MDA-MB-231 cells. Further, we have identified proof-of-concept compounds which are dual (substrate competitive) EZH2-EHMT1/2 inhibitors.
Combined inhibition of EZH2 and EHMT2 is more effective at inducing gene re-expression and inhibiting tumour cell growth than single HKMT inhibition
Next, the effects on cell viability of GSK343 and UNC0638 used alone or in combination were examined (Fig. 1c). Treatment alone with GSK343 showed no significant reduction in cell viability up to 15 μM, while UNC0638 sole treatment caused a dose dependant reduction in cell viability, with a calculated IC50 of 9 μM. When the cells were treated with both compounds in combination, a marked increase in growth inhibition was observed when compared to single agent treatment using UNC0638 or GSK343 (Fig. 1c). This is particularly apparent at a 5 μM concentration of both compounds, where alone they have no significant effect on reducing cell viability, while in combination, they markedly reduce cell viability to >50 % (p < 0.01).
Analogues of an EHMT2-specific inhibitor can upregulate EZH2 silenced genes
Both EZH2 and EHMT1/2 belong to the SET-domain superfamily , the catalytic SET-domain being responsible for the methylation of the targeted lysine residues. BIX-01294 has previously been shown, both structurally and biochemically, to bind to the substrate (histone)-binding pocket of EHMT1/2 . Since protein recognition motifs for histone binding at repressive sites are similar  and EHMT2 has been shown to be able to methylate H3K27, in addition to its more common H3K9 target , it is likely that there are common aspects to the histone substrate binding pockets of the repressive HKMTs EZH2 and EHMT1/2. We therefore felt it would be feasible to use quinazoline template of BIX-01294 in the discovery of dual (substrate competitive) EZH2-EHMT1/2 inhibitors.
RT-PCR data for a single dose of a panel of HKMT inhibitor compounds
Chromatin immunoprecipitation (ChIP) experiments were carried out on treated MDA-MB-231 cells to verify that the detected gene upregulation is indeed due to chromatin remodelling (Fig. 3b). We tested the silencing marks H3K9me3 and H3K27me3 as well as the activating marks H3K4me3, H3K4me2, H3K27ac and H3K9ac. All three compounds showed a clear decrease in repressive chromatin marks (H3K27me3, H3K9me3), and at least in some instances, an increase in permissive marks, at two target genes (Fig. 3b). This is consistent with the compounds having dual HKMT inhibitory activity in removing both H3K9me and H3K27me marks, while allowing activating marks to be established at these loci.
Genome-wide changes in gene expression
The same enrichment tests were repeated using target gene sets identified in an EZH2 siRNA knockdown study in another breast cancer cell line, MCF-7 . Almost no enrichment was observed of this gene set in MDA-MB-231 cells after treatment with any of the compounds (HKMTI-1-005, GSK343 and UNC0638) (Fig. 4a), suggesting that EZH2 has cell-type-specific targets. To investigate this further, we undertook a meta-analysis to identify consensus target genes based on 18 independent EZH2 siRNA studies (details of the meta-analysis are provided in ‘Methods’ also see Additional file 4: Table S5). Encouragingly, treatment of MDA-MB-231 cells with HKMTI-1-005 resulted in highly significant upregulation of these consensus EZH2-repressed genes (Fig. 4a). This suggests that key EZH2 target genes that are conserved across a wide range of cell lines are re-expressed upon treatment with our dual HKMT inhibitor. Furthermore, this identifies generally applicable pharmacodynamic biomarkers of EZH2 inhibitors across cell types.
Compound-induced changes in H3K9me and H3K27me levels in cells
Cell growth IC50 of HKMT inhibitors in MDA-MB-231 cells
Cell growth assay (IC50 μM)
Cell viability (MTT)
In order to identify specific pathways being transcriptionally modulated, the microarray data was analysed for enrichment of pathways belonging to each pathway listed on the ConsensusPathDB (CPDB) database . The Benjamini-Hochberg-adjusted  enrichment p value estimates for each treatment is given in Additional file 5: Table S6. Interestingly, genes belonging to the pathway ‘Apoptosis’ displayed a highly significant systematic shift towards upregulation on treatment with our hit compound(s) at 24 h (p < 1E-4) but not the selective EZH2 (GSK343) or EHMT2 (UNC0638) inhibitor compound (p = 0.42 and p = 0.30, respectively). Consistent with induction of apoptosis-related genes, hit compound HKMTI-1-005 induces apoptosis in MDA-MB-231 cells in a dose-dependent manner, as measured by Caspase 3/7 activity (Additional file 6: Figure S2).
Cell growth inhibition induced by HKMT inhibitors
Cell growth IC50 inhibitors in a panel of cell lines
(A) Lymphoma cell lines
(B) Breast cancer cell lines
The anti-proliferative effect of HKMTI-1-005 on a small panel of breast cancer cell lines was determined, with IC50 values in the range 2–10 μM (Table 3). All of the cancer breast cell lines examined were found to be more sensitive to HKMTI-1-005 compared to a normal breast epithelial cell line MCF10a. The breast cancer cell line BT-474, which is the cell line most sensitive to HKMTI-1-005 treatment, has the highest relative expression of EZH2, as detected by Western analysis (data not shown).
Hit compounds directly inhibit EZH2 and EHMT1/2 and are substrate competitive inhibitors
We have previously reported the EHMT2 IC50 of HKMTI-1-005, HKMTI-1-011 and HKMTI-1-022 to be 0.10, 3.19 and 0.47 μM, respectively . This data was generated using a scintillation proximity assay (SPA) which monitors the transfer of a tritium-labelled methyl group from [3H]S-adenosyl-L-methionine (SAM) to a biotinylated-H3 (1–25) peptide substrate, mediated by EHMT2. A comparable PRC2 enzymatic assay was employed here to assess biochemical inhibitory activity of our hits against EZH2. A trimeric PRC2 complex (EZH2:EED:SUZ12) was employed in this assay, along with a biotinylated-H3 (21–44) peptide substrate. This revealed HKMTI-1-005, HKMTI-1-011 and HKMTI-1-022 to have PRC2 IC50 values of 24, 12 and 16 μM, respectively, under these assay conditions (see Additional file 7: Figure S3). Since the peptide substrates used in these assays are poor models for the complex and dynamic structure of the chromatin substrate in cells and since the only minimal number of PRC2 proteins (EZH2:EED:SUZ12) required for enzymatically active EZH2 were employed in the PRC2 assay, care should be taken in the over interpretation of this in vitro inhibitory data. Nonetheless, we note that both the EHMT2 and PRC2 biochemical potency is comparable to the inhibitory concentrations employed in our cell-based assays.
Perhaps more importantly, in accordance with our design rationale, mechanism of inhibition studies on representative hit HKMTI-1-005 revealed it to have a well-defined, peptide substrate competitive mechanism of action (see Additional file 8: Figure S4), in contrast to all known EZH2 inhibitory chemotypes. Broad screening of our compound library against PRC2 using this assay revealed the IC50 values obtained for all actives to be dependent on peptide substrate concentration (data not shown), further confirming a substrate competitive inhibitory mode for this chemotype.
Finally, a methyltransferase selectivity screen was carried out for the hits on a panel of enzymes including 11 HKMTs, 3 protein arginine methyltransferases (PRMTs) and 1 DNA methyl-transferase (DNMT) (Additional file 9: Figure S5). None of the hits had any significant inhibitory activity against these 15 other methyltransferase targets (up to 100 μM), confirming them to be selective for EZH2 and EHMT1/2. Taken together, these data reveal our hit compounds to be dual EZH2 and EHMT1/2 inhibitors with a substrate competitive mechanism of action.
It is widely accepted that the installation, maintenance and functional output of epigenetic modifications occur in concert via combinatorial sets of modifications. Therefore, removal of a single specific repressive mark may not alone be sufficient for reversal of gene silencing. Elimination of multiple repressive methylation marks may instead be required to re-express a wider spectrum of genes. Given the complexities of epigenetic regulation and cross-talk between epigenetic regulators, the discovery of inhibitors of epigenetic processes that lead to reversal of epigenetic silencing may be more suited to cell-based methods measuring reactivation of a panel of target genes, rather than cell-free assays that use purified components. Through the use of a breast cancer (MDA-MB-231) cell assay based on the re-expression of epigenetically silenced genes, we report the identification of hit compounds that phenocopy the effects of dual EZH2/EHMT2 pharmacological inhibition and dual SiRNA gene knockdown.
The recently reported specific EZH2 inhibitors are all co-factor competitive, the majority of which have converged to a common chemotype (Fig. 2) [10–16]. Conversely, the dual EZH2/EHMT2 inhibitors we here report are substrate competitive. Not only do these represent the first inhibitors uniquely targeting the substrate binding site of EZH2 but also confirm our original hypothesis that the histone-binding sites of certain HKMTs are similar , and it is therefore possible to discover dual inhibitors targeting this supposedly divergent pocket. Indeed, the results herein suggest that there are common aspects to the histone-binding pockets of the repressive HKMTs, EZH2 and EHMT1/2, different from other HKMTs. Indeed, our selectivity data suggest EZH2 and EHMT1/2 to be the sole HKMT targets of our hit compounds, as does our cell-based data. It is interesting that small changes to the chemical structure of these molecules endow our hits with dual activity, something not observed for the structurally related UNC0638. Indeed, quinazoline EHMT2 inhibitors UNC0638   and UNC0642  have been previously shown not to significantly inhibit EZH2 in biochemical assays.
Amplification or overexpression of EZH2 has been observed in a wide range of tumour types [3–8]. Furthermore, it has been proposed that epigenetic dysregulation can be a contributing factor to acquired drug resistance [7, 8, 41]. In cancers, the specific signalling mechanisms that lead to rapid tumour cell proliferation or evasion of drug-induced apoptosis may vary from cell to cell. One of the appeals of epigenetic therapies in cancer is that, rather than trying to target each individual signalling aberration, the target is the means of acquiring aberrant signalling. Therefore, it is hoped that such therapies may fare better in a heterogeneous tumour environment than drugs targeting specific signalling proteins. In this light, we highlight the observation that a set of EZH2 target genes derived from siRNA knockdown in MDA-MB-231 cells were systematically upregulated following treatment of MDA-MB-231 cells with HKMTI-1-005, but not a set of EZH2 targets were identified from siRNA knockdown in MCF7 cells. This suggests that the compounds are able to elicit a transcriptional response that is specific to a particular cell line, and thus, represent a means of tailoring the response to the targets that are specifically epigenetically repressed in the cancer cells to be treated. However, this fact additionally suggests that it may be difficult to find generally appropriate pharmacodynamics biomarkers indicative of a cellular response to treatment with the compounds. To address this, we carried out a meta-analysis to identify genes with a consistent upregulation following EZH2 knockdown via siRNA across a panel of 18 cell lines. These genes may reflect useful biomarkers for extending the drug-screening process into a wider range of cancer cell lines. We also highlight a potential biomarker of dual EZH2/EHMT2 inhibition in MDA-MB-231 cells, the gene SPINK1, which shows no effect following transfection of siRNA targeting EZH2, nor siRNA targeting EHMT2, but shows increased expression following dual transfection of siRNAs targeting EZH2 and EHMT2. SPINK1 expression is increased upon treatment with HKMT-I-005, HKMT-I-011 and HKMT-I-022 at multiple doses and time points. SPINK1, a potent protease inhibitor, is also known as pancreatic secretory trypsin inhibitor (PSTI), and mutations in SPINK1 have been associated with chronic pancreatitis .
Genome-wide expression analysis revealed that genes upregulated upon treatment with HKMTI-1-005 were more enriched for genes silenced by EZH2 than treatment with either the specific EHMT2 inhibitor UNC0638 or the specific EZH2 inhibitor GSK343. It was interesting to note that the EHMT2 inhibitor UNC0638 seemed to be as effective as the specific EZH2 inhibitor GSK343 in terms of specific upregulation of genes silenced by EZH2. This could in part be explained by the fact that EHMT2 has the capacity to methylate H3K27 [26, 27] and that reversal of epigenetic silencing of certain EZH2 targets is dependent on inhibition of EHMT2 . Alternatively, it could be due to differences in the kinetics of the inhibitors that act through different mechanisms, and the fact that genome-wide expression analysis was only carried out within a limited time window.
We also note that the effects observed on gene expression, chromatin marks and global levels of H3K27me3 and H3K9me3 occur within 24–72 h, while some previously reported EZH2 inhibitors only show pharmacodynamic effects at later time points [10, 12, 14–16]. There may be many reasons for these differences, including the mechanism of action of the dual inhibitors, as well as their effects on mRNA levels of EZH2 and the H3K27 demethylase JMJD3. However, it should be noted that the kinetics of effects on gene expression we observe with the dual inhibitors are similar to the kinetics of effects on gene expression we observe with double siRNA knockdown of EZH2 and EHMT2. Also, we observe less effect of GSK343 on H3K27me3 after a 48-h treatment of MDA-MB-231 cells compared to previous studies in ovarian cell line models after 3 days of treatment . This lower activity on H3K27me3 could be cell specific or due to a shorter treatment period. The wealth of cellular data accumulated for our hit compounds, HKMTI-1-005 in particular, argue for direct effects on cells at the target H3K27me and H3K9me modifications at doses of drug less than or equivalent to growth inhibitory doses. Such data includes the specific expression of EZH2 target genes, global histone methylation changes by Western analysis and local chromatin changes on responsive genes. We also note the increased sensitivity of the mutant EZH2 DB lymphoma cell line to HKMTI-1-005, in accordance with an EZH2 inhibitory effect. Such cellular biological effects are observed at doses of hit compounds less than the in vitro biochemical IC50 detected for EZH2. We would argue that the cellular activity is a consequence of dual HKMT activity and so extrapolating from single enzyme IC50 values is difficult. Furthermore, since the in vitro biochemical EZH2 activity assay conditions used the minimal number of proteins (EZH2:EED:SUZ12) and a simple peptide substrate, rather than the complex (and dynamic) in vivo target of chromatin, care should be taken in drawing quantitative comparisons with cell-based data.
The hit compounds reported herein represent starting points for the further optimisation of dual EZH2/EHMT2 inhibitors. Indeed, recent reports suggest it is possible to improve the in vivo profile of this compound class . While this scaffold has been extensively pursued for selective EHMT1/2 inhibition, further studies are needed to confirm whether it is possible to simultaneously increase potency against both EZH2 and EHMT1/2 and whether it is possible to engineer EHMT1/2 activity out of this scaffold to identify a selective substrate competitive EZH2 inhibitor. Nonetheless, it will continue to be important to ‘repurpose’ existing HKMT inhibitor chemotypes, in light of the low number of validated HKMT inhibitory chemotypes currently available .
Many cancers show aberrant silencing of gene expression and overexpression of histone methyltransferases, including EZH2 and EHMT1/2. We have shown that combined inhibition of EHMT1/2 and EZH2 increases growth inhibition in tumour cells over inhibition of only EHMT1/2 or EZH2 and results in re-expression of silenced genes. We report the first dual EZH2-EHMT1/2 substrate competitive inhibitors and show that they may have greater activity in tumour cells that overexpress wild type EZH2.
qRT-PCR measurements for cell-based screening
Following compound treatment of MDA-MB-231 for 48 h (in 6-well plates), media was removed and 1.5 ml of TRIzol (Invitrogen) was added directly to lyse cells and RNA was isolated according to the manufacturer’s instructions. Reverse transcription was done using the SuperScript III First-Strand Synthesis System (Invitrogen) according to the manufacturer’s instructions. Each measurement was done in triplicate, and the List of Primers can be found in Additional file 2: Table S1. For normalisation, we have used GAPDH and RNA pol II. Experiments were also done with the ‘Fast Sybr Green Cell-to-CT™-Kit’ according to the manufacturer’s instructions (Applied Biosystem). Fifteen thousand cells per 96 well were plated, and after 24 h, treated with compounds at various concentrations.
SiRNA experiments were carried out on the MDA-MB-231 cell line using Qiagen reagents, according to the manufacturer’s instructions. In brief, cells were seeded at a density of 1 × 105 cells/6 cm well and treated for 48 h with siRNAs given in Additional file 2: Table S2.
Chromatin immunoprecipitation (ChIP-PCR) assay
ChIP was accomplished using Dynabeads Protein A (Invitrogen) according to , except that following the Chelex-DNA purification, an additional purification with QIAquick PCR Purification Kit (Qiagen) was carried out; here, the ChIP-products were eluted in 50 μl and for subsequent qPCR measurements (as described above). The list of Primers can be found in Additional file 2: Table S3. Results were calculated as a fold increase of the no-antibody control and then normalised to GAPDH (active marks) and beta-globin (inactive marks).
Cell viability assay
For clonogenic assays, cells were treated with compound for 48 h, media were removed and replaced with fresh media for a further 12 days to measure colony formation. For cell viability assays, lymphoma cells were plated at 20,000 cells in 200 μl per well in U-bottom 96 well plates in RPMI medium +20 % FCS. Cells were re-suspended 48 h later, diluted 10-fold in PBS + propidium iodide (PI) and the concentration of PI negative cells was counted using an Attune flow cytometer with autosampler. Breast cancer cells were seeded at a density of 10,000 cells/well in a sterile 96 clear-well plate with 150 μl of DMEM (+10 % FCS and 2 mM L-Glutamine). Each compound treatment was performed in triplicate for 72 h at concentrations of 100 nM and 1, 5, 10 and 50 μM in 100 μl of full-medium. After 72 h, 20 μl of MTT solution (3 mg of MTT Formazan, Sigma/1 ml PBS) was added to the medium and incubated for 4 h at 37 °C in a CO2 incubator. The MTT-product was solubilised with 100 μl DMSO, and for 1 h, incubated in the dark at room temperature. The optical density was read at 570 nm with PHERAstar.
MDA-MB-231 cells seeded in 6 well plates at a cell density of 3 x 105 were treated with HKMTI-1-005 (1–7.5uM) for 48 h. Following lysis in Triton Extraction Buffer (TEB, PBS containing 0.5 % Triton ×100 (v/v), 1/1000 protease inhibitor), nuclei were re-suspended in 0.2 N HCL at a density of 4 x 107 nuclei per millilitre and incubated overnight at 4 °C to acid-extract the histones, before being centrifuged at 6500g for 10 min at 4 °C. Protein concentration was determined using the Bradford assay. H3K27me3, H3K9me3, H3K9me2, H3K9me and total H3 protein expression levels in the histone extract samples were determined using Western blot analysis using H3K27me3 (1:1000; Abcam), H3K9me3 (1:1000; Abcam), H3K9me (1:1000) and H3 (1:2000; Abcam) antibodies. After washing, the membrane was incubated with a horseradish peroxidase-labelled secondary antibody (1 h, room temperature). The membrane was incubated for 1 min with 5 mL of Pierce ECL Western blotting substrate (Thermo Scientific). Images were captured using Konica Minolta SRX101A Tabletop X-Ray film processor.
Gene expression microarrays
Agilent 80k two-colour microarrays were used to profile gene expression changes induced by treatment with drug compounds in MDA-MB-231 cells, both at 24 and 48 h. In the initial microarray experiment, three replicates were used for each drug, time combination, and in the validation study, four replicates were used. A separate untreated control sample was used for comparison with each replicate. Sample labelling, array hybridization and scanning were performed by Oxford Gene Technologies, according to manufacturer’s instructions. Feature extracted files were imported into GeneSpring (Agilent), and data was normalised to produce log2 ratios of treated/untreated for each replicate of each drug, time combination.
Normalised log2 gene expression ratios were analysed using LIMMA  to obtain empirical Bayes-moderated t statistics for differential expression across the replicates for each drug treatment. After multiple testing adjustments by the Benjamini-Hochberg method, p < 0.1 was used to denote significant differential expression in the initial microarray experiment and p < 0.05 in the validation experiment.
A list of EZH2 targets in MDA-MB-231 cells were taken from . Statistical significance of systematic upregulation or downregulation of these targets was evaluated using the ‘GeneSetTest’ method from the Bioconductor package limma. The same method was used to evaluate systematic up- or downregulation of pathways as annotated in ConsensusPathDB . Further analysis was performed using DAVID  for exploration of functional annotation enrichments.
Identification of a set of consensus EZH2-suppressed genes via meta-analysis
A meta-analysis of 18 microarray experiments were carried out as described in ‘Supplementary Methods’, resulting in the list of consensus EZH2 target genes given in Additional file 3: Table S4.
We would like to thank Alun Passey for the help with the experiments during review of the manuscript. We would like to acknowledge the Ovarian Cancer Action and Cancer Research UK for the funding (grant C21484/A6944, C536/A13086). IG acknowledges PhD studentship from the Imperial Cancer Research UK Centre. SS acknowledges the European Commission for a Marie Curie International Incoming Fellowship (Agreement No. 299857). NS was supported by a Royal Thai Government Scholarship and the EPSRC-funded Institute of Chemical Biology Doctoral Training Centre. JC acknowledges support from the ARC. The SGC is a registered charity (number 1097737) that receives funds from AbbVie, Boehringer Ingelheim, the Canada Foundation for Innovation, the Canadian Institutes for Health Research, Genome Canada through the Ontario Genomics Institute (OGI-055), GlaxoSmithKline, Janssen, Lilly Canada, the Novartis Research Foundation, the Ontario Ministry of Economic Development and Innovation, Pfizer, Takeda and the Wellcome Trust [092809/Z/10/Z].
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- 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.View ArticlePubMedGoogle Scholar
- Kleer CG, Cao Q, Varambally S, Shen R, Ota I, Tomlins SA, et al. EZH2 is a marker of aggressive breast cancer and promotes neoplastic transformation of breast epithelial cells. Proc Natl Acad Sci. 2003;100:11606–11.PubMed CentralView ArticlePubMedGoogle Scholar
- Bachmann IM, Halvorsen OJ, Collett K, Stefansson IM, Straume O, Haukaas SA, et al. EZH2 expression is associated with high proliferation rate and aggressive tumor subgroups in cutaneous melanoma and cancers of the endometrium, prostate, and breast. J Clin Oncol. 2006;24:268–73.View ArticlePubMedGoogle Scholar
- Matsukawa Y, Semba S, Kato H, Ito A, Yanagihara K, Yokozaki H. Expression of the enhancer of zeste homolog 2 is correlated with poor prognosis in human gastric cancer. Cancer Sci. 2006;97:484–91.View ArticlePubMedGoogle Scholar
- Varambally S, Dhanasekaran SM, Zhou M, Barrette TR, Kumar-Sinha C, Sanda MG, et al. The polycomb group protein EZH2 is involved in progression of prostate cancer. Nature. 2002;419:624–9.View ArticlePubMedGoogle Scholar
- Collett K, Eide GE, Arnes J, Stefansson IM, Eide J, Braaten A, et al. Expression of enhancer of zeste homologue 2 is significantly associated with increased tumor cell proliferation and is a marker of aggressive breast cancer. Clin Cancer Res. 2006;12:1168–74.View ArticlePubMedGoogle Scholar
- Yehiely F, Moyano JV, Evans JR, Nielsen TO, Cryns VL. Deconstructing the molecular portrait of basal-like breast cancer. Trends Mol Med. 2006;12:537–44.View ArticlePubMedGoogle Scholar
- Bracken AP, Pasini D, Capra M, Prosperini E, Colli E, Helin K. EZH2 is downstream of the pRB-E2F pathway, essential for proliferation and amplified in cancer. EMBO J. 2003;22:5323–35.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Cai Q, Godwin AK, Zhang R. Enhancer of zeste homolog 2 promotes the proliferation and invasion of epithelial ovarian cancer cells. Mol Cancer Res. 2010;8:1610–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Verma SK, Tian X, LaFrance LV, Duquenne C, Suarez DP, Newlander KA, et al. Identification of potent, selective, cell-active inhibitors of the histone lysine methyltransferase EZH2. ACS Med Chem Lett. 2012;3:1091–6.PubMed CentralView ArticlePubMedGoogle Scholar
- McCabe MT, Ott HM, Ganji G, Korenchuk S, Thompson C, Van Aller GS, et al. EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations. Nature. 2012;492:108–12.View ArticlePubMedGoogle Scholar
- Knutson SK, Wigle TJ, Warholic NM, Sneeringer CJ, Allain CJ, Klaus CR, et al. A selective inhibitor of EZH2 blocks H3K27 methylation and kills mutant lymphoma cells. Nat Chem Biol. 2012;8:890–6.PubMedGoogle Scholar
- Qi W, Chan H, Teng L, Li L, Chuai S, Zhang R, et al. Selective inhibition of Ezh2 by a small molecule inhibitor blocks tumor cells proliferation. Proc Natl Acad Sci. 2012;109:21360–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Konze KD, Ma A, Li F, Barsyte-Lovejoy D, Parton T, MacNevin CJ, et al. An orally bioavailable chemical probe of the lysine methyltransferases EZH2 and EZH1. ACS Chem Biol. 2013;8:1324–34.PubMed CentralView ArticlePubMedGoogle Scholar
- Knutson SK, Warholic NM, Wigle TJ, Klaus CR, Allain CJ, Raimondi A, et al. Durable tumor regression in genetically altered malignant rhabdoid tumors by inhibition of methyltransferase EZH2. Proc Natl Acad Sci. 2013;110:7922–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Garapaty-Rao S, Nasveschuk C, Gagnon A, Chan Eric Y, Sandy P, Busby J, et al. Identification of EZH2 and EZH1 small molecule inhibitors with selective impact on diffuse large B cell lymphoma cell growth. Chem Biol. 2013;20:1329–39.View ArticlePubMedGoogle Scholar
- Yoo CB, Jones PA. Epigenetic therapy of cancer: past, present and future. Nat Rev Drug Discov. 2006;5:37–50.View ArticlePubMedGoogle Scholar
- Cho HS, Kelly JD, Hayami S, Toyokawa G, Takawa M, Yoshimatsu M, et al. Enhanced expression of EHMT2 is involved in the proliferation of cancer cells through negative regulation of SIAH1. Neoplasia. 2011;13:676.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding J, Li T, Wang X, Zhao E, Choi J-H, Yang L, et al. The histone H3 methyltransferase G9A epigenetically activates the serine-glycine synthesis pathway to sustain cancer cell survival and proliferation. Cell Metab. 2013;18:896–907.View ArticlePubMedGoogle Scholar
- Liu F, Chen X, Allali-Hassani A, Quinn AM, Wasney GA, Dong A, et al. Discovery of a 2,4-diamino-7-aminoalkoxyquinazoline as a potent and selective inhibitor of histone lysine methyltransferase G9a. J Med Chem. 2009;52:7950–3.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu F, Chen X, Allali-Hassani A, Quinn AM, Wigle TJ, Wasney GA, et al. Protein lysine methyltransferase G9a inhibitors: design, synthesis, and structure activity relationships of 2,4-diamino-7-aminoalkoxy-quinazolines. J Med Chem. 2010;53:5844–57.PubMed CentralView ArticlePubMedGoogle Scholar
- Vedadi M, Barsyte-Lovejoy D, Liu F, Rival-Gervier S, Allali-Hassani A, Labrie V, et al. A chemical probe selectively inhibits G9a and GLP methyltransferase activity in cells. Nat Chem Biol. 2011;7:566–74.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang Y, Ganesh T, Horton JR, Spannhoff A, Liu J, Sun A, et al. Adding a lysine mimic in the design of potent inhibitors of histone lysine methyltransferases. J Mol Biol. 2010;400:1–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu F, Barsyte-Lovejoy D, Allali-Hassani A, He Y, Herold JM, Chen X, et al. Optimization of cellular activity of G9a inhibitors 7-aminoalkoxy-quinazolines. J Med Chem. 2011;54:6139–50.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu F, Barsyte-Lovejoy D, Li F, Xiong Y, Korboukh V, Huang X-P, et al. Discovery of an in vivo chemical probe of the lysine methyltransferases G9a and GLP. J Med Chem. 2013;56:8931–42.View ArticlePubMedGoogle Scholar
- Tachibana M, Sugimoto K, Fukushima T, Shinkai Y. SET domain-containing protein, G9a, is a novel lysine-preferring mammalian histone methyltransferase with hyperactivity and specific selectivity to lysines 9 and 27 of histone H3. J Biol Chem. 2001;276:25309–17.View ArticlePubMedGoogle Scholar
- Wu H, Chen X, Xiong J, Li Y, Li H, Ding X, et al. Histone methyltransferase G9a contributes to H3K27 methylation in vivo. Cell Res. 2011;21:365–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoo KH, Hennighausen L. EZH2 methyltransferase and H3K27 methylation in breast cancer. Int J Biol Sci. 2012;8:59–65.PubMed CentralView ArticlePubMedGoogle Scholar
- Mozzetta C, Pontis J, Fritsch L, Robin P, Portoso M, Proux C, et al. The histone H3 lysine 9 methyltransferases G9a and GLP regulate polycomb repressive complex 2-mediated gene silencing. Mol Cell. 2014;53:277–89.View ArticlePubMedGoogle Scholar
- Tan J, Yang X, Zhuang L, Jiang X, Chen W, Lee PL, et al. Pharmacologic disruption of Polycomb-repressive complex 2-mediated gene repression selectively induces apoptosis in cancer cells. Genes Dev. 2007;21:1050–63.PubMed CentralView ArticlePubMedGoogle Scholar
- Dillon S, Zhang X, Trievel R, Cheng X. The SET-domain protein superfamily: protein lysine methyltransferases. Genome Biol. 2005;6:227.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang Y, Zhang X, Horton JR, Upadhyay AK, Spannhoff A, Liu J, et al. Structural basis for G9a-like protein lysine methyltransferase inhibition by BIX-01294. Nat Struct Mol Biol. 2009;16:312–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Yun M, Wu J, Workman JL, Li B. Readers of histone modifications. Cell Res. 2011;21:564–78.PubMed CentralView ArticlePubMedGoogle Scholar
- Foletta V, White L, Larsen A, Léger B, Russell A. The role and regulation of MAFbx/atrogin-1 and MuRF1 in skeletal muscle atrophy. Pflugers Arch - Eur J Physiol. 2011;461:325–35.View ArticleGoogle Scholar
- Lee Shuet T, Li Z, Wu Z, Aau M, Guan P, Karuturi RKM, et al. Context-specific regulation of NF-κB target gene expression by EZH2 in breast cancers. Mol Cell. 2011;43:798–810.View ArticlePubMedGoogle Scholar
- Kamburov A, Pentchev K, Galicka H, Wierling C, Lehrach H, Herwig R. ConsensusPathDB: toward a more complete picture of cell biology. Nucleic Acids Res. 2011;39:D712–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57:289–300.Google Scholar
- Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, et al. COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 2011;39:D945–50.PubMed CentralView ArticlePubMedGoogle Scholar
- Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483:603–307.PubMed CentralView ArticlePubMedGoogle Scholar
- Srimongkolpithak N, Sundriyal S, Li F, Vedadi M, Fuchter MJ. Identification of 2,4-diamino-6,7-dimethoxyquinoline derivatives as G9a inhibitors. Med Chem Commun. 2014.Google Scholar
- Brown R, Curry E, Magnani L, Wilhelm-Benartzi CS, Borley J. Poised epigenetic states and acquired drug resistance in cancer. Nat Rev Cancer. 2014;14(11)):747–53. doi:10.1038/nrc3819.View ArticlePubMedGoogle Scholar
- Witt H, Luck W, Hennies HC, Claszen M, Kage A, Lasz U, et al. Mutations in the gene encoding the serine protease inhibitor, Kazal type 1 are associated with chronic pancreatitis. Nat Genet. 2000;25:213–6.View ArticlePubMedGoogle Scholar
- Amatangelo MD, Garipov A, Li H, Conejo-Garcia JR, Speicher DW, Zhang R. Three-dimensional culture sensitizes epithelial ovarian cancer cells to EZH2 methyltransferase inhibition. Cell Cycle. 2013;12:2113–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson JD, Denisenko O, Bomsztyk K. Protocol for the fast chromatin immunoprecipitation (ChIP) method. Nat Protocols. 2006;1:179–85.View ArticlePubMedGoogle Scholar
- Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3.Google Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols. 2008;4:44–57.View ArticleGoogle Scholar