- Department of Otolaryngology, Head and Neck Surgery, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
Objectives: Given the high postoperative recurrence of locally advanced laryngeal squamous cell carcinoma (LSCC) and American Joint Committee on Cancer (AJCC) staging system prediction limitations, this study aims to construct and validate a postoperative recurrence-free survival (RFS) prediction model using multimodal feature fusion and explore data integration strategies to enhance prediction efficacy.
Methods: Data from 278 patients diagnosed with locally advanced LSCC between 2013 and 2024 were collected retrospectively. These data were then separated into a training dataset (n = 196) and a validation dataset (n = 82), using a near 7:3 allocation strategy. By integrating clinicopathological features, preoperative blood markers, and enhanced computed tomography imaging data, we constructed clinicopathological (Clinic-score), radiomics (Rad-score), and two fusion models: feature-level (FF-Model) and decision-level (DF-Model). Model performance was evaluated using the concordance index, time-dependent area under the receiver operating characteristic curve, calibration curve, and decision curve analyses. Improvement in model discriminative ability was assessed using continuous net reclassification improvement (cNRI) and integrated discrimination improvement (IDI).
Results: At 24.5 months median follow-up, 95 patients (34.2%) experienced recurrence. In the validation set, the DF-Model significantly outperformed the FF-Model, Rad-score and Clinic-score models, and AJCC stages. Additionally, the DF-Model demonstrated superior calibration and clinical utility, better prediction of 1-year, 3-year, and 5-year RFS through cNRI/IDI analysis, and excellent risk stratification across datasets, AJCC stages, and tumor locations.
Conclusion: The multimodal prediction DF-Model effectively integrates multi-source heterogeneous information, significantly improving the prediction accuracy of postoperative RFS in locally advanced LSCC, outperforming the FF-Model, single-modal models, and AJCC staging system, and demonstrating its potential clinical translational value.
1 Introduction
Laryngeal squamous cell carcinoma (LSCC), one of the most common malignant tumors of the head and neck region, is increasing in incidence among males annually, accounting for most of the approximately 180,000 new cases of laryngeal cancer globally each year (1, 2). Notably, approximately 43.1–65% of patients are diagnosed with locally advanced stage (III, IVa, and IVb) disease at initial presentation, resulting in a five-year disease-free survival rate of only 50–65% (3–7). Despite the widespread application of comprehensive treatment strategies, including surgery, radiotherapy, and chemotherapy, approximately 30–40% of patients experience local tumor recurrence or distant metastatic spread after surgery, which can substantially compromise long-term survival and quality of life (4, 7, 8).
Currently, clinical practice relies primarily on the TNM staging system of the American Joint Committee on Cancer (AJCC) for prognostic assessment. However, its reliance on anatomical criteria overlooks critical biological and systemic factors—such as tumor heterogeneity, host systemic status, and treatment quality—that significantly influence outcomes. Consequently, its predictive accuracy for recurrence, as measured by the concordance index (C-index), is suboptimal (C-index<0.65), limiting the utility of precise individualized risk stratification (9, 10). In addition, systemic inflammatory indicators (such as the neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR), and lymphocyte-to-monocyte ratio (LMR)) and nutritional prognostic indices (such as the prognostic nutritional index (PNI)) have attracted considerable attention; however, the predictive efficacy of single markers is limited and susceptible to various factors (11, 12). Novel molecular markers (such as epidermal growth factor receptor overexpression (13), WRAP53β (14), sex hormone receptors estrogen receptor-β and progesterone receptor (15), and p53 mutation (16)) have shown significant prognostic value; however, their clinical application requires further validation.
With advancements in artificial intelligence (AI) and precision medicine, the synergistic use of diverse data modalities—such as clinical, imaging, and molecular profiles—combined with AI algorithms—has shown great promise in enhancing prognosis prediction. In areas such as high-grade serous ovarian cancer (17), glioblastoma (18), thyroid cancer (19), renal cell carcinoma (20), breast cancer (21), and colorectal cancer (22), the area under the receiver operating characteristic (ROC) curve (AUC) values and C-indices (AUC > 0.8; C-index > 0.7) of multimodal data fusion models generally outperform traditional methods. However, research on predicting the risk of postoperative recurrence in locally advanced LSCC still has limitations, most of which are confined to application scenarios of a single type or single-modality data source, such as relying solely on clinicopathological features (23) or radiomics parameters (24–26), with AUC values generally<0.8, making it difficult to meet the need for precise stratification of patients.
Therefore, this study aimed to develop a comprehensive prediction model for postoperative recurrence-free survival (RFS) in locally advanced LSCC by synthesizing preoperative data from multiple domains—clinicopathological, laboratory, and radiomic—thereby leveraging the strengths of multimodal fusion. By comparing the efficacy of feature- and decision-level fusion strategies, this study sought to achieve precise stratification of postoperative recurrence risk, thereby providing data support for the formulation of personalized treatment plans (such as adjuvant radiotherapy dose adjustment) and dynamic prognosis management. We compared feature- and decision-level fusion because they represent fundamentally different integration strategies: the former concatenates features early (potentially capturing interactions but risking overfitting), while the latter combines model predictions, preserving modality-specific patterns and enhancing interpretability. This comparison helps identify the optimal architecture for clinical prediction in LSCC.
2 Materials and methods
2.1 Study participants
Following the principles of the TRIPOD-AI checklist (27), we conducted a retrospective cohort study of patients with locally advanced LSCC at the First Affiliated Hospital of Guangxi Medical University (Jan 2013–Jan 2024). Inclusion was determined by histopathologic staging per AJCC 8th edition (Stage III–IVb) (9), availability of more than six months of postoperative monitoring, documented results for standard blood counts as well as liver and kidney function indicators, and qualified preoperative contrast-enhanced multiplex spiral computed tomography (CT) scan images. The exclusion criteria were incomplete medical records; history of surgery, radiotherapy, or chemotherapy prior to surgery; presence of active infections, chronic inflammation, hematological diseases, or autoimmune diseases; and history of other malignancies.
The sample size was calculated using the Events Per Variable (EPV) method: N = (EPV × number of predictor variables/recurrence rate) × (1 + efficiency rate). With 7 expected variables, a 30% 5-year recurrence rate, EPV = 10, and 10% inefficiency, a minimum of 256 patients was required. After applying criteria, 278 patients were included, meeting statistical power needs. The study adhered to the Declaration of Helsinki and was approved by the Institutional Ethics Committee (No. 2025-E0564). Informed consent was waived for this retrospective study, and all data were anonymized. The study workflow is shown in Figure 1.

Figure 1. Methodological framework of the study design. (A) Multi-Scale Feature Fusion-Based Model; (B) Radiomic Model.
2.2 Acquisition of clinical variables
Clinical data were extracted from electronic health records, including age, sex, comorbidities (e.g., hypertension), smoking/alcohol history, preoperative blood count, liver/kidney function, and postoperative pathological features (tumor location, margin status, vascular invasion, lymph node metastasis, differentiation, and TNM stage). Based on preoperative peripheral blood test results, we calculated the following systemic inflammation-related biomarkers: NLR, PLR, LMR, Systemic Immune-Inflammation Index (SII) = platelet count × NLR, PNI = lymphocyte count × 5 + albumin, Advanced Lung Cancer Inflammation Index (ALI) = body mass index (BMI) × albumin/NLR, and Systemic Inflammation Response Index (SIRI) = monocyte count × NLR. Blood tests used the Beckman Coulter/LH 780 and Werfen/ACL TOP 750LAS, with collection between 6:00 and 10:00 AM. Smoking history was defined as >1 cigarette/day for >6 months pre-surgery; drinking history as ≥72g alcohol/week for >6 months pre-surgery. The AJCC 8th edition (2017) criteria were applied for pathological staging.
2.3 Image acquisition and preprocessing
Contrast-enhanced CT scans (skull base to supraclavicular) were performed. Venous phase images were acquired 60–90 s after injecting non-ionic contrast (iopromide/iohexol, 350 mg/ml, 1.0–2.0 ml/kg, 3–5 ml/s). Scanner parameters are in Supplementary Table 1. Image preprocessing followed the Evaluation of Radiomics Research (CLEAR) checklist (28): original DICOM images were resampled to 1×1×1 mm³ isotropic resolution, gray-level discretized (bin width: 25 HU), intensity range limited to 1–500 HU, and normalized using Z-score. Window width/level was set to 350/40 HU to enhance tumor boundaries. The tumor volume of interest (VOI) was outlined at the beginning of the study by randomly selecting CT images from 35 patients. The delineation of tumor boundaries on axial images was conducted independently by two head and neck specialists, who had accumulated 10 and 13 years of clinical expertise, using the ITK-SNAP platform (www.itksnap.org) without knowledge of other clinical data (Figure 1B). One week later, the more experienced rater repeated the segmentation. Inter-rater agreement was assessed using the intraclass correlation coefficient (ICC). Features with ICC > 0.75 were retained for analysis.
2.4 Prediction of outcomes and follow-up
RFS was defined as time from surgery to the first recurrence or metastasis or the end of follow-up (Jan 31, 2025), with time points set at 1, 3, and 5 years postoperatively. Follow-up used digital communication (WeChat), telephonic contact, and clinic visits (quarterly for years 1–3, semi-annually for years 4–5, annually thereafter), with<5% loss to follow-up.
2.5 Model development and validation
2.5.1 Dataset partitioning
To ensure reproducibility, we split the dataset into two subsets — a training set (n = 196) and a validation set (n = 82) — using a 7:3 partitioning strategy with a fixed random seed (seed = 42). The training set was used to train and optimize the model, whereas the validation set was used to evaluate the generalization ability of the model.
2.5.2 Development of clinical-pathological model (Clinic-score)
Based on 35 baseline clinicopathological characteristics and peripheral blood markers of patients in the training set, variables were screened using univariate Cox regression (p< 0.05), and the optimal variable combination was selected in conjunction with the Akaike Information Criterion (AIC) (29) to construct a multivariate Cox regression model, which was defined as the clinical-pathological model (Clinic-score). The results are presented as a nomogram, and a web-based calculator was developed.
2.5.3 Development of radiomics models (Rad-score)
Radiomics features were extracted from preprocessed tumor VOIs in the training cohort using PyRadiomics, encompassing ten feature classes: first-order statistics, 3D morphological features, gray-level co-occurrence matrix (GLCM), gray-level size zone matrix (GLSZM), gray-level run length matrix (GLRLM), neighborhood gray-tone difference matrix (NGTDM), Hessian-based features, fractal features, and topological features. Univariate Cox regression (P< 0.05) identified prognostic features, which were further refined using the least absolute shrinkage and selection operator (LASSO) with 10-fold cross-validation to minimize overfitting. The optimal regularization parameter (λ) was selected by minimizing cross-validated partial likelihood deviance, retaining only features with non-zero coefficients. A multivariable Cox model was then fitted to these features to compute the radiomics score (Rad-score):
Where is the value of the i-th radiomics feature, is its corresponding regression coefficient from the multivariable Cox model, and k is the number of selected features.
2.5.4 Development of fusion models
We compared feature-level and decision-level fusion strategies to determine the optimal approach for integrating multimodal data (clinical, blood, and radiomics). Feature-level fusion combines raw features, while decision-level fusion integrates model outputs, offering potentially greater robustness.
Feature-level fusion model (FF-Model)
1. Variable selection: This model involves the direct concatenation of clinicopathological features, peripheral blood markers, and radiomics features, which are selected through univariate Cox regression during the construction of the radiomics model. Subsequently, a LASSO-Cox regression approach incorporating 100 rounds of 10-fold cross-validation is applied to select relevant predictors.
2. Fusion strategy: Early fusion — raw features from different modalities are concatenated into a single vector.
3. Model construction: A multivariable survival model is constructed based on Cox regression. The resulting linear predictor is defined as the FF-Score.
Decision-level fusion model (DF-Model)
1. Variable selection: This model uses the outputs of two pre-trained submodels as input variables:
(i) Clinic-score: derived from the clinical-pathological model (Section 2.5.2);
(ii) Rad-score: derived from the radiomics model (Section 2.5.3).
2. Fusion strategy: Late fusion — predictions (scores) from separate submodels are combined at the decision level.
Model construction: The Clinic-score and Rad-score are combined in a multivariable Cox proportional hazards model as covariates. The model is defined as: . Where h0(t) is the baseline hazard, and β1, β2 are the maximum likelihood estimates for the respective regression coefficients.
2.5.5 Model evaluation and comparison
Model performance was evaluated based on discrimination, calibration, and clinical utility. Discrimination was assessed using the C-statistic and time-dependent AUC at 1, 3, and 5 years for RFS. Calibration was evaluated using calibration plots. Clinical utility was assessed via decision curve analysis (DCA), comparing models against each other and the AJCC 8th Edition TNM staging system. Additionally, the degree of improvement in the fusion model was assessed using continuous net reclassification improvement (cNRI) and integrated discrimination improvement (IDI). Finally, the optimal cutoff for the best-performing model was determined using X-tile software (30), stratified by AJCC stage and tumor location, to validate its risk stratification and generalizability across subgroups.
2.6 Statistical methods
Categorical variables (such as tumor location), are summarized with counts and proportions, while group differences were assessed using either the chi-square test or Fisher’s exact test. For continuous variables (such as NLR), the data are expressed as mean ± standard deviation (mean ± SD) or median and interquartile range (IQR), contingent upon whether the data adhered to a normal distribution. An independent samples t-test or Mann–Whitney U test was used to evaluate differences between groups. The statistical processing was carried out using R version 4.2.3 and Python 3.6, with key functions drawn from these packages: glmnet 4.1.8, pROC 1.18.5, rms 6.7.1, dplyr 1.1.4, survival 3.7.0, and timeROC 0.4. Statistical significance was set at p< 0.05.
3 Results
3.1 Patient characteristics
A total of 278 patients with locally advanced LSCC were enrolled in this study, including 106 (38.1%) patients with stage III, 170 (61.1%) with stage IVa, and two (0.7%) with stage IVb disease. The median follow-up duration was 24.5 (8.3–54.8) months. As of the follow-up cutoff date of January 31, 2025, 95 (34.2%) patients experienced recurrence or distant metastasis. The median time from the last treatment to recurrence or metastasis was 10 (6–21) months. The cumulative RFS rates at 1, 3, and 5 years were 79.5%, 70.1%, and 67.3%, respectively. The cumulative recurrence rate was 35.2% (69/196) in the training set and 31.7% (26/82) in the validation set, with no significant differences between the recurrence rates of the two sets (p > 0.05). Except for BMI (Z = 2.204, p = 0.028) and smoking history (χ² = 3.884, p = 0.049), there were no significant differences in the distribution of the remaining variables between the training and validation sets (all p > 0.05, Supplementary Table 2).
3.2 Clinic-score
Univariate Cox regression analysis identified 13 potential predictive variables (including tumor location and AJCC stage; all p< 0.05; Table 1). Based on the AIC, seven optimal variables were determined to construct the multivariate Cox regression model (Supplementary Table 3). Additionally, the model was visualized as a nomogram and deployed a web calculator at https://huangxiaoying.shinyapps.io/dynnomapp (Figure 2).

Table 1. Univariate Cox regression analysis of postoperative RFS in patients with locally advanced LSCC.

Figure 2. A clinicopathological model for predicting RFS after surgery in locally advanced LSCC. (A) Nomogram, (B) Web-based calculator.
3.3 Rad-score
Through PyRadiomics, 3232 radiomics features were extracted, with 2658 (82.24%) retained after ICC screening, including 10 categories of features such as first-order statistics (360) and Hessian matrix (780) (Supplementary Figure 1). Univariate Cox regression analysis further screened 413 significant features (p< 0.05), and LASSO ultimately retained seven non-zero coefficient features, including one first-order statistical feature, one three-dimensional shape feature, three GLCM, one GLRLM, and one GLSZM (Figure 3, Table 2). The Rad-score was computed as a linear combination of these features (Equation 1):

Figure 3. Selection of non-zero coefficient radiomics features using the least absolute shrinkage and selection operator (LASSO) regression model. (A) LASSO regularization path diagram; (B) C-index coefficient plot using 10-fold cross-validation; (C) 7 selected radiomics features and their weight coefficients; (D) correlation clustering heatmap of 7 radiomics features.
3.4 Fusion model construction
3.4.1 FF-Model
To develop the FF-Model, we integrated 13 clinicopathological variables with 413 radiomic features, resulting in a 426-dimensional dataset. Subsequently, LASSO-Cox regression identified 16 significant predictors, including PNI, tumor location, surgical margin status, contralateral cervical lymph node metastasis, and 12 radiomics features (Supplementary Figure 2, Supplementary Table 4). Multivariate Cox regression identified seven independent predictors of RFS (p< 0.05), among which were positive surgical margins and contralateral lymph node involvement (Figure 4).
3.4.2 DF-Model
The DF-Model formula based on multifactorial Cox regression is as follows:
Both the Clinic-score and Rad-score significantly impacted the postoperative RFS of locally advanced LSCC (p< 0.001; Figures 5, 6).
3.5 Model performance comparison
3.5.1 Discrimination
In the training set, the C-index of the DF-Model was 0.847 (95% CI: 0.811–0.884), which was significantly higher than the Clinic-score (0.723; p< 0. 001), Rad-score (0.828; p = 0.099), and AJCC stage (0.608; p< 0.001). The C-index of the FF-Model was 0.878 (95% CI: 0.838–0.917), which was significantly higher than that of the DF-Model (p = 0.024). In the validation set, the DF-Model achieved a C-index of 0.826 (95% CI: 0.763–0.889), showing a statistically significant improvement over the FF-Model (0.741; p = 0.047), Rad-score (0.734; p = 0.033), Clinic-score (0.723; p = 0.002), and AJCC stages (0.58; p< 0.001, Table 3).
The ROC analysis showed no significant difference between the AUC values of the FF-Model and DF-Model in the prediction of 1/3/5-year RFS in the training set (p > 0.05; Figures 7A–C). In the validation set, the AUC of the DF-Model was significantly higher than that of the FF-Model for predicting 3-year and 5-year RFS (all p = 0.022). no statistically significant difference was observed in the predictive performance of 1-year RFS across the models (p = 0.206; Figures 7D–F).

Figure 7. ROC comparison of four models for RFS prediction. (A–C) Training set; (D–F) Validation set.
3.5.2 Calibration
Calibration curve analysis demonstrated that in predicting 1-year, 3-year, and 5-year RFS, the calibration curves of the DF-Model in both the training and validation sets were closer to the ideal diagonal line than those of the FF-Model, Rad-score, and Clinic-score (Figure 8), indicating superior calibration performance over the other models.
3.5.3 Clinical utility
DCA showed that both fusion models provided greater clinical utility than the single-modality models and the AJCC staging system for predicting 1-, 3-, and 5-year RFS (Figure 9).

Figure 9. Decision curve analysis of five models for RFS prediction. (A–C) Training set; (D–F) Validation set.
In the training set, the FF-Model yielded the highest net benefit across most threshold probabilities. In the validation set, the DF-Model generally outperformed other models, particularly for 3- and 5-year predictions. Although the FF-Model showed slightly higher net benefit at some thresholds for 1-year RFS, the DF-Model remained superior to all single-modality models.
3.5.4 cNRI and IDI tests
The DF-Model exhibited better discriminative ability than the single-modal and the FF-Model in both the cNRI and IDI tests. In the 1-year prediction, the cNRI increased by 32.6% (p = 0.040), 51.3% (p = 0.020), and 23.4% (p< 0.001) compared to the FF-Model, Rad-score, and Clinic-score, respectively. Although the IDI gain did not reach statistical significance (p = 0.079), it transcended the minimal clinically important difference (MCID = 5%) (31) and thus, might have potential value. In the 3-year prediction, the cNRI showed significant improvement compared to the Rad-score (51.5%; p = 0.020) and Clinic-score (70.1%; p< 0.001). Additionally, the IDI exhibited superior performance compared to the Clinic-score (27.3%; p = 0.02) and Rad-score (19.1%; p< 0.001). Regarding the 5-year prediction, an even more robust performance was observed, with the cNRI increasing to 67.0% (p< 0.001) compared to the Clinic-score, and 47.8% (p = 0.020) compared to the Rad-score. The IDI exhibited a similar trend (Table 4).

Table 4. Comparison of the DF-Model with the FF-Model, Rad-score, and Clinic-score using IDI and cNRI metrics.
3.5.5 Subgroup analysis
X-tile software was used to determine the cutoff value of the DF-Model, which was 81.1, to divide the training set into low-risk (<81.1) and high-risk (≥81.1) subgroups. Kaplan-Meier analysis showed significantly worse recurrence-free survival (RFS) in the high-risk group compared to the low-risk group in both the training and validation cohorts (all p< 0.001; Figure 10). The prognostic value of the DF-Model remained significant across subgroups defined by AJCC stage (III vs. IV) and tumor location (glottic vs. non-glottic) (all p< 0.001; Figure 11).

Figure 10. Survival analysis using the fusion model threshold (81.1). (A) Training set; (B) Validation set.

Figure 11. Subgroup survival analysis using the DF-Model threshold (81.1). (A) Stage III; (B) Stage IVa/IVb; (C) Non-glottic; (D) Glottic.
4 Discussion
4.1 State of the art in multimodal data fusion applications
Locally advanced LSCC has a high postoperative recurrence rate, and the traditional AJCC staging system is not sufficient for precise prognosis prediction, making it difficult to develop personalized treatment plans. In this case, multimodal data fusion provides a different method for improving prediction effectiveness. Different data modalities complementarily reflect the biological behavior of tumors. Macroscopic invasiveness (such as rhabdomyolysis invasion and lymph node metastasis burden) and preoperative blood markers (such as PNI) characterize the host’s systemic inflammatory and nutritional status, whereas radiomic features (such as wavelet transform texture and gray entropy) quantify the heterogeneity of the tumor microenvironment. Studies have attempted to integrate multimodal data to improve predictive efficacy in head and neck squamous cell carcinoma (HNSCC). For example, Tseng et al. (32) integrated clinical, pathological, and genetic variation data to construct an elastic net Cox model to predict the survival risk in patients with oral cancer. Wang et al. (33) integrated radiomic and pathomic features to develop a Particle Swarm Optimization–Support Vector Machine model aimed at assessing the response to neoadjuvant chemotherapy among individuals diagnosed with nasopharyngeal carcinoma. Yin et al. (34) integrated multiple immuno-omics data to develop a Computational Model for Predicting Immunotherapy Response model for identifying patient populations sensitive to immunotherapy and chemotherapy. Cavalieri et al. (35) established a large-scale multiomics database, BD2Decide, to provide a rich data resource for head and neck cancer research. Additionally, studies have enhanced model performance by integrating radiomic features from different imaging modalities. For instance, Tomita et al. (36) utilized radiomic data obtained from multiple magnetic resonance imaging (MRI) scans to develop a deep learning model designed to predict the likelihood of 2-year progression-free survival (PFS) in patients with laryngeal and hypopharyngeal malignancies. Lin et al. (37) constructed a multimodal radiomics scoring model based on preoperative MRI images of patients with advanced sinonasal squamous cell carcinoma to predict early disease recurrence. Huynh et al. (38) evaluated the performance of traditional radiomics versus convolutional neural network (CNN) models in forecasting patient survival outcomes in HNSCC, and found that the predictive capability of CNN models improved when combined with clinical and radiomic features. Li et al. (39) developed multiple prediction models and demonstrated that the DL-Model achieved higher AUC, accuracy, and specificity in the validation set. Wang et al. (40) also showed that the DL-Model based on multimodal radiomic features achieved the highest AUC values (0.89–0.90) across all datasets, along with optimal sensitivity (82–88%) and specificity (79–85%).
4.2 Core findings and comparison with other studies
In the field of laryngeal cancer prognostic prediction, previous studies have mainly focused on single-modality data or small patient cohorts. Zhong et al. (25) used positron emission tomography (PET)-CT metabolic parameters to construct a random forest model to predict disease progression; Choi et al. (41) combined radiomic scores with clinical variables, which increased the C-index of survival prediction to 0.958; Lin et al. (42) found that intra-tumoral/peri-tumor features during the mid-radiation process performed better than pre-radiation models in prediction; Agarwal et al. (43) used pre-treatment CT image data of 60 patients treated with chemoradiotherapy and found that the entropy of medium-filtered texture features is an independent predictor. However, these studies have issues such as incomplete selection of regions of interest and some cases of hypopharyngeal cancer, which makes it difficult to apply these findings elsewhere. Al-Ibraheem et al. (44) used PET/CT data from 68 patients with laryngeal cancer to build a Cox model demonstrating that total lesion glycolysis and metabolic tumor volume have independent predictive values for PFS. Nakajo et al. (45) developed an RSF model using PET/CT features of 49 patients to predict PFS. Chen et al. (24) built a radiomics nomogram with a C-index of 0.913 using CT data from 136 cases of LSCC; however, these studies included early stage laryngeal cancer, which could add some selective bias. Rajgor et al. (26) studied 72 patients with advanced laryngeal cancer, using shape compactness and gray-level zone length matrix−gray-level nonuniformity modeling to predict 5-year survival, achieving a C-index of 0.759, outperforming the clinical model’s C-index of 0.655. Nevertheless, this investigation was conducted at a single institution and involved a limited number of participants; furthermore, its findings have yet to undergo external validation.
In contrast, this study provided an explicit definition and included 278 individuals with locally advanced LSCC, targeting the prediction of prognosis within this high-risk group of patients. Additionally, we included clinicopathological features such as surgical margin status and the ratio of positive lymph nodes, pre-PNI, and CT radiomics features to capture complementary information. The preoperative PNI indicates the patient’s nutritional and immune status, striated muscle invasion represents tumor aggressiveness, the number of positive nodes represents the extent of tumor metastasis, margin status reflects operation radicality, and radiomics features further quantify tumor microenvironment heterogeneity; therefore, the model enables a comprehensive, multidimensional evaluation of the “tumor characteristics- host status- treatment” framework. This represents a shift beyond the conventional AJCC staging system, which was not designed to predict recurrence at the individual patient level and primarily reflects anatomical extent rather than underlying biological behavior.
The radiomic model incorporated features reflecting diverse aspects of tumor phenotype. Tumor size was represented by the three-dimensional maximum diameter (original_shape_Maximum3DDiameter), with larger values generally associated with more advanced disease and poorer prognosis (47). Tumor intensity distribution asymmetry was captured by first-order skewness (gradient_firstorder_Skewness), where high positive values indicate a more aggressive and heterogeneous tumor phenotype (43). Tumor heterogeneity was further quantified using wavelet-based texture features (e.g., wavelet_LLL_glszm_ZoneEntropy), which measure gray-level inhomogeneity; higher values suggest greater internal complexity and are linked to adverse outcomes (43). Notably, four of the seven radiomic features in the final model were derived from wavelet transformation, underscoring the importance of capturing intratumoral heterogeneity in prognostic prediction (46).
In the present study, the DF-Model achieved a C-index of 0.826 (95% CI: 0.763–0.889) in validation, exceeding the FF-Model (0.741), Rad-score (0.734), and Clinic-score (0.657). Furthermore, it exhibited superior performance in calibration curve analysis and DCA. The cNRI and IDI metrics indicated that the DF-Model demonstrated superior performance in forecasting RFS at 1-, 3-, and 5-year intervals compared to alternative models. In particular, the cNRI enhancement for the 1-year prediction relative to the FF-Model reached 32.6% (p< 0.05). Although the improvements in the 3–5-year predictions did not achieve statistical significance, they exceeded the MCID and may hold clinical value in high-risk recurrence cases. The improving effect of the DF-Model compared with the FF-Model diminishes over time, aligning with the clinical observation that the likelihood of recurrence in locally advanced LSCC decreases with time.
Compared with AJCC staging (training set C-index: 0.608; validation set C-index: 0.58), the two fusion models constructed in this study, the single Rad-score and the single Clinic-score, demonstrated significant advantages in terms of discriminative ability and clinical utility. The DF-Model exhibited robust risk stratification capabilities across different datasets, AJCC staging, and tumor location subgroups (log-rank p< 0.001). For instance, in the stage III subgroup, the DF-Model achieved a significant difference in 3-year RFS (high-risk group: 31% vs. low-risk group: 89%), which could not be achieved by AJCC staging alone. This result indicates that the DF-Model transcends the limitations of anatomical staging and can identify truly low-risk groups among patients with advanced disease. Moreover, DCA results indicate that the DF-Model yields a greater net benefit over most of the clinically relevant threshold probability ranges. However, in the 1-year RFS prediction (threshold probability 0.18–0.4), the FF-Model is slightly superior, possibly because short-term recurrence risk is more dependent on the quality of treatment (such as surgical margin status), whereas medium- and long-term recurrence is driven by tumor heterogeneity and the immune microenvironment.
The enhanced risk stratification capability of the DF-Model holds significant clinical interpretability and direct implications for postoperative management. By accurately identifying high-risk patients beyond conventional staging, clinicians can consider intensifying adjuvant therapy—such as adding chemotherapy to radiotherapy or extending radiation fields—for those most likely to benefit, potentially improving survival outcomes. Conversely, low-risk patients identified by the DF-Model may be candidates for de-escalated treatment or less frequent follow-up, thereby reducing the burden of overtreatment, minimizing long-term toxicities (e.g., dysphagia, xerostomia, voice deterioration), and preserving quality of life. This ability to refine risk assessment within the same AJCC stage (e.g., distinguishing high- from low-risk Stage III patients) enables a more biologically driven, personalized approach to postoperative care, moving beyond anatomical staging alone.
4.3 Efficacy difference mechanism of multimodal fusion strategy
The differences in the results of the two fusion strategies can be attributed to their inherent characteristics. At the feature level, feature-level fusion is obtained by combining raw features directly before modeling, so the “information dilution effect” and “interaction noise interference” occur (48). The substantial number of radiomic features (413) compared to clinical feature data (13) suggests a high likelihood of overshadowing certain clinical information, such as the positive ratio of lymphatic metastasis, during the process of dimensionality reduction, which can adversely affect model performance. This phenomenon is evident in the results obtained from applying the FF-Model to the validation set, which demonstrated a decrease in the C-index by Δ = 0.137. These findings challenge the assertion that “fusion is always better than single-modality.”
The DF-Model combines the prediction probabilities or scores of individual modality models to build an ensemble model that preserves independent prediction information (49). In the DF-Model, both the Clinic-score and Rad-score are independent predictors of RFS. Analysis of β-values indicated that 80% of the DF-Model’s reliance is on the Rad-score, suggesting that tumor microenvironment heterogeneity significantly influences recurrence, whereas clinicopathological features act as supplementary factors. Patients misclassified as low-risk by the clinical model yet flagged as high-risk by the radiomic model can be accurately recognized by the DF-Model. These patients may benefit from intensive postoperative imaging follow-up, radiotherapy, and immunotherapy. For low-risk patients, the aim is to prevent overtreatment, conserve medical resources, and alleviate the burden of the illness. In addition, the DF-Model eliminates the demand for extensive data preparation by integrating the probable outputs from the Clinic-score and Rad-score, thereby decreasing technical prerequisites and enhancing applicability in resource-limited places, such as primary hospitals. Patients with lower preoperative PNI levels may have improved outcomes with improved nutrition (such as a protein-rich diet supplemented with ω-3 fatty acids). Further studies are required to verify this hypothesis.
4.4 Limitations and future directions
Several limitations should be acknowledged in this research, such as a limited participant count (n = 278), reliance on single-center data, and the absence of external validation to assess generalizability across different clinical environments. The strategies for feature selection and integration may require further refinement, as there is a risk of losing important features. Both feature-level and decision-level fusion may not fully optimize data utilization, leading to predictive uncertainty. Promoting the model is difficult because the model’s complexity and the cost associated with CT radiomics feature extraction hinder its clinical translation.
To facilitate the transition from research to clinical practice, essential next steps should prioritize multicenter prospective validation to confirm generalizability, followed by prospective clinical trials designed to evaluate whether model-guided management leads to tangible improvements in patient outcomes—such as recurrence, survival, and quality of life. Concurrently, integrating the model with established molecular biomarkers, including PD-L1 and HPV status, could refine risk stratification, while the development of an accessible web-based calculator would support real-time, bedside decision-making. Beyond these immediate translational priorities, future research should explore advanced feature engineering methods such as transformer networks, the integration of multiomics data—including genomics and proteomics—to enhance predictive accuracy, dynamic risk updating based on postoperative follow-up, and strategies for clinician education and seamless integration into routine clinical workflows.
5 Conclusion
The multimodal feature fusion model developed according to the decision-level fusion strategy (DF-Model) enhances the prediction of postoperative RFS in locally advanced LSCC. This improvement is primarily attributed to the integration of multimodal features derived from clinical pathology, PNI, and CT radiomics for the construction of the model. Its performance surpassed those of the FF-Model, single-modality model, and traditional AJCC staging. The model displayed robust risk stratification capabilities across various AJCC stages and subgroups according to tumor location. Further investigations are encouraged to emphasize multicenter validation, enhancement of predictive algorithms, and integration into clinical workflows in order to support real-world deployment and individualized therapeutic strategies.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
Ethics approval was obtained from the host institution’s review committee (NO.2025-E0564), and all patients signed an informed consent form before the examination.
Author contributions
FZ: Writing – original draft, Formal analysis, Visualization, Resources, Methodology, Conceptualization, Data curation. XH: Supervision, Project administration, Methodology, Writing – review & editing, Conceptualization. JML: Software, Methodology, Data curation, Writing – original draft. JH: Methodology, Data curation, Writing – original draft. JXL: Writing – original draft, Methodology, Data curation. GC: Resources, Writing – original draft, Data curation. ZZ: Supervision, Writing – review & editing, Project administration, Methodology, Conceptualization.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Guangxi Medical University Undergraduate Innovation and Entrepreneurship Training Program (No. X202510598314), the Guangxi Medical University Clinical Discipline Special Project for Educational Reform (No. 2025LCJG05), and the Guangxi Medical University Undergraduate Education and Teaching Reform Project (No. 2025XJGY18).
Acknowledgments
Although this study is retrospective in nature, the authors would like to express their sincere gratitude to all participants for their contributions. We also thank Professor Su Jiping from the First Affiliated Hospital of Guangxi Medical University for his valuable support and guidance throughout this research. We would like to thank Editage (www.editage.cn) for English language editing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1685737/full#supplementary-material
References
1. Chen B, Zhan Z, Fang W, Zheng Y, Yu S, Huang J, et al. Long-term trends and future projections of larynx cancer burden in China: a comprehensive analysis from 1990 to 2030 using GBD data. Sci Rep. (2024) 14:26523. doi: 10.1038/s41598-024-77797-6
2. Huang J, Chan SC, Ko S, Lok V, Zhang L, Lin X, et al. Updated disease distributions, risk factors, and trends of laryngeal cancer: a global analysis of cancer registries. Int J Surg. (2024) 110:810–9. doi: 10.1097/JS9.0000000000000902
3. Groome PA, O'Sullivan B, Irish JC, Rothwell DM, Schulze K, Warde PR, et al. Management and outcome differences in supraglottic cancer between Ontario, Canada, and the Surveillance, Epidemiology, and End Results areas of the United States. J Clin Oncol. (2003) 21:496–505. doi: 10.1200/JCO.2003.10.106
4. Haapaniemi A, Väisänen J, Atula T, Alho OP, Mäkitie A, and Koivunen P. Predictive factors and treatment outcome of laryngeal carcinoma recurrence. Head Neck. (2017) 39:555–63. doi: 10.1002/hed.24642
5. Li MM, Zhao S, Eskander A, Rygalski C, Brock G, Parikh AS, et al. Stage migration and survival trends in laryngeal cancer. Ann Surg Oncol. (2021) 28:7300–9. doi: 10.1245/s10434-021-10318-1
6. Lin Z, Lin H, Chen Y, Xu Y, Chen X, Fan H, et al. Long-term survival trend after primary total laryngectomy for patients with locally advanced laryngeal carcinoma. J Cancer. (2021) 12:1220–30. doi: 10.7150/jca.50404
7. Đokanović D, Gajanin R, Gojković Z, Marošević G, Sladojević I, Gajanin V, et al. Clinicopathological characteristics, treatment patterns, and outcomes in patients with laryngeal cancer. Curr Oncol. (2023) 30:4289–300. doi: 10.3390/curroncol30040327
8. Brandstorp-Boesen J, Sørum Falk R, Boysen M, and Brøndbo K. Impact of stage, management and recurrence on survival rates in laryngeal cancer. PloS One. (2017) 12:e0179371. doi: 10.1371/journal.pone.0179371
9. Huang SH and O'Sullivan B. Overview of the 8th edition TNM classification for head and neck cancer. Curr Treat Options Oncol. (2017) 18:40. doi: 10.1007/s11864-017-0484-y
10. Lydiatt WM, Patel SG, O'Sullivan B, Brandwein MS, Ridge JA, Migliacci JC, et al. Head and Neck cancers-major changes in the American Joint Committee on cancer eighth edition cancer staging manual. CA Cancer J Clin. (2017) 67:122–37. doi: 10.3322/caac.21389
11. Fu Y, Liu W, OuYang D, Yang A, and Zhang Q. Preoperative neutrophil-to-lymphocyte ratio predicts long-term survival in patients undergoing total laryngectomy with advanced laryngeal squamous cell carcinoma: A single-center retrospective study. Med (Baltimore). (2016) 95:e2689. doi: 10.1097/MD.0000000000002689
12. Chen H, Song S, Zhang L, Dong W, Chen X, and Zhou H. Preoperative platelet-lymphocyte ratio predicts recurrence of laryngeal squamous cell carcinoma. Future Oncol. (2020) 16:209–17. doi: 10.2217/fon-2019-0527
13. Kontić M, Milovanović J, Čolović Z, Poljak NK, Šundov Ž, Sučić A, et al. Epidermal growth factor receptor (EGFR) expression in patients with laryngeal squamous cell carcinoma. Eur Arch Otorhinolaryngol. (2015) 272:401–5. doi: 10.1007/s00405-014-3323-9
14. Tiefenböck-Hansson K, Haapaniemi A, Farnebo L, Palmgren B, Tarkkanen J, Farnebo M, et al. WRAP53β, survivin and p16INK4a expression as potential predictors of radiotherapy/chemoradiotherapy response in T2N0-T3N0 glottic laryngeal cancer. Oncol Rep. (2017) 38:2062–8. doi: 10.3892/or.2017.5898
15. Atef A, El-Rashidy MA, Elzayat S, and Kabel AM. The prognostic value of sex hormone receptors expression in laryngeal carcinoma. Tissue Cell. (2019) 57:84–9. doi: 10.1016/j.tice.2019.02.007
16. Yuen AP, Lam KY, Choy JT, Ho WK, and Wei WI. The clinicopathological significance of bcl-2 expression in the surgical treatment of laryngeal carcinoma. Clin Otolaryngol Allied Sci. (2001) 26:129–33. doi: 10.1046/j.1365-2273.2001.00441.x
17. Boehm KM, Aherne EA, Ellenson L, Nikolovski I, Alghamdi M, Vázquez-García I, et al. Multimodal data integration using machine learning improves risk stratification of high-grade serous ovarian cancer. Nat Cancer. (2022) 3:723–33. doi: 10.1038/s43018-022-00388-9
18. Beig N, Bera K, Prasanna P, Antunes J, Correa R, Singh S, et al. Radiogenomic-based survival risk stratification of tumor habitat on gd-T1w MRI is associated with biological processes in glioblastoma. Clin Cancer Res. (2020) 26:1866–76. doi: 10.1158/1078-0432.CCR-19-2556
19. Yu Y, Ouyang W, Huang Y, Huang H, Wang Z, Jia X, et al. AI-Based multimodal Multi-tasks analysis reveals tumor molecular heterogeneity, predicts preoperative lymph node metastasis and prognosis in papillary thyroid carcinoma: A retrospective study. Int J Surg. (2024) 111(1):839–56. doi: 10.1097/JS9.0000000000001875
20. Schulz S, Woerl AC, Jungmann F, Glasner C, Stenzel P, Strobl S, et al. Multimodal deep learning for prognosis prediction in renal cancer. Front Oncol. (2021) 11:788740. doi: 10.3389/fonc.2021.788740
21. Mondol RK, Millar E, Sowmya A, and Meijering E. BioFusionNet: deep learning-based survival risk stratification in ER+ Breast cancer through multifeature and multimodal data fusion. IEEE J BioMed Health Inform. (2024) 28:5290–302. doi: 10.1109/JBHI.2024.3418341
22. Zhou S, Sun D, Mao W, Liu Y, Cen W, Ye L, et al. Deep radiomics-based fusion model for prediction of bevacizumab treatment response and outcome in patients with colorectal cancer liver metastases: a multicentre cohort study. EClinicalMedicine. (2023) 65:102271. doi: 10.1016/j.eclinm.2023.102271
23. Santos TS, Estêvão R, Antunes L, Certal V, Silva JC, and Monteiro E. Clinical and histopathological prognostic factors in locoregional advanced laryngeal cancer. J Laryngol Otol. (2016) 130:948–53. doi: 10.1017/S002221511600880X
24. Chen L, Wang H, Zeng H, Zhang Y, and Ma X. Evaluation of CT-based radiomics signature and nomogram as prognostic markers in patients with laryngeal squamous cell carcinoma. Cancer Imaging. (2020) 20:28. doi: 10.1186/s40644-020-00310-5
25. Zhong J, Frood R, Brown P, Nelstrop H, Prestwich R, McDermott G, et al. Machine learning-based FDG PET-CT radiomics for outcome prediction in larynx and hypopharynx squamous cell carcinoma. Clin Radiol. (2021) 76:78.e9–78.e17. doi: 10.1016/j.crad.2020.08.030
26. Rajgor AD, Kui C, McQueen A, Cowley J, Gillespie C, Mill A, et al. Computed tomography-based radiomic markers are independent prognosticators of survival in advanced laryngeal cancer: a pilot study. J Laryngol Otol. (2024) 138:685–91. doi: 10.1017/S0022215123002372
27. Collins GS, Moons K, Dhiman P, Riley RD, Beam AL, Van Calster B, et al. TRIPOD+AI statement: updated guidance for reporting clinical prediction models that use regression or machine learning methods. BMJ. (2024) 385:e078378. doi: 10.1136/bmj-2023-078378
28. Kocak B, Baessler B, Bakas S, Cuocolo R, Fedorov A, Maier-Hein L, et al. CheckList for EvaluAtion of Radiomics research (CLEAR): a step-by-step reporting guideline for authors and reviewers endorsed by ESR and EuSoMII. Insights Imaging. (2023) 14:75. doi: 10.1186/s13244-023-01415-8
29. Ikemoto K, Takahashi K, Ozawa T, and Isobe H. Akaike's information criterion for stoichiometry inference of supramolecular complexes. Angew Chem Int Ed Engl. (2023) 62:e202219059. doi: 10.1002/anie.202219059
30. Milana F, Polidoro MA, Soldani C, Franceschini B, Famularo S, Di Tommaso L, et al. Unveiling the prognostic role of blood inflammatory indexes in a retrospective cohort of patients undergoing liver resection for intrahepatic cholangiocarcinoma. Int J Surg. (2024) 110:7088–96. doi: 10.1097/JS9.0000000000001924
31. Sedaghat AR. Understanding the minimal clinically important difference (MCID) of patient-reported outcome measures. Otolaryngol Head Neck Surg. (2019) 161:551–60. doi: 10.1177/0194599819852604
32. Tseng YJ, Wang HY, Lin TW, Lu JJ, Hsieh CH, and Liao CT. Development of a machine learning model for survival risk stratification of patients with advanced oral cancer. JAMA Netw Open. (2020) 3:e2011768. doi: 10.1001/jamanetworkopen.2020.11768
33. Wang Y, Zhang H, Wang H, Hu Y, Wen Z, Deng H, et al. Development of a neoadjuvant chemotherapy efficacy prediction model for nasopharyngeal carcinoma integrating magnetic resonance radiomics and pathomics: a multi-center retrospective study. BMC Cancer. (2024) 24:1501. doi: 10.1186/s12885-024-13235-0
34. Yin J, Xu L, Wang S, Zhang L, Zhang Y, Zhai Z, et al. Integrating immune multi-omics and machine learning to improve prognosis, immune landscape, and sensitivity to first- and second-line treatments for head and neck squamous cell carcinoma. Sci Rep. (2024) 14:31454. doi: 10.1038/s41598-024-83184-y
35. Cavalieri S, De Cecco L, Brakenhoff RH, Serafini MS, Canevari S, Rossi S, et al. Development of a multiomics database for personalized prognostic forecasting in head and neck cancer: The Big Data to Decide EU Project. Head Neck. (2021) 43:601–12. doi: 10.1002/hed.26515
36. Tomita H, Kobayashi T, Takaya E, Mishiro S, Hirahara D, Fujikawa A, et al. Deep learning approach of diffusion-weighted imaging as an outcome predictor in laryngeal and hypopharyngeal cancer patients with radiotherapy-related curative treatment: a preliminary study. Eur Radiol. (2022) 32:5353–61. doi: 10.1007/s00330-022-08630-9
37. Lin M, Lin N, Yu S, Sha Y, Zeng Y, Liu A, et al. Automated prediction of early recurrence in advanced sinonasal squamous cell carcinoma with deep learning and multi-parametric MRI-based radiomics nomogram. Acad Radiol. (2023) 30:2201–11. doi: 10.1016/j.acra.2022.11.013
38. Huynh BN, Groendahl AR, Tomic O, Liland KH, Knudtsen IS, Hoebers F, et al. Head and neck cancer treatment outcome prediction: a comparison between machine learning with conventional radiomics features and deep learning radiomics. Front Med (Lausanne). (2023) 10:1217037. doi: 10.3389/fmed.2023.1217037
39. Li W, Li Y, Wang L, Yang M, Iikubo M, Huang N, et al. Evaluating fusion models for predicting occult lymph node metastasis in tongue squamous cell carcinoma. Eur Radiol. (2025) 35(9):5228–38. doi: 10.1007/s00330-025-11473-9
40. Wang W, Liang H, Zhang Z, Xu C, Wei D, Li W, et al. Comparing three-dimensional and two-dimensional deep-learning, radiomics, and fusion models for predicting occult lymph node metastasis in laryngeal squamous cell carcinoma based on CT imaging: a multicentre, retrospective, diagnostic study. EClinicalMedicine. (2024) 67:102385. doi: 10.1016/j.eclinm.2023.102385
41. Choi JH, Choi JY, Woo SK, Moon JE, Lim CH, Park SB, et al. Prognostic value of radiomic analysis using pre- and post-treatment (18)F-FDG-PET/CT in patients with laryngeal cancer and hypopharyngeal cancer. J Pers Med. (2024) 14:71. doi: 10.3390/jpm14010071
42. Lin CH, Yan JL, Yap WK, Kang CJ, Chang YC, Tsai TY, et al. Prognostic value of interim CT-based peritumoral and intratumoral radiomics in laryngeal and hypopharyngeal cancer patients undergoing definitive radiotherapy. Radiother Oncol. (2023) 189:109938. doi: 10.1016/j.radonc.2023.109938
43. Agarwal JP, Sinha S, Goda JS, Joshi K, Mhatre R, Kannan S, et al. Tumor radiomic features complement clinico-radiological factors in predicting long-term local control and laryngectomy free survival in locally advanced laryngo-pharyngeal cancers. Br J Radiol. (2020) 93:20190857. doi: 10.1259/bjr.20190857
44. Al-Ibraheem A, Abdlkadir AS, Al-Adhami D, Hejleh TA, Mansour A, Mohamad I, et al. The prognostic and diagnostic value of [(18)F]FDG PET/CT in untreated laryngeal carcinoma. J Clin Med. (2023) 12:3514. doi: 10.3390/jcm12103514
45. Nakajo M, Nagano H, Jinguji M, Kamimura Y, Masuda K, Takumi K, et al. The usefulness of machine-learning-based evaluation of clinical and pretreatment 18F-FDG-PET/CT radiomic features for predicting prognosis in patients with laryngeal cancer. Br J Radiol. (2023) 96:20220772. doi: 10.1259/bjr.20220772
46. Ganeshan B and Miles KA. Quantifying tumour heterogeneity with CT. Cancer Imaging. (2013) 13:140–9. doi: 10.1102/1470-7330.2013.0015
47. Zhang H, Graham CM, Elci O, Griswold ME, Zhang X, Khan MA, et al. Locally advanced squamous cell carcinoma of the head and neck: CT texture and histogram analysis allow independent prediction of overall survival in patients treated with induction chemotherapy. Radiology. (2013) 269:801–9. doi: 10.1148/radiol.13130110
48. Mohsen F, Ali H, El Hajj N, and Shah Z. Artificial intelligence-based methods for fusion of electronic health records and imaging data. Sci Rep. (2022) 12:17981. doi: 10.1038/s41598-022-22514-4
49. Wang R, Dai W, Gong J, Huang M, Hu T, Li H, et al. Development of a novel combined nomogram model integrating deep learning-pathomics, radiomics and immunoscore to predict postoperative outcome of colorectal cancer lung metastasis patients. J Hematol Oncol. (2022) 15:11. doi: 10.1186/s13045-022-01225-3
Keywords: laryngeal squamous cell carcinoma, locally advanced, multimodal features, recurrence-free survival, decision-level fusion
Citation: Zhao F, Huang X, Li J, He J, Liu J, Chen G and Zhang Z (2025) Development and validation of a multimodal feature fusion-based model for predicting postoperative recurrence-free survival in locally advanced laryngeal squamous cell carcinoma. Front. Oncol. 15:1685737. doi: 10.3389/fonc.2025.1685737
Received: 14 August 2025; Accepted: 12 September 2025;
Published: 25 September 2025.
Edited by:
Mihai Dumitru, Carol Davila University of Medicine and Pharmacy, RomaniaReviewed by:
Claire Massagee Lanier, Wake Forest Baptist Medical Center, United StatesBogdan Cobzeanu, Grigore T. Popa University of Medicine and Pharmacy, Romania
Copyright © 2025 Zhao, Huang, Li, He, Liu, Chen and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhe Zhang, emhhbmd6aGVAZ3htdS5lZHUuY24=
†These authors have contributed equally to this work and share first authorship