Abstract
BACKGROUND: Heat shock proteins (HSPs) are widely involved in tumor occurrence and development and are prognostic markers for multiple tumors. However, the role of HSPs in clear cell renal cell carcinoma (ccRCC) remains unclear. METHODS: We used Cytoscape to identify hub genes in the ccRCC single-cell sequencing data set from the Gene Expression Omnibus (GEO) data repository. We identified subtypes, C1 and C2, of The Cancer Genome Atlas (TCGA) patients based on the expression of hub genes using unsupervised consensus clustering. Principal component analysis (PCA) was used to verify the clustering differences, and Kaplan-Meier (K-M) estimate was used to verify the survival differences between C1 and C2 patients. We used TIMER 2.0 and CIBERSORT to evaluate the immune cell infiltration of HSP genes and C1 and C2 patients. The R package "pRRophetic" was used to evaluate the sensitivity in C1 and C2 patients to the four first-line treatment drugs. RESULTS: We identified six hub genes (HSP90AA1, HSPH1, HSPA1B, HSPA8, and HSPA1A) encoding HSP, five of which were significantly downregulated in TCGA group, and four had a protective effect on prognosis (p <0.05). Survival analysis showed that C1 patients had a better overall survival (p <0.001). TIMER 2.0 analysis showed that three HSP genes were significantly correlated with the infiltration of CD4+ T cells and CD4+ Th1 cells (|cor|>0.5, p<0.001). CIBERSORT showed significant differences in multiple infiltrating immune cells between C1 and C2 patients. Meanwhile, the expression of PD1 was significantly lower in C1 patients than in C2 patients, and the expression of PDL1 is the another way around. Drug sensitivity analysis showed that C1 patients were more sensitive to sorafenib, pazopanib, and axitinib (p <0.001). CONCLUSION: Our research revealed two molecular subtypes of ccRCC based on 6 HSP genes, and revealed significant differences between the two subtypes in terms of clinical prognosis, immune infiltration, and drug sensitivity.