Abstract
Introduction
The neuropathology of multiple sclerosis (MS) is characterized by multifocal inflammatory demyelination and neuroaxonal injury.1,2 Whereas remyelination is an important mechanism of repair after acute inflammatory demyelination, 3 clinical recovery is facilitated by adaptive functional reorganization.4-7 Although the occurrence and amount of lesional damage may trigger functional reorganization,6,8 the integrity of the extralesional brain tissue explains part of the potential for network plasticity underlying clinical recovery.9-11 An increasing burden of pathology limiting the availability of such tissue has been proposed to lead to the progression of disability in MS through a failure of brain plasticity.12,13
Our previous work, however, suggests that adaptive brain plasticity may be substantially preserved across levels of disability and damage in MS. 14 We showed that practice-related improvements in visuomotor performance of MS patients were present even in patients with higher MRI lesion volumes and greater disabilities and could not be distinguished from those of healthy controls. 14 This and related evidence challenge the hypothesis that increasing MS pathology leads to disability progression through a general failure of adaptive brain plasticity.14,15
Here, we investigate the relationship between microstructural damage and brain plasticity in patients with MS. We probe brain plasticity using functional magnetic resonance imaging (fMRI) in patients undergoing short- (minutes) and longer-term (days) practice of a visuomotor task. We first test the relationship between microstructural pathology, as reflected in brain MRI, and improvements in visuomotor performance with practice. We then characterize brain functional plasticity underlying these visuomotor improvements. We contrast patients with healthy controls to test the hypothesis that the capacity for performance improvement with practice is preserved in patients 14 despite high levels of microstructural damage. However, we expect that functional patterns characterizing the systems-level mechanisms of brain plasticity underlying this preserved capacity in patients differ from those of healthy controls. 16
Because motor practice induces functional changes of cognitive systems similar to those occurring during recovery after brain injury, 17 our investigation provides insights into mechanisms supporting motor recovery in MS. As visuomotor integration constitutes the basis for reaching, grasping, and locomotion, 18 our findings also inform about important functional substrates of rehabilitative interventions.
Methods
Participants and Study Design
A total of 35 right-handed individuals (23 MS patients and 12 healthy controls) participated in this study approved by the Oxfordshire Ethics Committee. Only patients who had been free of relapses for the 6 months prior to screening for the study were included. Patients were asked to keep their medications constant during the study.
Two phases of motor practice were assessed in this study: short-term practice was assessed over 12 minutes of the visuomotor task during the first functional MRI session (fMRI1); longer-term practice was tested over days of daily home practice with the same task. Participants were asked to train for at least 15 days, although shorter practice periods were not an exclusion criterion. Compliance was checked at the end of the study using the output files automatically saved on the practice laptop. At the end of the home training, participants underwent a second scanning session (fMRI2). In the separate scanning session, participants underwent a structural brain MRI scan.
Data Acquisition
Clinical measures
At baseline, patients and controls underwent a behavioral assessment of limb and cognitive functions using the 9-hole peg test (9-HP), the 25-foot walk, and the Paced Auditory Serial Addition Test of the Multiple Sclerosis Functional Composite. 19 In patients, clinical disability was quantified by the Expanded Disability Status Scale score by a trained assessor (VT). 20
Visuomotor performance
Visuomotor performance was evaluated using an isometric visuomotor tracking task.14,21 Participants were asked to track the vertical movements of a computer-controlled bar (the target bar) by altering the height of an adjacent bar controlled through varying the grip pressure applied to a handheld plastic rod. The target bar moved in a smoothly oscillatory fashion with the amplitude changing in a repetitive fashion to form the “sequence.” Continuous feedback on performance was provided by the color of the participant’s bar changing from red to blue when the difference between the heights of the 2 bars fell below 10 mm. Participants alternated 38-s sequence blocks with 38-s rest blocks, during which they made no response while watching a random sequence of movements generated to simulate their performance. A total of 7 blocks of each condition were performed during the fMRI sessions. Participants undertook a practice session before the start of the experiment. After fMRI1, they underwent a 13-minute daily training session at home with 10 blocks of sequence alternated with 10 blocks of rest.
Functional MRI
Data acquisition was performed on 3-T Varian Inova MRI system (Oxford, UK) using multislice gradient-echo EPI sequence [repetition time (TR) = 3000 ms; echo time (TE) = 30 ms; 43 axial slices (3.5-mm thick) providing whole-brain coverage; field of view (FOV), 256 × 192 mm2; matrix 64 × 64].
Structural MRI
The structural scans were performed on a 1.5-T Siemens Sonata scanner. T1-weighted images (T1-WIs) were acquired for brain volume calculation (TR = 12 ms, TE = 5.65 ms, flip angle = 19°, with elliptical sampling of k-space, giving a voxel size of 1 × 1 × 1 mm). Diffusion-WIs (DWIs) were acquired using echo planar imaging [60 axial slices (2.5-mm thick), matrix size 128 × 104, FOV = 320×260 mm, giving a voxel size of 2.5 × 2.5 × 2.5 mm]. Diffusion weighting was isotropically distributed along 60 directions using a
Data Analysis
Measures of visuomotor performance
The distance between the target bar and the pressure-sensing bar constituted the tracking error. The 90th percentile of the tracking error (p90) for sequence represented the summary measure of error across each block for each participant. The mean of p90 tracking errors across all blocks in each scanning session represented the overall tracking error for sequence. Short-term practice was quantified by a reduction of the p90 over 7 blocks of sequence during fMRI1. We used the slope of a linear fit to the individual tracking errors over 7 blocks of sequence as a measure of the rate of short-term practice-related improvements. The mean daily tracking error of the home sequence practice was used to quantify the longer-term practice-related improvements, the slope over days being the rate of longer-term practice-related improvements. All behavioral measures are reported in arbitrary units as mean ± standard error (SE).
Changes in mean tracking error were tested using repeated-measures analysis of variance (ANOVA) with group (controls and patients) as a between-subject factor and either block (1 to 7, for short-term data) or day (1 to 15, for longer-term data) as a within-subject factor.
Structural MRI analysis
Analysis was carried out using tools from the FMRIB Software Library (FSL) (www.fmrib.ox.ac.uk/fsl) to characterize the amount and distribution of damage in patients.
Details on lesion volume calculation are reported in the supplementary material (available at http://nnr.sagepub.com/supplemental). Brain parenchyma volumes were measured on T1-WIs using SIENAX,22,23 which measures normalized brain volume (NBV).24,25
All DWIs were corrected for head motion and eddy currents. Mean fractional anisotropy (FA) images were created using DTIFit.
26
With Tract-Based Spatial Statistics,
27
individual FA maps were nonlinearly aligned to a common FA template. The across-subject mean FA image was calculated and used to generate a “skeleton” of white matter (WM) tracts, which was thresholded at FA > 0.2. Individual-participant maximum FA values nearest to the mean FA skeleton were perpendicularly projected onto this skeleton for statistical analysis. Mean FA values were calculated for the participants within the skeleton. To test for voxelwise differences in FA values, we carried out permutation-based nonparametric testing.
28
Results were considered significant for
T1-WIs were analyzed with FSL-VBM30,31 incorporating nonlinear registration and correction for local expansion or contraction. To test for local differences between gray-matter (GM) volumes in patients versus controls, a voxelwise general linear model (GLM) was applied using permutation-based nonparametric testing, forming clusters with threshold-free cluster enhancement and testing clusters for significance at
Analysis of the behavioral-structural relationships
Spearman correlation tested the relationship between the 9-HP test and overall tracking error as well as the relationship between measures of performance and structural MRI measures. The relationship between performance changes with short- and longer-term practice and structural MRI measures was tested in a partial correlation, which included initial performance as a covariate. Two-tailed unpaired
For all the statistical tests, differences were considered significant at the
fMRI of motor practice
Analysis was carried out using tools from FSL. First-level (step 1) data preprocessing included motion correction, brain extraction, spatial smoothing (Gaussian kernel of 8-mm full width at half maximum), high-pass temporal (available at http://nnr.sagepub.com/supplemental) (150-s cutoff), and correction for field inhomogeneities through fieldmap-based EPI unwarping using FEAT of FSL. 26 Nonlinear registration from high-resolution T1 structural to MNI standard space was carried out. The time series was analyzed using a GLM approach with local autocorrelation correction. The canonical gamma variate hemodynamic response function was used. Two principal explanatory variables (along with their temporal derivatives) specified (1) the onset and duration of sequence task periods to identify the mean effect associated with the task and (2) a linear trend over the course of the scanning session in signal change associated with the sequence task, which was orthogonalized with respect to the mean effect of the task. Head motion parameters were added to the model as confound regressors.
Higher-level (group) analyses were conducted
32
with automatic outlier deweighting.
33
GM partial volume information based on the individual structural images was added to the model as a voxel-dependent covariate. Group
A mean within-session group effect associated with the sequence task as well as with sequence linear change was identified (step 2). Next, a mean difference between sessions was tested to identify changes associated with home practice. This was achieved by using a fixed-effects, between-session, within-subject analysis (step 3), followed by a mixed-effects, between-session, group-level analysis (step 4).
The relationship between variation in behavioral measures and variation in BOLD signal change was tested in correlation analyses. To identify short-term BOLD behavioral relationships, individual-participant outputs from analysis of the first fMRI session in step 1 were tested for correlation with de-meaned individual measures of (1) p90 tracking error during the first block of sequence, (2) slopes of tracking error over 7 blocks of sequence, and, in a separate design, (3) baseline overall motor performance, that is, mean of p90 tracking error, across 7 blocks of fMRI1. To identify longer-term BOLD-behavioral relationships, individual-participant outputs quantifying BOLD change between sessions from step 3 were tested for correlation with de-meaned individual measures of (1) mean performance during the first day of home practice and (2) individual slopes of tracking error changes over the home practice period. Number of days of practice was included as a confound regressor.
Brain functional activations were labeled using the Harvard-Oxford Structural Atlas, the Juelich Histological Atlas, and the Oxford Thalamic Connectivity Probability Atlas (http://www.fmrib.ox.ac.uk/fsl/data/atlas-descriptions.html).
fMRI of control tasks
To assess any nonspecific brain BOLD signal changes over time, participants performed control tasks. During the visual checkerboard task, they were asked to fixate on a cross in the middle of the screen. No motor response was required. During the motor task, they were asked to tap their right index finger at a frequency of 1 Hz. The visual checkerboard was used to cue the motor response. Both tasks used a 30-s block design (5 blocks of “on” alternated with 5 blocks of “off”). Analysis was performed with the same general schema as for the visuomotor task. However, high-pass temporal filtering with a 60-s cutoff was used to remove low-frequency drifts.
Results
Baseline Characteristics
The MS patients (19 relapsing-remitting and 4 secondary progressive) had lower manual dexterity but were well matched to the healthy controls for demographic and cognitive measures (Table 1). The Expanded Disability Status Scale functional systems had the following median (range) scores: pyramidal = 3 (0-4); sensory = 2 (0-4); cerebellar = 2 (0-4); brainstem = 1 (0-4); bladder/bowel = 1 (0-3); visual = 1 (0-4); mental = 1 (0-2).
Demographic, Clinical, and Structural Characteristics of Controls and Patients a
Abbreviations: MS, multiple sclerosis; EDSS, Expanded Disability Status Scale; 9-HPT, 9-hole peg test; 25-FW, 25-foot walk; PASAT, Paced Auditory Serial Addition Test; T2-LV, T2-hyperintense lesion volume; T1-LV, T1-hypointense lesion volume; NBV, normalized brain volume; FA, fractional anisotropy; MD, mean diffusivity; L1, axial diffusivity; L2L3, radial diffusivity; GM, gray matter.
Values are reported as mean ± standard error and the significance was determined using an unpaired
Pearson χ2.
Mann-Whitney
Mean of 2 consecutive trials.
Supplementary Figure S1 shows the distribution of WM lesions in the patient group. Patients had smaller NBVs than controls, and there was widespread microstructural damage with reduced GM volume and WM FA (Table 1, Figure 1).

A. Between-group differences in gray-matter (GM) volume (
Visuomotor Performance and Its Changes With Practice
At baseline, patients had higher overall tracking error for the visuomotor task (177.1 ± 11.1 units) than did the controls (134.6 ± 8.2 units;

(A), Individual measures of overall tracking error at baseline (fMRI1) and after longer-term practice (fMRI2); (B), individual slope of short-term tracking error during the first scanning session (fMRI1); and (C) individual slopes of longer-term tracking performance error with home practice in controls (triangle) and in patients (circle). Color figure available online at nnr.sagepub.com.
Both groups improved their performance during the first practice session (
The duration of practice was shorter than 15 days in 2 patients and in 2 controls. The length of home practice was similar for patients (median, 17 days; range, 9-24 days; only 1 patient with <2 weeks of practice) and controls (median, 17 days; range, 13-21 days; only 1 control with <2 weeks of practice) as were rates of performance improvement over the practice period (controls, −1.1 ± 0.2 units/d; patients, −1.3 ± 0.2 units/d;
Visuomotor Performance and MRI Measures of Brain Pathology
At baseline, the mean tracking error during block 1 of sequence correlated with the NBV and mean FA within the major WM tracts in patients (NBV: ρ = −0.60,
Visuomotor Performance and fMRI Measures of Brain Plasticity
Relating initial performance to brain activity at baseline
Execution of the visuomotor task and the control visual and motor tasks were associated with activation of bilateral cortical, subcortical, and cerebellar regions both in controls and in patients (Figure 3A). There was no significant between-group difference in either the visuomotor or the control task.

BOLD signal changes associated with (A) baseline sequence-related activation (sequence vs rest contrast); (B), performance improvements for short-term (during fMRI1); and (C) long-term practice (contrasting fMRI1 > fMRI2) for controls (upper panel) and for patients (lower panel);
Within the control group, higher mean tracking error correlated with lower task-related activation in posterior parietal, occipital, and parahippocampal cortices during the first practice session (fMRI1; Table 2). No significant correlation was identified for the patient group, although a positive correlation between overall motor performance and BOLD signal in the left hippocampus, right superior parietal lobule, and right primary motor cortex was suggested using a voxelwise threshold of
Correlations Between Task-Related Brain Activation (Sequence vs Rest) and (1) Baseline Overall Motor Performance, (2) Short-term Practice-Related Improvements, and (3) Longer-term Practice-Related Improvements in Visuomotor Tracking Error a
Abbreviations: SPL, superior parietal lobule; IPL, intraparietal lobule; PMd, lateral premotor cortex dorsal; S1, primary somatosensory cortex; M1, primary motor cortex; V, visual cortex; R, right; L, left.
Localization of clusters is in Montreal Neurological Institute (MNI) Standard Brain Space.

A. Relationship between individual variation in overall tracking error at baseline and BOLD percentage signal change during fMRI1 (
Changes in brain activity with visuomotor practice
Reductions in BOLD signal with practice were observed in the prefrontal, premotor, motor, and parietal cortices in controls over the course of fMRI1. Additional reductions in activation were found in the temporal and occipital cortices and in the cerebellum in the patients (Figure 3B). Increases in BOLD signal change were not observed in the first practice session for either group.
After longer-term practice (fMRI1 vs fMRI2), there were significant decreases in visuomotor task–related activation in the occipital cortex for controls and in both the occipital and parietal cortices for patients (Figure 3C). We did not find increases in BOLD signal change between sessions for either group. There were no significant between-session changes in BOLD signal for the control visual or motor tasks for either group.
Relating performance improvements to changes in brain activity
Brain regions relevant for visuomotor improvements differed between patients and controls. Greater performance improvements with short-term practice (fMRI1) were associated with smaller task-related mean BOLD signal in sensorimotor, premotor, cingulate, temporal, and parahippocampal cortices for the patient group (Table 2). In many of these same regions, the mean BOLD signal was more strongly correlated with performance improvements in patients than in controls (Table 2, Figure 4B). We found no correlation between the mean BOLD signal and performance improvements for controls using fully corrected statistics. However, a voxelwise threshold of
Greater improvements in performance after longer-term practice (fMRI1 vs fMRI2) were associated with smaller changes in BOLD signal in the occipital cortex for the control group (Table 2). This relationship was not found for patients, resulting in a significant between-group difference (Table 2, Figure 4C). However, a post hoc analysis using a voxelwise threshold of
Discussion
Our results demonstrate that performance improvements in MS patients can occur despite considerable brain damage and disability. Baseline patterns of activation for the visuomotor and control tasks suggest adaptive plasticity in patients. The comparison of activation changes with visuomotor task practice for the 2 study groups provides direct evidence for functional plasticity that is both qualitatively and quantitatively different between patients and controls, extending previous neurophysiological observations.15,16
Relating Brain Damage to Brain Plasticity in MS
Along with previous behavioral evidence for a substantial preservation of brain plasticity with increasing burden of disability, 14 our results challenge the hypothesis that accumulation of greater burdens of microstructural damage leads to clinical progression in MS through an exhaustion of plastic reserve. They suggest that the concept of the progression of disability in MS should be extended from the simple focus on outright and general failure of adaptive plasticity to the consequences of more specific functional impairments arising from pathology in critical pathways 6 combined with those of chronic limb disuse. 37 Although systems-level adaptive plasticity may compensate for damage, 8 this compensatory potential can be exceeded with progression of neuroaxonal loss affecting functional “bottlenecks” (eg, the major afferent optic nerves for vision, the efferent corticospinal tracts for movement). Substantial evidence has highlighted that the integrity of long pathways correlates functionally with measures of impairment38-40 and behaviorally with measures of disability.6,41 Therefore, although we show that brain plasticity is largely preserved even at higher burdens of microstructural damage, the ultimate behavioral impact of systems-level adaptive plasticity may be limited by the extent of chronic damage to functionally relevant long pathways.
Learned disuse of the limb 42 with disease progression may also contribute to progression of disability in MS. We found a significant correlation of WM integrity or NBV with performance at baseline but not with performance improvements in the patients. Simple relationships between lesion size and prognosis for recovery with neurorehabilitation have also been weak for patients after stroke. 43 Learned disuse has provided a compelling interpretation of these results and motivated successful neurorehabilitation approaches based on forced use of the affected limb.44,45 Our observations suggest that learned disuse, rather than pathology alone, could contribute to disability in MS through maladaptive mechanisms, 37 highlighting a specific new focus for neurorehabilitation interventions.44,46
Characterizing Adaptive Brain Plasticity in MS
Patients showed different functional patterns underlying performance improvements with practice compared with controls. Short-term practice-related improvements were associated with signal changes in regions for sensorimotor control and integrative aspects of task encoding 47 : lower levels of activity in these regions were found in patients who improved their performance more with practice. Evidence that plasticity in the primary sensorimotor cortex mediates the establishment of practice-related and task-specific motor memories comes from a broad range of studies in both animal models48,49 and healthy humans. 50 Previous studies in MS have suggested that recruitment of sensorimotor areas ipsilateral to the hand moved may either limit the behavioral impact of brain pathology8,16,51 or may represent a functional marker of more severe microstructural damage that is reflected in impaired performance.52-55 The localized relationship between greater visuomotor improvements and lower levels of activity in sensorimotor and premotor cortices, which also show reduced activity with practice, provides direct evidence for an adaptive role of the ipsilateral sensorimotor areas, whose activity is no longer required when performance has improved with practice. These results support the rationale for therapeutic modulation of the sensorimotor activity ipsilateral to the hand moved as part of neurorehabilitation strategies.56-58
Brain activation changes after home practice could be consistently related to performance improvements only for the controls, for whom a sustained activation in the primary visual cortex correlated with larger improvements in performance. This relationship may reflect an aspect of visual memory related to directed attention for the visuomotor task.59-61 Although an exploratory analysis suggests that longer-term performance changes may be associated with changes in activity in the parietal and lateral occipital cortices, the failure to identify consistent correlations in the patients is striking. In the control experiment, the visual response to the checkerboard was indistinguishable between patients and controls, demonstrating that the potential for activation of the visual cortex was not impaired in patients. These results could reflect smaller adaptive functional changes, previously proposed as a disease trait associated with loss of functional connectivity. 12 Alternatively, there may be a greater heterogeneity of patterns between the patients arising from differences in individual pathologies and the functional systems that are predominantly impaired as well as from differences in the individual potential for improvement revealed or enhanced by practice. Differences in patterns of activation with a simple motor task across different stages of MS suggested the latter hypothesis. 62 In our study, indirect evidence supporting this notion comes from the identification of practice-related signal reductions in a wider range of areas in patients than in controls.
There were no significant between-session changes in BOLD signal for the control visual or motor tasks for either group. This suggests that practice of the visuomotor task over this short period does not influence low-level visual or motor responses. This observation and the localization of functional changes observed with practice of the task together emphasize that the major modulatory effects of practice of this task are expressed in higher-order visuomotor control regions rather than in primary sensory and motor regions.
Potential Limitations of the Study
When comparing participants over a broad range of disabilities, the choice of an optimal task to probe visuomotor learning is problematic. We developed an isometric grip task to allow acquisition of a motor action sequence without need for fine finger control. We previously tested this task behaviorally in healthy controls 21 as well as in a separate cohort of neurological patients.14,63 We have previously demonstrated a relationship between task performance and increasing MS disability and lesional damage, 14 suggesting that the task reflects clinically relevant aspects of MS pathology. The impact of differences in impairment of hand function in the patients was limited by patient selection and study design: All patients were able to grip well and perform the 9-HP test, and each practice session had an initial calibration phase.
The home practice design potentially introduces inconsistencies in the way the participants practice. To limit this, each participant was trained in the use of the equipment before the home practice sessions started. We confirmed good compliance with the practice schedule for all participants by examination of the electronic records of each practice session.
There are also potential analysis confounds intrinsic to imaging-based comparisons of pathological and normal brains, which vary nonlinearly in size and shape. To limit type I error caused by GM volume loss in patients, our analysis approach included GM as a covariate, allowing us to control for differences in local anatomy, which can contribute to differences in BOLD signal generation. We also controlled the quality of our imaging preprocessing by ensuring corrections of distortions through field mapping before nonlinear registration of serial images. We took into account the potential confound of motion by including it as a regressor in the GLM analysis.
Finally, there is the potential for differential variability of fMRI responses between patients and controls. However, our previous studies of intersession variability have confirmed good reproducibility for fMRI measures both in controls and in patients. 64 Along with the findings from the control visual and motor tasks, this supports our interpretation of the fMRI changes as reflecting brain plasticity specific to the task.
Conclusions
The results of this study refute the hypothesis of a general failure of adaptive brain plasticity in disabled MS patients with diffuse microstructural damage. They suggest that the concept of the progression of disability in MS should be extended from the simple focus on general failure of adaptive plasticity to the combined consequences of functional impairments arising from pathology in critical pathways and from chronic limb disuse. In this view, our study not only provides a novel framework for neurorehabilitation strategies in MS,44,45 but it also encourages efforts to develop repair strategies 65 and neuroprosthetic approaches 66 for functionally critical pathways in disabled MS patients.
Footnotes
Acknowledgements
The authors wish to thank Dr Zsigmond T. Kincses and Mr Iain Wilson for their technical support with the experimental equipment. We are also grateful to Dr Mark Jenkinson for his support with correction of field distortions on functional data. We thank the people with MS who participated in the study.
Declaration of Conflicting Interests
The authors have no conflicts of interest to declare. PMM is a full-time employee of GlaxoSmithKline, which is engaged in the development of new drugs for multiple sclerosis.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by a grant of the Multiple Sclerosis Society UK (Grant ref. 829/05), a Training Fellowship of the Multiple Sclerosis Society Italy (VT), and the Oxford Biomedical Research Centre funded by the UK National Institute for Health Research (VT). The Wellcome Trust funds HJ-B. RGW is grateful to the UK Medical Research Council.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
