2.1 Patients
Our study protocol was approved by the Ethics Review Committee of Chiba University (Approval reference No. 711 and 982). All procedures were carried out in accordance with the relevant guidelines and regulations. Written informed consent was obtained from participants and their parent and/or legal guardian. Overall, 169 patients were recruited at our institute for this study. Patients with a family history of PD, those with PARK2 mutations, and those treated with deep brain stimulation were excluded. Patients who did not undergo three-dimensional T1-weighted imaging on MRI or patients who showed brain tumors, hydrocephalus, traumatic lesions, and cerebrovascular disease on MRI were also excluded.
In total, 103 participants (61 men and 42 women; mean age ± standard deviation (SD): 68.1 ± 10.5 years; mean disease duration 7.4 ± 4.8 years) participated in this study. All patients fulfilled the PD criteria defined by the Movement Disorder Society Clinical Diagnostic Criteria for Parkinson's disease 16. The patients were evaluated using standard clinical tests, including i) disease duration, ii) disease severity according to the Unified Parkinson's Disease Rating Scale (UPDRS) part Ⅲ 17, iii) the Mini-Mental State Examination scores (MMSE) 18, and iv) medication use evaluated as levodopa-equivalent dose (LED) 19. 108 HS (48 men and 60 women; mean age ± SD: 65.9 ± 10.1 years) without any mental or neurological disorders were also recruited.
2.2 DNA
DNA was extracted from peripheral blood mononuclear cells from blood samples using a MagNA Pure Compact System (Roche Diagnostics, Penzberg, Germany). Genotyping for polymorphisms in DRD2 rs1800497 was performed by directly sequencing regions using High-Resolution Melting (HRM) analysis. Primers for HRM analysis were made using LightCycler Probe Design Software 2.0 (Roche), and HRM polymerase chain reaction (HRM-PCR) was performed using LightCycler 480 Instrument (Roche). HRM curve analysis was conducted with LightCycler 480 Gene Scanning Software (version 1.5) 20. PCR products were purified on a BigDye Xterminator (Life Technologies, Carlsbad, CA) and read by Applied Biosystems 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA). The A1 carriers were defined as having an AA genotype (15.5%) and AG genotype (59.2%), and the A2 homo carriers were defined as having a GG genotype (25.2%) 21. The Hardy–Weinberg equilibrium of the rs1800497 distribution had a p-value of >0.01. We did not perform DNA screening tests on HS.
2.3 MRI
All participants underwent MRI scans using a 1.5-T scanner (Siemens Medical System, Erlangen, Germany) or a 3-T scanner (Philips Medical Systems, Best, Netherlands). In this study, three-dimensional T1-weighted imaging was acquired by the following parameters: TR = 7 ms, TE = 3 ms, flip angle = 15°, field of view = 24 cm, section thickness = 1.4 mm, and resolution = 0.9 × 0.9 × 1.0 mm³ by the 1.5-T scanner, and TR = 8 ms, TE = 3 ms, flip angle = 15°, field of view = 22 cm, section thickness = 1.0 mm, and resolution = 0.9 × 0.9 × 0.7 mm³ by the 3-T scanner.
Statistical Parametric Mapping (SPM12) was used for image processing for VBM. The 3D-T1 images were segmented into gray matter, white matter, and cerebrospinal fluid images, and the Diffeomorphic Anatomical Registration Through Exponential Lie Algebra (DARTEL) toolbox on SPM12 was used for modulation 22. By incorporating the DARTEL template created by the affine transformation, the grey matter segments were warped into standard Montreal Neurological Institute space. Normalization of grey matter volumes was performed using Jacobian determinants from the registration step. A 9-mm Gaussian kernel was used to smoothen the modulated grey matter and white matter images.
2.4 Statistics
Analysis of variance (ANOVA) was performed to evaluate differences in age, sex, and total brain volume (TBV) between the A1 carriers, A2 homo carriers, and HS. Chi-square tests were used to compare the proportions of categorical variables, such as sex, between these groups. The Wilcoxon rank sum test was used as a non-parametric test to compare the averages of continuous variables, such as age and TBV, between A1 carriers and A2 homo carriers.
SPM12 software was used to perform the VBM analyses. A full-factorial analysis of covariance (ANCOVA) was used to evaluate morphological differences in grey matter (GM) between groups on a voxel-by-voxel basis throughout the whole brain. Age, gender, MRI scanners, and TBV were added as covariates of no interest to eliminate the effects of confounding factors. The initial voxel threshold was set at 0.001 uncorrected. The cluster-level and peak-level thresholds were set at P < 0.05 corrected for family-wise error (FWE). Region of interest (ROI) analyses were also conducted to calculate the caudate volumes. The caudate ROIs were obtained using automated anatomical labeling from the PickAtlas toolbox 23. According to a previous study 24, ROI measurements should be conducted in consideration of the intracranial volume and age. Therefore, the caudate volume ratios (caudate volume divided by TBV) were calculated 25. The target ROI was also determined as a sphere with a 6-mm radius from the opposite site of the peak coordinate obtained from the results of VBM analyses. ANOVA was used to compare group differences. We used R version 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria) for statistical analyses except for SPM12. We chose a p-value of <0.05 to indicate significant difference.