Abstract
Objectives:
To investigate the potential of intratumoral and peritumoral radiomics derived from CT to preoperatively predict tumor deposits (TDs) in patients with advanced gastric cancer (AGC).
Methods:
In this retrospective investigation, a total of 374 patients from two medical centers were recruited and divided into training (n = 186), validation (n = 80), and test (n = 108) cohorts. Intratumoral and peritumoral radiomics models were developed utilizing radiomics features derived from the corresponding 3D regions of interest (ROIs). A combined radiomics model integrating intratumoral and peritumoral features was further constructed through feature-level concatenation. Additionally, an ensemble model was established via the integration of this combined radiomics model with selected independent clinical prognostic factors. All models were evaluated using the area under the receiver operating characteristic curve (AUC), calibration curves, and decision curve analysis (DCA). Finally, the Shapley Additive Explanations (SHAP) method and nomogram were employed to elucidate the predictive mechanisms of the three radiomics models (intratumoral, peritumoral, and combined) and the ensemble model.
Results:
The combined intratumoral-peritumoral radiomics model showed higher AUC than the standalone intratumoral and peritumoral models across all cohorts (training: 0.874 vs. 0.751 vs. 0.830; validation: 0.846 vs. 0.720 vs. 0.713; test: 0.842 vs. 0.701 vs. 0.675). Moreover, the ensemble model yielded the highest AUCs across all cohorts (0.925, 0.865, 0.878 for training, validation, and test cohorts, respectively).
Conclusion:
Both intratumoral and peritumoral radiomics offer meaningful information regarding TDs, while the CT-based ensemble model holds the capacity to preoperatively predict TDs in AGC patients.
Introduction
Gastric cancer (GC) remains a significant public health challenge globally. It is the fifth most common malignancy and the fourth leading cause of cancer-related death worldwide (). We usually choose the appropriate treatment and assess prognosis based on TNM staging system. The TNM staging system, established by the American Joint Committee on Cancer (AJCC) and the International Union Against Cancer (UICC), is the most widely used tool for GC staging. However, an increasing number of studies have revealed the limitations of the TNM staging system (–). Patients with the same stage receiving identical treatments show significant variations in prognosis. Recurrence occurs even in patients who undergo radical resection for gastric cancer without lymph node metastasis. This indicates that in addition to tumor depth, gastric regional lymph node involvement, and distant metastasis, other factors also influence treatment selection and prognosis, such as tumor differentiation, lymphovascular invasion (LVI), perineural invasion (PNI), HER2 expression, MSI status, and tumor deposits (TDs) (–).
As a typical histopathological feature of colorectal tumors, TDs were first identified by Gabriel W.B. in 1935 (). With the advancement in TDs research, their clinical definition has been increasingly standardized. The AJCC 8th edition TNM staging system defines colorectal cancer TDs as discrete tumor nodules within the primary tumor’s lymphatic drainage area that lack identifiable lymph node tissue, vascular structures, or neural architecture (). Their presence indicates a higher risk of tumor dissemination and recurrence, leading to a poorer prognosis (–). Studies have shown that GC patients with TDs have a significantly lower 5-year survival rate than those without TDs (, ). However, TDs’ diagnosis mainly depends on postoperative pathology. Currently, there’s a lack of preoperative and noninvasive diagnostic tools for TDs. Preoperative diagnosis of TDs may assist clinicians in refining preoperative staging and devising perioperative treatment plans, with the potential to enhance patients’ prognoses.
Preoperative prediction of tumor deposits (TDs) in gastric cancer remains a significant clinical challenge. Previous studies have focused on assessment of TDs status using clinical characteristics: Anup et al. ()indicated that TDs are more prevalent in Borrmann IV gastric cancer and are positively associated with deeper tumor invasion and an increased likelihood of lymph node metastasis. Li et al. ()further reported that TDs are closely correlated with tumor diameter, Borrmann classification, tumor differentiation, pathological T stage (pT), pathological N stage (pN), pathological TNM stage (pTNM), perineural invasion, and vascular invasion. However, these studies rely on postoperative pathological features for the prediction of TDs, with limited practical utility for preoperative accurate prediction of TDs.
Notably, TDs are typically small in volume, and conventional computed tomography (CT) or magnetic resonance imaging (MRI) is constrained by scan slice thickness and spatial resolution, often failing to delineate the morphological characteristics of TDs adequately. Radiomics, a core technology in precision medicine, involves extracting high-throughput quantitative features from medical images and integrating them with machine learning algorithms to facilitate disease diagnosis, staging, treatment response prediction, and prognosis assessment (, ). Its application in gastric cancer has emerged as a research hotspot in recent years (, ).
Most prior radiomics studies have focused on intratumoral features derived from the tumor parenchyma. In contrast, peritumoral radiomics targets high-throughput quantitative characteristics of peritumoral tissues (including adjacent parenchyma, stroma, blood vessels, and fat spaces) (). These features indirectly reflect the biological properties of the tumor microenvironment (TME), such as inflammatory responses, angiogenesis, fibrosis, and invasive potential, thereby better capturing the tumor-host interaction (). As such, peritumoral radiomics serves as a key “imaging surrogate marker” for evaluating tumor invasiveness, metastatic potential, prognosis, and treatment response (–). However, the utility of both intratumoral and peritumoral radiomic features derived from CT in predicting TDs in patients with gastric cancer remains incompletely elucidated.
Therefore, the present study aimed to investigate the predictive value of CT-based intratumoral and peritumoral radiomics for TDs in AGC patients and to develop robust, interpretable models with potential clinical applicability.
Materials and methods
Patients and clinical characteristics
The study protocol was reviewed and approved by the Institutional Review Board of our hospital. Written informed consent was waived, given the retrospective observational design of the research. Pathologically confirmed gastric cancer patients were consecutively enrolled from two medical centers. Inclusion criteria included (1): pathologically confirmed advanced gastric cancer (AGC) and undergoing gastrectomy combined with lymph node dissection (2); complete clinical and postoperative pathological data (3); availability of non-enhanced and dual-phase enhanced computed tomography (CT) scans performed within 1 week prior to treatment. Exclusion criteria included (1): inadequate gastric distension (2); poor visualization of lesions on CT images (3); tumor diameter < 5 mm (rendering the region of interest insufficient) (4); poor image quality. Between December 2015 and March 2020, a total of 913 patients were enrolled from Center I, of whom 859 were eligible. This eligible cohort consisted of 133 TDs-positive and 726 TDs-negative patients. To mitigate potential data imbalance, 133 patients were randomly selected from the TDs-negative group and matched 1:1 with the TDs-positive group. In addition, in accordance with the inclusion and exclusion criteria, a total of 108 consecutive patients were enrolled from Center II between March 2019 and March 2020 as the test cohort, of whom 28 were TDs-positive. Patients from Center I were further randomly divided into a training cohort(n=186) and a validation cohort(n=80) at a ratio of 7:3. The detailed patient inclusion and exclusion process is illustrated in Figure 1.
Figure 1
Clinical data collected in the present study comprised sex, age, tumor location, tumor thickness, maximum tumor diameter, carcinoembryonic antigen (CEA), carbohydrate antigen 12-5 (CA12-5), carbohydrate antigen 19-9 (CA19-9), and carbohydrate antigen 72-4 (CA72-4). Tumor location was categorized as cardia, body, antrum, or involvement of ≥2/3 of the gastric wall. Tumor thickness was defined as the maximal perpendicular distance from the lesion surface to the deepest infiltration point, measured in the optimal plane (axial, sagittal, or coronal). Maximum tumor diameter was additionally measured at the lesion’s maximal cross-section. CEA, CA12-5, CA19-9, and CA72–4 were dichotomized into two categories (normal vs. abnormal) based on respective reference values (CEA ≤ 3.4 μg/L, CA12-5 ≤ 35 μg/L, CA19-9 ≤ 27 μg/L, CA72-4 ≤ 6.7 μg/L).
CT image acquisition
The workflow of this study is illustrated in Figure 2. All patients underwent scanning using a 64-slice spiral CT scanner (Discovery CT 750 HD; GE Healthcare, Waukesha, WI, USA; or Siemens Sensation 64 CT; Siemens Healthcare, Forchheim, Germany). Patients were required to fast for at least 6 hours before scanning and received oral administration of 500–1000 ml of water to distend the gastric cavity prior to the examination. All patients were placed in the supine position, and the scanning range for both plain and contrast-enhanced CT scans covered from the top of the diaphragm to the upper edge of the pubic symphysis. Plain CT scanning was performed first. For contrast-enhanced scanning, a high-pressure injector (Urich REF XD 2060-Touch, Ulrich Medical) was used to inject the contrast agent iohexol (Shanghai Bracco Sine Pharmaceutical Co., Ltd., iodine concentration 370 mg/mL) via the antecubital vein at an injection flow rate of 3.5 mL/s and a dosage of 1.5 mL/kg. Arterial phase scanning was initiated 25–30 seconds after contrast agent injection, and venous phase scanning was started at 60–70 seconds. The scanning parameters were as follows: tube voltage of 120 kV, rotation time of 0.5−0.6 s, tube current of 250−640 mA, reconstruction section thickness of 1.25 mm and 5 mm.
Figure 2
Image preprocessing and segmentation
Venous-phase images were extracted from the picture archiving and communication system (PACS) in Digital Imaging and Communications in Medicine (DICOM) format. To mitigate variability introduced by differences in scanning protocols or equipment on quantitative radiomics features, images were resampled to a voxel size of 1×1×1 mm³. Intratumoral three-dimensional (3D) regions of interest (ROIs) for gastric cancer were manually delineated slice-by-slice on axial images along the tumor boundary by YY using ITK-SNAP software (version 3.8.0; http://www.itksnap.org), with the delineation process supervised by a radiologist with 5 years of clinical experience. Following completion of intratumoral ROI outlining, a 2-mm peripheral ring was automatically generated using the SimpleITK package in Python software (version 3.6). Notably, regions outside the peritumoral area were manually excluded. To ensure intra-observer reliability, the ROIs of 30 randomly selected patients were re-delineated by YY one month later. For inter-observer reliability assessment, ZP independently performed ROI segmentation and feature extraction on the medical images of the same 30 patients. The intra-observer and inter-observer correlation coefficients (ICCs) were used to evaluate the consistency of radiomics features, with ICCs value > 0.80 indicating high consistency.
Radiomics feature extraction
Following delineation of all ROIs, the PyRadiomics package (Version 3.6) was employed to extract intratumoral, peritumoral, and combined intratumoral-peritumoral radiomics features. Specifically, image preprocessing and transformation approaches included the original image (without additional processing), wavelet transform, Laplacian of Gaussian (LoG) transform, 2D local binary pattern (LBP2D), 3D local binary pattern (LBP3D), square transformation, square root transformation, logarithmic transformation, exponential transformation, and gradient transformation. Radiomics features encompassed first-order features, shape-based features, gray-level co-occurrence matrix (GLCM), gray-level size zone matrix (GLSZM), gray-level run length matrix (GLRLM), gray-level dependence matrix (GLDM), and neighborhood gray-tone difference matrix (NGTDM). All selected features are well-established classic radiomics descriptors.
Feature selection criteria and dimension reduction
First, the consistency of the extracted radiomics features was assessed, and only those with ICCs ≥ 0.80 were retained. All eligible features were then normalized via Z-score normalization. A three-step sequential procedure was adopted for optimal feature selection: initially, the Mann–Whitney U test was employed to screen features with a P-value < 0.05; next, Spearman’s rank correlation analysis was performed to exclude features with a correlation coefficient > 0.9; subsequently, tenfold cross-validation was combined with the least absolute shrinkage and selection operator (LASSO) regression to identify optimal features with non-zero coefficients, which was determined by the optimal penalty parameter (λ).
Clinical, radiomics and ensemble model development
Univariable logistic regression analysis was conducted for clinical characteristics in the training cohort, including sex, age, tumor location, tumor thickness, maximum tumor diameter, CEA, CA12-5, CA19-9, and CA72-4. Statistically significant clinical characteristics (P-value < 0.05) were subsequently included in a multivariable logistic regression model to identify independent clinical predictors, which were used for the construction of the clinical model.
Using the selected radiomic features, logistic regression was used to construct intratumoral, peritumoral, and combined intratumoral-peritumoral models. The predicted probabilities derived from the combined intratumoral-peritumoral model, along with the previously identified independent clinical predictors, were integrated into a multivariable logistic regression model to establish the ensemble model. This ensemble model was visualized using a nomogram.
Statistical analysis
Statistical analyses and model construction/evaluation were performed using Python (version 3.6). Continuous variables with a normal distribution were presented as mean ± standard deviation (SD), while those with a non-normal distribution were expressed as median (interquartile range). For continuous variables, Student’s t-test or Mann–Whitney U test was employed based on their distribution characteristics. For categorical variables, the chi-square test or Fisher’s exact test was used, as appropriate.
The predictive performance of the models was evaluated using the area under the curve (AUC), accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). Calibration curves were constructed to visually compare predicted probabilities with observed outcomes. Decision curve analysis (DCA) was conducted to quantify the clinical net benefit of each model. Differences in AUCs between models were compared using the DeLong test. The SHapley Additive exPlanations (SHAP) method was applied to interpret the radiomics models by quantifying the contribution of each feature to the model’s predictions. A two-tailed P < 0.05 was considered statistically significant.
Results
Clinical characteristics
A total of 374 patients with AGC were enrolled in the present study, including 161 TDs-positive patients and 213 TDs-negative patients. All patients were further divided into a training cohort (n=186), a validation cohort (n=80), and a test cohort (n=108). Detailed demographic and clinical characteristics of the patients are summarized in Table 1.
Table 1
| Characteristics | Total, (n=374) | Training, (n=186) | Validation, (n = 80) | Test, (n = 108) | p |
|---|---|---|---|---|---|
| Age(years) | 60.38 ± 10.45 | 61.42 ± 10.16 | 59.58 ± 10.56 | 59.18 ± 10.78 | 0.154 |
| Longitude(mm) | 64.35 ± 23.94 | 63.74 ± 20.37 | 62.73 ± 23.99 | 66.58 ± 29.12 | 0.492 |
| Thickness(mm) | 20.94 ± 8.93 | 21.24 ± 8.91 | 20.98 ± 9.50 | 20.40 ± 8.59 | 0.741 |
| Sex | 0.523 | ||||
| Male | 286(76.47) | 140(75.27) | 65(81.25) | 81(75.00) | |
| Female | 88(23.53) | 46(24.73) | 15(18.75) | 27(25.00) | |
| CEA | 0.071 | ||||
| Normal | 270(72.19) | 128(68.82) | 55(68.75) | 87(80.56) | |
| Abnormal | 104(27.81) | 58(31.18) | 25(31.25) | 21(19.44) | |
| CA19-9 | 0.131 | ||||
| Normal | 294(78.61) | 140(75.27) | 62(77.50) | 92(85.19) | |
| Abnormal | 80(21.39) | 46(24.73) | 18(22.50) | 16(14.81) | |
| CA72-4 | 0.054 | ||||
| Normal | 272(72.73) | 135(72.58) | 51(63.75) | 86(79.63) | |
| Abnormal | 102(27.27) | 51(27.42) | 29(36.25) | 22(20.37) | |
| CA12-5 | 0.867 | ||||
| Normal | 333(89.04) | 164(88.17) | 72(90.00) | 97(89.81) | |
| Abnormal | 41(10.96) | 22(11.83) | 8(10.00) | 11(10.19) | |
| Location | 0.239 | ||||
| Cardia | 186(49.73) | 91(48.92) | 47(58.75) | 48(44.44) | |
| Body | 33(8.82) | 14(7.53) | 7(8.75) | 12(11.11) | |
| Antrum | 131(35.03) | 69(37.10) | 19(23.75) | 43(39.81) | |
| Whole | 24(6.42) | 12(6.45) | 7(8.75) | 5(4.63) |
Baseline demographic and clinical characteristics of patients.
CEA, carcinoembryonic antigen; CA19–9, carbohydrate antigen 19-9; CA72–4, carbohydrate antigen 72-4; CA12–5, carbohydrate antigen 12-5.
Clinical model
In the training set, univariate logistic regression analysis revealed that CEA, CA12-5, CA19-9, and CA72–4 exhibited statistically significant disparities between patients with positive TDs and those with negative TDs. Subsequently, multivariate logistic regression analysis identified CEA and CA19–9 as independent predictive factors for TDs status (Table 2). Clinical prediction models were further constructed based on these two serological biomarkers (CEA and CA19-9).
Table 2
| Variable | Univariable analysis | Multivariable analysis | ||
|---|---|---|---|---|
| Log (OR) | p | Log (OR) | p | |
| Gender | 0.928(0.773,1.115) | 0.505 | ||
| Location | 0.955(0.859,1.063) | 0.478 | ||
| Thickness | 0.997(0.987,1.008) | 0.682 | ||
| Age | 0.999(0.995,1.003) | 0.59 | ||
| Longitude | 0.999(0.995,1.002) | 0.608 | ||
| CA12-5 | 2.667(1.213,5.859) | 0.04 | 1.206(0.501,2.904) | 0.726 |
| CA72-4 | 2.923(1.723,4.958) | 0.001 | 1.124(0.577,2.190) | 0.773 |
| CEA | 4.800(2.710,8.499) | <0.001 | 2.770(1.358,5.652) | 0.019 |
| CA19-9 | 5.572(2.835,10.946) | <0.001 | 2.772(1.264,6.08) | 0.033 |
Univariable and multivariable analyses of the clinical characteristics in the training cohort (n=186).
CEA, carcinoembryonic antigen; CA19–9, carbohydrate antigen 19-9; CA72–4, carbohydrate antigen 72-4; CA12–5, carbohydrate antigen 12-5.
Feature selection and model construction
A total of 1,834 radiomic features were extracted from intratumoral and peritumoral regions based on both original and processed images. For the combined intratumoral-peritumoral analysis, 3,668 radiomic features were extracted. Following consistency assessment using the ICCs, 1,178, 1,356, and 2,309 features were retained for the intratumoral, peritumoral, and combined intratumoral-peritumoral models, respectively (ICCs ≥ 0.80). Subsequent to the three-step feature selection process, 12, 16, and 27 optimal radiomic features were identified for inclusion in the intratumoral, peritumoral, and combined intratumoral-peritumoral models, respectively. Detailed information regarding these selected radiomic features is provided in Supplementary Table 1. These final feature sets were then employed to construct the respective intratumoral, peritumoral, and combined intratumoral-peritumoral models using logistic regression. Finally, the ensemble model was established by fusing the predictive outputs of the combined intratumoral-peritumoral model with the previously identified independent clinical predictors.
Evaluation of model performance
A total of five models were constructed in this study. The area under the curve (AUC), accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) of each model for predicting TDs status in AGC are presented in Table 3. Receiver operating characteristic (ROC) curves for all models are illustrated in Figure 3. Pairwise comparisons of model performance using the DeLong test are summarized in Table 4.
Table 3
| Cohort | Models | Accuracy | Sensitivity | Specificity | PPV | NPV | AUC (95% CI) |
|---|---|---|---|---|---|---|---|
| Training | Clinic | 0.742 | 0.640 | 0.835 | 0.781 | 0.717 | 0.764(0.703,0.825) |
| Intratumoral | 0.704 | 0.787 | 0.629 | 0.66 | 0.762 | 0.751(0.682,0.820) | |
| Peritumoral | 0.753 | 0.820 | 0.691 | 0.709 | 0.807 | 0.830(0.772,0.889) | |
| Combined | 0.806 | 0.820 | 0.794 | 0.785 | 0.828 | 0.874(0.825,0.924) | |
| Ensemble | 0.86 | 0.888 | 0.835 | 0.832 | 0.89 | 0.925(0.888,0.961) | |
| Validation | Clinic | 0.725 | 0.636 | 0.833 | 0.824 | 0.652 | 0.735(0.636,0.834) |
| Intratumoral | 0.713 | 0.705 | 0.722 | 0.756 | 0.667 | 0.720(0.599,0.828) | |
| Peritumoral | 0.700 | 0.727 | 0.667 | 0.727 | 0.667 | 0.713(0.611,0.838) | |
| Combined | 0.825 | 0.818 | 0.833 | 0.857 | 0.789 | 0.846(0.753,0.939) | |
| Ensemble | 0.813 | 0.773 | 0.861 | 0.872 | 0.756 | 0.865(0.782,0.948) | |
| Test | Clinic | 0.741 | 0.536 | 0.812 | 0.500 | 0.833 | 0.694(0.585,0.802) |
| Intratumoral | 0.704 | 0.679 | 0.712 | 0.452 | 0.864 | 0.701(0.589,0.813) | |
| Peritumoral | 0.676 | 0.643 | 0.687 | 0.419 | 0.846 | 0.675(0.559,0.791) | |
| Combined | 0.806 | 0.714 | 0.837 | 0.606 | 0.893 | 0.842(0.760,0.925) | |
| Ensemble | 0.833 | 0.786 | 0.850 | 0.647 | 0.919 | 0.878(0.807,0.949) |
Predictive efficacy of different models for TDs status in advanced gastric cancer across the training, validation, and test cohorts.
Figure 3
Table 4
| Comparison Model | Clinic | Intratumoral | Peritumoral | Combined |
|---|---|---|---|---|
| Training | ||||
| Intratumoral | 0.786 | |||
| Peritumoral | 0.122 | 0.043 | ||
| Combined | 0.007 | <0.001 | 0.010 | |
| Ensemble | <0.001 | <0.001 | <0.001 | 0.004 |
| Validation | ||||
| Intratumoral | 0.842 | |||
| Peritumoral | 0.780 | 0.920 | ||
| Combined | 0.079 | 0.018 | 0.012 | |
| Ensemble | 0.005 | 0.004 | 0.009 | 0.373 |
| Test | ||||
| Intratumoral | 0.921 | |||
| Peritumoral | 0.820 | 0.758 | ||
| Combined | 0.031 | 0.045 | 0.034 | |
| Ensemble | <0.001 | 0.008 | 0.006 | 0.185 |
Statistical p-values from DeLong tests for pairwise comparisons of all models across the training, validation, and test cohorts.
For the training cohort, the ensemble model exhibited an AUC of 0.925 (95% confidence interval [CI] 0.888–0.961), which was statistically significantly superior to the clinical model (p<0.001), intratumoral model (p<0.001), peritumoral model (p<0.001), and combined intratumoral-peritumoral model (p=0.004). In the validation cohort, the ensemble model achieved an AUC of 0.865 (95%CI 0.782–0.948), demonstrating significant superiority over the clinical model (p=0.005), intratumoral model (p=0.004), and peritumoral model (p=0.009). Regarding the test cohort, the ensemble model yielded an AUC of 0.878 (95%CI 0.807–0.949), which was significantly higher than the clinical model (p<0.001), intratumoral model (p=0.008), and peritumoral model (p=0.006). Additionally, the combined intratumoral-peritumoral model outperformed both the intratumoral and peritumoral models across all cohorts, with statistical significance.
Decision curve analysis (DCA) revealed that the ensemble model yielded a higher net clinical benefit compared with the other models across various threshold probabilities (Figure 4). Calibration curve analysis demonstrated that the ensemble model exhibited superior predictive accuracy relative to both the individual radiomics models and the clinical model (Figure 5). SHAP (SHapley Additive exPlanations) analysis results for the constructed radiomics models are depicted in Figure 6. An illustrative example demonstrating the interpretability of the proposed radiomics model is presented in Figure 7.
Figure 4
Figure 5
Figure 6
Figure 7
Discussion
In this study, we constructed and validated a series of radiomics models to investigate the potential of CT-derived radiomics features in providing predictive value for the preoperative assessment of TDs status in patients with AGC. Notably, the ensemble model—integrating the output of the combined intratumoral-peritumoral radiomics model with independent clinical prognostic factors—achieved the most robust performance. The reliability and generalizability of this ensemble model were fully corroborated by its consistently favorable performance across the training, validation and test cohorts.
Traditionally, the prognosis assessment of gastric cancer patients has primarily centered on pT stage, pN stage, and the absence of distant metastasis (, ). However, numerous studies have highlighted the inadequacy of the 8th edition AJCC TNM staging system in accurately stratifying outcomes for TDs-positive patients. For instance, Sun et al. () demonstrated that when stratified by pN category, significant differences in survival were observed between patients with and without TDs for those in pN0/pT1-3, pN1/pT3, pN2/pT1–3 and pN3/pT2–3 category. Mounting evidence has further confirmed that TDs represent an independent prognostic risk factor for GC: Zhou et al. () reported a 3-year disease-free survival (DFS) rate of 61.5% in the TDs-negative group, significantly higher than the 40.8% observed in the TDs-positive group (HR 1.78, 95% CI 1.44–2.21; P < 0.001). Beyond prognostic stratification, preoperative identification of TDs status also plays a pivotal role in optimizing surgical strategies and enhancing curative efficacy. Gong et al. () noted that TDs are frequently localized in the mesogastric region, a site that cannot be fully resected via conventional radical gastrectomy with D2 lymph node dissection. Consequently, D2 lymphadenectomy combined with complete mesogastric excision (D2+CME) is recommended for TDs-positive GC patients to achieve adequate tumor clearance. Based on the aforementioned findings, the precise preoperative diagnosis of TDs in GC patients holds profound clinical significance. Its value permeates key aspects including clinical decision-making optimization, individualized treatment formulation, accurate prognostic assessment, and improved reliability of tumor TNM staging—all of which exert a crucial impact on enhancing treatment outcomes and prolonging long-term survival for patients. Therefore, the development of an efficient and non-invasive method for preoperative TDs detection in GC patients addresses an urgent clinical unmet need and carries substantial practical value.
Radiomics has enabled the shift from “qualitative description” to “quantitative analysis” by quantifying latent features in gastric cancer imaging, demonstrating tremendous application value in preoperative diagnosis, precise staging, treatment response prediction, prognostic evaluation, and other aspects (, ). However, early studies mainly focused on intratumoral radiomics features (, , ). Notably, tumor invasion and metastasis are not isolated events but depend on the “remodeling” of the peritumoral microenvironment. Peritumoral radiomics features essentially represent quantitative imaging surrogates of pathophysiological changes in the tumor microenvironment (, ). Their biological significance lies in overcoming the limitation that intratumoral features only reflect the intrinsic characteristics of tumor cells, thereby more comprehensively revealing tumor invasiveness and the host’s antitumor response. Recent studies showed that integrating intratumoral and peritumoral radiomics features improves the performance of models in predicting tumor metastasis in colorectal cancer (CRC) and cervical cancer, with AUC values reaching approximately 0.80 in independent test cohorts (, ). Therefore, in our study, we incorporated both intratumoral and peritumoral radiomics features, and the predictive performance of the peritumoral radiomics model was not inferior to that of the intratumoral radiomics model. When the intratumoral and peritumoral radiomics features were combined, the predictive ability of the integrated model further improved, with an AUC of 0.878 in the training cohort, 0.838 in the validation cohort, and 0.853 in the test cohort. Thus, the additional value of peritumoral radiomics features in predicting the status of tumor deposits in gastric cancer was confirmed. It is important to note that this study adopted a 3D region of interest (3D ROI) segmentation strategy. Its core advantage lies in the comprehensive coverage of the entire tumor volume information, enabling more accurate quantification of tumor heterogeneity features (). This approach effectively avoids the prevalent local sampling bias and subjective judgment bias associated with 2D ROI segmentation, thereby laying a solid methodological foundation for the development of high-performance radiomics prediction models.
Previous studies have shown that the presence of TDs is associated with multiple clinical and postoperative pathological features, which found that TDs are closely correlated with tumor diameter, Borrmann classification, tumor differentiation, pT, pN, pTNM, nerve invasion, and vascular invasion (, ). Given that the aim of our study was to develop a preoperative model for predicting TDs status in patients with gastric cancer, only clinical features were included in the analysis. Univariate and multivariate logistic regression analyses identified CEA and CA19–9 as independent predictive factors for TD status. Finally, an ensemble model was established by fusing the predictive outputs of the combined intratumoral-peritumoral radiomics model with the identified independent clinical predictors. Among all constructed models, the ensemble model achieved the highest area under the curve (AUC) values in the training (0.931, 95%CI 0.896-0.966), validation (0.880, 95%CI 0.807-0.983), and test cohorts (0.883, 95%CI 0.795-0.970).
Model interpretability is pivotal for the clinical translation of machine learning-based radiomics models. To address this critical requirement, the SHAP method was utilized to conduct interpretability analysis of the constructed radiomics models in the present study (). Specifically, this approach quantifies the direction and magnitude of each feature’s contribution to individual patient-specific predictive outcomes, thereby facilitating comprehensive comprehension of the model’s decision-making mechanism among both clinicians and patients. Furthermore, a nomogram was developed to visualize the ensemble model, enabling clinicians to intuitively interpret predictive results and quantify the relative importance of each input variable. Collectively, the integration of these two complementary interpretability strategies significantly enhances the clinical applicability, trustworthiness, and acceptability of our predictive models.
Our study also has certain limitations. First, the retrospective design of this study may introduce inherent biases. Therefore, we plan to conduct multicenter prospective studies in the future to validate our findings. Second, our study used a fixed 2-mm peritumoral region for analysis, a predefined boundary that fails to fully reflect the biological heterogeneity of gastric cancer invasion and the anatomical adjacency of surrounding abdominal organs to the stomach. Future studies should optimize peritumoral delimitation by integrating radiological and pathological gold standards—for example, stratifying the peritumoral region by radial distances (2 mm, 5 mm, 10 mm) from the tumor margin and validating the optimal boundary with postoperative gastric cancer pathological sections. This optimized delimitation can exclude confounding from adjacent normal tissues and enhance the clinical relevance of extracted peritumoral radiomic features (). Finally, we did not explore the correlation between the number of TDs and clinical or radiomic features. In future research, we intend to develop a preoperative predictive model for assessing the number of TDs, which is expected to provide more precise guidance for the individualized treatment of gastric cancer patients ().
In conclusion, both intratumoral and peritumoral radiomic features can yield valuable insights into the status of TDs in patients with AGC. The ensemble model, integrating these two radiomic feature sets with clinical characteristics, exhibits superior performance and holds significant promise for guiding precision medicine in the management of AGC patients.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Ethics Committee of the First Affiliated Hospital of Zhengzhou University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
YY: Methodology, Data curation, Conceptualization, Writing – review & editing, Writing – original draft. PZ: Investigation, Methodology, Writing – original draft. MY: Writing – original draft, Investigation. YC: Writing – original draft, Investigation. YF: Writing – original draft, Investigation. JG: Conceptualization, Supervision, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study has received funding by National Natural Science Foundation of China (No. 81971615).
Acknowledgments
We would like to express our thanks to all of the patients who participated in this study and each of the researchers and technicians who made this work possible.
Conflict of interest
The author(s) declared that this work 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) declared that generative AI was not 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.2026.1763646/full#supplementary-material
Abbreviations
TDs, Tumor deposits; AGC, Advanced gastric cancer; AUC, Under the receiver operating characteristic curve; LASSO, The Least Absolute Shrinkage and Selection Operator; ROC, Receiver operating characteristic; CEA, Carcinoembryonic antigen; CA19-9, Carbohydrate antigen 19-9; CA72-4, Carbohydrate antigen 72-4; CA12-5, Carbohydrate antigen 12-5; ROI, Regions of interest; SHAP, Shapley Additive Explanations.
References
1
MorganEArnoldMCamargoMCGiniAKunzmannATMatsudaTet al. The current and future incidence and mortality of gastric cancer in 185 countries, 2020-40: A population-based modelling study. EClinicalMedicine. (2022) 47:101404. doi: 10.1016/j.eclinm.2022.101404
2
FangCWangWDengJYSunZSeeruttunSRWangZNet al. Proposal and validation of a modified staging system to improve the prognosis predictive performance of the 8th AJCC/UICC pTNM staging system for gastric adenocarcinoma: a multicenter study with external validation. Cancer Commun (Lond). (2018) 38:67. doi: 10.1186/s40880-018-0337-5
3
LinJXDesiderioJLinJPWangWTuRHLiPet al. Multicenter validation study of the american joint commission on cancer (8th edition) for gastric cancer: proposal for a simplified and improved TNM staging system. J Cancer. (2020) 11:3483–91. doi: 10.7150/jca.36891
4
ChoiWHKimMJParkJHKangJGSeoSIKimHYet al. Lymphatic invasion might be considered as an upstaging factor in N0 and N1 gastric cancer. J Clin Med. (2020) 9:1275. doi: 10.3390/jcm9051275
5
ZhouMYangWZouWYangJZhouCZhangZet al. Prognostic significance of tumor deposits in radically resected gastric cancer: a retrospective study of a cohort of 1915 Chinese individuals. World J Surg Oncol. (2022) 20:304. doi: 10.1186/s12957-022-02773-1
6
FengFLiuJWangFZhengGWangQLiuSet al. Prognostic value of differentiation status in gastric cancer. BMC Cancer. (2018) 18:865. doi: 10.1186/s12885-018-4780-0
7
HwangJEHongJYKimJEShimHJBaeWKHwangECet al. Prognostic significance of the concomitant existence of lymphovascular and perineural invasion in locally advanced gastric cancer patients who underwent curative gastrectomy and adjuvant chemotherapy. Jpn J Clin Oncol. (2015) 45:541–6. doi: 10.1093/jjco/hyv031
8
BiliciASekerMUstaaliogluBBKefeliUYildirimEYavuzerDet al. Prognostic significance of perineural invasion in patients with gastric cancer who underwent curative resection. Ann Surg Oncol. (2010) 17:2037–44. doi: 10.1245/s10434-010-1027-y
9
GabrielWBDukesCBusseyHJR. Lymphatic spread in cancer of the rectum. Br J Surgery. (1935) 23:395–413. doi: 10.1002/bjs.1800239017
10
AminMBGreeneFLEdgeSBComptonCCGershenwaldJEBrooklandRKet al. The Eighth Edition AJCC Cancer Staging Manual: Continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin. (2017) 67:93–9. doi: 10.3322/caac.21388
11
ChenHTangZChenLLiHWangXLiuFet al. Evaluation of the impact of tumor deposits on prognosis in gastric cancer and a proposal for their incorporation into the AJCC staging system. Eur J Surg Oncol. (2018) 44:1990–6. doi: 10.1016/j.ejso.2018.10.062
12
AnupSLuJZhengCHLiPXieJWWangJBet al. Prognostic significance of perigastric tumor deposits in patients with primary gastric cancer. BMC Surg. (2017) 17:84. doi: 10.1186/s12893-017-0280-4
13
Graham MartínezCKnijnNVerheijMNagtegaalIDvan der PostRS. aNijmegen, the Netherlands: Radboud University (Radboud Universiteit) (2019) p. 809–16.
14
WenquanLYuhuaLJianxinCHongqingXKechengZJiyangLet al. Tumor deposit serves as a prognostic marker in gastric cancer: A propensity score-matched analysis comparing survival outcomes. Cancer Med. (2020) 9:3268–77. doi: 10.1002/cam4.2963
15
LiZLiXGuanSZhuGLinHWuHet al. Correlation analysis between tumor deposit and clinicopathologic characteristics and prognosis of gastric cancer: A multicenter retrospective study. Ann Surg Oncol. (2024) 31:5984–96. doi: 10.1245/s10434-024-15508-1
16
LambinPRios-VelazquezELeijenaarRCarvalhoSvan StiphoutRGGrantonPet al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer. (2012) 48:441–6. doi: 10.1016/j.ejca.2011.11.036
17
GilliesRJKinahanPEHricakH. Radiomics: images are more than pictures, they are data. Radiology. (2016) 278:563–77. doi: 10.1148/radiol.2015151169
18
WangWPengYFengXZhaoYSeeruttunSRZhangJet al. Development and validation of a computed tomography-based radiomics signature to predict response to neoadjuvant chemotherapy for locally advanced gastric cancer. JAMA Netw Open. (2021) 4:e2121143. doi: 10.1001/jamanetworkopen.2021.21143
19
HuangWXiongWTangLChenCYuanQZhangCet al. Non-invasive CT imaging biomarker to predict immunotherapy response in gastric cancer: a multicenter study. J Immunother Cancer. (2023) 11:e007807. doi: 10.1136/jitc-2023-007807
20
LinXLiuZZhouKLiYHuangGZhangHet al. Intratumoral and peritumoral PET/CT-based radiomics for non-invasively and dynamically predicting immunotherapy response in NSCLC. Br J Cancer. (2025) 132:558–68. doi: 10.1038/s41416-025-02948-z
21
LiuHWuJXuXWangHZhangCYinSet al. Peritumoral TIGIT(+)CD20(+) B cell infiltration indicates poor prognosis but favorable adjuvant chemotherapeutic response in gastric cancer. Int Immunopharmacol. (2022) 108:108735. doi: 10.1016/j.intimp.2022.108735
22
MaoNShiYLianCWangZZhangKXieHet al. Intratumoral and peritumoral radiomics for preoperative prediction of neoadjuvant chemotherapy effect in breast cancer based on contrast-enhanced spectral mammography. Eur Radiol. (2022) 32:3207–19. doi: 10.1007/s00330-021-08414-7
23
LinPXieWLiYZhangCWuHWanHet al. Intratumoral and peritumoral radiomics of MRIs predicts pathologic complete response to neoadjuvant chemoimmunotherapy in patients with head and neck squamous cell carcinoma. J Immunother Cancer. (2024) 12:e009616. doi: 10.1136/jitc-2024-009616
24
WuLCenCOuyangDZhangLLiXWuHet al. Interpretable machine learning model for predicting early recurrence of pancreatic cancer: integrating intratumoral and peritumoral radiomics with body composition. Int J Surg. (2025) 111:8198–211. doi: 10.1097/JS9.0000000000003078
25
BrouwerNPMNagtegaalID. Tumor deposits improve staging in colon cancer: what are the next steps? Ann Oncol. (2021) 32:1209–11. doi: 10.1016/j.annonc.2021.08.1751
26
LordACD’SouzaNShawARokanZMoranBAbulafiMet al. MRI-diagnosed tumor deposits and EMVI status have superior prognostic accuracy to current clinical TNM staging in rectal cancer. Ann Surg. (2022) 276:334–44. doi: 10.1097/SLA.0000000000004499
27
SunZWangZNXuYYZhuGLHuangBJXuYet al. Prognostic significance of tumor deposits in gastric cancer patients who underwent radical surgery. Surgery. (2012) 151:871–81. doi: 10.1016/j.surg.2011.12.027
28
XieDShenJLiuLCaoBWangYQinJet al. Complete mesogastric excision for locally advanced gastric cancer: short-term outcomes of a randomized clinical trial. Cell Rep Med. (2021) 2:100217. doi: 10.1016/j.xcrm.2021.100217
29
ChenYXuWLiYLLiuWSahBKWangLet al. CT-based radiomics showing generalization to predict tumor regression grade for advanced gastric cancer treated with neoadjuvant chemotherapy. Front Oncol. (2022) 12:758863. doi: 10.3389/fonc.2022.758863
30
SongRCuiYRenJZhangJYangZLiDet al. CT-based radiomics analysis in the prediction of response to neoadjuvant chemotherapy in locally advanced gastric cancer: A dual-center study. Radiother Oncol. (2022) 171:155–63. doi: 10.1016/j.radonc.2022.04.023
31
LiuCLiLChenXHuangCWangRLiuYet al. Intratumoral and peritumoral radiomics predict pathological response after neoadjuvant chemotherapy against advanced gastric cancer. Insights Imaging. (2024) 15:23. doi: 10.1186/s13244-023-01584-6
32
WangXWeiMChenYJiaJZhangYDaiYet al. Intratumoral and peritumoral MRI-based radiomics for predicting extrapelvic peritoneal metastasis in epithelial ovarian cancer. Insights Imaging. (2024) 15:281. doi: 10.1186/s13244-024-01855-w
33
ZhangZWanXLeiXWuYZhangJAiYet al. Intra- and peri-tumoral MRI radiomics features for preoperative lymph node metastasis prediction in early-stage cervical cancer. Insights Imaging. (2023) 14:65. doi: 10.1186/s13244-023-01405-w
34
LiHChenXLLiuHLuTLiZL. MRI-based multiregional radiomics for predicting lymph nodes status and prognosis in patients with resectable rectal cancer. Front Oncol. (2022) 12:1087882. doi: 10.3389/fonc.2022.1087882
35
YanRQinSXuJZhaoWXinPXingXet al. A comparison of 2D and 3D magnetic resonance imaging-based intratumoral and peritumoral radiomics models for the prognostic prediction of endometrial cancer: a pilot study. Cancer Imaging. (2024) 24:100. doi: 10.1186/s40644-024-00743-2
36
LundbergSMNairBVavilalaMSHoribeMEissesMJAdamsTet al. Explainable machine-learning predictions for the prevention of hypoxaemia during surgery. Nat BioMed Eng. (2018) 2:749–60. doi: 10.1038/s41551-018-0304-0
37
HuangZMoSWuHKongYLuoHLiGet al. Optimizing breast cancer diagnosis with photoacoustic imaging: An analysis of intratumoral and peritumoral radiomics. Photoacoustics. (2024) 38:100606. doi: 10.1016/j.pacs.2024.100606
38
XuRZhangYZhaoJChenKWangZ. Prognostic value of tumor deposits in lymph node-negative gastric cancer: A propensity score matching study. Eur J Surg Oncol. (2023) 49:845–52. doi: 10.1016/j.ejso.2022.12.004
Summary
Keywords
CT, gastric cancer, preoperative prediction, radiomics, tumor deposits
Citation
Yao Y, Zhan P, Yuan M, Chen Y, Fan Y and Gao J (2026) Preoperative prediction of tumor deposits in advanced gastric cancer using intratumoral and peritumoral CT radiomics: development and validation of an ensemble model. Front. Oncol. 16:1763646. doi: 10.3389/fonc.2026.1763646
Received
09 December 2025
Revised
26 January 2026
Accepted
18 February 2026
Published
04 March 2026
Volume
16 - 2026
Edited by
Arvind Mukundan, National Chung Cheng University, Taiwan
Reviewed by
Wujie Chen, University of Chinese Academy of Sciences, China
Xiaotang Yang, Shanxi Provincial Cancer Hospital, China
Updates
Copyright
© 2026 Yao, Zhan, Yuan, Chen, Fan and Gao.
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: Jianbo Gao, gjb372001@163.com
Disclaimer
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.