Skip to main content
  • Correspondence
  • Open access
  • Published:

Immune-related gene index predicts metastasis for prostate cancer patients undergoing radical radiotherapy

Abstract

In this study, we established a novel immunologic gene prognostic index (IGPI) to predict metastasis and provided new insights into tumor immune microenvironment (TIME) for PCa patients receiving radical radiotherapy. GBP2 and IGF1 were independent factors associated with metastasis-free survival. IGPI score was calculated based on GBP2 and IGF1 and this score was an independent risk factor for PCa patients undergoing radical radiotherapy. Patients with higher IGPI score were at higher risk of metastasis and biochemical recurrence, which were externally validated in the TCGA database and other GEO datasets. IGPI score had demonstrated moderate diagnostic ability of radiation resistance (AUC: 0.889). This score increased with the augment of Gleason score and T stage, as well as biochemical recurrence. Using EPIC, ESTIMATE and immunophenoscore (IPS) algorithms, cancer associated fibroblasts (CAFs), macrophages, stromal score, and estimate score were significantly higher in patients with metastasis group compared to their counterpart. Besides, for CAFs, macrophages, stromal score, and estimate score, patients with higher scores were at higher risk of metastasis, and the HRs were 3.65, 4.01, 4.27, and 3.78, respectively. IGPI score was highly positively associated with stromal score (coefficient: 0.39), immune score (coefficient: 0.43), estimate score (coefficient: 0.45), CAFs (coefficient: 0.42) and macrophages (coefficient: 0.42), while showing the opposite relationship with tumor purity (coefficient: − 0.45). In conclusion, we found that IGPI based on GBP2 and IGF1 might serve as a biomarker predicting metastasis for PCa patients. Besides, the current data further highlight the importance of CAFs in the metastatic process of PCa.

Dear editor,


Among genitourinary cancers, prostate cancer (PCa) has seen the greatest increase in global incident cases over the last 30 years [1]. There has been growing recognition of the vital links between PCa and immunotherapy, including checkpoint inhibitors, cytokines, and therapeutic cancer vaccines [2]. We previously developed an immune-related gene signature to predict progression and immune infiltration for PCa patients undergoing radical prostatectomy [3]. However, radical radiotherapy, as another preferred treatment for intermediate or high-risk localized patients, is rarely studied. About 30% of patients might undergo recurrence as a result of metastasis disease, radioresistant clones and unsatisfactory treatments [4]. In this study, we remedied these problems through establishing a novel immunologic gene prognostic index (IGPI) to predict metastasis and provides new insights into tumor immune microenvironment (TIME) for PCa patients receiving radical radiotherapy. We provided a full-text article in the Additional file 1 for the detailed methods and materials used in this study.

182 differentially expressed genes (Fig. 1A), 4519 tumor-related genes through weighted gene co-expression network analysis (Fig. 1B) and 2660 immune-related genes were used to identify 21 candidate genes (Fig. 1C). Univariate and multivariate COX regression analyses indicated that GBP2 and IGF1 were independent factors associated with metastasis-free survival (Fig. 1D). IGPI score was calculated based on GBP2 and IGF1 and this score was an independent risk factor for PCa patients undergoing radical radiotherapy (Fig. 1E). IGPI score presented certain diagnostic accuracy for metastasis (Fig. 1F). We divided PCa patients according to the median value of IGPI score and patients with higher IGPI score were at higher risk of metastasis (HR: 9.57; 95%CI: 3.23–28.33; Fig. 1G) and biochemical recurrence (BCR) (HR: 3.83; 95%CI: 2.22–6.61; Fig. 1H). In the GSE21034 [5], IGPI score also presented certain value of diagnostic accuracy (AUC: 0.802; 95%CI: 0.702–0.902; Fig. 1I). In addition, patients with higher IGPI score had higher risk of metastasis (HR: 3.85; 95%CI: 2.20–6.74; Fig. 1J) in the GSE134051 [6] and progression (HR: 1.62; 95%CI: 1.08–2.44; Fig. 1K) in the TCGA database. IGPI score had demonstrated moderate diagnostic ability of radiation resistance (AUC: 0.889) in the GSE53902 [7]. This score increased with the augment of Gleason score (Fig. 1M) and T stage (Fig. 1N), as well as BCR (Fig. 1O).

Fig. 1
figure 1

Identification of IGPI and its clinical applications. A volcano plot showing differentially expressed genes between tumor and normal samples; B modules and phenotype showing the tumor-related genes; C Venn diagram showing the intersection of differentially expressed genes, tumor-related genes and immune-related genes; D univariate and multivariate COX regression analysis of candidate genes; E COX regression analysis of clinical and genomic parameters; F time-dependent ROC curve of IGPI score; G analysis of metastasis free survival between high and low risk groups in the GSE116918; H analysis of biochemical recurrence free survival between high and low risk groups in the GSE116918; I ROC curve showing the diagnostic ability of metastasis in the GSE21034; J analysis of metastasis free survival between high and low risk groups in the GSE134051; K analysis of progression free survival between high and low risk groups in the TCGA database; L ROC curve showing the diagnostic ability of radiation resistance in the GSE53902; M comparison between Gleason score and IGPI score; N comparison between T stage and IGPI score; O comparison between biochemical recurrence and IGPI score. IGPI = immunologic gene prognostic index; ROC = receiver operating characteristic. ns, p ≥ 0.05; *p < 0.05; **p < 0.01; ***p < 0.001

Gene set enrichment analysis indicated that high risk group was positively associated with apoptosis, circadian rhythm, cell cycle, T cell receptor signaling pathway, chemokine signaling pathway, mismatch repair, homologous recombination, ECM receptor interaction, and so on, while it was negatively related to arachidonic acid metabolism, fatty acid metabolism, tyrosine metabolism, glutathione metabolism, and drug metabolism cytochrome P450 (Fig. 2A). Furthermore, we established the competing endogenous RNA network and identified the function axis of PART1/has-miR-6885-3p/IGF1 and GBP2 (Fig. 2B). Using EPIC, ESTIMATE and immunophenoscore (IPS) algorithms [8], cancer associated fibroblasts (CAFs), macrophages, stromal score, and estimate score were significantly higher in patients with metastasis group compared to their counterpart (Fig. 2C). Besides, for CAFs, macrophages, stromal score, and estimate score, patients with higher scores were at higher risk of metastasis, and the HRs were 3.65 (1.56–8.54), 4.01 (1.63–9.86), 4.27 (1.84–9.90), and 3.78 (1.58–9.01), respectively (Fig. 2C). However, tumor purity was significantly lower in metastasis group, and patients with higher score were less prone to metastasis compared to those with lower score (HR: 0.26; 95%CI: 0.11–0.63; Fig. 2C). IGPI score was highly positively associated with stromal score (coefficient: 0.39), immune score (coefficient: 0.43), estimate score (coefficient: 0.45), CAFs (coefficient: 0.42) and macrophages (coefficient: 0.42), while showing the opposite relationship with tumor purity (coefficient: − 0.45) (Fig. 2D).

Fig. 2
figure 2

Function enrichment analysis, interaction networks and tumor immune microenvironment analysis. A Gene set enrichment analysis; B competing endogenous RNA network; C TIME analysis; D heatmap presenting the correlations among immunophenoscore, stromal score, immune score, estimate score, tumor purity, cancer associated fibroblasts, macrophages, and IGPI score; E heatmap of clinical information, tumor immune microenvironment parameters, IGPI score, and metastasis-related differentially expressed genes. IGPI = immunologic gene prognostic index

Experimental evidence has implicated GBP2 and IGF1 in the metastasis of other tumors, such as breast cancer [9]. We found that immune-related pathways and ECM receptor interaction might participate in the metastasis of PCa. CAFs, as the most important stromal cells, mediate ECM deposition and remodeling in the reactive stroma and increase tissue stiffness and induce mechanical stress [10]. Prostatic CAFs could induce tumorigenesis in normal human prostatic epithelial cell in vitro via the secretion of CXCL12, and this mechanism was found to be dependent on the presence of TGF-β in vivo in a mouse model, which was consistent with our findings [11, 12].

Conclusion

We found that IGPI based on GBP2 and IGF1 might serve as a biomarker predicting metastasis for PCa patients.

Availability of data and materials

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/additional material.

References

  1. Zi H, He SH, Leng XY, Xu XF, Huang Q, Weng H, et al. Global, regional, and national burden of kidney, bladder, and prostate cancers and their attributable risk factors, 1990–2019. Mil Med Res. 2021;8(1):60.

    Google Scholar 

  2. Wurz GT, Kao CJ, DeGregorio MW. Novel cancer antigens for personalized immunotherapies: latest evidence and clinical potential. Ther Adv Med Oncol. 2016;8(1):4–31.

    Article  CAS  Google Scholar 

  3. Feng D, Zhang F, Li D, Shi X, Xiong Q, Wei Q, et al. Developing an immune-related gene prognostic index associated with progression and providing new insights into the tumor immune microenvironment of prostate cancer. Immunology. 2022;166(2):197–209.

    Article  CAS  Google Scholar 

  4. Jain S, Lyons CA, Walker SM, McQuaid S, Hynes SO, Mitchell DM, et al. Validation of a Metastatic Assay using biopsies to improve risk stratification in patients with prostate cancer treated with radical radiation therapy. Ann Oncol. 2018;29(1):215–22.

    Article  CAS  Google Scholar 

  5. Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, et al. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010;18(1):11–22.

    Article  CAS  Google Scholar 

  6. Friedrich M, Wiedemann K, Reiche K, Puppel SH, Pfeifer G, Zipfel I, et al. The role of lncRNAs TAPIR-1 and -2 as diagnostic markers and potential therapeutic targets in prostate cancer. Cancers (Basel). 2020;12(5):1122.

    Article  CAS  Google Scholar 

  7. Li XR, Zhou KQ, Yin Z, Gao YL, Yang X. Knockdown of FBP1 enhances radiosensitivity in prostate cancer cells by activating autophagy. Neoplasma. 2020;67(5):982–91.

    Article  CAS  Google Scholar 

  8. Zeng D, Ye Z, Shen R, Yu G, Wu J, Xiong Y, et al. IOBR: multi-omics immuno-oncology biological research to decode tumor microenvironment and signatures. Front Immunol. 2021;12: 687975.

    Article  CAS  Google Scholar 

  9. Godoy P, Cadenas C, Hellwig B, Marchan R, Stewart J, Reif R, et al. Interferon-inducible guanylate binding protein (GBP2) is associated with better prognosis in breast cancer and indicates an efficient T cell response. Breast Cancer. 2014;21(4):491–9.

    Article  Google Scholar 

  10. Bonollo F, Thalmann GN, Kruithof-de Julio M, Karkampouna S. The role of cancer-associated fibroblasts in prostate cancer tumorigenesis. Cancers (Basel). 2020;12(7):1887.

    Article  CAS  Google Scholar 

  11. Fiard G, Stavrinides V, Chambers ES, Heavey S, Freeman A, Ball R, et al. Cellular senescence as a possible link between prostate diseases of the ageing male. Nat Rev Urol. 2021;18(10):597–610.

    Article  Google Scholar 

  12. Feng D, Xiong Q, Wei Q, Yang L. Cellular landscape of tumour microenvironment in prostate cancer. Immunology. 2022. https://doi.org/10.1111/imm.13456.

    Article  Google Scholar 

Download references

Acknowledgements

The results showed here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga.

Funding

This program was supported by the National Key Research and Development Program of China (2021YFC2009303), Project of Health Commission of Sichuan Province (21PJ041) and the Key Research and Development Support Plan of Chengdu Science and Technology Bureau (2022-YF05-01568-SN), the National Natural Science Foundation of China (Grant Nos. 81974099, 82170785, 81974098, 82170784), programs from Science and Technology Department of Sichuan Province (Grant Nos. 2021YFH0172), Young Investigator Award of Sichuan University 2017 (Grant No. 2017SCU04A17), Technology Innovation Research and Development Project of Chengdu Science and Technology Bureau (2019-YF05-00296-SN), Sichuan University–Panzhihua science and technology cooperation special fund (2020CDPZH-4). The funders had no role in study design, data collection or analysis, preparation of the manuscript, or the decision to publish.

Author information

Authors and Affiliations

Authors

Contributions

DCF proposed the project, conducted data analysis, interpreted the data, and wrote the manuscript; WZZ, ZHW, and XS conducted data analysis, interpreted the data; LY and PH supervised the project, and interpreted the data; All authors reviewed and edited the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Lu Yang or Ping Han.

Ethics declarations

Ethics approval and consent to participate

The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Consent for publication

Not applicable.

Competing interests

The authors have no conflicts of interest to declare.

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.

The full-text of this correspondence.

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

Feng, D., Zhu, W., Shi, X. et al. Immune-related gene index predicts metastasis for prostate cancer patients undergoing radical radiotherapy. Exp Hematol Oncol 12, 8 (2023). https://doi.org/10.1186/s40164-022-00367-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40164-022-00367-x