Risk factors of regional lymph node (RLN) metastasis among patients with bone sarcoma and survival of patients with RLN-positive bone sarcoma.

BACKGROUND: Regional lymph node metastasis (RLNM) has been reported to be a prognostic factor for poor survival outcomes of bone sarcoma. However, studies about risk factors for RLNM of bone sarcoma are extremely rare, and the outcome of such patients remains to be explored. We aimed to identify risk factors for RLNM of bone sarcoma and conduct survival analysis for patients with bone sarcoma with RLNM.
METHODS: A total of 10,641 patients confirmed of malignant bone sarcomas from 1983 to 2014 were identified from the Surveillance, Epidemiology, and End Results (SEER) database, with 311 being regional lymph node positive. Logistic regression analysis was used to identify risk factors for RLNM, while the Cox proportional hazards model and the Fine and Gray's regression model were used for survival analysis.
RESULTS: The proportion of RLNM was 6.0% in Ewing sarcoma, 2.5% in osteosarcoma and 1.1% in chondrosarcoma. Other bone tumors together had a RLNM rate of 4.2%. Risk factors identified by the logistic regression analysis for RLNM were male patients, primary tumor site, tumor type and size. The multivariate Cox regression analysis suggested age, race, distant metastasis, tumor type and surgical treatment to be prognostic factors for the overall survival of patients with RLNM. Taking non-cancer-specific death as a competing risk, however, we found only age between 30-60 years [sub-distribution hazard ratio (SHR), 1.528, 95% CI, 1.028-2.271; P=0.02], distant metastasis (SHR, 2.418, 95% CI, 1.682-3.474; P<0.001) and surgery treatment (SHR, 0.493, 95% CI, 0.339-0.718; P<0.001) remained significant for the cancer-specific survival in the Fine and Gray's regression model.
CONCLUSIONS: Predictive factors for RLNM of bone sarcoma are sex, tumor site, type and size. In the presence of RKNM, only age, distant metastasis and surgery treatment are prognostic factors for the outcome of patients with bone sarcoma. 2021 Annals of Translational Medicine. All rights reserved.

Introduction

Malignant bone tumors are mainly sarcomas belonging to the mesenchymal tumor family, with osteosarcoma, Ewing sarcoma and chondrosarcoma being the three most common types (1). The survival of patients with bone sarcoma is largely dependent on the status of metastasis, which is generally considered to be achieved via hematogenous routes, with the lungs being the most common metastatic site (2). RLNM is relatively rare in bone sarcoma. Osteosarcoma, the most common type, has a RLNM rate of 1.4–2.3% in several studies (3,4). For Ewing sarcoma and chondrosarcoma, the reported rate was 6.1% and 1.6% for primary tumor on the bone (5-7), while for primary tumor on extra-skeletal sites, the rate of RLNM is higher (31.1% for Ewing sarcoma and 3.1% for chondrosarcoma) (8,9).

Despite the rare incidence, RLNM deserves great attention because it has been identified to be associated with significantly poor outcomes of patients with bone tumor (5,9,10). For osteosarcoma, one animal study has suggested that RLNM is a prognostic factor for unfavorable outcomes (11). In human beings, the reported 5-year overall survival of patients with RLNM is 10.9%, which is significantly lower than that (54.3%) of patients without RLNM (10). For Ewing sarcoma, the second most common bone sarcoma in children, RLNM is more common when the sarcoma originates from extra-skeletal sties, and the 5-year survival is decreased from 60.3% to 45.9% in the presence of RLNM (5). Chondrosarcoma mainly affects old people and is less malignant compared to the above-mentioned two kinds of sarcoma. Over 80% of patients with chondrosarcoma can survive 10 years after treatment (12). However, patients may also suffer local recurrence and death within 2 years after treatment if regional lymph nodes are involved (7). These findings indicate that clinical attention should be paid to the status of RLNM.

Preoperative identification of patients at a high risk of RLNM can help surgeons take related measures to optimize the treatment for these patients. However, it is still difficult to recognize such high-risk patients. Factors associated with increased risks for RLNM of bone sarcoma are still not clear, and prognostic factors for the overall and cancer-specific survival also remain to be identified. Based on these considerations, we conducted this study to deeply explore risk factors for RLNM and performed survival analyses of bone sarcoma with RLNM, which may help clinical surgeons better deal with this disease.

We present the following article in accordance with the STROBE reporting checklist (available at http://dx.doi.org/10.21037/atm-20-4681).

Methods

Study population

This retrospective study was conducted by using data from the Surveillance Epidemiology and End Results (SEER) database, which collects the clinical information of patients diagnosed with various cancers and covers approximately 30% of the total US population. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013), and individual consent for this retrospective analysis was waived as the patient data were extracted from the public SEER database. For our study, a total of 10,641 patients diagnosed with malignant bone sarcoma from 1983 to 2014 were finally enrolled from the database, with 311 having RLNM. Data of all patients were downloaded using the SEER*Stat software (version 8.3.5, National Cancer Institute, Bethesda, MD, USA). The diagnosis of each bone sarcoma was confirmed by positive histological examinations, while those diagnosed by positive exfoliative cytology, direct visualization or radiology were excluded. The final population to be investigated in this study was selected based on the following inclusion criteria: (I) malignant tumors diagnosed between 1983 and 2014, with the bone being the primary site; (II) confirmation of tumor by positive histology from living patients, not from biopsy or other dead cases; (III) known months of survival and complete follow-up information; (IV) known status of regional lymph node metastasis (RLNM); (V) known types of tumor and status of distant metastasis. A total of 10,641 patients meeting these criteria were finally included. Details of inclusion and exclusion processes are demonstrated in .

Figure 1

A flowchart demonstrating the inclusion criteria and exclusion process to the final study cohort.

Variable definition

Tumor type was determined according to the “ICD-O-3 Hist/behav” field in the SEER database, while RLNM information was obtained from the “EOD 4-nodes (1983–1987)”, “EOD 10-nodes (1988–2003)” and “CS-lymph nodes (2004–2015)” fields for different periods of time. The treatment in this study was mainly local management such as radiotherapy and surgery. Surgery treatment was specific to cancer-directed surgeries, including tumor destruction, radical resection, limb salvage and amputation. Those without surgery treatment or undergoing non-cancer-directed procedures such as incisional, needle or aspiration biopsy were considered to be no-surgery-treated cases. Because chemotherapy was documented as “yes” and “no/unknown” in SEER database, we did not know whether patients labeled as “unknown” received chemotherapy or not. In order to avoid possible bias in survival analysis, chemotherapy was excluded in the final multivariate analysis. Variables analyzed in our study included age, gender, race, tumor size, primary site, metastasis and treatment, most of which have been reported to be associated with the prognosis of bone sarcoma (3,13). Age was grouped into <30, 30–60, and >60 years. Tumor size was divided into four groups: <6, 6–12, >12 cm and unknown size. Metastatic information was obtained according to the SEER historic stage A [1973-2015] field. The primary outcome of interest was overall survival and cancer-specific survival. Overall survival was defined as survival span from first diagnosis to death of any cause, while cancer-specific survival was time from first diagnosis to death attributed to the primary bone sarcoma.

Statistical analysis

For baseline characteristics, categorical variates are presented as frequencies and percentages, whereas continuous variates are presented as median and interquartile range (IQR). Analysis of between-group difference in the baseline level was conducted by using the Wilcoxon rank-sum tests for continuous variates and the Chi-square test for categorical variates. In the survival analysis, age and tumor size were grouped as categorical variates. Univariate analysis was performed using the Kaplan-Meier curve and the log-rank test. Taking non-cancer-specific death as a competing risk to cancer-specific death, we plotted cumulative incidence function curves for each bone tumor and performed the Gray’s test to compare the difference in the incidence rate of death (14). Multivariate analysis was performed using the Cox proportional hazards regression model and the Fine and Gray’s regression model (14). Patients who were alive at last follow-up were considered as censored cases in the Cox analysis, while in the Fine and Gray’s regression model, non-cancer-specific death was considered as the competing risk, and sub-distribution hazards ratio (SHR) was used to represent each variable’s contribution to cancer-specific death. Statistical significance was defined as a P value less than 0.05. All statistical analyses were performed in the R software (version 3.6.3). In the software, the “survival” package was used for the Kaplan-Meier curve and the Cox regression analysis, whereas the “timereg” and “cmprsk” packages were used for the competing risk analysis.

Results

Demographics and clinical features of patients with or without RLNM

The proportion of RLNM was 2.9% in the total cohort and was significantly different among different sarcomas (, P<0.001), with Ewing sarcoma having the highest (6.0%) rate, followed by osteosarcoma (2.5%), while chondrosarcoma had the lowest rate of RLNM (1.1%). Other tumors together had a RLNM rate of 4.2%. The total cohort was divided into a RLNM and a non-RLNM group. Baseline characteristics of patients in the two groups are presented in . Demographic characteristics such as age and race were similar between patients with and without RLNM, while male patients were more likely to have RLNM (65% vs. 56.2%, P=0.003). RLNM tended to occur in bone sarcoma with a larger size [median (IQR), 9 cm (5–14 cm) vs. 7 cm (5–11 cm); P<0.001], and the median survival in the RLNM group was significantly shorter [median (IQR), 17 months (7–44 months) vs. 54 months (21–114 months); P<0.001] compared to the non-RLNM group. The two groups also differed in the distribution of primary tumor sites, tumor types and treatments, with all p-values being less than 0.001. Distant metastasis was present in over half (51.1%) of patients in the RLNM group, which was significantly higher than that in the non-RLNM group (15%, P<0.001).

Figure 2

Rate of regional lymph node involvement for different bone sarcomas. Other sarcomas were calculated together due to low incidence. The difference of RLNM rate among different sarcomas was analyzed by chi-square test. RLNM, regional lymph node metastasis.

Table 1

Characteristics of enrolled patients divided into non-RLNM and RLNM groups

Characteristics	Total	Non-RLNM†	RLNM	P
N (%)	10,641	10,330 (97.1%)	311 (2.9%)
Age, median [IQR‡]	38 [17, 58]	38 [17, 58]	34 [16, 59]	0.366
Race (%)				0.503
Black	1,127 (10.6)	1,088 (10.5)	39 (12.5)
White	8,620 (81.0)	8,375 (81.1)	245 (78.8)
Others§	894 (8.4)	867 (8.4)	27 (8.7)
Sex (%)				0.003
Female	4,634 (43.5)	4,525 (43.8)	109 (35.0)
Male	6,007 (56.5)	5,805 (56.2)	202 (65.0)
Primary site (%)				<0.001
Spine and pelvis	2,389 (22.6)	2,305 (22.5)	84 (27.6)
Extremities	5,679 (53.8)	5,556 (54.2)	123 (40.5)
Skull, face and mandible	2,492 (23.6)	2,395 (23.4)	97 (31.9)
Grade (%)				<0.001
Well/moderately differentiated	2,980 (28.0)	2,946 (28.5)	34 (10.9)
Poorly differentiated/undifferentiated	3,907 (36.7)	3,773 (36.5)	134 (43.1)
Unknown	3,754 (35.3)	3,611 (35.0)	143 (46.0)
Tumor type (%)				<0.001
Osteosarcoma	3,675 (34.5)	3,583 (34.7)	92 (29.6)
Chondrosarcoma	3,128 (29.4)	3,095 (30.0)	33 (10.6)
Ewing sarcoma	1,352 (12.7)	1,271 (12.3)	81 (26.0)
Other sarcomas	2,486 (23.4)	2,381 (23.0)	105 (33.8)
Distant metastasis (%)				<0.001
No	8,928 (83.9)	8776 (85.0)	152 (48.9)
Yes	1,713 (16.1)	1554 (15.0)	159 (51.1)
Size, median [IQR]	7 [5, 11]	7 [5, 11]	9 [5, 14]	<0.001
Radiation (%)				<0.001
No	8,961 (84.2)	8,730 (84.5)	231 (74.3)
Yes	1,680 (15.8)	1,600 (15.5)	80 (25.7)
Chemotherapy (%)				<0.001
No	5,545 (52.1)	5,454 (52.8)	91 (29.3)
Yes	5,096 (47.9)	4,876 (47.2)	220 (70.7)
Surgery (%)				<0.001
No	1,793 (16.9)	1,664 (16.2)	129 (41.7)
Yes	8,789 (83.1)	8,609 (83.8)	180 (58.3)
Survival, median (IQR), months	52 (21,112)	54 (21,114)	17 (7,44)	<0.001

†, regional lymph node metastasis. ‡, interquartile range. §, American Indian/Alaska Native, Asian/Pacific Islander. RLNM, regional lymph node metastasis.

Risk factors for RLNM of bone sarcoma

Univariate and multivariate logistic regression analyses were used to identify risk factors for RLNM. The factors with a P value less than 0.1 in the univariate analysis were included in the final multivariate analysis to adjust for potential cofounding factors (). The result indicated that male patients [odds ratio (OR): 1.334, 95% CI, 1.046–1.712; P=0.022], primary tumor site on head and face bones (OR: 1.778, 95% CI, 1.287–2.461; P<0.001), Ewing sarcoma (OR: 2.228, 95% CI, 1.556–3.185; P<0.001), other sarcomas (OR: 1.744, 95% CI, 1.234–2.469; P=0.002) and tumor size >12 cm (OR: 2.240, 95% CI, 1.545–3.249; P<0.001) were associated with a higher risk of RLNM, while chondrosarcoma had a lower risk of RLNM compared to osteosarcoma (OR: 0.550, 95% CI, 0.337–0.879; P=0.014).

Table 2

Logistic regression analysis for risk factors of regional lymph node metastasis

Variables	Univariate analysis			Multivariate analysis
	OR (95% CI)†	P		OR (95% CI)	P
Age, years
<30	Reference			Reference
30–60	0.759 (0.581, 0.988)	0.042		1.089 (0.794, 1.486)	0.594
>60	0.941 (0.704, 1.247)	0.675		1.247 (0.877, 1.764)	0.215
Race
Black	Reference
White	0.816 (0.586, 1.167)	0.247
Others‡	0.869 (0.523, 1.424)	0.58
Sex
Female	Reference			Reference
Male	1.445 (1.143, 1.835)	0.002		1.334 (1.046, 1.712)	0.022
Primary site
Spine and pelvis	Reference			Reference
Extremities	0.607 (0.459, 0.807)	0.001		0.807 (0.592, 1.104)	0.177
Skull, face and mandible	1.111 (0.826, 1.499)	0.487		1.778 (1.287, 2.461)	<0.001
Tumor type
Osteosarcoma	Reference			Reference
Chondrosarcoma	0.415 (0.274, 0.613)	<0.001		0.478 (0.289, 0.775)	0.003
Ewing sarcoma	2.482 (1.826, 3.367)	<0.001		2.228 (1.556, 3.185)	<0.001
Other sarcomas	1.717 (1.293, 2.286)	<0.001		1.744 (1.234, 2.469)	0.002
Distant metastasis
No	Reference			Reference
Yes	5.907 (4.699, 7.429)	<0.001		4.554 (3.569, 5.813)	<0.001
Grade
Well or moderately differentiated	Reference			Reference
Poorly differentiated or undifferentiated	3.077 (2.132, 4.566)	<0.001		1.466 (0.948, 2.317)	0.093
Unknown	3.431 (2.384, 5.081)	<0.001		1.132 (0.727, 1.805)	0.591
Tumor size, cm
<6	Reference			Reference
6–12	1.217 (0.895, 1.658)	0.211		1.230 (0.879, 1.726)	0.228
>12	2.395 (1.736, 3.303)	<0.001		2.240 (1.545, 3.249)	<0.001
Unknown	1.227 (0.880, 1.707)	0.225		1.070 (0.748, 1.528)	0.709

†, OR, odds ratio; CI, confidence interval. ‡, American Indian/Alaska Native, Asian/Pacific Islander.

Survival and outcome for bone sarcoma patients with RLNM

Of the total cohort, the 311 patients with RLNM were subsequently included for survival analysis. The 5-year overall survival rate for osteosarcoma, chondrosarcoma and Ewing sarcoma with RLNM was 25.6% (95% CI, 16.9–38.7%), 30.0% (95% CI, 16.9–53.2%) and 45.3% (95% CI, 34.5–59.4%), respectively. Log-rank tests suggested that Ewing sarcoma was associated with a better overall survival compared to osteosarcoma (P<0.001) or chondrosarcoma (P=0.003), while no difference was observed between chondrosarcoma and osteosarcoma (OS: P=0.9) (). We also plotted cumulative incidence function curves considering non-cancer-specific death as the competing risk to the cancer-specific death. As shown in , the cumulative incidence of cancer-specific death in Ewing sarcoma was lower than that in osteosarcoma (Gray’s test, P=0.033) but not in chondrosarcoma (Gray’s test, P=0.181). Unexpectedly, the incidence of the competing risk of Ewing sarcoma was also lower than that of osteosarcoma (Gray’s test, P=0.004) and chondrosarcoma (Gray’s test, P=0.017), suggesting that the difference of the overall survival between Ewing sarcoma and the other two sarcomas partially came from non-cancer-specific death.

Figure 3

Kaplan-Meier curves and cumulative incidence function curves of bone sarcomas with regional lymph node involvement. Overall survival between osteosarcoma and Ewing sarcoma (A), chondrosarcoma and Ewing sarcoma (B), osteosarcoma and chondrosarcoma (C). The cumulative incidence of cancer-specific death and competing risks between osteosarcoma and Ewing sarcoma (D), chondrosarcoma and Ewing sarcoma (E), osteosarcoma and chondrosarcoma (F).

To adjust for potential confounding factors, multivariate Cox regressions for overall and cancer-specific survival were conducted by including all possible prognostic factors. The results suggested that age, race, distant metastasis and surgery were independent prognostic factors for both the overall and the cancer-specific survival (). Compared to Ewing sarcoma, chondrosarcoma was associated with a decreased overall survival (HR, 1.540, 95% CI, 1.011–2.345; P=0.044) but not cancer-specific survival (HR, 1.377, 95% CI, 0.877–2.163; P=0.165), and similar results were observed for other sarcomas. For osteosarcoma, there was no difference in the overall (HR, 1.627, 95% CI, 0.914–2.896; P=0.098) and cancer-specific survival (HR, 1.664, 95% CI, 0.884–3.135; P=0.115) compared to Ewing sarcoma.

Table 3

Multivariate Cox analysis of overall and cancer-specific survival for 311 patients with RLNM

Variables	Overall survival			Cancer-specific survival
	HR (95% CI)†	P		HR (95% CI)	P
Age, years
<30	Reference			Reference
30–60	2.067 (1.439, 2.971)	<0.001		1.881 (1.259, 2.809)	0.002
>60	2.133 (1.462, 3.112)	<0.001		1.474 (0.948, 2.293)	0.085
Sex
Female	Reference			Reference
Male	1.277 (0.953, 1.713)	0.102		1.255 (0.896, 1.757)	0.187
Race
Black	Reference			Reference
White	0.622 (0.410, 0.946)	0.026		0.654 (0.404, 1.060)	0.085
Others	0.330 (0.171, 0.636)	0.001		0.337 (0.160, 0.707)	0.004
Grade
Low grade‡	Reference			Reference
High grade§	1.522 (0.952, 2.432)	0.079		1.705 (0.975, 2.981)	0.061
Unknown	0.637 (0.395, 1.027)	0.064		0.620 (0.350, 1.101)	0.103
Primary site
Spine and pelvis	Reference			Reference
Extremities	1.244 (0.864, 1.791)	0.24		1.351 (0.902, 2.025)	0.145
Skull, face and mandible	1.262 (0.851, 1.872)	0.247		1.160 (0.734, 1.832)	0.525
Distant metastasis
No	Reference			Reference
Yes	2.018 (1.477, 2.756)	<0.001		2.448 (1.702, 3.522)	<0.001
Tumor size
<6 cm	Reference			Reference
6–12 cm	1.207 (0.796, 1.831)	0.376		1.121 (0.693, 1.815)	0.641
>12 cm	1.160 (0.755, 1.783)	0.498		1.143 (0.704, 1.855)	0.589
Unknown	1.228 (0.796, 1.892)	0.353		1.010 (0.603, 1.693)	0.969
Tumor type
Ewing sarcoma	Reference			Reference
Chondrosarcoma	1.540 (1.011, 2.345)	0.044		1.377 (0.877, 2.163)	0.165
Osteosarcoma	1.627 (0.914, 2.896)	0.098		1.664 (0.884, 3.135)	0.115
Other sarcomas	1.832 (1.187, 2.828)	0.006		1.544 (0.961, 2.481)	0.072
Radiotherapy
No	Reference			Reference
Yes	0.876 (0.604, 1.270)	0.485		1.028 (0.674, 1.567)	0.9
Surgery
No	Reference			Reference
Yes	0.443 (0.307, 0.638)	<0.001		0.383 (0.253, 0.581)	<0.001

†, HR, hazard ratio; CI, confidence interval. ‡, well or moderately differentiated. §, poorly differentiated or undifferentiated. RLNM, regional lymph node metastasis.

In the traditional Cox proportional hazards regression analysis for cancer-specific survival, patients dying of other causes are usually censured, leading to discrepancies because death of other causes precludes death of the primary tumor (15). Considering death not caused by specific bone sarcoma as the competing risk, we further conducted the Fine and Gray’s regression analysis (). Significant factors associated with cancer-specific death after controlling for competing risks were only age between 30–60 years (SHR, 1.528, 95% CI, 1.028–2.271; P=0.02), distant metastasis (SHR, 2.418, 95% CI, 1.682–3.474; P<0.001) and surgery treatment (SHR, 0.493, 95% CI, 0.339–0.718; P<0.001). There was no difference in the cancer-specific death among different bone sarcomas with RLNM (chondrosarcoma: SHR, 1.603, 95% CI, 0.866–2.969, P=0.13; osteosarcoma: SHR, 1.260, 95% CI, 0.816–1.945, P=0.13; other sarcomas: SHR, 1.382, 95% CI, 0.875–2.184, P=0.17). Interestingly, tumor size was not an independent prognostic factor in all multivariate analyses.

Table 4

Multiple competing risk regression analysis of patients with RLNM

Variables	Sub-distribution hazard ratio (95% CI†)	P
Age, years
<30	Reference
30–60	1.528 (1.028, 2.271)	0.02
>60	0.958 (0.610, 1.505)	0.91
Sex
Female	Reference
Male	1.095 (0.778, 1.542)	0.67
Race
Black	Reference
White	0.953 (0.544, 1.668)	0.87
Others	0.511 (0.232, 1.125)	0.08
Grade
Low grade‡	Reference
High grade§	1.257 (0.763, 2.073)	0.32
Unknown	0.597 (0.359, 0.994)	0.47
Primary site
Spine and pelvis	Reference
Extremities	1.445 (0.976, 2.141)	0.07
Skull, face and mandible	1.086 (0.670, 1.762)	0.96
Distant metastasis
No	Reference
Yes	2.418 (1.682, 3.474)	<0.001
Tumor size, cm
<6	Reference
6–12	1.225 (0.732, 2.048)	0.44
>12	1.246 (0.757, 2.050)	0.39
Unknown		0.96
Tumor type
Ewing sarcoma	Reference
Chondrosarcoma	1.603 (0.866, 2.969)	0.13
Osteosarcoma	1.260 (0.816, 1.945)	0.3
Other sarcomas	1.382 (0.875, 2.184)	0.17
Radiotherapy
No	Reference
Yes	1.113 (0.746, 1.660)	0.6
Surgery
No	Reference
Yes	0.493 (0.339, 0.718)	<0.001

†, CI: confidence interval. ‡, well or moderately differentiated. §, poorly differentiated or undifferentiated. RLNM, regional lymph node metastasis.

To further confirm the credibility of the competing risk regression model, we conducted a series of sensitivity analyses, which included the following analyses: (I) analysis excluding patients without data of tumor size; (II) analysis excluding patients without tumor grade data; (III) analysis excluding patients with a survival less than 3 or 5 months to preclude the bias from surgery or other-treatment-related deaths; (IV) analysis including patients aged less than 30 years or older than 50 years to separately analyze young and old patients (Tables S1-S3). The results remained stable among these sensitivity analyses, except for that the survival of chondrosarcoma was poorer than Ewing sarcoma in patients younger than 30 years (SHR, 3.373, 95% CI, 1.05–10.77; P=0.04). Comprehensively, these results indicated that RLNM poses a similar cancer-specific outcome to patients with different bone sarcomas, and that patients with RLNM have a higher risk of death in the presence of distant metastasis. However, they can benefit from primary-cancer-directed surgeries.

Discussion

Bone sarcomas are generally considered to metastasize via the hematogenous routes, mostly to the lungs and the bone (2,16), while involvement of the lymph system is rare. Despite the rare incidence, the outcome of patients with RLNM is poorer than that in those without. Identifying patients at a high risk of RLNM is of great clinical significance, because it helps surgeons improve the treatment for such patients. However, to the best of our knowledge, there is no study available to explore risk factors for RLNM of bone sarcoma with relatively large size of sample. In addition, prognostic factors for patients with bone sarcoma with RLNM are also unclear. Considering this, we mined data from the SEER database and conducted this study, which is the first to provide a deep understanding of risk factors for RLNM and the survival outcomes of patients with bone sarcoma with RLNM.

The proportion of RLNM was significantly different among different bone sarcomas. Ewing sarcoma was found to have a higher rate of RLNM than osteosarcoma and other sarcomas, and chondrosarcoma was the least likely to have RLNM. This is consistent with previously reported rates of RLNM of different bone sarcomas (1.4–2.3% for osteosarcoma, 1.6% for chondrosarcoma and 6.1% for Ewing sarcoma) from separate single centers (3,4,6,7). The rate of RLNM for each sarcoma in our study is also similar to previous studies using data from the same database but including patients confirmed of bone sarcoma in different year spans (5,9,10), which further adds to the credibility that the rate in our study represents the actual rate of RLNM for different bone sarcomas.

Using the multivariate logistic regression model, we found tumor type was a risk factor for RLNM. Ewing sarcoma and other bone tumors were more likely to metastasize to regional lymph nodes, while chondrosarcoma was associated with a lower risk of RLNM. This is in agreement with the different tendencies of distant metastasis for different bone sarcomas. The rate of distant metastasis at first diagnosis has been reported to be about 20% for Ewing sarcoma, followed by 15% for osteosarcoma (17,18). For chondrosarcoma, this rate decreases to less than 10% (19), possibly due to the relatively low grade of malignancy and the less likelihood of metastasis of this sarcoma (20,21). In the multivariate logistic analysis, distant metastasis was also suggestive of an increased risk for RLNM. These two types of metastasis might together reflect the divergent traits of metastasis of different bone sarcomas. Another important factor predictive of RLNM was larger tumor size (>12 cm), which has been reported to be a risk factor for RLNM of osteosarcoma (22). This reminds clinical surgeons that more attention should be paid to RLNM when resecting tumor with a large size.

Although Ewing sarcoma was more likely to have RLNM in this study, patients with such tumor had a better overall survival compared to other tumors with RLNM. However, such difference was eliminated in the cancer-specific survival regarding non-cancer-specific deaths as the competing risk, which suggests a comparable lethality of different tumors in the presence of RLNM. This further arouses concerns about the role of non-cancer-specific death in patient’s survival. These deaths, including those attributed to cardiovascular events, diabetes or suicides, have been reported to be associated with the outcome of bone sarcoma or other cancers (23-25). Older age (>60 years) was associated with a poor overall survival, but not cancer-specific survival for bone sarcoma patients with RLNM, which might also be attributed to the fact that old patients are at a higher risks of non-cancer-caused death (26).

Old-age-related non-cancer-specific death may also explain the different overall survival but not cancer-specific survival among different sarcomas, as osteosarcoma and chondrosarcoma have a relatively higher proportion of old patients compared to Ewing sarcoma, which affects mainly young patients (20). Actually, using the Gray’s test, we found the cumulative incidence of non-cancer-specific death in Ewing sarcoma was significantly lower than that in osteosarcoma and chondrosarcoma. This difference is probably the cause of the different overall survival for patients with RLNM of different sarcomas. This reminds us that in addition to the treatment of tumor, preventing non-cancer-specific death in old patients is also a favorable strategy to improve the survival of patients with bone sarcoma.

The result of the multivariate competing risk regression model suggested that the cumulative incidence of cancer-specific death was not different in patients with different kinds of bone sarcoma. Chondrosarcoma is a less malignant tumor compared to other bone sarcomas. However, according to this study, RLNM confers equally poor outcomes to patients with chondrosarcoma compared to those with osteosarcoma and Ewing sarcoma. Tumor size, a risk factor for RLNM, was also not associated with the outcome of patients in the presence of RLNM. Small tumor is generally considered easy to resect and these patients are more likely to have a good outcome. However, our study shows that small tumor poses a comparable poor outcome to large tumor in the presence of RLNM. Cancer-directed surgery is associated with a decreased risk of cancer-specific death in the multivariate risk regression model, indicating that patients can benefit from surgery treatment. Thus, for patients at a high risk of RLNM, surgical treatment is recommended. In all, these findings remind doctors that the status of RLNM should be considered when resecting bone tumors, even for those with a small size or with a low-grade malignancy, to avoid potential poor outcomes associated with RLNM.

There were several limitations in this study. First, as lymph node examination is not a routine practice in bone tumor surgeries, the status of lymph node, the most important information for this study, was missing for a portion of patients. These patients were excluded from this study, causing difficulties in calculating the actual rate of RLNM. However, the rate of RLNM for different bone sarcomas analyzed in our study is similar to that of available single-center studies (3,4). This convinced us that the rate would not change much when those excluded patients were included. Second, information of some key factors associated with the survival of bone sarcoma, such as tumor size and grade (13,27), was not available for all patients included in the final analysis, which caused potential bias to the result. Considering this, we conducted related sensitivity analyses by excluding those patients missing information of tumor size or grade, and the results were in consistence with our main analysis. Third, although not effective for chondrosarcoma (28), chemotherapy is important for the treatment of osteosarcoma and Ewing sarcoma (29,30). However, in the SEER database, chemotherapy information is unclearly documented, with those without chemotherapy and with unknown information of chemotherapy being labeled together as “no/unknown” in one group. We were unable to tell the patients truly without chemotherapy from other patients in the group, and chemotherapy was not included in our analysis. Consequently, we could not determine the role of chemotherapy on patients with bone sarcoma in the presence of RLNM. Further studies are needed to address the role of chemotherapy for such patients.

In conclusion, our study is the first to provide a deep understanding, including risk factors and survival analysis, for bone sarcoma with RLNM. As RLNM is associated with a significantly poor outcome for all bone sarcomas, more attention should be paid by clinical surgeons to the status of RLNM, especially for Ewing sarcoma, which has the highest rate of RLNM. Although chondrosarcoma is less malignant than other bone sarcomas, we should be cautious about the status of RLNM, which confers an equally poor survival compared to other bone sarcomas. It should be admitted that the rate of RLNM in our and other studies is based on non-routine clinical practices. Further studies are needed to explore the actual rate of RLNM for different sarcomas with intended regional lymph node examination. Considering that RLNM poses a significantly poorer outcome to patients with bone sarcoma, preoperative detection of RLNM and improvement of survival for these patients by surgical, radiological or chemotherapeutic strategies also deserve more research. Finally, non-cancer-specific death, which contributes to different overall-survival outcomes among different bone sarcomas with RLNM, may be focused on to improve the outcome of cancer patients.

The article’s supplementary files as