Skip to main content

The expression profiles and prognostic values of HSP70s in hepatocellular carcinoma

Abstract

Background

The HSP70 family of heat shock protein plays a critical role in protein synthesis and transport to maintain protein homeostasis. Several studies have indicated that HSP70s are related to the development and occurrence of various cancers.

Methods

The relationship between the overall survival rate of hepatocellular carcinoma patients and the expression of 14 HSP70s from multiple databases, such as TCGA, ONCOMINE, cBioPortal was investigated. Western Blot and PCR were used to evaluate HSPA4 and HSPA14 expressions in various HCC cells to identify suitable cell lines for further experiments .Wound-healing assays, Transwell assays and EdU assays were used to verify the effects of HSPA4 and HSPA14 on the function of hepatocellular carcinoma cells, and statistical analysis was performed.

Results

Hepatocellular carcinoma tissues significantly expressed the 14 HSP70s compared to the normal samples. Besides, the high HSPA1A, HSPA1B, HSPA4, HSPA5, HSPA8, HSPA13, and HSPA14 expressions were inversely associated with the overall survival rate of patients, tumor grade, and cancer stage. A PPI regulatory network was constructed using the 14 HSP70s proteins with HSPA5 and HSPA8 at the network center. Univariate and multivariate analyses showed that HSPA4 and HSPA14 could be independent risk factors for the prognosis of hepatocellular carcinoma patients. Cell experiments have also confirmed that reducing HSPA4 and HSPA14 expressions can inhibit the invasion, metastasis, and proliferation of hepatocellular carcinoma cells.

Conclusions

Therefore, the HSP70s significantly influence the occurrence and development of hepatocellular carcinoma. For instance, HSPA4 and HSPA14 can be novel therapeutic targets and prognostic biomarkers for hepatocellular carcinoma.

Introduction

Hepatocellular carcinoma (HCC) is a common and deadly cancer. About 800,000 new cases and over 700,000 HCC-associated deaths are reported yearly [1]. Surgical treatment, including hepatectomy and liver transplantation (LT), is the main HCC treatment. Moreover, sorafenib is the only systematic therapy for HCC [2]. Despite significant progress in hepatocellular carcinoma treatment in the past, Its prognosis remains poor odue to its high invasiveness and metastasis rate. Moreover, the molecular characteristics of HCC are unknown. Understanding the potential pathogenesis and etiology of HCC can help find advanced biomarkers for HCC treatment and diagnosis, thus improving the prognosis and individualized patient treatment.

Increased conformational unstable proteins cause continuous proliferation and genetic instability of cancer cells. Also, excessive protein synthesis and metabolism lead to a stress state in cancer cells that must be permanently compensated. Therefore, tumor cells depend on the maintenance of protein homeostasis via the molecular chaperones. HSP70s are ubiquitous molecular chaperones essential for folding and remodeling of all cellular proteins. HSP70s help misfolded proteins to fold and refold correctly. HSP70s are essential in all protein synthesis and degradation stages, thus important in maintaining protein homeostasis [3]. Cancer cells can release HSP70s, which influence malignant characteristics via receptor-mediated signal transduction [4]. Therefore, HSP70s can influence HCC occurrence and development via several mechanisms. Some studies have revealed that the HSP70 expression is upregulated in some human malignant tumors, such as liver, breast, pancreatic, colorectal, and prostate cancers [5,6,7,8,9]. Therefore, a comprehensive study of various HSP70s in HCC is needed to clarify the molecular mechanisms involved in HCC development, providing potential therapeutic targets and prognostic markers of HCC.

Presently, 14 HSP70s members have been identified (HSPA1A, HSPA1B, HSPA1L, HSPA2, HSPA4, HSPA4L, HSPA5, HSPA6, HSPA8, HSPA9, HSPA12B, HSPA12A, HSPA13 and HSPA14). HSP70s have critical roles in tumor cells [10, 11]. Some studies have found abnormal expression of some HSP70s and their prognostic value. For instance, HSPA1A promotes hepatoma cell proliferation through TLR2 and TLR4 signaling pathways[12], and HSPA1L/HIF-1 α/GP78 significantly influences tumor progression [13]. HSPA1B (HSP70-2) down-regulation can reduce cell proliferation and tumor growth [14]. HSPA2 overexpression is associated with tumor angiogenesis, invasive progression, and poor prognosis [15, 16]. Previous studies showed that HSPA5 (GRP78) increases FAT10 expression via NF-κB pathway to promote HCC proliferation [17]. HSPA5 also accelerates HCC invasion and metastasis by activating the Wnt pathway via the LRP6 [18]. Similarly, HSPA4 [19], HSPA8 [20, 21], HSPA9 [22, 23], and HSPA12B [24, 25] are associated with tumor occurrence and development. However, the role and relationship of various HSP70s in HCC occurrence and development are unknown. Several DNA and RNA data have been produced due to the continuous progress in microarray technology. Several gene expressions or copy number variations published in the online database can be analyzed to test whether the expression of various HSP70 proteins is related to clinical parameters, such as the OS of HCC patients, and explore the function of HSP70s in HCC occurrence and development.

Materials and methods

ONCOMINE database

ONCOMINE database (www.oncomine.org) is a comprehensive online cancer microarray database for DNA or RNA sequence analysis that facilitates discovery from the gene-wide expression analyses [26]. Transcriptional expressions of 14 HSP70s in various cancer tissues and their corresponding adjacent normal control samples were obtained from ONCOMINE database.

UALCAN

UALCAN (http://ualcan.path.uab.edu) is an interactive web resource based on level 3 RNA-seq and clinical data of 31 cancer types from the TCGA database [27]. UALCAN was used to analyze the mRNA expressions of 14 HSPA70s in HCC tissues and their relationship with clinicopathologic parameters. Students’ t-test was used to compare the transcriptional expression difference. P < 0.01 was considered statistically significant.

cBioPortal

cBioPortal (www.cbioportal.org) is an online open-access website resource. It is used to explore, visualize, and analyze multidimensional cancer genomics data [28]. The genomic profiles of 14 HSP70s, including mutations, putative copy-number alterations from GISTIC, and mRNA Expression z-Scores (RNASeq V2 RSEM) with a z-score threshold ± 1.8 were analyzed.

The Cancer Genome Atlas (TCGA) database

TCGA [29] is a coordinated and comprehensive tool that ameliorates diagnostic methods and treatment criteria for cancer prevention. Over 30 human tumor sequencing and pathological data can be analyzed on TCGA. Transcriptional group RNA sequencing data of HCC samples, including the transcriptional group data and corresponding clinical information of 377 HCC patients, were downloaded from the TCGA (https://cancergenome.nih.gov/). A total of 89 patients with less than three months overall survival, or inadequate clinical information, or no R0 surgery in the TCGA database were not excluded to avoid the impact of non-related confounding factors. The data of the remaining 288 patients were used for univariate and multivariate analyses. The results are shown in Table 1.

Table 1 Basic characteristics of 288 HCC patients

Functional enrichment analysis

STRING v11 database (https://string-db.org/) was used to identify 30 frequently neighbor genes of HSP70s. DAVID database (https://david.ncifcrf.gov/summary.jsp) was then used to analyze the function of HSP70s and the 30 genes significantly related to HSP70s, such as cellular components (CC), biological processes (BP), molecular functions (MF), and KEGG analysis.

Cell culture and transfection

HCC human cell lines, including HCCLM3, Huh-7, SMCC7721, and MHCC97H, were sourced from the Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences (Shanghai, China). Cell Bank confirmed the cells using short tandem repeat profiling. The Cells were cultured as previously described. Lipofectamine™ 3000 Transfection Reagent (Invitrogen, Carlsbad, USA) was used for transfection [17].

Western blot analysis and Real-time quantitative polymerase chain reaction (qRT-PCR)

Western blotting and qRT-PCR were performed as earlier mentioned [30]. The specific primers used for PCR amplification are shown in Additional file 1: Table S1.

Wound-healing assays

Use the marker pen to draw a horizontal line evenly on the back of the 6-well plate, about every 0.5–1 cm, across the well. The cells were cultured in 6-well plates to 80–90%. Scratched the horizontal line perpendicular to the back with the pipette tip. Washed the cells gently with PBS for 3 times and added serum-free medium. Continue culturing in incubator at 37 ℃. Observed and took pictures at the time point of 0 h, 24 h, 48 h [30].

Transwell assays

5 × 104 cells were suspended in serum-free medium and placed in the upper chamber. After 24 h of culture, the unmigrated cells on the upper surface of the membrane were carefully removed with cotton swabs, the cells on the lower surface of the membrane were fixed and stained with 0.1% crystal violet, observed and photographed [18].

EdU assay

The cells were treated with 150 µl 5-ethynyl-20-deoxyuridine (EdU) (Ribobio, Guangzhou, China) for 2 h. Washed the cell with PBS for three times. After cells were treated with 2 mg/ml glycine and 0.5 % TritionX-100 in turn,the cells were stained with 300 µl of 1×Apollo® reaction cocktail for 30 min. Then, the cells were stained with 300 µl of Hoechst 33,342 (5 µg/ml) for 30 min. Washed the cells with PBS after each step. Observed and took pcitures  [31].

Statistical methods

R v3.6.1(https://www.r-project.org/) and SPSS 26.0(SPSS Inc., Chicago, IL) were used for all statistical analyses. Univariate Cox algorithm was used to analyze and evaluate the effects of clinical parameters and HSP70s mRNA expression on the survival of HCC patients. Those with P < 0.1 were used for follow-up analysis. A Multivariate Cox algorithm was used to analyze the relationship between HSP70s mRNA expression and patient survival. P < 0.05 was considered statistically significant. Each experiment was repeated thrice. All values were expressed as an average of ± SD. T-test was used to analyze all experimental results, * P < 0.05,** P < 0.01 and *** P < 0.001.

Results

HSP70 expressions in the HCC and normal tissues

ONCOMINE and TCGA databases were used to analyze the mRNA expression of HSP70s to explore the different prognosis and potential therapeutic value of HSP70s in HCC patients. ONCOMINE database datasets showed that the mRNA expression of HSPA4/4L/13/14 was significantly higher in liver cancer tissues. Roessler Liver 2 dataset showed significant HSPA4/13/14 overexpression in HCC tissues compared with the normal tissues with multiple changes of 2.202 times (p = 1.9E−72), 2.119 times (p = 9.4E−49), and 2.099 times (p = 7.76E−58), respectively. Wurmbach dataset observed that the mRNA expressions of HSPA4 and HSPA4L increased 2.053 times (p = 3.49E−06) and 5.204 times (p = 5.26E−07), respectively in HCC samples. Roessler dataset revealed that the mRNA expression of HSPA4 increased 2.492 times in HCC tissues (p = 1.21E−08) (Fig. 1; Table 2). The expression profile data of LIHC patients in the TCGA database were also analyzed. The mRNA expression level of HSP70s was significantly upregulated in HCC tissues (all p < 0.05) (Fig. 2).

Table 2 There was a significant change in the expression of HSP70s at the transcriptional level of (ONCOMINE) between hepatocellular carcinoma and normal liver tissues
Fig. 1
figure1

The transcriptional level of HSP70s in Oncomine dataset of different types of cancer

Fig. 2
figure2

HSP70s is differentially expressed between HCC and normal tissues. The mRNA of all HSP70s are higher than that of normal tissues (p < 0.05)

Prognostic value of differentially expressed HSP70s in hepatocellular carcinoma patients

The correlation between HSP70 expressions in HCC patients and the overall survival of HCC patients in the TCGA database was further analyzed. The mRNA expression of nearly half of HSP70s was associated with the prognosis of HCC patients. HSPA1A (p = 0.048), HSPA1B (p = 0.038), HSPA4 (p = 0.019), HSPA5 (p = 0.016), HSPA8 (p = 0.011), HSPA13 (p = 0.018) and HSPA14 (p = 0.00069)were associated with shorter OS. However, HSPA1L, HSPA2, HSPA4L, HSPA6, HSPA9, HSPA12A, and HSPA12B mRNA expressions were not associated with the prognosis of HCCC patients (Fig. 3). Therefore, the expression levels of HSPA1A, HSPA1B, HSPA4, HSPA5, HSPA8, HSPA13, and HSPA14 are significantly associated with the prognosis of HCC patients and can be used as potential biomarkers to predict prognosis.

Fig. 3
figure3

Prognostic value of different HSP70s mRNA expression in hepatocellular carcinoma. The mRNA levels of HSPA1A (A), HSPA1B (B), HSPA4 (E), HSPA5 (G), HSPA8 (I), HSPA13 (M), HSPA14 (N) were negatively correlated with the shorter OS in patients with HCC. However, the expression levels of HSPA1L (C), HSPA2 (D), HSPA4L (F), HSPA6 (H), HSPA9 (J), HSPA12A (K) and HSPA12B (L) mRNA were not related to the prognosis of patients with HCC (p > 0.05)

Relationship between HSP70 expressions and clinicopathological parameters in HCC patients

The relationship between HSP70expressions (novel prognostic markers) and clinicopathological parameters of HCC patients, including individual tumor grade and cancer stage, was further analyzed. The expression of seven HSP70s with prognostic value was significantly associated with the individual cancer stage and tumor grade. Furthermore, the HSP70 mRNA expression levels were higher at the high cancer stage (Fig. 4). The tumor grade increased with the increased mRNA expression. HSPA4 and HSPA14 expressions were significantly associated with various tumor grades (Fig. 5). Therefore, the mRNA expression of seven HSP70s with prognostic value is positively correlated with clinicopathological parameters.

Fig. 4
figure4

The relationship between the level of 7 HSP70s and neoplasm histologic grade of patients with HCC. With the increase of mRNA expression, their tumor grade tended to be higher. HSPA4 and HSPA14 were most significant among different histologic grade

Fig. 5
figure5

The relationship between the levels of 7 HSP70s and neoplasm disease stage of patients with HCC. With the increase of the mRNA expression, their tumor grade tended to be higher. HSPA4 and HSPA14 were most significant among different disease stage

HSP70s correlation analysis and gene mutation

The interaction between 14 HSP70s was also analyzed by assessing HSP70expressions to understand their relationship. HSPA1A and HSPA1B, HSPA4 and HSPA8, HSPA4 and HSPA9, HSPA6 and HSPA1A, HSPA6 and HSPA1B, HSPA13 and HSPA14 were significantly correlated (Fig. 6A). The genetic variation of HSP70s was then analyzed. HCC patients had a high mutation rate of HSP70s. A total of 202 of the 345 sequenced HCC patients showed genetic variations (58.55%). The mutation rates of HSPA1L, HSPA4, HSPA6 and HSPA9 were more than 10%. High and Multiple Alterations were the most common mutations (Fig. 6B).

Fig. 6
figure6

Analysis of correlation among HSP70s, gene mutation and functional enrichment analysis of HSP70. A HSPA1A and HSPA1B, HSPA4 and HSPA8, HSPA4 and HSPA9, HSPA6 and HSPA1A, HSPA6 and HSPA1B, HSPA13 and HSPA14 had very significant correlation (correlation coefficient > 0.3). B The mutation rate of HSP70s gene was 58.55%. the mutation rates of HSPA1L, HSPA4, HSPA6 and HSPA9 are all higher than 10%, mRNA High and Amplification are the most common mutation types. Thirty neighbor genes significantly related to HSP70s were found through STRING database. Using Cytoscape 3.7.1 to display PPI results through inter-gene connectivity, we found that HSPA8 and HSPA5 played a central role in this network. Then, through the analysis, annotation, visualization and comprehensive discovery of GO (E) and KEGG (F) analysis in DAVID database

Construction of an interaction network of HSP70s and 30 related genes and prediction of their functions and pathways

HSP70 correlation analysis showed that 30 neighbor genes were significantly related to HSP70s in the STRING database (Fig. 6C). Cytoscape 3.7.1 was used to visualize the PPI results via the inter-gene connectivity (Fig. 6D). HSPA8 and HSPA5 played a central part in this network, and the HSP90 and HSP70 families were closely related. Analysis, annotation, visualization, and comprehensive discovery of GO and KEGG enrichment in the DAVID database were used to predict the function of HSP70s and explore 30 adjacent genes. The 30 genes were significantly related to HSP70s and their potential molecular mechanisms. “Protein folding”, “U2-type catalytic step2 spliceosome” and “Unfolded Protein binding” are the most common biological terms in biological processes, molecular functions, and cellular components. Most genes were upregulated in Spliceosome and Protein processing pathways in the endoplasmic reticulum (Fig. 6F).

Independent prognostic value of HSP70 expressions in HCC patients

Independent prognostic value evaluation of HSP70 mRNA expression and OS in HCC patients showed HSPA1A, HSPA1B, HSPA4, HSPA5, HSPA8, HSPA13, and HSPA14 were significantly associated with HCC prognosis. Univariate analysis indicated that the TNM stage, Pathological stage, Child pugh classification grade with higher Tumor status, Pharmaceutical Therapy were associated with shorter OS in HCC patients. The higher mRNA expression of HSPA4L, HSPA4, HSPA8, HSPA12A, and HSPA14 was associated with shorter OS in HCC patients (Fig. 7A). Further multivariate analysis revealed that the higher mRNA expression of HSPA4 and HSPA14 was associated with shorter OS in HCC patients (Fig. 7B). Therefore, the expression levels of HSPA4 and HSPA14 are independent prognostic factors of OS in HCC patients.

Fig. 7
figure7

Univariate and multivariate Cox regression analyses of the association between clinicopathological factors (including HSP70s) and overall survival of patients

Interfering with the expression of HSPA4 or HSPA14 inhibits the proliferation, invasion and metastasis of hepatoma cells

The high HSPA4 expression was associated with vascular invasion, while the high HSPA14 expression with the higher T stage of the tumor (Fig. 8A and B). Therefore, the high HSPA4 and HSPA14 expressions are associated with tumor proliferation, invasion and metastasis, indicating that their high expressions could influence the prognosis of HCC patients by accelerating tumor proliferation, invasion, and metastasis of tumor. Cell function experiments were used to confirm the conjecture. Western Blot and qPCR were used to evaluate HSPA4 and HSPA14 expressions in normal liver cells (HL-7702) and HCC cells. The results revealed that HCCLM3 cells demonstrate relatively high levels of HSPA4 and HSPA14 expression (Fig. 8C and D). Therefore, HSPA4 and HSPA14 expressions in HCCLM3 cells were knocked down via interference fragments. PCR and Western blotting results revealed that HSPA4 and HSPA14 expressions significantly decreased in HCCLM3 cells (Fig. 8E–H). The NC group was the negative control group. A scratch experiment was then conducted using HCCLM3 cells. The scratch spacing significantly decreased in the NC group after 48 h, but there was no significant change in HSPA4 and HSPA14 knockdown groups, indicating that the migration ability of hepatoma cells decreased (Fig. 9A and B), similar to the migration and invasion assays. (Fig. 9C and D) Besides, Edu proliferation assays suggested that cell proliferation ability was significantly decreased after HSPA4 and HSPA14 expression interferences(Fig. 10A and B).

Fig. 8
figure8

Expression of HSPA4 and HSPA14 in normal liver cells and HCC cells. A The expression of HSPA4 was positively correlated with vascular invasion. B The expression level of HSPA14 was positively correlated with higher T stage. C and D The results of Western Blot and qPCR revealed that HCCLM3 cells demonstrate relatively high levels of HSPA4 and HSPA14 expression. E and F The expression of HSPA4 in siHSPA4 is less than that in NC group, and the expression of HSPA14 was less than that in NC group(p < 0.01), the NC group was the negative control group

Fig. 9
figure9

Knockdown of HSPA4 and HSPA14 expression inhibits HCC invasion and metastasis in vitro. The distance of scratche in siHSPA4 group is no significant change but that in NC group decreased (A). while that in siHSPA14 group is no significant change, but in NC group decreased (B). The number of cells in siHSPA4 group (C) and siHSPA14 group (D) were significantly less than that in NC group in invasion and migration assay (p < 0.01), the NC group was the negative control group

Fig. 10
figure10

Knockdown of HSPA4 and HSPA14 expression inhibits HCC proliferation in vitro. The number of proliferative cells in siHSPA4 group was less than that in NC group (A); The number of proliferative cells in siHSPA14 group was less than that in NC group (B), the NC group was the negative control group

Discussion

HCC is a serious malignant tumor in the world due to its complex molecular and cellular heterogeneity. Besides, HCC incidence continues to increase [32]. Over 200 genes related to HCC proliferation, invasion and metastasis have been reported. However, the specific prognostic biomarkers and therapeutic targets are insufficient [33]. Therefore, the screening of HCC molecular biological markers could improve prognosis and reduce mortality.

HSP70s are conserved and inducible heat shock proteins essential in protein homeostasis [34]. HSP70s inhibit non-specific protein aggregation and help proteins to obtain normal structure. HSP70s also have anti-apoptotic properties and are involved in many processes such as tumor cell proliferation, invasion, metastasis, and death [35]. HSP70s helps cells maintain protein homeostasis. The functional enrichment analysis showed that HSP70s are mainly involved in protein folding and unfolded protein reactions. Related HCC research showed that some HSP70s play a crucial role. However, it is necessary to determine the several roles of HSP70s in HCC. This is the first study to systematically investigate the expression, mutation, and prognostic value of various HSP70s in HCC.

Bioinformatics analysis indicated that 14 HSP70s are significantly increased in HCC. The high expression levels of HSPA1A, HSPA1B, HSPA4, HSPA5, HSPA8, HSPA13, and HSPA14 were significantly negatively correlated with overall survival. However, their high expression was positively associated with higher tumor grade and cancer stage. Some studies have shown that HSPA1A can enhance the apoptosis resistance of H22 cells to promote tumor growth [12] by activating NF- κB, and HSPA1A of extracellular to the endogenous ligand of TLR4 and TLR2 [12]. Garg et al. [36] did IHC analysis and RT-PCR and showed that HSPA1B (HSP70-2) is overexpressed in 86% of tumor tissues. Also, down-regulated HSP70-2 expression substantially inhibited cell growth, colony formation, migration, and invasion in vitro and inhibited tumor growth in vivo. Nirmala Jagadish et al. [14] also used RT-PCR and IHC to detect the high HSP 70-2 expressions in 78% of CRC patients. The decreased HSP 70-2 expression causes growth inhibition, colony formation, migration, and CRC cells invasion. Moreover, reduced HSP 70-2 expressions significantly decreased the tumor growth in the COLO205 human transplanted tumor model in mice. Little research has been conducted on the relationship between HSPA13 (STcH) and tumors. Zongguo Yang et al. [37] indicated that HSPA13 is overexpressed in tumor tissues. While HSPA1A, HSPA1B, and HSPA13 were negatively correlated with the overall survival of patients, they were highly expressed in HCC. The higher HSPA1A, HSPA1B, and HSPA13 expressions were positively correlated with higher cancer stage and tumor grade.

In this study, HSPA5 and HSPA8 were highly expressed in HCC, and their high expression was negatively associated with the overall survival of patients. HSPA5 and HSPA8 were in the center of the PPI network, suggesting that they are involved in the occurrence and development of hepatocellular carcinoma via several mechanisms. Therefore, HSPA5 and HSPA8 can be prognostic biomarkers in HCC diagnosis and treatment. A previous study showed that HSPA5 regulates FAT10 expression by modulating the NF- κ B pathway and directly activating the NF- κB pathway. Besides, the study showed that FAT10 is an anti-tumor suppressor p53 gene [17], and HSPA5 regulates HCC invasion and metastasis via the HSPA5-LRP6-HOXB9 axis [18]. Some studies have shown that HSPA8 (HSC 70, HSP 71, HSP 71, or HSP 73) is overexpressed in endometrial carcinoma. siRNA of HSPA8 transfection into endometrial cancer cells significantly inhibited cell proliferation and growth and promoted apoptosis. HSPA8 is key in the occurrence and development of endometrial carcinoma [21]. The above studies have shown that HSP70s promote the occurrence and development of tumor cells, and their high expression is associated with the poor prognosis of tumor patients.

Finally, whether HSP70s can be used as an independent risk factor to predict the prognosis of hepatocellular carcinoma patients was analyzed. Univariate and multivariate analyses showed that could predict the prognosis of HCC patients. High HSPA4 and HSPA14 expressions were also associated with vascular infiltration and higher T stage in HCC patients. HSPA4 and HSPA14 expression interference showed inhibited Cell proliferation, invasion, and metastasis. Zongguo Yang et al. [37] analyzed HSPA4 expressions in 220 HCC patients and found that HSPA4 is overexpressed in tumor tissues. They also showed that HSPA4 is associated with early HCC recurrence. HSPA4 overexpression is associated with the poor prognosis of breast cancer patients [19]. HSPA4 is a potential cancer stem cell (CSC) marker for gastric cancer [38]. While some studies have shown that HSPA14 (HSP70L1) is associated with the overall survival and recurrence of cervical cancer [39], siRNA-mediated inhibition of HSPA14 can reduce the migration, invasion, and transformation activity of NBS-1 overexpressed cells [40]. However, the mechanism of HSPA4 and HSPA14 effect on the occurrence and development of hepatocellular carcinoma is unknown and should be studied.

In conclusion, the differential expression of HSP70s in HCC and its potential diagnostic and prognostic value showed that HSP70s are highly expressed in HCC and are associated with poor prognosis. Further analysis showed that HSPA1A, HSPA1B, HSPA4, HSPA5, HSPA8, HSPA13, and HSPA14 could be used as prognosis biomarkers, and HSPA4 HSPA14 as independent risk factors to determine the prognosis of HCC patients. Besides, knockdowned HSPA4 or HSPA14 can inhibit HCC proliferation, invasion, and metastasis, indicating HSP70s promote HCC occurrence and development. Moreover, HSPA4 and HSPA14 can be novel biomarkers and therapeutic targets for HCC patients.

Availability of data and materials

The transcriptome and clinicopathological data set were obtained from TCGA (https://www.cancer.gov/tcga).

References

  1. 1.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. 2.

    Ingle PV, Samsudin SZ, Chan PQ, Ng MK, Heng LX, Yap SC, Chai AS, Wong AS. Development and novel therapeutics in hepatocellular carcinoma: a review. Ther Clin Risk Manag. 2016;12:445–55.

    CAS  Article  Google Scholar 

  3. 3.

    Rosenzweig R, Nillegoda NB, Mayer MP, Bukau B. The Hsp70 chaperone network. Nat Rev Mol Cell Biol. 2019;20(11):665–80.

    CAS  Article  Google Scholar 

  4. 4.

    Calderwood SK, Gong J. Heat shock proteins promote cancer: it’s a protection racket. Trends Biochem Sci. 2016;41(4):311–23.

    CAS  Article  Google Scholar 

  5. 5.

    Wang Z, Gou W, Liu M, Sang W, Chu H, Zhang W. Expression of P53 and HSP70 in chronic hepatitis, liver cirrhosis, and early and advanced hepatocellular carcinoma tissues and their diagnostic value in hepatocellular carcinoma: an immunohistochemical study. Med Sci Monit. 2015;21:3209–15.

    CAS  Article  Google Scholar 

  6. 6.

    Miao W, Fan M, Huang M, Li JJ, Wang Y. Targeted profiling of heat shock proteome in radioresistant breast cancer cells. Chem Res Toxicol. 2019;32(2):326–32.

    CAS  Article  Google Scholar 

  7. 7.

    Giri B, Sethi V, Modi S, Garg B, Banerjee S, Saluja A, Dudeja V. “Heat shock protein 70 in pancreatic diseases: friend or foe.” J Surg Oncol. 2017;116(1):114–22.

    Article  Google Scholar 

  8. 8.

    Moradi-Marjaneh R, Paseban M, Moradi Marjaneh M. Hsp70 inhibitors: implications for the treatment of colorectal cancer. IUBMB Life. 2019;71(12):1834–45.

    CAS  Article  Google Scholar 

  9. 9.

    Kumar S, Gurshaney S, Adagunodo Y, Gage E, Qadri S, Sharma M, Malik S, Manne U, Singh U, Singh R, et al. Hsp70 and gama-Semino protein as possible prognostic marker of prostate cancer. Front Bioscience (Landmark edition). 2018;23:1987–2000.

    CAS  Article  Google Scholar 

  10. 10.

    Cho W, Jin X, Pang J, Wang Y, Mivechi NF, Moskophidis D. The molecular chaperone heat shock protein 70 controls liver cancer initiation and progression by regulating adaptive DNA damage and mitogen-activated protein kinase/extracellular signal-regulated kinase signaling pathways. Mol Cell Biol. 2019;39(9):e00391.

    CAS  Article  Google Scholar 

  11. 11.

    Xiao H, Wang B, Xiong HX, Guan JF, Wang J, Tan T, Lin K, Zou SB, Hu ZG, Wang K. A novel prognostic index of hepatocellular carcinoma based on immunogenomic landscape analysis. J Cell Physiol. 2021;236(4):2572–91.

    CAS  Article  Google Scholar 

  12. 12.

    Wu FH, Yuan Y, Li D, Liao SJ, Yan B, Wei JJ, Zhou YH, Zhu JH, Zhang GM, Feng ZH. Extracellular HSPA1A promotes the growth of hepatocarcinoma by augmenting tumor cell proliferation and apoptosis-resistance. Cancer Lett. 2012;317(2):157–64.

    CAS  Article  Google Scholar 

  13. 13.

    Lee JH, Han YS, Yoon YM, Yun CW, Yun SP, Kim SM, Kwon HY, Jeong D, Baek MJ, Lee HJ, et al. Role of HSPA1L as a cellular prion protein stabilizer in tumor progression via HIF-1alpha/GP78 axis. Oncogene. 2017;36(47):6555–67.

    CAS  Article  Google Scholar 

  14. 14.

    Jagadish N, Parashar D, Gupta N, Agarwal S, Suri V, Kumar R, Suri V, Sadasukhi TC, Gupta A, Ansari AS, et al. Heat shock protein 70 – 2 (HSP70-2) is a novel therapeutic target for colorectal cancer and is associated with tumor growth. BMC Cancer. 2016;16:561.

    Article  Google Scholar 

  15. 15.

    Zhai LL, Xie Q, Zhou CH, Huang DW, Tang ZG, Ju TF. Overexpressed HSPA2 correlates with tumor angiogenesis and unfavorable prognosis in pancreatic carcinoma. Pancreatology. 2017;17(3):457–63.

    CAS  Article  Google Scholar 

  16. 16.

    Zhang H, Gao H, Liu C, Kong Y, Wang C, Zhang H. Expression and clinical significance of HSPA2 in pancreatic ductal adenocarcinoma. Diagn Pathol. 2015;10:13.

    Article  Google Scholar 

  17. 17.

    Luo C, Xiong H, Chen L, Liu X, Zou S, Guan J, Wang K. GRP78 Promotes Hepatocellular Carcinoma proliferation by increasing FAT10 expression through the NF-kappaB pathway. Exp Cell Res. 2018;365(1):1–11.

    CAS  Article  Google Scholar 

  18. 18.

    Xiong H, Xiao H, Luo C, Chen L, Liu X, Hu Z, Zou S, Guan J, Yang D, Wang K. GRP78 activates the Wnt/HOXB9 pathway to promote invasion and metastasis of hepatocellular carcinoma by chaperoning LRP6. Exp Cell Res. 2019;383(1):111493.

    CAS  Article  Google Scholar 

  19. 19.

    Gu Y, Liu Y, Fu L, Zhai L, Zhu J, Han Y, Jiang Y, Zhang Y, Zhang P, Jiang Z, et al. Tumor-educated B cells selectively promote breast cancer lymph node metastasis by HSPA4-targeting IgG. Nat Med. 2019;25(2):312–22.

    CAS  Article  Google Scholar 

  20. 20.

    Khosla R, Hemati H, Rastogi A, Ramakrishna G, Sarin SK, Trehanpati N. miR-26b-5p helps in EpCAM + cancer stem cells maintenance via HSC71/HSPA8 and augments malignant features in HCC. Liver Int. 2019;39(9):1692–703.

    CAS  Article  Google Scholar 

  21. 21.

    Shan N, Zhou W, Zhang S, Zhang Y. Identification of HSPA8 as a candidate biomarker for endometrial carcinoma by using iTRAQ-based proteomic analysis. Onco Targets Ther. 2016;9:2169–79.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Starenki D, Hong SK, Lloyd RV, Park JI. Mortalin (GRP75/HSPA9) upregulation promotes survival and proliferation of medullary thyroid carcinoma cells. Oncogene. 2015;34(35):4624–34.

    CAS  Article  Google Scholar 

  23. 23.

    Cui X, Li Z, Piao J, Li J, Li L, Lin Z, Jin A. Mortalin expression in pancreatic cancer and its clinical and prognostic significance. Hum Pathol. 2017;64:171–8.

    CAS  Article  Google Scholar 

  24. 24.

    Ma H, Lu T, Zhang X, Li C, Xiong J, Huang L, Liu P, Li Y, Liu L, Ding Z. HSPA12B: a novel facilitator of lung tumor growth. Oncotarget. 2015;6(12):9924–36.

    Article  Google Scholar 

  25. 25.

    Chen W, Liu X, Yuan S, Qiao T. HSPA12B overexpression induces cisplatin resistance in non-small-cell lung cancer by regulating the PI3K/Akt/NF-kappaB signaling pathway. Oncol Lett. 2018;15(3):3883–9.

    PubMed  PubMed Central  Google Scholar 

  26. 26.

    Rhodes D, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan A. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia (New York NY). 2004;6(1):1–6.

    CAS  Article  Google Scholar 

  27. 27.

    Chandrashekar D, Bashel B, Balasubramanya S, Creighton C, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    CAS  Article  Google Scholar 

  28. 28.

    Gao J, Aksoy B, Dogrusoz U, Dresdner G, Gross B, Sumer S, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.

    Article  Google Scholar 

  29. 29.

    Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge. Contemporary oncology (Poznan Poland). 2015;19:A68–77.

    Google Scholar 

  30. 30.

    Yuan R, Wang K, Hu J, Yan C, Li M, Yu X, Liu X, Lei J, Guo W, Wu L, et al. Ubiquitin-like protein FAT10 promotes the invasion and metastasis of hepatocellular carcinoma by modifying beta-catenin degradation. Cancer Res. 2014;74(18):5287–300.

    CAS  Article  Google Scholar 

  31. 31.

    Liu X, Chen L, Ge J, Yan C, Huang Z, Hu J, Wen C, Li M, Huang D, Qiu Y, et al. The ubiquitin-like protein FAT10 stabilizes eEF1A1 expression to promote tumor proliferation in a complex manner. Cancer Res. 2016;76(16):4897–907.

    CAS  Article  Google Scholar 

  32. 32.

    Bruix J, Gores GJ, Mazzaferro V. Hepatocellular carcinoma: clinical frontiers and perspectives. Gut. 2014;63(5):844–55.

    CAS  Article  Google Scholar 

  33. 33.

    Craig AJ, Villanueva A. Liver capsule: Molecular-based signatures in hepatocellular carcinoma. Hepatology. 2016;63(6):2018.

    Article  Google Scholar 

  34. 34.

    Elmallah MIY, Cordonnier M, Vautrot V, Chanteloup G, Garrido C, Gobbo J. Membrane-anchored heat-shock protein 70 (Hsp70) in cancer. Cancer Lett. 2020;469:134–41.

    CAS  Article  Google Scholar 

  35. 35.

    Lianos GD, Alexiou GA, Mangano A, Mangano A, Rausei S, Boni L, Dionigi G, Roukos DH. The role of heat shock proteins in cancer. Cancer Lett. 2015;360(2):114–8.

    CAS  Article  Google Scholar 

  36. 36.

    Garg M, Kanojia D, Saini S, Suri S, Gupta A, Surolia A, Suri A. Germ cell-specific heat shock protein 70 – 2 is expressed in cervical carcinoma and is involved in the growth, migration, and invasion of cervical cells. Cancer. 2010;116(16):3785–96.

    CAS  Article  Google Scholar 

  37. 37.

    Yang Z, Zhuang L, Szatmary P, Wen L, Sun H, Lu Y, Xu Q, Chen X. Upregulation of heat shock proteins (HSPA12A, HSP90B1, HSPA4, HSPA5 and HSPA6) in tumour tissues is associated with poor outcomes from HBV-related early-stage hepatocellular carcinoma. Int J Med Sci. 2015;12(3):256–63.

    CAS  Article  Google Scholar 

  38. 38.

    Morisaki T, Yashiro M, Kakehashi A, Inagaki A, Kinoshita H, Fukuoka T, Kasashima H, Masuda G, Sakurai K, Kubo N, et al. Comparative proteomics analysis of gastric cancer stem cells. PLoS One. 2014;9(11):e110736.

    Article  Google Scholar 

  39. 39.

    Hu YX, Zheng MJ, Zhang WC, Li X, Gou R, Nie X, Liu Q, Hao YY, Liu JJ, Lin B. Systematic profiling of alternative splicing signature reveals prognostic predictor for cervical cancer. J Transl Med. 2019;17(1):379.

    CAS  Article  Google Scholar 

  40. 40.

    Yang MH, Chiang WC, Chou TY, Chang SY, Chen PM, Teng SC, Wu KJ. Increased NBS1 expression is a marker of aggressive head and neck cancer and overexpression of NBS1 contributes to transformation. Clin Cancer Res. 2006;12(2):507–15.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

We are grateful to all of the reviewers for their comments.

Funding

The National Natural Science Foundation of China (82060454) ,the key research and development program of Jiangxi Province of China(20203BBGL73143) supported this study.

Author information

Affiliations

Authors

Contributions

BW, TL, XH, and KW conceived and designed the study. BW, TL, XH, ZHC, LFC, JFG, and KW performed experiments. BW, TL, XH, and KW contributed to Materials and methods. ZHC, WC,RFY, YX, ZGH, and ZD analyzed the data. BW, TL, XH, and KW wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kai Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors are aware of and agree to the content of the paper and their being listed as a co-author of the paper.

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 specific primers used for PCR amplification.

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

Wang, B., Lan, T., Xiao, H. et al. The expression profiles and prognostic values of HSP70s in hepatocellular carcinoma. Cancer Cell Int 21, 286 (2021). https://doi.org/10.1186/s12935-021-01987-9

Download citation

Keywords

  • HSP70 family
  • Hepatocellular carcinoma
  • Bioinformatics
  • HSPA4
  • HSPA14