Molecular subtype identification and prognosis stratification based on golgi apparatus-related genes in head and neck squamous cell carcinoma

Abstract Background Abnormal dynamics of the Golgi apparatus reshape the tumor microenvironment and immune landscape, playing a crucial role in the prognosis and treatment response of cancer. This study aims to investigate the potential role of Golgi apparatus-related genes (GARGs) in the heterogeneity and prognosis of head and neck squamous cell carcinoma (HNSCC). Methods Transcriptional data and corresponding clinical information of HNSCC were obtained from public databases for differential expression analysis, consensus clustering, survival analysis, immune infiltration analysis, immune therapy response assessment, gene set enrichment analysis, and drug sensitivity analysis. Multiple machine learning algorithms were employed to construct a prognostic model based on GARGs. A nomogram was used to integrate and visualize the multi-gene model with clinical pathological features. Results A total of 321 GARGs that were differentially expressed were identified, out of which 69 were associated with the prognosis of HNSCC. Based on these prognostic genes, two molecular subtypes of HNSCC were identified, which showed significant differences in prognosis. Additionally, a risk signature consisting of 28 GARGs was constructed and demonstrated good performance for assessing the prognosis of HNSCC. This signature divided HNSCC into the high-risk and low-risk groups with significant differences in multiple clinicopathological characteristics, including survival outcome, grade, T stage, chemotherapy. Immune response-related pathways were significantly activated in the high-risk group with better prognosis. There were significant differences in chemotherapy drug sensitivity and immune therapy response between the high-risk and low-risk groups, with the low-risk group being more suitable for receiving immunotherapy. Riskscore, age, grade, and radiotherapy were independent prognostic factors for HNSCC and were used to construct a nomogram, which had good clinical applicability. Conclusions We successfully identified molecular subtypes and prognostic signature of HNSCC that are derived from GARGs, which can be used for the assessment of HNSCC prognosis and treatment responses.