Yang-teng Fan
Ching-yi Wu2,3†
Ho-ling Liu
Keh-chung Lin
Yao-liang Chen- 1School of Occupational Therapy, College of Medicine, National Taiwan University and Division of Occupational Therapy, Department of Physical Medicine and Rehabilitation, National Taiwan University Hospital, Taipei, Taiwan
- 2Department of Occupational Therapy and Graduate Institute of Behavioral Sciences, College of Medicine, Chang Gung University, Taoyuan, Taiwan
- 3Healthy Aging Research Center, Chang Gung University, Taoyuan, Taiwan
- 4Department of Imaging Physics, Division of Diagnostic Imaging, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
- 5Department of Medical Imaging and Radiological Sciences, Chang Gung University, Taoyuan, Taiwan
- 6Department of Physical Medicine and Rehabilitation, Division of Occupational Therapy, National Taiwan University Hospital, Taipei, Taiwan
- 7Department of Diagnostic Radiology, Chang Gung Memorial Hospital, Keelung, Taiwan
- 8MRI Center, Chang Gung Memorial Hospital, Taoyuan, Taiwan
Most neuroimaging research in stroke rehabilitation mainly focuses on the neural mechanisms underlying the natural history of post-stroke recovery. However, connectivity mapping from resting-state fMRI is well suited for different neurological conditions and provides a promising method to explore plastic changes for treatment-induced recovery from stroke. We examined the changes in resting-state functional connectivity (RS-FC) of the ipsilesional primary motor cortex (M1) in 10 post-acute stroke patients before and immediately after 4 weeks of robot-assisted bilateral arm therapy (RBAT). Motor performance, functional use of the affected arm, and daily function improved in all participants. Reduced interhemispheric RS-FC between the ipsilesional and contralesional M1 (M1-M1) and the contralesional-lateralized connections were noted before treatment. In contrast, greater M1-M1 functional connectivity and disturbed resting-state networks were observed after RBAT relative to pre-treatment. Increased changes in M1-M1 RS-FC after RBAT were coupled with better motor and functional improvements. Mediation analysis showed the pre-to-post difference in M1-M1 RS-FC was a significant mediator for the relationship between motor and functional recovery. These results show neuroplastic changes and functional recoveries induced by RBAT in post-acute stroke survivors and suggest that interhemispheric functional connectivity in the motor cortex may be a neurobiological marker for recovery after stroke rehabilitation.
Introduction
Resting-state functional magnetic resonance imaging (RS-fMRI) has emerged as a powerful tool for evaluating intrinsic brain connectivity and regional interactions during wakeful rest (Friston, 1994; Raichle and Mintun, 2006; Van Essen et al., 2012). The advantage of RS-fMRI is that it can be performed without task demand or external input and is therefore particularly useful for investigating connectivity mapping in patients with different levels of neurological impairment (Park et al., 2011).
Resting state functional connectivity (RS-FC) reflects the strength of temporal synchrony of blood oxygen level—dependent signals between spatially remote brain regions (Biswal et al., 1995; Gusnard and Raichle, 2001) and has been applied in stroke survivors (Rehme and Grefkes, 2013). Widespread changes in network functional connectivity take place immediately after stroke and have important implications for recovery (Carter et al., 2012; Grefkes and Ward, 2014). RS-FC research shows that the functional connectivity between the ipsilesional and the contralesional primary sensorimotor cortex is significantly diminished at the early stage of stroke (Wang et al., 2010; Park et al., 2011; Golestani et al., 2013). Decreased functional connectivity with the ipsilesional primary motor cortex (M1) was also found in other brain regions such as the bilateral supplementary motor area (SMA), bilateral secondary somatosensory cortex, bilateral cerebellum, bilateral thalamus, contralesional premotor cortex, and contralesional posterior parietal cortex (Carter et al., 2010; Wang et al., 2010).
During the recovery process after stroke, the resting state interhemispheric connectivity between the bilateral M1 increases (Wang et al., 2010; Park et al., 2011; Golestani et al., 2013), which is associated with motor improvements (Carter et al., 2010; Varkuti et al., 2013; Xu et al., 2014) and can predict better upper limb (UL) motor improvements in the next 6 months after stroke (Park et al., 2011). Seminal research on treatment-induced RS-fMRI found that increases in functional connectivity between the sensorimotor areas were correlated with gains in the Fugl-Meyer Assessment (FMA) score and improvements in the performance of activities of daily living (ADL) (James et al., 2009; Varkuti et al., 2013; Young et al., 2014). These imply that RS-FC changes may modulate motor and functional recovery after stroke and provide preliminary neural evidence for movement and functional improvements after a rehabilitation intervention. Notably, the brain is primed to neurological recovery in the first 3 months after stroke (Kwakkel et al., 2004; Murphy and Corbett, 2009), and most of the existing treatment-induced RS-FC studies have paid scant attention to the post-acute stage of stroke.
Robot-assisted bilateral arm therapy (RBAT), one of the main approaches to bilateral arm training, has gained increasing popularity in stroke neurorehabilitation and shows promising results in improving recovery of post-stroke sensorimotor functions (Hesse et al., 2005; Lin et al., 2010; Basteris et al., 2014; De Santis et al., 2014; Hughes et al., 2015). Enhanced motor function of stroke patients after receiving RBAT was reflected by results on the FMA and the Wolf Motor Function Test (WMFT) at the end of treatment and/or at the 3-month follow-up (Liao et al., 2012; Wu et al., 2013; Hsieh et al., 2014). A few studies have examined the changes in interhemispheric ipsilesional and contralesional M1 (M1-M1) functional connectivity but focused on the administration of unilateral robot-assisted therapy. Results showed that interhemispheric M1-M1 connectivity was increased in patients receiving unilateral robot-assisted therapy compared with pre-treatment and that a greater improvement in the pre-to-post difference in M1-M1 RS-FC was associated with greater gains in motor function (Sergi et al., 2011; Saleh et al., 2012).
These findings on unilateral robot-assisted therapy imply that interhemispheric motor cortex connectivity during the resting state may be a potential marker of stroke recovery after rehabilitative therapies (Chen and Schlaug, 2013). Nevertheless, the connectivity patterns after unilateral robot-assisted therapy may be different from RBAT because these two approaches induce different effects on stroke recovery (Van Delden et al., 2012; Wu et al., 2013). To our best knowledge, no study has looked at the effect of bilateral robot-assisted training on the functional connectivity features of the involved neuronal network during the resting state. There are only task-based neuroimaging studies showing that bilateral training facilitates excitability of transcallosal projections from the ipsilesional to the contralesional M1 (Stinear et al., 2014) and induces greater activation in the ipsilesional motor-related areas (Waller et al., 2014) than unilateral training. Illuminating the functional connectivity patterns in more regions and at the resting state post-stroke RBAT is of crucial importance, because these patterns might be potential markers of stroke rehabilitation and influence outcomes after stroke.
Previous research reveals that the recovery of the UL motor impairments is associated with the performance of ADL in stroke survivors (Dromerick et al., 2006; Wei et al., 2011) and that motor scales, such as the FMA and the WMFT, have good predictive validity with the Functional Independence Measure (FIM) at post-stroke rehabilitation (Hsieh et al., 2009). In addition, brain-imaging research has found that interhemispheric functional coherence in the resting state, particularly between the primary motor cortices, affects cortical reorganization and functional recovery after stroke (Rehme and Grefkes, 2013; Grefkes and Ward, 2014). Although these findings suggest that the change of M1-M1 RS-FC is linked with motor and functional outcomes during stroke recovery, the possibility that these neurological and behavioral variables may be interrelated, especially post-stroke RBAT, has not been tested. Thus, if our data show associations between the changes in clinical outcomes and pre-to-post difference in RS-FC after intervention, we would use a mediation analysis, a post hoc strategy, to further explore whether intrinsic brain connections represent a mediator between motor and functional recovery.
In the present study, we used seed-voxel correlation mapping (Horwitz et al., 1998) to investigate the changes of interhemispheric RS-FC in stroke survivors receiving RBAT. In addition to clinical measures for motor and functional outcomes, the relationships between RS-FC changes and motor and daily functions were also examined. We hypothesized that such interhemispheric connectivity would increase, that motor and functional performance would improve from pre-treatment to the end of treatment, and that the increase in RS-FC would correlate with improvements in motor performance and functional independence. Further, we would conduct a mediation analysis to test whether the pre-to-post difference in interhemispheric RS-FC is a significant mediator for an association between motor and ADL outcomes.
Materials and Methods
This study was approved by the local Ethics Committee and conducted in accordance with the Declaration of Helsinki.
Participants
The study enrolled 10 participants (8 men; mean age, 52.7 ± 6.5 years) with post-acute stroke in the right (n = 7) or left (n = 3) hemisphere. All participants gave written informed consent for the study. The diagnosis of stroke was clinically confirmed by computed tomography scanning. The average time from stroke onset was 46.8 (standard deviation, 20.11) days.
The inclusion criteria were: (1) first episode of stroke in cortical or subcortical regions; (2) time since stroke of less than 3 months and more than 2 weeks; (3) no serious cognitive impairment (Mini Mental State Exam score > 24) (Teng and Chui, 1987); and (4) initial motor part of the FMA-upper limb (FMA-UL) score ranging from 24–52 (Fugl-Meyer et al., 1975). We chose participants with mild and moderate limb paresis because they may have beneficial results with robot-assisted training (Housman et al., 2009; Hsieh et al., 2012). Exclusion criteria were: (1) aphasia that might interfere with understanding instructions; (2) chronic inflammatory, autoimmune, and hematologic disorders; (3) intake of anti-inflammatory drugs; (4) other major health problems or poor physical condition that might limit participation; and (5) current participation in any other research. Participants’ background information is reported in Table 1.
Table 1. Demographic and clinical characteristics of right-handed study participants.
Procedure
Eligible participants received RBAT. Motor performance, functional use of the affected arm, ADL, and the RS-fMRI were assessed at pre-treatment and at the end of 4 weeks of RBAT by the same blinded rater. All assessments were done within 3 days of each timestamp. The blinded raters were trained to properly administer these measures, and 2 certified occupational therapists assessed the raters’ competence.
Intervention Protocols
The RBAT was implemented for 90 min/day, 5 days/week, for 4 weeks. All participants received 5 min of tone normalization for the arm at the beginning of therapy. Training was administrated during regularly scheduled occupational therapy, and all other routine interdisciplinary rehabilitation that did not focus on UL training was continued as usual.
Participants received 70 min of RBAT. They practiced with the Bi-Manu-Track (Reha-Stim Co, Berlin, Germany), which includes wrist flexion-extension and forearm pronation-supination movements. During the Bi-Manu-Track training, participants used their nonparetic and paretic hands under active and passive modes.
• In mode 1, both arms were guided passively by the device.
• In mode 2, the paretic arm guided the nonparetic arm and/or was guided by the nonparetic arm, depending on the participant’s level of motor severity.
• In mode 3, the nonparetic arm had to overcome continual resistance through the entire movement, and the paretic arm had to overcome only the initial resistance, which was set by the therapist according to the resistance against which the participant performed the voluntary movement with maximal force.
The robot was equipped with a computer game to provide instant visual movement feedback and to increase motivation.
The RBAT was followed by 15 min of functional task practice, which included various unilateral tasks and/or bilateral tasks such as picking up coins and opening a jar with one hand stabilizing while the other hand manipulated.
Clinical Measures
We used changes in motor impairment and motor function on the FMA-UL and the WMFT to evaluate motor recovery. The FMA is a widely used quantitative measure of sensorimotor conditions in stroke survivors. Each item was rated on a 3-point ordinal scale, with 2 points indicating performed completely, 1 point was given for partial performance, and 0 indicated cannot perform. The 33 UL items were used to measure the movement and reflexes of the shoulder, elbow, forearm, wrist, and hand, and coordination (Appendix 1). The FMA motor subscale shows high reliability, validity, and responsiveness for stroke survivors (Hsieh et al., 2009).
Post-stroke UL motor function was assessed with the WMFT (Wolf et al., 1989). Participants were scored for functional ability on a 6-point ordinal scale as they executed 15 activities that included gross and fine motor tasks. The functional ability scores (WMFT-FAS) of the paretic arm are reported in Appendix 2. We used a 6-point ordinal scale, where 0 indicated “does not attempt with the involved arm” and five indicated “arm does participate; movement appears to be normal.” The clinometric properties of the WMFT have been ascertained in stroke (Wolf et al., 2005).
The FIM consisted of 18 items grouped into six subscales measuring self-care, sphincter control, transfer, locomotion, communication, and social cognition ability (Appendix 3) (Hamilton et al., 1987). Each item was rated from 1–7 according to the required level of assistance to perform the tasks (e.g., 1, complete assistance; 7, complete independence). A higher score on any subscale represented less disability. The FIM has good interrater reliability, construct validity, and discriminant validity (Hamilton et al., 1994).
RS-fMRI Data Acquisition
MRIs were acquired using a 3.0-Tesla TIM Trio MRI scanner (Siemens, Erlangen, Germany). Tight but comfortable foam padding was used to minimize head movement, and earplugs were used to reduce scanner noise. RS-fMRI data were obtained using a single-shot gradient-echo sequence, echo planar imaging sequence (repetition time, 2000 ms; echo time, 30 ms; field of view, 220 mm; flip angle, 90°; matrix, 64 × 64; slice thickness, 4 mm; 36 slices/slab covering the entire brain, 180 volumes; and acquisition time, 6 min and 10 s). RS-fMRI data were acquired twice during a 4-week period. Participants were instructed to keep their eyes closed, to remain awake, to remain motionless, and not to think of anything in particular during the RS-fMRI scan.
Pre-processing for RS-fMRI Data
Before the data pre-processing, we flipped the imaging data from left to right along the midsagittal line for the three participants who had lesions on the left hemisphere. For all participants, the right side corresponded to the ipsilesional hemisphere. The RS-fMRI data were preprocessed using Statistical Parametric Mapping 8 software1. The beginning 10 volumes from each participant were discarded to allow the signal to reach T1 equilibrium and participants to adapt to the scanning noise. The remaining 170 volumes were first processed for slice timing and realigned to the middle volume to correct for interscan head motions. No participants had a maximum displacement of >2 mm or a maximum rotation of >2.0°. The remaining data set was spatially normalized to the Montreal Neurological Institute echo planar imaging template and resampled into 3 × 3 × 3 mm3 voxels. Thereafter, nuisance variables, including the averaged signals of the ventricular, white matter, and the whole brain, and Friston 24 regressors (Friston et al., 1996) were regressed out from the fMRI data.
Next, a band-pass frequency filter (0.01–0.08 Hz) was applied to reduce low-frequency drift and high-frequency noise. Finally, the filtered blood oxygen level—dependent images were spatially smoothed using an isotropic Gaussian kernel of 8 mm full width at half maximum.
RS-FC Analysis
RS-FC analyses were conducted using the Data Processing Assistant for Resting-State fMRI toolbox (Yan and Zang, 2010). We defined two regions of interest (ROIs; ipsilesional and contralesional M1) individually for each participant using a combination of anatomical and functional criteria (Fan et al., 2014). The anatomical M1 was defined to include voxels covering approximately the caudal half of the precentral gyrus along the anterior wall of the central sulcus, based on the Harvard–Oxford Atlas in the standard space. We overlaid each participant’s statistical parameter map for post-treatment vs. pre-treatment fMRI scans on his or her high-resolution anatomical scan and chose all active voxels within a radius of 10 mm around particular anatomical landmark.
The defined seed masks were used as the ROIs to perform the ROI-based RS-FC analysis. For each participant, the Pearson correlation coefficient between the mean time series of the ipsilesional and contralesional M1 was computed. A Fisher r-to-Z transform was used on the correlation coefficient value of each participant to improve the normality of the correlation coefficient.
In addition to the ROI approach, a whole-brain voxel-by-voxel analysis was performed for each participant using the average time course from the ipsilesional M1 as the signal of interest in a general linear model. A within-group voxel-based analysis across participants was then conducted using a general linear mixed model after transformation of the image data into the standardized space of the Montreal Neurological Institute. This analysis generated whole-brain connectivity maps for pre-treatment and post-treatment as well as a difference map between pre-treatment and post-treatment. The individual’s sex and age were entered as covariates of no interest. All maps were generated at a family-wise error-corrected p value of 0.05 and an extent threshold of κ > 10 voxels.
Statistical Analysis
One sample t-tests against zero were performed on the whole-brain connectivity maps to detect the brain regions showing significant functional connectivity with the ipsilesional M1 before and after treatment. Paired t-tests were then conducted to identify changes on whole-brain connectivity maps between pre-treatment and post-treatment. A paired Wilcoxon test was used to identify differences between pre-treatment and post-treatment scores on the FMA-UL, WMFT-FAS, FIM, and M1-M1 connectivity values. Differences (%) between pre-treatment and post-treatment scores [(post-treatment—pre-treatment)/(post-treatment + pre-treatment)] on the FMA-UL, WMFT-FAS, FIM, and M1-M1 functional connectivity were used to estimate the changes in motor and functional outcomes as well as interhemispheric RS-FC from pre-treatment to post-treatment. The Spearman correlation test was performed to examine relationships between the RS-FC and motor and functional recovery.
To assess the relationship between the changes of interhemispheric RS-FC and motor performance, motor function, and ADL recovery after RBAT, correlation analyses were performed to identify whether the pre-to-post difference in M1-M1 connectivity was significantly associated with changes in the FMA-UL score, WMFT-FAS score, and FIM total score.
Results
Clinical Measures
The results of the FMA-UL, WMFT-FAS, and FIM are presented in Table 1. All participants had substantial deficits in motor performance, functional use of the ULs, and daily function before treatment.
The results showed that there were significant differences between pre-treatment and post-treatment at the corrected level of significance (p < 0.017) on all clinical measures. The paired Wilcoxon test on the FMA-UL total scores revealed that participants showed significant improvements in levels of motor impairment from pre-treatment to the end of RBAT (Z = 2.82, p = 0.005). Moreover, the WMFT-FAS and FIM data indicated that eligible participants had better motor function (Z = 2.81, p = 0.005) and functional independence (Z = 2.80, p = 0.005) after RBAT relative to pre-treatment.
Functional Connectivity Results
The paired Wilcoxon test on the value of the M1-M1 RS-FC showed that participants had significantly increased M1-M1 functional connectivity from pre-treatment to the end of RBAT (Z = 2.80, p = 0.005).
A one sample t-test showed that for the ipsilesional M1 pre-treatment, participants had positive RS-FC with the bilateral middle frontal gyrus, bilateral cerebellum, bilateral inferior frontal gyrus, bilateral thalamus, ipsilesional angular gyrus, ipsilesional posterior cingulate cortex, ipsilesional superior frontal gyrus, contralesional M1, contralesional caudate nucleus, and contralesional precuneus. Moreover, negative RS-FC was observed before treatment between the ipsilesional M1 and the bilateral middle temporal gyrus, ipsilesional somatosensory cortex ipsilesional SMA, ipsilesional insula, ipsilesional superior parietal lobule, and contralesional M1 (Figure 1A and Table 2). Upon completion of RBAT, positive RS-FC with the ipsilesional M1 was seen in the bilateral somatosensory cortex (SI/SII), bilateral posterior cingulate cortex, bilateral cerebellum, bilateral thalamus, ipsilesional SMA, ipsilesional middle temporal gyrus, contralesional M1, contralesional inferior frontal gyrus, contralesional caudate nucleus, contralesional medial prefrontal cortex, contralesional anterior cingulate cortex (ACC), and contralesional middle frontal gyrus. However, participants had negative RS-FC between the ipsilesional M1 and the ipsilesional inferior frontal gyrus, ipsilesional middle frontal gyrus, ipsilesional superior frontal gyrus, contralesional temporal pole, contralesional inferior temporal gyrus, and contralesional insula after RBAT (Figure 1B and Table 2).
Figure 1. Positive resting-state functional connectivity results of one-sample t-test for (A) pre-treatment and (B) post-treatment.
Table 2. Regions showing significant positive and negative functional connectivity before and after treatment.
Figure 2 shows the maps exhibiting significant differences in RS-FC between pre-treatment and post-treatment. These brain regions are summarized in Table 3. When compared with post-treatment, greater RS-FC of the ipsilesional M1 with contralesional-lateralized brain regions was observed before treatment (Figure 2A). In contrast, increases in RS-FC were observed between the ipsilesional M1 seed and bilateral medial prefrontal cortex, bilateral M1, bilateral cerebellum, bilateral superior temporal gyrus, ipsilesional middle temporal gyrus, ipsilesional inferior parietal lobule (IPL), ipsilesional SMA, ipsilesional posterior cingulate cortex, ipsilesional SI/SII, ipsilesional caudate nucleus, contralesional ACC, contralesional insula, and contralesional middle occipital gyrus after RBAT relative to pre-treatment (Figure 2B).
Figure 2. Resting-state functional connectivity results of paired t-tests between (A) pre-treatment vs. post-treatment and (B) post-treatment vs. pre-treatment.
Table 3. Regions showing the significant differences in resting-state functional connectivity between pre-treatment and post-treatment.
Correlation of the RS-FC with Motor and Functional Recovery
Spearman correlation analysis showed that the pre-to-post difference in M1-M1 RS-FC was significantly positively correlated with changes in the WMFT-FAS score (R = 0.79, p = 0.006) and FIM total score (R = 0.92, p < 0.001). These indicated that participants with increased M1-M1 RS-FC after the intervention had greater gains in functional use of the affected arm and daily function. However, the relations between the pre-to-post difference M1-M1 connectivity and the changes of FMA-UL score were not significant (R = 0.55, p = 0.09).
Mediation Analysis Results
On the basis of a standard three-variable path model with a bootstrap test for the statistical significance of the product a × b, a single-level version of the mediation path model was used to get further insight of linkage between the clinical measures and RS-FC. Matlab coding implementing mediation analyses, developed by Wager et al. (2009) is freely available at2. In all participants, the change of interhemispheric M1-M1 functional connectivity from pre-treatment to post-treatment was a significant mediator in predicting the WMFT-FIM relation. The increased change in M1-M1 connectivity was associated with greater improvements in functional use of the affected arm and daily function after the intervention (a = 1.27, standard error = 0.61, p = 0.044; b = 0.17, standard error = 0.057, p = 0.021; a × b = 0.21, Z = 2.03, p = 0.042; Figure 3).
Figure 3. Mediation analysis results. (A) Path diagram shows the relationships between change scores on the WMFT-FAS and FIM from pre-treatment and the pre-to-post difference in interhemispheric M1-M1 RS-FC. The predictor region in the WMFT-FAS score is shown at the left, which predicts the M1-M1 RS-FC. This is the a path. The lines are labeled with path coefficients, and standard errors are shown in parentheses. The mediator factor (M1-M1 RS-FC) connection to the outcome (FIM total score) is the b path. This is calculated controlling for the WMFT-FAS and for the mediator factor, as is standard in mediation models. *p < 0.05, two-tailed. (B) Partial regression scatterplots for the (left panel) WMFT-FAS–M1-M1 RS-FC and for the (center panel) M1-M1 RS-FC–FIM total score relation. The right panel shows an example of a bootstrapped mediation effect (path a × b) for the M1-M1 RS-FC.
Discussion
This is the first study to demonstrate that RBAT may facilitate motor and functional recovery in participants with post-acute stroke by modulating the functional connectivity among the motor-related areas. All participants had greater gains in motor performance, motor function, and in daily function from pre-treatment to the end of RBAT. Increased M1-M1 RS-FC was shown after RBAT relative to pre-treatment. In addition, participants receiving RBAT had increased the intrahemispheric and interhemispheric RS-FC of the ipsilesional M1 for connections to the regions implicated in sensorimotor functions, spatial processing, and attention. Further, a greater improvement in the pre-to-post difference in M1-M1 connectivity was coupled with better motor and functional gains. Mediation analysis also confirmed that such interhemispheric M1-M1 interaction was a significant mediator for the relationship between the changes of WMFT-FAS scores and FIM total scores.
In agreement with previous research, motor performance and motor function were greatly improved in all participants after RBAT (Liao et al., 2012; Wu et al., 2013; Basteris et al., 2014). A significant improvement was also found for functional outcomes after RBAT compared with pre-treatment. Gains of more than 10 points were recorded on scores of the FMA in all participants, and these differences reached minimal clinically meaningful values (Arya et al., 2011). Compared with pre-treatment, each participant exhibited significantly greater improvement on UL functional use, as defined by the WMFT-FAS after treatment. These results indicated that greater motor improvements of RBAT may induce efficient UL functions. The result on ADL performance did not corroborate the findings of some previous studies (Lo et al., 2010), possibly because of the content of the robotic protocols (Basteris et al., 2014).
Our findings coincide with previous connectivity research on stroke recovery that M1-M1 RS-FC follows an evolution in the pattern of initial decrease and gradual restoration of the level of nearly normal (Carter et al., 2010; Wang et al., 2010; Park et al., 2011; Golestani et al., 2013). For ipsilesional M1 pre-treatment, participants had greater functional connectivity with the contralesional nonprimary sensorimotor regions compared with post-treatment. These results indicate that stroke lesions may induce changes in the contralesional nonprimary sensorimotor networks for supporting the imbalance of RS-FC between motor cortices during the early stage of stroke (Puh et al., 2007). After RBAT, participants had a greater interhemispheric M1-M1 connectivity compared with pre-treatment. This enhanced connectivity pattern may reflect the importance of proper equilibrium across hemispheres for optimal function (Bloom and Hynd, 2005; Manson et al., 2008).
Recent evidence also reveals that repetitive bimanual movements can increase the excitability of the ipsilesional M1 to inhibit the contralesional M1, which promotes a rebalancing of mutual transcallosal connections and enhances motor and functional recovery (Simkins et al., 2013; Stinear et al., 2014). Here, a greater improvement in the pre-to-post difference in M1-M1 RS-FC was coupled with better motor and functional gains after RBAT was completed. These results lend support to the notion that bilateral arm training can remodel functional connections between sensorimotor cortices and remediate imbalances between the two hemispheres for stroke recovery and that RS-FC can demonstrate this plasticity.
This study extended previous treatment-induced RS-FC studies on unilateral robotic therapy (Sergi et al., 2011; Saleh et al., 2012; Varkuti et al., 2013) and showed changes in functional connectivity between M1 and other motor-related areas after stroke and robotic training. First, the disturbed intracortical and intercortical functional connectivity of the ipsilesional M1 was noted, as in previous RS-fMRI studies (Wang et al., 2010; Park et al., 2011; Rehme et al., 2012; Golestani et al., 2013).
Second, we found increased functional connectivity between the SMA, middle temporal gyrus, superior temporal gyrus, medial prefrontal cortex, caudate nucleus, and cerebellum with the ipsilesional M1 from pre-treatment to the end of RBAT. The connectivity patterns within these networks may reflect the adaptive processes that occur during the reorganization of motor control and increased recruitment of regions disconnected by stroke (Bosnell et al., 2011). Moreover, increased pre-to-post differences in RS-FC in the ipsilesional M1-ACC and ipsilesional M1-IPL were also demonstrated after RBAT relative to pre-treatment. The IPL is important for spatial processing, motor intention, and visual-motor integration (Desmurget et al., 2009), and the ACC is engaged in attentional processes related to self-monitored movement in right hemispheric stroke(Hanlon et al., 2005). A strengthened association of these two regions has beneficial effects on motor and functional performances because it involves brain networks for processing of motor control, movement initiation, attention, and spatial perception (Varkuti et al., 2013).
We observed the increased pre-to-post differences in M1-M1 connectivity in response to the RBAT correlated positively with gains in the WMFT-FAS score and FIM total score. These findings support those from studies that showed similar relationships between the interhemispheric RS-FC with motor and functional outcomes from stroke (James et al., 2009; Carter et al., 2010; Hsieh et al., 2014; Urbin et al., 2014; Young et al., 2014). The results of our mediation analysis further indicated that a greater improvement in the pre-to-post difference in M1-M1 functional connectivity mediates the changes of functional use of the affected arm and daily function. Motor and functional recovery may be enhanced with reactivation and reorganization of the sensorimotor networks along with therapy (Calautti and Baron, 2003). Structural and functional imaging studies also suggest that disruption of the intraregional and interregional motor cortices are related to motor impairment and limited functional recovery (Granziera et al., 2012; Dijkhuizen et al., 2014). Importantly, the degree of functional recovery after stroke was associated with reinstatement of interhemispheric neuronal signal synchronization and normalization of cortical network organization (van Meer et al., 2012). The current study, combining correlation and mediation analyses, illuminates how treatment-induced changes in RS-FC in the primary motor cortices influence motor and functional outcomes after stroke. Specifically, RBAT may produce a restoration of the interregional neural connections for modulating later recovery of stroke.
This study has some limitations that warrant consideration. First, the sample size is rather small, and the eligible participants had heterogeneous patient-related (e.g., age, sex) and lesion-related (e.g., lesion size, stroke subtype, and lesion location) characteristics. Such interindividual variation might have differential effects on behavioral and neurological outcomes. Second, the study lacked a control group such as post-stroke patients receiving other interventions. Third, the study only examined the effectiveness at the end of treatment and did not explore the effects after some delay or longer-term treatment effects. Fourth, no task-evoked fMRI was performed as a means to determine the relationship between brain activity during motor tasks and RS-FC. Furthermore, studies involving acute or post-acute stroke patients should be implemented, where spontaneous recovery will need to be considered as a potential confounding factor. A specific study design, such as delayed treatment, should be applied to clarify possible treatment-induced effects. Future studies need to include a larger number of stroke survivors, have a small degree of interindividual variation, recruit control groups for comparison, and focus on the maintenance of therapy results over time.
Conclusion
We first examined RS-FC after bilateral robotic training, not only in M1-M1 but also in M1 with other motor-related areas. We also first used mediation analysis to capture the relations among intrinsic brain connections, motor performance, and functional independence after stroke rehabilitation. Increases in interhemispheric RS-FC between motor cortices were associated with better motor gains and improvements in ADLs performance after RBAT. These findings suggest that functional connectivity maps of the motor network may provide a prognostic value for motor and functional recovery after stroke and for monitoring the efficacy of rehabilitative therapies.
Author Contributions
FY-T, conducted the experiments, analyzed the data, interpreted the results and statistical analysis, and drafted the manuscript. WC-Y, developed the experimental setup, contributed to the study conception, interpreted the results, and edited the manuscript. LH-L, contributed to the fMRI experiment and assisted in interpreting the results. LK-C, contributed to the study conception and study design and interpreted the results. WY-Y, conducted the fMRI experiments and assisted in interpreting the results. CY-L, contributed to monitoring the participants and conducted the fMRI experiments. All authors have read and approved of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This study was supported in part by the National Health Research Institutes (NHRI-EX104–10403PI), the Ministry of Science and Technology (102–2314-B-002–154-MY2, 102–2628-B-182-005-MY3, and 103–2314-B-182–004-MY3), Healthy Ageing Research Center at Chang Gung University (EMRPD1E1711), Chang Gung Memorial Hospital (CMRPD 1C0403), and Excellent Research Projects of National Taiwan University in Taiwan. None of the authors have any conflicts of interest.
Supplementary Material
The Supplementary Material for this article can be found onlineat: http://journal.frontiersin.org/article/10.3389/fnhum.2015.00546/abstract
Footnotes
References
Arya, K. N., Verma, R., and Garg, R. K. (2011). Estimating the minimal clinically important difference of an upper extremity recovery measure in subacute stroke patients. Top Stroke Rehabil. 18(Suppl. 1), 599–610. doi: 10.1310/tsr18s01-599
Basteris, A., Nijenhuis, S. M., Stienen, A. H., Buurke, J. H., Prange, G. B., and Amirabdollahian, F. (2014). Training modalities in robot-mediated upper limb rehabilitation in stroke: a framework for classification based on a systematic review. J. Neuroeng. Rehabil. 11:111. doi: 10.1186/1743-0003-11-111
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
Bloom, J. S., and Hynd, G. W. (2005). The role of the corpus callosum in interhemispheric transfer of information: excitation or inhibition? Neuropsychol. Rev. 15, 59–71. doi: 10.1007/s11065-005-6252-y
Bosnell, R. A., Kincses, T., Stagg, C. J., Tomassini, V., Kischka, U., Jbabdi, S., et al. (2011). Motor practice promotes increased activity in brain regions structurally disconnected after subcortical stroke. Neurorehabil. Neural Repair. 25, 607–616. doi: 10.1177/1545968311405675
Calautti, C., and Baron, J. C. (2003). Functional neuroimaging studies of motor recovery after stroke in adults: a review. Stroke. 34, 1553–1566. doi: 10.1161/01.str.0000071761.36075.a6
Carter, A. R., Astafiev, S. V., Lang, C. E., Connor, L. T., Rengachary, J., Strube, M. J., et al. (2010). Resting interhemispheric functional magnetic resonance imaging connectivity predicts performance after stroke. Ann. Neurol. 67, 365–375. doi: 10.1002/ana.21905
Carter, A. R., Shulman, G. L., and Corbetta, M. (2012). Why use a connectivity-based approach to study stroke and recovery of function? Neuroimage. 62, 2271–2280. doi: 10.1016/j.neuroimage.2012.02.070
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
De Santis, D., Zenzeri, J., Casadio, M., Masia, L., Riva, A., Morasso, P., et al. (2014). Robot-assisted training of the kinesthetic sense: enhancing proprioception after stroke. Front. Hum. Neurosci. 8:1037. doi: 10.3389/fnhum.2014.01037
Desmurget, M., Reilly, K. T., Richard, N., Szathmari, A., Mottolese, C., and Sirigu, A. (2009). Movement intention after parietal cortex stimulation in humans. Science. 324, 811–813. doi: 10.1126/science.1169896
Dijkhuizen, R. M., Zaharchuk, G., and Otte, W. M. (2014). Assessment and modulation of resting-state neural networks after stroke. Curr. Opin. Neurol. 27, 637–643. doi: 10.1097/wco.0000000000000150
Dromerick, A. W., Lang, C. E., Birkenmeier, R., Hahn, M. G., Sahrmann, S. A., and Edwards, D. F. (2006). Relationships between upper-limb functional limitation and self-reported disability 3 months after stroke. J. Rehabil. Res. Dev. 43, 401–408. doi: 10.1682/jrrd.2005.04.0075
Fan, Y. T., Chen, C., Chen, S. C., Decety, J., and Cheng, Y. (2014). Empathic arousal and social understanding in individuals with autism: evidence from fMRI and ERP measurements. Soc. Cogn. Affect. Neurosci. 9, 1203–1213. doi: 10.1093/scan/nst101
Friston, K. J. (1994). Functional and effective connectivity in neuroimaging: A synthesis. Hum. Brain Mapp. 2, 56–78. doi: 10.1002/hbm.460020107
Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S., and Turner, R. (1996). Movement-related effects in fMRI time-series. Magn. Reson. Med. 35, 346–355. doi: 10.1002/mrm.1910350312
Fugl-Meyer, A. R., Jaasko, L., Leyman, I., Olsson, S., and Steglind, S. (1975). The post-stroke hemiplegic patient. 1. a method for evaluation of physical performance. Scand. J. Rehabil. Med. 7, 13–31.
Golestani, A. M., Tymchuk, S., Demchuk, A., and Goodyear, B. G. (2013). Longitudinal evaluation of resting-state FMRI after acute stroke with hemiparesis. Neurorehabil. Neural Repair. 27, 153–163. doi: 10.1177/1545968312457827
Granziera, C., Daducci, A., Meskaldji, D. E., Roche, A., Maeder, P., Michel, P., et al. (2012). A new early and automated MRI-based predictor of motor improvement after stroke. Neurology. 79, 39–46. doi: 10.1212/wnl.0b013e31825f25e7
Grefkes, C., and Ward, N. S. (2014). Cortical reorganization after stroke: how much and how functional? Neuroscientist. 20, 56–70. doi: 10.1177/1073858413491147
Gusnard, D. A., and Raichle, M. E. (2001). Searching for a baseline: functional imaging and the resting human brain. Nat. Rev. Neurosci. 2, 685–694. doi: 10.1038/35094500
Hamilton, B. B., Granger, C. V., Sherwin, F. S., Zielezny, M., and Tashman, J. S. (1987). “A uniform national data system for medical rehabilitation,” in Rehabilitation outcomes: Analysis and measurement, ed. M. J. Fuhrer (Baltimore, USA: Paul H. Brookes).
Hamilton, B. B., Laughlin, J. A., Fiedler, R. C., and Granger, C. V. (1994). Interrater reliability of the 7-level functional independence measure (FIM). Scand. J. Rehabil. Med. 26, 115–119.
Hanlon, C. A., Buffington, A. L., and McKeown, M. J. (2005). New brain networks are active after right MCA stroke when moving the ipsilesional arm. Neurology. 64, 114–120. doi: 10.1212/01.wnl.0000148726.45458.a9
Hesse, S., Werner, C., Pohl, M., Rueckriem, S., Mehrholz, J., and Lingnau, M. L. (2005). Computerized arm training improves the motor control of the severely affected arm after stroke: a single-blinded randomized trial in two centers. Stroke. 36, 1960–1966. doi: 10.1161/01.str.0000177865.37334.ce
Horwitz, B., Rumsey, J. M., and Donohue, B. C. (1998). Functional connectivity of the angular gyrus in normal reading and dyslexia. Proc. Natl. Acad. Sci. U. S. A. 95, 8939–8944. doi: 10.1073/pnas.95.15.8939
Housman, S. J., Scott, K. M., and Reinkensmeyer, D. J. (2009). A randomized controlled trial of gravity-supported, computer-enhanced arm exercise for individuals with severe hemiparesis. Neurorehabil. Neural Repair. 23, 505–514. doi: 10.1177/1545968308331148
Hsieh, Y. W., Lin, K. C., Horng, Y. S., Wu, C. Y., Wu, T. C., and Ku, F. L. (2014). Sequential combination of robot-assisted therapy and constraint-induced therapy in stroke rehabilitation: a randomized controlled trial. J. Neurol. 261, 1037–1045. doi: 10.1007/s00415-014-7345-4
Hsieh, Y. W., Wu, C. Y., Lin, K. C., Chang, Y. F., Chen, C. L., and Liu, J. S. (2009). Responsiveness and validity of three outcome measures of motor function after stroke rehabilitation. Stroke. 40, 1386–1391. doi: 10.1161/strokeaha.108.530584
Hsieh, Y. W., Wu, C. Y., Lin, K. C., Yao, G., Wu, K. Y., and & Chang, Y. J. (2012). Dose-response relationship of robot-assisted stroke motor rehabilitation the impact of initial motor status. Stroke. 43, 2729–2734. doi: 10.1161/strokeaha.112.658807
Hughes, C. M. L., Tommasino, P., Budhota, A., and Campolo, D. (2015). Upper extremity proprioception in healthy aging and stroke populations and the effects of therapist-and robot-based rehabilitation therapies on proprioceptive function. Front. Hum. Neurosci. 9:120. doi: 10.3389/fnhum.2015.00120
James, G. A., Lu, Z. L., VanMeter, J. W., Sathian, K., Hu, X. P., and Butler, A. J. (2009). Changes in resting state effective connectivity in the motor network following rehabilitation of upper extremity poststroke paresis. Top Stroke Rehabil. 16, 270–281. doi: 10.1310/tsr1604-270
Kwakkel, G., Kollen, B., and Lindeman, E. (2004). Understanding the pattern of functional recovery after stroke: facts and theories. Restorat. Neurol. Neurosci. 22, 281–299.
Liao, W. W., Wu, C. Y., Hsieh, Y. W., Lin, K. C., and Chang, W. Y. (2012). Effects of robot-assisted upper limb rehabilitation on daily function and real-world arm activity in patients with chronic stroke: a randomized controlled trial. Clin. Rehabil. 26, 111–120. doi: 10.1177/0269215511416383
Lin, K. C., Chen, Y. A., Chen, C. L., Wu, C. Y., and Chang, Y. F. (2010). The effects of bilateral arm training on motor control and functional performance in chronic stroke: a randomized controlled study. Neurorehabil. Neural. Repair. 24, 42–51. doi: 10.1177/1545968309345268
Lo, A. C., Guarino, P. D., Richards, L. G., Haselkorn, J. K., Wittenberg, G. F., Federman, D. G., et al. (2010). Robot-assisted therapy for long-term upper-limb impairment after stroke. N. Engl. J. Med. 362, 1772–1783. doi: 10.1056/NEJMoa0911341
Manson, S. C., Wegner, C., Filippi, M., Barkhof, F., Beckmann, C., Ciccarelli, O., et al. (2008). Impairment of movement-associated brain deactivation in multiple sclerosis: further evidence for a functional pathology of interhemispheric neuronal inhibition. Exp. Brain. Res. 187, 25–31. doi: 10.1007/s00221-008-1276-1
Murphy, T. H., and Corbett, D. (2009). Plasticity during stroke recovery: from synapse to behavior. Nat. Rev. Neurosci. 10, 861–872. doi: 10.1038/nrn2735
Park, C. H., Chang, W. H., Ohn, S. H., Kim, S. T., Bang, O. Y., Pascual-Leone, A., et al. (2011). Longitudinal changes of resting-state functional connectivity during motor recovery after stroke. Stroke. 42, 1357–1362. doi: 10.1161/strokeaha.110.596155
Puh, U., Vovk, A., Sevsek, F., and Suput, D. (2007). Increased cognitive load during simple and complex motor tasks in acute stage after stroke. Int. J. Psychophysiol. 63, 173–180. doi: 10.1016/j.ijpsycho.2006.03.011
Raichle, M. E., and Mintun, M. A. (2006). Brain work and brain imaging. Annu. Rev. Neurosci. 29, 449–476. doi: 10.1146/annurev.neuro.29.051605.112819
Rehme, A. K., and Grefkes, C. (2013). Cerebral network disorders after stroke: evidence from imaging-based connectivity analyses of active and resting brain states in humans. J. Physiol. 591, 17–31. doi: 10.1113/jphysiol.2012.243469
Rehme, A. K., Eickhoff, S. B., Rottschy, C., Fink, G. R., and Grefkes, C. (2012). Activation likelihood estimation meta-analysis of motor-related neural activity after stroke. Neuroimage. 59, 2771–2782. doi: 10.1016/j.neuroimage.2011.10.023
Saleh, S., Adamovich, S. V., and Tunik, E. (2012). Resting state functional connectivity and task-related effective connectivity changes after upper extremity rehabilitation: a pilot study. Conf. Proc. IEEE. Eng. Med. Biol. Soc. 2012. 4559–4562. doi: 10.1109/embc.2012.6346981
Sergi, F., Krebs, H. I., Groissier, B., Rykman, A., Guglielmelli, E., Volpe, B. T., et al. (2011). Predicting efficacy of robot-aided rehabilitation in chronic stroke patients using an MRI-compatible robotic device. Conf. Proc. IEEE. Eng. Med. Biol. Soc. 2011. 7470–7473. doi: 10.1109/iembs.2011.6091843
Simkins, M., Kim, H., Abrams, G., Byl, N., and Rosen, J. (2013). Robotic unilateral and bilateral upper-limb movement training for stroke survivors afflicted by chronic hemiparesis. IEEE. Int. Conf. Rehabil. Robot. 3:6650506. doi: 10.1109/icorr.2013.6650506
Stinear, C. M., Petoe, M. A., Anwar, S., Barber, P. A., and Byblow, W. D. (2014). Bilateral priming accelerates recovery of upper limb function after stroke: a randomized controlled trial. Stroke. 45, 205–210. doi: 10.1161/strokeaha.113.003537
Teng, E. L., and Chui, H. C. (1987). The Modified Mini-Mental State (3MS) examination. J. Clin. Psychiatry. 48, 314–318. doi: 10.1007/springerreference_183334
Urbin, M. A., Hong, X., Lang, C. E., and Carter, A. R. (2014). Resting-state functional connectivity and its association with multiple domains of upper-extremity function in chronic stroke. Neurorehabil. Neural. Repair. 28, 761–769. doi: 10.1177/1545968314522349
Van Delden, A. E. Q., Peper, C. E., Beek, P. J., and Kwakkel, G. (2012). Unilateral versus bilateral upper limb exercise therapy after stroke: a systematic review. J. Rehabil. Medicine. 44, 106–117. doi: 10.2340/16501977-0928
Van Essen, D. C., Ugurbil, K., Auerbach, E., Barch, D., Behrens, T. E., Bucholz, R., et al. (2012). The Human Connectome Project: a data acquisition perspective. Neuroimage. 62, 2222–2231. doi: 10.1016/j.neuroimage.2012.02.018
van Meer, M. P., Otte, W. M., van der Marel, K., Nijboer, C. H., Kavelaars, A., van der Sprenkel, J. W., et al. (2012). Extent of bilateral neuronal network reorganization and functional recovery in relation to stroke severity. J. Neurosci. 32 4495–4507. doi: 10.1523/jneurosci.3662-11.2012
Varkuti, B., Guan, C., Pan, Y., Phua, K. S., Ang, K. K., Kuah, C. W., 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
Wager, T. D., Waugh, C. E., Lindquist, M., Noll, D. C., Fredrickson, B. L., and Taylor, S. F. (2009). Brain mediators of cardiovascular responses to social threat: part I: Reciprocal dorsal and ventral sub-regions of the medial prefrontal cortex and heart-rate reactivity. Neuroimage. 47, 821–835. doi: 10.1016/j.neuroimage.2009.05.043
Waller, S., Whitall, J., Jenkins, T., Magder, L. S., Hanley, D. F., Goldberg, A., et al. (2014). Sequencing bilateral and unilateral task-oriented training versus task oriented training alone to improve arm function in individuals with chronic stroke. BMC. Neurol. 14:236. doi: 10.1186/s12883-014-0236-6
Wang, L., Yu, C., Chen, H., Qin, W., He, Y., Fan, F., et al. (2010). Dynamic functional reorganization of the motor execution network after stroke. Brain. 133, 1224–1238. doi: 10.1093/brain/awq043
Wei, X. J., Tong, K. Y., and Hu, X. L. (2011). The responsiveness and correlation between Fugl-Meyer Assessment, Motor Status Scale and the Action Research Arm Test in chronic stroke with upper-extremity rehabilitation robotic training. Int. J. Rehabil. Res. 34, 349–356. doi: 10.1097/mrr.0b013e32834d330a
Wolf, S. L., Lecraw, D. E., Barton, L. A., and Jann, B. B. (1989). Forced use of hemiplegic upper extremities to reverse the effect of learned nonuse among chronic stroke and head-injured patients. Exp. Neurol. 104, 125–132. doi: 10.1016/s0014-4886(89)80005-6
Wolf, S. L., Thompson, P. A., Morris, D. M., Rose, D. K., Winstein, C. J., Taub, E., et al. (2005). The EXCITE trial: attributes of the Wolf Motor Function Test in patients with subacute stroke. Neurorehabil. Neural. Repair. 19, 194–205. doi: 10.1177/1545968305276663
Wu, C. Y., Yang, C. L., Chen, M. D., Lin, K. C., and Wu, L. L. (2013). Unilateral versus bilateral robot-assisted rehabilitation on arm-trunk control and functions post stroke: a randomized controlled trial. J. Neuroeng. Rehabil. 10:35. doi: 10.1186/1743-0003-10-35
Xu, H., Qin, W., Chen, H., Jiang, L., Li, K., and Yu, C. (2014). Contribution of the resting-state functional connectivity of the contralesional primary sensorimotor cortex to motor recovery after subcortical stroke. PLoS One 9:e84729. doi: 10.1371/journal.pone.0084729
Yan, C.-G., and Zang, Y.-F. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4, 13. doi: 10.3389/fnsys.2010.00013
Keywords: functional connectivity, stroke, rehabilitation, robot, fMRI, motor cortex
Citation: Fan Y-t, Wu C-y, Liu H-l, Lin K-c, Wai Y-y and Chen Y-l (2015) Neuroplastic changes in resting-state functional connectivity after stroke rehabilitation. Front. Hum. Neurosci. 9:546. doi: 10.3389/fnhum.2015.00546
Received: 23 April 2015; Accepted: 17 September 2015;
Published: 06 October 2015.
Edited by:
Jean-Claude Baron, University of Cambridge, UKReviewed by:
Martin Lotze, University of Greifswald, GermanyErin Burke Quinlan, University of California, Irvine, USA
Copyright © 2015 Fan, Wu, Liu, Lin, Wai and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and 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: Keh-chung Lin, School of Occupational Therapy, College of Medicine, National Taiwan University and Division of Occupational Therapy, Department of Physical Medicine and Rehabilitation, National Taiwan University Hospital, 17 F4, Xu Zhou Road, Taipei, Taiwan, a2VoY2h1bmdsaW5AbnR1LmVkdS50dw==
† These authors have contributed equally to this work.