Skip to main content

Phosphorylation of EZH2 differs HER2-positive breast cancer invasiveness in a site-specific manner

A Correction to this article was published on 01 December 2023

This article has been updated

Abstract

HER2-positive breast cancer (BC) invasiveness and drug-resistance issue is the critical treatment obstacle recently. We investigated the total and phosphorylated status EZH2 expression in database and BC tissue microarray. We demonstrated for the first time that EZH2 is distributed both in cytoplasm and nucleus of breast cancer cells in a phosphorylation site-specific manner. High expressed-EZH2 cases more frequently had an advanced clinical stage (lymph node metastasis) and aggressive features than EZH2-low cases, potentially indicating the high risk of HER2-positive BC (p < 0.05). Notably, highly expressed phosphorylated EZH2 is differently located in cytoplasm or nucleus in a site-specific manner in breast cancer cells. Nucleus-located pEZH2-S21 is expressed in invasive and lymph node metastatic HER2-positive BC cases (p = 0.144, p = 0.001). Cytoplasmic pEZH2-T487 is correlated with HER2 positive status (p = 0.014).

In conclusion, high expression of nucleus-located EZH2 might be a predictor of invasive BC. Activation of phosphorylated EZH2-S21 site in nucleus would be a potential predictor of HER2-positve BC and poor efficacy of HER2-target therapy. These results point to a PRC2-independent non-epigenetic mechanism and therapeutic strategy of EZH2 in HER2-positive BC.

Peer Review reports

Introduction

Human epidermal growth factor receptor 2 (HER2)-positive breast cancer (BC) accounts for approximately 15–20% of all types of BC patients. Over the years, HER2-positive BC is associated with an increased risk of recurrence and metastasis and poor outcome because of its more aggressive biological behaviors [1]. Although there are emerging dramatic improvements of anti-HER2 target treatment in the recent years, HER2-positive BC invasiveness and drug-resistance still remain to be a tricky potential molecular issue in current era [2]. Therefore, it is an unmet need for seeking for the precise molecular mechanism of HER2-positive BC aggressive biological behavior.

Enhancer of Zeste Homolog 2(EZH2) is a histone methyltransferase that catalyzes the histone H3 lysine 27 methylation. It is one of the core subunits of the polycomb-repressive complex 2(PRC2), which is consisted of several key components—EED, SUZ12 and RpAp46/48 [3, 4]. Numerous studies have shown that EZH2 is involved in a wide range of biological processes, such as tumor development, growth [5], metastasis [6, 7], apoptosis [8], angiogenesis [9], stem cell renewal/ maintenance [10], and treatment-resistance [11,12,13]. Indeed, most of these functional integrities of EZH2 are based on epigenetic mechanism. Therefore, EZH2 is considered as a transcriptional silencer in tumors, which has been well established. Nevertheless, minority is focused on the non-epigenetic role in breast cancer, especially the posttranscriptional modification (PTM)---phosphorylation. In comparison to an epigenetic repressor of tumor-related genes in nucleus, little is known about the function of cytoplasmic EZH2 in tumor cells, breast cancer cells particularly.

Apart from epigenetic mechanism and methylation, phosphorylation of EZH2 function in breast tumors has been increasingly taken consideration and revealed in recent years [14]. Accordingly, the subcellular location of EZH2 in tumor cell has also been concerned by researchers. For instance, cytoplasmic pEZH2 T367 expression and subcellular localization could stratify metaplastic carcinoma subtypes [15]. In addition, pEZH2-T367 induces EZH2 cytoplasmic localization to promote breast cancer metastasis [16]. Similarly, it is found that high cytoplasmic EZH2 expression is significantly associated with short overall survival in cholangiocarcinoma [17]. Evidence shows that EZH2 in the cytoplasm is closely coupled to cancer stem cell properties [18], thus inhibition of cytoplasmic EZH2 induces antitumor activity in NSCLC [19]. As aforementioned, the extranuclear function of EZH2 in tumorigenesis has been gradually noticed and validated in solid tumors, less has been explored in HER2-positive BC concerning EZH2, especially in aggressiveness and drug-resistance.

Therefore, we examined the non-traditional biological function of phosphorylated EZH2 in breast cancer. In this work, we take the lead to reveal a unique and non-classical role of EZH2 phosphorylation in HER2-positive BC, particularly in a site-specific manner with different cell components.

Results

EZH2 is highly expressed in malignancies and correlated with poor prognosis of breast cancer

We explore the expression of EZH2 gene in pan-cancer tissues as well as normal tissues in The Cancer Genome Atlas (TCGA) RNAseq data set (TCGA). Data shows that EZH2 expression mRNA level in pan-cancer including breast cancer tissue is higher than the corresponding normal tissue in TCGA dataset (Fig. 1A-B). Consistently, the EZH2 gene transcript from Gene expression profiling interactive analysis (GEPIA2) dataset and TNMplot database also provide that higher level in breast cancer and metastasis specimens of EZH2 than normal mammary glands epithelium specimens (Fig. 1C-D). To explore the impact of EZH2 in breast cancer on BC patients’ survival, we analysis the correlation between EZH2 mRNA level and survival prognosis in different subtypes of BC patients in K-M plotter database (http://kmplot.com/analysis). Data shows that Patients with lower EZH2 expression had longer overall survival (OS), recurrence-free survival (RFS) and distant metastasis-free survival (DMFS) in total breast cancer dataset (HR = 1.56, p < 0.05; HR = 1.73, p < 0.05; HR = 1.71, p < 0.05). Of note, patients with higher EZH2 expression level had poor OS for luminal A and HER2-positive BCs (HR = 2.14, p < 0.05; HR = 1.75, p = 0.055) (Fig. 1E).

Fig. 1
figure 1

EZH2 expression in tumor and normal tissues and its prognosis of BC in database. A.B EZH2 expression level in pancancers and normal samples in TCGA dataset. C. EZH2 transcript in pancancers in GEPIA2 database. D. EZH2 mRNA level in TNMplot database. E. BC patients’prognosis of EZH2 expression in K-M plotter

EZH2 exerts cytoplasm-nucleus sublocation in breast cancer tissues, which associated with invasiveness

Two tissue microarray (TMA) slides of breast cancer specimen totally containing 16 benign tissues, 7 DCIS tissues and 90 breast cancer tissues (113 cases in total) were collected and produced by our department and lab from the Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou. The total specimen samples were classified as normal mammary and benign lesion, Ductal carcinoma in situ (DCIS), Invasive ductal carcinoma (IDC) and other invasive type breast carcinoma for further investigation. Besides, we stratified breast tumors into luminal-type, HER2-positive and Triple negative breast cancer (TNBC) subtypes for the following study.

It could be observed that EZH2 was highly expressed in both cytoplasm and nucleus of breast cancer tissues. In benign lesion, EZH2 staining is almost negative or weak in IHC assay. In DCIS tissues, 7 cases were almost positive for EZH2 expression as cytoplasm location (Fig. 2A, B). Strikingly, we found that the sublocation of EZH2 staining was associated with lymph node metastasis in HER2-positive invasive BC tissues—82.6% (19/23) cases were positive for EZH2 expression in LN negative IDC specimens, as cytoplasm-located staining while 53.3% (8/15) cases were positive for EZH2 expression in LN positive IDC specimens as nucleus-located staining (Fig. 2C, p < 0.05). In TCGA database, Fig. 2D show that EZH2 expression in breast cancer is significantly higher than normal breast tissue. Simultaneously, EZH2 expression level is statistically higher in aggressive subtypes of BC—TNBC and HER2-positive subtype (p < 0.05, Fig. 2E). To determine EZH2 expression in metastatic breast cancer, we make comparison of EZH2 in each breast cancer N stage group (N0:516 cases, N1:362cases, N2:120cases, N3:77cases) with normal mammary gland tissue (n = 114 cases) based on TCGA database. Consistently, data shows that EZH2 expression level is statistical significantly higher in high N stage of BC than normal breast epithelium group (p < 0.05, Fig. 2F). Consequently, the cytoplasm-nucleus sublocation of EZH2 in BC tissues is closely related to BC invasiveness, HER2-positive BC, in particular.

Fig. 2
figure 2

EZH2 expression in BC tissues and its clinical significance of BC patients. (A) Representative IHC staining of EZH2 expression in paracancerous breast tissue, DCIS, lymph node negative-IDC and lymph node metastasis IDC tissue as low and high expression with different sublocations. (The immunohistochemical images were at 5,20,63x magnification) (B) EZH2 expression IHC score in paracancerous breast tissue, DCIS and other invasive type tissues by statistical analysis. (C) Nucleus-cytoplasm location of EZH2 in LN metastasis BC tissue. D,E,F. EZH2 expression in BC tissue, different subclasses and from TCGA database

Phosphorylation of EZH2 T487 and S21 site display obvious importance in HER2-positive BCs

All phosphorylated sites of EZH2 are displayed in schematic diagram of Fig. 3A from PhosphoSitePlus database. The sequence of each S, T, Y amino acid phosphoSite of EZH2 is listed from PhosphoELM database (http://phospho.elm.eu.org). We analysed pEZH2-T487 site of breast cancer in Clinical Proteomic Tumor Analysis Consortium (CPTAC) database. It is found that pEZH2-T487 expression level is increased in breast cancer specimen (n = 125) than normal tissues (n = 18), especially in luminal (n = 64) and TNBC (n = 16) subtype (Fig. 3B). Moreover, pEZH2-T487 is positively associated with late BC tumor stage—pEZH2-T487 level is gradually elevated in stage I, II and III in CPTAC BC samples than normal tissues (Fig. 3B). In addition, it suggested that higher level of pEZH2-T487 is presented in mTOR altered BC tissues, which might imply underlying signaling pathway mechanisms uncovered (Fig. 3B).

Fig. 3
figure 3

Phosphorylated EZH2 (S21 and T487) profile in breast cancer tissues in database and clinical specimens. (A) All phosphorylated sites and sequences of EZH2 are displayed in schematic diagram. (B) Phosphorylated EZH2-T487 level in normal and tumor tissue, different subtypes of BC, different tumor stages and mTOR pathway alteration in CPTAC dataset. (C) phosphorylated EZH2 (S21 and T487) in anti-HER2 treatment resistance BC recurrence tissue

Early in 2018, the author and her colleges have revealed that transcription factor STAT3 enhances lincRNA HOTAIR transcription by interacting with pEZH2-Serine21, instead of pEZH2-T487 in head and neck squamous cell carcinoma (HNSCC), thus promoting tumor progression and cisplatin resistance [20]. To validate the phosphorylation status of EZH2 in invasive metastatic BC, we examined pEZH2-T487 and pEZH2-S21 expression by immunofluorescence two individual HER2-positive BC chest recurrence case tissue after multi-line anti-HER2 treatment. With one accord, images show that pEZH2-T487 is almost intensively staining in cytoplasm while pEZH2-S21 is almost intensively staining in nucleus in the anti-HER2 treatment recurrence BC tissues (Fig. 3C). Subsequently, we examined pEZH2-S21, pEZH2-T487 and HER2 staining in BC specimens by IF and IHC assay. The clinicopathological characteristics of BC patients in TMA were collected and analyzed statistically using Pearson’s chi-square test. In Table 1, data shows that cytoplasmic pEZH2-T487 is correlated with HER2 positive status of breast cancer tissues (p = 0.014). Meanwhile, nucleus-located pEZH2-S21 is expressed in invasive and lymph node metastatic HER2-positive BCs (p = 0.001, Table 2). Representative HE and HER2 Immunohistochemical images of benign lesion of breast (case No.60, No.58 on TMA), HER2-negative IDC (case No.25, No.6 on TMA) and HER2-positive IDC (case No.15, No.54 IDC on TMA) are displayed in Fig. 4A and B, respectively. Consistently, it is suggested that both pEZH2-S21(p = 0.059) and pEZH2-T487 (p < 0.0001) expression level are increased in IDC than benign lesion, DCIS and other invasive type tissue (Fig. 4E, F). Remarkably, pEZH2-S21(p < 0.05) and pEZH2-T487 (p < 0.01) are strictly correlated with HER2 status, respectively (Fig. 4E, F). Accordingly, we observed an increased EZH2 phosphorylation levels of T487 and S21 site in breast cancer tissues. Interestingly, the T487 site of phosEZH2 and S21 site of phosEZH2 are localized differently in subcellular segment of tumor tissues—cytoplasm and nucleus, respectively (Fig. 4C, D). The intensity of pEZH2-T487 is significantly stronger in HER2-positive IDC as cytoplasmic staining than in HER2-negative IDC, as well as in normal breast tissue (Fig. 4D). Meanwhile, the intensity of pEZH2-S21 is significantly stronger in HER2-positive IDC as nuclear staining than in HER2-negative IDC, as well as in normal breast tissue (Fig. 4C). This phenomenon of phosphorylated EZH2 cytoplasm-nucleus sublocation difference in site-specific manner is of statistical significance in our TMA data (p < 0.05).

Table 1 pEZH2-T487 in breast cancer tissue
Table 2 pEZH2-S21 in breast cancer tissue
Fig. 4
figure 4

Phosphorylation of EZH2 (Serine21, Tyrosine487 site) correlates with invasive and HER2-positive BC in different subcellular components. A.B Representative HE and HER2 IHC staining images from human breast benign lesion, HER2-negative IDC and HER2-positive IDC. C. Representative pEZH2-S21 IF staining images from human breast benign lesion, HER2-negative IDC and HER2-positive IDC. D. Representative pEZH2-T487 IF staining images from human breast benign lesion, HER2-negative IDC and HER2-positive IDC. E.F Expression level of pEZH2-S21 and pEZH2-T487 in breast benign lesion, DCIS, IDC and other type tissue by statistical analysis. Expression level of pEZH2-S21 in HER2-negative and HER2-positive BC tissue by statistical analysis

Regarding to the target genes of EZH2 in breast cancer, we performed Bioinformatics analysis based on TCGA-BRCA using STRING, and metascape (http://metascape.org) database, which illustrated in Supplementary Fig. 1. Genes such as DBF4 (R = 0.72, p < 0.01), CDCA8 (R = 0.75, p < 0.01), RAD54L (R = 0.71, p < 0.01) and STMN1 (R = 0.79, p < 0.01) are correlated with EZH2 expression, which have been displayed in sup. Figure 1 C. Ezh2 has also involved in several signaling pathways, such as mitosis cell cycle, DNA-damage and cellular response, etc. using Enriched Ontology Clusters (sup. Figure 1D, E) .

More specifically, we employed a highly aggressive breast cancer cell line (BT549) treated with EZH2 inhibitor (Dznep) for RNA sequencing assays, attempting to identify its downstream target genes in BC in vitro. Data was analyzed by Gene ontology (GO) enrichment and Kyoto Ency- clopedia of Genes and Genomes (KEGG) pathway. GO analysis for differentially expressed genes in EZH2 inhibition group and control group. Differentially expressed genes (DEGs) driven by EZH2 is obtained as gene correlation network figure exhibited as “pathway in cancer” and “metabolic pathway” panel, as representatives (Supplementary Fig. 2).

The downstream target genes of EZH2 would provide molecular clues for our further investigation.

Taken together, our findings demonstrate for the first time that the expression and subcellular localization of phosphorylated EZH2 differs HER2-positive breast cancer invasiveness in a site-specific manner, suggesting a correlation with lymph node metastasis. Therefore, further clinical diagnostic and treatment implications would be provided for HER2-positive BCs.

Discussion

The methyltransferase EZH2 play vital role in chromatin conformation regulation and gene epigenetic transcription [21]. EZH2 has been implicated to have critical role in multiple carcinogenesis and tumor development and progression, including breast cancer [22]. Although various studies have been investigated in EZH2 behavior in malignancies, few literatures have been focused on the cytoplasm-nucleus sublocation of EZH2 in tumor cells. Interestingly, we initially report and uncover a novel evidence that phosphorylated EZH2 behaves as an actionable target in HER2-positive BC in a site-specific manner. In this study, we present significant evidence to reveal that phosphorylated EZH2 is differently located in cytoplasm and nucleus of breast cancer cells in tumor tissue specimen in a site-specific manner, which suggesting pEZH2 is critical for identify invasiveness and metastasis in HER2-positive BC.

The biological function of phosphorylated EZH2 in tumor cells have been increasingly raising concern recently. Data shows that it is multifaceted and context-dependent. For instance, phosphorylation of EZH2 at T311 inhibits PRC2 methyltransferase activity to suppress tumorigenesis [23]. In natural killer/T-cell lymphoma, EZH2 S220 phosphorylation is increased by MELK to mediate its sensitivity to bortezomib [24]. CDK2-mediated phosphorylation of EZH2 at T416 activates EZH2 to silence target genes—ERα gene (ESR1), leading to in high-grade serous ovarian carcinoma (HGSOC) and TNBC [25]. Meanwhile, CDK2-mediated EZH2 phosphorylation(pT416) drives tumorigenesis— converts the luminal breast cancer to TNBC [26]. EZH2 phosphorylation at T372 reduces ovarian cancer cell proliferation, migration and tumor formation [27]. Calcium/Calmodulin dependent protein kinase II alpha (CAMK2A), a key calcium signaling molecule, phosphorylates EZH2 at T487 with suppression of EZH2 activity to support Tumor initiating cells (TIC) maintenance in lung adenocarcinoma [28]. Phosphorylation of EZH2 at T435 and T487 by CDK1 and CDK2 epigenetically silences target genes in breast cancers [29, 30]. Mechanically, phosphorylation of Thr-487 site, which mediated by CDK1 is necessary for EZH2 ubiquitination and subsequent degradation by ubiquitin-proteasome pathway [31]. Hence, degradation of EZH2 could attenuate breast cancer invasion and metastasis in MCF-7 and MDA-MB-231 cells [32]. Apart from EZH2 degradation mediated by pEZH2-T487 in luminal and TNBC, the biological function of cytoplasm-located pEZH2-T487 in HER2-positive BC progression and metastasis remains unexplored. In our study, we elucidate that cytoplasmic pEZH2-T487 is correlated with HER2 positive status of breast cancer tissues, which might be a potential regulatory target of HER2-positive BC.

As for serine21 of phosphorylated EZH2, views differ in the modulation mechanism between epigenetic and non-epigenetic characteristics. AKT-mediated EZH2 phosphorylation at serine21 drives hepatocarcinogenesis and associates with tumor recurrence and poor survival [33]. In the literature, pS21-EZH2 is functionally important in transactivating Bcl6 and the T follicular helper (TFH) cell program, revealing a critical contribution by EZH2 phosphorylation to a non-cancerous, physiological process [34]. Besides, EZH2 activity resulting from posttranslational phosphorylation at the serine-21 site is responsible for the increased enrichment of H3K27me3 at the RECK promoter region, which is related to tumor metastasis and angiogenesis [35]. In a recent study, dephosphorylates EZH2 at S21 could activate its methyltransferase function. Hence, EZH2-S21 phosphorylated-controlled Epithelial–Mesenchymal Transition (EMT) genes are upregulated during fibrotic disease progression of human eyes—3 major extracellular matrix (ECM) target genes: Col1A1, MMP17, and POSTN [36]. However, pEZH2-S21 biological behavior in HER2-positive BCs has not yet been further explored and elucidated. In the present study, we initially suggest that pEZH2-S21, which preferably located in nucleus of breast cancer, is intensively activated in invasive and lymph node metastatic HER2-positive BCs.

It is reported that EZH2-mediated histone modification modulates the response to HER2-targeted therapies via the protein phosphatase 2 A (PP2A) regulatory subunit PPP2R2B regulation [37]. In 2022, researchers demonstrated that inhibition of EZH2 and HDAC could reactivate IFI16-directed immune response of HER2-positive BC patients, changing “cold tumor” into “hot tumor”, and sensitizing anti-HER2 treatment responses [38]. Accordingly, we would speculate that cytoplasic-phosEZH2 and nuclear-phosEZH2 might serve its biological functions separately in invasive BCs, HER2-positive BC particularly.

In conclusion, our study first put forward the clinical significance of phosphorylated EZH2 in HER2-positive BCs. The expression and subcellular localization of phosphorylated EZH2 differs HER2-positive breast cancer invasiveness in a site-specific manner, which associated with metastasis. Further clinical potential therapeutic strategy would be explored in HER2-positive BC to improve the efficacy of target therapy.

Methods

Tissue specimen, tissue microarray (TMA) construction and reagents

Breast tumor tissues were collected from our department who underwent breast surgery. In total, 113 patients were recruited from the Second Affiliated Hospital of Zhejiang University between 2020 and 2021. All patients had signed informed consent of the study. Tissue microarray construction was performed with breast benign lesion, paracancerous tissue and breast tumor samples. All these were embedded into paraffin TMA blocks with a manual tissue arrayer. Two TMA cores of 1.0 mm in diameter were sampled from a cohort of 113 breast samples (16 benign tissues, 7 DCIS tissues and 90 breast cancer tissues). Sections of 4 μm were cut and stained with H&E. All the clinical data of these 113 patients has been collected, reviewed and analyzed statistically to testify the clinical importance of sub-location of (p)EZH2 in BC patients, which exhibited in Tables 1 and 2. Three randomly chosen fields were acquired per condition with identical microscope settings and images were scanned measured using Caseviewer software.

Bioinformatics database analysis

Gene expression and survival prognosis analysis

The GEPIA2 database (http://gepia2.cancer-pku.cn) could analysis gene expression profiles from TCGA dataset and Genotype-Tissue Expression (GTEx) projects. EZH2 copy number alteration in pancancers and normal samples from TCGA and GEPIA2 webserver could be obtained by Boxplot [39]. Survival analysis of breast cancer OS, RFS and DMFS with EZH2 expression level were detected by Kaplan-Meier plots. We identified EZH2 expression levels in different subtypes of BC prognosis based on K-M plotter database.

The sequence of phosphorylated EZH2 and schematic diagram is displayed from PhosphoELM and PhosphoSitePlus website database. The association of pEZH2-T487 in breast cancer with clinical features is conducted from CPTAC dataset.

The R package, clusterProfiler, was then used to perform Gene Ontology (GO) and Kyoto Ency- clopedia of Genes and Genomes (KEGG) [40,41,42] analyses of cancer-related biological processes or pathways affected by EZH2.

Immunohistochemical (IHC) and immunofluorescence (IF)staining and scoring analysis

They were and paraffin-embedded and sectioned, followed by deparaffinization, rehydration and antigen retrieval. Primary antibodies specific for EZH2, pEZH2-S21, pEZH2-T487 and HER2 were incubated while a negative staining was used as a control. TMA slides were examined and scored by two individual trained pathologists using semi-quantitative IHC and IF scoring system. The IHC score is mainly classified into intensive, moderate and negative based on the staining intensity and fraction of stained cells [43, 44]. Images were acquired with a Zeiss fluorescence microscope. Fields were randomly chosen with the DAPI filter. All screen settings were the same for each condition.

Statistical analysis

Statistical analysis was performed by SPSS and GraphPad software. To evaluate the significance of phosphorylated EZH2 association with clinicopathological data, Pearson’s Chi-square test and t test were conducted to examine clinical characteristic significance. P < 0.05 was considered as statistically significant.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data Availability

Gene expression data described in this study are deposited in TCGA and GTE database by GEPIA2 (http://gepia2.cancer-pku.cn) website. Survival analysis of EZH2 expression in breast cancer is exported from K-M plotter (https://kmplot.com/analysis/) database. The phosphorylated sites of EZH2 are analyzed from PhosphoSitePlus (http://phosphosite.org) database, PhosphoELM (http://phospho.elm.eu.org) and UALCAN (https://ualcan.path.uab.edu) Clinical Proteomic Tumor Analysis Consortium (CPTAC) database.

Additional relevant data that support the findings of this study are available from the corresponding authors upon reasonable request.

Change history

References

  1. Exman P, Tolaney SM. HER2-positive metastatic breast cancer: a comprehensive review. Clin Adv Hematol Oncol. 2021;19(1):40–50.

    PubMed  Google Scholar 

  2. Mishra A, Hourigan D, Lindsay AJ. Inhibition of the endosomal recycling pathway downregulates HER2 activation and overcomes resistance to tyrosine kinase inhibitors in HER2-positive breast cancer. Cancer Lett. 2022;529:153–67.

    CAS  PubMed  Google Scholar 

  3. Sun S, Yu F, Xu D, Zheng H, Li M. EZH2, a prominent orchestrator of genetic and epigenetic regulation of solid tumor microenvironment and immunotherapy. Biochim Biophys Acta Rev Cancer. 2022;1877(2):188700.

    CAS  PubMed  Google Scholar 

  4. Boyer LA, Plath K, Zeitlinger J, Brambrink T, Medeiros LA, Lee TI, Levine SS, Wernig M, Tajonar A, Ray MK, et al. Polycomb complexes repress developmental regulators in murine embryonic stem cells. Nature. 2006;441(7091):349–53.

    CAS  PubMed  Google Scholar 

  5. Zhang H, Qi J, Reyes JM, Li L, Rao PK, Li F, Lin CY, Perry JA, Lawlor MA, Federation A, et al. Oncogenic deregulation of EZH2 as an opportunity for targeted therapy in Lung Cancer. Cancer Discov. 2016;6(9):1006–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Hu FF, Chen H, Duan Y, Lan B, Liu CJ, Hu H, Dong X, Zhang Q, Cheng YM, Liu M, et al. CBX2 and EZH2 cooperatively promote the growth and metastasis of lung adenocarcinoma. Mol Ther Nucleic Acids. 2022;27:670–84.

    CAS  PubMed  Google Scholar 

  7. Ding S, Wang X, Lv D, Tao Y, Liu S, Chen C, Huang Z, Zheng S, Wei Y, Kang T, et al. EBF3 reactivation by inhibiting the EGR1/EZH2/HDAC9 complex promotes metastasis via transcriptionally enhancing vimentin in nasopharyngeal carcinoma. Cancer Lett. 2022;527:49–65.

    CAS  PubMed  Google Scholar 

  8. Zhou X, Ren Y, Kong L, Cai G, Sun S, Song W, Wang Y, Jin R, Qi L, Mei M, et al. Targeting EZH2 regulates tumor growth and apoptosis through modulating mitochondria dependent cell-death pathway in HNSCC. Oncotarget. 2015;6(32):33720–32.

    PubMed  PubMed Central  Google Scholar 

  9. Lu C, Han HD, Mangala LS, Ali-Fehmi R, Newton CS, Ozbun L, Armaiz-Pena GN, Hu W, Stone RL, Munkarah A, et al. Regulation of tumor angiogenesis by EZH2. Cancer Cell. 2010;18(2):185–97.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Suvà ML, Riggi N, Janiszewska M, Radovanovic I, Provero P, Stehle JC, Baumer K, Le Bitoux MA, Marino D, Cironi L, et al. EZH2 is essential for glioblastoma cancer stem cell maintenance. Cancer Res. 2009;69(24):9211–8.

    PubMed  Google Scholar 

  11. Xu C, Hao K, Hu H, Sheng Z, Yan J, Wang Q, Yu L. Expression of the enhancer of zeste homolog 2 in biopsy specimen predicts chemoresistance and survival in advanced non-small cell lung cancer receiving first-line platinum-based chemotherapy. Lung Cancer. 2014;86(2):268–73.

    PubMed  Google Scholar 

  12. Chen F, Liu J, Song X, DuCote TJ, Byrd AL, Wang C, Brainson CF. EZH2 inhibition confers PIK3CA-driven lung tumors enhanced sensitivity to PI3K inhibition. Cancer Lett. 2022;524:151–60.

    CAS  PubMed  Google Scholar 

  13. Zhang L, Li HT, Shereda R, Lu Q, Weisenberger DJ, O’Connell C, Machida K, An W, Lenz HJ, El-Khoueiry A, et al. DNMT and EZH2 inhibitors synergize to activate therapeutic targets in hepatocellular carcinoma. Cancer Lett. 2022;548:215899.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Kim E, Kim M, Woo DH, Shin Y, Shin J, Chang N, Oh YT, Kim H, Rheey J, Nakano I, et al. Phosphorylation of EZH2 activates STAT3 signaling via STAT3 methylation and promotes tumorigenicity of glioblastoma stem-like cells. Cancer Cell. 2013;23(6):839–52.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. McMullen ER, Skala SL, Gonzalez ME, Djomehri S, Chandrashekar DS, Varambally S, Kleer CG. Subcellular localization of EZH2 phosphorylated at T367 stratifies metaplastic breast carcinoma subtypes. Breast Cancer. 2021;28(2):496–505.

    PubMed  Google Scholar 

  16. Anwar T, Arellano-Garcia C, Ropa J, Chen YC, Kim HS, Yoon E, Grigsby S, Basrur V, Nesvizhskii AI, Muntean A, et al. p38-mediated phosphorylation at T367 induces EZH2 cytoplasmic localization to promote breast cancer metastasis. Nat Commun. 2018;9(1):2801.

    PubMed  PubMed Central  Google Scholar 

  17. Wasenang W, Puapairoj A, Settasatian C, Proungvitaya S, Limpaiboon T. Overexpression of polycomb repressive complex 2 key components EZH2/SUZ12/EED as an unfavorable prognostic marker in cholangiocarcinoma. Pathol Res Pract. 2019;215(7):152451.

    CAS  PubMed  Google Scholar 

  18. Venkatesan N, Wong JF, Tan KP, Chung HH, Yau YH, Cukuroglu E, Allahverdi A, Nordenskiold L, Goke J, Geifman-Shochat S, et al. EZH2 promotes neoplastic transformation through VAV interaction-dependent extranuclear mechanisms. Oncogene. 2018;37(4):461–77.

    CAS  PubMed  Google Scholar 

  19. Tripathi BK, AnDerman MF, Bhargava D, Boccuzzi L, Qian X, Wang D, Durkin ME, Papageorge AG, de Miguel FJ, Politi K, et al. Inhibition of cytoplasmic EZH2 induces antitumor activity through stabilization of the DLC1 tumor suppressor protein. Nat Commun. 2021;12(1):6941.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Sun S, Wu Y, Guo W, Yu F, Kong L, Ren Y, Wang Y, Yao X, Jing C, Zhang C, et al. STAT3/HOTAIR Signaling Axis regulates HNSCC Growth in an EZH2-dependent manner. Clin Cancer Res. 2018;24(11):2665–77.

    CAS  PubMed  Google Scholar 

  21. Kim KH, Roberts CW. Targeting EZH2 in cancer. Nat Med. 2016;22(2):128–34.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Yin H, Wang Y, Wu Y, Zhang X, Zhang X, Liu J, Wang T, Fan J, Sun J, Yang A, et al. EZH2-mediated epigenetic silencing of miR-29/miR-30 targets LOXL4 and contributes to Tumorigenesis, Metastasis, and Immune Microenvironment remodeling in breast Cancer. Theranostics. 2020;10(19):8494–512.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Wan L, Xu K, Wei Y, Zhang J, Han T, Fry C, Zhang Z, Wang YV, Huang L, Yuan M, et al. Phosphorylation of EZH2 by AMPK suppresses PRC2 methyltransferase activity and oncogenic function. Mol Cell. 2018;69(2):279–91. e275.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Li B, Yan J, Phyu T, Fan S, Chung TH, Mustafa N, Lin B, Wang L, Eichhorn PJA, Goh BC, et al. MELK mediates the stability of EZH2 through site-specific phosphorylation in extranodal natural killer/T-cell lymphoma. Blood. 2019;134(23):2046–58.

    PubMed  Google Scholar 

  25. Han Y, Wei Y, Yao J, Chu YY, Li CW, Hsu JL, Nie L, Hung MC. Inhibition of CDK2 reduces EZH2 phosphorylation and reactivates ERalpha expression in high-grade serous ovarian carcinoma. Am J Cancer Res. 2020;10(4):1194–206.

    PubMed  PubMed Central  Google Scholar 

  26. Nie L, Wei Y, Zhang F, Hsu YH, Chan LC, Xia W, Ke B, Zhu C, Deng R, Tang J, et al. CDK2-mediated site-specific phosphorylation of EZH2 drives and maintains triple-negative breast cancer. Nat Commun. 2019;10(1):5114.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Ozes AR, Pulliam N, Ertosun MG, Yilmaz O, Tang J, Copuroglu E, Matei D, Ozes ON, Nephew KP. Protein kinase A-mediated phosphorylation regulates STAT3 activation and oncogenic EZH2 activity. Oncogene. 2018;37(26):3589–600.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Wang SQ, Liu J, Qin J, Zhu Y, Tin VP, Yam JWP, Wong MP, Xiao ZJ. CAMK2A supported tumor initiating cells of lung adenocarcinoma by upregulating SOX2 through EZH2 phosphorylation. Cell Death Dis. 2020;11(6):410.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Wei Y, Chen YH, Li LY, Lang J, Yeh SP, Shi B, Yang CC, Yang JY, Lin CY, Lai CC, et al. CDK1-dependent phosphorylation of EZH2 suppresses methylation of H3K27 and promotes osteogenic differentiation of human mesenchymal stem cells. Nat Cell Biol. 2011;13(1):87–94.

    CAS  PubMed  Google Scholar 

  30. Chen S, Bohrer LR, Rai AN, Pan Y, Gan L, Zhou X, Bagchi A, Simon JA, Huang H. Cyclin-dependent kinases regulate epigenetic gene silencing through phosphorylation of EZH2. Nat Cell Biol. 2010;12(11):1108–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Wu SC, Zhang Y. Cyclin-dependent kinase 1 (CDK1)-mediated phosphorylation of enhancer of zeste 2 (Ezh2) regulates its stability. J Biol Chem. 2011;286(32):28511–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Li Z, Hou P, Fan D, Dong M, Ma M, Li H, Yao R, Li Y, Wang G, Geng P, et al. The degradation of EZH2 mediated by lncRNA ANCR attenuated the invasion and metastasis of breast cancer. Cell Death Differ. 2017;24(1):59–71.

    PubMed  Google Scholar 

  33. Feng H, Yu Z, Tian Y, Lee YY, Li MS, Go MY, Cheung YS, Lai PB, Chan AM, To KF, et al. A CCRK-EZH2 epigenetic circuitry drives hepatocarcinogenesis and associates with tumor recurrence and poor survival of patients. J Hepatol. 2015;62(5):1100–11.

    CAS  PubMed  Google Scholar 

  34. Li F, Zeng Z, Xing S, Gullicksrud JA, Shan Q, Choi J, Badovinac VP, Crotty S, Peng W, Xue HH. Ezh2 programs T(FH) differentiation by integrating phosphorylation-dependent activation of Bcl6 and polycomb-dependent repression of p19Arf. Nat Commun. 2018;9(1):5452.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Ning S, Ma X. Dephosphorylation-induced EZH2 activation mediated RECK downregulation by ERK1/2 signaling. J Cell Physiol. 2019;234(10):19010–8.

    CAS  PubMed  Google Scholar 

  36. Zhang L, Wang L, Hu XB, Hou M, Xiao Y, Xiang JW, Xie J, Chen ZG, Yang TH, Nie Q, et al. MYPT1/PP1-Mediated EZH2 dephosphorylation at S21 promotes epithelial-mesenchymal transition in fibrosis through control of multiple families of genes. Adv Sci (Weinh). 2022;9(14):e2105539.

    PubMed  Google Scholar 

  37. Bao Y, Oguz G, Lee WC, Lee PL, Ghosh K, Li J, Wang P, Lobie PE, Ehmsen S, Ditzel HJ, et al. EZH2-mediated PP2A inactivation confers resistance to HER2-targeted breast cancer therapy. Nat Commun. 2020;11(1):5878.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Ong LT, Lee WC, Ma S, Oguz G, Niu Z, Bao Y, Yusuf M, Lee PL, Goh JY, Wang P, et al. IFI16-dependent STING signaling is a crucial regulator of anti-HER2 immune response in HER2 + breast cancer. Proc Natl Acad Sci U S A. 2022;119(31):e2201376119.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.

    CAS  PubMed  Google Scholar 

  43. Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson A, Kampf C, Sjostedt E, Asplund A et al. Proteomics. Tissue-based map of the human proteome. Science 2015, 347(6220):1260419.

  44. Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, Alm T, Asplund A, Bjork L, Breckels LM et al. A subcellular map of the human proteome. Science 2017, 356(6340).

Download references

Acknowledgements

We thank Dr. Chenghui Yang for the Bioinformatics assay of the Supplementary part.

Funding

This work was supported by grants from Zhejiang Provincial Natural Science Foundation of China under Grant (LQ19H160043, to Yu), the National Natural Science Foundation of China (81802633, to Sun; 82302901, to Li), China Postdoctoral Science Foundation Grant (2019M652115, to Sun) and Natural Science Foundation of Zhejiang Province (LQ20H160051, to Zhang).

Author information

Authors and Affiliations

Authors

Contributions

All authors made substantial contributions to the intellectual content of this work via conception, design, data acquisition, analysis, writing and/or interpretation the manuscript. All authors provided their final approval of the manuscript.

Corresponding author

Correspondence to Shanshan Sun.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Second Affiliated Hospital of Zhejiang University School of Medicine and was conducted in accordance with the Declaration of Helsinki. Written informed consent for participation was obtained from the patient in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original version of this article was revised: The author affiliations have been corrected.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1: Supplementary Fig. 1

. EZH2-related genes in tumors provided by bioinformatics analysis based on TCGA-BRCA using STRING, and metascape database

Supplementary Material 2: Supplementary Fig. 2

. EZH2 downstream target genes in breast cancer cell examined by RNA sequencing assay and identified using GO enrichment and KEGG pathway, shown as gene correlation network

Supplementary Material 3

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, F., Li, L., Zhang, M. et al. Phosphorylation of EZH2 differs HER2-positive breast cancer invasiveness in a site-specific manner. BMC Cancer 23, 948 (2023). https://doi.org/10.1186/s12885-023-11450-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11450-9

Keywords