- 1Department of Radiology, Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education), Peking University Cancer Hospital and Institute, Beijing, China
- 2Center for Data Science, Peking University, Beijing, China
- 3Department of Gastrointestinal Oncology, Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education), Peking University Cancer Hospital and Institute, Beijing, China
- 4National Engineering Laboratory for Big Data Analysis and Applications, Peking University, Beijing, China
- 5Beijing International Center for Mathematical Research (BICMR), Peking University, Beijing, China
Introduction: To explore the prognostic value of CT-based delta radiomics in predicting the prognosis of patients with stage IV gastric cancer treated with immune checkpoint inhibitors (ICI).
Materials and methods: Forty-two patients with stage IV gastric cancer, who had received ICI monotherapy, were enrolled in this retrospective study. Baseline and first follow-up CT scans were analyzed. Intratumoral and peritumoral regions of interest (ROI) were contoured, enabling the extraction of 192 features from each ROI. The intraclass correlation coefficients were used to select features with high stability. The least absolute shrinkage and selection operator was used to select features with high weights for predicting patient prognosis. Kaplan–Meier analysis and log-rank test were performed to explore the association between features and progression free survival (PFS). Cox regression analyses were used to identify predictors for PFS. The C-index was used to assess the prediction performance of features.
Results: Two radiomics features of ΔVintra_ZV and postVperi_Sphericity were identified from intratumoral and peritumoral regions, respectively. The Kaplan–Meier analysis revealed significant differences in PFS between patients with low and high feature value (ΔVintra_ZV: P=0.000; postVperi_Sphericity: P=0.012), and the multivariable cox analysis demonstrated that ΔVintra_ZV was independent predictor for PFS (HR, 1.911; 95% CI: 1.163–3.142; P=0.011), with C-index of 0.705.
Conclusions: Based on CT scans at baseline and first follow-up, the delta radiomics features could efficiently predict the PFS of gastric cancer patients treated with ICI therapy.
Introduction
Gastric cancer is one of the malignancies with high mortality rate (1). Despite significant efforts to develop innovative treatment techniques based on cytotoxic chemotherapy, targeted therapy, and radiotherapy, a significant proportion of gastric cancer patients will still demonstrate poor response to conventional therapies or even fast progression after treatment (2). Immune checkpoint inhibitors (ICIs) have revolutionized the treatment of a variety of malignancies, including gastric cancer (3). More specifically, several large multicenter clinical trials demonstrated a significant and durable survival benefit in refractory gastric cancer patients who received ICI therapy, with a duration of response ranging from 8.4 to 9.5 months (3, 4). However, treatment response varied significantly as 60% of patients derived no benefit from ICI therapy, and 21% of patients even showed hyperprogression during treatment (2, 3, 5). Therefore, there is an urgent need for the introduction of precise biomarkers that can predict the response of gastric cancer to ICI at the early treatment stage.
Several predictive tumor biomarkers from biopsy tissue samples could indicate ICI treatment response and prognosis of gastric cancer patients (6), such as positivity of programmed death-ligand 1 (PD-L1), mismatch repair deficiency (dMMR), and Epstein–Barr virus (EBV) (6, 7). However, most patients tested negative for the above-mentioned biomarkers, e.g., 86% of patients with PD-L1 combined positive scores (CPS)<1, 78.4-92.5% of patients with mismatch repair proficiency (pMMR), and 91% of patients with negative EBV according to previous studies (3, 5, 8, 9). Some of these biomarkers predicted that patients with poor treatment response could still respond well to ICI therapy, with reported objective response rates of approximately 6.4-10.9% for PD-L1 CPS < 1, 12.3% for pMMR, and 16.4% for negative EBV, respectively (3, 6, 7). Besides, not all laboratories have the available resources to perform complex immunohistochemistry protocols that are necessary to identify or evaluate potential tumor biomarkers, hindering their subsequent application in clinical practice (10, 11). In addition, given the spatial heterogeneity of gastric cancer, biopsy samples may not always be evaluated appropriately.
Computed tomography (CT) has been widely and routinely used in clinical practice, yet traditional unidimensional measurements made both RECIST and iRECIST criteria no longer meet the needs of the ICI response evaluation and hindered the realization of the precision medicine (12). Radiomics is a useful tool to mine data from radiographic images, such as tumor texture characteristics, which may not be detectable by ‘naked-eye’ inspection (13). Several studies have verified that radiomics could predict response to neoadjuvant chemotherapy and palliative chemotherapy in patients with gastric cancer, with an area under the curve (AUC) of 0.74-0.82 (14–16). Recently, one study explored the response prediction performance of baseline CT radiomics in patients treated with immunotherapy combined with chemotherapy and showed promising results, with an AUC over 0.7 (17). To the best of our knowledge, the prognostic value of radiomics features in patients with gastric cancer treated with ICI monotherapy has not been elucidated. Therefore, this study aimed to use delta radiomics to extract information from CT scans (at baseline and first follow-up) and predict the survival of patients with stage IV gastric cancer treated by ICI.
Materials and methods
Patients
This study was performed in line with the principles of the Declaration of Helsinki. Written informed consent of this retrospective study was waived. Data from 101 consecutive patients with stage IV gastric cancer who had received anti-programmed cell death protein 1/programmed cell death ligand 1 (PD-1/PD-L1) antibody alone or in combination with anti-cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) antibody were collected in the Peking University Cancer Hospital, Beijing, China, between 2016 and 2020. Inclusion criteria were as follows: (a) histologically confirmed gastric adenocarcinoma; (b) patients treated with ICI monotherapy (anti-PD-1/PD-L1 alone or in combination with anti-CTLA-4 antibodies); (c) availability of baseline enhanced abdominal/pelvic CT scans performed <30 days before ICI treatment; (d) availability of the first follow-up enhanced abdominal/pelvic CT scans two to three cycles after ICI treatment initiation. Exclusion criteria were as follows: (a) patients with primary gastric surgical treatment (n=53) (b) patients with other synchronous or metachronous malignant neoplasms (n=4); (c) thickness of primary gastric lesions <10mm on CT (n=1); (d) CT images with obvious artifacts (n=1). Finally, 42 patients were included in our study. The following clinicopathological data were retrospectively collected from patients’ medical records: age, gender, Eastern Cooperative Oncology Group performance status score (ECOG PS), treatment regimen, treatment cycles of ICI before first follow-up, Lauren subtype, degree of differentiation, PD-L1 status, MMR status, EBV status, peritoneal metastasis, hepatic metastasis, the number of metastatic sites. We registered patients with PD-L1 CPS ≥ 1 as PD-L1 positive cases (4). The flowchart is shown in Figure 1.
Figure 1 Study design for the evaluation of response prediction to immune checkpoint inhibitors in stage IV gastric cancer patients based on radiomics features. preAintra: intratumoral regions on baseline atrial phase; preVintra: intratumoral regions on baseline venous phase, preAperi: peritumoral regions on baseline atrial phase; preVperi: peritumoral regions on baseline venous phase; postAintra: intratumoral regions on follow-up atrial phase; postVintra: intratumoral regions on follow-up venous phase; postAperi: peritumoral regions on follow-up atrial phase; postVperi: peritumoral regions on follow-up venous phase.
Treatment regimens and follow-up protocol
There were 29 patients received anti-PD-1/PD-L1 treatment alone, including seven patients received GLS-010 (zimberelimab) (240mg d1 d15 Q28d), six patients received CS1003 (nofazinlimab) (200mg d1 Q21d), four patients received toripalimab (3mg/kg d1 d15 Q28d), three patients received BGB-A317 (tislelizumab) (200mg d1 Q21d), three patients received pembrolizumab (200mg d1 Q14d), two patients received atezolizumab (1200mg d1 Q21d), two patients received MSB2311 (20mg/kg d1 Q21d), one patients received LZM009 (432mg d1 Q28d) and one patients received Sintilimab (200mg d1 Q21d). There were 13 patients received anti-PD-1 in combination with anti-CTLA-4 treatment, including seven patients received Sintilimab + IBI310 (Sintilimab 200mg d1 Q21d, IBI310 68mg d1 Q42d) and six patients received Nivolumab + Ipilimumab (Nivolumab 1mg/kg d1, d22 Q42d, Ipilimumab 3mg/kg d1 d22 Q42d).
All patients conducted follow-up every two to three cycles of ICI treatment, including enhanced abdominal/pelvic CT scans until the resistance to ICI therapy. PFS was defined as the time from the start of ICI treatment to disease progression, death from any cause, or the cutoff date of November 12, 2021. Patients without any progression or death at the end of the follow-up period were censored.
CT examination
All patients underwent abdominal/pelvic contrast-enhanced CT examinations after fasting for more than eight hours. 10 mg anisodamine (654-2, Hangzhou Minsheng Pharma) was administered intramuscularly to reduce gastrointestinal motility before CT examination. Next, 6g gas-producing crystals with 10ml warm water were given orally shortly before the examination. All patients underwent a quick respiratory training session to prevent potential respiratory artifacts. The CT scanner was either the LightSpeed 64 VCT or the Discovery CT750 HD, with a peak tube voltage of 120 kVp, an automatic tube current-time product, a collimation thickness of 64 x 0.625 mm, a helical pitch of 0.984:1, 5-mm scanning thickness, and 0.625-mm reconstructed thickness. Patients were scanned in the supine position, and scan coverage started from the diaphragmatic dome until 2cm below the lower margin of symphysis ossium pubis. All patients were injected with nonionic contrast material through the antecubital vein at a rate of 3.5ml/s (1.5ml/kg of body weight, iohexol 300mg I/ml, Omnipaque, GE Healthcare). Arterial and venous phase scanning were performed at 40s and 70s, respectively, following contrast media injection.
Image analysis and segmentation
Baseline and first follow-up CT scans in arterial and venous phases were analyzed by two radiologists with 20 and 3 years of experience in gastrointestinal CT interpretation, respectively (TL and LJZ). Both radiologists were blinded to the clinical and histopathological information. However, they did know the anatomical location of gastric cancer. Two intratumoral regions of interest (ROI) were manually contoured–one ROI for the arterial phase and another ROI for the venous phase–on the largest area of the gastric lesions (axial plane) using the ITK-SNAP software (v.3.6.0, http://www.itksnap.org). To capture peritumoral information, the slice image was uniformly interpolated to 0.6 mm per pixel, and a peripheral ring was then created automatically by dilating the tumor boundaries by 7 pixels (4.2mm) on the outside and shrinking by 7 pixels (4.2mm) on the inside (18). Secondly, the modification was conducted manually on the pre-modified peripheral ring to exclude the gastric cavity and the area covering the surrounding organs and large vessels (Figure 2).
Figure 2 Schematic illustration of the steps followed to draw intratumoral and peritumoral ROIs. (A, D): the largest axial sections of gastric lesions were manually contoured (intratumoral region of ROI). (B, E): the pre-modified peripheral ring was automatically generated. (C, F): modifying was conducted manually on the pre-modified peripheral ring to exclude the gastric cavity and the part covering the surrounding organs (pancreas) and large vessels (in C: solid line: peritumoral region of ROI; dotted line: the part of ROI been deleted).
Feature extraction and selection
First, we uniformly resampled the CT image and its ROI annotation so that the spacing parameters in the x-, y-, and z-axis were 0.6, 0.6, and 5.0, respectively. The texture features were extracted from each ROI of each patient’s CT image using the open-source python platform Pyradiomics (version 2.1.2, https://pyradiomics.readthedocs.io/en/latest/#). We extracted a total of 192 features for each ROI, including eight shape features, 36 first-order statistics, 46 gray level co-occurrence matrices, 32 gray level run length matrices, 32 gray level size zone matrices, 28 gray level dependence matrices, and ten neighboring gray tone difference matrices. Eight sets of radiomics features were derived from intratumoral and peritumoral regions at baseline arterial and venous phases, and follow-up CT scans, including features from intratumoral regions at baseline atrial phase (preAintra), intratumoral regions at baseline venous phase (preVintra), peritumoral regions at baseline atrial phase (preAperi), peritumoral regions at baseline venous phase (preVperi), intratumoral regions at follow-up atrial phase (postAintra), intratumoral regions at follow-up venous phase (postVintra), peritumoral regions at follow-up atrial phase (postAperi), and peritumoral regions at follow-up venous phase (postVperi). Calculate the changes between baseline and follow-up features by subtracting the values of CT features of follow-up and baseline, which provided four corresponding sets of delta features (ΔAintra, ΔVintra, ΔAperi, and ΔVperi).
All radiomics features were standardized by subtracting the mean value and dividing by the standard deviation. Intraclass correlation coefficients (ICCs) based on a multiple-raters, two-way random-effects model were calculated to assess the stability and reproducibility of radiomic features within groups. To ensure reliability for all twelve sets of radiomics features, we only reserved radiomic features with ICC estimates > 0.80, and further selection was then conducted in the data obtained by TL. Furthermore, we used the Cox proportional hazards regression method with the least absolute shrinkage and the selection operator (LASSO) penalty with four-fold cross-validation to select the most useful predictive features from intratumoral and peritumoral regions, respectively (19). Since the total patient number was limited, the most significant nonzero feature in intratumoral and peritumoral regions was selected to avoid overfitting.
Statistical analysis
Continuous variables were presented as the mean with standard deviation (SD) or median with interquartile ranges (IQR) based on normal distribution or not. Categorical variables were shown as numbers with percentages. PFS was estimated using the Kaplan–Meier method, and the log-rank test was employed to compare differences in survival probability. The Cox proportional hazards model was used for univariate and multivariate analyses. P values less than 0.10 in univariate analysis were subsequently included in the multivariate analyses where enter feature selection was used. Harrell’s concordance index (C-index) was calculated to evaluate prognostic ability. Statistical analysis was conducted using R software (R 4.0.4, The R Foundation for Statistical Computing, Vienna, Austria). All statistical tests were two-sided, and a value of P<0.05 was considered significant.
Results
Patient characteristics
A total of 42 patients were included in this study. The median follow-up time and the median time for PFS were 736 (IQR: 656, 1266) and 133 (IQR: 61, 483) days, respectively. The patients’ clinicopathological data are summarized in Table 1. Univariate analysis revealed that age and Lauren type were associated with PFS. In contrast, other clinicopathological characteristics were not found to have a prognostic impact. The K–M analysis showed that older patients (>62 years, median value) had more prolonged PFS compared to younger patients [median PFS time: younger patients, 92 (IQR: 45, 165) days; older patients, 483 (IQR: 73, not reached) days; P=0.001]. Patients with intestinal-type gastric cancer showed more prolonged PFS than patients with a different Lauren type [median PFS time: intestinal type, 195 (IQR: 100, 649) days, reference; diffuse type, 63 (IQR: 45, 92) days, P =0.003; mixed type, 127 (IQR: 54, 134) days, P=0.087].
Radiomics feature selection
A three-step radiomics feature selection procedure was applied. In the first step, 2304 radiomics features were extracted from twelve sets of features. Consequently, 99 features were further enrolled with ICC>0.80 as a reliability standard, including 70 intratumoral features (preAintra: 14; preVintra: 25; postAintra: 8; postVintra:13; ΔAintra: 2, ΔVintra: 8) and 29 peritumoral features (preAperi: 5; preVperi: 8; postAperi: 7; postVperi: 8; ΔAperi: 0, ΔVperi: 1). The third step involved the selection of features with the highest coefficient based on the Lasso COX method, which included ΔVintra_original_glszm_Zone Variance (ΔVintra_ZV) from the intratumoral regions and postVperi_original_shape_Sphericity (postVperi_Sphericity) from the peritumoral regions.
Radiomics feature analysis
The optimal cut-off values were -0.09 and 0.88 for ΔVintra_ZV and postVperi_Sphericity determined by X-tile (version 3.6.1), respectively. The K–M analysis suggested that the PFS of stage IV gastric cancer patients with a high ΔVintra_ZV value (> cutoff value) was worse than that of patients with a low value (≤ cutoff value), with a median PFS of 402 vs. 64 days (P=0.000, log-rank test). The PFS of stage IV gastric cancer patients with a low postVperi_Sphericity value was worse than that of patients with a high value, with a median PFS of 100 vs. 589 days (P=0.012, log-rank test) (Figure 3). We performed additional analyses within subgroups of gastric cancer patients who had either not been tested or had already tested negative for biomarkers, including PD-L1, MMR, and EBV. Our findings revealed that ΔVintra_ ZV and postVperi Sphericity could stratify patients in all three subgroups according to their PFS (Figures 4A–F).
Figure 3 Kaplan–Meier estimates of progression-free survival in all patients according to (A) ΔVintra_ZV (ΔVintra_original_glszm_Zone Variance); (B) postVperi _Sphericity (postVperi_original_shape_Sphericity).
Figure 4 Kaplan–Meier estimates of progression-free survival in patients who had either not been tested or had already tested negative for biomarkers. (A, D) PD-L1, (B, E) MMR, and (C, F) EBV according to radiomics features. ΔVintra_ZV, ΔVintra_original_glszm_Zone Variance; postVperi_Sphericity, postVperi_original_shape_Sphericity.
ΔVintra_ZV and postVperi_Sphericity were both significant in univariate analysis (hazard ratio [HR], 2.320; 95% confidence interval [CI]: 1.478–3.641, P=0.000; HR, 0.601; 95% CI: 0.410–0.882, P=0.009). After controlling for age, Lauren type, peritoneal metastasis, and number of metastatic sites, ΔVintra_ZV were still independent predictor of survival (HR, 1.911; 95% CI: 1.163–3.142; P=0.011). However, postVperi_Sphericity had no association with PFS (HR, 0.690; 95% CI: 0.421–1.132; P=0.142). ΔVintra_ ZV and postVperi_Sphericity yielded a C-index of 0.705 (95% CI: 0.625–0.785) and 0.632 (95% CI: 0.528–0.736), respectively.
Discussion
This study initially explored the relationship between delta radiomics with the prognosis of patients with stage IV gastric cancer receiving ICI. Our findings revealed that ΔVintra_ZV and postVperi_Sphericity from the intratumoral and peritumoral regions, respectively, could classified patients with survival outcomes and ΔVintra_ZV was the independent predictor for PFS.
Previous studies have reported that on-treated tumor samples of patients with effective ICI response showed increased immune cell abundance and a low percentage of tumor cells (20, 21). A previous study explored the association between radiomics features in pan-cancer and CD8 cell abundance within the tumor. The relatively homogeneous tumors were associated with increased pre-existing CD8+ cell infiltration and better prognosis; in contrast, tumors composed of highly proliferating tumor cells exhibited a more heterogeneous radiomics texture (22). In our study, the low ΔVintra_ZV score indicated that the texture of gastric lesions changed from heterogeneous to homogenous and thus were more likely to be observed in patients with prolonged survival after ICI treatment. We hypothesized that this change may indicate immune cell infiltration and good tumor response to ICI therapy, whereas texture changes towards non-uniformity may indicate a high extent of tumor cell proliferation and resistance to ICI. Similar to our results, Basler et al. also suggested that the changes in the CT texture of metastatic melanoma from heterogeneity to homogeneity during ICI treatment are more likely to represent pseudoprogression, whereas changes from homogeneity to heterogeneity may indicate true tumor progression (23). Accordingly, patients with pseudoprogression showed longer survival compared with patients with true progression (23).
Given the dynamic change of tumor-immune interactions, biomarkers capable of tracking tumor evolution during the treatment course may provide more information on patients’ prognoses. A previous histological study showed that early on-treatment samples were more predictive of the response to ICI compared to the mere assessment of baseline samples (20). Although biopsies provide a method to capture the dynamic change of tumors, invasive re-biopsy may not be frequently conducted in real-world clinical practice. In our study, both predictors, ΔVintra_ ZV and postVperi_ Sphericity, incorporated post-treatment CT texture features and could predict the response of patients receiving ICI. Consistently, Khorrami et al. developed radiomics models to predict the ICI response and OS of patients with non–small cell lung cancer (NSCLC). The results showed that the performance of models combining baseline and follow-up features was better than the baseline radiomics model alone (21). Similar results also have been reported in patients treated with radiation therapy and chemotherapy (24, 25).
Khorrami et al. have shown that the ICI response prediction performance of combined radiomics from intra- and peritumoral regions in NSCLC was superior to radiomics from the intratumoral region alone (21). The authors also found that the density of immune infiltration in surgical specimens after ICI was correlated with peritumoral delta radiomics (21). In previous articles, the association between peritumoral radiomics and pathological characteristics of gastric cancer were also studied, but the prediction value of peritumoral radiomics were different (18, 26, 27). Some large-scale studies showed peritumoral radiomics features were one of the important factors to determine the tumor immune microenvironment of gastric cancer and had the prognosis predicting value, while another large-scale study showed peritumoral features may be inapplicable for predicting the Lauren classification of gastric cancer (18, 26, 27). We noticed that the peritumoral ROI in their studies were all a peripheral ring, the same with the peritumoral ROI used in lung cancer (18, 21, 26, 27). However, unlike lung cancer which is surrounded by consistent pulmonary tissue, gastric cancer is usually surrounded by air in the stomach cavity, fat tissue of peritoneum and adjacent organs. We suppose the peritumoral ROI covering air, fat, gastric cancer, and even other organs may influence the precision of information from radiomics features of the peritumoral region, although thickness of ROI around the tumor used in previous studies were smaller than ours. Therefore, in our research, we put effort into modifying the automatically generated peripheral ring, especially deleting the adjacent organs and air covered by the automatically generated ROI. However, unfortunately, such procedure increased the interobserver variability, and only 29 peritumoral features had ICC > 0.80 (70 intratumoral features had ICC > 0.80). Moreover, it was labor-intense to modify the peritumoral ROI of all patients. Taken together, we believe that further research is needed to explore the appropriate method for extracting information from the peritumoral region of gastric cancer.
In our study, patients with a high postVperi_Sphericity score demonstrated a trend towards a more promising survival outcome compared to patients with a low score. Sphericity measured the roundness of the shape and a larger value meant that the shape of ROI resembled a circle (28). Given that all ROIs of the peritumoral area appeared long and narrow, high score of postVperi_ Sphericity could be considered in two aspects, the larger the width and the shorter the length of the ROI. The width of the pre-modified peripheral ring was consistent among different patients (9.4mm in total) when first generated automatically. In patients with low visceral adipose tissue, the pre-modified peripheral ring may cover adjacent organs and thus should be manually modified, contributing to a smaller width. Poor nutritional status, including low visceral fat, has been associated with worse survival outcomes in patients treated with ICI therapy (29, 30). The length of ROI could be regarded as the maximum tumor extension on stomach. Maximum tumor diameter has been proved to be a negative factor for prognosis of patients with gastric cancer (31). Therefore, we considered that a low sphericity score may reflect poor nutritional status and high tumor burden and indicate worse survival after ICI treatment.
Our study has some limitations. First, the sample size of this retrospective study was relatively small. However, the data obtained from patients treated with ICI monotherapy were informative and of great value for assessing response after ICI treatment. In contrast, a combination regimen, such as ICI and chemotherapy, may cause confounding factors. Our study should be considered exploratory. Second, histological biomarker data were unavailable from all patients in this study. Since not all hospitals have accredited laboratories to carry out complex immunohistochemistry protocols, it is worthwhile to investigate the predictive value of radiomics features in patients who have not been tested or have already tested negative for biomarkers to provide a method of selecting appropriate treatment. Third, pathology confirmation of immune cell infiltration from post-treatment samples was absent. Future studies should aim to evaluate the relationship between radiomics features and immune cell infiltration in post-treatment gastric cancer samples.
Conclusions
Given the complexity of the intrinsic biological pathway of the tumor microenvironment, current biomarkers alone, including PD-L1, dMMR, and EBV status, cannot predict patient prognosis completely. Radiomics features complement these widely accepted histological biomarkers and can be considered candidate biomarkers that can reflect tumor phenotype and provide longitudinal surveillance. Radiomics features have the potential to be used as cost-effective screening tools that can be applied in clinical practice when administering ICI treatment to patients with gastric cancer.
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 Peking Cancer Hospital. The ethics committee waived the requirement of written informed consent for participation.
Author contributions
JL and ZC contributed equally to this work and share first authorship; BD, XZ, LT, and LS contributed equally to this work and share corresponding authorship. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by Beijing Natural Science Foundation (No. Z180001; Z200015); National Natural Science Foundation of China (No. 81801778, 12090022, 11831002, 91959205); Science Foundation of Peking University Cancer Hospital (JC202301); PKU-Baidu Foundation (No. 2020BD027); Pilot Project of Public Welfare Development Reform of Beijing-based Medical Research Institutes under Grant 2019-1.
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.
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:209–49. doi: 10.3322/caac.21660
2. Shitara K, Özgüroğlu M, Bang YJ, Di Bartolomeo M, Mandalà M, Ryu MH, et al. Pembrolizumab versus paclitaxel for previously treated, advanced gastric or gastro-oesophageal junction cancer (KEYNOTE-061): a randomised, open-label, controlled, phase 3 trial. Lancet (2018) 392:123–33. doi: 10.1016/S0140-6736(18)31257-1
3. Kang YK, Boku N, Satoh T, Ryu MH, Chao Y, Kato K, et al. Nivolumab in patients with advanced gastric or gastro-oesophageal junction cancer refractory to, or intolerant of, at least two previous chemotherapy regimens (ONO-4538-12, ATTRACTION-2): A randomised, double-blind, placebo-controlled, phase 3 trial. Lancet (2017) 390:2461–71. doi: 10.1016/S0140-6736(17)31827-5
4. Fuchs CS, Doi T, Jang RW, Muro K, Satoh T, Machado M, et al. Safety and efficacy of pembrolizumab monotherapy in patients with previously treated advanced gastric and gastroesophageal junction cancer: Phase 2 clinical KEYNOTE-059 trial. JAMA Oncol (2018) 4:e180013. doi: 10.1001/jamaoncol.2018.0013
5. Sasaki A, Nakamura Y, Mishima S, Kawazoe A, Kuboki Y, Bando H, et al. Predictive factors for hyperprogressive disease during nivolumab as anti-PD1 treatment in patients with advanced gastric cancer. Gastric Cancer (2019) 22:793–802. doi: 10.1007/s10120-018-00922-8
6. Hagi T, Kurokawa Y, Kawabata R, Omori T, Matsuyama J, Fujitani K, et al. Multicentre biomarker cohort study on the efficacy of nivolumab treatment for gastric cancer. Br J Cancer (2020) 123:965–72. doi: 10.1038/s41416-020-0975-7
7. Kim ST, Cristescu R, Bass AJ, Kim KM, Odegaard JI, Kim K, et al. Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer. Nat Med (2018) 24:1449–58. doi: 10.1038/s41591-018-0101-z
8. Eso Y, Shimizu T, Takeda H, Takai A, Marusawa H. Microsatellite instability and immune checkpoint inhibitors: Toward precision medicine against gastrointestinal and hepatobiliary cancers. J Gastroenterol (2020) 55:15–26. doi: 10.1007/s00535-019-01620-7
9. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of gastric adenocarcinoma. Nature (2014) 513:202–9. doi: 10.1038/nature13480
10. Munari E, Rossi G, Zamboni G, Lunardi G, Marconi M, Sommaggio M, et al. PD-L1 assays 22C3 and SP263 are not interchangeable in non-small cell lung cancer when considering clinically relevant cutoffs: An interclone evaluation by differently trained pathologists. Am J Surg Pathol (2018) 42:1384–9. doi: 10.1097/PAS.0000000000001105
11. Park Y, Koh J, Na HY, Kwak Y, Lee KW, Ahn SH, et al. PD-L1 testing in gastric cancer by the combined positive score of the 22C3 PharmDx and SP263 assay with clinically relevant cut-offs. Cancer Res Treat (2020) 52:661–70. doi: 10.4143/crt.2019.718
12. Seymour L, Bogaerts J, Perrone A, Ford R, Schwartz LH, Mandrekar S, et al. iRECIST: guidelines for response criteria for use in trials testing immunotherapeutics. Lancet Oncol (2017) 18:e143–52. doi: 10.1016/S1470-2045(17)30074-8
13. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images are more than pictures, they are data. Radiology (2016) 278:563–77. doi: 10.1148/radiol.2015151169
14. Xu Q, Sun Z, Li X, Ye C, Zhou C, Zhang L, et al. Advanced gastric cancer: CT radiomics prediction and early detection of downstaging with neoadjuvant chemotherapy. Eur Radiol (2021) 31:8765–74. doi: 10.1007/s00330-021-07962-2
15. Park KJ, Lee JL, Yoon SK, Heo C, Park BW, Kim JK. Radiomics-based prediction model for outcomes of PD-1/PD-L1 immunotherapy in metastatic urothelial carcinoma. Eur Radiol (2020) 30:5392–403. doi: 10.1007/s00330-020-06847-0
16. Yoon SH, Kim YH, Lee YJ, Park J, Kim JW, Lee HS, et al. Tumor heterogeneity in human epidermal growth factor receptor 2 (HER2)-positive advanced gastric cancer assessed by CT texture analysis: Association with survival after trastuzumab treatment. PloS One (2016) 11:e0161278. doi: 10.1371/journal.pone.0161278
17. Liang Z, Huang A, Wang L, Bi J, Kuang B, Xiao Y, et 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
18. Jiang Y, Wang H, Wu J, Chen C, Yuan Q, Huang W, et al. Noninvasive imaging evaluation of tumor immune microenvironment to predict outcomes in gastric cancer. Ann Oncol (2020) 31:760–8. doi: 10.1016/j.annonc.2020.03.295
19. Wada T, Yokota H, Horikoshi T, Starkey J, Hattori S, Hashiba J, et al. Diagnostic performance and inter-operator variability of apparent diffusion coefficient analysis for differentiating pleomorphic adenoma and carcinoma ex pleomorphic adenoma: Comparing one-point measurement and whole-tumor measurement including radiomics approach. Jpn J Radiol (2020) 38:207–14. doi: 10.1007/s11604-019-00908-1
20. Helmink BA, Reddy SM, Gao J, Zhang S, Basar R, Thakur R, et al. B cells and tertiary lymphoid structures promote immunotherapy response. Nature (2020) 577:549–55. doi: 10.1038/s41586-019-1922-8
21. Khorrami M, Prasanna P, Gupta A, Patil P, Velu PD, Thawani R, et al. Changes in CT radiomic features associated with lymphocyte distribution predict overall survival and response to immunotherapy in non-small cell lung cancer. Cancer Immunol Res (2020) 8:108–19. doi: 10.1158/2326-6066.CIR-19-0476
22. Sun R, Limkin EJ, Vakalopoulou M, Dercle L, Champiat S, Han SR, et al. A radiomics approach to assess tumour-infiltrating CD8 cells and response to anti-PD-1 or anti-PD-L1 immunotherapy: An imaging biomarker, retrospective multicohort study. Lancet Oncol (2018) 19:1180–91. doi: 10.1016/S1470-2045(18)30413-3
23. Basler L, Gabryś HS, Hogan SA, Pavic M, Bogowicz M, Vuong D, et al. Radiomics, tumor volume, and blood biomarkers for early prediction of pseudoprogression in patients with metastatic melanoma treated with immune checkpoint inhibition. Clin Cancer Res (2020) 26:4414–25. doi: 10.1158/1078-0432.CCR-20-0020
24. Lu L, Dercle L, Zhao B, Schwartz LH. Deep learning for the prediction of early on-treatment response in metastatic colorectal cancer from serial medical imaging. Nat Commun (2021) 12:6654. doi: 10.1038/s41467-021-26990-6
25. Xu Y, Hosny A, Zeleznik R, Parmar C, Coroller T, Franco I, et al. Deep learning predicts lung cancer treatment response from serial medical imaging. Clin Cancer Res (2019) 25:3266–75. doi: 10.1158/1078-0432.CCR-18-2495
26. Huang W, Jiang Y, Xiong W, Sun Z, Chen C, Yuan Q, et al. Noninvasive imaging of the tumor immune microenvironment correlates with response to immunotherapy in gastric cancer. Nat Commun (2022) 13:5095. doi: 10.1038/s41467-022-32816-w
27. Wang XX, Ding Y, Wang SW, Dong D, Li HL, Chen J, et al. Intratumoral and peritumoral radiomics analysis for preoperative Lauren classification in gastric cancer. Cancer Imaging (2020) 20:83. doi: 10.1186/s40644-020-00358-3
28. Lorensen WE, Cline HE. Marching cubes: A high resolution 3D surface construction algorithm. ACM SIGGRAPH Comput Graph (1987) 21:163–9. doi: 10.1145/37402.37422
29. Kim YY, Lee J, Jeong WK, Kim ST, Kim JH, Hong JY, et al. Prognostic significance of sarcopenia in microsatellite-stable gastric cancer patients treated with programmed death-1 inhibitors. Gastric Cancer (2021) 24:457–66. doi: 10.1007/s10120-020-01124-x
30. Martini DJ, Shabto JM, Goyal S, Liu Y, Olsen TA, Evans ST, et al. Body composition as an independent predictive and prognostic biomarker in advanced urothelial carcinoma patients treated with immune checkpoint inhibitors. Oncologist (2021) 26:1017–25. doi: 10.1002/onco.13922
31. Kunisaki C, Makino H, Takagawa R, Oshima T, Nagano Y, Kosaka T, et al. Tumor diameter as a prognostic factor in patients with gastric cancer. Ann Surg Oncol (2008) 15(7):1959–67. doi: 10.1245/s10434-008-9884-3
Glossary
Keywords: immunotherapy, gastric cancer, prognosis, radiomics, computed tomography
Citation: Li J, Chen Z, Chen Y, Zhao J, He M, Li X, Zhang L, Dong B, Zhang X, Tang L and Shen L (2023) CT-based delta radiomics in predicting the prognosis of stage IV gastric cancer to immune checkpoint inhibitors. Front. Oncol. 12:1059874. doi: 10.3389/fonc.2022.1059874
Received: 02 October 2022; Accepted: 29 November 2022;
Published: 04 January 2023.
Edited by:
Damiano Caruso, Sapienza University of Rome, ItalyReviewed by:
Paolo Tini, Siena University Hospital, ItalyHuan Zhang, Shanghai Jiao Tong University, China
Copyright © 2023 Li, Chen, Chen, Zhao, He, Li, Zhang, Dong, Zhang, Tang and Shen. 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: Bin Dong, dongbin@math.pku.edu.cn; Xiaotian Zhang, zhangxiaotianmed@163.com; Lei Tang, tangl@bjcancer.org; Lin Shen, shenlin@bjmu.edu.cn
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share corresponding authorship