Skip to main content
Download PDF

Analysis on EZH2: mechanism identification of related CeRNA and its immunoassay in hepatocellular carcinoma

BMC Medical Genomics volume 16, Article number: 201 (2023) Cite this article

Abstract

Objective

To screen the possible potential signaling pathways related to enhancer of zeste homolog 2 (EZH2) based on ceRNA mechanism, and to analyze the correlation between E2H2 and depths of various immune cell infiltration depths. The relationship between different immune checkpoints were also analyzed.

Methods

First, the expression of EZH2 in pan-cancer (18 malignancies) was analyzed with the TCGA database. Hepatocellular carcinoma (HCC) tissues of 374 cases and normal tissues of 50 cases were analyzed in terms of the differential expression, overall survival (OS) and progression-free-survival (PFS). Then, we conducted GO and KEGG enrichment analysis on target gene. We also analyzed mRNA-miRNA and MicroRNA (miRNA)- long non-coding RNA (lncRNA) correlation with starbase databse, so as to determine the potential ceRNA mechanism associated with EZH2. Finally, immunoassay and drug-sensitivity analysis of EZH2 was performed.

Results

Seven potential EZH2-related ceRNA pathways were screened out, namely lncRNA: Small Nucleolar RNA Host Gene 1 (SNHG1), SNHG 3, and SNHG 6-miR-101-3p-EZH2; and lncRNA: Long Intergenic Non-Protein Coding RNA 1978 (LINC01978), SNHG12, Ring Finger Protein 216 Pseudogene 1 (RNF216P1), and Coiled-coil Domain Containing 18 Antisense RNA 1 (CCDC18-AS1)-let-7c-5p—EZH2. Finally, 4 potential EZH2-related ceRNA pathways were identified through qPCR.According to immune correlation analysis, EZH2 may be positively correlated with T cells follicular helper, T cells Cluster of differentiation (CD)4 memory activated, Macrophages M0, and B cells memory (P < 0.05, cof > 0.2); while be negatively correlated with T cells CD4 + memory resting (P < 0.05, cof < -0.2). And EZH2 is positively correlated with Programmed Cell Death 1 (PDCD1) (R = 0.22), CD274 (R = 0.3) and Cytotoxic T-Lymphocyte Associated Protein 4 (CTLA4) (R = 0.23). According to drug sensitivity analysis, patients in the high expression group were more susceptible to the effects of various drugs including Sorafenib, 5-Fluorouracil, Doxorubicin, Etoposide, Paclitaxel, and Vinorelbine than those with low expression.

Conclusion

This study revealed seven potential pathways of Enhancer of Zeste Homolog 2 (EZH2)-related ceRNA mechanisms: lncRNA (SNHG3, 6) -Mir-101-3P-ezh2; lncRNA (SNHG12, RNF216P1)-let-7c-5p—EZH2. We also analyzed the immunity and drug sensitivity of EZH2. Our study proves that EZH2 still has great research prospects in HCC.

Peer Review reports

Introduction

Primary liver cancer is the third leading cause of cancer death, and the incidence of HCC accounts for 80% of primary liver cancer. The major pathogenic causes of HCC are hepatitis (e.g., hepatitis B virus [HBV], hepatitis C virus [HCV], etc.), alcoholism, smoking, obesity, and congenital inheritance. There were approximately 906,000 new cases of liver cancer and 830,000 deaths worldwide in 2020 [1]. The advent of immunosuppressant has extended the window period for effective HCC treatment to a certain extent, but this does not mean the ultimate solution for HCC treatment. Perhaps only by further exploring effective biological therapeutic targets and on this basis combined with research in more tumor-related fields, including immunology, is the right way out for oncology research.

LncRNAs have been misunderstood over the past few decades. Since hardly no lncRNAs are expressed as proteins, they were once considered junk. It was not until this century that the role of lncRNA was gradually explored by researchers and became particularly popular in the recent decade. Even though lncRNAs do not express proteins at last, they still play an important role in many life activities such as Dosage compensation effect, epigenetic regulation, cell cycle regulation and cell differentiation regulation [2]. Abnormal regulation of the above processes will cause changes in proto-oncogenes and tumor suppressor genes, and eventually lead to the “alienation” of normal tissues.

As the core component of the Polycomb repression complex 2 (PRC2), EZH2 conventionally catalyzes H3K27me3 as a histone methyltransferase, thereby repressing target gene transcription [3]. EZH2 is highly expressed in most solid tumors, including HCC, and plays a critical role in cell cycle regulation [4], PD-L1 expression [5], and chemo sensitivity [6]. This also stimulates our interest in further study of EZH2. Therefore, we took EZH2 as the research target, and analyzed miRNAs and lncRNAs that potentially correlated with ceRNA. The correlation between EZH2 and immune infiltration levels of related immune cells was also analyzed. We hope to find new therapeutic targets to provide ideas and references for future research. Our research workflow is shown below (Fig. 1).

Fig. 1
figure 1

Research workflow

Full size image

Methods

EZH2-related differential analysis, prognostic analysis and protein expression

The R language package including the ggpubr was utilized to analyze the differential expression of E2F1 in 18 malignant tumors and plot the differential expression of EZH2 in pan-cancer. According to the magnitude of the P value of EZH2 in different cancer types, we divided samples into groups with P < 0.001, P < 0.01, and P < 0.05. The differential expression analysis of the target gene EZH2 was analyzed and visualized by the “limma”, “ggplot2” and “ggpubr” packages, and the prognosis-related survival curves were drawn by the “survival” and “survminer” packages. (P < 0.05 was defined as significant statistical difference.) Protein expression data were obtained through the HPA online database (https://www.proteinatlas.org/).

Analysis of prognostic factors

The correlation between target gene expression and each prognostic factor was analyzed and visualized by “limma” and “ggpubr” packages. Heatmaps for each clinical prognostic factor were drawn with the “limma” and “ComplexHeatmap” packages. ROC curve, calibration curve and nomogram were obtained by “survival”, “survminer”, “timeROC”, “regplot”, and “rms” analysis. Univariate and multivariate COX analysis and forest plots were obtained by the “survival” package. P < 0.05 was defined as statistically significant difference.

GO, KEGG and GSEA analysis

GO, KEGG and GSEA enrichment analysis were conducted and visualized by R language programs “clusterProfiler”, “org.Hs.eg.db”, “enrichplot”, “ggplot2”, “circlize”, “RColorBrewer”, “dplyr”, and “ComplexHeatmap” with 0.05 as P value filter. P < 0.05 was defined as statistically significant difference.

Construction of potential signaling pathways based on ceRNA mechanism and analysis of miRNA and lncRNA survival

The mRNA-miRNA, miRNA-LncRNA interaction data were downloaded from the starbase database (https://rnasysu.com/encori/index.php). Among them, programNum  ≥ 2 was used as one of the mRNA-miRNA screening criteria. Correlation coefficient plot (R > 0.2 was defined as positive correlation, while R < -0.2 as negative correlation), differential expression (P < 0.01 was defined as statistical difference) and survival curve (P < 0.05 was defined as statistical difference) were screened and plotted by the R language package including “limma”, “reshape2”, “ggpubr”, “ggExtra”, “survival” and “survminer”. The diagram of potential EZH2-related ceRNA mechanism was drawned with BioRender software.

Correlation analysis of EZH2 and immune cell infiltration

The correlation between EZH2 and various immune cells, immune checkpoints and immunotherapy analysis were visualized by R language packages including “limma”, “reshape2”, “ggplot2”, “ggpubr”, “vioplot”, “ggExtra” and “corrplot”. The data of immunotherapy analysis were obtained from the TCIA database (http://tcia.at/). The analysis of Programmed Cell Death Protein 1 (PDCD1), CD274 and CTLA4 that are related with EZH2 was conducted by the TIMER 2.0 GEPIA database (http://gepia.cancer-pku.cn), and the P value was calculated by Spearman statistical method. Among them, positive correlation was defined as P < 0.05 and Rp > 0.2; negative correlation was defined as P < 0.05 and R < -0.2; and P > 0.05 was defined as no significance.

Drug sensitivity evaluation

We calculated the half maximal inhibitory concentration (IC50) of drugs by the R package with “pRRophetic” and its dependencies “car, ridge preprocessCore, genefilter, and sva”. The package included the information on the effects of 138 drugs. Boxplots were drawn using the R package “ggplot2”, and P < 0.05 was defined as statistically significant difference.

RNA extraction, cDNA synthesis, and quantitative real-time PCR

RNA extraction and cDNA synthesis were performed using Trizol reagent (TaKaRa, Otsu, Japan) and the PrimeScript RT Reagent Kit with gDNA Eraser (Perfect Real Time) (TaKaRa, Otsu, Japan). Quantitative real-time PCR was performed using SYBR® Fast qPCR Mix (TAKARA), The relative expression of gene was normalized to GAPDH or U6 using 2–ΔΔCt methods, respectively.

Results

Expression of EZH2 in pan-cancer and EZH2- related survival analysis

With the data of various cancer types downloaded from TCGA through UCSC Xena, which were analyzed by R language in terms of the expression of EZH2 in pan-cancer, it was found that EZH2 of 17 malignant tumor tissues (including HCC) was significantly differentially expressed compared with normal tissues (P < 0.001), as shown in Fig. 2a. Then, with the R package, we analyzed the differential expressions between HCC tumor group (374 cases) and normal group (50 cases), as shown in Fig. 2b and c. We also analyzed the disparity between EZH2 low-risk group and high-risk group in terms of disease free survival (DFS) and overall survival (OS). The results showed that EZH2 was differentially expressed between the HCC group and normal groups (TCGA normal and GTEx data) (P < 0.01), and the EZH2 low-risk group was significantly outperformed high-risk group in terms of both DFS (P < 0.001) and OS (P < 0.001), as shown in Fig. 2d and e. In addition, according to immunohistochemical staining test, EZH2 is significantly more stained in HCC as shown in Human Protein Atlas (HPA) (Fig. 2f).

Fig. 2
figure 2

a Differential expression of EZH2 in pan-cancer between tumor groups and normal groups b Differential expression of EZH2 in HCC between tumor group and normal group c Pairwise differential analysis of EZH2 in HCC d Overall survival curve e Overall survival curve f Protein expression between tumor group and normal group in HPA database (T: tumor, N: normal) (* p < 0.05; ** p < 0.01; *** p < 0.001)

Full size image

Analysis of prognostic factors of EZH2 in HCC

Clinical correlation analysis showed that the expression of EZH2 was significantly different in T1 compared with in T2 and T3 (P < 0.01) (Fig. 3a); there was a significant difference between stageI, stageII, and stageIII (P < 0.01) (Fig. 3b); both grade1 and grade2 were significantly different from the other grades (Fig. 3c). And the classifications of T, stage and grade showed significant difference between the high and low expression groups (P < 0.001) (Fig. 3d). In addition, we drew the receiver operating characteristic (ROC) curve based on the target gene EZH2, which showed that the area under curve (AUC) at 1, 3, and 5 years were 0.724, 0.659, and 0.608 (Fig. 3e) respectively. We then built a gene correlation nomogram to evaluate the possibility that EZH2 can predict survival time in HCC (Fig. 3f), and validated the feasibility of the prediction method by a calibration curve (Fig. 3g). Finally, we concluded that EZH2 expression and HCC stage were independent risk factors for prognosis through univariate and multivariate COX regression analysis (Fig. 3h, i).

Fig. 3
figure 3

Prognostic factors: a T b Stage c Grade d Heat map e ROC curve f Nomogram g Calibration curve h Prognostic related risk factors identified by univariate COX regression analysis i Prognostic related risk factors identified by multivariate COX regression analysis (* p < 0.05; ** p < 0.01; *** p < 0.001)

Full size image

EZH2-related Gene Set Enrichment Analysis (GSEA)

To elucidate the biological function of EZH2 in HCC, we analyzed the DEGs between the EZH2 low- and high-expressing groups according to the median EZH2 expression value (P < 0.001). GO analysis showed that: the top five enriched BPs were mainly related to biological processes such as chromosome segregation and DNA replication; the top five CCs mainly occurred in chromosomal, centromeric region and kinetochore; the top five MFs were mainly related to DNA-dependent ATPase activity, gated channel activity, DNA helicase activity, ion channel activity and catalytic activity, acting on DNA correlation (P < 0.001) (Fig. 4a, c). The GO-related GSEA pathway analysis showed that: the most enriched functions and pathways of patients with high EZH2 expression were immunoglobulin complex; while the low expression of EZH2 was associated with the catabolic process of cellular amino acid, fatty acid, monocarboxylic acid and organic acid (Fig. 3b). Meanwhile, in the KEGG enrichment analysis, the top 5 were mainly enriched in cell cycle, DNA replication, neuroactive ligand-receptor interaction, nicotine addiction and homotogous recombination (P < 0.001) (Fig. 4d). And low EZH2 expression was associated with complement and coagulation cascades, cytochrome P450 mediated drug metabolism, fatty acid metabolism, peroxisome and primary bile acid biosynthesis (Fig. 4e).

Fig. 4
figure 4

a GO enrichment analysis b GESA of GO c Circle diagram of GO enrichment analysis d GESA of KEGG e KEGG enrichment analysis

Full size image

Screening of signaling pathway of EZH2-related ceRNA mechanism

miRNAs that may have mRNA-miRNA interaction with EZH2 were searched out in the starbase database. And along with the search results with the original data of the TCGA database, we analyzed the correlation between these miRNAs and EZH2 (cor < -0.2, LogFC < 0, P < 0.05). Further screening was performed by examining the differential expression (P < 0.05) and survival curves (P < 0.05) of these miRNAs. Finally, two miRNAs were identified that may be involved in the ceRNA mechanism targeting EZH2, namely miR-101-3p and let-7c-5p (Fig. 4a, b). The lncRNAs that were related to the above two miRNAs, namely miR-101-3p and let-7c-5p, were screened out, especially ones that were positively correlated with EZH2 and meanwhile negatively correlated with the above two miRNAs. The screening was conducted with starbase database and by R language. Then the above lncRNAs were further analyzed and screened for their differential expression (P < 0.05) and survival curve (P < 0.05) (Fig. 5). Finally, a total of 7 signaling pathways that may have ceRNA mechanisms targeting EZH2 genes were found, namely lncRNA (SNHG1, SNHG3, SNHG6)-miR-101-3p-EZH2 (Fig. 6a), and lncRNA (LINC01978, SNHG12, RNF216P1, CCDC18-AS1)-let-7c-5p-EZH2. (Fig. 6b), We then verified the target gene, and the results showed that: The expression level of SNHG3 in hepatoma cell line HepG2 was significantly higher than that of normal liver cell line LO2 (P < 0.05), and the expression levels of SNHG6, SNHG12 and RNF216P1 in HepG2 and Huh7 hepatoma cell lines were significantly higher than that of LO2 (P < 0.05). The expression levels of miR-101-3p and let-7c-5p in hepatoma cell lines HepG2 and Huh7 were significantly lower than those of LO2 (P < 0.05). Among the screened pathways, lncRNA (SNHG3, SNHG6)-miR-101-3p-EZH2, and lncRNA (SNHG12, RNF216P1,)-let-7c-5p-EZH2 were preliminatively confirmed. It may be a potential EZH2-related signaling pathway with ceRNA mechanism (Fig. 7). To summarize these potential mechanisms, we drew a conceptual diagram of the EZH2-related ceRNA mechanism (Fig. 8).

Fig. 5
figure 5

Correlation analysis, difference analysis and survival curve between two miRNAs and the target gene EZH2 a let-7c-5p b miR-101-3p

Full size image
Fig. 6
figure 6

a Correlation, difference analysis and survival curve of miR-101-3p- EZH2-related lncRNAs with miRNA and EZH2, respectively; b Correlation, difference analysis and survival curves of let-7c-5p-ezh2-related lncRNAs with miRNA and EZH2, respectively

Full size image
Fig. 7
figure 7

The verification of qPCR results of target gene

Full size image
Fig. 8
figure 8

Diagram of EZH2-related ceRNA mechanism (hypothesis)

Full size image

Analysis on EZH2 immune cells and checkpoint correlation

With the TIMER 2.0 database, we analyzed the correlation of infiltration among different types of immune cells, and the results showed that: EZH2 was probably positively correlated with T cells follicular helper, T cells CD4 memory activated, Macrophages M0, and B cells memory (P < 0.05, cof > 0.2); while it may be negatively correlated with T cells CD4 + memory resting (Fig. 9a, b). We also analyzed the correlation between EZH2 and various immune checkpoints (Fig. 9c). Then another analysis on the three immune checkpoints, namely PDCD1, CD274 and CTLA4, was performed with the GEPIA database. It showed that EZH2 was positively correlated with these three immune checkpoint (R > 0.2), as shown in Fig. 9d. In CTLA4, the results of PD-1 immunotherapy analysis showed that the EZH2 low expression group was significantly different from the high expression group when receiving anti-CTLA4 and anti-PD-1 combined with anti-CTLA4 treatment (P < 0.05) (Fig. 9e).

Fig. 9
figure 9

Analysis of the correlation between EZH2 and various immune cells. a lollipop chart b Boxplot c Correlation analysis between EZH2 and various immune checkpoints d Correlation of EZH2 with PDCD1, CD274 and CTLA4 in GEPIA database e The difference violin plot of PD-1, CTLA4 and the combined regimen in EZH2 high and low expression groups (* p < 0.05; ** p < 0.01; *** p < 0.001)

Full size image

To further investigate the possibility of individualized treatment of EZH2 in HCC patients, we examined the relationship between risk score and IC50 of various drugs used in clinical treatment of HCC, including Sorafenib, 5-Fluorouracil, Doxorubicin, Etoposide, Paclitaxel, and Vinorelbine. In addition, we analyzed and listed the sensitivity of potential therapeutic agents. It showed that patients in the high-expression group appeared to be more susceptible to most drugs than those with low expression (Fig. 10).

Fig. 10
figure 10

The ic50 of EZH2 high and low expression groups in various HCC treatment drugs (P < 0.001)

Full size image

Discussion

As a member of the family of polycomb group genes (PcGs) [3], EZH2 is located on chromosome 7q35, covering a total of 20 exons and encoding 746 amino acid residues. It is composed of five domains, including 1) EED-interaction domain (EID); 2) Domain I and Domain II, 3) cysteine-rich domain (CXC domain); 4) C-terminal suppressor of variegation 39; and 5) enhancer of zeste and trithorax domain (SET domain). The histone methyltransferase activity of EZH2 is mainly maintained by the SET domain, which is also the central hub of the EZH2-related PRC2-dependent H3K27 methylation action mode. Gene silencing is mediated by PRC2, which is one of the PcG protein core complexes, mainly by regulating chromatin structure [7]. Under normal conditions, EZH2 mainly functions through three basic mechanisms: PRC2-dependent H3K27 methylation [3, 8], PRC2-dependent non-histone protein methylation [9] and PRC2-independent gene transactivation [10]. However, the abnormal expressed EZH2 will lead to the occurrence, development, invasion and metastasis of tumors. In recent years, a series of adverse chain reactions of malignant tumor progression due to EZH2 overexpression were reported more frequently, including malignant tumors like HCC [11], gastric cancer [12], breast cancer [13], esophageal cancer [14] and endometrial carcinoma [15]. Recent reports on HCC found a strong positive correlation between EZH2 and cycle-related genes such as CCND1, CDK2, CDK4, CCNB1, and CCNB2 [4]. And regulating the genes related to the above cells’ cycle results in HCC proliferation. Not only that, studies have shown that EZH2 can also inhibit the expression of the immune checkpoint inhibitor PD-L1 by directly upregulating the levels of CD274 and IFN regulatory factor 1 (IRF1) promoter H3K27me3 in liver cancer cells [5]. In addition, EZH2 also represses the tumor suppressor gene Atonal bHLH transcription factor 8 (ATOH8) by controlling DNA methylation and H3K27 methylation of its promoter, which in turn increases stem cell properties and reduces HCC chemo sensitivity [6]. The above findings all provide an important rationale for EZH2 as a therapeutic target for HCC, which motivates our in-depth research.

In terms of EZH2-related immune analysis, based on the TIMER database, different algorithms were chosen for the correlation of immune infiltration levels against different types of immune cells, so as to ensure the greatest accuracy and effectiveness. The results showed that: EZH2 may be positively correlated with T cells follicular helper, T cells CD4 memory activated, Macrophages M0, and B cells memory (P < 0.05, cof > 0.2); while EZH2 may be negatively correlated with T cells CD4 + memory resting (P < 0.05, cof < -0.2). Previous studies found that activated CD4 + T cells from HCCs stimulate CXCL10 production by macrophages, CXCL10 bind CXC chemokine receptor 3 on B cells and cause them to become IgG-producing plasma cells. IgG activated Fc receptor in macrophages to produce cytokines that reduce the anti-tumor immune response. Activation of CD4 + T cells in HCC stimulated macrophage production of CXCL10, and in turn, bind CXC chemokine receptor 3 on B cells and made them IgG-producing plasma cells. IgG activated Fc receptors in macrophages to produce cytokines that reduced antitumor immune responses [16]. This pathway witnessed the increased expression of DNA methyltransferase 1 and EZH2, which resulted from HCC and hepatoma cells. And reports about T cells follicular helper showed that EZH2 instructed the early commitment to TFH cell differentiation by stabilizing a cluster of TFH cell lineage-associated genes, including Bcl6 [17]. It was also confirmed the role of EZH2 in governing TFH cell differentiation by integrating phosphorylation-dependent Bcl6 activation and H3K27me3-dependent repression of p19Arf. Therefore, EZH2 regulates TH differentiation in a cell-type-specific manner [18]. In addition, a study of lymphoma pointed out that dysregulation of the GC, which result from constitutively active EZH2, activates lymphoma formation and identifies EZH2 as a possible therapeutic target for NHL and other GC-derived B-cell diseases. EZH2 maintains activation-induced cytidine deaminase function and prevents terminal differentiation of B cells in the germinal center, resulting in antibody diversification and affinity maturation. Dysregulation of the GC response by constitutively active EZH2 promotes lymphoma formation [19]. In addition to the above-mentioned reports on immune cells that were significantly different in our analysis, it was also found that phenotypic analysis of human CD8 + EZH2 + cells confirmed its greater effector capacity and the ability to reduce apoptosis sensitivity [20]. In the melanoma B16-tumor-bearing mice model of transferred Ezh2fl/flCd4Cre Pmel CD8 + T cells, the tumor suppressor effect of these cells was significantly reduced when EZH2 was insufficient [21]. EZH2 decreased CD8 + T cell suppression by inhibiting Tregs function. For instance, it promoted the frequency of CD8 + T cells producing ifn γ and reduced inhibitory marker expression on tumor-infiltrating T cells [22]. It has been reported that CAFs inhibited the expression of VASH1 by affecting the expression of VEGF and then upregulating the expression of EZH2 in HCC, thereby promoting the proliferation and angiogenesis of human umbilical vein endothelial cells (HUVECs). And this phenomenon has been further confirmed in vivo and in vitro experiments [23]. Furthermore, genetic deletion of EZH2 increased the expression of IL-15R, CD122 and NKG2D-activating receptors, resulting in enhanced NK cell production in hematopoietic stem and progenitor cells (HSPCs). And inhibition of EZH2 also improved mature NK cell function [24].

In this study, based on bioinformatics, we analyzed lncRNAs and miRNAs with potential ceRNA mechanisms targeting EZH2, namely miR-101-3p and three lncRNAs (SNHG1, 3, 6) that may interact with it, as well as hsa-let-7c-5p and 4 lncRNAs (LINC01978, SNHG12, RNF216P1, CCDC18-AS1) that may interact with it. Among the above ncRNAs, there are many miRNAs and lncRNAs that have been found closely related to the occurrence and development of tumors. In terms of miRNAs, a recent report showed that the expression of miR-101-3p in plasma exosomes from medulloblastoma (MB) patients was significantly higher than that of healthy controls, and these exosomal miRNAs could be delivered to tumor cells via exosomes. According to subsequent in vitro functional analysis, treatment of MB cells with the corresponding mimics of miR-101-3p significantly inhibited proliferation, colony formation, migration and invasion of tumor cells, and promoted apoptosis. Meanwhile, miR- 101-3p also targeted EZH2, thereby enhancing tumor suppression [25]. And it has been confirmed in several researches that SNHG6 inhibited apoptosis by regulating EZH2 expression via the sponging of MiR-101-3p in esophageal squamous-cell carcinoma [26, 27]. In addition, miR-101-3p also had coordinated or independent effects on the occurrence and development of malignant tumors such as breast cancer, [28] renal cell carcinoma [29] and oral squamous cell carcinoma [30]. Studies have shown that let-7c-5p enhanced 5-FU exposure (by inhibiting ABCC5/MRP5 expression) and co-targeted thymidylate synthase with 5-FU (let-7c reduced protein expression, 5-FU made the enzyme irreversibly inactivated). And this finally resulted in a strong synergistic effect in inhibiting the viability of hepatoma cells [31]. In terms of LncRNAs, SNHG1 repressed the transcription of Cyclin Dependent Kinase Inhibitor 1A (CDKN1A) and CDKN2B in the nucleus by enhancing EZH2-mediated H3K27me3 in the promoters of CDKN1A and CDKN2B, leading to de-repression of the cell cycle, thereby promoting HCC cell growth, cell cycle progression, metastasis and epithelial–mesenchymal transition (EMT) [32]. It has also been reported that SNHG3 promoted CD151 by inhibiting miR-128 signaling, resulting in HCC invasion, EMT, and sorafenib resistance [33]. SNHG12 regulated the SNHG12/miR-199a/b-5p pathway through a ceRNA mechanism, thereby reducing the inhibitory effect on Mixed Lineage Kinase 3 (MLK3, its expression was upregulated in tumor tissues and was associated with tumor progression and metastasis), thus enhanced expression of MLK3 and its downstream effectors in the NF-κB pathway, and finally promoted the occurrence and development of HCC [34]. Recent research reports suggest that CCDC18-AS1 was a lncRNA involved in pyroptosis which was related to the immune landscape of COAD patients [35].

In summary, 4 EZH2-related potential signaling pathways were screened out based on the ceRNA mechanism in multiple public databases and combined with R language analysis. In addition, we analyzed EZH2 drug sensitivity and the correlations between EZH2 and various immune cells as well as immune checkpoints were analyzed, which provided clear ideas and directions for future basic experimental research. However, it should be pointed out that although the above-mentioned signaling pathways and immune-related analysis have been verified in the previous literature, some of the results reached in this article were still in the stage of research and verification. We hope that this article can provide new research directions and biological targets for future HCC research.

Conclusion

In conclusion, this study revealed four potential pathways of EZH2-related ceRNA mechanisms: lncRNA [SNHG3, SNHG6] -Mir-101-3P-ezh2; lncRNA (SNHG12, RNF216P1)-let-7c-5p—EZH2. We also analyzed the immunity and drug sensitivity of EZH2. At present, EZH2 still has great potential and prospect for the treatment of HCC. By further understanding its functional scope, we expect to find more potential therapeutic and diagnostic options for EZH2 in HCC.

Availability of data and materials

The raw data used in this experiment come from TCGA website (https://portal.gdc.cancer.gov/), starbase (https://rnasysu.com/encori/index.php), TIMER2.0 (http://timer.cistrome.org/), TCIA database (http://tcia.at/), KEGG pathway database (www.kegg.jp/kegg/kegg1.html), HPA online database (https://www.proteinatlas.org/) and GEPIA (http://gepia.cancer-pku.cn) (The relation of the three immune checkpoints associated with EZH2 by clicking on the link to the GEPIA home page, typing “EZH2” in the box, and then clicking “GoPIA”. Then click “correlation” and enter “EZH2” in “Gene A” and “CTLA4”, “PDCD1” and “CD274” respectively in “GENE B”. Then select “LIHC Tumor” in “TCGA Tumor” and finally click “Plot”). We would like to express our sincere thanks to the collectors and producers of the research data.

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660. Epub 2021 Feb 4. PMID: 33538338.

    Article  PubMed  Google Scholar 

  2. Ferrè F, Colantoni A, Helmer-Citterich M. Revealing protein-lncRNA interaction. Brief Bioinform. 2016;17(1):106–16. https://doi.org/10.1093/bib/bbv031. Epub 2015 Jun 2. PMID: 26041786; PMCID: PMC4719072.

    Article  CAS  PubMed  Google Scholar 

  3. Duan R, Du W, Guo W. EZH2: a novel target for cancer treatment. J Hematol Oncol. 2020;13(1):104. https://doi.org/10.1186/s13045-020-00937-8. PMID:32723346;PMCID:PMC7385862.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Zhang DY, Sun QC, Zou XJ, Song Y, Li WW, Guo ZQ, Liu SS, Liu L, Wu DH. Long noncoding RNA UPK1A-AS1 indicates poor prognosis of hepatocellular carcinoma and promotes cell proliferation through interaction with EZH2. J Exp Clin Cancer Res. 2020;39(1):229. https://doi.org/10.1186/s13046-020-01748-y. PMID:33121524;PMCID:PMC7596966.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Xiao G, Jin LL, Liu CQ, Wang YC, Meng YM, Zhou ZG, Chen J, Yu XJ, Zhang YJ, Xu J, Zheng L. EZH2 negatively regulates PD-L1 expression in hepatocellular carcinoma. J Immunother Cancer. 2019;7(1):300. https://doi.org/10.1186/s40425-019-0784-9. PMID:31727135;PMCID:PMC6854886.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Zhang Y, Tang B, Song J, Yu S, Li Y, Su H, He S. Lnc-PDZD7 contributes to stemness properties and chemosensitivity in hepatocellular carcinoma through EZH2-mediated ATOH8 transcriptional repression. J Exp Clin Cancer Res. 2019;38(1):92. https://doi.org/10.1186/s13046-019-1106-2. PMID:30786928;PMCID:PMC6381703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Pasini D, Di Croce L. Emerging roles for Polycomb proteins in cancer. Curr Opin Genet Dev. 2016;36:50–8. https://doi.org/10.1016/j.gde.2016.03.013. Epub 2016 May 2. PMID: 27151431.

    Article  CAS  PubMed  Google Scholar 

  8. Wang R, Xin M, Li Y, Zhang P, Zhang M. The functions of histone modification enzymes in cancer. Curr Protein Pept Sci. 2016;17(5):438–45.

    Article  CAS  PubMed  Google Scholar 

  9. He A, Shen X, Ma Q, Cao J, von Gise A, Zhou P, Wang G, Marquez VE, Orkin SH, Pu WT. PRC2 directly methylates GATA4 and represses its transcriptional activity. Genes Dev. 2012;26(1):37–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Xu L, Lin J, Deng W, Luo W, Huang Y, Liu CQ, Zhang FP, Qin YF, Wong PP, Liu C. EZH2 facilitates BMI1-dependent hepatocarcinogenesis through epigenetically silencing microRNA-200c. Oncogenesis. 2020;9(11):101. https://doi.org/10.1038/s41389-020-00284-w. PMID:33168810;PMCID:PMC7652937.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Pan YM, Wang CG, Zhu M, Xing R, Cui JT, Li WM, Yu DD, Wang SB, Zhu W, Ye YJ, Wu Y, Wang S, Lu YY. STAT3 signaling drives EZH2 transcriptional activation and mediates poor prognosis in gastric cancer. Mol Cancer. 2016;15(1):79. https://doi.org/10.1186/s12943-016-0561-z. PMID:27938379;PMCID:PMC5148878.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Yin H, Wang Y, Wu Y, Zhang X, Zhang X, Liu J, Wang T, Fan J, Sun J, Yang A, Zhang R. 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. https://doi.org/10.7150/thno.44849. PMID:32754259;PMCID:PMC7392008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Qiu BQ, Lin XH, Ye XD, et al. Long non-coding RNA PSMA3-AS1 promotes malignant phenotypes of esophageal cancer by modulating the miR-101/EZH2 axis as a ceRNA. Aging (Albany NY). 2020;12(2):1843–56.

    Article  CAS  PubMed  Google Scholar 

  15. Zhang C, Shao S, Zhang Y, Wang L, Liu J, Fang F, Li P, Wang B. LncRNA PCAT1 promotes metastasis of endometrial carcinoma through epigenetical downregulation of E-cadherin associated with methyltransferase EZH2. Life Sci. 2020;243:117295. https://doi.org/10.1016/j.lfs.2020.117295. Epub 2020 Jan 9. PMID: 31927050.

    Article  CAS  PubMed  Google Scholar 

  16. Wei Y, Lao XM, Xiao X, Wang XY, Wu ZJ, Zeng QH, Wu CY, Wu RQ, Chen ZX, Zheng L, Li B, Kuang DM. Plasma cell polarization to the immunoglobulin G phenotype in hepatocellular carcinomas involves epigenetic alterations and promotes hepatoma progression in mice. Gastroenterology. 2019;156(6):1890-1904.e16. https://doi.org/10.1053/j.gastro.2019.01.250. Epub 2019 Jan 31. PMID: 30711627.

    Article  CAS  PubMed  Google Scholar 

  17. Chen X, Cao G, Wu J, Wang X, Pan Z, Gao J, Tian Q, Xu L, Li Z, Hao Y, Huang Q, Wang P, Xiao M, Xie L, Tang S, Liu Z, Hu L, Tang J, He R, Wang L, Zhou X, Wu Y, Chen M, Sun B, Zhu B, Huang J, Ye L. The histone methyltransferase EZH2 primes the early differentiation of follicular helper T cells during acute viral infection. Cell Mol Immunol. 2020;17(3):247–60. https://doi.org/10.1038/s41423-019-0219-z. Epub 2019 Mar 6. PMID: 30842630; PMCID: PMC7052164.

    Article  CAS  PubMed  Google Scholar 

  18. Li F, et al. Ezh2 programs TFH differentiation by integrating phosphorylation-dependent activation of Bcl6 and polycomb-dependent repression of p19Arf. Nat Commun. 2018;9:5452.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Caganova M, Carrisi C, Varano G, Mainoldi F, Zanardi F, Germain PL, George L, Alberghini F, Ferrarini L, Talukder AK, Ponzoni M, Testa G, Nojima T, Doglioni C, Kitamura D, Toellner KM, Su IH, Casola S. Germinal center dysregulation by histone methyltransferase EZH2 promotes lymphomagenesis. J Clin Invest. 2013;123(12):5009–22. https://doi.org/10.1172/JCI70626. Epub 2013 Nov 8. Erratum in: J Clin Invest. 2014;124(4):1869. PMID: 24200695; PMCID: PMC3859423.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Stairiker CJ, Thomas GD, Salek-Ardakani S. EZH2 as a regulator of CD8+ T cell fate and function. Front Immunol. 2020;11:593203. https://doi.org/10.3389/fimmu.2020.593203. PMID: 33117406; PMCID: PMC7574680.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. He S, Liu Y, Meng L, Sun H, Wang Y, Ji Y, et al. Ezh2 phosphorylation state determines its capacity to maintain CD8(+) T memory precursors for antitumor immunity. Nat Commun. 2017;8:2125. https://doi.org/10.1038/s41467-017-02187-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Goswami S, Apostolou I, Zhang J, Skepner J, Anandhan S, Zhang X, et al. Modulation of EZH2 expression in T cells improves efficacy of anti-CTLA-4 therapy. J Clin Invest. 2018;128:3813–8. https://doi.org/10.1172/JCI99760.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Huang B, Huang M, Li Q. Cancer-Associated Fibroblasts Promote Angiogenesis of Hepatocellular Carcinoma by VEGF-Mediated EZH2/VASH1 Pathway. Technol Cancer Res Treat. 2019;18:1533033819879905. https://doi.org/10.1177/1533033819879905. PMID:31757187;PMCID:PMC6876164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Yin J, Leavenworth JW, Li Y, Luo Q, Xie H, Liu X, Huang S, Yan H, Fu Z, Zhang LY, Zhang L, Hao J, Wu X, Deng X, Roberts CW, Orkin SH, Cantor H, Wang X. Ezh2 regulates differentiation and function of natural killer cells through histone methyltransferase activity. Proc Natl Acad Sci U S A. 2015;112(52):15988–93. https://doi.org/10.1073/pnas.1521740112. Epub 2015 Dec 14. PMID: 26668377; PMCID: PMC4702963.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Xue P, Huang S, Han X, Zhang C, Yang L, Xiao W, Fu J, Li H, Zhou Y. Exosomal miR-101–3p and miR-423–5p inhibit medulloblastoma tumorigenesis through targeting FOXP4 and EZH2. Cell Death Differ. 2022;29(1):82–95. https://doi.org/10.1038/s41418-021-00838-4. Epub 2021 Jul 22. PMID: 34294888; PMCID: PMC8738741.

    Article  CAS  PubMed  Google Scholar 

  26. Wang J, Yang X, Li R, Zhang R, Hu D, Zhang Y, Gao L. LncRNA SNHG6 inhibits apoptosis by regulating EZH2 Expression via the sponging of MiR-101-3p in esophageal squamous-cell carcinoma. Onco Targets Ther. 2020;13:11411–20. https://doi.org/10.2147/OTT.S275135. PMID:33192074;PMCID:PMC7656962.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhang Y, Li R, Ding X, He M, Zhang R. Long noncoding RNA SNHG6 promotes oesophageal squamous cell carcinoma by downregulating the miR-101–3p/EZH2 pathway. J Biochem Mol Toxicol. 2021;36:e22959. https://doi.org/10.1002/jbt.22959. Epub ahead of print. PMID: 34766670.

    Article  CAS  PubMed  Google Scholar 

  28. Jiang H, Li L, Zhang J, Wan Z, Wang Y, Hou J, Yu Y. MiR-101–3p and Syn-Cal14.1a synergy in suppressing EZH2-Induced progression of breast cancer. Onco Targets Ther. 2020;13:9599–609. https://doi.org/10.2147/OTT.S264600. PMID: 33061442; PMCID: PMC7532305.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Dong Y, Gao Y, Xie T, Liu H, Zhan X, Xu Y. miR-101-3p Serves as a Tumor Suppressor for Renal Cell Carcinoma and Inhibits Its Invasion and Metastasis by Targeting EZH2. Biomed Res Int. 2021;2021:9950749. https://doi.org/10.1155/2021/9950749. PMID:34307682;PMCID:PMC8282380.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wu K, Jiang Y, Zhou W, Zhang B, Li Y, Xie F, Zhang J, Wang X, Yan M, Xu Q, Ren Z, Chen W, Cao W. Long noncoding RNA RC3H2 facilitates cell proliferation and invasion by targeting MicroRNA-101–3p/EZH2 Axis in OSCC. Mol Ther Nucleic Acids. 2020;20:97–110. https://doi.org/10.1016/j.omtn.2020.02.006. Epub 2020 Feb 13. PMID: 32163895; PMCID: PMC7066035.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Jilek JL, Tu MJ, Zhang C, Yu AM. Pharmacokinetic and pharmacodynamic factors contribute to synergism between Let-7c-5p and 5-Fluorouracil in inhibiting hepatocellular carcinoma cell viability. Drug Metab Dispos. 2020;48(12):1257–63. https://doi.org/10.1124/dmd.120.000207. Epub 2020 Oct 13. PMID: 33051247; PMCID: PMC7684025.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Li B, Li A, You Z, Xu J, Zhu S. Epigenetic silencing of CDKN1A and CDKN2B by SNHG1 promotes the cell cycle, migration and epithelial-mesenchymal transition progression of hepatocellular carcinoma. Cell Death Dis. 2020;11(10):823. https://doi.org/10.1038/s41419-020-03031-6. PMID:33009370;PMCID:PMC7532449.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhang PF, Wang F, Wu J, Wu Y, Huang W, Liu D, Huang XY, Zhang XM, Ke AW. LncRNA SNHG3 induces EMT and sorafenib resistance by modulating the miR-128/CD151 pathway in hepatocellular carcinoma. J Cell Physiol. 2019;234(3):2788–94. https://doi.org/10.1002/jcp.27095. Epub 2018 Aug 21. PMID: 30132868.

    Article  CAS  PubMed  Google Scholar 

  34. Lan T, Ma W, Hong Z, Wu L, Chen X, Yuan Y. Long non-coding RNA small nucleolar RNA host gene 12 (SNHG12) promotes tumorigenesis and metastasis by targeting miR-199a/b-5p in hepatocellular carcinoma. J Exp Clin Cancer Res. 2017;36(1):11. https://doi.org/10.1186/s13046-016-0486-9. PMID:28073380;PMCID:PMC5223416.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Xia F, Yan Y, Shen C. A prognostic pyroptosis-related lncRNAs risk model correlates with the immune microenvironment in colon adenocarcinoma. Front Cell Dev Biol. 2021;9:811734. https://doi.org/10.3389/fcell.2021.811734. PMID: 34966747; PMCID: PMC8710686.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

No funding.

Author information

Authors and Affiliations

  1. Department of Hepatobiliary and Pancreatic Surgery, Affiliated Cancer Hospital of Harbin Medical University, Harbin, Heilongjiang Province, 150040, China

    Haoran Zhao, Haishi Liu & Chao Zhan

  2. Department of Oncology, Beidahuang Industry Group General Hospital, No. 235 Hashuang Road, Harbin, Heilongjiang Province, 150088, China

    Wenli Kang, Yingchun Man & Tong Qu

Authors
  1. Haoran Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Haishi Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Wenli Kang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Chao Zhan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yingchun Man
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Tong Qu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

TQ and YM contributed to study design and supervision; CZ contributed to data collection, HZ and HL contributed to data analysis and article writing; TQ and YM revised the manuscript. WK finished the experimental part. All authors read and approved the final manuscript, and all authors reviewed the manuscript.

Corresponding authors

Correspondence to Yingchun Man or Tong Qu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

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.

Supplementary Information

Additional file 1: Supplement 1.

lncRNA-miRNA-mRNA pathway related parameters.

Additional file 2: Supplement 2.

Other drug susceptibility evaluation.

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

Zhao, H., Liu, H., Kang, W. et al. Analysis on EZH2: mechanism identification of related CeRNA and its immunoassay in hepatocellular carcinoma. BMC Med Genomics 16, 201 (2023). https://doi.org/10.1186/s12920-023-01594-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-023-01594-9

Keywords

BMC Medical Genomics

ISSN: 1755-8794

Contact us