Improved prediction of outcome in Parkinson's disease using radiomics analysis of longitudinal DAT SPECT images

Individualized diagnosis of Parkinson's disease based on multivariate magnetic resonance imaging radiomics and clinical indexes

Objective

To explore MRI-based radiomics models, integrating clinical characteristics, for differential diagnosis of Parkinson's disease (PD) to evaluate their diagnostic performance.

Methods

A total of 256 participants [153 PD, 103 healthy controls (HCs)] from the First Affiliated Hospital of Wenzhou Medical Hospital, were enrolled as the training set, and 120 subjects (74 PD, 46 HCs) from the PPMI dataset served as the test set. Radiomics features were extracted from structural MRI (T1WI and T2-FLair). Support Vector Machine (SVM) classifiers were developed using MRI radiomics data from both monomodal and multimodal radiomics models. The clinical-radiomics model was constructed by integrating clinical variables, including UPDRS, Hoehn-Yahr stage, age, sex, and MMSE scores. Receiver operating characteristic (ROC) curves were generated to evaluate the performance of the models. Decision curve analysis (DCA) was performed to access the clinical usefulness of the models.

Results

In the training set, the T2-FLair and T1WI radiomics model achieved an AUC of 0.896 (95% CI, 0.812-0.900) and 0.899 (95% CI, 0.818-0.908), respectively. The double-sequence radiomics model demonstrated superior diagnostic performance, with an AUC of 0.965 (95% CI, 0.885-0.978) in the training set and an AUC of 0.852 (95% CI, 0.748-0.910) in the test set. The integrated clinical-radiomics model showed enhanced diagnostic accuracy, with AUC = 0.983 (95% CI, 0.897-0.996) in the training set and AUC = 0.837 (95% CI, 0.786-0.902) in the test set. Rad-scores derived from the radiomics model were significantly correlated with diagnostic outcomes (P < 0.001). DCA confirmed the substantial clinical usefulness of the clinical-radiomics integrated model.

Conclusion

The integrated clinical-radiomics model offered superior diagnostic performance compared to models based relying solely on imaging or clinical data, underscoring its potential as a non-invasive and effective tool in routine clinical practice for the early diagnosis of PD.

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.

Impact of Harmonization on MRI Radiomics Feature Variability Across Preprocessing Methods for Parkinson's Disease Motor Subtype Classification

This study aimed to assess the reproducibility of MRI-derived radiomic features across multiple preprocessing methods for classifying Parkinson's disease (PD) motor subtypes and to evaluate the impact of ComBat harmonization on feature stability and machine learning performance. T1-weighted MRI scans from 140 PD patients (70 tremor-dominant and 70 postural instability gait difficulty) and 70 healthy controls were obtained from the Parkinson's Progression Markers Initiative (PPMI) database, acquired using different scanner models. Radiomic features were extracted from 16 brain regions using various preprocessing pipelines. ComBat harmonization was applied using a combined batch variable incorporating both scanner models and preprocessing methods. Intraclass correlation coefficients (ICC) and Kruskal-Wallis tests assessed feature reproducibility before and after harmonization. Feature selection was performed using Linear Support Vector Classifier with L1 regularization. Support vector machine classifiers were used for PD subtype classification. ComBat harmonization significantly improved feature reproducibility across all feature groups. The percentage of features showing excellent robustness (ICC ≥ 0.90) increased from 40.2 to 56.3% after harmonization. First-order statistic features showed the highest robustness, with 71.11% demonstrating excellent ICC after harmonization. The proportion of features significantly affected by preprocessing methods was reduced following harmonization. Classification accuracy improved dramatically, from a range of 34-75% before harmonization to 89-96% after harmonization across all preprocessing methods. AUC values similarly increased from 0.28-0.87 to 0.95-0.99 after harmonization. ComBat harmonization significantly enhanced the reproducibility of radiomic features across preprocessing methods and improved PD motor subtype classification performance. This study highlights the importance of harmonization in radiomics research for PD and suggests potential clinical applications in personalized treatment planning.

Rest Tremor in Parkinson's Disease Is Associated with Ipsilateral Striatal Dopamine Transporter Binding

BACKGROUND

The cardinal motor symptoms of Parkinson's disease (PD) include rigidity, bradykinesia, and rest tremor. Rigidity and bradykinesia correlate with contralateral nigrostriatal degeneration and striatal dopamine deficit, but association between striatal dopamine function and rest tremor has remained unclear.

OBJECTIVE

The aim of this study was to investigate the possible link between dopamine function and rest tremor using Parkinson's Progression Markers Initiative dataset, the largest prospective neuroimaging cohort of patients with PD.

METHODS

Clinical, [123I]N-ω-fluoropropyl-2β-carbomethoxy-3β-(4-iodophenyl)nortropane ([123I]FP-CIT) single photon emission computed tomography (SPECT), and structural magnetic resonance imaging data from 354 early PD patients and 166 healthy controls were included in this study. We employed a novel approach allowing nonlinear registration of individual scans accurately to a standard space and voxelwise analyses of the association between motor symptoms and striatal dopamine transporter (DAT) binding.

RESULTS

Severity of both rigidity and bradykinesia was negatively associated with contralateral striatal DAT binding (PFWE < 0.05 [FWE, family-wise error corrected]). However, rest tremor amplitude was positively associated with increased ipsilateral DAT binding (PFWE < 0.05). The association between rest tremor and binding remained the same controlling for Hoehn & Yahr stage, Movement Disorder Society Unified Parkinson's Disease Rating Scale (MDS-UPDRS) part III score, bradykinesia-rigidity score, or motor phenotype. The association between rest tremor and binding was independent of bradykinesia-rigidity and replicated using 2-year follow-up data (PFWE < 0.05).

CONCLUSION

In agreement with the existing literature, we did not find a consistent association between rest tremor and contralateral dopamine defect. However, our results demonstrate a link between rest tremor and increased or less decreased ipsilateral DAT binding. Our findings provide novel information about the association between dopaminergic function and parkinsonian rest tremor. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Prospective study of dual-phase 99mTc-MIBI SPECT/CT nomogram for differentiating non-small cell lung cancer from benign pulmonary lesions

OBJECTIVES

To establish and validate a technetium 99m sestamibi (99mTc-MIBI) single-photon emission computed tomography/computed tomography (SPECT/CT) nomogram for predicting non-small cell lung cancer (NSCLC). Comparing the diagnostic performance of early and delayed SPECT/CT nomogram, and compare the diagnostic performance of SPECT/CT radiomics models with single SPECT and CT radiomics models.

METHODS

This prospective study included 119 lesions (NSCLC: n = 92, benign pulmonary lesions: n = 27) from 103 patients (mean age: 59.68 ± 8.94 years). Patients underwent dual-phase 99mTc-MIBI SPECT/CT imaging. They were divided into the training (n = 83) and validation (n = 36) cohorts. Logistic regression, support vector machine, random forest, and light-gradient boosting machine were applied to train and determine the optimal machine learning model. Then, combining radiomics score and clinical factors, establish nomograms for diagnosing NSCLC.

RESULT

CYFRA21-1 was selected for constructing the clinical model. In early imaging, the areas under the curve (AUCs) of the clinical model, radiomics model, and nomogram were 0.571, 0.830, and 0.875, respectively. The nomogram performed better than the clinical model and similarly to the radiomics model (P=0.020, P=0.216), and there are no statistically significant differences in the predictive performance between the radiomics model and the clinical model (P=0.103). In delayed imaging, the AUC was 0.643, 0.888, and 0.893, respectively. The predictive performance of the nomogram was superior compared to the clinical model and comparable to the radiomics model (P=0.042, P=0.480), and the radiomics model also demonstrated superior diagnostic performance compared to the clinical model (P=0.049). Compared to early SPECT/CT results, the AUC values of the nomogram and radiomics models in the delayed phase were higher, although no statistical differences were found (P=0.831, P=0.568). In delayed imaging, the AUC of the radiomics models for CT and SPECT was 0.696 and 0.768, respectively, SPECT/CT radiomics exhibited significant differences compared with CT and SPECT alone (P=0.042, P=0.038).

CONCLUSION

Dual-phase 99mTc-MIBI SPECT/CT nomograms and radiomics models can effectively predict NSCLC, providing an economically and non-invasive imaging method for diagnosing NSCLC, moreover, these findings provide a basis for early diagnosis and treatment strategies in NSCLC patients. Delayed-phase SPECT/CT imaging may offer greater practical value than early-phase imaging for diagnosing NSCLC. However, this novel approach necessitates further validation in larger, multi-center cohorts.

CLINICAL RELEVANCE

Radiomics nomogram based on SPECT/CT for discriminating NSCLC from benign lung lesions helps to aid early diagnosis and guide treatment.

KEY POINTS

Nomograms, based on dual-phase SPECT/CT, was constructed to discriminate between non-small cell lung cancer and benign lesions. SPECT/CT radiomics model has better predictive performance than SPECT and CT radiomics model.

Multi-modality radiomics of conventional T1 weighted and diffusion tensor imaging for differentiating Parkinson's disease motor subtypes in early-stages

This study aimed to develop and validate a multi-modality radiomics approach using T1-weighted and diffusion tensor imaging (DTI) to differentiate Parkinson's disease (PD) motor subtypes, specifically tremor-dominant (TD) and postural instability gait difficulty (PIGD), in early disease stages. We analyzed T1-weighted and DTI scans from 140 early-stage PD patients (70 TD, 70 PIGD) and 70 healthy controls from the Parkinson's Progression Markers Initiative database. Radiomics features were extracted from 16 brain regions of interest. After harmonization and feature selection, four machine learning classifiers were trained and evaluated for both three-class (HC vs TD vs PIGD) and binary (TD vs PIGD) classification tasks. The light gradient boosting machine (LGBM) classifier demonstrated the best overall performance. For the three-class classification, LGBM achieved an accuracy of 85% and an area under the receiver operating characteristic curve (AUC) of 0.94 using combined T1 and DTI features. In the binary classification task, LGBM reached an accuracy of 95% and AUC of 0.95. Key discriminative features were identified in the Thalamus, Amygdala, Hippocampus, and Substantia Nigra for the three-group classification, and in the Pallidum, Amygdala, Hippocampus, and Accumbens for binary classification. The combined T1 + DTI approach consistently outperformed single-modality classifications, with DTI alone showing particularly low performance (AUC 0.55-0.62) in binary classification. The high accuracy and AUC values suggest that this approach could significantly improve early diagnosis and subtyping of PD. These findings have important implications for clinical management, potentially enabling more personalized treatment strategies based on early, accurate subtype identification.

Automatic substantia nigra segmentation with Swin-Unet in susceptibility- and T2-weighted imaging: application to Parkinson disease diagnosis

Background

Accurately distinguishing between Parkinson disease (PD) and healthy controls (HCs) through reliable imaging method is crucial for appropriate therapeutic intervention. However, PD diagnosis is hindered by the subjective nature of the evaluation. We aimed to develop an automatic deep-learning method that can segment the substantia nigra areas on susceptibility-weighted imaging (SWI) and T2-weighted imaging (T2WI) and further differentiate patients with PD from HCs using a machine learning algorithm.

Methods

Magnetic resonance imaging (MRI) data from 83 patients with PD and 83 age- and sex-matched HCs were obtained on the same 3.0-T MRI scanner. A deep learning method with Swin-Unet was developed to segment volumes of interest (VOIs) on SWI and then map the VOIs on SWI to the corresponding T2WI; features were then extracted from the VOIs on SWI and T2WI. Three machine learning models were developed and compared to differentiate those with PD from HCs.

Results

Swin-Unet achieved a better Dice coefficient than did U-Net in SWI segmentation (0.832 vs. 0.712). Machine learning models outperformed visual analysis (P>0.05), and logistic regression (LR) achieved the best performance [area under the curve (AUC) ≥0.819] and the most stable (relative standard deviations in AUC ≤0.05). The test results showed that the AUC of the LR model based on SWI segmentation was 0.894 while that of the LR model based on T2WI segmentation was 0.876. There was no significant difference in VOIs based on manual labeling or automatic segmentation across T2WI, SWI, or a combination of the two (P>0.05). The AUCs of the LR model based on automatic segmentation were close to those of the model based on manual labeling (P>0.05).

Conclusions

Our approach could provide a powerful and useful method for automatically and rapidly diagnosing PD in the clinic with only T2WI.

Voxel-based texture similarity networks reveal individual variability and correlate with biological ontologies

The human brain is organized as a complex, hierarchical network. However, the structural covariance patterns among brain regions and the underlying biological substrates of such covariance networks remain to be clarified. The present study proposed a novel individualized structural covariance network termed voxel-based texture similarity networks (vTSNs) based on 76 refined voxel-based textural features derived from structural magnetic resonance images. Validated in three independent longitudinal healthy cohorts (40, 23, and 60 healthy participants, respectively) with two common brain atlases, we found that the vTSN could robustly resolve inter-subject variability with high test-retest reliability. In contrast to the regional-based texture similarity networks (rTSNs) that calculate radiomic features based on region-of-interest information, vTSNs had higher inter- and intra-subject variability ratios and test-retest reliability in connectivity strength and network topological properties. Moreover, the Spearman correlation indicated a stronger association of the gene expression similarity network (GESN) with vTSNs than with rTSNs (vTSN: r = 0.600, rTSN: r = 0.433, z = 39.784, P < 0.001). Hierarchical clustering identified 3 vTSN subnets with differential association patterns with 13 coexpression modules, 16 neurotransmitters, 7 electrophysiology, 4 metabolism, and 2 large-scale structural and 4 functional organization maps. Moreover, these subnets had unique biological hierarchical organization from the subcortex-limbic system to the ventral neocortex and then to the dorsal neocortex. Based on 424 unrelated, qualified healthy subjects from the Human Connectome Project, we found that vTSNs could sensitively represent sex differences, especially for connections in the subcortex-limbic system and between the subcortex-limbic system and the ventral neocortex. Moreover, a multivariate variance component model revealed that vTSNs could explain a significant proportion of inter-subject behavioral variance in cognition (80.0 %) and motor functions (63.4 %). Finally, using 494 healthy adults (aged 19-80 years old) from the Southwest University Adult Lifespan Dataset, the Spearman correlation identified a significant association between aging and vTSN strength, especially within the subcortex-limbic system and between the subcortex-limbic system and the dorsal neocortex. In summary, our proposed vTSN is robust in uncovering individual variability and neurobiological brain processes, which can serve as biologically plausible measures for linking biological processes and human behavior.

Machine learning for predicting cognitive decline within five years in Parkinson's disease: Comparing cognitive assessment scales with DAT SPECT and clinical biomarkers

OBJECTIVE

Parkinson's disease (PD) is an age-related neurodegenerative condition characterized mostly by motor symptoms. Although a wide range of non-motor symptoms (NMS) are frequently experienced by PD patients. One of the important and common NMS is cognitive impairment, which is measured using different cognitive scales. Monitoring cognitive impairment and its decline in PD is essential for patient care and management. In this study, our goal is to identify the most effective cognitive scale in predicting cognitive decline over a 5-year timeframe initializing clinical biomarkers and DAT SPECT.

METHODS

Machine Learning has previously shown superior performance in image and clinical data classification and detection. In this study, we propose to use machine learning with different types of data, such as DAT SPECT and clinical biomarkers, to predict PD-CD based on various cognitive scales. We collected 330 DAT SPECT images and their clinical data in baseline, years 2,3,4, and 5 from Parkinson's Progression Markers Initiative (PPMI). We then designed a 3D Autoencoder to extract deep radiomic features (DF) from DAT SPECT images, and we then concatenated it with 17 clinical features (CF) to predict cognitive decline based on Montreal Cognitive Assessment (MoCA) and The Movement Disorder Society-Unified Parkinson's Disease Rating Scale (MDS-UPDRS-I).

RESULTS

The utilization of MoCA as a cognitive decline scale yielded better performance in various years compared to MDS-UPDRS-I. In year 4, the application of the deep radiomic feature resulted in the highest achievement, with a cross-validation AUC of 89.28, utilizing the gradient boosting classifier. For the MDS-UPDRS-I scale, the highest achievement was obtained by utilizing the deep radiomic feature, resulting in a cross-validation AUC of 81.34 with the random forest classifier.

CONCLUSIONS

The study findings indicate that the MoCA scale may be a more effective predictor of cognitive decline within 5 years compared to MDS-UPDRS-I. Furthermore, deep radiomic features had better performance compared to sole clinical biomarkers or clinical and deep radiomic combined. These results suggest that using the MoCA score and deep radiomic features extracted from DAT SPECT could be a promising approach for identifying individuals at risk for cognitive decline in four years. Future research is needed to validate these findings and explore their utility in clinical practice.

Radiomics incorporating deep features for predicting Parkinson's disease in 123I-Ioflupane SPECT

PURPOSE

123I-Ioflupane SPECT is an effective tool for the diagnosis and progression assessment of Parkinson's disease (PD). Radiomics and deep learning (DL) can be used to track and analyze the underlying image texture and features to predict the Hoehn-Yahr stages (HYS) of PD. In this study, we aim to predict HYS at year 0 and year 4 after the first diagnosis with combined imaging, radiomics and DL-based features using 123I-Ioflupane SPECT images at year 0.

METHODS

In this study, 161 subjects from the Parkinson's Progressive Marker Initiative database underwent baseline 3T MRI and 123I-Ioflupane SPECT, with HYS assessment at years 0 and 4 after first diagnosis. Conventional imaging features (IF) and radiomic features (RaF) for striatum uptakes were extracted from SPECT images using MRI- and SPECT-based (SPECT-V and SPECT-T) segmentations respectively. A 2D DenseNet was used to predict HYS of PD, and simultaneously generate deep features (DF). The random forest algorithm was applied to develop models based on DF, RaF, IF and combined features to predict HYS (stage 0, 1 and 2) at year 0 and (stage 0, 1 and ≥ 2) at year 4, respectively. Model predictive accuracy and receiver operating characteristic (ROC) analysis were assessed for various prediction models.

RESULTS

For the diagnostic accuracy at year 0, DL (0.696) outperformed most models, except DF + IF in SPECT-V (0.704), significantly superior based on paired t-test. For year 4, accuracy of DF + RaF model in MRI-based method is the highest (0.835), significantly better than DF + IF, IF + RaF, RaF and IF models. And DL (0.820) surpassed models in both SPECT-based methods. The area under the ROC curve (AUC) highlighted DF + RaF model (0.854) in MRI-based method at year 0 and DF + RaF model (0.869) in SPECT-T method at year 4, outperforming DL models, respectively. And then, there was no significant differences between SPECT-based and MRI-based segmentation methods except for the imaging feature models.

CONCLUSION

The combination of radiomic and deep features enhances the prediction accuracy of PD HYS compared to only radiomics or DL. This suggests the potential for further advancements in predictive model performance for PD HYS at year 0 and year 4 after first diagnosis using 123I-Ioflupane SPECT images at year 0, thereby facilitating early diagnosis and treatment for PD patients. No significant difference was observed in radiomics results obtained between MRI- and SPECT-based striatum segmentations for radiomic and deep features.

Multimodal data fusion reveals functional and neurochemical correlates of Parkinson's disease

BACKGROUND

Neurotransmitter deficits and spatial associations among neurotransmitter distribution, brain activity, and clinical features in Parkinson's disease (PD) remain unclear. Better understanding of neurotransmitter impairments in PD may provide potential therapeutic targets. Therefore, we aimed to investigate the spatial relationship between PD-related patterns and neurotransmitter deficits.

METHODS

We included 59 patients with PD and 41 age- and sex-matched healthy controls (HCs). The voxel-wise mean amplitude of the low-frequency fluctuation (mALFF) was calculated and compared between the two groups. The JuSpace toolbox was used to test whether spatial patterns of mALFF alterations in patients with PD were associated with specific neurotransmitter receptor/transporter densities.

RESULTS

Compared to HCs, patients with PD showed reduced mALFF in the sensorimotor- and visual-related regions. In addition, mALFF alteration patterns were significantly associated with the spatial distribution of the serotonergic, dopaminergic, noradrenergic, glutamatergic, cannabinoid, and acetylcholinergic neurotransmitter systems (p < 0.05, false discovery rate-corrected).

CONCLUSIONS

Our results revealed abnormal brain activity patterns and specific neurotransmitter deficits in patients with PD, which may provide new insights into the mechanisms and potential targets for pharmacotherapy.

Robustness of radiomic features in 123I-ioflupane-dopamine transporter single-photon emission computer tomography scan

Radiomic features are usually used to predict target variables such as the absence or presence of a disease, treatment response, or time to symptom progression. One of the potential clinical applications is in patients with Parkinson's disease. Robust radiomic features for this specific imaging method have not yet been identified, which is necessary for proper feature selection. Thus, we are assessing the robustness of radiomic features in dopamine transporter imaging (DaT). For this study, we made an anthropomorphic head phantom with tissue heterogeneity using a personal 3D printer (polylactide 82% infill); the bone was subsequently reproduced with plaster. A surgical cotton ball with radiotracer (123I-ioflupane) was inserted. Scans were performed on the two-detector hybrid camera with acquisition parameters corresponding to international guidelines for DaT single photon emission tomography (SPECT). Reconstruction of SPECT was performed on a clinical workstation with iterative algorithms. Open-source LifeX software was used to extract 134 radiomic features. Statistical analysis was made in RStudio using the intraclass correlation coefficient (ICC) and coefficient of variation (COV). Overall, radiomic features in different reconstruction parameters showed a moderate reproducibility rate (ICC = 0.636, p <0.01). Assessment of ICC and COV within CT attenuation correction (CTAC) and non-attenuation correction (NAC) groups and within particular feature classes showed an excellent reproducibility rate (ICC > 0.9, p < 0.01), except for an intensity-based NAC group, where radiomic features showed a good repeatability rate (ICC = 0.893, p <0.01). By our results, CTAC becomes the main threat to feature stability. However, many radiomic features were sensitive to the selected reconstruction algorithm irrespectively to the attenuation correction. Radiomic features extracted from DaT-SPECT showed moderate to excellent reproducibility rates. These results make them suitable for clinical practice and human studies, but awareness of feature selection should be held, as some radiomic features are more robust than others.

Neuroimaging Biomarkers in Parkinson's Disease

Idiopathic Parkinson's disease (PD) is a neurodegenerative disorder that affects multiple systems in the body and is characterized by a variety of motor and non-motor (e.g., psychiatric, autonomic) symptoms. As the fastest growing neurological disorder expected to affect over 12 million people globally by 2040 (Dorsey, Bloem JAMA Neurol 75(1):9-10. https://doi.org/10.1001/jamaneurol.2017.3299 . PMID: 29131880, 2018), PD poses an enormous individual and public health burden. Currently, there are no therapies that can slow down the disease progression in PD, and existing therapies are limited to symptomatic treatment. Importantly, people in the prodromal phase who are at high risk of developing PD can now be identified, which makes disease prevention an achievable goal. An in-depth understanding of the pathological processes in PD is crucial for prevention and treatment development. Advanced multimodal neuroimaging techniques provide unique biomarkers that can further our understanding of PD at multiple levels ranging from neurotransmitters to neural networks. These neuroimaging biomarkers also have value in clinical application, for example, in the differential diagnosis of PD. As the field continues to advance, neuroimaging biomarkers are expected to become more specific, more widely accessible, and can be readily incorporated into translational research for treatment development in PD.

Radiotherapy modification based on artificial intelligence and radiomics applied to (18F)-fluorodeoxyglucose positron emission tomography/computed tomography

Over the last decades, the refinement of radiation therapy techniques has been associated with an increasing interest for individualized radiation therapy with the aim of increasing or maintaining tumor control and reducing radiation toxicity. Developments in artificial intelligence (AI), particularly machine learning and deep learning, in imaging sciences, including nuclear medecine, have led to significant enthusiasm for the concept of "rapid learning health system". AI combined with radiomics applied to (18F)-fluorodeoxyglucose positron emission tomography/computed tomography ([18F]-FDG PET/CT) offers a unique opportunity for the development of predictive models that can help stratify each patient's risk and guide treatment decisions for optimal outcomes and quality of life of patients treated with radiation therapy. Here we present an overview of the current contribution of AI and radiomics-based machine learning models applied to (18F)-FDG PET/CT in the management of cancer treated by radiation therapy.

Standardized Radiomics Analysis of Clinical Myocardial Perfusion Stress SPECT Images to Identify Coronary Artery Calcification

PURPOSE

Myocardial perfusion (MP) stress single-photon emission computed tomography (SPECT) is an established diagnostic test for patients suspected of coronary artery disease (CAD). Meanwhile, coronary artery calcification (CAC) scoring obtained from diagnostic CT is a highly sensitive test, offering incremental diagnostic information in identifying patients with significant CAD yet normal MP stress SPECT (MPSS) scans. However, after decades of wide utilization of MPSS, CAC is not commonly reimbursed (e.g. by the CMS), nor widely deployed in community settings. We studied the potential of complementary information deduced from the radiomics analysis of normal MPSS scans in predicting the CAC score.

METHODS

We collected data from 428 patients with normal (non-ischemic) MPSS (99mTc-sestamibi; consensus reading). A nuclear medicine physician verified iteratively reconstructed images (attenuation-corrected) to be free from fixed perfusion defects and artifactual attenuation. Three-dimensional images were automatically segmented into four regions of interest (ROIs), including myocardium and three vascular segments (left anterior descending [LAD]-left circumference [LCX]-right coronary artery [RCA]). We used our software package, standardized environment for radiomics analysis (SERA), to extract 487 radiomic features in compliance with the image biomarker standardization initiative (IBSI). Isotropic cubic voxels were discretized using fixed bin-number discretization (eight schemes). We first performed blind-to-outcome feature selection focusing on a priori usefulness, dynamic range, and redundancy of features. Subsequently, we performed univariate and multivariate machine learning analyses to predict CAC scores from i) selected radiomic features, ii) 10 clinical features, and iii) combined radiomics + clinical features. Univariate analysis invoked Spearman correlation with Benjamini-Hotchberg false-discovery correction. The multivariate analysis incorporated stepwise linear regression, where we randomly selected a 15% test set and divided the other 85% of data into 70% training and 30% validation sets. Training started from a constant (intercept) model, iteratively adding/removing features (stepwise regression), invoking the Akaike information criterion (AIC) to discourage overfitting. Validation was run similarly, except that the training output model was used as the initial model. We randomized training/validation sets 20 times, selecting the best model using log-likelihood for evaluation in the test set. Assessment in the test set was performed thoroughly by running the entire operation 50 times, subsequently employing Fisher's method to verify the significance of independent tests.

RESULTS

Unsupervised feature selection significantly reduced 8×487 features to 56. In univariate analysis, no feature survived the false-discovery rate (FDR) to directly correlate with CAC scores. Applying Fisher's method to the multivariate regression results demonstrated combining radiomics with the clinical features to enhance the significance of the prediction model across all cardiac segments.  Conclusions: Our standardized and statistically robust multivariate analysis demonstrated significant prediction of the CAC score for all cardiac segments when combining MPSS radiomic features with clinical features, suggesting radiomics analysis can add diagnostic or prognostic value to standard MPSS for wide clinical usage.

Prediction of Cognitive Decline in Parkinson's Disease Using Clinical and DAT SPECT Imaging Features, and Hybrid Machine Learning Systems

BACKGROUND

We aimed to predict Montreal Cognitive Assessment (MoCA) scores in Parkinson's disease patients at year 4 using handcrafted radiomics (RF), deep (DF), and clinical (CF) features at year 0 (baseline) applied to hybrid machine learning systems (HMLSs).

METHODS

297 patients were selected from the Parkinson's Progressive Marker Initiative (PPMI) database. The standardized SERA radiomics software and a 3D encoder were employed to extract RFs and DFs from single-photon emission computed tomography (DAT-SPECT) images, respectively. The patients with MoCA scores over 26 were indicated as normal; otherwise, scores under 26 were indicated as abnormal. Moreover, we applied different combinations of feature sets to HMLSs, including the Analysis of Variance (ANOVA) feature selection, which was linked with eight classifiers, including Multi-Layer Perceptron (MLP), K-Neighbors Classifier (KNN), Extra Trees Classifier (ETC), and others. We employed 80% of the patients to select the best model in a 5-fold cross-validation process, and the remaining 20% were employed for hold-out testing.

RESULTS

For the sole usage of RFs and DFs, ANOVA and MLP resulted in averaged accuracies of 59 ± 3% and 65 ± 4% for 5-fold cross-validation, respectively, with hold-out testing accuracies of 59 ± 1% and 56 ± 2%, respectively. For sole CFs, a higher performance of 77 ± 8% for 5-fold cross-validation and a hold-out testing performance of 82 + 2% were obtained from ANOVA and ETC. RF+DF obtained a performance of 64 ± 7%, with a hold-out testing performance of 59 ± 2% through ANOVA and XGBC. Usage of CF+RF, CF+DF, and RF+DF+CF enabled the highest averaged accuracies of 78 ± 7%, 78 ± 9%, and 76 ± 8% for 5-fold cross-validation, and hold-out testing accuracies of 81 ± 2%, 82 ± 2%, and 83 ± 4%, respectively.

CONCLUSIONS

We demonstrated that CFs vitally contribute to predictive performance, and combining them with appropriate imaging features and HMLSs can result in the best prediction performance.

Myocardial Perfusion SPECT Imaging Radiomic Features and Machine Learning Algorithms for Cardiac Contractile Pattern Recognition

A U-shaped contraction pattern was shown to be associated with a better Cardiac resynchronization therapy (CRT) response. The main goal of this study is to automatically recognize left ventricular contractile patterns using machine learning algorithms trained on conventional quantitative features (ConQuaFea) and radiomic features extracted from Gated single-photon emission computed tomography myocardial perfusion imaging (GSPECT MPI). Among 98 patients with standard resting GSPECT MPI included in this study, 29 received CRT therapy and 69 did not (also had CRT inclusion criteria but did not receive treatment yet at the time of data collection, or refused treatment). A total of 69 non-CRT patients were employed for training, and the 29 were employed for testing. The models were built utilizing features from three distinct feature sets (ConQuaFea, radiomics, and ConQuaFea + radiomics (combined)), which were chosen using Recursive feature elimination (RFE) feature selection (FS), and then trained using seven different machine learning (ML) classifiers. In addition, CRT outcome prediction was assessed by different treatment inclusion criteria as the study's final phase. The MLP classifier had the highest performance among ConQuaFea models (AUC, SEN, SPE = 0.80, 0.85, 0.76). RF achieved the best performance in terms of AUC, SEN, and SPE with values of 0.65, 0.62, and 0.68, respectively, among radiomic models. GB and RF approaches achieved the best AUC, SEN, and SPE values of 0.78, 0.92, and 0.63 and 0.74, 0.93, and 0.56, respectively, among the combined models. A promising outcome was obtained when using radiomic and ConQuaFea from GSPECT MPI to detect left ventricular contractile patterns by machine learning.

Radiomics insight into the neurodegenerative "hot" brain: A narrative review from the nuclear medicine perspective

The application of radiomics for non-oncologic diseases is currently emerging. Despite its relative infancy state, the evidence highlights the potential of radiomics approaches to serve as neuroimaging biomarkers in the field of the neurodegenerative brain. This systematic review presents the last progress and potential application of radiomics in the field of neurodegenerative nuclear imaging applied to positron-emission tomography (PET) and single-photon emission computed tomography (SPECT) by focusing mainly on the two most common neurodegenerative disorders, Alzheimer's (AD) and Parkinson's disease (PD). A comprehensive review of the current literature was performed using the PubMed and Web of Science databases up to November 2022. The final collection of eighteen relevant publications was grouped as AD-related and PD-related. The main efforts in the field of AD dealt with radiomics-based early diagnosis of preclinical AD and the prediction of MCI to AD conversion, meanwhile, in the setting of PD, the radiomics techniques have been used in the attempt to improve the assessment of PD diagnosis, the differential diagnosis between PD and other parkinsonism, severity assessment, and outcome prediction. Although limited evidence with relatively small cohort studies, it seems that radiomics-based analysis using nuclear medicine tools, mainly [18F]Fluorodeoxyglucose (FDG) and β-amyloid (Aβ) PET, and dopamine transporter (DAT) SPECT, can be used for computer-aided diagnoses in AD-continuum and parkinsonian disorders. Combining nuclear radiomics analysis with clinical factors and introducing a multimodality approach can significantly improve classification and prediction efficiency in neurodegenerative disorders.

Relevance of 18F-DOPA visual and semi-quantitative PET metrics for the diagnostic of Parkinson disease in clinical practice: a machine learning-based inference study

PURPOSE

To decipher the relevance of visual and semi-quantitative 6-fluoro-(18F)-L-DOPA (¹⁸F-DOPA) interpretation methods for the diagnostic of idiopathic Parkinson disease (IPD) in hybrid positron emission tomography (PET) and magnetic resonance imaging.

MATERIAL AND METHODS

A total of 110 consecutive patients (48 IPD and 62 controls) with 11 months of median clinical follow-up (reference standard) were included. A composite visual assessment from five independent nuclear imaging readers, together with striatal standard uptake value (SUV) to occipital SUV ratio, striatal gradients and putamen asymmetry-based semi-quantitative PET metrics automatically extracted used to train machine learning models to classify IPD versus controls. Using a ratio of 70/30 for training and testing sets, respectively, five classification models-k-NN, LogRegression, support vector machine, random forest and gradient boosting-were trained by using 100 times repeated nested cross-validation procedures. From the best model on average, the contribution of PET parameters was deciphered using the Shapley additive explanations method (SHAP). Cross-validated receiver operating characteristic curves (cv-ROC) of the most contributive PET parameters were finally estimated and compared.

RESULTS

The best machine learning model (k-NN) provided final cv-ROC of 0.81. According to SHAP analyses, visual PET metric was the most important contributor to the model overall performance, followed by the minimum between left and right striatal to occipital SUV ratio. The 10-time cv-ROC curves of visual, min SUVr or both showed quite similar performance (mean area under the ROC of 0.81, 0.81 and 0.79, respectively, for visual, min SUVr or both).

CONCLUSION

Visual expert analysis remains the most relevant parameter to predict IPD diagnosis at 11 months of median clinical follow-up in ¹⁸F-FDOPA. The min SUV ratio appears interesting in the perspective of simple semi-automated diagnostic workflows.

Introduction to radiomics for a clinical audience

Radiomics is a rapidly developing field of research focused on the extraction of quantitative features from medical images, thus converting these digital images into minable, high-dimensional data, which offer unique biological information that can enhance our understanding of disease processes and provide clinical decision support. To date, most radiomics research has been focused on oncological applications; however, it is increasingly being used in a raft of other diseases. This review gives an overview of radiomics for a clinical audience, including the radiomics pipeline and the common pitfalls associated with each stage. Key studies in oncology are presented with a focus on both those that use radiomics analysis alone and those that integrate its use with other multimodal data streams. Importantly, clinical applications outside oncology are also presented. Finally, we conclude by offering a vision for radiomics research in the future, including how it might impact our practice as radiologists.

Application of novel hybrid machine learning systems and radiomics features for non-motor outcome prediction in Parkinson's disease

Objectives.Parkinson's disease (PD) is a complex neurodegenerative disorder, affecting 2%-3% of the elderly population. Montreal Cognitive Assessment (MoCA), a rapid nonmotor screening test, assesses different cognitive dysfunctionality aspects. Early MoCA prediction may facilitate better temporal therapy and disease control. Radiomics features (RF), in addition to clinical features (CF), are indicated to increase clinical diagnoses, etc, bridging between medical imaging procedures and personalized medicine. We investigate the effect of RFs, CFs, and conventional imaging features (CIF) to enhance prediction performance using hybrid machine learning systems (HMLS).Methods.We selected 210 patients with 981 features (CFs, CIFs, and RFs) from the Parkinson's Progression-Markers-Initiative database. We generated 4 datasets, namely using (i), (ii) year-0 (D1) or year-1 (D2) features, (iii) longitudinal data (D3, putting datasets in years 0 and 1 longitudinally next to each other), and (iv) timeless data (D4, effectively doubling dataset size by listing both datasets from years 0 and 1 separately). First, we directly applied 23 predictor algorithms (PA) to the datasets to predict year-4 MoCA, which PD patients this year have a higher dementia risk. Subsequently, HMLSs, including 14 attribute extraction and 10 feature selection algorithms followed by PAs were employed to enhance prediction performances. 80% of all datapoints were utilized to select the best model based on minimum mean absolute error (MAE) resulting from 5-fold cross-validation. Subsequently, the remaining 20% was used for hold-out testing of the selected models.Results.When applying PAs without ASAs/FEAs to datasets (MoCA outcome range: [11,30]), Adaboost achieved an MAE of 1.74 ± 0.29 on D4 with a hold-out testing performance of 1.71. When employing HMLSs, D4 + Minimum_Redundancy_Maximum_Relevance (MRMR)+K_Nearest_Neighbor Regressor achieved the highest performance of 1.05 ± 0.25 with a hold-out testing performance of 0.57.Conclusion.Our study shows the importance of using larger datasets (timeless), and utilizing optimized HMLSs, for significantly improved prediction of MoCA in PD patients.

Joint EANM/SNMMI guideline on radiomics in nuclear medicine : Jointly supported by the EANM Physics Committee and the SNMMI Physics, Instrumentation and Data Sciences Council

PURPOSE

The purpose of this guideline is to provide comprehensive information on best practices for robust radiomics analyses for both hand-crafted and deep learning-based approaches.

METHODS

In a cooperative effort between the EANM and SNMMI, we agreed upon current best practices and recommendations for relevant aspects of radiomics analyses, including study design, quality assurance, data collection, impact of acquisition and reconstruction, detection and segmentation, feature standardization and implementation, as well as appropriate modelling schemes, model evaluation, and interpretation. We also offer an outlook for future perspectives.

CONCLUSION

Radiomics is a very quickly evolving field of research. The present guideline focused on established findings as well as recommendations based on the state of the art. Though this guideline recognizes both hand-crafted and deep learning-based radiomics approaches, it primarily focuses on the former as this field is more mature. This guideline will be updated once more studies and results have contributed to improved consensus regarding the application of deep learning methods for radiomics. Although methodological recommendations in the present document are valid for most medical image modalities, we focus here on nuclear medicine, and specific recommendations when necessary are made for PET/CT, PET/MR, and quantitative SPECT.

Predicting hepatocellular carcinoma early recurrence after ablation based on magnetic resonance imaging radiomics nomogram

BACKGROUND

The aim of this study is to investigate a model for predicting the early recurrence of hepatocellular carcinoma (HCC) after ablation.

METHODS

A total of 181 patients with HCC after ablation (train group was 119 cases; validation group was 62 cases) were enrolled. The cases of early recurrence in the set of train and validation were 63 and 31, respectively. Radiomics features were extracted from the enhanced magnetic resonance imaging scanning, including pre-contrast injection, arterial phase, late arterial phase, portal venous phase, and delayed phase. The least absolute shrinkage and selection operator cox proportional hazards regression after univariate and multivariate analysis was used to screen radiomics features and build integrated models. The nomograms predicting recurrence and survival of patients of HCC after ablation were established based on the clinical, imaging, and radiomics features. The area under the curve (AUC) of the receiver operating characteristic curve and C-index for the train and validation group was used to evaluate model efficacy.

RESULTS

Four radiomics features were selected out of 34 texture features to formulate the rad-score. Multivariate analyses suggested that the rad-score, number of lesions, integrity of the capsule, pathological type, and alpha-fetoprotein were independent influencing factors. The AUC of predicting early recurrence at 1, 2, and 3 years in the train group was 0.79 (95% CI: 0.72-0.88), 0.72 (95% CI: 0.63-0.82), and 0.71 (95% CI: 0.61-0.83), respectively. The AUC of predicting early recurrence at 1, 2, and 3 years in the validation group was 0.72 (95% CI: 0.58-0.84), 0.61 (95% CI: 0.45-0.78) and 0.64 (95% CI: 0.40-0.87).

CONCLUSION

The model for early recurrence of HCC after ablation based on the clinical, imaging, and radiomics features presented good predictive performance. This may facilitate the early treatment of patients.

Use of deep learning-based radiomics to differentiate Parkinson's disease patients from normal controls: a study based on [18F]FDG PET imaging

OBJECTIVES

We proposed a novel deep learning-based radiomics (DLR) model to diagnose Parkinson's disease (PD) based on [¹⁸F]fluorodeoxyglucose (FDG) PET images.

METHODS

In this two-center study, 255 normal controls (NCs) and 103 PD patients were enrolled from Huashan Hospital, China; 26 NCs and 22 PD patients were enrolled as a separate test group from Wuxi 904 Hospital, China. The proposed DLR model consisted of a convolutional neural network-based feature encoder and a support vector machine (SVM) model-based classifier. The DLR model was trained and validated in the Huashan cohort and tested in the Wuxi cohort, and accuracy, sensitivity, specificity and receiver operator characteristic (ROC) curve graphs were used to describe the model's performance. Comparative experiments were performed based on four other models including the scale model, radiomics model, standard uptake value ratio (SUVR) model and DLR model.

RESULTS

The DLR model demonstrated superiority in differentiating PD patients and NCs in comparison to other models, with an accuracy of 95.17% [90.35%, 98.13%] (95% confidence intervals, CI) in the Huashan cohort. Moreover, the DLR model also demonstrated greater performance in diagnosing PD early than routine methods, with an accuracy of 85.58% [78.60%, 91.57%] in the Huashan cohort.

CONCLUSIONS

We developed a DLR model based on [¹⁸F]FDG PET images that showed good performance in the noninvasive, individualized prediction of PD and was superior to traditional handcrafted methods. This model has the potential to guide and facilitate clinical diagnosis and contribute to the development of precision treatment.

KEY POINTS

The DLR method on [¹⁸F]FDG PET images helps clinicians to diagnose PD and PD subgroups from normal controls. A prospective two-center study showed that the DLR method provides greater diagnostic accuracy.

Neuromelanin and T2*-MRI for the assessment of genetically at-risk, prodromal, and symptomatic Parkinson's disease

MRI was suggested as a promising method for the diagnosis and assessment of Parkinson's Disease (PD). We aimed to assess the sensitivity of neuromelanin-MRI and T₂* with radiomics analysis for detecting PD, identifying individuals at risk, and evaluating genotype-related differences. Patients with PD and non-manifesting (NM) participants [NM-carriers (NMC) and NM-non-carriers (NMNC)], underwent MRI and DAT-SPECT. Imaging-based metrics included 48 neuromelanin and T₂* radiomics features and DAT-SPECT specific-binding-ratios (SBR), were extracted from several brain regions. Imaging values were assessed for their correlations with age, differences between groups, and correlations with the MDS-likelihood-ratio (LR) score. Several machine learning classifiers were evaluated for group classification. A total of 127 participants were included: 46 patients with PD (62.3 ± 10.0 years) [15:LRRK2-PD, 16:GBA-PD, and 15:idiopathic-PD (iPD)], 47 NMC (51.5 ± 8.3 years) [24:LRRK2-NMC and 23:GBA-NMC], and 34 NMNC (53.5 ± 10.6 years). No significant correlations were detected between imaging parameters and age. Thirteen MRI-based parameters and radiomics features demonstrated significant differences between PD and NMNC groups. Support-Vector-Machine (SVM) classifier achieved the highest performance (AUC = 0.77). Significant correlations were detected between LR scores and two radiomic features. The classifier successfully identified two out of three NMC who converted to PD. Genotype-related differences were detected based on radiomic features. SBR values showed high sensitivity in all analyses. In conclusion, neuromelanin and T₂* MRI demonstrated differences between groups and can be used for the assessment of individuals at-risk in cases when DAT-SPECT can't be performed. Combining neuromelanin and T₂*-MRI provides insights into the pathophysiology underlying PD, and suggests that iron accumulation precedes neuromelanin depletion during the prodromal phase.

A Pipeline for the Implementation and Visualization of Explainable Machine Learning for Medical Imaging Using Radiomics Features

Machine learning (ML) models have been shown to predict the presence of clinical factors from medical imaging with remarkable accuracy. However, these complex models can be difficult to interpret and are often criticized as “black boxes”. Prediction models that provide no insight into how their predictions are obtained are difficult to trust for making important clinical decisions, such as medical diagnoses or treatment. Explainable machine learning (XML) methods, such as Shapley values, have made it possible to explain the behavior of ML algorithms and to identify which predictors contribute most to a prediction. Incorporating XML methods into medical software tools has the potential to increase trust in ML-powered predictions and aid physicians in making medical decisions. Specifically, in the field of medical imaging analysis the most used methods for explaining deep learning-based model predictions are saliency maps that highlight important areas of an image. However, they do not provide a straightforward interpretation of which qualities of an image area are important. Here, we describe a novel pipeline for XML imaging that uses radiomics data and Shapley values as tools to explain outcome predictions from complex prediction models built with medical imaging with well-defined predictors. We present a visualization of XML imaging results in a clinician-focused dashboard that can be generalized to various settings. We demonstrate the use of this workflow for developing and explaining a prediction model using MRI data from glioma patients to predict a genetic mutation.

Dopamine transporter single-photon emission computed tomography-derived radiomics signature for detecting Parkinson's disease

BACKGROUND

We hypothesised that the radiomics signature, which includes texture information of dopamine transporter single-photon emission computed tomography (DAT-SPECT) images for Parkinson's disease (PD), may assist semi-quantitative indices. Herein, we constructed a radiomics signature using DAT-SPECT-derived radiomics features that effectively discriminated PD from healthy individuals and evaluated its classification performance.

RESULTS

We analysed 413 cases of both normal control (NC, n = 101) and PD (n = 312) groups from the Parkinson's Progression Markers Initiative database. Data were divided into the training and two test datasets with different SPECT manufacturers. DAT-SPECT images were spatially normalised to the Montreal Neurologic Institute space. We calculated 930 radiomics features, including intensity- and texture-based features in the caudate, putamen, and pallidum volumes of interest. The striatum uptake ratios (SURs) of the caudate, putamen, and pallidum were also calculated as conventional semi-quantification indices. The least absolute shrinkage and selection operator was used for feature selection and construction of the radiomics signature. The four classification models were constructed using a radiomics signature and/or semi-quantitative indicator. Furthermore, we compared the classification performance of the semi-quantitative indicator alone and the combination with the radiomics signature for the classification models. The receiver operating characteristics (ROC) analysis was used to evaluate the classification performance. The classification performance of SUR_putamen was higher than that of other semi-quantitative indicators. The radiomics signature resulted in a slightly increased area under the ROC curve (AUC) compared to SUR_putamen in each test dataset. When combined with SUR_putamen and radiomics signature, all classification models showed slightly higher AUCs than that of SUR_putamen alone.

CONCLUSION

We constructed a DAT-SPECT image-derived radiomics signature. Performance analysis showed that the current radiomics signature would be helpful for the diagnosis of PD and has the potential to provide robust diagnostic performance.

Development and validation of a radiomic model for the diagnosis of dopaminergic denervation on [18F]FDOPA PET/CT

Purpose

FDOPA PET shows good performance for the diagnosis of striatal dopaminergic denervation, making it a valuable tool for the differential diagnosis of Parkinsonism. Textural features are image biomarkers that could potentially improve the early diagnosis and monitoring of neurodegenerative parkinsonian syndromes. We explored the performances of textural features for binary classification of FDOPA scans.

Methods

We used two FDOPA PET datasets: 443 scans for feature selection, and 100 scans from a different PET/CT system for model testing. Scans were labelled according to expert interpretation (dopaminergic denervation versus no dopaminergic denervation). We built LASSO logistic regression models using 43 biomarkers including 32 textural features. Clinical data were also collected using a shortened UPDRS scale.

Results

The model built from the clinical data alone had a mean area under the receiver operating characteristics (AUROC) of 63.91. Conventional imaging features reached a maximum score of 93.47 but the addition of textural features significantly improved the AUROC to 95.73 (p < 0.001), and 96.10 (p < 0.001) when limiting the model to the top three features: GLCM_Correlation, Skewness and Compacity. Testing the model on the external dataset yielded an AUROC of 96.00, with 95% sensitivity and 97% specificity. GLCM_Correlation was one of the most independent features on correlation analysis, and systematically had the heaviest weight in the classification model.

Conclusion

A simple model with three radiomic features can identify pathologic FDOPA PET scans with excellent sensitivity and specificity. Textural features show promise for the diagnosis of parkinsonian syndromes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00259-022-05816-7.

The Disease Modification Conundrum in Parkinson’s Disease: Failures and Hopes

In the last half-century, Parkinson’s disease (PD) has played a historical role in demonstrating our ability to translate preclinical scientific advances in pathology and pharmacology into highly effective clinical therapies. Yet, as highly efficacious symptomatic treatments were successfully developed and adopted in clinical practice, PD remained a progressive disease without a cure. In contrast with the success story of symptomatic therapies, the lack of translation of disease-modifying interventions effective in preclinical models into clinical success has continued to accumulate failures in the past two decades. The ability to stop, prevent or mitigate progression in PD remains the “holy grail” in PD science at the present time. The large number of high-quality disease modification clinical trials in the past two decades with its lessons learned, as well as the growing knowledge of PD molecular pathology should enable us to have a deeper understanding of the reasons for past failures and what we need to do to reach better outcomes. Periodic reviews and mini-reviews of the unsolved disease modification conundrum in PD are important, considering how this field is rapidly evolving along with our views and understanding of the possible explanations.

Longitudinal clustering analysis and prediction of Parkinson's disease progression using radiomics and hybrid machine learning

BACKGROUND

We employed machine learning approaches to (I) determine distinct progression trajectories in Parkinson's disease (PD) (unsupervised clustering task), and (II) predict progression trajectories (supervised prediction task), from early (years 0 and 1) data, making use of clinical and imaging features.

METHODS

We studied PD-subjects derived from longitudinal datasets (years 0, 1, 2 & 4; Parkinson's Progressive Marker Initiative). We extracted and analyzed 981 features, including motor, non-motor, and radiomics features extracted for each region-of-interest (ROIs: left/right caudate and putamen) using our standardized standardized environment for radiomics analysis (SERA) radiomics software. Segmentation of ROIs on dopamine transposer - single photon emission computed tomography (DAT SPECT) images were performed via magnetic resonance images (MRI). After performing cross-sectional clustering on 885 subjects (original dataset) to identify disease subtypes, we identified optimal longitudinal trajectories using hybrid machine learning systems (HMLS), including principal component analysis (PCA) + K-Means algorithms (KMA) followed by Bayesian information criterion (BIC), Calinski-Harabatz criterion (CHC), and elbow criterion (EC). Subsequently, prediction of the identified trajectories from early year data was performed using multiple HMLSs including 16 Dimension Reduction Algorithms (DRA) and 10 classification algorithms.

RESULTS

We identified 3 distinct progression trajectories. Hotelling's t squared test (HTST) showed that the identified trajectories were distinct. The trajectories included those with (I, II) disease escalation (2 trajectories, 27% and 38% of patients) and (III) stable disease (1 trajectory, 35% of patients). For trajectory prediction from early year data, HMLSs including the stochastic neighbor embedding algorithm (SNEA, as a DRA) as well as locally linear embedding algorithm (LLEA, as a DRA), linked with the new probabilistic neural network classifier (NPNNC, as a classifier), resulted in accuracies of 78.4% and 79.2% respectively, while other HMLSs such as SNEA + Lib_SVM (library for support vector machines) and t_SNE (t-distributed stochastic neighbor embedding) + NPNNC resulted in 76.5% and 76.1% respectively.

CONCLUSIONS

This study moves beyond cross-sectional PD subtyping to clustering of longitudinal disease trajectories. We conclude that combining medical information with SPECT-based radiomics features, and optimal utilization of HMLSs, can identify distinct disease trajectories in PD patients, and enable effective prediction of disease trajectories from early year data.

Single-Photon Emission Computed Tomography/Computed Tomography Image-Based Radiomics for Discriminating Vertebral Bone Metastases From Benign Bone Lesions in Patients With Tumors

Purpose: The purpose of this study was to investigate the feasibility of Single-Photon Emission Computed Tomography/Computed Tomography (SPECT/CT) image-based radiomics in differentiating bone metastases from benign bone lesions in patients with tumors. Methods: A total of 192 lesions from 132 patients (134 in the training group, 58 in the validation group) diagnosed with vertebral bone metastases or benign bone lesions were enrolled. All images were evaluated and diagnosed independently by two physicians with more than 20 years of diagnostic experience for qualitative classification, the images were imported into MaZda software in Bitmap (BMP) format for feature extraction. All radiomics features were selected by least absolute shrinkage and selection operator (LASSO) regression and 10-fold cross-validation algorithms after the process of normalization and correlation analysis. Based on these selected features, two models were established: The CT model and SPECT model (radiomics features were derived from CT and SPECT images, respectively). In addition, a combination model (ComModel) combined CT and SPECT features was developed in order to better evaluate the predictive performance of radiomics models. Subsequently, the diagnostic performance between each model was separately evaluated by a confusion matrix. Results: There were 12, 13, and 18 features contained within the CT, SPECT, and ComModel, respectively. The constructed radiomics models based on SPECT/CT images to discriminate between bone metastases and benign bone lesions not only had high diagnostic efficacy in the training group (AUC of 0.894, 0.914, 0.951 for CT model, SPECT model, and ComModel, respectively), but also performed well in the validation group (AUC; 0.844, 0.871, 0.926). The AUC value of the human experts was 0.849 and 0.839 in the training and validation groups, respectively. Furthermore, both SPECT model and ComModel show higher classification performance than human experts in the training group (P = 0.021 and P = 0.001, respectively) and the validation group (P = 0.037 and P = 0.007, respectively). All models showed better diagnostic accuracy than human experts in the training group and the validation group. Conclusion: Radiomics derived from SPECT/CT images could effectively discriminate between bone metastases and benign bone lesions. This technique may be a new non-invasive way to help prevent unnecessary delays in diagnosis and a potential contribution in disease staging and treatment planning.

Cardiac SPECT radiomic features repeatability and reproducibility: A multi-scanner phantom study

BACKGROUND

The aim of this work was to assess the robustness of cardiac SPECT radiomic features against changes in imaging settings, including acquisition, and reconstruction parameters.

METHODS

Four commercial SPECT and SPECT/CT cameras were used to acquire images of a static cardiac phantom mimicking typical myorcardial perfusion imaging using 185 MBq of 99mTc. The effects of different image acquisition and reconstruction parameters, including number of views, view matrix size, attenuation correction, as well as image reconstruction related parameters (algorithm, number of iterations, number of subsets, type of post-reconstruction filter, and its associated parameters, including filter order and cut-off frequency) were studied. In total, 5,063 transverse views were reconstructed by varying the aforementioned factors. Eighty-seven radiomic features including first-, second-, and high-order textures were extracted from these images. To assess reproducibility and repeatability, the coefficient of variation (COV), as a widely adopted metric, was measured for each of the radiomic features over the different imaging settings.

RESULTS

The Inverse Difference Moment Normalized (IDMN) and Inverse Difference Normalized (IDN) features from the Gray Level Co-occurrence Matrix (GLCM), Run Percentage (RP) from the Gray Level Co-occurrence Matrix (GLRLM), Zone Entropy (ZE) from the Gray Level Size Zone Matrix (GLSZM), and Dependence Entropy (DE) from the Gray Level Dependence Matrix (GLDM) feature sets were the only features that exhibited high reproducibility (COV ≤ 5%) against changes in all imaging settings. In addition, Large Area Low Gray Level Emphasis (LALGLE), Small Area Low Gray Level Emphasis (SALGLE) and Low Gray Level Zone Emphasis (LGLZE) from GLSZM, and Small Dependence Low Gray Level Emphasis (SDLGLE) from GLDM feature sets turned out to be less reproducible (COV > 20%) against changes in imaging settings. The GLRLM (31.88%) and GLDM feature set (54.2%) had the highest (COV < 5%) and lowest (COV > 20%) number of the reproducible features, respectively. Matrix size had the largest impact on feature variability as most of the features were not repeatable when matrix size was modified with 82.8% of them having a COV > 20%.

CONCLUSION

The repeatability and reproducibility of SPECT/CT cardiac radiomic features under different imaging settings is feature-dependent. Different image acquisition and reconstruction protocols have variable effects on radiomic features. The radiomic features exhibiting low COV are potential candidates for future clinical studies.

Dopamine Transporter Imaging, Current Status of a Potential Biomarker: A Comprehensive Review

A major goal of current clinical research in Parkinson's disease (PD) is the validation and standardization of biomarkers enabling early diagnosis, predicting outcomes, understanding PD pathophysiology, and demonstrating target engagement in clinical trials. Molecular imaging with specific dopamine-related tracers offers a practical indirect imaging biomarker of PD, serving as a powerful tool to assess the status of presynaptic nigrostriatal terminals. In this review we provide an update on the dopamine transporter (DAT) imaging in PD and translate recent findings to potentially valuable clinical practice applications. The role of DAT imaging as diagnostic, preclinical and predictive biomarker is discussed, especially in view of recent evidence questioning the incontrovertible correlation between striatal DAT binding and nigral cell or axon counts.

Texture Analysis and its Applications in Biomedical Imaging: A Survey

Texture analysis describes a variety of image analysis techniques that quantify the variation in intensity and pattern. This paper provides an overview of several texture analysis approaches addressing the rationale supporting them, their advantages, drawbacks, and applications. This surveys emphasis is in collecting and categorising over five decades of active research on texture analysis. Brief descriptions of different approaches are presented along with application examples. From a broad range of texture analysis applications, this surveys final focus is on biomedical image analysis. An up-to-date list of biological tissues and organs in which disorders produce texture changes that may be used to spot disease onset and progression is provided. Finally, the role of texture analysis methods as biomarkers of disease is summarised.

Convolutional neural network based attenuation correction for123I-FP-CIT SPECT with focused striatum imaging

SPECT imaging with¹²³I-FP-CIT is used for diagnosis of neurodegenerative disorders like Parkinson's disease. Attenuation correction (AC) can be useful for quantitative analysis of¹²³I-FP-CIT SPECT. Ideally, AC would be performed based on attenuation maps (μ-maps) derived from perfectly registered CT scans. Suchμ-maps, however, are most times not available and possible errors in image registration can induce quantitative inaccuracies in AC corrected SPECT images. Earlier, we showed that a convolutional neural network (CNN) based approach allows to estimate SPECT-alignedμ-maps for full brain perfusion imaging using only emission data. Here we investigate the feasibility of similar CNN methods for axially focused¹²³I-FP-CIT scans. We tested our approach on a high-resolution multi-pinhole prototype clinical SPECT system in a Monte Carlo simulation study. Three CNNs that estimateμ-maps in a voxel-wise, patch-wise and image-wise manner were investigated. As the added value of AC on clinical¹²³I-FP-CIT scans is still debatable, the impact of AC was also reported to check in which cases CNN based AC could be beneficial. AC using the ground truthμ-maps (GT-AC) and CNN estimatedμ-maps (CNN-AC) were compared with the case when no AC was done (No-AC). Results show that the effect of using GT-AC versus CNN-AC or No-AC on striatal shape and symmetry is minimal. Specific binding ratios (SBRs) from localized regions show a deviation from GT-AC≤2.5% for all three CNN-ACs while No-AC systematically underestimates SBRs by 13.1%. A strong correlation (r≥0.99) was obtained between GT-AC based SBRs and SBRs from CNN-ACs and No-AC. Absolute quantification (in kBq ml^-1) shows a deviation from GT-AC within 2.2% for all three CNN-ACs and of 71.7% for No-AC. To conclude, all three CNNs show comparable performance in accurateμ-map estimation and¹²³I-FP-CIT quantification. CNN-estimatedμ-map can be a promising substitute for CT-basedμ-map.

An interpretable multiparametric radiomics model for the diagnosis of schizophrenia using magnetic resonance imaging of the corpus callosum

There is a growing need to develop novel strategies for the diagnosis of schizophrenia using neuroimaging biomarkers. We investigated the robustness of the diagnostic model for schizophrenia using radiomic features from T1-weighted and diffusion tensor images of the corpus callosum (CC). A total of 165 participants [86 schizophrenia and 79 healthy controls (HCs)] were allocated to training (N = 115) and test (N = 50) sets. Radiomic features of the CC subregions were extracted from T1-weighted, apparent diffusion coefficient (ADC), and fractional anisotropy (FA) images (N = 1605). Following feature selection, various combinations of classifiers were trained, and Bayesian optimization was adopted in the best performing classifier. Discrimination, calibration, and clinical utility of the model were assessed. An online calculator was constructed to offer the probability of having schizophrenia. SHapley Additive exPlanations (SHAP) was applied to explore the interpretability of the model. We identified 30 radiomic features to differentiate participants with schizophrenia from HCs. The Bayesian optimized model achieved the highest performance, with an area under the curve (AUC), accuracy, sensitivity, and specificity of 0.89 (95% confidence interval: 0.81-0.98), 80.0, 83.3, and 76.9%, respectively, in the test set. The final model offers clinical probability in an online calculator. The model explanation by SHAP suggested that second-order features from the posterior CC were highly associated with the risk of schizophrenia. The multiparametric radiomics model focusing on the CC shows its robustness for the diagnosis of schizophrenia. Radiomic features could be a potential source of biomarkers that support the biomarker-based diagnosis of schizophrenia and improve the understanding of its neurobiology.

Oropharyngeal cancer patient stratification using random forest based-learning over high-dimensional radiomic features

To improve risk prediction for oropharyngeal cancer (OPC) patients using cluster analysis on the radiomic features extracted from pre-treatment Computed Tomography (CT) scans. 553 OPC Patients randomly split into training (80%) and validation (20%), were classified into 2 or 3 risk groups by applying hierarchical clustering over the co-occurrence matrix obtained from a random survival forest (RSF) trained over 301 radiomic features. The cluster label was included together with other clinical data to train an ensemble model using five predictive models (Cox, random forest, RSF, logistic regression, and logistic-elastic net). Ensemble performance was evaluated over the independent test set for both recurrence free survival (RFS) and overall survival (OS). The Kaplan-Meier curves for OS stratified by cluster label show significant differences for both training and testing (p val < 0.0001). When compared to the models trained using clinical data only, the inclusion of the cluster label improves AUC test performance from .62 to .79 and from .66 to .80 for OS and RFS, respectively. The extraction of a single feature, namely a cluster label, to represent the high-dimensional radiomic feature space reduces the dimensionality and sparsity of the data. Moreover, inclusion of the cluster label improves model performance compared to clinical data only and offers comparable performance to the models including raw radiomic features.

Feature selection and machine learning methods for optimal identification and prediction of subtypes in Parkinson's disease

OBJECTIVES

The present work focuses on assessment of Parkinson's disease (PD), including both PD subtype identification (unsupervised task) and prediction (supervised task). We specifically investigate optimal feature selection and machine learning algorithms for these tasks.

METHODS

We selected 885 PD subjects as derived from longitudinal datasets (years 0-4; Parkinson's Progressive Marker Initiative), and investigated 981 features including motor, non-motor, and imaging features (SPECT-based radiomics features extracted using our standardized SERA software). Two different hybrid machine learning systems (HMLS) were constructed and applied to the data in order to select optimal combinations in both tasks: (i) identification of subtypes in PD (unsupervised-clustering), and (ii) prediction of these subtypes in year 4 (supervised-classification). From the original data based on years 0 (baseline) and 1, we created new datasets as inputs to the prediction task: (i,ii) CSD0 and CSD01: cross-sectional datasets from year 0 only and both years 0 & 1, respectively; (iii) TD01: timeless dataset from both years 0 & 1. In addition, PD subtype in year 4 was considered as outcome. Finally, high score features were derived via ensemble voting based on their prioritizations from feature selector algorithms (FSAs).

RESULTS

In clustering task, the most optimal combinations (out of 981) were selected by individual FSAs to enable high correlation compared to using all features (arriving at 547). In prediction task, we were able to select optimal combinations, resulting in an accuracy >90% only for timeless dataset (TD01); there, we were able to select the most optimal combination using 77 features, directly selected by FSAs. In both tasks, however, using combination of only high score features from ensemble voting did not enable acceptable performances, showing optimal feature selection via individual FSAs to be more effective.

CONCLUSION

Combining non-imaging information with SPECT-based radiomics features, and optimal utilization of HMLSs, can enable robust identification of subtypes as well as appropriate prediction of these subtypes in PD patients. Moreover, use of timeless dataset, beyond cross-sectional datasets, enabled predictive accuracies over 90%. Overall, we showed that radiomics features extracted from SPECT images are important in clustering as well as prediction of PD subtypes.

A three-stage, deep learning, ensemble approach for prognosis in patients with Parkinson's disease

BACKGROUND

Diagnosis of Parkinson's disease (PD) is informed by the presence of progressive motor and non-motor symptoms and by imaging dopamine transporter with [123I]ioflupane (DaTscan). Deep learning and ensemble methods have recently shown promise in medical image analysis. Therefore, this study aimed to develop a three-stage, deep learning, ensemble approach for prognosis in patients with PD.

METHODS

Retrospective data of 198 patients with PD were retrieved from the Parkinson's Progression Markers Initiative database and randomly partitioned into the training, validation, and test sets with 118, 40, and 40 patients, respectively. The first and second stages of the approach extracted features from DaTscan and clinical measures of motor symptoms, respectively. The third stage trained an ensemble of deep neural networks on different subsets of the extracted features to predict patient outcome 4 years after initial baseline screening. The approach was evaluated by assessing mean absolute percentage error (MAPE), mean absolute error (MAE), Pearson's correlation coefficient, and bias between the predicted and observed motor outcome scores. The approach was compared to individual networks given different data subsets as inputs.

RESULTS

The ensemble approach yielded a MAPE of 18.36%, MAE of 4.70, a Pearson's correlation coefficient of 0.84, and had no significant bias indicating accurate outcome prediction. The approach outperformed individual networks not given DaTscan imaging or clinical measures of motor symptoms as inputs, respectively.

CONCLUSION

The approach showed promise for longitudinal prognostication in PD and demonstrated the synergy of imaging and non-imaging information for the prediction task.

Radiomics approach in the neurodegenerative brain

As claimed by Robert Gilles et al., "Images are more than pictures, they are data". This statement refers to the power of imaging to provide large amounts of quantitative features for improving diagnosis, prognosis and therapy response. The conversion of digital medical images into high-dimensional mineable data is called radiomics. Radiomics analysis is based on data-characterisation algorithms which have the potential to uncover disease heterogeneity characteristics that might escape from the expert evaluation. This method has been widely applied in oncology and genetic fields, while the literature on neurodegenerative disorders is in its relative infancy. Here, we provide a preliminary evaluation of the main results reached applying radiomics analyses on well-established MRI features of patients with Alzheimer's Disease and Parkinson's disease.

Aberrant functional connectivity and activity in Parkinson's disease and comorbidity with depression based on radiomic analysis

INTRODUCTION

The current diagnosis of Parkinson's disease (PD) comorbidity with depression (DPD) largely depends on clinical evaluation. However, the modality may tend to lack precision in detecting PD with depression. A radiomic approach that combines functional connectivity and activity with clinical scores has the potential to achieve accurate and differential diagnosis between PD and DPD.

METHODS

In this study, we aimed to employ the radiomic approach to extract large-scale features of functional connectivity and activity for differentiating among DPD, PD with no depression (NDPD), and healthy controls (HC). We extracted 6,557 features of five types from all subjects including clinical characteristics, resting-state functional connectivity (RSFC), amplitude of low-frequency fluctuation (ALFF), regional homogeneity (ReHo), and voxel-mirrored homotopic connectivity (VMHC). Lasso, random forest, and support vector machine (SVM) were implemented for feature selection and dimension reduction based on the training sets, and the prediction performance for different methods in the testing sets was compared.

RESULTS

The results showed that nineteen features were selected for the group of DPD versus HC, 34 features were selected for the group of NDPD versus HC, and 17 features were retained for the group of DPD versus NDPD. In the testing sets, Lasso prediction achieved the accuracies of 0.95, 0.96, and 0.85 for distinguishing between DPD and HC, NDPD and HC, and DPD and NDPD, respectively. Random forest achieved the accuracies of 0.90, 0.82, and 0.90 for distinguishing between DPD and HC, NDPD and HC, and DPD and NDPD, respectively, while SVM yielded the accuracies of 1, 0.86 and 0.65 for distinguishing between DPD and HC, NDPD and HC, and DPD and NDPD, respectively.

CONCLUSIONS

By identifying aberrant functional connectivity and activity as potential biomarkers, the radiomic approach facilitates a deeper understanding and provides new insights into the pathophysiology of DPD to support the clinical diagnosis with high prediction accuracy.

Multivariate radiomics models based on 18F-FDG hybrid PET/MRI for distinguishing between Parkinson's disease and multiple system atrophy

PURPOSE

To construct multivariate radiomics models using hybrid ¹⁸F-FDG PET/MRI for distinguishing between Parkinson's disease (PD) and multiple system atrophy (MSA).

METHODS

Ninety patients (60 with PD and 30 with MSA) were randomized to training and test sets in a 7:3 ratio. All patients underwent ¹⁸F-fluorodeoxyglucose (¹⁸F-FDG) PET/MRI to simultaneously obtain metabolic images (¹⁸F-FDG), structural MRI images (T1-weighted imaging (T1WI), T2-weighted imaging (T2WI) and T2-weighted fluid-attenuated inversion recovery (T2/FLAIR)) and functional MRI images (susceptibility-weighted imaging (SWI) and apparent diffusion coefficient). Using PET and five MRI sequences, we extracted 1172 radiomics features from the putamina and caudate nuclei. The radiomics signatures were constructed with the least absolute shrinkage and selection operator algorithm in the training set, with progressive optimization through single-sequence and double-sequence radiomics models. Multivariable logistic regression analysis was used to develop a clinical-radiomics model, combining the optimal multi-sequence radiomics signature with clinical characteristics and SUV values. The diagnostic performance of the models was assessed by receiver operating characteristic and decision curve analysis (DCA).

RESULTS

The radiomics signatures showed favourable diagnostic efficacy. The optimal model comprised structural (T1WI), functional (SWI) and metabolic (¹⁸F-FDG) sequences (Radscore_{FDG_T1WI_SWI}) with the area under curves (AUCs) of the training and test sets of 0.971 and 0.957, respectively. The integrated model, incorporating Radscore_{FDG_T1WI_SWI}, three clinical symptoms (disease duration, dysarthria and autonomic failure) and SUV_max, demonstrated satisfactory calibration and discrimination in the training and test sets (0.993 and 0.994, respectively). DCA indicated the highest clinical benefit of the clinical-radiomics integrated model.

CONCLUSIONS

The radiomics signature with metabolic, structural and functional information provided by hybrid ¹⁸F-FDG PET/MRI may achieve promising diagnostic efficacy for distinguishing between PD and MSA. The clinical-radiomics integrated model performed best.

Improved motor outcome prediction in Parkinson's disease applying deep learning to DaTscan SPECT images

PURPOSE

Dopamine transporter (DAT) SPECT imaging is routinely used in the diagnosis of Parkinson's disease (PD). Our previous efforts demonstrated the use of DAT SPECT images in a data-driven manner by improving prediction of PD clinical assessment outcome using radiomic features. In this work, we develop a convolutional neural network (CNN) based technique to predict clinical motor function evaluation scores directly from longitudinal DAT SPECT images and non-imaging clinical measures.

PROCEDURES

Data of 252 subjects from the Parkinson's Progression Markers Initiative (PPMI) database were used in this work. The motor part (III) score of the unified Parkinson's disease rating scale (UPDRS) at year 4 was selected as outcome, and the DAT SPECT images and UPDRS_III scores acquired at year 0 and year 1 were used as input data. The specified inputs and outputs were used to develop a CNN based regression method for prediction. Ten-fold cross-validation was used to test the trained network and the absolute difference between predicted and actual scores was used as the performance metric. Prediction using inputs with and without DAT images was evaluated.

RESULTS

Using only UPDRS_III scores at year 0 and year 1, the prediction yielded an average difference of 7.6 ± 6.1 between the predicted and actual year 4 motor scores (range [5, 77]). The average difference was reduced to 6.0 ± 4.8 when longitudinal DAT SPECT images were included, which was determined to be statistically significant via a two-sample t-test, and demonstrates the benefit of including images.

CONCLUSIONS

This study shows that adding DAT SPECT images to UPDRS_III scores as inputs to deep-learning based prediction significantly improves the outcome. Without requiring segmentation and feature extraction, the CNN based prediction method allows easier and more universial application.

Artificial Intelligence and the Medical Physicist: Welcome to the Machine

Artificial intelligence (AI) is a branch of computer science dedicated to giving machines or computers the ability to perform human-like cognitive functions, such as learning, problem-solving, and decision making. Since it is showing superior performance than well-trained human beings in many areas, such as image classification, object detection, speech recognition, and decision-making, AI is expected to change profoundly every area of science, including healthcare and the clinical application of physics to healthcare, referred to as medical physics. As a result, the Italian Association of Medical Physics (AIFM) has created the “AI for Medical Physics” (AI4MP) group with the aims of coordinating the efforts, facilitating the communication, and sharing of the knowledge on AI of the medical physicists (MPs) in Italy. The purpose of this review is to summarize the main applications of AI in medical physics, describe the skills of the MPs in research and clinical applications of AI, and define the major challenges of AI in healthcare.

Radiomics and Prostate MRI: Current Role and Future Applications

Multiparametric prostate magnetic resonance imaging (mpMRI) is widely used as a triage test for men at a risk of prostate cancer. However, the traditional role of mpMRI was confined to prostate cancer staging. Radiomics is the quantitative extraction and analysis of minable data from medical images; it is emerging as a promising tool to detect and categorize prostate lesions. In this paper we review the role of radiomics applied to prostate mpMRI in detection and localization of prostate cancer, prediction of Gleason score and PI-RADS classification, prediction of extracapsular extension and of biochemical recurrence. We also provide a future perspective of artificial intelligence (machine learning and deep learning) applied to the field of prostate cancer.

Texture-based markers from structural imaging correlate with motor handicap in Parkinson's disease

There is a growing need for surrogate biomarkers for Parkinson’s disease (PD). Structural analysis using magnetic resonance imaging with T1-weighted sequences has the potential to quantify histopathological changes. Degeneration is typically measured by the volume and shape of morphological changes. However, these changes appear late in the disease, preventing their use as surrogate markers. We investigated texture changes in 108 individuals, divided into three groups, matched in terms of sex and age: (1) healthy controls (n = 32); (2) patients with early-stage PD (n = 39); and (3) patients with late-stage PD and severe L-dopa-related complications (n = 37). All patients were assessed in off-treatment conditions. Statistical analysis of first- and second-order texture features was conducted in the substantia nigra, striatum, thalamus and sub-thalamic nucleus. Regions of interest volumetry and voxel-based morphometry were performed for comparison. Significantly different texture features were observed between the three populations, with some showing a gradual linear progression between the groups. The volumetric changes in the two PD patient groups were not significantly different. Texture features were significantly associated with clinical scores for motor handicap. These results suggest that texture features, measured in the nigrostriatal pathway at PD diagnosis, may be useful in predicting clinical progression of motor handicap.

Robust identification of Parkinson's disease subtypes using radiomics and hybrid machine learning

OBJECTIVES

It is important to subdivide Parkinson's disease (PD) into subtypes, enabling potentially earlier disease recognition and tailored treatment strategies. We aimed to identify reproducible PD subtypes robust to variations in the number of patients and features.

METHODS

We applied multiple feature-reduction and cluster-analysis methods to cross-sectional and timeless data, extracted from longitudinal datasets (years 0, 1, 2 & 4; Parkinson's Progressive Marker Initiative; 885 PD/163 healthy-control visits; 35 datasets with combinations of non-imaging, conventional-imaging, and radiomics features from DAT-SPECT images). Hybrid machine-learning systems were constructed invoking 16 feature-reduction algorithms, 8 clustering algorithms, and 16 classifiers (C-index clustering evaluation used on each trajectory). We subsequently performed: i) identification of optimal subtypes, ii) multiple independent tests to assess reproducibility, iii) further confirmation by a statistical approach, iv) test of reproducibility to the size of the samples.

RESULTS

When using no radiomics features, the clusters were not robust to variations in features, whereas, utilizing radiomics information enabled consistent generation of clusters through ensemble analysis of trajectories. We arrived at 3 distinct subtypes, confirmed using the training and testing process of k-means, as well as Hotelling's T2 test. The 3 identified PD subtypes were 1) mild; 2) intermediate; and 3) severe, especially in terms of dopaminergic deficit (imaging), with some escalating motor and non-motor manifestations.

CONCLUSION

Appropriate hybrid systems and independent statistical tests enable robust identification of 3 distinct PD subtypes. This was assisted by utilizing radiomics features from SPECT images (segmented using MRI). The PD subtypes provided were robust to the number of the subjects, and features.

Radiomics-based machine learning model to predict risk of death within 5-years in clear cell renal cell carcinoma patients

PURPOSE

The aim of this study was to develop radiomics-based machine learning models based on extracted radiomic features and clinical information to predict the risk of death within 5 years for prognosis of clear cell renal cell carcinoma (ccRCC) patients.

METHODS

According to image quality and clinical data availability, we eventually selected 70 ccRCC patients that underwent CT scans. Manual volume-of-interest (VOI) segmentation of each image was performed by an experienced radiologist using the 3D slicer software package. Prior to feature extraction, image pre-processing was performed on CT images to extract different image features, including wavelet, Laplacian of Gaussian, and resampling of the intensity values to 32, 64 and 128 bin levels. Overall, 2544 3D radiomics features were extracted from each VOI for each patient. Minimum Redundancy Maximum Relevance (MRMR) algorithm was used as feature selector. Four classification algorithms were used, including Generalized Linear Model (GLM), Support Vector Machine (SVM), K-nearest Neighbor (KNN) and XGBoost. We used the Bootstrap resampling method to create validation sets. Area under the receiver operating characteristic (ROC) curve (AUROC), accuracy, sensitivity, and specificity were used to assess the performance of the classification models.

RESULTS

The best single performance among 8 different models was achieved by the XGBoost model using a combination of radiomic features and clinical information (AUROC, accuracy, sensitivity, and specificity with 95% confidence interval were 0.95-0.98, 0.93-0.98, 0.93-0.96 and ~1.0, respectively).

CONCLUSIONS

We developed a robust radiomics-based classifier that is capable of accurately predicting overall survival of RCC patients for prognosis of ccRCC patients. This signature may help identifying high-risk patients who require additional treatment and follow up regimens.