Skip to main content
Download PDF
Download ePub
  • Research article
  • Open access
  • Published:

Intratumoral and peritumoral MRI-based radiomics prediction of histopathological grade in soft tissue sarcomas: a two-center study

Cancer Imaging volume 23, Article number: 103 (2023) Cite this article

Abstract

Objectives

This study aims to develop a model based on intratumoral and peritumoral radiomics from fat-suppressed T2-weighted(FS-T2WI) images to predict the histopathological grade of soft tissue sarcoma (STS).

Methods

This retrospective study included 160 patients with STS from two centers, of which 82 were low-grade and 78were high-grade. Radiomics features were extracted and selected from the region of tumor mass volume (TMV) and peritumoral tumor volume (PTV) respectively. The TMV, PTV, and combined(TM-PTV) radiomics models were established in the training cohort (n = 111)for the prediction of histopathological grade. Finally, a radiomics nomogram was constructed by combining the TM-PTV radiomics signature (Rad-score) and the selected clinical-MRI predictor. The ROC and calibration curves were used to determine the performance of the TMV, PTV, and TM-PTV models in the training and validation cohort (n = 49). The decision curve analysis (DCA) and calibration curves were used to investigate the clinical usefulness and calibration of the nomogram, respectively.

Results

The TMV model, PTV model, and TM-PTV model had AUCs of 0.835, 0.879, and 0.917 in the training cohort and 0.811, 0.756, 0.896 in the validation cohort. The nomogram, including the TM-PTV signatures and peritumoral hyperintensity, achieved good calibration and discrimination with a C-index of 0.948 (95% CI, 0.906 to 0.990) in the training cohort and 0.921 (95% CI, 0.840 to 0.995) in the validation cohort. Decision curve analysis demonstrated the clinical usefulness of the nomogram.

Conclusion

The proposed model based on intratumoral and peritumoral radiomics showed good performance in distinguishing low-grade from high-grade STSs.

Key points

• Peritumoral and intratumoral radiomics models based on FS-T2WI images had a good predictive ability for the histopathological grading of STSs.

• Peritumoral hyperintensity and radiomic scores were independently associated with histopathological grading of STSs.

• The nomogram integrating clinic features and radiomics scores demonstrated good performance for predicting the histopathological grading of STSs.

Soft tissue sarcomas (STSs) are highly heterogeneous malignant tumors originating from mesenchymal tissue. According to the histological grading method used by the French Federation of Cancer Centers Sarcoma Group(FNCLCC) [1], STSs are divided into classes I–III based on how aggressively malignant cells manifest themselves, where grade I is a low grade and grade II, III, is a high grade. Different histopathological grades of treatment strategies and prognoses have variances [2]. Preoperative adjuvant chemoradiotherapy is recommended for high-grade STSs, which can help improve the survival rate of patients. Correspondingly, the side effects of preoperative chemoradiotherapy can be avoided if patients with low-grade STS are identified before surgery. At present, the preoperative histopathological grading diagnosis of STSs mainly relies on the core needle biopsy [3]. However, the results are easily influenced by the sampling site, the size, and the makeup of the lesion, making it impossible to do an overall assessment of the lesion [4].

The most common non-invasive technique for STS preoperative diagnosis and evaluation is magnetic resonance imaging (MRI). Although experienced radiologists can easily identify tumors from MRI, tumor heterogeneity makes grading STSs challenging [5]. Recently, there has been promising progress in the grading of STSs based on radiomic features of MRI. Zhang et al. [6] employed radiomics-based features to establish different diagnostic models to identify the grading of STSs and found the model using the support vector machine (SVM) classifier method performed best. Peeken et al. [7] established independent and combined radiomics models based on FS-T2WI and T1WI-enhanced sequences to predict STSs histopathological grades and found that the radiomics model based on FS-T2WI had the highest predictive performance for high-grade STSs. Yan et al. [6, 8] constructed a radiomics nomogram method to predict high-grade STSs, conducted model development and validation with 180 cases of STSs in two centers, and found that the nomogram based on radiomics features, T staging, and MR boundaries was superior to a single radiomics model or a clinical feature model. Although previous studies have achieved great success in differentiating pathological grades, extracting more valuable radiomics features to boost prediction accuracy remains a challenge.

However, the existing radiomics-based approaches focus on the intratumoral area and ignore the role of the peritumoral environment in STSs grading. Endothelial cells, fibroblasts, immune cells, and other cell types as well as extracellular components make up the peritumoral area, also known as the tumor microenvironment [9, 10]. The microenvironment determines many aspects of tumor behavior, including tumor progression, treatment response, and metastasis [11]. White et al. have shown that satellite-like single or clustered tumor cells that are not visible on imaging can be found beyond the tumor margins in a population of patients undergoing surgery [12]. Clinical evidence suggested the heterogeneity of STSs is not limited to tumor margins but also involves peritumoral regions [13]. As a result, the tumor’s peritumoral environment is also promising and may offer important data for the clinical evaluation of tumor invasive biological behavior. Recently, there have been some studies combining intratumoral and peritumoral radiomic features to determine the histopathological classification of clear cell renal cell carcinoma [14, 15], identify benign and malignant nodules in the lung [16, 17], predict lymph node metastasis and distant metastasis in lung cancer [18], and predict the risk of breast cancer response to chemotherapy [19], which achieved good performance.

The purpose of this study was to investigate, using a two-center dataset, the capacity of intratumoral and peritumoral radiomics signatures based on MRI to noninvasively predict STS histopathological grade.

Materials and methods

Patients and MRI morphologic characteristics

This two-center retrospective study’s ethical approval was provided by the two institutional review boards, and the requirement for informed consent was waived. Between June 2016 and July 2022, 160 patients with STS confirmed by pathology and met the inclusion criteria were retrospectively collected. Inclusion criteria:(1) Patients underwent surgical resection; (2)STS was diagnosed by histopathology; (3) Axial FS-T2WI MRI scans ≤ 2 weeks before surgery; Exclusion criteria: (1) Incomplete clinical or imaging data; (2) MRI image quality is poor, signal-to-noise ratio ≤ 1.0; (3) Development of other subsequent tumors; (4) The patient has received prior treatment, such as chemotherapy, radiation therapy, or needle biopsy.

In total, 160 STS patients were analyzed. Clinical-MRI characteristics included age, gender, location, and MRI morphological features. All images were independently reviewed by two radiologists with more than 5 years of skeletal muscle MRI experience while remaining blind to the clinical and histopathological data. Decisions on MRI findings were made through team negotiation. According to Zhao et al. [20], the following MRI morphological features were selected for comparison: (1) size (maximum diameter of tumor, < 5 cm or ≥ 5 cm); (2) margin (well- or poorly-defined); (3) signal intensity (homogeneous or heterogeneous, > 30% of the whole volume was considered heterogeneous); (4) peritumoral hyperintensity. All these MRI features were labeled as dichotomous variables and recorded using Yes or No.

The final histopathological results of the 160 STS patients were shown in Table 1. The FNCLCC system assigns a score to the tumor based on its mitotic index, differentiation, and amount of necrosis, and the tumor grade was calculated by adding these three scores. According to their FNCLCC tumor grade, the patients were divided into two groups: low-grade (N = 82) and high-grade (N = 78). The workflow was shown in Fig. 1.

Table 1 The pathologic data of the 160 STS patients
Full size table
Fig. 1
figure 1

(Top) Flow chart of patient enrollment. (Bottom) Work flow of the radiomics implementation

Full size image

MRI Acquisition

All 160 patients underwent FS-T2WI with Siemens Verio3.0 T, Siemens Aera 1.5T (Siemens Medical AG, Erlangen, Germany), or Philips Achia1.5t, Philips Achieva3.0T (Philips Medical Systems, Best, The Netherlands), with adapted position and coils depending on tumor size and location. The scan parameters listed below were used: TR: 2640–5000 ms; TE:30-102ms; slice gap:1 mm; slice thickness:3-4 mm; matrix:320 × 320; The field of view ranges from 200 × 200mm2 to 400 × 400mm2.

Image segmentation and extraction

For image segmentation, all FS T2WI sequence images from patients were uploaded into 3D slicer (version 4.10.2, https://www.slicer.org/, Accessed 8 June 2023). In FS-T2WI images, tumor mass volume (TMV) VOIs were delineated within the margins of tumor masses, encompassing necrotic, cystic change, and hemorrhagic areas but omitting peritumoral edema. The TMV VOIs were then used as a template to construct the corresponding peritumoral tumor volume (PTV) VOIs. The PTV VOIs were generated automatically by uniformly dilating the tumor’s boundary by 10 mm in three dimensions, and adjacent air and bone were manually removed (Fig. 2). The segmentation process was independently performed by two readers (Reader1 and Reader 2) with more than five years of experience, blinded to clinical information and histopathological results. Reader 1 segmented 40 random cases to assess intra-observer reliability two weeks later. Additionally, Reader 2 completed the same 40 random cases to assess inter-observer reliability. Intra- and inter-class Dice coefficients were calculated to assess the stability of delineated VOIs. Features extracted from VOIs with ICCs greater than 0.75 were retained for subsequent investigation.

Fig. 2
figure 2

Example of delineated ROI on FS-T2WI mapping. A 43-years-old woman with pleomorphic sarcoma. A The TMV region is marked in green. B The PTV region is marked in red, and the air region beyond the human tissue has been removed. 1037 feature values were extracted from each of the two disjoint regions through the Slicer-radiomics extension package of 3Dslicer (Feature types, and extraction methods are included in Supplementary Material 1–2)

Full size image

Preprocessing procedures were used to reduce the bias of the features and to counteract the intensity inhomogeneity caused by different imaging protocols before radiomics feature extraction. All VOIs were normalized and resampled to the same resolution (1 mm×1 mm×1 mm) to eliminate data heterogeneity. The limitation of dynamics to µ ± 3σ (µ gray level mean, σ standard deviation) was used to minimize the influence of contrast and brightness variation [21].

Radiomics features were extracted via the Slicer-Radiomics extension in 3D Slicer which enables processing and extraction of radiomic features from medical image data using a large panel of engineered hard-coded feature algorithms by accessing PyRadiomics (https://github.com/AIM-Harvard/pyradiomics, Accessed 8 June 2023) [22]. The detailed operation of extracting features is shown in the Supplementary material (M1-2) 0.1037 radiomics features were extracted from each VOI of TMV and PTV, including first-order statistics(first-order), shape-based(3D) features, shape-based (2D) features, grey-level cooccurrence matrix (GLCM), grey-level run length matrix (GLRLM), grey-level size zone matrix (GLSZM), neighboring grey tone difference matrix (NGTDM), grey-level dependence matrix (GLDM), and wavelet decomposition features. Before further analysis, all the extracted radiomics features of TMV, and PTV were normalized by Z score transformation [23] and ComBat compensation [24] to eliminate the differences in the value scales of the data and remove the batch effects derived from multiple sources of variability caused by different scanners and protocols.

Feature selection

Feature selection was conducted using python software (version 3.8.8, https://www.python.org/, Accessed 8 June 2023), which is mainly implemented by calling the scikit-learn library, a widely used Python library for machine learning and data science [25]. In this step “levene,” “ttest” and “LassoCV” function will be used to select features. A two-step feature selection methodology was performed for the training - validation cohort. Firstly, A t-test was used to filter out the features that were significantly different between the low-grade and high-grade groups. Secondly, the least absolute shrinkage and selection operator (LASSO) method was applied to select the most powerful features in the training - validation cohort and selected non-zero coefficients based on 10 cross-validation. All codes and additional details can be found online (https://github.com/mystic1602/radiomics, Accessed 8 June 2023).

Model construction, rad-score building

To assess the feasibility and promise of the FS-T2WI-based peritumoral radiomics signature for detecting low- and high-grade STSs, the following 3 types of radiomic signatures were extracted: (1) radiomics signatures from TMV features; (2) radiomics signatures from PTV features; (3) radiomics signatures from the merged features of TMV and PTV (TM-PTV).TMV and PTV radiomics signatures were created utilizing the same approach described in the “Radiomics features extraction” and “Feature selection” subsections. TM-PTV features were created by combining TVM and PTV features, and the statistically significant features were chosen using the approach described in the “Feature selection” part.

Prediction models of 3 types of radiomics signatures were created using logistic regression, and three types of radiomics signatures were fed into the the GridSearchCV to establish an ideal parameter configuration [26]. In the external test cohort, their predictive performance was assessed utilizing the area under the curve (AUC) of receiver operating characteristic (ROC) curve analysis. The AUC of each model was evaluated first, and the best model was picked for further investigation. The Rad-score was then calculated using a LASSO logistic regression model based on the best type of radiomics signature.

Development and validation of Nomogram

Univariate and multivariate logistic regression analyses were used to select clinical features and the Rad-score, and a nomogram was constructed based on the independent risk factors in the multivariate study. The model’s discriminative capacity was evaluated using Harrell’s concordance (C-index) with confidence intervals of 95% for both cohorts. The calibration curve was plotted to investigate the model’s predictive accuracy. To assess clinical usefulness, decision curve analysis (DCA) was used to calculate the net benefit of the nomogram model in training and validation groups.

Statistical analysis

Statistical analyses were performed by GraphPad Prism (version 9.4, https://www.graphpad.com/, Accessed 8 June 2023), and R software (version 3.6.2, http://www.Rproject.org, Accessed 8 June2023). When comparing clinical data, The t-test or Mann-Whitney U test was used for continuous variables and Fisher’s exact test for categorical variables. The plot nomograms and calibration curves using the “RMS” software package and the DCA curve were drawn using the “RMDA” software package. Two-sided P < 0.05 was considered statistically significant for all tests.

Results

Clinical data: patient and MRI morphological characteristics

One hundred sixty patients were recruited for this study. In the training cohort, there were 51 high-grade (grade II-III) patients and 60 low-grade (grade I) patients. In the validation cohort, there were 27 high-grade (grade II-III) patients and 22 low-grade (grade I) patients. The detailed clinical and MRI morphologic characteristics of the study population were listed in Table 2. There was no significant difference in age, sex, tumor location, peritumoral hyperintensity, size, and margin between the two cohorts. There were significant differences in signal uniformity between the two cohorts (p < 0.05).

Table 2 Demographic data of patients in the training and validation cohorts
Full size table

Feature selection

The intra-class Dice coefficient of TMV was 0.908 ± 0.032, and the inter- was 0.876 ± 0.057; the intra-class Dice coefficient of PTV was 0.842 ± 0.060, and the inter- was 0.816 ± 0.050 (Supplementary material F1). Features for TMV and PTV were reduced to 845 and 688 respectively after excluding features with ICC less than 0.75. The reserved features were analyzed using a t-test (p < 0.05) to identify features with significant between-group differences between low- and high-grade STSs. Following statistical analysis-based feature selection, three radiomics feature subsets were obtained: (1) 72 significant TMV features, (2) 322 significant PTV features, and (3) 394 significant TMV and PTV features (TM-PTV). Finally, LASSO was used to pick all of the significant features in each feature subset, 8, 14, and 18 discriminative features were chosen, respectively, to create a peritumoral radiomics signature for malignancy grading from TMV, PTV, and TM-PTV. The details of these features were shown in Supplementary material M3-4 and F2.

Performance of radiomics signatures

The AUC of the TM-PTV model was higher than that of the TMV model or PTV model in both the training and validation cohorts. All of the models in the training group had similar performance. However, in the validation cohort, the sensitive, specificity, accuracy, PPV and NPV of the TMV-PTV model is higher than other models (Fig. 3; Table 3).

Fig. 3
figure 3

The ROCs and radar chart performance of Model_TMV, Model_PTV, and Model_TM-PTV in training (A, C) and validation (B, D) cohorts, respectively. AUC: area under the ROC curve. NPV: negative predictive value; PPV: positive predictive value

Full size image
Table 3 Performance of LASSO models in train and validation cohort
Full size table

Rad-score building, identification of Independent risk factors

The TM-PTV model displayed better performance than the other two models. As a result, a LASSO logistic regression model was used to calculate the radiomics score (Rad-score) based on TM-PTV. Supplementary material M4 and F3 illustrate the details of the Rad-score computation formula. A logistic regression analysis with backward stepwise selection identified the Rad-score and peritumoral hyperintensity as independent predictors (Table 4), which were then used to create a personalized prediction nomogram.

Table 4 Result of Univariate Logistic Regression and Multivariable Logistic Regression
Full size table

Development and validation of nomogram

A model containing the above independent predictors was shown in a nomogram (Fig. 4). The model showed a favorable C-index of 0.948 (95% CI, 0.906 to 0.990) in the training cohort and0.921 (95% CI, 0.840 to 0.995) in the validation cohort. Calibration curves of radiomics nomograms used to predict the histopathological grade of STSs showed good agreement between predicted and observed outcomes for both cohorts (Fig. 5). The Hosmer-Lemeshow test showed no statistically significant difference between the calibration curve and the ideal curve (training cohort: χ2 = 8.275, p = 0.407, validation cohort: χ2 = 9.790, p = 0.280), indicating no deviation from a perfect fit.

Fig. 4
figure 4

Nomogram used to distinguish between low- and high-grade levels in soft tissue sarcomas. Development of the radiomics-clinical nomogram, which includes the Rad-score within the second row, and “Peritumoral hyperintensity” within the third row are summed to give the “total points,” which are marked on the “Total points” row. The probability of high grad is read off from the scale in the last row by vertically drawing a line from the total points

Full size image
Fig. 5
figure 5

Calibration curve of the nomogram in the training cohort (A) and validation cohort (B). The blue line represents the perfect prediction of the ideal model, and the red line represents the performance of the model. The Grey line represents corrected predictive performance. The closer the red line and the blue line are, it means that the predicted results are in good agreement with the actual results, and the prediction ability is better

Full size image

The DCA of the TM-PTV model (based on Rad-score) and nomogram model is shown in Fig. 6. In the training cohort, using the nomogram to predict histopathological grade added more benefits than using the TM-PTV model for threshold probabilities from 4 to 90%. In the validation cohort, using the nomogram to predict histopathological grade added more benefit than using the TM-PTV model for threshold probabilities ranging from 8 to 64% and greater than 75%.

Fig. 6
figure 6

Decision curve analysis (DCA). The y axis represents the net benefit, which was determined by calculating the difference between the expected benefit and the expected harm associated with each proposed model [net benefit = true-positive rate (TPR) – (false-positive rate (FPR)× weighting factor), where the weighting factor = threshold probability/ (1-threshold probability)]. The gray line represents the assumption that all tumors were histopathological high-grade (the treat-all scheme). The black line represents the assumption that all tumors were histopathological low-grade expression (the treat-none scheme). A DCA in the training cohort. For threshold probabilities from 4 to 90%, using the nomogram to predict the histopathological grade added more benefit than using the radiomics model. B DCA in the validation cohort. for threshold probabilities ranging from 8–64% and greater than 75%, using the nomogram to predict the histopathological grade added more benefit than using the TM-PTV model

Full size image

The DCA revealed that the net benefit of the nomogram was higher than that of the TM-PTV model, suggesting that the predicting strategy based on nomogram has better clinical utility.

Discussion

A radiomics nomogram was developed based on the combination of intra-peritumoral radiomic features and peritumoral hyperintensity, which successfully differentiated low-grade and high-grade STSs, and its recognition performance was higher than the TM-PTV model, and showed good calibration in both groups, suggesting that it may be a promising tool for clinical strategy. To the best of our knowledge, this is the first study to develop a model based on intratumoral and peritumoral radiomics to predict the histopathological grade of STSs. The present study found peritumoral hyperintensity and Rad-score were independent risk factors for predicting histopathological grade. Peritumoral hyperintensity on FS-T2WI may be caused by surrounding tissue edema, inflammatory cell infiltration, angiogenesis, and other changes [27]. It has been widely recognized as a significantly poor prognostic factor for STS [28].

Previous studies have reported on an MRI-based radiomics model for predicting STS histopathological grade [6, 8]. Zhang et al. [6] developed an MRI-based radiological model with an AUC of 0.92 in predicting histopathological grade. Yan et al. [8] developed an MRI-based radiomics nomogram for predicting the grade and demonstrated good calibration and good clinical utility. These two available studies were all based on TMV models. However, Significant clinical evidence suggests the heterogeneity of STSs is not limited to tumor margins but also involves peritumoral regions [13]. White et al. have demonstrated in a population of patients undergoing surgery, satellite-like single or clustered tumor cells that are not visible on imaging can be found beyond the tumor margins [12].In our study, the TMV, PTV, and TM-PTV models demonstrated outstanding discrimination in both the training and validation cohorts. Compared with the TMV model, the combination of the TMV and PTV(PT-TMV) model significantly improved the AUC, accuracy, sensitivity, specificity, PPV, and NPV of the predictive model. This result indicated that the combination of peritumoral features provides more information about the tumor microenvironment, which can reflect the biological behavior of the tumor better.

Our study used a 10 mm extension from the lesion, which is based on the current standard for STSs surgical margins. To protect critical neurovascular structures or bones, a margin aiming at 10 mm is the minimum appropriate width to be considered acceptable according to National Comprehensive Cancer Network (NCCN) guidelines [29]. Previous peritumoral radiomics studies of lung cancer [30] and glioma [14] have shown that the closer the peritumoral region is to the intratumoral region, the more information it contains. As the expansion distance increases, more normal soft tissue is incorporated into the region of interest (ROI), resulting in a smaller difference in peritumoral tissue heterogeneity. Radiomics features extracted from a 10 mm peritumoral ROI were most likely to provide important information for predicting histopathological grade for STS. In this study, the radiomics model was constructed based on a single FS-T2WI sequence, which was currently the most widely used sequence for STS radiomics. Peeken’s study found the FS-T2WI-based model showed better reproducibility compared with the T1WI model and the combined T1 and T2WI models [7]. Importantly, the present study focused on the efficacy of radiomics models of the tumor itself and the tumor microenvironment. However, most lesions have unclear borders on T1WI, and image segmentation is difficult, which may affect the division of intratumoral and peritumoral ROIs. In addition, extracting a large number of features from many sequences tends to increase the risk of model overfitting [31].

Radiomics features are sensitive to all the acquisition conditions including MR protocols, scanners, and MR adjustments. To lower bias and variance, We performed several preprocessing methods, In particular, the ComBat compensation,which is easy to apply and is suitable for retrospective information analysis [24], eliminates batch effects due to multi-source variation caused by different scanners and protocols in multi-center radiation analysis while preserving the excellent properties of its texture patterns and has been used in previous reports to improve reproducibility between different centers [32]. Therefore, the key to addressing the heterogeneity of acquisition conditions may lie in adequate preprocessing and consistent scanning parameters.

When dealing with radiomics, the segmentation method has a large impact on the reproducibility and reliability of the radiomics signature [33]. The current image segmentation methods include manual, semiautomatic, and fully automatic. Manual segmentation is prone to interobserver variability, which may hinder the reproducibility of radiomics analyses. Automatic segmentation techniques based on deep learning are the mainstream of the current research field and are currently being considered for clinical trials, showing improvements in image classification prediction and recognition tasks [34]. While results from automatic image segmentation are promising, errors due to contrast blur and biased fields are common and often require manual correction to ensure accuracy [35]. We tried to use U-net neural network for fully automatic segmentation of lesions, but because STS comes from various types of tissues, shapes, and positions are not fixed, the automatic segmentation of lesions was not effective. Therefore, this study still used the traditional manual segmentation method. Although both the intra-observer and inter-observer ICC coefficients were high, this was a hugely time-consuming process, which was a limitation of this study. In the future, in-depth research was expected to break through the limitations of image segmentation. This study has several limitations that need to be addressed. First, our study was retrospective, so despite our strict criteria, there is a potential selection bias. Second, our data came from two institutions using similar but different scanners and protocols. Therefore, the resampling methodology and the combat compensation method were adopted to reduce the difference in image specifications to improve the stability of features and different models. Third, multimodal MR radiomics may have better potential, such as DWI. Finally,this study did not include all histopathological types of STS because our study was retrospective and limited by the number of patients. Therefore, the sample needs to be enlarged in future research to improve the generalization ability of the model.

In conclusion, peritumoral radiomics features can provide complementary information to intratumoral regions to predict the histopathological grade of STS. Such quantitative radiomics prognostic models of STS may potentially be useful for precision medicine.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request. All codes and additional details can be found online (https://github.com/mystic1602/radiomics, Accessed 8 June 2023).

Abbreviations

AUC:

Area Under Curve

DCA:

Decision Curve Analysis

DWI:

Diffusion Weighted Imaging

FNCLCC:

French Federation of Cancer Centers Sarcoma Group

FST2WI:

Fat-Saturated T2-Weighted

ICC:

Intraclass Correlation Coefficient

LASSO:

Least Absolute Shrinkage and Selection Operator

MRI:

Magnetic Resonance Imaging

PTM:

Peritumoral Tumor Volume

ROC:

Receiver Operating Characteristic Curve

ROI:

Region of Interest

STSs:

Soft Tissue Sarcomas

SVM:

Support Vector Machine

TM-PTV:

Combined TMV and PTV

TMV:

Tumor Mass Volume

VOI:

Volume of Interest

References

  1. Coindre JM. Grading of soft tissue sarcomas: review and update. Arch Pathol Lab Med. 2006;130(10):1448–53. https://doi.org/10.5858/2006-130-1448-GOSTSR.

    Article  PubMed  Google Scholar 

  2. Italiano A, Le Cesne A, Mendiboure J, Blay JY, Piperno-Neumann S, Chevreau C, et al. Prognostic factors and impact of adjuvant treatments on local and metastatic relapse of soft-tissue sarcoma patients in the competing risks setting. Cancer. 2014;120(21):3361–9. https://doi.org/10.1002/cncr28885.

    Article  PubMed  Google Scholar 

  3. Schneider N, Strauss DC, Smith MJ, Miah AB, Zaidi S, Benson C, et al. The adequacy of Core Biopsy in the Assessment of smooth muscle Neoplasms of soft tissues: implications for treatment and prognosis. Am J Surg Pathol. 2017;41(7):923–31. https://doi.org/10.1097/PAS.0000000000000867.

    Article  PubMed  Google Scholar 

  4. Wang H, Chen H, Duan S, Hao D, Liu J. Radiomics and Machine Learning with Multiparametric Preoperative MRI May accurately predict the histopathological grades of soft tissue sarcomas. J Magn Reson Imaging. 2020;51(3):791–7. https://doi.org/10.1002/jmri.26901.

    Article  PubMed  Google Scholar 

  5. Fang S, Yang Y, Xu N, Tu Y, Yin Z, Zhang Y, et al. An update in imaging evaluation of histopathological Grade of Soft tissue sarcomas using structural and Quantitative Imaging and Radiomics. J Magn Reson Imaging. 2022;55(5):1357–75. https://doi.org/10.1002/jmri27954.

    Article  PubMed  Google Scholar 

  6. Zhang Y, Zhu Y, Shi X, Tao J, Cui J, Dai Y, et al. Soft tissue sarcomas: preoperative predictive histopathological grading based on Radiomics of MRI. Acad Radiol. 2019;26(9):1262–8. https://doi.org/10.1016/jacra201809025.

    Article  PubMed  Google Scholar 

  7. Peeken JC, Spraker MB, Knebel C, Dapper H, Pfeiffer D, Devecka M, et al. Tumor grading of soft tissue sarcomas using MRI-based radiomics. EBioMedicine. 2019;48:332–40. https://doi.org/10.1016/jebiom201908059.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Yan R, Hao D, Li J, Liu J, Hou F, Chen H, et al. Magnetic resonance imaging-based Radiomics Nomogram for Prediction of the histopathological Grade of Soft tissue sarcomas: a two-Center Study. J Magn Reson Imaging. 2021;53(6):1683–96. https://doi.org/10.1002/jmri27532.

    Article  PubMed  Google Scholar 

  9. Wu T, Dai Y. Tumor microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8. https://doi.org/10.1016/jcanlet201601043.

    Article  CAS  PubMed  Google Scholar 

  10. Hirata E, Sahai E. Tumor Microenvironment and Differential responses to Therapy. Cold Spring Harb Perspect Med. 2017;7(7):a026781. https://doi.org/10.1101/cshperspect.a026781.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Semenza GL, Ruvolo PP. Introduction to Tumor microenvironment regulation of cancer cell survival, Metastasis, inflammation, and immune surveillance. Biochim Biophys Acta. 2016;1863(3):379–81. https://doi.org/10.1016/j.bbamcr.2015.12.015.

    Article  CAS  PubMed  Google Scholar 

  12. White LM, Wunder JS, Bell RS, O’Sullivan B, Catton C, Ferguson P, et al. Histologic assessment of peritumoral edema in soft tissue sarcoma. Int J Radiat Oncol Biol Phys. 2005;61(5):1439–45. https://doi.org/10.1016/jijrobp200408036.

    Article  PubMed  Google Scholar 

  13. Hoefkens F, Dehandschutter C, Somville J, Meijnders P, Van Gestel D. Soft tissue sarcoma of the extremities: pending questions on Surgery and radiotherapy. Radiat Oncol. 2016;11(1):136. https://doi.org/10.1186/s13014-016-0668-9.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Cheng J, Liu J, Yue H, Bai H, Pan Y, Wang J. Prediction of Glioma Grade using Intratumoral and Peritumoral Radiomic features from multiparametric MRI images. IEEE/ACM Trans Comput Biol Bioinform. 2022;19(2):1084–95. https://doi.org/10.1109/TCBB20203033538.

    Article  PubMed  Google Scholar 

  15. Zhou Z, Qian X, Hu J, Ma X, Zhou S, Dai Y, et al. CT-based peritumoral radiomics signatures for malignancy grading of clear cell renal cell carcinoma. Abdom Radiol (NY). 2021;46(6):2690–8. https://doi.org/10.1007/s00261-020-02890-z.

    Article  PubMed  Google Scholar 

  16. Perez-Morales J, Tunali I, Stringfield O, Eschrich SA, Balagurunathan Y, Gillies RJ, et al. Peritumoral and intratumoral radiomic features predict survival outcomes among patients diagnosed in Lung cancer screening. Sci Rep. 2020;10(1):10528. https://doi.org/10.1038/s41598-020-67378-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Das SK, Fang KW, Xu L, Li B, Zhang X, Yang HF. Integrative nomogram of intratumoral, peritumoral, and lymph node radiomic features for prediction of lymph node Metastasis in cT1N0M0 lung adenocarcinomas. Sci Rep. 2021;11(1):10829. https://doi.org/10.1038/s41598-021-90367-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Dou TH, Coroller TP, van Griethuysen JJM, Mak RH, Aerts H. Peritumoral radiomics features predict distant Metastasis in locally advanced NSCLC. PLoS ONE. 2018;13(11): e0206108. https://doi.org/10.1371/journal.pone.0206108.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Braman NM, Etesami M, Prasanna P, Dubchuk C, Gilmore H, Tiwari P, et al. Intratumoral and peritumoral radiomics for the pretreatment prediction of pathological complete response to neoadjuvant chemotherapy based on breast DCE-MRI. Breast Cancer Res. 2017;19(1):57. https://doi.org/10.1186/s13058-017-0846-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhao F, Ahlawat S, Farahani SJ, Weber KL, Montgomery EA, Carrino JA, et al. Can MR imaging be used to predict Tumor grade in soft-tissue sarcoma? Radiology. 2014;272(1):192–201. https://doi.org/10.1148/radiol.14131871.

    Article  PubMed  Google Scholar 

  21. Meyer HJ, Renatus K, Hohn AK, Hamerla G, Schopow N, Fakler J, et al. Texture analysis parameters derived from T1-and T2-weighted magnetic resonance images can reflect Ki67 index in soft tissue sarcoma. Surg Oncol. 2019;30:92–7. https://doi.org/10.1016/j.suronc.2019.06.006.

    Article  PubMed  Google Scholar 

  22. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic phenotype. Cancer Res. 2017;77(21):e104. https://doi.org/10.1158/0008-5472.CAN-17-0339.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zheng Z, Gu Z, Xu F, Maskey N, He Y, Yan Y, et al. Magnetic resonance imaging-based radiomics signature for preoperative prediction of Ki67 expression in Bladder cancer. Cancer Imaging. 2021;21(1):65. https://doi.org/10.1186/s40644-021-00433-3.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Orlhac F, Lecler A, Savatovski J, Goya-Outi J, Nioche C, Charbonneau F, et al. How can we combat multicenter variability in MR radiomics? Validation of a correction procedure. Eur Radiol. 2021;31(4):2272–80. https://doi.org/10.1007/s00330-020-07284-9.

    Article  PubMed  Google Scholar 

  25. Ahn S. Building and analyzing machine learning-based warfarin dose prediction models using scikit-learn. Transl Clin Pharmacol. 2022;30(4):172–81. https://doi.org/10.12793/tcp.2022.30.e22.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Alhakeem ZM, Jebur YM, Henedy SN, Imran H, Bernardo LFA, Hussein HM. Prediction of ecofriendly concrete compressive strength using gradient boosting regression Tree Combined with GridSearchCV Hyperparameter-optimization techniques. Mater (Basel). 2022;15(21):7432. https://doi.org/10.3390/ma15217432.

    Article  CAS  Google Scholar 

  27. Liu QY, Li HG, Chen JY, Liang BL. [Correlation of MRI features to histopathologic grade of soft tissue sarcoma]. Ai Zheng. 2008;27(8):856–60.

    PubMed  Google Scholar 

  28. Gruber L, Gruber H, Luger AK, Glodny B, Henninger B, Loizides A. Diagnostic hierarchy of radiological features in soft tissue tumours and proposition of a simple diagnostic algorithm to estimate malignant potential of an unknown mass. Eur J Radiol. 2017;95:102–10. https://doi.org/10.1016/j.ejrad.2017.07.020.

    Article  PubMed  Google Scholar 

  29. von Mehren M, Randall RL, Benjamin RS, Boles S, Bui MM, Casper ES, et al. Soft tissue sarcoma, version 2.2014. J Natl Compr Canc Netw. 2014;12(4):473–83. https://doi.org/10.6004/jnccn.2014.0053.

    Article  Google Scholar 

  30. Tang X, Huang H, Du P, Wang L, Yin H, Xu X. Intratumoral and peritumoral CT-based radiomics strategy reveals distinct subtypes of non-small-cell Lung cancer. J Cancer Res Clin Oncol. 2022;148(9):2247–60. https://doi.org/10.1007/s00432-022-04015-z.

    Article  PubMed  Google Scholar 

  31. Fusco R, Granata V, Grazzini G, Pradella S, Borgheresi A, Bruno A, et al. Radiomics in medical imaging: pitfalls and challenges in clinical management. Jpn J Radiol. 2022;40(9):919–29. https://doi.org/10.1007/s11604-022-01271-4.

    Article  PubMed  Google Scholar 

  32. Fortin JP, Cullen N, Sheline YI, Taylor WD, Aselcioglu I, Cook PA, et al. Harmonization of cortical thickness measurements across scanners and sites. Neuroimage. 2018;167:104–20. https://doi.org/10.1016/j.neuroimage.2017.11.024.

    Article  PubMed  Google Scholar 

  33. Heye T, Merkle EM, Reiner CS, Davenport MS, Horvath JJ, Feuerlein S, et al. Reproducibility of dynamic contrast-enhanced MR imaging. Part II. Comparison of intra- and interobserver variability with manual region of interest placement versus semiautomatic lesion segmentation and histogram analysis. Radiology. 2013;266(3):812–21. https://doi.org/10.1148/radiol.12120255.

    Article  PubMed  Google Scholar 

  34. Zhu HT, Zhang XY, Shi YJ, Li XT, Sun YS. Automatic segmentation of rectal Tumor on diffusion-weighted images by deep learning with U-Net. J Appl Clin Med Phys. 2021;22(9):324–31. https://doi.org/10.1002/acm2.13381.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Clendenen TV, Zeleniuch-Jacquotte A, Moy L, Pike MC, Rusinek H, Kim S. Comparison of 3-point Dixon imaging and fuzzy C-means clustering methods for breast density measurement. J Magn Reson Imaging. 2013;38(2):474–81. https://doi.org/10.1002/jmri.24002.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the plastic surgery team of Sichuan Provincial People’s Hospital and the First People’s Hospital of Yibin for collecting and providing the MRI images.

Funding

The authors state that this work has not received any funding.

Author information

Author notes

  1. Liyuan Zhang and Yang Yang contributed equally to this work and should be considered co-first authors.

Authors and Affiliations

  1. Department of Plastic Surgery, Sichuan Provincial People’s Hospital, School of Medicine,University of Electronic Science and Technology of China, Chengdu, 610000, People’s Republic of China

    Liyuan Zhang, Junnan Deng, Zhen Cai & Wei Cui

  2. Department of Radiology, Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu, 610000, People’s Republic of China

    Yang Yang & Mingyue Tang

  3. Department of Plastic Surgery, The First People’s Hospital of Yibin, Yibin, 644000, People’s Republic of China

    Ting Wang

  4. Sichuan College of Traditional Chinese Medicine, Mianyang, 621000, People’s Republic of China

    Xi Chen

Authors
  1. Liyuan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yang Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ting Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xi Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mingyue Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Junnan Deng
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Zhen Cai
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Wei Cui
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Wei Cui and Zhen Cai participated in the design of the study and the initial syntheses of results. Liyuan Zhang and Yang Yang prepared figures and wrote the manuscript text. Xi Chen,Mingyue Tang and Junnan Deng made statistical analysis. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zhen Cai or Wei Cui.

Ethics declarations

Ethics approval and consent to participate

The Ethics Committee of Sichuan Provincial People’s Hospital and the First People’s Hospital of Yibin approved this study.

Consent for publication

Not applicable.

Competing interests

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Material 1.

Detailed information about extracted radiomics features. Material 2. Image Types were used to extract features from. Material 3. Features screened by the LASSO method. Material 4. Calculation formula of Rad-score. Figure 1. The Dice coefficient of TMV and PTV. Figure 2. Detailed information about feature extract by python. Figure 3. Rad-score for each patient in each cohort.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, L., Yang, Y., Wang, T. et al. Intratumoral and peritumoral MRI-based radiomics prediction of histopathological grade in soft tissue sarcomas: a two-center study. Cancer Imaging 23, 103 (2023). https://doi.org/10.1186/s40644-023-00622-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40644-023-00622-2

Keywords

Cancer Imaging

ISSN: 1470-7330

Contact us