Abstract
BACKGROUND: Survival outcomes among patients with colorectal cancer (CRC) often differ despite identical disease stages, partly due to variations in nutritional and immune status. Malnutrition can impair immune defense, exacerbate inflammatory responses, and influence tumor progression, ultimately contributing to a poorer prognosis. However, current clinical prognostic systems rarely integrate nutritional immune indicators with tumor biomarkers, limiting the application of nutritional intervention in CRC management. This study aimed to develop a nutritional immune risk score (NIRS) model to improve long-term prognostic evaluation in patients with CRC. METHODS: In this retrospective study, 892 inpatients with primary CRC who underwent curative resection in 2017 were included and followed until 2023. Unsupervised learning was applied to nutritional and tumor biomarkers for feature extraction and patient stratification. K-means clustering was used to identify subgroups, and principal component analysis was used to derive composite features, which were then used to construct the NIRS model for long-term prognostic assessment. RESULTS: Four variables-prognostic nutritional index (PNI), carcinoembryonic antigen (CEA), carbohydrate antigen 19-9 (CA19-9), and carbohydrate antigen 72-4 (CA72-4)-were selected for model construction. The final model was defined as: NIRS = 0.572 × PNI - 0.101 × CEA - 0.412 × CA19-9 - 0.028 × CA72-4. Using an optimal cutoff value of 21.34, patients were stratified into a low-risk group and a high-risk group. The Kaplan-Meier analysis showed that patients in the low-risk group had significantly better overall survival than those in the high-risk group (p < 0.001). Multivariable Cox regression analysis indicated that the high-risk group had a 1.72-fold higher mortality risk than the low-risk group (HR = 1.72, 95% CI: 1.34-2.21, p < 0.001). In addition, PNI was negatively correlated with maximum tumor diameter in both survivors and non-survivors (survivors: r = -0.434, p < 0.001; non-survivors: r = -0.214, p < 0.001). Locally estimated scatterplot smoothing (LOESS) analysis further demonstrated that among patients with PNI ≥ 50, survivors had smaller tumors than non-survivors, whereas the opposite pattern was observed among patients with PNI < 50. CONCLUSION: We developed a novel NIRS for long-term prognostic assessment in patients with CRC. The NIRS model demonstrated robust risk stratification and potential clinical utility. PNI may serve as a complementary factor to refine risk classification, and its interaction with maximum tumor diameter may improve the sensitivity and precision of prognostic assessment across different nutritional immune states.