Abstract
Immunogenic cell death (ICD) is a specialized form of cell death that triggers antitumor immune responses. In tumors, ICD promotes the release of tumor-associated and tumor-specific antigens, thereby reshaping the immune microenvironment, restoring antitumor immunity, and facilitating tumor eradication. However, the regulatory mechanisms of ICD and its immunological effects vary across tumor types, and a comprehensive understanding remains limited. We systematically analyzed the expression of 34 ICD-related regulatory genes across 33 tumor types. Differential expression at the RNA, copy number variation (CNV), and DNA methylation levels was assessed in relation to clinical features. Associations between patient survival and RNA expression, CNVs, single-nucleotide variations (SNVs), and methylation were evaluated. Patients were stratified into immunological subtypes and further divided into high- and low-risk groups based on optimal prognostic models built using a machine learning framework. We explored the relationships between ICD-related genes and immune cell infiltration, stemness, heterogeneity, immune scores, immune checkpoint and regulatory genes, and subtype-specific expression patterns. Moreover, we examined the influence of immunotherapy and anticancer immune responses, applied three machine learning algorithms to identify prognostic biomarkers, and performed drug prediction and molecular docking analyses to nominate therapeutic targets. ICD-related genes were predominantly overexpressed in ESCA, GBM, KIRC, LGG, PAAD, and STAD. RNA expression of most ICD-related genes was associated with poor prognosis, while DNA methylation of these genes showed significant survival correlations in LGG and UVM. Prognostic models were successfully established for 18 cancer types, revealing intrinsic immune regulatory mechanisms of ICD-related genes. Machine learning identified several key prognostic biomarkers across cancers, among which NT5E emerged as a predictive biomarker in head and neck squamous cell carcinoma (HNSC), mediating tumor-immune interactions through multiple ligand-receptor pairs. This study provides a comprehensive view of ICD-related genes across cancers, identifies NT5E as a potential biomarker in HNSC, and highlights novel targets for predicting immunotherapy response and improving clinical outcomes in cancer patients.