Abstract
Purpose. Distal metastasis in chondrosarcoma is associated with a poor prognosis. The aim of this study was to develop and validate columnar maps to predict the risk of distal metastasis in patients with chondrosarcoma, thereby contributing to clinical diagnosis and treatment. Methods. Data from chondrosarcoma patients obtained from the Surveillance, Epidemiology, and End Results (SEER) database from 2010 to 2016 were then screened by univariate and multifactorial logistic regressions to construct a predictive distal metastasis risk. The model discrimination of nomogram was assessed by calibration plots, while the predictive accuracy and clinical values were measured by decision curve analysis (DCA). In addition, predictive column line plots were validated in an in-house test set. Results. A total of 1,290 patients were included and randomized in a 7 to 3 ratio into a training group (n = 906) and a test group (n = 384). After logistic regression analysis, the significant variables were gender, tumor pathological grade, laterality, primary tumor stage, regional lymph node metastasis, surgical treatment, and chemotherapy. Calibration curves showed agreement between column line graph predictions and actual observations, while DCA showed the clinical utility of the nomogram. In addition, ROC showed good discrimination and calibration in the training (AUC = 0.937, 95% CI 0.919–0.952) and validation groups (AUC = 0.91, 95% CI 0.877–0.937). Conclusions. The nomogram for distal metastasis risk in patients with chondrosarcoma can effectively predict the individualized risk of distal metastasis and provide clinicians with enlightening information to optimize treatment options.
1. Introduction
Chondrosarcomas (CSs) are a group of heterogeneous bone malignancy with diverse histopathological and clinical features characterized by the production of a cartilaginous stroma [1], which are the second most frequent type of bone malignancy after osteosarcoma, accounting for approximately 20% of all types of bone malignancy. They are commonly found in male adults within flat bones; the pelvis and femur are two common sites of involvement, although any bone may be affected. Most chondrosarcomas exhibit indolent, with approximately 90% below intermediate grade, and stably behave and rarely metastasize [1]. Approximately 8% of the patients with chondrosarcoma developed distant metastasis [2]. Recently, there are no common therapies such as radiotherapy or chemotherapy regimens, and targeted drugs are still under basic research and clinical trials [3]. Fortunately, CS is insensitive to chemotherapy and radiotherapy, after complete surgical resection at present. The prognosis of combined pre- and postsurgical radiotherapy and chemotherapy is significantly better than that of surgical treatment alone [3]. Some studies have shown that distant metastases occur in 8–38% of patients with chondrosarcoma, which can greatly affect the execution of surgery and make complete resection of the tumor extremely difficult [4]. Pulmonary metastases with local recurrence are the most common reason for death in CS. Prospective research treatment of chondrosarcoma is, therefore, necessary [5], and clinicians treating patients with chondrosarcoma must determine the likelihood of metastasis [6], and it is necessary to identify risk factors for distal metastases, requiring an expansion of treatment options and approaches to improve clinical outcomes. For the current studies, predicting survival in individual patients remains difficult [7].
Since chondrosarcomas with lower incidence are relatively rare, studies assessing prognostic factors are difficult to carry out while requiring a large sample size. The National Cancer Institute’s SEER program is a comprehensive source of population-based cancer incidence and survival information from the United States, collecting 18 population-based cancer registries representing approximately 27.8% of the total U.S. population [8]. This study is based on the SEER database, which does not require patient authorization, to investigate risk factors for distal metastasis of chondrosarcoma at the time of initial diagnosis, considering that this database has been providing site-specific data on metastatic tumors since 2010, by collecting information on the demographic and clinical characteristics of chondrosarcoma.
Nomogram is commonly used to generate the likelihood of clinical events through complex computational formulas [4, 9]. With nomograms, clinicians can assess the risk of clinical events, personalize individual treatment alternatives, and optimize treatment strategy. Considering the important role of distal metastasis in the chondrosarcoma prognosis, this research study aimed to assess patients at high risk of distal metastasis from chondrosarcoma through the use of a nomogram.
2. Methods
2.1. Data Sources and Inclusion Criteria
This study was based on the clinical data of patients with chondrosarcoma from the SEER database. Patients with a pathological diagnosis of chondrosarcoma in the SEER database were retrieved through the SEERStat software, along with the third edition of the International Classification of Oncology (ICDO-3), morphology code (9220) used to identify chondrosarcoma. The data in this study included patients diagnosed with chondrosarcoma from 2010 to 2016. Exclusion criteria were as follows: (1) patients with no positive pathology; (2) patients with unknown survival; (3) tumors that were not the first occurrence; (4) more than one primary tumor; (5) distal metastasis information unknown; and (6) regional lymphatic fluid information was incomplete for lymph node metastases.
Data were extracted from the SEER database including age, sex, race, primary site, survival duration, laterality, tumor pathology grade, primary tumor stage, surgical treatment, radiotherapy, chemotherapy, and lymph node metastasis. These data were derived from the variable “CS site-specific factor 6.” In addition, less than 20 tumor sites were categorized as “other.”
2.2. Nomogram Construction, Validation, and Clinical Application
Patients with chondrosarcoma who met the inclusion criteria were randomly divided into training and validation groups in a ratio of 7 to 3. Subsequently, the following variables were selected for the study: age, race, gender, laterality, survival time, tumor pathological grade, primary tumor metastasis, primary site, surgical treatment, radiotherapy, chemotherapy, and lymph node metastasis. Univariate and multivariate binary logistic regressions were applied to identify independent risk factors with a forward stepwise regression method. Nomogram was constructed based on the results of logistic regression analysis. Calibration plot chart of clinical prediction model (calibration plot) and ROC curve were plotted, and ROC was used to estimate the predictive performance of column line plot. The larger area under the ROC curve (AUC) means the better the discriminatory ability or prognostic accuracy of the variable. In addition, decision curve analysis (DCA) plots the net benefit (NB) under a range of reasonable risk thresholds that were consistent with clinical practice and were used to assess the clinical utility of column line plots in decision-making.
2.3. Statistical Methods and Software
Continuous variables were expressed as mean ± standard deviation (SD), and categorical variables were expressed as proportions. Continuous and categorical variables were compared by t-test and chi-square test of SPSS, respectively. IBM SPSS Statistics version 26.0 (SPSS Inc., Chicago, IL, USA) and R software version 4.0.5 (http://www.r-project.org) performed the above statistical methods, and several R packages (including regplot, rms, rmda, and pROC) were applied to plot graphs, such as nomogram, calibration plot, DCA plot, and ROC curve. KM curves were plotted by GraphPad Prim 8.0. All values were bivariate, values of were considered statistically significant, and confidence intervals (CIs) were expressed at the 95% confidence level.
3. Results
3.1. Results of Single-Factor and Multifactor Logistic Regression
A total of 1,290 patients were included in the statistics, and by the univariate and multifactorial logistic regressions, the extracted variables were first subjected to univariate logistic regression analysis, which showed that age, survival time, tumor pathology grading classification, gender, laterality, primary tumor stage, surgical treatment, and chemotherapy were prognostic factors affecting distal metastasis (). Further multifactorial logistic risk was performed, resulting in seven factors as independent prognostic factors for distal metastasis (Table 1), such as gender (female: dominance ratio (OR) 0.368, 95% CI (0.150–0.901), ), tumor pathological grade (GIII: OR = 6.921, 1.410–33.983, ; unclear tumor grading: OR = 6.216, 1.385–27.891, ), laterality (right: OR = 2.674, 1.029–6.945, ), primary tumor stage (T2 : OR = 4.698, 1.732–12.738, ; T3 : OR = 59.117, 8.152–428.709, ; TX : OR = 4.074, 1.221–13.594, ), regional lymph node metastasis (positive: OR = 9.168, 1.572–53.669, ; unknown: 6.743, 1.874–24.269, ), surgical treatment (OR = 0.237, 0.101–0.555, ), and chemotherapy (OR = 19.188, 7.554–48.740, ). The specific results are shown in Table 1 as single- and multifactor logistic regression analysis.
3.2. Baseline Characteristics of Patients with Nonmetastatic and Metastatic Chondrosarcoma
The 1,290 patients were grouped by M0 and M1 and are then summarized in Table 2. Age, tumor pathological grade, laterality, primary tumor stage, regional lymph node metastasis, surgical treatment, chemotherapy, and survival time were shown to be prognostic factors affecting distal metastasis ().
3.3. Population Baseline Characteristics
The 1,290 patients were randomized into a training group (n = 906, 7 to 3 ratio) and a validation group (n = 384), which are then summarized in Table 3. There were no statistically significant differences between the training and validation groups ().
3.4. Construction and Validation of Nomogram for Distal Metastasis of Chondrosarcoma
The results of univariate and multivariate logistic regressions were used to construct the nomogram of distal metastasis (Figure 1(a)). The primary tumor stage was the best predictor, followed by chemotherapy, tumor pathological grade, lymph node metastasis, laterality, surgical treatment, and gender. The calibration plots of the nomogram (Figures 1(b) and 1(c)) showed that the apparent curves, in both training and validation groups, showed good consistency. The AUC values of nomogram were 0.937 (95% CI 0.919–0.952) and 0.91 (95% CI 0.877–0.937) in the training and validation groups, respectively (Figures 2(a) and 2(b)). According to the ROC curves in the training set, the values of nomogram were more important than other variables, including tumor pathological grade (AUC = 0.733, 95% CI 0.703 to 0.762), laterality (0.591, 0.558 to 0.623), primary tumor stage (AUC = 0.608, 95% CI 0.575 to 0.640), gender (AUC = 0.635, 95% CI 0.603 to 0.666), and surgical treatment (0.718, 95% CI 0.688 to 0.748). Similarly, nomogram values were higher in the test set than in the univariate, as shown in Table 4.

(a)

(b)

(c)

(a)

(b)
3.5. Clinical Usefulness of the Distal Metastasis Nomogram
Kaplan-Meier survival curves for overall survival (OS) were plotted for the 1,290 patients enrolled in the study (Figure 3(a)), and there was a significant difference between the Kaplan-Meier survival curves for the two groups (), suggesting that patients with chondrosarcoma who are expected to develop distal metastases would have a significant survival disadvantage. Subsequently, the DCA curve shown (Figure 3(b)) with a threshold of approximately 0–0.7 was the maximum gain for distal metastasis.

(a)

(b)

(c)
4. Discussion
This study is a study to analyze the risk of distal metastasis in chondrosarcoma based on the SEER database. According to the study, approximately 8% of patients with chondrosarcoma develop distal metastases [10]. Since patients with chondrosarcoma combined with distal metastases have a poor prognosis, it is necessary to identify relevant factors to identify patients with chondrosarcoma at high risk for distal metastases [11]. The results of this study showed that patients with tumor pathological grade 3, i.e., hypofractionated tumor and tumor with an unclear grade, laterality to the right, higher stage of primary tumor, lymph node metastasis, male, nonsurgical treatment, and chemotherapy alone were at higher risk of distal metastasis.
This study found that tumors with grade 3 pathology versus unclear grade had a higher risk of distal metastasis. In both univariate and multifactorial logistic regression analyses, tumor pathology grade 3 and unclear grade were associated with the risk of distal metastasis, with OR values of 6.921 and 6.216, respectively, representing a 6.921- and 6.216-fold higher risks of distal metastasis in hypofractionated tumors and tumors with unclear grade than in highly differentiated tumors. It has also been shown that the main prognostic factor for patients with chondrosarcoma is the grade of the tumor, with increased pathological grade suggesting a poor prognosis [12]. Approximately 85% of these tumors are low grade, and overall survival is favorable [13], due to the fact that low-grade chondrosarcomas have abundant cartilage stroma, low cell density, are easily locally confined, and have a good prognosis after surgical resection, whereas high-grade tumors have little cartilage stroma, high cell density, and are prone to metastasis, leading to a poor prognosis [1]. Therefore, high-grade tumors have been considered as an independent risk factor for metastasis and death [14], which is also consistent with the findings of this study. Regarding laterality at diagnosis, this study concluded that the more right-sided tendency at diagnosis, the greater the risk of developing metastasis. The results of logistics analysis in Table 1 show that with right-sided tendency, the risk of developing metastasis is 2.674 times greater than with left-sided tendency. Chondrosarcoma is often a lateral growth, and there are fewer related studies. This study found that right-sided laterality has a higher risk of developing distal metastasis, and the reasons for this need further study. It was found that the higher the primary tumor stage, the higher the risk of distal metastasis. In the univariate and multifactorial logistic regression analyses, T2, T3, and TX were associated with the risk of distal metastasis, with OR values of 4.698, 59.117, and 4.074, respectively, representing that the risk of developing distal metastasis will be higher when the primary tumor is limited to the bone cortex, exceeds the bone cortex, or cannot be determined 4.698, 59.117, and 4.074 times, with the highest risk of distal metastasis at T3 stage. Possible reasons for this phenomenon are that as the tumor volume increases, the depth of infiltration and the extent of collateral tissue involvement increase, and the tumor becomes more aggressive and also predicts a higher degree of malignancy, thus increasing the likelihood of metastasis [15]. It has been suggested that larger chondrosarcomas may be a predictor of poor survival expectations [16], which is also consistent with the findings of this study.
It is also noteworthy that this study found studies showing that inability to undergo surgical resection and chemotherapy alone has a higher risk of metastasis, with an OR of 19.188 in univariate and multifactorial logistic regression analyses, representing an inability to resect and a 19.188-fold higher risk of distal metastasis after chemotherapy alone. A related study found chemotherapy to be an important risk factor, and the results indicated that patients receiving chemotherapy were more likely to have higher tumor grade, larger tumor size, and greater tumor extent [17], consistent with the findings of this study. This may be related to the ineffective delivery of chemotherapeutic agents, with some studies showing that tumors poorly respond to chemotherapy, and no significant improvement in disease-free survival or overall survival was seen compared to cohorts that did not receive chemotherapy [18]. The role of chemotherapy in the treatment of patients with localized and advanced chondrosarcoma is unclear. Although its use in conventional chondrosarcoma has been largely ineffective, recent data suggest that it may play a role in certain chondrosarcoma subtypes, particularly dedifferentiated and mesenchymal variants, and reports suggest that chemotherapy may provide benefit for this particular subtype [19]. Therefore, further research is needed on the role of chemotherapy in the treatment of chondrosarcoma. In this study, we found that the OR of the population treated with surgery compared to the population not treated with surgery showed that the risk of distal metastases was only 0.2 times higher in patients who underwent surgery. Surgical resection is the primary treatment for both primary and metastatic chondrosarcoma [20]. The goal of resection is to remove the primary tumor with clear margins to limit recurrence and metastasis. Currently, 90% to 95% of patients with osteosarcoma of the extremities successfully avoid amputation through limb-preserving surgery [21]. It has also been demonstrated in relevant studies that high-grade chondrosarcoma (grade 2 and above) is best treated with extensive resection [22]. Approximately half of the patients have a good or very good prognosis [23], which is consistent with the findings of this study. This study showed that male patients with chondrosarcoma had a higher and statistically different risk of developing distal metastases compared to female patients. OR values showed that female patients had only 0.4 times the risk of developing LM compared to male patients. A related study found that men were an independent risk factor for survival in patients with chondrosarcoma [24]. Considering that men have a higher risk of distal metastasis, it may, therefore, affect survival expectations.
According to the results of the regression analysis in Table 1 logistics, patients with lymphatic metastases had an approximately 9-fold higher risk of distal metastasis than patients without lymphatic metastases (OR = 9.168), and patients with unknown lymphatic metastases had an approximately 7-fold higher risk of distal metastasis than patients without lymphatic metastases (OR = 6.743). The prevalence of regional lymph node involvement in patients with chondrosarcoma is 1.3% due to the lack of lymphatic vessels in the bone, which rarely metastasize through lymph nodes. Therefore, lymph node spread in chondrosarcoma is extremely rare, and lymph node metastasis, although extremely rare in primary osteochondrosarcoma, has been shown to have an overall 5-year survival rate of 28% and 77% in patients with and without regional lymph node metastasis, respectively. Patients with chondrosarcoma with lymph node metastases had a worse prognosis than those who did not report regional lymph node metastases [25, 26], consistent with the findings of this study. This study also found that patients with regional lymph node metastases had different tumor characteristics compared to those without regional lymph node metastases. Most importantly, the finding of regional lymph node involvement independently indicates a lower survival rate in patients with chondrosarcoma; this may be important when planning treatment or advising patients on their condition, so clinicians should more carefully examine patients with chondrosarcoma who have lymph node metastases [27]. Therefore, factors related to lymphatic metastasis and distal metastasis still need to continue to be studied in depth.
Nomogram is a quantitative mathematic tool to assess risk and benefit and has been widely used in the medical field for clinical decision-making in a variety of diseases [28, 29]. In previous studies, several nomograms have been developed and validated to predict specific survival and overall survival in chondrosarcoma [30]. However, a nomogram for predicting distal metastasis has not been reported. In this study, 1,290 chondrosarcoma cases were obtained from the SEER database, and patient prognostic factors (i.e., gender, tumor pathology grading classification, laterality, primary tumor stage, regional lymph node metastasis, surgical treatment, and chemotherapy) based on seven of the logistics regression analyses were used to establish a nomogram for predicting the distal nomogram for distal metastasis that showed better diagnostic efficiency compared to other individual variables as evidenced by calibration plots and ROC curves (Figures 1 and 2, Table 4). This all proves the value of the use of nomograms in this study, which can be further applied and improved in clinical work, and clinicians can choose better medical tests and optimize treatment plans with the help of nomograms.
This study still has limitations. First, this study is a retrospective analysis, and the data are biased and indeed lack systematic and prospective data. Also, as a single-center study, even though it was divided into training and validation groups, it still lacks external validation from other institutions, which may lead to overfitting of the nomogram for predicting distal metastasis.
5. Conclusions
A large population-based cohort from the SEER dataset was screened for this study and statistically analyzed to conclude that age, gender, survival time, tumor pathological grade, laterality, primary tumor stage, and whether or not surgical treatment and chemotherapy were prognostic factors affecting distal metastasis; higher or unclear tumor pathological grade, laterality to the right, higher primary tumor stage, male, lymph node metastasis, chemotherapy, and not the higher tumor pathology grade or rightward, higher primary tumor stage, male, lymph node metastasis, chemotherapy, and no surgical treatment were independent risk factors for distal metastasis. A nomogram was further constructed based on the results of statistical analysis to predict distal metastasis in patients with chondrosarcoma. Based on the results of internal validation, DCA curves and clinical impact maps, the nomogram in this study can effectively predict the individualized risk of distal metastasis.
Data Availability
SEER database within the article is a public dataset.
Ethical Approval
The SEER database is a comprehensive data source developed based on population data and annually updated since its launch in 1973. It is public and identifiably accessible that data analysis is treated as nonhuman subjects by the Office for Human Research Protections. As such, no institutional review board approval and informed consent were required.
Conflicts of Interest
All authors declare that they have no conflicts of interest in this paper.
Authors’ Contributions
Wenle Li, Rong Li, and Wanying Li have equally contributed to this work. CLY, RL, QL, and HWP jointly carried out the entire research design. WLL, RL, WYL, CX, and MMM participated in the research and collected and analyzed data. WLL, RL, WYL, and BW drafted the manuscript. CLY and QL provided expert consultation and advice. All authors conceived this research, participated in its design and coordination, and helped polish the language. All authors reviewed the final version of the manuscript.