Skip to main content

LincRNA ZNF529-AS1 inhibits hepatocellular carcinoma via FBXO31 and predicts the prognosis of hepatocellular carcinoma patients

Abstract

Background

Invasion and metastasis of hepatocellular carcinoma (HCC) is still an important reason for poor prognosis. LincRNA ZNF529-AS1 is a recently identified tumour-associated molecule that is differentially expressed in a variety of tumours, but its role in HCC is still unclear. This study investigated the expression and function of ZNF529-AS1 in HCC and explored the prognostic significance of ZNF529-AS1 in HCC.

Methods

Based on HCC information in TCGA and other databases, the relationship between the expression of ZNF529-AS1 and clinicopathological characteristics of HCC was analysed by the Wilcoxon signed-rank test and logistic regression. The relationship between ZNF529-AS1 and HCC prognosis was evaluated by Kaplan‒Meier and Cox regression analyses. The cellular function and signalling pathways involved in ZNF529-AS1 were analysed by GO and KEGG enrichment analysis. The relationship between ZNF529-AS1 and immunological signatures in the HCC tumour microenvironment was analysed by the ssGSEA algorithm and CIBERSORT algorithm. HCC cell invasion and migration were investigated by the Transwell assay. Gene and protein expression were detected by PCR and western blot analysis, respectively.

Results

ZNF529-AS1 was differentially expressed in various types of tumours and was highly expressed in HCC. The expression of ZNF529-AS1 was closely correlated with the age, sex, T stage, M stage and pathological grade of HCC patients. Univariate and multivariate analyses showed that ZNF529-AS1 was significantly associated with poor prognosis of HCC patients and could be an independent prognostic indicator of HCC. Immunological analysis showed that the expression of ZNF529-AS1 was correlated with the abundance and immune function of various immune cells. Knockdown of ZNF529-AS1 in HCC cells inhibited cell invasion and migration and inhibited the expression of FBXO31.

Conclusion

ZNF529-AS1 could be a new prognostic marker for HCC. FBXO31 may be the downstream target of ZNF529-AS1 in HCC.

Peer Review reports

Introduction

Hepatocellular carcinoma (HCC) is one of the most common malignant tumours in the world [1, 2]. The morbidity and mortality of HCC rank sixth and third, respectively [3]. The incidence of HCC in China is higher than that in Western countries [4, 5]. At present, the preferred method for the treatment of HCC is surgical resection for BCLC stage 0-B disease [6]. However, most patients have HCC metastasis at the time of treatment, and some patients have a high recurrence and metastasis rate after surgery, resulting in poor prognosis [7, 8]. In recent years, immunotherapy has been gradually applied by clinicians and has shown efficacy in the clinic [9], but some patients are not sensitive to immunotherapy. Therefore, it is necessary to identify key molecules related to the immune infiltrate of HCC, which provide better predictive targets for the prognosis and immunotherapy sensitivity of HCC.

ZNF529 is a member of the zinc finger protein (ZNF) family. The ZNF family plays a role through sequence-specific binding to downstream target molecules such as DNA, RNA, and DNA‒RNA, or binding to other zinc finger proteins (ZNFs) [10]. Many studies have shown that the abnormal expression of ZNFs plays an important role in the occurrence and development of HCC [11, 12]. The ZNF family also plays an important role in cancer prognosis, immune infiltration, the tumour microenvironment, epigenetics, and drug sensitivity [13]. ZNF529-AS1 (ZNF529-antisense 1) is the antisense strand of ZNF529 and plays the opposite role to ZNF529. We speculate that ZNF529-AS1 may play a regulatory role in the occurrence and development of HCC.

In this study, we first analysed the correlation of ZNF529-AS1 with clinicopathological characteristics, histological characteristics, survival prognosis, immune infiltration and sensitivity to chemotherapy in HCC patients in the database. In addition, we investigated the biological role of ZNF529-AS1 in regulating HCC cell invasion via FBXO31 in vitro.

Materials and methods

Patient datasets

We obtained RNA-seq data in TPM format from UCSC XENA (https://xenabrowser.net/datapages/), which were uniformly processed by the Toil program of TCGA and GTEx. The downloaded data were used to compare the expression level between 33 cancerous tissues and matched peri-cancerous tissues. We extracted the corresponding normal tissue data in LIHC and GTEx of TCGA, compared RNA-seq data in TPM format and performed logarithmic transformation on the expression between samples. We collected RNA-seq data and clinical information from a total of 424 HCC patients, including 50 paired peri-cancerous tissues for subsequent analysis. A total of 122 immunomodulatory factors, including MHC, receptors, chemokines, and immunostimulatory factors, were collected from the study by Charoentong et al. [14]. In addition, immune checkpoint molecules were retrieved from the study by Auslander et al. [15]. Spearman correlation analysis was used to explore the correlation between ZNF529-AS1 and the above factors in cancer.

Survival analysis

The data of 374 HCC patients were collected and ranked according to the expression level of ZNF529-AS1 from low to high, with the median between the 187th and 188th patients. A value higher than this median was regarded as the high expression group; otherwise, it was regarded as the low expression group. The Kaplan‒Meier method [16] was used to analyse the relationship between ZNF529-AS1 expression and overall survival (OS), disease-specific survival (DSS) and progression-free interval (PFI) in HCC patients. The correlation between ZNF529-AS1 expression and DFS of HCC was analysed based on the GEPIA database (http://gepia.cancer-pku.cn/) [17]. Hazard ratios (HRs) and 95% confidence intervals were calculated using univariate survival analysis.

Univariate and multivariate logistic regression analysis

We used univariate Cox regression analysis to calculate the relationship between the ZNF529-AS1 expression level and HCC prognosis. Multivariate analysis was used to evaluate whether ZNF529-AS1 was an independent prognostic factor for the survival of HCC patients. p < 0.05 was considered statistically significant.

ZNF529-AS1-related functional enrichment analysis and drug sensitivity correlation analysis

Differentially expressed genes between the ZNF529-AS1 high expression and low expression groups were determined using the ‘limma’ R package [18]. The coexpression heatmap of ZNF529-AS1 and the top 20 differentially expressed genes was plotted according to Pearson correlation. Gene set enrichment analysis (GSEA) and GO and KEGG enrichment analyses were performed using the ‘clusterProfiler’ R package [19]. The differentially expressed genes between the ZNF529-AS1 high expression and low expression groups were analysed using the ‘DEseq2’ R package. All differentially expressed genes were pooled for GSEA, and the genes with |log2FC|> 1 and adjusted p value < 0.05 were selected for GO and KEGG enrichment analysis. The relationship between ZNF529-AS1 expression and drug sensitivity was analysed using the CellMiner database (http://discover.nci.nih.gov/cellminer/) [20, 21]. Data processing and plotting were performed through the “impute”, “limma” and “ggpubr" packages in R.

Comprehensive analysis of the tumour microenvironment (TME) and immune cell infiltration

We obtained the gene set related to the cancer-immunity cycle from the website developed by Xu et al. (http://biocc.hrbmu.edu.cn/TIP/) [22]. We obtained the data of a set of gene signatures positively correlated with the clinical response to an anti-PD-L1 drug (atezolizumab) from the study by Mariathasan et al. [23]. GSVA was used to calculate the enrichment scores of gene signatures positively correlated with the cancer-immunity cycle and immunotherapy. We also collected gene sets related to immune function to analyse the activity of related pathways.

The stromal score, immune score and ESTIMATEScore of each HCC sample were estimated using the ESTIMATE algorithm (https://bioinformatics.mdanderson.org/estimate/). Single-sample gene set enrichment analysis (ssGSEA) was used to quantify the relative activity of various immune cells and immune functions in the TME, which was normalized to a unit distribution from 0 to 1. The Pearson correlation of ZNF529-AS1 with these immune cells and immune functions was analysed. A total of 22 inhibitory immune checkpoints with therapeutic potential were collected from the study by Auslander et al. [24], and their differences between the ZNF529-AS1 high expression and low expression groups were analysed.

Cell culture and transfection

The HCC cells Hep3B and HUH7 and normal hepatocyte MHAs were obtained from the Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences (Cat. No.: SCSP-5045, SCSP-526). The small interfering RNA (siRNA) of ZNF529-AS1 was purchased from Hippo Biotechnology Co., Ltd. Cells were transfected with siRNA as previously reported [25].

Quantitative real-time PCR (qRT‒PCR)

The expression levels of ZNF529-AS and FBXO31 in HCC were determined by q-RT‒PCR according to protocols in a previous report [26]. The primers were as follows: ZNF529-AS1 forward 5′–3′: GACTCGCCAGGCTTAACTCA, reverse 5′–3′: GAGCTTCAGGGCTCTTTCGT; FBXO31, forward 5′–3′: TGGAGGACATCTTCCACGAG, reverse 5′–3′: GCATTGTCGTACTGACTGG; and GAPDH, forward 5′–3′: GGAGCGAGATCCCTCCAAAAT reverse 5′–3′: GGCTGTTGTCATACTTCTCATGG.

Transwell assay

The invasion and migration assays of HCC cells were carried out 48 h after HCC cells were transfected with ZNF529-AS1. The details were the same as previously reported [27].

Western blotting

Western blot analysis was performed using previously reported protocols [28]. The materials used included primary antibodies against FBXO31 (1:1000, Proteintech, USA) and GAPDH (1:500, Proteintech, USA) and a horseradish peroxidase (HRP)-linked secondary antibody (Cell Signalling, USA).

Statistical analysis

All analyses were performed using the R programming language, and the data from different groups were compared by the Kruskal‒Wallis test, T test or Wilcoxon test. If not specified, p < 0.05 was considered statistically significant.

Results

Significant correlation of ZNF529-AS1 expression with immune factors in HCC

Differential expression of ZNF529-AS1 was found in various tumour tissues and their peri-cancerous tissues (Fig. 1A, B), including breast cancer, HCC and lung adenocarcinoma. In HCC, ZNF529-AS1 expression was higher in cancerous tissues than in peri-cancerous tissues (Fig. 1C). The area under the curve (AUC) of ZNF529-AS1 was 0.916, indicating that ZNF529-AS1 may be a potential biomarker in HCC tissues (Fig. 1D). ZNF529-AS1 was significantly associated with MHC molecules, chemokines, receptors and immunostimulatory factors of cancer immune cells (Fig. 1E–H). These results show that the function of the tumour-associated molecule ZNF529-AS1 is associated with immune infiltration.

Fig. 1
figure 1

Expression and immunological characteristics of ZNF529-AS1 in different cancers. A The Wilcoxon rank sum test was used to analyse the differential expression of ZNF529-AS1 in tumour tissues and adjacent tissues. B The Wilcoxon rank sum test was used to analyse the differential expression of ZNF529-AS1 in normal adjacent tissues of GTEx-bound TCGA and tumour tissues of TCGA. C Differential expression of ZNF529-AS1 in tumour samples and matched adjacent samples. D ROC curve shows the efficiency of the ZNF529-AS1 expression level to distinguish HCC tissues from nontumor tissues. The X-axis denotes the false-positive rate, and the Y-axis denotes the true positive rate. EH Heatmaps show the relationship between ZNF529-AS1 and H MHC molecules, F chemokines, G receptors and H immune checkpoints. For each relationship, the upper-left triangle is coloured to denote the p value, and the lower-right triangle is coloured to denote the correlation coefficient

Correlation analysis between ZNF529-AS1 expression and clinical characteristics of HCC

The clinicopathological characteristics of HCC patients are shown in Table 1. The clinical and gene expression characteristics of 374 HCC patients were collected from the TCGA-LIHC database. According to the mean value of ZNF529-AS1 expression, HCC patients were divided into two groups: a high expression group (n = 187) and a low expression group (n = 187). The relationship between the ZNF529-AS1 expression level and the clinicopathological characteristics of HCC was analysed by the Wilcoxon rank sum test. By analysing the relationship between ZNF529-AS1 expression and clinical parameters, we found that ZNF529-AS1 expression was correlated with T stage (p < 0.05) (Fig. 2A), pathological stage (p < 0.05) (Fig. 2B), histological grade (p < 0.001) (Fig. 2C), age (p = 0.01) (Fig. 2D), sex (p = 0.01) (Fig. 2E), weight (p < 0.001) (Fig. 2F), BMI (p = 0.03) (Fig. 2G), AFP (p < 0.001) (Fig. 2H), vascular invasion (p < 0.01) (Fig. 2I), and Child‒Pugh grade (p < 0.01) (Fig. 2J) but not with other clinical characteristics (Additional file 1: Fig. S1). The relationship between the ZNF529-AS1 expression level and the clinicopathological characteristics of HCC patients was analysed by logistic regression (Table 2).

Table 1 Correlation between ZNF529-AS1 expression and clinicopathological characteristics in HCC
Fig. 2
figure 2

The relationship between ZNF529-AS1 expression and clinical characteristics of HCC patients: A T stage, B pathological stage, C histologic grade, D age, E sex, F weight, G BMI, H AFP, I vascular invasion, and J Child‒Pugh grade

Table 2 ZNF529-AS1 expression associated with clinicopathologic characteristics (logistic regression)

Correlation of high ZNF529-AS1 expression with poor prognosis of HCC

Survival analysis of the TCGA-LIHC dataset showed that high ZNF529-AS1 expression in patients was significantly correlated with OS (p < 0.01) (Fig. 3A), DSS (p < 0.05) (Fig. 3B), and PFI (p < 0.05) (Fig. 3C) (poor prognosis). In addition, the analysis of the GIPIA2 database showed that high ZNF529-AS1 expression in patients was significantly correlated with low DFS (p < 0.05) (Fig. 3D).

Fig. 3
figure 3

Kaplan‒Meier survival analysis of different ZNF529-AS1 expression states in LIHC from the TCGA dataset. A Overall survival; B disease-specific survival; C progression-free interval; D disease-free survival

Univariate Cox analysis showed that T stage (p < 0.001), M stage (p < 0.05), pathological grade (p < 0.001) and ZNF529-AS1 expression (p < 0.01) were significantly correlated with OS in HCC patients (Fig. 4A). To further confirm the prognostic factors of HCC patients, multivariate risk analysis was performed by Cox regression, and the results showed that M stage (p < 0.05), histological grade (p < 0.05) and ZNF529-AS1 expression level (p = 0.013) were independent prognostic factors of OS (Fig. 4B). In addition, by constructing a Cox univariate regression model based on the DSS and PFI data of HCC patients and using Cox regression for multivariate risk analysis, we confirmed that ZNF529-AS1 expression was an independent prognostic factor of DSS (p < 0.05) and PFI (p < 0.01) in HCC patients (Additional file 1: Fig. S2). To further expand the clinical predictive value of ZNF529-AS1, we constructed a nomogram by combining multiple clinicopathological characteristics (Fig. 4C), and the nomograms of 1-, 3-, and 5-year survival probabilities were calibrated. The results showed that there was strong consistency between the prediction and the actual results (Fig. 4D).

Fig. 4
figure 4

ZNF529-AS1 is an independent prognostic indicator of HCC and quantitatively predicts the survival probability of 1-, 3-, and 5-year OS in HCC patients. Univariate (A) and multivariate (B) COX regression analyses of ZNF529-AS1 and different clinical characteristics. C Nomograms predicting the 1-, 3-, and 5-year overall survival probabilities of HCC patients. D Calibration plot of the nomograms to predict the 1-, 3-, and 5-year overall survival probabilities of HCC patients

Negative correlation of ZNF529-AS1 with HCC chemosensitivity

The correlation between ZNF529-AS1 and HCC chemosensitivity was explored through the CellMiner database. We identified 16 antitumour drugs that were significantly associated with ZNF529-AS1 expression. Among them, 14 drug sensitivities were significantly negatively correlated with ZNF529-AS1 expression, and 2 drug sensitivities were significantly positively correlated with ZNF529-AS1 (Fig. 5). The results showed that the expression level of ZNF529-AS1 was negatively correlated with sensitivity to most anticancer drugs (such as sunitinib, fluorouracil, and trametinib). Therefore, high expression of ZNF529-AS1 may affect chemotherapy efficacy against HCC.

Fig. 5
figure 5

Correlation analysis between the ZNF529-AS1 expression level and drug sensitivity. The horizontal axis denotes gene expression, and the vertical axis denotes drug sensitivity

Functional enrichment of ZNF529-AS1-related genes in HCC

To elucidate the underlying molecular mechanism of the role of ZNF529-AS1 in HCC, the “DESeq2” R package was used to evaluate differentially expressed genes (DEGs) between the ZNF529-AS1 high expression and low expression groups. A total of 187 samples with high ZNF529-AS1 expression and 187 samples with low ZNF529-AS1 expression were analysed. A total of 3034 DEGs were identified, including 2,455 upregulated and 579 downregulated DEGs, which were statistically significant between the ZNF529-AS1 high expression and low expression groups (|log2FC|> 1, p < 0.05) (Fig. 6A). The heatmap showed the top 10 upregulated DEGs and the top 10 downregulated DEGs between the ZNF529-AS1 high expression and low expression groups (Fig. 6B). To predict the functional enrichment information between high expression and low expression of 4987 DEGs, the clusterProfiler package was used to perform GO and KEGG functional enrichment analysis (Additional file 2: Tables S1 and S2).

Fig. 6
figure 6

Functional enrichment analysis of ZNF529-AS1-related genes in HCC. A Heatmap of the differentially expressed genes, including 2455 upregulated genes and 579 downregulated genes; B heatmap of 20 differentially expressed genes, including 10 upregulated genes and 10 downregulated genes; C abundance of ZNF529-AS1-related genes associated with biological process; D KEGG pathway annotation

Thresholds of p. adj < 0.05 and FDR < 0.25 were used to screen significantly enriched items. Biological process (BP) association included signal release, response to xenobiotic stimulus, and drug catabolic process; cellular component (CC) included synaptic membrane, postsynaptic membrane, and integral component of synaptic membrane; molecular function (MF) included receptor ligand activity, passive transmembrane transporter activity, and channel activity (Fig. 6C). KEGG pathways included neuroactive ligand‒receptor interaction, bile secretion, and retinol metabolism (Fig. 6D).

The relationship between ZNF529-AS1 expression and immune infiltration in HCC

We systematically analysed the immune cells in the tumour and found that the expression of ZNF529-AS1 was positively correlated with the abundances of Th2 cells, NK CD56 bright cells and TFH cells and was negatively correlated with the abundances of Th17 cells, CD cells, and neutrophils (Fig. 7A). By comparing the proportions of immune cells in the ZNF529-AS1 high-expression group and low-expression group in the TCGA dataset, we found that the levels of Th2 (p < 0.001), TFH (p < 0.01), NK CD56 bright cells (p < 0.01), and cytotoxic cells (p < 0.001) in the high ZNF529-AS1 expression group were higher, while the levels of Treg (p < 0.01), Th17 (p < 0.001), Tgd (p < 0.01), Tcm (p < 0.001), NK cells (p < 0.05), NK CD56 bright cells (p < 0.01) and DCs (p < 0.001) in the low ZNF529-AS1 expression group were higher (Fig. 7B). We further analysed the correlations between ZNF529-AS1 expression and immune cell markers, including B cells, CD8+ T cells, M1/M2 macrophages, tumour-associated macrophages (TAMs), neutrophils, natural killer cells and dendritic cells, in ccRCC using the GEPIA database. In addition, some different types of functional T cells, such as Th1, Th2, Th9, Th17, Th22, Tfh, exhausted T cells and Treg cells (Fig. 7C), were also detected. The expression of immune checkpoint proteins can be deregulated in tumours and is an important immune resistance mechanism. Therefore, we analysed the differences in immune checkpoint genes between the ZNF529-AS1 high expression and ZNF529-AS1 low expression groups. Compared with the low expression group, 11 immune checkpoint genes, including PD-1, CTAL4 and LAG-3, were significantly upregulated, and 2 immune checkpoint genes were significantly downregulated in the high expression group (Fig. 7D). Finally, the TME composition of HCC samples was analysed by ESTIMATE. The results showed that compared with the high ZNF529-AS1 expression group, the low ZNF529-AS1 expression group had a higher immune score (p = 0.38) (Fig. 7E), stromal score (p < 0.01) (Fig. 7F) and ESTIMATE score (p = 0.063) (Fig. 7G) and lower tumour purity (p = 0.063) (Fig. 7H), indicating that ZNF529-AS1 had a good ability to predict the TME components of HCC patients.

Fig. 7
figure 7

Immune infiltration analysis of ZNF529-AS1 in the HCC microenvironment. A Correlation analysis between ZNF529-AS1 and immune cells; B Wilcoxon test analysis of TIICs of HCC patients with high expression and low expression of ZNF529-AS1; C and D differential analysis of immune function and immune checkpoints between HCC patients with high expression and low expression of ZNF529-AS1; differences in the E immune score, F stromal score, G ESTIMATE score, and H tumour purity between the ZNF529-AS1 high expression and low expression groups (*p < 0.05; **p < 0.01; ***p < 0.001)

Correlation analysis of ZNF529-AS1 in immunotherapy

Given the important role of immunotherapy in HCC, we analysed the correlation between the ZNF529-AS1 expression level and the pathway enrichment scores predicted by immunotherapy as well as key steps of the cancer-immunity cycle (Fig. 8C, D). The results showed that in the high ZNF529-AS1 expression group, risk scores were positively correlated with most immune pathways, including DNA replication, cell cycle, Fanconi anaemia pathway, homologous recombination, mismatch repair, nucleotide excision repair, oocyte meiosis, viral carcinogenesis, base excision repair, p53 signalling pathway, spliceosome, progesterone-mediated oocyte maturation, microRNAs in cancer, pyrimidine and RNA degradation (Fig. 8A, B). In addition, the high expression of ZNF529-AS1 was correlated with MDSC recruitment and monocyte recruitment at Step 4, a key step in the cancer-immunity cycle (Additional file 1: Fig. S3).

Fig. 8
figure 8

Evaluation of ZNF529-AS1 for predicting treatment response in HCC patients. A Difference analysis between ZNF529-AS1 and clinical response to HCC immunotherapy; B correlation between the high expression group and low expression group in the enrichment score of the immunotherapy prediction path; C correlation between ZNF529-AS1 and the enrichment score of the immunotherapy prediction pathways; D correlations between ZNF529-AS1 and various steps of the cancer-immunity cycle

ZNF529-AS1 regulated HCC cell invasion by targeting FBXO31

The expression level of ZNF529-AS1 was relatively high in HCC tissues and HCC cells (Fig. 9A). Transwell assays showed that low ZNF529-AS1 expression inhibited the invasion and migration of tumour cells (Fig. 9B). After knockdown of ZNF529-AS1, the expression level of the downstream target gene FBXO31 decreased (Fig. 9C, D). It was inferred that FBXO31 may be a downstream target of ZNF529-AS1.

Fig. 9
figure 9

ZNF529-AS1 affects the invasion and migration of hepatocellular carcinoma. A Differential expression of ZNF529-AS1 in hepatocellular carcinoma cells and tissue was determined by qPCR assay. B The effect of ZNF529-AS1 on invasion and migration was determined by the Transwell assay. C and D Effects on downstream genes were determined by qPCR and WB assays after knocking down ZNF529-AS1. (*p < 0.05; **p < 0.01; ***p < 0.001)

Discussion

HCC is a common malignant tumour of the digestive tract and has an insidious onset and high postoperative recurrence and metastasis rates, leading to high mortality [29, 30]. As an important treatment method for cancer, immunotherapy has become the main force in the treatment of advanced HCC. However, many patients are not sensitive to immunotherapy, requiring us to identify potential prognostic and therapeutic targets.

In this study, we identified a novel tumour-associated molecule, ZNF529-AS1, in the TCGA database. ZNF529-AS1 is differentially expressed in various tumours, including HCC. At present, there are few studies on ZNF529-AS1. According to previous reports, ZNF family members play a regulatory role in tumorigenesis. For example, ZNF212 and ZNF367 can regulate the apoptosis of ovarian tumour cells and the proliferation of colorectal cancer cells, respectively [31, 32]. In the ZNF protein family, ZNF 471, ZNF1 and ZNF20 are involved in tumour progression through various pathways, including prognostic features, immune infiltration, TME, and tumour-related signalling pathways [33,34,35]. ZNF529-AS1 is the antisense strand of ZNF529, and its partial sequence is structurally opposite to that of the ZNF family, so it may have the opposite tumour regulation effect. We found that ZNF529-AS1 is upregulated in a variety of cancers in a TCGA database analysis. We analysed the role of ZNF529-AS1 in HCC and found that ZNF529-AS1 is closely correlated with T stage, M stage and pathological stage of HCC, indicating that ZNF529-AS1 may be potentially associated with the occurrence and development of HCC. By comparing the prognosis of HCC patients with high expression and low expression of ZNF529-AS1, we found a significant difference between the two groups (the higher the expression level was, the shorter the survival time), and ZNF529-AS1 was an independent factor affecting the prognosis of HCC patients. However, in our study, we used data from the TCGA and GTEx databases to analyse our results. Nonetheless, there are still some shortcomings, because there may be some asymmetry in the data, such as different TNM classification systems and different time courses, such as different TNM classification systems.

We speculate that ZNF529-AS1 may be a potential therapeutic target for HCC. By immunological correlation analysis, we found that ZNF529-AS1 is negatively correlated with immune-related molecules. ZNF529-AS1 inhibits the killing effect of immune cells on tumours and promotes the occurrence and development of HCC in vivo. Moreover, high expression of ZNF529-AS1 can affect the expression of CD80, PD-1 and CD276 in immune checkpoint cells, thereby promoting the occurrence of HCC. ZNF529-AS1 plays an important role in the immune system of HCC, not only affecting the cellular immune process in HCC but also affecting the immune checkpoint activity of HCC. ZNF529-AS1 promotes the occurrence and development of HCC through this dual role. At the same time, through targeted drug sensitivity analysis, we found that the expression of ZNF529-AS1 was negatively correlated with sensitivity to multiple targeted drugs, indicating that ZNF529-AS1 may be a potential therapeutic target for HCC. Knockdown of ZNF529-AS1 expression may improve the sensitivity to targeted drug therapy for HCC, promote the activation of the body’s immune system, inhibit the occurrence and development of HCC, enhance clinical efficacy, and improve patient prognosis. Additionally, this depends on tumour heterogeneity. Some patients were not sensitive to ZNF529-AS1.

The zinc lipoprotein family has an important role in the immune infiltration of tumours. In this study, we made a correlation prediction for ZNF529-AS1 in the immune infiltration of tumours, and the results showed that ZNF529-AS1 could affect the immune infiltration of hepatocellular carcinoma. Moreover, we correlated the functional expression of FBXO31 with ZNF529-AS1 by predicting that FBXO31, as a member of the forkhead box protein family, was reported to play an important role in immune infiltration of tumours[36], such as in hepatocellular carcinoma, lung cancer, breast cancer and colorectal cancer[37,38,39,40], in which forkhead box proteins play an important role in the function of immune infiltration. Second, we also investigated the function of ZNF529-AS1 at the cellular level. qRT‒PCR showed that ZNF529-AS1 was relatively highly expressed in HCC cells and tissues. Moreover, interfering with ZNF529-AS1 can significantly inhibit the invasion and migration of HCC cells. Next, to identify the downstream targets of ZNF529-AS1, we predicted that ZNF529-AS1 can regulate hsa-miR-561-5p by using the STARBASE database. hsa-miR-561-5p can be closely related to FBXO31 based on the relationship in the database, so we speculate that ZNF529-AS1 can regulate the expression of FBXO31. To verify the relationship between ZNF529-AS1 and FBXO31. We also verified the correlation between ZNF529-AS1 and FBXO31 by using the GEPIA database, and the results showed that ZNF529-AS1 was correlated with FBXO31 (p = 0.042). Therefore, we speculate that the expression of FBXO31 is correlated with ZNF529-AS1 and that FBXO31 is a downstream target of ZNF529-AS1 (Additional file 1: Fig. S4A–C). The expression of FBXO31 was downregulated after interfering with ZNF529-AS1 in HCC cells. ZNF529-AS1 may affect the invasion and migration of HCC cells by regulating FBXO31 expression. The regulatory mechanism of ZNF529-AS1 in HCC needs to be further explored to better understand its role in HCC.

Summary

As a tumour-associated molecule, ZNF529-AS1 is highly expressed in HCC and is closely associated with the prognosis and immune infiltration of HCC. ZNF529-AS1 can affect the invasion and migration of HCC cells by regulating the expression of the downstream target FBXO31. ZNF529-AS1 may be a potential therapeutic target for HCC. The specific molecular mechanism of ZNF529-AS1 in HCC needs to be further studied.

Availability of data and materials

Simulated datasets and the code used to produce the results presented in this paper are available at https://www.scidb.cn/s/VF7fmi.

References

  1. Shangguan D, Meng L, Cao ZC, Xiao Z, Fang X, Li Y, Cardona D, Witek RP, Liu C, Tan W. Identification of liver cancer-specific aptamers using whole live cells. Anal Chem. 2008;80:721–8. https://doi.org/10.1021/ac701962v.

    Article  CAS  PubMed  Google Scholar 

  2. Sumazin P, Chen Y, Treviño LR, Sarabia SF, Hampton OA, Patel K, Mistretta T-A, Zorman B, Thompson P, Heczey A, et al. Genomic analysis of hepatoblastoma identifies distinct molecular and prognostic subgroups. Hepatology (Baltimore, MD). 2017;65:104–21. https://doi.org/10.1002/hep.28888.

    Article  CAS  PubMed  Google Scholar 

  3. 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:209–49. https://doi.org/10.3322/caac.21660.

    Article  PubMed  Google Scholar 

  4. Wallace MC, Preen D, Jeffrey GP, Adams LA. The evolving epidemiology of hepatocellular carcinoma: a global perspective. Expert Rev Gastroenterol Hepatol. 2015;9:765–79. https://doi.org/10.1586/17474124.2015.1028363.

    Article  CAS  PubMed  Google Scholar 

  5. Ioannou GN. Epidemiology and risk-stratification of NAFLD-associated HCC. J Hepatol. 2021;75:1476–84. https://doi.org/10.1016/j.jhep.2021.08.012.

    Article  PubMed  Google Scholar 

  6. Reig M, Forner A, Rimola J, et al. BCLC strategy for prognosis prediction and treatment recommendation: the 2022 update. J Hepatol. 2022;76(3):681–93. https://doi.org/10.1016/j.jhep.2021.11.018.

    Article  PubMed  Google Scholar 

  7. Xu X-F, Xing H, Han J, Li Z-L, Lau W-Y, Zhou Y-H, Gu W-M, Wang H, Chen T-H, Zeng Y-Y, et al. Risk factors, patterns, and outcomes of late recurrence after liver resection for hepatocellular carcinoma: a multicenter study from China. JAMA Surg. 2019;154:209–17. https://doi.org/10.1001/jamasurg.2018.4334.

    Article  PubMed  Google Scholar 

  8. Verna EC, Patel YA, Aggarwal A, Desai AP, Frenette C, Pillai AA, Salgia R, Seetharam A, Sharma P, Sherman C, et al. Liver transplantation for hepatocellular carcinoma: management after the transplant. Am J Transplant. 2020;20:333–47. https://doi.org/10.1111/ajt.15697.

    Article  PubMed  Google Scholar 

  9. Zhang Y, Zhang X, Kuang M, Yu J. Emerging insights on immunotherapy in liver cancer. Antioxid Redox Signal. 2022. https://doi.org/10.1089/ars.2022.0047.

    Article  PubMed  Google Scholar 

  10. Guo HJ, Wang LJ, Wang C, Guo DZ, Xu BH, Guo XQ, Li H. Identification of an Apis cerana zinc finger protein 41 gene and its involvement in the oxidative stress response. Arch Insect Biochem Physiol. 2021;108:e21830. https://doi.org/10.1002/arch.21830.

    Article  CAS  PubMed  Google Scholar 

  11. He L, Fan X, Li Y, Chen M, Cui B, Chen G, Dai Y, Zhou D, Hu X, Lin H. Overexpression of zinc finger protein 384 (ZNF 384), a poor prognostic predictor, promotes cell growth by upregulating the expression of Cyclin D1 in hepatocellular carcinoma. Cell Death Dis. 2019;10:444. https://doi.org/10.1038/s41419-019-1681-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Xie W, Qiao X, Shang L, Dou J, Yang X, Qiao S, Wu Y. Knockdown of ZNF233 suppresses hepatocellular carcinoma cell proliferation and tumorigenesis. Gene. 2018;679:179–85. https://doi.org/10.1016/j.gene.2018.08.070.

    Article  CAS  PubMed  Google Scholar 

  13. Fan L, Wang J, Deng P, Wang Y, Zhang A, Yang M, Zeng G. Foxhead box D1 promotes the partial epithelial-to-mesenchymal transition of laryngeal squamous cell carcinoma cells via transcriptionally activating the expression of zinc finger protein 532. Bioengineered. 2022;13:3057–69. https://doi.org/10.1080/21655979.2021.2024978.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sokolov A, Carlin DE, Paull EO, Baertsch R, Stuart JM. Pathway-based genomics prediction using generalized elastic net. PLoS Comput Biol. 2016;12:e1004790. https://doi.org/10.1371/journal.pcbi.1004790.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Auslander N, Zhang G, Lee JS, Frederick DT, Miao B, Moll T, Tian T, Wei Z, Madan S, Sullivan RJ, et al. Robust prediction of response to immune checkpoint blockade therapy in metastatic melanoma. Nat Med. 2018;24:1545–9. https://doi.org/10.1038/s41591-018-0157-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Liu J, Lichtenberg T, Hoadley KA, Poisson LM, Lazar AJ, Cherniack AD, Kovatich AJ, Benz CC, Levine DA, Lee AV, et al. An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell. 2018;173:400–16. https://doi.org/10.1016/j.cell.2018.02.052.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:E90-W102. https://doi.org/10.1093/nar/gkx247.

    Article  CAS  Google Scholar 

  18. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:1–25.

    Article  Google Scholar 

  19. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7. https://doi.org/10.1089/omi.2011.0118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Shankavaram UT, Varma S, Kane D, Sunshine M, Chary KK, Reinhold WC, Pommier Y, Weinstein JN. Cell miner: a relational database and query tool for the NCI-60 cancer cell lines. BMC Genom. 2009;10:277. https://doi.org/10.1186/1471-2164-10-277.

    Article  CAS  Google Scholar 

  21. Reinhold WC, Sunshine M, Liu H, Varma S, Kohn KW, Morris J, Doroshow J, Pommier Y. Cell miner: a web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set. Cancer Res. 2012;72:3499–511. https://doi.org/10.1158/0008-5472.CAN-12-1370.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Xu L, Deng C, Pang B, Zhang X, Liu W, Liao G, Yuan H, Cheng P, Li F, Long Z, et al. TIP: a web server for resolving tumor immunophenotype profiling. Cancer Res. 2018;78:6575–80. https://doi.org/10.1158/0008-5472.CAN-18-0689.

    Article  CAS  PubMed  Google Scholar 

  23. Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, Kadel EE, Koeppen H, Astarita JL, Cubas R, et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554:544–8. https://doi.org/10.1038/nature25501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Auslander N, Zhang G, Lee JS, Frederick DT, Miao B, Moll T, Tian T, Wei Z, Madan S, Sullivan RJ, et al. Publisher Correction: robust prediction of response to immune checkpoint blockade therapy in metastatic melanoma. Nat Med. 1942;2018:24. https://doi.org/10.1038/s41591-018-0247-8.

    Article  CAS  Google Scholar 

  25. Zhang D, Ma X, Sun W, Cui P, Lu Z. Down-regulated FSTL5 promotes cell proliferation and survival by affecting Wnt/β-catenin signaling in hepatocellular carcinoma. Int J Clin Exp Pathol. 2015;8:3386–94.

    PubMed  PubMed Central  Google Scholar 

  26. Zhang DY, Ma SS, Sun W-L, Lv XCH, Lu Z. KIF4A as a novel prognostic biomarker in cholangiocarcinoma. Medicine. 2021;100:e26130. https://doi.org/10.1097/MD.0000000000026130.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhang D-Y, Lei J-S, Sun W-L, Wang D-D, Lu Z. Follistatin Like 5 (FSTL5) inhibits epithelial to mesenchymal transition in hepatocellular carcinoma. Chin Med J. 2020;133:1798–804. https://doi.org/10.1097/CM9.0000000000000847.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang D-Y, Sun W-L, Ma X, Zhang P, Wu W, Wu H, Zhou S, Lu Z. Up-regulated FSTL5 inhibits invasion of hepatocellular carcinoma through the Wnt/β-catenin/YAP pathway. Int J Clin Exp Pathol. 2017;10:10325–33.

    PubMed  PubMed Central  Google Scholar 

  29. Rutherford MJ, Arnold M, Bardot A, Ferlay J, De P, Tervonen H, Little A, Bucher O, St Jacques N, Gavin A, et al. Comparison of liver cancer incidence and survival by subtypes across seven high-income countries. Int J Cancer. 2021;149:2020–31. https://doi.org/10.1002/ijc.33767.

    Article  CAS  PubMed  Google Scholar 

  30. Bao X, Chen L, Liu Y, Sheng H, Wang K, Luo Y, Qin T, Liu Y, Qiu Y. Treatment of liver cancer: role of the traditional mongolian medicine. Evid Based Complement Alternat Med. 2022;2022:6535977. https://doi.org/10.1155/2022/6535977.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Qu D, Zou X, Liu Z. Propofol modulates glycolysis reprogramming of ovarian tumor via restraining circular RNA-zinc finger RNA-binding protein/microRNA-212-5p/superoxide dismutase 2 axis. Bioengineered. 2022;13:11881–92. https://doi.org/10.1080/21655979.2022.2063649.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Lei T, Gao Y, Duan Y, Cui C, Zhang L, Si M. Inhibition of zinc finger protein 367 exerts a tumor suppressive role in colorectal cancer by affecting the activation of oncogenic YAP1 signaling. Environ Toxicol. 2021;36:2278–90. https://doi.org/10.1002/tox.23341.

    Article  CAS  PubMed  Google Scholar 

  33. Cao L, Wang S, Zhang Y, Wong KC, Nakatsu G, Wang X, Wong S, Ji J, Yu J. Correction to: Zinc-finger protein 471 suppresses gastric cancer through transcriptionally repressing downstream oncogenic PLS3 and TFAP2A. Oncogene. 2022;41:3298–9. https://doi.org/10.1038/s41388-022-02216-x.

    Article  CAS  PubMed  Google Scholar 

  34. Cui SB, Wang TX, Liu ZW, Yan JY, Zhang K. Zinc finger protein A20 regulates the development and progression of osteoarthritis by affecting the activity of NF-κB p65. Immunopharmacol Immunotoxicol. 2021;43:713–23. https://doi.org/10.1080/08923973.2021.1970764.

    Article  CAS  PubMed  Google Scholar 

  35. Piskor EM, Ross J, Möröy T, Kosan C. Myc-interacting zinc finger protein 1 (Miz-1) is essential to maintain homeostasis and immunocompetence of the B cell lineage. Biology (Basel). 2022;11:504. https://doi.org/10.3390/biology11040504.

    Article  CAS  PubMed  Google Scholar 

  36. Bach DH, Long NP, Luu TT, Anh NH, Kwon SW, Lee SK. The dominant role of forkhead box proteins in cancer. Int J Mol Sci. 2018;19:3279. https://doi.org/10.3390/ijms19103279.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Hu G, Yan Z, Zhang C, Cheng M, Yan Y, Wang Y, et al. FOXM1 promotes hepatocellular carcinoma progression by regulating KIF4A expression. J Exp Clin Cancer Res. 2019;38:188. https://doi.org/10.1186/s13046-019-1202-3.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Giannos P, Kechagias KS, Gal A. Identification of prognostic gene biomarkers in non-small cell lung cancer progression by integrated bioinformatics analysis. Biology (Basel). 2021;10:1200. https://doi.org/10.3390/biology10111200.

    Article  CAS  PubMed  Google Scholar 

  39. Fan D, Zeng C, Wang S, Han J, Zhu L, Zhao H, et al. Forkhead box P3 promotes breast cancer cell apoptosis by regulating programmed cell death 4 expression. Oncol Lett. 2020;20:292. https://doi.org/10.3892/ol.2020.12155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Tang H, Zheng J, Bai X, Yue KL, Liang JH, Li DY, et al. Forkhead box Q1 is critical to angiogenesis and macrophage recruitment of colorectal cancer. Front Oncol. 2020;10:564298. https://doi.org/10.3389/fonc.2020.564298.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to the Natural Science Research Program of Anhui University and the Bengbu Medical College for funding and to our teachers and peers for their contributions and efforts in completing the thesis.

Funding

This work was supported by the Anhui University Natural Science Research Project (Dengyong Zhang, Grant Number [KJ2021A0731]).

Author information

Authors and Affiliations

Authors

Contributions

ZL and D-YZ contributed to the study conception and design. Material preparation, data assembly and analysis were performed by YM, S-SM, GZ, WS and ZL. The first draft of the manuscript was written by D-YZ, YM and WS. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zheng Lu or Dengyong Zhang.

Ethics declarations

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

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.

The supplementary figure 1-figure 4.

Additional file 2.

Supplementary table 1-table 2.

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

Ma, Y., Sun, Wl., Ma, S.S. et al. LincRNA ZNF529-AS1 inhibits hepatocellular carcinoma via FBXO31 and predicts the prognosis of hepatocellular carcinoma patients. BMC Bioinformatics 24, 54 (2023). https://doi.org/10.1186/s12859-023-05189-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12859-023-05189-0

Keywords