Development and validation of a Surveillance, Epidemiology, and End Results (SEER)-based prognostic nomogram for predicting survival in gastric cancer with multi-organ metastases.

Background: Nomogram can be used to accurately predict the prognosis of patients and guide treatment according to the individual situation of patients. This study is to investigate the independent prognostic factors for multi-organ metastases in gastric cancer (GC) patients, and construct and validate prognostic nomograms for overall survival (OS) and cancer-specific survival (CSS).
Methods: The clinical data of GC patients with multi-organ metastases from 2010 to 2018 were extracted from the Surveillance, Epidemiology, and End Results (SEER) database. The independent prognostic factors affecting the OS and CSS of the patients were screened using univariate and multivariate Cox's proportional hazards model and the Fine-Gray competing risk model. Corresponding nomogram models were constructed to predict the OS and CSS of the patients. The reliability and accuracy of the prediction model were evaluated by consistency index (C-index), area under receiver operating characteristic (ROC) curve (AUC) and calibration curve.
Results: A total of 1,386 patients were included and randomly divided into a training group (972 cases) and a validation group (414 cases) in a 7:3 ratio. Cox proportional hazards analysis showed that age [P<0.001, hazard ratio (HR) =1.29 (1.11-1.49)], race (P=0.018, HR =0.79 (0.65-0.96)], metastases [P=0.036, HR =1.96 (1.05-3.67)], tumor size [P=0.045, HR =1.35 (1.01-1.82)], degree of differentiation [P=0.002, HR =1.99 (1.30-3.06)] and metastasis surgery (P=0.005, HR =0.52 (0.33-0.82)] were independent prognostic factors for OS in GC patients with multi-organ metastases. The Fine-Gray competing risk analysis showed that age [P=0.006, HR =1.23 (1.06-1.42)], histological type [P=0.037, HR =1.53 (1.03-2.27)], metastases [P=0.009, HR =2.02 (1.19-3.41)], tumor size [P=0.028, HR =1.33 (1.03-1.70)], degree of differentiation [P=0.009, HR =1.65 (1.13-2.40)] and metastasis surgery [P=0.001, HR =0.50 (0.32-0.76)] were independent prognostic factors for CSS in GC patients with multi-organ metastases. The above factors were used to construct nomogram models for predicting OS and CSS. Both C-index and AUC of the training group and the validation group showed that the models had an acceptable predictive performance. The calibration curve showed that the predicted and ideal curves fit well, indicating that the constructed models were well-calibrated. Conclusions: Using data from the SEER database, this study established and validated nomogram models for OS and CSS in GC patients with multi-organ metastases, to help clinicians formulate accurate and individualized treatment plans. 2022 Translational Cancer Research. All rights reserved.

Introduction

Gastric cancer (GC) is one of the most common malignancies of the digestive system worldwide and was responsible for more than one million new cases and approximately 783,000 deaths in 2018 (1,2). It has high rates of fatality, recurrence, and distant metastasis (3). Due to its atypical early symptoms, most patients are already at an advanced stage when first diagnosed (4), and advanced GC patients are prone to distant metastases, where treatment options are limited. Several studies have shown that approximately 40% of newly diagnosed GC patients have synchronized distant metastases at the time of initial diagnosis (5,6). A recent study showed that compared with T1N+M1 GC patients, superficial GC patients with distant metastases without lymph node metastases had a worse prognosis because of more aggressive behaviors (7). The most common distant metastatic organs of GC are liver, lung and brain, etc. GC patients with multi-organ metastases have a lower survival rate and a worse prognosis, with a median survival time of only 6 months, which brings a huge burden to the family and society (8-10). Some previous studies showed that age, gender, race, primary tumor site, histological type, grade, T stage and N stage were prognostic factors for GC patients with distant metastases (11,12). To date, researchers have yet to explore the prognostic factors in GC patients with multi-organ metastases, so this study has clinical value, which can help improve the survival rate of patients and reduce the mortality rate (13).

A nomogram is a model that predicts the probability of a patient’s clinical events based on multivariate regression analysis, which can quickly and intuitively predict the prognosis of patients, and is widely used in tumor-related research (14). A key advantage of nomograms is their ability to estimate personalized risk based on patient and disease characteristics (15). By establishing a direct evaluation system, the predictive nomogram may also help patients and physicians make the optimal decisions for disease management (16). These nomograms will be part of the Cancer Survival Query System, a resource for clinicians and patients. The Surveillance, Epidemiology, and End Results (SEER) database registered by the National Cancer Institute is a recognized resource for information about the survival of cancer patients. The SEER website provides information including survival statistics divided by race, age, gender, disease stage and cancer site.

Currently, there has been published nomogram models for the prognosis of GC based on SEER (17), but researchers have not yet constructed a prognostic model for GC patients with multi-organ metastases. Therefore, in this study, we used the data extracted from SEER database to investigate prognostic factors of GC patients with multi-organ metastasis. Then we established a survival prediction model for such patients, to provide clinicians with a practical individualized prediction tool to optimize treatment plans and improve patient prognosis. As far as we know, this is the first predictive model to explore the prognosis of GC patients with multi-organ metastatic. We present the following article in accordance with the TRIPOD reporting checklist (available at https://tcr.amegroups.com/article/view/10.21037/tcr-21-2569/coif).

Methods

Data source

This study analyzed data obtained from searching the SEER database {SEER Research Data, 18 Registries, Nov 2020 Sub [2000-2018]}. All the data were downloaded using SEER*Stat 8.3.9 software (www.seer.cancer.gov). Data released by the SEER database has been granted patients’ informed consent and is freely available, so an ethical review and obtainment of informed consent were not required. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Inclusion and exclusion criteria

The inclusion criteria for our study were as follows: (I) patients diagnosed with GC between 2010 and 2018, inclusive; (II) age over 18 years; (III) pathological or imaging diagnosis of multi-organ metastasis (the SEER database only contains information about liver, lung, bone and brain metastases, so this study defined multi-organ metastasis as two or more metastases from those particular organs); (IV) a single primary tumor; (V) patient records clearly showing survival time and survival status. The exclusion criteria were as follows: (I) incomplete clinical or pathological data, such as age, race, histological grade, T stage, and surgical status; (II) patients diagnosed with multiple primary cancers; (III) metastases outside the liver, lung, bone, or brain, such as peritoneal metastasis or distant lymph node metastasis; (IV) a single metastatic site. After being screened according to these inclusion and exclusion criteria, the qualifying patients were randomly divided into a training group (n=972) and a validation group (n=414).

Factors included in the analysis

The research factors included in the analysis included age, gender, race, primary site, histological type, number of lymph nodes detected, metastases, tumor size, degree of differentiation, T stage, N stage, primary tumor surgery, metastasis surgery, survival time, survival status, and cause of death. Age was changed from a continuous variable to a categorical variable using X-tile software and divided into 2 groups: <70 and ≥70 years. Tumors were also classified into 2 size-related groups: <90 and ≥90 mm. Information about tumor size was drawn from the pathology report, operative report, imaging/radiographic techniques, and physical examination. The race parameter included black, white, and other races (American Indian/AK Native and Asian/Pacific Islander). The degree of tumor differentiation was divided into 4 groups, or grades, by the SEER database: well-differentiated (grade I), moderately differentiated (grade II), poorly differentiated (grade III), and undifferentiated (grade IV). We combined poorly differentiated and undifferentiated tumors into a single group (“grade III or IV”). Following the seventh and eighth edition of the American Joint Committee on Cancer (AJCC) guidelines, T stages were divided into T0–1, T2, T3, T4, and TX, and N stages were divided into N0, N1, N2, N3, and NX. The most common histological type of GC is adenocarcinoma; mucinous adenocarcinoma and signet-ring cell carcinoma are the most malignant, yet clinically common. Therefore, we divided the histological types into 4 groups, following the third edition of the International Classification of Diseases for Oncology (ICD-O-3): adenocarcinoma, mucinous adenocarcinoma, signet-ring cell carcinoma, and other types. The primary site code for the stomach (C16.0–C16.9) was based on the ICD-O-3 code. As only liver, lung, bone, and brain metastases are included in the SEER database, we classified the metastases into 11 groups according to the various combinations of 2 or more organs. Primary and metastasis surgeries were categorized as having been performed (yes) or not performed (no). Our main observation indicators were overall survival (OS) and cancer-specific survival (CSS). In this study, the most recent patient information was released in December 2018. Therefore, OS was defined as the number of months from GC diagnosis to either death from any cause or the end of follow-up (31 December, 2018). We defined CSS as the time from GC diagnosis to death from GC. The follow-up time for this study was 107 months.

Construction and verification of nomograms

After univariate and multivariate analysis, independent prognostic factors were determined using the forward stepwise selection method. In the univariate analysis, the variables were further analyzed using the Cox proportional hazards regression model when P<0.05. A nomogram was then constructed to predict the 15-, 19-, and 23-month OS for patients with multi-organ metastases of GC. Based on the Kaplan-Meier and log-rank test methods, the OS survival curves of the main variables were drawn. In addition, the Fine-Gray competing risk model was used to screen out independent risk factors affecting CSS in patients with GC and multi-organ metastasis. A nomogram was constructed to predict the 15-, 19-, and 23-month CSS for patients with multi-organ metastases of GC to evaluate the outcome after eliminating the impact of competing risk events. Based on the cumulative risk model and Fine-Gray test methods, the cumulative incidence curve of each main variable was drawn.

The performance of OS and CSS nomograms was evaluated by discrimination and calibration. Discrimination was evaluated using Harrell’s consistency index (C-index). Discrimination of the nomogram models was further evaluated using the area under the receiver operating characteristic (ROC) curve (AUC). Calibration plots between the predicted and actual values were drawn to evaluate the prediction efficiency of the nomograms.

Statistical analysis

In this study, statistical analysis was performed using the software packages X-tile 3.6.1, SPSS 26.0, and R 4.1.1. The test level was α=0.05, and a P value <0.05 indicated a statistically significant difference.

Results

Baseline characteristics of the training and validation group

After applying the inclusion and exclusion criteria, we collected data from a total of 1,386 GC patients with multi-organ metastasis from 2010 to 2018 and randomly grouped them into a training group (972 patients) and a verification group (414 patients), at a ratio of 7:3 (). The cutoff age value was 70 years (, low subset = blue, high subset = gray) and the cutoff value for tumor size was 90 mm (, low subset = blue, high subset = gray). There were no statistical differences in the basic characteristics of the 2 patient groups (P>0.05), and they were considered comparable ().

Figure 1

The inclusion flowchart of recruited GC patients with multi-organ metastases. GC, gastric cancer.

Figure 2

Determination of best-cutoff points of age and tumor size variables by the X-tile software.

Table 1

Baseline characteristics of the training and validation group

Variables	Training group (n=972)	Validation group (n=414)	Total (n=1,386)	P
Age, n (%)				0.708
<70 years	665 (68.4)	279 (67.4)	944 (68.1)
≥70 years	307 (31.6)	135 (32.6)	442 (31.9)
Sex, n (%)				0.938
Female	272 (28.0)	115 (27.8)	387 (27.9)
Male	700 (72.0)	299 (72.2)	999 (72.1)
Race, n (%)				0.230
Black	134 (13.8)	44 (10.6)	178 (12.8)
White	711 (73.1)	309 (74.6)	1,020 (73.6)
Other*	127 (13.1)	61 (14.7)	188 (13.6)
Primary site*, n (%)				0.763
Cardia	488 (50.2)	217 (52.4)	705 (50.9)
Fundus	59 (6.1)	26 (6.3)	85 (6.1)
Body	69 (7.1)	24 (5.8)	93 (6.7)
Antrum	86 (8.8)	34 (8.2)	120 (8.7)
Pylorus	10 (1.0)	3 (0.7)	13 (0.9)
Lesser	49 (5.0)	21 (5.1)	70 (5.1)
Greater	30 (3.1)	6 (1.4)	36 (2.6)
Overlapping	66 (6.8)	28 (6.8)	94 (6.8)
Stomach, NOS	115 (11.8)	55 (13.3)	170 (12.3)
Histological type*, n (%)				0.870
Adenocarcinoma	750 (77.2)	315 (76.1)	1,065 (76.8)
Mucinous adenocarcinoma	10 (1.0)	5 (1.2)	15 (1.1)
Signet ring cell carcinoma	106 (10.9)	43 (10.4)	149 (10.8)
Other	106 (10.9)	51 (12.3)	157 (12.3)
Regional nodes examined, n (%)				0.882
≤15	933 (96.0)	397 (95.9)	1,330 (96.0)
>15	3 (0.3)	2 (0.5)	5 (0.4)
Unknown	36 (3.7)	15 (3.6)	51 (3.7)
Metastases, n (%)				0.410
Liver + Lung	448 (46.1)	205 (49.5)	653 (47.1)
Liver + Bone	229 (23.6)	101 (24.4)	330 (23.8)
Liver + Brain	20 (2.1)	2 (0.5)	22 (1.6)
Lung + Bone	97 (10.0)	35 (8.5)	132 (9.5)
Lung + Brain	11 (1.1)	5 (1.2)	16 (1.2)
Bone + Brain	13 (1.3)	2 (0.5)	15 (1.1)
Liver + Lung + Bone	106 (10.9)	48 (11.6)	154 (11.1)
Liver + Lung + Brain	19 (2.0)	5 (1.2)	24 (1.7)
Liver + Bone + Brain	12 (1.2)	6 (1.4)	18 (1.3)
Lung + Bone + Brain	7 (0.7)	1 (0.2)	8 (0.6)
Liver + Lung + Bone + Brain	10 (1.0)	4 (1.0)	14 (1.0)
Tumor size, n (%)				0.197
<90 mm	356 (36.6)	149 (36.0)	505 (36.4)
≥90 mm	57 (5.9)	15 (3.6)	72 (5.2)
Unknown	559 (57.5)	250 (60.4)	809 (58.4)
Grade, n (%)				0.582
Grade I	28 (2.9)	9 (2.2)	37 (2.7)
Grade II	282 (29.0)	129 (31.2)	411 (29.7)
Grade III, Grade IV	662 (68.1)	276 (66.7)	938 (67.7)
Stage T, n (%)				0.561
T0-1	179 (18.4)	66 (15.9)	245 (17.7)
T2	28 (2.9)	13 (3.1)	41 (3.0)
T3	104 (10.7)	53 (12.8)	157 (11.3)
T4	167 (17.2)	79 (19.1)	246 (17.7)
TX	494 (50.8)	203 (49.0)	697 (50.3)
Stage N, n (%)
N0	292 (30.0)	108 (26.1)	400 (28.9)	0.137
N1	402 (41.4)	150 (36.2)	552 (39.8)	0.074
N2	49 (5.0)	24 (5.8)	73 (5.3)	0.564
N3	46 (4.7)	25 (6.0)	71 (5.1)	0.313
NX	183 (18.8)	107 (25.8)	290 (20.9)	0.003
Surg prim site, n (%)				0.397
No	943 (97.0)	398 (96.1)	1,341 (96.8)
Yes	29 (3.0)	16 (3.9)	45 (3.2)
Surg dis site, n (%)				0.413
No	946 (97.3)	406 (98.1)	1,352 (97.5)
Yes	26 (2.7)	8 (1.9)	34 (2.5)

Other*: American Indian/Alaska Native, Asian/Pacific Islander; Primary site*: C16.0-Cardia, NOS, C16.1-Fundus of stomach, C16.2-Body of stomach, C16.3-Gastric antrum, C16.4-Pylorus, C16.5-Lesser curvature of stomach NOS, C16.6-Greater curvature of stomach NOS, C16.8-Overlapping lesion of stomach, C16.9-Stomach, NOS; Histological type*: Adenocarcinoma: Adenocarcinoma other than mucinous adenocarcinoma and signet-ring cell carcinoma; Other: Squamous cell carcinoma, Carcinoid tumor, Adenosquamous carcinoma, Neuroendocrine carcinoma et al.

Prognostic factors affecting patient OS and CSS

Prognostic factors affecting the OS of patients

According to univariate analysis, age, race, metastasis, tumor size, degree of differentiation and metastasis surgery were prognostic factors affecting OS in GC patients with multi-organ metastasis (all P<0.05). Survival curves for the main variables were drawn based on the Kaplan-Meier and log-rank tests (). The results of multivariate analysis showed that being aged ≥70 years, being of black race, having liver+bone+brain metastasis, having a tumor size ≥90 mm, having a degree of differentiation of grade III or IV, and not undergoing metastasis surgery were independent prognostic factors for OS in GC patients with multi-organ metastasis (all P<0.05) ().

Figure 3

Survival curve of GC patients with multi-organ metastasis by age (A), race (B), metastasis (C), tumor size (D), grade (E) and metastasis surgery (F). GC, gastric cancer.

Table 2

Univariate and multivariate analysis of OS in gastric cancer patients with multi-organ metastases

Variables	Univariate analysis			Multivariate analysis
	HR (95% CI)	P value		HR (95% CI)	P value
Age
<70 years	Reference			Reference
≥70 years	1.25 (1.09–1.44)	0.002*		1.29 (1.11–1.49)	<0.001*
Sex
Female	Reference			–	–
Male	1.00 (0.87–1.16)	0.972		–	–
Race
Black	Reference			Reference
White	0.82 (0.67–0.99)	0.038*		0.79 (0.65–0.96)	0.018*
Other	0.85 (0.65–1.09)	0.203		0.82 (0.63–1.07)	0.146
Primary site
Cardia	Reference			–	–
Fundus	1.10 (0.83–1.46)	0.519		–	–
Body	1.06 (0.81–1.40)	0.671		–	–
Antrum	1.14 (0.89–1.46)	0.304		–	–
Pylorus	0.86 (0.45–1.67)	0.661		–	–
Lesser	0.86 (0.62–1.19)	0.353		–	–
Greater	1.29 (0.88–1.90)	0.187		–	–
Overlapping	1.20 (0.92–1.57)	0.170		–	–
Stomach, NOS	1.03 (0.84–1.28)	0.755		–	–
Histological type
Adenocarcinoma	Reference			–	–
Mucinous adenocarcinoma	1.32 (0.71–2.47)	0.380		–	–
Signet ring cell carcinoma	1.19 (0.96–1.48)	0.112		–	–
Other	1.12 (0.90–1.38)	0.316		–	–
Regional nodes examined
≤15	Reference			–	–
>15	0.81 (0.26–2.53)	0.723		–	–
Unknown	1.36 (0.97–1.92)	0.077		–	–
Metastases
Liver + Lung	Reference			Reference
Liver + Bone	0.89 (0.75–1.05)	0.162		0.85 (0.71–1.00)	0.054
Liver + Brain	1.17 (0.73–1.87)	0.525		1.25 (0.78–2.02)	0.357
Lung + Bone	1.09 (0.87–1.38)	0.447		0.98 (0.77–1.24)	0.844
Lung + Brain	0.75 (0.38–1.44)	0.383		0.74 (0.38–1.44)	0.374
Bone + Brain	1.38 (0.78–2.46)	0.271		1.22 (0.68–2.18)	0.503
Liver + Lung + Bone	0.94 (0.75–1.18)	0.603		0.93 (0.74–1.17)	0.540
Liver + Lung + Brain	1.26 (0.78–2.01)	0.345		1.58 (0.98–2.55)	0.062
Liver + Bone + Brain	1.37 (0.75–2.51)	0.303		1.96 (1.05–3.67)	0.036*
Lung + Bone + Brain	0.25 (0.09–0.67)	0.006*		0.41 (0.15–1.12)	0.081
Liver + Lung + Bone + Brain	1.47 (0.76–2.85)	0.252		1.35 (0.69–2.62)	0.379
Tumor size
<90 mm	Reference			Reference
≥90 mm	1.40 (1.05–1.88)	0.023*		1.35 (1.01–1.82)	0.045*
Unknown	1.14 (0.99–1.32)	0.061		1.08 (0.94–1.25)	0.275
Grade
Grade I	Reference			Reference
Grade II	1.44 (0.93–2.23)	0.103		1.25 (0.80–1.94)	0.324
Grade III, Grade IV	2.09 (1.37–3.21)	0.001*		1.99 (1.30–3.06)	0.002*
Stage T
T0-1	Reference			–	–
T2	0.98 (0.64–1.51)	0.937		–	–
T3	0.92 (0.71–1.18)	0.495		–	–
T4	1.05 (0.84–1.31)	0.658		–	–
TX	1.05 (0.88–1.26)	0.561		–	–
Stage N
N0	Reference			–	–
N1	0.98 (0.84–1.15)	0.818		–	–
N2	0.94 (0.69–1.27)	0.670		–	–
N3	1.23 (0.90–1.69)	0.195		–	–
NX	1.13 (0.92–1.39)	0.242		–	–
Surg prim site
No	Reference			–	–
Yes	0.78 (0.53–1.17)	0.232		–	–
Surg dis site
No	Reference			Reference
Yes	0.51 (0.33–0.78)	0.002*		0.52 (0.33–0.82)	0.005*

*, P<0.05. OS, overall survival; HR, hazard ratio; CI, confidence interval.

Prognostic factors affecting the CSS of patients

In this study, the Fine-Gray competing risk model was used to analyze CSS in GC patients with multi-organ metastases. In survival analysis, a competing risk is an event that precludes the primary event of interest and can easily lead to the miscalculation of the survival function (18). These events are said to be competing. In our study, death from other causes precluded the outcome of death from GC. Among the 972 patients in the training group, 880 died, of whom 846 were cancer-specific deaths and 34 were competing deaths. The results of stepwise screening based on the Fine-Gray competitive risk model and multivariate analysis showed that independent prognostic factors related to patient CSS were age, histological type, metastasis, tumor size, degree of differentiation, and metastasis surgery (all P<0.05) (). We also drew the cumulative incidence curve of the main variables ().

Table 3

Multivariate analysis of CSS in gastric cancer patients with multi-organ metastases

Variables	HR (95% CI)	P value
Age
<70 years	Reference
≥70 years	1.23 (1.06–1.42)	0.006*
Histological type
Adenocarcinoma	Reference
Mucinous adenocarcinoma	1.53 (1.03–2.27)	0.037*
Signet ring cell carcinoma	1.06 (0.86–1.31)	0.580
Other	0.93 (0.74–1.18)	0.560
Metastases
Liver + Lung	Reference
Liver + Bone	0.92 (0.79–1.09)	0.340
Liver + Brain	1.30 (0.91–1.84)	0.150
Lung + Bone	0.91 (0.71–1.16)	0.440
Lung + Brain	0.82 (0.47–1.41)	0.460
Bone + Brain	1.31 (0.88–1.96)	0.190
Liver + Lung + Bone	0.87 (0.69–1.09)	0.220
Liver + Lung + Brain	1.40 (0.96–2.03)	0.083
Liver + Bone + Brain	2.02 (1.19–3.41)	0.009*
Lung + Bone + Brain	0.53 (0.28–1.03)	0.060
Liver + Lung + Bone + Brain	1.44 (0.82–2.55)	0.210
Tumor size
<90 mm	Reference
≥90 mm	1.33 (1.03–1.70)	0.028*
Unknown	1.02 (0.89–1.18)	0.730
Grade
Grade I	Reference
Grade II	1.14 (0.77–1.68)	0.500
Grade III, Grade IV	1.65 (1.13–2.40)	0.009*
Surg dis site
No	Reference
Yes	0.50 (0.32–0.76)	0.001*

*, P<0.05. CSS, cancer-specific survival; HR, hazard ratio; CI, confidence interval.

Figure 4

Cumulative incidence curve of age (A), histological type (B), metastasis (C), tumor size (D), grade (E) and metastasis surgery (F) of GC patients with multi-organ metastases. 1, cancer-specific death; 2, competitive death; GC, gastric cancer.

Construction and verification of the nomogram prediction model

Construction of nomogram prediction models

Based on the selected independent prognostic factors affecting patient OS and CSS (), we constructed nomogram models to predict patient OS and CSS (). In the prediction models, the prognostic factors scores () were added together to obtain a total score, and the value corresponding to the total score could be used to predict a patient’s 15-, 19-, and 23-month OS and CSS.

Figure 5

Forest plot of all variables with hazard ratios in GC patients with multi-organ metastasis with OS (A) and CSS (B). GC, gastric cancer; OS, overall survival; CSS, cancer-specific survival.

Figure 6

Establishment of nomograms regarding both OS (A) and CSS (B). OS, overall survival; CSS, cancer-specific survival.

Table 4

Scores of prognostic factors in the OS and CSS nomograms

Variables	OS nomogram	CSS nomogram
Age
<70 years	0	0
≥70 years	16	18
Race
Black	15	–
White	0	–
Other	3	–
Histological type
Adenocarcinoma	–	2
Mucinous adenocarcinoma	–	31
Signet ring cell carcinoma	–	7
Other	–	0
Metastases
Liver + Lung	57	59
Liver + Bone	46	51
Liver + Brain	71	74
Lung + Bone	56	56
Lung + Brain	38	50
Bone + Brain	70	80
Liver + Lung + Bone	53	51
Liver + Lung + Brain	86	82
Liver + Bone + Brain	100	100
Lung + Bone + Brain	0	0
Liver + Lung + Bone + Brain	76	83
Tumor size
<90 mm	0	0
≥90 mm	19	21
Unknown	5	3
Grade
Grade I	0	0
Grade II	14	14
Grade III, Grade IV	44	43
Surg dis site
No	42	51
Yes	0	0

OS, overall survival; CSS, cancer-specific survival.

Verification of the nomogram prediction models

The bootstrap method was used to repeatedly sample 1,000 times to verify the modeling effect of the nomogram. The C-indices of the OS nomogram prediction model in the training and validation group were 0.608 [95% confidence interval (CI): 0.584–0.632] and 0.615 (95% CI: 0.582–0.648), respectively. We drew the ROC curve of the training group’s nomogram to predict the OS of patients at 15, 19, and 23 months, and calculated their AUC values to be 0.585, 0.573, 0.558, respectively (). At the same time, the AUC values of the verification group at 15, 19, and 23 months were 0.602, 0.719, and 0.719, respectively (). In addition, the prediction curve and the ideal curve in the calibration diagrams for the 2 groups fitted well, indicating that the model had good accuracy ().

Figure 7

ROC of nomogram using OS of the training group (A-C) and the validation group (D-F) in 15-, 19- and 23-month. TP, true positive; FP, false positive; ROC, receiver operating characteristic; OS, overall survival.

Figure 8

Evaluation of calibration plots based on OS of the training group (A-C) and the validation group (D-F) in 15-, 19- and 23-month. OS, overall survival.

The C-indices of the CSS nomogram prediction model for the training and validation groups were 0.605 (95% CI: 0.581–0.629) and 0.646 (95% CI: 0.622–0.670), respectively. The AUC values of the 15-, 19- and 23-month ROC curves of the CSS for the training group were 0.668, 0.679 and 0.667, respectively (), while in the validation group they were 0.781, 0.807 and 0.790, respectively (). In addition, the calibration curves of the 2 groups also showed a good fit between the predicted curve and the ideal curve, indicating that the model was well calibrated ().

Figure 9

ROC of nomogram and evaluation of calibration plots based on CSS of the training group (A,C) and the validation group (B,D) in 15-, 19- and 23-month. TP, true positive; FP, false positive; ROC, receiver operating characteristic; CSS, cancer-specific survival.

Discussion

Although GC is one of the most common malignant tumors, in the past decade, continuous improvements in diagnosis and treatment level have led to a decline in the incidence and mortality of GC (19). However, as the population ages, the number of cases of GC may increase. In addition, the onset of GC is insidious, leading to advanced disease stages at first diagnosis. Patients with advanced GC are prone to metastasis, leading to increased patient mortality (20). Therefore, it is of great significance to explore the prognostic factors for GC patients with multi-organ metastases.

A nomogram is a visual statistical graph used to predict the prognosis of various diseases. It can score each independent prognostic factor based on the results of multivariate analysis, and the sum of the points of each factor corresponds to the incidence of the endpoint event to predict the probability of the patient’s endpoint event (21). To date, researchers have constructed a nomogram model for predicting the prognosis of advanced GC, and confirmed that the model has good predictive ability (22,23). In the present study, a nomogram model was constructed to predict OS and CSS in GC patients with multi-organ metastasis. The model’s C-indices and AUC values were approximately 0.6 and 0.7, respectively. Considering that multi-organ metastatic tumors are affected by confounding variables, the predictive ability of the model constructed in this study could be considered relatively good. To our knowledge, this was the first study to highlight the potential value of nomograms using a population dataset of GC patients with multi-organ metastasis. By including the patient’s personal clinical and pathological information, the model constructed herein could briefly and intuitively predict the patient’s OS and CSS.

This study showed that the older the GC patients with multi-organ metastasis, the lower their OS and CSS. A previous study also showed that age was an independent prognostic factor for many cancer patients, including those with advanced GC (24). There are 2 possible factors for the low survival rate of elderly patients. First, elderly patients are prone to severe complications and have low immunity (25). Secondly, because younger patients are generally in better physical condition, they can easily tolerate extensive lymphadenectomy and standardized chemotherapy, while older patients cannot tolerate these treatments so easily (26). The results of this study indicated that tumor metastasis was an independent prognostic factor affecting OS and CSS in GC patients with multi-organ metastasis. Among the patients studied, those with liver, lung, bone and brain metastases had the least favorable prognoses, with a median survival time of only 1 month, consistent with the results of previous studies in hepatocellular carcinoma (27). This could be because patients with systemic metastases have almost passed the radical surgery window, and such patients have a heavier tumor burden and achieve poor chemotherapy effects, which leads to a very low survival rate. The nomogram of this study showed that race was an independent prognostic factor in OS for GC patients with multi-organ metastasis, which could increase the impact weight by 15 points. The poor prognosis of black people might be related to genetic differences, economic status, living habits and less frequent communication between doctors and nurses (28). In addition, this study also found that patients with mucinous adenocarcinoma had lower CSS compared with patients with other histological types. The impact weight in these cases was increased by 31 points, which might be due to the fact that patients with mucinous adenocarcinoma had developed to an advanced stage at the time of diagnosis and had a higher degree of differentiation (29).

Alshehri et al. (24) retrospectively analyzed the prognostic factors of 2,005 advanced GC patients and found that a tumor size >50 mm was an independent prognostic factor affecting OS rates. Similarly, the present study found that GC patients with multi-organ metastases with a tumor size of 90 mm or more had lower OS and CSS rates. This may be due to the larger tumors easily invading the serosal layer, which may lead to a later stage and poorer prognosis. our study used X-tile software to calculate that the optimum tumor size cut-off value affecting the prognosis of GC patients with multi-organ metastases was 90 mm, which could provide a certain reference value for future screening of patients (30). The degree of differentiation represents the relative malignancy of the tumor. The worse the differentiation, the higher the malignancy of the tumor. The present study also concluded that poor differentiation was a prognostic factor affecting patient OS and CSS, reaching an OS nomogram model score of 44 points and CSS nomogram model score of 43 points. The vast majority of patients (97.3%) in this study did not receive surgical metastasis treatment, which may be due to the fact that patients with 2 or more organ metastases rarely meet the requirements for radical resection. Our analysis found that the prognosis of patients who did not undergo metastasis surgery was unfavorable, with the OS nomogram showing an increase of 42 points in impact weight, and the CSS nomogram showing an increase of 51 points in impact weight. Previous studies have shown that metastasis surgery could prolong patient survival and improve prognosis (31,32). Therefore, when GC patients with multi-organ metastasis are in a physical condition that could tolerate surgery and the residual metastases function well, surgical treatment could be recommended to prolong the patient’s survival (33).

This study had several limitations. First, although the number of cases included in the SEER database was large and of high quality, it did not include factors affecting the prognosis of GC, such as smoking history, drinking history, family history, history of Helicobacter pylori infection and tumor markers, which caused bias. Secondly, since the SEER database has not provided relevant patient radiotherapy and chemotherapy information since 2015; considering that this information is affected by the subjective willingness of patients, the variables of radiotherapy and chemotherapy were not included in this study, which may also have caused bias. Thirdly, due to the incomplete clinicopathological information of some GC patients with multi-organ metastasis, we eliminated them from the study according to our inclusion and exclusion criteria, which may also have caused bias. Finally, this study constructed nomograms based on a retrospective analysis, and the level of research evidence was low, so our findings need further verification through prospective studies.

Conclusions

This study established a reliable prognostic risk assessment model for GC patients with multi-organ metastasis. These patients can find a total score based on the nomogram, and then evaluate their future survival, which has a strong clinical guiding role. For patients with high scores, clinicians can make treatment interventions and treatment plan adjustments in a timelier manner to improve the prognosis. We have confirmed the reliability and accuracy of this model through verification, and proved that the nomogram has a better prediction, and it can provide more information for clinical work.