Abstract
We used a network approach to study the effects of anti-parkinsonian treatment on motor sequence learning in humans. Eight Parkinson's disease (PD) patients with bilateral subthalamic nucleus (STN) deep brain stimulation underwent H215O positron emission tomography (PET) imaging to measure regional cerebral blood flow (rCBF) while they performed kinematically matched sequence learning and movement tasks at baseline and during stimulation. Network analysis revealed a significant learning-related spatial covariance pattern characterized by consistent increases in subject expression during stimulation (p = 0.008, permutation test). The network was associated with increased activity in the lateral cerebellum, dorsal premotor cortex, and parahippocampal gyrus, with covarying reductions in the supplementary motor area (SMA) and orbitofrontal cortex. Stimulation-mediated increases in network activity correlated with concurrent improvement in learning performance (p < 0.02). To determine whether similar changes occurred during dopaminergic pharmacotherapy, we studied the subjects during an intravenous levodopa infusion titrated to achieve a motor response equivalent to stimulation. Despite consistent improvement in motor ratings during infusion, levodopa did not alter learning performance or network activity. Analysis of learning-related rCBF in network regions revealed improvement in baseline abnormalities with STN stimulation but not levodopa. These effects were most pronounced in the SMA. In this region, a consistent rCBF response to stimulation was observed across subjects and trials (p = 0.01), although the levodopa response was not significant. These findings link the cognitive treatment response in PD to changes in the activity of a specific cerebello-premotor cortical network. Selective modulation of overactive SMA–STN projection pathways may underlie the improvement in learning found with stimulation.
Introduction
Parkinson's disease (PD) can affect cognitive functioning, even at its earliest clinical stages (Williams-Gray et al., 2007; Aarsland et al., 2009; Kehagia et al., 2010). Mild cognitive deficits in early PD typically involve visuomotor processing, working memory, and aspects of executive performance (Levin and Katzen, 2005; Aarsland et al., 2010). Motor sequence learning, which encompasses elements of these cognitive functions, has been studied extensively in human subjects (Doyon, 2008; Orban et al., 2010), and mild impairment of task performance has been noted in early-stage PD patients (Nakamura et al., 2001; Ghilardi et al., 2003; Carbon and Eidelberg, 2006; Wilkinson et al., 2009; Carbon et al., 2010b). In this vein, the effects of anti-parkinsonian interventions on cognitive functioning have been found generally to be independent of motor benefit (Cools et al., 2006; Frank and Claus, 2006; Kehagia et al., 2010). Dopamine replacement for the treatment of PD motor symptoms is associated with subtle changes in sequence learning performance assessed at the group mean level (Feigin et al., 2003; Ghilardi et al., 2007). Indeed, recent evidence suggests that individual differences in the cognitive response to levodopa relate to baseline performance capacity (Argyelan et al., 2008; Mattis et al., 2011). This contrasts with the consistent improvement in sequence learning performance observed in medication refractory PD patients treated with internal globus pallidus (GPi) deep brain stimulation (DBS) (Fukuda et al., 2002; Carbon et al., 2003). Although the subthalamic nucleus (STN) DBS procedure that is in routine use has been found to have minimal impact on executive functioning (Parsons et al., 2006; Williams et al., 2011), little is known of its specific effects on motor sequence learning performance. Although levodopa treatment and STN stimulation have both been found to have similar effects on the expression of PD motor networks (Asanuma et al., 2006), it is not known how these interventions affect the neural pathways that mediate cognitive functioning in this disorder.
In the current study, we used a combined psychophysical and functional imaging strategy to evaluate the changes in sequence learning performance and associated network activity that occur during STN stimulation, which were then compared with those observed during levodopa treatment. PD patients and healthy control subjects were scanned with 15O-labeled water (H215O) positron emission tomography (PET) while performing kinematically equivalent motor sequence learning and execution tasks (Nakamura et al., 2001; Ghilardi et al., 2003). The PD subjects performed these tasks in the untreated baseline condition, during STN stimulation, and again while receiving an intravenous levodopa infusion titrated to a comparable degree of clinical improvement. A within-group functional connectivity model (Habeck et al., 2005; Moeller and Habeck, 2006) was applied to identify a specific learning-related network in the imaging data acquired in the baseline and on-stimulation conditions. We compared the degree of network modulation achieved by STN stimulation and levodopa infusion in the same subjects. Additionally, treatment-mediated differences in network activity were correlated with concurrent changes in learning performance recorded during imaging.
Materials and Methods
Subjects
We studied eight right-handed PD subjects (six men and two women; age, 55.8 ± 8.5 years, mean ± SD; Hoehn and Yahr Stage 2.5 ± 0.4) with bilaterally implanted STN DBS electrodes. The demographic features of these subjects are presented in Table 1. The participants exhibited no evidence of dementia or depression on neuropsychological testing conducted in the baseline untreated condition [estimated intelligence quotient (Nelson, 1982; Blair and Spreen, 1989), 121 ± 7; Beck Depression Inventory (Beck, 1987), 5 ± 3; Brief Test of Attention (Schretlen, 1989), 7 ± 1; Hopkins Verbal Learning Test Trial 1–5, 20 ± 6; delayed recall, 4 ± 4 (Brandt, 1991; Benedict et al., 1998); Symbol Digit Test (Smith, 1982), 33 ± 8; controlled oral word association (Benton et al., 1983), 44 ± 15]. There was no evidence of atrophy or other structural brain abnormalities on routine magnetic resonance imaging (MRI).
The subjects were scanned with H215O PET in a single imaging session while they performed two kinematically matched reaching tasks (described below) in each of the following experimental conditions: (1) OFF: a baseline condition in which anti-parkinsonian medications were discontinued for at least 12 h beforehand and with stimulators turned off for at least 1 h; (2) DBS: a treatment condition in which stimulation on the two sides (amplitude, 3.29 ± 0.41 V; pulse width, 67.5 ± 12.5 μs; frequency, 176.7 ± 20.4 Hz) was optimized to achieve a clinically meaningful improvement in motor signs, defined as a reduction on Part III of the Unified Parkinson's Disease Rating Scale (UPDRS) (Fahn et al., 1987) of at least five points or 20% of baseline; and (3) LD: a treatment condition in which the subjects received a steady-state intravenous infusion of levodopa without stimulation (Feigin et al., 2003; Hirano et al., 2008). In this condition, the levodopa infusion rate (mean of 1.13 ± 0.41 mg · kg−1 · h−1) was individually titrated to produce an improvement in motor UPDRS ratings equivalent to that achieved by stimulation without dyskinesia. Eight healthy volunteer subjects (six men and two women; age, 52.7 ± 13.7 years) served as controls for these experiments. These subjects underwent H215O PET while performing the same set of tasks in a single imaging session.
Behavioral testing
Tasks.
In each experimental condition, the subjects were scanned during the performance of a motor sequence learning task (LEARN) and a kinematically equivalent motor execution task (MOVE). The characteristics of these reaching tasks have been described in detail previously (Ghilardi et al., 2003; Argyelan et al., 2008; Ghilardi et al., 2009). Briefly, in performing the two reaching tasks, 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. In LEARN, the eight targets appeared in an unknown, pseudorandom order without repeating elements that was repeated over the 90 s trial block (SEQ). Subjects were informed that a sequence was to be presented. They were instructed to learn the sequence order while reaching for the targets, to anticipate successive targets, and to reach each target in synchrony with the tone. Sequences were different for each trial run. In MOVE, the eight targets appeared in a predictable counterclockwise (CCW) order. To reach the target in synchrony with the tone, subjects had to initiate movement before the target appeared. In addition to the reaching tasks, the subjects were also scanned in a non-movement sensory-matched resting state (REST) in which they passively observed randomly appearing targets (Ghilardi et al., 2000).
To evaluate the effects of STN stimulation on sequence learning and associated brain networks, we designed the imaging experiments so that the order of the trials was randomized and counterbalanced across tasks (LEARN, MOVE) and treatment conditions (OFF, DBS). Because of the long-duration effect of levodopa on motor symptoms (Hauser and Holford, 2002), the infusion studies (LD) followed the DBS and OFF trials. Thus, the possibility existed that the observed differences in learning across experimental conditions, if any, were influenced by treatment order, with practice improving and/or fatigue worsening task performance during levodopa infusion relative to STN stimulation. To address this potential confound, we measured the random reaction time outside the scanner before each LEARN trial. To quantify this measure of perceptual motor processing speed, targets were presented in pseudorandomized, nonrepeating, and unpredictable order, with the subjects required to reach each target as soon as possible (Nakamura et al., 2001). Random reaction time values were measured separately in each of the three experimental conditions. As described below, this measure was used as an empiric criterion for target anticipation during learning performance in each condition (Ghilardi et al., 2003, 2009). The random reaction time measurements were also compared across conditions to determine whether order-dependent changes in motor processing speed (an index of general attentiveness) were evident in the data.
The order of task trials was randomized and counterbalanced within each of the three experimental conditions. Of note, the pace of the movements in LEARN and MOVE, represented by the intertone interval, was empirically adapted to the subject's motor ability in the baseline OFF state. This allowed us to fix the rate of movement across tasks (LEARN, MOVE) and conditions (OFF, DBS, LD). The intertone interval was 1.5 s in five of the eight PD subjects, 2.0 s in two of the subjects, and 1.83 s in one of the subjects. Similarly, we held movement size constant by adjusting the gain to a target extent of 1 cm for all subjects, tasks, and experimental conditions.
Performance measures.
To quantify learning performance in each LEARN run, we computed the number of correct movements initiated below a threshold determined by the random reaction time (Nakamura et al., 2001; Ghilardi et al., 2003, 2009). This threshold was determined outside the scanner and separately quantified for each experimental condition. Movements initiated below this reaction time reflect anticipation and successful retrieval of previously acquired targets. The number of correctly anticipated targets in each cycle were summed defining a composite retrieval index that quantified task performance in each LEARN trial (Ghilardi et al., 2003; Argyelan et al., 2008; Ghilardi et al., 2009). The effects of experimental condition (OFF, DBS, LD) on motor disability (UPDRS motor ratings), perceptual motor processing speed (random reaction time), and learning performance (retrieval index) were assessed in separate one-way repeated-measures analysis of variance (RMANOVA) models; pairwise contrasts were assessed post hoc using Tukey's HSD tests. Additionally, psychophysical performance measures in each condition were compared with control values using two-tailed Student's t tests. Contrasts were considered significant at p < 0.05.
Imaging studies
Positron emission tomography.
All subjects were scanned with H215O PET at baseline (OFF) and in each of the two treatment conditions (DBS, LD) in a single 1 d imaging session (Argyelan et al., 2008). The subjects fasted overnight before imaging; anti-parkinsonian medications were discontinued at least 12 h beforehand. PET imaging was performed using the GE Advance tomograph (GE Healthcare) in 3D mode. In the baseline and treatment conditions, the subjects performed the SEQ and CCW tasks in randomized order using the dominant right arm. Four of the subjects performed the tasks twice in the three experimental conditions. Because of fatigue, the remaining subjects performed the tasks once per condition. Thus, 12 complete scan sets (i.e., pairs of SEQ and CCW scans from individual subjects studied in the OFF and DBS conditions) were available for pattern derivation (see below).
Ethical permission for these studies was obtained from the Institutional Review Board of The Feinstein Institute for Medical Research. Written consent was obtained from each subject after detailed explanation of the procedures.
Network analysis.
Imaging data processing was performed using SPM5 (Institute of Neurology, London, UK). Standard preprocessing (realignment, spatial normalizing, smoothing with an isotropic Gaussian kernel (10 mm) in all directions) was conducted. To identify a specific functional brain network associated with STN stimulation effects on motor sequence learning, we applied a within-subject network modeling approach to the imaging data termed Ordinal Trends/Canonical Variates Analysis (OrT/CVA) (Habeck et al., 2005; Moeller and Habeck, 2006). The OrT/CVA algorithm is designed to identify specific spatial covariance patterns that increase in their expression across experimental conditions on a subject-by-subject basis while the relationships between brain regions remain constant (Carbon et al., 2010a; Mure et al., 2011). OrT/CVA differs from routine univariate analysis in that it requires that pattern expression values exhibit an “ordinal trend”: the property of consistent change across conditions at the individual subject level (i.e., network activity) is required to increase monotonically in all or most of the subjects. As in groupwise spatial covariance analysis (Eidelberg, 2009; Habeck and Stern, 2010; Spetsieris and Eidelberg, 2011), large-scale networks are described in terms of the voxel loadings (“region weights”) on each of the relevant principal component (PC) topographies. Likewise, the expression of a given pattern in each scan is quantified by a specific network activity value (“subject score”), the PC scalar multiplier for the subject in each task and experimental condition. The significance of networks resulting from OrT/CVA is assessed using nonparametric tests. In pattern derivation datasets, permutation tests of the relevant subject scores are used to confirm that the observed monotonic changes in pattern expression across conditions did not occur by chance. The reliability of the voxel loadings comprising the network topography itself is assessed using bootstrap resampling procedures (Efron and Tibshirani, 1994).
In the current study, OrT/CVA was used to characterize the effects of STN stimulation on learning-related cerebral blood flow at the network level. The computational algorithm was applied to the LEARN scans of the PD patients acquired in the OFF and DBS conditions. To minimize confounds stemming from the concurrent effects of anti-parkinsonian treatment on movement (as opposed to learning), we restricted the search for learning-related network topographies to the portion of the subject × voxel space that was independent of motor activation. This was accomplished by defining a subspace relating to motor activation in the current patient and control sample. Specifically, in a previous study (Carbon et al., 2010a), we used OrT/CVA to identify a significant reproducible movement-related spatial covariance pattern in motor execution (MOVE) and rest (REST) H215O PET scan pairs from 18 healthy volunteer subjects. We then quantified the expression of this normal movement-related activation pattern (NMRP) in the MOVE and REST scans acquired in the individual PD patients in each of the three treatment conditions and the healthy control subjects that comprised the current dataset. Analysis of the resulting subject scores disclosed consistent MOVE > REST increases in each group and condition (p < 0.003, paired Student's t tests) with minimal exceptions. Having demonstrated that the individual motor activation responses were present within the NMRP subspace, we conducted the OrT/CVA search for learning-related patterns in the subspace orthogonal to that network, i.e., that which was independent of movement-related subject × voxel variability in the data.
The resulting OrT/CVA networks were considered significant if a consistent monotonic trend (p < 0.01, permutation test, 1000 iterations) was demonstrable for pattern expression values quantified in OFF/DBS scan pairs acquired during learning performance (i.e., increasing subject scores for DBS vs OFF computed in the LEARN scans). The corresponding spatial topographies were displayed at a voxel weight threshold of Z = 2.81 (p = 0.005) with an extent cutoff of 50 voxels. Clusters present at this threshold were considered to represent significant regional “nodes” of the network if they were additionally found to be reliable by bootstrap resampling, as defined by an inverse coefficient of variation >1.64 (p < 0.05, 1000 iterations).
After network identification in OFF/DBS LEARN scan pairs, pattern expression was prospectively quantified in the MOVE scans of the PD patients acquired in these two experimental conditions and in the healthy control subjects. Network activity was additionally measured in the LEARN and MOVE scans of the PD patients acquired during levodopa infusion. As with network identification, prospective computation of pattern expression values in individual subjects/scans was conducted in the space orthogonal to the NMRP motor activation subspace (see above). Network values computed for each task/condition were z-transformed with respect to the corresponding healthy control values so that the normal (NL) group had a mean of 0 and an SD of 1. All OrT/CVA computations were performed using software freely available at http://www.nitrc.org/projects/gcva_pca.
Data analysis.
Network scores computed in the baseline LEARN and MOVE scans of the PD patients were compared with those from the healthy control subjects. Group differences were evaluated using 2 × 2 RMANOVA with group (PD, NL) and task (LEARN, MOVE) as the between- and within-subject factors. For each subject, we additionally computed a learning-specific measure of network activity. This was defined as the difference in pattern expression computed in the LEARN and MOVE scans from that individual (i.e., the LEARN − MOVE differences in subject scores). The resulting value was correlated with the corresponding learning performance measure from the same subject. Correlation analysis was performed separately in each group by calculating Pearson's correlation coefficients. Network activity and task performance values from individual subjects with multiple LEARN runs were averaged for this analysis.
Similar statistical analyses were performed to quantify subject differences in network activity occurring during treatment and to correlate these measures with concurrent changes in task performance. Stimulation-mediated changes in pattern expression measured in LEARN and MOVE scans were assessed using paired Student's t tests. Likewise, network values measured during stimulation (DBS) were compared with the corresponding normal values using Student's t tests. Analogously, network values computed in LEARN and MOVE scans acquired during stimulation were compared with corresponding values from treatment scans acquired during levodopa infusion (LD). Network values computed during levodopa treatment were also compared with those from the healthy control subjects. For each intervention, we also correlated learning-specific network modulation during treatment with concurrent changes in task performance. As above, this was done by computing Pearson's product moment correlation coefficients.
Nodal analysis
After characterizing the effects of treatment at the network level, we evaluated the changes that occurred concurrently at each of the network nodes. This was done by measuring regional cerebral blood flow (rCBF) post hoc in volumes of interest (VOIs) corresponding to the significant clusters that defined the network topography (see above). Each node was represented by a spherical VOI (radius of 4 mm) centered at the peak voxel of the corresponding cluster. Only nodes with significant changes in rCBF with STN stimulation (p < 0.05, paired Student's t test) were selected for additional analysis. In these VOIs, rCBF values measured in the OFF, DBS, and LD conditions in the PD patients were compared with corresponding values from the healthy control subjects using Student's t tests. Differences in rCBF values measured in the DBS and LD conditions were assessed separately for each node using paired Student's t tests. The results were considered significant at p < 0.05.
Finally, we used the rCBF data from each of the VOIs to explore the relative influence of the individual nodes on overall network behavior during treatment. Specifically, for each node, we determined the degree to which the effects of STN stimulation on rCBF (an index of local synaptic activity) were consistent across subjects and trials. This was done separately for each region by computing the number of trials (in proportion to the total) that violated the ordinal trend observed for the network as a whole. We hypothesized that network nodes exhibiting uniform rCBF responses to stimulation (i.e., with significantly fewer individual trials that violated the network-level trend) had closer functional connection to the DBS target than those that exhibited less consistent (i.e., more random) responses at the individual trial level. For the positive nodes (i.e., those characterized by network-related increases in activity during stimulation), trials exhibiting concurrent treatment-mediated declines in rCBF were considered to be at variance with the network trend and were classified as exceptions. Analogously, for the negative nodes (i.e., those with network-related decreases in activity during stimulation), trials with increasing rCBF during treatment were considered inconsistent with the network-level response and were counted as exceptions. Operationally, rCBF responses to stimulation at a given node were viewed as consistent if the frequency of exceptions in the individual trial data was lower than expected by chance. This possibility was evaluated using the binomial test and was considered significant for p < 0.05 after Bonferroni's correction for multiple regional comparisons.
All statistical analyses were conducted using SPSS software, and the results were considered significant at p < 0.05.
Results
Behavioral measures
Mean ± SD values for the clinical motor ratings, random reaction times, and learning performance measures for the PD patients in the OFF, DBS, and LD conditions, and for the healthy control subjects are presented in Table 1. UPDRS motor ratings (Fig. 1A) differed significantly across the three experimental conditions (F(2,14) = 24.6, p < 0.001, one-way RMANOVA), with improvement in clinical motor ratings during anti-parkinsonian treatment with either STN stimulation or levodopa infusion (p < 0.001 relative to baseline, post hoc tests). Between-condition differences (F(2,13) = 8.2, p < 0.005) were also observed for random reaction time values (Fig. 1B), with significant reductions in this measure in both treatment conditions (p < 0.05 relative to baseline, post hoc tests). Sequence learning performance (Fig. 1C) also differed across conditions (F(2,14) = 7.1, p < 0.01). However, in contrast to the motor ratings and reaction time measures, sequence learning differed for the two interventions, with significantly better performance during STN stimulation relative to both levodopa and the untreated baseline condition (p = 0.01, post hoc tests). Of note, there was no correlation between stimulation-related improvement in learning performance and the DBS parameters used during imaging (p > 0.34). In contrast to STN stimulation, learning performance measured during levodopa infusion did not differ from baseline values (p = 0.99). Indeed, target retrieval was abnormally low both at baseline and during levodopa treatment (p < 0.02 compared with healthy control values, Student's t tests), whereas learning performance measured during STN stimulation did not differ significantly from normal (p = 0.45).
Network characterization
A significant functional network was identified by OrT/CVA of on- and off-stimulation scan pairs acquired in the PD patients during sequence learning. This spatial covariance pattern (Table 2) was characterized by stimulation-related increases (Fig. 2A, red) in the right lateral cerebellum (lobule VII) and parahippocampal gyrus (BA 36) and in the left dorsal premotor (BA 6) and inferior parietal (BA 40) regions. These changes were associated with stimulation-related reductions medially (Fig. 2A, blue) in the orbitofrontal cortex (BA 11) and caudal supplementary motor area (SMA).
Network activity values (Fig. 2B) measured in LEARN scans acquired at baseline (OFF) and during stimulation (DBS) exhibited a significant monotonic trend across conditions (p = 0.008, permutation test) with increasing pattern expression in 11 of 12 trials. Pattern expression was computed prospectively in the LEARN and MOVE scans from the healthy control subjects. In this group, network values computed in LEARN (Fig. 2C, left) were found to be higher than in MOVE (p = 0.02, paired Student's t test). In contrast, PD patients scanned at baseline (Fig. 2C, right) exhibited reduced network activity in LEARN relative to MOVE (p < 0.005). Indeed, learning-specific network responses (i.e., difference between LEARN and MOVE subject scores) for the PD patients in the untreated condition and the healthy control subjects were significantly different between the two groups (F(1,16) = 11.2, p < 0.005; group × task interaction effects). Network activity values computed in the resting condition in the baseline PD scans were also lower than the corresponding values from the healthy control subjects (p < 0.05). In both the control and PD groups (Fig. 2D,E), network activity correlated with concurrent measurements of task performance obtained during each LEARN trial (NL, r = 0.84, p < 0.01; PD (baseline), r = 0.74, p = 0.03, Pearson's product moment correlations).
Changes in network activity during treatment
STN stimulation was associated with significant increases in network activity (Fig. 3A,B) measured during the sequence learning and motor execution tasks (p < 0.002, paired Student's t tests for separate DBS vs OFF comparisons of network values computed in LEARN and MOVE scans). Nonetheless, stimulation-mediated changes in network activity were greater during learning than during movement (task × condition interaction: F(1,20) = 4.8, p = 0.04, 2 × 2 RMANOVA). Indeed, individual differences in the degree of network modulation that occurred with stimulation [(LEARN − MOVE)DBS − (LEARN − MOVE)OFF] (Fig. 3C) correlated with concurrent changes in learning performance observed in the same subjects (r = 0.78, p = 0.02, Pearson's correlation). The stimulation-mediated network changes did not, however, correlate with differences in the DBS parameters used during imaging (p > 0.26).
Despite comparable degrees of improvement in clinical motor ratings, levodopa infusion and STN stimulation differed in their impact on learning-related network activity. Network activity measured during learning (Fig. 3A) was greater for STN stimulation than levodopa treatment (p = 0.03, paired Student's t test for DBS vs LD comparison of network values measured in LEARN scans). Significant between-intervention differences in network modulation (Fig. 3B) were not evident in the non-learning tasks (MOVE, p = 0.25; REST, p = 0.49). A correlation of marginal significance (r = 0.68, p = 0.06) was noted between levodopa-mediated differences in network activity measured during learning and concurrent changes in task performance (Fig. 3D).
Nodal analysis
We next examined the regional correlates of the observed network changes. First, we examined rCBF values in regions with increases in learning-related activity during stimulation (Fig. 2A, red–yellow; Table 2). Significant effects of treatment on LEARN scan rCBF (p < 0.05, paired Student's t tests) were evident in the left cerebellar hemisphere (lobule VII) and in the right premotor cortex (Fig. 4, top panels). Baseline rCBF values measured in both these areas were lower than normal in the PD patients, regardless of task (cerebellum: LEARN, p < 0.01; MOVE, p < 0.05; REST, p < 0.05; premotor cortex: LEARN, p < 0.01; MOVE, p < 0.05; REST, p < 0.01, Student's t tests). In both regions, rCBF values were increased by stimulation during learning (DBS > OFF: p < 0.05, paired Student's t tests for LEARN scan comparisons) but not during motor execution (cerebellum, p = 0.054; premotor cortex, p = 0.18 for MOVE scan comparisons) or rest (p > 0.55 for REST scan comparisons). rCBF values in these regions did not differ significantly across the two interventions (LEARN, p > 0.09; MOVE, p > 0.37; REST, p > 0.56, paired Student's t tests for comparisons of DBS and LD).
Similarly, we examined the rCBF data from regions with reductions in learning-related activity during stimulation (Fig. 2A, blue–green; Table 2). Significant treatment effects on rCBF values measured during learning were evident in the orbitofrontal cortex (BA 11) and the SMA (Fig. 5A,B). In both regions, rCBF values measured in LEARN scans were abnormally elevated in the baseline condition (orbitofrontal cortex, p < 0.05; SMA, p < 0.01, Student's t tests). In the orbitofrontal cortex, baseline rCBF measured during movement and rest did not differ from control values (p > 0.65). In this region, stimulation was associated with reduced rCBF during learning (p < 0.04, paired Student's t test). Significant changes in rCBF during stimulation were not present in this region when measured during movement or rest (p > 0.67). Likewise, orbitofrontal rCBF values measured during each of the three tasks did not differ for the two interventions (p > 0.51, paired Student's t tests for DBS vs LD comparisons).
In contrast, in the SMA, abnormal elevations in baseline rCBF values were present during learning and also during the performance of the non-learning movement and rest tasks (MOVE, p < 0.05; REST, p < 0.001, Student's t tests for comparisons with normal values). STN stimulation resulted in significant reductions in SMA values measured during LEARN (p < 0.02) and REST (p < 0.01); a trend-level reduction (p = 0.1) was discerned during MOVE. Indeed, in this region, significant between-intervention differences were evident during learning and rest (LEARN, p < 0.02; REST, p < 0.003, paired Student's t tests for comparisons of DBS and LD), with lower rCBF during STN stimulation compared with levodopa. A borderline difference between interventions (p = 0.07) was noted in SMA for rCBF values measured during movement. Thus, significant baseline rCBF elevations were evident in the SMA, independent of task. These abnormalities were modulated toward normal by anti-parkinsonian treatment with STN stimulation but not levodopa.
Last, we assessed the consistency of the stimulation effects on rCBF values measured at each node with the concurrent changes observed at the network level (see Materials and Methods) (Fig. 6). Of the four nodal regions with significant stimulation-mediated rCBF changes, only the SMA exhibited local responses that consistently paralleled the network trend. Specifically, rCBF values in this region were remarkable for the presence of only a single violation in the 12 stimulation trials that were performed during learning (p = 0.01, binomial test, corrected for multiple comparisons). This contrasted with the other nodes (the cerebellum, the premotor cortex, and the orbitofrontal cortex) in which four to five exceptions were noted per region, which did not differ from the number expected by chance (p > 0.48).
Despite the presence of significant rCBF changes during treatment in nodal regions, these regional measures did not correlate with concurrent changes in learning performance (p > 0.16).
Discussion
Sequence learning network associated with STN stimulation
Electrophysiological and functional imaging studies have consistently linked the learning of motor sequences to the activation of cortico-striato-pallido-thalamo-cortical and cerebello-thalamo-cortical neural pathways (Doyon et al., 2009; Orban et al., 2010). Indeed, we have found previously that sequence learning performance was associated with the activation of discrete brain networks in healthy subjects and in unmedicated PD patients (Nakamura et al., 2001; Carbon et al., 2003; Carbon and Eidelberg, 2006). Moreover, changes in learning network activity during the treatment of PD motor symptoms with GPi stimulation or levodopa administration correlated with concurrent recordings of task performance during imaging (Carbon and Eidelberg, 2002; Carbon et al., 2003). That said, these previous studies did not address the possibility of distinct treatment-specific network effects on learning pathways.
In the present study, we used OrT/CVA, a multivariate algorithm designed to detect and quantify functional brain networks with monotonically changing subject expression across experimental conditions on an individual case basis (Habeck et al., 2005). This analytical approach differed from that used in our previous studies of learning-related activation patterns in PD patients (Nakamura et al., 2001; Carbon et al., 2003). We initially applied spatial covariance analysis to the subtracted images (LEARN − MOVE) from healthy subjects (or from untreated PD patients) to identify patterns correlating with learning performance. The resulting learning-related patterns were then used prospectively to quantify the effects of treatment on network activity in individual subjects. This strategy, however, relied on the premise that therapeutic interventions such as DBS act by modulating the activity of normal brain networks (Carbon and Eidelberg, 2002) as opposed to inducing a distinctive set of network changes related specifically to treatment. To address the latter possibility, OrT/CVA was used in the current study to indentify a specific functional network characterized by stimulation-mediated changes. A number of steps were taken to ensure that the resulting network was learning related and not a manifestation of nonspecific stimulation effects on rCBF. First, the network identification procedure (as well as the prospective single scan computations) was conducted in a subspace of the mean centered subject × voxel covariance matrix that was independent of (orthogonal to) that originally specified in relation to motor activation in normal individuals (Carbon et al., 2010a) and that was demonstrated in the current cohort of PD and healthy control subjects. Furthermore, network activity values computed in the baseline scans of the PD patients and in the healthy control subjects correlated significantly with concurrent measurements of learning performance. Last, network values computed in the on-stimulation condition did not correlate with individual differences in the DBS tuning parameters, nor did the latter influence the slope of the observed network-performance correlation line. Thus, although the network was identified using scan pairs obtained off- and on-stimulation, the relevant topography was specified by applying OrT/CVA to images acquired during learning and therefore related preferentially to this cognitive process.
Anatomical and functional basis for the learning-related network
Many of the areas identified as nodes of the learning-related network identified in this study have been described previously in univariate analyses of PET and fMRI activation data from cohorts of healthy subjects performing a variety of sequence learning tasks (for review, see Carbon and Eidelberg, 2006; Doyon et al., 2009; Wilkinson et al., 2009). Typically, the cerebellum, cingulate, premotor cortex, and inferior parietal area are activated in the early phases of learning. In contrast, the striatum, SMA, prefrontal cortex, and precuneus are activated during the later stages of this cognitive process. In the current study, the patients were scanned comparatively early in the learning process, which accords with the areas of increased activity featured by the network. Overall, these nodal regions exhibited abnormally low neural activity in the baseline condition, regardless of the task performed. Nevertheless, the normalization of rCBF that was observed in these areas during stimulation was specific for sequence learning and paralleled the cognitive response to treatment.
At two network nodes, STN stimulation influenced neural activity by lowering abnormal baseline elevations of neural activity. In the orbitofrontal cortex, neural activity was significantly elevated in the baseline scans of the PD patients acquired during learning but not during motor execution or rest. Accordingly, STN stimulation lowered rCBF in this region when measured during learning but effected no change in movement and rest scans from the same subjects. These findings accord with the notion that the orbitofrontal cortex can be deployed during sequence learning under pathological conditions (Feigin et al., 2006). In this vein, the elevation of learning-related neural activity observed in this region in untreated PD patients may be compensatory. Indeed, normalization of rCBF in this region by either STN stimulation or levodopa is consistent with this interpretation.
In addition to the orbitofrontal cortex, abnormal baseline rCBF elevations in the PD patients were evident in the SMA proper. Given the prominent contribution of increased SMA activity to the abnormal resting metabolic network characteristic of PD (Eidelberg, 2009; Tang et al., 2010), it is not surprising that the baseline rCBF increases observed in this region were task independent. Moreover, STN stimulation improved these regional abnormalities in a similar task-independent manner, with significant reductions in rCBF during learning and rest and a trend in that direction during movement.
The dominant role of the SMA in mediating the changes in network activity seen during STN stimulation is underscored by the uniform rCBF treatment responses that characterized this brain region. These findings are consistent with previous data from human subjects and experimental animal models. A recent resting-state fMRI study demonstrated a significant increase in the functional connectivity of pathways linking the STN and cortical motor regions (including SMA) in PD patients (Baudrexel et al., 2011). Indeed, abnormal elevations in the metabolic activity of the STN and motor cortex (M1/SMA) have been observed during PD progression (Huang et al., 2007), in parallel with increases in the PD-related metabolic network mentioned above (Lin et al., 2008; Eidelberg, 2009; Tang et al., 2010). These changes are likely attributable to overactivity of the “hyperdirect” M1/SMA–STN projection system demonstrated in animal models (Dejean et al., 2008; Nachev et al., 2008). Oversynchronization of high-frequency oscillatory activity in this loop may account for the abnormally high cerebral blood flow and metabolic activity seen in these regions in untreated PD patients. In this context, STN stimulation lowered these values toward normal, with concomitant improvement in learning performance as well as generalized motor functioning. The current data are also compatible with the notion that high-frequency STN DBS exerts its therapeutic effect through antidromic modulation of cortical–STN axonal projections (Gradinaru et al., 2009). In the rodent, localized modulation of motor cortical activity may be needed for the therapeutic effects of STN stimulation to occur. Nonetheless, in human PD patients, clinical outcome is likely to be contingent on changes in activity of the network as a whole. Indeed, concurrent normalization of learning-related rCBF abnormalities at the other network nodes (e.g., cerebellum, premotor cortex, and inferior parietal area) argues in favor of a combination of stimulation-mediated functional effects at multiple sites.
One may speculate as to why analogous therapeutic effects on learning performance and associated network activity were not observed when levodopa was administered to the same subjects. Prospective quantification of pattern expression in scans acquired during levodopa infusion suggested that the learning-related network was indeed deployed in this treatment condition. Nonetheless, on average, network activity was significantly lower with levodopa than the motorically equivalent stimulation condition. In this vein, subject scores for the network measured in LEARN scans acquired during levodopa infusion were on average 1 SD below the normal mean, whereas corresponding values measured during STN stimulation were ∼1 SD above normal. Of note, there was no correlation between clinical descriptors, such as subject age, symptom duration, and motor severity ratings, and treatment-mediated changes in learning performance and associated network activity.
Analogous trends were also seen at the regional level: STN stimulation was relatively more effective than levodopa in normalizing local neural activity in network areas. This difference was most pronounced at the critical SMA node. It is likely that modulation of abnormal baseline activity in this region by STN stimulation—and the accompanying improvement in learning performance—stems from the selectivity of this intervention for the axons of the hyperdirect M1/SMA–STN pathway. That functional changes in this pathway were seen in response to cognitive as well as motor demands suggests a potential role for the SMA as a discrete treatment target. Indeed, therapeutic interventions directed specifically at this brain region may benefit PD patients with medication-refractory symptoms, including commonly encountered cognitive deficits in sequence learning and skill acquisition (Ghilardi et al., 2003; Carbon and Eidelberg, 2006; Doyon, 2008). Importantly, the findings also lend support to the idea of conceptualizing the therapeutic target as the entire functional network as a single entity, as opposed to each of its regional elements in isolation.
Footnotes
This work was supported by National Institutes of Health/National Institute of Neurological Disorders and Stroke Grants RO1 NS 35069 and P50 NS 071675 (D.E.). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NINDS or the NIH. We thank Dr. Thomas Chaly for radiochemistry support and Claude Margouleff for technical assistance. Special thanks to Toni Fitzpatrick for valuable help with copyediting.
- Correspondence should be addressed to Dr. David Eidelberg, Center for Neurosciences, The Feinstein Institute for Medical Research, 350 Community Drive, Manhasset, NY 11030. david1{at}nshs.edu