A Cytokine-Related Gene Signature for Pan-Cancer Prognostic Stratification and Malignant Phenotype Characterization

1. Introduction

Cytokines are central mediators of intercellular communication within the tumor microenvironment and play critical roles in regulating inflammation, immune responses, tissue remodeling, and cancer progression [1,2,3]. Dysregulated cytokine signaling has been widely implicated in tumor initiation, malignant transformation, immune evasion, and therapeutic resistance across diverse cancer types [1,4,5,6]. Rather than acting as isolated factors, cytokines function through complex and dynamic networks that integrate signals from tumor cells, immune cells, stromal components, and the extracellular matrix [7]. This intricate cytokine landscape contributes to substantial intertumoral and intratumoral heterogeneity, posing major challenges for accurate prognosis assessment and precision oncology [8,9].

Accumulating evidence suggests that cytokine-driven transcriptional programs are closely associated with aggressive tumor phenotypes, including enhanced proliferation, epithelial-mesenchymal transition, metastatic potential, and remodeling of the immune microenvironment [10,11,12]. Specific cytokines such as IL6, IL10, and IL1 family members have been linked to poor clinical outcomes in multiple malignancies, while others may exert context-dependent or even protective effects [13,14]. However, most prior studies have focused on individual cytokines or limited signaling axes, which may fail to capture the coordinated and system-level impact of cytokine networks on cancer biology. A comprehensive evaluation of cytokine-related gene expression patterns may therefore provide more robust and biologically meaningful prognostic information than single-gene biomarkers [15,16,17].

With the increasing availability of large-scale transcriptomic datasets, integrative computational approaches enable systematic dissection of cytokine-related signatures across cancer types. Regularized regression models, such as Least Absolute Shrinkage and Selection Operator Cox regression, allow the construction of parsimonious prognostic signatures while minimizing overfitting and multicollinearity [18,19,20]. When combined with pan-cancer analyses and independent validation cohorts, these strategies facilitate the identification of conserved cytokine-associated risk patterns and their relevance to patient survival. Moreover, linking cytokine-related signatures to hallmark malignant processes, including cell cycle progression, EMT, angiogenesis, and metabolic reprogramming, may provide mechanistic insights into how inflammatory signaling shapes tumor behavior [21,22,23].

In this study, we systematically investigated cytokine-related genes across pan-cancer cohorts to develop and validate a robust prognostic risk score. By integrating transcriptomic data with survival outcomes, we constructed a cytokine-related gene signature that stratifies patients into distinct risk groups with significantly different clinical prognoses. We further evaluated the distribution and prognostic performance of this signature across cancer types, assessed its clinical utility through nomogram modeling, and explored its biological relevance by correlating the risk score with malignant phenotypes and functional enrichment patterns. Together, our findings provide a comprehensive framework for understanding the prognostic and biological significance of cytokine dysregulation in cancer and highlight its potential value for pan-cancer risk stratification and personalized clinical management.

2. Results

2.1. Identification and Construction of a Cytokine-Related Prognostic Signature Using LASSO Cox Regression

To systematically characterize the prognostic relevance of cytokine-related genes across cancers, we applied LASSO Cox proportional hazards regression to the training cohort to construct a stable and interpretable prognostic signature. As illustrated in Figure 1A, coefficient trajectories demonstrated progressive shrinkage with increasing penalization, indicating effective control of multicollinearity among cytokine-associated candidates.

Using the one-standard-error criterion, an optimal lambda value of 9e-04 was selected, resulting in a parsimonious model comprising 16 cytokine-related genes with non-zero coefficients, including PANX1, FRMD8, GAPDH, MMP8, IL10, IL1A, NOD2, ZC3H12A, CLEC5A, S100A12, IL6, FERMT1, FN1, CX3CL1, CLEC9A, and ARG2. These genes collectively represent multiple facets of cytokine-driven tumor biology, encompassing pro-inflammatory signaling, innate immune sensing, leukocyte recruitment, extracellular matrix remodeling, and metabolic adaptation. Subsequent univariate Cox regression analysis revealed marked heterogeneity in prognostic effects among individual genes (Figure 1B), with PANX1, FRMD8, GAPDH, and MMP8 showing hazard ratios greater than 1, consistent with their reported roles in sustaining inflammatory microenvironments and facilitating tumor invasion, whereas ARG2, CLEC9A, and CX3CL1 exhibited protective associations, suggesting potential involvement in immune surveillance or anti-tumor immune modulation. Importantly, the coexistence of both risk-associated and protective cytokine-related genes within the same signature underscores the complex, context-dependent nature of cytokine signaling in cancer progression. To further evaluate redundancy and coordinated regulation within the model, pairwise Pearson correlation analysis was performed.

As shown in Figure 1C, moderate positive correlations were observed among subsets of inflammatory mediators such as IL6, IL10, CLEC5A, and NOD2, reflecting shared activation of innate immune and cytokine signaling pathways, while overall correlation coefficients remained below levels indicative of problematic collinearity, supporting the robustness and stability of the model. Collectively, these results demonstrate that the 16-gene cytokine-related signature integrates diverse yet complementary inflammatory and immune-regulatory signals, providing a biologically coherent framework for downstream prognostic stratification and functional interpretation in pan-cancer analyses.

2.2. Pan-Cancer Distribution and Prognostic Relevance of the Cytokine-Related Risk Score

Following construction of the 16-gene cytokine-related prognostic signature, we systematically investigated the distribution of the cytokine-related risk score across cancer types and its prognostic implications at the pan-cancer level. As shown in Figure 2A, the risk score exhibited pronounced heterogeneity among different tumor entities. Cancers characterized by chronic inflammation and immune activation, including head and neck squamous cell carcinoma, lung squamous cell carcinoma, bladder cancer, and multiple gastrointestinal malignancies, generally displayed higher risk score distributions, whereas endocrine-related and relatively indolent tumors, such as thyroid carcinoma and prostate adenocarcinoma, showed lower scores. This pattern suggests that cytokine-driven transcriptional programs are strongly shaped by tissue context and tumor-specific immune microenvironments. To further evaluate the prognostic impact of the cytokine-related risk score, univariate Cox regression analyses were conducted across cancer types for progression-free interval (PFI), disease-specific survival (DSS), and overall survival (OS). As illustrated by the bubble plots in Figure 2A, elevated risk scores were significantly associated with adverse survival outcomes in a wide spectrum of malignancies, with bubble size reflecting statistical significance and color intensity indicating effect magnitude. In the TCGA training cohort, forest plot analysis demonstrated that the cytokine-related risk score functioned predominantly as a hazard factor for survival in most cancer types (Figure 2B), with hazard ratios exceeding 1 for OS in multiple tumors, supporting its unfavorable prognostic role. Importantly, these findings were independently validated in the TCGA test cohort, where the majority of cancer types again showed consistent hazard effects associated with higher risk scores (Figure 2C), confirming the robustness and generalizability of the prognostic signal. To assess the clinical discriminative ability of the signature at the cohort level, patients in the TCGA training cohort were stratified into high-risk and low-risk groups according to the median cytokine-related risk score. Kaplan-Meier survival analyses revealed that high-risk patients experienced significantly worse DSS, OS, and PFI compared with low-risk patients (p < 0.001, Figure 2D).

This prognostic stratification was reproducibly observed in the TCGA test cohort, where high-risk patients consistently exhibited inferior survival outcomes across all endpoints (Figure 2E). Collectively, these results demonstrate that the cytokine-related risk score captures clinically relevant inflammatory and immune dysregulation signals and provides stable prognostic stratification across cancer types and independent cohorts, underscoring its potential value as a pan-cancer prognostic biomarker.

2.3. Cancer-Specific Prognostic Performance of the Cytokine-Related Risk Score in Representative Tumor Types

To further delineate the cancer-type-specific prognostic value of the cytokine-related risk score, Kaplan-Meier survival analyses were performed in multiple representative malignancies using disease-specific survival (DSS), overall survival (OS), and progression-free interval (PFI) as clinical endpoints. As shown in Figure 3A, patients with bladder urothelial carcinoma (BLCA) stratified into the high-risk group exhibited significantly poorer DSS and OS compared with those in the low-risk group, whereas the difference in PFI did not reach statistical significance. These findings suggest that cytokine-driven transcriptional programs in BLCA are more closely linked to mortality risk rather than early disease progression. In pancreatic adenocarcinoma (PAAD), a highly aggressive and inflammation-associated malignancy, the cytokine-related risk score demonstrated strong and consistent prognostic power across all three endpoints (Figure 3B). High-risk patients showed markedly reduced DSS, OS, and PFI (p < 0.001), highlighting the critical role of dysregulated cytokine signaling in driving rapid disease progression and poor outcomes in PAAD. In lung adenocarcinoma (LUAD), elevated risk scores were also associated with significantly worse DSS, OS, and PFI (Figure 3C), indicating that cytokine-related inflammatory and immune-regulatory pathways contribute to both survival and disease progression in this tumor type. By contrast, in cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), the prognostic impact of the cytokine-related risk score appeared endpoint-dependent (Figure 3D).

While high-risk patients experienced significantly inferior DSS and OS, no significant difference was observed for PFI, suggesting that cytokine-related dysregulation in CESC may preferentially influence long-term survival rather than short-term progression dynamics. Collectively, these cancer-specific analyses demonstrate that although the cytokine-related risk score consistently stratifies patient survival in multiple tumor types, its association with disease progression varies across cancers, reflecting underlying biological heterogeneity and context-dependent roles of cytokine signaling within distinct tumor microenvironments.

2.4. Development and Validation of a Cytokine-Related Nomogram for Individualized Survival Prediction

To enhance the clinical applicability of the cytokine-related risk score and provide an intuitive tool for individualized survival prediction, we integrated the risk score with key clinical variables to construct a comprehensive prognostic nomogram. As shown in Figure 4A, the nomogram incorporates the cytokine-related risk score, patient age, and cancer type, with each variable assigned a weighted point value proportional to its contribution to overall survival risk. By summing the total points for an individual patient, the nomogram enables direct estimation of 1-, 3-, 5-, and 10-year overall survival probabilities, thereby facilitating personalized risk assessment in routine clinical settings. Calibration analysis demonstrated a high degree of concordance between nomogram-predicted and observed survival outcomes. As illustrated in Figure 4B, the calibration curve for overall survival showed close alignment with the 45-degree reference line, indicating that the nomogram provided accurate survival probability estimates with minimal systematic bias. To further quantify predictive performance, time-dependent receiver operating characteristic analyses were performed. The ROC curves at 5 years (Figure 4C) revealed that the nomogram consistently outperformed the cytokine-related risk score alone in both the TCGA training and test cohorts, achieving higher area under the curve values, thereby confirming the added prognostic value of integrating clinical factors with molecular information. This advantage was maintained across time, as demonstrated by the dynamic AUC profiles in Figure 4D, where the nomogram exhibited superior and more stable discrimination ability over long-term follow-up compared with the risk score alone. Beyond statistical performance, the potential clinical utility of the nomogram was evaluated using decision curve analysis.

As shown in Figure 4E and Figure 4F, the nomogram yielded a higher standardized net benefit than the cytokine-related risk score across a wide range of threshold probabilities in both the training and test cohorts. Importantly, the nomogram consistently outperformed default strategies of treating all or no patients, indicating that its application could improve clinical decision-making by optimizing the balance between benefit and harm. Collectively, these results demonstrate that the cytokine-related nomogram provides accurate, robust, and clinically meaningful survival prediction, supporting its potential role as a practical decision-support tool for risk stratification and personalized management across diverse cancer types.

2.5. Association Between the Cytokine-Related Risk Score and Malignant Biological Features Across Pan-Cancer

To explore the biological implications underlying the cytokine-related risk score, we systematically examined its relationship with key malignant features, including angiogenesis, epithelial-mesenchymal transition (EMT), and cell cycle activity, using pathway-level scores across the pan-cancer cohort. As shown in Figure 5A-C, at the global pan-cancer level, the cytokine-related risk score exhibited a significant positive correlation with EMT score and cell cycle score, indicating that tumors with higher cytokine-related risk scores tend to display enhanced mesenchymal characteristics and proliferative capacity. In contrast, a weak but statistically significant negative correlation was observed between the risk score and angiogenesis score, suggesting that cytokine-driven risk may not uniformly reflect angiogenic activity and may instead capture alternative inflammatory or immune-mediated malignant programs.

To further delineate tumor-specific patterns, stratified correlation analyses were performed across individual cancer types. As illustrated in Figure 5D, the relationship between the cytokine-related risk score and angiogenesis score varied markedly across cancers. Positive correlations were observed in selected tumor types such as COAD, DLBC, and READ, whereas weak or negative correlations predominated in other cancers, highlighting substantial inter-tumoral heterogeneity in cytokine-angiogenesis coupling. In contrast, the association between the risk score and EMT score was consistently positive across the majority of cancer types (Figure 5E), including BLCA, BRCA, COAD, KIRC, LGG, LUAD, and STAD, supporting the notion that cytokine-related transcriptional programs are broadly linked to EMT activation and invasive phenotypes. Similarly, correlation analyses between the risk score and cell cycle score revealed widespread positive associations across multiple malignancies (Figure 5F), particularly in highly proliferative tumors such as LGG, LUAD, PRAD, and THCA, indicating that elevated cytokine-related risk is closely aligned with enhanced cell cycle progression and tumor growth dynamics. Collectively, these findings demonstrate that the cytokine-related risk score is not merely a prognostic indicator but also reflects fundamental malignant biological behaviors, particularly EMT activation and proliferative signaling, while its relationship with angiogenesis appears to be context-dependent and tumor-specific. These results provide mechanistic insights into how cytokine dysregulation may drive aggressive tumor phenotypes across diverse cancer types.

2.6. Functional Characterization of Cytokine-Related Risk Groups Reveals Distinct Biological Programs

To further elucidate the biological mechanisms underlying the prognostic stratification defined by the cytokine-related risk score, we performed differential expression analysis between high-risk and low-risk groups followed by comprehensive functional enrichment analyses. Gene Ontology and KEGG pathway enrichment results demonstrated that genes upregulated in the high-risk group were strongly associated with fundamental malignant processes (Figure 6A). At the biological process level, high-risk tumors showed significant enrichment of cell cycle-related terms, including nuclear division, mitotic nuclear division, chromosome segregation, and sister chromatid segregation, indicating enhanced proliferative activity. Consistently, cellular component analysis highlighted enrichment in chromosomal regions and condensed chromosomes, further supporting active mitosis and genomic instability in cytokine-high tumors. Molecular function terms were predominantly related to extracellular matrix structural constituents, integrin binding, proteoglycan binding, and microtubule binding, suggesting coordinated activation of extracellular matrix remodeling and cytoskeletal reorganization. Notably, KEGG pathway analysis revealed enrichment of ECM-receptor interaction and IL-17 signaling pathway, implicating cytokine-driven inflammatory signaling in promoting invasive and proliferative tumor phenotypes. In contrast, genes upregulated in the low-risk group exhibited a markedly different functional profile (Figure 6B). Biological processes were enriched for transmembrane transport-related activities, including anion and carboxylic acid transport, vascular transport, and transport across the blood-brain barrier, reflecting a more differentiated and metabolically regulated cellular state. Cellular component terms were dominated by plasma membrane-associated structures, such as apical, basal, and basolateral plasma membranes, as well as neuronal cell bodies, suggesting preserved cellular polarity and tissue-specific functions. Molecular function enrichment further emphasized active and organic anion transmembrane transporter activity. KEGG pathways enriched in the low-risk group included adrenergic signaling in cardiomyocytes, thyroid hormone synthesis, bile secretion, and one-carbon pool by folate, indicating engagement of endocrine, metabolic, and homeostatic signaling rather than aggressive malignant programs.

Collectively, these results demonstrate that high cytokine-related risk tumors are characterized by hyperproliferation, extracellular matrix remodeling, and inflammatory signaling, whereas low-risk tumors preferentially maintain transport, metabolic, and differentiated cellular functions. This functional divergence provides mechanistic insight into how cytokine dysregulation contributes to tumor aggressiveness and poor clinical outcomes.

2.7. Immunohistochemical Validation of PANX1 and FRMD8 Expression in Normal and Malignant Tissues

To further support the transcriptomic findings and validate the protein-level expression of key genes within the cytokine-related prognostic signature, immunohistochemical staining data for representative signature genes were examined in normal and tumor tissues. As shown in Figure 7A-D, PANX1 protein expression was generally absent or undetectable in corresponding normal tissues, including breast, pancreas, lung, and cervix, whereas moderate cytoplasmic staining was consistently observed in the matched cancer tissues. This differential expression pattern suggests that PANX1 is upregulated during malignant transformation across multiple tissue types and supports its role as a risk-associated gene within the cytokine-related signature. Consistent with these observations, FRMD8 exhibited low or undetectable expression in normal breast and lung tissues, while markedly increased staining intensity was observed in breast cancer, pancreatic cancer, lung cancer, and endometrial cancer specimens (Figure 7E-H). Notably, FRMD8 showed high and diffuse cytoplasmic expression in several tumor types, indicating broad activation at the protein level in malignant tissues. Together, these immunohistochemical results corroborate the transcriptomic analyses and demonstrate that PANX1 and FRMD8 are upregulated at the protein level in diverse cancers, providing additional biological and translational support for their inclusion in the cytokine-related prognostic signature.

3. Discussion

In the present study, we systematically characterized cytokine-related transcriptional programs across pan-cancer cohorts and developed a robust prognostic signature based on cytokine-associated genes. By integrating large-scale transcriptomic data with survival outcomes, we identified a 16-gene cytokine-related risk score derived from an initial pool of 131 cytokine-related genes curated from GSEA and MSigDB. This signature consistently stratified patients into distinct risk groups with significantly different clinical outcomes across multiple cancer types and independent cohorts, underscoring the broad prognostic relevance of cytokine dysregulation in cancer. Unlike single-cytokine biomarkers, the composite nature of this signature captures the coordinated and context-dependent effects of inflammatory and immune signaling networks within the tumor microenvironment [24].

The cytokine-related risk score demonstrated pronounced heterogeneity across tumor types, reflecting tissue-specific immune landscapes and inflammatory states. Tumors characterized by chronic inflammation or strong immune cell infiltration tended to exhibit higher risk scores and poorer prognoses, whereas hormonally regulated or relatively indolent cancers displayed lower scores and attenuated prognostic effects. Importantly, the risk score acted predominantly as a hazard factor in both training and test cohorts, indicating robust generalizability [25]. Cancer-specific analyses further revealed that while the signature consistently predicted survival, its association with disease progression varied across malignancies, suggesting that cytokine-driven programs may differentially influence tumor growth, invasion, and metastatic dissemination depending on biological context [26].

Beyond its prognostic value, the cytokine-related risk score was closely linked to fundamental malignant phenotypes [27]. Correlation analyses revealed strong associations between elevated risk scores and enhanced epithelial-mesenchymal transition and cell cycle activity across most cancer types, supporting a role for cytokine dysregulation in promoting invasive and proliferative tumor behavior [28]. Functional enrichment analyses further demonstrated that high-risk tumors were characterized by activation of cell cycle progression, chromosome segregation, extracellular matrix remodeling, and inflammatory signaling pathways, including IL-17-related processes [29]. In contrast, low-risk tumors preferentially enriched pathways related to transmembrane transport, metabolic regulation, and tissue-specific homeostatic functions [1]. Together, these findings provide mechanistic insight into how cytokine-related transcriptional programs may drive tumor aggressiveness and adverse clinical outcomes [3].

To facilitate clinical translation, we integrated the cytokine-related risk score with clinical variables to construct a prognostic nomogram, which showed improved predictive accuracy and clinical net benefit compared with the molecular score alone [30]. This integrated model enables individualized survival estimation and highlights the potential utility of cytokine-based signatures as decision-support tools in precision oncology [31]. Nevertheless, several limitations should be acknowledged. First, the study is based on retrospective transcriptomic data, and prospective validation in independent clinical cohorts is warranted [32]. Second, although the signature reflects coordinated cytokine activity, it does not directly capture post-transcriptional regulation or protein-level cytokine dynamics within the tumor microenvironment [33]. Finally, functional experiments are required to delineate the causal roles of individual genes within the signature [34]. Despite these limitations, our findings provide a comprehensive framework linking cytokine dysregulation to pan-cancer prognosis and tumor biology, and they support further exploration of cytokine-related signatures as biomarkers and potential therapeutic targets in cancer.

4. Materials and Methods

4.1. Data Collection and Preprocessing

Transcriptomic expression profiles and corresponding clinical annotations were obtained from The Cancer Genome Atlas pan-cancer cohort through the UCSC Xena platform. Only primary tumor samples with complete survival information were included. Overall survival, disease-specific survival, and progression-free interval were defined according to TCGA standardized clinical endpoints [35]. Gene expression matrices were log2-transformed and subsequently standardized by z-score normalization within each dataset to reduce technical variability. Samples were randomly divided into a training cohort and an internal test cohort at a ratio of 7:3 for model development and validation [36].

4.2. Cytokine-Related Gene Curation Based on GSEA and MSigDB

To comprehensively define cytokine-related genes, we systematically curated gene sets associated with cytokine signaling and cytokine-mediated immune responses from Gene Set Enrichment Analysis resources and the Molecular Signatures Database (MSigDB; URL: https://www.gsea-msigdb.org/gsea/msigdb/human/collections.jsp). Specifically, cytokine-related gene sets were retrieved from the MSigDB collections, including Hallmark, C2 curated gene sets, and Gene Ontology biological process categories related to cytokine activity, cytokine-mediated signaling pathways, inflammatory response, interleukin signaling, and immune regulation [37]. After removing redundant entries and consolidating overlapping genes across gene sets, a total of 131 cytokine-related genes were identified and used as the initial candidate pool for downstream prognostic modeling [38,39].

4.3. Survival Analysis and Pan-Cancer Validation

Patients were stratified into high-risk and low-risk groups according to the median cytokine-related risk score within each cohort. Kaplan-Meier survival analyses were performed to compare survival outcomes between groups, and statistical significance was assessed using the log-rank test. Univariate Cox regression analyses were conducted across cancer types to evaluate the prognostic impact of the risk score on overall survival, disease-specific survival, and progression-free interval [24]. The robustness and generalizability of the signature were independently validated in the TCGA test cohort [40].

4.4. Nomogram Construction and Performance Evaluation

To improve clinical applicability, a prognostic nomogram was developed by integrating the cytokine-related risk score with key clinical variables, including age and cancer type, using multivariable Cox regression analysis. Calibration curves were generated to assess agreement between predicted and observed survival probabilities [30]. Time-dependent receiver operating characteristic analyses were used to evaluate predictive accuracy over time, and dynamic area under the curve values were compared between the risk score alone and the nomogram model. Decision curve analysis was performed to quantify the clinical net benefit of each model across a range of threshold probabilities [41,42].

4.5. Functional Enrichment and Pathway Analysis

Differentially expressed genes between high-risk and low-risk groups were identified using predefined statistical thresholds. Gene Ontology enrichment analyses covering biological process, cellular component, and molecular function categories, as well as Kyoto Encyclopedia of Genes and Genomes pathway analyses, were performed using the clusterProfiler package in R [43]. In addition, pathway-level scores for angiogenesis, epithelial-mesenchymal transition, and cell cycle activity were calculated, and Pearson correlation analyses were conducted to assess associations between these malignant features and the cytokine-related risk score at both pan-cancer and cancer-specific levels [44].

4.6. Immunohistochemistry (IHC) Analysis

To validate the protein-level expression of representative cytokine-related signature genes, immunohistochemistry analyses were performed using publicly available tissue microarray-based staining data. Representative IHC images for PANX1 and FRMD8 in normal and corresponding tumor tissues were obtained from a curated human protein expression database [45]. Multiple tissue types, including breast, pancreas, lung, cervix, and endometrium, were examined to assess expression patterns across different malignancies. IHC staining intensity was annotated according to the database-defined scoring system and categorized as not available (NA), low, medium, or high based on staining strength and the proportion of positively stained cells. Protein expression levels in tumor tissues were qualitatively compared with those in matched normal tissues to evaluate differential expression associated with malignant transformation. All IHC images were independently reviewed to ensure consistency of staining patterns across tissue types [46]. These analyses were used to corroborate transcriptomic findings and to provide translational evidence supporting the involvement of PANX1 and FRMD8 in cytokine-related tumor biology.

4.7. Immune Infiltration Analysis

All statistical analyses were performed using R (4.4.1) software. Two-sided p values less than 0.05 were considered statistically significant unless otherwise specified. Continuous variables were summarized using appropriate descriptive statistics, and correlations were evaluated using Pearson correlation coefficients. All analyses followed reproducible research principles.

5. Conclusions

In summary, this study systematically delineates the prognostic, biological, and translational significance of cytokine-related gene expression across pan-cancer cohorts. By integrating transcriptomic profiles with survival data, we established and validated a robust 16-gene cytokine-related prognostic signature that consistently stratifies patients into distinct risk groups with significantly different clinical outcomes across multiple cancer types and independent cohorts. The cytokine-related risk score was strongly associated with key malignant phenotypes, including enhanced epithelial-mesenchymal transition, cell cycle activation, and extracellular matrix remodeling, providing mechanistic insight into its unfavorable prognostic impact. Importantly, immunohistochemical validation further confirmed the upregulation of representative risk-associated genes, PANX1 and FRMD8, at the protein level in multiple tumor tissues compared with corresponding normal tissues, reinforcing the biological credibility and translational relevance of the signature. Moreover, integration of the cytokine-related risk score with clinical variables into a nomogram improved individualized survival prediction and demonstrated potential clinical utility. Collectively, these findings highlight cytokine-related transcriptional programs as critical determinants of tumor aggressiveness and patient prognosis and support their further development as clinically relevant biomarkers for pan-cancer risk stratification and precision oncology.