Identification of ISG15 and ZFP36 as novel hypoxia- and immune-related gene signatures contributing to a new perspective for the treatment of prostate cancer by bioinformatics and experimental verification

Abstract

Background: Prostatic cancer (PCa) is one of the most common malignant tumors in men worldwide. Emerging evidence indicates significance of hypoxia and immunity in PCa invasion and metastasis. This study aimed to develop a hypoxia- and immune-related gene risk signature and explore the molecular mechanisms to formulate a better prognostic tool for PCa patients. Methods: The hypoxia and immune scores of all PCa patients in The Cancer Genome Atlas (TCGA) dataset were calculated via the maximally selected rank statistics method and the ESTIMATE algorithm. From common genes identified overlapping hypoxia- and immune-related differentially expressed genes (DE-HRGs and DE-IRGs), a hypoxia- and immune-related gene risk signature was developed utilizing univariate and multivariate Cox regression analyses, and validated in the Memorial Sloan Kettering Cancer Centre (MSKCC) database. The immune cell infiltration level of PCa samples were evaluated with ssGSEA algorithm. Differential expression of prognostic genes was evidenced by immunohistochemistry and western blot (WB) in paired PCa samples. Expression levels of these genes and their variations under regular and hypoxic conditions were examined in cell lines. The functional effects of the prognostic gene on PCa cells were examined by wound healing and transwell assays. Results: A hypoxia- and immune-related gene risk signature constructed by ISG15 and ZFP36 displays significant predictive potency, with higher risk score representing worse survival. A nomogram based on independent prognostic factors including the risk score and Gleason score exhibited excellent clinical value in the survival prediction of PCa. Infiltration levels of eosinophils, neutrophils, Tcm, Tem, TFH, Th1 cells, and Th17 cells were significantly lower in the high-risk group. Conversely, aDC, pDC, T helper cells, and Tregs were significantly higher. Additionally, the two prognostic genes were closely correlated with the tumor-infiltrating immune cell subset in PCa progression. RT-qPCR and WB presented higher and lower expression of ISG15 and ZFP36 in PCa cells, respectively. They were correspondingly increased and decreased in PCa cells under hypoxic conditions. Wound healing and transwell assays showed that over-expression of ISG15 promoted the migration and invasion of PCa cells. Conclusion: Our study identified a novel hypoxia- and immune-related gene signature, contributing a new perspective to the treatment of PCa