Research Article | Open Access
Xinjie Wu, Yanlei Wang, Wei Sun, Mingsheng Tan, "Prognostic Factors and a Nomogram Predicting Overall Survival in Patients with Limb Chondrosarcomas: A Population-Based Study", BioMed Research International, vol. 2021, Article ID 4510423, 17 pages, 2021. https://doi.org/10.1155/2021/4510423
Prognostic Factors and a Nomogram Predicting Overall Survival in Patients with Limb Chondrosarcomas: A Population-Based Study
Introduction. We aimed to develop and validate a nomogram for predicting the overall survival of patients with limb chondrosarcomas. Methods. The Surveillance, Epidemiology, and End Results (SEER) program database was used to identify patients diagnosed with chondrosarcomas, from which data was extracted from 18 registries in the United States between 1973 and 2016. A total of 813 patients were selected from the database. Univariate and multivariate analyses were performed using Cox proportional hazards regression models on the training group to identify independent prognostic factors and construct a nomogram to predict the 3- and 5-year survival probability of patients with limb chondrosarcomas. The predictive values were compared using concordance indexes (-indexes) and calibration plots. Results. All 813 patients were randomly divided into a training group () and a validation group (). After univariate and multivariate Cox regression, a nomogram was constructed based on a new model containing the predictive variables of age, site, grade, tumor size, histology, stage, and use of surgery, radiotherapy, or chemotherapy. The prediction model provided excellent -indexes (0.86 and 0.77 in the training and validation groups, respectively). The good discrimination and calibration of the nomograms were demonstrated for both the training and validation groups. Conclusions. The nomograms precisely and individually predict the overall survival of patients with limb chondrosarcomas and could assist personalized prognostic evaluation and individualized clinical decision-making.
Chondrosarcoma is the most common bone sarcoma in aging populations and accounts for 30% of all malignant bone tumors [1, 2]. Central chondrosarcoma is the most common type and accounts for 70% of all chondrosarcomas that involve the proximal femur and proximal humerus . It is generally thought that, because of their extracellular matrix, low percentage of dividing cells, and poor vascularity, chondrosarcomas are chemo- and radiotherapy-resistant. Surgical resection is the most widely accepted treatment modality for chondrosarcoma. For selected patients, however, surgery alone is clearly an inadequate treatment [4, 5]. Rare histologic chondrosarcoma subtypes, many with characteristic clinicopathologic features, share the challenges for care of the common subtypes. The disease demands concerted effort from researchers to undertake collaborative research if outcomes are to improve.
Survival is multifactorial; therefore, no single factor can accurately predict survival outcomes for patients with limb chondrosarcoma. Thus, it would be desirable to establish a statistical prediction model that could integrate all individual prognostic factors to precisely predict the survival of limb chondrosarcoma patients. Nomograms are reliable and convenient statistical tools used to predict the overall probability of a specific outcome in various diseases by incorporating numerous prognostic factors [6, 7]. However, to the best of our knowledge, no previous research has reported a comprehensive nomogram that predicts the overall survival of patients with limb chondrosarcoma. We therefore aim to develop and validate a nomogram for predicting the overall survival of patients with limb chondrosarcomas.
2.1. Data Source and Inclusion Criteria
The present study analyzed the deidentified data obtained from the SEER 18 registry via SEERStat software (version 8.3.5; NCI, Bethesda, USA). The SEER Program is one of the largest registry sources of cancer information supported by the National Cancer Institute of the United States and covers approximately 30% of the U.S. population .
Inclusion criteria were as follows: (1) diagnosis of chondrosarcoma as primary malignancy, (2) positive histological confirmation of chondrosarcoma, (3) site limited to the limb, (4) known survival months after diagnosis and cause of death, and (5) complete follow-up. Exclusion criteria were as follows: (1) chondrosarcoma not the first cancer; (2) clear cell chondrosarcoma or mesenchymal chondrosarcoma; (3) unknown or blank information about tumor size, grade, stage, laterality, and total number of tumors; and (4) unknown use of surgery, radiotherapy, or chemotherapy.
2.2. Prognostic Factors
Data on patient age, sex, race, year of diagnosis, primary site, tumor size, histology, grade, total number of tumors, laterality, stage, use of surgery, use of radiation, use of chemotherapy, and survival time until death or the time of the last follow-up were collected from the SEER database. The age at diagnosis was categorized as younger than 40, 40-59, and 60 years and older . Race was categorized as White, Black, and others (American Indian/Alaskan Native, Asian/Pacific Islander). Histologic subtypes were categorized as conventional, myxoid, and dedifferentiated chondrosarcoma. The pathologic grade was classified according to four categories: grades I, II, III, and IV. Tumor size was classified into three groups (<8, 8-13, and >13 cm) [10, 11]. The total number of tumors was categorized into two groups ( and ). The tumor stage was categorized as localized, regional, and distant. As described in the 2018 version of the Summary Stage Manual provided by SEER (https://seer.cancer.gov/tools/ssm/), localized tumors were defined as a tumor confined to the periosteum. Regional tumors were defined as a tumor that has extended beyond the periosteum without distant metastasis.
2.3. Statistical Analysis
All statistical analyses were completed using R software (version 3.3.0; http://www.r-project.org/). Continuous variables that conformed to a normal distribution were expressed as the values and compared with Student’s -test. Categorical variables were expressed as frequencies and proportions and compared with the chi-squared and Fisher exact tests. The entire database of patients was randomly divided into a training group and a validation group to construct and validate the nomograms via the R caret package (https://cran.r-project.org/package=caret). Survival curves were created using Kaplan-Meier methods for selected variables performed for overall survival and compared with the log-rank test. Univariate and multivariate analyses were performed using Cox proportional hazards regression models to determine the factors associated with survival. Based on the predictive model with identified prognostic factors, a nomogram was constructed to predict the 3- and 5-year survival rates of patients. The nomogram was validated by a concordance index (-index) and calibration curve both internally (training group) and externally (validation group). For all statistical analyses, a value of <0.05 was considered significant.
3.1. Patient Characteristics
A total of 5196 patients diagnosed with chondrosarcoma were registered in the SEER database from 1973 to 2016, of whom 2644 (50.9%) had primary tumors in the limb. According to the inclusion and exclusion criteria, 813 patients were eligible for and were eventually enrolled in the present study and randomly assigned into the training group (, for construction and internal validation of the nomograms) or the validation group (, for external validation of the nomograms) (Figure 1).
Patients characteristics are shown in Table 1. Among these 813 patients, 452 (55.6%) were male and 361 (44.4%) were female. The mean age was 51.51 years (range, 8-90 years), and the median survival time was 51 months (range, 1-155 months). Of the whole population, the majority of patients were White ( (86.3%)), and 668 (82.2%) had conventional chondrosarcoma. The major disease sites were the long bones of the upper and lower limb and associated joints ( (34.6%) and (56.5%), respectively).
Note: LU: long bones: upper limb, scapula, and associated joints; SU: short bones: upper limb, scapula, and associated joints; LL: long bones of lower limb and associated joints; SL: short bones of lower limb and associated joints.
Grade II was the most common of the pathologic grades of chondrosarcoma ( (39.1%)), followed by grade I ( (36.0%)), grade III ( (13.4%)), and grade IV ( (11.4%)). The laterality of tumors was slightly more right than left ( (51.5%) and (48.5%), respectively). Localized and regional were the most common tumor stages, accounting for 61.1% and 31.1%, respectively. The tumor size was less than 8 cm in 471 patients (57.9%), between 8 and 13 cm in 193 (23.7%), and more than 8 cm in 149 (18.3%). The total number of tumors equaling one accounts for 93% (). Moreover, surgery was performed in most patients (, (94.8%)). However, only a small number of patients received chemotherapy and radiotherapy ( (8.6%) and (5.9%), respectively). Among the patients who underwent surgery, 497 of them (61.1%) had localized disease, 253 of them (31.1%) had regional disease, and 63 patients (7.7%) had distant disease. There were no significant differences for variables of patient characteristics between the training and validation groups.
3.2. Prognostic Factors of Overall Survival
According to the univariate analysis, overall survival was significantly related to age, site, histology, grade, stage, use of surgery, radiotherapy or chemotherapy, and tumor size (Table 2). Significant factors identified by univariate analysis were further explored in multivariate analysis, which demonstrated that age, grade, stage, and use of surgery were the independent prognostic factors (Figure 2).
Note: HR: hazard ratio; CI: confidence intervals; LU: long bones: upper limb, scapula, and associated joints; SU: short bones: upper limb, scapula, and associated joints; LL: long bones of lower limb and associated joints; SL: short bones of lower limb and associated joints.
3.3. Kaplan-Meier Curve Analyses
According to the Kaplan-Meier survival and log-rank tests, older age, long bones of lower limbs and associated joints, histologic subtype of dedifferentiated, higher grade, tumor stage of distant, no surgery, use of radiotherapy or chemotherapy, and larger tumor size resulted in worse rates of survival. However, gender, race, laterality, and total number of tumors were not significant for survival (Figures 3 and 4). Notably, the curves were crossed, suggesting that there might be multiple factors of interference.
Based on the result of the Cox regression, we built and interpreted a nomogram for predicting the 3- and 5-year overall survival (Figure 5). The point scale at the top of each nomogram was used first to give each prognostic variable a score; then, the scale at the bottom of each nomogram was used (adding up the scores of all variables) to predict the 3- and 5-year survival probability. There were parallel lines below the figure with linear relationship scales with each other. The nomogram revealed that tumor stage “distant” contributed the highest risk to prognosis, followed by years, grade IV, long bones of lower limb and associated joints, histology subtype of dedifferentiated, no use of surgery, tumor cm, and use of radiotherapy or chemotherapy.
The nomograms were validated internally and externally. In the training group for internal validation, the -index for overall survival prediction was 0.86 (95% confidence interval (CI), 0.81-0.88). In the validation group for external validation, the -index for overall survival prediction was 0.77 (95% CI, 0.73-0.79). Internal and external calibration plots for 3- and 5-year overall survival showed excellent agreement between nomogram predictions and observed outcomes (Figure 6).
Various prognostic factors influence the survival outcome of patients with cancer, and the ability of a single prognostic factor has limited capability to predict individual survival probability. Moreover, relying merely on traditional staging systems is not enough to accurately assess the cancer prognosis [11–13]. Recently, nomograms have been widely used as reliable and convenient tools for the individualized prediction of a patient’s survival outcome. However, to date, there are no studies that report nomograms that predict survival specifically for patients with limb chondrosarcomas. Therefore, utilizing the SEER database which covers approximately 28% of the overall U.S. population, we developed a comprehensive nomogram to predict 3- and 5-year overall survival for patients with limb chondrosarcoma.
The majority of included patients were older than 40 years of age, and the prognosis for survival worsened with increased age, according to the Kaplan-Meier curves and log-rank analyses. Similar results were also observed in previous studies including chondrosarcoma patients in other sites [14–16]. van Praag et al. found that patients older than 50 years of age had worse outcomes . The role of gender in chondrosarcoma remains controversial [18, 19]. In the present study, the ratio of sex was approximately 1.25 to 1. However, we did not identify sex as an independent risk factor in the Cox regression analysis.
In the present study, we did not find a significant difference for the laterality of chondrosarcoma, and the ratio of left to right was 10 to 9. The primary site of chondrosarcoma was more likely to occur in the long bones of the extremities and associated joints, which accounted for 91.0%. We found that patients with chondrosarcoma in long bones had a worse prognosis. However, this should be interpreted with caution because tumors that occurred in other sites had relatively lower incidence. In addition, the Kaplan-Meier curve for the primary site of the tumor was crossed, indicating that there may be multiple factors of interference.
Pathology plays an important role in the prognosis of patients with cancer. Conventional chondrosarcoma of bone constitutes approximately 85% of all chondrosarcomas . In our study, the conventional subtype accounted for 82.2% of all chondrosarcomas. Grade I chondrosarcomas are lowly cellular, with an abundant hyaline cartilage matrix, and rarely metastasize . In contrast, grade III chondrosarcomas are highly cellular, with a myxoid matrix and mitoses, and with metastases developing in 70% of patients . In the present study, we also found that the grade played a crucial role in the survival of patients through Cox regression and the Kaplan-Meier curve and log-rank test. Notably, grading is also subject to interobserver variability, which is worrisome due to the fact that treatment for grade I and II chondrosarcomas may differ.
In addition to conventional subtypes, other rare subtypes constituted 10-15% of all chondrosarcomas . In the present study, rare subtypes accounted for approximately 17% of the chondrosarcomas. Dedifferentiated chondrosarcoma, which accounts for approximately 13% of chondrosarcoma, is diagnosed when a high-grade noncartilaginous tumor borders low-grade disease . Detection of skeletal and/or lung metastases at diagnosis or within 12 to 18 months is common and associated with a dismal prognosis . It was reported that the 5-year survival rate of dedifferentiated chondrosarcoma was less than 20% . In our research, we also found that the dedifferentiated subtype had a worse prognosis compared to conventional and myxoid subtypes.
We demonstrated that the tumor stage was an independent prognostic factor for survival, and the localized stage was related to a significant survival benefit compared to the regional or distant stage. Previous studies reported similar results for chondrosarcoma in other sites or the whole skeletal system [9, 10, 26].
Optimal surgical treatment remains under debate, including systemic chemotherapy and radiotherapy. In line with previous studies, we found that surgery confers a significant benefit in the overall survival of patients with limb chondrosarcoma. Wide, en bloc excision was the preferred surgical treatment for intermediate- and high-grade chondrosarcoma patients. Considering that chondrosarcomas grow slowly, with a relatively low fraction of dividing cells, and radiotherapy acts on dividing cells, chondrogenic tumors are considered relatively radiotherapy-resistant . Radiotherapy can be considered in two situations: after incomplete resection, aiming at maximal local control (curative), and in situations where resection is not feasible or would cause unacceptable morbidity (palliative). Chemotherapy is generally not effective in chondrosarcoma, especially in the most frequently observed conventional type and the rare clear cell variant. Previous studies reported that expression of the multidrug-resistance 1 gene, P-glycoprotein, may contribute to chemotherapy resistance in vivo and in vitro [28, 29]. In our study, we also found that chemotherapy and radiotherapy had limited effects for improving the prognosis of patients.
Based on the independent predictors of overall survival that we identified, we constructed a nomogram to predict 3- and 5-year overall survival probability for patients with limb chondrosarcoma. Our nomogram can be easily used in the clinical practice to estimate a patient’s prognosis. For example, a 55-year-old man was diagnosed with grade II dedifferentiated chondrosarcoma in the long bones of his lower limb and associated joints, with tumor size of 6.0 cm, along with the use of surgery and no use of chemotherapy and radiotherapy and the stage of localized. Totaling the points for this patient, we see that he would have 187.5 points in the overall survival nomogram. The results for estimated 3-year overall survival and cancer-specific survival rates would be 80% and 75%, respectively, according to the nomogram. In our study, we also investigated the performance of the nomogram through the -index and calibration plots internally and externally. They all demonstrated excellent consistency between the prediction model and actual outcomes.
Notably, there are some limitations in the present study. First, the nomogram was developed based on retrospective data. It is inevitable that some patient information was unclear or blanked. Second, due to the unavailability of some variables such as local recurrence in the database, we only considered 3- and 5-year survival as endpoints. However, local recurrence at certain time points could also be evaluated as one of the endpoints. Besides, there was no clear distinction between the type of surgery performed or chemotherapy or radiotherapy regimen. Finally, the nomogram was constructed and validated via the same database. Although we randomly divided patients into the training and validation groups to validate the nomogram internally and externally, it is more reliable to validate the nomogram prospectively or using another data source.
With data from a large population-based group, we were the first to develop and validate a nomogram to provide individualized estimates of rates of 3- and 5-year overall survival in patients with limb chondrosarcomas. The nomograms showed a high degree of accuracy and applicability. Using the nomogram as an effective and convenient tool, clinical practitioners could make personalized survival evaluations and identify patients with a high risk of mortality.
The dataset generated and analyzed during the current study is available in the National Cancer Institute’s Surveillance, Epidemiology, and End Results (SEER) repository.
Analysis of the deidentified data from the Surveillance, Epidemiology, and End Results (SEER) program was exempt from medical ethics review. All procedures performed in studies involving human participants were in accordance with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Conflicts of Interest
The authors declare that there are no conflicts of interest in this work.
Xinjie Wu and Yanlei Wang participated in concept development, data generation, quality control of data, data analysis and interpretation, and writing of the manuscript. Wei Sun and Mingsheng Tan were responsible for data analysis and participated in the interpretation and presentation of the data. Xinjie Wu provided input into data interpretation. Wei Sun and Xinjie Wu were involved in concept development, quality control of data, and data analysis and interpretation. All authors have read and approved the final version of the submitted manuscript. Xinjie Wu and Yanlei Wang contributed equally to this study and should be considered co-first authors.
This study was funded by the National Natural Science Foundation of China (grant nos. 82072524 and 81871830), the Natural Science Foundation of Beijing Municipality (grant no. 7182146), the Fundamental Research Funds for the Central Universities, and the Research Projects on Biomedical Transformation of the China-Japan Friendship Hospital (grant no. PYBZ1828).
- J. S. Whelan and L. E. Davis, “Osteosarcoma, chondrosarcoma, and chordoma,” JCO., vol. 36, no. 2, pp. 188–193, 2018.
- R. L. Jones, D. Katz, E. T. Loggers, D. Davidson, E. T. Rodler, and S. M. Pollack, “Clinical benefit of antiangiogenic therapy in advanced and metastatic chondrosarcoma,” Medical Oncology, vol. 34, no. 10, p. 167, 2017.
- E. Nazeri, M. Gouran Savadkoohi, K. Majidzadeh-A, and R. Esmaeili, “Chondrosarcoma: an overview of clinical behavior, molecular mechanisms mediated drug resistance and potential therapeutic targets,” Critical Reviews in Oncology/Hematology., vol. 131, pp. 102–109, 2018.
- K. D. A. Bernstein, N. Liebsch, Y.-L. Chen et al., “Clinical outcomes for patients after surgery and radiation therapy for mesenchymal chondrosarcomas,” Journal of Surgical Oncology, vol. 114, no. 8, pp. 982–986, 2016.
- Y. Tsuda, S. Evans, J. D. Stevenson et al., “Is the width of a surgical margin associated with the outcome of disease in patients with peripheral chondrosarcoma of the pelvis? A multicenter study,” Clinical Orthopaedics and Related Research, vol. 477, no. 11, pp. 2432–2440, 2019.
- S.-S. Feng, H. Li, F. Fan et al., “Clinical characteristics and disease-specific prognostic nomogram for primary gliosarcoma: a SEER population-based analysis,” Scientific Reports, vol. 9, no. 1, p. 10744, 2019.
- X.-Q. Zheng, J.-F. Huang, J.-L. Lin et al., “Incidence, prognostic factors, and a nomogram of lung cancer with bone metastasis at initial diagnosis: a population-based study,” Transl Lung Cancer Res., vol. 8, no. 4, pp. 367–379, 2019.
- K. M. Doll, A. Rademaker, and J. A. Sosa, “Practical guide to surgical data sets: Surveillance, Epidemiology, and End Results (SEER) database,” JAMA Surgery, vol. 153, no. 6, p. 588, 2018.
- L. Chen, C. Long, J. Liu, X. Duan, and Z. Xiang, “Prognostic nomograms to predict overall survival and cancer-specific survival in patients with pelvic chondrosarcoma,” Cancer Medicine, vol. 8, no. 12, pp. 5438–5449, 2019.
- K. Song, J. Song, X. Shi et al., “Development and Validation of Nomograms Predicting Overall and Cancer-Specific Survival of Spinal Chondrosarcoma Patients,” Spine, vol. 43, no. 21, pp. E1281–E1289, 2018.
- J. Zhang, Z. Pan, F. Zhao et al., “Development and validation of a nomogram containing the prognostic determinants of chondrosarcoma based on the Surveillance, Epidemiology, and End Results database,” International Journal of Clinical Oncology, vol. 24, no. 11, pp. 1459–1467, 2019.
- J. Zhang, Z. Pan, J. Yang, X. Yan, Y. Li, and J. Lyu, “A nomogram for determining the disease-specific survival in Ewing sarcoma: a population study,” BMC Cancer, vol. 19, no. 1, p. 667, 2019.
- K. Song, J. Song, F. Chen, K. Lin, X. Ma, and J. Jiang, “Prognostic nomograms for predicting overall and cancer-specific survival of high-grade osteosarcoma patients,” Journal of Bone Oncology., vol. 13, pp. 106–113, 2018.
- B. J. Miller, “CORR Insights®: survival in mesenchymal chondrosarcoma varies based on age and tumor location: a survival analysis of the SEER database,” Clinical Orthopaedics and Related Research, vol. 475, no. 3, pp. 806-807, 2017.
- B. A. Schneiderman, S. A. Kliethermes, and L. M. Nystrom, “Survival in mesenchymal chondrosarcoma varies based on age and tumor location: a survival analysis of the SEER database,” Clinical Orthopaedics and Related Research, vol. 475, no. 3, pp. 799–805, 2017.
- L.-E. Bohman, M. Koch, R. L. Bailey, M. Alonso-Basanta, and J. Y. K. Lee, “Skull base chordoma and chondrosarcoma: influence of clinical and demographic factors on prognosis: a SEER analysis,” World Neurosurgery, vol. 82, no. 5, pp. 806–814, 2014.
- V. M. van Praag, A. J. Rueten-Budde, V. Ho et al., “Incidence, outcomes and prognostic factors during 25 years of treatment of chondrosarcomas,” Surgical Oncology, vol. 27, no. 3, pp. 402–408, 2018.
- Z. Nie, Q. Lu, and H. Peng, “Prognostic factors for patients with chondrosarcoma: a survival analysis based on the Surveillance, Epidemiology, and End Results (SEER) database (1973–2012),” Journal of Bone Oncology., vol. 13, pp. 55–61, 2018.
- Z. Gao, F. Ren, H. Song et al., “Marital status and survival of patients with chondrosarcoma: a population-based analysis,” Medical Science Monitor, vol. 24, pp. 6638–6648, 2018.
- H. Gelderblom, P. C. W. Hogendoorn, S. D. Dijkstra et al., “The clinical approach towards chondrosarcoma,” The Oncologist., vol. 13, no. 3, pp. 320–329, 2008.
- X. Ma, Y. Dong, C. Zhang, and B. Zeng, “Recurrence analysis in 66 cases with grade I and grade II chondrosarcomas in the extremities,” Orthopaedic Surgery., vol. 1, no. 2, pp. 132–136, 2009.
- C. Cipriano, A. M. Griffin, P. C. Ferguson, and J. S. Wunder, “Developing an evidence-based followup schedule for bone sarcomas based on local recurrence and metastatic progression,” Clinical Orthopaedics and Related Research, vol. 475, no. 3, pp. 830–838, 2017.
- R. J. Grimer, G. Gosheger, A. Taminiau et al., “Dedifferentiated chondrosarcoma: prognostic factors and outcome from a European group,” European Journal of Cancer., vol. 43, no. 14, pp. 2060–2065, 2007.
- R. Capanna, F. Bertoni, G. Bettelli et al., “Dedifferentiated chondrosarcoma,” The Journal of Bone and Joint Surgery. American Volume, vol. 70, no. 1, pp. 60–69, 1988.
- E. L. Staals, P. Bacchini, and F. Bertoni, “Dedifferentiated central chondrosarcoma,” Cancer, vol. 106, no. 12, pp. 2682–2691, 2006.
- A. Y. Giuffrida, J. E. Burgueno, L. G. Koniaris, J. C. Gutierrez, R. Duncan, and S. P. Scully, “Chondrosarcoma in the United States (1973 to 2003): an analysis of 2890 cases from the SEER database,” The Journal of Bone and Joint Surgery. American Volume, vol. 91, no. 5, pp. 1063–1072, 2009.
- D. J. Indelicato, R. L. Rotondo, D. Begosh-Mayne et al., “A prospective outcomes study of proton therapy for chordomas and chondrosarcomas of the spine,” International Journal of Radiation Oncology • Biology • Physics, vol. 95, no. 1, pp. 297–303, 2016.
- J. J. Wyman, A. M. Hornstein, P. A. Meitner et al., “Multidrug resistance-1 and p-glycoprotein in human chondrosarcoma cell lines: expression correlates with decreased intracellular doxorubicin and in vitro chemoresistance,” Journal of Orthopaedic Research, vol. 17, no. 6, pp. 935–940, 1999.
- C. Guo and X. Jin, “Chemoprotection effect of multidrug resistance 1 (MDR1) gene transfer to hematopoietic progenitor cells and engrafted in mice with cancer allows intensified chemotherapy,” Cancer Investigation, vol. 24, no. 7, pp. 659–668, 2006.
Copyright © 2021 Xinjie Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.