Neurophysiological signatures of hand motor response to dual-transcranial direct current stimulation in subacute stroke: a TMS and MEG study
Journal of NeuroEngineering and Rehabilitation volume 17, Article number: 72 (2020)
Dual transcranial direct current stimulation (tDCS) to the bilateral primary motor cortices (M1s) has potential benefits in chronic stroke, but its effects in subacute stroke, when behavioural effects might be expected to be greater, have been relatively unexplored. Here, we examined the neurophysiological effects and the factors influencing responsiveness of dual-tDCS in subacute stroke survivors.
We conducted a randomized sham-controlled crossover study in 18 survivors with first-ever, unilateral subcortical ischaemic stroke 2–4 weeks after stroke onset and 14 matched healthy controls. Participants had real dual-tDCS (with an ipsilesional [right for controls] M1 anode and a contralesional M1 [left for controls] cathode; 2 mA for 20mins) and sham dual-tDCS on separate days, with concurrent paretic [left for controls] hand exercise. Using transcranial magnetic stimulation (TMS) and magnetoencephalography (MEG), we recorded motor evoked potentials (MEPs), the ipsilateral silent period (iSP), short-interval intracortical inhibition, and finger movement-related cortical oscillations before and immediately after tDCS.
Stroke survivors had decreased excitability in ipsilesional M1 with a relatively excessive transcallosal inhibition from the contralesional to ipsilesional hemisphere at baseline compared with controls, as quantified by decreased MEPs and increased iSP duration. Dual-tDCS led to increased MEPs and decreased iSP duration in ipsilesional M1. The magnitude of the tDCS-induced MEP increase in stroke survivors was predicted by baseline contralesional-to-ipsilesional transcallosal inhibition (iSP) ratio. Baseline post-movement synchronization in α-band activity in ipsilesional M1 was decreased after stroke compared with controls, and its tDCS-induced increase correlated with upper limb score in stroke survivors. No significant adverse effects were observed during or after dual-tDCS.
Task-concurrent dual-tDCS in subacute stroke can safely and effectively modulate bilateral M1 excitability and inter-hemispheric imbalance and also movement-related α-activity.
Transcranial direct current stimulation (tDCS) has been demonstrated to non-invasively modulate cortical excitability in the primary motor cortex (M1) in both controls and stroke survivors [1, 2]. Motor evoked potentials (MEPs) typically increase following unilateral anodal tDCS and decrease following unilateral cathodal tDCS, with effects outlasting stimulation by minutes to hours . This polarity-specific modulation has been suggested as a putative way to promote post-stroke motor recovery, either by enhancing ipsilesional M1 excitability with anodal tDCS or decreasing contralesional M1 excitability with cathodal tDCS . Some tDCS studies have suggested that these approaches may be promising, but the quality of evidence for tDCS in promoting post-stroke motor recovery is still low to moderate , partially due to heterogeneity in study designs and stroke survivor profiles across the literature, and small sample sizes within studies . To successfully translate tDCS effects into clinical benefits, patient selection based on neurophysiological status for a given effective tDCS montage may be required [6,7,8,9].
Dual, or bi-hemispheric, tDCS involves concurrent anodal stimulation to one M1 and cathodal stimulation to the other M1. The effects of dual tDCS have been studied in both healthy controls and stroke survivors in a number of different ways. Dual-tDCS has been shown to have a significant beneficial effect in dexterity  and motor learning [11,12,13,14] in controls. The effects of dual-tDCS on motor performance were better or at least equal to the unilateral anodal stimulation effect in most of these studies [10, 11, 13, 14]. Consistent with these behavioural findings, increased MEPs in the anode-targeted M1 and decreased MEPs in the cathode-targeted M1 after dual-tDCS have been demonstrated in most studies in healthy controls [10, 14,15,16,17], suggesting that dual-tDCS may result in additive effects of unilateral stimulation, although these effects are not entirely consistent .
The effects of dual-tDCS post stroke, when an ipsilesional anode and a contralesional cathode have been paired with concurrent rehabilitation, are much more varied across studies, with some studies demonstrating significantly enhanced motor score, dexterity or grip strength [19,20,21,22,23] in chronic stroke survivors, but other studies showing no effect [24,25,26]. Only one previous study has explored the immediate post-stimulation effects of dual-tDCS on cortical excitability, which showed no significant changes in MEPs or transcallosal inhibition in six subacute stroke survivors . Two studies have combined repetitive dual-tDCS with rehabilitation in chronic stroke survivors and demonstrated an increase in ipsilesional MEPs compared to sham [19, 20]. Sensorimotor dynamic activity, recorded using electroencephalography (EEG) or magnetoencephalography (MEG), are also altered following stroke [28,29,30,31], and changes in these measures are correlated with motor function [30, 31]. The event-related desynchronization (ERD) before and during movement has been associated with motor preparation and execution, while event-related synchronization (ERS) after cessation of movement has been suggested to reflect motor deactivation . Dual-tDCS has been shown to alter motor imagery-related hemispheric lateralization in sensorimotor rhythms in controls  and after stroke . However, the effect of dual-tDCS on movement-related dynamic activity has not yet been examined.
We hypothesized that task-concurrent dual-tDCS could enhance ipsilesional corticospinal excitability, modulate movement-related brain oscillations, and also rebalance the asymmetry of interhemispheric interactions after subacute stroke. Potential for recovery is most prominent early after stroke and declines gradually after 6 months (chronic stage) [35, 36]. Most rehabilitation intervention studies, including tDCS studies, recruit patients at the chronic stage of recovery , but it is likely that rehabilitative effects would be stronger if applied earlier. In addition, recent studies have also demonstrated that different adaptive mechanisms and bi-hemispheric interactions may exist during the subacute and chronic stages [38, 39]. Studying effects in subacute stroke is therefore important if we are to fully explore the potential of dual-tDCS for stroke recovery. Cortical and subcortical infarctions have highly different outcomes, network reorganization and tDCS responsiveness . To minimize variations in tDCS substrates (the cortex) and responsiveness in this study, we enrolled stroke survivors with relatively homogenous subcortical infarctions at 2–4 weeks post-stroke to investigate the effects of dual-tDCS using transcranial magnetic stimulation (TMS) and MEG in a randomized, sham-controlled, crossover study design to compare the effects of dual-tDCS in subacute stroke survivors and age- and gender-matched controls, as well as explore factors predicting responsiveness to dual-tDCS.
Materials and methods
Subjects and study design
We recruited patients aged 20–80 years with first-ever, unilateral, ischaemic stroke with mild to moderate hand paresis (Medical Research Council motor score 3–4) 2–4 weeks after stroke. Exclusion criteria were: major neurological or medical comorbidities or absent MEPs in the paretic hand. Age- and gender-matched healthy controls were also recruited. All participants gave their informed consent to participate in the study, which was approved by the institutional review board of Taipei Veterans General Hospital.
Twenty-one stroke survivors fulfilled the inclusion/exclusion criteria and 3 of them dropped out (14%) after recruitment because of unwillingness to complete the full assessment. The remaining 18 stroke survivors were enrolled a median of 23 days (range:14–28) after stroke. Survivors had a median age of 63 (range:31–76 years) and equal gender distribution (Table 1). They had mild to moderate hand paresis as defined by the Medical Research Council (MRC) score 3–4 of the extensor carpi radialis (ECR) muscle for being suitable for TMS measures from the ECR and MEG measures of finger lifting task. Their median scores of Fugl-Meyer assessment-upper extremity scale (FMA-UE) and Action Research Arm Test (ARAT) were 60 (range:24–64) and 56 (range:3–57), respectively. All stroke survivors were right-handed , and 10 of them had left hemiparesis. Lesions were mainly located in the deep nuclei and white matter of the middle cerebral artery territory. Four subjects had a pons infarction (Fig. 1a). 14 age (median:59, range:37–72) and gender (6 male) matched right-handed controls were also enrolled [Age: U = 78, p = 0.07; Gender χ2(1) = 0.69, p = 0.74].
All subjects participated in four experimental sessions on separate days: two real and two sham dual-tDCS sessions, which were combined with either TMS or MEG recordings (i.e. TMS + real tDCS, TMS + sham tDCS, MEG + real tDCS, MEG + sham tDCS). The order of the sessions was counterbalanced across the groups and there was at least 72 h between sessions (Fig. 1b).
Transcranial direct current stimulation
A tDCS stimulator (Eldith, UK) with 5x5cm conductive rubber electrodes and saline-soaked sponges was used. M1 representations were localized in each participant from ECR hotspot derived from single-pulse TMS. The anode was placed over the ipsilesional [right for controls] M1 [anode-targeted M1; M1Anode], and the cathode over the contralesional [left for controls] M1 [cathode-targeted M1; M1Cathode] (Fig. 1c). Impedance was kept below 5 kΩ. For real stimulation, a 2 mA current (current density 0.08 mA/cm2) was applied for 20 min, including 30 s ramp-up and ramp-down. This tDCS dosage was safe according to the safety guideline and used by previous dual-tDCS studies [6, 41]. For sham stimulation, the stimulator was turned off immediately after the initial ramp-up period.
During the tDCS, subjects performed full range self-paced extension of the paretic [left for controls] ECR, with wrist extension over 3 s followed by 3 s of wrist flexion back to neutral position. This was repeated for 4 min followed by 1 min of rest for 4 cycles (Fig. 1b and d). The Adverse Effects Questionnaire  was administered after each tDCS session.
Transcranial magnetic stimulation
In the TMS sessions, MEPs, the ipsilateral silent period (iSP), and short-interval intracortical inhibition (SICI) were recorded from both M1s before (i.e. baseline) and immediately after real or sham tDCS (0, 15, 30 min) Both resting motor threshold (rMT) and active motor threshold (aMT) were measured once before either real or sham tDCS stimulation as the baseline values. The fixed evaluation order of each TMS session was rMT → MEP → SICI→iSP → aMT of the non-paretic hand [right for controls)] and then the paretic hand [left for controls]. In line with the literatures, we observed the after-effects of a single-session tDCS lasted for about 30–40 min . Therefore, the TMS measurements were completed by 40 min after tDCS to be within an estimated effective window. We utilized a Magstim Rapid2 stimulator (Magstim, UK) with a 90 mm diameter double cone coil and a synchronized surface electromyography recording system (MedelecSynergy; VIASYS HealthCare, UK) for rMT, aMT, MEP amplitude, and iSP measurements. Two Magstim 2002 stimulators were used for SICI. The rMT was defined as the minimum intensity required to elicit MEPs with a peak-to-peak amplitude greater than 50 μV in 5 out of 10 trials in the relaxed ECR . The aMT was defined as the minimum intensity required to evoke MEPs with a peak-to-peak amplitude greater than 200 μV in 5 out of 10 trials, while the subject maintained at ~ 20% of maximum contraction. The strength of muscle contraction was visually monitored and maintained by concurrent surface EMG recording (Fig. 1d). MEP amplitude was quantified as the mean peak-to-peak amplitude from 12 successive TMS pulses evoked every 5 s at an intensity of 120% baseline rMT. The iSP was calculated from the mean of 5 trials and was defined as the length of disrupted ongoing muscular activity of the ECR after an ipsilateral M1 TMS pulse at an intensity of 150% rMT. The participant was asked to maintain steady wrist extension with maximal strength monitored by surface EMG during the test period. The onset point of iSP was defined as disrupted ongoing muscular activity below the mean EMG amplitude of the baseline before the TMS pulse. The offset point of iSP was defined as EMG amplitude recovery to 50% of the baseline before the TMS pulse [44, 45]. The iSP ratio of contralesional / ipsilesional [left/right for controls] iSP was used to examine the interhemispheric transcallosal inhibition imbalance between contralesional to ipsilesional [left to right for controls] hemisphere. SICI was quantified using a subthreshold conditioning stimulus (80% of the baseline rMT) and suprathreshold test stimulus (120% of the baseline rMT), with individualized inter-stimulus intervals (ISI) between 1 and 5 ms to prevent ceiling or floor effects. ISI remained constant within subject across the experiment. SICI was calculated as the percentage of conditioned MEP amplitude relative to unconditioned MEPs for each trial and was obtained from the mean of 5 trials with a 5 s interval [46, 47].
Magnetoencephalography acquisition and analysis
In the MEG sessions, movement-related cortical oscillations were recorded before (i.e. baseline), and immediately after real and sham tDCS. The MEG measurements were completed by 40 min after tDCS to be within an estimated effective window. A Neuromag Vectorview MEG (Elekta, Helsinki, Finland) was used to acquire electrophysiological data while performing a simple motor task (see  for more details). Briefly, subjects performed a self-paced, unilateral index finger lifting task every 7 s with the paretic [left for controls] hand, keeping other muscles relaxed and avoiding excessive blinking. The acquisition lasted for up to 15 min to collect around 50 adequate trials as determined by online MEG and electrooculography (EOG) monitoring with a pre-defined rejection threshold (4000 fT for MEG and 200 mV for EOG).
MEG data were acquired with a 500 Hz sampling rate and analysed using Brainstorm . MEG epochs comprised the − 3 s to 3 s relative to movement onset, identified by an optic detection pad. Trials were discarded if the interval between two consecutive finger lifts was shorter than 6 s, or the corresponding EOG or MEG data were noisy, as defined by visual inspection and aforementioned EOG or MEG pre-defined rejection threshold. MEG channels were co-registered with the individual-MRI. MEG data were then filtered into α-band (8–12 Hz) and β-band (16–30 Hz), and projected to source level by minimum norm estimate . The event-related desynchronization/synchronization (ERD/S) were then calculated as follows: ERD/ERS % = A-R/R × 100 (A: the power within the frequency band of interest during the active period of the event; R: the mean power of the reference period) . The peak amplitude of ERD [or ERS] was determined as the minimum [or maximum] between − 2 to 2 s [or 0 to 3 s] relative to movement onset for each individual . This procedure resulted in four MEG measures of interest, i.e. α-ERD, α-ERS, β-ERD, β-ERS peak amplitude.
SPSS version 24.0.0 (Chicago, USA) was used for statistical analyses. Data were transformed where necessary to achieve a normal distribution. Outliers were identified if studentized residuals were greater than ±3 and removed from analysis.
Repeated measures analysis of variance (ANOVAs) were used to compare the baseline TMS (rMT, aMT, MEPs, iSP and SICI) and MEG measures (α-ERD, α-ERS, β-ERD, β-ERS) within groups (baseline session 1, baseline session 2), while Mann-Whitney U tests were used to compare these baseline TMS and MEG measures between groups (stroke survivors, healthy controls).
To investigate the effect of tDCS, TMS measures were normalized to baseline (post / baseline) and normalized values were used for all analyses of change. MEG measures were normalized as follows: (post - baseline) / (post + baseline). ANOVA with a post hoc Bonferroni correction for multiple comparisons was used to compare the normalized TMS and MEG measures. For TMS, the ISI was used as the covariate in SICI analysis.
To define tDCS responsiveness, post-tDCS MEPs changes were categorized into three groups : High: if the mean at two or more time points (0, 15 mins, 30 mins) was greater than 15% for anodal tDCS or smaller than 15% for cathodal tDCS; Medium: between 10 and 15%; Low: < 10%. Stepwise multiple regression analysis with backward elimination was used to identify the independent factors for the tDCS response.
Dual-tDCS is well tolerated
The only self-reported discomfort was tingling during the initial 1–2 min of tDCS stimulation (four subjects, 14%). No seizures or other severe adverse events were reported.
Stroke survivors had decreased baseline excitability in ipsilesional M1 compared with controls
We first wanted to investigate whether there were any differences between stroke survivors and controls in baseline TMS measures, i.e. before tDCS was applied. There were no significant intra-individual variabilities in the baseline TMS measures before real versus sham stimulation in healthy controls and stroke survivors, respectively (all p > 0.1, supplementary Table A.1). Therefore, we used the mean of the baseline measures for subsequent analyses. Compared to controls, stroke survivors had increased rMT (51 ± 21% versus 38 ± 5% maximum stimulator output, U = 63, p = 0.017, r = 0.42), increased aMT (40% ± 12% versus 30 ± 6% maximum stimulator output, U = 59, p = 0.009, r = 0.45), and decreased ipsilesional M1 excitability (MEP amplitude at 120% rMT 282 ± 233 μV versus 753 ± 374 μV, U = 27, p < 0.001, r = 0.66). Stroke survivors also showed increased contralesional-to-ipsilesional transcallosal inhibition (iSP: 81 ± 13 ms versus 58 ± 4 ms, U = 1, p < 0.001, r = 0.18), but comparable ipsilesional SICI (30 ± 12% versus 27 ± 9%, U = 122, p = 0.88, r = 21.6; Fig. 2).
Dual-tDCS leads to significant changes in TMS metrics in stroke survivors and in controls
We then wanted to explore the effects of dual-tDCS on our TMS metrics. We first wished to investigate whether stroke survivors and controls responded differently to dual-tDCS. We therefore ran two mixed-design 2 × 2 × 3 ANOVAs, with group (patients, controls) as the between-subjects factor, stimulation (real, sham) and time (0, 15, 30 min) as within-subjects factors and normalized TMS measures (MEP, iSP, SICI) as the dependent variable for each hemisphere separately. There was no significant group difference between stroke survivors and controls in response to tDCS in either M1Anode or M1Cathode in any of the TMS metrics (all p’s > 0.1).
We were interested in the effects of dual-tDCS on TMS measures in stroke survivors to understand group-specific features (Fig. 3). In each group, a set of 2 × 3 repeated measures ANOVAs with factors of stimulation (real, sham) and time (0, 15 min, 30 min) were performed for post-tDCS changes of each normalized TMS measure (MEP, iSP, SICI) from M1Anode and M1Cathode separately.
In stroke survivors, there were significant main effects of stimulation [M1Anode: F (1,16) = 12.7, p = 0.003, ηp2 = 0.44; M1Cathode F (1,17) = 11.23, p = 0.004, ηp2 = 0.40], with an increase MEP amplitude in M1Anode and a decrease in MEP amplitude in M1Cathode after real stimulation relative to sham. There were no significant main effects of time or interactions between stimulation and time (see Additional file 1, Table A.2 with full statistics).
In terms of iSP there were significant main effects of stimulation [M1Anode: F (1,15) = 13.7, p = 0.002, ηp2 = 0.48; M1Cathode: F (1,15) = 12.5, p = 0.003, ηp2 = 0.46], reflecting a shorter iSP in M1Anode and a longer iSP in M1Cathode after real stimulation relative to sham. There were no significant main effects of time, and no significant interactions (Table A.2).
There were no main effects or interactions in either M1 in terms of SICI change (all p’s > 0.1, Table A.2). In line with the non-significant effects of group between stroke survivors and controls, the results in controls were highly similar to those in stroke survivors (Fig. 3, Table A.3).
Stroke survivors exhibit smaller baseline ipsilesional α- and β- ERSs than controls
We then wished to see whether the changes in neural excitability as revealed by TMS was reflected in changes in frequency-specific neural activity as measured using MEG. MEG data were acquired from 12 healthy controls and 11 stroke survivors. Initially, we wanted to investigate whether there were any differences in MEG metrics (α-ERD, α-ERS, β-ERD, β-ERS) before tDCS, between stroke survivors and controls (Fig. 4a). Since there were no differences in any baseline MEG metric between the two baseline measures in either group (all p’s > 0.1), the mean of the two baseline measures were used for subsequent analyses. The mixed model ANOVAs and statistics details were summarized in Table A.4. The baseline ERS peak amplitudes in the ipsilesional [right for controls] hemisphere was significantly smaller in stroke survivors than controls in both the α- [U = 27, p = 0.015, r = 0.51] and β-frequency bands [U = 17, p = 0.002, r = 0.64]. Baseline ERS peak amplitude in the contralesional [left for controls] hemisphere was significantly smaller in stroke survivors than controls in the β-band [U = 20, p = 0.005, r = 0.59], but not in the α-band [U = 32, p = 0.036, r = 0.44, α = 0.025].
Dual-tDCS leads to an increase in α-ERS in M1Anode in stroke survivors
Next, we wished to investigate the neural response to dual-tDCS. We first wished to examine whether there was a significant difference between stroke survivors and controls in response to stimulation. We therefore performed a set of 2 × 2x2 mixed-design ANOVAs with group (patients, controls) as a between-subjects factor, stimulation (real, sham) and frequency band (α, β) as within-subjects factors, and normalized MEG measure as dependent variable for M1Anode-ERD, M1Anode-ERS, M1Cathode-ERD and M1Cathode-ERS separately. These demonstrated that there were no significant differences between controls and stroke survivors in response to dual-tDCS in either hemisphere (all p’s > 0.1).
We then wanted to explore the group-specific effects of dual-tDCS in stroke survivors alone. Two 2 × 2 × 2 repeated measures ANOVAs with factors of stimulation (real, sham), frequency band (α, β) and time window (ERD, ERS) were performed for M1Anode and M1Cathode separately. In M1Anode, there was a significant main effect of stimulation [F (1,10) = 6.04, p = 0.034, ηp2 = 0.38], but not of frequency band or time window. There was a significant interaction between stimulation and time window [F (1,11) = 5.44, p = 0.042, ηp2 = 0.35]. Therefore, two 2x2 repeated measures ANOVAs with factors of stimulation (real, sham) and frequency band (α, β) were performed for M1Anode-ERD and M1Anode-ERS separately. No significant differences were observed in M1Anode-ERD (all p’s > 0.1). However, there was a significant main effect of stimulation in M1Anode-ERS, [F (1,10) = 10.88, p = 0.008, ηp2 = 0.52], but no significant main effect of frequency band nor stimulation by frequency band interaction. Post-hoc testing revealed that α-ERS in M1Anode significantly increased after stimulation compared with sham (t (10) = 3.19, p = 0.01, d = 0.96). There was also a trend of increase in β-ERS in M1Anode after stimulation compared with sham [t (10)=2.42, p = 0.036, d = 0.73, α = 0.025; Fig. 4b and c, full statistics in Table A.5).
The equivalent analyses in M1Cathode revealed no significant effects of stimulation (Table A.5). There were no significant effects of stimulation in controls in any MEG metric studied (all p’s > 0.1).
Dual-tDCS-induced change in α-ERS in ipsilesional M1 correlates with behaviour
This study was performed to understand the neurophysiological underpinnings of tDCS post-stroke. We therefore wished to explore the relationship between the observed stimulation-related increased α-ERS in the ipsilesional M1Anode, to address the hypothesis that tDCS-induced changes in neurophysiological measures may relate to functional status. We demonstrated that the stimulation-related change in α-ERS (i.e. change with real stimulation compared to sham) was positively correlated with concurrent paretic upper limb motor scores (r = 0.904, p < 0.001 for ARAT; r = 0.744, p = 0.009 for FMA-UE, Fig. 4d), such that stroke survivors with better function had larger stimulation-related α-ERS change in ipsilesional M1. No significant correlations were demonstrated between tDCS-induced change in any TMS measure and tDCS-induced change in any other MEG metrics (all p’s > 0.1).
Responsiveness to dual-tDCS was predicted by interhemispheric inhibition imbalance, as reflected by iSP ratio
There was no statistically difference in response rate between stroke survivors and controls. In M1Anode 93% of controls were “high” responders and 7% “medium”, 61% of stroke survivors were “high” responders and 28% “medium” [χ2(2) = 4.40, p = 0.07]. In M1Cathode the rates were 86 and 7% in controls and 67 and 5% in stroke survivors [χ2(2) = 2.89, p = 0.24]. 79% of controls and 44% of stroke survivors were high responders in both hemispheres [χ2(1) = 3.80, p = 0.075]. Only 3 stroke survivors were “low” responders in both hemispheres, all of whom had relatively severe hand paresis (cases 3, 8, 15).
A summary schematic of the neural changes elicited by dual-tDCS in both stroke survivors and controls is given in Fig. 5a. Given the variability in response to tDCS, we believed it would be important to explore the physiological underpinnings of this variability if possible. We therefore ran a multiple regression in stroke survivors to predict MEP change (change of MEP amplitude 15 min after real tDCS stimulation relative to sham) from six independent variables: age, gender, baseline FMA-UE score, baseline MEP ratio (baseline MEP amplitude of ipsilesional hemisphere (IH) relative to contralesional hemisphere (CH)), baseline iSP ratio (CH → IH baseline iSP / IH → CH baseline iSP), and baseline α-ERS ratio (baseline M1 ERS peak amplitude of IH relative to CH). After eliminating the least significant variables (p = 0.59 for age, p = 0.99 for gender, p = 0.57 for FMA-UE, p = 0.52 for baseline MEPs ratio, and p = 0.41 for baseline α-ERS ratio), the greatest explanatory power for the responsiveness to tDCS was achieved with the iSP ratio (R2 = -0.253, p = 0.047). In this study, stroke survivors showed significantly higher iSP ratio than controls [Stroke survivors 1.22 (0.99–1.95) [Median (Range)]; Controls 1.02 (0.99–1.95); U = 29.0, p = 0.001, r = 0.63)]. Greater baseline iSP ratios predicted a smaller response to tDCS (Fig. 5b). Similar analysis applied in healthy controls using the same independent variables (except baseline FMA-UE) revealed no significant predictors.
To our knowledge, this is the first study to demonstrate the changes of motor cortical excitability, transcallosal inhibition and neural oscillations in subacute stroke survivors after M1-M1 dual-tDCS stimulation. We demonstrated that dual-tDCS led to an increase in cortical excitability and a decrease in transcallosal inhibition, as reflected by a shorter iSP, in the anode-targeted, ipsilesional M1. We also demonstrated a decrease in excitability and an increase in transcallosal inhibition, in the cathode-targeted, contralesional M1. The interhemispheric balance in iSP ratio predicted an individual’s response to subsequent dual-tDCS. In addition, we demonstrated an increase in movement-related ERS in the α-band. The increase in α-ERS was correlated with the stroke survivors’ concurrent ARAT and FMA-UE.
Increased interhemispheric inhibition to ipsilesional M1 in stroke survivors was significantly modulated by dual-tDCS and predicted response to stimulation
In line with previous literature [53,54,55], we demonstrated a significant decrease in excitability in the ipsilesional hemisphere and a relatively excessive transcallosal inhibition from the contralesional to ipsilesional hemisphere at baseline, compared with controls, as evidenced by a decreased MEP amplitude and an increased iSP duration. The iSP ratio reflects the degree of asymmetry interhemispheric inhibition, with higher ratio representing more transcallosal inhibition from contralesional [left in controls] M1 toward ipsilesional [right in controls] M1. In addition, this interhemispheric inhibition asymmetry predicted response to dual-tDCS: the greater the inhibition from the contralesional towards ipsilesional M1, the less the ipsilesional M1 excitability could be modulated. This finding is consistent with the hypothesis that excessive transcallosal inhibition may further inhibit ipsilesional hemisphere excitability  and, therefore, decrease neuromodulatory response. This hypothesis is supported by animal studies, highlights the role of post-stroke inhibition in motor recovery . The iSP ratio was higher at baseline in stroke survivors than controls and predicted response to tDCS in stroke survivors, but not in controls. It may be, therefore, that while the underlying physiological processes this metric represents is modulated by tDCS in health and disease, it has a functionally important role in representing interhemispheric imbalance. However, we may have to consider the individual difference when applying brain stimulation based on the interhermispheric imbalance hypothesis [7, 38].
The ipsilesional SICI has been suggested to decrease within 3 months post-stroke compared to contralesional and healthy control SICI in a recent meta-analysis by McDonnell et.al . However, we found comparable ipsilesional SICI between subacute (2–4 weeks post-stroke) patients and healthy controls (Fig. 2d). The inconsistency may arise from different stroke severity (our patients had milder paresis than those with significant SICI decrease) and other stroke profiles (our patients were more homogenous as subacute and subcortical infarction) . Further studies are required to determine ipsilesional SICI changes in subacute stroke. The iSP reflects interhemispheric inhibition mediated by γ-Aminobutyric acid (GABA)-A and GABA-B receptors [53, 57], while SICI reflects intra-cortical inhibition mediated by GABA-A activity. Our stroke survivors showed more interhemispheric inhibition than controls, and this was modulated to a greater extent in stroke survivors than controls by dual-tDCS. Both GABA-A and GABA-B circuits have been reported to be modulated by unilateral tDCS [58,59,60], but none has been explored in dual-tDCS studies.
Initially decreased ipsilesional M1 α-ERS in stroke survivors increased after tDCS and correlated with motor function
The way in which dual-tDCS modulates sensorimotor rhythms in subacute stroke survivors has not yet been investigated from the literature. Here, we demonstrated that stroke survivors had decreased ERS in both the α- and β-bands in ipsilesional M1 compared with controls. In addition, tDCS-induced change in α-ERS was positively correlated with motor function (FMA-UE and ARAT). Since ERS is thought to represent post-movement cortical inhibition or fine movement control [61, 62] and increased during stroke recovery , tDCS may promote stroke motor recovery by modulating ipsilesional M1 circuit interneuronal activity following movement to obtain better motor control.
Some previous studies have demonstrated an increase in α-ERD amplitude during motor imagery after dual- or unilateral- anodal tDCS in chronic stroke [34, 64] and healthy adults [33, 65]. It is not clear why ERD changes were not observed here. One possibility may be derived from the differences between motor imagery and motor execution tasks (we adopted the latter only), which include different cortical excitability states, intracortical circuits and task complexity [66, 67].
Stroke survivors and controls showed similar high responsiveness to dual-tDCS
In our stroke survivors, the response rate and overall magnitude of the modulation effects in MEP were similar to controls (Fig. 3). Few tDCS studies have addressed responsiveness, but our cohort were much more responsive than those in a previous study using unilateral anodal tDCS . One previous dual-tDCS study  in six subacute stroke survivors without healthy controls demonstrated no MEP changes. Further studies are required to compare responsiveness across different tDCS montages and stages of post-stroke motor recovery.
There are some limitations to this study. First, we enrolled a homogenous group of stroke survivors for more consistent dual-tDCS effects in subacute stroke, which may limit the interpretation of our results in other stroke groups. Second, while we observed robust neurophysiological effects of a single session of tDCS on upper limbs, we did not seek to quantify its behavioural effects or concurrent influences on lower limbs. It is therefore difficult to fully describe the potential clinical significance of our findings, though as discussed above, changes in pathophysiological markers may be seen as robust proxy markers of function.
Our study demonstrated that dual-tDCS modulates neurophysiological markers in post-stroke survivors. These findings add to our understanding of the effects of tDCS in this understudied population. Our results, particularly that inter-hemispheric inhibition asymmetry may influence the effects of dual-tDCS, are important for the design of future clinical studies, suggesting that baseline electrophysiological measures should be therefore considered in translational applications of tDCS in stroke recovery.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Active motor threshold
Analysis of variance
Action research arm test
Extensor carpi radialis
Fugl-Meyer assessment-upper extremity scale
Ipsilateral silent period
Primary motor cortex
Motor evoked potentials
Medical Research Council scale of muscle strength
Modified Rankin Scale
National Institutes of Health Stroke Scale
Resting motor threshold
Short-interval intracortical inhibition
Transcranial direct current stimulation
Transcranial magnetic stimulation
Kubis N. Non-invasive brain stimulation to enhance post-stroke recovery. Front Neural Circuits. 2016;10:56.
Stagg CJ, Nitsche MA. Physiological basis of transcranial direct current stimulation. Neuroscientist. 2011;17(1):37–53.
Nitsche MA, Paulus W. Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation. J Physiol. 2000;527(Pt 3):633–9.
Takeuchi N, Izumi S. Noninvasive brain stimulation for motor recovery after stroke: mechanisms and future views. Stroke Res Treat. 2012;2012:584727.
Elsner B, Kugler J, Pohl M, Mehrholz J. Transcranial direct current stimulation (tDCS) for improving activities of daily living, and physical and cognitive functioning, in people after stroke. Cochrane Database Syst Rev. 2016;3:CD009645.
Calderon MAF, Miralles AN, Pimienta MJ, Estella JMG, Ledesma MJS. Analysis of the factors related to the effectiveness of transcranial current stimulation in upper limb motor function recovery after stroke: a systematic review. J Med Syst. 2019;43(3):69.
Di Pino G, Pellegrino G, Assenza G, Capone F, Ferreri F, Formica D, et al. Modulation of brain plasticity in stroke: a novel model for neurorehabilitation. Nat Rev Neurol. 2014;10(10):597–608.
Jo JY, Lee A, Kim MS, Park E, Chang WH, Shin Y-I, et al. Prediction of motor recovery using quantitative parameters of motor evoked potential in patients with stroke. Ann Rehabil Med. 2016;40(5):806.
Labruna L, Jamil A, Fresnoza S, Batsikadze G, Kuo M-F, Vanderschelden B, et al. Efficacy of anodal Transcranial direct current stimulation is related to sensitivity to Transcranial magnetic stimulation. Brain Stimulation. 2016;9(1):8–15.
Kidgell DJ, Goodwill AM, Frazer AK, Daly RM. Induction of cortical plasticity and improved motor performance following unilateral and bilateral transcranial direct current stimulation of the primary motor cortex. BMC Neurosci. 2013;14:64.
Karok S, Fletcher D, Witney AG. Task-specificity of unilateral anodal and dual-M1 tDCS effects on motor learning. Neuropsychologia. 2017;94:84–95.
Koyama S, Tanaka S, Tanabe S, Sadato N. Dual-hemisphere transcranial direct current stimulation over primary motor cortex enhances consolidation of a ballistic thumb movement. Neurosci Lett. 2015;588:49–53.
Vines BW, Cerruti C, Schlaug G. Dual-hemisphere tDCS facilitates greater improvements for healthy subjects’ non-dominant hand compared to uni-hemisphere stimulation. BMC Neurosci. 2008;9:103.
Karok S, Witney AG. Enhanced motor learning following task-concurrent dual transcranial direct current stimulation. PLoS One. 2013;8(12):e85693.
Tazoe T, Endoh T, Kitamura T, Ogata T. Polarity specific effects of transcranial direct current stimulation on interhemispheric inhibition. PLoS One. 2014;9(12):e114244.
Mordillo-Mateos L, Turpin-Fenoll L, Millan-Pascual J, Nunez-Perez N, Panyavin I, Gomez-Arguelles JM, et al. Effects of simultaneous bilateral tDCS of the human motor cortex. Brain Stimul. 2012;5(3):214–22.
McCambridge AB, Stinear JW, Byblow WD. ‘I-wave’ recruitment determines response to tDCS in the upper limb, but only so far. Brain Stimul. 2015;8(6):1124–9.
O'Shea J, Boudrias MH, Stagg CJ, Bachtiar V, Kischka U, Blicher JU, et al. Predicting behavioural response to TDCS in chronic motor stroke. NeuroImage. 2014;85:924–33.
Goodwill AM, Teo WP, Morgan P, Daly RM, Kidgell DJ. Bihemispheric-tDCS and upper limb rehabilitation improves retention of motor function in chronic stroke: a pilot study. Front Hum Neurosci. 2016;10:258.
Bolognini N, Vallar G, Casati C, Latif LA, El-Nazer R, Williams J, et al. Neurophysiological and behavioral effects of tDCS combined with constraint-induced movement therapy in poststroke patients. Neurorehabilitat Neural Repair. 2011;25(9):819–29.
Lindenberg R, Renga V, Zhu LL, Nair D, Schlaug G. Bihemispheric brain stimulation facilitates motor recovery in chronic stroke patients. Neurology. 2010;75(24):2176–84.
Lefebvre S, Dricot L, Laloux P, Gradkowski W, Desfontaines P, Evrard F, et al. Neural substrates underlying stimulation-enhanced motor skill learning after stroke. Brain. 2015;138(Pt 1):149–63.
Williams JA, Pascual-Leone A, Fregni F. Interhemispheric modulation induced by cortical stimulation and motor training. Phys Ther. 2010;90(3):398–410.
Straudi S, Fregni F, Martinuzzi C, Pavarelli C, Salvioli S, Basaglia N. tDCS and robotics on upper limb stroke rehabilitation: effect modification by stroke duration and type of stroke. Biomed Res Int. 2016;2016:5068127.
Del Felice A, Daloli V, Masiero S, Manganotti P. Contralesional cathodal versus dual transcranial direct current stimulation for decreasing upper limb spasticity in chronic stroke individuals: a clinical and neurophysiological study. J Stroke Cerebrovasc Dis. 2016;25(12):2932–41.
Dehem S, Gilliaux M, Lejeune T, Delaunois E, Mbonda P, Vandermeeren Y, et al. Effectiveness of a single session of dual-transcranial direct current stimulation in combination with upper limb robotic-assisted rehabilitation in chronic stroke patients: a randomized, double-blind, cross-over study. Int J Rehabil Res. 2018;41(2):138–45.
Wang QM, Cui H, Han SJ, Black-Schaffer R, Volz MS, Lee Y-T, et al. Combination of transcranial direct current stimulation and methylphenidate in subacute stroke. Neurosci Lett. 2014;569:6–11.
Tangwiriyasakul C, Verhagen R, WLC R, van MJAM P. Temporal evolution of event-related desynchronization in acute stroke: A pilot study. Clinical Neurophysiology. 2013;125:1112–20.
Gerloff C, Bushara K, Sailer A, Wassermann EM, Chen R, Matsuoka T, et al. Multimodal imaging of brain reorganization in motor areas of the contralesional hemisphere of well recovered patients after capsular stroke. Brain. 2006;129(Pt 3):791–808.
Shiner CT, Tang H, Johnson BW, McNulty PA. Cortical beta oscillations and motor thresholds differ across the spectrum of post-stroke motor impairment, a preliminary MEG and TMS study. Brain Res. 1629;2015:26–37.
Bartur G, Pratt H, Soroker N. Changes in Mu and Beta Amplitude of the EEG during Upper Limb Movement Correlate with Motor Impairment and Structural Damage in Subacute Stroke. Clinical Neurophysiol. 2019;130(9):1644–51.
Pfurtscheller G. Lopes da Silva FH. Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin Neurophysiol. 1999;110(11):1842–57.
Zich C, Harty S, Kranczioch C, Mansfield KL, Sella F, Debener S, et al. Modulating hemispheric lateralization by brain stimulation yields gain in mental and physical activity. Sci Rep. 2017;7(1):13430.
Ang KK, Guan C, Phua KS, Wang C, Zhao L, Teo WP, et al. Facilitating effects of transcranial direct current stimulation on motor imagery brain-computer interface with robotic feedback for stroke rehabilitation. Arch Phys Med Rehabil. 2015;96(3 Suppl):S79–87.
Cramer SC. Repairing the human brain after stroke: I. mechanisms of spontaneous recovery. Ann Neurol. 2008;63(3):272–87.
Lee S-H. Rehabilitation in subacute and chronic stage after stroke. 2017. In: Stroke revisited : diagnosis and treatment of ischemic stroke. Springer, Singapore. 1st edition. [351-360]. Available from: doi: https://doi.org/10.1007/978-981-10-1424-6.
Stinear C, Ackerley S, Byblow W. Rehabilitation is initiated early after stroke, but most motor rehabilitation trials are not: a systematic review. Stroke. 2013;44(7):2039–45.
Xu J, Branscheidt M, Schambra H, Steiner L, Widmer M, Diedrichsen J, et al. Rethinking interhemispheric imbalance as a target for stroke neurorehabilitation. Ann Neurol. 2019;85(4):502–13.
Stinear CM, Petoe MA, Byblow WD. Primary motor cortex excitability during recovery after stroke: implications for Neuromodulation. Brain Stimul. 2015;8(6):1183–90.
Verdino M, Dingman S. Two measures of laterality in handedness: the Edinburgh handedness inventory and the Purdue pegboard test of manual dexterity. Percept Mot Skills. 1998;86(2):476–8.
Bikson M, Grossman P, Thomas C, Zannou AL, Jiang J, Adnan T, et al. Safety of Transcranial direct current stimulation: evidence based update 2016. Brain Stimul. 2016;9(5):641–61.
Brunoni AR, Amadera J, Berbel B, Volz MS, Rizzerio BG, Fregni F. A systematic review on reporting and assessment of adverse effects associated with transcranial direct current stimulation. Int J Neuropsychopharmacol. 2011;14(8):1133–45.
Groppa S, Oliviero A, Eisen A, Quartarone A, Cohen LG, Mall V, et al. A practical guide to diagnostic transcranial magnetic stimulation: report of an IFCN committee. Clin Neurophysiol. 2012;123(5):858–82.
Lang N, Nitsche MA, Paulus W, Rothwell JC, Lemon RN. Effects of transcranial direct current stimulation over the human motor cortex on corticospinal and transcallosal excitability. Exp Brain Res. 2004;156(4):439–43.
Yang YR, Tseng CY, Chiou SY, Liao KK, Cheng SJ, Lai KL, et al. Combination of rTMS and treadmill training modulates corticomotor inhibition and improves walking in Parkinson disease: a randomized trial. Neurorehabil Neural Repair. 2013;27(1):79–86.
Chiou SY, Wang RY, Liao KK, Wu YT, Lu CF, Yang YR. Co-activation of primary motor cortex ipsilateral to muscles contracting in a unilateral motor task. Clin Neurophysiol. 2013;124(7):1353–63.
Kujirai T, Caramia MD, Rothwell JC, Day BL, Thompson PD, Ferbert A, et al. Corticocortical inhibition in human motor cortex. J Physiol. 1993;471:501–19.
Hsu YF, Huang YZ, Lin YY, Tang CW, Liao KK, Lee PL, et al. Intermittent theta burst stimulation over ipsilesional primary motor cortex of subacute ischemic stroke patients: a pilot study. Brain Stimul. 2012;6(2):166–74.
Tadel F, Baillet S, Mosher JC, Pantazis D, Leahy RM. Brainstorm: a user-friendly application for MEG/EEG analysis. Comput Intell Neurosci. 2011;2011:879716.
Hämäläinen M, Ilmoniemi R. Interpreting magnetic fields of the brain: minimum norm estimates. Med Biol Eng Comput. 1994;32(1):35–42.
Pfurtscheller G, Da Silva FHL. Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin Neurophysiol. 1999;110(11):1842–57.
Nettekoven C, Volz LJ, Leimbach M, Pool EM, Rehme AK, Eickhoff SB, et al. Inter-individual variability in cortical excitability and motor network connectivity following multiple blocks of rTMS. Neuroimage. 2015;118:209–18.
Palmer JA, Wheaton LA, Gray WA, Saltao da Silva MA, Wolf SL, Borich MR. Role of Interhemispheric cortical interactions in Poststroke motor function. Neurorehabil Neural Repair. 2019;33(9):762–74.
McDonnell MN, Stinear CM. TMS measures of motor cortex function after stroke: a meta-analysis. Brain Stimul. 2017;10(4):721–34.
Takechi U, Matsunaga K, Nakanishi R, Yamanaga H, Murayama N, Mafune K, et al. Longitudinal changes of motor cortical excitability and transcallosal inhibition after subcortical stroke. Clin Neurophysiol. 2014;125(10):2055–69.
Clarkson AN, Huang BS, Macisaac SE, Mody I, Carmichael ST. Reducing excessive GABA-mediated tonic inhibition promotes functional recovery after stroke. Nature. 2010;468(7321):305–9.
Ziemann U, Reis J, Schwenkreis P, Rosanova M, Strafella A, Badawy R, et al. TMS and drugs revisited 2014. Clin Neurophysiol. 2015;126(10):1847–68.
Tremblay S, Beaule V, Lepage JF, Theoret H. Anodal transcranial direct current stimulation modulates GABAB-related intracortical inhibition in the M1 of healthy individuals. Neuroreport. 2013;24(1):46–50.
Antal A, Terney D, Kuhnl S, Paulus W. Anodal transcranial direct current stimulation of the motor cortex ameliorates chronic pain and reduces short intracortical inhibition. J Pain Symptom Manag. 2010;39(5):890–903.
Hummel F, Celnik P, Giraux P, Floel A, Wu WH, Gerloff C, et al. Effects of non-invasive cortical stimulation on skilled motor function in chronic stroke. Brain. 2005;128(Pt 3):490–9.
Pfurtscheller G. Induced oscillations in the alpha band: functional meaning. Epilepsia. 2003;44(Suppl. 12):2–8.
Hummel F, Andres F, Altenmuller E, Dichgans J, Gerloff C. Inhibitory control of acquired motor programmes in the human brain. Brain. 2002;125(Pt 2):404–20.
Laaksonen K, Kirveskari E, Mäkelä JP, Kaste M, Mustanoja S, Nummenmaa L, et al. Effect of afferent input on motor cortex excitability during stroke recovery. Clin Neurophysiol. 2012;123(12):2429–36.
Kasashima Y, Fujiwara T, Matsushika Y, Tsuji T, Hase K, Ushiyama J, et al. Modulation of event-related desynchronization during motor imagery with transcranial direct current stimulation (tDCS) in patients with chronic hemiparetic stroke. Exp Brain Res. 2012;221(3):263–8.
Vecchio F, Pellicciari MC, Miraglia F, Brignani D, Miniussi C, Rossini PM. Effects of transcranial direct current stimulation on the functional coupling of the sensorimotor cortical network. Neuroimage. 2016;140:50–6.
Meng HJ, Pi YL, Liu K, Cao N, Wang YQ, Wu Y, et al. Differences between motor execution and motor imagery of grasping movements in the motor cortical excitatory circuit. PeerJ. 2018;6:e5588.
Avanzino L, Gueugneau N, Bisio A, Ruggeri P, Papaxanthis C, Bove M. Motor cortical plasticity induced by motor learning through mental practice. Front Behav Neurosci. 2015;9:105.
Nuzum ND, Hendy AM, Russell AP, Teo WP. Measures to predict the individual variability of Corticospinal responses following Transcranial direct current stimulation. Front Hum Neurosci. 2016;10:487.
We gratefully thank Wen-Yung Sheng for providing statistical advice. We thank the staff of the Clinical Research Core Laboratory and the Stroke Registry of the Taipei Veterans General Hospital for providing experimental facilities and collecting registry data, respectively. We also thank the Pervasive Artificial intelligence Research Labs in Taiwan.
This work was supported by the Taiwan Ministry of Science and Technology [grant numbers 106–2314-B-075-019-MY3, 108–2634-F-008-003]; the Taipei Veterans General Hospital [grant numbers V107C-064, V108C-052, V109C-034]; the Taiwan Ministry of Education [Academic Strategic Alliance Project between Taiwan and Oxford University]; and the National Yang-Ming University [Brain Research Center, from The Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Taiwan Ministry of Education]. CJS holds a Sir Henry Dale Fellowship, funded by the Wellcome Trust and the Royal Society (102584/Z/13/Z). The work was supported by the NIHR Oxford Health Biomedical Research Centre. The Wellcome Centre for Integrative Neuroimaging is supported by core funding from the Wellcome Trust (203139/Z/16/Z).
Ethics approval and consent to participate
The study was approved by the ethics committee of Taipei Veterans General Hospital. All subjects provided written informed consent before participation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kuo, IJ., Tang, CW., Tsai, YA. et al. Neurophysiological signatures of hand motor response to dual-transcranial direct current stimulation in subacute stroke: a TMS and MEG study. J NeuroEngineering Rehabil 17, 72 (2020). https://doi.org/10.1186/s12984-020-00706-1