Patient-specific multi-modal modeling uncovers neurotransmitter receptor involvement in motor and non-motor axes of Parkinson’s disease ========================================================================================================================================= * Ahmed Faraz Khan * Quadri Adewale * Sue-Jin Lin * Tobias R. Baumeister * Yashar Zeighami * Felix Carbonell * Nicola Palomero-Gallagher * Yasser Iturria-Medina ## Abstract Multi-systemic neurodegeneration in Parkinson’s disease (PD) is increasingly acknowledged, involving several neurotransmitter systems beyond the classical dopaminergic circuit and resulting in heterogeneous motor and non-motor symptoms. Nevertheless, the mechanistic basis of neuropathological and symptomatic heterogeneity remains unclear. Here, we use patient-specific generative brain modeling to identify neurotransmitter receptor-mediated mechanisms involved in PD progression. Combining receptor maps with longitudinal neuroimaging (PPMI data), we detect a diverse set of receptors influencing gray matter atrophy, microstructural degeneration, and dendrite loss in PD. Importantly, identified receptor mechanisms correlate with symptomatic variability along two distinct axes, representing motor/psychomotor symptoms with large GABAergic contributions, and cholinergically-driven visuospatial dysfunction. Furthermore, we map cortical and subcortical regions where receptors exert significant influence on neurodegeneration. Our work constitutes the first personalized causal model linking the progression of multi-factorial brain reorganization in PD across spatial scales, including molecular systems, accumulation of neuropathology in macroscopic brain regions, and clinical phenotypes. Keywords * neurotransmitter receptors * multimodal neuroimaging * Parkinson’s disease * whole-brain computational model * personalized medicine ## Introduction Parkinson’s disease (PD) has traditionally been primarily associated with a nigrostriatal dopamine deficit resulting in the characteristic motor symptoms of tremor, rigidity, and bradykinesia. However, the involvement of other brain circuits is now widely recognized [1], and the majority of patients also present numerous non-motor symptoms such as dementia, depression, sleep disorders, or apathy [2]. For this multi-system disease with significant inter-patient heterogeneity in pathology, symptoms and treatment response [3] [4] [5], consistent links between genetic, neuropathological and clinical subtypes remain elusive [6]. With no cure [7], symptomatic pharmacological treatment (e.g. levodopa) is at best partially effective [8] and may result in undesired side effects with chronic administration [9]. Given that diagnostic accuracy in untreated or medication non-responder PD patients is as low as 26% [10], an improved understanding of biological mechanisms and potential therapeutic targets underlying PD heterogeneity is imperative to bridging the treatment gap in PD [11] [12] [13]. Neurotransmission underlies many disease-related mechanisms as well as pharmacological response [8] [14]. Regional variability in neurotransmitter receptor gene expression correlates with altered macroscopic interactions such as neurovascular decoupling [15] and structural-functional decoupling [16]. Many non-dopaminergic nuclei are affected in PD [17], [18], with specific neurotransmitter systems linked to symptoms such as cholinergic freezing of gait and dementia [19], serotonergic depression and tremor [20], and adrenergic postural symptoms [21]. The dual syndrome hypothesis of PD [18] proposes a dichotomy between dopamine-mediated fronto-striatal executive impairment and a cholinergically-mediated prodromal visuospatial dementia. To better characterize the role of neurotransmission in mediating neurodegenerative brain reorganization, an integrative model linking multiple receptor systems, macroscopic brain reorganization and clinical symptoms would be essential. However, we are limited by the absence of whole-brain spatial distribution maps of neurotransmitter receptors in PD patients [8]. On the other hand, neuroimaging supports the multi-factorial and heterogeneous view of PD [22]. Various modalities are routinely used to support differential diagnosis [23] [11] [24] and evaluate treatment effects [25]. Multi-modal modeling of neuroimaging alterations can elucidate the temporal ordering, disease trajectories, and interactions of various pathologies in neurodegeneration [26] [27], and link these macroscopic observations with underlying genetic and transcriptomic determinants [28]. Multifactorial causal modeling (MCM) is a mechanistic modeling approach that is able to identify contributions of interacting factors to longitudinal changes [29], which can be used in a personalized medicine context to design optimal therapeutic interventions [30]. Combining multi-modal neuroimaging with spatial distribution templates of 15 neurotransmitter receptors from post-mortem autoradiography [31] in an MCM-based approach significantly improved the explanation of degenerative changes in individual patients’ neuroimaging data, and linked specific receptor-pathology interactions to clinical symptoms in Alzheimer’s disease (AD) [32]. Furthermore, this approach was able to estimate individualized receptor alterations based on inter-subject differences in receptor-neuroimaging interactions. Here, we extend previous molecular-phenotypic PD characterizations in four fundamental ways: i) combining spatial distribution templates of 15 key neurotransmitter receptors derived from post-mortem autoradiography [31] with longitudinal neuroimaging data in a personalized modeling framework that infers patient-specific receptor-mediated alterations (N=71, PPMI data), ii) demonstrating the improved ability of receptor-enriched multifactorial causal modeling (re-MCM) to explain imaging-measured neurodegeneration and identify consistent mechanistic changes across patients, iii) characterizing inter-patient receptor-based heterogeneity, specifically linking mechanistic alterations to two main axes of motor, cognitive and psychiatric symptoms, iv) quantitatively mapping brain regions with high receptor influence on PD neurodegeneration. This work represents the first attempt to integrate a broad range of neurotransmitter receptors, multi-modal neuroimaging and clinical symptoms for the identification of the molecular pathways underlying personalized treatment needs in PD. ## Results ### Model-based approach to inferring personalized neurotransmitter receptor alterations To characterize neurotransmitter receptor contributions to the multifaceted neurodegenerative processes of PD, we fit receptor-informed individualized generative computational models to the longitudinal alterations of 6 biological factors. Each biological factor is associated with neurodegeneration in PD, namely atrophy, dysregulated functional activity, dopaminergic deficiency directed and microstructural damage, and dendrite loss, represented by the neuroimaging-derived measures of gray matter density (GM), fractional amplitude of low frequency fluctuations (fALFF), dopamine transporter SPECT (DAT-SPECT), fractional anisotropy (FA), mean diffusivity (MD), and t1/t2 ratio [33] [34]. Neuroimaging data was acquired over multiple imaging scans for N=71 PD patients (PPMI data, *Methods: Data description and processing*). In addition, regional densities for 15 neurotransmitter receptors (from glutamatergic, GABAergic, cholinergic, adrenergic, serotonergic, and dopaminergic families) were derived from averaged templates (*Methods: Data description and processing: Receptor densities and brain parcellation),* and anatomical connectivity was estimated from the high-resolution Human Connectome Project template (HCP-1065; *Methods: Anatomical connectivity estimation*). The *neurotransmitter receptor-enriched multifactorial causal model* (re-MCM; Fig. 1) decomposes the spatiotemporal evolution of pathology of multiple biological factors into localized receptor- and network-mediated effects. Model parameters explicitly represent distinct biological mechanisms, namely i) direct and ii) receptor-mediated pairwise interactions between imaging-derived biological factors (dopaminergic deficiency, functional activity, microstructural damage, dendrite density, and atrophy), iii) effects of local neurotransmitter receptor densities on factor-specific longitudinal deterioration, and iv) spreading of pathology to and from anatomically-connected regions. Notice that, in the absence of true personalized longitudinal receptor imaging, model weights of specific receptor-mediated biological mechanisms compensate to fit individualized trajectories of neurodegeneration. Thus, inter-subject variability in model weights serves as a proxy for the corresponding receptor densities or receptor-pathology interactions. Specifically, i) model fit being significantly improved by the inclusion of healthy aged receptor templates validates their application to this clinical population, ii) biological mechanisms that are statistically stable across subjects represent mechanistic pathways shared by all PD patients in our cohort, iii) inter-patient co-variability between biological mechanisms and clinical symptoms represents overlapping disease processes (Fig. 1b), and iv) inter-region variability in the model fit improvement due to receptor templates can identify regions differentially affected by neurotransmitter receptor alterations in PD (Fig. 1c). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F1) Figure 1: Neurotransmitter receptor-enriched multifactorial causal modeling. **a)** Each patients’ longitudinal pathological progression is decomposed into local effects due to: i) direct influence of every imaging-derived biological factor (e.g., atrophy on resting state functional activity), ii) receptor density distribution (e.g., D1 receptor density on DAT loss), and iii) receptor-pathology interactions (e.g., D1 receptors × DAT interactions on functional activity), in addition to iv) network-mediated inter-region propagation. Combining this data across (NROI=95) brain regions and multiple visits results in a multivariate regression problem to identify the patient-specific parameters {α}. **b)** Decomposing the covariance matrix of patients’ model-derived biological mechanism weights and clinical scores (specifically, the rates of decline of composite clinical scores; *Methods: Clinical scores*) identifies multivariate axes of receptor-factor interactions that are robustly correlated with the severity of combinations of clinical symptoms in PD (*Methods: Biological parameters and relationship with cognition)*. **c)** The regional contributions of receptor interactions to neurobiological changes are estimated by a feature importance analysis. We fit individualized models for every biological factor with and without each receptor map, and performed permutation tests on the improvement in regional model residuals due to the inclusion of receptor maps. The resulting z-scores are the regional influence of receptors on each target biological factor model. ### Neurotransmitter receptor maps significantly improve the explainability of multi-factorial brain reorganization in PD Before proceeding to identify relevant model-derived biological mechanisms in PD, we first aimed to validate that re-MCM robustly fits patient-specific neuroimaging data. For each of the 6 biological factors and all subjects (N=71), we calculated the coefficient of determination (R2) as a measure of the data variance explained. On average, re-MCM explained 74% ± 18% of the variance in rate of pathology accumulation (Fig. 2a), although model fit varied by biological factor, with dopaminergic (DAT-SPECT; 80% ± 13%) and dendrite loss (t1/t2 ratio; 80% ± 12%) being explained better than gray matter atrophy (GM; 58% ± 14%), or microstructural damage (MD; 70% ± 14%, and FA; 0.74 ± 0.13). For validation, we repeated the model-fitting without receptor-pathology interactions or direct local receptor density effects. On average, neuroimaging-only models without receptor data explained 52% ± 20% of the variance in neuroimaging rate of change (Fig. 2b). On average, the inclusion of receptor templates improves the data variance explained by 42.3%. Dopaminergic loss (DAT-SPECT) was the least improved by the addition of receptor maps, with imaging-only models explaining 60% ± 17%, a drop of 20% of variance on average compared to the full re-MCM. On the other hand, gray matter atrophy (GM: 22% ± 17% variance explained without receptor maps) was the most reliant on receptor data. While DAT-SPECT scans themselves already image the density of presynaptic dopaminergic transporters, gray matter atrophy models benefit more from regional differentiation based on receptor expression. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F2) Figure 2: Contribution of receptor distributions to explaining multimodal brain reorganization in PD. The improvement in modeling the accumulation of pathology was evaluated in terms of i) the additional explanatory power due to receptor information, and ii) the significance of true receptor maps compared to null distributions. The histograms show the distribution of the coefficient of determination (R2) of N=71 individual models of longitudinal neuroimaging changes including **(a)** and excluding **(b)** receptor predictors. Notably, including receptor terms improves model fit for all biological factors, although to varying extents. **(c)** Subject-wise F-tests between models with and without receptor maps (113 and 8 parameters, respectively) show proportions of subjects for whom the F-statistic is above the critical threshold (red dotted line). This critical threshold corresponds to a statistically significant (P<0.05) improvement due to the receptor terms in the re-MCM model, accounting for the increase in adjustable model parameters. Furthermore, to validate the benefit of the receptor templates over randomized null maps, re-MCM models were fit with 1000 spatially-shuffled receptor maps for each subject. The p-value of the model fit (R2) using true receptor templates compared to the distribution of R2 of models using randomized templates was calculated for each subject. **(d)** Proportion of subjects for whom the true receptor maps resulted in a statistically significant improvement in model fit (P<0.05; red dotted line). Figure 2c presents the improvement in each participant’s model fit due to receptor mechanisms, compared to the restricted, neuroimaging-only models. Accounting for the increased model size from 8 to 113 parameters, the F-statistics of 80.3% (MD) to 100% (DAT-SPECT) of patients is significant (p<0.05 red dotted line in Fig. 2c). We then performed a permutation test for the significance of the informativeness of receptor maps, by randomly shuffling each receptor map across brain regions 1000 times, and fitting the re-MCM with each set of permuted maps. The resulting distribution of model fit (R2) was used to calculate significance levels for re-MCM with true receptor data from Figure 2a. For each biological factor, we plotted the number of subjects with significantly better model fit (p<0.05) compared to the null distribution in Fig. 2d. Notably, nearly all patients’ biological factor models are significantly improved by the inclusion of receptor maps, with the exception of undirected microstructural damage (MD; 67.6% or 48 subjects). Across all participants, Fisher’s method gives χ2 statistics in the range of 800< χ2 < 2300 (depending on the biological factor), corresponding to a near-zero combined P-value. These analyses validate the use of averaged receptor templates in patient-specific PD models. ### Identifying stable neurobiological mechanisms and receptor-pathology interactions in PD We proceeded to identify biological mechanisms consistently involved in dopaminergic, functional and structural brain alterations in PD. For this, 99% confidence intervals for each re-MCM parameter across all patients were calculated and used to identify stable predictors. Since all predictors were standardized before data fitting, model weights are the relative effect sizes of different biological mechanisms on the rate of change of their target biological factor over the course of PD progression. Specifically, these neurobiological mechanisms are i) direct effects of local pathology, ii) direct effects of local receptor densities, iii) local receptor-pathology interactions, and iv) network propagation of pathology (*Methods: Receptor-Enriched Multifactorial Causal Model*). Figure 3a shows the relative effective sizes of stable biological mechanisms. The most influential stable predictors of each biological factor’s rate of change are the direct effects of local alterations to the same modality. Propagation of pathology along the structural connectome is also a minor yet stable predictor for all data modalities except functional activity (fALFF) and directed microstructural damage (FA), with a much lower effect than the local evolution of neurodegeneration. Notably, from Figure 3b, functional brain alterations (fALFF) do not appear to drive structural alterations (GM and MD), instead interacting bidirectionally with dendritic density (t1/t2). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F3) Figure 3: Receptor-mediated interactions explaining longitudinal neurodegeneration in PD. **a)** Statistically stable biological mechanisms in PD show significant receptor-mediated contributions. The angle of each outer sector is proportional to the mean weight of each stable (99% confidence interval) re-MCM model weight across the PD patients. The inner sectors represent the 6 modeled biological factors. Within each factor, the intermediate sectors represent the neurotransmitter system involved, while the outer sector consists of the specific two-way receptor-pathology interactions or direct predictor terms in the model. Notably, biological factors may appear as both model predictors (outer sector) and targets (inner sector). **b)** Effect size (number of chords) of statistically stable interactions between any pair of biological factors modeled in PD. Nevertheless, local interactions between imaging-based biological factors, whether direct or receptor-mediated, constitute a significant driver of PD neurodegeneration in all cases, and form a complex network with potentially bidirectional influences (Fig. 3b). While comparatively smaller for functional activity, dopaminergic transporter density and directed microstructural integrity (FA), receptor-mediated interactions constitute approximately half the model effects for gray matter atrophy (GM), overall microstructural integrity (MD) and dendrite density (t1/t2). We observed that a relatively sparse set of receptors is involved in stable interactions for each biological factor (Fig. 4). The muscarinic M2 and nicotinic α4β2 receptors contribute significantly to gray matter atrophy, neuronal activity dysfunction, and dopaminergic loss. The Bz site is also prominently associated with neuronal activity dysfunction and dopaminergic loss. The serotonergic 5HT2 receptor is involved in functional and undirected microstructural alterations, while glutamatergic effects are marked by NMDA affecting gray matter atrophy, AMPA and kainate affecting directed microstructure and kainate affecting dendrite density, respectively. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F4) Figure 4: Receptors mediating degenerative alterations to different macroscopic biological factors in PD. The combined statistically stable model effects of each receptor type on each biological factor are shown. The muscarinic M2 and nicotinic α4β2 receptors contribute significantly to gray matter density, neuronal activity and dopamine transporter alterations. The Bz site is prominently associated with activity and dopamine transporter alterations. The serotonergic 5HT2 receptor is involved in functional and microstructural (MD) alterations, while glutamatergic effects are marked by NMDA affecting gray matter atrophy, AMPA and kainate affecting directed microstructural damage (FA) and kainate affecting dendrite density (t1/t2), respectively. Notably, the D1 receptor distribution is relatively homogeneous and not marginally informative in the presence of DAT imaging. Generally, the dopaminergic, cholinergic, serotonergic, glutamatergic and GABAergic systems broadly affect (micro-)structural alterations (GM, MD and t1/t2). Serotonergic mechanisms are most associated with undirected microstructural alterations (MD), and secondarily dysfunctional neural activity (fALFF). Cholinergic receptors are prominent predictors of atrophy, microstructural damage and loss of dendrites (GM, MD and t1/t2), with minor influence on functional activity and dopaminergic transporter density. Glutamatergic receptors have a moderate influence across structural modalities (GM, MD, FA and t1/t2). GABAergic influence is minor yet stable across functional (fALFF and SPECT) and (micro-)structural (MD and t1/t2) modalities. Adrenergic and dopaminergic receptors are the least involved in stable neurobiological mechanisms, with α2 adrenergic receptor modulating directed microstructural damage (FA), and the D1 dopaminergic receptors mediating the effect of atrophy on microstructure (MD). For atrophy (GM), functional activity (fALFF) and microstructure (MD) models, the direct effects of specific receptor density maps reflect local susceptibility to neurodegeneration. The densities of the muscarinic M2 and nicotinic α4β2 cholinergic receptors help explain inter-region variability in the rate of gray matter atrophy, while M2 and the serotonergic 5HT2 receptor densities are stable predictors of both altered activity (fALFF) and microstructural damage (MD). ### Two axes of receptor-pathology alterations underlie clinical symptoms in PD To link model-derived receptor-mediated neurobiological mechanisms with clinical presentation in PD, we identified shared axes of covariance between re-MCM-derived biological mechanisms and motor, non-motor, cognitive and psychiatric symptoms (*Methods, Clinical scores*). Singular value decomposition (SVD) across all patients (N=71) was used to identify multivariate and overlapping relationships between identified biological parameters and clinical symptoms (*Statistical analysis: Covariance of biological mechanisms with clinical symptoms*). Two SVD components were relevant based on permutation tests, explaining 48.4% (P<0.001, FWE-corrected) and 13.2% (P<0.07, FWE-corrected) of the population co-variance, respectively. Projections of biological mechanisms and clinical scores to these components show moderate to high correlations of r=0.70 (P<10-12; Fig. 5a) and 0.86 (P<10-20; Fig. 5b). ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F5) Figure 5: Two axes of covariance between biological mechanisms and symptom severity in PD. **a)** Based on a permutation analysis, two latent SVD components were significant or near-significant, explaining 48.4.0% (P<0.001, FWE-corrected) and 13.2% (P<0.07, FWE-corrected) of the covariance respectively. **a,b)** High correlations of r=0.70 (P<0.001) and 0.86 (P<0.001), between the projections of statistically stable biological mechanisms and rates of clinical decline onto the latent space were observed. **c,d)** Bootstrap ratios of each clinical assessment to the two latent components, providing a relative ranking of motor, nom-motor, psychiatric and cognitive domains. These saliences are proportional to the contribution of each term relative to every other term, for example showing that MDS-UPDRS scores, UPDRS, SDM and HVLT scores are the top contributors to PC1. Interestingly, the first component (PC1; Fig. 5c) largely corresponds to variance of the MDS-UPDRS Parts 1-3 scores (composed of cognitive, psychiatric and motor aspects of daily living, as well as a motor exam), and SDM (assessing attention, perceptual speed, motor speed, and visual scanning [35]). On the other hand, the second component (PC2; Fig. 5d) is associated with the BJLOT (visuospatial judgment), LNS (working memory), STAIAD (anxiety) and the GDS (depression in older adults). The statistically stable biological mechanisms contributing to each axis are summarized in Figure 6. Both components show that inter-subject symptom variability is associated with multiple receptor-mediated biological mechanisms and neuropathological changes. PC1 is largely driven by GABAergic alterations (explaining 5.97% of the total covariance via this component), although glutamatergic (4.85%), cholinergic (4.77%), and serotonergic (3.77%) alterations are also prominent. PC2 is instead associated primarily with cholinergic alterations (1.74%), although GABAergic (1.24%) and glutamatergic (1.19%) alterations also play a role. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F6) Figure 6: Distinct combinations of receptor-mediated interactions are associated with the two axes of clinical symptoms. Biological mechanisms correlated with clinical severity in PD via the **a)** motor/psychomotor and **b)** visuospatial/memory/psychiatric axes are plotted. Representing the effects of receptor densities, local pathology, receptor-pathology interactions, and network propagation of pathological factors, combinations of patient-specific mechanisms co-vary with specific clinical symptoms. For each mechanism, the angle is proportional to the percentage of mechanistic-clinical covariance explained. The outer sector contains the specific mechanisms, while the middle sector is grouped by receptor families and the inner sector by target biological factor. While the local (regional) evolution of pathology in each considered biological factor and its network propagation are prominent stable predictors of PD neurodegeneration (Fig. 3), the influence of these mechanisms does not co-vary significantly with symptom severity. Instead, we find a broad array of receptors with clinical effects along both PCs, as shown in Fig 5. For example, the primarily motor symptoms of PC1 are associated with inter-subject variability in glutamatergic and GABAergic interactions affecting functional activity (fALFF), microstructural integrity (MD) and dendrite density (t1/t2). In contrast, the more psychiatric and visuospatial dysfunction of PC2 is associated more with inter-subject variability in cholinergic interactions affecting microstructure (MD) and dendritic density (t1/2), as well as a larger adrenergic component. ### Obtaining PD receptor influence maps Finally, we inferred the degree of receptor influence on multi-modal PD neurodegeneration at different brain regions. For each receptor, we fit individualized, single receptor-enriched models, and compared their ability to explain the accumulation of pathology at each brain region with restricted, neuroimaging-only models (see *Statistical analysis: Regional analysis*). At each brain region, we studentized residuals across all patients, with each residual representing the unexplained pathology in a region at a given imaging visit. For all regions, we computed the Wilcoxon rank sum statistics of the population residuals from the two models, and repeated the model-fitting procedure with 1000 randomly shuffled receptor maps to obtain a null distribution of Wilcoxon statistics. We used this permutation test to filter brain regions with significant residual improvements (P<0.05) over the null distributions. In Figure 7, we summarize the receptor influence maps for the top 4 receptor-pathology pathways (Fig. 3a): 5HT2 and M2 on microstructural alterations (MD), α4β2 on gray matter atrophy (GM), and kainate on dendrite density (t1/t2). Receptor influence maps for all biological factors are presented in Supplementary Figures S1-S6. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/22/2023.03.22.23287578/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2023/03/22/2023.03.22.23287578/F7) Figure 7: Model-derived maps of receptor influence on PD neurodegeneration. We compared the (a) receptor density and (b) influence maps. Influence maps represent the population-wide improvement in modeling accuracy at each region due to the inclusion of receptor density maps and receptor-pathology interactions as model predictors for each PD patient. Among other regions, the 5HT2 receptor most prominently influences microstructure (MD) in the anterior and medial thalamus, left posterior cingulate region (Brodmann area 31), anterior prefrontal cortex, left primarily motor cortex, right premotor cortex and supplementary motor area (Brodmann area 6). The muscarinic M2 receptor influences microstructural alterations in the somatosensory cortex, left distal visual area V3d, right primary motor cortex, left hippocampus (CA), right primary somatosensory cortex (Brodmann area 2), lateral prefrontal cortex (Brodmann areas 46 - left and 47 - right), and entorhinal cortex (Brodmann areas 36-right and 37-left). The nicotinic α4β2 receptor influences gray matter atrophy in the (left and right) thalamus, primary somatosensory cortex (Brodmann area 2), right temporal inferior parietal area, left caudate nucleus and entorhinal cortex (left Brodmann region 22). Kainate influences dendrite density in a broad set of regions, focused on the thalamus, visual areas (V1, V2 and the ventral parts of V3 and V4 in the right hemisphere, and V1 and ventral V4 in the left hemisphere), and prefrontal areas. Across biological factors, glutamatergic receptors contribute significantly to explaining neurodegeneration in fronto-temporal regions (Supplementary Fig. S1). Particularly, both AMPA and kainate receptors contribute strongly to most factors (with the exception of dopamine transporter loss) in frontal regions. The influences of GABAA receptors, GABAB receptors and the benzodiazepine binding site (Bz site) generally follow their distribution (Supplementary Fig. S2), peaking at visual, visual-parietal and fronto-temporal areas, respectively. Notably, dendrite loss is most pronounced at subcortical and fronto-temporal regions for all GABAeric receptors. ## Discussion The complex pathophysiology of PD involves multiple difficult-to-map neurotransmitter systems, and the selective vulnerability of various non-dopaminergic nuclei [4]. We apply a novel personalized brain modeling approach that combines longitudinal neuroimaging data and clinical assessments with averaged spatial receptor templates, to infer the previously uncharacterized roles of receptor-mediated interactions in PD neurodegeneration and symptomatic heterogeneity. This constitutes the first attempt at a causal model linking the progression of multi-modal neurodegeneration in PD across spatial scales, including molecular mechanisms, macroscopic accumulation of neuropathology, and clinical phenotype. Although dopaminergic neuroimaging is common [36], the expense of PET imaging and the lack of suitable *in vivo* radioligands have impeded the study of other receptor alterations in a PD population. Our method circumvents this limitation by inferring the importance of receptor interactions in individualized models of brain reorganization. While we used autoradiography-derived templates of receptor density, receptor gene expression may be used as a proxy [37]. For example, the Allen Human Brain Atlas ([http://human.brain-map.org](http://human.brain-map.org)) gene expression template has been used to identify transcriptomic pathways mediating neurodegeneration in Alzheimer’s disease (AD) [38]. However, several translational and trafficking steps separate gene expression and synaptically integrated receptors. Although receptor densities and gene expression are correlated for selected receptor subunit genes and across certain cytoarchitectonically-defined regions [39], this is not universally true [40]. Low correlations are also observed between gene expression and *in vivo* PET imaging of dopamine transporters [41]. On the other hand, our averaged receptor templates are correlated with neurobiological processes such as drug-induced cerebral blood flow changes [42]. Furthermore, *in vitro* autoradiography allows access to a broader class of receptors at a sub-millimeter resolution (as low as 0.3mm slice width per receptor [31]) compared to PET with its theoretical bound of ~2mm spatial resolution [43]. Future work will extend the presented results with voxel-scale whole brain receptor maps rather than macroscopically averaged values. We incorporated several neuroimaging-derived measures sensitive to PD progression [44], from structural MRI-based gray matter density (GM) and dendrite density (t1/t2 ratio), diffusion-based measures of microstructural integrity (MD and FA) [45], functional neuronal activity (fALFF) and presynaptic dopamine transporter availability (DAT-SPECT). Resting-state fMRI-derived metrics such as fALFF can distinguish PD patients from controls [46], with fALFF being able to explain up to 25% of variability in MDS-UPDRS scores [47]. While initially proposed as a quantitative measure of demyelination from routine MRI scans, t1/t2 ratio has since been demonstrated to have a stronger correlation with dendritic density [33] [34], particularly relevant to synaptic integrity and receptor activity. Furthermore, our flexible modeling approach can be extended to incorporate other relevant modalities. Although receptor maps were averaged from 4 neurologically healthy aged brains, earlier work has demonstrated their informativeness to explaining the accumulation of neuropathology in other cohorts, namely healthy aged subjects, mildly cognitively impaired subjects and AD patients from the Alzheimer’s disease Neuroimaging Initiative (ADNI) [32]. Extending this validation to the PPMI cohort, we note an approximately 42.3% improvement in the explanation of neuropathology accumulation in receptor-enriched models. These improvements are statistically significant for well over 90% of subjects (P<0.05 in both F-tests and permutation tests; Fig. 2c,d) for all biological factors with the exception of undirected microstructural damage (MD). For a third of all subjects, the improvement in model fit of undirected microstructure was not significantly better than permuted null distributions of receptors. While diffusion MRI can be sensitive to aspects of gray matter microstructure [48] [49], it is less accurate than in white matter due to the heterogeneity of tissues and their (lack of) organization [50]. Yet, despite the limitation of partial volume effects in gray matter ROIs [51], receptor-enriched models fit longitudinal alterations to microstructure reasonably well (average r2 = 0.70 for undirected MD and r2 = 0.74 for directed FA; Fig. 2). Differential neurotransmitter and receptor expression may underpin the selective vulnerability of several neuronal populations, from the dopaminergic substantia nigra to the adrenergic locus coeruleus and serotonergic raphe nuclei, and their cortical projections [52]. Furthermore, PD neurodegeneration may alter both the spatial distributions as well as functional interactions of specific dopaminergic and non-dopaminergic receptors, with symptomatic consequences [14]. In our mechanistic modeling framework, each model weight is interpretable as the importance of specific neurobiological mechanism. Receptors contribute to neurodegeneration in re-MCM either as i) direct effects representing regional susceptibility to neurodegeneration based on receptor expression, or ii) receptor-mediated interactions involving a source and target biological factor. Additionally, biological factors have i) local effects on themselves and other factors, and ii) intra-factor network effects due to propagation of pathology along the structural connectome. Lacking inter-subject variability in receptor data, our model compensates by assigning weights differently across subjects. Consistent trends in model weights reflect the importance of the corresponding neurobiological mechanism across the PD population, while co-variability with symptoms suggests clinical relevance. First, we identified specific mechanisms affecting neurodegeneration across the PD cohort (Fig. 3). We observed a complex network of interactions between biological factors, with distinct receptor profiles affecting each factor. The large contributions of receptor-mediated inter-factor interactions (Fig. 3a) supports the multi-system view of PD. Fewer receptors are statistically stable predictors of longitudinal changes to functional activity (fALFF), directed microstructural damage (FA) and dopaminergic neurotransmission (SPECT), while gray matter atrophy (GM), dendrite density (t1/t2 ratio) and undirected microstructural changes (MD) show greater influence from a more diverse set of receptors. Notably, the D1 receptor map is not a stable predictor of DAT alterations. While presynaptic DAT density and postsynaptic dopaminergic receptor distributions are strongly related under normal conditions, they may be affected differently by disorders. For example, while D2 receptor availability is reduced in alcoholism, DAT availability is preserved [53]. In PD, DAT-SPECT and receptor PET imaging have distinct clinical interpretations [54], and increased dopamine turnover early at symptom onset has implicated presynaptic mechanisms at this disease stage [55]. Furthermore, healthy aged D1 receptor expression is relatively uninformative as it is comparatively homogeneous across cortical regions (Supplementary Fig. S6) and likely redundant to the model in the presence of individualized DAT imaging. Network degeneration hypotheses of PD pathogenesis implicate various mechanisms from the propagation of neurotoxic alpha-synuclein [56] to the structural and functional neurodegeneration following striatal denervation [57]. We note that propagation is only a small contributor to the accumulation of pathology, and is dwarfed by local effects in our models (Fig. 3a). These findings may potentially reflect distinct disease phases. Our cohort was composed entirely of PD patients, for whom propagative, disease seeding processes may have already occurred, and neurodegeneration may now be driven by local effects. Furthermore, white matter tractography may not completely capture the connectivity between our cyto- and receptor-architectonically defined regions. A more complete treatment may consider vascular connectivity as well [29] [30], which may also be a substrate for pathology propagation. We find notable glutamatergic effects on multiple (micro-)structural factors (Supplementary Fig. S1): gray matter atrophy (NMDA), directed microstructural damage (AMPA and kainate), and dendrite density (kainate and NMDA). As NMDA and AMPA receptors are postsynaptic targets of glutamate, these mechanisms likely reflect the structural consequences of excitotoxicity and cell death [58]. On the other hand, kainate is believed to modulate synaptic transmission and plasticity [59], which may affect dendritic density. In our models, NMDA receptor influence is focused on occipital and temporal regions, AMPA influence is highest in frontal regions, and kainate influences mainly dendrite loss in both frontal and occipital regions. Among glutamatergic receptors, influence on microstructure of the motor cortex (MD, FA and t1/t2) is prominent, although it is more limited for atrophy or functional alterations. The notable stable roles of GABAergic receptors suggest their involvement via altered neuronal activity inhibition, interaction with the dopaminergic system, and potential regional vulnerability to microstructural degradation or dendrite loss. Inter-subject variability along the primary, mainly motor axis correlates with GABAergic mechanisms affecting microstructure (MD and t1/t2) and functional activity. Furthermore, a magnetic resonance spectroscopy (MRS) study found reduced levels of GABA in the visual cortex of PD patients [60], consistent with the regions of maximal influence of GABAA and GABAB receptors in our model. Various non-dopaminergic neurotransmitter systems have been associated with specific symptoms in PD, including cholinergic memory defects, adrenergic impairment of attention, and serotonin-driven depression [61] and visual hallucinations [62] [63]. Comparing model-derived receptor mechanisms and clinical assessments across PD patients, we observe two main axes of co-variability. The primary component represents motor/psychomotor symptoms associated prominently with GABAergic mechanisms, with secondary contributions from glutamatergic, cholinergic, and serotonergic systems (Supplementary Table S7). The secondary component is defined by visuospatial, memory and psychiatric symptoms, with the cholinergic system being the dominant receptor family. Mechanisms affecting neuronal activity dysfunction are more prominent in the primary component, while those affecting dendrite loss are greater in the secondary component. Notably, receptor mechanisms affecting microstructure contribute strongly to both axes. The secondary component is consistent with the cholinergically-driven visuospatial aspect of the dual-syndrome hypothesis of PD [18]. Stable cholinergic mechanisms are also present for every biological factor except directed microstructure, most notably the muscarinic M2 and nicotinic α4β2 contributions to gray matter atrophy, and muscarinic contributions to undirected microstructure. Specifically, we note prominent muscarinic M2 and nicotinic α4β2 receptor influences (on MD and GM, respectively) on the primary somatosensory cortex, a site of reduced activation in PD [64]. Our model suggests that nicotinic and muscarinic cholinergic systems strongly affect PD symptoms along specific pathways primarily involving dendritic density and degradation of microstructure (Fig. 6b). While typically associated with cognitive impairment and dementia in PD, cholinergic degeneration is also linked to depressive mood, apathy, olfaction, sleep disorder, and postural and gait disorder [65]. Epidemiological studies of smokers suggest a neuroprotective role for nicotinic receptors [66], which experience widespread decrease in PD [67]. The cholinergic and dopaminergic systems interact at biochemical, circuit and functional levels [61], tightly coupled by nicotinic receptors expressed on striatal dopaminergic neurons and acetylcholine [61] [68] modulate dopaminergic neurotransmission. An imbalance of cholinergic and dopaminergic neurotransmission may thus underlie PD cognitive dysfunction [61]. Our results suggest that cholinergic receptor distributions contribute to both motor and non-motor axes, albeit via distinct (i.e. functional/dopaminergic and microstructural) pathways (Fig. 6a,b). In addition to mediating inter-factor interactions, dysfunctional interactions between receptors may also be involved in neurodegeneration. Neurotransmitter release is regulated by presynaptic auto- and hetero-receptors [69], which in PD is potentially impaired in the dopaminergic system [70] and in GABAergic inhibition of the motor cortex [71]. Where possible, concurrent receptor or transporter imaging in a PD cohort would help clarify the role of neurotransmission balance in neurodegeneration. The lack of healthy control subjects currently limits our model’s ability to disentangle mechanisms of healthy ageing from PD-specific neurodegeneration. As longitudinal data collection progresses in large cohorts, model-derived mechanisms may help differentiate mechanisms distinct to PD and its (genetic or clinical) subtypes, Parkinson-plus syndromes, other neurodegenerative diseases and healthy ageing. Despite the prevalence of PD, the causes of this neurodegenerative condition remain unknown, and treatment is limited to symptomatic therapy complicated by individual variability in clinical presentation, side effects and treatment response [72]. Our work sheds light on the complex, especially non-dopaminergic neurotransmitter receptor-mediated mechanisms underlying brain reorganization and symptomatic variability in PD. Since neurotransmitter receptors are clinically efficacious drug targets [73], future work will explore the use of our personalized modeling approach to design personalized receptor-based therapy. ## Methods ### Ethics Statement Neuroimaging and clinical data in this study was acquired through the multi-center Parkinson’s Progression Markers Initiative (PPMI; [ppmi-info.org](http://ppmi-info.org)). Following good clinical practices, study subjects and/or authorized representatives gave written informed consent at the time of enrollment for sample collection and completed questionnaires approved by each participating site Institutional Review Board (IRB). The authors obtained approval from the PPMI for data use and publication, see documents [https://www.ppmi-info.org/documents/ppmi-data-use-agreement.pdf](https://www.ppmi-info.org/documents/ppmi-data-use-agreement.pdf) and [https://www.ppmi-info.org/documents/ppmi-publication-policy.pdf](https://www.ppmi-info.org/documents/ppmi-publication-policy.pdf), respectively. ### Data description and processing #### Study participants This study used longitudinal data from N=71 participants from the PPMI from 12 international sites. Demographic information is summarized in Supplementary Table S1. The inclusion criterion was the presence of at least 3 different imaging modalities (i.e. structural MRI, resting functional MRI, diffusion MRI and/or dopamine SPECT) over at least 3 visits at the time of our analysis. #### Structural MRI acquisition/processing Brain structural T1- and T2-weighted 3D images were acquired for all N=71 subjects. A detailed description of acquisition details can be found from the PPMI procedures manuals at [http://www.ppmi-info.org/](http://www.ppmi-info.org/). T1- and T2-weighted images from 3T scanners were acquired as a 3D sequence with a slice thickness of 1.5 mm or less, under three different views: axial, sagittal and coronal. All images underwent non-uniformity correction using the N3 algorithm [74]. Next, they were segmented into gray matter probabilistic maps using SPM12 (fil.ion.ucl.ac.uk/spm). Gray matter segmentations were standardized to MNI space [75] using the DARTEL tool [76]. Each map was modulated in order to preserve the total amount of signal/tissue. Mean gray matter density [76] values were calculated for the regions described in *Methods: Data description and processing: Receptor densities and brain parcellation*. #### Resting fMRI acquisition/processing Resting-state functional images were obtained using an echo-planar imaging sequence on 3T MRI scanners for N=71 subjects. For a detailed description of acquisition protocols, please see [http://www.ppmi-info.org](http://www.ppmi-info.org). Acquisition parameters were: 140 time points, repetition time (TR)=2400 ms, echo time (TE)=25 ms, flip angle=80°, number of slices=40, slice thickness=3.3 mm, in plane resolution=3.3 mm and in plane matrix=68×66. Pre-processing steps included: 1) motion correction, 2) slice timing correction, 3) alignment to the structural T1 image, and 4) spatial normalization to MNI space using the registration parameters obtained for the structural T1 image with the nearest acquisition date, and 5) signal filtering to keep only low frequency fluctuations (0.01–0.08 Hz) [77]. For each brain region, our model requires a local (i.e. intra-regional, non-network) measure of functional activity, in order to maintain mechanistic interpretability and to prevent data leakage of network information into local model terms (described further in *Receptor-Enriched Multifactorial Causal Model*). Due to its high correlation with glucose metabolism [78] and disease progression in PD [46], we calculated regional fractional amplitude of low-frequency fluctuation (fALFF) [79] as a measure of functional integrity. #### Diffusion MRI acquisition/processing Diffusion MRI (dMRI) images were acquired using standardized protocol on 3T MRI machines from 32 different international sites. Diffusion-weighted images were acquired along 64 uniformly distributed directions using a b-value of 1000 s/mm2 and a single b = 0 image. Single shot echo-planar imaging (EPI) sequence was used (116 × 116 matrix, 2 mm isotropic resolution, TR/TE 900/88 ms, and twofold acceleration). An anatomical T1-weighted 1 mm3 MPRAGE image was also acquired. Each patient underwent two baseline acquisitions and a further two one year later. More information on the dMRI acquisition and processing can be found online at [http://www.ppmi-info.org/](http://www.ppmi-info.org/). Preprocessing steps included: 1) motion and eddy current correction [80], 2) EPI distortion correction, 2) alignment of the T1-weighted image to the b0 image based on mutual information, 3) calculation of the deformation field between the diffusion and T1-weighted images, 4) calculation of the voxelwise diffusion tensors, 5) alignment to the structural T1 image, and 6) spatial normalization to MNI space [75] using the registration parameters obtained for the structural T1 image with the nearest acquisition date, and 6) calculation of mean values of summary metrics (FA and MD) for each considered brain region. #### Dopamine SPECT acquisition/processing A 111-185 MBq (3-5 mCi) bolus injection of I-123 FB-CIT was administered to each participant (N=71), and the SPECT scan was performed 4 hours post-injection. Raw projection data was acquired as a 128×128 matrix and the SPECT image was reconstructed. Attenuation correction and Gaussian blurring with a 3D 6mm filter were applied. The reconstructed and corrected SPECT images were normalized and registered to MNI space [75], and average values were calculated for all considered regions of interest. #### Receptor densities and brain parcellation In-vitro quantitative receptor autoradiography was applied to measure the densities of 15 receptors in 57 cytoarchitectonically defined cortical areas spread throughout the brain [81]. These receptors span major neurotransmitter systems, and show significant regional variability across the brain. Brains were obtained through the body donor programme of the University of Düsseldorf. Donors (three male and one female; between 67 and 77 years of age) had no history of neurological or psychiatric diseases, or long-term drug treatments. Causes of death were non-neurological in each case. Each hemisphere was sliced into 3 cm slabs, shock frozen at −40C, and stored at −80C. Receptors for the neurotransmitters glutamate (AMPA, NMDA, kainate), GABA (GABAA, GABAA-associated benzodiazepine binding sites, GABAB), acetylcholine (muscarinic M1, M2, M3, nicotinic α4β2), noradrenaline (α1, α2), serotonin (5-HT1A, 5-HT2), and dopamine (D1) were labeled according to previously published binding protocols consisting of pre-incubation, main incubation and rinsing steps [81]. The ligands used are summarized in Supplementary Table S3. Receptor densities were quantified by densitometric analysis of the ensuing autoradiographs, and areas were identified by cytoarchitectonic analysis in sections neigbouring those processed for receptor autoradiography, and which had been used for the visualization of cell bodies [82]. A brain parcellation was then defined with the aid of the Anatomy Toolbox [83] using 57 regions of interest for which receptor densities were available [31]. This parcellation was based on areas identified by cortical cytoarchitecture, as well as other cyto- and receptor-architectonically defined regions with receptor measurements (regions are summarized in Supplementary Table S4). These 57 regions were mirrored across left and right hemispheres for a total of 114 brain regions in our parcellation. For each receptor, regional densities were normalized using the mean and standard deviation across all brain regions. The structural T1 images of the Jülich [83], Brodmann [84], AAL3 [85] and DISTAL [86] brain parcellations were registered to the MNI ICBM152 T1 template using FSL 5.0’s FLIRT affine registration tool [87], and the obtained transformations were used to project the corresponding parcellations to the MNI ICBM152 space (using nearest neighbor interpolation to conserve original parcellation values). In the MNI ICBM152 space, voxels corresponding to the cytoarchitectonically-defined regions from [31] were identified from the regions in the Anatomy Toolbox, with the remaining Brodmann regions filled in using the Brodmann brain atlas. Supplementary Table S4 summarizes the ROI maps used to create the Brain atlas for regions with receptor data. The resulting parcellation of 114 brain regions in the common template space was then quality controlled, and small regions under 50 voxels were excluded. The resulting atlas with 155 bilateral brain regions (95 of which had receptor data) was used to extract whole-brain multi-modal neuroimaging data and estimate the diffusion-based connectivity matrix, as described in *Methods: Multimodal neuroimaging data and Methods: Anatomical connectivity estimation*. #### Anatomical connectivity estimation The connectivity matrix was constructed using DSI Studio ([http://dsi-studio.labsolver.org](http://dsi-studio.labsolver.org)). A group average template was constructed from a total of 1065 subjects [88]. A multi-shell diffusion scheme was used, and the b-values were 990, 1985 and 2980 s/mm2. The number of diffusion sampling directions were 90, 90, and 90, respectively. The in-plane resolution was 1.25 mm. The slice thickness was 1.25 mm. The diffusion data were reconstructed in the MNI space using q-space diffeomorphic reconstruction [89] to obtain the spin distribution function [90]. A diffusion sampling length ratio of 2.5 was used, and the output resolution was 1 mm. The restricted diffusion was quantified using restricted diffusion imaging [91]. A deterministic fiber tracking algorithm [92] was used. A seeding region was placed at whole brain. The QA threshold was 0.159581. The angular threshold was randomly selected from 15 degrees to 90 degrees. The step size was randomly selected from 0.5 voxel to 1.5 voxels. The fiber trajectories were smoothed by averaging the propagation direction with a percentage of the previous direction. The percentage was randomly selected from 0% to 95%. Tracks with length shorter than 30 or longer than 300 mm were discarded. A total of 100000 tracts were calculated. A custom brain atlas based on cytoarchitectonic regions with neurotransmitter receptor data [31] was used as the brain parcellation, as described in *Methods*: *Data description and processing: Receptor densities and brain parcellation*, and the connectivity matrix was calculated by using count of the connecting tracks. ### Multimodal neuroimaging data After pre-processing PPMI neuroimaging data for all 6 modalities, data harmonization was performed using ComBat [93] to correct for site and scanner effects. After extracting harmonized neuroimaging data for the cytoarchitectonically defined atlas described in *Methods: Data description and processing: Receptor densities and brain parcellation*, subjects lacking sufficient longitudinal or multimodal data were discarded. The disqualification criteria were i) fewer than 4 imaging modalities with data, or ii) fewer than 3 longitudinal samples for all modalities. For the remaining subjects, missing neuroimaging modalities (primarily FA, MD and t1/t2 ratios) at each visit were imputed using trimmed scores regression. Finally, a total of N=71 subjects were left with all 6 neuroimaging modalities with an average of 3.59 (± 0.50) time points. We used the mean and variance of each neuroimaging modality across all regions to calculate z-scores of neuroimaging data across all subjects. Please see Supplementary Table S1 for demographic characteristics. ### Clinical scores We used multiple composite scores derived from the PPMI clinical (motor, non-motor, psychiatric, cognitive, etc.) testing battery, namely the Benton Judgment of Line Orientation Test (BJLOT [94]), Geriatric Depression Scale (GDS [95]), Hopkins Verbal Learning Test (HVLT [96]), Letter Number Sequencing (LNS [97]), Movement Disorders Society – Unified Parkinson’s Disease Rating Scale (MDS-UPDRS [98]) Parts 1 (non-motor aspects of daily living; NP1), 2 (motor aspects of daily living; NP2), and 3 (motor exam; NP3), the Montreal Cognitive Assessment (MoCA [99]), semantic fluency (SF), State-Trait Anxiety Inventory for Adults (STAIAD [100]), and Symbol Digit Modalities (SDM [101]) tests. Protocols for deriving each score are described in the respective PPMI protocols documentation. We calculated cognitive decline as the linear best fit rate of change of each cognitive score with respect to examination date. Thus, for each patient, symptomatic decline was represented by a set of 11 rates of change. Average numbers of longitudinal evaluations per clinical score are summarized in Supplementary Table S2. ### Receptor-Enriched Multifactorial Causal Model (re-MCM) Multifactorial causal modeling is a generalized framework [29] [32] that treats the brain as a dynamical system of ROIs characterized by multiple interacting neuroimaging-quantified biological factors. Pathology may develop over time in each factor, affecting other factors locally and propagating to neigbouring regions via anatomical connections. We introduce the receptor-enriched multifactorial causal model (re-MCM), in which the local densities of various neurotransmitter receptors mediate interactions between biological factors at each brain region. In this work, the biological factors are gray matter density, neuronal activity, presynaptic dopamine, demyelination/dendritic density and two measures of white matter integrity, derived from structural T1 MRI, resting state functional MRI (rs-fMRI), DAT-SPECT, T1/T2 ratio, FA and MD, respectively. For any given subject and at a particular brain region *i*, the level of pathology of each biological factor *m* is represented by a single variable *Sm,i*, calculated as the deviation from the neuroimaging signal at the baseline visit. The temporal evolution of pathology *Sm,i* in modality *m* at brain region *i* is given by following differential equation: ![Formula][1] The functions f and *g* govern the global biological factor dynamics that are consistent across all brain regions. The local component *f*(***S****,*i*(*t*), ***R****,*i*) is the cumulative effect of all biological factors on factor *m* within region *i* mediated by ***R****,*i*, composed of local densities *rk,i* of a receptor *k* at a region *i*. The propagation term *g*(*sm*,*(*t*), *Ci*↔*) represents the net spreading of pathology in factor *m* along anatomical connections *Ci*↔* of the region *i*. Since the inter-visit interval of approximately 6 months is significantly shorter than the temporal scale of neurodegeneration, we assume a locally linear, time-invariant dynamical system: ![Formula][2] where *Ci*↔j is the directed anatomical connectivity from region *i* to *j*, and ![Graphic][3] the local rate of change of neuroimaging data for successive longitudinal samples at times *t′* and *t*: ![Formula][4] Local effects include i) direct factor-factor effects, ii) interaction terms mediated by *N*rec = 15 receptor types, and iii) direct receptor effects on the biological factor rate of change ![Graphic][5] (the second term in Equation 2). The first term in Equation 2 is thus expanded: ![Formula][6] The propagation term assumes symmetric connectivity *Cj*↔i between regions *i* and *j*, using a template connectivity matrix for all subjects, as described in *Anatomical connectivity estimation*, so we define the propagation component as: ![Formula][7] Thus, for each subject, the evolution of pathology in each biological factor m at region *i* is described by: ![Formula][8] Each model contains a set of *N*params = *N*fac × (1 + *N*rec) + *N*rec + 1 = 113 parameters ![Graphic][9] for subject *x* and factor *m* (or 678 total parameters per subject), each with a distinct neurobiological interpretation (e.g. the effect of reduced white matter integrity on gray matter atrophy mediated by glutamatergic receptor density). We perform linear regression, using the terms in Equation 6 as predictors with longitudinal PPMI neuroimaging samples *Sm,i*(*t*) and receptor maps **R**, to fit parameters ![Graphic][10] for each subject *x* and modality *m*. Separate regression models were built for i) each of the N=71 qualifying subjects, and ii) each of the 6 neuroimaging factors. These subjects were drawn from the PPMI dataset with at least 3 recorded neuroimaging modalities, and at least 3 longitudinal samples for at least one modality. We then calculate the coefficient of determination (*R*2) for each model to evaluate model fit, summarized in Figure 2. With the data vector **y** with elements ![Graphic][11], and model predictions ![Graphic][12] with ![Graphic][13], the coefficient of determination is ![Formula][14] where < **y***m* > is the mean of neuroimaging data for modality *m* across all brain regions and longitudinal samples. ### Statistical analysis #### Model fit For each subject and neuroimaging modality, we evaluated the quality of model fit by calculating the coefficient of determination (R2). Secondly, to evaluate the improvement in model fit due to receptor and receptor-mediated interaction terms while accounting for the difference in model size for each subject, we used F-tests (p<0.05) to compare the model fit of the full, receptor-neuroimaging interaction models (113 parameters per modality) with restricted, neuroimaging-only (8 parameters per modality) models. Finally, we evaluated the significance of the improvement in model fit (R2) due to actual receptor distributions with a permutation test using 1000 iterations of randomly permuted receptor maps (with receptor densities shuffled across regions independently for each receptor type), calculating the p-value of the true receptor data model R2 compared to the distribution from the permutations. #### Covariance of biological mechanisms with clinical symptoms To identify multivariate links between receptor-mediated biological mechanisms and to clinical symptoms in PD, we performed a data-driven cross-covariance analysis. Using singular value decomposition (SVD) to factorize the population covariance matrix between re-MCM parameters and clinical assessments (summarized in *Methods: Clinical scores*) to its eigenvectors, we identify multivariate axes of co-varying features. Different axes represent orthogonal disease processes affecting symptom severity. Permutation tests and bootstrapping ensure the statistical significance of the axes and the stability of identified mechanisms and symptoms, respectively. The algorithm is summarized as follows 1. We performed SVD on the cross-covariance matrix between all 678 re-MCM parameters and rates of clinical decline for N=71 PD patients, adjusted for covariates (baseline age, education and gender). SVD simultaneously reduces the dimensionality of features, and ranks them by their contribution to each axis. The cross-covariance matrix C = XY′ of the z-scores of re-MCM parameters *X* and the z-scores of the clinical decline rates *Y* is decomposed as ![Formula][15] where *U* and *V* are orthonormal matrices of spatial loadings for the parameters and clinical scores, respectively, and *S* is a diagonal matrix of singular values {*s*1,…,*s*7}. 2. We then performed permutation tests by shuffling the mapping between subjects’ re-MCM parameters and clinical scores, and repeating Step 1 for 1000 iterations, to evaluate the significance of SVD components. We performed a Procrustes transformation to align the axes of singular components in order to compare components from permuted iterations. We retained only those significant (p<0.05 with respect to the permuted distribution) singular components. 3. To discard non-stable re-MCM parameters and clinical assessments in each axis, we performed 1000 iterations of bootstrapping on the parameters *X* and clinical scores *Y*. To compare permuted iterations, we performed a Procrustes transformation to align the axes of singular components. We discarded the parameters with non-stable 95% confidence intervals. 4. For the remaining stable re-MCM parameters and clinical scores, and significant SVD components, we computed the variance explained per parameter *j* along each axis *i*: ![Formula][16] ### Regional influence To infer the spatial patterns of receptor involvement in neurodegeneration, we examined the improvement in neuroimaging models due to the inclusion of each receptor map. For each biological factor *m*, receptor *k* and brain region *i*, we fit a restricted, single-receptor version of the model ![Formula][17] where the longitudinal rate of change of each factor is predicted by its network propagation, direct factor effects, the local density of a single receptor *k*, and factor interactions with the density of only receptor *k*. We compare this model with a restricted, neuroimaging-only model excluding receptor density and interactions: ![Formula][18] To generate brain maps representing receptor influence on neuroimaging changes, 1. for each subject, we fit the single receptor and neuroimaging-only models for all biological factors and receptors, and studentize the residuals across regions and time points, 2. we combine the studentized residuals corresponding to each region across subjects and time points, and calculate the Wilcoxon rank sum statistic ![Graphic][19] between studentized residuals of the two models, 3. we compute a null distribution of the Wilcoxon statistic by repeating Steps 1-2 with 1000 randomly permuted receptor maps per imaging modality and receptor, 4. to estimate the significance of the Wilcoxon maps of each receptor across all 6 imaging modalities, we calculate the z-scores ![Graphic][20] of the Wilcoxon statistic ![Graphic][21] to its null distribution. ## Data and code availability The three datasets used in this study are available from the PPMI database (neuroimaging and cognitive evaluations; [https://www.ppmi-info.org/](https://www.ppmi-info.org/)), the HCP database (tractography template for connectivity estimation; [http://www.humanconnectomeproject.org/](http://www.humanconnectomeproject.org/)), and receptor autoradiography data published in [31]. We anticipate that the re-MCM method will be released soon as part of our available and open-access, user-friendly software [102] ([https://www.neuropm-lab.com/neuropm-box.html](https://www.neuropm-lab.com/neuropm-box.html)). ## Supporting information Supplementary [[supplements/287578_file03.docx]](pending:yes) ## Data Availability The three datasets used in this study are available from the PPMI database (neuroimaging and cognitive evaluations; [https://www.ppmi-info.org/](https://www.ppmi-info.org/)), the HCP database (tractography template for connectivity estimation; [http://www.humanconnectomeproject.org/](http://www.humanconnectomeproject.org/)), and receptor autoradiography data published in [Zilles and Palomero-Gallagher, 2017]. [https://www.ppmi-info.org/](https://www.ppmi-info.org/) [http://www.humanconnectomeproject.org/](http://www.humanconnectomeproject.org/) ## Funding This research was undertaken thanks in part to funding from: the Parkinson’s Canada graduate training award to AFK, the *Canada First Research Excellence Fund*, awarded to McGill University for the *Healthy Brains for Healthy Lives Initiative*, the Canada Research Chair tier-2, *Fonds de la recherche en santé du Québec* (FRQS) Junior 1 Scholarship, Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant, and Weston Brain Institute awards to YIM, the *Brain Canada Foundation* and *Health Canada* support to the McConnell Brain Imaging Center at the Montreal Neurological Institute, and the *European Union’s Horizon 2020 Framework Programme for Research and Innovation* under the Specific Grant Agreements 785907 (Human Brain Project SGA2) and 945539 (Human Brain Project SGA3) awarded to NPG and KZ. Multimodal imaging and clinical data collection and sharing for this project was funded by PPMI. A public-private partnership, PPMI is funded by the Michael J. Fox Foundation for Parkinson’s Research and funding partners, including AbbVie, Allergan, Amathus Therapeutics, Avid Radiopharmaceuticals, Biogen, BioLegend, Bristol Myers Squibb, Celgene, Denali Therapeutics, GE Healthcare, Genentech, GlaxoSmithKline plc., Golub Capital, Handl Therapeutics, Insitro, Janssen Neuroscience, Eli Lilly and Company, Lundbeck, Merck Sharp & Dohme Corp., Meso Scale Discovery, Neurocrine Biosciences, Pfizer Inc., Piramal Group, Prevail Therapeutics, Roche, Sanofi Genzyme, Servier Laboratories, Takeda Pharmaceutical Company Limited, Teva Pharmaceutical Industries Ltd., UCB, Verily Life Sciences, and Voyager Therapeutics Inc. ## Competing Interests The authors report no competing interests. ## Supplementary material Supplementary material is available online. ## Acknowledgements The authors would like to acknowledge the integral role of the late Prof. Dr. Karl Zilles in collecting the receptor autoradiography data. ## Abbreviations PD : Parkinson’s disease re-MCM : receptor-enriched multifactorial causal model ROI : region(s) of interest SVD : singular value decomposition PPMI : Parkinson’s Progression Marker Initiative. * Received March 22, 2023. * Revision received March 22, 2023. * Accepted March 22, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1]. G. Alexander, “Biology of Parkinson’s disease: pathogenesis and pathophysiology of a multisystem neurodegenerative disorder,” Dialogues in clinical neuroscience, vol. 6, no. 3, pp. 259–280, 2004. 2. [2]. J. Han, Y. Ahn, W. Kim, C. Shin, S. Jeong, Y. Song, Y. Bae and J. Kim, “Psychiatric manifestation in patients with Parkinson’s disease,” Journal of Korean medical science, vol. 33, no. 47, p. e300, 2018. 3. [3]. K. Jellinger, “Neuropathology of sporadic Parkinson’s disease: evaluation and changes of concepts,” Movement disorders, vol. 27, no. 1, 2012. 4. [4]. N. Titova, S. Lewis, C. Padmakumar and K. Chaudhuri, “Parkinson’s: a syndrome rather than a disease?,” Journal of Neural Transmission, vol. 124, no. 8, 2017. 5. [5]. A. Sauerbier, M. Qamar, T. Rajah and K. Chaudhuri, “New concepts in the pathogenesis and presentation of Parkinson’s disease,” Clinical Medicine, vol. 16, no. 4, p. 365, 2016. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImNsaW5tZWRpY2luZSI7czo1OiJyZXNpZCI7czo4OiIxNi80LzM2NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 6. [6]. R. von Coelln and L. Shulman, “Clinical subtypes and genetic heterogeneity: of lumping and splitting in Parkinson disease,” Current opinion in neurology, vol. 29, no. 6, pp. 727–734, 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/WCO.0000000000000384&link_type=DOI) 7. [7]. O. W. and J. Schulz, “Current and experimental treatments of Parkinson disease: A guide for neuroscientists,” Journal of neurochemistry, vol. 139, pp. 325–337, 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jnc.13750&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27577098&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 8. [8]. L. Brichta, P. Greengard and M. Flajolet, “Advances in the pharmacological treatment of Parkinson’s disease: targeting neurotransmitter systems,” Trends in neurosciences, vol. 36, no. 9, pp. 543–54, 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tins.2013.06.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23876424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325044000006&link_type=ISI) 9. [9]. A. Espay, F. Morgante, A. Merola, A. Fasano, L. Marsili, S. Fox, E. Bezard, B. Picconi, P. Calabresi and A. Lang, “Levodopa-induced dyskinesia in Parkinson disease: current and evolving concepts,” Annals of Neurology, vol. 84, no. 6, pp. 797–811, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.25364&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30357892&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 10. [10]. C. Adler, T. Beach, J. Hentz, H. Shill, J. Caviness, E. Driver-Dunckley, M. Sabbagh, L. Sue, S. Jacobson and C. D. B. Belden, “Low clinical diagnostic accuracy of early vs. advance Parkinson disease (Clinicopathologic study),” Neurology, vol. 83, p. 406–412, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000000641&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24975862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 11. [11]. G. Pagano, F. Niccolini and M. Politis, “Imaging in Parkinson’s disease,” Clinical Medicine, vol. 16, no. 4, p. 371, 2016. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImNsaW5tZWRpY2luZSI7czo1OiJyZXNpZCI7czo4OiIxNi80LzM3MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 12. [12]. J. Ellis and M. Fell, “Current approaches to the treatment of Parkinson’s Disease,” Bioorganic & Medicinal Chemistry Letters, vol. 27, no. 18, pp. 4247–55, 2017. 13. [13]. P. A. LeWitt and K. Chaudhuri, “Unmet needs in Parkinson disease: Motor and non-motor,” Parkinsonism & Related Disorders, vol. 80, pp. S7–S12, 2020. 14. [14]. Y. Xu, J. Yan, P. Zhou, J. Li, H. Gao, Y. Xia and Q. Wang, “Neurotransmitter receptors and cognitive dysfunction in alzheimer’s disease and parkinson’s disease,” Progress in neurobiology, vol. 97, no. 1, p. 1–13, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pneurobio.2012.02.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22387368&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 15. [15]. S. Shang, H. Zhang, Y. Feng, J. Wu, W. Dou, Y. Chen and X. Yin, “Region-Specific neurovascular decoupling associated with cognitive decline in parkinson’s disease,” Frontiers in aging neuroscience, vol. 13, p. 770528, 2021. 16. [16]. A. Zarkali, P. McColgan, L. Leyland, A. Lees, G. Rees and R. Weil, “Organisational and neuromodulatory underpinnings of structural-functional connectivity decoupling in patients with Parkinson’s disease,” Communications biology, vol. 4, no. 1, pp. 1–13, 2021. 17. [17]. G. M. Halliday, J. B. Leverenz, J. S. Schneider and C. H. Adler, “The neurobiological basis of cognitive impairment in Parkinson’s disease,” Movement Disorders, vol. 29, no. 5, pp. 634–650, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.25857&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24757112&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 18. [18]. A. A. Kehagia, R. A. Barker and T. W. Robbins, “Cognitive impairment in Parkinson’s disease: the dual syndrome hypothesis,” Neurodegenerative diseases, vol. 11, no. 2, pp. 79–92, 2013. 19. [19]. N. Bohnen, A. Yarnall, R. Weil, E. Moro, M. Moehle, P. Borghammer, M. Bedard and R. Albin, “Cholinergic system changes in Parkinson’s disease: emerging therapeutic approaches,” The Lancet Neurology, vol. 21, no. 4, pp. 381–392, 2022. 20. [20]. M. Politis and F. Niccolini, “Serotonin in Parkinson’s Disease,” Behavioural Brain Research, vol. 277, pp. 136–145, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbr.2014.07.037&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25086269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 21. [21]. Y. Grimbergen, J. Langston, R. Roos and B. Bloem, “Postural instability in Parkinson’s disease: the adrenergic hypothesis and the locus coeruleus,” Expert review of neurotherapeutics, vol. 9, no. 2, pp. 279–290, 2009. 22. [22]. C. Weingarten, M. Sundman, P. Hickey and N. Chen, “Neuroimaging of Parkinson’s disease: Expanding views,” Neuroscience & Biobehavioral Reviews, vol. 59, pp. 16–52, 2015. 23. [23]. N. Bidesi, I. Vang Andersen, A. Windhorst, V. Shalgunov and M. Herth, “The role of neuroimaging in Parkinson’s disease,” Journal of Neurochemistry, vol. 159, no. 4, pp. 660–89, 2021. 24. [24]. C. Atkinson-Clement, S. Pinto, E. A. and O. Coulon, “Diffusion tensor imaging in Parkinson’s disease: review and meta-analysis,” NeuroImage: Clinical, vol. 16, pp. 98–110, 2017. 25. [25]. S. Lin, R. Rodriguez-Rojas, T. Baumeister, C. Lenglos, J. Pineda-Pardo, J. Máñez-Miró, M. Del Alamo, R. Martinez-Fernandez, J. Obeso and Y. Iturria-Medina, “Neuroimaging signatures predicting motor improvement to focused ultrasound subthalamotomy in Parkinson’s disease,” npj Parkinson’s Disease, vol. 8, no. 1, p. 70, 2022. 26. [26]. C. R. Jack Jr., D. S. Knopman, W. J. Jagust, R. C. Petersen, M. W. Weiner, P. S. Aisen, L. M. Shaw, P. Vemuri, H. J. Wiste, S. D. Weigand, T. G. Lesnick, V. S. Pankratz, M. C. Donohue and J. Q. Trojanowski, “Tracking pathophysiological processes in Alzheimer’s disease: an updated hypothetical model of dynamic biomarkers,” The Lancet Neurology, vol. 12, no. 2, pp. 207–216, 2013. 27. [27]. Y. Iturria-Medina, R. C. Sotero, P. J. Toussaint, J. M. Mateos-Pérez, A. C. Evans and t. A. D. N. Initiative, “Early role of vascular dysregulation on late-onset Alzheimer’s disease based on multifactorial data-driven analysis,” Nature Communications, vol. 7, no. 1, pp. 1–14, 2016. 28. [28]. C. Lenglos, S. Lin, Y. Zeighami, T. Baumeister, F. Carbonell and Y. Iturria-Medina, “Multivariate genomic and transcriptomic determinants of imaging-derived personalized therapeutic needs in Parkinson’s disease,” Scientific Reports, vol. 12, no. 1, p. 5483, 2022. 29. [29]. Y. Iturria-Medina, F. M. Carbonell, R. C. Sotero, F. Chouinard-Decorte and A. C. Evans, “Multifactorial causal model of brain (dis)organization and therapeutic intervention: Application to Alzheimer’s disease,” NeuroImage, vol. 152, pp. 60–77, 2017. 30. [30]. Y. Iturria-Medina, F. M. Carbonell and A. C. Evans, “Multimodal imaging-based therapeutic fingerprints for optimizing personalized interventions: Application to neurodegeneration,” NeuroImage, vol. 179, pp. 40–50, 2018. 31. [31]. K. Zilles and N. and Palomero-Gallagher, “Multiple transmitter receptors in regions and layers of the human cerebral cortex,” Frontiers in neuroanatomy, vol. 11, p. 78, 2017. 32. [32]. A. Khan, Q. Adewale, T. Baumeister, F. Carbonell, K. Zilles, N. Palomero-Gallagher and Y. Iturria-Medina, “Personalized brain models identify neurotransmitter receptor changes in alzheimer’s disease,” Brain, 2021. 33. [33]. M. Mühlau, “T1/T2-weighted ratio is a surrogate marker of demyelination in multiple sclerosis: No,” Multiple Sclerosis Journal, vol. 28, no. 3, pp. 355–356, 2022. 34. [34]. R. Righart, V. Biberacher, L. Jonkman, R. Klaver, P. Schmidt, D. Buck, A. Berthele, J. Kirschke, C. Zimmer, B. Hemmer, J. Geurts and M. Mühlau, “Cortical pathology in multiple sclerosis detected by the T 1/T 2-weighted ratio from routine magnetic resonance imaging,” Annals of neurology, vol. 82, no. 4, pp. 519–529., 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.25020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28833433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 35. [35]. K. Kiely, P. Butterworth, N. Watson and M. Wooden, “The Symbol Digit Modalities Test: Normative data from a large nationally representative sample of Australian,” Archives of Clinical Neuropsychology, vol. 29, no. 8, pp. 767–775, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/arclin/acu055&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25352087&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 36. [36]. G. Pagano, F. Niccolini and M. Politis, “Imaging in Parkinson’s disease,” Clinical Medicine, vol. 16, no. 4, p. 371, 2016. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImNsaW5tZWRpY2luZSI7czo1OiJyZXNpZCI7czo4OiIxNi80LzM3MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 37. [37]. G. Ballentine, S. Friedman and D. Bzdok, “Trips and neurotransmitters: Discovering principled patterns across 6850 hallucinogenic experiences,” Science Advances, vol. 8, no. 11, 16 March 2022. 38. [38]. Q. Adewale, A. Khan, F. Carbonell and Y. Iturria-Medina, “Integrated transcriptomic and neuroimaging brain model decodes biological mechanisms in aging and Alzheimer’s disease,” Elife, vol. 10, p. e62589, 2021. 39. [39]. D. Zachlod, S. Bludau, S. Cichon, N. Palomero-Gallagher and K. Amunts, “Combined analysis of cytoarchitectonic, molecular and transcriptomic patterns reveal differences in brain organization across human functional brain systems,” NeuroImage, p. 119286, 2022. 40. [40]. J. Hansen, R. Markello, L. Tuominen, M. Nørgaard, E. Kuzmin, N. Palomero-Gallagher, A. Dagher and B. Misic, “Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain,” bioRxiv, 2021. 41. [41]. K. Pak, S. Seo, M. Lee, H. Im, K. Kim and I. Kim, “Limited power of dopamine transporter mRNA mapping for predicting dopamine transporter availability,” Synapse, vol. 76, no. 5-6, p. e22226, 2022. 42. [42]. J. Dukart, Š. Holiga, C. Chatham, P. Hawkins, A. Forsyth, R. McMillan, J. Myers, A. Lingford-Hughes, D. Nutt, E. Merlo-Pich and C. Risterucci, “Cerebral blood flow predicts differential neurotransmitter activity,” Scientific reports, vol. 8, no. 1, pp. 1–11, 2018. 43. [43]. W. Moses, “Fundamental limits of spatial resolution in PET.,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 648, pp. S236–S240, 2011. 44. [44]. E. Mak, L. Su, G. Williams and J. O’Brien, “Neuroimaging correlates of cognitive impairment and dementia in Parkinson’s disease,” Parkinsonism & related disorders, vol. 21, no. 8, pp. 862–870, 2015. 45. [45]. W. Zhan, G. Kang, G. Glass, Y. Zhang, C. Shirley, R. Millin, K. Possin, M. Nezamzadeh, M. Weiner, W. Marks Jr and N. Schuff, “Regional alterations of brain microstructure in Parkinson’s disease using diffusion tensor imaging,” Movement disorders, vol. 27, no. 1, pp. 90–97, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.23917&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21850668&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 46. [46]. Y. Tang, L. Meng, C. Wan, Z. Liu, W. Liao, X. Yan, X. Wang, B. Tang and J. Guo, “Identifying the presence of Parkinson’s disease using low-frequency fluctuations in BOLD signals,” Neuroscience letters, vol. 645, pp. 1–6, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neulet.2017.02.056&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28249785&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 47. [47]. K. Nguyen, V. Raval, A. Treacher, C. Mellema, F. Yu, M. Pinho, R. Subramaniam, D. J. R.B. and A. Montillo, “Predicting Parkinson’s disease trajectory using clinical and neuroimaging baseline measures,” Parkinsonism & related disorders, vol. 85, pp. 44–51, 2021. 48. [48]. Y. Assaf, “Imaging laminar structures in the gray matter with diffusion MRI,” Neuroimage, vol. 197, pp. 677–688, 2019. 49. [49]. P. Weston, I. Simpson, N. Ryan, S. Ourselin and N. Fox, “Diffusion imaging changes in grey matter in Alzheimer’s disease: a potential marker of early neurodegeneration,” Alzheimer’s research & therapy, vol. 7, no. 1, pp. 1–8, 2015. 50. [50]. M. Cercignani, M. Inglese, E. Pagani, G. Comi and M. Filippi, “Mean diffusivity and fractional anisotropy histograms of patients with multiple sclerosis,” American Journal of Neuroradiology, vol. 22, no. 5, pp. 952–958, 2001. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo4OiIyMi81Lzk1MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 51. [51]. M. Bergamino, E. Keeling, V. Mishra, A. Stokes and R. Walsh, “Assessing white matter pathology in early-stage Parkinson disease using diffusion MRI: a systematic review,” Frontiers in neurology, vol. 11, p. 314, 2020. 52. [52]. H. Fu, J. Hardy and K. Duff, “Selective vulnerability in neurodegenerative diseases,” Nature Neuroscience, vol. 21, no. 10, pp. 1350–1358, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-018-0221-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30250262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 53. [53]. N. Volkow, G. Wang, J. Fowler, J. Logan, R. Hitzemann, Y. Ding, N. Pappas, C. Shea and K. Piscani, “Decreases in dopamine receptors but not in dopamine transporters in alcoholics,” Alcoholism: Clinical and Experimental Research, vol. 20, no. 9, pp. 1594–1598, 1996. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1530-0277.1996.tb05936.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8986209&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VZ86300016&link_type=ISI) 54. [54]. S. Thobois, S. Prange, C. Scheiber and E. Broussolle, “What a neurologist should know about PET and SPECT functional imaging for parkinsonism: A practical perspective,” Parkinsonism & related disorders, vol. 59, pp. 93–100, 2019. 55. [55]. M. Bu, M. Farrer and H. Khoshbouei, “Dynamic control of the dopamine transporter in neurotransmission and homeostasis,” NPJ Parkinson’s disease, vol. 7, no. 1, pp. 1–11, 2021. 56. [56]. Y. Zeighami, M. Ulla, Y. Iturria-Medina, M. Dadar, Y. Zhang, K. Larcher, V. Fonov, A. Evans, D. Collins and A. Dagher, “Network structure of brain atrophy in de novo Parkinson’s disease,” eLife, vol. 4, p. e08440, 7 Sep 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.08440&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26344547&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 57. [57]. M. Ruppert, A. Greuel, M. Tahmasian, F. Schwartz, S. Stürmer, F. Maier, J. Hammes, M. Tittgemeyer, L. Timmermann, T. Van Eimeren and A. Drzezga, “Network degeneration in Parkinson’s disease: multimodal imaging of nigro-striato-cortical dysfunction,” Brain, vol. 143, no. 3, pp. 944–959, 1 Mar 2020. 58. [58]. W. Caudle and J. Zhang, “Glutamate, excitotoxicity, and programmed cell death in Parkinson disease,” Experimental neurology, vol. 220, no. 2, pp. 230–233, 2009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.expneurol.2009.09.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19815009&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272428700005&link_type=ISI) 59. [59]. Z. Zhang, S. Zhang, P. Fu, Z. Zhang, K. Lin, J. Ko and K. Yung, “Roles of glutamate receptors in Parkinson’s disease,” International journal of molecular sciences, vol. 20, no. 18, p. 4391, 2019. 60. [60]. M. Firbank, J. Parikh, N. Murphy, A. Killen, C. Allan, D. Collerton, A. Blamire and J. Taylor, “Reduced occipital GABA in Parkinson disease with visual hallucinations,” Neurology, vol. 91, no. 7, pp. e675–e685, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000006007&link_type=DOI) 61. [61]. P. Calabresi, B. Picconi, L. Parnetti and M. Di Filippo, “A convergent model for cognitive dysfunctions in Parkinson’s disease: the critical dopamine–acetylcholine synaptic balance,” The Lancet Neurology, vol. 5, no. 11, pp. 974–983, 2006. 62. [62]. S. Vegas-Suarez, E. Paredes-Rodriguez, A. Aristieta, J. Lafuente, C. Miguelez and L. Ugedo, “Dysfunction of serotonergic neurons in Parkinson’s disease and dyskinesia,” International Review of Neurobiology, vol. 146, pp. 259–279, 2019. 63. [63]. B. Ballanger, A. Strafella, T. van Eimeren, M. Zurowski, P. Rusjan, S. Houle and S. Fox, “Serotonin 2A receptors and visual hallucinations in Parkinson disease,” Archives of neurology, vol. 67, no. 4, pp. 416–421, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archneurol.2010.35&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20385906&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276559500007&link_type=ISI) 64. [64]. A. Nelson, T. Hoque, C. Gunraj and R. Chen, “Altered somatosensory processing in Parkinson’s disease and modulation by dopaminergic medications,” Parkinsonism & Related Disorders, vol. 53, pp. 76–81, 2018. 65. [65]. M. Müller and N. Bohnen, “Cholinergic dysfunction in Parkinson’s disease,” Current neurology and neuroscience reports, vol. 13, no. 9, pp. 1–9, 2013. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 66. [66]. M. Quik and J. Kulak, “Nicotine and nicotinic receptors; relevance to Parkinson’s disease,” Neurotoxicology, vol. 23, no. 4-5, pp. 581–594, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0161-813X(02)00036-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12428730&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178985400017&link_type=ISI) 67. [67]. M. Fujita, M. Ichise, S. Zoghbi, J. Liow, S. Ghose, D. Vines, J. Sangare, J. Lu, V. Cropley, H. Iida and K. Kim, “Widespread decrease of nicotinic acetylcholine receptors in Parkinson’s disease,” Annals of neurology, vol. 59, no. 1, pp. 174–177, 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.20688&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16374823&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234447100020&link_type=ISI) 68. [68]. D. Lester, T. Rogers and C. Blaha, “Acetylcholine–dopamine interactions in the pathophysiology and treatment of CNS disorders,” CNS neuroscience & therapeutics, vol. 16, no. 3, pp. 137–162, 2010. 69. [69]. E. Schlicker and T. Feuerstein, “Human presynaptic receptors,” Pharmacology & Therapeutics, vol. 172, pp. 1–21, 2017. 70. [70]. R. de la Fuente-Fernández, M. Schulzer, E. Mak, D. Calne and A. Stoessl, “Presynaptic mechanisms of motor fluctuations in Parkinson’s disease: a probabilistic model,” Brain, vol. 127, no. 4, pp. 888–899, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awh102&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14960500&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220485500022&link_type=ISI) 71. [71]. J. Chu, A. Wagle-Shukla, C. Gunraj, A. Lang and R. Chen, “Impaired presynaptic inhibition in the motor cortex in Parkinson disease,” Neurology, vol. 72, no. 9, pp. 842–849, 2009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/01.wnl.0000343881.27524.e8&link_type=DOI) 72. [72]. K. Jamebozorgi, E. Taghizadeh, D. Rostami, H. Pormasoumi, G. Barreto, S. Hayat and A. Sahebkar, “Cellular and molecular aspects of Parkinson treatment: future therapeutic perspectives,” Molecular Neurobiology, pp. 4799–4811, 2019. 73. [73]. L. Brichta, P. Greengard and M. Flajolet, “Advances in the pharmacological treatment of Parkinson’s disease: targeting neurotransmitter systems,” Trends in neurosciences, vol. 36, no. 9, pp. 543–554, 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tins.2013.06.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23876424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325044000006&link_type=ISI) 74. [74]. J. G. Sled, A. P. Zijdenbos and A. C. Evans, “A nonparametric method for automatic correction of intensity nonuniformity in MRI data,” IEEE transactions on medical imaging, vol. 17, no. 1, p. 87–97, 1998. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/42.668698&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9617910&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000073646700008&link_type=ISI) 75. [75]. A. C. Evans, M. Kamber, D. Collins and D. MacDonald, “An MRI-based probabilistic atlas of neuroanatomy,” Magnetic resonance scanning and epilepsy, pp. 263–274, 1994. 76. [76]. J. Ashburner, “A fast diffeomorphic image registration algorithm,” NeuroImage, vol. 38, no. 1, pp. 95–113, 2007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2007.07.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17761438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249992900011&link_type=ISI) 77. [77]. C. Yan and Y. Zang, “DPARSF: A matlab toolbox for “pipeline” data analysis of resting-state fMRI,” Frontiers in systems neuroscience, vol. 4, p. 13, 2010. 78. [78]. S. E. C. A. P. S. C. C. Aiello M, A. Prinster, E. Nicolai, M. Salvatore, J. Baron and M. Quarantelli, “Relationship between simultaneously acquired resting-state regional cerebral glucose metabolism and functional MRI: a PET/MR hybrid scanner study,” Neuroimage, no. 113, pp. 111–121, 1 June 2015. 79. [79]. Q.-H. Zou, C.-Z. Zhu, Y. Yang, X.-N. Zuo, X.-Y. Long, Q.-J. Cao, Y.-F. Wang and Y.-F. Zang, “An improved approach to detection of amplitude of low-frequency fluctuation (ALFF) for resting-state fMRI: Fractional ALFF,” Journal of neuroscience methods, vol. 172, no. 1, pp. 137–141, 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2008.04.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18501969&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257485800020&link_type=ISI) 80. [80]. G. Rohde, A. Barnett, P. Basser, S. Marenco and C. Pierpaoli, “Rohde, G.K., Barnett, A.S., Basser, P.J., Marenco, S. and Pierpaoli, C., 2004. Comprehensive approach for correction of motion and distortion in diffusion-weighted MRI. Magnetic Resonance in Medicine:,” Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine, vol. 51, no. 1, pp. 103–114, 2004. 81. [81]. N. Palomero-Gallagher and K. Zilles, “Cyto-and receptor architectonic mapping of the human brain,” Handbook of clinical neurology, vol. 150, pp. 355–387, 2018. 82. [82]. B. Merker, “Silver staining of cell bodies by means of physical development,” Journal of neuroscience methods, vol. 9, no. 3, pp. 235–241, 1983. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0165-0270(83)90086-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6198563&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1983RX24700005&link_type=ISI) 83. [83]. S. Eickhoff, T. Paus, S. Caspers, M. Grosbras, A. Evans, K. Zilles and K. Amunts, “Assignment of functional activations to probabilistic cytoarchitectonic areas revisited,” NeuroImage, vol. 36, no. 3, pp. 511–521, 2007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2007.03.060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17499520&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247634400003&link_type=ISI) 84. [84]. K. Brodmann, Vergleichende Lokalisationslehre der Grosshirnrinde in ihren Prinzipien dargestellt auf Grund des Zellenbaues, Leipzig: Barth JA, 1909. 85. [85]. E. Rolls, C. Huang, C. Lin, J. Feng and M. Joliot, “Automated anatomical labelling atlas 3,” Neuroimage, vol. 206, p. 116189, 1 February 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2019.116189&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 86. [86]. S. Ewert, P. Plettig, N. Li, M. Chakravarty, D. Collins, T. Herrington, A. Kühn and A. Horn, “Toward defining deep brain stimulation targets in MNI space: a subcortical atlas based on multimodal MRI, histology and structural connectivity,” Neuroimage, vol. 170, pp. 271–282, 15 April 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.05.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28536045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 87. [87]. M. Jenkinson, C. Beckmann, T. Behrens, M. Woolrich and S. Smith, “FSL,” vol. 62, pp. 782–90, 2012. 88. [88]. F. Yeh, S. Panesar, D. Fernandes, A. Meola, M. Yoshino, J. Fernandez-Miranda, J. Vettel and T. Verstynen, “Population-averaged atlas of the macroscale human structural connectome and its network topology,” NeuroImage, vol. 178, pp. 57–68, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2018.05.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29758339&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 89. [89]. F. Yeh and W. I. Tseng, “NTU-90: A high angular resolution brain atlas constructed by q-space diffeomorphic reconstruction,” NeuroImage, vol. 58, no. 1, pp. 91–99, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2011.06.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21704171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000293548500012&link_type=ISI) 90. [90]. F. Yeh, V. J. Wedeen and W.-Y. I. Tseng, “Generalized q-sampling imaging,” IEEE transactions on medical imaging, vol. 29, no. 9, pp. 1626–1635, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TMI.2010.2045126&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20304721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283222600003&link_type=ISI) 91. [91]. F. Yeh, L. Liu, T. K. Hitchens and Y. L. Wu, “Mapping immune cell infiltration using restricted diffusion MRI,” Magnetic resonance in medicine, vol. 77, no. 2, pp. 603–612, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.26143&link_type=DOI) 92. [92]. F. Yeh, T. D. Verstynen, Y. Wang, J. C. Fernández-Miranda and W. I. Tseng, “Deterministic diffusion fiber tracking improved by quantitative anisotropy,” PloS One, pp. vol. 8, no. 11, 2013. 93. [93]. J. P. D. T. B. W. T. E. M. R. K. R. D. Fortin, T. Satterthwaite, R. Gur, R. Gur and R. Schultz, “Harmonization of multi-site diffusion tensor imaging data,” Neuroimage, vol. 161, pp. 149–170, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.08.047&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 94. [94]. C. Qualls, N. Bliwise and A. Stringer, “Short forms of the Benton judgment of line orientation test: Development and psychometric properties,” Archives of Clinical Neuropsychology, vol. 15, no. 2, pp. 159–163, 2000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/arclin/15.2.159&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14590559&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000084900600006&link_type=ISI) 95. [95]. J. Yesavage, “Geriatric depression scale,” Psychopharmacol Bull, vol. 24, no. 4, pp. 709–711, 1988. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3249773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1988T349900028&link_type=ISI) 96. [96]. J. Brandt, “The Hopkins Verbal Learning Test: Development of a new memory test with six equivalent forms,” The clinical neuropsychologist, vol. 5, no. 2, pp. 125–142, 1991. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/13854049108403297&link_type=DOI) 97. [97]. D. Wechsler, “Wechsler adult intelligence scale,” Archives of Clinical Neuropsychology, 1955. 98. [98]. C. Goetz, B. Tilley, S. Shaftman, G. Stebbins, S. Fahn, P. Martinez-Martin, W. Poewe, C. Sampaio, M. Stern, R. Dodel and B. Dubois, “Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results,” Movement disorders: official journal of the Movement Disorder Society, vol. 23, no. 15, pp. 2129–2170, 2008. 99. [99]. Z. Nasreddine, N. Phillips, V. Bédirian, S. Charbonneau, V. Whitehead, I. Collin, J. Cummings and H. Chertkow, “The Montreal Cognitive Assessment, MoCA: a brief screening tool for mild cognitive impairment,” Journal of the American Geriatrics Society, vol. 53, no. 4, pp. 695–699, 2005. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1532-5415.2005.53221.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15817019&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227899200021&link_type=ISI) 100.[100]. C. Spielberger, “State-trait anxiety inventory for adults,” 1983. 101.[101]. A. Smith, “Symbol digit modalities test,” Western psychological services, pp. 1–22, 1973. 102.[102]. Y. Iturria-Medina, F. Carbonell, A. Assadi, Q. Adewale, A. Khan, R. Baumeister and L. Sanchez-Rodriguez, “NeuroPM toolbox: integrating Molecular, Neuroimaging and Clinical data for Characterizing Neuropathological Progression and Individual Therapeutic Needs,” medRvix DOI: 10.1101/2020.09.24.20200964, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wOS4yNC4yMDIwMDk2NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjIvMjAyMy4wMy4yMi4yMzI4NzU3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 103.[103]. K. Roy, Computational modeling of drugs against Alzheimer’s disease, Springer, 2018. 104.[104]. A. Bejanin, D. Schonhaut, R. La Joie, J. Kramer, S. Baker, N. Sosa, N. Ayakta, A. Cantwell, M. Janabi, M. Lauriola and J. O’Neil, “Tau pathology and neurodegeneration contribute to cognitive impairment in Alzheimer’s disease,” Brain, vol. 140, no. 12, pp. 3286–3300, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awx243&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29053874&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 105.[105]. N. A. Bishop, T. Lu and B. A. Yankner, “Neural mechanisms of ageing and cognitive decline,” Nature, vol. 464, no. 7288, p. 529–535, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature08983&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20336135&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275974200036&link_type=ISI) 106.[106]. M. Busche, S. Wegmann, S. Dujardin, C. Commins, J. Schiantarelli, N. Klickstein, T. Kamath, G. Carlson, I. Nelken and B. Hyman, “Tau impairs neural circuits, dominating amyloid-β effects, in Alzheimer models in vivo,” Nature neuroscience, vol. 22, no. 1, pp. 57–64, 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-018-0289-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30559471&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 107.[107]. D. Butterfield and C. Pocernich, “The glutamatergic system and Alzheimer’s disease,” CNS drugs, vol. 17, no. 9, pp. 641–652, 2003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2165/00023210-200317090-00004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12828500&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184376000004&link_type=ISI) 108.[108]. J.-K. Choi, Y. I. Chen, E. Hamel and B. G. Jenkins, “Brain hemodynamic changes mediated by dopamine receptors: Role of the cerebral microvasculature in dopamine-mediated neurovascular coupling,” Neuroimage, vol. 30, no. 3, pp. 700–712, 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2005.10.029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16459104&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236894800003&link_type=ISI) 109.[109]. P. T. Francis, M. J. Ramírez and M. K. Lai, “Neurochemical basis for symptomatic treatment of alzheimer’s disease,” Neuropharmacology, vol. 59, no. 4-5, p. 221–229, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuropharm.2010.02.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20156462&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 110.[110]. M. Grothe, L. Zaborszky, M. Atienza, E. Gil-Neciga, R. Rodriguez-Romero, S. Teipel, K. Amunts, A. Suarez-Gonzalez and J. Cantero, “Reduction of basal forebrain cholinergic system parallels cognitive impairment in patients at high risk of developing Alzheimer’s disease,” Cerebral Cortex, vol. 20, no. 7, pp. 1685–1695, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhp232&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19889714&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278690800016&link_type=ISI) 111.[111]. S. Guntupalli, J. Widagdo and V. Anggono, “Amyloid-β-induced dysregulation of AMPA receptor trafficking,” Neural plasticity, vol. 2016, 2016. 112.[112]. H. Hampel, M.-M. M. A. C. Cuello, A. S. Khachaturian, A. Vergallo, M. R. Farlow, P. J. Snyder, E. Giacobini and Z. S. Khachaturian, “Revisiting the cholinergic hypothesis in alzheimer’s disease: Emerging evidence from translational and clinical research,” Alzheimer’s Dementia, vol. 6, no. 1, pp. 2–15, 2017. 113.[113]. J. J. Harris, R. Jolivet and D. Attwell, “Synaptic energy use and supply,” Neuron, vol. 75, no. 5, pp. 762–777, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2012.08.019&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22958818&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308684300005&link_type=ISI) 114.[114]. W. Heiss and K. Herholz, “Brain receptor imaging,” Journal of Nuclear Medicine, vol. 47, no. 2, pp. 302–312, 2006. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam51bWVkIjtzOjU6InJlc2lkIjtzOjg6IjQ3LzIvMzAyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjIvMjAyMy4wMy4yMi4yMzI4NzU3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 115.[115]. M. V. D. Heuvel and H. Pol, “Exploring the brain network: a review on resting-state fMRI functional connectivity,” European neuropsychopharmacology, vol. 20, no. 8, pp. 519–534, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.euroneuro.2010.03.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20471808&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 116.[116]. M. R. Hynd, H. L. Scott and P. R. Dodd, “Glutamate-mediated excitotoxicity and neurodegeneration in Alzheimer’s disease,” Neurochemistry international, vol. 45, no. 5, pp. 583–595, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuint.2004.03.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15234100&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223645400001&link_type=ISI) 117.[117]. C. Iadecola, “Neurovascular regulation in the normal brain and in Alzheimer’s disease,” Nature Reviews Neuroscience, vol. 5, no. 5, pp. 347–360, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn1387&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15100718&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221158600011&link_type=ISI) 118.[118]. L. M. Ittner and J. Götz, “Amyloid-β and tau—a toxic pas de deux in alzheimer’s disease,” Nature Reviews Neuroscience, vol. 12, no. 2, p. 67–72, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn2967&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21193853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 119.[119]. C. R. Jack Jr, D. S. Knopman, W. J. Jagust, L. M. Shaw, P. S. Aisen, M. W. Weiner, R. C. Petersen and J. Q. Trojanowski, “Hypothetical model of dynamic biomarkers of the alzheimer’s pathological cascade,” The Lancet Neurology, vol. 9, no. 1, p. 119–128, 2010. 120.[120]. W. Jagust, D. Bandy, K. Chen, N. Foster, S. Landau, C. Mathis, J. Price, E. Reiman, D. Skovronsky, R. Koeppe and ADNI, “The Alzheimer’s Disease Neuroimaging Initiative positron emission tomography core,” Alzheimer’s & Dementia, vol. 6, no. 3, pp. 221–229, 2010. 121.[121]. R. Kandimalla and P. H. Reddy, “Therapeutics of neurotransmitters in alzheimer’s disease,” Journal of Alzheimer’s Disease, vol. 57, no. 4, pp. 1049–1069, 2017. 122.[122]. S. Kaur, G. DasGupta and S. Singh, “Altered Neurochemistry in Alzheimer’s Disease: Targeting Neurotransmitter Receptor Mechanisms and Therapeutic Strategy,” Neurophysiology, pp. 1–17, 2019. 123.[123]. K. Kosik, “Personalized medicine for effective alzheimer disease treatment,” JAMA Neurology, vol. 72, no. 5, pp. 497–498, 2015. 124.[124]. B. Lam, M. Masellis, M. Freedman, D. T. Stuss and S. E. Black, “Clinical, imaging, and pathological heterogeneity of the alzheimer’s disease syndrome,” Alzheimer’s research & therapy, vol. 5, no. 1, p. 1, 2013. 125.[125]. C. Lau and R. Zukin, “NMDA receptor trafficking in synaptic plasticity and neuropsychiatric disorders,” Nature Reviews Neuroscience, vol. 8, no. 6, pp. 413–426, 2007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn2153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17514195&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246658800011&link_type=ISI) 126.[126]. J. Levenga, P. Krishnamurthy, H. Rajamohamedsait, H. Wong, T. Franke, P. Cain, E. Sigurdsson and C. Hoeffer, “Tau pathology induces loss of GABAergic interneurons leading to altered synaptic plasticity and behavioral impairments,” Acta neuropathologica communiucations, vol. 1, no. 1, p. 34, 2013. 127.[127]. C. Maes, L. Hermans, L. Pauwels, S. Chalavi, I. Leunissen, O. Levin, K. Cuypers, R. Peeters, S. Sunaert, D. Mantini and N. Puts, “Age-related differences in GABA levels are driven by bulk tissue changes,” Human Brain Mapping, vol. 39, no. 9, pp. 3652–3662, 2018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.24201&link_type=DOI) 128.[128]. C. M. McCann, J. Tapia, H. Kim, J. S. Coggan and J. Lichtman, “Rapid and modifiable neurotransmitter receptor dynamics at a neuronal synapse in vivo,” Nature neuroscience, vol. 11, no. 7, p. 807, 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.2145&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18568021&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257089400018&link_type=ISI) 129.[129]. J. W. Mink, R. J. Blumenschine and D. B. Adams, “Ratio of central nervous system to body metabolism in vertebrates: Its constancy and functional basis,” American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, vol. 231, no. 3, pp. R203–R212, 1981. 130.[130]. W. Moses, “Fundamental limits of spatial resolution in PET,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 648, pp. S236–S240, 2011. 131.[131]. F. Nasrallah, J. Griffin, V. Balcar and C. Rae, “Understanding your inhibitions: modulation of brain cortical metabolism by GABAB receptors,” Journal of Cerebral Blood Flow & Metabolism, vol. 27, no. 8, pp. 1510–1520, 2007. 132.[132]. J. Pavia, M. De Ceballos and F. De La Cuesta, “Alzheimer’s disease: relationship between muscarinic cholinergic receptors, β-amyloid and tau proteins,” Fundamental & clinical pharmacology, vol. 12, no. 5, pp. 473–481, 1998. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1472-8206.1998.tb00975.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9794144&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076061800001&link_type=ISI) 133.[133]. R. C. Petersen, P. S. Aisen, L. A. Beckett, M. C. Donohue, A. C. Gamst, D. J. Harvey, C. R. Jack, W. J. Jagust, L. M. Shaw, A. W. Toga, J. Q. Trojanowski and M. W. Weiner, “Alzheimer’s disease neuroimaging initiative (ADNI): Clinical characterization,” Neurology, vol. 74, no. 3, p. 201–209, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3181cb3e25&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20042704&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 134.[134]. A. Prakash, J. Kalra, V. Mani, K. Ramasamy and A. B. A. Majeed, “Pharmacological approaches for alzheimer’s disease: Neurotransmitter as drug targets,” Expert review of neurotherapeutics, vol. 15, no. 1, pp. 53–71, 2015. 135.[135]. N. J. Schork, “Personalized medicine: Time for one-person trials,” Nature, vol. 520, no. 7549, pp. 609–611, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/520609a&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25925459&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 136.[136]. S. Verma, A. Kumar, T. Tripathi and A. Kumar, “Muscarinic and nicotinic acetylcholine receptor agonists: current scenario in Alzheimer’s disease therapy,” Journal of Pharmacy and Pharmacology, vol. 70, no. 8, pp. 985–993, 2018.. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jphp.12919&link_type=DOI) 137.[137]. R. Wang and P. H. Reddy, “Role of glutamate and nmda receptors in Alzheimer’s disease,” Journal of Alzheimer’s Disease, vol. 57, no. 4, pp. 1041–1048, 2017. 138.[138]. Y. Wang, Q. Ren, W. Gong, D. Wu, X. Tang, X. Li, F. Wu, F. Bai, L. Xu and Z. Zhang, “Escitalopram attenuates β-amyloid-induced tau hyperphosphorylation in primary hippocampal neurons through the 5-HT1A receptor mediated Akt/GSK-3β pathway,” Oncotarget, vol. 7, no. 12, p. 13328, 2016. 139.[139]. R. Whittington, L. Virág, M. Gratuze, H. Lewkowitz-Shpuntoff, M. Cheheltanan, F. Petry, I. Poitras, F. Morin and E. Planel, “Administration of the benzodiazepine midazolam increases tau phosphorylation in the mouse brain,” Neurobiology of aging, vol. 75, pp. 11–24, 2019. 140.[140]. J. Wu, S. Hussaini, I. Bastille, G. Rodriguez, A. Mrejeru, K. Rilett, D. Sanders, C. Cook, H. Fu, R. Boonen and M. Herman, “Neuronal activity enhances tau propagation and tau pathology in vivo,” Nature neuroscience, vol. 19, no. 8, pp. 1085–10, 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.4328&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27322420&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 141.[141]. L. Yang, Y. Yan, Y. Wang, X. Hu, J. Lu, P. Chan, T. Yan and Y. Han, “Gradual disturbances of the amplitude of low-frequency fluctuations (ALFF) and fractional ALFF in Alzheimer spectrum,” Frontiers in neuroscience, vol. 12, p. 975, 2018. 142.[142]. F. Yeh, T. D. Verstynen, Y. Wang, J. C. Fernández-Miranda and W.-Y. I. Tseng, “Deterministic diffusion fiber tracking improved by quantitative anisotropy,” PloS One, vol. 8, no. 11, p. e80713, 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0080713&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24348913&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 143.[143]. Y. Zang, T. Jiang, Y. Lu, Y. He and L. Tian, “Regional homogeneity approach to fMRI data,” NeuroImage, vol. 22, no. 1, pp. 394–400, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2003.12.030&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15110032&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221190200041&link_type=ISI) 144.[144]. F. Zhang, M. Gannon, Y. Chen, S. Yan, S. Zhang, W. Feng, J. Tao, B. Sha, Z. Liu, T. Saito and T. Saido, “β-amyloid redirects norepinephrine signaling to activate the pathogenic GSK3β/tau cascade,” Science Translational Medicine, vol. 12, no. 526, 2020. 145.[145]. P. Whitehouse and K.-S. Au, “Neurotransmitter receptor alterations in Alzheimer’s disease,” Senile Dementia of the Alzheimer Type, pp. 175–182, 1985. 146.[146]. A. Folch-Fortuny, F. Arteaga and A. Ferrer, “Missing data imputation toolbox for MATLAB,” Chemometrics and Intelligent Laboratory Systems, vol. 154, pp. 93–100, 2016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chemolab.2016.03.019&link_type=DOI) 147.[147]. E. Chang, M. Savage, D. Flood, J. Thomas, R. Levy, V. Mahadomrongkul, T. Shirao, C. Aoki and P. Huerta, “AMPA receptor downscaling at the onset of Alzheimer’s disease pathology in double knockin mice,” Proceedings of the National Academy of Sciences, vol. 103, no. 9, pp. 3410–3415, 2006. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTAzLzkvMzQxMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 148.[148]. R. Yasuda, M. Ikonomovic, R. Sheffield, R. Rubin, B. Wolfe and D. Armstrong, “Reduction of AMPA-selective glutamate receptor subunits in the entorhinal cortex of patients with Alzheimer’s disease pathology: a biochemical study,” Brain research, vol. 678, no. 1-2, pp. 161–167, 1995. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0006-8993(95)00178-S&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7542540&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QW66200018&link_type=ISI) 149.[149]. T. Carter, R. Rissman, A. Mishizen-Eberz, B. Wolfe, R. Hamilton, S. Gandy and D. Armstrong, “Differential preservation of AMPA receptor subunits in the hippocampi of Alzheimer’s disease patients according to Braak stage,” Experimental neurology, vol. 187, no. 2, pp. 299–309, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.expneurol.2003.12.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15144856&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221602200009&link_type=ISI) 150.[150]. E. Miller, P. Teravskis, B. Dummer, X. Zhao, R. Huganir and D. Liao, “Tau phosphorylation and tau mislocalization mediate soluble Aβ oligomer-induced AMPA glutamate receptor signaling deficits,” European Journal of Neuroscience, vol. 39, no. 7, pp. 1214–1224, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ejn.12507&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24713000&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 151.[151]. A. Limon, J. Reyes-Ruiz and R. Miledi, “Loss of functional GABAA receptors in the Alzheimer diseased brain,” Proceedings of the National Academy of Sciences, vol. 109, no. 25, pp. 10071–10076, 2012. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzI1LzEwMDcxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjIvMjAyMy4wMy4yMi4yMzI4NzU3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 152.[152]. Y. Li, H. Sun, Z. Chen, H. Xu, G. Bu and H. Zheng, “Implications of GABAergic neurotransmission in Alzheimer’s disease,” Frontiers in aging neuroscience, vol. 8, p. 31, 2016. 153.[153]. G. Bischof, F. Jessen, Fliessbach, D. J. K., J. Hammes, B. Neumaier, O. Onur, G. Fink, J. Kukolja, A. Drzezga and T. van Eimeren, “Impact of tau and amyloid burden on glucose metabolism in Alzheimer’s disease,” Annals of clinical and translational neurology, vol. 3, no. 12, pp. 934–939, 2016. 154.[154]. N. Menkes-Caspi, H. Yamin, V. Kellner, T. Spires-Jones, D. Cohen and E. Stern, “Pathological tau disrupts ongoing network activity,” Neuron, vol. 85, no. 5, pp. 959–966, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2015.01.025&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25704951&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 155.[155]. A. Bryant, M. Hu, B. Carlyle, S. Arnold, M. Frosch, S. Das, B. Hyman and R. Bennett, “Cerebrovascular Senescence Is Associated With Tau Pathology in Alzheimer’s Disease,” Frontiers in neurology, vol. 11, p. 1058, 2020. 156.[156]. K. Bangen, A. Clark, E. Edmonds, N. Evangelista, M. Werhane, K. Thomas, L. Locano, M. Tran, Z. Zlatar, D. Nation and M. Bondi, “Cerebral Blood Flow and Amyloid-β Interact to Affect Memory Performance in Cognitively Normal Older Adults,” Frontiers in Aging Neuroscience, vol. 9, p. 181, 2017. 157.[157]. C. Iadecola, “Neurovascular regulation in the normal brain and in Alzheimer’s disease,” Nature Reviews Neuroscience, vol. 5, no. 5, pp. 347–360, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn1387&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15100718&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221158600011&link_type=ISI) 158.[158]. E. Planel, T. Miyasaka, T. Launey, D. Chui, K. Tanemura, S. Sato, O. Murayama, K. Ishiguro, Y. Tatebayashi and A. Takashima, “Alterations in glucose metabolism induce hypothermia leading to tau hyperphosphorylation through differential inhibition of kinase and phosphatase activities: implications for Alzheimer’s disease,” Journal of Neuroscience, vol. 24, no. 10, pp. 2401–2411, 2004. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIyNC8xMC8yNDAxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjIvMjAyMy4wMy4yMi4yMzI4NzU3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 159.[159]. R. Desikan, C. Fan, Y. Wang, A. Schork, H. Cabral, L. Cupples, W. Thompson, L. Besser, W. Kukull and D. C. C. Holland, “Genetic assessment of age-associated Alzheimer disease risk: Development and validation of a polygenic hazard score,” PLoS medicine, vol. 14, no. 3, p. e1002258, 21 March 2017. 160.[160]. T. Aso, G. Sugihara, T. Murai, S. Ubukata, S. Urayama, T. Ueno, G. Fujimoto, D. Thuy, H. Fukuyama and K. Ueda, “A venous mechanism of ventriculomegaly shared between traumatic brain injury and normal ageing,” Brain, vol. 143, no. 6, pp. 1843–1856, June 2020. 161.[161].Alzheimer’s Disease Neuroimaging Initiative, “ADNI2 Procedures Manual,” July 2008. [Online]. Available: [https://adni.loni.usc.edu/wp-content/uploads/2008/07/adni2-procedures-manual.pdf](https://adni.loni.usc.edu/wp-content/uploads/2008/07/adni2-procedures-manual.pdf). [Accessed May 2021]. 162.[162]. L. Gibbons, A. Carle, R. Scott-Mackin, D. Harvey, S. Mukherjee, P. Insel, S. Curtis, A. Gross, R. Jones, D. Mungas, M. Weiner, P. Crane and ADNI, “Composite measures of executive function and memory: ADNI\_EF and ADNI\_Mem,” 23 October 2015. [Online]. [Accessed May 2021]. 163.[163]. S. Choi, S. Mukherjee, L. Gibbons, R. Sanders, R. Jones, D. Tommet, J. Mez, E. Trittschuh, A. Saykin, M. Lamar and L. Rabin, “Development and validation of language and visuospatial composite scores in ADNI,” Alzheimer’s & Dementia: Translational Research & Clinical Interventions, vol. 6, no. 1, p. e12072, 2020. 164.[164]. F. Mora, G. Segovia and A. del Arco, “Aging, plasticity and environmental enrichment: structural changes and neurotransmitter dynamics in several areas of the brain,” Brain research reviews, vol. 55, no. 1, pp. 78–88, 2007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.brainresrev.2007.03.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17561265&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249510600007&link_type=ISI) 165.[165]. C. Honey, O. Sporns, L. Cammoun, X. Gigandet, J. Thiran, R. Meuli and P. Hagmann, “Predicting human resting-state functional connectivity from structural connectivity,” Proceedings of the National Academy of Sciences, vol. 106, no. 6, pp. 2035–2040, 2009. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTA2LzYvMjAzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 166.[166]. S. Bhat, U. R. Acharya, Y. Hagiwara, N. Dadmehr and H. Adeli, “Parkinson’s disease: Cause factors, measurable indicators, and early diagnosis,” Computers in biology and medicine, vol. 102, pp. 234–241, 2018. 167.[167]. K. Marek, D. Jennings, S. Lasch, A. Siderowf, C. Tanner, T. Simuni, C. Coffey, K. Kieburtz, E. Flagg, S. Chowdhury, W. Poewe and et al., “The Parkinson progression marker initiative (PPMI),” Progress in neurobiology, vol. 95, no. 4, pp. 629–635, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pneurobio.2011.09.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21930184&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 168.[168]. M. Aiello, E. Salvatore, A. Cachia, S. Pappatà, C. Cavaliere, A. Prinster, E. Nicolai, M. Salvatore, J. Baron and M. Quarantelli, “Relationship between simultaneously acquired resting-state regional cerebral glucose metabolism and functional MRI: a PET/MR hybrid scanner study,” Neuroimage, no. 113, pp. 111–121, 1 June 2015. 169.[169]. F. de Vos, M. Koini, T. Schouten, S. Seiler, J. van der Grond, A. Lechner, R. Schmidt, M. de Rooij and S. Rombouts, “A comprehensive analysis of resting state fMRI measures to classify individual patients with Alzheimer’s disease.,” Neuroimage, vol. 15, no. 167, pp. 62–72, 15 February 2018. 170.[170]. W. Dong, D. Fong, J. Yoon, E. Wan, L. Bedford, E. Tang and C. Lam, “Generative adversarial networks for imputing missing data for big data clinical research,” BMC Medical Research Methodology, vol. 21, no. 1, pp. 1–10, 2021. 171.[171]. J. Yoon, J. Jordon and M. Schaar, “GAIN: Missing data imputation using generative adversarial nets,” in International Conference on Machine Learning, Stockholm, 2018. 172.[172]. K. Jellinger, “Dementia with Lewy bodies and Parkinson’s disease-dementia: current concepts and controversies,” Journal of neural transmission, vol. 125, no. 4, pp. 615–650, 2018. 173.[173]. V. Kotagal, R. Albin, M. Müller, R. Koeppe, R. Chervin, K. Frey and N. Bohnen, “Symptoms of rapid eye movement sleep behavior disorder are associated with cholinergic denervation in Parkinson disease,” Annals of Neurology, vol. 71, no. 4, pp. 560–8, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.22691&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22522445&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 174.[174]. D. Miller and J. O’Callaghan, “Biomarkers of Parkinson’s disease: present and future,” Metabolism, vol. 64, no. 3, pp. S40–6, 2015. 175.[175]. G. T. E. R. Z. M. M. R. S. L. L. H. S. D. V. K. H. G. T. J. Oliveira L, “Alpha-synuclein research: defining strategic moves in the battle against Parkinson’s disease,” npj Parkinson’s Disease, vol. 7, no. 1, pp. 1–23, 2021. 176.[176]. D. Sulzer and R. Edwards, “The physiological role of α-synuclein and its relationship to Parkinson’s Disease,” Journal of neurochemistry, vol. 150, no. 5, pp. 475–86, 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jnc.14810&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31269263&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 177.[177]. C. Atkinson-Clement, S. Pinto, A. Eusebio and O. Coulon, “Diffusion tensor imaging in Parkinson’s disease: review and meta-analysis,” Neuroimage: Clinical, vol. 16, pp. 98–110, 2017. 178.[178]. E. Sykova, “Extrasynaptic volume transmission and diffusion parameters of the extracellular space,” Neuroscience, vol. 129, no. 4, pp. 861–876, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroscience.2004.06.077&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15561404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225726000003&link_type=ISI) 179.[179]. S. Song, S. Sun, M. Ramsbottom, C. Chang, J. Russell and A. Cross, “Dysmyelination revealed through MRI as increased radial (but unchanged axial) diffusion of water,” Neuroimage, vol. 17, no. 3, pp. 1429–1436, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.2002.1267&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12414282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179012800031&link_type=ISI) 180.[180]. A. Sauerbier, P. Jenner, A. Todorova and K. Chaudhuri, “Non motor subtypes and Parkinson’s disease,” Parkinsonism & related disorders, vol. 22, pp. S41–S46, 2016. 181.[181]. H. Braak, K. Del Tredici, U. Rüb, R. De Vos, E. Steur and E. Braak, “Staging of brain pathology related to sporadic Parkinson’s disease,” Neurobiology of aging, vol. 24, no. 2, pp. 197–211, 2003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0197-4580(02)00065-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12498954&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180616100001&link_type=ISI) 182.[182]. R. Righart, V. Biberacher, L. Jonkman, R. Klaver, P. Schmidt, D. Buck, A. Berthele, J. Kirschke, C. Zimmer, B. Hemmer, J. Geurts and M. Mühlau, “Cortical pathology in multiple sclerosis detected by the T 1/T 2-weighted ratio from routine magnetic resonance imaging,” Annals of neurology, vol. 82, no. 4, pp. 519–529, 2017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.25020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28833433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 183.[183]. A. Goulas, J. Changeux, K. Wagstyl, K. Amunts, N. Palomero-Gallagher and C. Hilgetag, “The natural axis of transmitter receptor distribution in the human cerebral cortex,” Proceedings of the National Academy of Sciences, vol. 118, no. 3, p. e2020574118, 2021. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxNzoiMTE4LzMvZTIwMjA1NzQxMTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMy8yMi8yMDIzLjAzLjIyLjIzMjg3NTc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 184.[184]. A. Butt, R. Fern and C. Matute, “Neurotransmitter signaling in white matter,” Glia, vol. 62, no. 11, pp. 1762–1779, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/glia.22674&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24753049&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 185.[185]. D. Jones, T. Knösche and R. Turner, “White matter integrity, fiber count, and other fallacies: the do’s and don’ts of diffusion MRI,” Neuroimage, vol. 73, pp. 239–254, 2013. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.06.081&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22846632&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317084500024&link_type=ISI) 186.[186]. J. Yoon, J. Jordon and M. Schaar, “GAIN: Missing data imputation using generative adversarial nets,” in Proceedings of International Conference on Machine Learning, 2018. 187.[187]. I. Rektor, A. Svátková, L. Vojtíšek, I. Zikmundová, J. Vaníček, A. Király and N. Szabó, “White matter alterations in Parkinson’s disease with normal cognition precede grey matter atrophy,” PloS one, vol. 13, no. 1, p. e0187939, 2018. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 188.[188]. M. Politis, K. Wu, C. Loane, L. Kiferle, S. Molloy, D. Brooks and P. Piccini, “Staging of serotonergic dysfunction in Parkinson’s disease: an in vivo 11C-DASB PET study,” Neurobiology of disease, vol. 40, no. 1, pp. 216–221, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nbd.2010.05.028&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20594979&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 189.[189]. R. Albin, M. Müller, N. Bohnen, C. Spino, M. Sarter, R. Koeppe, A. Szpara, K. Kim, C. Lustig and W. Dauer, “α4β2* Nicotinic Cholinergic Receptor Target Engagement in Parkinson Disease Gait–Balance Disorders,” Annals of Neurology, vol. 90, no. 1, pp. 130–142, 2021. 190.[190]. S. Roshanbin, M. Xiong, G. Hultqvist, L. Söderberg, O. Zachrisson, S. Meier, S. Ekmark-Lewén, J. Bergström, M. Ingelsson, D. Sehlin and S. Syvänen, “In vivo imaging of alpha-synuclein with antibody-based PET,” Neuropharmacology, vol. 208, p. 1089, 2022. 191.[191]. S. Carriedo, H. Yin and J. Weiss, “Motor neurons are selectively vulnerable to AMPA/kainate receptor-mediated injury in vitro,” Journal of Neuroscience, vol. 16, no. 13, pp. 4069–4079, 1996. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIxNi8xMy80MDY5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjIvMjAyMy4wMy4yMi4yMzI4NzU3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 192.[192]. S. Rao, S. Banack, P. Cox and J. Weiss, “BMAA selectively injures motor neurons via AMPA/kainate receptor activation,” Experimental neurology, vol. 201, no. 1, pp. 244–252, 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.expneurol.2006.04.017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16764863&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240152100026&link_type=ISI) 193.[193]. P. Seeman and H. Niznik, “Dopamine receptors and transporters in Parkinson’s disease and schizophrenia,” The FASEB Journal, vol. 4, no. 10, pp. 2737–2744, 1990. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2197154&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DP78300009&link_type=ISI) 194.[194]. B. Noudoost and T. Moore, “Control of visual cortical signals by prefrontal dopamine,” Nature, vol. 474, no. 7351, pp. 372–375, 2011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09995&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21572439&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291647100044&link_type=ISI) 195.[195]. F. Zhou, K. Lesch and D. Murphy, “Serotonin uptake into dopamine neurons via dopamine transporters: a compensatory alternative,” Brain research, vol. 942, no. 1-2, pp. 109–119, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-8993(02)02709-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12031859&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176397900013&link_type=ISI) 196.[196]. E. Foran and D. Trotti, “Glutamate transporters and the excitotoxic path to motor neuron degeneration in amyotrophic lateral sclerosis,” Antioxidants & redox signaling, vol. 11, no. 7, pp. 1587–1602, 2009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/ars.2009.2444&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19413484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266373100009&link_type=ISI) 197.[197]. G. Ambrosi, S. Cerri and F. Blandini, “A further update on the role of excitotoxicity in the pathogenesis of Parkinson’s disease,” Journal of neural transmission, vol. 121, no. 8, pp. 849–859, 2014. 198.[198]. F. Benedetti, I. Bollettini, S. Poletti, C. Locatelli, C. Lorenzi, A. Pirovano, E. Smeraldi and C. Colombo, “White matter microstructure in bipolar disorder is influenced by the serotonin transporter gene polymorphism 5-HTTLPR,” Genes, Brain and Behaviour, vol. 14, no. 3, pp. 238–250, 2015. 199.[199]. G. Alexopoulos, C. Murphy, F. Gunning-Dixon, C. Glatt, V. Latoussakis, R. Kelly Jr, D. Kanellopoulos, S. Klimstra, K. Lim, R. Young and M. Hoptman, “Serotonin transporter polymorphisms, microstructural white matter abnormalities and remission of geriatric depression,” Journal of affective disorders, vol. 119, no. 1-3, pp. 132–41, 2009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2009.03.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19375170&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 200.[200]. C. Lochner, J. Fouché, S. du Plessis, B. Spottiswoode, S. Seedat, N. Fineberg, S. Chamberlain and D. Stein, “Evidence for fractional anisotropy and mean diffusivity white matter abnormalities in the internal capsule and cingulum in patients with obsessive–compulsive disorder,” Journal of Psychiatry and Neuroscience, vol. 37, no. 3, pp. 193–199, 2012. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianBuIjtzOjU6InJlc2lkIjtzOjg6IjM3LzMvMTkzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjIvMjAyMy4wMy4yMi4yMzI4NzU3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 201.[201]. G. Hardingham and H. Bading, “Synaptic versus extrasynaptic NMDA receptor signalling: implications for neurodegenerative disorders,” Nature Reviews Neuroscience, vol. 11, no. 10, pp. 682–696, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn2911&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20842175&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281928500002&link_type=ISI) 202.[202]. A. Butt, R. Fern and C. Matute, “Neurotransmitter signaling in white matter,” Glia, vol. 62, no. 11, pp. 1762–1779, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/glia.22674&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24753049&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 203.[203]. N. Bohnen, D. Kaufer, L. Ivanco, B. Lopresti, R. Koeppe, J. Davis, C. Mathis, R. Moore and S. DeKosky, “Cortical cholinergic function is more severely affected in parkinsonian dementia than in Alzheimer disease: an in vivo positron emission tomographic study,” Archives of neurology, vol. 60, no. 12, pp. 1745–1748, 2003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archneur.60.12.1745&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14676050&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187178700012&link_type=ISI) 204.[204]. P. Tiraboschi, L. Hansen, M. Alford, M. Sabbagh, B. Schoos, E. Masliah, L. Thal and J. Corey-Bloom, “Cholinergic dysfunction in diseases with Lewy bodies,” Neurology, vol. 54, pp. 407–411, 2000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.54.2.407&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10668703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 205.[205]. G. Pagano, F. Niccolini, P. Fusar-Poli and M. Politis, “Serotonin transporter in Parkinson’s disease: A meta-analysis of positron emission tomography studies,” Annals of neurology, vol. 81, no. 2, pp. 171–180, 2017. 206.[206]. M. van den Hurk, S. Lau, M. Marchetto, J. Mertens, S. Stern, O. Corti, A. Brice, B. Winner, J. Winkler, F. Gage and C. Bardy, “Druggable transcriptomic pathways revealed in Parkinson’s patient-derived midbrain neurons,” npj Parkinson’s Disease, vol. 8, no. 1, pp. 1–18, 2022. 207.[207]. M. Asahina, T. Suhara, H. Shinotoh, O. Inoue, K. Suzuki and T. Hattori, “Brain muscarinic receptors in progressive supranuclear palsy and Parkinson’s disease: a positron emission tomographic study,” Journal of Neurology, Neurosurgery & Psychiatry, vol. 65, no. 2, pp. 155–163, 1998. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI2NS8yLzE1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzIyLzIwMjMuMDMuMjIuMjMyODc1NzguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 208.[208]. N. Bohnen, M. Müller, V. Kotagal, R. Koeppe, M. Kilbourn, S. Gilman, R. Albin and K. Frey, “Heterogeneity of cholinergic denervation in Parkinson’s disease without dementia,” Journal of Cerebral Blood Flow & Metabolism, vol. 32, no. 8, pp. 1609–1617, 2012. 209.[209]. S. Pimlott, M. Piggott, J. Owens, E. Greally, J. Court, E. Jaros, R. Perry, E. Perry and D. Wyper, “Nicotinic acetylcholine receptor distribution in Alzheimer’s disease, dementia with Lewy bodies, Parkinson’s disease, and vascular dementia: in vitro binding study using 5-[125I]-A-85380,” Neuropsychopharmacology, vol. 29, no. 1, pp. 108–116, 2004. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.npp.1300302&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12955099&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187993600013&link_type=ISI) 210.[210]. S. van der Zee, P. Kanel, M. Gerritsen, J. Boertien, A. Slomp, M. Müller, N. Bohnen, J. Spikman and T. van Laar, “Altered Cholinergic Innervation in De Novo Parkinson’s Disease with and Without Cognitive Impairment,” Movement Disorders, vol. 37, no. 4, pp. 713–723, 2022. 211.[211]. B. Chen and M. Rice, “Synaptic regulation of somatodendritic dopamine release by glutamate and GABA differs between substantia nigra and ventral tegmental area,” Journal of neurochemistry, vol. 81, no. 1, pp. 158–169, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1471-4159.2002.00811.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12067228&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174784200016&link_type=ISI) 212.[212]. K. Kawabata and H. Tachibana, “Evaluation of benzodiazepine receptor in the cerebral cortex of Parkinson’s disease using 123I-iomazenil SPECT,” Nihon rinsho. Japanese Journal of Clinical Medicine, vol. 55, no. 1, pp. 244–248, 1997. 213.[213]. K. Jellinger, “Is Braak staging valid for all types of Parkinson’s disease?,” Journal of Neural Transmission, vol. 126, no. 4, pp. 423–431, 2019. 214.[214]. J. Klein, C. Eggers, E. Kalbe, S. Weisenbach, C. Hohmann, S. Vollmar, S. Baudrexel, N. Diederich, W. Heiss and R. Hilker, “Neurotransmitter changes in dementia with Lewy bodies and Parkinson disease dementia in vivo,” Neurology, vol. 74, no. 11, pp. 885–892, 2010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3181d55f61&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20181924&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F22%2F2023.03.22.23287578.atom) 215.[215]. S. Pandya, Y. Zeighami, B. Freeze, M. Dadar, D. Collins, A. Dagher and A. Raj, “Predictive model of spread of Parkinson’s pathology using network diffusion,” NeuroImage, vol. 192, pp. 178–194, 15 May 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2019.03.001&link_type=DOI) [1]: /embed/graphic-8.gif [2]: /embed/graphic-9.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/graphic-10.gif [5]: /embed/inline-graphic-2.gif [6]: /embed/graphic-11.gif [7]: /embed/graphic-12.gif [8]: /embed/graphic-13.gif [9]: /embed/inline-graphic-3.gif [10]: /embed/inline-graphic-4.gif [11]: /embed/inline-graphic-5.gif [12]: /embed/inline-graphic-6.gif [13]: /embed/inline-graphic-7.gif [14]: /embed/graphic-14.gif [15]: /embed/graphic-15.gif [16]: /embed/graphic-16.gif [17]: /embed/graphic-17.gif [18]: /embed/graphic-18.gif [19]: /embed/inline-graphic-8.gif [20]: /embed/inline-graphic-9.gif [21]: /embed/inline-graphic-10.gif