Abstract
Background:
Endoscopic resection has been introduced as an alternative treatment for superficial adenocarcinoma of the esophagogastric junction (AEG), but is limited by positive nodal status. We aimed to investigate the predictors of lymph node metastasis (LNM) in patients with Siewert type II T1 AEG.
Methods:
The Surveillance, Epidemiology, and End Results (SEER) database was used to identify eligible patients with Siewert type II T1 AEG. The prevalence of LNM was assessed. Logistic regression analysis with multivariable adjustment was used to determine predictors of LNM. We also performed Cox regression analysis to examine the prognostic value of LNM, which was further confirmed by competing risk analysis and cumulative incidence function (CIF).
Results:
In total, 2651 patients with T1 AEG were included, with a median age of 69 years and a median follow-up of 28 months. The overall prevalence of LNM was 17.2% in T1 AEG. When stratified by tumor invasion depth, the prevalence of LNM was 8.5% for intramucosal tumors and 22.6% for submucosal tumors. Adjusted logistic regression analysis showed that age, sex, tumor grade, tumor size and tumor infiltration depth were independent predictors of LNM in T1 AEG. Multivariate Cox regression analysis revealed that positive nodal status was significantly associated with worse overall survival and cancer-specific survival (CSS). Subgroup analysis consistently demonstrated that patients with LNM had significantly poorer CSS than those without LNM in most subgroups. Finally, the CIF was calculated, showing that patients with LNM had a significantly higher cancer-specific death rate than those without LNM.
Conclusions:
This population-based study identified age, sex, tumor grade, tumor infiltration depth and tumor size as independent predictors of LNM in T1 AEG. Considering the high prevalence of LNM in T1 AEG, endoscopic resection for curative aims may only be introduced in patients without high risks of LNM.
Introduction
Although adenocarcinoma of the esophagogastric junction (AEG) is uncommon, its incidence has been rapidly increasing over time globally. 1 -4 The incidence of AEG increased by approximately 2.5-fold from the early 1970s to the early 1990s according to the statistics from the National Cancer Institute’s Surveillance, Epidemiology, and End Results (SEER) program. 2 Similarly, a Japanese cohort of consecutive patients with gastric adenocarcinoma revealed that the overall proportion of AEG increased from 2.3% (1962-1965) to 10.0% (2001-2005). 3 The survival of AEG patients is generally poor and might vary greatly depending on regional lymph node involvement and distant metastasis. 5
AEG is commonly considered as a separate tumor entirety of digestive tract cancer. 6,7 Due to its special anatomical location, the classification of AEG has been historically complicated. Siewert classified AEG into 3 subgroups based on the anatomical location of the tumor epicenter relative to the esophagogastric junction (EGJ): Siewert type I (5 to 1 cm above the EGJ), type II (1 cm above to 2 cm below the EGJ), and type III (2 to 5 cm below the EGJ). 8,9 Among the 3 subtypes, Siewert type II is generally considered as the true cardia carcinoma arising from EGJ. 9,10
Despite the overall poor survival of patients with AEG, the prognosis of patients with superficial lesions is relatively favorable if curative resection is performed. Superficial AEG has been traditionally managed with surgical resection in most cases, 11 mainly including radical esophagectomy and lymphadenectomy. 12 Although radical surgery is conventionally linked with secure long-term outcomes, it also has several drawbacks. 13 First, the surgical procedure may increase the risk of overtreatment in mucosal or submucosal tumor lesions without high risks of local recurrence and distant metastasis. Additionally, the in-hospital mortality rate of esophagectomy is reported to be as high as 5.0%, 5 and operative resection of the gastric cardia and postoperative complications can diminish the quality of life of patients. 14
Endoscopic resection, a minimally invasive technique, has been increasingly propagated as a reliable treatment option for superficial AEG when properly adopted based on rigorous indication criteria. 15,16 In a retrospective study enrolling 53 patients with superficial AEG who underwent endoscopic submucosal dissection (ESD), Yamada et al reported that the cause-specific survival rate was 100%, without recurrence or metastasis among patients after curative resection (median follow-up: 6.1 years). 17 Another retrospective Japanese study from 13 centers revealed similar findings, showing that the 5-year cause-specific survival rate was 100% among superficial AEG patients with a low risk of lymph node metastasis (LNM). 15 Endoscopic resection can eliminate superficial cancer that is confined to the primary site, while it cannot be used to cure cancers with regional lymph node involvement or even distant metastasis. 18 Among patients with superficial AEG who were treated with endoscopic mucosal resection (EMR) or ESD, the 5-year overall survival (OS) was 93.9% in patients with a low risk of LNM, which sharply dropped to approximately 80% among those with a high risk of LNM. 15 Therefore, the accurate prediction of lymph node involvement is an essential prerequisite for the success of endoscopic resection, which is also of great significance in pretreatment decision making.
The SEER database, an authoritative source of cancer data in the US, records and reports cancer incidence and survival data by covering approximately 28% of the total US population. 19,20 By providing information on patient-specific and tumor-specific characteristics, the SEER database is particularly useful for studying uncommon malignant tumors.
In the present study, the SEER database was used to assess the prevalence of LNM in T1 AEG (mucosal and submucosal tumor lesions) and to identify the predictors of LNM in T1 AEG.
Materials and Methods
Patient Selection
We performed this retrospective study by retrieving relevant data from the SEER database. Although detailed information on the Siewert classification of AEG (type I, II or III) was not directly available in the SEER database, we were still able to specifically identify Siewert type II AEG according to 2 parameters. Cancers simultaneously satisfying 2 conditions (“TNM 7/CS v0204 + Schema” encoded 28 (Esophagus GE Junction) and “Primary Site-Labeled” encoded 160 (Cardia, NOS)) were extracted and classified as Siewert type II AEG. 2,21,22
We downloaded the data of patients diagnosed with AEG from 2004 to 2015 from the SEER database. For this analysis, the inclusion criteria were as follows: 1) age at diagnosis of 18 years or older (in consideration of the extremely small proportion of patients under 18 years and the large proportion of older patients); 2) pathologically diagnosed T1 Siewert Type II AEG; 3) available lymph node status; 4) active follow-up; and 5) first or the only 1 primary malignancy. Patients were excluded if they had in situ cancer. Moreover, patients with distant metastasis and those with survival times less than 1 month were also excluded. The details of patient selection are depicted in Figure 1.

Flowchart of patient selection.
Patient demographics (age, sex, race, year of diagnosis and marital status), tumor characteristics (tumor grade, tumor size, T stage, N stage and number of lymph nodes examined), treatment regimens and patient survival were collected from the SEER database for subsequent analysis. Since the SEER database is publicly available and the data are de-identified, the requirement for approval was waived by the local ethics committee.
Statistical Analysis
Eligible patients were divided into N-negative and N-positive groups according to their regional lymph node status. The Chi-square test or Fisher’s exact test was used to test the independence of the clinicopathological categorical variables. Predictors of LNM were assessed and identified by unadjusted and adjusted logistic regression models as well as backward logistic regression model. Odds ratios (ORs) along with 95% confidence intervals (CIs) were calculated. Afterward, a nomogram model was generated based on the independent LNM predictors identified from the adjusted logistic regression analysis. The performance of the nomogram-based prediction of LNM risk was evaluated by a calibration curve. A receiver operating characteristic (ROC) curve was plotted to assess the predictive accuracy of the nomogram model.
In this study, the primary endpoints included OS and cancer-specific survival (CSS). The former was defined as the duration from the cancer diagnosis to death from any cause, while the latter referred to the period between the date of diagnosis and the date of death attributed to this type of cancer. Survival curves for both OS and CSS were generated by the Kaplan-Meier method. The difference between survival curves was evaluated by the log-rank test. We further applied a Cox regression model to identify independent prognostic factors for OS and CSS. Finally, in consideration of both oncological and non-oncological risks among tumor patients, competing risk analysis was performed, and the cumulative incidence function (CIF) was calculated. 23,24
SPSS version 26.0 (SPSS Inc., Chicago, IL, USA) and R software for Mac version 3.6.1 (The R Foundation for Statistical Computing, Vienna, Austria) were used to analyze data and to plot figures. The level of statistical significance was set at 2-sided P values < 0.05.
Results
Demographic and Clinicopathological Characteristics of Patients
The detailed process of patient selection was shown in Figure 1. Among the 21482 patients with Siewert II AEG diagnosed between 2004 and 2015, 2651 eligible patients were finally enrolled based on the inclusion and exclusion criteria. The median age was 69 years [interquartile range (IQR), 60-78], and the median follow-up was 28 months (IQR, 11-63). Most patients were male (76.8%) and white (88.8%). Overall, 457 of 2651 patients (17.2%) had LNM. Table 1 summarized the patient demographics and clinicopathological characteristics.
Clinicopathological Characteristics of Patients.
Abbreviation: IQR, interquartile range.
Independent Predictors of Lymph Node Metastasis
Adjusted multivariable logistic regression was performed to identify the risk factors for LNM. The results showed that age, sex, tumor grade, depth of tumor invasion and tumor size were significant predictive factors for LNM in T1 AEG (Table 2). Interestingly, a decreased LNM risk was detected in older patients [OR = 0.64 (age: 66-80 years), OR = 0.40 (age: over 80 years), both P < 0.05]. The risk of LNM was attenuated in female patients [OR = 0.59 (0.44-0.79), P < 0.001]. Regarding oncological factors, patients with moderately-differentiated [OR = 1.62 (1.02-2.70), P = 0.049] and poorly-differentiated [OR = 3.10 (1.96-5.11), P < 0.001] AEG had a higher risk of LNM than those with well-differentiated lesions. Tumor invasion depth was also significantly associated with LNM risk. Patients with submucosal tumors had a 2.28-fold higher risk of LNM than those with mucosal lesions. Compared with patients with small tumor lesions sized ≤ 1 cm, the risk of LNM was significantly increased in those with tumor sizes exceeding 1 cm [OR = 2.22 (1.1-2 cm), OR = 4.48 (2.1-4 cm), OR = 5.81 (> 4 cm), all P < 0.05]. Moreover, the backward logistic regression model robustly showed that age, sex, tumor grade, tumor invasion and tumor size were independent predictors of LNM in T1 AEG (Table 2).
Logistic Regression Analysis of Risk Factors for Lymph Node Metastasis in T1 AEG.
Abbreviations: AEG, adenocarcinoma of the esophagogastric junction; OR, odds ratio; 95% CI, 95% confidence intervals.
To better visualize and present the risk factors of LNM, we further constructed a nomogram model (Supplementary Figure 1A), which could be used to estimate the numerical probability for a specific individual by integrating these parameters. To assess the performance of the nomogram model, a calibration curve was constructed and showed good agreement between the nomogram-predicted risks and the actual risks of LNM (C index: 0.742) (Supplementary Figure 1B). In addition, to assess the predictive capacity of the nomogram for predicting LNM risk, a ROC curve was plotted. As shown in Supplementary Figure 1C, the area under the curve (AUC) of the ROC curve was 0.742 (95% CI: 0.718-0.765).
Lymph Node Metastasis and Patient Survival
Unadjusted and adjusted Cox regression models were adopted to investigate the prognostic significance of LNM. Our findings showed that age, marital status, tumor grade, tumor size, lymph node status, local treatment and radiation were significant prognostic factors for both OS (Table 3) and CSS (Table 4) in patients with T1 AEG. Tumor invasion was significantly associated with CSS but not OS in the adjusted Cox regression analysis. Kaplan-Meier curves were further plotted to depict the survival in patients stratified by lymph node status. As shown in Figure 2A and B, OS and CSS rates were significantly decreased in patients with positive LNM compared with those without LNM (both P < 0.0001).
Cox Regression Analysis of Prognostic Factors for Overall Survival in T1 AEG.
Abbreviations: AEG, adenocarcinoma of the esophagogastric junction; HR, hazard ratio; 95% CI, 95% confidence intervals.
Cox Regression Analysis of Prognostic Factors for Cancer-Specific Survival in T1 AEG.
Abbreviations: AEG, adenocarcinoma of the esophagogastric junction; HR, hazard ratio; 95% CI, 95% confidence intervals.

Prognostic value of lymph node metastasis in T1 AEG. Kaplan Meier curves for (A), Overall survival and (B), cancer-specific survival in patients stratified by lymph node status. (C) Forest plot for subgroup analysis. The patients were divided into subgroups according to different clinicopathological characteristics. Patients with Lymph Node Metastasis (LNM) had significantly worse CSS than those without LNM in most subgroups.
We further performed survival analysis in subsets of patients according to different clinicopathological features. As shown in Figure 2C, patients with positive lymph node status had significantly poorer CSS than those without LNM in most subgroups. Therefore, subgroup analysis further demonstrated that LNM was an independent prognostic factor for patients with T1 AEG.
Competing Risk Analysis
The long-term survival outcomes of cancer patients are affected by oncological factors and non-oncological factors. During follow-up, patients might die from other causes, such as cardiovascular disease and car accidents, before the occurrence of cancer-specific death. 24,25 To accurately reveal the prognostic value of LNM in T1 AEG, a competing risk model was applied for a direct and exact interpretation of the effects of risk factors on the cause-specific cumulative incidence of death. 26 Multivariate analysis showed that age over 80 years [subdistribution hazard ratio (SHR) = 1.51, P = 0.006], unmarried status (SHR = 1.26, P < 0.001), poorly differentiated/undifferentiated tumor grade (SHR = 1.31, P = 0.045), submucosal lesion (SHR = 1.23, P = 0.037), tumor size over 2 cm, positive nodal status (SHR = 1.51, P < 0.001) and the administration of local treatment were all significant prognostic factors for CSS (Table 5). Finally, the CIF was calculated to elucidate the probability of cancer-specific death and death attributable to other causes. 27 The results showed that patients with LNM had a significantly higher cancer-specific death rate than those without LNM (P < 0.01) (Figure 3).

Cumulative incidence function for cancer-specific death in T1 AEG. The red curve indicates the cancer-specific death in patients with Lymph Node Metastasis (LNM), and the black curve suggests the cancer-specific death in patients without LNM.
Competing Risk Analysis for Cancer-Specific Death.
Abbreviations: SHR, subdistribution hazard ratio; 95% CI, 95% confidence intervals.
Discussion
Endoscopic resection of early-stage AEG based on rigorous indication criteria and complete resection of the tumor is advantageous. In addition to the comparable long-term clinical outcomes between endoscopic resection and surgical resection, 15,17,28,29 endoscopic treatment has been widely introduced due to the dramatically decreased postoperative morbidity and significantly increased quality of life.
Our population-based analysis revealed that the overall risk of LNM in patients with T1 AEG was relatively high (17.2%). In addition, the prevalence of LNM in patients with intramucosal cancer was 8.5% (85 out of 998), which sharply rose to 22.6% in those with submucosal tumors (182 out of 807). Consistent with our findings, by analyzing 453 patients with T1 AEG who underwent surgical resection between 2004 and 2010, Dubecz et al 12 previously reported that the prevalence of LNM was 9.5% for T1a tumors and 22.9% for T1b tumors. Similarly, the proportion of LNM was 22.2% in patients with T1 adenocarcinoma of the esophagus and EGJ, according to an Australian study from 1985 to 2003. 5
In consideration of the decisive role of LNM in choosing endoscopic or surgical resection, we further examined the predictors of LNM in patients with T1 AEG. Age, sex, tumor grade, depth of tumor infiltration and tumor size were identified as significant predictors of lymph node involvement in T1 AEG. In a previous population-based study, Dubecz et al investigated the predictors of LNM in patients with pT1 carcinoma of the esophagus and the gastric cardia who underwent surgical resection. 12 Similarly, their study revealed that tumor infiltration of the submucosa, large tumor size (exceeding 1 cm) and poor tumor differentiation were independent predictors of LNM.
In our study, we showed that poorly differentiated or undifferentiated tumor grades increased the LNM risk by 3 times, compared to well-differentiated tumor grade. In a Chinese cohort involving 393 AEG patients who underwent radical resection and lymphadenectomy, tumor differentiation was also an independent influencing factor for LNM. 30 Poor differentiation indicates higher tumor heterogeneity, resulting in more aggressive biological characteristics compared to well and moderately differentiated tumors. 7 Tumor size has also been considered to be associated with LNM risk in AEG. Gross tumor volume detected by multidetector computed tomography is an independent risk factor for LNM. 31 In addition, gross tumor volume can also be used to differentiate negative lymph nodes from positive lymph nodes, 31 which can also assist preoperative clinical decision making. Tumor infiltration depth is another important predictor of lymph node positivity in superficial esophageal carcinoma. 12,32 As expected, we also found that tumor infiltration depth was an independent predictor of LNM in T1 AEG.
Intriguingly, we found that the LNM risk was significantly decreased in older patients. Specifically, the risk of LNM in patients aged 66-80 years and over 80 years significantly dropped to 0.64 and 0.40 (both P < 0.05), respectively, compared to that in patients aged 50 years or under. In a single-center study involving 137 AEG patients, the rate of LNM was slightly but not significantly higher in younger patients (71.9%) than in older patients (64.8%). 33 The insignificant statistical outcomes might be due to the relatively small sample size. Similarly, several studies have reported decreased LNM risk in older populations with other malignancies of the digestive system, 34,35 but not in AEG. Although relevant studies have indicated that the overall poor prognosis of young patients might be due to a more aggressive biological process of the tumor, 36 the underlying causes for the decreased LNM risk in older patients still remain elusive. Further studies are warranted to explore the intrinsic associations between age and LNM in AEG. Our findings indicate for the first time that endoscopic resection might be the optimal option for low-risk older T1 AEG patients who are at high risk of perioperative and postoperative complications if surgically treated. In addition, we also showed that female patients had a significantly lower risk of LNM than male patients (OR = 0.59, P < 0.001). The incidence of AEG shows a male predominance. 37 Siewert and Stein reported that the male-to-female ratio was 4.9 in Siewert type II AEG. 38 Similarly, the male-to-female ratio was 3.1 in patients with Siewert type II T1 AEG in our study. Although there are few studies investigating the male predominance and its possible role on patient prognosis in AEG, several studies have shown longer survival in females than in males with esophageal cancer, 39 -42 which is caused by both sex itself (sex hormones and reproductive factors) and other extrinsic risk factors for mortality. 39 Thus, it is intriguing to examine the possible underlying causes of sex differences in the incidence and prognosis of AEG in the future.
In the survival analysis, we found that age, marital status, tumor grade, tumor size, lymph node status, local treatment and radiation were significant prognostic factors for both OS and CSS. The cancer-specific death rate was 35.8% in patients without LNM, which dramatically increased to 59.5% in nodal-positive patients (Table 1). The above findings indicate that lymph node status definitely plays a vital role in the outcomes attributable to cancer.
A distinctive feature of follow-up in cancer patients is that their survival is threatened by both oncological and non-oncological factors. For instance, a patient who dies of a car accident is no longer at risk of death attributable to cancer. These non-oncological events are called competing events. 25 In our study, we performed competing risk analysis and calculated the CIF to accurately determine the prognostic significance of lymph node status. Consequently, lymph node involvement was robustly associated with poor survival in T1 AEG.
We investigated the predictors of LNM in T1 AEG by enrolling 2651 eligible patients from the SEER database. With a median follow-up of 28 months and a relatively large sample size, our findings can achieve a high degree of statistical power. However, certain limitations still exist. In addition to the studied predictors of LNM, lymphovascular invasion has been reported as an independent predictor of LNM. 43 However, we could not assess these factors due to the limited data available in the SEER database.
Conclusions
In summary, in this population-based study, we report that the prevalence of LNM is as high as 17.2% among patients with T1 AEG. Age, sex, tumor grade, depth of tumor infiltration and tumor size are independent predictors of LNM in T1 AEG. Despite the obvious advantage of endoscopic resection, decision making in the treatment of T1 AEG should be cautiously approached according to the individualized conditions of patients. Endoscopic resection might be reasonably performed among T1 nodal-negative AEG patients with well-differentiated, small-size and mucosa-confined lesions, especially in the elder population, for curative aims.
Supplemental Material
Supplemental Material, sj-tif-1-ccx-10.1177_10732748211026668 - Predictors of Lymph Node Metastasis in Siewert Type II T1 Adenocarcinoma of the Esophagogastric Junction: A Population-Based Study
Supplemental Material, sj-tif-1-ccx-10.1177_10732748211026668 for Predictors of Lymph Node Metastasis in Siewert Type II T1 Adenocarcinoma of the Esophagogastric Junction: A Population-Based Study by Liubo Chen, Kejun Tang, Sihan Wang, Dongdong Chen and Kefeng Ding in Cancer Control
Footnotes
Abbreviations
AEG, adenocarcinoma of esophagogastric junction; CI, confidence interval; CIF, cumulative incidence function; CSS, cancer-specific survival; EGJ, esophagogastric junction; EMR, endoscopic mucosal resection; ESD, endoscopic submucosal dissection; IQR, interquartile range; LNM, lymph node metastasis; OR, odd ratio; OS, overall survival; SEER, Surveillance, Epidemiology, and End Results; SHR, subdistribution hazard ratio.
Authors’ Note
Liubo Chen, MD, Kejun Tang, PhD, and Sihan Wang, PhD contributed equally to this work. Since SEER database is publicly available and de-identified, the requirement for approval was waived by the local ethics committee in this retrospective study.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by Key Technology Research and Development Program of Zhejiang Province (No.2017C03017), National Natural Science Foundation of China (No. 81802750) and Zhejiang Provincial Natural Science Foundation of China (No. LQ19H160043).
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
