Original Research ARTICLE
Neural substrates underlying motor skill learning in chronic hemiparetic stroke patients
- 1Institute of Neuroscience, Université Catholique de Louvain, Brussels, Belgium
- 2Neurology Department, CHU Dinant-Godinne UCL Namur, Université Catholique de Louvain, Yvoir, Belgium
- 3Imagilys, Brussels, Belgium
- 4Faculty of Electronics and Information Technology, Institute of Radioelectronics, Warsaw University of Technology, Warsaw, Poland
- 5Neurology Department, Site Saint-Joseph, CHC, Liège, Belgium
- 6Neurology Department, Clinique Saint-Pierre, Ottignies, Belgium
- 7Service de Neurologie, Unité Neuro-Vasculaire, Cliniques Universitaires Saint Luc UCL, Université Catholique de Louvain, Brussels, Belgium
- 8Scientific Support Unit, CHU Dinant-Godinne UCL Namur, Université Catholique de Louvain, Yvoir, Belgium
- 9Louvain Bionics, Université Catholique de Louvain, Louvain-la-Neuve, Belgium
Motor skill learning is critical in post-stroke motor recovery, but little is known about its underlying neural substrates. Recently, using a new visuomotor skill learning paradigm involving a speed/accuracy trade-off in healthy individuals we identified three subpopulations based on their behavioral trajectories: fitters (in whom improvement in speed or accuracy coincided with deterioration in the other parameter), shifters (in whom speed and/or accuracy improved without degradation of the other parameter), and non-learners. We aimed to identify the neural substrates underlying the first stages of motor skill learning in chronic hemiparetic stroke patients and to determine whether specific neural substrates were recruited in shifters versus fitters. During functional magnetic resonance imaging (fMRI), 23 patients learned the visuomotor skill with their paretic upper limb. In the whole-group analysis, correlation between activation and motor skill learning was restricted to the dorsal prefrontal cortex of the damaged hemisphere (DLPFCdamh: r = −0.82) and the dorsal premotor cortex (PMddamh: r = 0.70); the correlations was much lesser (−0.16 < r > 0.25) in the other regions of interest. In a subgroup analysis, significant activation was restricted to bilateral posterior parietal cortices of the fitters and did not correlate with motor skill learning. Conversely, in shifters significant activation occurred in the primary sensorimotor cortexdamh and supplementary motor areadamh and in bilateral PMd where activation changes correlated significantly with motor skill learning (r = 0.91). Finally, resting-state activity acquired before learning showed a higher functional connectivity in the salience network of shifters compared with fitters (qFDR < 0.05). These data suggest a neuroplastic compensatory reorganization of brain activity underlying the first stages of motor skill learning with the paretic upper limb in chronic hemiparetic stroke patients, with a key role of bilateral PMd.
Stroke is a devastating disorder that causes life-long upper limb hemiparesis in 30–70% of survivors (Lai et al., 2002; Kwakkel et al., 2003). The biochemical mechanisms triggered by acute stroke (e.g., edema resolution, inflammation, up- and down-regulation of neurotransmitters) play a prominent role in early recovery (Kreisel et al., 2006; Carey and Seitz, 2007). Beyond these biochemical cascades, recovery of motor function also relies on neuroplastic reconfiguration of the cortical motor network and its descending projections, which support transfer of impaired functions toward undamaged areas of the brain (Feydy et al., 2002; Johansen-Berg et al., 2002; Lotze et al., 2006; Lindenberg et al., 2010; Schulz et al., 2012). Although this neuroplastic reorganization may reflect a simple re-routing of information flow through pre-existing, undamaged pathways, stroke patients must learn how to recruit these neuronal resources. To some extent, recovering from hemiparesis might be conceptualized as a particular form of motor skill learning, in other words, learning to use the reconfigured motor network to optimize planning, execution and movement control of the paretic upper limb. Indeed, the idea that motor skill learning plays a central role in post-stroke motor recovery is becoming a major focus in neurorehabilitation (Matthews et al., 2004; Krakauer, 2006; Dipietro et al., 2012; Kitago and Krakauer, 2013).
The neural substrates of motor skill learning are relatively well elucidated in healthy individuals. Functional magnetic resonance imaging (fMRI) studies demonstrated that motor skill learning relies on a network encompassing the primary motor cortex (M1), supplementary motor area (SMA), premotor cortex (PM), dorsolateral prefrontal cortex (DLPFC), cerebellum and basal ganglia (Ghilardi et al., 2000; Halsband and Lange, 2006; Debas et al., 2010; Hardwick et al., 2013). Recently, the definition of motor skill learning has been refined to a training-induced acquisition and improvement of motor performance (i.e., skills), persisting over time and characterized by a shift of the speed/accuracy trade-off (SAT), automatisation and reduction of performance variability (Reis et al., 2009; Dayan and Cohen, 2011; Krakauer and Mazzoni, 2011). Using a visuomotor skill learning paradigm involving a SAT, we demonstrated that different that different behavioral trajectories can be observed in healthy individuals during the first stages of motor skill learning (Lefebvre et al., 2012). The first one was characterized either by improvement in both speed and accuracy or by improvement of one parameter without a concomitant worsening of the other one; this resulted in a shift of the SAT, which suggests a rapid and successful motor skill learning. The second behavioral trajectory was characterized by opposite changes in speed and accuracy over time, resulting in less efficient motor skill learning. E.g., when speed improved, accuracy worsened, resulting in slight improvement in the SAT. Finally, the third behavioral trajectory was characterized by a deterioration of both speed and accuracy or a lack of any improvement. According to their behavioral trajectory, the subjects were refereed as shifters, fitters and non-learners respectively. These different behavioral trajectories were observed despite identical instructions and experimental conditions, and they were associated with specific brain activation patterns. Specifically, in the efficient shifters, activation was found in the M1, cerebellum and the SMA where the activation changes correlated with performance improvement, suggesting that the SMA plays a key role in early motor skill learning involving a SAT. In the less efficient fitters, there was only a non-significant correlation in the cerebellum (Lefebvre et al., 2012).
In stroke patients, functional brain imaging has been used extensively to explore the reorganization of the network controlling the paretic arm or hand. Grossly, early after stroke, this reorganized network is characterized by compensatory recruitment of the undamaged hemisphere, especially the motor and premotor areas (Feydy et al., 2002; Tombari et al., 2004; Jaillard et al., 2005; Ward and Frackowiak, 2006) and/or widespread activation in the damaged hemisphere with extensive activation of the somatosensory and premotor areas (Tombari et al., 2004; Ward et al., 2006). Over time, motor recovery is associated with a shift of activation back toward the damaged hemisphere (Pineiro et al., 2001; Jaillard et al., 2005; Favre et al., 2014; Grefkes and Ward, 2014) and a progressive recruitment of the cerebellum ipsilateral to the paretic hand (Small et al., 2002). In addition, changes in brain connectivity have been associated with motor function recovery after stroke (Schaechter et al., 2009; Jiang et al., 2013; Grefkes and Ward, 2014). Early after stroke, both anatomical and functional connectivity (FC) decrease within the damaged hemisphere; over time, motor function recovery is associated with gradual recovery of connectivity (Pannek et al., 2009; Westlake et al., 2012; Golestani et al., 2013). Thus, the more similar the reorganized motor network becomes to that of healthy individuals, the better the recovery. Nevertheless, the undamaged hemisphere may still play a vicarious role in recovered motor control of the paretic hand (Johansen-Berg et al., 2002; Werhahn et al., 2003; Tombari et al., 2004; Lotze et al., 2006; Bestmann et al., 2010; Grefkes and Ward, 2014). It has also been suggested that the resting-state FC correlates with the motor recovery potential (Park et al., 2011, 2014; Yin et al., 2012; Golestani et al., 2013; Ovadia-Caro et al., 2013; Dacosta-Aguayo et al., 2014).
Since functional reorganization occurs in the network supporting motor recovery of the paretic upper limb after stroke, it seems logical that similar neuroplasticity should occur in the network underlying motor skill learning. However, despite extensive fMRI studies of the functional neuroanatomy of motor skill learning in healthy individuals (Ghilardi et al., 2000; Halsband and Lange, 2006; Debas et al., 2010; Lefebvre et al., 2012; Hardwick et al., 2013), very few studies have assessed stroke patients. Using a region of interest (ROI) approach, one study with 10 chronic stroke patients performing visuomotor tracking with the paretic hand showed a bilaterally reorganized pattern with a predominance in the undamaged hemisphere during the pre-training fMRI session (Carey et al., 2002). After training, activation was partially transferred back toward the damaged hemisphere, suggesting functional reorganization (Carey et al., 2002). Another study using ROI showed decreased task-related fMRI activation in the contralesional M1 of nine chronic stroke patients after 3 days of training on a serial targeting task (Boyd et al., 2010). A recent fMRI study highlighted the differences in brain activation patterns between nine healthy individuals and nine chronic stroke patients during training over several days on an implicit sequential visuomotor tracking task (Meehan et al., 2011). Compared with healthy individuals, motor skill learning and retention in stroke patients relied on a reorganized network involving compensatory activations, especially in prefrontal attentional areas such as the DLPFC. Finally, during baseline performance of a sequential grip-force tracking task, 10 chronic stroke patients showed reduced fMRI activation in the damaged hemisphere compared with healthy controls (Bosnell et al., 2011). After repeated training, fMRI activation decreased in healthy controls but was maintained or increased in stroke patients.
These four studies involved relatively small cohorts of mostly high-functioning patients, typically with sub-cortical strokes, and they did not characterize motor skill learning through SAT. Instead, they compared fMRI activation related to motor performance pre- and post-training, and two used an ROI approach (Carey et al., 2002; Boyd et al., 2010). Since motor learning plays a key role in motor function recovery, better knowledge of the neurophysiology of motor skill learning after stroke should lead to the refinement of recovery models and translate into the development of specific neurorehabilitation methods based on the principles of motor learning.
Motor skill learning can be divided in two stages: a fast on-line learning process leading to large performance improvement over a single training session (i.e., early stages of motor skill learning as described is the present study), and a slower process involving smaller performance gains obtained through repeated training sessions (Dayan and Cohen, 2011).
This study aimed to specifically explore the early stages of motor skill learning with the paretic hand in chronic hemiparetic stroke patients, using an innovative motor skill learning paradigm with a SAT. This is a first step to understand the “recovery process” in stroke patients, whether residual motor learning aptitudes are present and which brain areas are (neuroplastically?) involved. A better knowledge about motor (skill) learning in stroke patients could help to refine neurorehabilitation protocols, in which motor learning is often imbedded as an implicit assumption but poorly recognized. The purposes of this study were: (i) to use random effect (RFX) analyses of whole-brain fMRI activation to identify the neural substrates underlying the first stages of motor skill learning involving a SAT in a larger cohort of chronic stroke patients using their paretic upper limb, (ii) to determine whether shifter and fitter stroke patients recruit specific neural substrates, and (iii) to determine whether resting-state FC acquired before training would predict the behavioral trajectory (shifter/fitter) during the first stage of motor skill learning and/or correlate with the amount of motor skill learning.
Material and Methods
The experimental protocol was approved by the local Ethical Committee (CHU Dinant-Godinne, UCL) and was conducted according to the recommendations of the Helsinki declaration. After providing written informed consent, 25 chronic stroke patients meeting the following criteria were selected: Inclusion: (i) being a chronic (>6 months) stroke patient aged 18-90 years, (ii) presenting a chronic motor deficit in the upper limb (i.e., chronic hemiparesis), (iii) having a vascular brain lesion demonstrated by cerebral imaging (Figure 1); Exclusion: (i) being unable to perform the task or to understand instructions, (ii) presence of intracranial metal, (iii) alcoholism, (iv) pregnancy, (v) cognitive impairment or psychiatric disorder, (vi) any contraindication to MRI. The following measures were assessed: disability with the modified Rankin Scale (mRS), level of impairment with the National Institutes of Health Stroke Scale (NHISS), residual dexterity with the Purdue Pegboard Test (PPT), maximal hand force (MaxHF) with a whole-hand Jamar dynamometer and manual ability with the ABILHAND scale (Penta et al., 2001) (Table 1). Patients #12, 13, and 15 had participated in a previous study at least 1 year before, exploring motor skill learning enhancement by transcranial direct current stimulation (tDCS) (Lefebvre et al., 2013). They were included in the present study because this circuit was different from those used in the previous experiment, so they learned from scratch a new circuit, as did the other patients. These three patients showed no significant difference from the other patients on (i) their new baseline motor performances and (ii) their new performance evolution (Table 2). This analysis was performed using the Crawford and Howell statistical test to compare an individual score to a small population (Crawford and Howell, 1998; Crawford and Garthwaite, 2002). Additional fMRI analyses did not demonstrate difference with the other patients (Supplementary Figure 1 and Supplementary Table 1).
Figure 1. Stroke localization and overlap. Upper panel: T1 magnetic resonance imaging (MRI) at the level of the main stroke injury. Lower panel: Lesion overlap in stroke patients symbolized by the color scale. Purple represents the stroke area of a single patient, green represents localization shared by half of the patients, and red indicates localization shared by all of the patients. For patients with lesions on the right side of the brain, the 3D-T1 MRI was flipped. The map of lesion localization and overlap was created with MRIcro 1.4. DamH: damaged hemisphere.
Prior to the fMRI activation runs, a resting-state run was acquired except in patients #1–4 and #17. The resting-state acquisition consisted of one single 6-min run during which the patients kept their eyes closed and had to avoid moving or falling asleep.
Then, the patients performed two consecutive activation runs of motor skill learning with their paretic upper limb, using a MR-compatible mouse (NAtA Technologies, Canada). Visual feedback was projected on a screen, which was viewed via a mirror placed on the head coil. As described in a previous study (Lefebvre et al., 2012), each run (duration 8 min 41 s; 172 volumes) contained a REST condition (fixation cross) and three experimental conditions: LEARNING, EASY, and REPLAY. LEARNING required performing the motor skill learning paradigm described below. EASY required moving the cursor back and forth between two horizontal or vertical targets, without speed or accuracy constraint. This condition was designed to explore brain activation related to simple movement execution under visual control. During REPLAY, a video clip of the last LEARNING performance was displayed, and the patients were instructed to follow the cursor's displacement with their eyes without moving the hand. The REPLAY condition was designed to isolate activation related to visual and oculomotor processes.
Each condition was presented four times (30 s each; 10 volumes) during each run, and conditions were separated by four volumes (12 s) of REST (Figure 2). Beforehand, patients performed one habituation run (4 × 30 s separated by four REST volumes) of a simple circuit square.
Figure 2. 1 simple square circuit for the habituation run (design not shown). 2 (green): LEARNING: subjects were instructed to move the cursor within the track as quickly and accurately as possible; 3 (red): REPLAY: subjects were instructed to follow the cursor displacement with their eyes while watching a video clip of their last LEARNING block, keeping the hands relaxed; 4 (blue): EASY: subjects were instructed to move the cursor between the two targets at a comfortable speed (50% trials with vertical movements, 50% with horizontal movements); 5: REST (gray): fixation cross.
Motor Skill Learning Paradigm
The motor skill learning paradigm consisted in moving a cursor controlled by an MR-compatible mouse along a complex circuit under visual feedback (Lefebvre et al., 2012, 2013). Patients were asked to complete as many laps as possible while keeping the cursor within the track (SAT) and to improve their performance over time.
Error was defined as the surface area between the actual trajectory and the ideal trajectory (i.e., midline of the track since the cursor had the same width as the track). Velocity was the first derivative of the position. Velocity and error during the LEARNING blocks were averaged as mean error and mean velocity over 3 s time bins [corresponding to repetition time (TR)]. To quantify motor skill learning as change of SAT, error and velocity were combined. Based on mean velocity and error, four indices were computed (for details see Lefebvre et al., 2012, 2013): error index (Pe) = constant error/patient mean error; velocity index (Pv) = patient mean velocity/constant velocity; Performance Index (PI) = Pv * Pe over bins of 3 s, then averaged for each LEARNING block; and finally Learning Index (LI) = [(PI − PI initial)/PI initial] * 100, which was computed for each LEARNING block as a percentage of the PI relative to the baseline performance during the first block of LEARNING (PI initial). Based on the LI changes over time observed previously in healthy individuals and chronic stroke patients (Lefebvre et al., 2012, 2013), motor skill learning could follow three behavioral trajectories, as described in the Introduction: shift, fit or no learning.
Student's t-tests were performed to determine whether baseline clinical characteristics (age, time since stroke, mRS, NIHSS, PPT, and ABILHAND) were related to the shift or fit behaviors. Similarly, Chi-square tests were calculated to investigate differences between subgroups (shifter or fitter stroke patients) based on stroke localization (cortical/subcortical), gender and whether the paretic hand was dominant or non-dominant. These comparisons were done with SigmaStat 3.5 (Systat Software San Jose, CA).
Imaging Acquisition Parameters
A 3T scanner with a 32-channel head coil (Siemens Verio, Germany) was used to acquire fMRI scans of brain activity with repeated single-shot echo-planar imaging, using the following parameters: TR = 3000 ms, echo time (TE) = 23 ms, flip angle (FA) = 90°, matrix size = 64 × 64, field of view (FOV) = 220 × 220 mm2, slice order descending and interleaved, slice thickness = 2 mm (no gap) and number of slices = 59 (whole brain), voxel size = 23.8 mm3 (3.4375 * 3.4375 * 2). Resting-state scans were acquired via repeated single-shot echo-planar imaging with the following parameters: TR = 1800 ms, TE = 23 ms, FA = 80°, matrix size = 64 × 64, FOV = 240 × 240 mm, slice order ascending and interleaved, slice thickness = 3.8 mm (no gap) and number of slices = 35 (whole brain), voxel size = 53.4 mm3 (3.75 * 3.75 * 3.8). The whole brain was acquired 200 times (200 volumes = 6 min). The whole brain was scanned 172 times per learning run and 60 times during the habituation run. A three-dimensional (3D) T1-weighted data set covering the whole brain was acquired (1 mm3, TR = 1600 ms, TE = 2.39 ms, FA = 9°, matrix size = 512 × 512, FOV = 256 × 256 mm3, 176 slices, slice thickness = 1 mm, no gap).
Pre-Processing and Statistical Analyses
fMRI data were analyzed using BrainVoyager QX (Version 22.214.171.1240, Brain Innovation, The Netherlands). For patients with stroke lesions on the right side of the brain, both 3D-T1 and functional data were flipped. The pre-processing of the functional data consisted of a slice time scan correction, temporal high-pass filtering (removing temporal frequencies below three cycles/run) and 3D motion correction for head movements using a rigid body algorithm. If a movement >2 mm (or degrees) was observed in one of the 6 directions, the related movement parameters were using as regressors of non interest to remove potential residuals motion artifacts. Functional data were analyzed using the General Linear Model (GLM), consisting of predictors based on specific experimental conditions, in which beta weights measure the potential contribution of the predictors in each voxel time course. Coregistrations using gradient-driven affine transformations with 9 alignment parameters (3 translations, 3 rotations and 3 scales) between functional runs and 3D-T1 weighted scans of each patient were performed automatically, and then manually corrected. All anatomical and functional volumes were spatially normalized in Talairach space (Talairach and Tournoux, 1988) to allow group analysis. The statistical maps obtained were overlaid on the 3D T1-weighted scans. Functional runs were smoothed in the spatial domain with a Gaussian filter of 5 mm. The resting-state scans were not smoothed.
First, a whole-group RFX analysis involving the 23 patients who achieved learning was constructed. The brain activation associated with each condition was explored at the whole-brain level using the following contrasts [LEARNING] (main condition); [EASY] (areas involved in lower aspects of movement control and execution); [REPLAY] (areas involved in visual and oculomotor activity minus REST activation); and [LEARNING - (REPLAY + EASY)]. This contrast was designed to explore the activation underlying the first stages of motor skill learning after subtracting the activation related to lower aspects of movement control and to visual-oculomotor activity (which were different between [REPLAY] and [EASY] since their respective visual feedbacks were different). Correlation analyses between the fMRI activation during [LEARNING] and the patient's clinical characteristic (the mRS score, the NIHSS score, the Abilhand score and the PPT score of the paretic hand) were performed.
Pearson correlation analyses were performed to identify the key area(s) where activation changes had the highest correlations with motor skill learning. For this analysis, performance (PI) values of each patient were averaged, and the correlation was performed between the eight beta weights [one averaged value for each block (mean beta for the whole group at each block)] of each ROIs with significant activation with [REPLAY], [LEARNING] and [LEARNING − (REPLAY + EASY)] and the eight PI values [one averaged value for each block (mean PI for the whole group at each block)]. The beta weights were extracted as the mean beta weight of the ROI. The correlations were done with SigmaStat 3.5.
For each patient (n = 20) an independent component analysis (ICA) extracted 40 components, using the “FastICA” algorithm implemented in BrainVoyager. Then, a self-organizing group-level ICA was used to classify the 40 components based on a clustering procedure that accepts only one component per subject in each cluster. Across these 40 clusters, five networks of interest were identified: the default-mode network (DMN) encompassing the posterior cingulate cortex, the medial prefrontal cortex, the lateral parietal cortex and the parahippocampal gyrus (Buckner et al., 2008); the somatomotor network including the primary and higher order motor and sensory areas (Biswal et al., 1995; Carter et al., 2010; Laird et al., 2011), the salience network including the anterior insula and the dorsal anterior cingulate (Seeley et al., 2007), the dorsal attentional network (DAN) including frontoparietal areas (Carter et al., 2010), and the visual network (VN) involving much of the occipital cortex (Beckmann et al., 2005).
The baseline PI values, the mRS score, the NIHSS score, the Abilhand score, the PPT score of the paretic hand and the amount of learning at the end of the session (LI) were used as covariates to assess whether the resting-state FC in the pre-mentioned networks correlated with the baseline motor performances or motor skill learning potential of the patients.
Subgroup (Shifters/Fitters) Analyses
A whole-brain 2 factors ANOVA using condition estimates (beta values) from the first-level whole-group RFX GLM analysis constructed with 46 runs (23 stroke patients participated in a single sessions with two runs) was performed to directly compare the activation associated with the first stages of motor skill learning between each patients' subgroup [first factor: fMRI conditions (LEARNING, EASY, and REPLAY), second factor: patients' subgroup (shifters/fitters)]. The age and time since stroke were used as additional covariates.
The mean beta weights of the ROIs involved in the first stages of motor skill learning (identified with [LEARNING - (REPLAY + EASY)] in the whole-group RFX analysis) were used to perform Student t-tests to compare respective activation in these ROIs between shifter and fitter stroke patients (one averaged beta weigh per condition for each patient).
Separate RFX subgroup analyses
Two additional whole-brain RFX were constructed separately for each subgroup (i.e., shifters and fitters). The brain activation associated specifically with the first stages of motor skill learning was explored at the whole-brain level using the [LEARNING − (REPLAY + EASY)] contrast.
Subgroup Pearson correlations
Pearson correlation analyses were performed to identify the key area(s) where activation changes had the highest correlations with motor skill learning separately for each subgroup. The correlations were performed between the eight beta weights of each ROIs with significant activation with [LEARNING − (REPLAY + EASY)] and the eight PI values separately for each subgroup. The beta weights are extracted as the mean beta weights of the ROIs.
Resting-state FC comparison between subgroups
Finally, a two factors whole-brain ANOVA was performed to compare the five pre-specified networks between the two subgroups [first factor: resting-state networks (DMN, somatomotor network, salience network, DAN, VN), second factor: subgroups (shifters, fitters)].
Of the 25 stroke patients, two (#2 and 6) were classified as non-learners since their motor performance deteriorated; they were excluded from further analyses. The 23 remaining patients achieved motor skill learning: nine were classified as shifters and 14 as fitters (Figure 3; Supplementary Figure 2). At the end of the second learning run, the performance of the shifters had improved significantly more (LI: 49 ± 30%; mean ± SD) than that of the fitters (LI: 13 ± 10%; p < 0.001). Raw PI scores are listed in Supplementary Table 2 and SAT for each patient are illustrated in Supplementary Figure 2.
Figure 3. Learning Index (LI) evolution across learning blocks. The eight LI values (mean ± SD) correspond to the LI during each learning block. Black squares: shifter stroke subgroup (n = 9; LI, 49 ± 30% at the last block); white triangles: fitters stroke subgroup (n = 14; LI: 13 ± 10%; p = 0.0005).
Age, time since stroke, mRS, NIHSS, PPT, and ABILHAND scores were not significantly different between the two subgroup (p = 0.43, 0.76, 0.66, 0.90, 0.29, and 0.73, respectively, Student's t-test). Finally, the two subgroups were also not different based on the type of stroke (cortical versus subcortical; p = 0.47, Chi-square), the gender (p = 0.96, Chi-square) or whether the paretic hand was dominant or non-dominant (p = 0.18, Chi-square).
Whole-group RFX analysis
As shown in Figure 4 and Table 3, whole-group RFX analysis revealed the areas activated during each condition [t(22) = 2.19; pUNCORRECTED < 0.04]. Areas activated during LEARNING were: M1 [Brodmann Area (BA) 4] in the damaged hemisphere (M1damH), the dorsal premotor cortex (PMddamH, BA 6), SMAdamH (BA 6), bilateral posterior parietal cortex (PPC, BA 7), bilateral primary somatosensory cortex (S1, BA 3), DLPFCdamH (BA 9) and bilateral visual cortex. Areas activated during EASY were: M1damH, PMddamH, bilateral SMA, bilateral PPC, inferior parietal cortex (IPCdamH, BA 40), S1damH, bilateral anterior cerebellum and bilateral visual cortex. As expected, the intense visual area activations during LEARNING were also observed in REPLAY, as well as activation in bilateral ventral premotor cortex (PMv), bilateral PMd, thalamusundamH, bilateral putamen, bilateral PPC and IPC. No significant correlation was observed between the fMRI activation during LEARNING and the baseline clinical characteristic of the patients.
Figure 4. Whole-group activation. BOLD activation for the 23 chronic stroke patients in the following contrasts, [LEARNING], [REPLAY], and [EASY] [RFX t(22) = 2.19; pUNCORRECTED < 0.04], and [LEARNING − (REPLAY + EASY)] [RFX t(22) = 2.10; pUNCORRECTED < 0.04]. DamH, damaged hemisphere.
Next, to focus on the areas involved in the first stages of motor skill learning, the [LEARNING − (REPLAY + EASY)] contrast was computed [t(22) = 2.10; pUNCORRECTED < 0.04], it showed significant activation in the M1damH (44 mm2) and PMddamH (42 mm2).
Pearson correlation analyses
At the whole-group level, correlation analyses between the PI and beta weights of each ROI activated in [LEARNING] in the whole-group RFX analysis, showed a statistically significant positive correlation in the PMddamH (r = 0.70, p = 0.05) and a negative correlation in the DLPFCdamH (r = −0.82, p = 0.01). Correlations in the other areas activated in [LEARNING] (i.e., M1damH, SMAdamH, bilateral S1, bilateral PPC, and bilateral visual cortical areas) were not statistically significant (r ranging between −0.16 and +0.25).
Correlation analyses performed based on [LEARNING − (REPLAY + EASY)] showed a statistically significant correlation exclusively in the PMddamH (r = 0.71, p = 0.048). It is worth noting that correlation analyses between the PI and beta weights of each ROI activated in [REPLAY] did not demonstrate any statistically significant correlation.
The resting-state FC was obtained for the 5 following networks at the whole-group level (n = 20): the DMN, the DAN, the VN, the somatomotor network, and the salience network [Figure 5, (qFDR) < 0.05, t(76)].
Figure 5. ICA resting-state networks at the whole-group level. Statistical maps are thresholded at (qFDR) < 0.05. DMN, default mode network; VN, Visual Network; DAN, dorsal attentional network. DamH, damaged hemisphere.
No significant correlation were observed in any of these five networks between the resting-state FC and the baseline clinical characteristic of the patients or the amount of LI at the end of the session.
Subgroup (Shifters/Fitters) Analysis
The whole-brain ANOVA revealed a significant effect of the two factors and their interaction [fMRI conditions: F(2, 42) = 6.27, q(FDR, False Discovery Rate) = 0.05; patients' subgroup: F(1, 21) = 6,27, p = 0.02; interaction: F(2, 42) = 6,27, p = 0.004].
The post-hoc analyses contrasted the patients' subgroup for each fMRI condition. There was no significant difference between the two patients' subgroup for [EASY] [t(44) = 3.00; q(FDR) = 0.05] and for [REPLAY] [t(44) = 3.00]. With [LEARNING], in the fitters subgroup 23 voxels in the posterior parietal cortex of the damaged hemisphere [PPCdamH, BA 7, (mean x = −17, mean y = −66, mean z = 61)] were significantly more activated; in the shifters subgroup 16 voxels in the cerebellum ipsilateral to the paretic hand (mean x = 38, mean y = −37, mean z = −27) were significantly more activated [t(44) = 3.00].
In each ROI significantly activated, Student t-tests were performed to compare the BOLD activation between shifter and fitter patients for [LEARNING − (REPLAY + EASY)]. These comparisons, summarized in Table 4 and in Supplementary Table 3 (for additional contrasts), showed a statistically significant differences exclusively in the PMddamH (shifters > fitters, pBonferroni−corrected = 0.02) and a non-significant trend in PMdundamH (shifters > fitters, pBonferroni−corrected = 0.15).
Table 4. Differential activation between “fitter” and “shifter” stroke patients (whole-group RFX analysis).
Separate RFX subgroup analyses
To compare the activation patterns in the shifter and fitter patients, separate RFX subgroup analyses were computed with [LEARNING − (REPLAY + EASY)] (Figure 6, Table 5). In the shifters subgroup [t(8) = 2.31; pUNCORRECTED < 0.05], five areas were activated: the SMA, bilateral PMd, M1damH, and S1damH. In the fitters subgroup [t(13) = 2.17; pUNCORRECTED < 0.05], only the bilateral PPC was significantly activated.
Figure 6. Activation in the shifters and fitters stroke subgroups. BOLD activation for [LEARNING − (REPLAY + EASY)] for shifter stroke patients [RFX t(8) = 2.31; pUNCORRECTED < 0.05] and fitter stroke patients [RFX t(13) = 2.17; pUNCORRECTED < 0.05]. Note that the activation in the shifter stroke subgroup is distributed in a sensorimotor/premotor network (SMA, bilateral PMd, M1damH, and S1damH), whereas significant activation in the fitters stroke subgroup is limited to the bilateral PPC. DamH, damaged hemisphere.
Analyses based on additional areas classically involved in motor skill learning but not revealed in the present study were performed (e.g., in healthy volunteers: the cerebellum ipsilateral to the working hand, (Lefebvre et al., 2012); in stroke patients: the DLPFCdamH, (Meehan et al., 2011); and the striatumdamH, Bosnell et al., 2011). No significant difference in the beta weights between shifters and fitters subgroups was found in the cerebellum [ROI from (Lefebvre et al., 2012): x = −21, y = −45, z = −21, 73 mm3; p = 0.39], nor in the DLPFCdamH (ROI from whole-group analysis of [LEARNING]: x = −40, y = 42, z = 30, 363 mm3; p = 0.24), nor in the striatumdamH (ROI sphere of 3 mm centered on x = −26, y = 7, z = 5). Ten patients were excluded from this analysis centered on the striatumdamH, since their stroke involved part of the striatum (#3, #8, #12, #13, #14, #17, #21, #22, #23, #24).
Subgroup Pearson correlation analyses
Based on the RFX in the shifters subgroup, the correlation analysis between the PI and beta weight of each ROI activated in [LEARNING − (REPLAY + EASY)] revealed a significant correlation exclusively in the PMddamH r = 0.91, pBonferonni−corrected = 0.01) and a non-significant trend in the PMdundamH (r = 0.79, pBonferonni−corrected = 0.1). In the fitters subgroup [t(13) = 2.17; pUNCORRECTED < 0.05], only the bilateral PPC was significantly activated but there was no significant correlation between the PI and beta weight changes.
In both the fitters and shifters subgroups, in each area where a significant correlation between the beta weight and the PI was found, there was also a positive correlation with the velocity's change and a negative correlation with the error's change. No areas showed correlations exclusively with error or velocity change.
Resting-State FC Comparison between Subgroups
A whole-group ANOVA was created to compare the FC between the two patients' subgroup [shifters (n = 9)/fitters (n = 11)]. At qFDR < 0.05 for the DMN, the VN, the DAN and the somatomotor network no statistically significant differences were observed between the two patients' subgroup. By contrast, there was a statistically significant difference for the salience network where the FC was significantly higher in the shifters compared with the fitters subgroup [qFDR < 0.05; t(76) = 4.78] especially in the anterior cingulate cortex (BA 24; 124 voxels; x = 0, y = 34, z = 7), insula [BA 13; 39 voxels; x = 35, y = 20, z = 5] and temporal lobe [BA 42; 26 voxels; x = 53, y = −13, z = 10] (Figure 7).
Figure 7. Resting-state FC Whole-brain ANOVA. (A) FC in the salience network in the Whole-group (qFDR < 0.05), (B) FC in the salience network in each subgroup (qFDR < 0.05). (C) The result of the whole-brain ANOVA (qFDR < 0.05), comparing the shifters with the fitters subgroups in the salience network, with higher FC in the shifter subgroup. DamH, damaged hemisphere; red circle, anterior cingulate cortex (BA 24); green circle, insula (BA 13); blue circle, temporal lobe (BA 42).
In a group of 23 chronic stroke patients, the first stages of motor skill learning with the paretic upper limb were associated with a bilateral fMRI activation pattern characterized by a positive correlation in the PMddamH and a negative correlation in the DLPFCdamH. After subtracting activation related to visual/visuomotor processes and lower aspects of motor control, correlation with motor skill learning was restricted to the PMddamH. The FC acquired with resting-state fMRI before the learning runs showed correlation neither with the baseline motor performance nor with the amount of motor skill learning. Comparison between the fitters and shifters stroke subgroups revealed two distinct brain activation patterns. In the less efficient fitters, significant fMRI activation was restricted to the bilateral PPC but did not correlate with motor skill learning. In contrast, in the more efficient shifters, fMRI activation encompassed the bilateral PMd, SMAdamH, M1damH and S1damH. The bilateral PMd, especially PMddamH, showed the most significant correlation between fMRI activation changes and the first stages of motor skill learning in shifters. In the resting-state salience network, the FC measured before learning was higher in the shifters subgroup; no significant difference was found in the other resting-state networks.
Learning a Motor Skill with the Paretic Upper Limb
A progressive worsening of motor performance was observed in two patients (8%). There were no obvious differences in demographic or stroke characteristic between these two non-learners and the other patients. We speculate that an excessive fatigue or a lack of motivation may have played a role, but these aspects were not assessed. Alternatively, in these individuals, strokes may have destroyed key areas that support motor skill learning. However, this seems unlikely for two reasons. First, their strokes were similar in extent and location to those of the 23 other patients. Second, to date, there is no compelling evidence that a specific brain injury may abolish motor skill learning with the paretic upper limb. Indeed, stroke patients remain able to learn new motor skills (Platz et al., 1994; Carey et al., 2002; Boyd et al., 2010; Meehan et al., 2011) and motor skill learning is impaired rather than abolished after a focal lesion of the basal ganglia (Vakil et al., 2000; Exner et al., 2001) or prefrontal cortex (Gomez Beldarrain et al., 2002).
Reorganized Motor Skill Learning Network after Stroke
In chronic hemiparetic stroke patients, the evolution of the beta weights across the two runs correlated positively with successful (early) motor skill learning in the PMddamH and negatively in the DLPFCdamH. In healthy individuals, we previously found a (positive) correlation only in the SMA (Lefebvre et al., 2012). The visual comparison of fMRI activations between non-age matched healthy individuals and hemiparetic stroke patients shown obvious pattern differences (see Supplementary Materials, and Supplementary Figure 3). Although the SMAdamH was also recruited in chronic hemiparetic stroke patients, the PMddamH was the key area correlating with the first stages of motor skill learning, suggesting a neuroplastic reorganization. Previous motor skill learning studies in stroke patients led to conflicting observations about PMd. After extensive tracking training, performance improvement was associated with an increased activation in PMddamH at the expense of PMdundamH (Carey et al., 2002). In another study, activation in PMdundamH correlated with motor sequence accuracy at retention (Meehan et al., 2011). The current study clearly showed that PMddamH plays a key role during the early stage of learning a motor skill involving a SAT with the paretic hand. This - likely compensatory - activation in PMddamH is coherent with the observations that PMd (i) is crucial for motor skill learning in healthy individuals (Kantak et al., 2012; Hardwick et al., 2013), (ii) is involved in motor function recovery after stroke (Carey et al., 2002; Fridman et al., 2004; Tombari et al., 2004; Lotze et al., 2006; Bestmann et al., 2010; Schulz et al., 2012), and (iii) is part of the reorganized network involved in stroke patients performing simple motor tasks (Johansen-Berg et al., 2002; Calautti et al., 2010).
A recent study in stroke patients showed strong compensatory activation of the bilateral DLPFC and prefrontal BA 46 in the undamaged hemisphere during motor sequence learning and of the prefrontal BA 8 in the undamaged hemisphere during the retention test (Meehan et al., 2011). In the current study, activation in the DLPFCdamH correlated negatively with performance improvement during the first stages of motor skill learning. This may suggest an early, transient recruitment of prefrontal attentional areas (here, DLPFCdamH) in parallel with rising activation of the (pre-)motor network in stroke patients.
During REPLAY, intense activation was observed in the visual areas, as well as widespread activation including the bilateral PMv, PMd, putamen, PPC, IPC, and the thalamusundamH. However, there was no correlation between fMRI activation in these areas and motor skill learning, demonstrating that REPLAY was not involved in learning this skill through rehearsal or memory retrieval.
During the simple motor condition EASY, the bilateral cerebellum was activated, consistently with its role (i) in motor control (Sakai et al., 1998; Miall et al., 2001) and (ii) in recovered motor function after stroke (Small et al., 2002). In contrast, during motor skill learning, there was neither significant activation nor correlation within the cerebellum. The heterogeneity of the strokes might have prevented reliable recruitment of the cerebellum. However, this seems unlikely since consistent cerebellar activation was observed during simple motor performance (EASY). It is worth noting that no change in cerebellar activation associated with motor learning was reported in previous whole-brain studies in subcortical stroke patients (Bosnell et al., 2011; Meehan et al., 2011). We speculate that, in contrast to healthy individuals, this might reflect a preferential or more consistent recruitment of the premotor and/or prefrontal areas rather than the cerebellum during the first stages of motor skill learning.
Post-stroke Neuroplasticity or Ageing-related Modifications?
Could the observed activation pattern reflect ageing-related reorganization, rather than post-stroke neuroplasticity? Previous studies demonstrated a modification of the motor control/execution network associated with ageing (Mattay et al., 2002; Heuninckx et al., 2008). It is thus logical to infer that ageing would also influence the motor skill learning network, even though such a reorganization has yet not been formally demonstrated (Daselaar et al., 2003). Furthermore, studies comparing fMRI activation between stroke patients and age-matched healthy individuals consistently found striking differences suggestive of post-stroke neuroplastic reorganization, both for motor skill learning (Carey et al., 2002; Bosnell et al., 2011; Meehan et al., 2011) and motor performance (Zemke et al., 2003; Ward and Frackowiak, 2006; Schaechter and Perdue, 2008). This clearly demonstrates that the impact of stroke on the reorganization of brain activation is considerably greater than the impact of ageing. In addition, an external multiple regression analysis (data not shown) combining the healthy individuals (Lefebvre et al., 2012) and stroke patients did not demonstrate a statistically significant correlation between age (18–82 years old) and fMRI activation in all the activated ROIs with [LEARNING] (M1, SMA, PMd, S1, PPC, thalamus, DLPFC and cerebellum). Therefore, we conclude that the observed differences in fMRI patterns between chronic stroke patients and healthy individuals (Lefebvre et al., 2012) predominantly reflect post-stroke neuroplastic reorganization rather than ageing. In other words, chronic stroke patients recruited a reorganized (likely compensatory) network to achieve motor skill learning.
Differential Brain Activation in Shifter and Fitter Stroke Patients
As in healthy individuals under similar experimental conditions (Lefebvre et al., 2012), the chronic stroke patients could be classified as shifters, fitters or non-learners. In a previous study using the same motor skill learning paradigm, more chronic stroke patients behaved as fitters rather than as shifters after 30 min of motor skill learning under control condition (sham transcranial direct current stimulation) (Lefebvre et al., 2013). These observations suggest that this motor skill learning paradigm based on a SAT allows a fine dissection of behavioral trajectories in healthy individuals and stroke patients. Demographic factors, stroke characteristics [time since stroke or stroke localization (cortical vs. subcortical lesion)] and level of impairment did not predict whether stroke patients would behave as shifters or fitters. It remains to be determined whether a genetic background, such as a polymorphism of the brain-derived neurotrophic factor (BDNF) (Bath and Lee, 2006; Kleim et al., 2006; McHughen et al., 2010), underlies these different behaviors.
In the shifter stroke subgroup, activation was found in the bilateral PMd, SMAdamH, M1damH, and S1damH; and a significant correlation between motor skill learning and fMRI activation changes was found exclusively in the bilateral PMd, especially PMddamH. This suggests that the shifter stroke patients recruited neuronal resources from both hemispheres but predominantly from the damaged hemisphere to achieve efficient motor skill learning. The stronger involvement of the bilateral PMd (especially PMddamH) compared with the SMAdamH (strongly involved in healthy shifters, see Supplementary Materials) likely reflects a neuroplastic reorganization of the motor skill network. Finally, the whole-brain ANOVA demonstrated a higher level of activation in the cerebellum ipsilateral to the paretic hand in the shifters compared with fitters, suggesting a more consistent recruitment of the cerebellum in the more efficient shifter stroke patients.
As previously observed in healthy fitters (Lefebvre et al., 2012), fitter stroke patients seemed unable to efficiently or consistently engage the motor learning network observed in shifter stroke patients. Here, significant activation was found exclusively in the bilateral PPC and did not correlate with motor skill learning. Furthermore, the whole-brain ANOVA confirmed that the PPCdamh was more activated in the fitters stroke subgroup compared with the shifters subgroup. The PPC is involved in the planning of visuomotor tasks (Desmurget et al., 1999; Torres et al., 2013). This may suggest that the fitter stroke patients were focused (or jammed) on the visuospatial aspects of the task, possibly reflecting a protracted exploratory phase.
The current fMRI results shed new light on a previous observation in chronic stroke patients. Compared with sham, tDCS applied bilaterally over M1 (dual-tDCS) enhanced motor skill learning in 100% of chronic stroke patients (n = 18) (Lefebvre et al., 2013). In sharp contrast, after sham dual-tDCS, 44% of the stroke patients (n = 8) showed performance worsening (non-learners). Given the spread of the direct current delivered through tDCS electrodes centered over M1, it is possible that the bilateral PMd were also stimulated. If bilateral PMd are indeed crucial for efficient motor skill learning with the paretic hand in chronic stroke patients, as suggested by the current fMRI observation, stimulation of PMd could partially explain the enhancement of motor skill learning and long-term retention by dual-tDCS in stroke patients (Lefebvre et al., 2013).
Baseline Resting-State FC
In contrast with previous studies which showed that resting-state FC correlates with residual motor function in chronic stroke (Carter et al., 2010; Yin et al., 2012), no correlation between the level of motor performance and the baseline FC in any of the five resting-state networks of interest was observed in the current study. As discussed in the limitations section, this might be related to the heterogeneity of our chronic stroke patients and the smaller sample (n = 20) for resting-state fMRI. The choice of the outcome measure (scales) may also be of key importance when looking for correlation between resting-state FC and residual functions in chronic stroke patients.
Previous studies suggested that resting-state-fMRI FC might also be of interest for predicting the recovery potential after stroke (Park et al., 2011, 2014; Golestani et al., 2013; Ovadia-Caro et al., 2013; Dacosta-Aguayo et al., 2014). In the current study, the only correlation between baseline resting-state FC and subsequent motor skill learning was observed in the salience network where FC was higher in the shifters subgroup compared with the fitters subgroup. No difference was observed in the others resting-state networks. Thus, a higher (i.e., supposed “better”) resting-state FC within the salience network before training may predispose chronic stroke patients to learn more efficiently a new motor skill with a SAT. Since the salience network is involved in the detection of pertinent stimuli to perform a task and in the coordination of the neural resources needed to achieve this task (Uddin, 2015), one may hypothesize that the salience network plays a key role in recruiting the attentional resources that facilitate learning. The anterior insula triggers attentional shift to salient event in the anterior cingulate cortex which in turn modulates the activation of PMd (Menon and Uddin, 2010). One could thus hypothesize that the higher connectivity in the anterior cingulate cortex and anterior insula in the shifter stroke patients could favor the early recruitment of PMd during motor skill learning. This observation needs to be confirmed in a larger stroke cohort. The lack of significant correlation within the other resting-state networks may be due to the already discussed limitation. Furthermore, it is also possible that this motor skill learning session was too short, compared with the longitudinal recovery driven by neurorehabilitation (Park et al., 2011, 2014; Golestani et al., 2013; Ovadia-Caro et al., 2013; Dacosta-Aguayo et al., 2014).
Limitations of the Study
This study had four main limitations. First, patients were relatively heterogeneous in terms of stroke localization (cortical, subcortical, brainstem) and etiology (large arteries, lacunar infarcts, intracerebral hemorrhages). This heterogeneity might have enhanced the inter-individual variability and dampened the significance of the activation and resting-state BOLD signals. However, this relative heterogeneity may also be considered to be a strength, since this cohort represents real-life hemiparetic stroke patients.
Second, for all the RFX analyses (excepted for the whole-brain ANOVA), a liberal threshold of pUNCORRECTED < 0.04 was applied. This threshold was chosen given the heterogeneity of the stroke population and the variability associated with large brain lesions. Due to this variability and to the technical difficulties associated with stroke patients, the level of statistical significance warranted in a population of healthy individuals is difficult to achieve, unless the recruitment is restricted to a subpopulation, typically patients with slight impairment and a single subcortical stroke. Nevertheless, the results described at group level were also confirmed at the individual level (individual correlations, paired t-test ROIs comparison) and were correlated with the behavioral performance evolution. Thus, even if these results have to been confirmed, they are both informative and relevant for the field.
Third, the results of the correlation analyses presented here are preliminary and need to be confirmed in a larger cohort of stroke patients. It has however to be noted that compared to the previous studies exploring motor skill learning in stroke patients (n = 9 or 10), a larger sample (n = 23; i.e., 25 − 2) was involved in the current study.
Fourth, a formal comparison with a group of age-matched control subjects has not been performed. However, as discussed previously, the reorganization of the motor skill learning network we observed in stroke patients is coherent with previous studies showing post-stroke compensatory reorganization of the motor control network in aged patients (Carey et al., 2002; Jaillard et al., 2005). Further studies are required to unveil the full dynamic evolution of motor skill learning in stroke patients compared with healthy controls.
Compared with previous studies based on ROI analyses and/or smaller cohorts of subcortical stroke patients, the current study represents several important advances. First, the use of whole-brain fMRI analysis unveiled a network dynamically engaged during the first stages of motor skill learning in chronic hemiparetic stroke patients. Second, the application of RFX analyses to a relatively large cohort of chronic patients with various types of stroke identified more general reorganization mechanisms which are likely shared by most hemiparetic chronic stroke patients. Third, careful behavioral dissection of this motor skill learning paradigm involving a SAT provided an unprecedented level of precision to investigate the fine neuronal mechanisms underlying the first stages of motor skill learning in stroke patients.
When learning a new visuomotor skill involving a SAT with the paretic upper limb, chronic stroke patients presented a reorganized brain activation pattern involving both hemispheres with a predominant recruitment of the damaged hemisphere. The key area underlying efficient motor skill learning was bilateral PMd and especially PMddamH, in contrast to the key role played by the SMA in healthy individuals. This suggests a neuroplastic, compensatory recruitment of additional areas during motor skill learning in chronic hemiparetic stroke patients. Surprisingly, the correlations between baseline resting-state FC and residual motor function or motor learning potential revealed only a higher resting-state FC in the salience network in the shifter stroke subgroup. Further explorations of FC as a biomarker of recovery potential are needed. A better understanding of the neural substrates underlying motor learning in stroke patients is a crucial step forward to design the next generation of neurorehabilitation paradigms.
The work of YV was supported by the following grants: Fonds de la Recherche Scientifique Médicale (FRSM) 3.4.525.08.F in 2008, 2010 and 2012, Fonds Spécial de Recherche (FSR) grant from the Université Catholique de Louvain (UCL) in 2008 and 2010, and Fondation Van Goethem-Brichant. The work of SL was supported by UCL FSR grants in 2008 and 2010, and by a grant from the Fondation Mont-Godinne 2012. The work of LD was supported by a grant from the Fondation LOUVAIN (UCL).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the stroke patients who participated in this experiment and their families, the collaborators from Siemens (Roger Demeure, Dimitri Haye, Jean-Pierre Bastin and Patrick Caluwaerts), and the radiology team from the CHU Dinant-Godinne UCL for flexibility and assistance. We are grateful to Dr J-J. Orban de Xivry and Pr Ph. Lefèvre for helpful discussion and comments about the study design. We thank G. Tincani (Hopital Saint Joseph, CHC, Liège) and A. Findik (Cliniques Universitaires Saint-Luc UCL, Brussels) for their help with patients screening and recruitment.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnhum.2015.00320/abstract
Beckmann, C. F., Deluca, M., Devlin, J. T., and Smith, S. M. (2005). Investigations into resting-state connectivity using independent component analysis. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 360, 1001–1013. doi: 10.1098/rstb.2005.1634
Bestmann, S., Swayne, O., Blankenburg, F., Ruff, C. C., Teo, J., Weiskopf, N., et al. (2010). The role of contralesional dorsal premotor cortex after stroke as studied with concurrent TMS-fMRI. J. Neurosci. 30, 11926–11937. doi: 10.1523/JNEUROSCI.5642-09.2010
Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409
Bosnell, R. A., Kincses, T., Stagg, C. J., Tomassini, V., Kischka, U., Jbabdi, S., et al. (2011). Motor practice promotes increased activity in brain regions structurally disconnected after subcortical stroke. Neurorehabil. Neural Repair 25, 607–616. doi: 10.1177/1545968311405675
Boyd, L. A., Vidoni, E. D., and Wessel, B. D. (2010). Motor learning after stroke: is skill acquisition a prerequisite for contralesional neuroplastic change? Neurosci. Lett. 482, 21–25. doi: 10.1016/j.neulet.2010.06.082
Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain's default network: anatomy, function, and relevance to disease. Ann. N.Y. Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011
Calautti, C., Jones, P. S., Naccarato, M., Sharma, N., Day, D. J., Bullmore, E. T., et al. (2010). The relationship between motor deficit and primary motor cortex hemispheric activation balance after stroke: longitudinal fMRI study. J. Neurol. Neurosurg. Psychiatry 81, 788–792. doi: 10.1136/jnnp.2009.190512
Carey, J. R., Kimberley, T. J., Lewis, S. M., Auerbach, E. J., Dorsey, L., Rundquist, P., et al. (2002). Analysis of fMRI and finger tracking training in subjects with chronic stroke. Brain 125, 773–788. doi: 10.1093/brain/awf091
Carey, L. M., and Seitz, R. J. (2007). Functional neuroimaging in stroke recovery and neurorehabilitation: conceptual issues and perspectives. Int. J. Stroke 2, 245–264. doi: 10.1111/j.1747-4949.2007.00164.x
Carter, A. R., Astafiev, S. V., Lang, C. E., Connor, L. T., Rengachary, J., Strube, M. J., et al. (2010). Resting interhemispheric functional magnetic resonance imaging connectivity predicts performance after stroke. Ann. Neurol. 67, 365–375. doi: 10.1002/ana.21905
Crawford, J. R., and Garthwaite, P. H. (2002). Investigation of the single case in neuropsychology: confidence limits on the abnormality of test scores and test score differences. Neuropsychologia 40, 1196–1208. doi: 10.1016/S0028-3932(01)00224-X
Dacosta-Aguayo, R., Grana, M., Savio, A., Fernandez-Andujar, M., Millan, M., Lopez-Cancio, E., et al. (2014). Prognostic value of changes in resting-state functional connectivity patterns in cognitive recovery after stroke: A 3T fMRI pilot study. Hum. Brain Mapp. 35, 3819–3831. doi: 10.1002/hbm.22439
Daselaar, S. M., Rombouts, S. A., Veltman, D. J., Raaijmakers, J. G., and Jonker, C. (2003). Similar network activated by young and old adults during the acquisition of a motor sequence. Neurobiol. Aging 24, 1013–1019. doi: 10.1016/S0197-4580(03)00030-7
Debas, K., Carrier, J., Orban, P., Barakat, M., Lungu, O., Vandewalle, G., et al. (2010). Brain plasticity related to the consolidation of motor sequence learning and motor adaptation. Proc. Natl. Acad. Sci. U.S.A. 107, 17839–17844. doi: 10.1073/pnas.1013176107
Desmurget, M., Epstein, C. M., Turner, R. S., Prablanc, C., Alexander, G. E., and Grafton, S. T. (1999). Role of the posterior parietal cortex in updating reaching movements to a visual target. Nat. Neurosci. 2, 563–567. doi: 10.1038/9219
Dipietro, L., Krebs, H. I., Volpe, B. T., Stein, J., Bever, C., Mernoff, S. T., et al. (2012). Learning, not adaptation, characterizes stroke motor recovery: evidence from kinematic changes induced by robot-assisted therapy in trained and untrained task in the same workspace. IEEE Trans. Neural Syst. Rehabil. Eng. 20, 48–57. doi: 10.1109/TNSRE.2011.2175008
Favre, I., Zeffiro, T. A., Detante, O., Krainik, A., Hommel, M., and Jaillard, A. (2014). Upper limb recovery after stroke is associated with ipsilesional primary motor cortical activity: a meta-analysis. Stroke 45, 1077–1083. doi: 10.1161/STROKEAHA.113.003168
Feydy, A., Carlier, R., Roby-Brami, A., Bussel, B., Cazalis, F., Pierot, L., et al. (2002). Longitudinal study of motor recovery after stroke: recruitment and focusing of brain activation. Stroke 33, 1610–1617. doi: 10.1161/01.STR.0000017100.68294.52
Fridman, E. A., Hanakawa, T., Chung, M., Hummel, F., Leiguarda, R. C., and Cohen, L. G. (2004). Reorganization of the human ipsilesional premotor cortex after stroke. Brain 127, 747–758. doi: 10.1093/brain/awh082
Ghilardi, M., Ghez, C., Dhawan, V., Moeller, J., Mentis, M., Nakamura, T., et al. (2000). Patterns of regional brain activation associated with different forms of motor learning. Brain Res. 871, 127–145. doi: 10.1016/S0006-8993(00)02365-9
Golestani, A. M., Tymchuk, S., Demchuk, A., and Goodyear, B. G. (2013). Longitudinal evaluation of resting-state FMRI after acute stroke with hemiparesis. Neurorehabil. Neural Repair 27, 153–163. doi: 10.1177/1545968312457827
Gomez Beldarrain, M., Grafman, J., Ruiz de Velasco, I., Pascual-Leone, A., and Garcia-Monco, C. (2002). Prefrontal lesions impair the implicit and explicit learning of sequences on visuomotor tasks. Exp. Brain Res. 142, 529–538. doi: 10.1007/s00221-001-0935-2
Hardwick, R. M., Rottschy, C., Miall, R. C., and Eickhoff, S. B. (2013). A quantitative meta-analysis and review of motor learning in the human brain. Neuroimage 67, 283–297. doi: 10.1016/j.neuroimage.2012.11.020
Heuninckx, S., Wenderoth, N., and Swinnen, S. P. (2008). Systems neuroplasticity in the aging brain: recruiting additional neural resources for successful motor performance in elderly persons. J. Neurosci. 28, 91–99. doi: 10.1523/JNEUROSCI.3300-07.2008
Jaillard, A., Martin, C. D., Garambois, K., Lebas, J. F., and Hommel, M. (2005). Vicarious function within the human primary motor cortex? A longitudinal fMRI stroke study. Brain 128, 1122–1138. doi: 10.1093/brain/awh456
Johansen-Berg, H., Rushworth, M. F., Bogdanovic, M. D., Kischka, U., Wimalaratna, S., and Matthews, P. M. (2002). The role of ipsilateral premotor cortex in hand movement after stroke. Proc. Natl. Acad. Sci. U.S.A. 99, 14518–14523. doi: 10.1073/pnas.222536799
Kantak, S. S., Stinear, J. W., Buch, E. R., and Cohen, L. G. (2012). Rewiring the brain: potential role of the premotor cortex in motor control, learning, and recovery of function following brain injury. Neurorehabil. Neural Repair 26, 282–292. doi: 10.1177/1545968311420845
Kleim, J. A., Chan, S., Pringle, E., Schallert, K., Procaccio, V., Jimenez, R., et al. (2006). BDNF val66met polymorphism is associated with modified experience-dependent plasticity in human motor cortex. Nat. Neurosci. 9, 735–737. doi: 10.1038/nn1699
Kwakkel, G., Kollen, B. J., Van Der Grond, J., and Prevo, A. J. (2003). Probability of regaining dexterity in the flaccid upper limb: impact of severity of paresis and time since onset in acute stroke. Stroke 34, 2181–2186. doi: 10.1161/01.STR.0000087172.16305.CD
Laird, A. R., Fox, P. M., Eickhoff, S. B., Turner, J. A., Ray, K. L., Mckay, D. R., et al. (2011). Behavioral interpretations of intrinsic connectivity networks. J. Cogn. Neurosci. 23, 4022–4037. doi: 10.1162/jocn_a_00077
Lefebvre, S., Dricot, L., Gradkowski, W., Laloux, P., and Vandermeeren, Y. (2012). Brain activations underlying different patterns of performance improvement during early motor skill learning. Neuroimage 62, 290–299. doi: 10.1016/j.neuroimage.2012.04.052
Lefebvre, S., Laloux, P., Peeters, A., Desfontaines, P., Jamart, J., and Vandermeeren, Y. (2013). Dual-tDCS enhances online motor skill learning and long-term retention in chronic stroke patients. Front. Hum. Neurosci. 6:343. doi: 10.3389/fnhum.2012.00343
Lindenberg, R., Renga, V., Zhu, L. L., Betzler, F., Alsop, D., and Schlaug, G. (2010). Structural integrity of corticospinal motor fibers predicts motor impairment in chronic stroke. Neurology 74, 280–287. doi: 10.1212/WNL.0b013e3181ccc6d9
Lotze, M., Markert, J., Sauseng, P., Hoppe, J., Plewnia, C., and Gerloff, C. (2006). The role of multiple contralesional motor areas for complex hand movements after internal capsular lesion. J. Neurosci. 26, 6096–6102. doi: 10.1523/JNEUROSCI.4564-05.2006
Mattay, V. S., Fera, F., Tessitore, A., Hariri, A. R., Das, S., Callicott, J. H., et al. (2002). Neurophysiological correlates of age-related changes in human motor function. Neurology 58, 630–635. doi: 10.1212/WNL.58.4.630
Matthews, P. M., Johansen-Berg, H., and Reddy, H. (2004). Non-invasive mapping of brain functions and brain recovery: applying lessons from cognitive neuroscience to neurorehabilitation. Restor. Neurol. Neurosci. 22, 245–260.
McHughen, S. A., Rodriguez, P. F., Kleim, J. A., Kleim, E. D., Marchal Crespo, L., Procaccio, V., et al. (2010). BDNF val66met polymorphism influences motor system function in the human brain. Cereb. Cortex 20, 1254–1262. doi: 10.1093/cercor/bhp189
Meehan, S. K., Randhawa, B., Wessel, B., and Boyd, L. A. (2011). Implicit sequence-specific motor learning after subcortical stroke is associated with increased prefrontal brain activations: an fMRI study. Hum. Brain Mapp. 32, 290–303. doi: 10.1002/hbm.21019
Ovadia-Caro, S., Villringer, K., Fiebach, J., Jungehulsing, G. J., Van Der Meer, E., Margulies, D. S., et al. (2013). Longitudinal effects of lesions on functional networks after stroke. J. Cereb. Blood Flow Metab. 33, 1279–1285. doi: 10.1038/jcbfm.2013.80
Pannek, K., Chalk, J. B., Finnigan, S., and Rose, S. E. (2009). Dynamic corticospinal white matter connectivity changes during stroke recovery: a diffusion tensor probabilistic tractography study. J. Magn. Reson. Imaging 29, 529–536. doi: 10.1002/jmri.21627
Park, C. H., Chang, W. H., Ohn, S. H., Kim, S. T., Bang, O. Y., Pascual-Leone, A., et al. (2011). Longitudinal changes of resting-state functional connectivity during motor recovery after stroke. Stroke 42, 1357–1362. doi: 10.1161/STROKEAHA.110.596155
Park, J. Y., Kim, Y. H., Chang, W. H., Park, C. H., Shin, Y. I., Kim, S. T., et al. (2014). Significance of longitudinal changes in the default-mode network for cognitive recovery after stroke. Eur. J. Neurosci. 40, 2715–2722. doi: 10.1111/ejn.12640
Penta, M., Tesio, L., Arnould, C., Zancan, A., and Thonnard, J. L. (2001). The ABILHAND questionnaire as a measure of manual ability in chronic stroke patients: rasch-based validation and relationship to upper limb impairment. Stroke 32, 1627–1634. doi: 10.1161/01.STR.32.7.1627
Pineiro, R., Pendlebury, S., Johansen-Berg, H., and Matthews, P. M. (2001). Functional MRI detects posterior shifts in primary sensorimotor cortex activation after stroke: evidence of local adaptive reorganization? Stroke 32, 1134–1139. doi: 10.1161/01.STR.32.5.1134
Reis, J., Schambra, H. M., Cohen, L. G., Buch, E. R., Fritsch, B., Zarahn, E., et al. (2009). Noninvasive cortical stimulation enhances motor skill acquisition over multiple days through an effect on consolidation. Proc. Natl. Acad. Sci. U.S.A. 106, 1590–1595. doi: 10.1073/pnas.0805413106
Schaechter, J. D., Fricker, Z. P., Perdue, K. L., Helmer, K. G., Vangel, M. G., Greve, D. N., et al. (2009). Microstructural status of ipsilesional and contralesional corticospinal tract correlates with motor skill in chronic stroke patients. Hum. Brain Mapp. 30, 3461–3474. doi: 10.1002/hbm.20770
Schaechter, J. D., and Perdue, K. L. (2008). Enhanced cortical activation in the contralesional hemisphere of chronic stroke patients in response to motor skill challenge. Cereb. Cortex 18, 638–647. doi: 10.1093/cercor/bhm096
Schulz, R., Park, C. H., Boudrias, M. H., Gerloff, C., Hummel, F. C., and Ward, N. S. (2012). Assessing the integrity of corticospinal pathways from primary and secondary cortical motor areas after stroke. Stroke 43, 2248–2251. doi: 10.1161/STROKEAHA.112.662619
Seeley, W. W., Menon, V., Schatzberg, A. F., Keller, J., Glover, G. H., Kenna, H., et al. (2007). Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356. doi: 10.1523/JNEUROSCI.5587-06.2007
Small, S. L., Hlustik, P., Noll, D. C., Genovese, C., and Solodkin, A. (2002). Cerebellar hemispheric activation ipsilateral to the paretic hand correlates with functional recovery after stroke. Brain 125, 1544–1557. doi: 10.1093/brain/awf148
Tombari, D., Loubinoux, I., Pariente, J., Gerdelat, A., Albucher, J. F., Tardy, J., et al. (2004). A longitudinal fMRI study: in recovering and then in clinically stable sub-cortical stroke patients. Neuroimage 23, 827–839. doi: 10.1016/j.neuroimage.2004.07.058
Torres, E. B., Quian Quiroga, R., Cui, H., and Buneo, C. A. (2013). Neural correlates of learning and trajectory planning in the posterior parietal cortex. Front. Integr. Neurosci. 7:39. doi: 10.3389/fnint.2013.00039
Vakil, E., Kahan, S., Huberman, M., and Osimani, A. (2000). Motor and non-motor sequence learning in patients with basal ganglia lesions: the case of serial reaction time (SRT). Neuropsychologia 38, 1–10. doi: 10.1016/S0028-3932(99)00058-5
Ward, N. S., Newton, J. M., Swayne, O. B., Lee, L., Thompson, A. J., Greenwood, R. J., et al. (2006). Motor system activation after subcortical stroke depends on corticospinal system integrity. Brain 129, 809–819. doi: 10.1093/brain/awl002
Werhahn, K. J., Conforto, A. B., Kadom, N., Hallett, M., and Cohen, L. G. (2003). Contribution of the ipsilateral motor cortex to recovery after chronic stroke. Ann. Neurol. 54, 464–472. doi: 10.1002/ana.10686
Westlake, K. P., Hinkley, L. B., Bucci, M., Guggisberg, A. G., Byl, N., Findlay, A. M., et al. (2012). Resting state alpha-band functional connectivity and recovery after stroke. Exp. Neurol. 237, 160–169. doi: 10.1016/j.expneurol.2012.06.020
Yin, D., Song, F., Xu, D., Peterson, B. S., Sun, L., Men, W., et al. (2012). Patterns in cortical connectivity for determining outcomes in hand function after subcortical stroke. PLoS ONE 7:e52727. doi: 10.1371/journal.pone.0052727
Keywords: motor skill learning, stroke, fMRI, neurorehabilitation, premotor cortex, resting-state fMRI, functional connectivity
Citation: Lefebvre S, Dricot L, Laloux P, Gradkowski W, Desfontaines P, Evrard F, Peeters A, Jamart J and Vandermeeren Y (2015) Neural substrates underlying motor skill learning in chronic hemiparetic stroke patients. Front. Hum. Neurosci. 9:320. doi: 10.3389/fnhum.2015.00320
Received: 05 February 2015; Accepted: 19 May 2015;
Published: 03 June 2015.
Edited by:Jean-Claude Baron, University of Cambridge, UK
Reviewed by:Pavel Lindberg, CNRS, INSERM, Université Paris Descartes, France
Assia Jaillard, University Hospital of Grenoble, France
Copyright © 2015 Lefebvre, Dricot, Laloux, Gradkowski, Desfontaines, Evrard, Peeters, Jamart and Vandermeeren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yves Vandermeeren, Department of Neurology, CHU Dinant-Godinne UCL Namur, Avenue Dr. G. Therasse, Yvoir, Belgium, firstname.lastname@example.org