- Research
- Open access
- Published:
ZMIZ2 facilitates hepatocellular carcinoma progression via LEF1 mediated activation of Wnt/β-catenin pathway
Experimental Hematology & Oncology volume 13, Article number: 5 (2024)
Abstract
Background
Hepatocellular carcinoma (HCC) is one of the most common malignancies with a high lethality rate. ZMIZ2 is a transcriptional co-activator implicated in various human diseases. However, the role and molecular mechanism of ZMIZ2 in HCC remains to be elucidated.
Methods
The expression and prognostic value of ZMIZ2 in HCC was excavated from public databases and explored by bioinformatic analysis. Then the expression of ZMIZ2 and related genes was further validated by quantitative RT-PCR, western blotting, and immunohistochemistry. Loss and gain-of-function experiments were performed in vitro and in vivo to investigate the function of ZMIZ2 in HCC. In addition, transcriptome sequencing and immunoprecipitation was conducted to explore the potential molecular mechanisms of ZMIZ2.
Results
ZMIZ2 was highly expressed in HCC and associated with poor prognosis. Silencing ZMIZ2 significantly inhibited HCC cell proliferation, cell cycle process, migration, and invasion in vitro, and also inhibited the progression of HCC in vivo. Additionally, ZMIZ2 expression was correlated with immune cell infiltration in HCC samples. Somatic mutation analysis showed that ZMIZ2 and TP53 mutations jointly affected the progression of HCC. Mechanistically, ZMIZ2 interacted with LEF1 to regulate malignant progression of HCC by activating the Wnt/β-catenin pathway.
Conclusion
ZMIZ2 was overexpressed in HCC and associated with poor prognosis. The overexpression of ZMIZ2 was corelated with malignant phenotype, and it facilitated HCC progression via LEF1-mediated activation of the Wnt/β-catenin pathway. Furthermore, ZMIZ2 could be served as a prognostic biomarker and a new therapeutic target for HCC.
Introduction
Hepatocellular carcinoma (HCC) is a worldwide common malignancy with high incidence rate and high mortality [1]. Aggressive is widely considered as the main driving force for HCC metastasis and recurrence [2], and the rapid disease progression caused by tumor recurrence is the primary cause of the mortality rate [3]. In addition, the majority of HCC patients have a poor prognosis due to insufficient early diagnosis and limited effective treatment options for advanced HCC [4, 5], which drives global researchers to reveal deeper into the molecular mechanism of hepatocarcinogenesis and explore new treatment strategies.
HCC is characterized by the decreased genomic stability and accumulation of somatic DNA alterations such as TP53 mutations [6]. Besides, sustained activation of oncogene signaling such as RHO GTPase related signaling [7], HGF/c-MET signaling [2], PI3K/AKT signaling [3], NRAS/ERK signaling [8], is a key driver of hepatocellular carcinoma (HCC) progression. Hyperactivation of the Wnt/β-catenin signaling is an early event in HCC pathogenesis [9, 10]. Approximately 50% of HCC exhibit activation of the Wnt/β-catenin signaling pathway due to mutations in CTNNB1, AXIN1 or APC [11]. Agents such as salinomycin and NVP-TNKS656 targeting the upstream molecules of Wnt/β-catenin signaling have been developed. However, these agents may be ineffective when the pathway is activated downstream. Blocking the interactions between β-catenin and its transcriptional partners with agents such as PRI-724, ICG001, and isoquercitrin is a promising therapeutic approach [12]. Therefore, extensive studies and the ongoing clinical trials highlight the promise for developing novel therapeutic strategies targeting the Wnt/β-catenin pathway.
Zinc finger MIZ type containing 2 (ZMIZ2) is a co-activator of transcription and an inhibitor of activated STAT-like proteins (PIAS) [13]. The protein also possesses an intrinsic transcription activation domain (TAD) that enhances the transcriptional activity of nuclear hormone receptors and other transcription factors [14,15,16,17]. While aberrant ZMIZ2 expression has related to several cancer types, such as prostate cancer [13], breast cancer [18] and colorectal cancer [19], its potential involvement in HCC remains unclear.
Lymphoid enhancer-binding factor-1 (LEF1) is a member of the LEF1/T-cell factor family and serves as a nuclear effector in the Wnt/β-catenin signaling pathway [20]. It is characterized by a highly conserved high-mobility group DNA-binding domain [20, 21]. As a transcription regulatory factor, LEF1 controls the expression of various genes and influences the progression of multiple malignancies [22,23,24,25]. Prior studies have shown that LEF1 promotes cell self-renewal ability, dedifferentiation, invasion, and drug resistance in HCC [22, 26,27,28]. As a key player in the Wnt/β-catenin pathway, LEF1 binds to accumulated β‐catenin in the nucleus and promotes the transcription of Wnt‐responsive genes, which in turn drives cell growth and malignant progression [29]. However, whether LEF1 interacts with ZMIZ2 to modulate the Wnt/β-catenin pathway is yet to be determined.
The TP53 gene is the most commonly mutated gene in cancer [30, 31]. It is reported that TP53 mutational inactivation causes β-catenin accumulation and activates the Wnt signaling pathway [32,33,34]. The relationship between the role of ZMIZ2 and TP53 mutations remains unclear. Transcriptional co-activators, such as YAP/TAZ and CBP/p300, have been implicated in the establishment of tumor microenvironment (TME) and immune escape [35,36,37]. However, the role of ZMIZ2 in HCC tumor immunity deserves further exploration.
In the present study, the expression and clinical significance of ZMIZ2 in HCC were explored by bioinformatics analysis and experimental validation. In particular, loss and gain-of-function study, transcriptome sequencing and co-immunoprecipitation experiments were conducted to explore the significant role of ZMIZ2 in activating the Wnt/β-catenin signaling pathways to promote HCC progression.
Materials and methods
Patients specimens and clinical data collection
Eight pairs of primary liver cancer tissues and adjacent non-tumor tissue samples were obtained from Zhongnan Hospital of Wuhan University. Informed consent was obtained from patients participating in this study. Furthermore, detailed clinical information of these HCC patients was obtained from the electronic medical record system of Zhongnan Hospital. This study was also approved by the ethics committee of Zhongnan Hospital (Approval No.2022011K).
Cell culture
Human immortalized liver cell line of L02 and liver cancer cell lines including Huh7, Hep3B, HCCLM3, HepG2, SK-hep1, MHCC-97 H were were routinely cultured in high glucose DMEM (Gibco, USA) or MEM (Gibco, USA) containing 10% fetal bovine serum (Gibco, USA). All cells were nurtured in a saturated humidity at 37 °C with 5% CO2.
RNA extraction, reverse transcription and quantitative real-time PCR
Total RNA was extracted from cells and tissues using Trizol Reagent (Invitrogen, USA). The Nanodrop 2000 spectrophotometer (ThermoFisher, USA) was used to measure the quantity and quality of RNA. cDNA was produced through a reverse transcription process with the TOYOBO ReverTra Ace Kit (Toyobo, Japan). Gene expression was accomplished using UltraSYBR Mixture (CWbio, China) and the CFX96 Touch Real-Time PCR system (Biorad, USA). The mRNA expression levels of ZMIZ2 and related genes were standardized to GAPDH. The relative quantification was calculated using the 2−△△CT method. Primers designed for RT-PCR in this research are listed in Additional file 2: Table S2.
EDU assay
All transfected with siRNA, plasmid cells were seeded into 96-well plates at 4 × 103 cells/well and incubated for 24 h. The proliferation ability was detected with an EDU assay kit (#C10310-1, RIBO, Guangzhou, China) following the manufacturer’s instructions. Images were acquired using a fluorescence microscope (Olympus, Shanghai, China), and the percentage of EDU positive cells was calculated using the formula EDU-positive cell count/total cell count × 100%.
Cell counting kit 8 assay
Transfected correspondingly siRNA, plasmid or control to Huh7, Hep3B, SK-hep1 and HepG2 cells for 24 h, Thereafter, 2000 and 4000 transfected cells were seeded into separate 96-well plates. Cell viability was assessed with the Cell Counting Kit-8 (CCK8) (#CK04, Dojindo, Japan) according to the manufacturer’s instructions. CCK-8 reagent was added at 0 h,24 h,48 h,72 h,96 h. After 2 h of incubation, absorbance reading was tested on a microtiter plate reader at 450 nm (ELX-800, BioTek, USA). Then the growth curves were plotted according to the OD values.
Luminescent cell viability assay
Cell transfection, counting, and seeding are the same as CCK8. The CellTiter-Glo® Luminescent Cell Viability Assay Kit (#G7570, Promega, USA) was used to determine the intracellular ATP and reflect cell viability. First, 100 µl CellTiter-Glo solution was directly mixed with the culture medium and incubated at room temperature for 10 min. Following this, the intracellular ATP luminescence was recorded using Victor Nivo® multimode plate reader (Perkin-Elmer, Midrand, South Africa). Add reaction solution at 0 h, 24 h, 48 h, 72 h, and 96 h to detect the luminescence signal.
Wound healing, migration and matrigel invasion assays
The 200 µl pipette tip was employed to create wounds when 90% of the tumor cells had merged following transfection. The previous cell culture medium was replaced by a medium containing 2% fetal bovine serum. At 0 h, 24 h, and 48 h, the wounds were noted by a microscope and photographed, respectively. For migration assays, 5 × 104 cells suspended in serum-free medium and 400 µl of cell suspension were seeded into the upper chamber. Then the lower chamber was filled with 600 µl medium containing 10% FBS, which acted as a chemoattractant. After 24 to 48 h, the cells were wiped off the upper chamber (6.5 mm. in diameter, 8.0 μm pore size, Corning, USA) with a cotton swab. The invaded cells to the underside of the membrane were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet (#C0121, Beyotime, Shanghai, China), then counted in three randomly selected microscopic fields (200×) per filter. For Matrigel invasion assays, we paved Matrigel (#356234, Corning, USA) in the surface of upper chamber and 1 × 105 cells were seeded. The remaining steps were same as migration assays.
Flow cytometry analysis
Huh7, Hep3B, SK-hep1 and HepG2 cells were transfected with required siRNA, plasmid, or control, then cell apoptosis assay was performed with the Annexin V-FITC/PI Apoptosis Detection Kit (#4101-2, BestBio, China). The cells were stained with PI to conduct cell cycle analysis using cell cycle staining kit (#CCS012, Multisciences, China) according to the manufacturer’s protocol. The cells were analyzed by flow cytometry (Cytoflex Beckman, China).
Immunohistochemistry (IHC)
HCC tissues and mouse tumor tissues were fixed by 4% paraformaldehyde at room temperature for 24 h. Embedded in paraffin, dewaxed and rehydrated, antigen extraction was performed using corresponding antigen retrieval solution at 100 °C for 20 min. Later, 3% hydrogen peroxide solution was applied to inhibit endogenous peroxidase activity for 15 min at room temperature. The slides were incubated with 5% normal goat serum in PBST for 1 h at room temperature to prevent nonspecific binding of antibodies. After that, samples are incubated with the designated primary antibody overnight at 4 °C. Each slide was incubated with the corresponding HRP-labeled secondary antibody, signal generated with DAB solution, and counterstained with hematoxylin after three TBST washings. Image J was used to analyze the intensity of IHC staining. The antibodies involved in this study are specified in Additional file 2: Table S1.
Animal experiment
Six-week-old male BALB/c nude mice were bought from Charles River Labs (Beijing, PR China). All animal experimental procedures were approved by the Experimental Animal Welfare Ethics Committee, Zhongnan Hospital of Wuhan University (ZN2022088). The six-week-old male mice were randomized into different groups. SK-hep1 cells (5 × 106 cells/mice) with stable knockdown ZMIZ2 and control SK-hep1 cells were subcutaneously injected into the nude mice. After tumor formation, tumor volume (mm3) was measured by Vernier calliper every other day and calculated by the formula (length × width × width)/2. All the mice were sacrificed after 37 days. The tumor samples were then fixed in 4% paraformaldehyde and embedded in paraffin for IHC analysis.
Coimmunoprecipitation (Co-IP)
Protein samples were first prepared, then LEF1 antibody (1:50, ABclonal, A0909) and IgG were added separately, and the tubes were gently vortexed and incubated overnight at 4 °C. The next day, 5 µl ProteinA/ProteinG (CST, 70024) Magnetic Beads were added to each tube and mixed gently at 4 °C for 3 h. Then the precipitate was centrifuged at 4 °C and washed with 1 × wash buffer. Next, the proteins were eluted from the beads using 1× SDS loading buffer and boiled at 95 °C for 10 min. Finally, western blot assays were performed for further verification.
Western blot
Cells and HCC samples were collected using RIPA buffer (P0013B, Beyotime, China) containing protease inhibitor (#HY-K0010, MCE, China) and phosphatase inhibitor (#HY-K0021&#HY-K0022, MCE, China). Proteins were centrifuged at 15,000 rpm for 20 min at 4 °C; Then supernatants were collected, and protein concentrations were calculated by the BCA Kit (P0012, Beyotime, China). RIPA buffer was mixed with protein loading buffer and protein was denatured at 100 °C for 5 min. An equal amount of proteins was separated using SDS-polyacrylamide gels and then transferred to polyvinylidene fluoride (PVDF) membranes (Millipore, USA). Membranes were blocked for 2 h with 5% skimmed milk and incubated with primary antibodies overnight at 4 °C. The membrane was incubated with a secondary antibody at room temperature for 1.5 h After three washes with TBST. Protein bands were obtained using an enhanced chemiluminescence (ECL) kit (SQ202, Epizyme, China) using GENESys (Synoptics Ltd, China). The antibodies involved in this study are prescribed in Additional file 2: Table S1.
Dataset analysis
Gene expression of ZMIZ2 of HCC and normal tissues were obtained from publicly available datasets, including The Cancer Genome Atlas (TCGA)-LIHC cohort and Gene Expression Omnibus (GEO) datasets (accession number: GSE54236, GSE55092 and GSE144269). The gene expression profiles and corresponding clinical information were downloaded from the UCSC Xena browser (https://xenabrowser.net) and GEO database (http://wwwncbinlmnih.gov/geo/), respectively.
Tumor somatic mutation and immune infiltration analysis
Somatic mutation data of the TCGA-LIHC cohort were downloaded from TCGA database (http://cancergenome.nih.gov/). The top 15 mutation genes were obtained and then compared between ZNIZ2 subgroups. According to the definition of the tumor mutation burden (TMB) in published literature [16], we evaluated the TMB score of patients in the TCGA-LIHC cohort. The relative abundance of 28 immune cell infiltration in each HCC sample was quantified via ssGSEA algorithm by previously described methods [38].
Functional annotation, gene set variation analysis (GSVA), and gene set enrichment analysis (GSEA)
To elucidate the role of ZMIZ2 in HCC, gene ontology (GO) functional annotation analysis was performed based on the top 200 differentially expressed genes (DEGs) (high-ZMIZ2 versus low-ZMIZ2). The DEGs were screened using “limma” R package. GSVA and GSEA was also conducted between ZMIZ2 subgroups via the “clusterProfler” and “GSVA” R packages, respectively. The gene set of “H.all.v7.4” was downloaded from the Molecular Signatures Database [39].
Statistical analysis
Data were represented as means ± standard deviation. Data analyses, including normality test, paired t-test, Spearman correlation test, Kaplan-Meier method, and univariate and multivariate analyses by Cox proportional hazard regression models, were performed using R (version 4.1.1), SPSS version 20.0 (IBM SPSS, Armonk, NY, USA) and GraphPad PRISM version 9.0 (GraphPad Software Inc., San Diego, CA, USA) software. Differences between data were considered statistically significant at P < 0.05.
Results
ZMIZ2 was overexpressed in HCC and high ZMIZ2 was correlated with poor clinical prognosis
A pan-cancer analysis of ZMIZ2 mRNA expression levels in 26 common human cancers showed that ZMIZ2 was overexpressed in various tumors compared with adjacent normal tissue (logFC > 0). Of note, ZMIZ2 mRNA expression levels were remarkably higher in HCC (n = 369) than that in normal tissues (n = 50) (logFC = 1.30, P < 0.001) (Fig. 1A). Univariate Cox analysis revealed that high ZMIZ2 expression was associated with poor prognosis in Liver hepatocellular carcinoma (LIHC) (HR = 1.26, p = 0.04) (Fig. 1B). Paired Student’s t-test analysis confirmed the ZMIZ2 mRNA expression pattern in TCGA-LIHC (Fig. 1C). In addition, analysis in three independent GEO HCC cohorts (GSE54236, GSE55092, and GSE144269) further confirmed that ZMIZ2 was significantly upregulated in HCC tissues (Fig. 1D). Furthermore, the mRNA and protein expression of ZMIZ2 was highly increased in several hepatoma cell lines compared with immortalized LO2 human hepatocytes (Figs. 1E and F and 2A). Higher protein levels were also found in eight HCC specimens than in adjacent normal tissues via western blot (Fig. 1G) and IHC (Fig. 1H). Regarding the potential indicating role of ZMIZ2 in clinical prognosis, Kaplan-Meier survival analysis showed that higher ZMIZ2 expression was associated with shorter overall survival times (Fig. 1I), and subsequent univariate Cox regression analysis indicated that ZMIZ2 expression was a risk factor for HCC prognosis (Table 1). The correlation between ZMIZ2 and clinical features of HCC patients in TCGA-LIHC was presented in Table 2. Overall, these results indicated that ZMIZ2 was upregulated in HCC and elevated expression of ZMIZ2 was associated with worse clinical outcomes.
ZMIZ2 promotes HCC cell proliferation
To investigate the role of ZMIZ2 in the development of HCC, we chose several HCC cells lines to perform related experiments. Among the six HCC cell lines, Huh7, Hep3B, MHCC97H and SK-hep1 cells showed the highest ZMIZ2 expression at mRNA levels, while HepG2 and LM3 cells showed relatively low expression levels of ZMIZ2 (Figs. 1E and 2A). Therefore, Huh7, Hep3B and SK-hep1 were selected for ZMIZ2 silencing, and HepG2 was chosen for ZMIZ2 overexpression. Three siRNAs (si-ZMIZ2#1, si-ZMIZ2#2, si-ZMIZ2#3) and plasmids (OE-ZMIZ2) were transfected into the cells. The efficiencies of ZMIZ2 silencing and overexpression were validated by qRT-PCR and Western blotting (Figs. 3A and B and 2B and C). The results showed that si-ZMIZ2#2 and si-ZMIZ2#3 had a better inhibition efficiency (Fig. 3A, B), which were chosen for subsequent experiments. CCK8 assay, Luminescent Cell Viability Assay and EDU assay revealed that ZMIZ2 knockdown significantly impaired HCC cell proliferation (Figs. 3C-E and 2D-F). In contrast, ZMIZ2 overexpression promoted the HCC cell proliferation (Fig. 3C-E).
ZMIZ2 promotes HCC cell migration and invasion in vitro
Previous studies have shown that ZMIZ2 plays a role in the invasion and migration of breast cancer cells [18]. Herein, the wound healing assays also revealed a significant reduction of migration in the siZMIZ2 group compared with the control group (Figs. 3F and 2G). The transwell assay confirmed that downregulation of ZMIZ2 suppressed the migration and invasion ability of Huh7, Hep3B and SK-hep1 cells (Figs. 3G and 2H). Conversely, overexpression of ZMIZ2 promoted the migration and invasion ability of HepG2 cells (Fig. 3F, G).
ZMIZ2 downregulation arrested cell cycle and promoted apoptosis of HCC cells
To determine the potential biological functions of ZMIZ2 in HCC, we performed GSEA and found that patients with high ZMIZ2 expression were significantly enriched in the cell cycle signaling and apoptosis pathways in the TCGA cohort (Additional file 1: Fig. S1). Subsequently, the flow cytometry analysis revealed the suppression of ZMIZ2 could significantly block the process of cell cycle and particularly induced cell cycle G0/G1 phase arrest in Huh7, Hep3B and SK-hep1 cells (Figs. 4A and 2I). Besides, the percentage of apoptosis cells was significantly increased by ZMIZ2 silencing (Figs. 4B and 2J). In contrast, ZMIZ2 overexpression showed opposite results (Fig. 4A, B). Western blot showed that ZMIZ2 silencing reduced the expression of Bcl-2, Cyclin B1, and Cyclin D1 but increased Bax, p21, p27, and p53 protein levels (Fig. 4C), indicating that ZMIZ2 was associated with cell cycle and apoptosis.
ZMIZ2 enhanced the oncogenesis of HCC cells in vivo
In order to explore the function of ZMIZ2 in vivo, we utilized lentivirus-mediated ZMIZ2 knockdown to establish sh-ZMIZ2 SK-hep1 cells. These cells were injected subcutaneously into the axilla of nude mice. We found the tumor growth rate was significantly slowed in the sh-ZMIZ2 group compared to the control group (Fig. 5A, B). The tumor size and weight were significantly reduced after ZMIZ2-knockdown (Fig. 5C, D), indicating that ZMIZ2 deficiency inhibited HCC tumor development in vivo. IHC staining also revealed that the ZMIZ2 knockdown group had lower levels of ZMIZ2, LEF1, and Ki-67 staining intensity compared to the normal tissues (Fig. 5E). Taken together, these findings suggest that ZMIZ2 plays a significant role in promoting oncogenesis and progression of HCC in vivo.
Different patterns of tumor somatic mutations in patients between ZMIZ2 expression subgroups
We next investigated the tumor somatic alterations between high and low ZMIZ2 expression groups. The top 15 variant mutations were identified in 358 HCC patients of the TCGA-LIHC cohort (Fig. 6A). We found no significant difference in tumor mutation burden score between the two subgroups (Fig. 6B). However, TP53 mutations (p < 0.0001) were more frequent in samples with high ZMIZ2 expression (Fig. 6A), and ZMIZ2 expression was higher in TP53-mutated samples (Fig. 6C). Patients without TP53 mutation had a better prognosis than those with TP53 mutation (Fig. 6D), which is consistent with previous studies [40, 41]. Notably, patients with TP53 mutation as well as high ZMIZ2 expression had the poorest prognosis among all subgroups (Fig. 6E). These results indicated that ZMIZ2 and TP53 mutations may have a potential cooperative effect on the development and prognosis of HCC.
Correlation between ZMIZ2 expression and immune cells infiltration in HCC
Next, we analyzed the relationship between ZMIZ2 expression and immune cells infiltration by the ssGSEA algorithm. As shown in Fig. 7A, there was a positive correlation between ZMIZ2 expression and the infiltration of CD4 T cells and dendritic cells (DCs). In contrast, ZMIZ2 expression was significantly negatively correlated with the infiltration of multiple anti-tumor immune cells, such as NK cells, CD8 T cells, and Th1 cells. Consistently, the infiltration levels of those anti-tumor immune cells were significantly decreased in the ZMIZ2 high expression groups, while the infiltration scores of CD4 T cells and DCs were higher (Additional file 1: Fig. S3A). Since these results were based on bioinformatic analysis, we further validated the association between immune cells infiltration and ZMIZ2 expression in HCC through IHC. As the IHC staining shown (Fig. 7C), CD8A, the marker of CD8 T cells was highly expressed in HCC tissues with lower ZMIZ2 expression. In contrast, the markers of CD4 T cells (CD4) and DC cells (CD11c) had the positive correlations with the ZMIZ2 expression. These data were consistent with the bioinformatics results. Furthermore, we found that decreased activated CD8 T cells, enriched CD4 T cells, and enriched DCs in HCC were associated with poor survival outcomes, severing as risk factors for HCC prognosis (Fig. 7B). Notably, patients with high ZMIZ2 expression together with one of these factors had relatively shorter OS times than other subgroups (Additional file 1: Fig. S3B). These results suggest that ZMIZ2 may influence the prognosis of HCC by affecting the infiltration of immune cells.
ZMIZ2 activated the Wnt/β-catenin signaling pathways through LEF1
Next, we investigated potential biological distinction between high- and low-ZMIZ2 expression samples. First, differential expression analysis was conducted between high and low ZMIZ2 expression samples separated by the median expression in the TCGA-LIHC cohort (Additional file 1: Fig. S2A). According to GO analysis, the top 200 differential expression genes (DEGs) were significantly enriched in genomic modification-related processes such as chromatin modification, histone modification, and RNA splicing. (Additional file 1: Fig. S2B). Additionally, GSVA analysis revealed an activation of tumor progression-related pathways, including DNA repair, PI3K/AKT/MTOR signaling, and Wnt/β-catenin signaling, in the high-ZMIZ2 group (Fig. 8A). To further investigate the role of ZMIZ2 in HCC, we performed transcriptome sequencing in ZMIZ2-silenced Huh7 cells and control cells, and visualized the DEGs in a heat map (Additional file 1: Fig. S2C). These DEGs (siZMIZ2 versus siNC; adjust. p < 0.05) were substantially enriched in the Wnt signaling pathway according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis (Fig. 8B), and the siNC samples had higher activation levels of Wnt/β-catenin signaling pathway as confirmed by the GSEA analysis (Additional file 1: Fig. S2D). Taken together, these results suggest that ZMIZ2 plays a crucial role in regulating HCC malignant progression via the Wnt/β-catenin signaling pathway.
Ten DEGs (siZMIZ2 versus siNC) were enriched in the Wnt/β-catenin pathway according to the enrichment analysis, as shown in the heat map (Fig. 8C). We conducted qRT-PCR analysis and observed that three out of the ten genes were the potential targets of ZMIZ2, including LEF1, NFATC2 and NFATC3 (Fig. 8D). Verified by public databases, LEF1 was highly expressed in HCC tissues, which showed a significant positive correlation with the ZMIZ2 expression (Fig. 8E). In contrast, NFATC2 and NFATC3 showed no significant difference between HCC and normal liver tissues (Additional file 1: Fig. S2E). After inhibiting ZMIZ2 expression, the mRNA and protein expression levels of LEF1 were significantly reduced (Fig. 8D, F), and the protein expression levels of positive regulators in the Wnt/β-catenin pathways (Survivin, c-Myc, and β-Catenin) were decreased, whereas the negative regulators (GSK3β and P-GSK3β) were increased (Fig. 8F). An opposite trend was observed in the ZMIZ2-overexpressed HCC cell lines (Fig. 8F). We performed the immunoprecipitation (IP) assay and validated the direct interaction between ZMIZ2 and LEF1 (Fig. 8G). Collectively, these results suggest that ZMIZ2 could activate the Wnt/β-catenin signaling pathways by directly interact with LEF1.
Discussion
HCC remains a major cause of cancer-related mortality worldwide, despite significant research efforts over decades of years [42, 43]. Therefore, identification of novel prognostic biomarkers and treatment targets of HCC is an urgent need. This study aimed to investigate the biological functions of ZMIZ2 in HCC, which has not been previously described. We found that upregulation of ZMIZ2 was prevalent in HCC and elevated ZMIZ2 expression was associated with worse clinical outcome. Furthermore, our results demonstrate that ZMIZ2 promotes the malignant progression of HCC. These results highlight the potential of ZMIZ2 as a novel prognostic biomarker and therapeutic target for HCC.
ZMIZ2 functions as a transcriptional co-activator that regulates the activity of multiple transcription factors, including p53, Smads and nuclear hormone receptors [14,15,16,17]. As a tumor suppressor gene, TP53 is frequently mutant and inactive in liver cancer [44], indicating its inactivation plays a crucial role in the initiation and development of HCC. Smads are downstream targets of the TGF-β signaling pathway, TGF-β1/Smad3 signaling has been demonstrated to promote hepatic gluconeogenesis [45]. ZMIZ2 enhances the transcriptional activity of androgen receptor(AR) [46], which is highly expressed in HCC and is associated with increased tumor recurrence and decreased overall survival [47,48,49,50]. These findings suggest that ZMIZ2 plays an important role in HCC development. In this study, we found that ZMIZ2 was highly expressed in HCC tissues and cells, and the expression of ZMIZ2 was a significant independent risk factors for HCC patient survival. In addition, we found that ZMIZ2 promoted the proliferation of primary HCC in vitro and vivo. ZMIZ2 knockdown significantly inhibited the invasion and migration abilities of HCC cells. These results were consistent with previous findings in the breast cancer cells [18] and colorectal cancers [19]. Zou et al. speculated that ZMIZ2 may promote the progression of TNBC by regulating CCL5 to activate the NOD-like receptor signaling pathway and toll-like receptor signaling pathway [18]. Zhu et al. demonstrated that depletion of ZMIZ2 distinctly attenuated colorectal tumorigenesis in mice. ZMIZ2 recruits the enzyme USP7, which deubiquitylates and stabilizes β-catenin, thereby facilitating colorectal tumorigenesis [19]. These data support our conclusion that ZMIZ2 was an oncogene in HCC, and thus it could be a potential prognostic biomarker for this disease.
The tumor suppressor gene TP53 is the molecular hub of multiple crucial signaling pathways, and recent results of whole genome and RNA sequencing analysis for 254 HCC samples also suggested that TP53 was the most frequently mutated gene in HCC [6]. Our research found that HCC patients with high ZMIZ2 expression had higher mutation rate of TP53 and poorer overall survival outcomes, suggesting that ZMIZ2 overexpression and TP53 mutations may jointly affect the development and prognosis of HCC. In HCC, TP53 also shows a high mutation rate and is closely associated with poor prognosis [51, 52], consistent with our analysis. A growing number of studies have indicated a correlation between tumor immune infiltration and cancer progression [53, 54]. Lin et al. found that circRanGAP1 accelerated the HCC progression by regulating the miR-27b-3p/NRAS/ERK axis and affecting the infiltration levels of tumor-associated macrophages [8]. Similarly, we observed a correlation between ZMIZ2 expression and immune infiltration in HCC, which suggested that ZMIZ2 may affect the infiltration of immune cells.
The Wnt/β-catenin signaling pathway regulates the proliferation and differentiation of stem cells and also participates in the expansion of cancer cells. Accumulating evidence has implicated critical roles for both ZMIZ proteins and Wnt/β-catenin signaling pathways in development and tumorigenesis. ZMIZ1 is a coactivator of Notch1 that induces TCF-1 to form a complex with β-catenin, leading to transcriptional activation of target genes via the classical Wnt pathway [55]. In the prostate cancer, ZMIZ2 interacts with β-catenin and improves β-catenin-mediated transcriptions [13]. A similar effect was found in our study. GSEA and transcriptome sequencing analysis showed that ZMIZ2 knockdown had important effects on genes that are mainly related to Wnt/β-catenin signaling pathway. At present, the potential interaction between ZMIZ2 and other members of the Wnt pathway remains unclear. LEF1, as a key component of the Wnt/β-catenin signaling, plays an important role in HCC tumorigenesis, differentiation, and cell survival. Overexpression of LEF1 has been frequently observed in primary liver cancer tissues, and it promotes stemness and poor differentiation of hepatocellular carcinoma by directly activating the NOTCH pathway [26]. LEF1 can increase self-renewal and invasiveness in HCC cells through transcriptional regulation of Oct4 and EMT regulators [22]. Our study suggests that ZMIZ2 interacts with LEF1, therefore activates the Wnt/β-catenin signaling pathway. However, the binding site between ZMIZ2 and LEF1 warrants further studies.
Conclusion
Our study suggests that ZMIZ2 is significantly upregulated in HCC tissues and is associated with poor prognosis in HCC patients. Mechanistically, ZMIZ2 promotes HCC cell proliferation, migration, and invasion by activating the Wnt/β-catenin signaling pathway through the upregulation of LEF1 (Fig. 9). Therefore, ZMIZ2 may serve as a potential prognostic biomarker and therapeutic target for HCC.
Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- AR:
-
Androgen receptor
- CCK8:
-
Cell counting kit-8
- DEGs:
-
Differential expression genes
- DFI:
-
Disease-free interval
- GEO:
-
Gene expression omnibus
- GO:
-
Gene ontology
- GSEA:
-
Gene set enrichment analysis
- GSVA:
-
Gene set variation analysis
- HCC:
-
Hepatocellular carcinoma
- IHC:
-
Immunohistochemistry
- LEF1:
-
Lymphoid enhancer-binding factor-1
- OS:
-
Overall survival
- PFI:
-
Progression-free interval
- TAD:
-
Transcription activation domain
- TCGA:
-
The cancer genome atlas
- TERT:
-
Telomere reverse transcriptase
- TMB:
-
Tumor mutation burden
- TME:
-
Tumor microenvironment
- ZMIZ2:
-
Zinc finger MIZ type containing 2
References
Vogel A, Meyer T, Sapisochin G, Salem R, Saborowski A. Hepatocellular carcinoma. Lancet. 2022;400(10360):1345–62.
Ma XL, Nie YY, Xie SH, Zheng H, Tong Y, Wang YC, et al. ASAP2 interrupts c-MET-CIN85 interaction to sustain HGF/c-MET-induced malignant potentials in hepatocellular carcinoma. Exp Hematol Oncol. 2023;12(1):38.
Wang W, Huang Q, Liao Z, Zhang H, Liu Y, Liu F, et al. ALKBH5 prevents hepatocellular carcinoma progression by post-transcriptional inhibition of PAQR4 in an m6A dependent manner. Exp Hematol Oncol. 2023;12(1):1.
Zhang H, Zhang W, Jiang L, Chen Y. Recent advances in systemic therapy for hepatocellular carcinoma. Biomark Res. 2022;10(1):3.
Chen E, He Y, Jiang J, Yi J, Zou Z, Song Q, et al. CDCA8 induced by NF-YA promotes hepatocellular carcinoma progression by regulating the MEK/ERK pathway. Exp Hematol Oncol. 2023;12(1):9.
Chang YS, Tu SJ, Chen HD, Chung CC, Hsu MH, Chou YP, et al. Whole genome and RNA sequencing analyses for 254 Taiwanese hepatocellular carcinomas. Biomark Res. 2023;11(1):68.
Wang T, Rao D, Yu C, Sheng J, Luo Y, Xia L, et al. RHO GTPase family in hepatocellular carcinoma. Exp Hematol Oncol. 2022;11(1):91.
Lin XH, Liu ZY, Zhang DY, Zhang S, Tang WQ, Li DP, et al. circRanGAP1/miR-27b-3p/NRAS Axis may promote the progression of hepatocellular Carcinoma. Exp Hematol Oncol. 2022;11(1):92.
Suzuki T, Yano H, Nakashima Y, Nakashima O, Kojiro M. Beta-catenin expression in hepatocellular carcinoma: a possible participation of beta-catenin in the dedifferentiation process. J Gastroenterol Hepatol. 2002;17(9):994–1000.
Fujie H, Moriya K, Shintani Y, Tsutsumi T, Takayama T, Makuuchi M, et al. Frequent beta-catenin aberration in human hepatocellular carcinoma. Hepatol Res. 2001;20(1):39–51.
Vilchez V, Turcios L, Marti F, Gedaly R. Targeting Wnt/beta-catenin pathway in hepatocellular carcinoma treatment. World J Gastroenterol. 2016;22(2):823–32.
Zhang Y, Wang X. Targeting the Wnt/beta-catenin signaling pathway in cancer. J Hematol Oncol. 2020;13(1):165.
Lee SH, Zhu C, Peng Y, Johnson DT, Lehmann L, Sun Z. Identification of a novel role of ZMIZ2 protein in regulating the activity of the Wnt/beta-catenin signaling pathway. J Biol Chem. 2013;288(50):35913–24.
Jackson PK. A new RING for SUMO: wrestling transcriptional responses into nuclear bodies with PIAS family E3 SUMO ligases. Genes Dev. 2001;15(23):3053–8.
Kotaja N, Karvonen U, Janne OA, Palvimo JJ. PIAS proteins modulate transcription factors by functioning as SUMO-1 ligases. Mol Cell Biol. 2002;22(14):5222–34.
Megidish T, Xu JH, Xu CW. Activation of p53 by protein inhibitor of activated Stat1 (PIAS1). J Biol Chem. 2002;277(10):8255–9.
Schmidt D, Muller S. PIAS/SUMO: new partners in transcriptional regulation. Cell Mol Life Sci. 2003;60(12):2561–74.
Zou X, Liu Y, Di J, Wei W, Watanabe N, Li J, et al. ZMIZ2 promotes the development of triple-receptor negative breast cancer. Cancer Cell Int. 2022;22(1):52.
Zhu Y, Gu L, Lin X, Cui K, Liu C, Lu B, et al. LINC00265 promotes colorectal tumorigenesis via ZMIZ2 and USP7-mediated stabilization of beta-catenin. Cell Death Differ. 2020;27(4):1316–27.
Behrens J, von Kries JP, Kuhl M, Bruhn L, Wedlich D, Grosschedl R, et al. Functional interaction of beta-catenin with the transcription factor LEF-1. Nature. 1996;382(6592):638–42.
Petropoulos K, Arseni N, Schessl C, Stadler CR, Rawat VP, Deshpande AJ, et al. A novel role for Lef-1, a central transcription mediator of wnt signaling, in leukemogenesis. J Exp Med. 2008;205(3):515–22.
Chen CL, Tsai YS, Huang YH, Liang YJ, Sun YY, Su CW, et al. Lymphoid enhancer factor 1 contributes to Hepatocellular Carcinoma Progression through Transcriptional Regulation of epithelial-mesenchymal transition regulators and stemness genes. Hepatol Commun. 2018;2(11):1392–407.
Zhang Z, Min L, Li H, Chen L, Zhao Y, Liu S, et al. Asporin represses gastric cancer apoptosis via activating LEF1-mediated gene transcription independent of beta-catenin. Oncogene. 2021;40(27):4552–66.
Chen X, Tu J, Liu C, Wang L, Yuan X. MicroRNA-621 functions as a metastasis suppressor in colorectal cancer by directly targeting LEF1 and suppressing Wnt/beta-catenin signaling. Life Sci. 2022;308:120941.
Zhou X, Li X, Wang R, Hua D, Sun C, Yu L, et al. Recruitment of LEF1 by Pontin chromatin modifier amplifies TGFBR2 transcription and activates TGFbeta/SMAD signalling during gliomagenesis. Cell Death Dis. 2022;13(9):818.
Fang S, Liu M, Li L, Zhang FF, Li Y, Yan Q, et al. Lymphoid enhancer-binding factor-1 promotes stemness and poor differentiation of hepatocellular carcinoma by directly activating the NOTCH pathway. Oncogene. 2019;38(21):4061–74.
Chen Y, Guo Y, Li Y, Yang J, Liu J, Wu Q, et al. miR–300 regulates tumor proliferation and metastasis by targeting lymphoid enhancer–binding factor 1 in hepatocellular carcinoma. Int J Oncol. 2019;54(4):1282–94.
Tian X, Li J, Ma ZM, Zhao C, Wan DF, Wen YM. Role of hepatitis B surface antigen in the development of hepatocellular carcinoma: regulation of lymphoid enhancer-binding factor 1. J Exp Clin Cancer Res. 2009;28(1):58.
Wada H, Sato Y, Fujimoto S, Okamoto K, Bando M, Kawaguchi T, et al. Resveratrol inhibits development of colorectal adenoma via suppression of LEF1; comprehensive analysis with connectivity map. Cancer Sci. 2022;113(12):4374–84.
Zehir A, Benayed R, Shah RH, Syed A, Middha S, Kim HR, et al. Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat Med. 2017;23(6):703–13.
Kandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, et al. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502(7471):333–9.
Sadot E, Geiger B, Oren M. Ben-Ze’ev A. Down-regulation of beta-catenin by activated p53. Mol Cell Biol. 2001;21(20):6768–81.
Cagatay T, Ozturk M. P53 mutation as a source of aberrant beta-catenin accumulation in cancer cells. Oncogene. 2002;21(52):7971–80.
Prowald A, Cronauer MV, von Klot C, Eilers T, Rinnab L, Herrmann T, et al. Modulation of beta-catenin-mediated TCF-signalling in prostate cancer cell lines by wild-type and mutant p53. Prostate. 2007;67(16):1751–60.
Pan Z, Tian Y, Cao C, Niu G. The emerging role of YAP/TAZ in Tumor Immunity. Mol Cancer Res. 2019;17(9):1777–86.
van Janse HJ, Azad T, Ling M, Hao Y, Snetsinger B, Khanal P, et al. The Hippo Pathway Component TAZ promotes Immune Evasion in Human Cancer through PD-L1. Cancer Res. 2018;78(6):1457–70.
Lu J, He X, Zhang L, Zhang R, Li W. Acetylation in Tumor Immune Evasion Regulation. Front Pharmacol. 2021;12:771588.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25.
Villanueva A, Hoshida Y. Depicting the role of TP53 in hepatocellular carcinoma progression. J Hepatol. 2011;55(3):724–5.
Zhan P, Ji YN. Prognostic significance of TP53 expression for patients with hepatocellular carcinoma: a meta-analysis. Hepatobiliary Surg Nutr. 2014;3(1):11–7.
Kotewall CN, Cheung TT. Optimizing hepatectomy for hepatocellular carcinoma in Asia-patient selection and special considerations. Transl Gastroenterol Hepatol. 2018;3:75.
Cao W, Chen HD, Yu YW, Li N, Chen WQ. Changing profiles of cancer burden worldwide and in China: a secondary analysis of the global cancer statistics 2020. Chin Med J (Engl). 2021;134(7):783–91.
Ueda H, Ullrich SJ, Gangemi JD, Kappel CA, Ngo L, Feitelson MA, et al. Functional inactivation but not structural mutation of p53 causes liver cancer. Nat Genet. 1995;9(1):41–7.
Shi X, Yang J, Deng S, Xu H, Wu D, Zeng Q, et al. TGF-beta signaling in the tumor metabolic microenvironment and targeted therapies. J Hematol Oncol. 2022;15(1):135.
Huang CY, Beliakoff J, Li X, Lee J, Li X, Sharma M, et al. hZimp7, a novel PIAS-like protein, enhances androgen receptor-mediated transcription and interacts with SWI/SNF-like BAF complexes. Mol Endocrinol. 2005;19(12):2915–29.
Qin Z, Liu X, Li Z, Wang G, Feng Z, Liu Y, et al. LncRNA LINC00667 aggravates the progression of hepatocellular carcinoma by regulating androgen receptor expression as a miRNA-130a-3p sponge. Cell Death Discov. 2021;7(1):387.
Ren QN, Zhang H, Sun CY, Zhou YF, Yang XF, Long JW, et al. Phosphorylation of androgen receptor by mTORC1 promotes liver steatosis and tumorigenesis. Hepatology. 2022;75(5):1123–38.
Montgomery EJ, Xing E, Campbell MJ, Li PK, Blachly JS, Tsung A, et al. Constitutively active androgen receptor in Hepatocellular Carcinoma. Int J Mol Sci. 2022;23(22):13768.
Oaks Z, Patel A, Huang N, Choudhary G, Winans T, Faludi T, et al. Cytosolic aldose metabolism contributes to progression from cirrhosis to hepatocarcinogenesis. Nat Metab. 2023;5(1):41–60.
Shen T, Li SF, Wang JL, Zhang T, Zhang S, Chen HT, et al. TP53 R249S mutation detected in circulating tumour DNA is associated with prognosis of hepatocellular carcinoma patients with or without hepatectomy. Liver Int. 2020;40(11):2834–47.
Woo HG, Wang XW, Budhu A, Kim YH, Kwon SM, Tang ZY, et al. Association of TP53 mutations with stem cell-like gene expression and survival of patients with hepatocellular carcinoma. Gastroenterology. 2011;140(3):1063–70.
Bayik D, Lathia JD. Cancer stem cell-immune cell crosstalk in tumour progression. Nat Rev Cancer. 2021;21(8):526–36.
Pena-Romero AC, Orenes-Pinero E. Dual effect of Immune cells within Tumour Microenvironment: Pro- and Anti-tumour effects and their triggers. Cancers (Basel). 2022;14(7):1681.
Bao X, Qin Y, Lu L, Zheng M. Transcriptional regulation of early T-Lymphocyte development in Thymus. Front Immunol. 2022;13:884569.
Acknowledgements
We sincerely thank the team members for their dedication to this project.
Funding
This study was sponsored by grants from the National Nature Science Foundation of China (No.82172983, 81670554, 82073095), the Wuhan Science and Technology Project (No.2020020601012208), science and technology innovation and cultivation fund of Zhongnan Hospital of Wuhan University (No. CXPY2020042) and subjects and platform construction in Zhongnan Hospital of Wuhan University (No. PTXM2020003).
Author information
Authors and Affiliations
Contributions
Conceived and designed the experiments: YC, XXH. Performed the experiments: YD, HK. Analysed data: YD, YMN. Contributed reagents/materials/analysis tools: Yuan Y, LK, CW, Yun Y, JHH, LRL. Wrote the paper: YD, YMN. YC and XXH reviewed the paper. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the ethics committee of Zhongnan Hospital of Wuhan University (Approval No.2022011 K), and all individuals provided informed consent for inclusion of their tissues in this study. Moreover, all animal experimental procedures were approved by the Experimental Animal Welfare Ethics Committee, Zhongnan Hospital of Wuhan University (Approval No. ZN2022088).
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.
Electronic supplementary material
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.
About this article
Cite this article
Ding, Y., Ning, Y., Kang, H. et al. ZMIZ2 facilitates hepatocellular carcinoma progression via LEF1 mediated activation of Wnt/β-catenin pathway. Exp Hematol Oncol 13, 5 (2024). https://doi.org/10.1186/s40164-024-00475-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40164-024-00475-w