ORIGINAL RESEARCH article

Front. Immunol., 31 March 2025

Sec. Cancer Immunity and Immunotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1504387

Nomogram based on CT imaging and clinical data to predict the efficacy of PD-1 inhibitors combined with chemotherapy in advanced gastric cancer

  • 1. Department of Radiology, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, Shandong, China

  • 2. School of Radiology, Shandong First Medical University, Taian, Shandong, China

  • 3. Department of General Education, Shandong First Medical University, Jinan, Shandong, China

Abstract

Background:

PD-1 inhibitors, in combination with chemotherapy, have become the first-line treatment option for patients with advanced metastatic gastric cancer. However, some patients still do not benefit from this treatment, highlighting an urgent need for simple and reliable markers to predict the efficacy of immunotherapy.

Methods:

Immunotherapy efficacy was evaluated using RECIST 1.1 and categorized into complete remission (CR), partial remission (PR), stable disease (SD), and disease progression (PD). Patients with CR, PR, and SD were classified as non-PD responders, while PD patients were categorized as PD responders. Clinical characteristics and CT imaging features of gastric cancer patients from two centers, before receiving PD-1 inhibitor combination chemotherapy, were retrospectively analyzed. A univariate logistic regression analysis was performed for each variable, and separate models for clinical and imaging characteristics, as well as a nomogram, were developed. Area under the curve (AUC), accuracy, sensitivity, specificity, and decision curve analysis (DCA) were used to evaluate all models.

Results:

Data from 272 patients (non-PD responders = 206, PD responders = 66) from Center 1 were collected for this study. Data from 76 patients (non-PD responders = 54, PD responders = 22) from Center 2 were used as an external validation cohort to verify the robustness of the models. We developed a clinical model, an imaging features model, and a nomogram. The nomogram, combining clinical and imaging features, demonstrated superior performance with an AUC of 0.904 (95% CI: 0.862–0.947) in the training set and an AUC of 0.801 (95% CI: 0.683-0.918) in the validation set, with sensitivity, specificity, and accuracy of 0.889, 0.682, and 0.829, respectively. Calibration curves further confirmed the agreement between actual results and predictions.

Conclusions:

A nomogram combining clinical features and CT imaging features before treatment was developed, which can effectively and simply predict the efficacy response of advanced gastric cancer patients treated with PD-1 inhibitors combined with chemotherapy. This tool can aid in optimizing treatment strategies in clinical practice.

1 Introduction

Gastric cancer ranks among the top five cancers worldwide in terms of incidence and mortality (), and remains a significant threat to human health. Due to its atypical early symptoms and rapid progression, most patients are diagnosed at advanced stages (). Recent advancements in immunotherapy have transformed cancer treatment, particularly with the introduction of immune checkpoint inhibitors (ICIs), such as anti-programmed cell death-1 (PD-1) and anti-programmed cell death ligand-1 (PD-L1) inhibitors. These agents have made significant strides in treating various malignancies, including advanced gastric cancer (AGC) (). Anti-PD-1 inhibitors have been officially approved for clinical use in multiple cancers (), and in recent years, chemotherapy combined with immunotherapy has become a standard regimen for many patients with advanced gastric cancer in clinical practice.

Although PD-1 inhibitors combined with chemotherapy have demonstrated improved efficacy compared to chemotherapy alone in some patients with advanced gastric cancer (), not all patients benefit from this approach in clinical practice (). Some patients continue to experience disease progression, immune-related adverse events (), and potentially excessive medical costs. Therefore, there is an urgent need for simple and reliable predictors of the efficacy of PD-1 inhibitors combined with chemotherapy to identify patients who are likely to benefit.

Current research predominantly focuses on biomarkers such as PD-L1 (programmed death ligand 1) expression, tumor mutational burden (TMB), microsatellite instability-high (MSI-H), and mismatch repair (MMR) status to predict the efficacy of immunotherapy in patients with advanced gastric cancer (–). However, many of these biomarker tests require invasive and costly biopsy sampling. In contrast, clinical features such as laboratory tests and imaging findings are more readily accessible and non-invasive. Tumor heterogeneity can complicate the assessment of overall tumor characteristics from biopsy specimens, while imaging tests can provide a comprehensive view of tumor characteristics. Wang et al. () demonstrated that changes in certain CT imaging features were associated with the degree of pathological remission in gastric cancer patients undergoing chemotherapy. Additionally, previous studies have explored the relationship between peripheral blood markers or innovative markers derived from them and the prognosis of gastric cancer patients undergoing immunotherapy (, ). However, most studies have only proposed associations between predictive factors and clinical outcomes without synthesizing multiple factors into a simple and accurate prediction model.

To more conveniently and accurately assess the efficacy of PD-1 inhibitor combination chemotherapy for gastric cancer, we retrospectively collected clinical characteristics and CT imaging features before treatment from patients with advanced gastric cancer treated with anti-PD-1 inhibitor combination chemotherapy. Our goal is to develop simple and effective predictive models to help identify patients who are likely to benefit from this treatment.

2 Materials and methods

2.1 Patient inclusion and exclusion criteria

The study was approved by the Ethics Committee of Shandong Provincial Hospital (Ethics number: NSFC: NO.2022-402). Informed consent was waived due to the retrospective nature of the study.

We retrospectively reviewed patients with advanced gastric cancer who were treated with anti-PD-1 inhibitors combined with chemotherapy from June 2019 to June 2023 at Shandong Provincial Hospital of Shandong First Medical University (Center 1) and from June 2019 to June 2023 at Shandong Cancer Hospital (Center 2). The inclusion criteria were as follows: (1) Pathologically confirmed gastric cancer. (2) Receipt of at least four cycles of anti-PD-1 inhibitor combined with chemotherapy. (3) Presence of at least one evaluable tumor lesion meeting the Response Evaluation Criteria in Solid Tumors (RECIST 1.1). (4) Availability of complete CT imaging before and after anti-PD-1 inhibitor combination chemotherapy. Exclusion criteria included: (1) Prior or concurrent surgical treatment of the stomach during anti-PD-1 inhibitor combination chemotherapy. (2) Inadequate gastric distension preventing accurate measurement of the lesion. Figure 1 shows the participant selection flowchart.

Figure 1

2.2 Treatment regimens and outcome assessment

All patients received first-line standard chemotherapy in conjunction with PD-1 inhibitors. The PD-1 inhibitors used in the study included camrelizumab, sintilimab, nivolumab, tislelizumab, and pembrolizumab, administered every three weeks. The primary chemotherapy regimens included SOX (S-1 + oxaliplatin), XELOX (capecitabine + oxaliplatin), and FOLFOX (leucovorin, fluorouracil, and oxaliplatin). The specific chemotherapy regimen was chosen based on the patient’s clinical condition and preference. Tumor response was evaluated approximately four weeks after the completion of four cycles of PD-1 inhibitor combination chemotherapy. The assessment was conducted by radiologists using computed tomography (CT), following the Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1. The response outcomes were categorized as complete remission (CR), partial remission (PR), stable disease (SD), or disease progression (PD). Based on these outcomes, patients were grouped as either non-PD responders (CR, PR, or SD) or PD responders.

2.3 Data collection

In this study, clinical characteristics of patients within two weeks prior to their first PD-1 inhibitor combined chemotherapy were obtained through the clinical electronic medical record system. The clinical characteristics include the patient’s gender, age at onset, body mass index (weight[kg]/height [m]²), treatment regimen, presence of distant metastasis, and certain laboratory indicators.

These laboratory indices included alpha-fetoprotein, carcinoembryonic antigen, carbohydrate antigen 125 (CA 125), carbohydrate antigen 19-9 (CA 19-9), carbohydrate antigen 72-4 (CA 72-4), platelet count, and the absolute values of lymphocytes (ALC), monocytes (AMC), and neutrophils (ANC). The inflammation indices included the neutrophil-to-lymphocyte ratio (NLR), calculated as ANC/ALC (absolute neutrophil count [×109/L]/absolute lymphocyte count [×109/L]), the monocyte-to-lymphocyte ratio (MLR), calculated as AMC/ALC (absolute monocyte count [×109/L]/absolute lymphocyte count [×109/L]), and the systemic immune-inflammation index (SII), calculated as (peripheral platelet count [×109/L] × absolute neutrophil count [×109/L])/absolute lymphocyte count [×109/L].

In this study, CT images of patients taken within two weeks before their first PD-1 inhibitor combination chemotherapy and after four cycles of treatment were retrieved from the Picture Archiving and Communication System (PACS). CT imaging characterization was performed by two radiologists with more than three years of experience in abdominal diagnostics. They independently reviewed all CT images using ITK-SNAP (RRID: SCR_017341) (http://www.itksnap.org) (). The CT images were segmented into lesions using ITK-SNAP, and the lesion areas were outlined as regions of interest (ROI). Disagreements in the assessments were resolved by a senior physician with 30 years of diagnostic radiology experience to reach a consensus. CT signatures included the following characteristics (, ): (1) Tumor size (maximum axial diameter and thickness), primary location (based on endoscopic findings: cardia, gastric body, or gastric antrum), and whether the tumor invades the visceral peritoneum; (2) External features: whether the perigastric fat is infiltrated or not and whether the tumor is surrounded by enlarged lymph nodes (defined as regional lymph nodes with a short diameter of 1 cm or more); (3) CT value parameters: The largest cross-section of the tumor was measured in the axial plane, avoiding necrotic and vascular areas, and its CT values were recorded in the plain, arterial, venous, and delayed phases. The final result was the average of measurements from two radiologists. Tumor arterial attenuation (HUA = HUA - HUU), venous attenuation (HUV = HUV - HUU), and arterial enhancement fraction (AEF) = (HUA - HUU)/(HUV - HUU) × 100% were also calculated () (where HUA represents arterial phase attenuation; HUV represents venous phase attenuation; HUA represents CT values in the arterial phase; HUV represents CT values in the venous phase; and HUU represents CT value of the lesion in the unenhanced phase) ().

2.4 Statistical analysis and nomogram development

Several analytical scales were designed to evaluate clinical and imaging characteristics of the patients. Statistical data were analyzed and visualized using IBM SPSS Statistics (RRID: SCR_019096) (version 25.0; SPSS, Inc., Chicago, IL, USA) and GraphPad Prism (RRID: SCR_002798) (version 9.0; GraphPad, San Diego, CA) for statistical analysis and graphing. Continuous variables with a normal distribution were expressed as mean ± standard deviation, while those not following a normal distribution were expressed as median (interquartile range). Categorical variables were expressed as frequencies and percentages. All statistical tests were two-sided, and a p-value < 0.05 was considered statistically significant.

First, a univariate logistic regression analysis was conducted for each variable. Variables with p < 0.05 were then integrated and modeled using R software (RRID: SCR_001905) (version 4.2.3, https://www.r-project.org/) for clinical and imaging features, respectively. Covariate diagnostics were performed for each model’s variables. Multifactor logistic regression analysis was then conducted by combining indicators with p < 0.05 from both models, and nomograms were developed.

In addition, the area under the curve (AUC) of the Receiver operating characteristic (ROC) of the subjects was plotted in this study to assess the differential diagnostic ability of each model. The DeLong test was used to compare the AUC values between the models. The Hosmer-Lemeshow goodness-of-fit test was used to assess the goodness-of-fit of the nomograms, and calibration curves were plotted to assess the agreement between predicted and actual results. In addition, the performance of all models was assessed by calculating sensitivity, specificity, and accuracy. Finally, clinical decision curve analysis (DCA) was performed in order to understand the patient benefit. The data from center 2 patients were used as an external validation cohort to verify the robustness, generalization ability of each model. There were no statistically significant differences in most of the characteristics of the patients in the two centers (Table 1).

Table 1

Training cohort (N=272)Validation cohort (N=76)
PD responders (N=66)non-PD responders (N=206)p-valuePD responders (N=22)non-PD responders (N=54)p-value
Age (%)0.4560.819
≤60 years31 (47.0)86 (41.7)10 (45.5)23 (42.6)
>60 years35 (53.0)120 (58.3)12 (54.5)31 (57.4)
Gender (%)0.2240.770
Male45 (68.2)156 (75.7)17 (77.3)40 (74.1)
Female21 (31.8)50 (24.3)5 (22.7)14 (25.9)
Immunotherapy Regimen (%)0.670.729
Camrelizumab36 (54.5)91 (44.2)11 (50.0)25 (46.3)
Sintilimab23 (34.8)88 (42.7)10 (45.5)24 (44.4)
Nivolumab3 (4.5)11 (5.3)0 (0.0)3 (5.6)
Tislelizumab3 (4.5)10 (4.9)1 (4.5)2 (3.7)
Pembrolizumab1 (1.5)6 (2.9)00
Tumor location (%)0.8970.847
Cardia22 (33.3)74 (35.9)8 (36.4)19 (35.2)
Body24 (36.4)69 (33.5)10 (45.5)22 (40.7)
Antrum20 (30.3)63 (30.6)4 (18.2)13 (24.1)
Lymph nodes with a short diameter>1 cm (%)0.8980.77
Absence19 (28.8)61 (29.6)6 (27.3)13 (24.1)
Presence47 (71.2)145 (70.4)16 (72.7)41 (75.9)
Obscuration of the perigastric fat space (%)0.4150.935
Absence9 (13.6)37 (18.0)4 (18.2)12 (22.2)
Presence57 (86.4)169 (82.0)18 (81.8)42 (77.8)
Serosal invasion (%)0.0020.387
Absence3 (4.5)43 (20.9)1 (4.5)8 (14.8)
Presence63 (95.5)163 (79.1)21 (95.5)46 (85.2)
Distant metastasis (%)<0.0010.002
Absence19 (28.8)145 (70.4)4 (18.2)31 (57.4)
Presence47 (71.2)61 (29.6)18 (81.8)23 (42.6)
BMI (kg/m2)23.22 (21.88, 23.73)23.26 (21.36, 24.61)0.66323.23 ± 3.0423.07 ± 3.280.843
Maximum tumor thickness (cm)2.21 (1.50, 2.98)2.10 (1.63, 2.70)0.8422.100 (1.825,2.675)2.200 (1.500-2.675)0.566
Maximum tumor diameter (cm)5.55 (4.03, 8.28)5.850 (4.00, 7.68)0.7186.38 ± 2.405.99 ± 2.860.573
CT values in the unenhanced phase (HU)36.850 (32.33, 39.10)36.850 (33.10, 40.28)0.84331.250 (30.13, 35.40)32.450 (28.73, 35.20)0.714
CT values in the arterial phase (HU)59.90 (47.50, 71.38)67.350 (54.80, 77.15)0.00154.050 (49.58, 61.18)55.300 (49.60, 62.68)0.457
CT values in the venous phase (HU)77.15 (67.73, 93.23)77.300 (65.58, 88.85),0.41966.750 (61.03, 76.20)66.450 (60.63, 74.90)0.909
CT values in the delayed phase (HU)76.80 (67.10, 87.60)76.050 (67.78, 87.83)0.88363.950 (60.50, 73.48)69.950 (64.18, 77.90)0.159
Tumor arterial attenuation (HU)21.55 (14.30, 31.48)30.100 (21.33, 37.98)<0.00120.350 (14.35, 28.80)23.100 (18.53, 29.88)0.186
Tumor venous attenuation on portal phase (HU)40.25 (32.33, 55.00)39.650 (31.45, 50.98)0.53335.200 (29.00, 45.25)36.600 (27.13, 44.20)0.766
AEF (%)55.12 (45.66, 63.75)71.700 (57.98, 83.61)<0.00162.10 (51.75, 68.35)70.740 (58.68, 75.28)0.004
CEA (ng/ml)7.96 (2.63, 16.85)3.98 (1.96, 15.10)0.03713.77 (3.91, 48.63)4.66 (2.00, 34.34)0.162
AFP (ng/ml)3.00 (1.90, 11.08)2.50 (1.60, 6.80)0.0953.68 (2.33, 18.42)3.67 (2.11, 9.26)0.532
CA125 (U/ml)25.17 (12.15, 77.75)14.40 (8.87, 31.70)0.00151.79 (20.88, 187.25)33.21 (18.75, 79.28)0.236
CA199 (U/ml)22.34 (12.15, 83.45)16.95 (10.25, 43.75)0.14620.25 (16.60, 126.68)22.10 (6.35, 110.13)0.287
CA724 (U/ml)14.90 (4.00, 140.50)5.66 (2.30, 21.80)0.00116.20 (9.70, 32.33)8.10 (3.53, 18.23)0.052
PLT (109/L)223.00 (150.75, 276.75)255.500 (199.00, 300.50)0.037259.27 ± 87.60244.30 ± 91.640.515
lymphocyte (109/L)1.26 (0.98, 1.71)1.525 (1.20, 1.86)0.0031.29 (1.09, 1.71)1.34 (1.12, 1.78)0.433
Neutrophilcount (109/L)3.58 (2.67, 4.90)3.26 (2.43, 4.34)0.1183.22 (2.81, 4.60)3.77 (2.79, 4.94)0.397
Monocyte (109/L)0.50 (0.35, 0.65)0.47 (0.36, 0.59)0.720.51 (0.41, 0.68)0.54 (0.42, 0.69)0.868
SII541.76 (350.56, 1044.94)517.41 (323.15, 893.94)0.302657.72 (411.83, 970.80)666.99 (361.34, 1008.44)0.805
MLR0.39 (0.29, 0.53)0.32 (0.24, 0.43)0.0020.36 (0.30, 0.53)0.41 (0.28, 0.50)0.779
NLR2.65 (1.86, 4.56)2.05 (1.53, 3.15)0.0022.69 (2.08,3.98)2.79 (2.09,4.01)0.963

Baseline characteristics of the training cohort and the external validation cohort.

3 Results

3.1 Patients

A total of 272 patients (206 non-PD responders and 66 PD responders) from Center 1 and 76 patients (54 non-PD responders and 22 PD responders) from Center 2 were included in this study. Patients from center 1 were used as a training set to build the model, and patients from center 2 were used as an external validation cohort to test the robustness of the model.

3.2 Clinical features

Each enrolled patient was classified as either a PD responder or non-responder based on treatment efficacy after receiving four cycles of PD-1 inhibitor combination chemotherapy. In the training cohort, no statistically significant differences were observed between the PD responder and non-PD responder groups in terms of age, gender, and body mass index (BMI) (p > 0.05). However, the presence of distant metastases was significantly different between the two groups (p < 0.05). Laboratory tests revealed no significant differences between the PD responder and non-PD responder groups in alpha-fetoprotein, carcinoembryonic antigen, carbohydrate antigen 125 (CA 125), carbohydrate antigen 19-9 (CA 19-9), platelet count, or absolute monocyte count (AMC). Conversely, significant differences were observed in carbohydrate antigen 72-4 (CA 72-4), absolute neutrophil count (ANC), absolute lymphocyte count (ALC), neutrophil-to-lymphocyte ratio (NLR), monocyte-to-lymphocyte ratio (MLR), and systemic immune-inflammation index (SII) (p < 0.05). Univariate and multivariate logistic regression analyses of clinical characteristics are presented in Table 2.

Table 2

VariableUnivariate analysisMultivariate analysis
OR (95%CI)P-valueOR (95%CI)P-value
Age (year)1.24 (0.71-2.16)0.456
Gender0.69 (0.37-1.26)0.226
BMI (kg/m2)1.03 (0.93-1.14)0.569
Distant metastasis0.17 (0.09-0.31)<0.0010.22 (0.11-0.44)<0.001
AFP (ng/ml)1.00 (1.00-1.00)0.161
CEA (ng/ml)1.00 (1.00-1.00)0.352
CA125 (U/ml)1.00 (1.00-1.00)0.277
CA199 (U/ml)1.00 (1.00-1.00)0.265
CA724 (U/ml)0.99 (0.99-0.99)<0.0010.99 (0.98-0.99)<0.001
PLT (109/L)1.00 (1.00-1.01)0.06
lymphocyte (109/L)2.47 (1.33-4.59),0.0040.73 (0.18-2.92)0.659
Monocyte (109/L)0.62 (0.16-2.47)0.501
Neutrophilcount (109/L)0.86 (0.75-0.99)0.0351.26 (0.78-2.04)0.353
NLR0.71 (0.60-0.85)<0.0010.43 (0.23-0.83)0.012
MLR0.06 (0.01-0.35)0.0021.72 (0.13-22.53)0.681
SII1.00 (1.00-1.00)0.0461.00 (1.00-1.00)0.057

Univariate and multivariate logistic regression analysis of clinical features.

3.3 Computed tomography imaging features

CT examination showed no statistically significant differences between the two groups regarding tumor location, maximum diameter, maximum thickness, whether the surrounding fat was infiltrated, or the presence of enlarged lymph nodes in the surrounding area (p > 0.05). However, the presence of serosal invasion (involving the visceral peritoneum) showed a statistically significant difference between the PD responder and non-PD responder groups (p < 0.05). There were no statistically significant differences in CT values or tumor venous attenuation values between the plain, venous, and delayed phases (p > 0.05). Notably, CT values, arterial attenuation values, and arterial enhancement fraction (AEF) of tumors in the arterial phase were significantly higher in the non-PD responder group compared to the PD responder group (p < 0.05). Univariate and multivariate logistic regression analyses of CT imaging features are shown in Table 3.

Table 3

VariableUnivariate analysisMultivariate analysis
OR (95%CI)P-valueOR (95%CI)P-value
Tumor location0.97 (0.69-1.36)0.897
Maximum tumor thickness (cm)0.86 (0.65-1.13)0.283
Maximum tumor diameter (cm)1.00 (0.91-1.11)0.965
Lymph nodes with a short diameter>1 cm0.96 (0.52-1.77)0.898
Obscuration of the perigastric fat space0.72 (0.33-1.59)0.416
Serosal invasion0.18 (0.05-0.60)0.0050.12 (0.03-0.43)0.001
CT values in the unenhanced phase (HU)0.99 (0.95-1.04)0.790
CT values in the arterial phase (HU)1.04 (1.02-1.06)0.0010.98 (0.92-1.04)0.548
CT values in the venous phase (HU)0.99 (0.98-1.01)0.439
CT values in the delayed phase (HU)1.00 (0.98-1.01)0.906
Tumor Arterial attenuation (HU)1.05 (1.03-1.08)<0.0011.00 (0.93-1.08)0.927
Tumor venous attenuation on portal phase (HU)0.99 (0.98-1.01)0.439
AEF (%)1.08 (1.06-1.11)<0.0011.10 (1.06-1.13)<0.001

Univariate and multivariate logistic regression analysis of imaging features.

3.4 Development of differentiation models for PD and non-PD patients

Comparative analyses of laboratory parameters and imaging features were conducted to develop the most effective model for optimal discrimination. Initially, all indicators were analyzed individually using univariate logistic regression, and those with p < 0.05 were included in subsequent multivariate logistic regression analysis. Clinical and imaging feature models were established separately, with covariate diagnosis performed for each. The clinical feature model included carbohydrate antigen 72-4 (CA 72-4), the neutrophil-to-lymphocyte ratio (NLR), and the presence of distant metastasis. The imaging feature models included the arterial enhancement fraction (AEF) and the presence or absence of visceral peritoneum infiltration. The area under the curve (AUC) for the clinical and imaging feature models in the training cohort were 0.792 and 0.829, respectively, while in the external validation cohort, the AUCs were 0.713 and 0.729, respectively. The ROC curves for the two models are shown in Figure 2.

Figure 2

3.5 Development and evaluation of the nomogram

Indicators with p < 0.05 from both the clinical feature model and the imaging feature model were selected to construct the nomogram. These included the neutrophil-to-lymphocyte ratio (NLR), CA 72-4, the presence of distant metastasis, the presence of visceral peritoneum infiltration, and the arterial enhancement fraction (AEF). The nomogram is illustrated in Figure 3.

Figure 3

The AUC for the nomogram was 0.904 in the training cohort and 0.801 in the external validation cohort (Figure 4). The DeLong test indicated a statistically significant difference between the nomogram and both the clinical and imaging feature models in both cohorts (p < 0.05). According to the Hosmer-Lemeshow goodness-of-fit test, the nomogram demonstrated good fit in both the training cohort (p = 0.236) and the external validation cohort (p = 0.191). The nomogram also showed strong agreement with the calibration curves for both cohorts (Figure 5). Additionally, decision curve analysis for the clinical feature model, imaging feature model, and nomogram is depicted in Figure 6, showing that the nomogram provided the highest net benefit. The accuracy, specificity, and sensitivity of all models are detailed in Table 4.

Figure 4

Figure 5

Figure 6

Table 4

CohortModelsAUCAccuracySpecificitySensitivity
Training cohortClinical model0.792 (95% CI:0.730-0.855)0.7170.8030.689
Imaging model0.829 (95% CI:0.775-0.883)0.7280.8330.694
Nomogram (clinical+imaging)0.904 (95% CI:0.862-0.947)0.8490.8180.859
External validation cohortClinical model0.713 (95% CI:0.593-0.833)0.6840.7270.667
Imaging model0.729 (95% CI:0.604-0.854)0.6970.6820.704
Nomogram (clinical+imaging)0.801 (95% CI:0.683-0.918)0.8290.6820.889

The performance of models in the training cohort and external validation cohort.

4 Discussion

In this study, we used clinical and CT imaging features before treatment to predict the early treatment response of gastric cancer patients receiving PD-1 inhibitors combined with chemotherapy. The nomogram that integrated both clinical and imaging features showed the best performance, with AUCs of 0.904 in the training set and 0.801 in the validation set. This performance surpassed both the clinical feature model and imaging feature model alone and slightly outperformed the radiomics model by Liang et al. (). In a previous study, Li et al. () developed a nomogram to predict complete pathological remission in gastric cancer patients undergoing immunotherapy combined with chemotherapy, using tumor diameter, clinical N stage, and the combined positive score (CPS) of PD-L1. In contrast, our nomogram incorporates readily available, non-invasive biomarkers, eliminating the need for immunohistochemical markers like PD-L1 CPS from biopsy samples, thereby simplifying data acquisition while maintaining high predictive accuracy.

Liang et al. () previously constructed a model using radiomics-extracted imaging features combined with two clinical risk factors to predict the efficacy of PD-1 inhibitors in gastric cancer patients, achieving an AUC of 0.865 in the training set and 0.778 in the internal validation set. While their model demonstrated strong predictive performance, our nomogram offers a more accessible approach by incorporating standard CT imaging features, simplifying the process of obtaining lesion data. Moreover, our nomogram’s AUC, sensitivity, and specificity are slightly higher than Liang’s model, and we further validated our nomogram externally, supporting its generalizability and robustness.

In this study, we initially developed a clinical characterization model, which achieved an area under the curve (AUC) of 0.792 in the training cohort. However, this model exhibited a lower sensitivity of 0.689, and the decision curve analysis (DCA) indicated a reduced patient benefit. Subsequently, we developed an imaging feature model based on CT imaging characteristics. This model demonstrated superior performance compared to the clinical feature model, with an AUC of 0.829, and improved accuracy and specificity. Nevertheless, the sensitivity of the CT imaging model was comparable to that of the clinical feature model, at 0.694. Finally, we integrated the statistically significant predictors (p < 0.05) from both models to construct a nomogram, which simplified the metrics while enhancing diagnostic efficacy. The nomogram exhibited excellent predictive performance, achieving an AUC of 0.904, surpassing both the CT imaging feature model and the clinical feature model. Additionally, the nomogram demonstrated sensitivity, specificity, and accuracy exceeding 0.80, and showed promising results in external validation with an AUC of 0.806.

We constructed the nomogram based on five relevant predictive factors (p < 0.05) identified from the clinical feature model and imaging feature model. The clinical model included three key factors: neutrophil-to-lymphocyte ratio (NLR), CA724, and the presence of distant metastases. NLR (neutrophils/lymphocytes) has been found to be positively correlated with the density of exhausted CD8 T immune cells in the tumor microenvironment (), and exhausted CD8 T cells represent a state of progressively diminished T-cell function (). High levels of serum carbohydrate antigen 72-4 (CA72-4) typically indicate the persistence and proliferation of tumor cells (), and previous studies have demonstrated that CA72-4 can predict chemotherapy efficacy in gastric cancer patients () as well as prognosis following immunosuppressive therapy (). The presence of tumor metastasis is considered a key factor influencing the effectiveness of PD-1 inhibitor therapy in gastric cancer patients (). Tumor cells and the tumor immune microenvironment (TME) within metastatic lesions may impact the local immune microenvironment of distant tissues via circulation (), thus affecting the response to immunotherapy.

Relevant predictors (p < 0.05) from the imaging characterization model used in our nomogram included arterial enhancement fraction (AEF) and the presence or absence of serosal invasion (visceral peritoneal). Our study revealed that tumors with a higher AEF on pre-treatment CT images were more likely to show a better response to treatment in patients with advanced gastric cancer undergoing PD-1 inhibitor combination chemotherapy. This observation may be related to the association between higher AEF values and elevated composite positive score (CPS) for PD-L1 expression (), with gastric cancer patients who have high CPS scores for PD-L1 potentially deriving greater benefit from immunotherapy (). Additionally, serosal invasion (visceral peritoneal) is a significant factor affecting the prognosis of gastric cancer patients (). Such invasion may contribute to the development of an immunosuppressive inflammatory tumor microenvironment through chronic inflammation (, ), which can subsequently impact the effectiveness of immunotherapy.

In our study, we developed a straightforward and effective nomogram using the predictors mentioned above, which are easily accessible and non-invasive. In contrast to individual biomarker predictions, our nomogram assists clinicians in forecasting treatment efficacy and crafting personalized treatment strategies. By quantifying each predictor as numerical probabilities, it enables more accurate early-stage predictions of treatment response for gastric cancer patients undergoing PD-1 inhibitor combination chemotherapy.

As with most retrospective studies, our research has several limitations. Firstly, the study population was relatively small, indicating the need for larger, multicenter prospective studies to confirm our findings. Secondly, although CT imaging is a crucial tool for assessing the efficacy of immunotherapy in gastric cancer, discrepancies may arise between CT assessments and pathological findings, which can affect the comprehensiveness and accuracy of the efficacy evaluation. Since this was a retrospective study, most patients did not undergo surgical resection and therefore the degree of their pathological remission could not be determined, further limiting our ability to comprehensively assess treatment efficacy. Therefore, in future studies, we plan to include more cases and use the degree of pathological remission as a criterion for efficacy assessment to enhance the accuracy and comprehensiveness of evaluating immunotherapy efficacy.

5 Conclusion

In summary, we developed and validated a nomogram utilizing clinical features before treatment (NLR, CA724, and the presence or absence of distant metastases) alongside CT imaging features before treatment (arterial enhancement fraction and serosal invasion). This nomogram serves as a simple and effective tool for early prediction of therapeutic efficacy in patients with advanced gastric cancer receiving PD-1 inhibitor-combined chemotherapy, thereby aiding clinicians in formulating personalized treatment plans.

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 authors.

Author contributions

YM: Writing – original draft, Conceptualization, Data curation, Investigation, Methodology, Validation. ZW: Writing – original draft, Methodology, Supervision. CQ: Writing – original draft, Methodology, Validation. MX: Writing – original draft. SW: Writing – original draft. KH: Writing – original draft. HX: Writing – review & editing. HW: Writing – review & editing.

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 Key Technology Research and Development Program of Shandong Province, China (2021SFGC0104) and the Key Technology Research and Development Program, China (2022YFC2408400).

Acknowledgments

We thank all authors for their contributions to the conceptualization and design of the study.

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.

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/fimmu.2025.1504387/full#supplementary-material

References

  • 1

    SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660

  • 2

    JohnstonFMBeckmanM. Updates on management of gastric cancer. Curr Oncol Rep. (2019) 21:67. doi: 10.1007/s11912-019-0820-4

  • 3

    TakeiSKawazoeAShitaraK. The new era of immunotherapy in gastric cancer. Cancers (Basel). (2022) 14. doi: 10.3390/cancers14041054

  • 4

    YiMZhengXNiuMZhuSGeHWuK. Combination strategies with pd-1/pd-L1 blockade: current advances and future directions. Mol Cancer. (2022) 21:28. doi: 10.1186/s12943-021-01489-2

  • 5

    JanjigianYYShitaraKMoehlerMGarridoMSalmanPShenLet al. First-line nivolumab plus chemotherapy versus chemotherapy alone for advanced gastric, gastro-oesophageal junction, and oesophageal adenocarcinoma (Checkmate 649): A randomised, open-label, phase 3 trial. Lancet. (2021) 398:27–40. doi: 10.1016/S0140-6736(21)00797-2

  • 6

    MareiHEHasanAPozzoliGCenciarelliC. Cancer immunotherapy with immune checkpoint inhibitors (Icis): potential, mechanisms of resistance, and strategies for reinvigorating T cell responsiveness when resistance is acquired. Cancer Cell Int. (2023) 23:64. doi: 10.1186/s12935-023-02902-0

  • 7

    TopalianSLHodiFSBrahmerJRGettingerSNSmithDCMcDermottDFet al. Safety, activity, and immune correlates of anti-pd-1 antibody in cancer. N Engl J Med. (2012) 366:2443–54. doi: 10.1056/NEJMoa1200690

  • 8

    ChoYALeeHKimDGKimHHaSYChoiYLet al. Pd-L1 expression is significantly associated with tumor mutation burden and microsatellite instability score. Cancers (Basel). (2021) 13. doi: 10.3390/cancers13184659

  • 9

    WangMRanXLeungWKawaleASaxenaSOuyangJet al. Atr inhibition induces synthetic lethality in mismatch repair-deficient cells and augments immunotherapy. Genes Dev. (2023) 37:929–43. doi: 10.1101/gad.351084.123

  • 10

    ChengYBuDZhangQSunRLyleSZhaoGet al. Genomic and transcriptomic profiling indicates the prognosis significance of mutational signature for tmb-high subtype in chinese patients with gastric cancer. J Adv Res. (2023) 51:121–34. doi: 10.1016/j.jare.2022.10.019

  • 11

    YangYShiZBaiRHuW. Heterogeneity of msi-H gastric cancer identifies a subtype with worse survival. J Med Genet. (2021) 58:12–9. doi: 10.1136/jmedgenet-2019-106609

  • 12

    WangZLLiYLLiXTTangLLiZYSunYS. Role of ct in the prediction of pathological complete response in gastric cancer after neoadjuvant chemotherapy. Abdom Radiol (NY). (2021) 46:3011–8. doi: 10.1007/s00261-021-02967-3

  • 13

    FuMZhangXShenFMaJLiZ. Prognostic value of peripheral blood neutrophil/lymphocyte ratio, platelet/lymphocyte ratio, pan-immune-inflammation value and systemic immune-inflammation index for the efficacy of immunotherapy in patients with advanced gastric cancer. Immunotherapy. (2024) 16(8):551–563. doi: 10.2217/imt-2024-0031

  • 14

    HouYLiXYangYShiHWangSGaoM. Serum cytokines and neutrophil-to-lymphocyte ratio as predictive biomarkers of benefit from pd-1 inhibitors in gastric cancer. Front Immunol. (2023) 14:1274431. doi: 10.3389/fimmu.2023.1274431

  • 15

    YushkevichPAPivenJHazlettHCSmithRGHoSGeeJCet al. User-guided 3d active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. (2006) 31:1116–28. doi: 10.1016/j.neuroimage.2006.01.015

  • 16

    XuMLiuSQiaoXLiLJiCZhouZ. Clinicopathological features and ct findings of papillary gastric adenocarcinoma. Abdom Radiol (NY). (2022) 47:3698–711. doi: 10.1007/s00261-022-03635-w

  • 17

    LiuYKimJQuFLiuSWangHBalagurunathanYet al. Ct features associated with epidermal growth factor receptor mutation status in patients with lung adenocarcinoma. Radiology. (2016) 280:271–80. doi: 10.1148/radiol.2016151455

  • 18

    KimKWLeeJMKlotzEParkHSLeeDHKimJYet al. Quantitative ct color mapping of the arterial enhancement fraction of the liver to detect hepatocellular carcinoma. Radiology. (2009) 250:425–34. doi: 10.1148/radiol.2501072196

  • 19

    ZhouYXuYKGengDWangJWChenXBSiYet al. Added value of arterial enhancement fraction derived from dual-energy computed tomography for preoperative diagnosis of cervical lymph node metastasis in papillary thyroid cancer: initial results. Eur Radiol. (2024) 34:1292–301. doi: 10.1007/s00330-023-10109-0

  • 20

    LiangZHuangAWangLBiJKuangBXiaoYet al. A radiomics model predicts the response of patients with advanced gastric cancer to pd-1 inhibitor treatment. Aging (Albany NY). (2022) 14:907–22. doi: 10.18632/aging.203850

  • 21

    LiLChenGChenEYStricklandMRZhaoWZhangJet al. Development and validation of a nomogram to predict pathological complete response in patients with locally advanced gastric adenocarcinoma treated with neoadjuvant chemotherapy in combination with pd-1 antibodies. J Gastrointest Oncol. (2023) 14:2373–83. doi: 10.21037/jgo-23-751

  • 22

    YuanJZhaoXLiYYaoQJiangLFengXet al. The association between blood indexes and immune cell concentrations in the primary tumor microenvironment predicting survival of immunotherapy in gastric cancer. Cancers (Basel). (2022) 14. doi: 10.3390/cancers14153608

  • 23

    DolinaJSVan-Braeckel-BudimirNThomasGDSalek-ArdakaniS. Cd8(+) T cell exhaustion in cancer. Front Immunol. (2021) 12:715234. doi: 10.3389/fimmu.2021.715234

  • 24

    XuYZhangPZhangKHuangC. The application of ca72-4 in the diagnosis, prognosis, and treatment of gastric cancer. Biochim Biophys Acta Rev Cancer. (2021) 1876:188634. doi: 10.1016/j.bbcan.2021.188634

  • 25

    ZhuYChenNChenMCuiXYangHZhuXet al. Circulating tumor cells: A surrogate to predict the effect of treatment and overall survival in gastric adenocarcinoma. Int J Biol Markers. (2021) 36:28–35. doi: 10.1177/1724600820981972

  • 26

    SunHChenLHuangRPanHZuoYZhaoRet al. Prognostic nutritional index for predicting the clinical outcomes of patients with gastric cancer who received immune checkpoint inhibitors. Front Nutr. (2022) 9:1038118. doi: 10.3389/fnut.2022.1038118

  • 27

    MishimaSKawazoeANakamuraYSasakiAKotaniDKubokiYet al. Clinicopathological and molecular features of responders to nivolumab for patients with advanced gastric cancer. J Immunother Cancer. (2019) 7:24. doi: 10.1186/s40425-019-0514-3

  • 28

    OliverAJKeamSPvon ScheidtBZankerDJHarrisonAJTantaloDGet al. Primary and metastatic breast tumors cross-talk to influence immunotherapy responses. Oncoimmunology. (2020) 9:1802979. doi: 10.1080/2162402X.2020.1802979

  • 29

    XuMRenTDengJYangJLuTXiHet al. Correlation of ct parameters and pd-L1 expression status in gastric cancer. Abdom Radiol (NY). (2024) 49:1320–9. doi: 10.1007/s00261-024-04200-3

  • 30

    YuHYLiCPHuangYHHsuSJWangYPHsiehYCet al. Microsatellite instability, epstein-barr virus, and programmed cell death ligand 1 as predictive markers for immunotherapy in gastric cancer. Cancers (Basel). (2022) 14. doi: 10.3390/cancers14010218

  • 31

    NakamuraKTajimaKKanamoriKYatabeKOgimiMHiguchiTet al. Impact of subclassification of serosal invasion on the survival of patients with T4a gastric cancer. In Vivo. (2022) 36:1923–9. doi: 10.21873/invivo.12913

  • 32

    YuraMYoshikawaTWadaTOtsukiSHayashiTYamagataYet al. The prognostic impact of macroscopic serosal change on resectable advanced gastric cancer. BMC Cancer. (2021) 21:1056. doi: 10.1186/s12885-021-08767-8

  • 33

    ZhaoHWuLYanGChenYZhouMWuYet al. Inflammation and tumor progression: signaling pathways and targeted intervention. Signal Transduct Target Ther. (2021) 6:263. doi: 10.1038/s41392-021-00658-5

Summary

Keywords

nomogram, gastric cancer, computed tomography, programmed cell death 1 inhibitors, chemo-immunotherapy

Citation

Ma Y, Wang Z, Qiu C, Xiao M, Wu S, Han K, Xu H and Wang H (2025) Nomogram based on CT imaging and clinical data to predict the efficacy of PD-1 inhibitors combined with chemotherapy in advanced gastric cancer. Front. Immunol. 16:1504387. doi: 10.3389/fimmu.2025.1504387

Received

30 September 2024

Accepted

12 March 2025

Published

31 March 2025

Volume

16 - 2025

Edited by

Abhishek Mahajan, The Clatterbridge Cancer Centre, United Kingdom

Reviewed by

Dmitry Aleksandrovich Zinovkin, Gomel State Medical University, Belarus

Junmeng Li, Henan Provincial People’s Hospital, China

Updates

Copyright

*Correspondence: Hui Xu, ; Haiyan Wang,

†These authors share first authorship

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics