High-Grade Soft-Tissue Sarcomas: Can Optimizing Dynamic Contrast-Enhanced MRI Postprocessing Improve Prognostic Radiomics Models?

Integration of pre-treatment computational radiomics, deep radiomics, and transcriptomics enhances soft-tissue sarcoma patient prognosis

Our objective was to capture subgroups of soft-tissue sarcoma (STS) using handcraft and deep radiomics approaches to understand their relationship with histopathology, gene-expression profiles, and metastatic relapse-free survival (MFS). We included all consecutive adults with newly diagnosed locally advanced STS (N = 225, 120 men, median age: 62 years) managed at our sarcoma reference center between 2008 and 2020, with contrast-enhanced baseline MRI. After MRI postprocessing, segmentation, and reproducibility assessment, 175 handcrafted radiomics features (h-RFs) were calculated. Convolutional autoencoder neural network (CAE) and half-supervised CAE (HSCAE) were trained in repeated cross-validation on representative contrast-enhanced slices to extract 1024 deep radiomics features (d-RFs). Gene-expression levels were calculated following RNA sequencing (RNAseq) of 110 untreated samples from the same cohort. Unsupervised classifications based on h-RFs, CAE, HSCAE, and RNAseq were built. The h-RFs, CAE, and HSCAE grouping were not associated with the transcriptomics groups but with prognostic radiological features known to correlate with lower survivals and higher grade and SARCULATOR groups (a validated prognostic clinical-histological nomogram). HSCAE and h-RF groups were also associated with MFS in multivariable Cox regressions. Combining HSCAE and transcriptomics groups significantly improved the prognostic performances compared to each group alone, according to the concordance index. The combined radiomic-transcriptomic group with worse MFS was characterized by the up-regulation of 707 genes and 292 genesets related to inflammation, hypoxia, apoptosis, and cell differentiation. Overall, subgroups of STS identified on pre-treatment MRI using handcrafted and deep radiomics were associated with meaningful clinical, histological, and radiological characteristics, and could strengthen the prognostic value of transcriptomics signatures.

The impact of radiomics in the management of soft tissue sarcoma

INTRODUCTION

Soft tissue sarcomas (STSs) are rare malignancies. Pre-therapeutic tumour grading and assessment are crucial in making treatment decisions. Radiomics is a high-throughput method for analysing imaging data, providing quantitative information beyond expert assessment. This review highlights the role of radiomic texture analysis in STSs evaluation.

MATERIALS AND METHODS

We conducted a systematic review according to the Systematic Reviews and Meta-Analyses (PRISMA) guidelines. A comprehensive search was conducted in PubMed/MEDLINE and Scopus using the search terms: 'radiomics [All Fields] AND ("soft tissue sarcoma" [All Fields] OR "soft tissue sarcomas" [All Fields])'. Only original articles, referring to humans, were included.

RESULTS

A preliminary search conducted on PubMed/MEDLINE and Scopus provided 74 and 93 studies respectively. Based on the previously described criteria, 49 papers were selected, with a publication range from July 2015 to June 2023. The main domains of interest were risk stratification, histological grading prediction, technical feasibility/reproductive aspects, treatment response.

CONCLUSIONS

With an increasing interest over the last years, the use of radiomics appears to have potential for assessing STSs from initial diagnosis to predicting treatment response. However, additional and extensive research is necessary to validate the effectiveness of radiomics parameters and to integrate them into a comprehensive decision support system.

Radiomics and artificial intelligence for soft-tissue sarcomas: Current status and perspectives

This article proposes a summary of the current status of the research regarding the use of radiomics and artificial intelligence to improve the radiological assessment of patients with soft tissue sarcomas (STS), a heterogeneous group of rare and ubiquitous mesenchymal malignancies. After a first part explaining the principle of radiomics approaches, from raw image post-processing to extraction of radiomics features mined with unsupervised and supervised machine-learning algorithms, and the current research involving deep learning algorithms in STS, especially convolutional neural networks, this review details their main research developments since the formalisation of 'radiomics' in oncologic imaging in 2010. This review focuses on CT and MRI and does not involve ultrasonography. Radiomics and deep radiomics have been successfully applied to develop predictive models to discriminate between benign soft-tissue tumors and STS, to predict the histologic grade (i.e., the most important prognostic marker of STS), the response to neoadjuvant chemotherapy and/or radiotherapy, and the patients' survivals and probability for presenting distant metastases. The main findings, limitations and expectations are discussed for each of these outcomes. Overall, after a first decade of publications emphasizing the potential of radiomics through retrospective proof-of-concept studies, almost all positive but with heterogeneous and often non-replicable methods, radiomics is now at a turning point in order to provide robust demonstrations of its clinical impact through open-science, independent databases, and application of good and standardized practices in radiomics such as those provided by the Image Biomarker Standardization Initiative, without forgetting innovative research paths involving other '-omics' data to better understand the relationships between imaging of STS, gene-expression profiles and tumor microenvironment.

Role of Imaging in Initial Prognostication of Locally Advanced Soft Tissue Sarcomas

BACKGROUND

Although imaging is central in the initial staging of patients with soft tissue sarcomas (STS), it remains underused and few radiological features are currently used in practice for prognostication and to help guide the best therapeutic strategy. Yet, several prognostic qualitative and quantitative characteristics from magnetic resonance imaging (MRI) and positron emission tomography (PET) have been identified over these last decades.

OBJECTIVE

After an overview of the current validated prognostic features based on baseline imaging and their integration into prognostic tools, such as nomograms used by clinicians, the aim of this review is to summarize more complex and innovative MRI, PET, and radiomics features, and to highlight their role to predict indirectly (through histologic grade) or directly the patients' outcomes.

Distinct patterns of the natural evolution of soft tissue sarcomas on pre-treatment MRIs captured with delta-radiomics correlate with gene expression profiles

OBJECTIVES

Radiomics of soft tissue sarcomas (STS) is assumed to correlate with histologic and molecular tumor features, but radiogenomics analyses are lacking. Our aim was to identify if distinct patterns of natural evolution of STS obtained from consecutive pre-treatment MRIs are associated with differential gene expression (DGE) profiling in a pathway analysis.

METHODS

All patients with newly diagnosed STS treated in a curative intent in our sarcoma reference center between 2008 and 2019 and with two available pre-treatment contrast-enhanced MRIs were included in this retrospective study. Radiomics features (RFs) were extracted from fat-sat contrast-enhanced T1-weighted imaging. Log ratio and relative change in RFs were calculated and used to determine grouping of samples based on a consensus hierarchical clustering. DGE and oncogenesis pathway analysis were performed in the delta-radiomics groups identified in order to detect associations between delta-radiomics patterns and transcriptomics features of STS. Secondarily, the prognostic value of the delta-radiomics groups was investigated.

RESULTS

Sixty-three patients were included (median age: 63 years, interquartile range: 52.5-70). The consensus clustering identified 3 reliable delta-radiomics patient groups (A, B, and C). On imaging, group B patients were characterized by increase in tumor heterogeneity, necrotic signal, infiltrative margins, peritumoral edema, and peritumoral enhancement before the treatment start (p value range: 0.0019-0.0244), and, molecularly, by downregulation of natural killer cell-mediated cytotoxicity genes and upregulation of Hedgehog and Hippo signaling pathways. Group A patients were characterized by morphological stability of pre-treatment MRI traits and no local relapse (log-rank p = 0.0277).

CONCLUSIONS

This study highlights radiomics and transcriptomics convergence in STS. Proliferation and immune response inhibition were hyper-activated in the STS that were the most evolving on consecutive imaging.

KEY POINTS

• Three consensual and stable delta-radiomics clusters were identified and captured the natural patterns of morphological evolution of STS on pre-treatment MRIs. • These 3 patterns were explainable and correlated with different well-known semantic radiological features with an ascending gradient of pejorative characteristics from the A group to C group to B group. • Gene expression profiling stressed distinct patterns of up/downregulated oncogenetic pathways in STS from B group in keeping with its most aggressive radiological evolution.

Impact of Dynamic Contrast Enhanced and Diffusion-Weighted MR Imaging on Detection of Early Local Recurrence of Soft Tissue Sarcoma

BACKGROUND

Diagnosis of residual or recurrent tumor in soft-tissue sarcomas (STS) is a differential diagnostic challenge since post-therapeutic changes impede diagnosis.

PURPOSE

To evaluate the diagnostic accuracy of quantitative dynamic contrast enhanced (DCE)-MRI and diffusion-weighted imaging (DWI) to detect local recurrence of STS of the limb.

STUDY TYPE

Prospective.

POPULATION

A totalof 64 consecutive patients with primary STS of the limbs were prospectively included 3-6 months after surgery between January 2016 and July 2021.

FIELD STRENGTH/SEQUENCE

A 1.5 T; axial DWI echo-planar imaging sequences and DCE-MRI using a 3D T1-weighted spoiled gradient-echo sequence.

ASSESSMENT

The quantitative DCE-MRI parameters relative plasma flow (rPF) and relative mean transit time (rMTT) were calculated and ADC mapping was used to quantify diffusion restriction. Regions of interest of tumor growth and postoperative changes were drawn in consensus by two experts for diffusion and perfusion analysis. An additional morphological assessment was done by three independent and blinded radiologists.

STATISTICAL TEST

Unpaired t-test, ROC-analysis, and a logistic regression model were applied. Interobserver reliability was calculated using Fleiss kappa statistics. A P value of 0.05 was considered statistically significant.

RESULTS

A total of 11 patients turned out to have local recurrence. rPF was significantly higher in cases of local recurrence when compared to cases without local recurrence (61.1-4.5) while rMTT was slightly and significantly lower in local recurrence. ROC-analysis showed an area under the curve (AUC) of 0.95 (SEM ± 0.05) for rPF while a three-factor multivariate logistic regression model showed a high diagnostic accuracy of rPF (R²  = 0.71). Compared with morphological assessment, rPF had a distinct higher specificity and true positive value in detection of LR.

DATA CONCLUSION

DCE-MRI is a promising additional method to differentiate local recurrence from benign postoperative changes in STS of the limb. Especially specificity in detection of LR is increased compared to morphological assessment.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

The diagnosis, classification, and treatment of sarcoma in this era of artificial intelligence and immunotherapy

Soft-tissue sarcomas (STS) represent a group of rare and heterogeneous tumors associated with several challenges, including incorrect or late diagnosis, the lack of clinical expertise, and limited therapeutic options. Digital pathology and radiomics represent transformative technologies that appear promising for improving the accuracy of cancer diagnosis, characterization and monitoring. Herein, we review the potential role of the application of digital pathology and radiomics in managing patients with STS. We have particularly described the main results and the limits of the studies using radiomics to refine diagnosis or predict the outcome of patients with soft-tissue sarcomas. We also discussed the current limitation of implementing radiomics in routine settings. Standard management approaches for STS have not improved since the early 1970s. Immunotherapy has revolutionized cancer treatment; nonetheless, immuno-oncology agents have not yet been approved for patients with STS. However, several lines of evidence indicate that immunotherapy may represent an efficient therapeutic strategy for this group of diseases. Thus, we emphasized the remarkable potential of immunotherapy in sarcoma treatment by focusing on recent data regarding the immune landscape of these tumors. We have particularly emphasized the fact that the development of immunotherapy for sarcomas is not an aspect of histology (except for alveolar soft-part sarcoma) but rather that of the tumor microenvironment. Future studies investigating immunotherapy strategies in sarcomas should incorporate at least the presence of tertiary lymphoid structures as a stratification factor in their design, besides including a strong translational program that will allow for a better understanding of the determinants involved in sensitivity and treatment resistance to immune-oncology agents.

Tumor cellularity beyond the visible in soft tissue sarcomas: Results of an ADC-based, single center, and preliminary radiomics study

Purpose

Soft tissue sarcomas represent approximately 1% of all malignancies, and diagnostic radiology plays a significant role in the overall management of this rare group of tumors. Recently, quantitative imaging and, in particular, radiomics demonstrated to provide significant novel information, for instance, in terms of prognosis and grading. The aim of this study was to evaluate the prognostic role of radiomic variables extracted from apparent diffusion coefficient (ADC) maps collected at diagnosis in patients with soft tissue sarcomas in terms of overall survival and metastatic spread as well as to assess the relationship between radiomics and the tumor grade.

Methods

Patients with histologically proven soft tissue sarcomas treated in our tertiary center from 2016 to 2019 who underwent an Magnetic Resonance (MR) scan at diagnosis including diffusion-weighted imaging were included in this retrospective institution review board–approved study. Each primary lesion was segmented using the b50 images; the volumetric region of interest was then applied on the ADC map. A total of 33 radiomic features were extracted, and highly correlating features were selected by factor analysis. In the case of feature/s showing statistically significant results, the diagnostic accuracy was computed. The Spearman correlation coefficient was used to evaluate the relationship between the tumor grade and radiomic features selected by factor analysis. All analyses were performed applying p<0.05 as a significant level.

Results

A total of 36 patients matched the inclusion criteria (15 women; mean age 58.9 ± 15 years old). The most frequent histotype was myxofibrosarcoma (16.6%), and most of the patients were affected by high-grade lesions (77.7%). Seven patients had pulmonary metastases, and, altogether, eight were deceased. Only the feature Imc1 turned out to be a predictor of metastatic spread (p=0.045 after Bonferroni correction) with 76.7% accuracy. The value -0.16 showed 73.3% sensitivity and 71.4% specificity, and patients with metastases showed lower values (mean Imc1 of metastatic patients -0.31). None of the examined variables was a predictor of the overall outcome (p>0.05, each). A moderate statistically significant correlation emerged only between Imc1 and the tumor grade (r=0.457, p=0.005).

Conclusions

In conclusion, the radiomic feature Imc1 acts as a predictor of metastatic spread in patients with soft tissue sarcomas and correlates with the tumor grade.

Soft Tissue Sarcomas: The Role of Quantitative MRI in Treatment Response Evaluation

BACKGROUND

Although curative surgery remains the cornerstone of the therapeutic strategy in patients with soft tissue sarcomas (STS), neoadjuvant radiotherapy and chemotherapy (NART and NACT, respectively) are increasingly used to improve operability, surgical margins and patient outcome. The best imaging modality for locoregional assessment of STS is MRI but these tumors are mostly evaluated in a qualitative manner.

OBJECTIVE

After an overview of the current standard of care regarding treatment for patients with locally advanced STS, this review aims to summarize the principles and limitations of (i) the current methods used to evaluate response to neoadjuvant treatment in clinical practice and clinical trials in STS (RECIST 1.1 and modified Choi criteria), (ii) quantitative MRI sequences (i.e., diffusion weighted imaging and dynamic contrast enhanced MRI), and (iii) texture analyses and (delta-) radiomics.

Multi-parametric MRI-based radiomics for the diagnosis of malignant soft-tissue tumor

PURPOSE

To develop and validate a multiparametric magnetic resonance imaging-based radiomics nomogram for differentiating malignant and benign soft-tissue tumors.

METHODS

A total of 91 patients with pathologically confirmed soft-tissue tumors were enrolled between January 2017 and October 2020. Forty-eight patients were consecutively enrolled between November 2020 and March 2022, as a time-independent cohort. All patients underwent contrast-enhanced T1-weighted and T2-weighted fat-suppression magnetic resonance scans at 3.0 T. Radiomics features were extracted and selected from the two modalities to develop the radiomics signature. Significant clinical/morphological characteristics were identified using a multivariate logistic regression analysis. The least absolute shrinkage and selection operator regression were applied to identify discriminative features. A clinical-radiomics nomogram was constructed based on clinical/morphological characteristics and radiomics features. Finally, the performance of the nomogram was validated using the receiver operating characteristic and decision curve analysis (DCA).

RESULTS

Six features were selected to establish the combined RS. Size, margin, and peritumoral edema were identified as the most important clinical and morphological factors, respectively. The radiomics signature outperformed the clinical model in terms of AUC and sensitivity. The nomogram integrating the combined RS, size, margin, and peritumoral edema achieved favorable predictive efficacy, generating AUCs of 0.954 (95% confidence interval [CI]: 0.907-1.000, Sen = 0.861, Spe = 0.917), 0.962 (95% CI: 0.901-1.000, Sen = 0.944, Spe = 0.923), and 0.935 (95% CI: 0.871-0.998, Sen = 0.815, Spe = 0.952) in the training (n = 60), validation (n = 31) and time-independent (n = 48) cohorts, respectively. The DCA curve indicated good clinical usefulness of the nomogram.

CONCLUSIONS

Our study demonstrated the clinical potential of multiparametric MRI-based radiomics in distinguishing malignant from benign soft-tissue tumors, which can be considered as a noninvasive tool for individual treatment management.

Radiomic features as biomarkers of soft tissue paediatric sarcomas: preliminary results of a PET/MR study

BACKGROUND

Pediatric soft tissue sarcomas are rare tumors with rhabdomyosarcoma being the most frequent histotype. Diagnostic imaging plays a significant role in the evaluation of this type of tumors. Thus, aim of this study was to assess the prognostic and diagnostic value of radiomic features extracted from axial T2w images of the primary lesion in children with soft tissue sarcomas examined by PET/MR for staging.

METHODS

Using an open source software, each lesion was segmented and 33 radiomic features then extracted. Factor and logistic regression analyses were applied to select highly correlating features and evaluate their prognostic role, respectively. Differences in radiomic, demographics, metabolic, and laboratory variables according to tumor grade and histotype were investigated by the Students' and Chi-square tests. In case of differences the diagnostic value of the variable/s was assessed by receiver operating curves.

RESULTS

Eighteen children (11 female; mean age 7.8 ± 4.6-year-old) matched the inclusion criteria. The factor analysis allowed the selection of five highly correlating features which, according to regression analysis, did not influence the outcome (p > 0.05, each). The feature lmc1 was significantly higher in low grade lesions (p = 0.045) and showed 70.4% accuracy in classifying high grade tumors while the feature variance was significantly lower in rhabdomyosarcomas (p = 0.008) and showed 83.3% accuracy for this histotype.

CONCLUSIONS

In conclusion, our preliminary results suggest that specific radiomic features may act as biomarkers of pediatric soft tissue sarcoma grade and histotype.

MRI-based delta-radiomics predicts pathologic complete response in high-grade soft-tissue sarcoma patients treated with neoadjuvant therapy

PURPOSE

In high-grade soft-tissue sarcomas (STS) the standard of care encompasses multimodal therapy regimens. While there is a growing body of evidence for prognostic pretreatment radiomic models, we hypothesized that temporal changes in radiomic features following neoadjuvant treatment ("delta-radiomics") may be able to predict the pathological complete response (pCR).

METHODS

MRI scans (T1-weighted with fat-saturation and contrast-enhancement (T1FSGd) and T2-weighted with fat-saturation (T2FS)) of patients with STS of the extremities and trunk treated with neoadjuvant therapy were gathered from two independent institutions (training: 103, external testing: 53 patients). pCR was defined as <5% viable cells. After segmentation and preprocessing, 105 radiomic features were extracted. Delta-radiomic features were calculated by subtraction of features derived from MRI scans obtained before and after neoadjuvant therapy. After feature reduction, machine learning modeling was performed in 100 iterations of 3-fold nested cross-validation. Delta-radiomic models were compared with single timepoint models in the testing cohort.

RESULTS

The combined delta-radiomic models achieved the best area under the receiver operating characteristic curve (AUC) of 0.75. Pre-therapeutic tumor volume was the best conventional predictor (AUC 0.70). The T2FS-based delta-radiomic model had the most balanced classification performance with a balanced accuracy of 0.69. Delta-radiomic models achieved better reproducibility than single timepoint radiomic models, RECIST or the peri-therapeutic volume change. Delta-radiomic models were significantly associated with survival in multivariate Cox regression.

CONCLUSION

This exploratory analysis demonstrated that MRI-based delta-radiomics improves prediction of pCR over tumor volume and RECIST. Delta-radiomics may one day function as a biomarker for personalized treatment adaptations.

Implementing a Machine Learning Strategy to Predict Pathologic Response in Patients With Soft Tissue Sarcomas Treated With Neoadjuvant Chemotherapy

PURPOSE

Neoadjuvant chemotherapy (NAC) has been increasingly used in patients with locally advanced high-risk soft tissue sarcomas in the past decade, but definition and prognostic impact of a good histologic response (GHR) are lacking. Our aim was to investigate which histologic feature from the post-NAC surgical specimen independently correlated with metastatic relapse-free survival (MFS) in combination with clinical, radiologic, and pathologic features using a machine learning approach.

METHODS

This retrospective study included 175 consecutive patients (median age: 59 years, 75 women) with resectable disease, treated with anthracycline-based NAC between 1989 and 2015 in our sarcoma reference center, and with quantitative histopathologic analysis of the surgical specimen. The outcome of interest was the MFS. A multimodel, multivariate survival analysis was used to define GHR. The added prognostic value of GHR was investigated through the comparisons with the standard model (including histologic grade, size, and depth) and SARCULATOR nomogram using concordance indices (c-index) and Monte-Carlo cross-validation.

RESULTS

Seventy-two patients (72 of 175, 41.1%) had a metastatic relapse. Stepwise Cox regression, random survival forests, and least absolute shrinkage and selection operator-penalized Cox regression all converged toward the same definition for GHR, ie, < 5% stainable tumor cells. The five-year MFS probability was 1 (95% CI, 1 to 1) in patients with GHR versus 0.73 (95% CI, 0.65 to 0.81) in patients without GHR (log-rank P = .0122). The final prognostic model incorporating the GHR was significantly better than the standard model and SARCULATOR (average c-index in testing sets = 0.72 [95% CI, 0.61 to 0.82] v 0.57 [95% CI, 0.44 to 0.70] and 0.54 [95% CI, 0.45 to 0.64], respectively; P = .0414 and .0091).

CONCLUSION

Histologic response to NAC improves the prediction of MFS in patients with soft tissue sarcoma and represents a possible end point in future studies exploring innovative regimens in the neoadjuvant setting.

Manual versus semiautomatic segmentation of soft-tissue sarcomas on magnetic resonance imaging: evaluation of similarity and comparison of segmentation times

Objective

To evaluate the degree of similarity between manual and semiautomatic segmentation of soft-tissue sarcomas on magnetic resonance imaging (MRI).

Materials and Methods

This was a retrospective study of 15 MRI examinations of patients with histopathologically confirmed soft-tissue sarcomas acquired before therapeutic intervention. Manual and semiautomatic segmentations were performed by three radiologists, working independently, using the software 3D Slicer. The Dice similarity coefficient (DSC) and the Hausdorff distance were calculated in order to evaluate the similarity between manual and semiautomatic segmentation. To compare the two modalities in terms of the tumor volumes obtained, we also calculated descriptive statistics and intraclass correlation coefficients (ICCs).

Results

In the comparison between manual and semiautomatic segmentation, the DSC values ranged from 0.871 to 0.973. The comparison of the volumes segmented by the two modalities resulted in ICCs between 0.9927 and 0.9990. The DSC values ranged from 0.849 to 0.979 for intraobserver variability and from 0.741 to 0.972 for interobserver variability. There was no significant difference between the semiautomatic and manual modalities in terms of the segmentation times (p > 0.05).

Conclusion

There appears to be a high degree of similarity between manual and semiautomatic segmentation, with no significant difference between the two modalities in terms of the time required for segmentation.

Development and External Validation of Deep-Learning-Based Tumor Grading Models in Soft-Tissue Sarcoma Patients Using MR Imaging

BACKGROUND

In patients with soft-tissue sarcomas, tumor grading constitutes a decisive factor to determine the best treatment decision. Tumor grading is obtained by pathological work-up after focal biopsies. Deep learning (DL)-based imaging analysis may pose an alternative way to characterize STS tissue. In this work, we sought to non-invasively differentiate tumor grading into low-grade (G1) and high-grade (G2/G3) STS using DL techniques based on MR-imaging.

METHODS

Contrast-enhanced T1-weighted fat-saturated (T1FSGd) MRI sequences and fat-saturated T2-weighted (T2FS) sequences were collected from two independent retrospective cohorts (training: 148 patients, testing: 158 patients). Tumor grading was determined following the French Federation of Cancer Centers Sarcoma Group in pre-therapeutic biopsies. DL models were developed using transfer learning based on the DenseNet 161 architecture.

RESULTS

The T1FSGd and T2FS-based DL models achieved area under the receiver operator characteristic curve (AUC) values of 0.75 and 0.76 on the test cohort, respectively. T1FSGd achieved the best F1-score of all models (0.90). The T2FS-based DL model was able to significantly risk-stratify for overall survival. Attention maps revealed relevant features within the tumor volume and in border regions.

CONCLUSIONS

MRI-based DL models are capable of predicting tumor grading with good reproducibility in external validation.

CT and MRI radiomics of bone and soft-tissue sarcomas: a systematic review of reproducibility and validation strategies

Background

Feature reproducibility and model validation are two main challenges of radiomics. This study aims to systematically review radiomic feature reproducibility and predictive model validation strategies in studies dealing with CT and MRI radiomics of bone and soft-tissue sarcomas. The ultimate goal is to promote achieving a consensus on these aspects in radiomic workflows and facilitate clinical transferability.

Results

Out of 278 identified papers, forty-nine papers published between 2008 and 2020 were included. They dealt with radiomics of bone (n = 12) or soft-tissue (n = 37) tumors. Eighteen (37%) studies included a feature reproducibility analysis. Inter-/intra-reader segmentation variability was the theme of reproducibility analysis in 16 (33%) investigations, outnumbering the analyses focused on image acquisition or post-processing (n = 2, 4%). The intraclass correlation coefficient was the most commonly used statistical method to assess reproducibility, which ranged from 0.6 and 0.9. At least one machine learning validation technique was used for model development in 25 (51%) papers, and K-fold cross-validation was the most commonly employed. A clinical validation of the model was reported in 19 (39%) papers. It was performed using a separate dataset from the primary institution (i.e., internal validation) in 14 (29%) studies and an independent dataset related to different scanners or from another institution (i.e., independent validation) in 5 (10%) studies.

Conclusions

The issues of radiomic feature reproducibility and model validation varied largely among the studies dealing with musculoskeletal sarcomas and should be addressed in future investigations to bring the field of radiomics from a preclinical research area to the clinical stage.

Intensity harmonization techniques influence radiomics features and radiomics-based predictions in sarcoma patients

Intensity harmonization techniques (IHT) are mandatory to homogenize multicentric MRIs before any quantitative analysis because signal intensities (SI) do not have standardized units. Radiomics combine quantification of tumors' radiological phenotype with machine-learning to improve predictive models, such as metastastic-relapse-free survival (MFS) for sarcoma patients. We post-processed the initial T2-weighted-imaging of 70 sarcoma patients by using 5 IHTs and extracting 45 radiomics features (RFs), namely: classical standardization (IHT_std), standardization per adipose tissue SIs (IHT_fat), histogram-matching with a patient histogram (IHT_HM.1), with the average histogram of the population (IHT_HM.All) and plus ComBat method (IHT_HM.All.C), which provided 5 radiomics datasets in addition to the original radiomics dataset without IHT (No-IHT). We found that using IHTs significantly influenced all RFs values (p-values: < 0.0001-0.02). Unsupervised clustering performed on each radiomics dataset showed that only clusters from the No-IHT, IHT_std, IHT_HM.All, and IHTHM.All.C datasets significantly correlated with MFS in multivariate Cox models (p = 0.02, 0.007, 0.004 and 0.02, respectively). We built radiomics-based supervised models to predict metastatic relapse at 2-years with a training set of 50 patients. The models performances varied markedly depending on the IHT in the validation set (range of AUROC from 0.688 with IHT_std to 0.823 with IHT_HM.1). Hence, the use of intensity harmonization and the related technique should be carefully detailed in radiomics post-processing pipelines as it can profoundly affect the reproducibility of analyses.

Systematic review of sarcomas radiomics studies: Bridging the gap between concepts and clinical applications?

OBJECTIVES

Sarcomas are a model for intra- and inter-tumoral heterogeneities making them particularly suitable for radiomics analyses. Our purposes were to review the aims, methods and results of radiomics studies involving sarcomas METHODS: Pubmed and Web of Sciences databases were searched for radiomics or textural studies involving bone, soft-tissues and visceral sarcomas until June 2020. Two radiologists evaluated their objectives, results and quality of their methods, imaging pre-processing and machine-learning workflow helped by the items of the Quality Assessment of Diagnostic Accuracy Studies (QUADAS-2), Image Biomarker Standardization Initiative (IBSI) and 'Radiomics Quality Score' (RQS). Statistical analyses included inter-reader agreements, correlations between methodological assessments, scientometrics indices, and their changes over years, and between RQS, number of patients and models performance.

RESULTS

Fifty-two studies were included involving: soft-tissue sarcomas (29/52, 55.8 %), bone sarcomas (15/52, 28.8 %), gynecological sarcomas (6/52, 11.5 %) and mixed sarcomas (2/52, 3.8 %), mostly imaged with MRI (36/52, 69.2 %), for a total of distinct patients. Median RQS was 4.5 (28.4 % of the maximum, range: -7 - 17). Performances of predictive models and number of patients negatively correlated (p = 0.027). None of the studies detailed all the items from the IBSI guidelines. There was a significant increase in studies' impact factors since the establishing of the RQS in 2017 (p = 0.038).

CONCLUSION

Although showing promising results, further efforts are needed to make sarcoma radiomics studies reproducible with an acceptable level of evidence. A better knowledge of the RQS and IBSI reporting guidelines could improve the quality of sarcoma radiomics studies and accelerate clinical applications.