- 1Bioengineering College of Chongqing University, Chongqing University Central Hospital (Chongqing Emergency Medical Center), Chongqing, China
- 2Department of Radiology, Chongqing Western Hospital, Chongqing, China
- 3Department of Medical Service, Yanzhuang Central Hospital of Gangcheng District, Chongqing, China
- 4Department of Research and Development, Shanghai United Imaging Intelligence, Co., Ltd. Shanghai, China
Background: Parkinson’s disease (PD) is the second most common neurodegenerative disease. An objective diagnosis method is urgently needed in clinical practice. In this study, deep learning and radiomics techniques were studied to automatically diagnose PD from healthy controls (HCs).
Methods: 155 PD patients and 154 HCs were randomly divided into a training set (246 patients) and a testing set (63 patients). The brain subregions identification and segmentation were automatically performed with a VB-net, and radiomics features of billateral thalamus, caudatum, putamen and pallidum were extracted. Five independent machine learning classifiers [Support Vector Machine (SVM), Stochastic gradient descent (SGD), random forest (RF), quadratic discriminant analysis (QDA) and decision tree (DT)] were trained on the training set, and validated on the testing. Delong test was used to compare the performance of different models.
Results: Our VB-net could automatically identify and segment the brain into 109 regions. 2,264 radiomics features were automatically extracted from the billateral thalamus, caudatum, putamen or pallidum of each patient. After four step of features dimensionality reduction, Delong tests showed that the SVM model based on combined features had the best performance, with AUCs of 0.988 (95% CI: 0.979 ~ 0.998, specificity = 91.1%, sensitivity =100%, accuracy = 89.4% and precision = 88.2%) and 0.976 (95% CI: 0.942 ~ 1.000, specificity = 100%, sensitivity = 87.1%, accuracy = 93.5% and precision = 88.6%) in the training set and testing set, respectively. Decision curve analysis showed that the clinical benefit of the line graph model was high.
Conclusion: The SVM model based on combined features could be used to diagnose PD with high accuracy. Our fully automatic model could rapidly process the MRI data and distinguish PD and HCs in one minute. It greatly improved the diagnostic efficiency and has a great potential value in clinical practice to help the early diagnosis of PD.
Introduction
Parkinson’s disease (PD) is a complex and progressive neurodegenerative disorder characterized by an insidious onset, high incidence, and significant disability rate (1, 2). It poses a serious threat to the physical and mental health, as well as the overall quality of life, of middle-aged and elderly individuals (3). Its primary clinical manifestations include motor symptoms such as tremors, muscle rigidity, bradykinesia, postural instability; non-motor symptoms such as sleep disorders and olfactory dysfunction; autonomic nervous system dysfunction; cognitive impairment; and psychiatric disturbances (4). Currently, more than 6 million people worldwide suffer from PD, and the number is expected to further increase, bringing a huge burden to families and the society (5).
The rapid and accurate diagnosis of PD is of great significance for targeted treatment, prevention of disease progression, improvement of quality of life and overall prognosis (6). At present, the diagnosis of PD still relies on subjective clinical symptoms (7). Objective diagnosis methods are urgently needed in clinical practice. As a medical imaging technique, MRI has the advantages of non-invasive, non-radiation exposure, and high-resolution capabilities, making it widely used in the diagnosis and staging of neurological diseases (8, 9). Previous studies have identified alterations in both the structure and function of the brain in individuals diagnosed with PD (10, 11). Vogt et al. (12) discovered that the cingulate cortex played a crucial role in identifying new biomarkers for patients with early PD. Nyberg et al. (13) found a significant increasing in the volume of bilateral hippocampal and right nucleus accumbens of PD patients. Kassubek et al. discovered a significant increasing in the gray matter volume of the ventral medial thalamic nucleus on the contralateral side of the tremor limb in patients with PD. Furthermore, they observed a positive correlation between changes in thalamic gray matter volume and tremor amplitude (14). However, these studies mainly focused on macroscopic changes of the brain of PD patients, but overlook the small structural changes.
Radiomics, proposed by Lambin et al. (15), is a new field of computer-aided imaging diagnosis that assists in diagnosing and differentiating diseases by quantifying subtle information in medical images that is difficult to evaluate with the naked eye. Tupe-Waghmare et al. (16) had extracted radiomic features from the subcortical structure, cerebellum, brainstem, and used a random forest machine learning model to identify PD and HCs with an accuracy of 70%. Tafuri et al. (17) used Freesurfer software to extract radiomics features from the subcortical nucleus and used SVM model to distinguish PD and HCs patients, achieving an AUC (area under the receiver operating characteristic curve) of 0.77. However, these methods mostly showed low accuracy and were time-consuming, taking approximately 4 h to process a patient. These shortcomings limited its clinical application. In our previous study, we had developed a CNN-based artificial intelligence model for the automatic segmentation and measurement of the whole brain (109 brain regions), and the entire segmentation and reconstruction process took less than half a minute (18). In this study, we used this network to quickly segment and analyze the radiomics characteristics of the cerebral cortex and neuclei, and try to establish artificial intelligence models to help distinguish PD and HCs.
Materials and methods
Participants
All PD patients included in this study were sourced from the PD Progression Initiative (PPMI) database. The PD subjects in the PPMI were newly diagnosed patients who were not receiving any medication. The clinical diagnostic criteria for PD were based on the Movement Disorder Society guidelines. For the most up-to-date information, please visit https://www.ppmi-info.org/. The T1-weighted MR images were obtained using a 3 T Tesla scanner manufactured by Siemens. The imaging parameters were as follows: Acquisition Type = 3D; Flip Angle = 9.0 degree; field of view (FOV) = 256 × 256; matrix = 256 × 256; TR = 2300.0 ms; TE = 3.0 ms; Slice Thickness = 1.0 mm; interslice gap = 0 mm. A total of 155 PD patients and 154 HCs were enrolled in the study. The HC subjects and PD patients were matched for age and sex. All patients were randomly assigned to a training set and a testing set in a ratio of 8:2.
Brain subregions segmentation
The brain subregions segmentation module was implemented using a deep learning algorithm based on a 3D VB-NET network. The data preprocessing module performed a series of operations, including rotation, resampling, resizing, skull stripping, image non-uniform correction, histogram matching, and gray-scale normalization, on the MRI images used for training and testing. All images are standardized to the size of 256*256*256*1 mm3 in the standard Cartesian LPI coordinate system, and the gray-scale range was within the interval (−1, 1). The network training module used an end-to-end deep convolutional neural network, taking each sample image and its corresponding brain substructure partition atlas as the training sample. The sample image was the network input, and the output label was the brain atlas correspondent to the sample image. The network parameters were adjusted according to the difference between the output brain division and the actual brain division, and the training continues until the network basically converges and the output label image was substantially consistent with the partition image corresponding to the sample. In the overall network training process, a coarse-to-fine cascading segmentation strategy was used, simplifying the complexity and difficulty of the brain segmentation problem by step decomposition, providing extra information to the lower level network by the upper level network to enhance the network’s segmentation performance, and achieving fine segmentation of the large brain region, medium brain region, and brain substructures on a stage-by-stage basis. The model was constructed based on 1,800 subjects and evaluation showed an averaged 0.92 Dice overlap with ground truth on 295 subjects. Detailed segmentation process information was descripted in our previously published literature (18).
The entire brain was automatically divided into 109 subregions, which included 22 subregions in the temporal lobe, 20 subregions in the frontal lobe, 12 subregions in the parietal lobe, 8 subregions in the occipital lobe, 8 subregions in the cingulate gyrus, 2 subregions in the insula, 12 subcortical neuclei, white matter structures, ventricles, cerebellum structures and other structures (Supplementary material 1). The segmentation process took less than half a minute for each patient. The flow chart of brain subregions segmentation was shown in Figure 1.
Features extraction and dimensionality reduction
The volume of 109 brain subregions of each patient was automatically extracted through the deep learning model. 2,264 radiomics features were automatically extracted from the billateral thalamus, caudatum, putamen or pallidum of each patient. There were 18 first-order statistics and 14 shape features, which reflect the shape and size of the region accurately. Texture features included 21 Gray Level Co-occurrence Matrix (GLCM) features, 16 Gray Level Run Length Matrix (GLRLM) features, 16 Gray Level Size Zone Matrix (GLSZM) features, 5 Neighbouring Gray Tone Difference Matrix (NGTDM) features, and 14 Gray Level Dependence Matrix (GLDM) features. The high level features were obtained through 24 filters (including Box Mean, additive Gaussian Noise, binomial blur, curvature flow, Box-Sigma, normalization, Laplace Sharpening, discrete Gaussian, mean, speck noise, recursive Gaussian, Shot Noise and LoG with sigma values of 0.5, 1, 1.5, and 2), as well as wavelet transformations (LLL, LLH, LHL, LHH, HLL, HLH, HHL, and HHH). All radiomic features were then normalized using z-score normalization, and the reproducibility of these features was assessed using a pipeline that adheres to the recommendations of the Image Biomarker Standardization Initiative.
The relief and least absolute shrinkage and selection operator (LASSO) method were used to select the most robust features. Hyperparameter for LASSO was evaluated using stratified 5-fold cross-validation-based grid search method on the training set. The parameters that provided the highest cross-validation AUC was selected.
Models building and evaluation
Based on the retained features, five independent machine learning classifiers, including Support Vector Machine (SVM), Stochastic gradient descent (SGD), random forest (RF), quadratic discriminant analysis (QDA) and decision tree (DT) algorithm were trained on the training set, and validated on the testing set. The performance of classifier models on the test subset was evaluated by the mean and 95% confidence intervals (CI) of the accuracy, sensitivity/recall, specificity, and precision based on a case probability cut-off value of 0.5, as well as the F-score metric and area under the receiver operating characteristic curve (AUC). The calibration curve was used to evaluate the calibration of the model, and DCA was used to evaluate the clinical applicability of the predictive model. The flow chart of our study was shown in Figure 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1303501/fmed-10-1303501-HTML-r1/image_m/fmed-10-1303501-g002.jpg)
Figure 2. The flow chart of brain subregions segmentation, features extraction, features dimension reduction and models construction.
Statistical analysis
Statistical analyses were performed with R (version 4.0.4) and Python software. Mann Whitney U test or Student’s t test was used for the continuous variables according to the test of normal distribution. The chi-square test was used to compare categorical variables. Statistical significance was indicated by a two-tailed p value <0.05. Delong test was used for the comparison of different models.
Results
A total of 155 PD patients were enrolled in the study, including 99 males and 56 females. Meanwhile, 154 HCs were included, including 99 males and 55 females. There was no significant difference in age and gender between the two groups. The demographic and clinical features of PD and HC were shown in Table 1.
In the training set, the volume of 109 brain subregions of each patient was automatically obtained. 100 radiomics features were selected from 4,528 radiomics features of the bilateral thalamus using the relief method, and then 16 optimal features were obtained using the LASSO method. According to the same method, 12, 15, and 10 optimal features were obtained from the caudatum, putamen, and pallidum, respectively. Finally, feature selection of relief method and LASSO were performed again for all the above brain regions together, to obtain the 13 optimal combined features (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1303501/fmed-10-1303501-HTML-r1/image_m/fmed-10-1303501-g003.jpg)
Figure 3. Lasso algorithm for features selection of PD patients and HCs (A). The Lasso path displayed coefficient profiles for radiomic features across the entire range of possible values (B). Rad_Score distribution in the training sets (C) and testing set (D). “0” group represented HCs. “1” group represented PD.
Twenty-five models were established based on algorithms of SVM, SGD, RF, QDA, DT and thalamus features, caudatum features, putamen features, pallidum feature and combined features. ROC curves of SVM, SGD, RF, QDA and DT models based on thalamus, caudatum, putamen, pallidum and combined features were shown in Table 2. It was found that the SVM model based on the combined features (5 features from pallidum, 4 features from putamen, 4 features from caudatum) showed the best performance, with AUCs of 0.988 (95% CI: 0.979 ~ 0.998, specificity = 91.1%, sensitivity =100%, accuracy = 89.4% and precision = 88.2%), 0.976 (95% CI: 0.942 ~ 1.000, specificity = 100%, sensitivity = 87.1%, accuracy = 93.5% and precision = 88.6%) in the training set and testing set, respectively (Figure 4). The calibration curve showed a good agreement between the actual and predicted probabilities of the samples (Figure 5A). Decision curve analysis showed that the clinical benefit of the line graph model was high (Figure 5B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1303501/fmed-10-1303501-HTML-r1/image_m/fmed-10-1303501-g004.jpg)
Figure 4. Receiver Operating Characteristic (ROC) curves for SVM, SGD, RF, QDA and DT models in training set (A) and testing set (B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1303501/fmed-10-1303501-HTML-r1/image_m/fmed-10-1303501-g005.jpg)
Figure 5. Calibration curves for nomogram goodness of fit (A). Decision curves of different models (B).
Discussion
It is well known that PD is a complex progressive neurodegenerative disease with high incidence rate (19). Developing an objective, accurate, and effective method to distinguish PD and HC has become an urgent issue in clinical practice. The changes in the cerebral cortex, deep nucleus, cerebellum and ventricles were found closely related to the occurrence of PD (20). Traditional imaging diagnostic methods were time-consuming, subjective and unable to detect these subtle changes in brain structure (21). Deep learning and radiomics were new fields of computer-aided imaging diagnosis (22), which overcomed the limitations of visual diagnosis and could quantify the subtle information in medical images, providing new hope for the rapid and accurate clinical diagnosis of PD.
In this study, a VB-net network evolved from U-Net network was developed to automatically identify and segment the cerebral cortex, cerebellum, ventricle, and 14 subcortical nuclei in 3D-T1-weighted imaging sequences. 2,264 radiomics features were automatically extracted from each region of the bilateral thalamus, caudatum, putamen, or pallidum. Five machine learning algorithms were applied to classify PD and HCs automatically and an independent verification group was set up to verify the performance of the model. We found among SVM, SGD, RF, QDA and DT classifiers, SVM classifier had the highest classification performance, with an AUC of 0.988 in the training set and 0.974 in the testing set. Our results demonstrated that the application of artificial intelligence technology to analyze raw T1-weighted MRI images could accurately differentiation PD patients from HCs. As a completely objective method, it did not rely on patients’ personal history or doctors’ clinical experience, making it more subjective. To the best of our knowledge, this study achieved the first fully automatic differentation of PD and HCs with high accuracy. Previously, Liu et al. (23) had used a radiomic model based on T2-weighted images of caudate nucleus and putamen to distinguish between PD and HCs, and the AUCs in the training set and testing set were 0.8767 (95% CI: 0.8066 ~ 0.9469) and 0.7143 (95% CI:0.5540 ~ 0.8746), respectively. Their results showed low accuracy and the manual brain segmentation took a long time, greatly reducing its work efficiency and clinical value. Adeli et al. (24) had utilized MRI and Single-Photon Emission Computed Tomography (SPECT) data from 538 subjects in the PPMI database to establish a classification framework, achieving a diagnostic accuracy of 97.5% for distinguishing between PD and HCs. However, the accuracy of their study mainly depended on the SPECT data, while the MRI data only made a small contribution to diagnostic ability. Additionally, the SPECT examination was expensive and involved radiopharmaceuticals that were not easily acceptable by patients (25). In our study, conventional MRI T1 weighted images were used, and all processes were fully automated, resulting in significant clinical significance. Our study was the first to use a deep learning model to automatically detect and segment brain subregions and extract valuable information. Our model could identify PD patients within 1 min and provided high accuracy. With this artificial intelligence system, doctors could quickly and accurately diagnose PD and perform early treatment timely, which could significantly improve prognosis. On the other hand, it could improve medical efficiency, shorten patient waiting time and reduce medical labor costs.
In the optimal model of this study, 13 radiomics features were retained, including 5 features from pallidum, 4 features from putamen and 4 features from caudatum. The 5 features from the pallidum included one Gray-level dependence matrix (GLDM) feature, one Gray-level run length matrix (GLRLM) feature, one first-order feature and two GLDM features. The pallidum receives afferent fibers from the subthalamic nucleus and is the primary region of the basal ganglia that emits efferent fibers (26). It could continuously release inhibitory neurotransmitter γ-aminobutyric acid (GABA), which may cause tremors (27). Hutchison et al. (28) had discovered the presence of neurons that exhibit tremor frequency activity ranging from 4 to 6 Hz in the pallidum, thus supporting the significant role of the pallidum in generating resting tremors in patients with PD. The first-order feature of it reflected changes in the shape and volume of the pallidum (29). Previous studies had demonstrated significant alterations in pallidum volume between PD patients and HCs (30, 31). The loss and degeneration of pallidus neurons, as well as the proliferation of glial cells, could lead to texture changes in MRI images, which were reflected in GLDM and GLRLM. The four features of the putamen were all first-order characteristics, including Mean Absolute Deviation, Kurtosis, Energy, and Minimum. Previous studies had indicated that the putamen appeared to be the first region affected in PD (32). The four features of the caudatum included one first-order feature, one GLCM feature and two GLRLM features. The caudatum is a key substructure of the unique basal ganglia circuit associated with emotional and psychomotor fatigue. It plays a crucial role in the pathological and physiological regulation of PD (33). Previous studies had proved that dopaminergic dysfunction could lead to damage in the caudatum of PD patients, resulting in morphological and pathological changes (34, 35). The first-order features reflected the asymmetry and flatness of the morphology of the caudate nucleus, while the features of GLCM and GLRLM reflected the roughness and heterogeneity caused by pathological changes in the internal structure of the caudate nucleus.
Our study had several limitations. Firstly, this was a retrospective cross-sectional study that did not reflect the dynamic changes of the brain in PD. A prospective longitudinal follow-up study was needed in the future. Secondly, the sample size of PD patients was relatively small and all data were obtained through 3 T Siemens scanners; therefore, the generalization of the model needed to be further verified. In the future, images from multiple research centers were required for independent external validation, so that it could be adapted to a wider range of clinical scenarios. Finally, although we had confirmed that radiomics features have high accuracy in distinguishing PD from HCs, the pathological basis behind them still needed to be further revealed.
Conclusion
In this study, we established an AI model which could distinguish PD and HCs accurately. It was fully automated and could quickly process the routine MRI data within one minute to obtain accurate results. Our research results have greatly improved the diagnostic efficiency and had a great potential value in clinical practice to help the early diagnosis of PD.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ppmi-info.org/.
Ethics statement
The studies involving humans were approved by The PPMI study was approved by an ethics standards committee on human experimentation at each institution. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
ChaL: Conceptualization, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing, Resources. DH: Data curation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing. FW: Methodology, Writing – original draft, Resources. YX: Formal Analysis, Methodology, Software, Supervision, Validation, Writing-original draft. FS: Formal Analysis, Methodology, Software, Supervision, Validation, Writing-original draft. MY: Writing – original draft, Methodology, Software, Validation. JZ: Data curation, Writing – review & editing. CP: Model validation and analysis of Rad_ Score values. JF: Conceptualization, Funding acquisition, Writing – original draft, Writing – review & editing, Data curation, Investigation, Methodology, Resources, Software, Validation. ChuL: Conceptualization, Funding acquisition, Writing – original draft, Writing – review & editing, Investigation, Methodology, Resources.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Fundamental Research Funds for the Central Universities (project no. 2022CDJYGRH-004); Natural Science Foundation Project of Chongqing (cstc2020jcyj-msxm0030, cstb2023nscq-bhx00074); and Chongqing Medical Scientific Research Project (Joint project of Chongqing Health Commission and Science and Technology Bureau) (2022QNXM013).
Conflict of interest
Authors YX and FS were employed by company Shanghai United Imaging Intelligence, Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2023.1303501/full#supplementary-material
References
1. Tolosa, E, Garrido, A, Scholz, SW, and Poewe, W. Challenges in the diagnosis of Parkinson's disease. Lancet Neurol. (2021) 20:385–97. doi: 10.1016/S1474-4422(21)00030-2
2. Albin, RL. Parkinson's disease: background, diagnosis, and initial management. Clin Geriatr Med. (2006) 22:735–51. v. doi: 10.1016/j.cger.2006.06.003
3. Qi, S, Yin, P, Wang, L, Qu, M, Kan, GL, Zhang, H, et al. Prevalence of Parkinson's disease: a community-based study in China. Mov Disord. (2021) 36:2940–4. doi: 10.1002/mds.28762
4. Bloem, BR, Okun, MS, and Klein, C. Parkinson's disease. Lancet (London, England). (2021) 397:2284–303. doi: 10.1016/S0140-6736(21)00218-X
5. Simon, DK, Tanner, CM, and Brundin, P. Parkinson disease epidemiology, pathology, genetics, and pathophysiology. Clin Geriatr Med. (2020) 36:1–12. doi: 10.1016/j.cger.2019.08.002
6. Mele, B, Van, S, Holroyd-Leduc, J, Ismail, Z, Pringsheim, T, and Goodarzi, Z. Diagnosis, treatment and management of apathy in Parkinson's disease: a scoping review. BMJ Open. (2020) 10:e037632. doi: 10.1136/bmjopen-2020-037632
7. Waller, S, Williams, L, Morales-Briceño, H, and Fung, VS. The initial diagnosis and management of Parkinson's disease. Aust J Gen Pract. (2021) 50:793–800. doi: 10.31128/AJGP-07-21-6087
8. Donzuso, G, Agosta, F, Canu, E, and Filippi, M. MRI of motor and nonmotor therapy-induced complications in Parkinson's disease. Mov Dis. (2020) 35:724–40. doi: 10.1002/mds.28025
9. Zhou, Y, Song, Z, Han, X, Li, H, and Tang, X. Prediction of Alzheimer's disease progression based on magnetic resonance imaging. ACS Chem Neurosci. (2021) 12:4209–23. doi: 10.1021/acschemneuro.1c00472
10. Baggio, HC, and Junqué, C. Functional MRI in Parkinson's disease cognitive impairment. Int Rev Neurobiol. (2019) 144:29–58. doi: 10.1016/bs.irn.2018.09.010
11. Xu, X, Han, Q, Lin, J, Wang, L, Wu, F, and Shang, H. Grey matter abnormalities in Parkinson's disease: a voxel-wise meta-analysis. Eur J Neurol. (2020) 27:653–9. doi: 10.1111/ene.14132
12. Vogt, BA. Cingulate cortex in Parkinson's disease. Handb Clin Neurol. (2019) 166:253–66. doi: 10.1016/B978-0-444-64196-0.00013-3
13. Nyberg, EM, Tanabe, J, Honce, JM, Krmpotich, T, Shelton, E, Hedeman, J, et al. Morphologic changes in the mesolimbic pathway in Parkinson's disease motor subtypes. Parkinsonism Relat Disord. (2015) 21:536–40. doi: 10.1016/j.parkreldis.2015.03.008
14. Kassubek, J, Juengling, FD, Hellwig, B, Spreer, J, and Lücking, CH. Thalamic gray matter changes in unilateral parkinsonian resting tremor: a voxel-based morphometric analysis of 3-dimensional magnetic resonance imaging. Neurosci Lett. (2002) 323:29–32. doi: 10.1016/s0304-3940(02)00111-8
15. Lambin, P, Rios-Velazquez, E, Leijenaar, R, Carvalho, S, van Stiphout, RG, Granton, P, et al. Radiomics: Extracting more information from medical images using advanced feature analysis. Eur J Cancer. (2012) 48:441–6. doi: 10.1016/j.ejca.2011.11.036
16. Tupe-Waghmare, P, Rajan, A, Prasad, S, Saini, J, Pal, PK, and Ingalhalikar, M. Radiomics on routine T1-weighted MRI can delineate Parkinson's disease from multiple system atrophy and progressive supranuclear palsy. Eur Radiol. (2021) 31:8218–27. doi: 10.1007/s00330-021-07979-7
17. Tafuri, B, Lombardi, A, Nigro, S, Urso, D, Monaco, A, Pantaleo, E, et al. The impact of harmonization on radiomic features in Parkinson's disease and healthy controls: a multicenter study. Front Neurosci. (2022) 16:1012287. doi: 10.3389/fnins.2022.1012287
18. Liu, S, Jie, C, Zheng, W, Cui, J, and Wang, Z. Investigation of underlying association between whole brain regions and Alzheimer's disease: a research based on an artificial intelligence model. Front Aging Neurosci. (2022) 14:872530. doi: 10.3389/fnagi.2022.872530
19. Tysnes, OB, and Storstein, A. Epidemiology of Parkinson's disease. J Neural Transm. (1996) 124:901–5. doi: 10.1007/s00702-017-1686-y
20. Li, R, Zou, T, Wang, X, Wang, H, Hu, X, Xie, F, et al. Basal ganglia atrophy-associated causal structural network degeneration in Parkinson's disease. Hum Brain Mapp. (2022) 43:1145–56. doi: 10.1002/hbm.25715
21. Cabreira, V, and Massano, J. Parkinson's disease: clinical review and update. Acta medica portuguesa. (2019) 32:661–70. doi: 10.20344/amp.11978
22. Loh, HW, Hong, W, Ooi, CP, Chakraborty, S, Barua, PD, Deo, RC, et al. Application of deep learning models for automated identification of Parkinson's disease: a review (2011-2021). Sensors (Basel). (2021) 21:7034. doi: 10.3390/s21217034
23. Liu, P, Wang, H, Zheng, S, Zhang, F, and Zhang, X. Parkinson's disease diagnosis using Neostriatum Radiomic features based on T2-weighted magnetic resonance imaging. Front Neurol. (2020) 11:248. doi: 10.3389/fneur.2020.00248
24. Adeli, E, Wu, G, Saghafi, B, An, L, Shi, F, and Shen, D. Kernel-based joint feature selection and max-margin classification for early diagnosis of Parkinson's disease. Sci Rep. (2017) 7:41069. doi: 10.1038/srep41069
25. Bidesi, NSR, Vang Andersen, I, Windhorst, AD, Shalgunov, V, and Herth, MM. The role of neuroimaging in Parkinson's disease. J Neurochem. (2021) 159:660–89. doi: 10.1111/jnc.15516
26. Parent, A, and Hazrati, LN. Functional anatomy of the basal ganglia. II. The place of subthalamic nucleus and external pallidum in basal ganglia circuitry. Brain research. Brain Res Rev. (1995) 20:128–54. doi: 10.1016/0165-0173(94)00008-d
27. Hua, S, Reich, SG, Zirh, AT, Perry, V, Dougherty, PM, and Lenz, FA. The role of the thalamus and basal ganglia in parkinsonian tremor. Mov Dis. (1998) 13:40–2. doi: 10.1002/mds.870131307
28. Hutchison, WD, Lozano, AM, Tasker, RR, Lang, AE, and Dostrovsky, JO. Identification and characterization of neurons with tremor-frequency activity in human globus pallidus. Exp Brain Res. (1997) 113:557–63. doi: 10.1007/pl00005606
29. Li, G, Zhai, G, Zhao, X, An, H, Spincemaille, P, Gillen, KM, et al. 3D texture analyses within the substantia nigra of Parkinson's disease patients on quantitative susceptibility maps and R2(∗) maps. NeuroImage. (2019) 188:465–72. doi: 10.1016/j.neuroimage.2018.12.041
30. Zhang, S, Zi, S, Xiong, S, Peng, H, Hu, K, and He, H. Apathy following bilateral deep brain stimulation of subthalamic nucleus and Globus pallidus internus in Parkinson's disease: a Meta-analysis. Parkinsons Dis. (2022) 2022:4204564. doi: 10.1155/2022/4204564
31. Santos, AN, Kherif, F, Melie-Garcia, L, Lutti, A, Chiappini, A, Rauschenbach, L, et al. Parkinson's disease may disrupt overlapping subthalamic nucleus and pallidal motor networks. NeuroImage Clin. (2023) 38:103432. doi: 10.1016/j.nicl.2023.103432
32. Sigirli, D, Ozdemir, ST, Erer, S, Sahin, I, Ercan, I, Ozpar, R, et al. Statistical shape analysis of putamen in early-onset Parkinson's disease. Clin Neurol Neurosurg. (2021) 209:106936. doi: 10.1016/j.clineuro.2021.106936
33. Schwab, NA, Tanner, JJ, Nguyen, PT, Schmalfuss, IM, Bowers, D, Okun, M, et al. Proof of principle: transformation approach alters caudate nucleus volume and structure-function associations. Brain Imaging Behav. (2015) 9:744–53. doi: 10.1007/s11682-014-9332-x
34. Shen, W, Zhai, S, and Surmeier, DJ. Striatal synaptic adaptations in Parkinson's disease. Neurobiol Dis. (2022) 167:105686. doi: 10.1016/j.nbd.2022.105686
Keywords: PD, radiomics, artificial intelligence, T1-weighted, MRI
Citation: Li C, Hui D, Wu F, Xia Y, Shi F, Yang M, Zhang J, Peng C, Feng J and Li C (2024) Automatic diagnosis of Parkinson’s disease using artificial intelligence base on routine T1-weighted MRI. Front. Med. 10:1303501. doi: 10.3389/fmed.2023.1303501
Edited by:
Minjeong Kim, University of North Carolina at Greensboro, United StatesReviewed by:
Cheng-Yu Wei, Chang Bing Show Chwan Memorial Hospital, TaiwanWeisheng Zhuang, Henan Provincial People's Hospital, China
Andrea Quattrone, University Magna Graecia, Italy
Copyright © 2024 Li, Hui, Wu, Xia, Shi, Yang, Zhang, Peng, Feng and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chuanming Li, licm@cqu.edu.cn; Junbang Feng, junbangfeng@163.com; Chao Peng, 523945695@qq.com
†These authors have contributed equally to this work