Robustness of radiomic features in magnetic resonance imaging for patients with glioblastoma: multi-center study

Impact of intensity standardisation and ComBat batch size on clinical-radiomic prognostic models performance in a multi-centre study of patients with glioblastoma

PURPOSE

To assess the effect of different intensity standardisation techniques (ISTs) and ComBat batch sizes on radiomics survival model performance and stability in a heterogenous, multi-centre cohort of patients with glioblastoma (GBM).

METHODS

Multi-centre pre-operative MRI acquired between 2014 and 2020 in patients with IDH-wildtype unifocal WHO grade 4 GBM were retrospectively evaluated. WhiteStripe (WS), Nyul histogram matching (HM), and Z-score (ZS) ISTs were applied before radiomic feature (RF) extraction. RFs were realigned using ComBat and minimum batch size (MBS) of 5, 10, or 15 patients. Cox proportional hazards models for overall survival (OS) prediction were produced using five different selection strategies and the impact of IST and MBS was evaluated using bootstrapping. Calibration, discrimination, relative explained variation, and model fit were assessed. Instability was evaluated using 95% confidence intervals (95% CIs), feature selection frequency and calibration curves across the bootstrap resamples.

RESULTS

One hundred ninety-five patients were included. Median OS = 13 (95% CI: 12-14) months. Twelve to fourteen unique MRI protocols were used per MRI sequence. HM and WS produced the highest relative increase in model discrimination, explained variation and model fit but IST choice did not greatly impact on stability, nor calibration. Larger ComBat batches improved discrimination, model fit, and explained variation but higher MBS (reduced sample size) reduced stability (across all performance metrics) and reduced calibration accuracy.

CONCLUSION

Heterogenous, real-world GBM data poses a challenge to the reproducibility of radiomics. ComBat generally improved model performance as MBS increased but reduced stability and calibration. HM and WS tended to improve model performance.

KEY POINTS

Question ComBat harmonisation of RFs and intensity standardisation of MRI have not been thoroughly evaluated in multicentre, heterogeneous GBM data. Findings The addition of ComBat and ISTs can improve discrimination, relative model fit, and explained variance but degrades the calibration and stability of survival models. Clinical relevance Radiomics risk prediction models in real-world, multicentre contexts could be improved by ComBat and ISTs, however, this degrades calibration and prediction stability and this must be thoroughly investigated before patients can be accurately separated into different risk groups.

T2-weighted MRI radiomics for the prediction of pediatric and young adult rhabdomyosarcoma alveolar subtype and distant metastasis: a pilot study

INTRODUCTION

Rhabdomyosarcomas are the most common soft tissue sarcoma in children. While treatment outcomes have improved, risk-based therapy classification relies on staging and tumor subtypes for therapeutic planning.

OBJECTIVE

This study investigated the utility of T2-weighted MR radiomics features and machine learning models in identifying the presence of distant metastasis and alveolar histological subtypes at baseline imaging in children diagnosed with rhabdomyosarcoma.

MATERIALS AND METHODS

This retrospective cross-sectional study utilized MRIs from 86 patients, 49 (median age (IQR) 59 months (37-161), alveolar subtype=15, distant metastasis=9) of whom had been imaged at outside imaging centers (training set); and 37 (median age 52 months (24-164), alveolar subtype=14, distant metastasis=8) of whom were imaged at our institution (holdout validation set). Radiomic features were extracted from T2-weighted images. We selected features that demonstrated intra-scan repeatability and used maximum relevance and minimum redundancy supervised feature selection to identify the 50 most important features. Lasso logistic regression and support vector machine (SVM) classifiers were trained to predict binary outcomes. The median of all predictions for a given patient was used as patient-level predictions. DeLong's test compared the area under the receiver operating characteristic curves (AUC). Cut-offs obtained by maximizing the Youden index were evaluated on an external validation set, and accuracy metrics were reported.

RESULTS

On the validation set, the Lasso and SVM classifiers obtained patient level AUCs of 0.76 (95% CI 0.59-0.94) and 0.73 (0.54-0.92), respectively, in predicting alveolar subtype, with the Lasso regressor obtaining 71.4% (41.9-91.6) sensitivity and 60.9% (38.5-80.3) specificity. When predicting the presence of distant metastasis, the Lasso and SVM classifier had AUCs of 0.81 (0.67-0.95) and 0.77 (0.58-0.97), respectively. There were no differences between model performance (P>0.05). A total of 12 and 18 features had nonzero coefficients in the Lasso regressors for predicting alveolar subtype and tumor metastasis, respectively.

CONCLUSION

MRI radiomics from baseline T2-weighted MRI demonstrated potential in predicting alveolar subtype and distant metastatic disease at presentation. Larger studies are needed to explore multinomial multiclass models for better prognostication of pediatric rhabdomyosarcomas.

Robust vs. Non-robust radiomic features: the quest for optimal machine learning models using phantom and clinical studies

PURPOSE

This study aimed to select robust features against lung motion in a phantom study and use them as input to feature selection algorithms and machine learning classifiers in a clinical study to predict the lymphovascular invasion (LVI) of non-small cell lung cancer (NSCLC). The results of robust features were also compared with conventional techniques without considering the robustness of radiomic features.

METHODS

An in-house developed lung phantom was developed with two 22mm lesion sizes based on a clinical study. A specific motor was built to simulate motion in two orthogonal directions. Lesions of both clinical and phantom studies were segmented using a Fuzzy C-means-based segmentation algorithm. After inducing motion and extracting 105 radiomic features in 4 feature sets, including shape, first-, second-, and higher-order statistics features from each region of interest (ROI) of the phantom image, statistical analyses were performed to select robust features against motion. Subsequently, these robust features and a total of 105 radiomic features were extracted from 126 clinical data. Various feature selection (FS) and multiple machine learning (ML) classifiers were implemented to predict the LVI of NSCLC, followed by comparing the results of predicting LVI using robust features with common conventional techniques not considering the robustness of radiomic features.

RESULTS

Our results demonstrated that selecting robust features as input to FS algorithms and ML classifiers surges the sensitivity, which has a gentle negative effect on the accuracy and the area under the curve (AUC) of predictions compared with commonly used methods in 12 of 15 outcomes. The top performance of the LVI prediction was achieved by the NB classifier and RFE FS without considering the robustness of radiomic features with 95% area under the curve of AUC, 67% accuracy, and 100% sensitivity. Moreover, the top performance of the LVI prediction using robust features belonged to the NB classifier and Boruta feature selection with 92% AUC, 86% accuracy, and 100% sensitivity.

CONCLUSION

Robustness over various influential factors is critical and should be considered in a radiomic study. Selecting robust features is a solution to overcome the low reproducibility of radiomic features. Although setting robust features against motion in a phantom study has a minor negative impact on the accuracy and AUC of LVI prediction, it boosts the sensitivity of prediction to a large extent.

Reproducibility and repeatability of 18F-(2S, 4R)-4-fluoroglutamine PET imaging in preclinical oncology models

INTRODUCTION

Measurement of repeatability and reproducibility (R&R) is necessary to realize the full potential of positron emission tomography (PET). Several studies have evaluated the reproducibility of PET using 18F-FDG, the most common PET tracer used in oncology, but similar studies using other PET tracers are scarce. Even fewer assess agreement and R&R with statistical methods designed explicitly for the task. 18F-(2S, 4R)-4-fluoro-glutamine (18F-Gln) is a PET tracer designed for imaging glutamine uptake and metabolism. This study illustrates high reproducibility and repeatability with 18F-Gln for in vivo research.

METHODS

Twenty mice bearing colorectal cancer cell line xenografts were injected with ~9 MBq of 18F-Gln and imaged in an Inveon microPET. Three individuals analyzed the tumor uptake of 18F-Gln using the same set of images, the same image analysis software, and the same analysis method. Scans were randomly re-ordered for a second repeatability measurement 6 months later. Statistical analyses were performed using the methods of Bland and Altman (B&A), Gauge Reproducibility and Repeatability (Gauge R&R), and Lin's Concordance Correlation Coefficient. A comprehensive equivalency test, designed to reject a null hypothesis of non-equivalence, was also conducted.

RESULTS

In a two-way random effects Gauge R&R model, variance among mice and their measurement variance were 0.5717 and 0.024. Reproducibility and repeatability accounted for 31% and 69% of the total measurement error, respectively. B&A repeatability coefficients for analysts 1, 2, and 3 were 0.16, 0.35, and 0.49. One-half B&A agreement limits between analysts 1 and 2, 1 and 3, and 2 and 3 were 0.27, 0.47, and 0.47, respectively. The mean square deviation and total deviation index were lowest for analysts 1 and 2, while coverage probabilities and coefficients of the individual agreement were highest. Finally, the definitive agreement inference hypothesis test for equivalency demonstrated that all three confidence intervals for the average difference of means from repeated measures lie within our a priori limits of equivalence (i.e. ± 0.5%ID/g).

CONCLUSIONS

Our data indicate high individual analyst and laboratory-level reproducibility and repeatability. The assessment of R&R using the appropriate methods is critical and should be adopted by the broader imaging community.

Radiomics score derived from T1-w/T2-w ratio image can predict motor symptom progression in Parkinson's disease

OBJECTIVES

To clarify the association between a radiomics score (Rad-score) derived from T1-weighted signal intensity to T2-weighted signal intensity (T1-w/T2-w) ratio images and the progression of motor symptoms in Parkinson's disease (PD).

MATERIALS AND METHODS

This retrospective study included patients with PD enrolled in the Parkinson's Progression Markers Initiative. The Movement Disorders Society-Unified Parkinson's Disease Rating Scale Part III score ≥ 33 and/or Hoehn and Yahr stage ≥ 3 indicated motor function decline. The Rad-score was constructed using radiomics features extracted from T1-w/T2-w ratio images. The Kaplan-Meier analysis and Cox regression analyses were used to assess the time differences in motor function decline between the high and low Rad-score groups.

RESULTS

A total of 171 patients with PD were divided into training (n = 101, mean age at baseline, 61.6 ± 9.3 years) and testing (n = 70, mean age at baseline, 61.6 ± 10 years). The patients in the high Rad-score group had a shorter time to motor function decline than those in the low Rad-score group in the training dataset (log-rank test, p < 0.001) and testing dataset (log-rank test, p < 0.001). The multivariate Cox regression using the Rad-score and clinical factors revealed a significant association between the Rad-score and motor function decline in the training dataset (HR = 2.368, 95%CI:1.423-3.943, p < 0.001) and testing dataset (HR = 2.931, 95%CI:1.472-5.837, p = 0.002).

CONCLUSION

Rad-scores based on radiomics features derived from T1-w/T2-w ratio images were associated with the progression of motor symptoms in PD.

CLINICAL RELEVANCE STATEMENT

The radiomics score derived from the T1-weighted/T2-weighted ratio images offers a predictive tool for assessing the progression of motor symptom in patients with PD.

KEY POINTS

Radiomics score derived from T1-weighted/T2-weighted ratio images is correlated with the motor symptoms of Parkinson's disease. A high radiomics score correlated with faster motor function decline in patients with Parkinson's disease. The proposed radiomics score offers predictive insight into the progression of motor symptoms of Parkinson's disease.

A Comprehensive Evaluation of Radiomic Features in Normal Brain Magnetic Resonance Imaging: Investigating Robustness and Region Variations

Background

Despite extensive research on various brain diseases, a few studies have focused on radiomic feature distribution in healthy brain images. The present study applied a novel radiomic framework to investigate the robustness and baseline values of radiomic features in normal brain magnetic resonance imaging (MRIs) regions.

Materials and Methods

Analyses were performed on T1 and T2 images including 276 normal brains and 14 healthy volunteers were scanned with three scanners using the same protocols. The images were divided into 1024 three-dimensional nonoverlap patches with the same pixel size. Seven patches located in the thalamus, putamen, hippocampus and brain stem were selected as volume of interest (VOI). Eighty-five radiomic features were generated. To investigate the variation of features across VOIs, the analysis of variance was performed and coefficient of variation (COV) and intraclass correlation coefficient (ICC) were explored to examine the features repeatability.

Results

Thalamus (right and left) and hippocampus (left) resulted in more stable features (COV ≤ 6%) in T1 and T2 images, respectively. The inter-scanner ICC analysis demonstrated the features of T2 sequences represented more repeatable results and the brain stem and thalamus (both T1 and T2) showed particularly high repeatability (higher ICC values). Robust results (ICC ≥ 0.9) were identified for energy and range features of the first order class and several textures features across different brain regions.

Conclusion

Our results indicated the baselines of the repeatable texture features in healthy brain structural MRI highlighting inter-scanner stability. According to the findings, MRI sequencing and VOI location impact feature robustness and should be considered in brain radiomic studies.

Magnetic resonance imaging radiomic features stability in brain metastases: Impact of image preprocessing, image-, and feature-level harmonization

Background and purpose

Magnetic resonance imaging (MRI) scans are highly sensitive to acquisition and reconstruction parameters which affect feature stability and model generalizability in radiomic research. This work aims to investigate the effect of image pre-processing and harmonization methods on the stability of brain MRI radiomic features and the prediction performance of radiomic models in patients with brain metastases (BMs).

Materials and methods

Two T1 contrast enhanced brain MRI data-sets were used in this study. The first contained 25 BMs patients with scans at two different time points and was used for features stability analysis. The effect of gray level discretization (GLD), intensity normalization (Z-score, Nyul, WhiteStripe, and in house-developed method named N-Peaks), and ComBat harmonization on features stability was investigated and features with intraclass correlation coefficient >0.8 were considered as stable. The second data-set containing 64 BMs patients was used for a classification task to investigate the informativeness of stable features and the effects of harmonization methods on radiomic model performance.

Results

Applying fixed bin number (FBN) GLD, resulted in higher number of stable features compare to fixed bin size (FBS) discretization (10 ± 5.5 % higher). `Harmonization in feature domain improved the stability for non-normalized and normalized images with Z-score and WhiteStripe methods. For the classification task, keeping the stable features resulted in good performance only for normalized images with N-Peaks along with FBS discretization.

Conclusions

To develop a robust MRI based radiomic model we recommend using an intensity normalization method based on a reference tissue (e.g N-Peaks) and then using FBS discretization.

Deeply 3D-T1-TFE hypointense voxels are characteristic of phase-rim lesions in multiple sclerosis

OBJECTIVES

The development of new drugs for the treatment of progressive multiple sclerosis (MS) highlights the need for new prognostic biomarkers. Phase-rim lesions (PRLs) have been proposed as markers of progressive disease but are difficult to identify and quantify. Previous studies have identified T1-hypointensity in PRLs. The aim of this study was to compare the intensity profiles of PRLs and non-PRL white-matter lesions (nPR-WMLs) on three-dimensional T1-weighted turbo field echo (3DT1TFE) MRI. We then evaluated the performance of a derived metric as a surrogate for PRLs as potential markers for risk of disease progression.

METHODS

This study enrolled a cohort of relapsing-remitting (n = 10) and secondary progressive MS (n = 10) patients for whom 3 T MRI was available. PRLs and nPR-WMLs were segmented, and voxel-wise normalized T1-intensity histograms were analyzed. The lesions were divided equally into training and test datasets, and the fifth-percentile (p5)-normalized T1-intensity of each lesion was compared between groups and used for classification prediction.

RESULTS

Voxel-wise histogram analysis showed a unimodal histogram for nPR-WMLs and a bimodal histogram for PRLs with a large peak in the hypointense limit. Lesion-wise analysis included 1075 nPR-WMLs and 39 PRLs. The p5 intensity of PRLs was significantly lower than that of nPR-WMLs. The T1 intensity-based PRL classifier had a sensitivity of 0.526 and specificity of 0.959.

CONCLUSIONS

Profound hypointensity on 3DT1TFE MRI is characteristic of PRLs and rare in other white-matter lesions. Given the widespread availability of T1-weighted imaging, this feature might serve as a surrogate biomarker for smoldering inflammation.

CLINICAL RELEVANCE STATEMENT

Quantitative analysis of 3DT1TFE may detect deeply hypointense voxels in multiple sclerosis lesions, which are highly specific to PRLs. This could serve as a specific indicator of smoldering inflammation in MS, aiding in early detection of disease progression.

KEY POINTS

• Phase-rim lesions (PRLs) in multiple sclerosis present a characteristic T1-hypointensity on 3DT1TFE MRI. • Intensity-normalized 3DT1TFE can be used to systematically identify and quantify these deeply hypointense foci. • Deep T1-hypointensity may act as an easily detectable, surrogate marker for PRLs.

MRI-based Machine Learning Radiomics Can Predict CSF1R Expression Level and Prognosis in High-grade Gliomas

The purpose of this study is to predict the mRNA expression of CSF1R in HGG non-invasively using MRI (magnetic resonance imaging) omics technology and to evaluate the correlation between the established radiomics model and prognosis. We investigated the predictive value of CSF1R in the Cancer Genome Atlas (TCGA) and The Cancer Imaging Archive (TCIA) database. The Support vector machine (SVM) and the Logistic regression (LR) algorithms were used to create a radiomics_score (Rad_score), respectively. The effectiveness and performance of the radiomics model was assessed in the training (n = 89) and tenfold cross-validation sets. We further analyzed the correlation between Rad_score and macrophage-related genes using Spearman correlation analysis. A radiomics nomogram combining the clinical factors and Rad_score was constructed to validate the radiomic signatures for individualized survival estimation and risk stratification. The results showed that CSF1R expression was markedly elevated in HGG tissues, which was related to worse prognosis. CSF1R expression was closely related to the abundance of infiltrating immune cells, such as macrophages. We identified nine features for establishing a radiomics model. The radiomics model predicting CSF1R achieved high AUC in training (0.768 in SVM and 0.792 in LR) and tenfold cross-validation sets (0.706 in SVM and 0.717 in LR). Rad_score was highly associated with tumor-related macrophage genes. A radiomics nomogram combining the Rad_score and clinical factors was constructed and revealed satisfactory performance. MRI-based Rad_score is a novel way to predict CSF1R expression and prognosis in high-grade glioma patients. The radiomics nomogram could optimize individualized survival estimation for HGG patients.

Role of MRI radiomics for the prediction of MYCN amplification in neuroblastomas

OBJECTIVES

We evaluate MR radiomics and develop machine learning-based classifiers to predict MYCN amplification in neuroblastomas.

METHODS

A total of 120 patients with neuroblastomas and baseline MR imaging examination available were identified of whom 74 (mean age ± standard deviation [SD] of 6 years and 2 months ± 4 years and 9 months; 43 females and 31 males, 14 MYCN amplified) underwent imaging at our institution. This was therefore used to develop radiomics models. The model was tested in a cohort of children with the same diagnosis but imaged elsewhere (n = 46, mean age ± SD: 5 years 11 months ± 3 years 9 months, 26 females and 14 MYCN amplified). Whole tumour volumes of interest were adopted to extract first-order histogram and second-order radiomics features. Interclass correlation coefficient and maximum relevance and minimum redundancy algorithm were applied for feature selection. Logistic regression, support vector machine, and random forest were employed as the classifiers. Receiver operating characteristic (ROC) analysis was performed to evaluate the diagnostic accuracy of the classifiers on the external test set.

RESULTS

The logistic regression model and the random forest both showed an AUC of 0.75. The support vector machine classifier obtained an AUC of 0.78 on the test set with a sensitivity of 64% and a specificity of 72%.

CONCLUSION

The study provides preliminary retrospective evidence demonstrating the feasibility of MRI radiomics in predicting MYCN amplification in neuroblastomas. Future studies are needed to explore the correlation between other imaging features and genetic markers and to develop multiclass predictive models.

KEY POINTS

• MYCN amplification in neuroblastomas is an important determinant of disease prognosis. • Radiomics analysis of pre-treatment MR examinations can be used to predict MYCN amplification in neuroblastomas. • Radiomics machine learning models showed good generalisability to external test set, demonstrating reproducibility of the computational models.

What Matters in Radiological Image Segmentation? Effect of Segmentation Errors on the Diagnostic Related Features

Segmentation is a crucial step in extracting the medical image features for clinical diagnosis. Though multiple metrics have been proposed to evaluate the segmentation performance, there is no clear study on how or to what extent the segmentation errors will affect the diagnostic related features used in clinical practice. Therefore, we proposed a segmentation robustness plot (SRP) to build the link between segmentation errors and clinical acceptance, where relative area under the curve (R-AUC) was designed to help clinicians to identify the robust diagnostic related image features. In experiments, we first selected representative radiological series from time series (cardiac first-pass perfusion) and spatial series (T2 weighted images on brain tumors) of magnetic resonance images, respectively. Then, dice similarity coefficient (DSC) and Hausdorff distance (HD), as the widely used evaluation metrics, were used to systematically control the degree of the segmentation errors. Finally, the differences between diagnostic related image features extracted from the ground truth and the derived segmentation were analyzed, using the statistical method large sample size T-test to calculate the corresponding p values. The results are denoted in the SRP, where the x-axis indicates the segmentation performance using the aforementioned evaluation metric, and the y-axis shows the severity of the corresponding feature changes, which are expressed in either the p values for a single case or the proportion of patients without significant change. The experimental results in SRP show that when DSC is above 0.95 and HD is below 3 mm, the segmentation errors will not change the features significantly in most cases. However, when segmentation gets worse, additional metrics are required for further analysis. In this way, the proposed SRP indicates the impact of the segmentation errors on the severity of the corresponding feature changes. By using SRP, one could easily define the acceptable segmentation errors in a challenge. Additionally, the R-AUC calculated from SRP provides an objective reference to help the selection of reliable features in image analysis.

The role of parametric feature maps to correct different volume of interest sizes: an in vivo liver MRI study

BACKGROUND

Different volume of interest (VOI) sizes influence radiomic features. This study examined if translating images into feature maps before feature sampling could compensate for these effects in liver magnetic resonance imaging (MRI).

METHODS

T1- and T2-weighted sequences from three different scanners (two 3-T scanners, one 1.5-T scanner) of 66 patients with normal abdominal MRI were included retrospectively. Three differently sized VOIs (10, 20, and 30 mm in diameter) were drawn in the liver parenchyma (right lobe), excluding adjacent structures. Ninety-three features were extracted conventionally using PyRadiomics. All images were also converted to 93 parametric feature maps using a pretested software. Agreement between the three VOI sizes was assessed with overall concordance correlation coefficients (OCCCs), while OCCCs > 0.85 were rated reproducible. OCCCs were calculated twice: for the VOI sizes of 10, 20, and 30 mm and for those of 20 and 30 mm.

RESULTS

When extracted from original images, only 4 out of the 93 features were reproducible across all VOI sizes in T1- and T2-weighted images. When the smallest VOI was excluded, 5 features (T1-weighted) and 7 features (T2-weighted) were reproducible. Extraction from parametric maps increased the number of reproducible features to 9 (T1- and T2-weighted) across all VOIs. Excluding the 10-mm VOI, reproducibility improved to 16 (T1-weighted) and 55 features (T2-weighted). The stability of all other features also increased in feature maps.

CONCLUSIONS

Translating images into parametric maps before feature extraction improves reproducibility across different VOI sizes in normal liver MRI.

RELEVANCE STATEMENT

The size of the segmented VOI influences the feature quantity of radiomics, while software-based conversion of images into parametric feature maps before feature sampling improves reproducibility across different VOI sizes in MRI of normal liver tissue.

KEY POINTS

• Parametric feature maps can compensate for different VOI sizes. • The effect seems dependent on the VOI sizes and the MRI sequence. • Feature maps can visualize features throughout the entire image stack.

Advancements in MRI-Based Radiomics and Artificial Intelligence for Prostate Cancer: A Comprehensive Review and Future Prospects

The use of multiparametric magnetic resonance imaging (mpMRI) has become a common technique used in guiding biopsy and developing treatment plans for prostate lesions. While this technique is effective, non-invasive methods such as radiomics have gained popularity for extracting imaging features to develop predictive models for clinical tasks. The aim is to minimize invasive processes for improved management of prostate cancer (PCa). This study reviews recent research progress in MRI-based radiomics for PCa, including the radiomics pipeline and potential factors affecting personalized diagnosis. The integration of artificial intelligence (AI) with medical imaging is also discussed, in line with the development trend of radiogenomics and multi-omics. The survey highlights the need for more data from multiple institutions to avoid bias and generalize the predictive model. The AI-based radiomics model is considered a promising clinical tool with good prospects for application.

Radiogenomics Analysis Linking Multiparametric MRI and Transcriptomics in Prostate Cancer

Prostate cancer (PCa) is a highly prevalent cancer type with a heterogeneous prognosis. An accurate assessment of tumor aggressiveness can pave the way for tailored treatment strategies, potentially leading to better outcomes. While tumor aggressiveness is typically assessed based on invasive methods (e.g., biopsy), radiogenomics, combining diagnostic imaging with genomic information can help uncover aggressive (imaging) phenotypes, which in turn can provide non-invasive advice on individualized treatment regimens. In this study, we carried out a parallel analysis on both imaging and transcriptomics data in order to identify features associated with clinically significant PCa (defined as an ISUP grade ≥ 3), subsequently evaluating the correlation between them. Textural imaging features were extracted from multi-parametric MRI sequences (T2W, DWI, and DCE) and combined with DCE-derived parametric pharmacokinetic maps obtained using magnetic resonance dispersion imaging (MRDI). A transcriptomic analysis was performed to derive functional features on transcription factors (TFs), and pathway activity from RNA sequencing data, here referred to as transcriptomic features. For both the imaging and transcriptomic features, different machine learning models were separately trained and optimized to classify tumors in either clinically insignificant or significant PCa. These models were validated in an independent cohort and model performance was used to isolate a subset of relevant imaging and transcriptomic features to be further investigated. A final set of 31 imaging features was correlated to 33 transcriptomic features obtained on the same tumors. Five significant correlations (p < 0.05) were found, of which, three had moderate strength (|r| ≥ 0.5). The strongest significant correlations were seen between a perfusion-based imaging feature-MRDI A median-and the activities of the TFs STAT6 (-0.64) and TFAP2A (-0.50). A higher-order T2W textural feature was also significantly correlated to the activity of the TF STAT6 (-0.58). STAT6 plays an important role in controlling cell proliferation and migration. Loss of the AP2alpha protein expression, quantified by TFAP2A, has been strongly associated with aggressiveness and progression in PCa. According to our findings, a combination of texture features extracted from T2W and DCE, as well as perfusion-based pharmacokinetic features, can be considered for the prediction of clinically significant PCa, with the pharmacokinetic MRDI A feature being the most correlated with the underlying transcriptomic information. These results highlight a link between quantitative imaging features and the underlying transcriptomic landscape of prostate tumors.

Deriving quantitative information from multiparametric MRI via Radiomics: Evaluation of the robustness and predictive value of radiomic features in the discrimination of low-grade versus high-grade gliomas with machine learning

PURPOSE

Analysis pipelines based on the computation of radiomic features on medical images are widely used exploration tools across a large variety of image modalities. This study aims to define a robust processing pipeline based on Radiomics and Machine Learning (ML) to analyze multiparametric Magnetic Resonance Imaging (MRI) data to discriminate between high-grade (HGG) and low-grade (LGG) gliomas.

METHODS

The dataset consists of 158 multiparametric MRI of patients with brain tumor publicly available on The Cancer Imaging Archive, preprocessed by the BraTS organization committee. Three different types of image intensity normalization algorithms were applied and 107 features were extracted for each tumor region, setting the intensity values according to different discretization levels. The predictive power of radiomic features in the LGG versus HGG categorization was evaluated by using random forest classifiers. The impact of the normalization techniques and of the different settings in the image discretization was studied in terms of the classification performances. A set of MRI-reliable features was defined selecting the features extracted according to the most appropriate normalization and discretization settings.

RESULTS

The results show that using MRI-reliable features improves the performance in glioma grade classification (AUC=0.93±0.05) with respect to the use of raw (AUC=0.88±0.08) and robust features (AUC=0.83±0.08), defined as those not depending on image normalization and intensity discretization.

CONCLUSIONS

These results confirm that image normalization and intensity discretization strongly impact the performance of ML classifiers based on radiomic features. Thus, special attention should be provided in the image preprocessing step before typical radiomic and ML analysis are carried out.

How Radiomics Can Improve Breast Cancer Diagnosis and Treatment

Recent technological advances in the field of artificial intelligence hold promise in addressing medical challenges in breast cancer care, such as early diagnosis, cancer subtype determination and molecular profiling, prediction of lymph node metastases, and prognostication of treatment response and probability of recurrence. Radiomics is a quantitative approach to medical imaging, which aims to enhance the existing data available to clinicians by means of advanced mathematical analysis using artificial intelligence. Various published studies from different fields in imaging have highlighted the potential of radiomics to enhance clinical decision making. In this review, we describe the evolution of AI in breast imaging and its frontiers, focusing on handcrafted and deep learning radiomics. We present a typical workflow of a radiomics analysis and a practical "how-to" guide. Finally, we summarize the methodology and implementation of radiomics in breast cancer, based on the most recent scientific literature to help researchers and clinicians gain fundamental knowledge of this emerging technology. Alongside this, we discuss the current limitations of radiomics and challenges of integration into clinical practice with conceptual consistency, data curation, technical reproducibility, adequate accuracy, and clinical translation. The incorporation of radiomics with clinical, histopathological, and genomic information will enable physicians to move forward to a higher level of personalized management of patients with breast cancer.

Impact of Wavelet Kernels on Predictive Capability of Radiomic Features: A Case Study on COVID-19 Chest X-ray Images

Radiomic analysis allows for the detection of imaging biomarkers supporting decision-making processes in clinical environments, from diagnosis to prognosis. Frequently, the original set of radiomic features is augmented by considering high-level features, such as wavelet transforms. However, several wavelets families (so called kernels) are able to generate different multi-resolution representations of the original image, and which of them produces more salient images is not yet clear. In this study, an in-depth analysis is performed by comparing different wavelet kernels and by evaluating their impact on predictive capabilities of radiomic models. A dataset composed of 1589 chest X-ray images was used for COVID-19 prognosis prediction as a case study. Random forest, support vector machine, and XGBoost were trained (on a subset of 1103 images) after a rigorous feature selection strategy to build-up the predictive models. Next, to evaluate the models generalization capability on unseen data, a test phase was performed (on a subset of 486 images). The experimental findings showed that Bior1.5, Coif1, Haar, and Sym2 kernels guarantee better and similar performance for all three machine learning models considered. Support vector machine and random forest showed comparable performance, and they were better than XGBoost. Additionally, random forest proved to be the most stable model, ensuring an appropriate balance between sensitivity and specificity.

MRI radiomics and potential applications to glioblastoma

MRI plays an important role in the evaluation of glioblastoma, both at initial diagnosis and follow up after treatment. Quantitative analysis via radiomics can augment the interpretation of MRI in terms of providing insights regarding the differential diagnosis, genotype, treatment response, and prognosis. The various MRI radiomic features of glioblastoma are reviewed in this article.