Skip to main content

Integrated analysis of ceRNA network reveals potential prognostic Hint1-related lncRNAs involved in hepatocellular carcinoma progression

Abstract

Background

Hint1 is a novel tumor suppressor gene, and inactivation of its expression is closely associated with the carcinogenesis of a variety of malignancies. The effects of Hint1 deficiency on the competing endogenous RNA (ceRNA) regulatory network in the context of HCC remains to be fully characterized. This study aims to explore Hint1-related hub lncRNAs in HCC and to establish a reliable prognostic model for HCC patients based on these hub lncRNAs.

Methods

lncRNA + mRNA microarray was used to identify differentially expressed (DE) lncRNAs and mRNAs in Huh7 cells before and after Hint1 knockdown. A Hint1-related ceRNA network was mapped by bioinformation technology. The DEmRNAs in the network were analyzed via GO and KEGG enrichment analyses. Hub DElncRNAs associated with HCC patient prognosis were then detected through univariate and multivariate Cox regression analyses and were incorporated into a prognostic model. The prognostic value of this model was then assessed through the use of Kaplan-Meier curves, time-related ROC analyses, and nomograms. We also utilized Kaplan-Meier curves to validate the relationship between hub lncRNAs and the overall survival (OS) of HCC patients. Finally, A Hint1-related core ceRNA network based on the hub DElncRNAs and DEmRNAs was mapped.

Results

We identified 417 differentially expressed DElncRNAs and 2096 DEmRNAs in Huh7 cells before and after Hint1 knockdown. Three hub DElncRNAs (LINC00324, SNHG3, and DIO3OS) in the Hint1-associated ceRNA network were screened out using univariate and multivariate Cox regression analyses. A hepatocellular carcinoma (HCC) prognostic risk-scoring model and nomogram were constructed using these three hub lncRNAs, and it was confirmed that the risk score of the model could be used as an independent predictor of HCC prognosis. A Hint1-related core ceRNA network based on the hub DElncRNAs and DEmRNAs was also mapped.

Conclusion

We constructed a reliable prognostic model for HCC patients based on three Hint1-related hub lncRNAs, and we believe these three hub lncRNAs may play critical roles in hepatocarcinogenesis, and progression.

Introduction

Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide [1]. More than 700,000 people are diagnosed with HCC each year, and about 600,000 people die due to HCC and related complications [2]. Although great progress has been made in the research and treatment of liver cancer and the 10-year survival rate after radical resection has been improved to about 30%, about 50–70% of patients with radical resection of small liver cancer still have recurrence and metastasis within 5 years [3]. Therefore, exploring the carcinogenesis and progression mechanism of HCC, finding an effective way to inhibit metastasis and recurrence of HCC, and investigating the predictive factors of postoperative HCC recurrence are key to further improving the survival rate of patients with liver cancer.

At one time, non-coding RNAs (ncRNAs) were considered transcriptional noise with no biological activity [4]. However, in recent years, it has become clear that ncRNAs play a critical role in many cellular processes, such as epigenetic regulation, chromosome remodeling, transcriptional regulation, and post-translation modification [5,6,7]. Its abnormal expression can lead to the disorder of important biological processes such as cell proliferation, apoptosis, and cell migration and invasion, which leads to malignant cell transformation and even tumorigenesis [8,9,10]. In these ncRNAs, microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) have attracted attention [4]. lncRNAs can regulate the inhibitory effects of miRNAs on downstream target genes through a mechanism involving competing endogenous RNA (ceRNA) and then participate in the carcinogenesis and progression of HCC [11,12,13].

Histidine triad nucleotide-binding protein 1 (HINT1) is a tumor suppressor, which belongs to the triple histidine superfamily [14, 15]. The epigenetic inactivation of the Hint1 gene is a common and potential early event in the carcinogenesis of HCC [16]. HINT1 protein can regulate the transcription of a variety of cancer-related genes by directly binding a variety of transcription factors, including AP1, MITF, and USF2, thus affecting the carcinogenesis and progression of tumors [14, 15, 17]. However, the effects of Hint1 deficiency on the ceRNA regulatory network of tumor cells need to be clarified.

Therefore, in this study, we used siRNA to knock down the expression of Hint1 in Huh7 cells, comprehensively analyzed differentially expressed (DE) lncRNAs and mRNAs before and after knockdown using bioinformatics technology, and predicted the possible ceRNA network. Then, combined with TCGA (The Cancer Genome Atlas) database, a DElncRNAs diagnosis model was constructed, and hub DElncRNAs that may be regulated by Hint1 were screened out. These hub DElncRNAs may help further clarify the molecular mechanism underlying the promotion of HCC progression by Hint1. We also expect these hub DElncRNAs to become promising biomarkers for predicting HCC progression and prognosis.

Materials and methods

siRNA transfection

Hint1 siRNA (sc-92005) and interference control siRNA (sc-37007) were purchased from Santa Cruz Biotechnology (CA, USA). Transfection was performed with lipofectamine RNAiMax according to the manufacturer’s instructions.

Total RNA extraction and quantitative real-time PCR

Total RNA was extracted as previously described [15]. Quantitative RT-PCR was performed using the SYBR green mix kit (Roche, Basel, Switzerland) and specific primers (Supplementary Table S1) in the ABI StepOne Plus real-time PCR system (Applied Biosystems, Foster City, CA, USA) in a volume of 20 μL. The relative quantification of mRNA levels was computed using the 2-ΔΔCT method.

Western blot

The relative levels of the target and control β-actin proteins were determined by Western blotting using cell lysates (30 μg/lane). Anti-HINT1 (ab124912, 1:1000) and anti-β-actin (ab8226, 1:500) were purchased from Abcam (Cambridge, MA, USA).

Gene expression microarray profiling

TRIzol reagent was used to extract the total RNA of Huh7 cells treated with siHINT1 and siCTRL. The cDNA obtained by reverse transcription was hybridized to Agilent human lncRNA + mRNA Array V4.0(Boao, Beijing, China), which was designed with four identical arrays per slide (4×180K format) with each array containing probes interrogating about 41,000 human lncRNAs and about 34,000 human mRNAs. Microarray hybridization and computational analysis was performed as previously described [18].

Construction of a ceRNA network

The possible target miRNA of DElncRNAs was predicted using the miRcode [19] database (http://www.mircode.org/). Three online databases including miRDB (http://mirdb.org), miRTarBase [20] (http://mirtarbase.mbc.nctu.edu.tw/php/index.php), and TargetScan [21] (http://www.ta-rgetscan.og) were used to search the targeted mRNAs of miRNAs. To improve the reliability of the results, only the miRNA-mRNA relational pairs found in all three databases were selected as candidate genes for the construction of ceRNA networks. Finally, Cytoscape software (version: 3.8.2) was used to visualize the ceRNA network based on the interactions among predicted DElncRNAs, predicted miRNAs, and predicted DEmRNAs [22].

KEGG pathway and GO enrichment analyses

GO and KEGG enrichment analyses of DEmRNAs within the ceRNA network were performed using the “enrichgo” and “enrichkegg” functions within the “cluster profiler” package of the R software (version: 4.0.4); both P and Q values ≤ 0.05 were considered significantly enriched terms.

Data from TCGA database

The RNA-sequencing (RNA-seq) data of clinical specimens and corresponding clinical data of hepatocellular carcinoma (LIHC) patients were downloaded from the TCGA database (https://portal.g-dc.cancer.gov/). The “limma” package of R software was used to normalize the RNA-seq data, and significantly differentially expressed genes with an average of less than 1 were excluded. Then, using the “createDataPartition” function in the R software “caret” package, the TCGA case cohort was randomly split (50:50%) into a training group and testing group for subsequent analysis.

Construction of prognosis prediction model and nomogram

Using the “coxph” and “survdiff” functions of the “survival” package of R software, the survival correlated univariate Cox regression analysis of ceRNA related DElncRNA was carried out, and significant candidate hub DElncRNAs were screened out applying P<0.05 as a cutoff. Then, a multivariate stepwise Cox proportional risk regression model was constructed using the “caret,” “glmnet”, and “survminer” packages of R software, and the risk score model of each sample was obtained according to the following formula:

$$\mathrm{Risk}\ \mathrm{score}=\sum \left(\upbeta \ast \mathrm{ExpDElncRNAs}\right)$$

where β represents the coefficient (coef) of gene in the multivariate stepwise Cox analysis., and Exp represented the gene expression level. Then the receiver operating characteristic curve (ROC) of the training and testing groups was analyzed using the “survivalroc” package of R software to confirm and evaluate the predictive ability of the above models. Finally, the nomogram function of the “RMS” package of the R software was used to generate the nomogram.

Validation of prognostic significance

The “Kaplan Meier plotter” [23] online tool (https://kmplot.com/analysis/) was used to verify the prognostic value of hub DElncRNAs in HCC patients.

Statistical analysis

Numerical data are presented as mean ± standard error of the mean (SEM). Differences between groups were determined using a Student’s t test with a significance set of P≤0.05. Pearson’s chi-square tests were used to analyze the correlation between risk score and clinical features; the median risk score was used as the cutoff, and P ≤ 0.05 was considered significant.

Results

Knockdown of Hint1 and microarray profiling

The study design is illustrated in Fig. 1. siRNA was used for knockdown of Hint1 expression in Huh7 cells. As shown in Fig. 2, the expression of Hint1 decreased significantly by treatment of siHINT1. After confirming this knockdown effect, we ran a lncRNA + mRNA expression microarray with total RNA extracted from Hint1-deleted Huh7 cells (Huh7 cells treated with siCTRL as control). Microarray raw data have been uploaded to the GEO (Gene Expression Omnibus) database (GEO Accession: GSE177624). After preliminary processing of the microarray data, we identified 417 DElncRNAs (231 were upregulated and 186 were downregulated, Fig. 2C) and 2096 DEmRNAs (899 were upregulated and 1197 were downregulated, Fig. 2D) with |FC|≥2 and FDR<0.05 as the cutoff values. We selected three lncRNAs with the most significant changes in expression, as well as mRNAs, for PCR verification. There was a good correlation with the microarray results (Fig. 2E–G).

Fig. 1
figure 1

Flow chart of Hint1-related ceRNA network analysis. DElncRNAs, differentially expressed RNAs; DEmRNAs, differentially expressed mRNAs; ceRNA, competing endogenous RNA; GO, Gene Ontology; KEGG, Kyoto Encyclopaedia of Genes and Genomes; TCGA, The Cancer Genome Atlas; LIHC, liver hepatocellular carcinoma; RNA-seq, RNA-sequencing; ROC, receiver operating characteristic curve

Fig. 2
figure 2

Microarray profiling and validation in Hint1-deficient Huh7 cells. A mRNA level of Hint1 in both siCTRL- and siHINT1-treated Huh7 cells determined by qRT -PCR (normalized to the mRNA level of β-ACTIN, data was presented as mean ± SEM; n = 3). B Protein level of HINT1 and β-ACTIN in both siCTRL- and siHINT1-treated Huh7 cells determined by western blot. Duplicate lncRNAs (C) and mRNAs (D) expression microarray analysis of Huh7 cells treated with either siCTRL or siHINT1 were subjected to log2-transformed gene expression followed by K-means clustering (fold-change (FC) ≥ 2, while false discovery rate (FDR) ≤ 0.05, up- and downregulated lncRNAs and mRNAs are indicated in red and green, respectively). E qRT-PCR was performed on twelve lncRNAs and mRNAs with the highest levels of elevation and decrease in the microarray analysis and normalized to mRNA level of β-ACTIN (FCs were computed in terms of percentage fold change of the siHINT1-treated transcriptome relative to the siCTRL-treated transcriptome, data was presented as mean ± SEM; n = 3); Mini heat map representing the selected twelve lncRNA (F) and mRNA (G) in microarray analysis. ***P ≤ 0.001

Prediction of ceRNA network regulated by Hint1

We first used the miRcode database to predict the potential miRNAs that may interact with DElncRNAs. We then utilized three online databases (miRDB, miRTarBase, and TargetScan) to predict the potential target mRNAs of miRNAs. Next, we intersected these predicted mRNAs with the DEmRNAs to obtain DEmRNAs that may be involved in the ceRNA network. We ultimately identified 135 predicted DElncRNAs and 185 predicted DEmRNAs that may have regulatory relationships, as well as 40 predicted miRNAs that may be involved in ceRNA networks. We visualized the final ceRNA network using Cytoscape software; the network included 360 nodes with 1738 interactions (Fig. S1).

GO and KEGG enrichment analyses

The 185 predicted DEmRNAs that may be involved in pertinent biological processes and pathways were further investigated. Using the “Cluster Profiler” package of the R software, we carried out GO and KEGG enrichment analyses, and used the “Enrich Plot” package to visually display the results. The top 10 of 44 significant GO terms and the top 10 of 17 significant KEGG pathways are listed in Table 1. Of these terms, “RNA polymerase II transcription regulator complex,” “intrinsic apoptotic signaling pathway,” and “misfolded protein binding” were the top GO terms in the cellular component (CC), biological process (BP), and molecular function (MF) categories, respectively. “Kaposi sarcoma-associated herpesvirus infection,” “microRNAs in cancer,” “measles,” “p53 signaling pathway,” and “Pertussis” were the first five KEGG pathways.

Table 1 Top 10 significant GO and KEGG terms identified by enrichment analysis

Screening of prognosis-related lncRNA

The RNA-seq data and their corresponding clinical data in the LIHC subset of the TCGA database were randomly divided into a training group (186 cases) and testing group (184 cases) for the following process. Considering that lncRNA was in the upstream part of the ceRNA network and it is the main effector of miRNAs and mRNAs, we carried out survival-related univariate Cox regression analysis of the 135 predicted DElncRNAs in the ceRNA network with the training group data. Six survival-related candidate hub DElncRNAs were screened out (Fig. S2A). Subsequently, these six candidate hub DElncRNAs were used to perform multiple stepwise Cox regression, to further investigate their potential effects on the survival time and clinical outcomes of patients. Finally, three hub DElncRNAs, LINC00324, SNHG3, and DIO3OS, were identified as survival-related independent predictors for HCC patients (Fig. S2B, Table 2).

Table 2 Three prognosis-associated hub DElncRNAs identified by stepwise multivariate Cox regression analysis in the training group

Construction and analysis of prognosis-related lncRNA risk-scoring model

Next, we used the three hub lncRNAs and their coefficient values obtained in multiple stepwise Cox regression analysis to construct a multivariate Cox proportional hazards regression model, and calculated a risk score for each patient in the training group according to the following formula:

$$Risk\;score=\left(-0.2320\ast ExpLINC00324\right)+\left(0.0782\ast ExpSNHG3\right)+\left(-0.6007\ast ExpDIO3OS\right)$$

Furthermore, the training group was divided into low-risk and high-risk subgroups according to the median risk score. Then the differences in overall survival (OS) between these two subgroups were compared, the high-risk subgroup had a significantly OS than the low-risk subgroup (P<0.001, Fig. 3A). To further evaluate the utility of the three hub DElncRNAs as prognostic biomarkers, time-dependent ROC analysis was performed. The area under the ROC curve (AUC) of the lncRNA risk scoring model was 0.775, indicating moderate diagnostic performance (Fig. 3B–D).

Fig. 3
figure 3

Risk score analysis of three-hub lncRNAs prognostic model in the training group of TCGA LIHC cohort (n = 186). A Survival curve for low- and high-risk subgroups. B ROC curves for forecasting OS based on risk score. C Expression heat map of the three-hub lncRNAs for low- and high-risk subgroups. D Risk score distribution, and the survival status for low- and high-risk subgroups

In addition, we evaluated whether the risk-scoring model was against the testing cohort. Similar to the results for the training group, the OS of patients with a high risk score in the testing group was worse than that of patients with a low risk score (Fig. 4A). And the model had good sensitivity and specificity, with an AUC of 0.722 (Fig. 4B–D).

Fig. 4
figure 4

Risk score analysis of three-hub lncRNAs prognostic model in the testing group of TCGA LIHC cohort (n = 184). A Survival curve for low- and high-risk subgroups. B ROC curves for forecasting OS based on risk score. C Expression heat map of the three-hub lncRNAs for low- and high-risk subgroups. D Risk score distribution, and the survival status for low- and high-risk subgroups

Meanwhile, we also analyzed the association between risk score and other clinicopathological parameters in the testing group. As was shown in Table S2, the level of risk score was correlated with tumor stage and tumor sizes.

Nomogram drawn based on the three hub DElncRNAs

To predict the prognosis of HCC patients more intuitively, the “rms” package of the R software was used to draw a nomogram based on the three hub DElncRNAs (Fig. 5). Using that nomogram, the estimated survival rate of patients at 1, 3, and 5 years could be obtained by drawing a vertical line between the total integral axis and each pre-posterior axis, which could help doctors make clinical decisions for HCC patients. In addition, the prognostic significance of different clinical features in the testing group was also evaluated using univariate Cox regression analysis. Tumor stage, tumor size, and risk score were correlated with OS in patients with HCC (Table 3); however, a multiple Cox regression analysis showed that risk score was the only independent prognostic factor related to OS (Table 3).

Fig. 5
figure 5

Nomogram for predicting 1-, 3-, and 5-year OS of LIHC patients in the TCGA cohort

Table 3 The prognostic value of different clinical parameters in the testing group

Validation of the prognostic value of hub DElncRNAs

To further evaluate the prognostic value of the three hub lncRNAs in HCC, the “Kaplan Meier-plotter” online tool was used to analyze the relationship between each hub lncRNA and OS and draw Kaplan-Meier curves. All three hub lncRNAs were closely related to OS in HCC patients (Fig. S3).

Core ceRNA network related to hub DElncRNAs

Finally, on the basis of the original ceRNA network, Cytoscape software was used to screen and draw the core ceRNA network based on the three hub DElncRNAs. As shown in Fig. 6, the core network was composed of 172 nodes and 251 interactions, involving 143 DEmRNAs and 26 miRNAs that may participate in the network. This provides an important reference for further study of the specific mechanism of the ceRNA network regulated by Hint1.

Fig. 6
figure 6

Hint1-related core ceRNA network created by hub DElncRNAs. Red diamonds represent 3 hub DElncRNAs (2.2% of all DElncRNAs); blue ellipses represent 143 DEmRNAs (77.3% of all DEmRNAs); yellow triangles represent 26 predicted miRNAs (65.0% of all predicted miRNAs)

Discussion

HCC is the sixth most common malignant tumor worldwide [24]. Although hepatectomy, liver transplantation, and other therapies can improve the prognosis of HCC patients to a certain extent, it is still the third leading cause of cancer-related death in the world due to its high rate of invasion, metastasis, and postoperative recurrence [25]. Therefore, it is important to elucidate the molecular mechanism of the occurrence and progression of HCC and to find an effective target for its diagnosis and treatment. In recent years, with the development of high-throughput screening and bioinformatics technology, increasing numbers of lncRNAs have been found to be involved in gene regulation, particularly in ceRNA networks [26]. Meanwhile, the important role of ceRNA regulatory networks in the oncogenesis and progression of malignant tumors has gradually been recognized [27]. For example, lncRNA HIF1A-AS1 [28] and ZEB1-AS1 [29] were recently proved to be associated with HCC progression via ceRNA pattern, while lncRNA TTN-AS1 was demonstrated to act as a tumor promoter in gallbladder carcinoma by sponging miR-107 and upregulating HMGA1 [30].

Hint1 is a newly discovered tumor suppressor gene. In 2003, Su et al. [31] first found that Hint1 may have a tumor inhibitory effect by establishing a gene knockdown mouse model. Subsequent studies have found that its expression is decreased in many human malignant tumors, including HCC, non-small-cell lung cancer, and colon cancer [31,32,33]. HINT1 protein can affect the transcription of downstream target genes by interacting with a variety of transcription factors, thereby regulating the proliferation, apoptosis, invasion, and migration of tumor cells [34]. Specifically, our team found that Hint1 can inhibit the proliferation of rectal cancer cells by directly binding to Posh-JNK2 complexes, inhibiting the transcriptional activity of AP-1, which is an important cancer transcription factor [15]. Motzik et al. [17] confirmed that HINT1 protein can promote apoptosis of melanoma cells by directly binding to the DNA binding protein MITF. Considering the role of Hint1 in the regulation of transcription, we speculated that it may indirectly regulate the ceRNA network by affecting the transcription of lncRNAs, thus influencing the progress of malignant tumors. Therefore, in this study, we first used a lncRNA + mRNA microarray to detect the DElncRNAs and DEmRNAs in Huh7 cells before and after Hint1 knockdown. This identified 417 DElncRNAs and 2096 DEmRNAs, from which we constructed the ceRNA network.

As the proteins translated by mRNA in ceRNA regulatory networks are the final functional implementers, GO and KEGG enrichment analyses were conducted on the DEmRNAs in the ceRNA network, to better understand how Hint1-related ceRNA networks may be involved in the carcinogenesis and progression of HCC. GO enrichment analysis showed that DEmRNAs in the network were mainly involved in “RNA polymerase II transcription regulator complex” terms, indicating that the transcriptional regulatory activity of Hint1 may partly be based on the ceRNA network. In addition, DEmRNAs were highly enriched in “intrinsic apoptotic signaling pathway” terms, which to some extent explains the effects of Hint1 on apoptosis of HCC cells reported by Hsieh et al. [35], namely, that Hint1 may partly regulate the apoptosis of HCC cells by regulating the ceRNA network. In KEGG enrichment analysis, the DEmRNAs were mainly involved “microRNAs in cancer,” “p53 signaling pathway,” and other cancer-related pathways, which further confirmed that Hint1-related ceRNA regulatory networks played an important role in HCC.

lncRNAs are located in the upstream of ceRNA network and its expression level may be closely related to the prognosis of patients with malignant tumors [6]. Therefore, we used univariate and multivariate stepwise Cox regression analyses to find hub DElncRNAs related to the prognosis of HCC using the TCGA database. LINC00324, SNHG3, and DIO3OS were identified as survival-related hub lncRNAs. Among them, DIO3OS has been reported to inhibit the proliferation and invasion of HCC cells by competitive binding to miR-328 [36]; our results are consistent with that finding. In addition, Zhang et al. [37] reported that SNHG3 can promote the EMT process of HCC cells through competitive binding of miR-128, and then release CD151. However, there are no reports suggesting that LINC00324 is involved in the occurrence or development of HCC, although some studies have reported that it is a risk factor for disease progression in non-small-cell lung cancer and osteosarcoma [38, 39]. The functional differences of LINC00324 between different tissues is worth further investigating. Furthermore, we established a prognostic risk-scoring model of HCC based on three hub DElncRNAs using multivariate stepwise Cox regression analysis. The DElncRNAs had good diagnostic ability and could screen out HCC patients with a poor prognosis. Subsequently, we mapped a nomogram to help us more intuitively predict the OS of HCC patients for 1, 3, and 5 years. Validation tests indicated that all of them were closely related to OS in HCC patients. Finally, we constructed an Hint1-related core ceRNA network based on the three hub DElncRNAs. Interestingly, in this core network, these three hub DElncRNAs regulated 77% of the DEmRNAs through 24 miRNAs, confirming their importance in the network. Among these 24 miRNAs, miR-23b-3p was reported to be involved in the epithelial-mesenchymal transition (EMT) of HCC by targeting EOMES mRNA [40], and miR-10a-5p has also been proved to be associated with migration and invasion in HCC cells [41]. At the same time, the apoptosis-related gene PIM1 [42], proliferation-related gene CCND1 [43], migration- and invasion-related gene BMPR2 [44], and tumor immunocyte infiltration-related gene IFNAR1 also appeared in the core network [45]. This indicates the diversity of hub lncRNAs in the regulation of the tumor malignant phenotype, and also confirms the significant role of the Hint1-related ceRNA network in malignant tumors. Overall, our prognostic model is based on three lncRNAs, which significantly reduces the cost of clinical gene sequencing and is more favorable for clinical application. In addition, our prognostic model provided a complementary perspective on individual tumors and developed an individual scoring system for HCC patients, and the nomogram based on this scoring system may be a promising tool for clinicians in the future.

However, there were some limitations to this study. First, due to the lack of a liver cancer lncRNA database, our prognosis model is only based on data from the TCGA database and has not been verified in a clinical cohort or other databases. Second, although we used several online databases to predict the miRNAs involved in the ceRNA network, miRNAs expression needs to be validated in both in vivo and in vitro experiments.

Conclusions

In conclusion, we systematically investigated the composition of the Hint1-related ceRNA network for HCC and the prognostic value of specific lncRNAs in HCC through a series of bioinformatics analyses. We identified three hub DElncRNAs related to patient prognosis, constructed an HCC prognostic risk score model with these three lncRNAs, and verified that the model could be used as an independent predictor of HCC prognosis. To the best of our knowledge, this is the first study to establish an HCC prognostic model for Hint1-related lncRNAs. Our results may contribute to further research into the molecular mechanisms underlying the promotion of HCC progression by Hint1, and we expect these three lncRNAs will become promising biomarkers for predicting HCC progression and prognosis.

Availability of data and materials

Microarray raw data have been uploaded to the GEO database ( https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE177624). The other data used or analyzed during the study are available from the corresponding author on reasonable request.

References

  1. Villanueva A. Hepatocellular Carcinoma. N Engl J Med. 2019;380:1450–62.

    CAS  PubMed  Google Scholar 

  2. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34.

    PubMed  Google Scholar 

  3. Chen Z, Xie H, Hu M, Huang T, Hu Y, Sang N, et al. Recent progress in treatment of hepatocellular carcinoma. Am J Cancer Res. 2020;10:2993–3036.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Anastasiadou E, Jacob LS, Slack FJ. Non-coding RNA networks in cancer. Nat Rev Cancer. 2018;18:5–18.

    CAS  PubMed  Google Scholar 

  5. Dossin F, Pinheiro I, Żylicz JJ, Roensch J, Collombet S, Le Saux A, et al. SPEN integrates transcriptional and epigenetic control of X-inactivation. Nature. 2020;578:455–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Mercer TR, Mattick JS. Structure and function of long noncoding RNAs in epigenetic regulation. Nat Struct Mol Biol. 2013;20:300–7.

    CAS  PubMed  Google Scholar 

  7. Ni W, Yao S, Zhou Y, Liu Y, Huang P, Zhou A, et al. Long noncoding RNA GAS5 inhibits progression of colorectal cancer by interacting with and triggering YAP phosphorylation and degradation and is negatively regulated by the m (6) A reader YTHDF3. Mol Cancer. 2019;18:143.

    PubMed  PubMed Central  Google Scholar 

  8. Luo H, Xu C, Le W, Ge B, Wang T. lncRNA CASC11 promotes cancer cell proliferation in bladder cancer through miRNA-150. J Cell Biochem. 2019;120:13487–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Wang K, Long B, Zhou LY, Liu F, Zhou QY, Liu CY, et al. CARL lncRNA inhibits anoxia-induced mitochondrial fission and apoptosis in cardiomyocytes by impairing miR-539-dependent PHB2 downregulation. Nat Commun. 2014;5:3596.

    PubMed  Google Scholar 

  10. Zhuang C, Ma Q, Zhuang C, Ye J, Zhang F, Gui Y. LncRNA GClnc1 promotes proliferation and invasion of bladder cancer through activation of MYC. FASEB J. 2019;33:11045–59.

    CAS  PubMed  Google Scholar 

  11. Han TS, Hur K, Cho HS, Ban HS. Epigenetic Associations between lncRNA/circRNA and miRNA in Hepatocellular Carcinoma. Cancers (Basel). 2020;12:2622.

    CAS  Google Scholar 

  12. Wang H, Huo X, Yang XR, He J, Cheng L, Wang N, et al. STAT3-mediated upregulation of lncRNA HOXD-AS1 as a ceRNA facilitates liver cancer metastasis by regulating SOX4. Mol Cancer. 2017;16:136.

    PubMed  PubMed Central  Google Scholar 

  13. Yan Y, Lu Y, Mao K, Zhang M, Liu H, Zhou Q, et al. Identification and validation of a prognostic four-genes signature for hepatocellular carcinoma: integrated ceRNA network analysis. Hepatol Int. 2019;13:618–30.

    PubMed  Google Scholar 

  14. Wang L, Li H, Zhang Y, Santella RM, Weinstein IB. HINT1 inhibits beta-catenin/TCF4, USF2 and NFkappaB activity in human hepatoma cells. Int J Cancer. 2009;124:1526–34.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Wang L, Zhang Y, Li H, Xu Z, Santella RM, Weinstein IB. Hint1 inhibits growth and activator protein-1 activity in human colon cancer cells. Cancer Res. 2007;67:4700–8.

    CAS  PubMed  Google Scholar 

  16. Zhang YJ, Li H, Wu HC, Shen J, Wang L, Yu MW, et al. Silencing of Hint1, a novel tumor suppressor gene, by promoter hypermethylation in hepatocellular carcinoma. Cancer Lett. 2009;275:277–84.

    CAS  PubMed  Google Scholar 

  17. Motzik A, Amir E, Erlich T, Wang J, Kim BG, Han JM, et al. Post-translational modification of HINT1 mediates activation of MITF transcriptional activity in human melanoma cells. Oncogene. 2017;36:4732–8.

    CAS  PubMed  Google Scholar 

  18. Peng QL, Zhang YM, Yang HB, Shu XM, Lu X, Wang GC. Transcriptomic profiling of long non-coding RNAs in dermatomyositis by microarray analysis. Sci Rep. 2016;6:32818.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. 2012;28:2062–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011;39:D163–9.

    CAS  PubMed  Google Scholar 

  21. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.

    PubMed Central  Google Scholar 

  22. Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48:D127–31.

    CAS  PubMed  Google Scholar 

  23. Menyhárt O, Nagy Á, Győrffy B. Determining consistent prognostic biomarkers of overall survival and vascular invasion in hepatocellular carcinoma. R Soc Open Sci. 2018;5:181006.

    PubMed  PubMed Central  Google Scholar 

  24. Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol. 2019;16:589–604.

    PubMed  PubMed Central  Google Scholar 

  25. Anwanwan D, Singh SK, Singh S, Saikam V, Singh R. Challenges in liver cancer and possible treatment approaches. Biochim Biophys Acta Rev Cancer. 1873;2020:188314.

    Google Scholar 

  26. Lou W, Ding B, Fu P. Pseudogene-Derived lncRNAs and Their miRNA Sponging Mechanism in Human Cancer. Front Cell Dev Biol. 2020;8:85.

    PubMed  PubMed Central  Google Scholar 

  27. Zhang C, Ke Y, Liu X, Wang X, Li Y, Zhou J, et al. High Expression of the Long Noncoding RNA SNHG15 in Cancer Tissue Samples Predicts an Unfavorable Prognosis of Cancer Patients: A Meta-Analysis. J Oncol. 2020;2020:3417036.

    PubMed  PubMed Central  Google Scholar 

  28. Hong F, Gao Y, Li Y, Zheng L, Xu F, Li X. Inhibition of HIF1A-AS1 promoted starvation-induced hepatocellular carcinoma cell apoptosis by reducing HIF-1α/mTOR-mediated autophagy. World J Surg Oncol. 2020;18:113.

    PubMed  PubMed Central  Google Scholar 

  29. Xue S, Lu F, Sun C, Zhao J, Zhen H, Li X. LncRNA ZEB1-AS1 regulates hepatocellular carcinoma progression by targeting miR-23c. World J Surg Oncol. 2021;19:121.

    PubMed  PubMed Central  Google Scholar 

  30. Lin Z, Li Y, Shao R, Hu Y, Gao H. LncRNA TTN-AS1 acts as a tumor promoter in gallbladder carcinoma by regulating miR-107/HMGA1 axis. World J Surg Oncol 2021; 19: 163.31. Su T, Suzui M, Wang L, Lin CS, Xing WQ, Weinstein IB. Deletion of histidine triad nucleotide-binding protein 1/PKC-interacting protein in mice enhances cell growth and carcinogenesis. Proc Natl Acad Sci U S A. 2003;100:7824–9.

    PubMed  PubMed Central  Google Scholar 

  31. Jacobson BA, Ahmad Z, Chen S, Waldusky G, Dillenburg M, Stoian E, et al. 4Ei-10 interdiction of oncogenic cap-mediated translation as therapy for non-small cell lung cancer. Invest New Drugs. 2021;39:636–43.

    CAS  PubMed  Google Scholar 

  32. Cen B, Deguchi A, Weinstein IB. Activation of protein kinase G Increases the expression of p21CIP1, p27KIP1, and histidine triad protein 1 through Sp1. Cancer Res. 2008;68:5355–62.

    CAS  PubMed  Google Scholar 

  33. Wu XS, Bao TH, Ke Y, Sun DY, Shi ZT, Tang HR, et al. Hint1 suppresses migration and invasion of hepatocellular carcinoma cells in vitro by modulating girdin activity. Tumour Biol. 2016;37:14711–9.

    CAS  PubMed  Google Scholar 

  34. Duan DD, Xie H, Shi HF, Huang WW, Ding F, Hong JK, et al. Hint1 Overexpression Inhibits the Cell Cycle and Induces Cell Apoptosis in Human Osteosarcoma Cells. Onco Targets Ther. 2020;13:8223–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Hsieh SY, Hsu CY, He JR, Liu CL, Lo SJ, Chen YC, et al. Identifying apoptosis-evasion proteins/pathways in human hepatoma cells via induction of cellular hormesis by UV irradiation. J Proteome Res. 2009;8:3977–86.

    CAS  PubMed  Google Scholar 

  36. Wang Z, Song L, Ye Y, Long LW, Noncoding RNA. DIO3OS Hinders Cell Malignant Behaviors of Hepatocellular Carcinoma Cells Through the microRNA-328/Hhip Axis. Cancer Manag Res. 2020;12:3903–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Zhang PF, Wang F, Wu J, Wu Y, Huang W, Liu D, et al. LncRNA SNHG3 induces EMT and sorafenib resistance by modulating the miR-128/CD151 pathway in hepatocellular carcinoma. J Cell Physiol. 2019;234:2788–94.

    CAS  PubMed  Google Scholar 

  38. Pan ZH, Guo XQ, Shan J, Luo SX. LINC00324 exerts tumor-promoting functions in lung adenocarcinoma via targeting miR-615-5p/AKT1 axis. Eur Rev Med Pharmacol Sci. 2018;22:8333–42.

    PubMed  Google Scholar 

  39. Wu S, Gu Z, Wu Y, Wu W, Mao B, Zhao S. LINC00324 accelerates the proliferation and migration of osteosarcoma through regulating WDR66. J Cell Physiol. 2020;235:339–48.

    CAS  PubMed  Google Scholar 

  40. Yan Z, Lijuan Y, Yinhang W, Yin J, Jiamin X, Wei W, et al. Screening and analysis of RNAs associated with activated memory CD4 and CD8 T cells in liver cancer. World J Surg Oncol. 2022;20(1):2.

    PubMed  PubMed Central  Google Scholar 

  41. Wu Y, Zhou Y, Huan L, Xu L, Shen M, Huang S, et al. LncRNA MIR22HG inhibits growth, migration and invasion through regulating the miR-10a-5p/NCOR2 axis in hepatocellular carcinoma cells. Cancer Sci. 2019;110(3):973–84.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Brasó-Maristany F, Filosto S, Catchpole S, Marlow R, Quist J, Francesch-Domenech E, et al. PIM1 kinase regulates cell death, tumor growth and chemotherapy response in triple-negative breast cancer. Nat Med. 2016;22:1303–13.

    PubMed  PubMed Central  Google Scholar 

  43. Gennaro VJ, Stanek TJ, Peck AR, Sun Y, Wang F, Qie S, et al. Control of CCND1 ubiquitylation by the catalytic SAGA subunit USP22 is essential for cell cycle progression through G1 in cancer cells. Proc Natl Acad Sci U S A. 2018;115:E9298–307.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Pickup MW, Hover LD, Polikowsky ER, Chytil A, Gorska AE, Novitskiy SV, et al. BMPR2 loss in fibroblasts promotes mammary carcinoma metastasis via increased inflammation. Mol Oncol. 2015;9:179–91.

    CAS  PubMed  Google Scholar 

  45. Lu C, Klement JD, Ibrahim ML, Xiao W, Redd PS, Nayak-Kapoor A, et al. Type I interferon suppresses tumor growth through activating the STAT3-granzyme B pathway in tumor-infiltrating cytotoxic T lymphocytes. J Immunother Cancer. 2019;7:157.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the grants from the National Natural Science Foundation of China [grant numbers 81660399, 81860423, 82103173]; the Innovative Research Team Project of Yunnan Province [grant number 2015HC033]; the Yunnan Provincial Academician Workstation of Xiaoping Chen [grant number 2017IC018]; the Breeding Program for Major Scientific and Technological Achievements of Kunming Medical University [grant number CGYP201607]; the Medical Leading Talent Project of Yunnan Province [grant number L201622]; and Yunnan Provincial Clinical Center of Hepato-biliary-pancreatic Diseases [no specific number] to L. W; and the Leading Academic and Technical Young and Mid-aged Program of Kunming Medical University [grant number 60118260108]; and the Educational Research and Educational Reform Program of Kunming Medical University [grand number 2019-JY-Z-12] to Y. K.

Author information

Authors and Affiliations

Authors

Contributions

C.Z., Y.K., X.L., X.W., and T.B. designed this study. C.Z., Y.K., and X.L. analyzed the data and wrote the manuscript. X.W., W.L., Y.H., and T.B. collected the data. L.W. revised the manuscript. All authors approved the final version for submission.

Corresponding author

Correspondence to Lin Wang.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of the Second Affiliated Hospital of Kunming Medical University, China. Signed written informed consents were obtained from the patients and/or guardians.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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: Table S1.

The sequences of primers for quantitative PCR (qPCR). Table S2. Correlation between risk score and other clinicopathological parameters in testing group.

Additional file 2: Fig

ure S1. Hint1-related ceRNA network. Blue ellipses represent 185 DEmRNAs; red diamonds represent 135 DElncRNAs; yellow triangles represent 40 predicted miRNAs.

Additional file 3: Figure S2.

(A) Univariate Cox regression analysis for identification of hub lncRNAs in the training group of TCGA LIHC cohort. (B) Multivariate Cox regression analysis to identify prognosis related hub lncRNAs in the training group of TCGA LIHC cohort.

Additional file 4: Figure S3.

The expression of (A) DIO3OS (B) LINC00324 and (C) SNHG3 were associated with the OS in LIHC patients in Kaplan Meier Plotter server.

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

Verify currency and authenticity via CrossMark

Cite this article

Zhang, C., Bao, T., Ke, Y. et al. Integrated analysis of ceRNA network reveals potential prognostic Hint1-related lncRNAs involved in hepatocellular carcinoma progression. World J Surg Onc 20, 67 (2022). https://doi.org/10.1186/s12957-022-02535-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12957-022-02535-z

Keywords

  • Hint1
  • ceRNA
  • Hepatocellular carcinoma
  • lncRNA
  • Risk score