Keywords
Predictors, brain connectivity, robotic, motor function
Predictors, brain connectivity, robotic, motor function
The development of new methods for rehabilitation of deficits after stroke has enabled research into the brain mechanisms of improved function after such therapy. This has been accomplished in Constraint Induced Therapy1,2, Bilateral Arm Training3 and in one form of robotic hand training4. Some common themes in these studies include: 1. Changes in motor task-related brain activation after therapy (although both positive and negative changes have been reported) and, 2. Expansion of shrunken motor maps as measured by transcranial magnetic stimulation (TMS)5,6. However, there remain ambiguities and even controversies regarding the effects of repetitive task practice on brain activity and whether modern, well-defined therapeutic methods differ in their brain effects.
Robotic rehabilitation has certain mechanistic advantages over other therapeutic methods7. Robotic therapy is a better option for more severely affected stroke patients who may not be able to practice certain movements without external assistance. In such patients the mechanisms of recovery may be qualitatively different and who have the most to gain with improved understanding of the mechanisms of improvement of any particular therapy. In addition, although robotic therapy is well defined by the algorithms it uses for training, the therapy is flexible enough to train patients in various types of movements.
We had the opportunity to perform a multi-center investigation of the brain mechanisms underlying robotic rehabilitation by studying a subset of participants in two multi-center VA studies that compared robotic rehabilitation to both an intensity-matched non-robotic therapy regiment and usual care. The hypotheses for this study related to both prognosis (e.g. greater cortical motor excitability and reduction in transcallosal inhibition will predict greater functional improvement) and treatment effects (e.g. intensive rehabilitation will more effectively increase the ability to activate multiple muscles through motor cortical activity, partly through reduced interhemispheric inhibition to the affected motor cortex.) It was also an opportunity to test a connectivity-based approach that has shown promise in studies of recovery of function8,9.
This was a substudy of two multi-center clinical trials whose methods and results have been published10–12. It was originally intended to enroll approximately 40 participants across four sites but due to regulatory and staffing at issues at two sites, 13 subjects across two sites were enrolled. Briefly, all participants were chronic hemiparetic stroke patients with a significant degree of impairment (Upper Extremity Fugl-Meyer scale 7–38.) Figure 1 shows the basic design of the substudy, with TMS and magnetic resonance imaging (MRI) measures bracketing a 6–12 week intervention. Clinical Trial Registration URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00372411 & NCT00333983.
Stimulation of the motor cortex responsible for upper-extremity impairment was performed using a MagStim 200 or 2002 Magnetic Stimulator (MagStim Ltd., Wales, UK and 70 mm D double circular coil. Motor evoked potentials (MEP) were recorded unilaterally by surface electrodes fixed over the biceps and extensor digitorum communus (EDC) muscles in bipolar montage with 3 cm spacing. Responses were amplified by a battery-powered surface electromiography (EMG) integrated electrode and amplifier (B&L Engineering, Tustin, CA or DelSys, Boston, MA), and fed into a personal computer through a multifunctional I/O board and LabView acquisition/analysis software (National Instruments, Austin, TX). A 100 ms period after stimulus was examined with time window adjusted to capture only the MEP. Amplitudes were measured peak-to-peak. Bandpass was 30–1000 Hz and digitization 2000 Hz.
The target muscle was ensured to be at rest during the entire procedure, through audio and visual monitoring of muscle activity. Motor threshold was determined using International Federation of Clinical Neurophysiology criteria13 except that a 25 µV limit was used because of the bipolar montage. The coil was localized on the frontoparietal region contralateral to the target muscle in the examined limb and moved until each muscle’s hot-spot, where the response threshold was the lowest, was found. Exact position of stimulation was recorded using a stereotactic system (BrainSight, Rogue Research, Montreal, QC, Canada) and guided by a 1 cm Cartesian coordinate system projected onto the subject’s own MRI. If a different hotspot was found on a subsequent visit, threshold and recruitment curves were obtained at both the original and new hotspot, but the original hotspot data were used for group analysis.
Recruitment curves were measured by stimulating at a range of intensities, from 10% below threshold, increasing in increments of 10% of threshold until the response plateaued or the maximum output of the stimulator was reached. Ten stimuli at each intensity were delivered.
The ipsilateral silent period14 was measured by stimulation of the unaffected (defined here as contralesional) cortex and voluntary activation of the affected arm. The maximum force that the subject could sustain in each target muscle was determined. For the biceps the weight was placed on the wrist and for the extensor digitorum communis (EDC) on the proximal interphalangial joints. The coil was placed over the hand knob (hand representation within M1) of the unaffected hemisphere. After the subject stabilized the weight against gravity, a TMS pulse 100% of the maximal stimulator output was delivered. The subject was allowed to rest for several seconds and the procedure was repeated two further times. The EMG signal was integrated, and a ratio of post-relative to pre-stimulation activity was computed, making the appropriate adjustment for the length of period.
Anatomical and Functional Magnetic Resonance Imaging was performed at each center on a Tim Trio 3T scanner (Siemens AG, Erlangen, Germany) equipped with an 8-channel receive-only head coil. Anatomical Imaging: These consisted of a high-resolution three-dimensional sagittal T1-weighted magnetization-prepared rapid gradient echo (MP-RAGE) image, and oblique proton density and T2-weighted images acquired with 2 mm slice thickness. Diffusion-tensor images were acquired using a single-shot echo-planar technique and 65 directions. A b value of 1000 s/mm2 was used with an average of six images acquired to increase signal-to-noise ratio. A fractional anisotropy (FA) map was created from these data. A 5 mm radius spherical Region of Interest (ROI) was centered on the posterior limb of the left and right internal capsules (IC) on the FA images and the mean, standard deviation, and ratio (affected/unaffected) were computed.
Functional Imaging. Two eyes-closed rest scans were obtained, each with 128 coronal blood oxygenation-level dependent (BOLD) weighted volumes (echo planar imaging; 3 sec TR, 30 ms TE, 4 mm slice thickness with no gap, flip angle = 90°, 36 axial slices, 1.8 × 1.8 mm2 inplane resolution, FOV = 23 cm.) These were separated in time by at least 5 minutes. A tape and cushion technique was employed to reduce head motion and remind the subjects of the need to keep their head as still as possible. We examined head motion parameters within the analysis and rejected runs with absolute head movement greater than 2 voxels. Images were corrected for head motion by realignment and Independent Component Analysis (ICA) was used to remove movement related signal15.
Region of Interest (ROI) resting state correlation analysis. ROI based analysis was performed without spatial normalization in AFNI16 and MATLAB (MathWorks Inc., Natick, MA). All of a participant’s resting state scans were corrected for slice timing and spatially registered to the first resting sate scan from their first session. The structural image was skull-stripped and also spatially registered to the subject’s first functional scan. A 6-mm FWHM Gaussian blur was then applied to all spatially registered EPI scans. Nine ROIs were selected manually identifying the following anatomical landmarks: medial part of the precentral gyrus, postcentral gyrus, cerebellar hemispheres, supramarginal gyrus, supplementary motor area (caudal supplementary motor area between medial precentral gyrus and a coronal plane through the anterior commisure17 and superior, middle, and inferior frontal gyri, Figure 2. Pairwise ROI correlations were computed on the time series for each ROI and Z-transformed.
A: The 11 regions of interest (ROIs) are shown on representative axial slices of an example brain MRI. The top slices show the cerebellar ROI, then the PMAv on the left bottom slice, and PMAd, M1, and superior parietal regions from anterior to posterior. The SMA is represented by a single midline ROI. B. Correlation matrix with correlations at baseline in two resting states scans in the same participant. C. Example correlation in a single slice with a affected side M1 ROI as the seed.
SAS (Cary, North Carolina) was used for all analysis. Pearson correlation was used to calculate recruitment curve slope and associations between variables. Mixed models (REML with compound symmetry) were used to analyze predictive factors such as presence of TMS responses and recruitment curve slope.
Fourteen subjects were enrolled between 2008 and 2010 at the Baltimore and North Florida/South Georgia VAMC. Thirteen were stroke patients, three of whom were randomized to intensive comparison therapy and ten of which were randomized to robot therapy, and one subject was a healthy control who received no therapy. The relationship between initial and follow-up Fugl-Meyer (FM) impairment score and presence of MEP is shown in Figure 3. While MEP were absent in most of the lower functioning participants initially, there were both low and high functioning participants with absent and present MEP. Controlling for the effects of baseline FM in a fixed effects analysis, presence of an MEP at baseline was associated with a mean 3.3 ± 6.2 S.E. (N.S.) higher change in FM across all post-baseline visits. (There were up to four post-baseline visits.) A biceps MEP was never present without an EDC MEP, but not vice versa and the predictive value of these two measures was approximately equal.
The baseline and change in Fugl-Meyer score are shown for each participant, grouped by whether there was initial motor evoked potential as measured by TMS. Change 1 is across the intervention; Change 2 is between the end of the intervention and the last outcome measurement (12 weeks). Negative changes are always shown below the baseline.
One of the hypotheses regarding recruitment curves was that steeper recruitment curves would correlate with better function. However, there were no significant correlations between either EDC or biceps recruitment curves and function, within the population in which recruitment curves could be evaluated (N=6). Almost all participants had very shallow recruitment curves, with one moderately affected individual (FM=35) being the exception (Figure 4). There was a non-significant trend for higher recruitment curve slope at baseline correlating with functional improvement in a mixed effects model that controlled for baseline FM.
A. Single participant (#1) recruitment curves in the EDC and the second baseline and two follow-up measurements. Stimulation strength is indicated on the x-axis as a percentage of resting motor threshold. B. The slope of the recruitment curve between 100 and 120% of resting motor threshold stimulation strength was extracted from each recruitment curve of each participant that had measurable recruitment curves that could be measured on the affected side EDC. The first two measurements are both baseline periods. Changes in the recruitment curve slope were not significant. C. Recruitment curve slope in biceps (otherwise, as in B).
A long-lasting stimulus artifact contaminated too many cases to allow group analysis. One example of change in silent period is noted in Figure 5. In this case the subject had an increase in voluntary activity after the intervention, despite the same amount of force requirement, and demonstrated a clearer iSP only after the intervention. However, in most other subjects, there was no visible iSP.
An example of ipsilateral silent period measured at baseline (left) and immediately after the 12 weeks intervention (right) in participant RB5. EMG is measured in mV in the EDC muscle contralateral to the TMS stimulator in the upper trace and ipsilateral in the lower trace, which is offset by 2 mV. Before the intervention there is little activation of the muscle and also no apparent silent period. There is much more activation of the muscle after the intervention and also a visible short ipsilateral silent period.
Predictive measures. Resting state analysis resulted in a correlation matrix for the chosen ROI. While there was not an age-matched control population, there were clear asymmetries in the correlation matrices, as within hemisphere connections were decreased on the affected side as compared to the unaffected, but with exceptions such as the parietal area (data not shown). We were particularly interested in exploring the changes of correlation of the affected side motor cortex (M1) with other brain areas. Correlation of the affected M1 with all frontal lobe motor regions increased over the course of treatment, but there was no change in correlation of M1 with the cerebellum (Figure 6). The change in the unaffected M1, SMA and the unaffected side superior parietal area were most striking.
The mean change and S.D. of the Z-transformed correlation coefficient of the affected primary motor cortex (AM1) with each of the other regions is shown. Note that all regions showed an increase in connectivity except the parietal regions. Region names include ‘A’ for affected side hemisphere (the side opposite to the affected hemisphere in the case of the cerebellum, ‘CER’) and ‘U’ for unaffected.
Correlative measures. Improvement in Fugl-Meyer correlated with a trend towards reduction in two pairwise correlations. Greater FM increase was associated with decrease in affected M1-unaffected M1 connectivity (r2 = 0.31, p = 0.07), and unaffected M1-affected superior parietal area (r2 = 0.34, p = 0.06.) All other changes in functional connectivity measures correlated less well with changes in connectivity.
FA of the affected internal capsule was correlated with baseline motor ability (r2 = 0.48, p < 0.01.) but did not predict motor recovery, although the trend was for greater FA to be associated with better recovery.
There were no significant differences in any outcome measure for robotic vs non-robotic comparison therapy.
The purpose of this study was to obtain feasibility data on the use of TMS and MRI to provide predictive and mechanistic information about the motor functional response to intensive arm rehabilitation. It was not expected to provide definitive results in a field that has been marked by inconsistency. Some of the lessons learned in this study are that the lack of TMS responses in a majority of the moderate-to-severe population limits the utility of TMS for measuring change, although when MEP are present this predicts a better response to intervention, as has been demonstrated previously24. The ipsilateral silent period, a measure of transcallosal inhibition that can be performed even when MEP cannot be elicited, has limitations as well, and was not useful in this particular study, partly for technical reasons. MRI measures of resting state connectivity were more revealing, demonstrating both the deficits and changes with therapy, although in a purely exploratory manner.
While there were technical limitations to use of silent periods, as shown in Figure 5, the silent period could become more apparent after therapy. Since the appearance of a silent period depends on cortical activation of an affected muscle, a silent period could appear to be absent if there is little such cortical activation. Increased corticomotor effectiveness can thus result in the appearance of a silent period, and give a misleading impression that more intercortical inhibition is related to better function. The role of intercortical inhibition in shaping motor function is complex, and interpretation of tests that require activity in interhemispheric networks needs to be sophisticated.
In the resting state connectivity analysis, all connectivity with the affected motor cortex was negatively correlated with impairment. The two exceptions were the superior parietal area, in which increased connectivity was correlated with impairment. There have been a number of reports of the role of the superior parietal area in recovery of function after stroke18–22. However, its correlation here would suggest an association with more impairment and a role in compensation in only the more severely affected individuals.
The correlation matrix for even a small number of ROIs is a large amount of data and causes a multiple comparison problem. We focused on connectivity with the affected motor cortex as being most relevant to recovery of function. Out of the ten other regions, correlation with eight of them increased over the course of therapy. The largest increases were in the contralesional superior parietal area, ipsilesional dorsal premotor area, and supplementary motor area. These regions have strong bilateral connections and are good candidates for brain regions that would be engaged by the practice involving visual motor feedback and proximal arm movements. There were no significant associations between change in interregional correlation and a change in motor function. The best correlation with recovery was in connections of the unaffected M1 with two areas in the affected hemisphere: affected M1 and affected superior parietal area. The fact that this was a negative correlation suggests that intensive unimanual therapy may be decreasing the importance of the unaffected primary motor area in movement of the affected side. But other changes, not measured in this study, may be related to recovery of function, and the measured network changes, whether or not they are a significant effect of the intervention, may not be necessary for recovery or may represent compensatory changes. Likely because of the relatively small size of the study, were not able to find significant differences between such changes in the two types of treatment, if any differences exist. It would be interesting to speculate that the superior parietal activity would be more involved in robotic rehabilitation, with its visuomotor component and SMA in the intensive comparison treatment that involved more self-initiated activity.
Subject numbers started at 1 in Baltimore, and at 50 in Florida. All had anterior circulation ischemic strokes except for #1 who had a thalamic hemorrhage. Therapy assignment included either robotic or intensive comparison (comp.) therapy. When length of therapy was 6 weeks, there was no therapy with the wrist robot, only planar and vertical robots.
Measurement of brain changes related to motor recovery in moderate-to-severely affected stroke patients is complicated by difficulties in measuring brain function noninvasively. But our study showed that simple MEP presence might be useful in predicting response to rehabilitation in chronic stroke, while resting state connectivity appears to be responsive to treatment, with increase in affected primary motor cortical connectivity to other frontal motor areas. Motor cortical functional connectivity with the superior parietal cortex may be marker for compensatory changes that do not respond to affected side intensive practice.
F1000Research: Dataset 1. Raw data for predictors and brain connectivity changes associated with arm and motor function improvement from intensive robotic practice in cronic strocke, 10.5256/f1000research.8603.d13317525
GFW and ACL conceived of the study, SRR and RPG did MRI protocol design and analysis, PDG and SY statistical analysis, GFW, LGR, LMJL and LGR performed TMS studies.
This study was supported by a VA Rehabilitation Research and Development Merit Review award to GF Wittenberg.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors would like to thank Jaime Lush for expert technical assistance and Jui Panda for preliminary analysis. Dr. Jodie Haselkorn and Dr. Skip Rodriguez were involved in setting up a third site but regulatory issues prevented enrolling participants.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Competing Interests: No competing interests were disclosed.
Competing Interests: No competing interests were disclosed.
Competing Interests: No competing interests were disclosed.
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | |||
---|---|---|---|
1 | 2 | 3 | |
Version 2 (revision) 28 Feb 17 |
read | ||
Version 1 31 Aug 16 |
read | read | read |
Click here to access the data.
Spreadsheet data files may not format correctly if your computer is using different default delimiters (symbols used to separate values into separate cells) - a spreadsheet created in one region is sometimes misinterpreted by computers in other regions. You can change the regional settings on your computer so that the spreadsheet can be interpreted correctly.
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)