Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 28 February 2023
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Editor's Challenge: Abhishek Mahajan - How Can Precision Oncology be Advanced with Validated Imaging-Based Nomograms? View all 8 articles

Development and validation of a MRI-based combined radiomics nomogram for differentiation in chondrosarcoma

Xiaofen Li,&#x;Xiaofen Li1,2†Min Lan&#x;Min Lan3†Xiaolian Wang&#x;Xiaolian Wang2†Jingkun ZhangJingkun Zhang4Lianggeng GongLianggeng Gong5Fengxiang LiaoFengxiang Liao6Huashan LinHuashan Lin7Shixiang DaiShixiang Dai2Bing Fan*Bing Fan2*Wentao Dong*Wentao Dong2*
  • 1Medical College of Nanchang University, Nanchang, China
  • 2Department of Radiology, Jiangxi Provincial People’s Hospital, The First Affiliated Hospital of Nanchang Medical College, Nanchang, China
  • 3Department of Orthopedics, Jiangxi Provincial People’s Hospital, The First Affiliated Hospital of Nanchang Medical College, Nanchang, China
  • 4Department of Radiology, The Affiliated Hospital of Jiangxi University of Traditional Chinese Medicine, Nanchang, China
  • 5Department of Medical Imaging Center, Second Affiliated Hospital of Nanchang University, Nanchang, China
  • 6Department of Nuclear Medicine, Jiangxi Provincial People’s Hospital, The First Affiliated Hospital of Nanchang Medical College, Nanchang, China
  • 7Department of Pharmaceutical Diagnosis, General Electric Healthcare, Changsha, China

Objective: This study aims to develop and validate the performance of an unenhanced magnetic resonance imaging (MRI)-based combined radiomics nomogram for discrimination between low-grade and high-grade in chondrosarcoma.

Methods: A total of 102 patients with 44 in low-grade and 58 in high-grade chondrosarcoma were enrolled and divided into training set (n=72) and validation set (n=30) with a 7:3 ratio in this retrospective study. The demographics and unenhanced MRI imaging characteristics of the patients were evaluated to develop a clinic-radiological factors model. Radiomics features were extracted from T1-weighted (T1WI) images to construct radiomics signature and calculate radiomics score (Rad-score). According to multivariate logistic regression analysis, a combined radiomics nomogram based on MRI was constructed by integrating radiomics signature and independent clinic-radiological features. The performance of the combined radiomics nomogram was evaluated in terms of calibration, discrimination, and clinical usefulness.

Results: Using multivariate logistic regression analysis, only one clinic-radiological feature (marrow edema OR=0.29, 95% CI=0.11-0.76, P=0.012) was found to be independent predictors of differentiation in chondrosarcoma. Combined with the above clinic-radiological predictor and the radiomics signature constructed by LASSO [least absolute shrinkage and selection operator], a combined radiomics nomogram based on MRI was constructed, and its predictive performance was better than that of clinic-radiological factors model and radiomics signature, with the AUC [area under the curve] of the training set and the validation set were 0.78 (95%CI =0.67-0.89) and 0.77 (95%CI =0.59-0.94), respectively. DCA [decision curve analysis] showed that combined radiomics nomogram has potential clinical application value.

Conclusion: The MRI-based combined radiomics nomogram is a noninvasive preoperative prediction tool that combines clinic-radiological feature and radiomics signature and shows good predictive effect in distinguishing low-grade and high-grade bone chondrosarcoma, which may help clinicians to make accurate treatment plans.

Introduction

Chondrosarcoma is a malignant bone tumor characterized by the production of cartilage matrix, which accounts for about 20% of all primary malignant bone tumors (1). After myeloma and osteosarcoma, chondrosarcoma is the third most common primary bone malignant tumor. Chondrosarcoma occurs in a wide range of ages, more males than females, mostly in long tubular bones, mostly in the proximal femur and humerus, and rarely in the intracranial, larynx, sacrum, and sternum. According to the 5th edition of WHO Classification of Bone Neoplasms in 2020, conventional chondrosarcoma is classified into grade I, grade II, and grade III (2, 3). Different grades of chondrosarcoma have widely different treatment methods and prognosis (4). The main clinical treatment is tumor resection to prevent recurrence and distant metastasis (5). Studies have confirmed that low-grade chondrosarcoma does not metastasize, and the 10-year survival rate of grade I chondrosarcoma is about 77-89% (6). Therefore, curettage of the lesion is safe and effective. About 10-33% of grade II chondrosarcomas and 70% of grade III chondrosarcomas can metastasize, with 10-year survival rates of about 53-59% and 36-38%, respectively, requiring extensive resection or amputation (7). Due to the heterogeneity of the tumor and biopsy of small sampling, and the gross results have a certain difference, can not accurately predict the tumor differentiation degree, affect the clinical risk of illness misjudgment leads to increased postoperative recurrence or metastasis, therefore, noninvasive accurate preoperative assessment of the precision medical treatment and prognosis in patients with chondrosarcoma is particularly important.

On conventional X-ray and CT, puncture, annular, semi-annular, or gravel calcification were seen in the solid portion of the chondrosarcoma, and scallop-like changes were seen in the adjacent cortical bone (8). Due to its high soft tissue resolution, MRI has been widely used in the diagnosis and differential diagnosis of chondrosarcoma, especially for the morphological characteristics of the tumor, such as marrow edema, soft tissue mass, periosteal reaction, etc. (9). However, these features extracted from traditional imaging have certain limitations in predicting the differentiation degree of chondrosarcoma (10).

In recent years, radiomics through different kinds of machine algorithm to extract the quantitative features of high flux, and mining related to tumor heterogeneity, or the degree of differentiation of biological information, has been widely applied in glioma, the preoperative differentiation degree of prediction and prognosis of gastric cancer, colorectal cancer evaluation (1114). Meanwhile, some researches focused on identification of chondrosarcoma and enchondroma radiomics study (15). The objective of our study was to develop and validate a combined radiomics nomogram that combines unenhanced MRI radiomic features with clinic-radiological features for preoperative differentiation of low-grade (grade I) and high-grade (grade II, III) chondrosarcomas.

Materials and methods

Patients

Institutional Review Board approval and a waiver for informed consent were obtained. This retrospective research included patients with low-to-high grade chondrosarcoma who underwent MRI, between June 2015 and January 2022, at our Hospital. The inclusion criteria were as follows: (1) primary low-grade and high-grade chondrosarcoma that underwent surgery, such as intralesional curettage or resection. (2) the definitive histological diagnosis was based on the evaluation of the surgical specimen and was regarded as a reference standard. (3) MRI was performed within three months before surgery and T1WI and T2WI sequences were provided. The exclusion criteria were as follows: (1) patients received radiotherapy or chemotherapy before MRI scans. (2) poor imaging quality making difficulties in segmentation. Finally, a total of 102 patients with low-grade (n=44, 17 male and 27 female) and high-grade (n=58, 31 male and 27 female) chondrosarcoma were included in the present study.

44 low-grade chondrosarcoma patients. These lesions were located in the rib (n = 5), sternum (n = 2), pelvis (n = 4), humerus (n=8), tibia (n=8), femur (n=10), radius (n=4) and scapula (n = 3), respectively.

58 high-grade chondrosarcoma patients. These lesions were located in the sternum (n=2), pelvis (n=16), humerus (n=15), tibia (n=3), femur (n=16), scapula (n = 3), calcaneus (n=1), sacrum (n=1) and patella (n=1), respectively.

Patients were randomly assigned to the training and validation sets in a ratio of approximately 7:3. The training and validation sets were stratified to maintain the same proportion of low-grade and high-grade tumors in the training and validation groups.

MRI image acquisition

All enrolled patients underwent 3.0T MRI (SIMENS MAGNETOM Skyra) with at least T1WI and T2WI MRI sequences. T1WI scanning parameters: TR619ms, TE12ms, F0V390mm, matrix 384×384, layer thickness 5mm, interval 1mm, reverse Angle 160°, echo chain length 3. T2WI scanning parameters: TR3400ms, TE85ms, FOV380mm, matrix 384×384, layer thickness 5mm, interval 1mm, reverse Angle 160°, echo chain length 15.

MRI characteristic evaluation

The MRI image were scrutinized by two radiologists with 5 years (doctor 1, D.W.) and 15 years (doctor 2, L.X.) of diagnostic bone imaging experience. Blinded to the clinic-pathologic data, the two doctors interpreted the following MRI features by consensus: “Longest diameter” were the longest diameter of the tumor on an axial MRI image; Periosteal reaction (present or not; “periosteal reaction” were fluid signal intensity on T2W sequences on the surface of the bone); Marrow edema (present or not; “marrow edema” were patch of hyperintensity around the tumor on the T2W sequences); Soft tissue mass: (present or not; “soft tissue mass” were intramedullary cartilage tumors extend directly through the damaged cortex and have the same radiographic features and well-defined peripheral margins as intramedullary tumors). If there was disagreement, the two radiologists needed to reach a consensus.

Construction of the clinic-radiological model

The clinical parameters of gender and age were retrieved from the electronic medical record system of our hospital. Two radiologists with 5 and 15 years of experience in bone imaging were blinded to imaging reports and pathological details, and imaging features including longest diameter, periosteal reaction, marrow edema and soft tissue mass were reviewed and reported. When constructing the clinic-radiological model, the above clinic-radiological characteristics were analyzed by univariate regression, and then the statistically significant characteristics in the univariate regression analysis were processed by multivariate regression model. Ultimately, features with P values less than 0.05 were used to establish the clinic-radiological model.

Tumor segmentation

All T1-weighted MRI images in DICOM format, original size and resolution were transferred to ITK-SNAP software (Version 3.8, www.itksnap.org) for three-dimensional (3D) region of interest (ROI) segmentation. To ensure accurate tumor boundaries, ROIs on all slices were carefully delineated manually by a radiologist (doctor 1, D.W.) with 5 years of experience in bone imaging, who was blinded to the pathological findings. To test the stability of features, doctor 1 and another radiologist with 15 years of experience (doctor 2, L.X.) underwent re-extraction of radiomic features from 40 randomly selected patients from the entire study set. Intra-class correlation coefficient (ICC) was calculated to evaluate the consistency and reproducibility of the features. Subsequent analyses included features of ICC>0.75 in intra-observer and inter-observer consistency analyses. In order to avoid local volume effect, the top layer and bottom layer were eliminated. Figure 1 shows two typical T1WI MRI images of chondrosarcoma segmentation on ITK-SNAP, one of which is low-grade (Figures 1A, B) and the other of which is high-grade (Figures 1C, D).

FIGURE 1
www.frontiersin.org

Figure 1 Two representative T1WI MRI segmentation on ITK-SNAP. (A, B) A patient with low-grade chondrosarcoma of the tibia. (C, D) A patient with high-grade chondrosarcoma of the tibia.

Radiomics feature extraction

Prior to radiomics feature extraction, to reduce feature variability, we performed the following image preprocessing steps, containing gray discretization, intensity normalization and voxel resampling. Then, radiomics features were extracted from T1WI MRI images through the open source PyRadiomics library. They are divided into four categories: size and morphology features, descriptors of image intensity histograms, descriptors of the relationship between image voxels and higher-order texture features extracted from filtered images.

Construction of the radiomics signature

In order to prevent signature overfitting, dimensionality of features is reduced before signature construction. Succinctly, radiomics features that met the inter-observer and intra-observer ICCs criteria greater than 0.75 and were significantly different between the two groups as assessed by one-way analysis of variance (ANOVA) were included in the LASSO regression model to select the most valuable features in the training set. Finally, the selected radiomics features were used to construct radiomics signature. Rad-score was calculated for each patient by linear combination of selected features and weighted by the respective LASSO coefficients.

Development and validate of combined radiomics nomogram

A combined radiomics nomogram was developed by combining significant variables of clinic-radiological characteristics and Rad-score. The calibration curve was used to evaluate the calibration of nomogram. Hosmer-Lemeshow test evaluated the goodness of fit of nomogram. The ROC curve of the training and validation sets and DeLong tests were used to evaluate the diagnostic performance of the clinic-radiological model, radiomics signature and combined radiomics nomogram in differentiating low-grade from high-grade in chondrosarcoma. DCA was performed by calculating the net benefits over a range of threshold probabilities across the entire set, in order to assess the clinical usefulness of the nomogram.

Statistical analysis

Statistical analysis was conducted using SPSS 19.0 and R software (Version 3.4.4; http://www.Rproject.org). The chi-square test was used for the analysis of categorical variables and the Mann-Whitney U test was used for the analysis of continuous variables. The clinic-radiological factor and Rad-score were correlated using Pearson’s correlation coefficient. The “glmnet (R)” package was used to perform LASSO regression. The “Regression Modeling Strategy (RMS)” package was used to construct the nomogram and calibration curves. The Hosmer–Lemeshow test was performed on the”generalhoslem” package. ROC curves were plotted using the “partial Receiver Operating Characteristic (pROC)” software package. Delong test was used to estimate the difference of AUC values between different models (clinic-radiological model, radiomics signature, combined radiomics nomogram). Risk Model Decision Analysis (RMDA) software package was used to analyze DCA. The significance level was set as two-sided p < 0.05.

Results

Clinic-radiological characteristics and construction of the clinic-radiological model

A total of 102 patients were enrolled in our study, and clinical and radiological data were collected. Table 1 summarizes the differences in clinical radiological variables between low-grade and high-grade chondrosarcoma patients in the training and validation sets. The training set including 72 patients (36 males and 36 females), of whom 31 were diagnosed with low-grade chondrosarcoma and 41 with high-grade chondrosarcoma. Patients with high-grade chondrosarcoma was only significantly different from those with low-grade chondrosarcoma in terms of marrow edema (p=0.02 in the training set). Univariate analysis showed that marrow edema and periosteal reaction served as the risk factors of low-grade and high-grade in chondrosarcoma. After multivariate logistic regression analysis, only marrow edema (OR=0.29, 95% CI=0.11-0.76, p=0.012) remained an independent clinic-radiological predictor (Table 2).

TABLE 1
www.frontiersin.org

Table 1 Clinic-radiological characteristics of chondrosarcoma patients in the training and validation sets.

TABLE 2
www.frontiersin.org

Table 2 Univariate and multivariate logistic regression analysis of the clinic-radiological features in predicting the high-grade and low-grade chondrosarcoma.

Feature extraction, selection, and radiomics signature building

A total of 1316 radiomics features were extracted from T1WI MRI images of each chondrosarcoma patient, among which 829 features were proved to have good inter-observer and intra-observer agreement, which ICCs achieve greater than 0.75. The significant difference between high-grade and low-grade chondrosarcoma in one-way ANOVA (P < 0.05) was enrolled into the LASSO logistic regression model to select the most valuable radiomics features (Figures 2A, B). Finally, four radiomics features were used to construct radiomics signature (Figure 2C).

FIGURE 2
www.frontiersin.org

Figure 2 Radiomics features were selected using the least absolute shrinkage and selection operator (LASSO) regression. (A) Selection of tuning parameter (λ) in LASSO model. Selected the optimal value with λ= 0.025 and log(λ) =-3.605. (B) LASSO coefficient profiles of radiomics features. The coefficient profiles corresponding to the selected logarithm (λ) values were generated by five-fold cross validation. (C) Selected radiomics features and their coefficients.

Development and validate of combined radiomics nomogram

The clinic-radiological factor (marrow edema) and Rad-score were incorporated into construction of the MRI-based combined radiomics nomogram (Figure 3A). After spearman correlation analysis, the result was that R values were all very small, and P values were all > 0.05, which indicated that there was no correlation between marrow edema and Rad-score, and they were used to construct the combined radiomics nomogram from different perspectives. The calibration curve (Figure 3B) showed a good agreement between the predicted and actual probabilities for predicting the low-grade and high-grade chondrosarcoma in the training and validation sets, and the Hosmer-Lemeshow test yielded a nonsignificant statistical difference (P =0.169 and 0.317).

FIGURE 3
www.frontiersin.org

Figure 3 (A) Integrating marrow edema and radiomics signature, the MRI-based nomogram was established. (B) For the calibration curve of the nomogram, the closer the fit between the red diagonal and the ideal dashed line, the higher the prediction accuracy of the nomogram.

Table 3 and Figure 4 present the clinic-radiological, radiomics signature and combined radiomics nomogram diagnostic performance and ROC curves in the training and validation sets, respectively. The AUC of the clinic-radiological model was 0.65 (95%CI 0.54-0.76) in the training set and 0.54 (95%CI 0.36-0.73) in the validation set, while the radiomics signature yielded an AUC value of 0.75 (95%CI 0.63-0.86) and 0.73 (95%CI 0.55-0.91) in both sets. The combined radiomics nomogram achieved the optimal discrimination in the training (AUC, 0.78; 95% CI, 0.67-0.89) and validation (AUC, 0.77; 95% CI, 0.59-0.94) sets, with accuracy of 0.750 and 0.700, sensitivity of 0.805 and 0.765, and specificity of 0.677 and 0.615, respectively. Using the Delong test (Table 4), there was a significant difference between the clinic-radiological and the MRI-based combined radiomics nomogram AUC in the training set (P = 0.013) and the validation set (P =0.037), but the difference did not reach statistical significance between nomogram and radiomics (P=0.337, 0.491) and between clinic-radiological and radiomics (P=0.230, 0.208), respectively, in the training set and validation set. The DCA was shown in Figure 5. The combined radiomics nomogram demonstrated a higher overall net benefit than radiomics model, indicating that the combined radiomics nomogram had an excellent clinical application value in distinguishing high-grade from low-grade chondrosarcoma.

TABLE 3
www.frontiersin.org

Table 3 Predictive performance of clinic-radiological model, radiomics signature, and combined radiomics nomogram.

FIGURE 4
www.frontiersin.org

Figure 4 ROC curves (AUC) of the three models in the training (left) and validation sets (right).

TABLE 4
www.frontiersin.org

Table 4 Comparison of the prediction with the combined radiomics nomogram, clinic-radiological model, and radiomics signature.

FIGURE 5
www.frontiersin.org

Figure 5 Decision curve analysis (DCA) for the three models. The DCA showed that more net benefits were obtained using the combined radiomics nomogram approach over most threshold probability ranges.

Discussion

In this retrospective study, we developed and validated a combined radiomics nomogram for noninvasive, individualized prediction of low-grade and high-grade chondrosarcoma. The combination of clinic-radiological features and T1WI MRI radiomics signature showed predictive efficiency (AUC=0.77, 95%CI =0.59-0.94) in distinguishing chondrosarcoma patients with satisfactory reproducibility and reliability.

Accurate grading of chondrosarcoma is an urgent need to select the most appropriate treatment. For low-grade lesions, conservative treatment is recommended, and for high-grade lesions, aggressive treatment is required (4, 16). However, preoperative biopsies may mistakenly lead to a reduction in chondrosarcoma grade because only small lesion areas are sampled, interobserver differences in tumor grade can occur even among professional bone pathologists (17, 18). Furthermore, inaccurate preoperative grading may lead to inadequate treatment, subsequent need for further surgery, and increased morbidity.

Before surgery, imaging plays a crucial role, and MRI is the preferred method (9). H. Douis et al. found that periosteal reaction and soft tissue mass on MRI may indicate a high degree of malignancy in patients with chondrosarcoma (19). Our qualitative image analysis was less consistent with this finding, only marrow edema was proved to be independent factors of high-grade chondrosarcomas (OR=0.29, 95% CI=0.11-0.76, P=0.012), probably due to the sample of his research and ours were not large and statistical bias may occur.

We only extract tumor texture features on the T1WI sequence, because the boundary of the tumor is easy to delineate on the T1WI sequence. T1WI sequence which is the most generalized and has the most stable image quality was adopted as the research sequence, and all T1WI images were from the same device. We only used T1WI sequence to have a radiomics model, which had a similar diagnostic efficacy, indicating that our study was very effective. In addition, T2WI sequence and diffusion-weighted MRI have been demonstrated unable to distinguish low-grade from high-grade chondrosarcomas (10, 20). Besides, contrast-enhanced MRI (CE-MRI) is helpful in diagnosis because of the rapid enhancement of high-grade tumor areas due to rich vascularized intralesional septations (21, 22). However, CE-MRI were not evaluated in our research because not all patients underwent CE-MRI scans.

At present, imaging techniques could be further improved to better grade and accurately diagnose chondrosarcoma before surgery (23). With the development of artificial intelligence, radiomics can provide valuable assistance by providing quantitative data to integrate existing qualitative image information (24). Recently, the radiomics studies of chondrosarcoma mainly focus on the differentiation of enchondroma from chondrosarcoma, chondrosarcoma recurrence or survival evaluation, limited studies have focused on radiomics and chondrosarcoma differentiation and grading. To our knowledge, Gitto S and his team did some study on chondrosarcoma differentiation (10), they evaluated 58 patients with low-grade and high-grade chondrosarcoma and observed that MRI-based volumetric texture analysis could distinguish these two lesions by many separate texture features. Recently, they also investigated the diagnostic performance of machine learning model based on MRI-based radiomics in distinguishing atypical cartilaginous tumour (ACT) from grade II chondrosarcoma in long bones (25). However, these studies established a machine learning model based only on radiomics features, ignoring the importance of clinical and radiological features. For computerized clinical decision support systems, radiomics-derived data are not a panacea. Our research summarizes the subtle differences between clinical radiologic features and MRI radiomics features in patients with chondrosarcoma and analyses them with nomogram. We chose to investigate several indicators in terms of clinical information that could independently predict the grade of differentiation.

Limitations of our study deserve consideration. First, this is a retrospective study, the performance of the training model is limited by the sample size and only comes from single-center tertiary hospitals, which lacks external validation. It is necessary to carry out large-sample and multicenter studies in the future, and the efficiency of the model should be further tested. Second, more clinical factors and MRI sequence need to incorporated into our nomogram, especially contrast-enhanced T1WI sequence. Future studies will focus on improving the predictive power by combining multiparametric MRI data with other clinical factors and validating the predictive nomogram on multicenter data.

In conclusion, our findings suggest that combined radiomics nomogram based on T1WI sequential MRI have great potential as biomarkers for differentiating low-grade from high-grade chondrosarcomas. As a non-invasive, preoperative method, the combined radiomics nomogram may helpful for the differentiation and clinical decision-making of chondrosarcoma patients.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by the ethics committee of the Jiangxi Provincial People’s Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

The authors made the following contributions: WD and XL made the conception for this research. Data collection and analysis were performed by WD, ML, XW, JZ, FL and SD. WD, ML and XW analyzed the data and drafted the article. XL, BF, LG and HL reviewed/edited the manuscript. All the authors critically revised the article for important intellectual content. All authors contributed to the article and approved the submitted version.

Conflict of interest

Author HL was employed by company General Electric Healthcare.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Thorkildsen J, Taksdal I, Bjerkehagen B, Haugland HK, Børge Johannesen T, Viset T, et al. Chondrosarcoma in Norway 1990-2013; an epidemiological and prognostic observational study of a complete national cohort. Acta Oncol (2019) 58(3):273–82. doi: 10.1080/0284186X.2018.1554260

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Anderson WJ, Doyle LA. Updates from the 2020 world health organization classification of soft tissue and bone tumours. Histopathology (2021) 78(5):644–57. doi: 10.1111/his.14265

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Choi JH, Ro JY. The 2020 WHO classification of tumors of bone: An updated review. Adv Anat Pathol (2021) 28(3):119–38. doi: 10.1097/PAP.0000000000000293

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Dierselhuis EF, Goulding KA, Stevens M, Jutte PC. Intralesional treatment versus wide resection for central low-grade chondrosarcoma of the long bones. Cochrane Database Syst Rev (2019) 3(3). doi: 10.1002/14651858.CD010778.pub2

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ferguson JL, Turner SP. Bone cancer: Diagnosis and treatment principles. Am Fam Phys (2018) 98(4):205–13.

Google Scholar

6. Amer KM, Munn M, Congiusta D, Abraham JA, Basu Mallick A. Survival and prognosis of chondrosarcoma subtypes: SEER database analysis. J Orthop Res (2020) 38(2):311–9. doi: 10.1002/jor.24463

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tap WD, Villalobos VM, Cote GM, Burris H, Janku F, Mir O, et al. Phase I study of the mutant IDH1 inhibitor ivosidenib: Safety and clinical activity in patients with advanced chondrosarcoma. J Clin Oncol (2020) 38(15):1693–701. doi: 10.1200/JCO.19.02492

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Jang BG, Huh KH, Kang JH, Kim JE, Yi WJ, Heo MS, et al. Imaging features of chondrosarcoma of the temporomandibular joint: Report of nine cases and literature review. Clin Radiol (2020) 75(11):878.e1–878.e12. doi: 10.1016/j.crad.2020.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sharif B, Rajakulasingam R, Sharifi S, O'Donnell P, Saifuddin A. MRI Features of low-grade and high-grade chondrosarcoma in enchondromatosis. Skeletal Radiol (2021) 50(8):1637–46. doi: 10.1007/s00256-021-03718-7

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Gitto S, Cuocolo R, Albano D, Chianca V, Messina C, Gambino A, et al. MRI Radiomics-based machine-learning classification of bone chondrosarcoma. Eur J Radiol (2020) 128:109043. doi: 10.1016/j.ejrad.2020.109043

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Mayerhoefer ME, Materka A, Langs G, Häggström I, Szczypiński P, Gibbs P, et al. Introduction to radiomics. J Nucl Med (2020) 61(4):488–95. doi: 10.2967/jnumed.118.222893

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Gore S, Chougule T, Jagtap J, Saini J, Ingalhalikar M. A review of radiomics and deep predictive modeling in glioma characterization. Acad Radiol (2021) 28(11):1599–621. doi: 10.1016/j.acra.2020.06.016

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sah BR, Owczarczyk K, Siddique M, Cook GJR, Goh V. Radiomics in esophageal and gastric cancer. Abdom Radiol (NY). (2019) 44(6):2048–58. doi: 10.1007/s00261-018-1724-8

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Wang R, Dai W, Gong J, Huang M, Hu T, Li H, et al. Development of a novel combined nomogram model integrating deep learning-pathomics, radiomics and immunoscore to predict postoperative outcome of colorectal cancer lung metastasis patients. J Hematol Oncol (2022) 15(1):11. doi: 10.1186/s13045-022-01225-3

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Pan J, Zhang K, Le H, Jiang Y, Li W, Geng Y, et al. Radiomics nomograms based on non-enhanced MRI and clinical risk factors for the differentiation of chondrosarcoma from enchondroma. J Magn Reson Imaging (2021) 54(4):1314–23. doi: 10.1002/jmri.27690

PubMed Abstract | CrossRef Full Text | Google Scholar

16. MacDonald IJ, Lin CY, Kuo SJ, Su CM, Tang CH. An update on current and future treatment options for chondrosarcoma. Expert Rev Anticancer Ther (2019) 19(9):773–86. doi: 10.1080/14737140.2019.1659731

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hodel S, Laux C, Farei-Campagna J, Götschi T, Bode-Lesniewska B, Müller DA. The impact of biopsy sampling errors and the quality of surgical margins on local recurrence and survival in chondrosarcoma. Cancer Manag Res (2018) 2018:3765–71. doi: 10.2147/CMAR.S178768

CrossRef Full Text | Google Scholar

18. Deckers C, de Leijer EM, Flucke U, de Rooy JWJ, Schreuder HWB, Dierselhuis EF, et al. Curettage and cryosurgery for enchondroma and atypical cartilaginous tumors of the long bones: Oncological results of a large series. J Surg Oncol (2021) 123(8):1821–7. doi: 10.1002/jso.26457

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Douis H, Singh L, Saifuddin A. MRI Differentiation of low-grade from high-grade appendicular chondrosarcoma. Eur Radiol (2014) 24(1):232–40. doi: 10.1007/s00330-013-3003-y

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Douis H, Jeys L, Grimer R, Vaiyapuri S, Davies AM. Is there a role for diffusion-weighted MRI (DWI) in the diagnosis of central cartilage tumors? Skeletal Radiol (2015) 44(7):963–9. doi: 10.1007/s00256-015-2123-7

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Geirnaerdt MJ, Hogendoorn PC, Bloem JL, Taminiau AH, van der Woude HJ. Cartilaginous tumors: Fast contrast-enhanced MR imaging. Radiology (2000) 214(2):539–46. doi: 10.1148/radiology.214.2.r00fe12539

PubMed Abstract | CrossRef Full Text | Google Scholar

22. De Coninck T, Jans L, Sys G, Huysse W, Verstraeten T, Forsyth R, et al. Dynamic contrast-enhanced MR imaging for differentiation between enchondroma and chondrosarcoma. Eur Radiol (2013) 23(11):3140–52. doi: 10.1007/s00330-013-2913-z

PubMed Abstract | CrossRef Full Text | Google Scholar

23. van de Sande MAJ, van der Wal RJP, Navas Cañete A, van Rijswijk CSP, Kroon HM, Dijkstra PDS, et al. Radiologic differentiation of enchondromas, atypical cartilaginous tumors, and high-grade chondrosarcomas-improving tumor-specific treatment: A paradigm in transit? Cancer (2019) 125(19):3288–91. doi: 10.1002/cncr.32404

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Tunali I, Gillies RJ, Schabath MB. Application of radiomics and artificial intelligence for lung cancer precision medicine. Cold Spring Harb Perspect Med (2021) 11(8):a039537. doi: 10.1101/cshperspect.a039537

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gitto S, Cuocolo R, van Langevelde K, van de Sande MAJ, Parafioriti A, Luzzati A, et al. MRI Radiomics-based machine learning classification of atypical cartilaginous tumour and grade II chondrosarcoma of long bones. EBioMedicine (2022) 75:103757. doi: 10.1016/j.ebiom.2021.103757

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: chondrosarcoma, radiomics, nomogram, magnetic resonance imaging, T1-weighted image, prediction tool

Citation: Li X, Lan M, Wang X, Zhang J, Gong L, Liao F, Lin H, Dai S, Fan B and Dong W (2023) Development and validation of a MRI-based combined radiomics nomogram for differentiation in chondrosarcoma. Front. Oncol. 13:1090229. doi: 10.3389/fonc.2023.1090229

Received: 05 November 2022; Accepted: 13 February 2023;
Published: 28 February 2023.

Edited by:

Manish Devendra Mair, University Hospitals of Leicester NHS Trust, United Kingdom

Reviewed by:

Emma D’Ippolito, University of Campania Luigi Vanvitelli, Italy
Zhenzhen Gao, Second Hospital of Jiaxing City, China

Copyright © 2023 Li, Lan, Wang, Zhang, Gong, Liao, Lin, Dai, Fan and Dong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wentao Dong, NDYyNDMyNzE2QHFxLmNvbQ==; Bing Fan, MjYxNzEzODFAcXEuY29t

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.