- Open Access
PRC2 targeting is a therapeutic strategy for EZ score defined high-risk multiple myeloma patients and overcome resistance to IMiDs
Clinical Epigenetics volume 10, Article number: 121 (2018)
Multiple myeloma (MM) is a malignant plasma cell disease with a poor survival, characterized by the accumulation of myeloma cells (MMCs) within the bone marrow. Epigenetic modifications in MM are associated not only with cancer development and progression, but also with drug resistance.
We identified a significant upregulation of the polycomb repressive complex 2 (PRC2) core genes in MM cells in association with proliferation. We used EPZ-6438, a specific small molecule inhibitor of EZH2 methyltransferase activity, to evaluate its effects on MM cells phenotype and gene expression prolile.
PRC2 targeting results in growth inhibition due to cell cycle arrest and apoptosis together with polycomb, DNA methylation, TP53, and RB1 target genes induction. Resistance to EZH2 inhibitor is mediated by DNA methylation of PRC2 target genes. We also demonstrate a synergistic effect of EPZ-6438 and lenalidomide, a conventional drug used for MM treatment, activating B cell transcription factors and tumor suppressor gene expression in concert with MYC repression. We establish a gene expression-based EZ score allowing to identify poor prognosis patients that could benefit from EZH2 inhibitor treatment.
These data suggest that PRC2 targeting in association with IMiDs could have a therapeutic interest in MM patients characterized by high EZ score values, reactivating B cell transcription factors, and tumor suppressor genes.
Epigenetic events are key mechanisms in the regulation of cell fate and cell identity. DNA methylation, miRNA-associated gene repression, and histone modifications have been implicated in numerous diseases, including cancers and represent new therapeutic targets .
Gene expression regulation through polycomb-induced histone modifications is a well-studied mechanism. The polycomb repressive complex 2 (PRC2) contains three core subunits: EED (embryonic ectoderm development), SUZ12 (suppressor of zeste 12 homolog), and EZH2 (enhancer of zeste homolog 2). PRC2 represses gene transcription through tri-methylation of lysine 27 of histone 3 (H3K27me3) by its catalytic subunit EZH2.
EZH2 deregulation has been described in many cancer types, including hematological malignancies. Its overexpression or gain of function mutations lead to abnormal H3K27me3 accumulation, repressing tumor suppressor genes, such as cell cycle inhibitors, apoptotic activators, and senescence and differentiation factors .
Multiple myeloma (MM) is a neoplasia characterized by the accumulation of clonal plasma cells within the bone marrow. Recent advances in treatment have led to an overall survival of intensively treated patients of 6–7 years and an event-free survival of 3–4 years . However, patients invariably relapse after multiple lines of treatment, with shortened intervals between relapses, and finally become resistant to all treatments, resulting in loss of clinical control over the disease. MM is a genetically and clinically heterogeneous disease. Genome sequencing studies have revealed considerable heterogeneity and genomic instability, a complex mutational landscape and a branching pattern of clonal evolution [3, 4]. Epigenetic marks such as DNA methylation or histone posttranslational modifications are also involved in MM pathophysiology and drug resistance [5, 6].
Global gene expression profiling indicated that, while EZH2 is upregulated, its target genes are downregulated in myeloma cells compared with normal plasma cells . In human MM cell lines (HMCL), EZH2 expression has been correlated with increased proliferation and an independence on growth factors . Inhibition of EZH2 expression and activity is associated with HMCL growth inhibition [9, 10] and decreased tumor load in a mouse model of MM [7, 11]. One study shows that this effect is related to epithelial tumor suppressor gene upregulation . However, the use of specific EZH2 inhibitors demonstrated that MM proliferation inhibition is time dependent and cell line specific, indicating that EZH2 does not play a universal and monotonous role in promoting MM . Furthermore, the first genome-wide profiling of H3K27me3 and H3K4me3 in MM patient samples was recently published, showing a unique epigenetic profile of primary MM cells compared to normal bone marrow plasma cells . EZH2 inhibition was associated with upregulation of microRNAs with potential tumor suppressor functions . More recently, EZH2 overexpression was reported to be associated with poor outcome and dysregulation of proliferation . These data underscore an oncogenic role of EZH2 in MM. EZH2 inhibitors are currently in phase 2 clinical development in relapsed or refractory non-Hodgkin lymphoma (NHL) and biomarkers are needed for patient selection since neither EZH2 mutations nor H3K27me3 levels are sufficient to predict NHL cell response to EZH2 inhibitors [1, 14].
Here, we identified that PRC2 core genes are overexpressed in MM cells in association with proliferation activation. Treatment of MM cells with EPZ-6438, a specific small molecule inhibitor of EZH2 methyltransferase activity, results in growth inhibition due to cell cycle arrest and apoptosis. Resistance to EZH2 inhibitor is mediated by DNA methylation of PRC2 target genes. We also observed a synergy between EPZ-6438 and lenalidomide, a conventional drug used for MM treatment. More interestingly, pretreatment of myeloma cells with EPZ-6438 significantly re-sensitizes drug-resistant MM cells to lenalidomide. EPZ-6438/lenalidomide combination induced a significant transcriptional reprogramming of MM cells targeting major B and plasma cell transcription factors in association with MYC repression. RNA sequencing combined with H3K27me3 ChIP analyses allowed us to build an EZ GEP-based score that is able to predict HMCL and primary MM cell sensitivity to EZH2 inhibitors.
Human myeloma cell lines (HMCLs)
XG human myeloma cell lines were obtained as previously described . JJN3 was kindly provided by Dr. Van Riet (Brussels, Belgium), JIM3 by Dr. MacLennan (Birmingham, UK), and MM1S by Dr. S. Rosen (Chicago, USA). AMO-1, LP1, L363, U266, OPM2, and SKMM2 were purchased from DSMZ (Braunsweig, Germany) and RPMI8226 from ATTC (Rockville, MD, USA). All HMCLs derived in our laboratory were cultured in the presence of recombinant IL-6. HMCLs were authenticated according to their short tandem repeat profiling and their gene expression profiling using Affymetrix U133 plus 2.0 microarrays deposited in the ArrayExpress public database under accession numbers E-TABM-937 and E-TABM-1088 .
Primary multiple myeloma cells
Bone marrow samples were collected after patients’ written informed consent in accordance with the Declaration of Helsinki and institutional research board approval from Heidelberg and Montpellier University Hospital. Bone marrows were collected from 206 patients treated with high-dose Melphalan (HDM) and autologous stem cell transplantation (ASCT), and this cohort is termed “Heidelberg-Montpellier” (HM) cohort . Patients’ MMCs were purified using anti-CD138 MACS microbeads (Miltenyi Biotec, Bergisch Gladbach, Germany) and their gene expression profile (GEP) obtained using Affymetrix U133 plus 2.0 microarrays as described . The CEL files and MAS5 files are available in the ArrayExpress public database (E-MTAB-372). The structural chromosomal aberrations, as well as numerical aberrations, were assayed by fluorescence in situ hybridization (iFISH). We also used publicly available Affymetrix GEP (Gene Expression Omnibus, accession number GSE2658) of a cohort of 345 purified MMC from previously untreated patients from the University of Arkansas for Medical Sciences (UAMS, Little Rock, AR, USA), termed in the following UAMS-TT2 cohort. These patients were treated with total therapy 2 including HDM and ASCT . We also used Affymetrix data from total therapy 3 cohort (UAMS-TT3; n = 158; E-TABM-1138)  of 188 relapsed MM patients subsequently treated with bortezomib (GSE9782) from the study by Mulligan et al. .
EZH2 inhibition in primary MM cells
BM of patients presenting with previously untreated MM (n = 7) at the university hospital of Montpellier was obtained after patients’ written informed consent in accordance with the Declaration of Helsinki and agreement of the Montpellier University Hospital Centre for Biological Resources (DC-2008-417). Mononuclear cells were treated with or without EPZ-6438 (370 nM or 1 μM) and/or lenalidomide (2 μM), and MMC cytotoxicity was evaluated using anti-CD138-phycoerythrin monoclonal antibody (Immunotech, Marseille, France) as described .
Cell growth assay
HMCLs were cultured for 15 days in RPMI 1640 medium, 10% FCS, and 2 ng/ml IL-6 (control medium) in the presence of EPZ-6438 (Selleckchem, Houston, TX, USA) and or decitabine. Cell concentration and viability were assessed using trypan blue dye exclusion test. The number of metabolic-active cells was also determined using intracellular ATP quantitation. For drug combination assay, HMCLs were cultured for 4 days in 96-well flat-bottom microtiter plates in RPMI 1640 medium, 10% FCS, and 2 ng/ml IL-6 (control medium) in the presence of lenalidomide. Cell growth was evaluated by quantifying intracellular ATP amount with a Cell Titer Glo Luminescent Assay (Promega, Madison, WI, USA) using a Centro LB 960 luminometer (Berthold Technologies, Bad Wildbad, Germany).
Global H3K27me3 and IKZF1 immunofluorescence
After deposition on slides using a Cytospin centrifuge, cells were fixed with 4% PFA, permeabilized with 0.5% Triton in PBS and saturated with 5% bovine milk in PBS. The rabbit anti-H3K27me3 (Active Motif, Rixensart, Belgium, #39156) and anti-IKZF1 (Santa Cruz Biotechnology, Heidelberg, Germany, H-100 sc-3039) antibodies were diluted 1/500 and 1/250 respectively in 5% bovine milk in PBS and deposited on cytospins for 60 min at room temperature. Slides were washed twice, and anti-rabbit alexa 555-conjugated antibodies (diluted 1/500 in 5% bovine milk in PBS) were added for 60 min at room temperature. Slides were washed and mounted with Vectashield and 1% DAPI. Images and fluorescence were captured with a ZEISS Axio Imager Z2 microscope (× 63 objective) (Oberkochen, Germany) and analyzed with Omero (omero.mri.cnrs.fr) server and ImageJ software.
HMCLs were cultured for 4 days without or with 1 μM of EPZ6438. RNA samples were collected as previously described. The RNA sequencing (RNA-seq) library preparation was done with 150 ng of input RNA using the Illumina TruSeq Stranded mRNA Library Prep Kit. Paired-end RNA-seq were performed with Illumina NextSeq sequencing instrument (Helixio, Clermont-Ferrand, France). RNA-seq read pairs were mapped to the reference human GRCh37 genome using the STAR aligner . All statistical analyses were performed with the statistics software R (version 3.2.3; available from https://www.r-project.org) and R packages developed by BioConductor project (available from https://www.bioconductor.org/) . The expression level of each gene was summarized and normalized using DESeq2 R/Bioconductor package . Differential expression analysis was performed using DESeq2 pipeline. P values were adjusted to control the global FDR across all comparisons with the default option of the DESeq2 package. Genes were considered differentially expressed if they had an adjusted P value of 0.05 and a fold change of 1.5. Pathway enrichment analyses were performed using online curated gene set collection on the Gene Set Enrichment Analysis software (http://software.broadinstitute.org/gsea/msigdb/index.jsp) [24, 25].
Gene expression profiling and statistical analyses
Gene expression data were normalized with the MAS5 algorithm and analyses processed with GenomicScape (http://www.genomicscape.com)  the R.2.10.1 and Bioconductor version 2.5 programs . Gene set expression analysis (GSEA) was used to identify genes and pathways differentially expressed between populations. Univariate and multivariate analysis of genes prognostic for patients’ survival was performed using the Cox proportional hazard model. Difference in overall survival between groups of patients was assayed with a log-rank test and survival curves plotted using the Kaplan-Meier method (Maxstat R package) . The EZ score was built using our previously published methodologies . EZ score is the sum of the Cox beta-coefficients of each of the 15 EPZ6438-deregulated genes with a prognostic value, weighted by ± 1 if the patient MMC signal for a given gene is above or below the Maxstat reference value of this gene (Table 1) .
Cell cycle analysis
HMCLs were cultured in 24-well, flat-bottomed microtiter plates at 105 cells per well in RPMI1640–10% FCS or X-VIVO 20 culture medium with or without IL-6 (3 ng/mL) and EPZ-6438 (Selleckchem). The cell cycle was assessed using DAPI staining (Sigma-Aldrich, Saint-Louis, MO, USA) and cells in the S phase using incubation with bromodeoxyuridine (BrdU) for 1 h and labeling with an anti-BrdU antibody (APC BrdU flow kit, BD Biosciences, San Jose, CA, USA) according to the manufacturer’s instructions.
Flow cytometry analysis
Cells were fixed for 10 min with Cytofix/Cytoperm (BD Biosciences, San Jose, CA, USA) at 4 °C. The overall expression of MYC, IKZF1, IRF4, and H3K27me3 was evaluated by incubating 105 cells with 5 μLof an alexa 647-conjugated mouse anti-H3K27me3 antibody (Cell Signaling, Danvers, MA, USA, 12158S), alexa 647-conjugated mouse anti-EZH2 antibody (BD Biosciences, 563491) PE-conjugated mouse anti-IKZF1 (BD Biosciences, 564476), rat anti-IRF4 (Biolegend, 646403), rabbit anti-MYC (Cell Signaling, #12189), or anti-Ki67 antibodies in phosphate-buffered saline (PBS) containing 2% FBS at 4 °C for 20 min.
For primary samples, cells were double stained with APC or PE-conjugated anti-CD138 (Beckman-Coulter, Brea, CA, USA). Flow cytometry analysis was done on a Fortessa flow cytometer (BD, Mountain View, CA, USA).
Study of apoptosis
HMCLs were cultured in 24-well, flat-bottomed microtiter plates at 105 cells per well in RPMI1640–10% FCS or X-VIVO 20 culture medium with or without IL-6 (3 ng/mL), EPZ-6438 (Selleckchem), and QVD. After 8 days of culture, cells were washed twice in PBS and apoptosis was assayed with PE-conjugated Annexin V labeling (BD Biosciences) using a Fortessa flow cytometer (BD).
H3K27me3-associated genes in MM patients
The list of H3K27me3 associated genes in MM patients was recovered from GEO GSE53215 and used for the development of the EZ score.
RNA was converted to cDNA using the Qiagen QuantiTect Reverse Transcription Kit (Qiagen, Hilden, Germany). The assays-on-demand primers and probes and the TaqMan Universal Master Mix were used according to the user’s manual (Biosystems, Courtaboeuf, France). The measurement of gene expression was performed using the Roche LC480 Sequence Detection System. For each primer, serial dilutions of a standard cDNA were amplified to create a standard curve, and values of unknown samples were estimated relative to this standard curve in order to assess PCR efficiency. Ct values were obtained for GAPDH and the respective genes of interest during log phase of the cycle. Gene expression was normalized to that of b2M (dCt = Ct gene of interest–Ct b2M) and compared with the values obtained for a known positive control using the following formula: 100/2ddCt where ddCt = dCt unknown–dCt positive control.
Cells were lysed in RIPA buffer (Cell Signaling Technology, Beverly, MA, USA) supplemented with 1 mM phenylmethylsulfonyl fluoride immediately before use. Lysates were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (10% gels) and transferred to nitrocellulose membranes using an iBlot® Gel Transfer Device (InVitrogen). Non-specific membrane sites were blocked by incubation at room temperature in 140 mM NaCl, 3 mM KCl, 25 mM Tris-HCl (pH 7.4), 0.1% Tween 20 (tris-buffered saline Tween-20), 5% non-fat milk for 2 h, and then immunoblotted with rabbit polyclonal antibodies against Ikaros (Santa Cruz Biotechnology, Dallas, TX, USA), cMyc (Cell Signaling Technology), or IRF4 (Santa Cruz Biotechnology, Dallas, TX, USA). As a control for protein loading, a mouse monoclonal anti-β-actin antibody (Sigma-Aldrich) was used. The primary antibodies were visualized with peroxidase-conjugated goat anti-rabbit (Sigma-Aldrich) or goat anti-mouse (Jackson ImmunoResearch, West Grove, PA, USA) antibodies and an enhanced chemiluminescence detection system. Western blots were quantified by densitometry using the NIH ImageJ software (National Institutes of Health, Bethesda, MD, USA), and protein levels were normalized according to those of β-actin.
Chromatin immunoprecipitation followed by sequencing
Cells were cross-linked in formaldehyde at a final concentration of 1% for 8 min. All experiments reagents were included in the Auto iDeal ChIP-seq kit for Histones (Diagenode, Liege, Belgium). Sonication was performed using a Bioruptor Plus sonication devise (Diagenode, Liege, Belgium) under the optimal conditions to shear cross-linked DNA to fragments if 100–300 base pairs in length. ChIP was conducted with the IPStar Compact Automated System (Diagenode, Liege, Belgium). ChIP were performed starting from 1 million cells per IP, using the indirect method in 200 μL final volume. Crossed-linked DNA was incubated 13 h with the antibody (H3K27me3 catalog number C1540195 or H3K4me3 catalog number C15410003) and 3 h with the beads. After 5 min washes, eluates were recovered and reverse cross-linked for 4 h at 65 °C. Samples were treated for 1 h with RNAse at 37 °C, prior to DNA purification with the Auto IPure kit v2 (Diagenode, C03010010). Libraries were performed using NEBNext Ultra Library Prep Kit for Illumina (New England Biolabs). Sequencing was performed with Illumina NextSeq500 technology (Helixio, Clermont-Ferrand, France) using the following parameters: single-read, 50 bp, 40 million reads.
DNA methylation analysis
Methylation analysis was performed using the Illumina Infinium HumanMethylation450 BeadChip array (HM450K, Illumina Inc.). The microarray raw intensities were preprocessed using the R/Bioconductor package minfi .The methylated CpGs overlapping with genes promoter region were extracted using the R package GenomicFeatures  implemented in Bioconductor. A CpG is defined as highly methylated when the methylation beta level is above 0.8.
PRC2 complex is overexpressed in malignant plasma cells in association with cell cycle deregulation
Using Affymetrix microarrays, we analyzed the expression of PRC2 core genes EZH2, SUZ12, and EED in normal bone marrow plasma cells (BMPCs, n = 5), primary myeloma cells from patients (MMCs, n = 206), and human myeloma cell lines (HMCLs, n = 26). PRC2 core genes are significantly overexpressed in MM cells (Fig. 1a) and EZH2 expression is significantly correlated with SUZ12 and EED expression (Additional file 1: Figure S1). At the opposite, PRC1 core genes were significantly downregulated in MM cells compared to normal plasma cells underlining polycomb complex deregulation in MM (Additional file 1: Figure S2).
Investigating the prognostic value of EZH2, SUZ12, and EED expression in independent cohorts of previously untreated MM patients, using Maxstat R algorithm , only EZH2 expression was found to be associated with MM patient’s outcome as recently reported  (Additional file 1: Figure S3). EZH2 is significantly overexpressed in patients with del17p and 1q21 gain (Additional file 1: Figure S4A and B). EZH2 expression is significantly correlated with MMC plasma cell labeling index (PCLI) in a cohort of 101 newly diagnosed patients (P < 0.005; Fig. 1b) underlining a link between PRC2 expression and deregulation of cell cycle in MM cells. Furthermore, GSEA analysis of patients with high EZH2 expression identified a significant enrichment for genes involved in cell cycle, upregulated in proliferating plasmablasts compared to mature BMPCs, and EZH2 targets (P < 0.0001) (Fig. 1c).
Altogether, these data underline that PRC2 complex is overexpressed in MMCs in association with a proliferative plasmablastic gene signature and a poor prognosis.
To investigate the therapeutic interest of PRC2 deregulation to target MM cells, XG1, XG12, XG19, LP1, XG25, and XG7 HMCLs were treated with clinically relevant doses of EZH2 inhibitor EPZ-6438 (370 nM and 1 μM) [14, 32]. EPZ-6438 treatment induced a significant decrease of global H3K27me3 in all the HMCLs tested (P < 0.01) (Additional file 1: Figures S5 and S6) and inhibited MM cell growth together with proliferation inhibition and apoptosis induction in 3 out of the 6 HMCLs tested (Fig. 2a and Additional file 1: Figure S7). The inhibitory effect appeared at day 6, suggesting that it is mediated by epigenetic reprogramming (Fig. 2a). EZH2 inhibitor induced apoptosis was partially rescued by the QVD pan-caspase inhibitor, suggesting a caspase-dependent mechanism (Additional file 1: Figure S7). LP1 and XG7 were more resistant to EZH2 inhibitor whereas XG25 HMCL was completely resistant. Primary MM cells co-cultured with their bone marrow microenvironment and recombinant IL-6 were also treated with EPZ-6438 as previously described . EZH2 targeting significantly reduced the median number of viable myeloma cells by 35% (P = 0.004) in 9/17 patients whereas MM cells of 8 patients were not significantly affected by the EZH2 inhibitor (Fig. 2b and Additional file 1: Figure S8). As described in HMCLs, EPZ-6438 induced a significant global H3K27me3 decrease in all the patients (Fig. 2c). The effect of EZH2 inhibitor was not correlated with EZH2 expression, H3K27me3 and levels in light of the 6 HMCLs and 17 primary MM samples tested. Moreover, UTX/JMJD3 demethylases mutation status did not seem to affect EPZ-6438 efficiency in the tested samples (Additional file 1: Figure S9 and Additional file 2: Table S1).
We therefore conclude that MM cell growth could be affected by PRC2 targeted therapy which presents a therapeutic interest in a subgroup of MM patients.
Genes deregulated by EZH2 targeted inhibition in myeloma cells
To better understand the molecular mechanisms associated with PRC2 deregulation in MM, six HMCLs were treated with 1 μM of EPZ-6438 for 4 days, and GEP were analyzed using RNA sequencing. Two hundred sixty-three genes are significantly upregulated in EPZ-6438-treated MMC compared to untreated cells (fold change ≥ 2, 1000 permutations, FDR < 0.05, Additional file 3: Table S2). Notably, no gene was significantly downregulated after EPZ-6438 treatment. EPZ-6438-regulated genes are significantly enriched in polycomb target genes, genes enriched in H3K27me3 histone mark, genes described to be associated with DNA methylation in MM, TP53, and RB1 target genes (Fig. 3a and Additional file 4: Table S3). Interestingly, a set of these genes is overexpressed in non-cycling mature normal plasma cells compared to proliferating plasmablasts (Fig. 3a and Additional file 4: Table S3). Among these 263 genes deregulated after EPZ-6438 treatment, 174 were also associated with H3K27me3 mark in XG7 HMCL including 160 bivalent genes characterized by H3K4me3 and H3K27me3 histone marks (Fig. 3b and Additional file 5: Table S4).
Resistance to EZH2 inhibitor is mediated by DNA methylation of PRC2 target genes
Since PRC2 target genes were associated with a significant enrichment of genes presenting DNA methylation in MM, we analyzed the methylation status at the promoter region of EZH2 inhibitor target genes using 450 k microarrays. Interestingly, 155 out the 263 EZH2 inhibitor target genes were associated with significant CpG methylation in their promoter (Fig. 3c). Furthermore, 111 of these 155 genes were also significantly upregulated after HMCL treatment with decitabine (ratio > 1.5) (Fig. 3d). A significant difference in CpG methylation status of EZH2 inhibitor target genes were identified between sensitive and resistant HMCLs (P < 0.01) (Fig. 3c). Interestingly, combination of a sublethal dose of decitabine (100 nM, IC10) with EPZ-6438 (1 μM) allowed to sensitize XG7 and XG25 EPZ-6438-resistant HMCLs. After 8 and 11 days of treatment, we observed a significant cell growth inhibition in both cell lines (67.4% and 62.4% respectively, P < 0.01 in XG7, and 72.2% and 85.2% respectively in XG25), compared to decitabine or EPZ-6438 used alone (P < 0.05) (Fig. 3e).
These data demonstrate that PRC2 target genes could be associated with DNA methylation underlying an overlap between epigenetic silencing mechanisms on potent key MM tumor suppressor genes.
PRC2 inhibition sensitizes myeloma cells to lenalidomide
We investigated whether EZH2 inhibition could enhance the anti-myeloma activity of melphalan, bortezomib, and lenalidomide treatment. EPZ-6438 pretreatment does not significantly enhance the effect of melphalan and bortezomib treatment on MMCs in vitro (data not shown). However, EPZ-6438 pretreatment significantly sensitized HMCLs to lenalidomide treatment. Treatment for 4 days with 1 μM of EPZ-6438 prior to lenalidomide treatment significantly enhanced toxicity on LP1 and XG19 HMCLs (P < 0.05) (Fig. 4a). Furthermore, EZH2 inhibition was able to overcome lenalidomide resistance in XG7 and XG25 cell lines (Fig. 4a). These data strongly suggest that combination of EZH2 inhibitor and IMiDs could be of therapeutic interest in MM.
These results were validated using primary MMCs from patients [6, 28]. The median percentage of viable MMCs was reduced of 43.4%, 26.9%, and 75.1% when cells were treated with lenalidomide, EPZ-6438, or the combination EPZ-6438/lenalidomide respectively (P = 0.003; n = 5) (Fig. 4b). Global H3K27me3 levels showed a tendency toward decrease under combination treatment, although the difference was not significant.
PRC2 targeting combined with Lenalidomide induced MM cell transcriptional reprogramming
Using RNA sequencing GEP data, we then compared MM cell lines treated by EPZ-6438, lenalidomide, or EPZ-6438/lenalidomide combination. Thirty-one genes were commonly upregulated by the EPZ-6438 and EPZ-6438/lenalidomide combination (Fig. 4c). Sixty-seven percent of the genes deregulated by lenalidomide in MM are also deregulated by the EPZ-6438/lenalidomide combination (Fig. 4c and Additional file 6: Table S6). Furthermore, EPZ-6438/lenalidomide induced a significantly higher deregulation of EPZ-6438 or lenalidomide target genes including PARP9, RGS1, DKK1, SEZ6L2, CAV2, ANTXR1, EMP1, and ANXA1 (ratio > 1.5, P < 0.05; Fig. 4d and Additional file 6: Table S6 and Additional file 7: Table S7).
Seventy-six percent of EPZ-6438/lenalidomide deregulated genes were not affected by EPZ-6438 or lenalidomide alone, suggesting that EZH2 inhibitor and IMiDs combination could modulate the expression of a specific set of genes. Among the 231 genes significantly upregulated uniquely after EPZ-6438/lenalidomide treatment, we found a significant enrichment of PRC2 and RB1 target genes, of genes downregulated by MYC, and genes silenced by DNA methylation (Fig. 4e and Additional file 8: Table S8). GSEA analysis of the 82 specifically downregulated genes by EPZ-6438/lenalidomide revealed a significant enrichment of MYC target genes in association with genes involved in proliferation and replicative stress response (Fig. 4f and Additional file 8: Table S8). Altogether, these data indicate that EPZ-6438/lenalidomide induced a significant upregulation of PRC2, RB1, and DNA methylation target genes in association with a significant downregulation of MYC target genes and cell proliferation gene program.
IMiDs promote Ikaros and Aiolos transcription factors binding to the E3 ubiquitin ligase cereblon (CRBN) leading to their ubiquitination and proteasomal degradation . Ikaros and Aiolos degradation is associated with a downregulation of interferon regulatory factor 4 (IRF4) and MYC which, in turn, reduces MMCs survival . Investigating the effect of EPZ-6438/lenalidomide combination in XG7 and LP1 MMCs, we found that Ikaros (IKZF1), IRF4, and MYC protein levels were significantly decreased by the combination treatment (65.5%, 63.9%, and 14.8% respectively) compared with lenalidomide (51.5%, 43% and 2.2%) or EPZ-6438 (45.2%, 38.7%, and 6.2%) alone (Fig. 5a). These data were validated by western blot and immunofluorescence (Additional file 1: Figures S10 and S11). Furthermore, EPZ-6438/lenalidomide combination strongly upregulated PAX5, BACH2, and BCL6 B cell transcription factors in association with downregulation of PRDM1 and IRF4 key plasma cell transcription factor (Fig. 5b). These data were validated by q-RT-PCR (Additional file 1: Figure S12). Interestingly, ChIP-seq analysis showed the bivalence of PAX5 gene. PAX5 expression could be directly regulated by EZH2 through H3K27me3 (Additional file 1: Figure S13). EPZ-6438/lenalidomide mediated transcriptional deregulation resulted in significant XG7 HMCL proliferation inhibition and quiescence induction characterized by a significant increase in the percentage of G0/G1 cells and ki67 negative cells together with E2F1 downregulation and CDKN1A induction (Fig. 5c–e).
EZ-GEP-based score allows to predict sensitivity of MMCs to EZH2 inhibitor treatment
Since the sensitivity of MMC to EZH2 inhibition is heterogeneous, we searched to define a biomarker allowing the identification of MM patients that could benefit from EZH2 inhibitor treatment. We built a GEP-based score (EZ score), using the prognostic information of 15 genes deregulated by EPZ-6438 and associated with H3K27me3 and a prognostic value in newly diagnosed MM patients (Fig. 6a and Additional file 9: TableS5). The EZ score is defined by the sum of the beta coefficients of the Cox model for each prognostic gene, weighted by − 1 according to the patient MMC signal above or below the probe set Maxstat value as previously described . EZ score levels in normal, premalignant, or malignant plasma cells are displayed in Fig. 6b. MMCs of patients had a significantly higher EZ score than normal BMPCs or plasma cells from MGUS patients (P < 0.01). HMCLs present an even higher EZ score compared with primary MMCs (P < 0.001; Fig. 6b), demonstrating an association between EZ score and the progression of the disease.
Using patient’s HM and UAMS-TT2 cohorts, the EZ score had prognostic value when used as a continuous variable or by splitting patients into two groups using the Maxstat R function . The EZ score split patients in a high-risk group (EZ score > 0.686) and a low-risk group (EZ score < 0.686) in the HM and UAMS-TT2 cohorts (P < 0.0001) (Fig. 6c).
We then sought to identify whether the EZ score could predict for the sensitivity of HMCLs to EPZ-6438. Using our large cohort of HMCLs , we analyzed the response of five HMCLs with the highest EZ score and five HMCLs with the lowest EZ score to EZH2 inhibitor treatment. The five HMCLs with the highest EZ score exhibit a significant 20-fold median higher sensitivity to EPZ-6438 compared with the five HMCLs with low EZ score (P = 0.04) (Fig. 6d).
To determine whether the EZ score could predict the sensitivity of primary MMCs to EPZ-6438, we analyzed the correlation between its toxicity on MMCs and the EZ score value in a panel of 14 patients. Primary MMCs were cultured together with their bone marrow environment, recombinant IL-6, and 1 μM of EPZ-6438 for 8 days. As identified in HMCLs, a significant correlation between EZ score and EZH2 inhibitor activity on primary MMCs was observed (r = − 0.68; P = 0.005) (Fig. 6e). A high EZ score value is associated with a higher toxicity of EZH2 inhibitor in MMCs. The EZ score allows the identification of a subgroup of MM patients with a poor outcome that could benefit from EZH2 inhibitor treatment.
Our results underline a significant deregulation of PRC1 and PRC2 complex in MM and extend recent studies which pointed out a role of EZH2 in MM biology [7,8,9,10,11, 13]. We demonstrated that PRC2 deregulation is associated with cell cycle deregulation and a proliferative plasmablastic gene expression signature. Our data suggest that PRC2 deregulation could support MM physiopathology, dissemination, and progression. EZH2 inhibitors have been reported to induce expression of miRNA and decrease expression of oncogenes [12, 35]. PRC2 targeting, using EZH2 inhibitor, represents a potent therapeutic strategy in a subgroup of MM patients resulting in a significant induction of polycomb target genes, genes associated with DNA methylation, TP53, and RB1 target genes. The analysis of genes deregulated by EED inhibitor or by knock down of other PRC2 component will be of interest to define the PRC2 specific target genes and the EZH2i unspecific targets. As reported in diffuse large B cell lymphoma , we found that H3K27me3 global levels are potently reduced after EPZ-6438 treatment in all HMCLs and primary MM samples tested. However, no link with EZH2 inhibitor sensitivity was identified. Even if UTX loss was recently reported to sensitize MM cell lines to EZH2 inhibition , neither EZH2 gene expression nor UTX mutation status was predictive of MM cell response to PRC2 targeting in our collection of MM cell lines representative of molecular heterogeneity or in the primary samples tested. However, analyses of a higher number of samples will be important to provide more significant conclusions. However, we identified a significant overlap between H3K27me3 and DNA methylation of EPZ-6438 target genes in MM cells resistant to EZH2 inhibitor. These two epigenetic repression systems are mechanistically linked. Indeed, EZH2, as part of PRC2 complex, is required for its target gene promoter methylation . These data underline that PRC2 target genes could comprise key MM tumor suppressor genes silenced by different epigenetic silencing. Interestingly, sublethal doses of DNMTi allow to sensitize EPZ-6438-resistant MM cell lines to EZH2 inhibitor. A major stake remains the identification of biomarkers that can quickly identify the subset patients who could benefit from EZH2 inhibitor treatment.
We developed an EZ score based on 15 EPZ-6438 target genes associated with H3K27me3 and prognostic value in primary MMCs. Our results demonstrate that the EZ score allows the identification of MM patients with a high EZ score value, an adverse prognosis and who could benefit from treatment with EZH2 inhibitors. The expression of 12 out of the 15 genes building the EZ score is associated with favorable outcome in MM. However, their biological functions in MM remain to be characterized.
Several of those genes are known tumor suppressors that are hypermethylated in solid cancers (CACNA1G, IRF6, ITGA9, and PPP2R2C) and used as biomarkers related to adverse outcome [38,39,40,41,42,43,44,45,46,47,48,49,50,51]. In MM, SORL1, a member of the low-density lipoprotein receptor family, is significantly hypermethylated at relapse compared to diagnosis in MM [52, 53]. SLC5A1 encodes a member of the sodium-dependent glucose transporter (SGLT) family. SLC5A1 is a good prognosis biomarker in cervical tumors. SLC5A1 activation through MAP17 increases ROS production .
FBLIM1 is a widely expressed protein presenting a role in cellular shape modulation, motility, and differentiation . Interestingly, FBLIM1 expression is reduced in breast cancer and has also been shown to sensitize glioma cells to cisplatin-induced apoptosis [56, 57]. Members of beta-arrestin such as ARRB1 participate in the sensitization/desensitization of G-protein-coupled receptors. The loss of ARRB1 is related to a poor outcome in NSCLC patients . ARRB1 can regulate apoptosis and DNA repair, and its overexpression induced DNA damage in NSCLC cells . THY1 (CD90) is a cell surface glycoprotein involved in cell adhesion and cell communication. In nasopharyngal carcinoma, THY1 is poorly expressed due to its promoter hypermethylation. THY1 overexpression in this cell type induces a G0/G1 cell cycle arrest . These data emphasize that EZH2 inhibitor induces transcriptional activation of potential MM tumor suppressor genes.
Given that monotherapy remains inefficient in MM, a particular association of drugs is usually used to treat MM disease. Unfortunately, despite rising advances in drug development and patient care, MM patients often relapse [61, 62]. By combining EPZ-6438 with MM conventional drugs, we uncovered its synergistic activity with lenalidomide. Furthermore, pre-treatment of MMCs with EZH2 inhibitor was able to overcome lenalidomide resistance (Fig. 4). Transcriptome analysis of cells treated with EPZ-6438/lenalidomide showed a combination-specific deregulated set of genes enriched in GSEA signature related to MYC targeting and genes involved in proliferation and replicative stress response. Furthermore, PRC2 inhibition and most significantly EPZ-6438/lenalidomide combination upregulated PAX5, BACH2, and BCL6 B cell transcription factors. Recently, EZH2i have been shown to induce BCL6 expression in MM cell lines in association with cell death . PAX5, BCL6, and BACH2 repress PRDM1 plasma cell transcription factor. PRDM1 is known to repress BCL6 and PAX5 [63,64,65,66], but epigenetic mechanisms may play a major role in B cell transcription factor silencing during plasma cell differentiation as underlined by the bivalent domain of PAX5, including H3K4me3 and H3K27me3 marks, identified in MM cells (Additional file 1: Figure S13). Accordingly, EPZ-6438/lenalidomide treatment resulted in downregulation of PRDM1 and IRF4 expression. By binding to cereblon (CRBN), the IMiD enhances Ikaros and Aiolos proteosomal degradation, thus leading to IRF4 and MYC downregulation . These two factors have been shown to be critical for MMCs survival and disease progression . Our data revealed that EPZ-6438 and lenalidomide combination interfere with MYC transcriptional activity and significantly decreased Ikaros, IRF4, and MYC expression compared with lenalidomide alone. Finally, we discovered in this study that EZH2 inhibition and IMiDs activity act synergistically to alter gene transcription in MM cells, specifically targeting MM oncogenes and cell cycle genes. The deregulation of B and plasma cell transcription factors mediated by PRC2 targeting and IMiDs results in a significant shift between proliferation and quiescence with a significant increase in the percentage of ki67 negative MM cells and a significant induction of quiescent cell features including repression of E2F1 and higher level of p21 Cdk inhibitor [67,68,69,70]. Since MYC activates E2F expression and represses p21 , this reprogramming may be explained by the dual targeting of MYC-IRF4 axis by EPZ-6438/lenalidomide combination.
Altogether, our data demonstrate a deregulation of PRC profiles in MM cells contributing to malignant phenotype. The transcriptional deregulation mediated by PRC2 represent a therapeutic target in MM and a synergy was identified with IMiDs. PRC2 targeting in association with IMiDs could have a therapeutic interest in high-risk MM patients characterized by high EZ score values, reactivating B cell transcription factors and tumor suppressor genes.
DNA methyltransferases inhibitor
Enhancer of zeste homolog 2
Human myeloma cell line
Multiple myeloma cells
Polycomb repressive complex
Herviou L, Cavalli G, Cartron G, Klein B, Moreaux J. EZH2 in normal hematopoiesis and hematological malignancies. Oncotarget. 2016;7:2284.
Barlogie B, Mitchell A, van Rhee F, Epstein J, Morgan GJ, Crowley J. Curing myeloma at last: defining criteria and providing the evidence. Blood. 2014;124:3043–51.
Lohr JG, Stojanov P, Carter SL, Cruz-Gordillo P, Lawrence MS, Auclair D, et al. Widespread genetic heterogeneity in multiple myeloma: implications for targeted therapy. Cancer Cell. 2014;25:91–101.
Bolli N, Avet-Loiseau H, Wedge DC, Van Loo P, Alexandrov LB, Martincorena I, et al. Heterogeneity of genomic evolution and mutational profiles in multiple myeloma. Nat Commun. 2014;5. https://doi.org/10.1038/ncomms3997.
Delmore JE, Issa GC, Lemieux ME, Rahl PB, Shi J, Jacobs HM, et al. BET Bromodomain inhibition as a therapeutic strategy to target c-Myc. Cell. 2011;146:904–17.
Moreaux J, Reme T, Leonard W, Veyrune J-L, Requirand G, Goldschmidt H, et al. Gene expression-based prediction of myeloma cell sensitivity to histone deacetylase inhibitors. Br J Cancer. 2013;109:676–85.
Kalushkova A, Fryknäs M, Lemaire M, Fristedt C, Agarwal P, Eriksson M, et al. Polycomb target genes are silenced in multiple myeloma. PLoS One. 2010;5:e11483.
Croonquist PA, Van Ness B. The polycomb group protein enhancer of zeste homolog 2 (EZH2) is an oncogene that influences myeloma cell growth and the mutant ras phenotype. Oncogene. 2005;24:6269–80.
Zhao F, Chen Y, Li R, Liu Y, Wen L, Zhang C. Triptolide alters histone H3K9 and H3K27 methylation state and induces G0/G1 arrest and caspase-dependent apoptosis in multiple myeloma in vitro. Toxicology. 2010;267:70–9.
Agarwal P, Alzrigat M, Párraga AA, Enroth S, Singh U, Ungerstedt J, et al. Genome-wide profiling of histone H3 lysine 27 and lysine 4 trimethylation in multiple myeloma reveals the importance of polycomb gene targeting and highlights EZH2 as a potential therapeutic target. Oncotarget. 2016;7:6809.
Hernando H, Gelato KA, Lesche R, Beckmann G, Koehr S, Otto S, et al. EZH2 inhibition blocks multiple myeloma cell growth through upregulation of epithelial tumor suppressor genes. Mol Cancer Ther. 2016;15:287–98.
Alzrigat M, Párraga AA, Agarwal P, Zureigat H, Österborg A, Nahi H, et al. EZH2 inhibition in multiple myeloma downregulates myeloma associated oncogenes and upregulates microRNAs with potential tumor suppressor functions. OncoTarget. 2016. https://www.researchgate.net/profile/Mohammad_Alzrigat/publication/312021585_EZH2_inhibition_in_multiple_myeloma_downregulates_myeloma_associated_oncogenes_and_upregulates_microRNAs_with_potential_tumor_suppressor_functions/links/58694b7008ae8fce4917d807.pdf. Accessed 6 Apr 2017.
Pawlyn C, Bright MD, Buros AF, Stein CK, Walters Z, Aronson LI, et al. Overexpression of EZH2 in multiple myeloma is associated with poor prognosis and dysregulation of cell cycle control. Blood Cancer J. 2017;7:e549.
Kurmasheva RT, Sammons M, Favours E, Wu J, Kurmashev D, Cosmopoulos K, et al. Initial testing (stage 1) of tazemetostat (EPZ-6438), a novel EZH2 inhibitor, by the pediatric preclinical testing program: Kurmasheva et al. Pediatr Blood Cancer 2016. doi:https://doi.org/10.1002/pbc.26218.
Moreaux J, Klein B, Bataille R, Descamps G, Maiga S, Hose D, et al. A high-risk signature for patients with multiple myeloma established from the molecular classification of human myeloma cell lines. Haematologica. 2011;96:574–82.
Hose D, Reme T, Hielscher T, Moreaux J, Messner T, Seckinger A, et al. Proliferation is a central independent prognostic factor and target for personalized and risk-adapted treatment in multiple myeloma. Haematologica. 2011;96:87–95.
Hose D, Rème T, Meissner T, Moreaux J, Seckinger A, Lewis J, et al. Inhibition of aurora kinases for tailored risk-adapted treatment of multiple myeloma. Blood. 2009;113:4331–40.
Barlogie B. Total therapy 2 without thalidomide in comparison with total therapy 1: role of intensified induction and posttransplantation consolidation therapies. Blood. 2006;107:2633–8.
Pineda-Roman M, Zangari M, van Rhee F, Anaissie E, Szymonifka J, Hoering A, et al. VTD combination therapy with bortezomib–thalidomide–dexamethasone is highly effective in advanced and refractory multiple myeloma. Leukemia. 2008;22:1419–27.
Mulligan G, Mitsiades C, Bryant B, Zhan F, Chng WJ, Roels S, et al. Gene expression profiling and correlation with outcome in clinical trials of the proteasome inhibitor bortezomib. Blood. 2007;109:3177–88.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:1.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15. https://doi.org/10.1186/s13059-014-0550-8.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27:1739–40.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102:15545–50.
Kassambara A, Rème T, Jourdan M, Fest T, Hose D, Tarte K, et al. GenomicScape: an easy-to-use web tool for gene expression data analysis. Application to investigate the molecular events in the differentiation of B cells into plasma cells. PLoS Comput Biol. 2015;11:e1004077.
Küffner R, Zach N, Norel R, Hawe J, Schoenfeld D, Wang L, et al. Crowdsourced analysis of clinical trial data to predict amyotrophic lateral sclerosis progression. Nat Biotechnol. 2014;33:51–7.
Moreaux J, Reme T, Leonard W, Veyrune J-L, Requirand G, Goldschmidt H, et al. Development of gene expression-based score to predict sensitivity of multiple myeloma cells to DNA methylation inhibitors. Mol Cancer Ther. 2012;11:2685–92.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.
Lawrence M, Huber W, Pagès H, Aboyoun P, Carlson M, Gentleman R, et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9. https://doi.org/10.1371/journal.pcbi.1003118.
Hothorn T, Lausen B. On the exact distribution of maximally selected rank statistics. Comput Stat Data Anal. 2003;43:121–37.
Cummin TE, Araf S, Du M, Barrans S, Bentley MA, Clipson A, et al. Prognostic significance and correlation to gene expression profile of EZH2 mutations in diffuse large b-cell lymphoma (DLBL) in 2 large prospective studies. Hematol Oncol. 2017;35:158–9.
Guirguis AA, Ebert BL. Lenalidomide: deciphering mechanisms of action in myeloma, myelodysplastic syndrome and beyond. Curr Opin Cell Biol. 2015;37:61–7.
Bjorklund CC, Lu L, Kang J, Hagner PR, Havens CG, Amatangelo M, et al. Rate of CRL4CRBN substrate Ikaros and Aiolos degradation underlies differential activity of lenalidomide and pomalidomide in multiple myeloma cells by regulation of c-Myc and IRF4. Blood Cancer J. 2015;5:e354.
Ezponda T, Dupéré-Richer D, Will CM, Small EC, Varghese N, Patel T, et al. UTX/KDM6A loss enhances the malignant phenotype of multiple myeloma and sensitizes cells to EZH2 inhibition. Cell Rep. 2017;21:628–40.
Garapaty-Rao S, Nasveschuk C, Gagnon A, Chan EY, 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.
Viré E, Brenner C, Deplus R, Blanchon L, Fraga M, Didelot C, et al. The Polycomb group protein EZH2 directly controls DNA methylation. Nature. 2005;439:871–4.
Toyota M, Ho C, Ohe-Toyota M, Baylin SB, Issa J-PJ. Inactivation of CACNA1G, a T-type calcium channel gene, by aberrant methylation of its 5′ CpG island in human tumors. Cancer Res. 1999;59:4535–41.
Cha Y, Kim K-J, Han S-W, Rhee YY, Bae JM, Wen X, et al. Adverse prognostic impact of the CpG island methylator phenotype in metastatic colorectal cancer. Br J Cancer. 2016;115:164–71.
Ohkubo T. T-type voltage-activated calcium channel Cav3.1, but not Cav3.2, is involved in the inhibition of proliferation and apoptosis in MCF-7 human breast cancer cells. Int J Oncol. 2012. https://doi.org/10.3892/ijo.2012.1422.
Cohen Y, Merhavi-Shoham E, Avraham RB, Frenkel S, Pe’er J, Goldenberg-Cohen N. Hypermethylation of CpG island loci of multiple tumor suppressor genes in retinoblastoma. Exp Eye Res. 2008;86:201–6.
Bailey CM, Abbott DE, Margaryan NV, Khalkhali-Ellis Z, Hendrix MJC. Interferon regulatory factor 6 promotes cell cycle arrest and is regulated by the proteasome in a cell cycle-dependent manner. Mol Cell Biol. 2008;28:2235–43.
Zengin T, Ekinci B, Kucukkose C, Yalcin-Ozuysal O. IRF6 is involved in the regulation of cell proliferation and transformation in MCF10A cells downstream of notch signaling. PLoS One. 2015;10:e0132757.
Mostovich LA, Prudnikova TY, Kondratov AG, Loginova D, Vavilov PV, Rykova VI, et al. Integrin alpha9 (ITGA9) expression and epigenetic silencing in human breast tumors. Cell Adhes Migr. 2011;5:395–401.
Senchenko VN, Kisseljova NP, Ivanova TA, Dmitriev AA, Krasnov GS, Kudryavtseva AV, et al. Novel tumor suppressor candidates on chromosome 3 revealed by NotI-microarrays in cervical cancer. Epigenetics. 2013;8:409–20.
Nawaz I, Hu L-F, Du Z-M, Moumad K, Ignatyev I, Pavlova TV, et al. Integrin α9 gene promoter is hypermethylated and downregulated in nasopharyngeal carcinoma. Oncotarget. 2015;6:31493.
Mumby M. PP2A: unveiling a reluctant tumor suppressor. Cell. 2007;130:21–4.
Bluemn EG, Spencer ES, Mecham B, Gordon RR, Coleman I, Lewinshtein D, et al. PPP2R2C loss promotes castration-resistance and is associated with increased prostate Cancer-specific mortality. Mol Cancer Res. 2013;11:568–78.
Bi D, Ning H, Liu S, Que X, Ding K. miR-1301 promotes prostate cancer proliferation through directly targeting PPP2R2C. Biomed Pharmacother. 2016;81:25–30.
Wu A-H, Huang Y, Zhang L-Z, Tian G, Liao Q-Z, Chen S-L. MiR-572 prompted cell proliferation of human ovarian cancer cells by suppressing PPP2R2C expression. Biomed Pharmacother. 2016;77:92–7.
Fan Y, Chen L, Wang J, Yao Q, Wan J. Over expression of PPP2R2C inhibits human glioma cells growth through the suppression of mTOR pathway. FEBS Lett. 2013;587:3892–7.
Sugita Y, Ohwada C, Kawaguchi T, Muto T, Tsukamoto S, Takeda Y, et al. Prognostic impact of serum soluble LR11 in newly diagnosed diffuse large B-cell lymphoma: a multicenter prospective analysis. Clin Chim Acta. 2016;463:47–52.
Krzeminski P, Corchete LA, García JL, López-Corral L, Fermiñán E, García EM, et al. Integrative analysis of DNA copy number, DNA methylation and gene expression in multiple myeloma reveals alterations related to relapse. Oncotarget. 2016. https://doi.org/10.18632/oncotarget.13025.
Guijarro MV, Leal JFM, Blanco-Aparicio C, Alonso S, Fominaya J, Lleonart M, et al. MAP17 enhances the malignant behavior of tumor cells through ROS increase. Carcinogenesis. 2007;28:2096–104.
Wu C. Migfilin and its binding partners: from cell biology to human diseases. J Cell Sci. 2005;118:659–64.
Fan J, Ou Y, Wu C, Yu C, Song Y, Zhan Q. Migfilin sensitizes cisplatin-induced apoptosis in human glioma cells in vitro. Acta Pharmacol Sin. 2012;33:1301–10.
Gkretsi V, Papanikolaou V, Zacharia LC, Athanassiou E, Wu C, Tsezou A. Mitogen-inducible Gene-2 (MIG2) and migfilin expression is reduced in samples of human breast cancer. Anticancer Res. 2013;33:1977–81.
Ma H, Wang L, Zhang T, Shen H, Du J. Loss of β-arrestin1 expression predicts unfavorable prognosis for non-small cell lung cancer patients. Tumor Biol. 2016;37:1341–7.
Shen H, Wang L, Zhang J, Dong W, Zhang T, Ni Y, et al. ARRB1 enhances the chemosensitivity of lung cancer through the mediation of DNA damage response. Oncol Rep. 2016. https://doi.org/10.3892/or.2016.5337.
Kumar A, Bhanja A, Bhattacharyya J, Jaganathan BG. Multiple roles of CD90 in cancer. Tumor Biol. 2016;37:11611–22.
Cavo M, Rajkumar SV, Palumbo A, Moreau P, Orlowski R, Blade J, et al. International myeloma working group consensus approach to the treatment of multiple myeloma patients who are candidates for autologous stem cell transplantation. Blood. 2011;117:6063–73.
Anderson KC. The 39th David a. Karnofsky lecture: bench-to-bedside translation of targeted therapies in multiple myeloma. J Clin Oncol. 2012;30:445–52.
Méndez A, Mendoza L. A network model to describe the terminal differentiation of B cells. PLoS Comput Biol. 2016;12:e1004696.
Boi M, Zucca E, Inghirami G, Bertoni F. PRDM1 /BLIMP1: a tumor suppressor gene in B and T cell lymphomas. Leuk Lymphoma. 2015;56:1223–8.
Alinikula J, Nera K-P, Junttila S, Lassila O. Alternate pathways for Bcl6-mediated regulation of B cell to plasma cell differentiation. Eur J Immunol. 2011;41:2404–13.
Inagaki Y, Hayakawa F, Hirano D, Kojima Y, Morishita T, Yasuda T, et al. PAX5 tyrosine phosphorylation by SYK co-operatively functions with its serine phosphorylation to cancel the PAX5-dependent repression of BLIMP1: a mechanism for antigen-triggered plasma cell differentiation. Biochem Biophys Res Commun. 2016;475:176–81.
Yao G. Modelling mammalian cellular quiescence. Interface Focus. 2014;4:20130074.
Wu S, Cetinkaya C, Munoz-Alonso MJ, von der Lehr N, Bahram F, Beuger V, et al. Myc represses differentiation-induced p21CIP1 expression via Miz-1-dependent interaction with the p21 core promoter. Oncogene. 2003;22:351–60.
Bretones G, Delgado MD, León J. Myc and cell cycle control. Biochim Biophys Acta BBA - Gene Regul Mech. 2015;1849:506–16.
Fecteau J-F, Corral LG, Ghia EM, Gaidarova S, Futalan D, Bharati IS, et al. Lenalidomide inhibits the proliferation of CLL cells via a cereblon/p21WAF1/Cip1-dependent mechanism independent of functional p53. Blood. 2014;124:1637–44.
We thank the Microarray Core Facility of IRMB, CHU Montpellier.
This work was supported by grants from French INCA (Institut National du Cancer) Institute (2012-109/087437 and PLBIO15-256), Languedoc Roussillon CRLR (R14026FF), LR-FEDER Hemodiag, Fondation de France (201400047510), ITMO Cancer (MM&TT), SIRIC Montpellier (INCa-DGOS-Inserm 6045), the German Federal Ministry of Education (BMBF) “CAMPSIMM” (01ES1103) and within the framework of the e:Med research and funding concept “CLIOMMICS” (01ZX1309), and the Deutsche Forschungsgemeinschaft (SFB/TRR79; subproject B1). LH is supported by a grant from Labex EpiGenMed.
Availability of data and materials
HMCLs gene expression profiling using Affymetrix U133 plus 2.0 microarrays are deposited in the ArrayExpress public database under accession numbers E-TABM-937 and E-TABM-1088 . Bone marrows were collected from 206 patients treated with high-dose Melphalan (HDM) and autologous stem  cell transplantation (ASCT), and this cohort is termed “Heidelberg-Montpellier” (HM) cohort . Patients’ MMCs were purified using anti-CD138 MACS microbeads (Miltenyi Biotec, Bergisch Gladbach, Germany) and their gene expression profile (GEP) obtained using Affymetrix U133 plus 2.0 microarrays as described . The CEL files and MAS5 files are available in the ArrayExpress public database (E-MTAB-372). The other datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
BM of patients presenting with previously untreated MM at the university hospital of Montpellier was obtained after patients’ written informed consent in accordance with the Declaration of Helsinki and agreement of the Montpellier University Hospital Centre for Biological Resources (DC-2008-417).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. EZH2 expression is correlated with EED and SUZ12 expression. Figure S2 PRC1 members are significantly downregulated in primary MM cells compared with normal bone marrow plasma cells. Figure S3 EZH2 is associated with a poor prognosis in MM. Figure S4 Association between EZH2 expression and patients’ genetic abnormalities. Figure S5 EPZ-6438 leads to H3K27me3 loss in HMCLs. Figure S6 EPZ-6438 leads to H3K27me3 loss in HMCLs. Figure S7 EPZ-6438-induced apoptosis is partially rescued by QVD pan caspase inhibitor. Figure S8 EZH2 inhibition induces mortality of primary MM cells from patients. Figure S9 Correlation between EZH2 expression/H3K27m3 staining and drug response to EZH2 inhibitor in HMCLs and patients. Figure S10 IKZF1 protein level decreases after HMCLs were treated with EPZ-6438 and lenalidomide combination. Figure S11 Lenalidomide targets protein levels after treatment. Figure S12 B cell transcription factors mRNA expression after treatment. Figure S13 PAX5 is a bivalent gene in XG7 HMCL. (PDF 21322 kb)
Table S1. HMCLs molecular characteristics (XLSX 15 kb)
Table S2. EPZ-6438 regulated genes in HMCLs (XLSX 32 kb)
Table S3. GSEA signature enrichment of the 264 EPZ-6438 target genes (XLSX 18 kb)
Table S4. EZH2i target genes are mostly bivalent in XG7 HMCLs (XLSX 10 kb)
Table S6. 67 Lenalidomide + combo upregulated genes (XLSX 17 kb)
Table S7. 31 EPZ-6438 + combo upregulated genes (XLSX 11 kb)
Table S8. Lenalidomide+EPZ-6438-regulated genes associated with GSEA signatures (XLSX 75 kb)
Table S5. H3K27me3-associated and EPZ-6438-regulated genes (XLSX 12 kb)
About this article
Cite this article
Herviou, L., Kassambara, A., Boireau, S. et al. PRC2 targeting is a therapeutic strategy for EZ score defined high-risk multiple myeloma patients and overcome resistance to IMiDs. Clin Epigenet 10, 121 (2018). https://doi.org/10.1186/s13148-018-0554-4
- Multiple myeloma
- Predictive score