
Citation: | Zhang Weifeng, Chen Han, Zhang Guoxin, Jin Guangfu. A nomogram for predicting lymph node metastasis in superficial esophageal squamous cell carcinoma[J]. The Journal of Biomedical Research, 2021, 35(5): 361-370. DOI: 10.7555/JBR.35.20210034 |
Superficial esophageal squamous cell carcinoma (SESCC) is defined as carcinoma with mucosal or submucosal invasion, regardless of regional lymph node (LN) metastasis[1]. The conventional curative treatment for SESCC is esophagectomy, although it is associated with considerable postoperative mortality[2]. Recently, endoscopic submucosal dissection (ESD), firstly reported by Japanese scholars to treat early gastrointestinal cancers[3], has become an alternative of esophagectomy for most non-metastatic SESCC cases of T1a stage. Several recent studies showed that the 3- and 5-year overall survival rates exceeded 90% in patients undergoing endoscopic resection with T1a-EP (epithelium) to T1b-SM1 (submucosal) cancer[4–5]. Thus, the indications of ESD for esophageal cancer have been gradually broadened considering the effectiveness and safety of endoscopic therapy[6–7]. According to the National Comprehensive Cancer Network (NCCN) guideline (Version 1, 2019), endoscopic resection is only applicable to T1a stage cancer while esophagectomy is indicated if the tumor infiltrates into the submucosal layer (T1b) or deeper[8]. But guidelines of the Japanese Esophageal Society (JES) and the Chinese Society of Clinical Oncology (CSCO) suggest if invasion depth is confined to the upper submucosal layer (T1b, SM1), ESD is still considered as the main strategy[9–10]. Currently, most Chinese patients with a fairly good general condition prefer minimally invasive endoscopic resection to surgical treatment, and therefore, ESD enjoys wider popularity in China[5]. However, ESD is deemed applicable only in patients without LN metastasis (LNM) because endoscopic treatment could not replace lymphadenectomy[11]. The LN status is not only a key factor to determine the training strategy, but also the most important prognostic factor in esophageal cancer[12]. Therefore, it is necessary for clinicians to establish a model for predicting the risk of LNM in SESCC patients so as to decide the therapeutic indication of endoscopic treatment.
Although such models have been reported in recent years, the methodology applied was usually suboptimal. Most previously established models were non-parametric[13–14] and/or insufficiently validated[13–17]. And some were established with a less-than-300 sample size[14–17]. The aim of this study is to establish an updated and validated nomogram for predicting LNM in patients with SESCC, which can not only help evaluate cancer prognosis but also provide guidance for intervention selection.
A retrospective cohort was established in patients receiving surgical resection and lymphadenectomy for SESCC between December 2009 and June 2018 in the First Affiliated Hospital of Nanjing Medical University (Nanjing, China). Inclusion criteria were: (1) histopathological diagnosis of esophageal squamous carcinoma on surgical specimens; (2) pT1 stage carcinoma (no tumor invasion beyond the submucosa); (3) surgical resection and at least two-field lymphadenectomy. Exclusion criteria were: (1) history of previous malignancies and anticancer therapies; (2) histopathological diagnosis of esophageal adenocarcinoma or other types of esophageal cancer; (3) mixed type of esophageal cancer; (4) tumor with undefined pathological origin and metastatic esophageal cancer; (5) esophagectomy after endoscopic resection; (6) aged <18 years; (7) perioperative mortality; and (8) no contrast-enhanced computed tomography (CT) scans, endoscopic ultrasonography (EUS) or fluorodeoxyglucose positron emission tomography (FDG-PET) before surgery. Following the same inclusion and exclusion criteria, an independent validation cohort of 203 patients was prospectively studied from June 2018 to January 2019. The study was approved by the Institutional Ethics Committee of the First Affiliated Hospital of Nanjing Medical University (the register number, 2018-SR-152). Informed consent of surgical risk and anesthesia services was obtained before surgery.
All participants undergoing endoscopic biopsy were assessed histopathologically to define esophageal cancer or high-grade intraepithelial neoplasia before surgery. And all must have at least one of the following radiographic results for preoperative LN status evaluation: contrast-enhanced CT scans, EUS or FDG-PET. Two experienced radiologists assessed the results of CT and FDG-PET scans, and two sonographers reviewed EUS results. Interpretations were given independently and disagreements were resolved by consensus. Positive radiographic results were defined as CT-, EUS-, or PET-reported positive LNM.
General clinical features, including age and sex, were collected. All the participants were divided into two age groups: 18 to 59 years and ≥60 years. The pathologic diagnosis of esophageal squamous cell carcinoma was performed by two experienced pathologists who assessed the specimens independently. In case of discordant outcome, the specimen was evaluated by a third independent expert pathologist. The diagnosis was determined only if at least two pathologists agreed.
Pathologic features, including tumor size, invasion depth, microscopic type, histologic grade, lymph-vessel invasion (LVI), and LN status (LNM, location of LNM, the total number of resected LN, and the number of metastatic LNs during surgery), were collected. The tumor size (the maximum diameter) was measured in two dimensions by Vernier's calipers. Tumor location is defined by the position of epicenter of tumor. If the statement of epicenter was unavailable, the following measurements were applied: (1) upper: 15 to 24 cm from incisors; (2) middle: 25 to 29 cm from incisors; (3) low: 30 to 40 cm from incisors. Histologic grade (G) was categorized as well-differentiated (G1), moderately differentiated (G2), and poorly differentiated (G3)[18]. The invasion depth was classified into four categories: epithelium (EP)/lamina propria mucosa (LPM), muscularis mucosa (MM), submucosal (SM)1, and SM2 or deeper. Macroscopic tumor type was classified by using the 2016 Japanese Classification of Esophageal Cancer (11th edition): superficial type, protruding type, ulcerative and localized type, ulcerative and infiltrative type, diffusely infiltrative type, and unclassified type[19].
Significant predictors for establishing nomogram were identified by IBM SPSS Statistics, Version 20.0 (SPSS, USA). The univariate logistic regression and chi-square analyses were applied to identify the significance of each variable in the training cohort. Variables significantly associated with LNM (P≤0.05) were identified as candidate predictors for multivariate logistic regression analysis. The optimal cutoff value was assessed by Youden index in the receiver operating characteristic (ROC) curve. The area under receiver operating characteristic curves (AUC) was calculated to evaluate the diagnostic accuracy and compared with DeLong's test. The net reclassification improvement (NRI) and the integrated discrimination improvement (IDI) were further applied to quantify the refinement in predictive accuracy[20].
Multivariate analysis was performed by using the "glmnet" package in R software (version 3.5.1, http://www.R-project.org/). The results of multivariate analysis were integrated into R software to formulate a nomogram by using the "rms" package. The validation of nomogram was measured by concordance index (C-index)[21–23], calibration curve, and decision curve analysis (DCA)[24–26]. NRI and IDI were calculated by using the "PredictABEL" package. The statistical significance level was two-sided and statistical significance was considered as P<0.05.
The establishment of the model consisted of five steps: (1) analyzing logistic regression for candidate predictor selection; (2) developing a parametric model (nomogram); (3) evaluating model performance by calibration plot, discrimination, and decision-curve analysis; (4) establishing a prospective cohort for external validation; (5) comparing different models by multiple statistical methods.
Fig. 1 shows the flowchart of generating the training and validation cohorts. A total of 8984 LNs were resected from 711 patients (510 [71.7%] males; median age, 61.2 years; range, 36 to 83) by surgery. The number of patients receiving CT, FDG-PET, and EUS is 711 (100%), 78 (11.0%), and 63 (8.9%), respectively. The baseline characteristics of both training and validation cohorts are shown in Supplementary Table 1 (available online). The overall incidence of LNM was 13.8% (112/711) in the training cohort and 17.2% (35/203) in the validation cohort, respectively. Of the 112 patients with metastatic LN in the training cohort, 1546 LNs were resected (average, 16.8; range 8 to 51) and 186 LNs were found positive for metastasis. The average tumor size was (1.50±1.2) cm. A total of 265 (37.3%) patients were defined as T1a stage and 446 (62.7%) patients were T1b stage. In T1a, the number of patients with EP/LPM and MM invasion was 110 and 155, respectively. In T1b, the number of patients with SM1 and SM2 or deeper invasion was 352 and 94, respectively.
The LNM rate of T1a and T1b tumors were 2.3% (6/265) and 23.8% (106/446), respectively. Positive radiographic results were as follows: 245 were CT positive (including 4 CT+PET/CT+EUS positive, and 9 CT+PET/CT positive), and 3 were PET/CT positive.
The sensitivity and specificity of positive radiographic results (CT-, FDG-PET- or EUS- reported positive LNM) were 59.0% and 59.1%, respectively. The incidences of the upper, middle, and lower thoracic esophagus cancer were 9.8%, 40.1%, and 50.1%, respectively. Supplementary Table 2 (available online) summarizes the LN status of the resected LNs according to regions and N Stage.
Table 1 summarizes the characteristics of patients according to lymph node status in both training and validation cohorts. In the training group, patients' age, sex, tumor location, and macroscopic type were not significantly associated with LNM. The optimal applicable cut-off value of tumor size was defined as 1.55 cm by calculating Youden's index. The area under the ROC curve was 0.70 (95% confidence interval [CI], 0.651–0.747) by using univariate logistic regression.
Variable | Training cohort | P-value | Validation cohort | P-value | ||
LNM+ (n [%]) | LNM− (n [%]) | LNM+ (n [%]) | LNM− (n [%]) | |||
Sex | 0.938 | 0.193 | ||||
Male | 80 (15.7) | 430 (84.3) | 30 (19.1) | 127 (80.9) | ||
Female | 32 (15.9) | 169 (84.1) | 5 (10.9) | 41 (89.1) | ||
Age | 0.596 | 0.548 | ||||
Media, range (years) | 61.2 (36–83) | 61.7 (42–79) | ||||
≥60 | 64 (16.4) | 326 (83.6) | 23 (16.2) | 119 (83.8) | ||
<60 | 48 (15.0) | 273 (85.0) | 12 (19.7) | 49 (80.3) | ||
Tumor location | 0.562 | 0.597 | ||||
Upper thoracic | 14 (20.0) | 56 (80.0) | 16 (18.0) | 73 (82.0) | ||
Middle thoracic | 45 (15.8) | 240 (84.2) | 13 (14.8) | 75 (85.2) | ||
Lower thoracic | 53 (14.9) | 303 (85.1) | 6 (23.1) | 20 (76.9) | ||
Tumor size | <0.001 | 0.007 | ||||
Mean, range (cm) | 1.50±1.20 (0.3–7.0) | 1.60±1.10 (0.5–3.5) | ||||
≥1.55 | 78 (26.7) | 214 (73.3) | 26 (23.9) | 83 (76.1) | ||
<1.55 | 34 (8.1) | 385 (91.9) | 9 (9.6) | 85 (90.4) | ||
Histologic grade | <0.001 | <0.001 | ||||
Well-differentiated | 2 (2.9) | 68 (97.1) | 0 (0.0) | 29 (100.0) | ||
Moderately-differentiated | 50 (11.3) | 391 (88.7) | 16 (11.7) | 121 (88.3) | ||
Poorly-differentiated | 60 (30.0) | 140 (70.0) | 19 (51.4) | 18 (48.6) | ||
Invasion depth | <0.001 | <0.001 | ||||
EP/LPM | 0 (0.0) | 110 (100.0) | 0 (0.0) | 20 (100.0) | ||
MM | 6 (3.9) | 149 (96.1) | 2 (6.7) | 28 (93.3) | ||
SM1 | 50 (14.2) | 302 (85.8) | 15 (12.4) | 106 (87.6) | ||
SM2 or deeper | 56 (59.6) | 38 (40.4) | 16 (53.3) | 14 (46.7) | ||
Lymph-vessel invasion | <0.001 | <0.001 | ||||
Positive | 28 (65.1) | 15 (34.9) | 17 (85.0) | 3 (15.0) | ||
Negative | 84 (12.6) | 584 (87.4) | 18 (9.8) | 165 (90.2) | ||
Pathological type | 0.312 | 0.411 | ||||
Protruding | 34 (19.3) | 142 (80.7) | 11 (24.4) | 34 (75.6) | ||
Superficial type | 29 (13.2) | 191 (86.8) | 13 (16.2) | 67 (83.8) | ||
Ulcerative and localized | 36 (16.5) | 182 (83.5) | 9 (19.1) | 38 (80.9) | ||
Infiltrative | 11 (16.4) | 56 (83.6) | 2 (8.0) | 23 (92.0) | ||
Diffusely infiltrative | 2 (6.7) | 28 (93.3) | 0 (0.0) | 6 (100.0) | ||
Radiographic results | <0.001 | 0.009 | ||||
Positive for LNM | 66 (21.2) | 245 (78.8) | 14 (29.8) | 33 (70.2) | ||
Negative for LNM | 46 (11.5) | 354 (88.5) | 21 (13.5) | 135 (86.5) | ||
ignificant predictors for establishing nomogram were identified by IBM SPSS Statistics. The chi-square analysis was applied to identify the significance of each variable in the training cohort. P-values were derived from the chi-square analyses between the single variables and LN status. *P<0.05. LNM: lymph node metastasis; LNM+: patient with LNM; LNM−: patient without LNM; EP: epithelium; LPM: lamina propria mucosa; MM: muscularis mucosa; SM: submucosal. |
The multivariate logistic regression were applied and confirmed that tumor size ≥1.55 cm (odds ratio [OR], 3.75; 95% CI, 2.265–6.2; P<0.001), submucosal invasion (OR, 9.03; 95% CI, 3.78-21.551; P<0.001), poorly-differentiated histologic type (OR 2.81; 95% CI, 1.731–4.5556; P<0.001), positive lymph-vessel invasion (OR, 11.64; 95% CI, 5.271–25.72; P<0.001), and positive radiographic results (OR, 2.3; 95% CI, 1.412–3.735; P=0.001) were significantly associated with LNM in SESCCs (Supplementary Fig. 1 and Supplementary Table 3, available online).
Next, estimated through stratification of significant risk factors, the probability of LNM with one, two, three, and more than three risk factors were 5.1%, 12.1%, 37.9%, and 78.7%, respectively (Supplementary Table 4, available online).
A nomogram was established with five significant candidate predictors: tumor size ≥1.55 cm, submucosal invasion, poorly-differentiated histologic type, positive lymph-vessel invasion, and positive radiographic results (Fig. 2).
Three key measurements were performed to evaluate model performance: calibration plot, discrimination, and decision-curve analysis. A new cohort consisting of 203 patients was formed for external validation. The calibration curve of nomogram presented good agreement in predicting probability and actual LNM rate in the training cohort, with a non-significant Hosmer-Lemeshow test statistic (P=0.476). Fig. 3 shows the results of calibration curve in predicting model. The calibration curves of nomogram were highly consistent with actual observation for LN status. ROC curves were further constructed and the area under curve exceeded 0.8 in all established cohorts (Fig. 4). Internal validation presented good discrimination with a C-index of 0.855 (95% CI, 0.820–0.890; sensitivity, 83.9%; specificity, 71.3%). In external validation, good calibration and discrimination were also achieved with a C-index of 0.820 (95% CI, 0.782–0.859; sensitivity, 82.9%; specificity, 94.6%). Cross-validation was conducted in the training cohort, which was randomly divided into two testing groups: 494 (70%) and 217 (30%). The AUC was 0.923 (95% CI, 0.846–0.998) in testing group 1 (70%) and 0.872 (95% CI, 0.835–0.909) in testing group 2 (30%). The predicted accuracy was 87.2% in testing group 1 and 84.3% in testing group 2.
To simplify the model and further validate model accuracy, another model was set up by subtracting predictive candidates: a simplified model integrating invasion depth, tumor size, and LVI. The calibration curves for the probability of LNM in the training and validation cohorts demonstrated the agreement between prediction and observation in the simplified model (Fig. 4A and B). The simplified model presented discrimination with a C-index of 0.916 (95% CI, 0.858–0.973; sensitivity, 78.6%; specificity, 75.0%) in training cohort while a C-index of 0.827 (95% CI, 0.748–0.906; sensitivity, 48.6%; specificity, 98.2%) in validation cohort (Fig. 4).
According to DeLong's test, the training model was superior to the simplified model (Z=3.654, P<0.001) in the training cohort. To further quantify the refinement in predictive accuracy, NRI and IDI were applied. The results showed that the training model was significantly improved in prediction performance when compared with the simplified one (NRI, 19.3%, P<0.001; IDI, 6.75%, P<0.001). Table 2 summarizes the results of the two models.
Variable | Traning model | Simplified model | |||||
β | OR (95% CI) | P-value | β | OR (95% CI) | P-value | ||
Tumor size | 1.321 | 3.75(2.265–6.200) | <0.001 | 1.335 | 3.80 (2.334–6.188) | <0.001 | |
Invasion depth | 2.200 | 9.03 (3.780–21.551) | <0.001 | 2.291 | 9.88 (4.177–23.388) | <0.001 | |
LV invasion | 2.455 | 11.64 (5.271–25.720) | <0.001 | 2.430 | 11.36 (5.425–23.779) | <0.001 | |
Histologic grade | 1.033 | 2.81 (1.731–4.556) | 0.015 | – | – | – | |
Radiographic results | 0.832 | 2.30 (1.412–3.735) | – | – | – | – | |
Constant | −5.207 | – | – | −4.492 | – | – | |
C-index | |||||||
Training cohort | 0.855 (0.820–0.890) | 0.820 (0.782–0.859) | |||||
Validation cohort | 0.916 (0.858–0.973) | 0.827 (0.748–0.906) | |||||
Sensitivity | 0.839 | 0.768 | |||||
Specificity | 0.713 | 0.750 | |||||
Multivariate analysis was performed by using the "glmnet" package in R software. β represents the regression coefficient. The validation of models was measured by concordance index (C-index), The statistical significance level was two-sided and statistical significance was considered as P<0.05. Data of histologic grade was modified by using a binary-tiered grading system in which well-differentiated and moderately differentiated tumors were combined into one category and poorly-differentiated ones as another category. SESCCs: superficial esophageal squamous cell carcinoma; OR: Odds ratio; CI: confidence interval; LV: lymph-vessel. |
Well-fitted DCA was achieved in the training group and better performance was found than that in another model (Fig.5). The results indicated that if the threshold probability of LNM ranged from 10% to 80%, the application of the nomogram could yield more benefit than either the treat-all-patients or the treat-none strategies. Net benefit between the training model overlapped at several points while net benefit of training model was superior to the simplified model almost within the whole range.
In the current study, a nomogram was developed and validated as an easy-to-use tool for clinicians to predict the possibility of LNM in SESCC patients. Our study showed the rates of LNM were 13.8% in the training cohort and 17.2% in the validation cohort, which was consistent with previously published data[27–28]. The sample size in the training cohort was twice larger than the previously published models, which ensured a higher precision of the present model[29–30].
T1 substage, LVI, tumor differentiation, and tumor size have been indicated as predictors of LNM in several studies[31], which were also verified in our study. The presence of LNM was observed in 28 of 43 (65.1%) and 17 of 20 (85%) of superficial esophageal cancer with vascular invasion in the training and validation cohorts, respectively, which was significantly more frequent than in cancers without lymph-vessel invasion in these two cohorts (12.6% and 9.8%, respectively). Some other studies also suggested that patients with lymph-vessel invasion were more likely to present LNM[32]. Meanwhile, tumor lymph-vessel invasion was considered as an independent risk factor for LNM of superficial esophageal cancer (OR=4.06).
The pathological type of well-differentiated tumors is commonly considered to be close to mature tissues with slow growth and low risk of LNM, while the poorly differentiated tumors are featured by rapid growth and high risk of LNM. Thus tumor differentiation is also considered an independent risk factor for LNM of superficial esophageal cancer[12]. In our study, the rates of LNM in patients with well-, moderately-, and poorly-differentiated superficial esophageal cancer were 2.9% (2/70), 11.3% (50/441), and 30% (60/200) in training cohort, and 0% (0/29), 11.7% (16/137), 51.4% (19/37) in validation cohort, respectively. In this study, only 2 cases with LNM were found in patients with well-differentiated tumors in the training cohort, and no LNM was found in the patients with well-differentiated tumors in validation cohort, which confirmed that tumor differentiation was an independent risk factor for lymph node metastasis (OR=3.2).
Submucosa is considered as the boundary layer of esophageal cancer cell metastasis[33]. The depth of tumor invasion was another independent risk factor associated with lymph node metastasis. Previous studies reported LNM rates in MM1, MM2, MM3, SM1, SM2 and SM3 were 0%, 1.5%–3.7%, 5.3%–30.8%, 8.7%–42.1%, 12.7%–40.7% and 28.4%–66.7%, respectively[29, 33–39]. Our results showed that the LNM rate of patients with superficial esophageal cancer confined to mucosa was much lower than that of tumors confined to submucosa (2.3% and 23.8%, respectively, P<0.001) in the training cohort, which was further confirmed in the validation cohort (4.0% and 21.6%, respectively, P<0.001). The current study also confirmed tumor infiltration depth as an independent risk factor for LNM (OR=8.71), which is in accordance with some other previous reports[15].
Wang et al found that the probability of lymph node metastasis and the prognosis are closely associated with tumor size[40]. In this study, tumor diameter showed a significant association with LNM. The rates of LNM were 8.1% and 26.7% of tumors with diameter <1.55 cm and ≥1.55 cm in the training cohort, respectively, and similar results were obtained in the validation cohort (9.6% and 23.9%, respectively).
We also included the preoperative radiographic result as a predictor in the model. These modalities, as the principal non-invasive examination to determine clinical nodal status by regional lymph node imaging before interventions, have similar performance in clinical nodal status assessment[41]. Therefore, we defined preoperative radiographic results as CT-, EUS- or FDG-PET-reported positive LNM. Although the sensitivity and specificity were only 59.0% and 59.1%, this radiographic predictor was still statistically significant (P<0.001) as an independent risk factor for LN status. Thus, we included this factor in the final model.
The previous studies usually classify tumors into two categories: low-grade (well-differentiated) tumors, and high-grade (moderately- and poorly- differentiated) ones[16–17]. We unexpectedly detected that the moderately-differentiated tumors were not statistically significant in multivariate logistic regression. Hence, we adopted a binary-tiered grading system: well-moderately differentiated (low grade) and poorly differentiated (high grade). This grading method improves precision of prognostic significance, especially for the poorly- and moderately-differentiated tumors.
To simplify the model, another model was set by subtracting the predictive candidates. Although calibration plot of the training model was not superior to the simplified model, the improved C-index and NRI/IDI analysis of the training model demonstrated a better performance in discrimination. In the external validation, adequate discrimination was achieved by using the comprehensive model (C-index, 0.860), which was then improved in the validation cohort (C-index, 0.920). Considering the comparable baseline characteristics between the two cohorts, the enhanced discrimination implied good performance of this model by external validation results. However, after validation, we paradoxically detected that when compared with the simplified model, the training model in both training and validation cohorts presented a relatively poor calibration (potential of over-estimating the risk) but good discrimination. Since calibration and discrimination were both significant references in model evaluation, it was difficult to tell which one is more important. It is reported that calibration is more crucial in prognostic settings while discrimination is more meaningful in diagnosis[42]. Hence, we gave preference to good discrimination because patients with true LNM would always have higher predicted risks than those without. Although the relatively poor calibration plot might overestimate the risk, it could still distinguish the patients with LNM from those without, which could result in a stricter and more indicative ESD procedure as such risks would not be under-estimated. Thus, the training model is more clinically useful than the other two. Next, a decision curve has been applied to evaluate the clinical applicability of the nomogram system. The results indicated the net benefit of the training model is almost superior to that of the simplified model within the whole range of 10%–80% (the threshold probability).
The limitations of our study include: 1) this is a single-center retrospective study, and multi-center study with external validation is expected in the future; 2) a selection bias may exist; 3) the molecular-biological or immunohistochemical characteristics are not included in our model; 4) there are differences in some predictive indicators between validation and training cohort, necessitating a larger sample size to further validate the efficacy of this model.
In conclusion, we have developed a comprehensive nomogram integrating both clinically and statistically significant factors. This model may facilitate the prediction of LNM in patients with SESCC.
This work was supported by the National Natural Science Foundation of China for studies on early-stage GI neoplasms (Grant No. 81072032 and No. 81470830).
[1] |
Wang GQ, Jiao GG, Chang FB, et al. Long-term results of operation for 420 patients with early squamous cell esophageal carcinoma discovered by screening[J]. Ann Thorac Surg, 2004, 77(5): 1740–1744. doi: 10.1016/j.athoracsur.2003.10.098
|
[2] |
Jacobs M, Macefield RC, Elbers RG, et al. Meta-analysis shows clinically relevant and long-lasting deterioration in health-related quality of life after esophageal cancer surgery[J]. Qual Life Res, 2014, 23(4): 1097–1115. doi: 10.1007/s11136-013-0545-z
|
[3] |
Hirao M, Masuda K, Asanuma T, et al. Endoscopic resection of early gastric cancer and other tumors with local injection of hypertonic saline-epinephrine[J]. Gastrointest Endosc, 1988, 34(3): 264–269. doi: 10.1016/S0016-5107(88)71327-9
|
[4] |
Berger A, Rahmi G, Perrod G, et al. Long-term follow-up after endoscopic resection for superficial esophageal squamous cell carcinoma: a multicenter Western study[J]. Endoscopy, 2019, 51(4): 298–306. doi: 10.1055/a-0732-5317
|
[5] |
Qi ZP, Chen T, Li B, et al. Endoscopic submucosal dissection for early esophageal cancer in elderly patients with relative indications for endoscopic treatment[J]. Endoscopy, 2018, 50(9): 839–845. doi: 10.1055/a-0577-2560
|
[6] |
Min YW, Lee H, Song BG, et al. Comparison of endoscopic submucosal dissection and surgery for superficial esophageal squamous cell carcinoma: a propensity score-matched analysis[J]. Gastrointest Endosc, 2018, 88(4): 624–633. doi: 10.1016/j.gie.2018.04.2360
|
[7] |
Tsujii Y, Nishida T, Nishiyama O, et al. Clinical outcomes of endoscopic submucosal dissection for superficial esophageal neoplasms: a multicenter retrospective cohort study[J]. Endoscopy, 2015, 47(9): 775–783. doi: 10.1055/s-0034-1391844
|
[8] |
Ajani JA, D'Amico TA, Almhanna K, et al. Esophageal and esophagogastric junction cancers, version 1.2015[J]. J Natl Compr Canc Netw, 2015, 13(2): 194–227. doi: 10.6004/jnccn.2015.0028
|
[9] |
Kuwano H, Nishimura Y, Oyama T, et al. Guidelines for diagnosis and treatment of carcinoma of the esophagus April 2012 edited by the Japan Esophageal Society[J]. Esophagus, 2015, 12(1): 1–30. doi: 10.1007/s10388-014-0465-1
|
[10] |
Wang FH, Shen L, Li J, et al. The Chinese society of clinical oncology (CSCO): clinical guidelines for the diagnosis and treatment of gastric cancer[J]. Cancer Commun, 2019, 39(1): 1–31. doi: 10.1186/s40880-019-0349-9
|
[11] |
Cho JW, Choi SC, Jang JY, et al. Lymph node metastases in esophageal carcinoma: an endoscopist's view[J]. Clin Endosc, 2014, 47(6): 523–529. doi: 10.5946/ce.2014.47.6.523
|
[12] |
Rice TW, Ishwaran H, Hofstetter WL, et al. Esophageal cancer: associations with (pN+) lymph node metastases[J]. Ann Surg, 2017, 265(1): 122–129. doi: 10.1097/SLA.0000000000001594
|
[13] |
Ancona E, Rampado S, Cassaro M, et al. Prediction of lymph node status in superficial esophageal carcinoma[J]. Ann Surg Oncol, 2008, 15(11): 3278–3288. doi: 10.1245/s10434-008-0065-1
|
[14] |
Ma DW, Jung DH, Kim JH, et al. Predicting lymph node metastasis for endoscopic resection of superficial esophageal squamous cell carcinoma[J]. J Thorac Cardiovasc Surg, 2019, 157(1): 397–402. doi: 10.1016/j.jtcvs.2018.07.034
|
[15] |
Zhou Y, Du J, Wang Y, et al. Prediction of lymph node metastatic status in superficial esophageal squamous cell carcinoma using an assessment model combining clinical characteristics and pathologic results: A retrospective cohort study[J]. Int J Surg, 2019, 66: 53–61. doi: 10.1016/j.ijsu.2019.04.014
|
[16] |
Zheng H, Tang H, Wang H, et al. Nomogram to predict lymph node metastasis in patients with early oesophageal squamous cell carcinoma[J]. Br J Surg, 2018, 105(11): 1464–1470. doi: 10.1002/bjs.10882
|
[17] |
Min BH, Yang JW, Min YW, et al. Nomogram for prediction of lymph node metastasis in patients with superficial esophageal squamous cell carcinoma[J]. J Gastroenterol Hepatol, 2020, 35(6): 1009–1015. doi: 10.1111/jgh.14915
|
[18] |
Rice TW, Patil DT, Blackstone EH. 8th edition AJCC/UICC staging of cancers of the esophagus and esophagogastric junction: application to clinical practice[J]. Ann Cardiothorac Surg, 2017, 6(2): 119–130. doi: 10.21037/acs.2017.03.14
|
[19] |
Japan Esophageal Society. Japanese classification of esophageal cancer, 11th edition: part II and III[J]. Esophagus, 2017, 14(1): 37–65. doi: 10.1007/s10388-016-0556-2
|
[20] |
Pencina MJ, D'Agostino RB Sr, D'Agostino RB Jr, et al. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond[J]. Stat Med, 2008, 27(2): 157–172. doi: 10.1002/sim.2929
|
[21] |
Friedman J, Hastie T, Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent.[J]. Journal of Statistical Software, 2010, 33(1): 1–22. https://pubmed.ncbi.nlm.nih.gov/20808728/
|
[22] |
S Kundu, YS Aulchenko, CM van Duijn, ACJW Janssens. PredictABEL: an R package for the assessment of risk prediction models[J]. Eur J Epidemiol, 2011, 26: 261–4. doi: 10.1007/s10654-011-9567-4
|
[23] |
Feng Z, McLerran D, Grizzle J. A comparison of statistical methods for clustered data analysis with Gaussian error[J]. Stat in Med, 1996, 15: 1793–1806. https://pubmed.ncbi.nlm.nih.gov/8870161/
|
[24] |
Lee YH, Bang H, Kim DJ. How to establish clinical prediction models[J]. Endocrinol Metab, 2016, 31(1): 38–44. doi: 10.3803/EnM.2016.31.1.38
|
[25] |
Steyerberg EW, Vergouwe Y. Towards better clinical prediction models: seven steps for development and an ABCD for validation[J]. Eur Heart J, 2014, 35(29): 1925–1931. doi: 10.1093/eurheartj/ehu207
|
[26] |
Alba AC, Agoritsas T, Walsh M, et al. Discrimination and calibration of clinical prediction models: users' guides to the medical literature[J]. JAMA, 2017, 318(14): 1377–1384. doi: 10.1001/jama.2017.12126
|
[27] |
Siewert JR, Stein HJ, Feith M, et al. Histologic tumor type is an independent prognostic parameter in esophageal cancer: lessons from more than 1,000 consecutive resections at a single center in the Western world[J]. Ann Surg, 2001, 234(3): 360–369. doi: 10.1097/00000658-200109000-00010
|
[28] |
Stein HJ, Feith M, Bruecher BLDM, et al. Early esophageal cancer: pattern of lymphatic spread and prognostic factors for long-term survival after surgical resection[J]. Ann Surg, 2005, 242(4): 566–573. doi: 10.1097/01.sla.0000184211.75970.85
|
[29] |
Aoyama J, Kawakubo H, Mayanagi S, et al. Discrepancy between the clinical and final pathological findings of lymph node metastasis in superficial esophageal cancer[J]. Ann Surg Oncol, 2019, 26(9): 2874–2881. doi: 10.1245/s10434-019-07498-2
|
[30] |
Xue L, Ren L, Zou S, et al. Erratum: parameters predicting lymph node metastasis in patients with superficial esophageal squamous cell carcinoma[J]. Mod Pathol, 2015, 28(1): 161. doi: 10.1038/modpathol.2014.133
|
[31] |
Suzuki G, Yamazaki H, Aibe N, et al. Endoscopic submucosal dissection followed by chemoradiotherapy for superficial esophageal cancer: choice of new approach[J]. Radiat Oncol, 2018, 13(1): 246. doi: 10.1186/s13014-018-1195-7
|
[32] |
Yamada T, Sugiyama H, Ochi D, et al. Risk factors for submucosal and lymphovascular invasion in gastric cancer looking indicative for endoscopic submucosal dissection[J]. Gastric Cancer, 2014, 17(4): 692–696. doi: 10.1007/s10120-013-0323-1
|
[33] |
Shen W, Shen Y, Tan L, et al. A nomogram for predicting lymph node metastasis in surgically resected T1 esophageal squamous cell carcinoma[J]. J Thorac Dis, 2018, 10(7): 4178–4185. doi: 10.21037/jtd.2018.06.51
|
[34] |
Ishihara R, Oyama T, Abe S, et al. Risk of metastasis in adenocarcinoma of the esophagus: a multicenter retrospective study in a Japanese population[J]. J Gastroenterol, 2017, 52(7): 800–808. doi: 10.1007/s00535-016-1275-0
|
[35] |
Mitobe J, Ikegami M, Urashima M, et al. Clinicopathological investigation of lymph node metastasis predictors in superficial esophageal squamous cell carcinoma with a focus on evaluation of lympho-vascular invasion[J]. Scand J Gastroenterol, 2013, 48(10): 1173–1182. doi: 10.3109/00365521.2013.832365
|
[36] |
Yachida T, Oda I, Abe S, et al. Risk of lymph node metastasis in patients with the superficial spreading type of esophageal squamous cell carcinoma[J]. Digestion, 2020, 101(3): 239–244. doi: 10.1159/000499017
|
[37] |
Xue L, Ren L, Zou S, et al. Parameters predicting lymph node metastasis in patients with superficial esophageal squamous cell carcinoma[J]. Mod Pathol, 2012, 25(10): 1364–1377. doi: 10.1038/modpathol.2012.89
|
[38] |
Zhou Y, Du J, LI H, et al. Clinicopathologic analysis of lymph node status in superficial esophageal squamous carcinoma[J]. World J Surg Oncol, 2016, 14(1): 259. doi: 10.1186/s12957-016-1016-0
|
[39] |
Jia R, Luan Q, Wang J, et al. Analysis of predictors for lymph node metastasis in patients with superficial esophageal carcinoma[J]. Gastroenterol Res Pract, 2016, 2016: 3797615. doi: 10.1155/2016/3797615
|
[40] |
Wang BY, Goan YG, Hsu PK, et al. Tumor length as a prognostic factor in esophageal squamous cell carcinoma[J]. Ann Thorac Surg, 2011, 91(3): 887–893. doi: 10.1016/j.athoracsur.2010.11.011
|
[41] |
Choi J, Kim SG, Kim JS, et al. Comparison of endoscopic ultrasonography (EUS), positron emission tomography (PET), and computed tomography (CT) in the preoperative locoregional staging of resectable esophageal cancer[J]. Surg Endosc, 2010, 24(6): 1380–1386. doi: 10.1007/s00464-009-0783-x
|
[42] |
Cook NR. Statistical evaluation of prognostic versus diagnostic models: beyond the ROC curve[J]. Clin Chem, 2008, 54(1): 17–23. doi: 10.1373/clinchem.2007.096529
|
1. | Wang J, Zhang X, Gan T, et al. Risk factors and a predictive nomogram for lymph node metastasis in superficial esophageal squamous cell carcinoma. World J Gastroenterol, 2023, 29(47): 6138-6147. DOI:10.3748/wjg.v29.i47.6138 |
Variable | Training cohort | P-value | Validation cohort | P-value | ||
LNM+ (n [%]) | LNM− (n [%]) | LNM+ (n [%]) | LNM− (n [%]) | |||
Sex | 0.938 | 0.193 | ||||
Male | 80 (15.7) | 430 (84.3) | 30 (19.1) | 127 (80.9) | ||
Female | 32 (15.9) | 169 (84.1) | 5 (10.9) | 41 (89.1) | ||
Age | 0.596 | 0.548 | ||||
Media, range (years) | 61.2 (36–83) | 61.7 (42–79) | ||||
≥60 | 64 (16.4) | 326 (83.6) | 23 (16.2) | 119 (83.8) | ||
<60 | 48 (15.0) | 273 (85.0) | 12 (19.7) | 49 (80.3) | ||
Tumor location | 0.562 | 0.597 | ||||
Upper thoracic | 14 (20.0) | 56 (80.0) | 16 (18.0) | 73 (82.0) | ||
Middle thoracic | 45 (15.8) | 240 (84.2) | 13 (14.8) | 75 (85.2) | ||
Lower thoracic | 53 (14.9) | 303 (85.1) | 6 (23.1) | 20 (76.9) | ||
Tumor size | <0.001 | 0.007 | ||||
Mean, range (cm) | 1.50±1.20 (0.3–7.0) | 1.60±1.10 (0.5–3.5) | ||||
≥1.55 | 78 (26.7) | 214 (73.3) | 26 (23.9) | 83 (76.1) | ||
<1.55 | 34 (8.1) | 385 (91.9) | 9 (9.6) | 85 (90.4) | ||
Histologic grade | <0.001 | <0.001 | ||||
Well-differentiated | 2 (2.9) | 68 (97.1) | 0 (0.0) | 29 (100.0) | ||
Moderately-differentiated | 50 (11.3) | 391 (88.7) | 16 (11.7) | 121 (88.3) | ||
Poorly-differentiated | 60 (30.0) | 140 (70.0) | 19 (51.4) | 18 (48.6) | ||
Invasion depth | <0.001 | <0.001 | ||||
EP/LPM | 0 (0.0) | 110 (100.0) | 0 (0.0) | 20 (100.0) | ||
MM | 6 (3.9) | 149 (96.1) | 2 (6.7) | 28 (93.3) | ||
SM1 | 50 (14.2) | 302 (85.8) | 15 (12.4) | 106 (87.6) | ||
SM2 or deeper | 56 (59.6) | 38 (40.4) | 16 (53.3) | 14 (46.7) | ||
Lymph-vessel invasion | <0.001 | <0.001 | ||||
Positive | 28 (65.1) | 15 (34.9) | 17 (85.0) | 3 (15.0) | ||
Negative | 84 (12.6) | 584 (87.4) | 18 (9.8) | 165 (90.2) | ||
Pathological type | 0.312 | 0.411 | ||||
Protruding | 34 (19.3) | 142 (80.7) | 11 (24.4) | 34 (75.6) | ||
Superficial type | 29 (13.2) | 191 (86.8) | 13 (16.2) | 67 (83.8) | ||
Ulcerative and localized | 36 (16.5) | 182 (83.5) | 9 (19.1) | 38 (80.9) | ||
Infiltrative | 11 (16.4) | 56 (83.6) | 2 (8.0) | 23 (92.0) | ||
Diffusely infiltrative | 2 (6.7) | 28 (93.3) | 0 (0.0) | 6 (100.0) | ||
Radiographic results | <0.001 | 0.009 | ||||
Positive for LNM | 66 (21.2) | 245 (78.8) | 14 (29.8) | 33 (70.2) | ||
Negative for LNM | 46 (11.5) | 354 (88.5) | 21 (13.5) | 135 (86.5) | ||
ignificant predictors for establishing nomogram were identified by IBM SPSS Statistics. The chi-square analysis was applied to identify the significance of each variable in the training cohort. P-values were derived from the chi-square analyses between the single variables and LN status. *P<0.05. LNM: lymph node metastasis; LNM+: patient with LNM; LNM−: patient without LNM; EP: epithelium; LPM: lamina propria mucosa; MM: muscularis mucosa; SM: submucosal. |
Variable | Traning model | Simplified model | |||||
β | OR (95% CI) | P-value | β | OR (95% CI) | P-value | ||
Tumor size | 1.321 | 3.75(2.265–6.200) | <0.001 | 1.335 | 3.80 (2.334–6.188) | <0.001 | |
Invasion depth | 2.200 | 9.03 (3.780–21.551) | <0.001 | 2.291 | 9.88 (4.177–23.388) | <0.001 | |
LV invasion | 2.455 | 11.64 (5.271–25.720) | <0.001 | 2.430 | 11.36 (5.425–23.779) | <0.001 | |
Histologic grade | 1.033 | 2.81 (1.731–4.556) | 0.015 | – | – | – | |
Radiographic results | 0.832 | 2.30 (1.412–3.735) | – | – | – | – | |
Constant | −5.207 | – | – | −4.492 | – | – | |
C-index | |||||||
Training cohort | 0.855 (0.820–0.890) | 0.820 (0.782–0.859) | |||||
Validation cohort | 0.916 (0.858–0.973) | 0.827 (0.748–0.906) | |||||
Sensitivity | 0.839 | 0.768 | |||||
Specificity | 0.713 | 0.750 | |||||
Multivariate analysis was performed by using the "glmnet" package in R software. β represents the regression coefficient. The validation of models was measured by concordance index (C-index), The statistical significance level was two-sided and statistical significance was considered as P<0.05. Data of histologic grade was modified by using a binary-tiered grading system in which well-differentiated and moderately differentiated tumors were combined into one category and poorly-differentiated ones as another category. SESCCs: superficial esophageal squamous cell carcinoma; OR: Odds ratio; CI: confidence interval; LV: lymph-vessel. |