- Research article
- Open Access
EZH2 negatively regulates PD-L1 expression in hepatocellular carcinoma
Journal for ImmunoTherapy of Cancer volume 7, Article number: 300 (2019)
Accumulating studies suggest that targeting epigenetic modifications could improve the efficacy of tumor immunotherapy; however, the mechanisms underlying this phenomenon remain largely unknown. Here, we investigated the ability of the epigenetic modifier, enhancer of zeste 2 polycomb repressive complex 2 subunit (EZH2), to regulate the expression of immune checkpoint inhibitor, programmed death-1 ligand 1 (PD-L1) in hepatocellular carcinoma (HCC).
Immunohistochemistry and multiplex immunofluorescence staining were performed to analyze the expression and correlation of EZH2 and PD-L1 in HCC tissues. Immunoblotting, quantitative real-time PCR, flow cytometry, chromatin immunoprecipitation, and dual-luciferase reporter gene assays were performed to evaluate the regulatory roles of EZH2 on PD-L1 expression.
In vitro cell experiments revealed that EZH2 negatively regulated the PD-L1 expression of hepatoma cell lines in IFNγ-dependent manner. Mechanistic studies demonstrated that EZH2 could suppress PD-L1 expression by upregulating the H3K27me3 levels on the promoters of CD274 (encoding PD-L1) and interferon regulatory factor 1 (IRF1), an essential transcription factor for PD-L1 expression, without affecting the activation of the IFNγ-signal transducer and activator of transcription 1 (STAT1) pathway. Clinical samples from HCC patients with immune-activated microenvironments showed negative correlations between EZH2 and PD-L1 expression in hepatoma cells. Multivariate Cox analysis demonstrated that the combination of EZH2 and PD-L1 was an independent prognostic factor for both OS and RFS for patients with HCC.
The epigenetic modificator EZH2 can suppress the expression of immune checkpoint inhibitor PD-L1 by directly upregulating the promoter H3K27me3 levels of CD274 and IRF1 in hepatoma cells, and might serve as a potential therapeutic target for combination of immunotherapy for immune-activated HCC.
Cancer immunotherapy has become an important approach to harness the immune system to fight against cancer cells [1,2,3]. The use of immune checkpoint blockade to induce and restore -immune activation has produced durable clinical responses in the treatment of various tumor types; however, the response rates are still low [4,5,6,7,8]. Increasing evidence has demonstrated the immune-regulatory properties of epigenetic modulators in some preclinical models and in patients with advanced tumors, thus suggesting a rationale for combining epigenetics and immunotherapy [9,10,11]. However, the underlying mechanisms by which epigenetic modifiers influence the efficacy of immunotherapy are poorly understood.
As a key component of polycomb repressor complex 2 (PRC2), enhancer of zeste 2 polycomb repressive complex 2 subunit (EZH2) mediates trimethylation on histone 3 lysine 27 (H3K27me3) and gene silencing, and is involved in various biological processes [12, 13]. Pathological activation of EZH2 histone methyltransferase (HMT) has been observed in tumor transcription programs , which influence cell growth , survival [16, 17], and metastasis [18, 19]. As a potential target for cancer therapy, the therapeutic effects of EZH2 inhibitors are generally interpreted as the consequence of direct reduction of tumor cells (TCs) . Recently, studies have revealed that EZH2 can also contribute to both local and systemic antitumor immunity . For example, EZH2 can affect CD8+ T cell-mediated anti-tumor immunity by regulating the expression of Th-1 chemokines [22, 23] or the infiltration of myeloid-derived suppressor cells (MDSCs) [24, 25]. These data suggested that EHZ2, as an important modifier, could link epigenetic regulation and immune function; however, the mechanisms underlying this phenomenon in human tumors require further explanation.
Lack of programmed death-1 ligand 1 (PD-L1) expression is an important mechanism of resistance to the anti-PD-1/PD-L1 therapies . Thus, many studies have attempted to determine the biological and clinical factors involved in PD-L1 expression [27,28,29,30,31]. Our recent study  revealed that expression pattern of PD-L1 in TCs or monocytes/macrophages (Mo/Mφs) is a predictive biomarker in patients with hepatocellular carcinoma (HCC). The active immune milieu is associated with PD-L1 expression on Mφs; whereas, tumor-expressing PD-L1 may be hampered by certain cell-intrinsic modifications . Given the role of EZH2 in tumor immunity, we speculated that it would have a regulatory effect on PD-L1 expression in HCC.
In the present study, we provided evidence from clinical samples and in vitro cellular experiment that hepatoma-intrinsic EZH2 represses the expression of PD-L1. The results showed that EZH2 regulates the expression of PD-L1 through the epigenetic machinery, and thus could serve as a potential therapeutic target in combination with anti-PD-L1 immunotherapy.
Patients and specimens
Formalin-fixed, paraffin-embedded tissue from 386 patients with pathologically confirmed hepatocellular carcinoma (HCC), who had all received resection of the tumors at the Sun Yat-sen University Cancer Center between 2006 and 2010, were enrolled as previously described . All samples were anonymously coded in accordance with local ethical guidelines (as requested by the Declaration of Helsinki), with written informed consent and using a protocol approved by the Review Board of Sun Yat-sen University Cancer Center. Overall survival (OS) was defined as the interval between surgery and death or between surgery and the last observation for the surviving patients. Relapse-free survival (RFS) was defined as the interval between surgery and the first of recurrence or death, or between surgery and the last observation for patients without recurrence. Tissues were used to construct a tissue microarray (TMA) as described previously . A total of 386 patients who had complete OS and RFS information were used for the survival analysis. The clinical characteristics of all the patients are summarized in Additional file 2: Table S1.
Immunostaining and image analysis
Immunostaining and image analysis were conducted according to our previous reports . In brief, TMA sections were dewaxed in xylene, rehydrated through a decreasing ethanol series, and then placed in 0.3% H2O2 to diminish the activity of endogenous peroxidase. The sections were then heated for antigen retrieval. Following incubation with rabbit anti-human EZH2 (BD Transduction Laboratories, BD Biosciences, San Jose, CA, USA), immunostaining was performed using the EnVision Detection System (DakoCytomation, Carpinteria, CA, USA) following the manufacturer’s instructions. Sections were counter-stained with hematoxylin. Image acquisition was performed using an Eclipse advanced research microscope (Nikon, Melville, NY, USA).
For multiplex immunofluorescence staining of EZH2, PD-L1 (clone: E1L3N™; Cell Signaling Technology, Danvers, MA, USA) and CD68 (DakoCytomation), Tyramide Signal Amplification (TSA) Plus Fluorescence Kits (PerkinElmer, Foster City, CA, USA) combined with immunohistochemistry (IHC) was used. To obtain multispectral images, the stained slides were scanned using the Vectra System (PerkinElmer). The definition of PD-L1 positive expression was the same as that described previously . For colocalization analysis, images were acquired using a laser confocal microscope (Olympus, Essex, UK) and analyzed using FV10-ASW Viewer software (Olympus).
The expression of EZH2 was determined by nuclear EZH2 expression on tumor cells and immunohistochemical scoring of EZH2 was analyzed using the Inform software (PerkinElmer) with the modified Histo-score (H-score), which involves assessing both the intensity of staining (graded as non-staining-0, weak-1, median-2 or strong-3) and the percentage of positive cells (Additional file 1: Figure S1). The range of possible scores was from 0 to 300, quantified by H-score. The correlation of EZH2 and PD-L1 expression was analyzed by χ2 test. The cutoff value for the H-score was set at 35 with the minimum P value to categorize the samples into EZH2 high or low groups.
The human hepatoma cell lines PLC/PRF/5, Huh7, and Hep3B used in this study were purchased from the American Type Culture Collection (Manassas, VA, USA). PLC/PRF/5 cells were cultured in Roswell Park Memorial Institute (RPMI) 1640 medium, and Huh7 and Hep3B cell lines were cultured in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% fetal bovine serum at 37 °C and 5% CO2. Hepatoma cells were treated with recombinant interferon gamma (IFNγ) (Sino Biological Inc.), DZNep (MedChemExpress, Monmouth Junction, NJ, USA), or GSK-126 (MedChemExpress) for different times and at different concentrations.
Monocytes were selected from peripheral blood mononuclear cells using anti-CD14 magnetic beads (Miltenyi Biotec, Bergisch Gladbach, Germany) as described previously .
RNA interference assay
Hepatoma cells were transfected with small interfering RNAs (siRNAs) using Lipofectamine® RNAiMAX Reagent (Invitrogen, Waltham, MA, USA). Reverse transfection was performed according to the manufacturer’s instruction manual. The sequences of the siRNAs are listed in Additional file 2: Table S2.
Cells were collected by 0.25% trypsin digestion, and incubated with Phycoerythrin (PE) conjugated PD-L1 or isotype antibodies (eBioscience, San Diego, CA, USA). The cells were then subjected to flow cytometry.
Quantitative real-time PCR (qPCR)
Total RNA was isolated from cultured cells using TRIZOL (Invitrogen). Reverse transcription and real-time PCR were then performed using 5× All-In-One RT MasterMix (Applied Biological Materials, Richmond, Canada) and a SYBR green real-time PCR kit (Toyobo, Osaka, Japan). Relative quantification was calculated according to the comparative Ct method with normalization to the expression of GAPDH (encoding glyceraldehyde-3-phosphate dehydrogenase). The primers used are listed in Additional file 2: Table S3.
Cells were washed in phosphate-buffered saline (PBS) and suspended in Radioimmunoprecipitation assay (RIPA) buffer (Pierce, Rockford, IL, USA). Supernatant protein concentrations were determined using a BCA protein assay kit (Pierce). Supernatant samples were resolved by 10% or 15% SDS–PAGE depending on the sizes of target proteins, transferred to Immobilon-P polyvinylidene fluoride (PVDF) membranes (Millipore, Billerica, MA, USA) using electroblotting, and then probed with primary antibodies. Membranes were then incubated with horseradish peroxidase-conjugated secondary antibodies. The signals from the immunoreactive proteins were detected using the ECL reagent (Millipore). The information about the antibodies is listed in Additional file 2: Table S4.
Dual-luciferase reporter assay
Huh7 and Hep3B cells pre-transfected with siRNAs, IFNγ, or not, were cotransfected with the pGL3-PD-L1 promoter-luc reporter or pGL3-basic control vectors. pRL-TK was used as an internal control. Cell lysates were harvested for the dual-luciferase assay, which was performed according to the manufacturer’s instructions (Promega, Madison, WI, USA). The primers used are listed in Additional file 2: Table S5.
Construction of Hep3B-shEZH2 cells
To generate a cell line with the stable knockdown of EZH2, lentiviral plasmids carrying a short hairpin RNA (shRNA) targeting EZH2 (VectorBuilder Inc., Shenandoah, TX, USA) were transfected to 293 T cells together with plasmids PMD2.G and pSPAX2 using Lipofectamine 3000 (Invitrogen). After 48 h, culture supernatants were collected, passed through 0.45-μm filters, and mixed with fresh media (1:1) and polybrene (8 μg/ml) to infect Hep3B cells. Cells infected with shEZH2 or control vectors were designated as Hep3B-shEZH2 and Hep3B-vector stable cell lines respectively, and were established using 1 μg/ml puromycin selection. The shRNA-targeted regions in EZH2 were at nt 784–804.
To construction of PD-L1 overexpression plasmid, EZH2-silenced Hep3B cells were first treated with IFNγ for 24 h. Then, mRNA was extracted from these cells, and reverse transcription PCR was performed to obtain cDNA, which was used as template to amplify the coding sequence (CDS) of CD274 (encoding PD-L1). After double restriction enzyme digestion (EcoR I/Kpn I), the CD274 CDS was inserted into plasmid p3 × flag-CMV-14. To obtain more effective expression, the KOZAK sequence was designed into primers, which were as follows:
Forward primer, GGCC GAATTC GCCGCCACC ATGAGGATATTTGCTGTCTTTATATTC;
Reverse primer, CTGA GGTACC TTACGTCTCCTCCAAATGTGTATCACTTTG.
The EZH2 and IRF1 overexpression plasmids were purchased from VectorBuilder company (https://www.vectorbuilder.cn/).
Chromatin immunoprecipitation ChIP
ChIP was performed by using SimpleChIP® Enzymatic Chromatin IP Kit (Cell Signaling Technology). Crosslinking was performed with 1% paraformaldehyde for 10 min. Micrococcal nuclease was added to digest the DNA to lengths of approximately 150–900 bp. The digested DNA solution was sonicated using a Qsonica Q700 sonicator (Qsonica, Newtown, CT, USA) for 1 min at an amplitude of 15%. Protein-DNA complexes were precipitated using specific antibodies against H3K27me3 (Cell Signaling Technology) and IgG control (Cell Signaling Technology). ChIP-enriched chromatin was used for RT-PCR with a SYBR green real-time PCR kit (Toyobo); the data were normalized to the input. The specific primers are listed in Additional file 2: Table S6.
Prediction of CpG Islands on the CD274 promoter and bisulfite sequencing PCR
The sequence of the human CD274 promoter was obtained from the EPD database (https://epd.vital-it.ch/index.php). MethPrimer software (http://www.urogene.org/methprimer/) was used to predict CpG islands and design bisulfite-specific primers for amplification and sequencing. Only one CpG island was predicted on the CD274 promoter (− 2000 bp to + 500 bp). The bisulfite-specific primers sequences were as follows:
Forward primer, ATTTGTTGTTTTGGGTAGAGGTG;
Reverse primer, TAACTCTACTACCCCCTAAACCATC.
Transcriptome profiling analysis
In current study, we used the same batch of transcriptome profiling data previously used , with differences in grouping. According to the status of immune activation which was defined by the expression level of PD-L1 on the infiltrated Mφs, the HCC tissues were divided into two groups: immune-activated and immune-suppressed  (12 cases in each group).
Gene with fold change (FC) more than two between groups and with a Student’s t tests p value < 0.05 was defined as a differentially expressed gene (DEG). And the expression data of all DEGs were analyzed to form heatmap by Funrich 3.1.3 software. The genes upregulated in immune-activated group were further analyzed for GO term enrichment by Funrich 3.1.3 software.
Differences in the means for continuous variables were compared using Student’s t test or analysis of variance, and differences in the proportions were tested using the χ2 test. Kaplan–Meier estimates were calculated and compared using the log-rank test. Cox proportional hazard regression models were applied to evaluate the prognostic variables for OS and RFS. IBM SPSS (version 21.0; IBM Corp., Armonk, NY, USA) statistics software was used for all statistical analyses. All data were analyzed using two-tailed tests unless otherwise specified, and P < 0.05 was considered statistically significant.
EZH2 negatively regulates the IFNγ-induced PD-L1 expression
To evaluate the potential role of EZH2 in regulating PD-L1, we first investigated its expression in HCC tissues (Fig. 1a). IHC staining showed that EZH2 was highly expressed on hepatoma cells in HCC tumors compared with that on parenchyma cells in the non-tumor region. We noted that a few stroma cells were also positive for EZH2. Considering that Mo/Mφs are the major PD-L1-expressing stroma cells in HCC tumors, we performed multiplex staining to analyze EZH2 expression on Mo/Mφs. However, EHZ2 was weakly detected on Mo/Mφs in HCC tumor tissues (Fig. 1b).
Recently, we observed that over 70% of HCC tissues were weak or negative for PD-L1 expression . To evaluate the effect of EZH2 on PD-L1 expression, we used different hepatoma cell lines treated with IFNγ, which is a potent PD-L1 induction factor in multiple tumors [34, 35]. As shown in Fig. 1c, hepatoma cells expressed high level of EZH2, and had only marginally increased PD-L1 expression upon IFNγ stimulation. In contrast, monocytes expressed a low level of EZH2, and showed significantly upregulated PD-L1 expression in response to IFNγ. We then performed RNAi assays to elucidate the regulatory effect of EZH2 on PD-L1. Immunoblotting assays showed that the effective EZH2-targeted siRNAs enhanced IFNγ-induced PD-L1 expression in hepatoma cell (Fig. 1d). Flow cytometry analyses showed that IFNγ-induced PD-L1 expression peaked at 18 to 24 h (Additional file 1: Figure S2) post stimulation. Furthermore, qPCR, immunoblotting, and flow cytometry analyses confirmed that IFNγ-induced PD-L1 expression was upregulated in a variety EZH2-silenced hepatoma cells (Fig. 1e-g). These data suggested that EZH2 could negatively regulate IFNγ-induced PD-L1 expression in hepatoma cells.
EZH2-mediated H3K27me3 on the CD274 promoter to control PD-L1 expression
EZH2-mediated H3K27me3 often leads to epigenetic silencing of target genes ; thus we first tested whether EZH2 could directly regulate the H3K27me3 levels on the CD274 (the gene encoding PD-L1) promoter to suppress IFNγ-induced PD-L1 expression. A concentration gradient of DZNep (an inhibitor of all Sadenosylmethionine (SAM)dependent enzymes, including EZH2)  and GSK126 (a selective inhibitor of EZH2 methyltransferase activity)  were applied to reprogram the epigenetic pathways in hepatoma cells. As expected, GSK126 significantly inhibited the level of H3K27me3 without affecting EZH2 expression, while DZNep simultaneously downregulated the expression of EZH2 (Additional file 1: Figure S3a-b). Inhibiting EZH2 by GSK126 or DZNep treatment effectively increased IFNγ-induced PD-L1 expression in hepatoma cells (Fig. 2a, Additional file 1: Figure S3c-e). ChIP-qPCR analysis revealed that H3K27me3 occupancy on the promoter of CD274 was significantly downregulated in Hep3B-shEZH2 cells compared with that of the control cells (Fig. 2b). These data suggested that H3K27me3 modification at the promoter level is involved in EZH2mediated PD-L1 repression.
It has been reported that EZH2 could also serve as a recruitment platform for the DNA methyltransferase, DNMT1 . From the MethPrimer online database, we predicted CpG islands that are enriched near the transcription start site of CD274  (Fig. 2c). However, these CpG islands are barely methylated in these hepatoma cell lines, as analyzed by bisulfite sequencing PCR (BSP) (Fig. 2d). Furthermore, The Cancer Genome Atlas (TCGA) data showed no correlation between mRNA and promoter methylation levels of CD274 in HCC tumor tissues (Fig. 2e).
It has been reported that the upregulated expression of key molecules such as, interferon gamma receptor 1 (IFNGR1), IFNGR2, Janus kinase 1 (JAK1), and JAK2, activates the IFNγ-STAT1 signaling in response to IFNγ stimulation. We evaluated whether their expression, and the activation of IFNγ-STAT1 signaling, were influenced by EZH2. The results showed that the mRNA levels of these molecules and STAT1 phosphorylation (Fig. 2f-g) were not affected by EZH2 silencing (Additional file 1: Figure S4a-b).
Taken together, these data demonstrated that EZH2 directly regulates the H3K27me3 levels, but not DNA methylation, of the CD274 promoter. Moreover, activation of IFNγ-STAT1 signaling is not influenced by EZH2 expression in hepatoma cells.
EZH2 inhibits PD-L1 expression by epigenetic silencing of IRF1 expression
To further explore the regulatory effect of EZH2 on the CD274 promoter, we constructed a promoter luciferase reporter plasmid containing different truncated versions of the CD274 promoter without H3K27me3 modification (Fig. 3a) and examined the transcription activity of these truncated promoters using dual-luciferase reporter assays. As shown in Fig. 3b, all the designed promoters exhibited similarly enhanced luciferase activity in EZH2-silenced cells when treated with IFNγ, indicating that EZH2 might regulate CD274 transcription by affecting the activities of certain transcription factors (TF) that bind to the P1 truncated promoter. Thirty-nine TFs were predicted to bind to the P1 promoter on the PROMO website [41, 42] and 469 genes who showed co-expression with CD274 with correlation coefficients more than 0.3 were screened out from the cBioportal website (HCC, TCGA, Provisional) [43, 44]. The Venn diagram analysis identified IRF1 as the only potential candidate gene that met both screening criteria (Fig. 3c-d, Additional file 3: Table S7).
Next, dual-luciferase reporter assays were performed to examine the effect of IRF1 on the transcriptional activity of CD274 promoter. The results showed that the luciferase activity of P1 promoter was enhanced by knockdown of EZH2, and this enhancement could be impaired by IRF1 silencing (Fig. 3e). Consistent with this, ectopic expression of IRF1 partially restored the luciferase activity of P1 promoter reduced by overexpression of EZH2 (Fig. 3f). Moreover, re-expression of IRF1 also partly rescued the reduced luciferase activity of P1 caused by knockdown of IRF1 in the EZH2-silenced hepatoma cells (Fig. 3g). We predicted two IRF1 response elements (IREs) on the P1 truncated promoter of CD274 using the JASPAR database (Fig. 3h). To investigate the activities of these response elements, we constructed P1 truncated promoter luciferase reporter plasmids that were deleted for IRE1, IRE2, or both, which were named as △IRE1, △IRE2, and △IREs, respectively (Fig. 3i). Luciferase assays identified that the transcriptional activity of △IRE1, △IRE2, and △IREs were weakened as compared with wild-type P1 in Huh7 and Hep3B cells after stimulation with IFNγ (Fig. 3j). These data indicated that the transcriptional regulation of CD274 by EZH2 is also influenced by the level of IRF1.
Based on the above results, we speculated that IRF1 is downstream of EZH2 to involve in PD-L1 regulation. To examine whether IRF1 is upregulated prior to PD-L1 in EZH2-silencd cells under IFNγ treatment, we evaluated the expression dynamics of IRF1 and PD-L1. As expected, IRF1 was induced and reached its maximum level prior to that of PD-L1 (Fig. 4a, Additional file 1: Figure S5a-b). Moreover, knockdown of IRF1 in EZH2-silenced hepatoma cells decreased the abundance of PD-L1 (Fig. 4b). These results suggested IRF1 should be involved in the EZH2-regulated expression of PD-L1.
Of note, we observed that the expression of IRF1 and PD-L1 were significantly reduced after the re-expression of EZH2 in the EZH2-silenced cells; however, PD-L1 expression was only slightly increased after further ectopic expression of IRF1 (Fig. 4c, Additional file 1: Figure S5c). Considering the epigenetic silencing of EZH2 on PD-L1, we speculated IRF1 has limited induction effect on PD-L1 expression in cells with high expression of EZH2. To further verify this, the induction effect of IRF1 on PD-L1 expression was evaluated in EZH2-silenced hepatoma cells. The results showed that knockdown of IRF1 significantly reduced the expression of PD-L1, and this reduction was rescued by re-expression of IRF1 (Fig. 4d, Additional file 1: Figure S5d).
We also analyzed the impact of EZH2/H3K27me3 axis inhibitors on IRF1 expression. The results showed that GSK126 and DZNep treatments promoted IFNγ-induced IRF1 expression (Fig. 4e, Additional file 1: Figure S5e). ChIP-qPCR analysis was then performed to reveal whether EZH2 could suppress the expression of IRF1 through epigenetic machinery. The results showed H3K27me3 occupancy on the IRF1 promoter was significantly downregulated in Hep3B-shEZH2 cells compared with that of the control (Fig. 4f). These results suggested that IRF1 expression was suppressed by the EZH2-H3K27me3 axis, which leads to inhibition of PD-L1 expression.
Taken together, these results demonstrated that IRF1 promotes the expression of PD-L1, which depends on the epigenetic modification levels of PD-L1 driven by EZH2.
Correlation between EZH2 and PD-L1 expression in HCC tissues
The above findings indicated the effect of EZH2 in regulating IFNγ-induced PD-L1 expression. Next, we analyzed the relationship between EZH2 and PD-L1 expression in HCC tumors. Our recent study showed that Mφ-PD-L1 expression was related to the activated tumor microenvironment . Transcriptome profiling and gene ontology biological analysis confirmed that Mφ-PD-L1+ HCC tumor samples displayed an immune-activated microenvironment and upregulated genes that were mainly involved in the IFNγ-mediated signaling pathway (Additional file 1: Figure S6a and b, Additional file 4: Table S8). Considering that IFNγ stimulation was demonstrated to be required for EZH2-mediated PD-L1 expression in the cell experiments, we analyzed the association between EZH2 and PD-L1 in samples with different microenvironments. Statistical analyses showed a significantly negative correlation between PD-L1 and EZH2 levels on TCs in immune-activated HCC tissues, but not in the total or immune-suppressed samples (Fig. 5a). Multiple immunofluorescence staining revealed that EZH2 protein was barely detected on either PD-L1+ Mφs or TCs (Fig. 5b). Collectively, these results suggested that EZH2 was negatively correlated with PD-L1 expression in the immune-activated HCC tumor microenvironment.
Survival analysis showed that the high levels of EZH2 in hepatoma cells indicated poor OS (P = 0.013) and RFS (P = 0.002) in HCC patients (Fig. 5c, f). Moreover, patients with low or high levels of EZH2 could be further divided into two groups by their PD-L1 expression on TCs. Kaplan–Meier curves showed that patients with PD-L1 expression on TCs had poor OS and RFS (Fig. 5d, OS P = 0.025; Fig. 5g, RFS P = 0.020) in EZH2 low-expression group but not in EZH2 high-expression group (Fig. 5e, OS P = 0.072; Fig. 5h, RFS P = 0.230). Multivariate Cox analysis demonstrated that the combination of EZH2 and PD-L1 was an independent prognostic factor for both OS (P = 0.01) and RFS (P = 0.001) in HCC patients (Table 1).
Accumulating evidence indicates that epigenetic inhibitors could improve the therapeutic efficacy of immune checkpoint blockade [21, 45]; however, the exact regulatory mechanism is not fully understood. In the present study, we reported that the epigenetic modifier EZH2 negatively regulated IFNγ-induced PD-L1 expression in hepatoma cells. Mechanistic studies demonstrated that EZH2 could suppress PD-L1 expression by upregulating the promoter H3K27me3 levels of CD274 (encoding PD-L1) and IRF1 without affecting activation of the IFNγ-STAT1 pathway. A negative correlation between EZH2 and PD-L1 expression on TCs was demonstrated in HCC tissues with an immune-activated microenvironment. Moreover, the combination of EZH2 and PD-L1 on TCs was an independent prognostic factor for OS and RFS in patients with HCC.
Our previous studies demonstrated the differentiated expression of PD-L1 on TCs and Mo/Mφs had opposite clinical impacts on HCC patients . Transcriptome profiling analysis showed that the immune-activated microenvironment was associated with PD-L1 expression on Mo/Mφs but not TCs, indicating that the immune-induced PD-L1 expression by TCs might be regulated by certain intrinsic factors. In the present study, we found that hepatoma cells expressed high levels of EZH2, which abrogated PD-L1 upregulation by IFNγ. In contrast, Mo/Mφs had low levels of EZH2 in HCC tumors and showed significantly upregulated PD-L1 expression after IFNγ stimulation. It should be noted that the criteria for evaluation of EZH2 positive expression were ambiguous in different studies [46, 47]. We determined the expression of EZH2 by nuclear expression on tumor cells, and quantified by H-score. The cutoff value was chosen by the minimum P value that divided patients with diverse clinical outcomes, and also highlighted the correlation between EZH2 and PD-L1 expression. The biological and clinical significance of this cutoff value should be further validated in other cohorts. Taken together, our data revealed that EZH2 acts as an intrinsic modifier that could influence PD-L1 expression in hepatoma cells.
Our mechanistic studies showed that EZH2 epigenetically silenced IFNγ-induced PD-L1 expression by upregulation of the H3K27me3 levels on the promoters of both CD274 and IRF1. Notably, EZH2 did not affect the activation of IFNγ-STAT1 signaling in hepatoma cells, as analyzed by qPCR and immunoblotting assays. This is different from other tumor models, such as in myc-driven prostate cancer, in which EZH2 knockdown restored IFNGR1 expression and further led to activation of IFN-JAK-STAT1 signaling . Recently, it has been reported that the expressions of EZH2 and PD-L1 were positively correlated in lung adenocarcinoma, and DNA methylation could be involved in regulating PD-L1 expression [49, 50]. However, our BSP analysis in hepatoma cells with low constitutive expression of PD-L1 showed almost no methylation on the CD274 promoter regions that were predicted to be CpG islands (Fig. 2d). These data indicated that the intrinsic modifiers and extrinsic inducers for PD-L1 could be varied in different tumor models, for instance, the inflamed liver microenvironments due to chronic viral infection could have great impact on PD-L1 expression.
EZH2 inhibitor intervention experiments showed that both GSK126 and DZNep increased IFNγ-induced PD-L1 expression. ChIP-qPCR analysis showed the downregulation of EZH2 decreased the H3K27me3 levels at the CD274 promoter. These findings suggested that EZH2 regulates PD-L1 expression in hepatoma cells partly by controlling the H3K27me3 levels on the CD274 promoter. However, we observed that GSK126, a direct HMT inhibitor, and DZNep, an indirect HMT inhibitor, displayed certain differences in terms of promoting PD-L1 expression. These results suggested that other mechanisms might be involved in DZNep’s activity, such as proteasomal degradation of PRC2 subunits, inhibition of other methylation reactions, or reactivation of thioredoxin-binding protein 2 (TXNIP), which causes disruption of PRC2 .
We also found that EZH2 could control the H3K27me3 level of the IRF1 promoter, which is prerequisite for IFNγ-induced upregulation of PD-L1 . Downregulated expression of EZH2 significantly inhibited IFNγ-induced upregulation of IRF1. Whereas, the rescue assays showed that downregulated expression of IRF1 significantly decreased the upregulated levels of CD274 promoter luciferase activity and protein expression caused by siEZH2 under IFNγ stimulation. In addition, in line with the results of Lee’s report , we identified two IRF1 response elements (IREs) on the CD274 promoter. We further identified that IRE1 was a more effective than IRE2 in IRF1-mediated transcription of CD274, which might reflect the higher number of IRF1 binding sites in IRE1, as shown by sequence analysis.
In addition, we also tested other inflammatory factors, such as IL-6 and TNFα, for their induction of PD-L1 expression on hepatoma cells. However, neither IL-6 nor TNFα could induce PD-L1 expression on hepatoma cells, with or without EZH2-silencing (Additional file 1: Figure S7a-d). Notably, we found TNFα synergistically promoted the expression of PD-L1 induced by IFNγ (Additional file 1: Figure S7e). This might reflect its role of promoting the stability of PD-L1 . We observed that compared with control cells, the degradation rates of ectopically expressed PD-L1 did not increase in EZH2-silenced hepatoma cells after treatment with the protein synthesis inhibitor cycloheximide (CHX), which suggested that the regulation by EZH2 of IFNγ-induced PD-L1 does not involve the protein stability of PD-L1 in hepatoma cells (Additional file 1: Figure S8).
Drugs targeting EZH2 has been shown to promote the secretion of Th1-type chemokine and subsequent local infiltration of CD8+ T cells in ovarian and colon cancer [22, 23]. A recent study demonstrated a promotion role of CCRK/EZH2-NF-κB/IL-6 axis in HCC by reshaping the balance of MDSCs and CD8+ T cells . Simultaneously blockade of CCRK and PD-L1 could inhibit MDSCs accumulation and engendered CD8+ T cell responses in tumor tissues, resulting in eradication of HCC. In this study, we observed that downregulating the EZH2 could enhance IFNγ-induced PD-L1 expression in hepatoma cells. Survival analysis showed that patients with PD-L1 expression on TCs had poor survival in the EZH2 low-expression group. These findings collectively suggest an important role of EZH2 in reformulating the tumor immune microenvironment.
In conclusion, the present study demonstrated that the epigenetic modifier EZH2 can suppress the expression of immune checkpoint inhibitor PD-L1 by directly upregulating the promoter H3K27me3 levels of CD274 and IRF1 in hepatoma cells and might serve as a potential therapeutic target for combination of cancer immunotherapy for immune-activated HCC.
Availability of data and materials
The datasets used for the current study are available from the corresponding author on reasonable request.
C-X-C motif chemokine ligand 10
Enhancer of zeste 2
Fluorescence activated cell sorting
Trimethylation on histone 3 lysine 27
IRF1 response element
Interferon regulatory factor 1
Myeloid-derived suppressor cells
Phosphate buffer saline
Programmed death-1 ligand 1
Quantitative real-time polymerase chain reaction
Sodium dodecyl sulfate-polyacrylamide gel electrophoresis
Small interfering RNA
The Cancer Genome Atlas
Tumor necrosis factor
Naoum GE, Morkos M, Kim B, Arafat W. Novel targeted therapies and immunotherapy for advanced thyroid cancers. Mol Cancer. 2018;17:51.
Mahoney KM, Rennert PD, Freeman GJ. Combination cancer immunotherapy and new immunomodulatory targets. Nat Rev Drug Discov. 2015;14:561–84.
Inarrairaegui M, Melero I, Sangro B. Immunotherapy of hepatocellular carcinoma: facts and hopes. Clin Cancer Res. 2018;24:1518–24.
Keenan TE, Burke KP, Van Allen EM. Genomic correlates of response to immune checkpoint blockade. Nat Med. 2019;25:389–402.
Topalian SL, Hodi FS, Brahmer JR, Gettinger SN, Smith DC, McDermott DF, Powderly JD, Carvajal RD, Sosman JA, Atkins MB, et al. Safety, activity, and immune correlates of anti-PD-1 antibody in cancer. N Engl J Med. 2012;366:2443–54.
Ribas A, Wolchok JD. Cancer immunotherapy using checkpoint blockade. Science. 2018;359:1350–5.
Burris HA, Infante JR, Ansell SM, Nemunaitis JJ, Weiss GR, Villalobos VM, Sikic BI, Taylor MH, Northfelt DW, Carson WE 3rd, et al. Safety and activity of Varlilumab, a novel and first-in-class agonist anti-CD27 antibody, in patients with advanced solid tumors. J Clin Oncol. 2017;35:2028–36.
Lesokhin AM, Ansell SM, Armand P, Scott EC, Halwani A, Gutierrez M, Millenson MM, Cohen AD, Schuster SJ, Lebovic D, et al. Nivolumab in patients with relapsed or refractory hematologic malignancy: preliminary results of a phase Ib study. J Clin Oncol. 2016;34:2698–704.
Duenas-Gonzalez A, Lizano M, Candelaria M, Cetina L, Arce C, Cervera E. Epigenetics of cervical cancer. An overview and therapeutic perspectives. Mol Cancer. 2005;4:38.
Weintraub K. Take two: combining immunotherapy with epigenetic drugs to tackle cancer. Nat Med. 2016;22:8–10.
Emran AA, Chatterjee A, Rodger EJ, Tiffen JC, Gallagher SJ, Eccles MR, Hersey P. Targeting DNA methylation and EZH2 activity to overcome melanoma resistance to immunotherapy. Trends Immunol. 2019;40:328–44.
Comet I, Riising EM, Leblanc B, Helin K. Maintaining cell identity: PRC2-mediated regulation of transcription and cancer. Nat Rev Cancer. 2016;16:803–10.
Crea F, Hurt EM, Mathews LA, Cabarcas SM, Sun L, Marquez VE, Danesi R, Farrar WL. Pharmacologic disruption of Polycomb repressive complex 2 inhibits tumorigenicity and tumor progression in prostate cancer. Mol Cancer. 2011;10:40.
Hou S, Clement RL, Diallo A, Blazar BR, Rudensky AY, Sharpe AH, Sage PT. FoxP3 and Ezh2 regulate Tfr cell suppressive function and transcriptional program. J Exp Med. 2019;216:605–20.
Mohammad F, Weissmann S, Leblanc B, Pandey DP, Hojfeldt JW, Comet I, Zheng C, Johansen JV, Rapin N, Porse BT, et al. EZH2 is a potential therapeutic target for H3K27M-mutant pediatric gliomas. Nat Med. 2017;23:483–92.
Jin X, Kim LJY, Wu Q, Wallace LC, Prager BC, Sanvoranart T, Gimple RC, Wang X, Mack SC, Miller TE, et al. Targeting glioma stem cells through combined BMI1 and EZH2 inhibition. Nat Med. 2017;23:1352–61.
Ho TH, Kapur P, Eckel-Passow JE, Christie A, Joseph RW, Serie DJ, Cheville JC, Thompson RH, Homayoun F, Panwar V, et al. Multicenter validation of enhancer of Zeste homolog 2 expression as an independent prognostic marker in localized clear cell renal cell carcinoma. J Clin Oncol. 2017;35:3706–13.
Tiwari N, Tiwari VK, Waldmeier L, Balwierz PJ, Arnold P, Pachkov M, Meyer-Schaller N, Schubeler D, van Nimwegen E, Christofori G. Sox4 is a master regulator of epithelial-mesenchymal transition by controlling Ezh2 expression and epigenetic reprogramming. Cancer Cell. 2013;23:768–83.
Li CH, To KF, Tong JH, Xiao Z, Xia T, Lai PB, Chow SC, Zhu YX, Chan SL, Marquez VE, Chen Y. Enhancer of zeste homolog 2 silences microRNA-218 in human pancreatic ductal adenocarcinoma cells by inducing formation of heterochromatin. Gastroenterology. 2013;144:1086–1097.e1089.
Hirukawa A, Smith HW, Zuo D, Dufour CR, Savage P, Bertos N, Johnson RM, Bui T, Bourque G, Basik M, et al. Targeting EZH2 reactivates a breast cancer subtype-specific anti-metastatic transcriptional program. Nat Commun. 2018;9:2547.
Ennishi D, Takata K, Beguelin W, Duns G, Mottok A, Farinha P, Bashashati A, Saberi S, Boyle M, Meissner B, et al. Molecular and genetic characterization of MHC deficiency identifies EZH2 as therapeutic target for enhancing immune recognition. Cancer Discov. 2019;9:546–63.
Peng D, Kryczek I, Nagarsheth N, Zhao L, Wei S, Wang W, Sun Y, Zhao E, Vatan L, Szeliga W, et al. Epigenetic silencing of TH1-type chemokines shapes tumour immunity and immunotherapy. Nature. 2015;527:249–53.
Nagarsheth N, Peng D, Kryczek I, Wu K, Li W, Zhao E, Zhao L, Wei S, Frankel T, Vatan L, et al. PRC2 epigenetically silences Th1-type chemokines to suppress effector T-cell trafficking in colon cancer. Cancer Res. 2016;76:275–82.
Zhou J, Liu M, Sun H, Feng Y, Xu L, Chan AWH, Tong JH, Wong J, Chong CCN, Lai PBS, et al. Hepatoma-intrinsic CCRK inhibition diminishes myeloid-derived suppressor cell immunosuppression and enhances immune-checkpoint blockade efficacy. Gut. 2018;67:931–44.
Huang S, Wang Z, Zhou J, Huang J, Zhou L, Luo J, Wan YY, Long H, Zhu B. EZH2 inhibitor GSK126 suppresses antitumor immunity by driving production of myeloid-derived suppressor cells. Cancer Res. 2019;79:2009–20.
Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, Sosman JA, McDermott DF, Powderly JD, Gettinger SN, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature. 2014;515:563–7.
Wei Y, Zhao Q, Gao Z, Lao XM, Lin WM, Chen DP, Mu M, Huang CX, Liu ZY, Li B, et al. The local immune landscape determines tumor PD-L1 heterogeneity and sensitivity to therapy. J Clin Invest. 2019;129:3347–60.
Chen DP, Ning WR, Jiang ZZ, Peng ZP, Zhu LY, Zhuang SM, Kuang DM, Zheng L, Wu Y. Glycolytic activation of peritumoral monocytes fosters immune privilege via the PFKFB3-PD-L1 axis in human hepatocellular carcinoma. J Hepatol. 2019;71:333–43.
Sabbatino F, Villani V, Yearley JH, Deshpande V, Cai L, Konstantinidis IT, Moon C, Nota S, Wang Y, Al-Sukaini A, et al. PD-L1 and HLA class I antigen expression and clinical course of the disease in intrahepatic Cholangiocarcinoma. Clin Cancer Res. 2016;22:470–8.
Sun C, Mezzadra R, Schumacher TN. Regulation and function of the PD-L1 checkpoint. Immunity. 2018;48:434–52.
Jiang X, Wang J, Deng X, Xiong F, Ge J, Xiang B, Wu X, Ma J, Zhou M, Li X, et al. Role of the tumor microenvironment in PD-L1/PD-1-mediated tumor immune escape. Mol Cancer. 2019;18:10.
Liu CQ, Xu J, Zhou ZG, Jin LL, Yu XJ, Xiao G, Lin J, Zhuang SM, Zhang YJ, Zheng L. Expression patterns of programmed death ligand 1 correlate with different microenvironments and patient prognosis in hepatocellular carcinoma. Br J Cancer. 2018;119:80–8.
Kuang DM, Peng C, Zhao Q, Wu Y, Chen MS, Zheng L. Activated monocytes in peritumoral stroma of hepatocellular carcinoma promote expansion of memory T helper 17 cells. Hepatology. 2010;51:154–64.
Moon JW, Kong SK, Kim BS, Kim HJ, Lim H, Noh K, Kim Y, Choi JW, Lee JH, Kim YS. IFNgamma induces PD-L1 overexpression by JAK2/STAT1/IRF-1 signaling in EBV-positive gastric carcinoma. Sci Rep. 2017;7:17810.
Mandai M, Hamanishi J, Abiko K, Matsumura N, Baba T, Konishi I. Dual faces of IFNgamma in cancer progression: a role of PD-L1 induction in the determination of pro- and antitumor immunity. Clin Cancer Res. 2016;22:2329–34.
Kim J, Shin Y, Lee S, Kim M, Punj V, Lu JF, Shin H, Kim K, Ulmer TS, Koh J, et al. Regulation of breast cancer-induced osteoclastogenesis by MacroH2A1.2 involving EZH2-mediated H3K27me3. Cell Rep. 2018;24:224–37.
Tan J, Yang X, Zhuang L, Jiang X, Chen W, Lee PL, Karuturi RK, Tan PB, Liu ET, Yu Q. Pharmacologic disruption of Polycomb-repressive complex 2-mediated gene repression selectively induces apoptosis in cancer cells. Genes Dev. 2007;21:1050–63.
McCabe MT, Ott HM, Ganji G, Korenchuk S, Thompson C, Van Aller GS, Liu Y, Graves AP, Della Pietra A 3rd, Diaz E, et al. EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations. Nature. 2012;492:108–12.
Vire E, Brenner C, Deplus R, Blanchon L, Fraga M, Didelot C, Morey L, Van Eynde A, Bernard D, Vanderwinden JM, et al. The Polycomb group protein EZH2 directly controls DNA methylation. Nature. 2006;439:871–4.
Li LC, Dahiya R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 2002;18:1427–31.
Farre D, Roset R, Huerta M, Adsuara JE, Rosello L, Alba MM, Messeguer X. Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 2003;31:3651–3.
Messeguer X, Escudero R, Farre D, Nunez O, Martinez J, Alba MM. PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics. 2002;18:333–4.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:pl1.
Loo Yau H, Ettayebi I, De Carvalho DD. The cancer epigenome: exploiting its vulnerabilities for immunotherapy. Trends Cell Biol. 2019;29:31–43.
Zhao Y, Ding L, Wang D, Ye Z, He Y, Ma L, Zhu R, Pan Y, Wu Q, Pang K, et al. EZH2 cooperates with gain-of-function p53 mutants to promote cancer growth and metastasis. EMBO J. 2019;38:e99599.
Cai MY, Tong ZT, Zheng F, Liao YJ, Wang Y, Rao HL, Chen YC, Wu QL, Liu YH, Guan XY, et al. EZH2 protein: a promising immunomarker for the detection of hepatocellular carcinomas in liver needle biopsies. Gut. 2011;60:967–76.
Wee ZN, Li Z, Lee PL, Lee ST, Lim YP, Yu Q. EZH2-mediated inactivation of IFN-gamma-JAK-STAT1 signaling is an effective therapeutic target in MYC-driven prostate cancer. Cell Rep. 2014;8:204–16.
Toyokawa G, Takada K, Tagawa T, Hamamoto R, Yamada Y, Shimokawa M, Oda Y, Maehara Y. A positive correlation between the EZH2 and PD-L1 expression in resected lung adenocarcinomas. Ann Thorac Surg. 2019;107:393–400.
Asgarova A, Asgarov K, Godet Y, Peixoto P, Nadaradjane A, Boyer-Guittaut M, Galaine J, Guenat D, Mougey V, Perrard J, et al. PD-L1 expression is regulated by both DNA methylation and NF-kB during EMT signaling in non-small cell lung carcinoma. Oncoimmunology. 2018;7:e1423170.
Zhou J, Bi C, Cheong LL, Mahara S, Liu SC, Tay KG, Koh TL, Yu Q, Chng WJ. The histone methyltransferase inhibitor, DZNep, up-regulates TXNIP, increases ROS production, and targets leukemia cells in AML. Blood. 2011;118:2830–9.
Lee SJ, Jang BC, Lee SW, Yang YI, Suh SI, Park YM, Oh S, Shin JG, Yao S, Chen L, Choi IH. Interferon regulatory factor-1 is prerequisite to the constitutive expression and IFN-gamma-induced upregulation of B7-H1 (CD274). FEBS Lett. 2006;580:755–62.
Lim SO, Li CW, Xia W, Cha JH, Chan LC, Wu Y, Chang SS, Lin WC, Hsu JM, Hsu YH, et al. Deubiquitination and stabilization of PD-L1 by CSN5. Cancer Cell. 2016;30:925–39.
This work was supported by the National 135 Major Project of China (2017YFA0505803 and 2018ZX10302205); the National Natural Science Foundation of China (81702818, 81772536 and 81572387); the China Postdoctoral Science Foundation (2016 M592586); and the Fundamental Research Funds for the Central Universities under Grant 171gjc32.
Ethics approval and consent to participate
All authors approved and directly participated in the planning, execution and/or analysis of the data presented herein. The research presented here has been performed in accordance with the Declaration of Helsinki and has been approved by the Review Board of Sun Yat-sen University Cancer Center. The patients were informed about the sample collection and had signed informed consent forms.
Consent for publication
The content of this manuscript has not been previously published and is not under consideration for publication elsewhere.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
The expression of EZH2 in HCC tissue microarray. Figure S2. Knockdown of EZH2 upregulates PD-L1 expression induced by IFNγ at different time points. Figure S3. EZH2 inhibitors promotes the IFNγ-induced PD-L1 expression.. Figure S4. The IFNγ-STAT1 signaling is not regulated by EZH2. Figure S5. The EZH2/H3K27me3 axis regulates the expression of IRF1. Figure S6. Transcriptome profiling and gene ontology biological analysis of immune-suppressed and immune-activated HCC tissues. Figure S7. Effect of IL-6 and TNFα on the expression of PD-L1 in hepatoma cells. Figure S8. EZH2 did not effect the protein stability of PD-L1.
Additional file 2: Table S1.
Patient characteristics. Table S2. siRNA target sequences of EZH2 and IRF1 gene. Table S3. Sequences of primers for quantitative real-time PCR. Table S4. Antibodies used in current study. Table S5. Sequences of primers for PD-L1 promoter luciferase reporter plasmids construction. Table S6. ChIP primers used in current study.
Additional file 3: Table S7.
The genes predicted in Fig. 3c.
Additional file 4: Table S8.
List of differentially expressed gene between immune-activated and immune-suppressed HCC tissues screened by transcriptome.
Rights and permissions
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.
About this article
Cite this article
Xiao, G., Jin, LL., Liu, CQ. et al. EZH2 negatively regulates PD-L1 expression in hepatocellular carcinoma. j. immunotherapy cancer 7, 300 (2019). https://doi.org/10.1186/s40425-019-0784-9