Abstract
We used a network approach to assess systems-level abnormalities in motor activation in humans with Parkinson's disease (PD). This was done by measuring the expression of the normal movement-related activation pattern (NMRP), a previously validated activation network deployed by healthy subjects during motor performance. In this study, NMRP expression was prospectively quantified in 15O-water PET scans from a PD patient cohort comprised of a longitudinal early-stage group (n = 12) scanned at baseline and at two or three follow-up visits two years apart, and a moderately advanced group scanned on and off treatment with either subthalamic nucleus deep brain stimulation (n = 14) or intravenous levodopa infusion (n = 14). For each subject and condition, we measured NMRP expression during both movement and rest. Resting expression of the abnormal PD-related metabolic covariance pattern was likewise determined in the same subjects. NMRP expression was abnormally elevated (p < 0.001) in PD patients scanned in the nonmovement rest state. By contrast, network activity measured during movement did not differ from normal (p = 0.34). In the longitudinal cohort, abnormal increases in resting NMRP expression were evident at the earliest clinical stages (p < 0.05), which progressed significantly over time (p = 0.003). Analogous network changes were present at baseline in the treatment cohort (p = 0.001). These abnormalities improved with subthalamic nucleus stimulation (p < 0.005) but not levodopa (p = 0.25). In both cohorts, the changes in NMRP expression that were observed did not correlate with concurrent PD-related metabolic covariance pattern measurements (p > 0.22). Thus, the resting state in PD is characterized by changes in the activity of normal as well as pathological brain networks.
Introduction
Network analysis of functional imaging data acquired in the rest state has provided valuable insights into the network abnormalities underlying Parkinson's disease (PD) and the modulation of these changes by treatment (Eidelberg, 2009). For example, the PD-related spatial covariance pattern (PDRP), an abnormal large-scale metabolic network associated with the akinetic-rigid manifestations of the disorder (Spetsieris and Eidelberg, 2011), has been linked to disease-related changes in basal ganglia output pathways (e.g., Lin et al., 2008; Mure et al., 2011; Niethammer and Eidelberg, 2012).
Although neurodegenerative disorders can be associated with the expression of pathological metabolic patterns in the rest state, the underlying disease process can also influence the activity of normal brain networks (i.e., those routinely deployed by healthy individuals during task performance) (e.g., Nakamura et al., 2001; Mentis et al., 2003b; compare Pievani et al., 2011). In the case of PD, the status of networks normally activated during the execution of simple movements becomes relevant. The akinetic-rigid manifestations of PD have been linked to failed suppression of cortical movement-related neural activity in the resting condition (Ridding et al., 1995; Berardelli et al., 1996; Van Der Werf and Paus, 2006; Van Der Werf et al., 2006). By the same token, treatment with either dopaminergic agents or with deep brain stimulation (DBS) can restore cortical inhibition to varying degrees (Chen et al., 2001; Pierantozzi et al., 2002; Lefaucheur, 2005; Fraix et al., 2008). It is not known, however, whether analogous network-level changes can be discerned in the PD rest state.
To address this issue, we measured the expression of a specific movement-related activation network that we have previously reported in healthy subjects (Carbon et al., 2010; Mure et al., 2012). Using a supervised multivariate approach (Habeck et al., 2005) to analyze 15O-water (H215O) PET scans from normal volunteers acquired during movement and rest, we identified a highly replicable activation network, which was termed the normal movement-related activation pattern (NMRP). Apart from delineating the spatial topography of this covariance pattern, the algorithm was used to quantify pattern expression on a prospective scan basis, yielding separate network activity measurements for the movement and rest states. In this context, failure to inhibit motor network activity in the absence of movement would be manifest by resting elevations in NMRP expression. To test this hypothesis, we measured network activity in PD patients scanned at rest and during movement. In addition to comparing these measures to corresponding healthy control values, we evaluated the network changes that occurred with advancing disease and in response to treatment. Last, we correlated the observed NMRP changes with concurrent PDRP measurements from the same subjects and assessed the relationship of the two networks to clinical outcome.
Materials and Methods
Subjects
We studied three groups of PD patients with akinetic-rigid predominant symptoms. A diagnosis of PD was made clinically according to United Kingdom Parkinson's Disease Society Brain Bank criteria (Hughes et al., 1992). Patients with tremor-dominant clinical manifestations (composite Unified Parkinson's Disease Rating Scale [UPDRS] limb tremor ratings ≥4 with at least one limb ≥2) (Isaias et al., 2010; Mure et al., 2011) were excluded as were subjects with Mini-Mental State Examination scores <27. Further exclusion criteria included severe hypertension, cardiovascular disease, diabetes mellitus, and past or current psychiatric history. Review of routine anatomical magnetic resonance imaging from the subjects revealed no evidence of atrophy or incidental structural brain abnormalities. Ethical permission for all studies was obtained from the Institutional Review Board of North Shore University Hospital. Written consent was obtained from each subject after detailed explanation of the procedures.
Group 1: longitudinal cohort.
The clinical features of the longitudinal cohort are presented in Table 1. This cohort consisted of 12 right-handed early-stage PD patients (age 60.2 ± 2.0 years, mean ± SD; off-state UPDRS motor ratings 9.0 ± 1.4) with predominantly akinetic-rigid clinical manifestations. These subjects underwent H215O PET at baseline (TP1) and again at a second time point (TP2) 2.4 ± 0.3 years later. Six of these subjects were rescanned at a third time point (TP3) 4.4 ± 0.3 years after baseline. All subjects in this group had unilateral motor signs (“hemiparkinsonism”) at baseline. At TP1, seven of the subjects exhibited clinical manifestations on the right side; the remainder had signs on the left side. At TP2, motor ratings increased to 14.4 ± 1.5 and to 18.4 ± 2.9 for those scanned at TP3. At baseline, five subjects were drug-naive, six were on levodopa/carbidopa, dopamine agonists, or both; the remaining subject was on deprenyl. By TP2, all subjects but one were on chronic oral dopaminergic treatment. In the treated patients, antiparkinsonian medications were discontinued 12 h before PET imaging. At each longitudinal time point, the subjects additionally underwent metabolic brain imaging with [18F]-fluorodeoxyglucose (FDG) PET for the measurement of PDRP expression in the rest state. Limited imaging data from these subjects have been reported previously (Carbon et al., 2007; Huang et al., 2007; Tang et al., 2010).
Demographic and clinical featuresa
Group 2: treatment cohort.
The clinical features of this cohort are presented in Table 1. The first treatment subgroup (Group 2A) consisted of 14 right-handed moderately advanced PD patients (age 58.6 ± 2.8 years; off-state motor UPDRS 25.5 ± 2.2). These patients were on chronic bilateral subthalamic nucleus (STN) stimulation for predominant akinesia-rigidity; they underwent H215O PET at baseline and again during stimulation. The baseline (OFF) condition was defined as 12 h after the cessation of medication and at least 1 h after stimulators were turned off. The treatment (DBS) condition was defined by adjusting the stimulation parameters on both sides to achieve an improvement in motor UPDRS ratings of at least 5 points and/or a reduction of 20% relative to the baseline condition. In this subgroup, treatment was associated with a change in motor UPDRS ratings of −7.4 ± 1.3 points (−20.2%, p < 0.001; paired Student's t test). Of these subjects, 10 additionally underwent FDG PET imaging to measure resting PDRP expression in the two treatment conditions. Limited imaging data from these subjects have been reported previously (Asanuma et al., 2006; Hirano et al., 2008; Mure et al., 2012).
The second treatment subgroup (Group 2B) consisted of 14 right-handed moderately advanced PD patients (age 57.1 ± 2.0 years; off-state motor UPDRS 25.4 ± 2.0) who were scanned in a baseline condition (OFF) after a 12 h medication washout and during an intravenous levodopa infusion (LD) that was titrated to achieve a stable treatment response with minimal or no dyskinesia. The details of the levodopa infusion protocol have been provided previously (Hirano et al., 2008; Mure et al., 2012). In this subgroup, treatment was associated with a change in motor UPDRS ratings of −8.9 ± 1.3 points (−35.1%, p < 0.001; paired Student's t test). Of these subjects, six also underwent FDG PET at baseline and during levodopa infusion. Limited imaging data from these subjects have been reported previously (Argyelan et al., 2008; Hirano et al., 2008; Mure et al., 2012).
Group 3: healthy control subjects.
This group consisted of eight right-handed normal volunteers (age 53.8 ± 4.7 years) who were scanned with H215O PET at a single time point. These individuals were not among the 18 healthy subjects whose scans were used originally for NMRP derivation (Carbon et al., 2010). Rather, H215O PET scans from the current normal control group were used to measure network activity prospectively on a subject-by-subject basis (see below). The resulting NMRP expression values were compared with the corresponding measures from the PD patients computed in the different experimental conditions.
Imaging procedures
PET imaging was performed in 3D mode using a GE Advance tomograph (General Electric). The 18-ring bismuth germanate scanner provided 35 image planes with an axial field of view of 14.5 cm and an intrinsic resolution of 4.2 mm (FWHM) in all directions. PET scans were conducted in a dimly lit room with minimal auditory stimulation with eyes open. Patients were positioned in the scanner using a stereoadapter with 3D laser alignment with reference to the orbitomeatal line; identical stereoadapter and laser settings were used in each imaging session.
Subjects were scanned with H215O PET in each experimental condition (time point for Group 1; treatment condition for Group 2; baseline session for Group 3) while performing a kinematically controlled motor execution task (MOVE) in which reaching movements were made toward a circular array of targets (Ghilardi et al., 2000). In MOVE, the subjects moved a cursor on a digitizing tablet with the dominant right hand. Movements were out and back from a central starting position to one of eight radial targets displayed on the screen. The eight targets appeared in a predictable counter-clockwise order. To reach the target in synchrony with the tone, subjects had to initiate movement before it appeared. For each subject, the extent and rate of the movements were adjusted to be kinematically equivalent across trials and conditions (e.g., Argyelan et al., 2008; Mure et al., 2012). The subjects were also scanned during NO MOVE, a nonmovement reference condition in which they passively observed randomly appearing targets. All longitudinal (Group 1) and levodopa infusion (Group 2B) subjects were scanned during both MOVE and NO MOVE. All 14 members of the STN stimulation cohort (Group 2A) were scanned without movement in the OFF and DBS conditions; eight were scanned in the two treatment conditions during both movement and rest.
Network analysis
NMRP.
We used Ordinal Trends Canonical Variates Analysis (OrT/CVA), a within-group network modeling approach (Habeck et al., 2005; Moeller and Habeck, 2006) (software available at http://groups.google.com/group/gCVA), to identify and validate the spatial covariance topography associated with motor activation in healthy subjects. OrT/CVA is a specialized form of supervised principal component analysis that searches for the presence of specific patterns of functional connectivity (i.e., large-scale brain networks) in serial imaging data acquired over multiple ordered experimental conditions. The analysis searches for spatial covariance patterns (i.e., the principal components) for which the associated subject scores (i.e., the principal component scalars) exhibit monotonically increasing expression across conditions, even as the functional relationship between the brain regions comprising the network remains constant (i.e., is group invariate). This approach is fundamentally different from routine mass-univariate analysis in which task-related changes in individual regions (clusters) are evaluated independently at the group mean level (Moeller and Habeck, 2006). OrT/CVA identifies brain networks that, if present in the data, exhibit consistent increases in subject expression across tasks or experimental conditions in all or most cases. This property is known as an ordinal trend. The presence of a significant ordinal trend in a given dataset is determined nonparametrically by permutation testing of the relevant subject scores to reject the possibility that the observed changes in pattern expression across conditions had occurred by chance. Likewise, the reliability of the regional contributions to the network topography (i.e., the voxel loadings on the pattern) is assessed by bootstrap estimation (Habeck and Stern, 2010). Importantly, the test-retest reliability of pattern expression within individual subjects is evaluated across repeat runs obtained in each condition.
In addition to quantifying the expression of the resulting pattern(s) in each subject and condition in the derivation cohort, OrT/CVA can also be used prospectively to measure network activity in new subjects/conditions on an individual scan basis. In the current study, we used the previously validated NMRP topography (Fig. 1A) identified by OrT/CVA of H215O PET data from 18 healthy subjects scanned during both MOVE and NO MOVE (Carbon et al., 2010; Mure et al., 2012). This pattern exhibited only minimal correspondence with the PDRP topography (Fig. 1B), sharing only 7.4% of the variance in their respective voxel weights (p = 0.176, adjusted; see below). Here, we quantified NMRP activity in the MOVE and NO MOVE scans (referred to as NMRP [MOVE] and NMRP [NO MOVE] values, respectively) from all groups, subjects, and experimental conditions. For the early-stage longitudinal cohort (Group 1), network values were computed in MOVE and NO MOVE scans at each time point. For the treatment cohort (Group 2), NMRP [MOVE] and NMRP [NO MOVE] values were computed in scans acquired at baseline and during treatment with either STN stimulation (Group 2A) or levodopa infusion (Group 2B). We also computed NMRP [MOVE] and NMRP [NO MOVE] values in the scans acquired in the new healthy volunteer cohort (Group 3). To standardize measurements of network activity across subjects, tasks, and experimental conditions, computed values were z-scored with respect to NMRP expression in the NO MOVE scans (n = 39) of the original derivation cohort such that the mean for this sample was 0 and the SD was 1.
Spatial covariance patterns: the NMRP and PDRP topographies. A, NMRP identified by network analysis of 39 MOVE and NO MOVE H215O PET scan pairs from 18 healthy volunteer subjects (Carbon et al., 2010; Mure et al., 2012). This spatial covariance pattern was characterized by activation of the left sensorimotor cortex (SMC), premotor cortex, and inferior parietal cortex. Additional contributions to this activation topography were found bilaterally in the cerebellar vermis and paramedian zone. The displayed voxel loadings on the pattern were shown to be reliable (p < 0.001) using bootstrap estimation procedures. B, PDRP identified by spatial covariance analysis of FDG PET metabolic images from a combined group of 33 PD patients and 33 controls (Ma et al., 2007; Eidelberg, 2009). This pattern was characterized by increased (red) pallidal, thalamic, and motor cortical metabolic activities associated with reduced (blue) lateral premotor and parieto-occipital cortical activities. The displayed voxel loadings on the pattern were shown to be reliable (p < 0.001) by bootstrap estimation. C, Voxelwise correlation of standardized regional loadings on the two spatial covariance patterns (Mure et al., 2011). The NMRP and PDRP exhibited minimal spatial correspondence, with <7.5% of the total variation in voxel weights shared by the two topographies. This correlation did not differ from chance (p = 0.176).
For validation, we demonstrated the presence of a significant ordinal trend (NMRP [MOVE] > NMRP [NO MOVE]) in the current dataset. This was done separately for each group and experimental condition by evaluating each MOVE and NO MOVE scan pair individually and counting the number of exceptions (violations) to the ordinal trend (i.e., scan pairs for which NMRP [MOVE] ≤ NMRP [NO MOVE]) in relation to the total number of trials. In each sample, the null hypothesis that the observed distribution of NMRP [MOVE] − NMRP [NO MOVE] differences was explainable by chance was examined using binomial tests, which were considered significant for p < 0.01 (Mure et al., 2012). To evaluate the stability of the NMRP measurements, we performed a test-retest analysis of the values obtained in repeat trials from the same subject. The test-retest reproducibility of NMRP [MOVE] and NMRP [NO MOVE] values were separately evaluated by computing intraclass correlations (ICC) for each measurement. Values were considered reliable for p < 0.05.
Last, to determine whether intercurrent tremor influences network expression, we separately compared NMRP [NO MOVE] and NMRP [MOVE] values measured in the atremulous PD subjects (n = 10), defined as those with tremor subscale ratings (UPDRS III item 20 and 21) of 0, with corresponding measurements from a subgroup of PD subjects (n = 10) with observable tremor who were matched for both age and composite UPDRS motor ratings (Table 2). Differences in NMRP expression measured in PD subjects with and without tremor, and in healthy control subjects (n = 8), were evaluated using one-way ANOVA and post hoc Bonferroni tests. Differences between groups were considered significant for p < 0.05.
Atremulous and tremulous patientsa
PDRP expression.
In addition to quantifying NMRP activity in MOVE and NO MOVE scans, we also measured resting PDRP expression in FDG PET scans from these subjects. For consistency with previously published measures of PDRP expression, the resulting values were z-scored with respect to the values computed in the FDG PET scans of the normal subjects (n = 33) used in the original pattern derivation (Ma et al., 2007; Spetsieris and Eidelberg, 2011). All single-case computations were performed using a fully automated voxel-based algorithm (software available at http://www.fil.ion.ucl.ac.uk/spm/ext/#SSM), blind to group and experimental condition.
Data analysis
Comparison of the NMRP and PDRP topographies.
To evaluate spatial similarities and differences between the NMRP and PDRP topographies, we compared the regional weights on the two spatial covariance patterns using a voxelwise correlational approach (Mure et al., 2011). Spatial autocorrelation can produce chance correlations within the network volumes (Liebhold and Sharov, 1998). To address this concern, we conducted a set of simulations to determine the number of times a correlation of the observed magnitude (i.e., r2 for the correlation of the NMRP and PDRP voxel weights) occurred by chance. To do this, we constructed 1000 pseudo-random volume pairs; each volume was comprised of 116 regions defined by the automated anatomical labeling algorithm (Tzourio-Mazoyer et al., 2002). Within a given volume, each region was assigned pseudo-random numbers (Gaussian distribution with mean of 0 and SD of 1). Gaussian noise (mean of 0 and SD of 0.05) was added to each volume and smoothed with a box filter of increasing kernel size (5 × 5 × 5 to 45 × 45 × 45 voxels) to produce an autocorrelation of similar degree to that observed empirically for the two patterns. (The spatial autocorrelation of each volume was estimated by the average of local Moran's I calculated within a test space defined by a common mask constructed to include both patterns.) In each iteration, the generated volumes were correlated; this simulation procedure was repeated 1000 times. We then calculated the rank of the r2 value that corresponded to the magnitude of the correlation (r2) that was directly observed between the two region weight vectors. In the correlation between the NMRP and PDRP voxel weights, 176 volume-pairs exhibited correlations of magnitude greater than the observed r2 value of 0.074, which corresponded to p = 0.176. Thus, the correlation of two network topographies was therefore not statistically significant.
Changes across groups and experimental conditions.
NMRP expression values were computed separately in the MOVE and NO MOVE scans obtained for each subject and experimental condition. Additionally, these values were used to compute ΔNMRP, the MOVE − NO MOVE difference in NMRP expression, corresponding to the network-level motor activation response for each subject and condition. Group differences in NMRP [MOVE], NMRP [NO MOVE], and ΔNMRP [MOVE − NO MOVE] were assessed separately by comparing the respective network values for the early-stage (Group 1 scans averaged across the three longitudinal time points) and the more advanced (Group 2 scans acquired in the baseline untreated condition) PD cohorts, and in the healthy control group (Group 3) using one-way ANOVA with post hoc Bonferroni tests. This analysis is repeated using subject age as a covariate. We additionally compared the group differences in NMRP [MOVE] expression with those for NMRP [NO MOVE] to determine whether a significant group × task interaction effect was present in the data. This was accomplished using 2 × 3 repeated-measures ANOVA.
Within-group analyses were performed to assess the changes in network expression that occurred over time (Group 1) and during treatment (Group 2). For Group 1 patients, longitudinal changes in pattern expression were evaluated for NMRP [MOVE], NMRP [NO MOVE], and ΔNMRP [MOVE − NO MOVE]. The time course of each network measure was evaluated using the general linear model (GLM) (McClullagh and Nelder, 1989). The trajectories of the values over time were compared for each pair of network activity measures by computing Bland-Altman within-subject correlation coefficients (Bland and Altman, 1995). Additionally, cross-sectional relationships between NMRP and PDRP expression values were separately examined for each longitudinal time point by computing Pearson product-moment correlation coefficients.
Treatment effects.
A similar approach was used to evaluate changes in NMRP expression during treatment. For each subgroup (Group 2A: STN stimulation; Group 2B: levodopa infusion), NMRP [MOVE], NMRP [NO MOVE], and ΔNMRP [MOVE − NO MOVE] values were computed in the baseline (OFF) and treatment (DBS or LD) conditions, and compared with corresponding healthy control values using Student's t tests. Treatment-mediated changes in each network measure were evaluated using paired Student's t tests. Correlations between NMRP and PDRP expression were assessed in the subjects who underwent imaging with both H215O and FDG PET; separate Pearson product-moment correlation coefficients were computed for scans acquired in the baseline off-state and during treatment. The relationship between treatment-mediated changes in the expression of the two networks was evaluated using correlation analysis. This analysis was limited to the subset of Group 2 subjects (10 STN DBS; 6 levodopa) who underwent imaging with both FDG and H215O PET in the baseline and treatment conditions.
Clinical correlates.
We used the GLM to determine whether the changes in UPDRS motor ratings over time (Group 1) or with treatment (Groups 2) correlated with concurrent measurements of NMRP [NO MOVE] and/or PDRP network activity. The following linear models were evaluated for the two groups separately: (1) UPDRS = B*subjects[2…n] + c; (2) UPDRS = B*subjects[2…n] + b1*PDRP + c; (3) UPDRS = B*subjects[2…n] + b1*NMRP + c; and (4) UPDRS = B*subjects[2…n] + b1*PDRP + b2*NMRP + c.
The models were tested hierarchically to determine whether clinical outcome was more accurately predicted using the two networks together as opposed to each separately. This approach provided quantitative information concerning the relative contributions of each network predictor to clinical severity as well as potential overlap between these variables. The Matlab 7.7.0 statistical toolbox was used for the GLM analyses. SPSS 13.0 for Windows was used for the remaining statistical analyses. Results were considered significant for p < 0.05.
Results
Quantification of NMRP expression in individual subjects: validation
We first verified that the NMRP was consistently activated by movement in each of the groups and experimental conditions analyzed in this study. To confirm the presence of an ordinal trend in NMRP expression across tasks, we quantified network activity in each MOVE/NO MOVE scan pair and determined the number of violations (i.e., trials in which NMRP [MOVE] ≤ NMRP [NO MOVE]) that were present in each of the samples (see Materials and Methods). Indeed, NMRP expression in MOVE exceeded NO MOVE in each group/condition (Fig. 2) with few, if any, violations (Group 1: TP1, 1 violation in 22 [1 of 22] trials, p < 0.001; TP2, 0 of 17, p < 0.001; TP3, 0 of 10, p = 0.002; Group 2A: OFF, 0 of 8; DBS, 0 of 8, p = 0.008; Group 2B: OFF, 0 of 20; LD, 2 of 20, p < 0.001; Group 3: 0 of 16, p < 0.001; binomial tests). We next assessed the within-subject stability of prospectively computed NMRP expression values. These measures exhibited excellent test-retest reliability in repeat trials conducted during MOVE (ICC = 0.89, p < 0.001; 55 test-retest scan pairs) and NO MOVE (ICC = 0.91, p < 0.001; 51 test-retest scan pairs). Because of the close correlation that was found in repeat scans obtained in a given subject and condition, the corresponding network values were averaged for further analysis.
Prospective computations of NMRP expression: confirmation of ordinal trend. Prospectively computed measures of NMRP expression in MOVE and NO MOVE scans from the following: (A) healthy volunteer subjects, (B) early-stage PD patients scanned at each time point, (C) moderate-stage PD patients scanned on and off STN stimulation, and (D) moderate-stage PD patients scanned on and off levodopa infusion. The prospective healthy control group exhibited an increase in NMRP expression in all MOVE/NO MOVE trials. In the PD group, the MOVE > NO MOVE ordinal trend was confirmed in each of the groups and experimental conditions (p < 0.01; binomial tests), with three violations of a total of 105 trials (see Results).
Last, because different normative samples were used to standardize measurements of NMRP and PDRP expression, we considered the possibility of bias in the comparisons of the corresponding subject scores. We therefore compared z-scored NMRP [NO MOVE] and resting PDRP expression values measured in the six Group 3 normatives who were scanned with both H215O and FDG PET. The resulting network values disclosed no evidence of bias (NMRP [NO MOVE]: 0.24 ± 0.26; PDRP: 0.22 ± 0.25, p = 0.93) from the different control samples used to standardize the two variables. These values were therefore directly compared without recalibration.
Differences in NMRP expression across groups
NMRP expression measured in the NO MOVE scans of early-stage and moderately advanced PD patients (Fig. 3A) differed from corresponding healthy control values (F(2,37) = 7.881, p = 0.001; one-way ANOVA), with abnormal elevations in both patient groups relative to healthy control values (early PD: p = 0.028; moderate PD: p = 0.001; post hoc Bonferroni tests). This group difference remained significant after adjustment for variation in subject age (F(2,36) = 7.184, p = 0.002). By contrast, the profile of pattern expression across groups for NMRP [MOVE] (Fig. 3B) differed from NMRP [NO MOVE], with no significant between-group difference in network activity measured during movement (F(2,33) = 1.116, p = 0.340). Of note, a significant group difference in ΔNMRP [MOVE − NO MOVE] values (Fig. 3C) was noted across the three groups (F(2,33) = 8.699, p = 0.001). Network-level activation responses were reduced in both patient groups relative to the healthy control subjects (early PD: p = 0.033; moderate PD: p = 0.001).
Abnormal NMRP expression in PD. A, NMRP expression measured in the NO MOVE condition was increased in early (EARLY) and in moderately advanced (MOD) PD (F(2,37) = 7.881, p = 0.001; one-way ANOVA), with significant elevations in both patient groups relative to normal (NL) control subjects. No difference in network activity was noted between the two PD groups (p = 0.821). B, Group differences in NMRP expression measured during motor execution (MOVE) were not significant (F(2,33) = 1.116, p = 0.340). C, A significant difference across groups was noted for ΔNMRP, the movement-specific activation response, defined as the difference between NMRP [MOVE] and NMRP [NO MOVE] for each subject (F(2,33) = 8.699, p = 0.001). ΔNMRP was reduced relative to control values for both PD groups (EARLY: p = 0.033, MOD: p = 0.001; post hoc Bonferroni test). NMRP expression values was z-scored with respect to the NO MOVE scans (n = 39) of the original derivation cohort (see Materials and Methods). *p < 0.05 (post hoc Bonferroni test). **p < 0.01 (post hoc Bonferroni test).
To address the possibility that the observed elevations in NMRP [NO MOVE] activity were caused by intercurrent tremor, we compared network activity in 10 atremulous PD subjects, with corresponding values from 10 age- and severity-matched patients with observable tremor and from the eight healthy control subjects (see Materials and Methods). This analysis (Fig. 4A) revealed a significant group difference in NMRP [NO MOVE] expression (F(2,25) = 4.793, p = 0.017), with increases of similar degree in tremulous and atremulous PD patients relative to healthy subjects (Tremor(+): p = 0.034, Tremor(−): p = 0.035). There was no significant difference between patients with and without tremor (p = 1.0). Group differences in NMRP expression measured during motor execution (Fig. 4B) were not significant (F(2,25) = 1.421, p = 0.260). However, a significant difference across the three groups (Fig. 4C) was present for the ΔNMRP [MOVE − NO MOVE] activation response (F(2,25) = 7.533, p = 0.003), with abnormally reduced activation responses in the tremulous PD subgroup (p = 0.002) and a trend-level reduction (p = 0.1) in the atremulous subgroup. No significant difference in ΔNMRP values was present, however, across the two patient subgroups (p = 0.293).
NMRP expression PD: absence of tremor effect. NMRP expression values computed in atremulous (Tremor(−)) PD patients were compared with corresponding values from their age- and severity-matched counterparts with observable tremor (Tremor(+)) as well as measurements in normal (NL) control subjects. A, Pattern expression in the nonmovement condition (NMRP [NO MOVE]) was increased to similar levels in tremulous and atremulous patients, with significant elevations in both PD subgroups relative to healthy control values (p < 0.04). B, Group differences in NMRP expression measured during motor execution (MOVE) were not significant (F(2,25) = 1.421, p = 0.260). C, A significant difference across groups was noted for ΔNMRP [MOVE − NO MOVE], the network-level motor activation response (F(2,25) = 7.533, p = 0.003). ΔNMRP was abnormally reduced in the tremulous PD subgroup (p = 0.002; post hoc Bonferroni test). A trend-level reduction (p = 0.1) in ΔNMRP was evident in the atremulous PD subgroup relative to healthy control subjects. The two PD subgroups, however, did not differ with respect to this measure (p = 0.293). *p < 0.05 (post hoc Bonferroni test). **p < 0.01 (post hoc Bonferroni test).
Longitudinal changes in NMRP expression
Analysis of the longitudinal NMRP [NO MOVE] data (Fig. 5A) revealed a significant increase in network activity over time (t(17) = 3.49, p = 0.003). Pattern expression values were elevated at each time point relative to healthy control values (TP1: t(18) = 2.11, p = 0.049; TP2: t(18) = 3.18, p < 0.005; TP3: t(12) = 2.79, p = 0.016; Student's t test). Longitudinal increases in NMRP [MOVE] values (Fig. 5B) were also present (t(17) = 4.25, p < 0.001), although pattern expression became significantly abnormal only at the last time point (TP1: t(18) = 0.735, p = 0.472; TP2: t(18) = 1.610, p = 0.125; TP3: t(12) = 3.68, p = 0.003). Longitudinal changes in NMRP [MOVE] and [NO MOVE] values were found to be intercorrelated (r = 0.70, p < 0.001; Bland-Altman correlation), and ΔNMRP [MOVE − NO MOVE] values (Fig. 5C) did not change significantly over the same time interval (t(17) = 0.002, p = 0.99).
Longitudinal changes in NMRP expression. A, NMRP expression measured in nonmovement (NO MOVE) scans was elevated relative to healthy control (NL) values at each of the longitudinal time points (TP1: t(18) = 2.108, p = 0.049; TP2: t(18) = 3.180, p = 0.005; TP3: t(12) = 2.786, p = 0.016). Within subjects, these values significantly increased with disease progression (t(17) = 3.487, p = 0.0028; general linear model). B, In the MOVE scans, by contrast, the NMRP elevations reached significance with respect to control values only at the final time point (TP1: t(18) = 0.735, p = 0.472; TP2: t(18) = 1.610, p = 0.125; TP3: t(12) = 3.676, p = 0.003). C, The ΔNMRP activation responses were significantly reduced at the first two time points (TP1: t(18) = −2.272, p = 0.036; TP2: t(18) = −2.274, p = 0.035; TP3: t(12) = −1.713, p = 0.112). NMRP expression in MOVE and NO MOVE increased in parallel over time (r = 0.70, p < 0.001; Bland-Altman correlation). Accordingly, the longitudinal change in ΔNMRP values was not significant (t(17) = 0.002, p = 0.998). *p < 0.05, compared with NL (Student's t test). **p < 0.01, compared with NL (Student's t test).
Changes in NMRP expression with treatment
NMRP expression in NO MOVE scans (Fig. 6) was reduced by STN stimulation (t(13) = 3.57, p = 0.003), although network activity during treatment remained elevated relative to healthy control values (OFF: t(20) = 3.71, p = 0.001; DBS: t(20) = 2.64, p = 0.016; Student's t test). The effects of STN stimulation on NMRP [MOVE] and ΔNMRP values (Fig. 6B,C) were not significant (p > 0.12). Despite producing a similar degree of motor benefit, levodopa treatment (Fig. 7A) did not alter NMRP [NO MOVE] expression (t(13) = 1.22, p = 0.25); abnormal elevations in these values were present both at baseline and with treatment (OFF: t(20) = 3.13, p = 0.005; LD: t(20) = 2.95, p = 0.008; Student's t test). Levodopa administration had no effect (t(13) = −0.536, p = 0.60) on NMRP [MOVE] values (Fig. 7B), although a marginal increase in ΔNMRP (Fig. 7C) was observed with this intervention (t(13) = 2.12, p = 0.054).
NMRP expression: effects of STN stimulation. A, NMRP expression measured in the NO MOVE scans of PD patients treated with bilateral STN stimulation was abnormally elevated at baseline and during treatment (OFF: t(20) = 3.710, p = 0.001; ON: t(20) = 2.640, p = 0.016). Therapeutic stimulation was associated with a significant reduction in NMRP [NO MOVE] values (t(13) = 3.571, p = 0.003). B, By contrast, NMRP expression measured during MOVE did not differ from control (NL) values in either treatment condition (OFF: t(14) = 1.667, p = 0.118; ON: t(14) = 0.634, p = 0.536). Indeed, STN stimulation did not significantly alter pattern expression in patients scanned while performing the MOVE task (t(7) = 1.799, p = 0.115). C, In this cohort, ΔNMRP was abnormally reduced at baseline (OFF: t(14) = −2.894, p = 0.012) and a trend-level decrease was present during stimulation (ON: t(14) = −2.138, p = 0.051). The treatment-mediated change in this measure was, however, not significant (t(7) = −0.849, p = 0.424). *p < 0.05, compared with NL (Student's t test). **p < 0.01, compared with NL (Student's t test).
NMRP expression: effects of levodopa infusion. A, NMRP expression measured in NO MOVE scans was abnormally elevated in both treatment conditions (OFF: t(20) = 3.133, p = 0.005; ON: t(20) = 2.951, p = 0.008). No significant effect of levodopa treatment was observed on this measure (t(13) = 1.217, p = 0.245). B, NMRP expression in MOVE scans did not differ from control (NL) values in either treatment condition (OFF: t(20) = 1.295, p = 0.210; ON: t(20) = 1.696, p = 0.105). Levodopa administration had no significant effect on this measure (t(13) = −0.536, p = 0.601). C, ΔNMRP was reduced relative to control values at baseline (OFF: t(20) = −3.654, p = 0.002) but not during treatment (ON: t(20) = −1.654, p = 0.114). Changes in ΔNMRP with treatment were not significant (t(13) = −2.116, p = 0.054). **p < 0.01, compared with NL (Student's t test).
Correlates of abnormal NMRP expression in PD
NMRP [NO MOVE] values measured in the baseline off-medication scans of the PD patients did not correlate with PDRP expression computed in metabolic images from the same subjects (r = 0.248, p = 0.203). Analysis of the longitudinal data (Fig. 8A) revealed that the changes in NMRP [NO MOVE] activity seen with disease progression were accompanied by concurrent increases in PDRP expression. Nonetheless, correlations between the two network measures were not significant at any of the longitudinal time points (TP1: r = 0.367, p = 0.241; TP2: r = 0.520, p = 0.083; TP3: r = 0.081, p = 0.878). Likewise, the changes in these variables observed over time were not significantly correlated (r = 0.235, p = 0.332; Bland-Altman correlation). The progressive increases in UPDRS motor ratings seen in the longitudinal cohort correlated with concurrent changes in the activity of both networks (NMRP [NO MOVE]: t(17) = 2.92, p = 0.01; PDRP: t(17) = 2.13, p = 0.04; GLM). Indeed, clinical deterioration in these subjects was better predicted using the two network measures together in a single model than by the PDRP alone (F(1,16) = 6.711, p = 0.019). The bivariate model provided borderline improvement in clinical prediction compared with the NMRP alone (F(1,16) = 3.174, p = 0.094).
Clinical correlates of changing network activity. A, In the progression cohort (n = 12), the observed longitudinal increases in UPDRS motor ratings correlated with concurrent changes in NMRP [NO MOVE] (t(17) = 2.92, p = 0.01) and PDRP (t(17) = 2.13, p = 0.04) expression. When changes in the two networks were entered together into a single statistical model, the prediction of clinical outcome was more accurate than for each separately. B, C, Similarly, treatment-mediated improvement in motor ratings correlated with NMRP [NO MOVE] values (t(15) = 2.67, p = 0.0176) and with resting PDRP expression (t(15) = 2.33, p = 0.0342). As above, the two networks were more predictive in combination than when applied separately (see Results).
For the treatment cohort, correlations between NMRP [NO MOVE] and PDRP expression were also not significant when assessed at baseline (OFF: r = 0.129, p = 0.634; Pearson's correlation) or in the treated condition (ON: r = 0.160, p = 0.555). The correlation between treatment-mediated changes in the two network measures was, likewise, not significant (r = −0.328, p = 0.215). Treatment-mediated improvement in motor ratings (Fig. 8B,C) correlated with changes in the activity of both networks (NMRP: t(15) = 2.67, p = 0.018; PDRP: t(15) = 2.33, p = 0.034). When subject expression values for the two networks were entered together into a single model, the prediction of clinical outcome was more accurate than when either measure was evaluated individually (p < 0.016).
Discussion
Increased NMRP expression in PD: deficient inhibition of resting network activity
In this study, we used the NMRP, a previously characterized spatial covariance pattern associated with motor activation in normal subjects, as a quantitative probe of network activity in PD patients scanned during movement and in a sensory-matched nonmovement condition. In earlier studies, we have found that normal task-related activation patterns can be altered in disease states, especially if cognitive demand is high (Nakamura et al., 2001; Mentis et al., 2003a; Carbon et al., 2008). Nonetheless, with regard to simple motor execution, PD patients activated the same network as healthy subjects, with minimal violations of the MOVE > NO MOVE ordinal trend. The validation steps that were undertaken to show that the NMRP was indeed deployed by the patients during movement, regardless of disease stage and treatment condition, and that pattern expression was stable in repeated trials, lend credence to the consistent elevations of motor network activity that were observed in the PD group. In this regard, the finding in PD of significant reductions in the ΔNMRP network activation response accords well with fMRI studies showing attenuated activation in the motor cortex of PD patients scanned during movement (Sabatini et al., 2000; Buhmann et al., 2003; Tessa et al., 2010). Surprisingly, however, when we explicitly quantified motor network activity in the MOVE and NO MOVE states, the disease-related changes were found to be greatest in the absence of movement.
The physiological basis for the abnormal NMRP elevations seen in PD is not known. It is tempting, however, to associate this phenomenon to pathologically increased synchronization of cortical-STN and corticostriatal discharges, as has been described consistently in relation to this disorder (Hammond et al., 2007; McCarthy et al., 2011; Moran et al., 2011). In particular, the akinetic-rigid manifestations of PD have been attributed to failure of the normal suppression of β-range oscillatory activity (11–30 Hz) in these pathways before movement onset (Tzagarakis et al., 2010; Zaidel et al., 2010; Crowell et al., 2012). Indeed, substantial evidence exists that antiparkinsonian interventions, such as STN stimulation and levodopa treatment, can improve symptoms of akinesia-rigidity by suppressing excessive synchrony in motor regions (Doyle et al., 2005; Zaidel et al., 2010; Weinberger and Dostrovsky, 2011). Along similar lines, the current findings are also compatible with the notion of deficient GABA-dependent cortical inhibition in PD (Cantello et al., 2002), which can also be ameliorated by treatment (Cunic et al., 2002; Fraix et al., 2008).
NMRP and PDRP: two rest state networks in one disease
Increases in NMRP activity with disease progression and concomitant reductions by treatment are reminiscent of similar changes in the expression of the abnormal PDRP network that have been described (Eidelberg, 2009; Niethammer and Eidelberg, 2012). Nonetheless, the current data indicate that the two networks behave independently with respect to these effects. In this vein, we note that the NMRP and PDRP topographies are spatially distinct (Fig. 1C), with <7.5% correspondence between the regional loadings on the two networks. This minimal correspondence was not significantly different from pseudo-randomly generated volume-pairs with similar level of autocorrelation (p = 0.176). Indeed, a similarly small nonsignificant relationship, accounting for <6% of the total subject score variance, was observed for the changes in the expression of these patterns that occurred over time and with treatment. The distinct nature of these two networks was highlighted by their clinical correlates. Specifically, individual differences in the progression of motor disability in the longitudinal cohort, and in the observed therapeutic response in the treatment cohort, were better explained using NMRP and PDRP expression values in concert, rather than by each network measure separately.
Although the correlations between changes in network activity and clinical motor ratings were similar for NMRP and PDRP, careful analysis revealed these relationships to be both independent and additive. This is particularly evident in the responses of the two networks to treatment. Although both STN DBS and levodopa provided a comparable degree of motor benefit, the former intervention was associated with relatively greater NMRP modulation than the latter. With respect to levodopa, one might consider the possibility that suppression of abnormal baseline NMRP activity during treatment was blunted by the concurrent induction of levodopa-induced dyskinesias. This is rather unlikely, however, given that levodopa-induced dyskinesias were evident during on-state H215O PET imaging in only three of the 14 Group 2B subjects. Nonetheless, during infusion, NMRP [NO MOVE] values measured in these individuals did not influence the result reported for the entire group. Indeed, when the 14 levodopa treatment cases were ordered according to descending on-state NMRP [NO MOVE] expression, the three levodopa-induced dyskinesia subjects ranked at or near the median values of the group. Similarly, measurements of NMRP [MOVE] and ΔNMRP [MOVE − NO MOVE] were also noninfluential for these individuals.
Conversely, levodopa treatment had relatively greater impact on resting PDRP expression than did STN stimulation. Overall, these findings accord with the consistent effects of STN stimulation on corticosubthalamic activity that have been described (e.g., Mure et al., 2012), with concomitant modulation of β-range synchronous oscillatory activity in motor regions. Although dopaminergic treatment can also restore deficient β-suppression in PD patients (Doyle et al., 2005; Weinberger and Dostrovsky, 2011), this intervention also reduces the abnormal synchronization of basal ganglia output activity that characterizes the pathological low-frequency oscillations seen in parkinsonism (Salenius et al., 2002; van Albada et al., 2009). This oscillatory system involves excessive synchronization of neurons in STN, internal globus pallidus, and ventrolateral thalamus, regions with salient contributions to the PDRP metabolic network (Lin et al., 2008; Eidelberg, 2009; Spetsieris and Eidelberg, 2011; Dhawan et al., 2012; Ma et al., 2012). Whether the observed responses of the two networks to different interventions reflect therapeutic effects on different oscillatory systems is unknown.
Interpretative caveats
A number of issues emerge in the interpretation of the current findings. From a practical standpoint, it is conceivable that, in PD patients, the abnormal elevations in NMRP activity seen in the absence of movement are simply a reflection of intercurrent tremor. This is unlikely, however, in that the PD participants were selected for predominant akinetic-rigid symptomatology and had minimal resting tremor during imaging. Moreover, we found that NMRP [NO MOVE] expression is elevated to very similar levels in severity-matched patient subgroups with and without clinically discernible tremor during the imaging session. Thus, it is improbable that tremor substantially impacted the current findings.
Another issue relates to the physiological interpretation of the NO MOVE state. In this condition, the subject visually attended to randomly appearing targets synchronized with auditory pacing tones (Ghilardi et al., 2000). Thus, NO MOVE belongs to a class of eyes-open rest states that require visual fixation. The metabolic scans, by contrast, were conducted in a less constrained eyes-open state in which visual fixation was not required. That said, recent fMRI studies have demonstrated similar patterns of rest state activity whether or not attention to a fixation target was required during the scans (Yan et al., 2009; Van Dijk et al., 2010). In this vein, we have observed nearly identical NMRP expression in PD patients scanned in NO MOVE, which requires visual fixation, and an unconstrained eyes-open nonmovement resting state (NO MOVE: 1.74 ± 0.46; eyes-open: 1.58 ± 0.29; t(12) = 0.266, p = 0.794). Moreover, the eyes-open rest state was not concatenated with movement blocks, as were the NO MOVE trials. Therefore, the abnormal increase in network activity that was observed is unlikely to represent a “hangover” effect of prior motor task performance or cue-induced motor imagery. This possibility is less likely in the current study wherein 10 min were allowed to elapse between trial runs (Nakamura et al., 2001). Indeed, between-trial intervals of this length have been found to be adequate for recovery to pretask baseline after an interceding task block (Barnes et al., 2009).
Although NMRP and PDRP expression was quantified in comparable behavioral conditions, the network measurements themselves were obtained using different scanning techniques. We have previously reported dissociation of network activity values computed in scans of cerebral blood flow and metabolism obtained during levodopa treatment (Hirano et al., 2008), a finding recently substantiated in an experimental rat model (Ohlin et al., 2012). This phenomenon would likely obscure comparison of the levodopa-mediated changes in the expression of the two networks that were observed. Given that the vasomotor effects of levodopa are manifest strictly in the blood flow scans (Hirano et al., 2008), we chose PDRP values computed in the FDG PET scans as the basis for this comparison. Nevertheless, very similar results were obtained when PDRP subject scores were computed in the NO MOVE scans acquired with H215O PET compared with when the network values were computed in the same subjects using their FDG PET scans. This was particularly striking for the progression cohort (main effect of time: F(1,17) = 18.67 and F(1,17) = 24.41 for PDRP [metabolism] and PDRP [blood flow], respectively, p < 0.001 for both measures). Thus, the data suggest that the differences in PDRP and NMRP expression observed in PD patients cannot be attributed to differences in the imaging technique that was used to measure these values.
In conclusion, the findings show that the activity of the NMRP, a large-scale brain network associated with motor functioning in healthy subjects, is elevated in PD patients scanned in the absence of movement. This distinct systems-level abnormality was evident at the earliest clinical stages of the disorder, increasing in magnitude with disease progression. Moreover, the observed changes in NMRP expression were modulated by STN stimulation but not levodopa treatment. Comparisons of the NMRP changes with those occurring concurrently in the expression of the PDRP, a pathological rest state network expressed selectively in PD patients, suggests that the two networks provide independent and additive predictions of clinical outcome. That said, the impact of these network changes on other topographic features of the PD rest state is not known and will be addressed in future studies.
Footnotes
This work was supported by the National Institutes of Health Grant R01 NS 35069 to D.E. and in part by the National Institute of Neurological Disorders and Stroke Grant P50 NS 071675 (Morris K. Udall Center of Excellence for Parkinson's Disease Research at the Feinstein Institute for Medical Research) to D.E. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Neurological Disorders and Stroke or the National Institutes of Health. The sponsor did not play a role in study design, collection, analysis and interpretation of data, writing of the report or in the decision to submit the paper for publication. We thank Dr. Thomas Chaly and Claude Margouleff for technical support and Yoon Young Choi and Toni Fitzpatrick for valuable assistance with manuscript preparation/editing.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. David Eidelberg, Center for Neurosciences, Feinstein Institute for Medical Research, 350 Community Drive, Manhasset, NY 11030. david1{at}nshs.edu
