Participants
All data used in this study were obtained from the Parkinson’s Progression Markers Initiative (PPMI) database (www.ppmi-info.org/data) accessed on 3 April 2021. The dataset contained 790 pre-processed 123I-FP-CIT DAT-SPECT images acquired at the screening stage. In this study, we included cases acquired with the same manufacturer’s SPECT system (SIEMENSE) and were continuously enrolled in PPMI from the 790 cases. Table 1 shows the characteristics of the patients in this study. There were 291 patients in both the normal control (NC) and PD groups. The NC and PD group comprised 70 (49 men and 21 women, 61.2 ± 9.6 years) and 221 (140 men and 81 women, 60.0 ± 11.0 years) individuals, respectively. All data were divided into the training and test sets at 7:3 so that the ratio of the NC and PD groups would be constant. The number of the training and test sets for the NC and PD groups were 49 and 155 and 21 and 66, respectively.
Table 1 Patient characteristics
Characteristics
|
NC
|
PD
|
Number of patients
|
70
|
221
|
Male/female
|
49/21
|
140/81
|
Age
|
61.2 ± 9.6
|
60.0 ± 11.0
|
Age of diagnosis
|
NA
|
60.7 ± 9.5
|
Hoehn and Yahr stage
|
0.01 ± 0.1
|
1.58 ± 0.5
|
MDS-UPDRS Ⅲ
|
1.0 ± 1.9
|
21.8 ± 9.3
|
MoCA
|
28.2 ± 1.1
|
27.3 ± 2.3
|
MDS-UPDRS Movement Disorder Society-Unified Parkinson’s Disease Rating Scale, MoCA Montreal Cognitive Assessment, NC normal control, PD Parkinson’s disease, NA not applicable
Reconstruction and spatial normalisation of SPECT images
Reconstructed DAT-SPECT images were downloaded from the PPMI website. As per PPMI documentation, pre-processing steps were performed at the Institute for Neurodegenerative Disorders and included the following steps: SPECT imaging and reconstruction: SPECT imaging was acquired at each imaging centre as per the PPMI imaging protocol and sent to the institute for neurodegenerative disorders for processing. SPECT raw projection data were imported to a HERMES (Hermes Medical Solutions, Stockholm, Sweden) system for iterative reconstruction. Iterative reconstruction was performed without filtering. The reconstructed files were transferred to the PMOD (PMOD Technologies, Zurich, Switzerland) for subsequent processing. Attenuation correction ellipses were drawn on the images, and a Chang 0 attenuation correction was applied to images utilising a site-specific μ that was empirically derived from phantom data acquired during site initiation for the trial. Once attenuation correction was completed, a standard Gaussian three-dimensional (3D) 6.0-mm filter was applied.
Then, the DAT-SPECT images were spatially normalised to Montreal Neurologic Institute (MNI) space using statistical parametric mapping (SPM12, Wellcome Trust Centre for Neuroimaging, London, UK) in MATLAB R2021a (version 9.10, The MathWorks, Inc. Massachusetts, USA). DAT-SPECT images were spatially normalised to the MNI-based template of 123I-FP-CIT [27, 28] using the old normalise function under identical conditions. After spatial normalisation, the radiological technologist with 13-year clinical experience visually assessed for misalignment between DAT-SPECT and the template. The pre-processed images were saved in the Neuroimaging Informatics Technology Initiative format using 91 × 109 × 91 isotropic voxels of 2 mm.
Calculation of radiomics features and semi-quantitative indicators
The automated anatomical labelling atlas (AAL) 3 [29] voxel of interest (VOI) template was used to calculate the radiomics features. The feature calculation VOIs were the caudate nucleus, putamen, and pallidum (Fig. 1). The pallidum was added to the VOIs to compensate for errors in spatial normalisation. Radiomics features were calculated using Standardized Environment for Radiomics Analysis (SERA) [30–32] worked on MATLAB. One hundred eighty-six image biomarker standardisation initiative-standardised features [24] were calculated using SERA, including 50 first-order features (statistical, histogram, and intensity histogram features) and higher-order136 3D features (Table 2). A total of 558 radiomics features were calculated for the caudate, putamen, and pallidum VOIs. We also calculated the ratio of the caudate to the putamen or pallidum of radiomics features. These totalled 930 radiomics features. Furthermore, the SUR of the caudate nucleus (SURcaudate), putamen (SURputamen), and pallidum (SURpallidum) were calculated as conventional semi-quantification indices. The SUR was calculated using the following formula [33]:

where Cstriatum is the average count of the caudate nucleus, putamen, or pallidum and Cbackground is the average count of the occipital lobe. In addition, the ratios of the caudate to the putamen or pallidum (CRputamen, CRpallidum) were calculated. All the semi-quantitative indices were compared between the NC and PD groups, and receiver operating characteristic (ROC) [34] analysis was performed.
Table 2 Number of radiomics features per region and their family names
Feature family
|
Number of features
|
Local intensity
|
2
|
Intensity-based statistics
|
18
|
Intensity histogram
|
23
|
Intensity-volume histogram
|
7
|
GLCM
|
50
|
GLRLM
|
32
|
GLZSM
|
16
|
GLDZM
|
16
|
NGTDM
|
5
|
NGLDM
|
17
|
Total
|
186
|
GLCM grey-level co-occurrence matrix, GLRLM grey-level run-length matrix, GLZSM grey-level zone size matrix, GLDZM grey-level distance zone matrix, NGTDM neighbourhood grey tone difference matrix, NGLDM neighbourhood grey-level dependence matrix
Radiomics feature selection and signature construction
The least absolute shrinkage and selection operator (LASSO) [35] function in MATLAB was used to select effective features from the radiomics features. All radiomics features were z-scored to mean 0 and standard deviation 1.0 before being inputted to LASSO. LASSO permits the estimation and selection of explanatory variables [36, 37], that is, radiomics features with non-zero coefficients. For the selection of radiomics features using LASSO, a 10-fold cross-validation test was conducted using the training set. Furthermore, the linear combination sum of features with non-zero coefficients when the largest lambda value such that the mean square error (MSE) is within one standard error of the minimum MSE was used as the radiomics signature. We compared the classification performance of the radiomics signature and semi-quantitative indicator that showed the highest classification performance.
Classification model construction with radiomics signature and semi-quantitative indicator
The classification models for the NC and PD groups were constructed using the radiomics signature and/or semi-quantitative indicator. The four classifiers used were support vector machine (SVM), k-nearest neighbour (KNN), linear discriminant analysis (LDA), and decision tree. The main parameters of each classifier were as follows: SVM (BoxConstraint = 1, KernelScale = 1, KernelFunction = polynomial [order = 3]), KNN (NumNeighbors = 1, Distance = minkowski, Exponent = 2), LDA (Gamma = 0), and decision tree (MinLeafSize = 1, MinParentSize = 10). The features used were radiomics signature alone, semi-quantitative indicator alone, and both. The training set was used to train the classifier, and the performance of each classification model was evaluated using the test set. Classification performance was evaluated using the area under the ROC curve (AUC).
Statistical analyses
The radiomics signature and SURs in the NC and PD groups were tested for significant differences using the Wilcoxon rank-sum test. ROC analysis was performed using semi-quantitative indicators and radiomics signature. We used the DeLong [38] test to examine the differences in AUCs, and for multiple comparisons, the Bonferroni correction was performed. The sensitivity, specificity, and accuracy of semi-quantitative indices and radiomics signature were calculated using the optimal cut-off values determined based on ROC analysis. Differences were considered statistically significant at P < 0.05. All statistical analyses were performed using RStudio (version 1.4.1106).