An artificial intelligence method to assess the tumor microenvironment with treatment outcomes for gastric cancer patients after gastrectomy
Journal of Translational Medicine volume 20, Article number: 100 (2022)
The tumor microenvironment (TME) plays an important role in the occurrence and development of gastric cancer (GC) and is widely used to assess the treatment outcomes of GC patients. Immunohistochemistry (IHC) and gene sequencing are the main analysis methods for the TME but are limited due to the subjectivity of observers, the high cost of equipment and the need for professional analysts.
The ImmunoScore (IS) was developed in the TCGA cohort and validated in GEO cohorts. The Radiomic ImmunoScore (RIS) was developed in the TCGA cohort and validated in the Nanfang cohort. A nomogram was developed and validated in the Nanfang cohort based on RIS and clinical features.
For IS, the area under the curves (AUCs) were 0.798 for 2-year overall survival (OS) and 0.873 for 4-year overall survival. For RIS, in the TCGA cohort, the AUCs distinguishing High-IS or Low-IS and predicting prognosis were 0.85 and 0.81, respectively; in the Nanfang cohort, the AUC predicting prognosis was 0.72. The nomogram performed better than the TNM staging system according to the ROC curve (all P < 0.01). Patients with TNM stage II and III in the High-nomogram group were more likely to benefit from adjuvant chemotherapy than Low-nomogram group patients.
The RIS and the nomogram can be used to assess the TME, prognosis and adjuvant chemotherapy benefit of GC patients after radical gastrectomy and are valuable additions to the current TNM staging system. High-nomogram GC patients may benefit more from adjuvant chemotherapy than Low-nomogram GC patients.
The tumor-node-metastasis (TNM) staging system is widely used in the classification of cancer; is based on the primary tumor (T), regional nodes (N) and metastasis (M); and helps in choosing surgery, chemotherapy, etc. However, recent studies have shown that the TNM staging system does not perfectly predict the prognosis of gastric cancer (GC) patients . Those who had the same TNM stage and received similar treatment had widely varied clinical outcomes [2, 3]. As increasing evidence has indicated the critical forecasting capability of immune infiltration in GC patients , having immune cell markers along with TNM staging may help.
There is reportedly an exact correlation between tumor and immunocyte infiltration of the tumor microenvironment (TME) , which might lead to different prognoses of GC patients. In addition, the evaluation of the immune microenvironment seems to be useful in predicting survival and chemotherapy response for patients . Therefore, specific analysis of an individual's tumor microenvironment and targeted selection of therapy regimens based on the assessment will be considerably helpful for patients.
Because the importance of immunocyte infiltration of the TME has been recognized, ways in which to analyze the immune microenvironment have been developed. To analyze the types and composition of immunocytes in the TME, biopsies were initially analyzed using immunohistochemistry (IHC) , which is quite accurate and realistic . However, to a certain extent, the interpretation of IHC results depends on the subjectivity of the observer. Another method of determining immunocyte infiltration of the TME uses indirect information about the composition of immunocytes that can be calculated by gene expression profiles , which was proven to be effective. Because the calculation process of the results is traceable and easy to repeat, the results might be more objective.
Although calculation by gene sequencing is already objective, it still has some limitations, including a high cost for patients and specialized gene sequencing equipment for hospitals. Medical imaging may be a viable solution to this problem. Medical images are more than pictures, but data , such as radiomics, have been developed and successfully applied in many fields. Furthermore, many scholars have found that some tumor imaging features are closely related to tumor gene expression, providing a new idea for survival prediction and clinical strategy adjustment by combining their analyses [9, 10]. For evaluating the tumor microenvironment of GC patients, radiomics enables the extraction of innumerable quantitative features from medical images with high-throughput computing  and the establishment of a model. Several studies have explored the association between imaging features and tumor-infiltrating lymphocytes, such as CD8 + cells .
Supervised machine learning (ML) is a type of artificial intelligence that is widely used and shows great potential in precision oncology [13, 14]. As a mature machine learning method, support vector machines (SVMs) have great potential in solving cancer prediction problems based on diverse and complex clinical data [1, 15].
In this study, we analyzed the relationship between the gene expression the and CT images of GC patients to assess immunocyte infiltration of the TME based on the SVM model and developed an excellent nomogram based on clinicopathological factors and the SVM model to predict survival and chemotherapy benefit in GC patients.
Methods and materials
Study design and patients
The design of the entire study is shown in Fig. 1. A total of 4 independent cohorts, which included 936 GC samples, were enrolled in this study. One cohort (n 1 = 44) was from The Cancer Genome Atlas (TCGA) datasets (https://portal.gdc.cancer.gov/) and 2 cohorts (GSE62254, n2 = 300, and GSE15459, n3 = 192) were from Gene Expression Omnibus (GEO) datasets (https://www.ncbi.nlm.nih.gov/geo/), and a cohort (n4 = 400) from Nanfang Hospital (Guangzhou, China), which were used to analyze the relationship between gene expression and immune microenvironment of tumors. We developed an SVM model to analyze the relationship between gene expression and CT images based on a TCGA cohort, aiming to learn the immune microenvironment from CT images. CT images were obtained from The Cancer Imaging Archive (TCIA, https://www.cancerimagingarchive.net/). One cohort that comprised 400 consecutive patients from Nanfang Hospital (Guangzhou, China) from October 2004 to September 2011 was used to develop the nomogram to predict survival and chemotherapy benefit for GC patients. All of these patients underwent partial or total radical gastrectomy. The inclusion criteria were histologically confirmed GC and patients who underwent standard unenhanced and contrast-enhanced abdominal CT < 30 days before surgical resection. Those who received preoperative chemotherapy, had severe surgical complications, did not have complete follow-up data over 3 years, or had other coexisting cancers were excluded. The entire study design is shown in Fig. 1.
The TNM staging was reclassified, and the postoperative chemotherapy regimen was taken based on the guidelines of the National Comprehensive Cancer Network (NCCN) . In the Nanfang Hospital cohort, up to 223 patients were treated with adjuvant chemotherapy, including 114 (51.1%) patients treated with the XELOX (capecitabine-oxaliplatin, capecitabine 1000 mg/m2/days day1–day14, oxaliplatin 130 mg/m2 once) regimen, 90 (40.4%) patients treated with the FOLFOX (fluorouracil-folinic acid-oxaliplatin, fluorouracil 400 mg/m2 at day1 and 2400–3600 mg/m2 for 46 h, folinic acid 400 mg/m2 once, oxaliplatin 85 mg/m2 once) regimen and 19 (8.5%) patients treated with 5-FU treatment alone or other regimens.
The follow-up duration was calculated from the day of surgery to the last follow-up. We defined the time to all-cause death as overall survival (OS) and defined the time to tumor recurrence at any site or death as disease-free survival (DFS). Ethical approval was obtained for this retrospective study, and informed consent was obtained from the patients.
Assessment of immunocyte infiltration of the TME and calculation of the ImmunoScore (IS)
To predict the prognosis of patients using gene information to assess immune infiltration, we included 44 cases from TCGA as the training cohort, 300 cases from GSE62254 and 192 cases from GSE15459 as 2 validation groups to establish and evaluate the role of the ImmunoScore (IS) in clinical prognosis.
To quantify the proportions of 22 types of infiltrating immunocytes of the TME, normalized RNA sequencing data were analyzed using the CIBERSORT algorithm . The immunocyte type fractions of patients in the training cohort using the CIBERSORT algorithm could be downloaded directly from a comprehensive resource  (http://timer.cistrome.org/). For each sample, the sum of all estimates of the immunocyte type component equaled 100%. The least absolute shrinkage and selection operator (LASSO) method was used to screen the most effective prognostic markers among 22 immune cell subsets, and the optimal values of the penalty parameter λ were determined by 20-fold cross-validations.
The ImmunoScore (IS) was calculated by a Cox proportional hazards model with the LASSO results. The output was expressed as coefficients for estimating the immunocytes of the TME. The calculation string formula was as follows:
where Xi represents the relative component of each immunocyte, βi is the regression coefficient for the predictor and α is often a constant. We divided GC patients into 2 groups (High-IS and Low-IS) based on the median IS in each cohort.
Exploration between ImmunoScore (IS) and immunocytes and immune checkpoint
Differential expression gene (DEG) analysis of the GSE15459 cohort was performed by applying the “limma” package at a corrected P < 0.05 and | logFC |≥ 1.5. The resulting data were used to generate volcano plots using R software (version 4.0.4). Data on the abundance of 22 types of immune infiltrates from the TIMER website were obtained. Using the “ggplot 2” package, we identified immunocytes with significant differences between the High-IS and Low-IS groups while analyzing the correlation between these cells and the IS using the “ggstatsplot” package.
CT image processing, feature extraction and development of the radiomic ImmunoScore (RIS)
To use CT information to obtain radiomics features associated with the IS of patients, we included 44 cases from TCGA as the training cohort and 400 cases from Nanfang Hospital as the validation cohort to establish a Radiomic ImmunoScore (RIS) to predict IS and estimate its connection with patient survival.
Using ITK-SNAP software (www.itksnap.org, version 3.8), two radiologists manually delineated the primary tumor on the CT images. The CT number, representing the absorptivity of tissue to X-rays, was set ranging from − 150 to − 50 Hounsfield units to exclude bone, muscle, blood vessels, and other intra-abdominal organs [19, 20]. The region of interest (ROI) was delineated by hand along the gastric wall at a distance of 1 mm from the stomach wall around the tumor lesion.
Based on extracted image features and ImmunoScore, the RIS model was developed by SVM, aiming to predict the ImmunoScore by CT images. The common image features included color, texture, shape and spatial relation. The software version we used was MATLAB 2018a, and the shape of the feature matrix extracted from a CT image was 1*9694. In the training cohort, the size of the original CT feature matrix was 44*9694, which made regression difficult. Therefore, the average weight of each feature was finally obtained after dimensional reduction using the Relief algorithm. The top 100 features were selected based on the feature weights. The AUC was regarded as the evaluation index. Six features were selected in the development of the Radiomic ImmunoScore (RIS) model when the AUC reached its peak (more details of this process are described in Additional file 1). The RIS of each cohort was also divided into High-RIS and Low-RIS groups according to the median.
Subsequently, the clinical significance of RIS was displayed intuitively by ROC, and the hazard ratios of each clinicopathological factor were analyzed between the High-RIS and Low-RIS groups.
Development and validation of the nomogram
We developed a comprehensive model, including RIS and other clinical features, and visualized it with a nomogram. To verify the accuracy and effectiveness of the nomogram, we randomly divided 400 patients from Nanfang Hospital into a training cohort and a validation cohort, with 200 cases in each group.
According to the influence degree of each factor on the outcome variable in the multifactor regression model, each value level of each variable was given a score, and the comprehensive nomogram obtained by adding was calculated by function to attain the prediction probability. Patients were then divided into High-nomogram and low-nomogram groups according to the median nomogram score.
Subsequently, the predictive ability of the nomogram for patients was demonstrated by a variety of evaluation methods, including ROC, survival curve, calibration plot decision curve analysis (DCA), net reclassification improvement (NRI) and risk factor association diagrams.
Performance of the nomogram in the identification of chemotherapy benefit
The association between the nomogram and adjuvant chemotherapy response was assessed in patients with stage II and III gastric cancer (n = 280, Nanfang cohort). All patients were divided into a chemotherapy cohort (n = 182) and a no chemotherapy cohort (n = 98). Both the total survival time and disease-free survival time were observed. Survival curves were generated with Kaplan–Meier analysis. Hazard ratios were generated with the Cox regression model.
All statistical tests were two-sided, and P < 0.050 was considered statistically significant. Statistical analysis was conducted with R software (version 4.0.4) and SPSS software (version 25.0). The packages for downloading data from TCGA and GEO are TCGAbiolinks and GEOquery. Correlations between the immune cell immune checkpoints and the Immunoscore were analyzed using Pearson’s correlation test. Survival curves were constructed by the Kaplan–Meier method and compared using the log rank test. The sensitivity and specificity of the survival prediction based on the 3 models were depicted by a time-dependent receiver operating characteristic (ROC) curve and a time-independent ROC curve, with quantification of the area under the ROC curve using the timeROC and pROC packages. The results of subgroup analysis were represented by a forest graph using the forestplot package, and the points of HR value in the graph were processed by log2 normalization. We performed Cox multivariate modeling using the rms package. After constructing the nomogram with the nomogramFormula package, we drew decision curves with the rmda package and risk factor correlation diagrams with the ggrisk package to measure the clinical value of the nomogram.
Up to 936 GC patients from 4 cohorts (TCGA, GSE62254, GSE15459 and Nanfang Hospital) were enrolled in this study. We collected and analyzed their clinical characteristics, which were summarized in Table 1: 44 patients in the TCGA cohort [38 (86.4%) were male, 6 (13.6%) were female and the mean age was 64.73 ± 9.33 years], 300 patients in the GSE62254 cohort [199 (66.3%) were male, 101 (33.7%) were female and the mean age was 61.94 ± 11.36 years], 192 patients in the GSE15459 cohort [123 (65.1%) were male, 69 (34.9%) were female and the mean age was 64.37 ± 13.24 years], 400 patients in the Nanfang Hospital cohort [276 (69.0%) were male, 124 (31.0%) were female and the mean age was 56.04 ± 10.87 years]. More information on tumors was analyzed and is listed in Table 1, including the depth of invasion, lymph node metastasis, metastasis, TNM stage, pathological type, venous invasion, lymphovascular invasion and perineural invasion.
Derivation and validation of the ImmunoScore (IS)
The immunocyte composition of the TME was calculated with RNA sequencing data from the TCGA cohort (Fig. 2A),the GSE62254 cohort (Additional file 3: Figure S1) and the GSE15459 cohort (Additional file 4: Figure S2). As a training cohort, 6 types of immunocytes (plasma cells, CD8 + T cells, memory resting CD4 T cells, Tregs, M2 macrophages, and M0 macrophages) were selected by LASSO Cox regression analysis based on the TCGA cohort (Fig. 2B, C) and calculated the immunoScore by a Cox proportional hazards model (Additional file 2). The prognostic accuracy of the IS was assessed through time-dependent ROC analysis (2-year AUC = 0.798, and 4-year AUC = 0.873, Fig. 2D). The same calculations were performed for GSE62254 and GSE15459 as the validation cohorts (the results are shown in Additional file 2, Additional file 10: Figure S3 and Additional file 11: Figure S4). The survival curves of High-IS and Low-IS groups were generated with Kaplan–Meier analysis [TCGA cohort: p < 0.001, HR = 4.55 (95% CI 1.92–10.82), GSE62254 cohort: p = 0.004, HR = 1.60 (95% CI 1.17–2.21), GSE15459 cohort: p = 0.015, HR = 1.64 (95% CI 1.09–2.47), Fig. 2E–G].
Relationship between the ImmunoScore (IS) and immune cells and immune checkpoint
Based on the gene expression data of the GSE62254 cohort, we obtained differentially expressed cells from the immune microenvironment (Fig. 2H). Notably, the infiltration of plasma cells, CD8 + T cells, activated CD4 + memory T cells, follicular helper T cells, activated NK cells and M1 macrophages significantly decreased with the higher immune risk score, while the infiltration of resting CD4 + memory T cells and M2 macrophages was significantly higher with the higher IS, which was likely related to the immunosuppressive reaction of the TME. Furthermore, we analyzed the relationship between IS and immune checkpoints, such as PD-1, PD-L1, and CTLA4 (more details are shown in Additional file 5: Table S1). In summary, IS has a strong relationship with immune cells and immune checkpoints.
Performance of the radiomic ImmunoScore (RIS) Model
As a training cohort from TCGA, radiomics signatures were compared with 44 cases based on the CT image data, which showed the best effect when 6 features were included (Fig. 3A, B). The ability of the RIS model to distinguish High-IS or Low-IS accurately and predict prognosis exhibited an AUC of 0.85 (Fig. 3D) and 0.81 (Fig. 3E), respectively. The ability of RIS to predict survival had an AUC of 0.72 (Fig. 3F) in the validation cohort from Nanfang Hospital. In the training cohort, overall survival curves were generated with Kaplan–Meier analysis, comparing High-RIS with Low-RIS [p < 0.001, HR = 4.55 (95% CI 1.91–10.82), Fig. 3G], and the result was similar in the validation cohort [p < 0.001, HR = 2.72 (95% CI 1.96–3.78), Fig. 3H]. Furthermore, in the forest plot of the Nanfang Hospital cohort, the significant Hazard Ratios were found for RIS in each subgroup, more specifically, in each age stage (< 50, 50–65, > 65), each differentiation stage (low, medium, high), and regardless of vascular invasion, nerve invasion and lymph node invasion (Fig. 3C).
Development, comparison and validation of the nomogram
A nomogram was developed by Cox regression based on the RIS and clinical features (Fig. 4A), and the Cox regression coefficient and corresponding score are shown in Additional file 6: Table S2. From the results of univariate survival analysis, RIS, tumor metastasis, lymph node positive detection rate, T stage and age were significant factors associated with OS and DFS (Additional file 7: Table S3, Additional file 8: Table S4). By calculating the total score and finding the corresponding prediction probability on the total point scale, the 1-year, 3-year and 5-year estimated survival probability could be obtained. Based on calibration plots for predicting 1-year, 3-year and 5-year survival (Fig. 4B–D), the predicted calibration curve was near the standard curve, or the 45-degree line, demonstrating that the derived nomogram had considerable prediction capability. Likewise, the decision curve analysis (DCA) depicted in Fig. 4E indicated that the nomogram had a higher net income than traditional TNM among a series of risk thresholds. The results of comparing the nomogram with the TNM staging system by Net Reclassification Improvement (NRI) are as follows: 1-year (NRI = 0.27, 95% CI = 0.07–0.44), 3-year (NRI = 0.03, 95% CI = − 0.17–0.18), and 5-year (NRI = 0.05, 95% CI = − 0.26–0.26), indicating that the nomogram had a relatively accurate prognosis. Simultaneously, the risk factor association diagram (Fig. 4F) showed the distribution difference of prognosis between High-nomogram and Low-nomogram groups distinguished by the nomogram. To predict the survival of GC patients, the performance of the nomogram was better than that of TNM staging (training cohort: nomogram AUC = 0.849 and TNM staging system AUC = 0.778, p = 0.007, Fig. 4G; validation cohort: nomogram AUC = 0.863 and TNM staging system AUC = 0.776, p = 0.001, Fig. 4J). The results of the time-dependent ROC curve analysis in the training cohort and validation cohort are shown in Fig. 4H, K. Comparing the High-nomogram and Low-nomogram groups, the survival curves of the training cohort and validation cohort are shown in Fig. 4I, L, which indicated that the prognosis of the Low-nomogram group was better than that of the High-nomogram group.
Identification of chemotherapy beneficiaries with a nomogram
In addition, we analyzed differences in chemotherapy benefit among TNM stage II and III GC patients with different groups based on the nomogram. Chemotherapy dose leads to survival benefit for GC patients (OS: p = 0.023, and DFS: p = 0.005), but the degree of benefit varied among patients in different nomogram groups (High-nomogram group and Low-nomogram group) (Fig. 5). In the High-nomogram group, GC patients could significantly benefit with chemotherapy regarding OS and DFS (OS: p = 0.004, and DFS: p < 0.001), while in the Low-nomogram group, chemotherapy did not provide a survival benefit to patients compared with no chemotherapy (OS: p = 0.21, and DFS: p = 0.47). Furthermore, we analyzed the interactions between chemotherapy and death or recurrence in the High-nomogram and Low-nomogram groups with Cox proportional hazards regression, and High-nomogram patients benefited more from chemotherapy [OS: HR 0.41 (0.26–0.64), p < 0.001; DFS: HR 0.54 (0.34–0.84), p < 0.001] than Low-nomogram patients (Additional file 9: Table S5). In summary, High-nomogram GC patients benefit more from chemotherapy than Low-nomogram GC patients.
Accurate assessment of the prognosis and beneficial choices for treating GC patients are essential for clinicians. Identifying an evidence-based and broadly utilized classification system for GC patients is challenging. Since its development, the AJCC TNM staging system has been widely used in the prognosis assessment of tumors and in guiding the selection of treatment regimens . However, as we learn more about the tumor microenvironment, we realize this staging system has inevitable limitations, including that it is difficult to explain the complex biological behavior of the tumor based on tumor invasion depth, lymph node metastasis, and distant organ metastasis only. In the same TNM stage, patients may have different outcomes, while patients in different TNM stages may have similar outcomes [22,23,24]. An increasing number of studies have found that the infiltration of immunocytes in the TME is closely related to the prognosis and treatment sensitivity of patients [25, 26]. Furthermore, immune checkpoint inhibitors have shown positive results in therapies for cancer, such as ipilimumab targeting CTLA4  and lambrolizumab targeting PD-1 .
The individualized analysis of the TME is a major advance in the history of fighting against cancer, and it has also led to the development of many methods for analyzing the immune microenvironment, such as single-cell sequencing (SNS) [29, 30] and immunohistochemistry (IHC) [1, 3]. The CIBERSORT algorithm is a method applied to calculate the proportion of immunocytes in the tumor microenvironment based on next-generation sequencing (NGS) data [17, 31]. In our study, we developed a method to calculate the immune score to assess the immunocyte infiltration of the TME based on NGS data from TCGA, and the method was validated in 2 cohorts from GEO.
We developed a further model (RIS) by SVM to assess the immune score of the TME, which established a bridge between CT images and NGS. The RIS was also validated in GC patients from the Nanfang Hospital cohort. The AUC value was significantly higher than that of the TNM staging system in assessing the prognosis of GC patients.
Using SNS, IHC and NGS may be effective in assessing immunocyte infiltration of the TME. However, these methods are costly for patients and complex for clinicians, especially in remote areas with underdeveloped medical resources, and the contradictions will become more obvious. Moreover, IHC staining of tissue slices to evaluate the immune environment might not be accurate due to the subjectivity of observers and the limitation of hysteresis.
Medical images are data , such as radiomics, which has been developed and widely used in learning more information about tumors [32, 33]. Because the characteristics of radiomics are complex and diverse, we used machine learning to identify the radiomics characteristics of patients with different immune scores rapidly and effectively. The SVM is a mature machine learning method that has demonstrated excellent performance in multifactor classification and survival prediction. Based on the RIS model, we can quickly obtain a patient's individual immune score from his/her CT images, which means that this RIS model can greatly improve the efficiency of clinicians in evaluating patients' TME and greatly reduce the high cost of sequencing for patients.
Prognosis is related not only to the TME of the tumor but also to the clinical characteristics of the patients. To improve the ability of the RIS model to evaluate the prognosis of patients and maximize the prognostic value of clinical characteristics, we combined the RIS and clinical characteristics by Cox regression and visualized it in a nomogram. Based on the nomogram, we can more easily and intuitively calculate the individual 1-year, 3-year and 5-year survival rates of patients.
Adjuvant chemotherapy is now considered a common treatment for stage II and III GC patients, improving their survival in two large-scale trials [34, 35] (the ACTS-GC trial and CLASSIC trial). However, some researchers noted that the OS benefit was moderate with adjuvant chemotherapy, meaning that not all GC patients could benefit from adjuvant chemotherapy . The identification of adjuvant chemotherapy benefits is still a challenging task for researchers . In our study, based on the nomogram, High-nomogram GC patients could benefit significantly more from adjuvant chemotherapy.
Due to the universality of CT scanning and the availability of clinical features, we can quickly and effectively obtain TME assessment, prognosis and adjuvant chemotherapy benefit of GC patients based on our RIS model and nomogram. The strategy that we have developed for evaluating the TME of GC patients is useful and easy to implement, especially for clinicians of grassroots community hospitals, in the management and treatment of GC patients.
At the same time, we have also noticed that many researchers developed novel theranostics trying to use in cancer treatment, such as nanomaterials [38, 39]. For the potential of nanomaterials to deliver, enhance, and finetune in treatment, it shows exciting effects in animal cancer models, which may be a kind of hopeful treatments in the future . However, it’s difficult to identify who would benefit in novel theranostics and easily assess the effectiveness of novel theranostics. Maybe, AI could be the solution to these bottleneck problems. For example, Tang et.al developed a AI model to analyse the distribution of intratumoral nanoparticles to optimize the cancer treatment . It's worth noting that antioxidation therapy is an important mechanism of nanomaterials in treating with cancer . It works by reducing tumor-promoting oxidative stress and remodeling the TME. Perhaps, our AI model could reflect the antitumor efficacy of nanomaterials by analyzing TME in real-time.
There were several limitations to this study. First, the number of GC patients enrolled in the TCGA cohort was small due to the limited number of cases with CT images. Second, due to the limitations of data acquisition, the data types contained in the different cohorts were not perfect, but this did not affect our RIS model and nomogram development. We hope that more data can be obtained to develop and validate our RIS model and nomogram. Third, the study was limited to collecting retrospective data, and we need further prospective research to validate our strategy in assessing TME, prognosis and adjuvant chemotherapy benefit.
In conclusion, the RIS model and the nomogram can be used to assess the TME, prognosis and adjuvant chemotherapy benefit of GC patients after radical gastrectomy and seems to be a valuable addition to the current TNM staging system. High-nomogram GC patients may benefit more from adjuvant chemotherapy than low-nomogram GC patients. By enrolling more GC patients for validation, we hope our strategy will provide more individualized assistance to clinicians and patients.
All data generated or analyzed during this study are included either in this paper or in the additional information. All readers can obtain the original code of the SVM model and the raw data of clinical trials in this study by contacting us.
Jiang Y, Xie J, Han Z, et al. Immunomarker support vector machine classifier for prediction of gastric cancer survival and adjuvant chemotherapeutic benefit. Clin Cancer Res. 2018;24(22):5574–84.
Roberto M, Botticelli A, Strigari L, et al. Prognosis of elderly gastric cancer patients after surgery: a nomogram to predict survival. Med Oncol. 2018;35(7):111.
Jiang Y, Zhang Q, Hu Y, et al. ImmunoScore signature: a prognostic and predictive tool in gastric cancer. Ann Surg. 2018;267(3):504–13.
Lazăr DC, Avram MF, Romoșan I, et al. Prognostic significance of tumor immune microenvironment and immunotherapy: novel insights and future perspectives in gastric cancer. World J Gastroenterol. 2018;24(32):3583–616.
Jiang Y, Xie J, Huang W, et al. Tumor immune microenvironment and chemosensitivity signature for predicting response to chemotherapy in gastric cancer. Cancer Immunol Res. 2019;7(12):2065–73.
Jiang Y, Wang H, Wu J, et al. Noninvasive imaging evaluation of tumor immune microenvironment to predict outcomes in gastric cancer. Ann Oncol. 2020;31:760.
Zeng D, Zhou R, Yu Y, et al. Gene expression profiles for a prognostic immunoscore in gastric cancer. Br J Surg. 2018;105(10):1338–48.
Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More than Pictures, They Are Data. Radiology. 2016;278(2):563–77.
Pope WB, et al. Non-invasive detection of 2-hydroxyglutarate and other metabolites in IDH1 mutant glioma patients using magnetic resonance spectroscopy. J Neurooncol. 2012;107(1):197–205.
Li H, et al. MR Imaging radiomics signatures for predicting the risk of breast cancer recurrence as given by research versions of MammaPrint, Oncotype DX, and PAM50 Gene Assays. Radiology. 2016;281(2):382–91.
Li Q, Ye ZX. Radiomics: the process and applications in tumor research. Zhonghua Zhong Liu Za Zhi. 2018;40(11):801–4.
Sun R, Limkin EJ, Vakalopoulou M, 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(9):1180–91.
Cuocolo R, Caruso M, Perillo T, et al. Machine Learning in oncology: a clinical appraisal. Cancer Lett. 2020;481:55–62.
Bera K, Schalper KA, Rimm DL, et al. Artificial intelligence in digital pathology - new tools for diagnosis and precision oncology. Nat Rev Clin Oncol. 2019;16(11):703–15.
Polano M, Chierici M, Dal Bo M, et al. A pan-cancer approach to predict responsiveness to immune checkpoint inhibitors by machine learning. Cancers (Basel). 2019;11:10.
National Comprehensive Cancer Network. Gastric Cancer. Version 2019. Accessed March 14, 2019. https://www.nccn.org/professionals/physician_gls/pdf/gastric.pdf
Newman AM, Liu CL, Green MR, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Sturm G, Finotello F, Petitprez F, et al. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019;35(14):i436–45.
Zanfardino M, Pane K, Mirabelli P, et al. TCGA-TCIA impact on radiogenomics cancer research: a systematic review. Int J Mol Sci. 2019;20:23.
Prior F, Smith K, Sharma A, et al. The public cancer radiology imaging collections of The Cancer Imaging Archive. Sci Data. 2017;4:170124.
Amin MB, Gress DM. AJCC Cancer Staging Manual. New York: Springer; 2017.
Amin MB, Greene FL, Edge SB, Compton CC, Gershenwald JE, Brookland RK, Meyer L, Gress DM, Byrd DR, Winchester DP. The eight edition AJCC cancer staging manual: continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin. 2017;67:93–9.
Fang C, Wang W, Deng JY, Sun Z, Seeruttun SR, Wang ZN, Xu HM, Liang H, Zhou ZW. Proposal and validation of a modified staging system to improve the prognosis predictive performance of the 8th AJCC/UICC pTNM staging system for gastric adenocarcinoma: a multicenter study with external validation. Cancer Commun. 2018;38:67.
Liu JY, Peng CW, Yang XJ, Huang CQ, Li Y. The prognosis role of AJCC/UICC 8 (th) edition staging system in gastric cancer, a retrospective analysis. Am J Transl Res. 2018;10:292–303.
Lei X, Lei Y, Li JK, et al. Immune cells within the tumor microenvironment: biological functions and roles in cancer immunotherapy. Cancer Lett. 2020;470:126–33.
Gong Z, Zhang J, Guo W. Tumor purity as a prognosis and immunotherapy relevant feature in gastric cancer. Cancer Med. 2020;9(23):9052–63.
Lynch TJ, Bondarenko I, Luft A, et al. Ipilimumab in combination with paclitaxel and carboplatin as first-line treatment in stage IIIB/IV non-small-cell lung cancer: results from a randomized, double-blind, multicenter phase II study. J Clin Oncol. 2012;30(17):2046–54.
Hamid O, Robert C, Daud A, et al. Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma. N Engl J Med. 2013;369(2):134–44.
Zhang L, Zhang Z. Recharacterizing tumor-infiltrating lymphocytes by single-cell RNA sequencing. Cancer Immunol Res. 2019;7(7):1040–6.
Yang KL, Sun Z, Bai CM, et al. Application of single-cell RNA sequencing in research on tumor immune microenvironment. Zhongguo Yi Xue Ke Xue Yuan Xue Bao. 2020;42(1):117–23.
Chen B, Khodadoust MS, Liu CL, et al. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol. 2018;1711:243–59.
Liu Z, Wang S, Dong D, et al. The applications of radiomics in precision diagnosis and treatment of oncology: opportunities and challenges. Theranostics. 2019;9(5):1303–22.
Lambin P, Leijenaar R, Deist TM, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. 2017;14(12):749–62.
Sasako M, Sakuramoto S, Katai H, Kinoshita T, Furukawa H, Yamaguchi T, Nashimoto A, Fujii M, Nakajima T, Ohashi Y. Five-year outcomes of a randomized phase III trial comparing adjuvant chemotherapy with S-1 versus surgery alone in stage II or III gastric cancer. J Clin Oncol. 2011;29:4387–93.
Bang YJ, Kim YW, Yang HK, Chung HC, Park YK, Lee KH, Lee KW, Kim YH, Noh SI, Cho JY, Mok YJ, Kim YH, Ji J, Yeh TS, Button P, Sirzén F, Noh SH. Adjuvant capecitabine and oxaliplatin for gastric cancer after D2 gastrectomy (CLASSIC): a phase 3 open-label, randomised controlled trial. Lancet. 2012;379:315–21.
Cheong JH, Yang HK, Kim H, Kim WH, Kim YW, Kook MC, Park YK, Kim HH, Lee HS, Lee KH, Gu MJ, Kim HY, Lee J, Choi SH, Hong S, Kim JW, Choi YY, Hyung WJ, Jang E, Kim H, Huh YM, Noh SH. Predictive test for chemotherapy response in resectable gastric cancer: a multi-cohort, retrospective analysis. Lancet Oncol. 2018;19:629–38.
Lee KG, Lee HJ, Oh SY, Yang JY, Ahn HS, Suh YS, Kong SH, Kim TY, Oh DY, Im SA, Lee KU, Kim WH, Bang YJ, Yang HK. Is there any role of adjuvant chemotherapy for T3N0M0 or T1N2M0 gastric cancer patients in stage II in the 7th TNM but stage I in the 6th TNM system? Ann Surg Oncol. 2016;23:1234–43.
Maleki DS, Sharifi S, Ahmadian E, Eftekhari A, Adibkia K, Lotfipour F. An update on calcium carbonate nanoparticles as cancer drug/gene delivery system. Exp Opin Drug Deliv. 2019;16:1.
Elham A, Solmaz MD, Simin S, Shahriar S, Rovshan K, Aziz E, Mohammad H. The potential of nanomaterials in theranostics of oral squamous cell carcinoma: recent progress. Trends Analyt Chem. 2019;116:23.
Bockamp E, Rosigkeit S, Siegl D, Schuppan D. Nano-enhanced cancer immunotherapy: immunology encounters nanotechnology. Cells. 2020;9:67.
Tang Y, Zhang J, He D, Miao W, Liu W, Li Y, Lu G, Wu F, Wang S. GANDA: A deep generative adversarial network conditionally generates intratumoral nanoparticles distribution pixels-to-pixels. J Control Release. 2021;336:336–43.
Zhang RX, Liu FF, Lip H, Liu J, Zhang Q, Wu XY. Pharmaceutical nanoformulation strategies to spatiotemporally manipulate oxidative stress for improving cancer therapies—exemplified by polyunsaturated fatty acids and other ROS-modulating agents. Drug Deliv Transl Res. 2022;88:67.
We are grateful to all of the patients and cooperators for participating in this retrospective study, especially for the data from TCGA and GEO.
This study was supported by Guangdong Provincial Key Laboratory of Precision Medicine for Gastrointestinal Cancer (2020B121201004), the Guangdong Natural Science Foundation General Project (2019A1515011520), Guangdong Natural Science Foundation Outstanding Youth Project (2021B1515020055) and Guangdong Provincial Major Talents Project (2019JC05Y361).
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki and comparable Chinese ethical standards, and we obtained written informed consent from all study participants.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: The process of features selecting at highest AUC.
: The immunoScore was calculated of TCGA,GSE62254 and GSE15459 cohorts.
: The immunocyte composition of the TME of GSE62254 cohort.
: The immunocyte composition of the TME of GSE15459 cohort.
: Correlation analyses between Immunoscore and gene expressions and immunocytes.
: Cox regression coefficients and nomogram score for the training cohort.
: Univariable analysis of the RIS, clinical features with overall survival in Nanfang cohort.
: Univariable analysis of the RIS, clinical features with disease-free survival in Nanfang cohort.
: Treatment Interaction with Nomogram for DFS and OS with Stage II and III Disease.
: Time-dependent ROC analysis of GSE62254 cohort.
: Time-dependent ROC analysis of GSE15459 cohort.
About this article
Cite this article
Chen, T., Li, X., Mao, Q. et al. An artificial intelligence method to assess the tumor microenvironment with treatment outcomes for gastric cancer patients after gastrectomy. J Transl Med 20, 100 (2022). https://doi.org/10.1186/s12967-022-03298-7