Skip to main content
  • Primary research
  • Open access
  • Published:

Pseudogene HSPA7 is a poor prognostic biomarker in Kidney Renal Clear Cell Carcinoma (KIRC) and correlated with immune infiltrates

Abstract

Background

Pseudogenes played important roles in tumorigenesis, while there are nearly no reports about the expression and roles of HSPA7 in the cancer.

Methods

Firstly, we used Logistic regression, the KS test, the GEPIA database, UALCAN database and qRT-PCR to analyze the expression level of HSPA7 in KIRC, then we used the Cox regression and the Kaplan–Meier curve to analyze the overall survival (OS) of KIRC patients with different Clinico-pathological parameters. Thirdly, we used the multivariate Cox analysis of influencing factors to compare the correlation between the HSPA7 expression level and the clinical parameters. Finally, we used multi-GSEA analysis and the Tumor Immunoassay Resource (TIMER) database to explore the functional role of HSPA7 in KIRC

Results

The HSPA7 is highly expressed in KIRC tumor tissues, and its expression is related to clinico-pathological features and survival in KIRC patients. GSEA analysis displayed the high expression of HSPA7 in KIRC were related to several tumor-related and immune-related pathways. With the TIMER database analysis we showed that HSPA7 levels were correlated with the CD4+ T cells, neutrophils and Dendritic Cell.

Conclusions

Our study showed that HSPA7 is very important in the tumor progression and may act as a poor prognostic biomarker for KIRC tumor by modulating immune infiltrating cells.

Introduction

The morbidity of renal cell carcinoma (RCC) is about 4.2% of all newly-appeared cancer cases, which make RCC become one of the most frequent malignances worldwide. According to a recent survey, there were about 73,820 new cases of RCC and 14,770 deaths occurred in United States in 2019 [1]. Kidney Renal Clear Cell Carcinoma (KIRC) is the most common kidney cancer subtype [2]. At present, the primary treatment of KIRC is surgery, while 30% of the patients who underwent surgery still experience metastasis [3], and the currently used drugs are not effective and have relatively great side effects. Early identification and diagnosis of KIRC patients can help with more precise clinical treatment. Therefore, it is urgent to discover new and reliable markers to predict the prognosis of patients.

Pseudogenes are non-coding genes lacking of protein-coding ability, and were once labeled as junk genes. However, there is growing evidence indicated that pseudogenes can influence the regulatory mechanisms of many human cancers and pseudogene expression is treated as a novel marker and used in a variety of cancer types to stratify patients subtypes [4, 5] and is therefore taken into account in cancer survival prognostic factors. For example, the pseudogene PRELID1P6 can promote glioma progression through the hnHNPH1-Akt/mTOR pathway [6]. OCT4 abnormally activated pseudogene 5 (OCT4-pg5) can enhance cell proliferation by competing with miR-145 in endometrial carcinoma via upregulating OCT4 expression [7]. High expression of the pseudogene ANXA2P2 has been found to be related to a worse prognosis pseudogene in hepatocellular carcinoma [8]. LDHAP5 was associated with the poor prognosis of ovarian serous cystadenocarcinoma [9]. The Pseudogene HSPA7 (HSP70B) belongs to the HSP70 family (HSPA), discovered in 1985 and encoded near the highly homologous HSPA6 (HSP70B′) on chromosome 1, although mRNA can be expressed after thermal stimulation, it cannot transcribe a functional protein [10]. Numerous investigations have shown that HSPA6 plays an important role in multiple human cancers, including esophageal cancer [11, 12], glioma [13], lung cancer [14], hepatocellular carcinoma [15] and leukemia [16]. However, little has been reported about the expression and role of HSPA7 in cancer. In this study, we reported that high expression of HSPA7 can indicate the poor prognosis of KIRC.

Our study examined the expression and prognostic value of HSPA7 in KIRC patients in the Cancer Genome Atlas (TCGA) and validated them in multiple independent cohorts. Moreover, GSEA [17] and Tumor Immunoassay Resource (TIMER) database [18] were used to assay the potential mechanisms of HSPA7 in KIRC. Our results implied that the functional role of HSPA7 in KIRC may through regulating immune cell infiltration.

Methods

Data mining and data collection

The KIRC data of TCGA consists of 72 normal tissues and 539 tumor samples, was acquired from the TCGA data portal (https://tcga-data.nci.nih.gov/tcga/). Clinical data pertaining to patients’ age, gender, survival, grade, stage, and recurred/progressed outcome were also acquired from the data portal. The dataset including mRNA expression counts and survival data with clinical information. The samples with missing expression data were excluded from the study. The dataset contains survival data with clinical information and mRNA expression counts. The samples with missing expression data were excluded from our study.

Data analysis

The R-3.6.2 project was used to analysis the acquired data. Firstly, we used the Logistic regression and the KS test to analyze the relation between the HSPA7 gene expression and Clinico-pathological features. Then we used the Cox regression and the Kaplan–Meier curve to analyze the overall survival of KIRC patients with different Clinico-pathological parameters from TCGA data. Finally, we used the multivariate Cox analysis of influencing factors to compare the correlation between the HSPA7 expression level and the clinical parameters, such as age, gender grade, stage, T classification, N classification, and M classification, related to survival. The Cutoff Finder.2 was used to determine the cut-off value of HSPA7 expression.

Gene set enrichment analysis (GSEA)

Gene Set Enrichment Analysis (GSEA) is a computational method that determines whether an a priori defined set of genes shows statistically significant between two biological expression states [17]. In our study, an ordered list of genes based on the pathways related to the HSPA7 expression level were generated by the GSEA, and then the significant differences between the high and low-level expression groups of HSPA7 were annotated. The multi-GSEA results and signaling pathway enrichment analysis of phenotypes and were ranked by normalized enrichment score (NES) and the nominal p-value.

Analysis of TIMER database

The TIMER (https://cistrome.shinyapps.io/timer/) database is designed for analysing immune cell infiltrates in multiple cancers. This database can estimate tumour immune infiltration by macrophages, dendritic cells, CD4/CD8+ T cells, neutrophils, and B cells [19]. We used the TIMER database to assess the HSPA7 different expression levels in particular tumours, and then we explored the correlation between HSPA7 expression level and the degree of infiltration in particular immune cell subsets. We further explored the differences in patient survival as a function of gene expression or immune cell infiltration by Kaplan–Meier curve analyses.

Analysis of GEPIA and UALCAN database

The GEPIA (http://gepia.cancer-pku.cn/index.html) database and UALCAN (http://ualcan.path.uab.edu) database can explore the association of mRNA expression level with overall survival (OS). We used these two database to explore the correlation between the HSPA7 expression and patient overall survival in KIRC.

RNA extraction and qRT-PCR analysis

A total of 20 primary KIRC cancer tissues was collected from patients who had undergone surgery at the First Affiliated Hospital of Nanjing Medical University and the Second Affiliated Hospital of Nanjing Medical University. The study was approved by the Ethics Committee of Nanjing Medical University (Nanjing, Jiangsu, PR China), and it was performed in compliance with the Declaration of Helsinki Principles. The clinical information of the 20 KIRC patients was shown in Additional file 1: Table S1). Written informed consent was obtained for all patient samples. RNA extraction and qRT-PCR of the KIRC cancer tissues were performed as the product manual described (Cat# R312-01, Cat# Q131-02, Vazyme, China). The primers used in this study are purchased from Generay (Shanghai, China) and listed as follows.

  • HSPA7-R: CATCCCCAAGGTGCAGAAGT;

  • HSPA7-F: ACCATCCTCTCCACCTCCTC;

  • GAPDH-R: GGGAGCCAAAAGGGTCAT;

  • GAPDH-F: GAGTCCTTCCACGATACCAA.

Results

Characteristics of the of the patients

537 patients’ clinical data were acquired from TCGA, including the age, gender, Histological grade, TNM classification of KIRC (Table 1).

Table 1 Clinical characteristics of TCGA KIRC patients (n = 537)

High HSPA7 mRNA expression in KIRC

First, we assessed the differences in HSPA7 expression between KIRC tumor tissues and adjacent tissues via differential expression scatter plots and paired difference analyses. We find that the expression level of HSPA7 was significantly higher in KIRC tumor tissues (p = 6.183e−35) and in paired cancer tissues (p = 3.311e−18) compared with adjacent tissues (Fig. 1A, B). Then, the expression level of HSPA7 in KIRC tumor tissues and adjacent tissues were verified by GEPIA [20] (Fig. 1C) database, UALCAN database (Fig. 1D) [21] and qRT-PCR analysis (Fig. 1E). The clinical data of 20 patients’ used in qRT-PCR were shown in Additional file 1: Table S1.

Fig. 1
figure 1

HSPA7 is overexpressed in KIRC. A HSPA7 mRNA expression in KIRC based on TCGA Data. B Paired difference analysis of HSPA7 mRNA expression in KIRC based on TCGA Data. C HSPA7 protein expression comparison between normal and tumor tissues obtained from the GEPIA web tool. D HSPA7 protein expression comparison between normal and tumor tissues obtained from the UALCAN web tool. E HSPA7 expression was inspected by qRT-PCR and normalized to GAPDH expression in human KIRC tissues compared with corresponding non-tumor tissues (n = 20) (*p < 0.05, log-rank test)

Correlation between HSPA7 expression level and clinico-pathological features in KIRC tumors

As the Table 2 shown the expression of HSPA7 was highly statistically significantly correlated with clinical stage (p = 0.044) and distant metastasis (positive vs. negative, p = 0.049).

Table 2 Correlation between the clinicopathologic characteristics and HSPA7 mRNA expressiona (logistic regression)

Correlation between KIRC patients survival and HSPA7 expression

To evaluate the effect of HSPA7 expression on KIRC patients survival, the log-rank test and Kaplan–Meier survival analysis were used to estimate the correlation between HSPA7 expression and KIRC patients prognosis. The patients with high HSPA7 expression level displayed relatively poor survival (p = 1.176e−04; Fig. 2A). The clinical subgroup analysis implied that the patients in Histological grade (G1–2 vs. G3–4), clinical stage (Stage I vs. Stage IV), M classification and T classification (the T1–2 vs. T3–4) with HSPA7 expression also had significantly poor overall survival (OS) (Fig. 2B–E), whereas not in the N classification (Fig. 2F). We performed the univariate analysis with the variables and listed in the Table 3. We also performed Multivariate analysis with the Cox proportional hazards model and the results implied that the expression of HSPA7 (HR = 1.304605, p = 0.005187) is a potential prognostic factor for KIRC patients (Table 4). Then we performed the forest plot analysis (Fig. 3), the outcome of KIRC patients are statistically significant correlation with age (p < 0.001), histological grade (p = 0.002), clinical stage (p = 0.019) and the expression of HSPA7 (p < 0.001). In conclusion, HSPA7 is a reliable and effective independent prognostic biomarker of KIRC patients.

Fig. 2
figure 2

Correlation of HSPA7 Expression with Survival and clinicopathologic characteristics. A Survival analysis of HSPA7 expression in terms of overall survival. Kaplan–Meier curves produced survival analysis. B Subgroup analysis of Histological grade (G1–2 vs. G3–4). C Subgroup analysis of clinical stage (stage I and stage IV). D Subgroup analysis of M classification (M0 and M1). E Subgroup analysis of T classification (T1–2 and T3–T4). F Subgroup analysis of N classification (N0/N1)

Table 3 Univariate cox regression of overall survival and clinicopathologic characteristics in TCGA KIRC patients
Table 4 Multivariate analyses of overall survival and clinicopathologic characteristics in TCGA KIRC patients
Fig. 3
figure 3

Forest map analysis of expression and clinicopathologic characteristics

HSPA7-related signaling pathway was identified by GSEA

The differentially regulated pathways between high and low expression of HSPA7 groups were identified by GSEA and then the activated signaling pathways in KIRC were founded. The results with significant differences in enrichment (FDR < 0.25, NOM p < 0.05) in the MSigDB gene set (c2.cp.kegg.v6.2.symbols.gmt) were selected based on the NES and listed in Table 5. Figure 4 showed that renin angiotensin system, primary immunodeficiency, O-glycan biosynthesis, JAK-STAT signaling pathway, hematopoietic cell lineage, intestinal immune network for IgA production, glycerophospholipid metabolism, cytokine-cytokine receptor interaction, cytosolic DNA sensing pathway, autoimmune thyroid disease and asthma. GSEA analysis displayed the high expression of HSPA7 in KIRC were related to several tumor- and immune-related pathways.

Table 5 GSEA identifies an HSPA7-related signaling pathway
Fig. 4
figure 4

Enrichment plots from multiple GSEA

HSPA7 expression correlated with immune cell infiltration in KIRC

Previous studies showed that lymph node metastasis and survival are independently predicted by the frequency of lymphocytes infiltrating in cancer patients. Also GSEA analysis displayed the high expression of HSPA7 in KIRC were related to immune-related pathways. Using TIMER database we investigated whether HSPA7 expression was correlated with six main infiltrating immune cells in KIRC. The result implied that expression of HSPA7 associated with CD4+ T cells (r = 0.395, p-value = 1.24e−18), Macrophage (r = 0.216, p-value = 3.97e−06), neutrophils (r = 0.335, p-value = 1.88e − 13) and Dendritic Cell (DC) (r = 0.212, p-value = 4.86e−06) (Fig. 5). The HSPA7 expression levels was also correlated with tumor purity (cor = 0.125, p-value = 6.98e−03). These results suggested that immune infiltration may serve as a important role in KIRC patient outcomes, and HSPA7 could modulate immune infiltrating cells into KIRC tissues.

Fig. 5
figure 5

HSPA7 expression is correlated with the level of immune infiltration in KIRC. A HSPA7 expression is correlated with the level of immune infiltration in KIRC. B Kaplan–Meier plots of immune infiltration and HSPA7 expression levels in KIRC

Discussion

Our study first reported that pseudogene HSPA7 was expressed highly in KIRC patients and can predict a poor prognosis. We showed that the up-regulated HSPA7 had statistical correlation with histological grade, clinical stage, M classification, T classification and overall survival in KIRC.

HSPA7 belongs to the heat shock protein 70 (HSP70) family, has long been considered as being a pseudogene which is transcribed in response to stress, but now suggest as a high homology to HSPA6 [22]. The HSP70 family is composed of about 13 members, including HSPA1L, HSPA2, HSPA5, HSPA6, HSPA7, HSPA8, HSPA12A, HSPA12B HSPA9, HSPA13 and HSPA14 [23, 24]. Accumulating data indicated that HSP70 family can play a causal role in cancer initiation. Evidence showed that HSPA1L can enhance cancer stem cell-like properties via regulating β-Catenin transcription and activating IGF1Rβ [25]. RNF144A interacted with HSPA2 can promote tumor growth and progression [26]. Down regulation of HSPA5 can promote ANXA1 and repress PSAT1 expression, which inhibiting the osteosarcoma cell proliferation and inducing cell apoptosis [27]. The expression of HSPA6 were found associated with the lung cancer [14], leukemia [16] and baldder cancer’s [28] migration, invasion and proliferation. HSPA8 could regulate the cell viability in pancreatic cancer cells [29] and serve as a molecular target in human hepatocellular carcinoma [30]. Overexpression of HSPA12A can suppresses renal carcinoma cell migration while promotes hepatocellular carcinoma growth [31]. Overexpression of HSPA12B can induce cisplatin resistance in non-small-cell lung cancer (NSCLC) [32]. HSPA9 is associating with survival and proliferation of thyroid carcinoma cells [33, 34]. Less information is available for HSPA7, HSPA13 and HSPA14 representing more distally related members of the HSP70 family. In our research we explored that highly-expressed HSPA7 is related to clinicopathological features of KIRC. Most importantly, univariate and multivariate Cox analyse demonstrated that HSPA7 expression is an independent prognostic indicator of KIRC survival and may be a promising biomarker for clinical applications. Through GSEA analysis, we found that the high expression of HSPA7 in KIRC may related to several immune pathways. HSPA7 expression was found to correlate with the degree of immune infiltration in KIRC through the TIMER database. Knowledge of the immune components has increased over the past decade. Several studies have reported that immune cells from infiltrating tumors are capable of acting as tumor suppressors or promoters in the tumor microenvironment. CD8+ T cells were reported to correlate with the improved survival of cancer patients [35, 36], while regulatory T cells and tumor-associated macrophages were correlated with the promotion of tumor development [37, 38]. Few studies have shown that the HSP70 family members can serve as immunes signature for prognosis of cancers [11]. And the role of Hsp70 in cell immune modulation has remained contentious, only several studies have shown that the HSP70 family members may related to the cell immune. For example, HSPA2 is related to the responses of bone marrow derived dendritic cells to LPS [39], HSPA8 is central at different key steps in the presentation of peptide antigens to CD4+ T cells, with a potential to regulate T and B cell activation and the final secretion of antibodies by plasma cells [40]. HSPA13 is critical for plasma cells development and may be a new target for eliminating pathologic plasma cells [41]. Our research showed that the expression of HSPA7 was significant correlated with macrophage, CD4+ T cells, neutrophils and dendritic cell infiltrating. With the subsequent Kaplan–Meier analysis we found that CD4+ T cells and macrophage cells can predict the KIRC patients prognosis.

Conclusions

In summary, we explored that the pseudogenes HSPA7 is highly expressed in KIRC tumors and is correlated with tumor survival and progression. We implied that the expression level of HSPA7 was moderately positively associated with degree of macrophage, neutrophil, CD4+ T cells and DC infiltration, and weakly positively correlated with the degree of B cells and CD8+ T cells infiltration in KIRC tumor tissues. The pseudogene are believed as therapeutic targets or potential prognostic markers for KIRC tumor patients, while the detailed mechanism of pseudogene affect the KIRC patients prognosis is still to be explored.

Availability of data and materials

The data generated or analyzed during this study are included in this article, or if absent are available from the corresponding author upon reasonable request.

Abbreviations

KIRC:

Kidney Renal Clear Cell Carcinoma

HSPA7:

Heat Shock Protein Family A (Hsp70) Member 7

TIMER:

Tumor Immunoassay Resource

RCC:

Renal cell carcinoma

TCGA:

The Cancer Genome Atlas

GSEA:

Gene Set Enrichment Analysis

NES:

Normalized enrichment score

OS:

Overall survival

DC:

Dendritic cell

References

  1. Siegel RL, Miller KD, Jemal A. CA: a cancer journal for clinicians. Cancer Stat. 2019;69(1):7–34.

    Google Scholar 

  2. Rini BI, Campbell SC, Escudier B. Renal cell carcinoma. Lancet. 2009;373(9669):1119–32.

    Article  CAS  PubMed  Google Scholar 

  3. Motzer RJ, Hutson TE, Cella D, Reeves J, Hawkins R, Guo J, Nathan P, Staehler M, de Souza P, Merchan JR, et al. Pazopanib versus sunitinib in metastatic renal-cell carcinoma. N Engl J Med. 2013;369(8):722–31.

    Article  CAS  PubMed  Google Scholar 

  4. Smerekanych S, Johnson TS, Huang K, Zhang Y. Pseudogene-gene functional networks are prognostic of patient survival in breast cancer. BMC Med Genom. 2020;13(Suppl 5):51.

    Article  CAS  Google Scholar 

  5. Poliseno L, Marranci A, Pandolfi PP. Pseudogenes in human cancer. Front Med (Lausanne). 2015;2:68.

    Google Scholar 

  6. Xi S, Cai H, Lu J, Zhang Y, Yu Y, Chen F, Huang Q, Wang F, Chen Z. The pseudogene PRELID1P6 promotes glioma progression via the hnHNPH1-Akt/mTOR axis. Oncogene. 2021;40(26):4453–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bai M, Yuan M, Liao H, Chen J, Xie B, Yan D, Xi X, Xu X, Zhang Z, Feng Y. OCT4 pseudogene 5 upregulates OCT4 expression to promote proliferation by competing with miR-145 in endometrial carcinoma. Oncol Rep. 2015;33(4):1745–52.

    Article  CAS  PubMed  Google Scholar 

  8. Wang QS, Shi LL, Sun F, Zhang YF, Chen RW, Yang SL, Hu JL. High expression of ANXA2 pseudogene ANXA2P2 promotes an aggressive phenotype in hepatocellular carcinoma. Dis Markers. 2019;2019:9267046.

    PubMed  PubMed Central  Google Scholar 

  9. Lin S, Meng Y, Cao C, Wu P, Gao P, Zhi W, Peng T, Wu P. Comprehensive analysis of LDHAP5 pseudogene expression and potential pathogenesis in ovarian serous cystadenocarcinoma. Cancer Cell Int. 2020;20:229.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Tutter M, Schug C, Schmohl KA, Urnauer S, Schwenk N, Petrini M, Lokerse WJM, Zach C, Ziegler S, Bartenstein P, et al. Effective control of tumor growth through spatial and temporal control of theranostic sodium iodide symporter (NIS) gene expression using a heat-inducible gene promoter in engineered mesenchymal stem cells. Theranostics. 2020;10(10):4490–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Guo X, Wang Y, Zhang H, Qin C, Cheng A, Liu J, Dai X, Wang Z. Identification of the prognostic value of immune-related genes in esophageal cancer. Front Genet. 2020;11:989.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Wang L, Wei Q, Zhang M, Chen L, Li Z, Zhou C, He M, Wei M, Zhao L. Identification of the prognostic value of immune gene signature and infiltrating immune cells for esophageal cancer patients. Int Immunopharmacol. 2020;87:106795.

    Article  CAS  PubMed  Google Scholar 

  13. Sun H, Zou HY, Cai XY, Zhou HF, Li XQ, Xie WJ, Xie WM, Du ZP, Xu LY, Li EM, et al. Network analyses of the differential expression of heat shock proteins in glioma. DNA Cell Biol. 2020;39(7):1228–42.

    Article  CAS  PubMed  Google Scholar 

  14. Wang L, Hou J, Wang J, Zhu Z, Zhang W, Zhang X, Shen H, Wang X. Regulatory roles of HSPA6 in Actinidia chinensis Planch. root extract (acRoots)-inhibited lung cancer proliferation. Clin Transl Med. 2020;10(2):e46.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Tabuchi Y, Ando H, Takasaki I, Feril LB Jr, Zhao QL, Ogawa R, Kudo N, Tachibana K, Kondo T. Identification of genes responsive to low intensity pulsed ultrasound in a human leukemia cell line Molt-4. Cancer Lett. 2007;246(1–2):149–56.

    Article  CAS  PubMed  Google Scholar 

  17. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, Jiang P, Shen H, Aster JC, Rodig S, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. 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(W1):W98–102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Brocchieri L, Conway de Macario E, Macario AJ. hsp70 genes in the human genome: conservation and differentiation patterns predict a wide array of overlapping and specialized functions. BMC Evol Biol. 2008;8:19.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Baquero-Perez B, Whitehouse A. Hsp70 isoforms are essential for the formation of kaposi’s sarcoma-associated herpesvirus replication and transcription compartments. PLoS Pathog. 2015;11(11):e1005274.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Waller SJ, Knighton LE, Crabtree LM, Perkins AL, Reitzel AM, Truman AW. Characterizing functional differences in sea anemone Hsp70 isoforms using budding yeast. Cell Stress Chaperones. 2018;23(5):933–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Choi SI, Lee JH, Kim RK, Jung U, Kahm YJ, Cho EW, Kim IG. HSPA1L enhances cancer stem cell-like properties by activating IGF1Rbeta and regulating beta-catenin transcription. Int J Mol Sci. 2020;21(18):6957.

    Article  CAS  PubMed Central  Google Scholar 

  26. Yang YL, Zhang Y, Li DD, Zhang FL, Liu HY, Liao XH, Xie HY, Lu Q, Zhang L, Hong Q, et al. RNF144A functions as a tumor suppressor in breast cancer through ubiquitin ligase activity-dependent regulation of stability and oncogenic functions of HSPA2. Cell Death Differ. 2020;27(3):1105–18.

    Article  CAS  PubMed  Google Scholar 

  27. Tang X, Luo L, Li Y, Wu H, Hu Q, Yue H, He X, Zou J, Min S. Therapeutic potential of targeting HSPA5 through dual regulation of two candidate prognostic biomarkers ANXA1 and PSAT1 in osteosarcoma. Aging. 2020;13:1212.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Shin SS, Song JH, Hwang B, Noh DH, Park SL, Kim WT, Park SS, Kim WJ, Moon SK. HSPA6 augments garlic extract-induced inhibition of proliferation, migration, and invasion of bladder cancer EJ cells; implication for cell cycle dysregulation, signaling pathway alteration, and transcription factor-associated MMP-9 regulation. PLoS ONE. 2017;12(2):e0171860.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Tian Y, Xu H, Farooq AA, Nie B, Chen X, Su S, Yuan R, Qiao G, Li C, Li X, et al. Maslinic acid induces autophagy by down-regulating HSPA8 in pancreatic cancer cells. Phytother Res. 2018;32(7):1320–31.

    Article  CAS  PubMed  Google Scholar 

  30. 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.

    Article  CAS  PubMed  Google Scholar 

  31. Min X, Zhang X, Li Y, Cao X, Cheng H, Li Y, Li C, Kong Q, Mao Q, Peng P, et al. HSPA12A unstabilizes CD147 to inhibit lactate export and migration in human renal cell carcinoma. Theranostics. 2020;10(19):8573–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. 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 

  33. Starenki D, Sosonkina N, Hong SK, Lloyd RV, Park JI. Mortalin (GRP75/HSPA9) promotes survival and proliferation of thyroid carcinoma cells. Int J Mol Sci. 2019;20(9):2069.

    Article  CAS  PubMed Central  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. 2013;14(10):1014–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Governa V, Trella E, Mele V, Tornillo L, Amicarella F, Cremonesi E, Muraro MG, Xu H, Droeser R, Daster SR, et al. The Interplay between neutrophils and CD8(+) t cells improves survival in human colorectal cancer. Clin Cancer Res. 2017;23(14):3847–58.

    Article  CAS  PubMed  Google Scholar 

  37. Nishikawa H, Sakaguchi S. Regulatory T cells in cancer immunotherapy. Curr Opin Immunol. 2014;27:1–7.

    Article  CAS  PubMed  Google Scholar 

  38. Noy R, Pollard JW. Tumor-associated macrophages: from mechanisms to therapy. Immunity. 2014;41(1):49–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Van Goor A, Slawinska A, Schmidt CJ, Lamont SJ. Distinct functional responses to stressors of bone marrow derived dendritic cells from diverse inbred chicken lines. Dev Comp Immunol. 2016;63:96–110.

    Article  PubMed  CAS  Google Scholar 

  40. Bonam SR, Ruff M, Muller S. HSPA8/HSC70 in immune disorders: a molecular rheostat that adjusts chaperone-mediated autophagy substrates. Cells. 2019;8(8):849.

    Article  CAS  PubMed Central  Google Scholar 

  41. He Y, Xu R, Zhai B, Fang Y, Hou C, Xing C, Xiao H, Chen G, Wang X, Ma N, et al. Hspa13 promotes plasma cell production and antibody secretion. Front Immunol. 2020;11:913.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

CD and RH collected and analyzed the data and were the major contributors to writing the original draft of the manuscript. ZD and JW designed the experiments and analyzed the data. ZD and JW reviewed and edited the manuscript. JZ coordinated the study. All authors have agreed to the published version of the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zhan Dong or Jun Wu.

Ethics declarations

Ethics approval and consent to participate

A total of 20 primary KIRC cancer tissues was collected from patients who had undergone surgery at the First Affiliated Hospital of Nanjing Medical University and the Second Affiliated Hospital of Nanjing Medical University. The study was approved by the Ethics Committee of Nanjing Medical University (Nanjing, Jiangsu, PR China), the approval number from the Ethical Committee is 2016-SRFA-011, and it was performed in compliance with the Declaration of Helsinki Principles. Written informed consent was obtained for all patient samples.

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.

Clinical characteristics of KIRC patients (n = 20).

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

Ding, C., He, R., Zhang, J. et al. Pseudogene HSPA7 is a poor prognostic biomarker in Kidney Renal Clear Cell Carcinoma (KIRC) and correlated with immune infiltrates. Cancer Cell Int 21, 435 (2021). https://doi.org/10.1186/s12935-021-02141-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-021-02141-1

Keywords