Original Research ARTICLE
Characterizing relationships of DTI, fMRI, and motor recovery in stroke rehabilitation utilizing brain-computer interface technology
- 1Department of Radiology, University of Wisconsin - Madison, Madison, WI, USA
- 2Department of Biomedical Engineering, University of Wisconsin - Madison, Madison, WI, USA
- 3Medical Scientist Training Program, University of Wisconsin School of Medicine and Public Health - Madison, Madison, WI, USA
- 4Neuroscience Training Program, University of Wisconsin School of Medicine and Public Health - Madison, Madison, WI, USA
- 5Departments of Orthopedics and Rehabilitation, University of Wisconsin - Madison, Madison, WI, USA
- 6Department of Kinesiology, University of Wisconsin - Madison, Madison, WI, USA
- 7Department of Medicine, University of Wisconsin - Madison, Madison, WI, USA
- 8Department of Neurology, University of Wisconsin - Madison, Madison, WI, USA
- 9Department of Neurosurgery, University of Wisconsin - Madison, Madison, WI, USA
- 10Department of Psychiatry, University of Wisconsin - Madison, Madison, WI, USA
- 11Department of Psychology, University of Wisconsin - Madison, Madison, WI, USA
The relationship of the structural integrity of white matter tracts and cortical activity to motor functional outcomes in stroke patients is of particular interest in understanding mechanisms of brain structural and functional changes while recovering from stroke. This study aims to probe these underlying mechanisms using diffusion tensor imaging (DTI) and fMRI measures. We examined the structural integrity of the posterior limb of the internal capsule (PLIC) using DTI and corticomotor activity using motor-task fMRI in stroke patients who completed up to 15 sessions of rehabilitation therapy using Brain-Computer Interface (BCI) technology. We hypothesized that (1) the structural integrity of PLIC and corticomotor activity are affected by stroke; (2) changes in structural integrity and corticomotor activity following BCI intervention are related to motor recovery; (3) there is a potential relationship between structural integrity and corticomotor activity. We found that (1) the ipsilesional PLIC showed significantly decreased fractional anisotropy (FA) values when compared to the contralesional PLIC; (2) lower ipsilesional PLIC-FA values were significantly associated with worse motor outcomes (i.e., ipsilesional PLIC-FA and motor outcomes were positively correlated.); (3) lower ipsilesional PLIC-FA values were significantly associated with greater ipsilesional corticomotor activity during impaired-finger-tapping-task fMRI (i.e., ipsilesional PLIC-FA and ipsilesional corticomotor activity were negatively correlated), with an overall bilateral pattern of corticomotor activity observed; and (4) baseline FA values predicted motor recovery assessed after BCI intervention. These findings suggest that (1) greater vs. lesser microstructural integrity of the ipsilesional PLIC may contribute toward better vs. poor motor recovery respectively in the stroke-affected limb and demand lesser vs. greater cortical activity respectively from the ipsilesional motor cortex; and that (2) PLIC-FA is a promising biomarker in tracking and predicting motor functional recovery in stroke patients receiving BCI intervention.
Studies have suggested that motor recovery after stroke is related to the structural remodeling of white matter tracts (Liu et al., 2008; Schaechter et al., 2009) and the reorganization of cortical activity (Dijkhuizen et al., 2001; Jaillard et al., 2005; Grefkes et al., 2008) in the ipsilesional and contralesional hemispheres. Little is known, however, about the relationship between the white matter structural integrity and functional cortical activity of the sensorimotor region and how these two factors interact with motor recovery in stroke patients. Therefore, a multimodal assessment of structure-function relationships may provide insights for examining factors influencing stroke recovery. Noninvasive brain imaging methods have been widely applied for understanding brain recovery following stroke. Diffusion tensor imaging (DTI) is one of these imaging methods, which allows for quantitative evaluations of the structural integrity of white matter tracts after a stroke (Werring et al., 2000; Stinear et al., 2007; Yu et al., 2009). DTI-derived measures have been shown as potential biomarkers used for tracking motor impairment (Schaechter et al., 2009; Lindenberg et al., 2010; Sterr et al., 2010; Yeo et al., 2010; Borich et al., 2012; Chen and Schlaug, 2013) and motor recovery (Liang et al., 2009) after stroke. DTI has also been investigated for its prognostic potential, with DTI measures assessed during the acute and sub-acute stages of stroke shown to predict motor impairments observed 1–7 months later (Cho et al., 2007; Koyama et al., 2012, 2013a,b; Groisser et al., 2014). A recent study shows evidence that DTI measures may be used as potentially predictive of individual recovery in stroke patients receiving newer neurorehabilitative therapies, such as transcranial direct current stimulation (Lindenberg et al., 2012). Besides DTI, fMRI is another non-invasive neuroimaging technique that has been used to gain better understanding of the processes of brain functional reorganization accompanying motor recovery after stroke (Calautti and Baron, 2003; Riecker et al., 2010; Garrison et al., 2013; Havsteen et al., 2013; Favre et al., 2014; Zhang et al., 2014). Measures derived from fMRI have also been shown as potential biomarkers to track recovery, with correlations between functional changes and fMRI measures demonstrated with treatments such as Brain-Computer Interface (BCI) therapy (Mukaino et al., 2014), constrain-induced movement therapy (Murayama et al., 2011; Kononen et al., 2012) and motor imagery therapy (Sun et al., 2013).
Fewer studies have taken a multimodal approach to characterize brain recovery after stroke by combining information from both DTI and fMRI measures, with a few attempts existing only as case studies (Jang et al., 2005; Caria et al., 2011). One recent study found that DTI-derived measures correlated more strongly with clinical outcomes than measures derived from fMRI (Qiu et al., 2011) and another study reported both DTI and fMRI derived measures correlated with motor outcomes (Chen and Schlaug, 2013). Correlations have also been identified between DTI and fMRI measures, with greater damage to white matter tracts showing an association with increased bilateral recruitment of motor areas and poorer motor performance in stroke patients (Wang et al., 2012), although there is also evidence that such correlations may be modulated after the functional electrical stimulation (FES) training (Wei et al., 2013).
In this study, one main goal is to investigate the relationship between DTI and fMRI measures and further investigate the relative contribution of each to the tracking and predicting of motor functional recovery in a group of stroke patients with persistent upper extremity impairment receiving BCI therapy. In the majority of stroke patients, the upper extremity is more severely involved than the lower limb, as most strokes occur in the territory of the middle cerebral artery (Shelton and Reding, 2001). Stroke that affects the posterior limb of the internal capsule (PLIC) has been reported to be significantly associated with poor recovery of isolated upper-limb movements (Shelton and Reding, 2001) and overall motor outcomes (Puig et al., 2011). Given the significance of PLIC involved in motor recovery, one specific aim of this study is to evaluate the stroke-induced changes in structural integrity of the PLIC using DTI fractional anisotropy (FA) and to investigate if these changes are related to motor recovery. Corticomotor activation is a fMRI way to examine the “integrity” of corticomotor functions. In this study corticomotor activity is evaluated using motor-task fMRI and quantified by counts of statistically significantly active voxels within the ipsilesional and contralesional motor cortices. Another specific aim of this study is to evaluate the changes in corticomotor activity, and to further examine if these changes are related to motor recovery. Combining both DTI and fMRI analysis, we examine the potential relationship between structural integrity of PLIC and functional integrity of motor cortex, and examine how this relationship interact with motor recovery in patients receiving BCI intervention.
Materials and Methods
A permuted-block design accounting for gender, stroke chronicity and severity of motor impairment was used to randomize patients to either a BCI intervention group or a crossover control group. Neuroimaging data and motor outcome assessments were acquired at four time points: before the start of intervention (i.e., pre-intervention), at the midpoint of intervention (i.e., mid-intervention), upon completion of intervention phase (i.e., immediately post-intervention), and 1 month following the last session of BCI intervention (i.e., 1-month-post-intervention). Patients in the BCI intervention group began to receive BCI intervention soon after recruitment. Patients in the control group first received three additional neuroimaging scans and motor outcome assessments during the control phase in which no BCI intervention was administered. These three additional assessments were acquired at intervals analogous to those administered during the BCI intervention phase. Upon completion the final-control neuroimaging and motor outcome assessment, these patients were crossed over to complete the BCI intervention phase. Table 1 illustrates the time frame of the study design. All current findings are based on neuroimaging and motor outcome measurements acquired from 9 patients during the BCI intervention phase.
Sixteen patients with persistent upper extremity motor impairment resulting from first-ever ischemic or hemorrhagic stroke were contacted regarding study participation in an on-going study investigating effects of EEG-BCI driven FES therapy of the impaired hand in stroke patients. This report is based on 9 patients who have completed the study (6 M, mean age of 61.9 years, chronicity of stroke range 2–23 months).
The inclusion criteria were: (1) ages 18 years and above; (2) no known neurologic, psychiatric or developmental disability; (3) persistent upper-extremity motor impairment resulting from ischemic or hemorrhagic stroke. The exclusion criteria were: (1) contraindications for MRI; (2) allergy to electrode gel, surgical tape and metals that would be used in BCI intervention; (3) under treatment for infectious disease or having apparent oral lesions or inflammation. This study was approved by the University of Wisconsin-Madison's Health Sciences Institutional Review Board. All patients provided written informed consent. Patient profiles are shown in Table 2.
Table 2. Patient profiles (age, gender, time since stroke, baseline NIHSS, baseline NIHSS-motor arm and stroke location).
BCI Intervention Procedures
All patients were administered up to 15 two-hour sessions of interventional BCI therapy (13.11 ± 2.20 sessions). Figure S1 illustrates a conceptual schematic of the system. A detailed description of the procedures followed during each session is provided in the Supplementary Material. These sessions took place over a period of up to 6 weeks with two to three intervention sessions per week.
Motor Functional Outcome Measures
All patients were assessed for clinical stroke severity in addition to neurologic examination at four time points throughout the intervention (Table 1). The neurologic deficit was evaluated on the basis of the severity of motor paresis using the National Institute of Health Stroke Scale (NIHSS) (Brott et al., 1989). All patients' motor function of the impaired arm was assessed using a neuropsychological battery which included objective measures such as the Action Research Arm Test (ARAT) (Carroll, 1965; Lang et al., 2006) and subjective measures such as the Stroke Impact Scale-Hand function domain (Duncan et al., 1999). The ARAT is a standardized measure of upper-limb functioning assessing grip, grasp, pinch and gross motor performance. Total ARAT scores ranged from 0 to 57. The Stroke Impact Scale (SIS) hand function subscale (SIS-Hand) was used to assess self-reported satisfaction with hand use and to evaluate the relationship between SIS-Hand scores and neuroimaging measures. Raw scores of SIS-Hand assessment were transformed using the following algorithm (Sullivan, 1995): Transformed scale = 100 × [(actual raw score − lowest possible raw score)/possible raw score]. The transformed scale of SIS-Hand score ranged from 0 to 100. All clinical assessments of the stroke-affected limb for each patient at each time point are shown in Table S3.
Neuroimaging data acquisition and processing are described in detail in the Supplementary Material. FA values were computed for the ipsilesional and contralesional PLIC. In addition, asymmetry indices between the ipsilesional and contralesional PLIC-FA (aFA) were calculated as aFA = (FAcontra − FAipsi)/(FAcontra + FAipsi) (Stinear et al., 2007; Schaechter et al., 2009; Lindenberg et al., 2010). This yields a value of aFA ranging from −1.0 to +1.0 with positive values indicating reduced FA in the ipsilesional PLIC, a value of 0 indicating symmetrical FA measurements from the two hemispheric PLIC, and a negative value indicating reduced FA in the contralesional PLIC.
All patients performed a block-design sequential finger tapping task during fMRI scans that consisted of alternating 20-s blocks of tapping vs. rest. Patients were cued to rest or to tap the fingers of one hand sequentially on a button box, using either visual or tactile (for visually impaired patients) cues. All patients underwent two fMRI scans using this paradigm –once when tapping with the impaired hand (passive tapping if unable to generate sufficient tapping independently) and again when tapping with the unimpaired hand. Patients were instructed to hold their heads still throughout the scans, and sufficient padding was provided to discourage head movement.
In the passive motor tasks, patients were assisted by the investigator in finger movements (flexion-extension) to complete the finger tapping tasks according to the experimental paradigm design.
Motor-Task Generated Active Voxels
A previously published mask of the cortical components of the motor network was used to identify statistically significantly active voxels in the motor cortex during finger tapping. This mask consisted of the cortical components of a previously identified motor network derived from an independent component analysis (ICA) of whole-brain resting-state fMRI (rs-fMRI) scans (Shirer et al., 2012). These independent components were visually selected based on previous reports and then thresholded independently and arbitrarily to generate distinct moderately sized ROIs in the cortex and subcortical gray matter (voxels = 25) (Shirer et al., 2012). It is worth noting that although rs-fMRI investigates synchronous activations between brain regions occurring in the absence of a task or stimulus, these synchronous activations have been observed in somatosensory, visual, attention and other higher-order brain areas and have shown close correspondence between the independent analyses of resting and activation brain dynamics (Biswal et al., 1995; Smith et al., 2009). In our study, only the cortical components from the motor network for each hemisphere were used for evaluation of corticomotor activity. Once fMRI data was processed (see Supplementary Material), this motor cortical mask was resampled into subject space and then applied to the functional data to identify those statistically significantly active voxels within the motor cortex using a threshold of t = 4 (p < 0.0001).
Considering the relatively small sample size (n = 9), we used non-parametric statistical tests for the analyses. The Wilcoxon signed-rank test was used to compare ipsilesional and contralesional PLIC FA values and to compare corticomotor fMRI activity between the two motor cortices. The Spearman rank correlation test was used for analyses of neuroimaging and motor outcome measurements. To take advantage of a longitudinal, repeated-measurement design of this study, we used generalized estimating equations (GEE) for regression analyses. GEE analyses use the generalized linear model to estimate more efficient and unbiased regression parameters relative to ordinary least squares regression (Ballinger, 2004). Most importantly, the GEE analyses take into account the dependency of repeated measurements from the same patient in the regression analysis. In addition, ANOVAs were used to examine how factors of time (pre-, mid-, immediately post-, and 1-month-post), PLIC (contralesional vs. ipsilesional side) and interaction between time and PLIC affecting DTI and fMRI measures. All statistical analyses were performed using RStudio (version 0.97.318). A p value less than or equal to 0.05 was considered statistically significant.
Patient Characteristics and Clinical Measures
Patient characteristics are summarized in Table 2. Average age was 61.89 years (SD = 9.25 years); average time from stroke onset was 13.11 months (SD = 7.90 months). There were no significant differences in terms of left or right hemisphere stroke (p = 0.14), cortical or non-cortical stroke (p = 0.36), or gender (p = 0.36). Standard clinical MRI was used to assess damage to PLIC by the neuroradiologist Dr. Prabhakaran. Six of the nine patients showed damage to PLIC due to stroke. Patient CT004 with a left middle cerebral artery (MCA) territory infarct showed minimal damage to PLIC. Patients CI003 with a small left frontal lobe infarct and CT003 with a right pontine infarct did not show damage to PLIC.
Clinical motor outcome measures are summarized in Table S3. The ARAT scores varied from zero, indicating no ability to perform, to a maximum of 57, indicating unimpaired performance. The SIS measure of hand function varied widely, with a value of zero indicating a patient reporting no ability to use the impaired hand, and higher positive values indicating decreasing levels of difficulty using the impaired hand in daily activities such as carrying heavy objects, turning a doorknob, opening a can or jar, tying a shoe lace and picking up a dime.
Relationship of PLIC-DTI Measures and Motor Outcomes
An ANOVA was computed to examine the main effects of time and PLIC (contralesional vs. ipsilesional side) as well as time × PLIC interaction (Figure 1A). Only PLIC factor was significant (p = 5.56e-07), and this was further validated with a Wilcoxon signed-rank test. Ipsilesional PLIC FA values were significantly lower when compared to the contralesional side (Wilcoxon signed-rank test: p < 0.05) except at time point 3 (immediately post-intervention) with a p-value equal to 0.055 trending toward significance (Figure 1B).
Figure 1. FA values compared between the ipsilesional and contralesional sides of the PLIC. (A) PLIC-FA changes across time compared between two hemispheres using ANOVA tests. Contra, contralesional PLIC; Ipsi, ipsilesional PLIC. (B) Boxplots showed significantly lower FA in ipsilesional PLIC when compared to contralesional PLIC (Time of pre-, mid-, post-, and 1-month-post intervention are indicated as time-point 1, 2, 3, and 4, respectively). White boxes represent the contralesional side and gray boxes represent the ipsilesional side of PLIC. (Wilcoxon signed-rank tests; all p-values are listed).
To assess the relationship between PLIC-FA values and motor outcome measures, a Spearman rank correlation test was first performed on the longitudinal data acquired from all patients and from all time points. The results suggested that higher ARAT scores and higher SIS-Hand scores were significantly correlated with higher FA values in ipsilesional PLIC (Figure 2). PLIC FA asymmetry (aFA) was negatively correlated with ARAT and SIS-Hand scores (Figure 2). A secondary statistical analysis, the GEE analysis, was further computed to control for the dependence of repeated measurements from each patient across time (Table 3). This analysis confirmed that the relationship observed between PLIC-FA and motor outcomes remained statistically significant. In addition, stroke severity as assessed by the NIHSS was significantly and negatively correlated with ipsilesional PLIC-FA values (Figure 3; GEE regression coefficient = −33.57, p-value = 3.69e-04).
Figure 2. Correlation analyses between ipsilesional FA, aFA and motor outcomes assessed in the impaired hand. Spearman rank correlation tests showed significant relationships between DTI and motor outcome measurements. (A) Ipsilesional PLIC-FA was positively correlated with SIS-hand function. (B) Ipsilesional PLIC-FA was positively correlated with ARAT. (C) PLIC-aFA was negatively correlated with SIS-hand function. (D) PLIC-aFA was negatively correlated with ARAT.
Figure 3. Correlation analyses on stroke severity (NIHSS scores) and structural integrity of PLIC (FA values). The Spearman rank correlation test (ρ = −0.695; p-value = 2.892e-05) and GEE regression analysis (GEE p-value = 3.691e-04) showed significant relationships between PLIC-FA values and patients' clinical severity assessments (NIHSS scores).
Prediction of Motor Function Recovery with Baseline Neuroimaging Measures
The linear regression analyses using a least-square fitting method revealed that baseline ipsilesional PLIC-FA correlated with the post-intervention ARAT and SIS-Hand scores (Figure 4; r-squared values > 0.7). Ipsilesional PLIC FA values measured at pre-intervention (baseline) were significantly and positively correlated with motor outcome scores measured immediately post- and 1-month-post intervention (p-value = 0.05). The same approach was applied to fMRI voxel counts, which did not reveal a predictive relationship with this fMRI measure on motor recovery.
Figure 4. Scatter plots of PLIC-FA predicting motor outcomes evaluated (A) immediately post-intervention (M3) and (B) 1-month-post intervention (M4). Linear regression analyses demonstrated that baseline PLIC-FA can predict motor outcomes. Pearson's correlation coefficients (ρ) and r-squared (R2) values are shown on the figures.
Relationship of Corticomotor Activity and Motor Outcomes
For the impaired finger tapping task, counts of active voxels were not significantly different between the two motor cortices (Wilcoxon signed-rank test: p = 0.28) (Figure 5A). Note, passive motor-task fMRI data was collected for patients CI001, CI002, CI004, CT001, and CT003 who were unable to perform motor tasks during fMRI scans. For the unimpaired finger tapping task, active voxel counts were significantly greater within contralesional motor cortex compared with the ipsilesional side (Wilcoxon signed-rank test: p = 0.038) (Figure 5B). An ANOVA was computed to examine the effects of time and PLIC (contralesional vs ipsilesional side) and the time × lesion interaction for both fMRI measures from impaired and unimpaired finger tapping. The results revealed no significant changes in corticomotor activity due to any of these factors (Figures 5C,D). However, the influence of PLIC trended toward significance (p = 0.064) for the unimpaired finger tapping task.
Figure 5. Comparison of active voxel counts during impaired and unimpaired finger tapping during fMRI scans. (A) Impaired finger tapping condition examining lesion factor, (B) Unimpaired finger tapping condition examining lesion factor, (C) Impaired finger tapping condition examining both lesion and time factors, (D) Unimpaired finger tapping condition examining both lesion and time factors.
Given the changes observed in corticomotor activity across time for impaired finger tapping compared to unimpaired finger tapping, another ANOVA was computed to test the effects of time (pre-, mid-, immediately post- and 1-month-post), PLIC (contralesional vs. ipsilesional), hand-impairment (impaired vs. unimpaired), and interaction between time and hand-impairment. We found that both the effect of hand-impairment and interaction between time and hand-impairment were significant (p-value = 0.005 and 0.015 respectively).
The Spearman rank correlation tests also demonstrated that fMRI measures (i.e., active voxel counts in ipsilesional motor cortex generated from impaired finger tapping task) were associated with motor outcomes. This negative association between motor outcomes and fMRI measurements may suggest that better motor outcomes after BCI-intervention are associated with a smaller number of active voxels within the ipsilesional motor cortex (Figure S2). However, this relationship was no longer significant after GEE regression analyses accounted for the dependence of repeated measurements from each patient (Table 3).
Relationship Between DTI and fMRI Measurements
The corticospinal pathway from the primary motor cortex through the PLIC approaches the midline of lower medulla oblongata and crosses to the contralateral side at the pyramidal decussation (Vulliemoz et al., 2005). A small percentage of fibers (10–25%), however, remain ipsilateral (Davidoff, 1990). A schematic illustration of the PLIC and motor cortex is shown in Figure 6. Taking this fact into account, we examined the DTI-fMRI relationship within each hemisphere for both impaired and unimpaired finger tapping tasks. The Spearman rank correlation test followed by GEE was used for correlation analyses.
Figure 6. A schematic illustration of the PLIC and motor cortex. A significant negative correlation was observed between ipsilesional PLIC-FA and ipsilesional motor activity.
For the impaired finger tapping, crossing fibers form the majority of ipsilesional PLIC, and may be affected by the small percentage of PLIC fibers on the contralesional side (Vulliemoz et al., 2005). We found that ipsilesional PLIC-FA negatively correlated with active voxel counts within the ipsilesional motor cortex (Figure 7). We did not, however, observe a significant relationship between contralesional PLIC-FA and voxel counts within the contralesional motor cortex. A similar approach was applied for measurements from unimpaired finger tapping. No significant relationship was found between PLIC-FA and voxel counts within ipsilesional or contralesional hemisphere during unimpaired finger tapping.
Figure 7. Relationship between ipsilesional corticomotor activity and ipsilesional PLIC-FA. Shown here are Spearman rank correlation coefficient (ρ) and corresponding p-value as well as GEE p-value.
Much of the focus of fMRI in stroke studies has been on whether its application provides a better understanding of brain functional reorganization accompanying motor recovery after stroke (Riecker et al., 2010; Garrison et al., 2013; Havsteen et al., 2013; Heiss and Kidwell, 2014). In a recent study, fMRI-derived measures have been correlated with movement recovery achieved with robot-assisted BCI therapy (Varkuti et al., 2013), while another study found that fMRI measures do not contribute significantly to the prediction of motor recovery (Zarahn et al., 2011). Another non-invasive MRI-based technique, DTI, has been widely used to evaluate the integrity of the white matter tracts after stroke. DTI-derived measures have been shown to be potential biomarkers for tracking motor impairment (Lindenberg et al., 2010; Chen and Schlaug, 2013) and motor recovery (Liang et al., 2009) after stroke. This study investigated the relationship between white matter integrity evaluated by DTI FA at the PLIC and corticomotor activity measured by motor-task fMRI, and further examined if these imaging measurements correlate with motor functional recovery in stroke patients receiving a BCI-facilitated intervention. Although our findings are preliminary and based on a moderate-size dataset, we observed consistent and robust results which are discussed here.
Structural Integrity of the PLIC vs. Motor Recovery
Previous human and animal studies have characterized changes in FA in the corticospinal system due to stroke (Liu et al., 2007; Kusano et al., 2009; Schaechter et al., 2009; Lindenberg et al., 2012). In our study, DTI analyses on 9 stroke patients with varying lesion locations and size of infarct affected the corticospinal system and yielded consistent observations, specifically, decreased FA in the ipsilesional PLIC compared to the contralesional side. This has been suggested as a characteristic of chronic white matter “Wallerian” degeneration (Yu et al., 2009; Lindenberg et al., 2012) and is thought to arise from the loss of tissue structural integrity (Liu et al., 2007).
In the current study, we found that higher ARAT and SIS-Hand scores were significantly correlated with higher FA values measured within ipsilesional PLIC after accounting for repeated measurements (Table 3). To account for potential changes in contralesional PLIC, FA asymmetry between the ipsilesional and contralesional sides were calculated (Stinear et al., 2007; Lindenberg et al., 2010) and then correlated with motor outcome measurements. The relationship between FA asymmetry and motor outcomes was significant even when accounting for repeated measurements. Lower or near zero aFA values indicate better preserved integrity of the ipsilesional PLIC and were correlated with better motor outcomes. Given our observation that the ipsilesional PLIC had decreased FA compared with the contralesional side, these findings suggest that the more the ipsilesional PLIC-FA resembled the contralesional PLIC, the greater the potential for functional recovery in stroke-affected limb.
Predictive Value of the PLIC-FA in Motor Recovery
One novel finding of this study is the predictive value of PLIC-FA in predicting recovery of motor function. Baseline FA values were significantly correlated with motor outcomes measured after intervention both immediately post-intervention and 1-month-post intervention (Figure 4). This suggests PLIC-FA may be useful as a novel biomarker to predict upper-limb motor functional recovery of stroke patients receiving BCI interventions.
For impaired finger tapping, we observed a bilateral pattern of corticomotor activity (Figure 5A), while unimpaired finger tapping produced significantly more lateralized activity within the contralesional motor cortex in comparison to ipsilesional motor cortex (Figure 5B). When comparing corticomotor activity for impaired and unimpaired finger tapping, we observed a significant difference due to the factor of hand-impairment (ANOVA test, p-value = 0.005). Furthermore, the pattern of changes in voxel counts (i.e., interaction between factors of time and hand-impairment) was also significant between the two tapping conditions (ANOVA test, p-value = 0.015). These findings may suggest that a potential bilateral pattern of corticomotor activity evolved during motor recovery in these stroke patients receiving BCI intervention.
Corticomotor Activity vs. Structural Integrity of the PLIC
Our results are that greater stroke severity is significantly correlated with compromised PLIC in the ipsilesional hemisphere (Figure 3). Furthermore, the ipsilesional PLIC-FA values are negatively correlated with ipsilesional corticomotor activity during impaired finger tapping (Figure 7). Thus, a greater burden of cortical activation is placed on the ipsilesional motor cortex for impaired finger tapping in patients with greater stroke severity with an overall bilateral pattern of motor cortical involvement. The BCI intervention may contribute to the increased utilization of ipsilesional and contralesional motor cortex with more bilateral activity seen during impaired finger tapping, rather than the lateralized activity seen normally during unimpaired finger tapping (Figures 2A,B). These factors may ultimately place demands on both ipsilesional and contralesional motor cortices during stroke recovery. This is novel in comparison to previous studies that suggest successful therapeutic intervention produces restoration of motor function mediated by re-lateralization of motor cortical activation (Ward et al., 2003a,b; Saur et al., 2006; Schlaug et al., 2008).
The small sample size (n = 9) and the heterogeneity of stroke patients (Table 2) were the primary limitations of this study. Four of the 9 patients exhibited no or little improvement in functional recovery as assessed by clinical behavioral performance. Those patients were severely impaired and minimally able to perform the designed intervention tasks, resulting in a floor effect in some outcome measurements. While changes in fMRI and DTI measurements were observed across time, an ANOVA did not show these changes to be significant, which may also be due to small sample size and high between-patient variance. It is also worth noting that our current findings are preliminary and are based on a moderate-size dataset.
Another limitation of our study is the combined analysis of passive and active finger tapping tasks performed by these patients. Passive vs. active tasks may have different effects on corticomotor activity which was further described and discussed in Supplementary Material.
Current DTI techniques remain limited in their ability to untangle the mix of PLIC fibers from the ipsilateral and contralateral hemispheres as they descend along the corticospinal pathway. Although we used task fMRI with each hand to investigate the structure-function relationship of the PLIC and the motor cortex, the current study design does not allow us to separate the mixture of white matter tracts from ipsilateral and contralateral corticospinal pathways within the PLIC which constrains our current findings. Future studies may therefore need to be done utilizing high resolution DTI and tractography techniques to further investigate the relationship between structural and functional changes in stroke patients.
Jie Song assisted with subject recruitment, data collection, data analysis, and writing. Brittany M. Young assisted with subject recruitment, data collection, data analysis, and writing. Zack Nigogosyan assisted in data collection and data analysis. Leo M. Walton assisted with data collection. Veena A. Nair assisted with subject recruitment, data collection, data analysis, and writing. Scott W. Grogan assisted with data collection and data analysis. Mitchell E. Tyler provided TDU hardware and expertise. Dorothy Farrar-Edwards assisted with study design and outcome measure selection and interpretation. Kristin E. Caldera assisted with subject recruitment. Justin A. Sattin assisted with study design and subject recruitment. Justin C. Williams is one of two lead PI's on this project and supervised the technical and engineering aspects of the work. Vivek Prabhakaran is one of two lead PI's on this project and supervised the neuroimaging aspects of this work.
Conflict of Interest Statement
There is one patent pending on the closed-loop neurofeedback device used for the BCI-facilitated intervention administered in this study (Pending U.S. Patent Application No. 12/715,090). This patent was filed jointly by the two lead investigators Justin C. Williams and Vivek Prabhakaran. 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.
The authors thank all patients and their families for their participation. Thanks also to study coordinator Jennifer Swartz and the MRI technicians at the Wisconsin Institutes for Medical Research. The authors also thank Wolfgang Gaggl for his generous advice in DTI analyses and interpretation, Alejandro Munoz Del Rio for his generous advice in statistical analyses, Hui-Chun (Ruby) Chen and colleagues in Occupational Therapy, Svyatoslav Vergun, Ricardo Pizarro, and Christian La for their assistance with the intervention. Portions of this work were previously presented in poster or e-poster form at the SFN Annual Meeting 2012, the International Stroke Conference 2013, the ISMRM Annual Meeting 2013, and as a platform talk at the RSNA Annual Meeting 2013. This work was supported by NIH grants RC1MH090912-01, K23NS086852, T32GM008692, UL1TR000427, and T32EB011434, by a Coulter Translational Research Award, an American Heart Association Postdoctoral Fellow Research Award, UW Milwaukee-Madison Intercampus Grants, by funding from the UW Graduate School, and by Grants from the Shapiro Foundation and Foundation of ASNR award.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fneng.2014.00031/abstract
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
Borich, M. R., Mang, C., and Boyd, L. A. (2012). Both projection and commissural pathways are disrupted in individuals with chronic stroke: investigating microstructural white matter correlates of motor recovery. BMC Neurosci. 13:107. doi: 10.1186/1471-2202-13-107
Brott, T., Adams, H. P. Jr., Olinger, C. P., Marler, J. R., Barsan, W. G., Biller, J., et al. (1989). Measurements of acute cerebral infarction: a clinical examination scale. Stroke 20, 864–870. doi: 10.1161/01.STR.20.7.864
Caria, A., Weber, C., Brotz, D., Ramos, A., Ticini, L. F., Gharabaghi, A., et al. (2011). Chronic stroke recovery after combined BCI training and physiotherapy: a case report. Psychophysiology 48. 578–582. doi: 10.1111/j.1469-8986.2010.01117.x
Chen, J. L., and Schlaug, G. (2013). Resting state interhemispheric motor connectivity and white matter integrity correlate with motor impairment in chronic stroke. Front. Neurol. 4:178. doi: 10.3389/fneur.2013.00178
Cho, S. H., Kim, D. G., Kim, D. S., Kim, Y. H., Lee, C. H., and Jang, S. H. (2007). Motor outcome according to the integrity of the corticospinal tract determined by diffusion tensor tractography in the early stage of corona radiata infarct. Neurosci. Lett. 426, 123–127. doi: 10.1016/j.neulet.2007.08.049
Dijkhuizen, R. M., Ren, J. M., Mandeville, J. B., Wu, O. N., Ozdag, F. M., Moskowitz, M. A., et al. (2001). Functional magnetic resonance imaging of reorganization in rat brain after stroke. Proc. Natl. Acad. Sci. U.S.A. 98, 12766–12771. doi: 10.1073/pnas.231235598
Duncan, P. W., Wallace, D., Lai, S. M., Johnson, D., Embretson, S., and Laster, L. J. (1999). The stroke impact scale version 2.0. evaluation of reliability, validity, and sensitivity to change. Stroke 30, 2131–2340. doi: 10.1161/01.STR.30.10.2131
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
Garrison, K. A., Aziz-Zadeh, L., Wong, S. W., Liew, S. L., and Winstein, C. J. (2013). Modulating the motor system by action observation after stroke. Stroke 44, 2247–2253. doi: 10.1161/STROKEAHA.113.001105
Grefkes, C., Nowak, D. A., Eickhoff, S. B., Dafotakis, M., Kust, J., Karbe, H., et al. (2008). Cortical connectivity after subcortical stroke assessed with functional magnetic resonance imaging. Ann. Neurol. 63, 236–246. doi: 10.1002/Ana.21228
Groisser, B. N., Copen, W. A., Singhal, A. B., Hirai, K. K., and Schaechter, J. D. (2014). Corticospinal tract diffusion abnormalities early after stroke predict motor outcome. Neurorehabil. Neural Repair. doi: 10.1177/1545968314521896. [Epub ahead of print].
Havsteen, I., Madsen, K. H., Christensen, H., Christensen, A., and Siebner, H. R. (2013). Diagnostic approach to functional recovery: functional magnetic resonance imaging after stroke. Front. Neurol. Neurosci. 32, 9–25. doi: 10.1159/000346408
Jang, S. H., You, S. H., Kwon, Y. H., Hallett, M., Lee, M. Y., and Ahn, S. H. (2005). Cortical reorganization associated lower extremity motor recovery as evidenced by functional MRI and diffusion tensor tractography in a stroke patient. Restor. Neurol. Neurosci. 23, 325–329.
Kononen, M., Tarkka, I. M., Niskanen, E., Pihlajamaki, M., Mervaala, E., Pitkanen, K., et al. (2012). Functional MRI and motor behavioral changes obtained with constraint-induced movement therapy in chronic stroke. Eur. J. Neurol. 19, 578–586. doi: 10.1111/j.1468-1331.2011.03572.x
Koyama, T., Marumoto, K., Miyake, H., and Domen, K. (2013a). Relationship between diffusion tensor fractional anisotropy and motor outcome in patients with hemiparesis after corona radiata infarct. J. Stroke Cerebrovasc. Dis. 22, 1355–1360. doi: 10.1016/j.jstrokecerebrovasdis.2013.02.017
Koyama, T., Marumoto, K., Miyake, H., Ohmura, T., and Domen, K. (2013b). Relationship between diffusion-tensor fractional anisotropy and long-term outcome in patients with hemiparesis after intracerebral hemorrhage. Neurorehabilitation 32, 87–94. doi: 10.3233/Nre-130825
Koyama, T., Tsuji, M., Miyake, H., Ohmura, T., and Domen, K. (2012). Motor outcome for patients with acute intracerebral hemorrhage predicted using diffusion tensor imaging: an application of ordinal logistic modeling. J. Stroke Cerebrovasc. Dis. 21, 704–711. doi: 10.1016/j.jstrokecerebrovasdis.2011.03.004
Kusano, Y., Seguchi, T., Horiuchi, T., Kakizawa, Y., Kobayashi, T., Tanaka, Y., et al. (2009). Prediction of functional outcome in acute cerebral hemorrhage using diffusion tensor imaging at 3t: a prospective study. Am. J. Neuroradiol. 30, 1561–1565. doi: 10.3174/Ajnr.A1639
Lang, C. E., Wagner, J. M., Dromerick, A. W., and Edwards, D. F. (2006). Measurement of upper-extremity function early after stroke: properties of the action research arm test. Arch. Phys. Med. Rehabil. 87, 1605–1610. doi: 10.1016/j.apmr.2006.09.003
Liang, Z., Zeng, J., Zhang, C., Liu, S., Ling, X., Wang, F., et al. (2009). Progression of pathological changes in the middle cerebellar peduncle by diffusion tensor imaging correlates with lesser motor gains after pontine infarction. Neurorehabil. Neural Repair 23, 692–698. doi: 10.1177/1545968308331142
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
Lindenberg, R., Zhu, L. L., Ruber, T., and Schlaug, G. (2012). “Predicting functional motor potential in chronic stroke patients using diffusion tensor imaging”. Hum. Brain Mapp. 33, 1040–1051. doi: 10.1002/hbm.21266
Liu, Y. R., D'Arceuil, H. E., Westmoreland, S., He, J. L., Duggan, M., Gonzalez, R. G., et al. (2007). Serial diffusion tensor MRI after transient and permanent cerebral ischemia in nonhuman primates. Stroke 38, 138–145. doi: 10.1161/01.Str.0000252127.07428.9c
Liu, Z. W., Li, Y., Zhang, X. G., Savant-Bhonsale, S., and Chopp, M. (2008). Contralesional axonal remodeling of the corticospinal system in adult rats after stroke and bone marrow stromal cell treatment. Stroke 39, 2571–2577. doi: 10.1161/Strokeaha.107.511659
Mukaino, M., Ono, T., Shindo, K., Fujiwara, T., Ota, T., Kimura, A., et al. (2014). Efficacy of brain-computer interface-driven neuromuscular electrical stimulation for chronic paresis after stroke. J. Rehabil. Med. 46, 378–382. doi: 10.2340/16501977-1785
Murayama, T., Numata, K., Kawakami, T., Tosaka, T., Oga, M., Oka, N., et al. (2011). Changes in the brain activation balance in motor-related areas after constraint-induced movement therapy; a longitudinal fMRI study. Brain Inj. 25, 1047–1057. doi: 10.3109/02699052.2011.607785
Puig, J., Pedraza, S., Blasco, G., Daunis-i-Estadella, J., Prados, F., Remollo, S., et al. (2011). Acute damage to the posterior limb of the internal capsule on diffusion tensor tractography as an early imaging predictor of motor outcome after stroke. Am. J. Neuroradiol. 32, 857–863. doi: 10.3174/Ajnr.A2400
Qiu, M. G., Darling, W. G., Morecraft, R. J., Ni, C. C., Rajendra, J., and Butler, A. J. (2011). White matter integrity is a stronger predictor of motor function than bold response in patients with stroke. Neurorehabil. Neural Repair 25, 275–284. doi: 10.1177/1545968310389183
Riecker, A., Groschel, K., Ackermann, H., Schnaudigel, S., Kassubek, J., and Kastrup, A. (2010). The role of the unaffected hemisphere in motor recovery after stroke. Hum. Brain Mapp. 31, 1017–1029. doi: 10.1002/hbm.20914
Saur, D., Lange, R., Baumgaertner, A., Schraknepper, V., Willmes, K., Rijntjes, M., et al. (2006). Dynamics of language reorganization after stroke. Brain 129(Pt 6), 1371–1384. doi: 10.1093/brain/awl090
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
Shirer, W. R., Ryali, S., Rykhlevskaia, E., Menon, V., and Greicius, M. D. (2012). Decoding subject-driven cognitive states with whole-brain connectivity patterns. Cereb. Cortex 22, 158–165. doi: 10.1093/cercor/bhr099
Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Beckmann, C. E., et al. (2009). Correspondence of the brain's functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi: 10.1073/pnas.0905267106
Sterr, A., Shen, S., Szameitat, A. J., and Herron, K. A. (2010). The role of corticospinal tract damage in chronic motor recovery and neurorehabilitation:a pilot study. Neurorehabil. Neural Repair 24, 413–419. doi: 10.1177/1545968309348310
Stinear, C. M., Barber, P. A., Smale, P. R., Coxon, J. P., Fleming, M. K., and Byblow, W. D. (2007). Functional potential in chronic stroke patients depends on corticospinal tract integrity. Brain 130(Pt 1), 170–180. doi: 10.1093/brain/awl333
Sun, L. M., Yin, D. Z., Zhu, Y. L., Fan, M. X., Zang, L. L., Hu, Y., et al. (2013). Cortical reorganization after motor imagery training in chronic stroke patients with severe motor impairment: a longitudinal fMRI study. Neuroradiology 55, 913–925. doi: 10.1007/s00234-013-1188-z
Varkuti, B., Guan, C., Pan, Y. Z., Phua, K. S., Ang, K. K., Sitaram, C. W. K., et al. (2013). Resting state changes in functional connectivity correlate with movement recovery for bci and robot-assisted upper-extremity training after stroke. Neurorehabil. Neural Repair 27, 53–62. doi: 10.1177/1545968312445910
Wang, L. E., Tittgemeyer, M., Imperati, D., Diekhoff, S., Ameli, M., Fink, G. R., et al. (2012). Degeneration of corpus callosum and recovery of motor function after stroke: a multimodal magnetic resonance imaging study. Hum. Brain Mapp. 33, 2941–2956. doi: 10.1002/Hbm.21417
Ward, N. S., Brown, M. M., Thompson, A. J., and Frackowiak, R. S. (2003a). Neural correlates of motor recovery after stroke: a longitudinal fMRI study. Brain 126(Pt 11), 2476–2496. doi: 10.1093/brain/awg245
Ward, N. S., Brown, M. M., Thompson, A. J., and Frackowiak, R. S. (2003b). Neural correlates of outcome after stroke: a cross-sectional fMRI study. Brain 126(Pt 6), 1430–1448. doi: 10.1093/brain/awg145
Wei, W. J., Bai, L. J., Wang, J., Dai, R. W., Tong, R. K. Y., Zhang, Y. M., et al. (2013). A longitudinal study of hand motor recovery after sub-acute stroke: a study combined fMRI with diffusion tensor imaging. PLoS ONE 8:154. doi: 10.1371/journal.pone.0064154
Werring, D. J., Toosy, A. T., Clark, C. A., Parker, G. J. M., Barker, G. J., Miller, D. H., et al. (2000). Diffusion tensor imaging can detect and quantify corticospinal tract degeneration after stroke. J. Neurol. Neurosurg. Psychiatry 69, 269–272. doi: 10.1136/jnnp.69.2.269
Yeo, S. S., Choi, B. Y., Chang, C. H., and Jang, S. H. (2010). Transpontine connection fibers between corticospinal tracts in hemiparetic patients with intracerebral hemorrhage. Eur. Neurol. 63, 154–158. doi: 10.1159/000281900
Yu, C., Zhu, C., Zhang, Y., Chen, H., Qin, W., Wang, M., et al. (2009). A longitudinal diffusion tensor imaging study on wallerian degeneration of corticospinal tract after motor pathway stroke. Neuroimage 47, 451–458. doi: 10.1016/j.neuroimage.2009.04.066
Zarahn, E., Alon, L., Ryan, S. L., Lazar, R. M., Vry, M. S., Weiller, C., et al. (2011). Prediction of motor recovery using initial impairment and fMRI 48 h poststroke. Cereb. Cortex 21, 2712–2721. doi: 10.1093/cercor/bhr047
Keywords: DTI, FA, fMRI, motor recovery, stroke rehabilitation, BCI
Citation: Song J, Young BM, Nigogosyan Z, Walton LM, Nair VA, Grogan SW, Tyler ME, Farrar-Edwards D, Caldera KE, Sattin JA, Williams JC and Prabhakaran V (2014) Characterizing relationships of DTI, fMRI, and motor recovery in stroke rehabilitation utilizing brain-computer interface technology. Front. Neuroeng. 7:31. doi: 10.3389/fneng.2014.00031
Received: 16 April 2014; Accepted: 08 July 2014;
Published online: 29 July 2014.
Edited by:Jaime A. Pineda, University of California, San Diego, USA
Copyright © 2014 Song, Young, Nigogosyan, Walton, Nair, Grogan, Tyler, Farrar-Edwards, Caldera, Sattin, Williams and Prabhakaran. 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: Jie Song, Wisconsin Institutes for Medical Research, 1111 Highland Avenue, Madison, WI 53705, USA e-mail: firstname.lastname@example.org;
Vivek Prabhakaran, Department of Radiology, 600 Highland Avenue, Madison, WI 53711, USA e-mail: email@example.com