- 1Department of Medical Oncology, The First Affiliated Hospital of Bengbu Medical College, Bengbu, China
- 2Department of Radiation Oncology, The First Affiliated Hospital of Bengbu Medical College, Bengbu, China
- 3Department of Clinical Laboratory, The First Affiliated Hospital of Bengbu Medical College, Bengbu, China
- 4Department of Surgical Oncology, The First Affiliated Hospital of Bengbu Medical College, Bengbu, China
- 5Department of Cancer Biology, Mayo Clinic, Jacksonville, FL, United States
- 6Department of Oncology, The First Affiliated Hospital of Anhui Medical University, Hefei, China
Purpose: To explore the predictive value of multiple immune-inflammatory biomarkers including serum VEGFA and systemic immune-inflammation index (SII) in HER2-negative advanced gastric cancer (AGC) and establish nomograms for predicting the first-line chemotherapeutic efficacy, progression-free survival (PFS) and overall survival (OS) of patients with this fatal disease.
Methods: From November 2017 to April 2022, 102 and 34 patients with a diagnosis of HER2-negative AGC at the First Affiliated Hospital of Bengbu Medical College were enrolled as development and validation cohorts, respectively. Univariate and multivariate analyses were performed to evaluate the clinical value of the candidate indicators. The variables were screened using LASSO regression analysis. Predictive models were developed using significant predictors and are displayed as nomograms.
Results: Baseline VEGFA expression was significantly higher in HER2-negative AGC patients than in nonneoplastic patients and was associated with malignant serous effusion and therapeutic efficacy (all p<0.001). Multivariate analysis indicated that VEGFA was an independent predictor for first-line therapeutic efficacy and PFS (both p<0.01) and SII was an independent predictor for first-line PFS and OS (both p<0.05) in HER2-negative AGC patients. The therapeutic efficacy model had an R2 of 0.37, a Brier score of 0.15, and a Harrell’s C-index of 0.82 in the development cohort and 0.90 in the validation cohort. The decision curve analysis indicated that the model added more net benefits than VEGFA assessment alone. The PFS/OS models had Harrell’s C-indexes of 0.71/0.69 in the development cohort and 0.71/0.62 in the validation cohort.
Conclusion: The established nomograms integrating serum VEGFA/SII and commonly available baseline characteristics provided satisfactory performance in predicting the therapeutic efficacy and prognosis of HER2-negative AGC patients.
1 Introduction
Gastric cancer (GC) is the fifth most diagnosed cancer and the fourth leading cause of cancer-related death worldwide (1). Currently, for human epidermal growth factor receptor 2 (HER2)-negative GC, chemotherapy remains the cornerstone of systemic therapy for patients with advanced disease despite the emergence of immune checkpoint inhibitors (ICIs) (2, 3). Due to the highly heterogeneous and malignant nature of GC (4), patient response to chemotherapy varies greatly from person to person, with survival varying widely; most patients with advanced disease die within 12 months of diagnosis (5). Therefore, it is critical to identify new approaches to predict the efficacy and prognosis of chemotherapy to make more appropriate treatment decisions and improve survival in patients with advanced gastric cancer (AGC).
The immune-inflammatory response is widely accepted to play important roles in tumor initiation, invasion, angiogenesis, and metastasis (6, 7). Angiogenesis is one of the hallmarks of malignancy and contributes significantly to the growth, invasion, and metastatic spread of cancer (8). As a major marker of angiogenesis, vascular endothelial growth factor A (VEGFA) expression is associated with the treatment response (9) and prognosis of various malignancies (10, 11). In addition, other immune-inflammatory biomarkers such as systemic immune-inflammation index (SII) (12), c-reactive protein/albumin ratio (CAR) (13); nutritional markers such as body mass index (BMI) (14), hemoglobin (Hb) (15), albumin/globulin ratio (AGR) (16), and lipoprotein cholesterol (17); and the classic tumor marker carcinoembryonic antigen (CEA) (18) were associated with GC or other malignancies. Since the predictive efficacy of a single indicator is limited, there is a need to rationally evaluate diverse indicators of treatment efficacy and prognosis and to select appropriate parameters to establish predictive models that can comprehensively assess treatment effectiveness and provide individualized treatment for patients with malignant tumors.
Currently, most of the predictive models for GC are diagnostic and prognostic models (19, 20) and are based on retrospective studies; some of them are based on gene expression and are difficult to be applied clinically (21, 22). This prospective, single-center, cohort study aimed to predict the effectiveness of first-line chemotherapy and the treatment outcome for HER2-negative AGC using clinicopathologic features and laboratory hematologic indicators to provide convenient and general predictive models for the treatment of AGC and to improve clinical decision making and patient−physician communication.
2 Materials and methods
2.1 Patient selection
From November 2017 to April 2022, we prospectively enrolled 102 and 34 patients who attended the First Affiliated Hospital of Bengbu Medical College for primary treatment of recurrent or metastatic GC as development and validation cohorts. The inclusion criteria were as follows: patients with histologically confirmed HER2-negative AGC; patients who had not received antitumor therapy against recurrent or metastatic GC; patients who were not suitable for or were unwilling to undergo surgery or radiation therapy; and patients with target lesions that could be evaluated for efficacy. The exclusion criteria were as follows: patients with a combination of other tumors or subtypes; patients with severe cardiac, hepatic, or renal diseases; and patients with severe bleeding or infectious diseases. Moreover, 108 patients with nonneoplastic disease without bleeding or infectious diseases were retrospectively enrolled as a control for the development cohort and 15 HER2-positive AGC patients treated with first-line trastuzumab combined with chemotherapy from September 2019 to April 2022 were retrospectively enrolled as a HER2-positive cohort. The present study was reported in accordance with the STROBE (23) and TRIPOD (24) guidelines as much as possible. This study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of the First Affiliated Hospital of Bengbu Medical College (BYYFY-2017KY09). Informed consent was obtained from all patients or their family members.
2.2 Bioinformatics analysis
VEGFA mRNA levels in several types of cancer, including GC, were analyzed using the Tumor Immune Estimation Resource (TIMER2.0) database with The Cancer Genome Atlas (TCGA) data (http://timer.cistrome.org/). The GDC dataset containing the gene expression profiles (HTSeq-Counts) of 27 pairs of GC tissues and their matched normal tissues was downloaded from The Cancer Genome Atlas stomach adenocarcinoma (TCGA-STAD, https://xenabrowser.net/). To ensure uniformity across the data, the RNA-seq profiles were transformed into TPM values, and the formula log2 (TPM+1) was used for normalization. GC samples from TCGA-STAD were selected for survival analysis. Using the Kaplan−Meier plotter database (25), the prognostic role of VEGFA mRNA expression in HER2-negative GC patients at stage IV was analyzed (Probe ID: 211527_at_x).
2.3 Predictive variables
Clinicopathological information such as sex, age, height, weight, Eastern Cooperative Oncology Group performance status (ECOG PS), diagnosis pattern (including recurrent or primary cancer), histology [including adenocarcinoma (AC) and signet-ring cell carcinoma (SRC)], position of the primary lesion, number and sites of metastases, comorbidities (such as malignant serous effusion), and laboratory indicators [such as VEGFA, CEA, total protein (TP), albumin (ALB), Hb, C-reactive protein (CRP), low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C)] were collected at baseline for each patient. Composite indicators were calculated based on the following formulas: AGR=ALB/(TP-ALB), BMI=height/weight2, CAR=CRP/ALB, LHR=LDL-C/HDL-C, and SII=platelet count × neutrophil count/lymphocyte count.
2.4 Detection of VEGFA and other laboratory indicators
Five milliliter samples of nonanticoagulated peripheral venous blood and 3 ml of EDTA-K2 anticoagulated peripheral venous blood were collected from each patient at baseline. The nonanticoagulated peripheral venous blood was centrifuged at 3000 rpm for 10 min to separate the serum, and the supernatant was quickly frozen and stored at -20°C. VEGFA detection was performed using a Weigao JR-1 Chemiluminescent Immunoassay Analyzer and Vascular Endothelial Growth Factor Assay Kit (chemiluminescence, Shandong Weigao Group Medical Polymer Co., Ltd., Weihai, China) according to the manufacturer’s instructions. A Sysmex XE-2100 automatic blood analyzer (Sysmex Corporation, Kobe, Japan) was used for hematology analysis. According to the standards of laboratory SOP documents and reagent instructions, CEA measured by chemiluminescence microparticle immunoassay was detected using an Abbott I-2000 chemiluminescence immunoanalyzer, and all reagents, calibration and quality control products were obtained from Abbott Laboratories Trading Co., Ltd. (Shanghai, China). ALB/TP (measured by colorimetric assay), LDL-C/HDL-C (measured by homogeneous enzyme colorimetric assay), and CRP (measured by immunoturbidimetric assay) were detected using a Roche Cobas 8000 C 701 automatic biochemical analyzer, and all reagents, calibration and quality control products were obtained from Roche Diagnostics Co., Ltd. (Shanghai, China). All the test technicians were professionally trained and did not have access to clinical data.
2.5 Therapeutic regimens
All patients received at least one cycle of first-line chemotherapy according to the NCCN guidelines for gastric cancer (version 5. 2017) (26). Two-drug regimens (fluorouracil or capecitabine/oxaliplatin or cisplatin; fluorouracil/irinotecan; paclitaxel or docetaxel/cisplatin) were mainly used in this study. Monotherapy (fluorouracil; capecitabine; paclitaxel or docetaxel) was administered to patients with a poor ECOG PS. The doses of the regimens were decided by physicians according to the actual condition of the patients.
2.6 Chemotherapeutic efficacy assessment
Chemotherapeutic efficacy was evaluated according to the Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1 (27) every 2 cycles or when necessary by 2 independent senior physicians. The optimal therapeutic efficacy was recorded as complete response (CR), partial response (PR), stable disease (SD) (all of which indicated effective treatment), or progressive disease (PD) (which indicated ineffective treatment). The following formulas were applied to determine the objective response rate (ORR) and the disease control rate (DCR): ORR = (CR + PR)/(CR + PR + SD + PD); DCR = (CR + PR + SD)/(CR + PR + SD + PD).
2.7 Follow-up and outcomes
Patients were followed-up by telephone or hospital review until disease progression, death, or loss to follow-up. PFS was calculated from the start of first-line chemotherapy until disease progression or death. Overall survival (OS) was calculated from the start of first-line chemotherapy until death. The outcome measures used in this study were chemotherapeutic efficacy (non-PD or PD), PFS, and OS.
2.8 Statistical analysis
R for Windows (version 4.2.0, https://www.r-project.org/) was used as the primary tool for data analysis and graphing. All statistical tests were two-sided, and P<0.05 was considered to indicate statistical significance. The age data followed a normal distribution and were reported as the means. Other continuous data did not follow a normal distribution and were reported as medians [interquartile range (IQR)]. Continuous variables of two groups of independent samples were compared using Student’s t test (normal distribution) or Wilcoxon test (nonnormal distribution). Continuous variables of more than two groups of independent samples were compared using the Kruskal−Wallis test (nonnormal distribution). Correlations between continuous variables were analyzed using Spearman’s correlation test. Cutoffs for Hb and AGR were determined based on clinical significance. The optimal cutoffs of VEGFA, CEA, SII, LHR, and CAR in therapeutic efficacy analyses were determined according to the maximum Youden’s index. Predictive variables associated with therapeutic efficacy were initially screened using univariate logistic regression analysis. Variables with P<0.05 were included in the multivariate logistic regression analysis, and odds ratios (ORs) and corresponding 95% confidence intervals (CIs) were calculated. A receiver operating characteristic (ROC) curve was used to evaluate the accuracy of the prediction of therapeutic efficacy, and the area under the curve (AUC) was calculated. Optimal cutoffs for VEGFA, CEA, SII, LHR and CAR in survival analyses were determined according to the method of maximally selected rank and statistics. Prognostic variables associated with PFS and OS were initially screened using a univariate Cox regression analysis. Variables with P<0.05 were included in the multivariate Cox regression analysis, and hazard ratios (HRs) and corresponding 95% confidence intervals (CIs) were calculated. Survival analyses were performed using the Kaplan−Meier method, and comparisons of survival between the two groups were performed using the log-rank test.
Least absolute shrinkage and selection operator (LASSO) logistic and Cox regression analyses were performed to screen for candidate variables associated with therapeutic efficacy and prognosis, respectively. The number of events per independent variable (EPV) of each model followed the 10EPV principle as much as possible (28). Four variables were selected to fit the predictive model of therapeutic efficacy using the logistic regression model and five variables were selected to fit the prognostic models of PFS and OS using the Cox regression model. The models’ overall performance was assessed using R2 or Brier scores, and discrimination was evaluated using ROC curves and Harrell’s C-index. The models were internally validated by a bootstrap resampling method (1000 random samples were drawn from the original dataset and replaced between each resample step) (29, 30). Calibration was assessed using a calibration curve or the Hosmer and Lemeshow goodness of fit test statistic. Decision curve analysis (DCA) was performed to estimate the clinical utility of the model predicting therapeutic efficacy by calculating the net benefit for a range of threshold probabilities. Finally, a nomogram and a dynamic nomogram for each model were established.
3 Results
3.1 Relationship between VEGFA expression and either GC development or GC patient prognosis according to the TCGA and Kaplan−Meier Plotter databases
The VEGFA mRNA level was significantly higher in GC tissues (n=415) than in normal gastric tissue (n=35) (p<0.001, Figure 1A) and matched tissues (n=27, p=0.003, Figure 1B) in the TCGA database. Survival analyses according to the mRNA expression levels of VEGFA in GC using the TCGA database showed that patients with high mRNA expression levels of VEGFA had a trend toward shorter OS (n=271, p=0.089, Figure 1C) in all GC populations, and high VEGFA expression was significantly associated with shorter PFS (n=100, HR: 2.08, 95% CI: 1.30-3.33, p=0.002, Figure 1D) and OS (n=102, HR: 1.85, 95% CI: 1.16-2.93, p=0.008, Figure 1E) in patients with HER2-negative GC at stage IV from Kaplan−Meier plotter database.
Figure 1 The mRNA expression level of VEGFA is significantly associated with the development of GC and the prognosis of GC patients. (A) The mRNA expression levels of VEGFA in cancerous tissues versus normal tissues among different types of cancer, including GC (indicated by red box). (B) The mRNA expression levels of VEGFA in cancerous versus corresponding normal tissues of GC patients. (C) Kaplan−Meier survival curves of OS for GC patients from TCGA database with different mRNA expression levels of VEGFA. (D, E) Kaplan−Meier survival curves of PFS and OS of HER2-negative AGC patients from the Kaplan−Meier Plotter database with different mRNA expression levels of VEGFA. **: p<0.01, ***: p<0.001.
3.2 Diagnostic value of baseline serum VEGFA between patients with HER2-negative AGC and nonneoplastic disease
The HER2-negative AGC patients (n=102) in the development cohort included 62 (60.8%) men and 40 (39.2%) women, and 71 (65.7%) men and 37 (34.3%) women were included in the nonneoplastic group (n=108). The mean age did not differ significantly between the two groups (62.7 vs. 59.7 years, p=0.084). The baseline serum VEGFA level was significantly higher in AGC patients than in patients with nonneoplastic disease (178.7 (106.5-226.7) vs. 71.7 (49.3-105.9) pg/ml, p<0.001, Figure 2A). The AUC of the ROC curve indicting the ability of baseline VEGFA to distinguish between HER2-negative AGC and nonneoplastic patients was 0.80 (95% CI: 0.74-0.86, p<0.001, Figure 2B). With a cutoff value of 106.3 pg/ml, VEGFA had a sensitivity of 75.5%, a specificity of 75.9%, a positive likelihood ratio (+LR) of 3.14, and a negative likelihood ratio (−LR) of 0.32.
Figure 2 Baseline serum VEGFA in patients with nonneoplastic disease versus HER2-negative AGC in the development cohort. (A) Baseline serum VEGFA levels were significantly higher in HER2-negative AGC patients (n=102) than in patients with nonneoplastic disease (n=108) (P<0.001). (B) The AUC of the ROC curve indicating the ability of baseline VEGFA to distinguish between patients with nonneoplastic disease and HER2-negative AGC patients was 0.80 (95% CI: 0.74-0.86, P<0.001).
3.3 Differences in baseline serum VEGFA, therapeutic efficacy and survival between HER2-negative and HER2-positive AGC patients
The median follow-up times of the development, validation and HER2-positive cohorts were 13.1 months (95% CI: 10.3-16.6 months), 20.9 months (95% CI: 13.4-NA months) and 18.2 months (95% CI: 13.4-NA months), respectively. There are 6 and 1 patients were lost to OS follow-up in the development and validation cohorts, respectively. The baseline characteristics of the HER2-positive (n=15) and HER2-negative cohorts (development cohort, n=102) and their relationships to baseline serum VEGFA levels are shown in Table S1, Figure 3A and Table 1, respectively. The HER2-negative AGC patients (n=102) in the development cohort included 62 (60.8%) men and 40 (39.2%) women, and 10 (66.7%) men and 5 (33.3%) women were included in the HER2-positive cohort (n=15). The mean age between the two groups did not show a significant difference (62.7 vs. 57.5 years, p=0.126). The baseline serum VEGFA levels were similar between HER2-negative and HER2-positive AGC patients (178.7 (106.5-226.7) vs. 160.5 (102.7-262.8) pg/ml, p=0.549, Figure 3B). The AUC of the ROC curve indicting the ability of baseline VEGFA to distinguish between patients with HER2-negative and HER2-positive AGC was 0.45 (95% CI: 0.29-0.62, p=0.565, Figure 3C). The ORR and DCR were significantly higher in the HER2-positive cohort than in the HER2-negative cohort (ORR: 0.67 vs. 0.23; DCR: 0.87 vs. 0.68). The median PFS (mPFS) and median OS (mOS) were significantly longer in the HER2-positive cohort than in the HER2-negative cohort (mPFS: 9.1 vs. 3.4 months, p=0.024, Figure 3D; mOS: 17.5 vs. 10.0 months, p=0.040, Figure 3E).
Figure 3 Baseline serum VEGFA and survival in HER2-negative AGC versus HER2-positive AGC patients. (A) The baseline characteristics and relationships between baseline serum VEGFA and the clinicopathological characteristics of HER2-positive AGC patients (histology plot is not shown). (B) Baseline serum VEGFA levels were comparable between HER2-negative AGC patients (n=102) and HER2-positive AGC patients (n=15) (P=0.549). (C) The AUC of the ROC curve indicating the ability of baseline VEGFA to distinguish between HER2-negative and HER2-positive AGC patients was 0.45 (95% CI: 0.29-0.62, P=0.565). (D, E) Kaplan−Meier survival curves of PFS and OS of HER2-negative and HER2-positive AGC patients.
Table 1 Relationship between baseline serum VEGFA and clinicopathological characteristics of HER2-negative AGC patients.
3.4 Relationships between baseline serum VEGFA and other clinicopathological characteristics of HER2-negative AGC patients in the development and validation cohorts
Considering that HER2-positive AGC patients show relatively better treatment responses and survival than HER2-negative AGC patients and that these patients represent a small proportion of the GC patient population according to epidemiological studies, we focused on HER2-negative AGC patients for an in-depth analysis. A flowchart of the study is shown in Figure S1A. The baseline characteristics and relationships between baseline serum VEGFA levels and the clinicopathological characteristics of patients in the development (n=102) and validation cohorts (n=34) are shown in Table 1. Regarding clinicopathological characteristics, baseline VEGFA was significantly associated with malignant serous effusion and therapeutic efficacy in the development cohort (both p<0.001, Figure S1B, Table 1) and only therapeutic efficacy in the validation cohort (p<0.001, Table 1). In addition, baseline VEGFA was not significantly correlated with any of the other assessed laboratory indicators in the development cohort (all p>0.05, Figure S1C).
3.5 Relationship between predictive variables and the therapeutic efficacy of first-line chemotherapy in HER2-negative AGC patients
In this study, 33 of 102 patients with therapeutic efficacy data were confirmed to have PD. As shown in Table S2, univariate analysis suggested that BMI (p=0.032), malignant serous effusion (p=0.007), VEGFA (p<0.001), CEA (p=0.031), SII (p=0.044) and CAR (p=0.013) were significantly associated with therapeutic efficacy in the development cohort. Multivariate analysis indicated that only VEGFA [≥180.2 vs. <180.2, odds ratio (OR): 6.95, 95% CI: 2.30-20.99, p=0.001] was an independent predictive factor for therapeutic efficacy in HER2-negative AGC patients (Table S2, Figure 4A).
Figure 4 Predictive variable screening and construction, evaluation and presentation of a predictive model of therapeutic efficacy in HER2-negative AGC patients. (A) Forest plot showing the multivariate logistic analysis of the predictive variables that were significantly associated with therapeutic efficacy in the univariate logistic analyses. (B, C) Predictive variables were screened using the LASSO regression model via 10-fold cross-validation based on minimum criteria. The dotted vertical line in C was plotted at the value selected using 10-fold cross-validation in C. (D) ROC curves of the predictive model for predicting the therapeutic efficacy of HER2-negative AGC patients in the development and validation cohorts. (E) Calibration plot of the predictive model for HER2-negative AGC in the development cohort. (F) Calibration plot of the predictive model for HER2-negative AGC in the validation cohort. (G) The DCA for the predictive model and VEGFA alone in the development cohort. (H) Nomogram to predict the probability of PD at first evaluation in patients with HER2-negative AGC.
3.6 Development of nomograms for predicting therapeutic efficacy in HER2-negative AGC patients with first-line chemotherapy
In all 102 HER2-negative AGC patients with therapeutic efficacy data, variables were screened according to the LASSO logistic regression model (Figures 4B, C), clinical significance, and recommendations from oncology experts, and BMI, VEGFA, CEA and CAR were selected to fit a logistic regression model for predicting the therapeutic efficacy of PD at first evaluation. The model had an R2 of 0.37 and a Brier score of 0.15. The discriminatory efficiency of the model was evaluated by ROC analysis, which indicated that the model achieved a Harrell’s C-index of 0.82 (95% CI: 0.73-0.91) before calibration and 0.79 after calibration with internal validation using bootstrap resampling (1,000 repetitions) in the development cohort and 0.90 (95% CI: 0.79-1.00) in the validation cohort (Figure 4D). The calibration curves of the model in the development and validation cohorts are illustrated in Figures 4E, F; the Hosmer−Lemeshow test statistics indicated nonsignificant differences between sample types in the development cohort (p=0.370). The DCA for the model is presented in Figure 4G, which indicated that the model adds more net benefit than VEGFA alone and the “treat all” or “treat none” strategies. Finally, the nomogram and web-based dynamic nomogram of the therapeutic efficacy predictive model are presented in Figure 4H and https://cancerprediction.shinyapps.io/efficacy_prediction_of_her2-agc/. The calculation formula for the model is presented on the webpage of the dynamic nomogram.
3.7 Relationship between predictive variables and PFS after first-line treatment in HER2-negative AGC patients
In this study, 83 of 102 patients with complete follow-up data reached the PFS endpoint. As shown in Table S3, univariate analysis suggested that the diagnosis pattern (p=0.033), histology (p=0.004), malignant serous effusion (p=0.003), VEGFA (p<0.001), SII (p=0.006) and CAR (p=0.045) were significantly associated with PFS in patients with HER2-negative AGC. Multivariate analysis indicated that diagnosis pattern (recurrent vs. primary, HR: 2.34, 95% CI: 1.43-3.85, p=0.001), histology (signet-ring cell carcinoma vs. adenocarcinoma, HR: 1.88, 95% CI: 1.07-3.3, p=0.027), VEGFA (≥179.1 vs. <179.1, HR: 2.49, 95% CI: 1.51-4.11, p<0.001) and SII (≥376.9 vs. <376.9, HR: 1.98, 95% CI: 1.10-3.56, p=0.023) were independent prognostic factors for PFS in HER2-negative AGC patients (Table S3, Figure 5A). Kaplan−Meier plots of patients stratified by each independent prognostic factor are shown in Figures 6A–D.
Figure 5 Predictive variable screening and construction, evaluation and presentation of the PFS predictive model in HER2-negative AGC patients. (A) Forest plot showing the multivariate Cox analysis of PFS with the predictive variables that were significantly associated with PFS in the univariate Cox analyses. (B, C) Predictive variables were screened using the LASSO regression model via 10-fold cross-validation based on minimum criteria. The dotted vertical line in C was plotted at the value selected using 10-fold cross-validation in (C). (D) Kaplan−Meier survival curves of PFS of HER2-negative AGC patients in the development and validation cohorts. (E) Time-dependent ROC curves of the PFS predictive model at 3 and 6 months for HER2-negative AGC in the development cohort. (F) Time-dependent ROC curves of the PFS predictive model at 3 and 6 months for HER2-negative AGC in the validation cohort. (G, H) Calibration plot of the PFS predictive model at 3 and 6 months for HER2-negative AGC in the development cohort. (I, J) Calibration plot of the PFS predictive model at 3 and 6 months for HER2-negative AGC in the validation cohort. (K) Baseline survival curves of PFS of HER2-negative AGC patients in the development and validation cohorts. (L) Nomogram to predict the probability of 3-month or 6-month PFS in patients with HER2-negative AGC.
Figure 6 Kaplan−Meier plots of patients stratified by each independent prognostic factor of first-line PFS and OS. (A–D) Kaplan−Meier plots of patients stratified by each independent prognostic factor of first-line PFS. (E, F) Kaplan−Meier plots of patients stratified by each independent prognostic factor of OS.
3.8 Development of nomograms for predicting PFS after first-line treatment in HER2-negative AGC patients
In all 102 HER2-negative AGC patients with complete follow-up data, variables were screened according to the LASSO Cox regression model (Figures 5B, C), clinical significance, and recommendations from oncology experts, and diagnosis pattern, histology, malignant serous effusion, VEGFA and SII were selected to fit a Cox regression model for predicting PFS that met the proportional hazards assumption. The mPFS was similar between the development and validation cohorts (p=0.730, Figure 5D). The model had Harrell’s C-indexes of 0.71 and 0.70 after calibration with internal validation using bootstrap resampling (1,000 repetitions) in the development cohort and 0.71 in the validation cohort. The time-dependent ROC curves and calibration curves of the 3-month and 6-month PFS prediction in the development (AUC: 0.85 and 0.75, respectively) and validation cohorts (AUC: 0.73 and 0.79, respectively) are illustrated in Figures 5E–J. The baseline survival curves of PFS were similar in both development and validation cohorts (Figure 5K). The nomogram and web-based dynamic nomogram of the predictive PFS model are presented in Figure 5L and https://cancerprediction.shinyapps.io/pfs_prediction_of_her2-agc/. The calculation formula for the PFS model is presented on the webpage of the dynamic nomogram.
3.9 Relationship between predictive variables and OS of HER2-negative AGC patients
In this study, 57 of 102 patients without loss to follow-up reached the OS endpoint. As shown in Table S4, univariate analysis suggested that histology (p<0.001), malignant serous effusion (p=0.007), VEGFA (p=0.008), SII (p=0.003) and CAR (p=0.010) were significantly associated with OS in patients with HER2-negative AGC. Multivariate analysis indicated that histology (signet-ring cell carcinoma vs. adenocarcinoma, HR: 2.93, 95% CI: 1.50-5.74, p=0.002) and SII (≥448.1 vs. <448.1, HR: 2.06, 95% CI: 1.03-4.12, p=0.040) were independent prognostic factors for OS in patients with HER2-negative AGC (Table S4; Figure 7A). Kaplan−Meier plots of patients stratified by each independent prognostic factor are shown in Figures 6E, F.
Figure 7 Predictive variable screening and construction, evaluation and presentation of the OS predictive model in HER2-negative AGC patients. (A) Forest plot showing the multivariate Cox analysis of OS with the predictive variables that were significantly associated with OS in the univariate Cox analyses. (B, C) Predictive variables were screened using the LASSO regression model via 10-fold cross-validation based on minimum criteria. The dotted vertical line in C was plotted at the value selected using 10-fold cross-validation in C. (D) Kaplan−Meier survival curves of OS of HER2-negative AGC patients in the development and validation cohorts. (E) Time-dependent ROC curves of the OS predictive model at 9 and 12 months for HER2-negative AGC in the development cohort. (F) Time-dependent ROC curves of the OS predictive model at 9 and 12 months for HER2-negative AGC in the validation cohort. (G, H) Calibration plot of the OS predictive model at 9 and 12 months for HER2-negative AGC in the development cohort. (I, J) Calibration plot of the OS predictive model at 9 and 12 months for HER2-negative AGC in the validation cohort. (K) Baseline survival curves of OS of HER2-negative AGC patients in the development and validation cohorts. (L) Nomogram to predict the probability of 9-month or 12-month OS in patients with HER2-negative AGC.
3.10 Development of nomograms for predicting the OS of HER2-negative AGC patients
In 102 HER2-negative AGC patients with complete follow-up data, variables were screened according to the LASSO Cox regression model (Figures 7B, C), clinical significance, and recommendations from oncology experts, and histology, malignant serous effusion, VEGFA, SII and CAR were selected to fit a Cox regression model for predicting OS that met the proportional hazards assumption. The mOS was similar between the development and validation cohorts (p=0.760, Figure 7D). The model had Harrell’s C-indexes of 0.69 and 0.66 after calibration with internal validation using bootstrap resampling (1,000 repetitions) in the development cohort and 0.62 in the validation cohort. The time-dependent ROC curves and calibration curves of the 9-month and 12-month OS prediction in the development (AUC: 0.69 and 0.78, respectively) and validation cohorts (AUC: 0.62 and 0.59, respectively) are illustrated in Figures 7E–J. The baseline survival curves of OS were similar in both development and validation cohorts (Figure 7K). The nomogram and web-based dynamic nomogram of the predictive OS model are presented in Figure 7L and https://cancerprediction.shinyapps.io/os_prediction_of_her2-agc/. The calculation formula for the OS model is presented on the webpage of the dynamic nomogram.
4 Discussion
HER2-negative GC, which lacks precision medical care, is the main type of GC and has a poor prognosis. Establishment of predictive models can help effectively stratify populations and optimize individualized chronic disease management strategies for this malignancy. In this study, the role of immune-inflammatory biomarkers including serum VEGFA and SII in the treatment response and survival of HER2-negative AGC patients was clarified by a cohort study, and predictive nomograms related to therapeutic efficacy and prognosis prediction were established. Internal and external validation has shown that these nomograms based on immune-inflammatory indicators perform well in predicting the treatment response and survival of HER2-negative AGC patients and will provide an important reference for the clinical management of this lethal disease.
VEGF is a family of pivotal growth factors and signaling molecules involved in angiogenesis. Binding to its receptor, VEGFR, activates a complex cascade of downstream signaling pathways that leads to neovascularization and vasodilation (31). Inhibition of VEGF signaling activity impairs these pathways, resulting in a reduction in tumor proliferation and invasion in GC (32). As a major member of the VEGF family, VEGFA has been reported to be associated with prognosis in various malignancies (10, 11) and is considered to be a useful biomarker of the progression and remission of diseases, including GC (33). Clinically, ramucirumab is the primary antiangiogenic agent used for AGC. Ramucirumab was confirmed as the preferred agent for second-line and beyond therapy in AGC by the REGARD (34) and RAINBOW (35) studies, but the failure of the RAINFALL (36) study prevented ramucirumab from entering first-line therapy. One reason for these findings may be that the combination of ramucirumab/fluoropyrimidine/cisplatin may not be the optimal regimen; however, the lack of a screened population of patients predicted to benefit from antiangiogenic therapy may also contribute to the failure of the study. To investigate the potential role of serum VEGFA in screening the population that may benefit from chemotherapy, we previously launched a prospective cohort study in highly malignant SCLC and revealed the unique ability of VEGFA to predict first-line chemotherapeutic efficacy and survival in patients with this kind of disease (37). In this previous study, we also found that the median VEGFA level decreased after 2 cycles of chemotherapy and returned to pretreatment levels at progression, but the changes in VEGFA were not associated with therapeutic efficacy and PFS (37). Therefore, we conducted a follow-up study of AGC patients and focused on investigating the role of VEGFA at baseline. The results showed that VEGFA was not only an independent predictor of first-line chemotherapy efficacy and corresponding PFS outcomes but also held the greatest weight in the therapeutic efficacy and PFS models. Despite its reduced value for OS prediction, it still plays a role in the OS predictive model. PFS reflects the survival benefit of stage-specific antitumor therapies. Since patients with AGC in China often receive 3 or more lines of antitumor therapy (38), including late-line chemotherapy and antiangiogenic targeted therapy, the benefit in terms of PFS after first-line therapy does not reliably reflect the benefit in terms of OS, which could be why baseline VEGFA did not predict OS in this study. Therefore, for HER2-negative AGC, VEGFA can be used as a biomarker for evaluating efficacy and short-term survival but is a weak predictor of long-term survival. VEGFA-based nomograms would help clinicians screen chemotherapy-resistant populations in advance and provide a basis for stratifying patients for future intensive therapy studies.
It is noted that the present study identified SII as an independent prognostic factor for PFS and OS after first-line chemotherapy in patients with HER2-negative AGC, and a model including the SII was effective in predicting patient OS. The SII considers the abundance of lymphocytes that play a major role in inducing cytotoxic cell death and suppressing cancer cell proliferation (39), neutrophils that are known to secrete a variety of cytokines (VEGF, IL-8, IL-16, etc.) that stimulate tumor cell growth and metastasis (40), and platelets that can promote cancer cell arrest at the endothelium and support the establishment of secondary lesions of cancer cells (41). Therefore, the elevation of the SII may reflect cancer-promoting activity in the tumor microenvironment (TME), which eventually leads to poor prognosis (42, 43). The predictive effect of the SII on PFS and OS in the present study supports the notion that the immune inflammatory response broadly affects antitumor therapy and ultimately affects the survival of patients both in the short and long term. Currently, immune checkpoint inhibition has become a new standard treatment of AGC. Results from KEYNOTE-063 (44) and JAVELIN Gastric 100 (45) studies have confirmed that ICIs monotherapies failed to outperform the benefit from chemotherapy in patients with AGC. In contrast, ICIs combined with chemotherapy explored in CheckMate-649 (46) and ORIENT-16 (47) studies significantly improved the response rate and survival and has become the first-line standard treatment for this kind of patients, which implied that chemotherapy remains the cornerstone of systemic treatment for this fatal disease. Chemotherapy has been validated as potent immunogenic cell death (ICD) inducer which can destruct cancer cells and help exposing immunostimulatory molecules like damage-associated molecular patterns (DAMPs) and cytokines like type I interferons (IFNs), and thus increase the sensitivity of tumor to ICIs (48). In addition, in our study focus on patients with chemotherapy, VEGFA, SII and CAR in the predictive models are also immune-inflammatory indicators. Therefore, the determination of the sensitivity of chemotherapy with our models may contribute to the judgment of the efficacy of immunotherapy and this needs to be validated in future study on patients receiving ICIs.
Most importantly, novel immune-inflammation-based clinical predictive nomograms in HER2-negative AGC were developed and validated in this study. Most of the previous prediction models were based on retrospective studies and mainly focused on diagnosis (49) or prediction of neoadjuvant chemotherapy response (50) or survival in GC (51). Models in the present study were based on a prospective cohort study with more reliable information and less recall bias than previous studies. Validation based on an independent cohort enhanced the generalizability of the study results. Moreover, this study not only established a satisfactory prediction model for OS but also innovatively developed predictive models for first-line chemotherapy efficacy and PFS, providing new approaches for the screening of primary chemotherapy-resistant populations. In addition, compared with previous models, this study added new blood indicators, such as VEGFA, SII, and other traditional markers, making the results more comprehensive, reliable, and clinically applicable.
This study had some limitations. As this was a single-center study, the sample size was relatively small, which limited the number of variables that could be included in the models, leading to a decrease in their predictive power. In addition, although we enrolled HER2-negative AGC patients for external validation, data from other institutions would allow a more robust evaluation of the universality of the predictive models, and this will be complemented by a subsequent validation study with multicenter data.
In summary, this study established novel and effective series of immune-inflammation-based predictive models that showed favorable performance in evaluating first-line chemotherapeutic efficacy and survival in HER2-negative AGC patients. Nomograms constructed using VEGFA, SII, and other accepted clinicopathological variables could be utilized as an efficient computational technique for predicting the clinical prognosis of this kind of patients. We considered that angiogenesis and immune inflammation are interlinked in the treatment of GC and collectively affect the prognosis of patients with AGC. In the future, the use of reliable predictive models to screen patients who may benefit or not benefit from specific treatment may lead to the broad implementation of individualized medicine.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Institutional Review Board of the First Affiliated Hospital of Bengbu Medical College. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Conceptualization, JY. Methodology, YY and YS. Software, YY and YS. Formal Analysis, YY, YS and HY. Investigation, JW, QC, YL, JL and YZ. Resources, ZZ, MW. Data Curation, JW. Writing-Original Draft Preparation, YY and YS. Writing-Review and Editing, BJ, JY. Supervision, JY, YY. Project Administration, JY, YY. Funding Acquisition, YY. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Natural Science Foundation of Anhui Province (No. 2008085MH238), the Excellent Young Talents Fund Program of Higher Education Institutions of Anhui Province (No. gxyq2022042), the 512 Talent Cultivation Plan of Bengbu Medical College (No. by51202208), and the Science Fund for Distinguished Young Scholars of the First Affiliated Hospital of Bengbu Medical College (No. 2019byyfyjq02).
Acknowledgments
We sincerely thank all the patients and medical staffs who participated in this 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.
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.2023.1185240/full#supplementary-material
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
2. Li K, Zhang A, Li X, Zhang H, Zhao L. Advances in clinical immunotherapy for gastric cancer. Biochim Biophys Acta Rev Cancer. (2021) 1876(2):188615. doi: 10.1016/j.bbcan.2021.188615
3. Joshi SS, Badgwell BD. Current treatment and recent progress in gastric cancer. CA Cancer J Clin (2021) 71(3):264–79. doi: 10.3322/caac.21657
4. Smyth EC, Nilsson M, Grabsch HI, van Grieken NC, Lordick F. Gastric cancer. Lancet (2020) 396(10251):635–48. doi: 10.1016/S0140-6736(20)31288-5
5. Wagner AD, Syn NL, Moehler M, Grothe W, Yong WP, Tai BC, et al. Chemotherapy for advanced gastric cancer. Cochrane Database Syst Rev (2017) 8:CD004064. doi: 10.1002/14651858.CD004064.pub4
6. Kiely M, Lord B, Ambs S. Immune response and inflammation in cancer health disparities. Trends Cancer. (2022) 8(4):316–27. doi: 10.1016/j.trecan.2021.11.010
7. Hu X, Chen Z, Wang Z, Xiao Q. Cancer evolution: Special focus on the immune aspect of cancer. Semin Cancer Biol (2022) 86(Pt 2):420–435. doi: 10.1016/j.semcancer.2022.05.006
8. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell (2011) 144(5):646–74. doi: 10.1016/j.cell.2011.02.013
9. Llovet JM. Focal gains of VEGFA: candidate predictors of sorafenib response in hepatocellular carcinoma. Cancer Cell (2014) 25(5):560–2. doi: 10.1016/j.ccr.2014.04.019
10. Dinami R, Porru M, Amoreo CA, Sperduti I, Mottolese M, Buglioni S, et al. TRF2 and VEGF-A: an unknown relationship with prognostic impact on survival of colorectal cancer patients. J Exp Clin Cancer Res (2020) 39(1):111. doi: 10.1186/s13046-020-01612-z
11. Feng H, Jin Z, Liang J, Zhao Q, Zhan L, Yang Z, et al. FOXK2 transcriptionally activating VEGFA induces apatinib resistance in anaplastic thyroid cancer through VEGFA/VEGFR1 pathway. Oncogene (2021) 40(42):6115–29. doi: 10.1038/s41388-021-01830-5
12. Hu B, Yang XR, Xu Y, Sun YF, Sun C, Guo W, et al. Systemic immune-inflammation index predicts prognosis of patients after curative resection for hepatocellular carcinoma. Clin Cancer Res (2014) 20(23):6212–22. doi: 10.1158/1078-0432.CCR-14-0442
13. Xu BB, Lu J, Zheng ZF, Xie JW, Wang JB, Lin JX, et al. The predictive value of the preoperative C-reactive protein-albumin ratio for early recurrence and chemotherapy benefit in patients with gastric cancer after radical gastrectomy: using randomized phase III trial data. Gastric Cancer. (2019) 22(5):1016–28. doi: 10.1007/s10120-019-00936-w
14. Shepshelovich D, Xu W, Lu L, Fares A, Yang P, Christiani D, et al. Body mass index (BMI), BMI change, and overall survival in patients with SCLC and NSCLC: A pooled analysis of the international lung cancer consortium. J Thorac Oncol (2019) 14(9):1594–607. doi: 10.1016/j.jtho.2019.05.031
15. Gou M, Zhang Y, Liu T, Qu T, Si H, Wang Z, et al. The prognostic value of pre-treatment hemoglobin (Hb) in patients with advanced or metastatic gastric cancer treated with immunotherapy. Front Oncol (2021) 11:655716. doi: 10.3389/fonc.2021.655716
16. Chi J, Xie Q, Jia J, Liu X, Sun J, Chen J, et al. Prognostic value of albumin/globulin ratio in survival and lymph node metastasis in patients with cancer: A systematic review and meta-analysis. J Cancer. (2018) 9(13):2341–8. doi: 10.7150/jca.24889
17. Pirro M, Ricciuti B, Rader DJ, Catapano AL, Sahebkar A, Banach M. High density lipoprotein cholesterol and cancer: Marker or causative? Prog Lipid Res (2018) 71:54–69. doi: 10.1016/j.plipres.2018.06.001
18. Kochi M, Fujii M, Kanamori N, Kaiga T, Kawakami T, Aizaki K, et al. Evaluation of serum CEA and CA19-9 levels as prognostic factors in patients with gastric cancer. Gastric Cancer. (2000) 3(4):177–86. doi: 10.1007/pl00011715
19. Tang D, Zhou J, Wang L, Ni M, Chen M, Hassan S, et al. A novel model based on deep convolutional neural network improves diagnostic accuracy of intramucosal gastric cancer (with video). Front Oncol (2021) 11:622827. doi: 10.3389/fonc.2021.622827
20. Chen D, Liu Z, Liu W, Fu M, Jiang W, Xu S, et al. Predicting postoperative peritoneal metastasis in gastric cancer with serosal invasion using a collagen nomogram. Nat Commun (2021) 12(1):179. doi: 10.1038/s41467-020-20429-0
21. Chen L, Ma G, Wang P, Dong Y, Liu Y, Zhao Z, et al. Establishment and verification of prognostic model for gastric cancer based on autophagy-related genes. Am J Cancer Res (2021) 11(4):1335–46.
22. Li J, Pu K, Li C, Wang Y, Zhou Y. A novel six-gene-based prognostic model predicts survival and clinical risk score for gastric cancer. Front Genet (2021) 12:615834. doi: 10.3389/fgene.2021.615834
23. Cuschieri S. The STROBE guidelines. Saudi J Anaesth. (2019) 13(Suppl 1):S31–4. doi: 10.4103/sja.SJA_543_18
24. Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ (2015) 350:g7594. doi: 10.1136/bmj.g7594
25. Szasz AM, Lanczky A, Nagy A, Forster S, Hark K, Green JE, et al. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients. Oncotarget (2016) 7(31):49322–33. doi: 10.18632/oncotarget.10337
26. Hohn A, Grune T. Lipofuscin: formation, effects and role of macroautophagy. Redox Biol (2013) 1:140–4. doi: 10.1016/j.redox.2013.01.006
27. Eisenhauer EA, Therasse P, Bogaerts J, Schwartz LH, Sargent D, Ford R, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Cancer. (2009) 45(2):228–47. doi: 10.1016/j.ejca.2008.10.026
28. Peduzzi P, Concato J, Kemper E, Holford TR, Feinstein AR. A simulation study of the number of events per variable in logistic regression analysis. J Clin Epidemiol. (1996) 49(12):1373–9. doi: 10.1016/s0895-4356(96)00236-3
29. Steyerberg EW, Bleeker SE, Moll HA, Grobbee DE, Moons KG. Internal and external validation of predictive models: a simulation study of bias and precision in small samples. J Clin Epidemiol. (2003) 56(5):441–7. doi: 10.1016/s0895-4356(03)00047-7
30. Steyerberg EW, Harrell FE Jr., Borsboom GJ, Eijkemans MJ, Vergouwe Y, Habbema JD. Internal validation of predictive models: efficiency of some procedures for logistic regression analysis. J Clin Epidemiol. (2001) 54(8):774–81. doi: 10.1016/s0895-4356(01)00341-9
31. Yang Y, Cao Y. The impact of VEGF on cancer metastasis and systemic disease. Semin Cancer Biol (2022) 86(Pt 3):251–261. doi: 10.1016/j.semcancer.2022.03.011
32. Park DJ, Thomas NJ, Yoon C, Yoon SS. Vascular endothelial growth factor a inhibition in gastric cancer. Gastric Cancer. (2015) 18(1):33–42. doi: 10.1007/s10120-014-0397-4
33. Macedo F, Ladeira K, Longatto-Filho A, Martins SF. Gastric cancer and angiogenesis: is VEGF a useful biomarker to assess progression and remission? J Gastric Cancer (2017) 17(1):1–10. doi: 10.5230/jgc.2017.17.e1
34. Fuchs CS, Tomasek J, Yong CJ, Dumitru F, Passalacqua R, Goswami C, et al. Ramucirumab monotherapy for previously treated advanced gastric or gastro-oesophageal junction adenocarcinoma (REGARD): an international, randomised, multicentre, placebo-controlled, phase 3 trial. Lancet (2014) 383(9911):31–9. doi: 10.1016/S0140-6736(13)61719-5
35. Wilke H, Muro K, Van Cutsem E, Oh SC, Bodoky G, Shimada Y, et al. Ramucirumab plus paclitaxel versus placebo plus paclitaxel in patients with previously treated advanced gastric or gastro-oesophageal junction adenocarcinoma (RAINBOW): a double-blind, randomised phase 3 trial. Lancet Oncol (2014) 15(11):1224–35. doi: 10.1016/S1470-2045(14)70420-6
36. Fuchs CS, Shitara K, Di Bartolomeo M, Lonardi S, Al-Batran SE, Van Cutsem E, et al. Ramucirumab with cisplatin and fluoropyrimidine as first-line therapy in patients with metastatic gastric or junctional adenocarcinoma (RAINFALL): a double-blind, randomised, placebo-controlled, phase 3 trial. Lancet Oncol (2019) 20(3):420–35. doi: 10.1016/S1470-2045(18)30791-5
37. Yao JH, Shao Y, Wang JJ, Li YL, Yang HQ, Liu J, et al. Evaluation of diagnostic and predictive values of the serum VEGF-A level and systemic immune-inflammation index in small cell lung cancer. J Cancer. (2021) 12(5):1356–64. doi: 10.7150/jca.51972
38. Qi CS, Cheng SY, Shen L. Current status and research progress of third-line treatment for patients with gastric cancer in China. Zhonghua Zhong Liu Za Zhi. (2020) 42(12):983–8. doi: 10.3760/cma.j.cn112152-20190910-00589
39. Loi S, Michiels S, Adams S, Loibl S, Budczies J, Denkert C, et al. The journey of tumor-infiltrating lymphocytes as a biomarker in breast cancer: clinical utility in an era of checkpoint inhibition. Ann Oncol (2021) 32(10):1236–44. doi: 10.1016/j.annonc.2021.07.007
40. Xiong S, Dong L, Cheng L. Neutrophils in cancer carcinogenesis and metastasis. J Hematol Oncol (2021) 14(1):173. doi: 10.1186/s13045-021-01187-y
41. Lazar S, Goldfinger LE. Platelets and extracellular vesicles and their cross talk with cancer. Blood (2021) 137(23):3192–200. doi: 10.1182/blood.2019004119
42. Aziz MH, Sideras K, Aziz NA, Mauff K, Haen R, Roos D, et al. The systemic-immune-inflammation index independently predicts survival and recurrence in resectable pancreatic cancer and its prognostic value depends on bilirubin levels: a retrospective multicenter cohort study. Ann Surg (2019) 270(1):139–46. doi: 10.1097/SLA.0000000000002660
43. Fest J, Ruiter R, Mulder M, Groot Koerkamp B, Ikram MA, Stricker BH, et al. The systemic immune-inflammation index is associated with an increased risk of incident cancer-A population-based cohort study. Int J Cancer. (2020) 146(3):692–8. doi: 10.1002/ijc.32303
44. Chung HC, Kang YK, Chen Z, Bai Y, Wan Ishak WZ, Shim BY, et al. Pembrolizumab versus paclitaxel for previously treated advanced gastric or gastroesophageal junction cancer (KEYNOTE-063): a randomized, open-label, phase 3 trial in Asian patients. Cancer (2022) 128(5):995–1003. doi: 10.1002/cncr.34019
45. Moehler M, Dvorkin M, Boku N, Ozguroglu M, Ryu MH, Muntean AS, et al. Phase III trial of avelumab maintenance after first-line induction chemotherapy versus continuation of chemotherapy in patients with gastric cancers: results from JAVELIN gastric 100. J Clin Oncol (2021) 39(9):966–77. doi: 10.1200/JCO.20.00892
46. Janjigian YY, Shitara K, Moehler M, Garrido M, Salman P, Shen L, et 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(10294):27–40. doi: 10.1016/S0140-6736(21)00797-2
47. Xu J, Jiang H, Pan Y, Gu K, Cang S, Han L, et al. Abstract CT078: First-line treatment with sintilimab (sin) vs placebo in combination with chemotherapy (chemo) in patients (pts) with unresectable gastric or gastroesophageal junction (G/GEJ) cancer: Final overall survival (OS) results from the randomized, phase III ORIENT-16 trial. Cancer Res (2023) 83(8_Supplement):CT078–8. doi: 10.1158/1538-7445.Am2023-ct078
48. Sprooten J, Laureano RS, Vanmeerbeek I, Govaerts J, Naulaerts S, Borras DM, et al. Trial watch: chemotherapy-induced immunogenic cell death in oncology. Oncoimmunology (2023) 12(1):2219591. doi: 10.1080/2162402X.2023.2219591
49. Zheng X, Wang R, Zhang X, Sun Y, Zhang H, Zhao Z, et al. A deep learning model and human-machine fusion for prediction of EBV-associated gastric cancer from histopathology. Nat Commun (2022) 13(1):2790. doi: 10.1038/s41467-022-30459-5
50. Nakauchi M, Court CM, Tang LH, Gonen M, Janjigian YY, Maron SB, et al. Validation of the memorial sloan kettering gastric cancer post-resection survival nomogram: does it stand the test of time? J Am Coll Surg (2022) 235(2):294–304. doi: 10.1097/XCS.0000000000000251
51. Sundar R, Barr Kumarakulasinghe N, Huak Chan Y, Yoshida K, Yoshikawa T, Miyagi Y, et al. Machine-learning model derived gene signature predictive of paclitaxel survival benefit in gastric cancer: results from the randomised phase III SAMIT trial. Gut (2021) 71(4):676–685. doi: 10.1136/gutjnl-2021-324060
Keywords: nomogram, predictive model, VEGFA, HER2-negative, advanced gastric cancer
Citation: Yang Y, Shao Y, Wang J, Cheng Q, Yang H, Li Y, Liu J, Zhou Y, Zhou Z, Wang M, Ji B and Yao J (2023) Development and validation of novel immune-inflammation-based clinical predictive nomograms in HER2-negative advanced gastric cancer. Front. Oncol. 13:1185240. doi: 10.3389/fonc.2023.1185240
Received: 15 March 2023; Accepted: 07 August 2023;
Published: 06 September 2023.
Edited by:
Dong Ren, UC Irvine Medical Center, United StatesReviewed by:
Jianhua Wan, Nanchang University, ChinaJun Xue, Huazhong University of Science and Technology, China
Copyright © 2023 Yang, Shao, Wang, Cheng, Yang, Li, Liu, Zhou, Zhou, Wang, Ji and Yao. 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: Jinghao Yao, bGVpaml5YXNpQG91dGxvb2suY29t
†These authors have contributed equally to this work and share first authorship