Construction and validation of an immune prognostic risk model in oral squamous cell carcinoma

doi:10.19439/j.sjos.2024.04.003

Abstract

Abstract: PURPOSE: To analyze the immune-related core genes differentially expressed in oral squamous cell carcinoma(OSCC) and construct an immune-related prognostic risk model for OSCC patients. METHODS: Weighted gene co-expression network analysis of RNA sequencing data from OSCC patients in the Cancer Genome Atlas (TCGA) database was conducted to identify immune-related modules and core genes. Core genes associated with immune prognosis were screened using univariate Cox regression analysis and survival analysis to construct an immune-related prognostic risk model for OSCC. The prognostic risk model's predictive ability was evaluated using Kaplan-Meier analysis, receiver operating characteristic curves, and external datasets from GSE41613. The expression of 8 immune prognostic core genes in tumor samples from OSCC patients was detected by real-time quantitative PCR assay(RT-qPCR), and the correlation between risk score and depth of invasion was assessed by calculating risk scores for OSCC patients. Statistical analysis was performed with SPSS 21.0 software package. RESULTS: Prognostic risk model for OSCC was successfully constructed based on 8 immune prognostic core genes(CSF2RA, CLEC4C, COL5A3, CTSG, EDNRA, GPC4, GUCY1A2, ANGPT2). The prognostic risk model demonstrated perfect predictive value validated using Kaplan-Meier analysis, receiver operating characteristic curve, and the GSE41613 dataset. The risk scores of OSCC patients calculated based on this model were positively correlated with the depth of invasion, indicating that the model have the ability to predict the potential risk of OSCC. CONCLUSIONS: An OSCC prognostic risk model is constructed based on the signatures of 8 immune prognostic core genes, which may effectively predict the prognosis of OSCC patients, providing an important reference for immune prevention of OSCC.

Key words: Oral squamous cell carcinoma, Weighted gene co-expression network analysis, Immune prognostic genes, Prognostic risk model

CLC Number:

R739.8

ZHAO Jiao, SUI Bai-yan, LIU Xin, RUAN Min. Construction and validation of an immune prognostic risk model in oral squamous cell carcinoma[J]. Shanghai Journal of Stomatology, 2024, 33(4): 345-353.

References

[1] Bray F, Ferlay J, Soerjomataram I, et al.Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA Cancer J Clin, 2018, 68(6): 394-424.
[2] Sung H, Ferlay J, Siegel RL, et al.Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA Cancer J Clin, 2021, 71(3): 209-249.
[3] Gupta B, Johnson NW, Kumar N.Global epidemiology of head and neck cancers: a continuing challenge[J]. Oncology, 2016, 91(1): 13-23.
[4] Siegel RL, Miller KD, Fuchs HE, et al.Cancer statistics, 2021[J]. CA Cancer J Clin, 2021,71(1): 7-33.
[5] Joyce JA, Pollard JW.Microenvironmental regulation of metastasis[J]. Nat Rev Cancer, 2009, 9(4): 239-252.
[6] Geng SK, Fu YP, Fu SM, et al.A tumor microenvironment-related risk model for predicting the prognosis and tumor immunity of breast cancer patients[J]. Front Immunol, 2022, 13: 927565.
[7] Wu J, Li L, Zhang HB, et al.A risk model developed based on tumor microenvironment predicts overall survival and associates with tumor immunity of patients with lung adenocarcinoma[J]. Oncogene, 2021, 40(26): 4413-4424.
[8] Zhang J, Ma C, Qin H, et al.Construction and validation of a metabolic-related genes prognostic model for oral squamous cell carcinoma based on bioinformatics[J]. BMC Med Genomics, 2022, 15(1): 269-284.
[9] Lu W, Yin C, Zhang T, et al.An oxidative stress-related prognostic signature for indicating the immune status of oral squamous cell carcinoma and guiding clinical treatment[J]. Front Genet, 2022, 13: 977902.
[10] Szklarczyk D, Morris JH, Cook H, et al.The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible[J]. Nucleic Acids Res, 2017, 45(D1): D362-d368.
[11] Franz M, Lopes CT, Huck G, et al.Cytoscape.js: a graph theory library for visualisation and analysis[J]. Bioinformatics, 2016, 32(2): 309-311.
[12] Autenshlyus A, Arkhipov S, Mikhailova E, et al.The relationship between cytokine production, CSF2RA, and IL1R2 expression in mammary adenocarcinoma, tumor histopathological parameters, and lymph node metastasis[J]. Technol Cancer Res Treat, 2019, 18: 1533033819883626.
[13] Kunze K, Frank M, Bodner J, et al.Differentiation of primary and metastatic tumours in synchronous multifocal colonic and bronchopulmonary adenocarcinoma by targeted next-generation sequencing[J]. Histopathology, 2014, 64(7): 1041-1043.
[14] Zhao DR, Liu SZ, Sun LL, et al.Glypican-4 gene polymorphism (rs1048369) and susceptibility to Epstein-Barr virus-associated and -negative gastric carcinoma[J]. Virus Res, 2016, 220: 52-56.
[15] Muendlein A, Heinzle C, Brandtner EM, et al.Circulating glypican-4 is a predictor of 24-month overall survival in metastatic breast cancer[J]. Oncol Res Treat, 2023, 46(4): 151-156.
[16] Huang TY, Wu CC, Su WT.Biological and cytoprotective effect of Piper kadsura ohwi against hydrogen-peroxide-induced oxidative stress in human SW1353 cells[J]. Molecules, 2021, 26(20): 6287-6310.
[17] Wu F, Du YY, Hou XJ, et al.A prognostic model for oral squamous cell carcinoma using 7 genes related to tumor mutational burden[J]. Bmc Oral Health, 2022, 22(1): 152-166.
[18] Truby LK, Kwee LC, Agarwal R, et al.Proteomic profiling identifies CLEC4C expression as a novel biomarker of primary graft dysfunction after heart transplantation[J]. J Heart Lung Transplant, 2021, 40(12): 1589-1598.
[19] Chen YS, Teng LC, Liu WH, et al.Identification of biological targets of therapeutic intervention for clear cell renal cell carcinoma based on bioinformatics approach[J]. Cancer Cell Int, 2016, 16: 16-24.
[20] Yang SF, Zou XM, Li JC, et al.Immunoregulation and clinical significance of neutrophils/NETs-ANGPT2 in tumor microenvironment of gastric cancer[J]. Front Immunol, 2022, 13: 1010434.
[21] Li C, Sun CJ, Fan JC, et al.Angiopoietin-2 expression is correlated with angiogenesis and overall survival in oral squamous cell carcinoma[J]. Med Oncol, 2013, 30(2): 571-581.
[22] Lee YJ, Jung ES, Choi J, et al.The EDN1/EDNRA/beta-arrestin axis promotes colorectal cancer progression by regulating STAT3 phosphorylation[J]. Int J Oncol, 2023, 62(1): 13-27.
[23] D'Cruz AK, Dhar H, Vaish R, et al. Depth of invasion in early oral cancers- is it an independent prognostic factor?[J]. Eur J Surg Oncol, 2021, 47(8): 1940-1946.
[24] Wang F, Tan RK, Feng K, et al.Magnetic resonance imaging-based radiomics features associated with depth of invasion predicted lymph node metastasis and prognosis in tongue cancer[J]. J Magn Reson Imaging, 2022, 56(1): 196-209.
[25] Faisal M, Abu Bakar M, Sarwar A, et al.Depth of invasion (DOI) as a predictor of cervical nodal metastasis and local recurrence in early stage squamous cell carcinoma of oral tongue (ESSCOT)[J]. PLoS One, 2018, 13(8): e0202632.