Alexander B. Remsik1,2,3*
Klevest Gjini1,4
Leroy Williams Jr.1,5,6
Peter L. E. van Kan2,7Shawna Gloe1Erik Bjorklund1,8
Cameron A. Rivera1Sophia Romero1 Brittany M. Young1,2,7,8,9Veena A. Nair1Kristin E. Caldera1
Justin C. Williams10,11Vivek Prabhakaran1,4,7,9,12,13- 1Department of Radiology, University of Wisconsin–Madison, Madison, WI, United States
- 2Department of Kinesiology, University of Wisconsin–Madison, Madison, WI, United States
- 3Institute for Clinical and Translational Research, University of Wisconsin–Madison, Madison, WI, United States
- 4Department of Neurology, University of Wisconsin–Madison, Madison, WI, United States
- 5Department of Educational Psychology, University of Wisconsin–Madison, Madison, WI, United States
- 6Center for Women’s Health Research, University of Wisconsin–Madison, Madison, WI, United States
- 7Neuroscience Training Program, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States
- 8Clinical Neuroengineering Training Program, University of Wisconsin–Madison, Madison, WI, United States
- 9Medical Scientist Training Program, University of Wisconsin School of Medicine and Public Health, Madison, WI, United States
- 10Department of Biomedical Engineering, University of Wisconsin–Madison, Madison, WI, United States
- 11Department of Neurological Surgery, University of Wisconsin–Madison, Madison, WI, United States
- 12Department of Psychiatry, University of Wisconsin–Madison, Madison, WI, United States
- 13Department of Psychology, University of Wisconsin–Madison, Madison, WI, United States
Stroke is a leading cause of acquired long-term upper extremity motor disability. Current standard of care trajectories fail to deliver sufficient motor rehabilitation to stroke survivors. Recent research suggests that use of brain-computer interface (BCI) devices improves motor function in stroke survivors, regardless of stroke severity and chronicity, and may induce and/or facilitate neuroplastic changes associated with motor rehabilitation. The present sub analyses of ongoing crossover-controlled trial NCT02098265 examine first whether, during movements of the affected hand compared to rest, ipsilesional Mu rhythm desynchronization of cerebral cortical sensorimotor areas [Brodmann’s areas (BA) 1-7] is localized and tracks with changes in grip force strength. Secondly, we test the hypothesis that BCI intervention results in changes in frequency-specific directional flow of information transmission (direct path functional connectivity) in BA 1-7 by measuring changes in isolated effective coherence (iCoh) between cerebral cortical sensorimotor areas thought to relate to electrophysiological signatures of motor actions and motor learning. A sample of 16 stroke survivors with right hemisphere lesions (left hand motor impairment), received a maximum of 18–30 h of BCI intervention. Electroencephalograms were recorded during intervention sessions while outcome measures of motor function and capacity were assessed at baseline and completion of intervention. Greater desynchronization of Mu rhythm, during movements of the impaired hand compared to rest, were primarily localized to ipsilesional sensorimotor cortices (BA 1-7). In addition, increased Mu desynchronization in the ipsilesional primary motor cortex, Post vs. Pre BCI intervention, correlated significantly with improvements in hand function as assessed by grip force measurements. Moreover, the results show a significant change in the direction of causal information flow, as measured by iCoh, toward the ipsilesional motor (BA 4) and ipsilesional premotor cortices (BA 6) during BCI intervention. Significant iCoh increases from ipsilesional BA 4 to ipsilesional BA 6 were observed in both Mu [8–12 Hz] and Beta [18–26 Hz] frequency ranges. In summary, the present results are indicative of improvements in motor capacity and behavior, and they are consistent with the view that BCI-FES intervention improves functional motor capacity of the ipsilesional hemisphere and the impaired hand.
Introduction
Stroke is a leading cause of acquired upper extremity (UE) motor disability and many survivors are left with persistent upper extremity motor impairments requiring rehabilitation (Benjamin et al., 2017, 2019). In the United States alone, on average, every 40 seconds someone suffers a stroke (Benjamin et al., 2019). About 6 months after stroke insult, approximately half of stroke survivors continue to suffer residual motor deficit (Benjamin et al., 2017). Stroke burden on the United States economy, by 2050, is expected to be above $2.2 trillion (Benjamin et al., 2017). Despite advances in acute stroke care, the estimated direct and indirect costs of stroke continue to escalate and are disproportionately associated with long-term care and rehabilitation (Benjamin et al., 2017, 2019).
Electroencephalogram-based brain-computer interface (BCI) intervention has been proposed as a novel intervention tool (Pfurtscheller and Berghold, 1989; Pfurtscheller et al., 1997, 2005; Neuper and Pfurtscheller, 2001; Wolpaw et al., 2002; Leuthardt et al., 2004; Schalk et al., 2004; Birbaumer, 2006; Birbaumer et al., 2006; Wilson et al., 2009; Prasad et al., 2010; Schalk and Mellinger, 2010; Bundy et al., 2012; Babaiasl et al., 2016; Irimia et al., 2016; Young et al., 2016; Mazrooyisebdani et al., 2018; Mohanty et al., 2018a, b; Remsik et al., 2018; Remsik A. B. et al., 2019) capable of enhancing motor recovery post-stroke. The neural mechanisms underlying BCI’s effect on motor rehabilitation, either through neural plasticity or otherwise, are not well understood. BCIs are a promising supplement to existing means of neurorehabilitation but may also function as tools that provide insight into the sensorimotor processes underlying motor function and motor learning in either healthy or stroke-lesioned brains. Damage to the input (afferent) or output (efferent) pathways of the sensorimotor system creates a demand for reorganization of existing neural network functions (with respect to completion of behavioral goals) (Nudo, 2011; Nudo and McNeal, 2013). BCIs may induce or facilitate neuroplasticity by strengthening such connections between brain areas (Song et al., 2014, 2015; Young et al., 2014c; Biasiucci et al., 2018; Remsik A. B. et al., 2019).
EEG activity recorded from sensorimotor cortices (BA 1-7) of each hemisphere desynchronizes with imagined and attempted movements, and preparation of movement. This phenomenon is known as event-related desynchronization (ERD) (Pfurtscheller and Berghold, 1989; Pfurtscheller et al., 1997, 2005; Pfurtscheller, 1999; Pfurtscheller and Lopes da Silva, 1999; Neuper and Pfurtscheller, 2001; Pineda, 2005; Neuper et al., 2006; Nam et al., 2011). Specific frequency bands are associated with components of event-related motor behaviors (Pfurtscheller et al., 1997; Neuper and Pfurtscheller, 2001; Pineda, 2005; Nam et al., 2011; Remsik A. B. et al., 2019). When healthy individuals plan and execute purposeful movements, Mu rhythms of the contralateral sensorimotor cortices are desynchronized and attenuated (Pfurtscheller et al., 1997) and increased presence of Beta rhythm ERD is associated with motor command output and control (Pfurtscheller and Lopes da Silva, 1999; Pineda, 2005; Nam et al., 2011). ERD and event related synchronicity (ERS) were key components in the use and development of early EEG-based BCI motor rehabilitations (Wolpaw et al., 1991, 2002; Pfurtscheller et al., 1997, 2005; Pfurtscheller, 1999; McFarland et al., 2000; Neuper and Pfurtscheller, 2001; Leuthardt et al., 2004; Schalk et al., 2004; Neuper et al., 2006). Thus, ERD/ERS provide measures of volitional movement-related brain activity that can be decoded by a BCI to control a device, such as a prosthetic limb, or an output command, such as a cursor on a computer screen (Pfurtscheller et al., 1997; Leuthardt et al., 2004; Neuper et al., 2006; Wilson et al., 2009).
In human subjects, at frequencies around 10–20 Hz, Mu and Beta sensorimotor rhythms (SMRs) are recorded exclusively over cortical sensorimotor areas (BA 1-7) (Pfurtscheller et al., 1997; Birbaumer, 2006; Birbaumer et al., 2006). Two basic strategies, motor imagery and attempted movement (Wolpaw et al., 1991, 2002; Wilson et al., 2009; Nam et al., 2011; Ortner et al., 2012; Wilson et al., 2012; Song et al., 2014, 2015; Young et al., 2014a,b,c,d, 2015, 2016; Nair et al., 2015; Irimia et al., 2016; Remsik et al., 2016, 2018; Mazrooyisebdani et al., 2018; Mohanty et al., 2018a; Remsik A. B. et al., 2019) as well as various therapeutic adjuvant approaches (e.g., BCI-FES) (Biasiucci et al., 2018) have been introduced for motor rehabilitation in stroke patients. Both approaches record SMRs as input signals (electrophysiological recordings by the EEG cap) to the BCI from overlapping neural architecture. This protocol was designed to utilize attempted hand movements during the intervention following the logic that a motor therapy designed to rehabilitate volitional motor function of the affected UE should utilize voluntary attempted movements of an impaired hand in a continuous effort to improve the participant’s UE capacity and performance.
The present sub analyses of ongoing crossover-controlled trial NCT02098265 examine first whether, during movements of the affected hand compared to rest, ipsilesional Mu rhythm desynchronization (Pfurtscheller and Berghold, 1989; Pfurtscheller et al., 1997; Neuper and Pfurtscheller, 2001; Pineda, 2005; Nam et al., 2011) of cerebral cortical sensorimotor areas (Brodmann’s areas (BA) 1-7) (Brodmann, 1909) is localized and tracks with changes in grip force strength. Secondly, we test the hypothesis that BCI intervention results in changes in frequency-specific directional flow of information transmission (direct path functional connectivity) in BA 1-7 by measuring changes in isolated effective coherence, iCoh (Pascual-Marqui et al., 2014; Kitaura et al., 2017), between cerebral cortical sensorimotor areas thought to relate to electrophysiological signatures of motor actions and motor learning. Brief reports of some of the results have been published (Remsik A. et al., 2019; Remsik A. et al., 2020; Remsik A. B. et al., 2020).
Materials and Methods
Participant Population
Participants were recruited from the greater Madison, Wisconsin, United States area as part of an on-going cross-over controlled prospective randomized rehabilitation study investigating interventional BCI in upper extremity (UE) motor function impairment resulting from stroke. This study is approved by the University of Wisconsin Health Sciences Institutional Review Board (Study ID 20150469). All subjects provided informed written consent upon enrollment.
Study Design and Procedures
The BCI system and intervention procedures were consistent with previously published works, Wilson et al. (2009, 2012) and Song et al. (2014) using BCI2000 software (Schalk et al., 2004) version 2 with in-house modifications for input from a 16-channel EEG cap and amplifier (Guger Technologies) and integration with a ball and target gaming visual display (Wilson et al., 2009) and functional electrical stimulation (FES) adjuvant (Popovic et al., 2009; Takahashi et al., 2012; Biasiucci et al., 2018). A more detailed and complete description of the session and run procedures used in data acquisition and analysis reported in this manuscript, as well as a description of the functional electrical stimulation protocol and justification for how this approach further improves the effects of BCI have been published previously (Remsik A. B. et al., 2019).
Sixteen right hemisphere stroke survivors (8 female, age = 62.5 ± 12.7 years [mean ± SD]) (Table 1) participated in 9–15 sessions (13.9 ± 1.28 [mean ± SD]; following removal of artifact EEG data: 11.1 ± 3.87 [mean ± SD]) with the BCI. Participants executed hand movements in response to visual cues displayed on a computer screen concordantly with corresponding audio instructions (e.g., Left, Right, Rest). One session (or run) of the cursor and target task consisted of 10 trials, or attempts, during which the participant attempted to drive a virtual cursor across the display screen into the target space using voluntary modulation of their sensorimotor rhythms (i.e., satisfaction of BCI classifier) before the trial timed out. BCI performance was measured as the number of successful trials (i.e., hits compared to misses) “out of a possible 10” attempts in any one run. This metric in included in Table 1 as “Average BCI Performance.”
Table 1. Demographics of n = 16 stroke survivors.
The BCI classifier was defined in “screening” sessions (Pre and Post BCI intervention). Screening sessions contained two runs, each consisting of 15 trials for rest, left-hand, and right-hand movements (i.e., 5 trials for each of the three conditions, the order of trials in a run was random). Only the rest and left-hand movement trials are considered here.
From the larger NCT02098265 cohort, only participants with left hand motor impairment were selected for these sub analyses in an effort to control for homogeneity of stroke-related upper extremity impairment, and language or communicative deficits that might interfere with comprehension of the BCI task requirements or task execution.
Preprocessing of the Scalp EEG Data
Various EEG-signal processing techniques were used to estimate and further evaluate the spectral perturbations recorded over the course of BCI intervention. EEG data analysis and statistics were consistent with methods detailed previously (Remsik A. B. et al., 2019).
EEG electrodes were positioned according to the standard 10-20 system, grounded to Fpz, and referenced to an ear electrode placed on the back of the participant’s right ear. Electrodes in the cap arrangement are highlight in purple in Figure 1 and contain electrodes F5, F6, FC1, FC2, C5, C3, Cz, C4, C6, CP1, CP2, P5, P4, Pz, P4, P6, A2. Inputs to the BCI electrodes over the sensorimotor cortices including, C3, C4, and Cz, were recorded in every session (Figure 1). EEG preprocessing steps include: (1) Digital filtering with a high-pass filter cutoff of 4 Hz, and a low-pass filter cutoff of 30 Hz. (2) Extraction and grouping of trials according to condition (rest, left/impaired hand movement, movement type (attempted), and the screening phase (Pre, Post intervention). This resulted in 10 trials for each of condition/movement/screening phase combinations. (3) Identification and repair of bad (noisy) channels (via spline interpolation), followed by the removal of the three trials with the highest variance. The channel was identified as bad (noisy, poor connection, etc.) if the variance was <10 or >250 μV2 in more than three trials. The units of the variance were those of the data squared: as the EEG data units were in μV, the variance units were squared μV. If more than four channels were identified as bad, the data for that session were removed from further analysis. At session level, this step resulted in 28 s of EEG data (7 trials × 4 s) for each condition/movement/screening phase combination set. (4) An average-reference montage was applied to the data (i.e., re-referencing from the original monopolar recordings). (5) Analyses with Fourier transforms computed using a multi-taper method (Thomson, 1982; Mitra and Pesaran, 1999) at a 0.25 Hz resolution. This resulted in estimates of absolute spectral power sampled for every 1 Hz bin during the interval of 4–30 Hz, and cross-spectral density. The trial length was 4 s and the resolution of Fourier transforms was 1/4 = 0.25 Hz.
Figure 1. BCI cap array. International 10–20 system for standardized EEG electrode locations on the head: C, central; P, parietal; T, temporal; F, frontal; Fp, pre-frontal; O, occipital. Electrode array and cap arrangement for all n = 16 participants is indicated in purple against the standard 10–20 electrode array head map. AFz (yellow) is the ground electrode. A2 (blue) electrode is the reference electrode, placed on the back of the right ear. Please note, electrode arrangement is designed for increased density around cortical sensorimotor areas.
Use of Brodmann Areas
To better illustrate the cerebral cortical sensorimotor areas represented by the co-registration of the 16-channel spectral EEG recordings to a three-dimensional head model space, we implemented ten regions of interest (ROIs). These ROIs are based on Brodmann’s areas 1–7 (Brodmann, 1909), and coincide with the electrode placement in the standardized EEG cap (Figure 2). Regions of interests were defined as follows: ROI #1 = left BA 1-3 (primary somatosensory cortices); ROI #2 = left BA 4 (primary motor cortex); ROI #3 = left BA 5 (somatosensory association area); ROI #4 = left BA 6 (premotor cortices); ROI #5 = left BA 7 (visuomotor coordination area); ROI #6 = right BA 1-3 (primary somatosensory cortices); ROI #7 = right BA 4 (primary motor cortex); ROI #8 = right BA 5 (somatosensory association area); ROI #9 = right BA 6 (premotor cortices); ROI #10 = right BA 7 (visuomotor coordination area).
Figure 2. EEG electrode locations (n = 16, blue spheres) are shown superimposed on a generic head model (left, middle, and right panels show top/bottom, back/front, and left/right views, respectively). BA, Brodmann’s areas. Regions of interests (n = 10) are defined as follows: ROI #1 = left BA 1-3 (primary somatosensory cortices, light purple); ROI #2 = left BA 4 (primary motor cortex, light blue); ROI #3 = left BA 5 (somatosensory association area, cyan); ROI #4 = left BA 6 (premotor cortices, light blue); ROI #5 = left BA 7 (visuomotor coordination area, cyan); ROI #6 = right BA 1-3 (primary somatosensory cortices, green); ROI #7 = right BA 4 (primary motor cortex, orange); ROI #8 = right BA 5 (somatosensory association area, red); ROI #9 = right BA 6 (premotor cortices, orange); ROI #10 = right BA 7 (visuomotor coordination area, red).
Source Localization (e/sLORETA)
Exact Low Resolution Electromagnetic Brain Tomography (eLORETA) is a weighted L2 minimum-norm distributed source localization algorithm used for the estimation of three-dimensional current density in the brain from the measured scalp EEG spectral data (Pascual-Marqui et al., 1994; Pascual-Marqui, 2007). In the eLORETA (unit: μA/mm2) implementation, computations are made in a realistic head model using the MNI152 template, with the three-dimensional solution space restricted to cortical gray matter, as determined by the probabilistic Talairach atlas. The solution space consists of gray matter of the hemispheres and hippocampus (6,239 voxels at 5 mm grid in the MNI coordinate system). The specific frequency band cross-spectra (frequency-domain) obtained from the average-referenced scalp potential data, were the inputs for source localization. In summary, the cortical current density estimates from the 16 EEG electrode signals are assigned to an atlas space to obtain regions-of-interest (ROIs) based signals representing the source signals reflected in the recorded scalp EEG signal characteristics.
Overall, segments of clean EEG data (left/impaired hand movements and rest), separately for Pre- and Post-BCI intervention screening sessions, were used for the computation of cross-spectra (Mu [8–12 Hz] and Beta [18–26 Hz] bands) (Remsik A. B. et al., 2019). Averages of cross-spectra (1 average per each subject (n = 16), separately for movements and rest trials) were computed, and then eLORETA estimates of Mu [8–12 Hz] band power at 6239 cortical locations/voxels were obtained, using subject-wise normalization of the estimates (which means dividing the eLORETA current density estimate values by the grand-average over all voxels and frequency bands evaluated). Finally, the derived (movement – rest) eLORETA normalized estimates were used in paired-sample statistics (using statistical non-parametric mapping approach (SnPM) as implemented in the s/eLORETA software; voxel-based comparison).
Source-localized ERD in cortical space was chosen as a means of representing scalp-recorded brain signal changes during attempted hand movements. Activity recorded on the scalp is a representation of multiple source generators within cortical space. Extrapolating from the scalp to cortical space provides the justification for examining Mu ERD, an accepted brain signal associated with planning and executing movement (Pineda, 2005), and offers a means for visualizing surface EEG recordings in three dimensions. In fact, estimated source signals are a better representation of the underlying cortical generators that produce the activity recorded by scalp electrodes (Yuan et al., 2008).
Isolated Effective Coherence
Isolated effective coherence (iCoh) is a metric for frequency-specific directional flow of information transmission and offers a means of assessing direct paths of intracortical causal information flow of oscillatory activity. iCoh is based on formulating a multivariate autoregressive model from time series measurements and calculating the corresponding partial coherences after setting all irrelevant connections to zero, according to Pascual-Marqui et al. (2014). From the spectral density matrix (bandwidth of 4–30 Hz, with spectral resolution of 0.25 Hz, including Mu [8–12 Hz] and Beta [18–26 Hz] ranges) obtained from estimated signals in the selected ROIs, the partial coherences between any pair of nodes/ROIs can be calculated. The t-statistics was performed for iCoh values between post and pre-intervention screening session trials of the impaired (left) hand movement with threshold set at p = 0.05 (t = 2.13, uncorrected for multiple comparisons).
Hand Grip Function
Hand grip strength was assessed with a dynamometer (Boissy et al., 1999). Participants were asked to squeeze the spring-loaded dynamometer lever as hard as possible (i.e., maximal single hand grasp) with their entire hand, and then release. Three trials were performed with the affected hand and the average of the three trials was recorded as a handgrip score in pounds (lb).
Statistical Analysis
Statistical differences between post- and pre-intervention cortical estimates were tested using paired-sample t-statistics as part of the statistical non-parametric mapping approach (SnPM) implemented in the s/eLORETA software (voxel-based). Paired-sample t-statistics were also used to test differences in the ROI estimates (power and connectivity). Finally, the Pearson correlation coefficient was used to quantify the correlation between change in hand grip strength and change in Mu rhythm desynchronization (Post vs. Pre BCI intervention). The p values reported are uncorrected for multiple comparisons.
Results
Localization of Mu ERD Changes Following Brain-Computer Interface Intervention
Exact Low Resolution Electromagnetic Brain Tomography (eLORETA) was used to estimate the three-dimensional current density in the brain from the measured scalp EEG spectral data. In the eLORETA implementation, computations are made in a realistic head model using the MNI152 template, with the three-dimensional solution space restricted to cortical gray matter, as determined by the probabilistic Talairach atlas.
The resulting distributions of ERD and ERS heatmaps of Figure 3 illustrate the distribution of ERD and ERS changes from Pre to Post BCI intervention in the n = 16 stroke survivors represented in Figure 4. Considering the distribution of ERS (red) and ERD (blue) voxels in Figure 3, it is apparent that as a result of BCI intervention, participants realized an increase of Mu ERD in ipsilesional hemisphere voxels of sensorimotor cortices during task performance (i.e., hand grasping) with the affected (left) hand compared to rest.
Figure 3. (A) Localization of Mu rhythm desynchronization of cerebral cortical sensorimotor areas during movement of the affected hand compared to rest. Computations are made in a realistic head model using the MNI152 template, with the three-dimensional solution space restricted to cortical gray matter, as determined by the probabilistic Talairach atlas. The specific frequency band cross-spectra (frequency-domain) obtained from the average-reference potential data were the inputs for source localization. The resulting distribution of ERD and ERS heatmaps illustrate the distribution of ERD and ERS changes from Pre to Post BCI intervention in the n = 16 stroke survivors. (B) Obtained eLORETA of Mu band power estimates at 6239 cortical locations (i.e., voxels) normalized across subjects.
Figure 4. eLORETA estimates PRE (blue bars) and POST (orange bars) BCI intervention. eLORETA estimates between the two conditions [Movement – Rest] were based on measured EEG spectral power data in the Mu band [8–12 Hz] for each of the 10 regions of interest (ROI #1-10). P values refer to statistical group comparisons of Post vs. Pre estimates using a ROI-based one-tailed t-test. The p values reported are uncorrected for multiple comparisons.
Figure 4 shows pre-BCI intervention (blue) and post-BCI intervention (orange) computed eLORETA cortical spectral power estimates for the Mu [8–12 Hz] band (i.e., 1 average per each of the 16 subjects included in these sub analyses, separately for movement (left) and rest trials) for each of the ten regions of interest and assumed underlying Brodmann’s Areas (Figure 2). ERD increases (Post vs. Pre BCI intervention) were observed in regions of interest ROI #6 through ROI #9. The largest increases in Mu ERD were observed for regions of interest ROI #7 and ROI #8, which correspond to ipsilesional primary motor cortex (BA 4) (paired-samples t-test: p = 0.11) and ipsilesional somatosensory association area (BA 5) (paired-samples t-test: p = 0.056), respectively (Figure 4).
Correlation Between Mu Desynchronization and Hand Grip Strength
Figure 5A plots change in hand grip strength as a function of change in Mu rhythm desynchronization (Post vs. Pre BCI intervention) in the primary motor cortex of the ipsilesional hemisphere (ROI #7). Each data point represents an individual subject. Please note that negative numbers in the change in Mu rhythm desynchronization (i.e., Post-Pre), represent higher desynchronization in the Post phase. That is, larger negative numbers subtracted from smaller negative numbers will yield negative numbers. Improved hand grip function is positively correlated with increased Mu rhythm desynchronization (Post-Pre) in the ipsilesional primary motor cortex (Pearson’s correlation coefficient r = 0.435, p = 0.046).
Figure 5. (A) Increase in Mu rhythm desynchronization in ROI #7 (i.e., ipsilesional primary motor cortex), Post vs. Pre BCI intervention, is significantly correlated with increase in hand grip strength (Pearson correlation coefficient r = 0.435, p = 0.046). Two data points overlap at (−0.20, 0.00). (B) Voxel-wise correlation of Hand Grip change with the Mu [8–12 Hz] rhythm desynchronization. The heat maps show colored r values (threshold set at r = 0.5). Scale: −0.5 (blue) to 0 (white) to + 0.5 (red). Note that stronger ipsilesional (R) Mu rhythm desynchronization (negative values) (Post-Pre) correlates with larger increase in hand grip strength.
A voxel-wise correlation approach (refer to Methods section) was used to represent this relationship in brain atlas space (Figure 5B). The heatmaps show colored r values (e.g., Pearson’s correlation coefficient values) and indicate that the correlation between improvement in hand grip strength and increase in Mu rhythm r values is strongest in voxels representing ipsilesional cerebral cortical sensorimotor areas.
In summary, the above results suggest that for the stroke-lesioned hemisphere, BCI intervention facilitates increased Mu desynchronization associated with movement rehabilitation of the impaired upper extremity.
Frequency-Specific Directional Flow of Information Transmission in BA 1-7
Direct paths of Pre to post BCI intervention frequency-specific flow of information transmission between cerebral cortical sensorimotor areas were analyzed by computing isolated effective coherence (iCoh) developed by Pascual-Marqui et al. (2014).
Figure 6 plots iCoh values as a function of spectral frequency, Pre (blue) and Post (red) BCI intervention, to and from the 10 regions of interest (ROI #1-10) for the affected (left) hand. Mean values (Pre, Post) were compared via paired-sample t-statistics, thresholded at t = 2.13 [df = 15], p = 0.05 uncorrected).
Figure 6. Isolated effective coherence (iCoh) – Frequency Specific Directional Flow of Information Transmission (Left Hand, Impaired Movement). Isolated effective coherence (iCoh), Pre (blue) and Post (red) BCI intervention. Vertical axis: 0 to 1. Horizontal axis: spectral frequency, 4–30 Hz. Columns are the senders, rows are receivers (i.e., FROM the top row of ROIs TO the horizontal columns of ROIs).
Figure 7 summarizes t-statistics for the iCoh data shown in Figure 6. Statistically significant increases and decreases in iCoh at relative frequencies from Pre to Post BCI intervention between pairs of ROIs are plotted in red and blue, respectively. Significant increases in iCoh from Pre to post BCI intervention were seen going from ROI #7 to ROI #9, suggesting more causal informational flow is going from ROI #7 (ipsilesional primary motor cortex) to ROI #9 (ipsilesional premotor cortices) Post compared to Pre BCI intervention. Significant decreases in iCoh were seen going from ROI #8 to ROI #5, suggesting less causal informational flow is going from ROI #8 (ipsilesional somatosensory association area) to ROI #5 (contralesional visuomotor coordination area) Post compared to Pre BCI intervention.
Figure 7. Plots of t-statistics for isolated effective coherence (iCoh) values Post vs. Pre BCI intervention of left (i.e., impaired hand) movement trials. Vertical axis: 0 to 1. Horizontal axis: spectral frequency: 4–30 Hz. Columns are the senders, rows are receivers (i.e., FROM the top row of ROIs TO the horizontal columns of ROIs). Red portions of the iCoh between two ROIs denote coherence increases at relative frequencies from Pre to Post, whereas blue portions denote significant decreases in iCoh from Pre to Post between two ROIs. Mean values (Pre, Post) compared via paired-sample t-statistics, with threshold set at t = 2.13 [df = 15], p = 0.05, uncorrected for multiple comparisons.
Signals in both Mu [8–12 Hz] and Beta [18–26 Hz] frequency ranges served as input commands in the present BCI design (Remsik A. B. et al., 2019). Figures 8, 9 summarize for Mu [8–12 Hz] and Beta [18–26 Hz] signal ranges, respectively, changes in functional connectivity of sensorimotor cerebral cortices during attempted movement of the impaired upper extremity as a result of BCI intervention. Each of the panels of Figures 8, 9 compares [Pre vs Post] iCoh values in Mu (Figure 8) and Beta (Figure 9) frequency ranges and represents a direct path of intracortical causal information flow of oscillatory activity between two ROIs. Statistically significant differences between Pre (blue bars) and Post (red bars) iCoh values between pairs of ROIs are represented by surrounding boxes with blue denoting significant decreases in iCoh, and red boxes signifying significant increases in iCoh.
Figure 8. Plots of iCoh between two given ROIs in the Mu frequency band [8–12 Hz], Pre and Post BCI intervention. Vertical axis: Coherence value –1 to 1. Horizontal axis: Pre (blue) and Post (red). Columns are senders, rows are receivers (i.e., FROM the top row of ROIs TO the horizontal columns of ROIs). Blue surrounding boxes indicate significant decrease in iCoh while red surrounding boxes denote significant increases in iCoh (Post vs. Pre).
Figure 9. Plots of iCoh between two given ROIs in the Beta frequency band [18–26 Hz], Pre and Post BCI intervention, Vertical axis: Coherence value –1 to 1. Horizontal axis: Pre (blue) and Post (red). Columns are senders, rows are receivers (i.e., FROM the top row of ROIs TO the horizontal columns of ROIs). Blue surrounding boxes indicate significant decrease in iCoh while red surrounding boxes denote significant increases in iCoh (Post vs. Pre).
In the Mu [8–12 Hz] frequency band, iCoh values significantly decreased Pre vs Post BCI intervention from ROI #4 (contralesional premotor area) to ROI #8 (ipsilesional somatosensory association area) (Figure 8). Mu band iCoh values significantly increased Pre vs Post BCI intervention from ROI #7 (ipsilesional primary motor) to ROI #9 (ipsilesional premotor cortices).
In the Beta [18–26 Hz] frequency band, iCoh significantly decreased Pre vs. Post BCI intervention from ROI #8 (ipsilesional somatosensory association area) to ROI #5 (contralesional visuo-motor coordination area) (Figure 9). Beta band iCoh values significantly increased Pre vs. Post BCI intervention from ROI #5 (contralesional visuo-motor coordination area) to ROI #9 (ipsilesional premotor cortices), from ROI #6 (ipsilesional primary somatosensory areas) to ROI #10 (ipsilesional visuo-motor coordination area), and from ROI #7 (ipsilesional primary motor cortex) to ROI #9 (ipsilesional premotor cortices) (Figure 9).
Discussion
We have investigated potential neural substrates underlying the effect of BCI intervention on motor recovery in stroke survivors. Here we report, for grasping movements of the affected hand compared to rest, significantly greater desynchronization of sensorimotor rhythms in the Mu frequency band [8–12 Hz] recorded via scalp electrodes in the ipsilesional hemisphere, indicating greater activation of the ipsilesional sensorimotor system following BCI intervention. Moreover, we demonstrate that the increased Mu desynchronization in the ipsilesional primary motor cortex, Post vs. Pre BCI intervention, correlates significantly with improvements in hand function as assessed by grip force measurements. Furthermore, analyses of frequency-specific directional flow of information transmission between cerebral cortical sensorimotor areas, deduced from measurements of isolated effective coherence (iCoh), revealed both intra- and interhemispheric changes as a result of BCI intervention, in particular, increased functional connectivity from ipsilesional primary motor to ipsilesional premotor cortices, and from contralesional visuo-motor coordination area to ipsilesional premotor areas, respectively.
A growing body of evidence supports the effectiveness of EEG-based BCIs on improvement of upper extremity motor function following stroke (for meta-analyses/reviews see Cervera et al., 2018; Bai et al., 2020; Simon et al., 2021). BCI paradigms utilizing FES and/or attempted voluntary movements of the hemiparetic hand are particularly promising interventions for the rehabilitation of stroke survivors (Ramos-Murguialday et al., 2013; Kim et al., 2016; Biasiucci et al., 2018; Nishimoto et al., 2018; Remsik et al., 2018; Tabernig et al., 2018) because they may induce and/or facilitate neuroplastic changes, at both structural and functional levels, that link movement intention with muscle contraction (Daly and Wolpaw, 2008; Ramos-Murguialday et al., 2013; Biasiucci et al., 2018; Bai et al., 2020).
Source Localization of Mu ERD in Sensorimotor Cortices
In a less tightly controlled analysis of available NCT02098265 participants at the time, our laboratory demonstrated (Remsik A. et al., 2019) that ipsilesional Mu [8–12 Hz] ERD increases, recorded over cerebral cortical sensorimotor areas by scalp electrodes C3/C4 (Figure 1), were correlated with motor recovery of the affected upper extremity in stroke survivors. Here, we confirm and expand upon our previous findings in three ways. First, we present a more detailed source mapping to cortical space of Mu rhythm desynchronization and show that motor recovery is predominantly associated with Mu ERD changes ipsilesionally, near the primary motor cortex (Brodmann Area (BA) 4) and somatosensory association area (BA 5) (Figure 3). Second, we show that greater Mu desynchronization in ipsilesional primary motor cortex (BA 4) is significantly correlated with improvement in hand grip strength, an objective, quantitative measure of UE function and capacity, following BCI-FES intervention (Figures 4, 5). Third, we show evidence of functional connectivity changes between ipsilesional sensorimotor cortices following BCI-FES intervention, which is consistent with the view that BCI-FES use supports and/or engages fundamental neural mechanics that link volitionally controlled EEG signal changes with changes in objective behavioral measures of motor function and capacity. Such mechanisms may be specific to motor recovery (Biasiucci et al., 2018; Bai et al., 2020) and/or may be signatures of motor learning in general (Yuan et al., 2008; Edelman et al., 2019).
Our results are in line with findings of other studies that have employed similar BCI paradigms in cohorts of stroke survivors with similar demographics (Daly et al., 2009; Ramos-Murguialday et al., 2013; Ang et al., 2014; Ono et al., 2014; Popovic, 2014; Hu et al., 2015; Kim et al., 2016; Biasiucci et al., 2018; Nishimoto et al., 2018; Tabernig et al., 2018; Bai et al., 2020; Simon et al., 2021). For example, Biasiucci et al. (2018) report that only their BCI-FES group exhibited significant functional improvement following BCI intervention. Our result that increased ipsilesional Mu desynchronization is significantly correlated with increased grip strength is in good agreement with the finding of Biasiucci et al. (2018) that only their BCI-FES group exhibited a significant and lasting increase of the strength of the muscle targeted by the FES (in their study: extensor digitorum communis, which elicits full extension of the wrist and fingers). Furthermore, our result that SMR desynchronization is largely ipsilesional, indicating greater activation of the ipsilesional motor system following BCI-FES intervention, is in keeping with previous reports in stroke survivors (Mihara et al., 2013; Li et al., 2014; Biasiucci et al., 2018; Bai et al., 2020). BCI-FES devices may interface with native underlying mechanisms of neuromotor plasticity and control that have been demonstrated to engage and facilitate brain changes and behavior changes indicative of clinical motor recovery.
Functional Connectivity and Motor Recovery
Importantly, although voxel-wise distribution of eLORETA output revealed a sensorimotor-centered focus of Mu ERD change associated with motor recovery, ipsilesional-focused Mu ERD change may not be the sole driver of motor recovery in stroke survivors using a BCI. Mu signal changes in cortical sensorimotor areas may be part of a broader integration of disparate ROIs in a neural network that includes both intra- and interhemispheric components (Figures 8, 9), although recovery of the primary motor region has been demonstrated to be essential to motor recovery (Grefkes and Fink, 2014). Therefore, we sought to understand whether the underlying neural network is specific to localized neural populations represented by the Brodmann’s area segmentations or, whether there is a larger geography of populations – functionally connected neural networks – that may be influenced by participation in the BCI intervention. We used discrete data in time series to represent directionality or sequencing of activity associated with these changes to test the hypothesis that BCI intervention results in changes in frequency-specific directional flow of information transmission (functional connectivity) in ipsilesional and/or contralesional sensorimotor areas (Ward et al., 2003; Rehme et al., 2011; Dubovik et al., 2012; Várkuti et al., 2013; Nicolo et al., 2015; Pichiorri et al., 2015; Pundik et al., 2015; Wu et al., 2015; Biasiucci et al., 2018).
The changes in functional connectivity reported here, particularly between ipsilesional primary motor cortex and ipsilesional premotor cortices, agree with previous studies that have employed comparable BCI designs (Wu et al., 2015; Biasiucci et al., 2018). Increased interactions between ipsilesional motor brain areas are associated with better recovery and motor performance (Dubovik et al., 2012; Nicolo et al., 2015; Wu et al., 2015). The change in effective connectivity between ipsilesional motor cortex and premotor cortices following BCI intervention was observed for both Mu and Beta frequency ranges in both the present results (Figures 8, 9) and those of Biasiucci et al. (2018), which is consistent with the view that the functional improvement following BCI intervention is due to facilitation and induction of neuroplastic changes associated with motor recovery as well as motor learning more broadly.
Brain-Computer Interface-Functional Electrical Stimulation and Motor Recovery
BCI-mediated FES leads to the most significant recovery of motor function following stroke (Bai et al., 2020; Simon et al., 2021). The inclusion of FES is thought to contribute to the clinical effect of a BCI not only through the somatosensory contribution of facilitated muscle stimulation but also through pairing of volitionally modified CNS efferent signals with stimulation (i.e., activation) of the impaired distal muscle (Bergquist et al., 2011). Biasiucci et al. (2018) have demonstrated through their BCI-FES versus sham FES experimental design that it is not the FES alone but rather the contingency between rich sensory feedback and suitable activation of cortical motor areas that may drive activity-dependent, Hebbian plasticity that may underlie motor recovery in BCI-FES interventions. Whereas our study design does not allow us to draw the same conclusions as Biasiucci et al. (2018) with respect to the clinical significance of the FES adjuvant, both our results and BCI-FES device are similar to those of Biasiucci et al. (2018). Therefore, it is likely that the clinically relevant functional gains demonstrated here are due to the same strict contingency of BCI-driven FES detailed by Biasiucci et al. (2018).
Future research is needed to better identify and track the genesis and progression of neuroplastic changes, and to determine the relative importance of the various intra- and interhemispheric network connectivity changes presented here. In addition, further research is necessary to discover the mechanistic origins of any such neuroplasticity, and how such mechanism(s) may improve rehabilitation strategies that enable caregivers to provide maximal benefit to patients (Bai et al., 2020; Simon et al., 2021).
Limitations
Localizing cortical neuronal signatures with EEG scalp electrodes is a process that neuroscience researchers and others have used for some time (Pascual-Marqui et al., 1994; Pascual-Marqui, 2007). It is generally believed that higher density electrode arrays in the EEG cap provide more reliable and accurate localization estimations of the underlying neurophysiological processes. The use of eLORETA in a 16-channel cap arrangement, in non-research settings, likely pushes the boundaries of this method’s practical limits as an estimation tool for source localization of motor-EEG brain signals. Co-registration of spectral EEG data to MNI atlas space is a first step toward addressing this issue. In this work, scalp electrodes, while limited in number (16 channels), specifically covered sensorimotor cortices with greater density than the conventional 10–20 electrode placement system. Arguments can also be made for the use of “individualized” head modeling over “standard” head models. The use and reference to Brodmann’s areas in this work are intended solely as representative labels of the brain segmentations proposed by Brodmann, and we assume that these areas are generally emblematic of the functional cortical brain areas recorded by the scalp EEG electrodes of this cap arrangement (Figure 1).
In future research, repeating or initiating a similar BCI study with high density or 19 channel standard electrode arrays might allow a more strict evaluation of the “boundaries” in channel count and “electrode density as spatial coverage over areas of interest” by comparing the source localization results from these findings to the aforementioned high- and low-density electrode setups (i.e., with “full” head coverage). With any channel count placed “strategically” with higher density than the basic n = 19 (i.e.,10-20 system) over local functional network areas of interest on the scalp, reliable recordings of neuro-electrophysiological functioning might be feasible.
Conclusion
The results of the present study are consistent with the view that EEG-based BCI-FES intervention enhances information flow between cerebral cortical sensorimotor areas involved in motor planning, motor execution, and motor learning, and as such aides in establishing BCI-FES intervention as an effective therapy for motor rehabilitation of stroke survivors.
BCI-FES intervention may help facilitate adaptive brain changes, such as increased ERD during movement of the impaired upper extremity in the stroke-lesioned brain, but it is the brain’s ability to adapt its functional connections (i.e., plasticity of the sensorimotor system) that may underlie the potential of BCI-FES intervention as a rehabilitation strategy.
In conclusion, this study not only helps to establish the efficacy of BCI-FES intervention as a therapy for stroke survivors but is also important for increasing our understanding how sensorimotor processing contributes to the transformations of plans for voluntary limb movements into muscle commands necessary for their execution in healthy individuals. Thus, the results have the potential to guide development of innovative strategies for motor rehabilitation and are also important for increasing understanding of motor control in general.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by University of Wisconsin Health Sciences Institutional Review Board. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
AR involved in subject recruitment, staff training, data collection, data processing, data analysis, interpreting data, and writing of the manuscript. KG involved in analysis, interpreting writing, and editing of the manuscript. LW involved in staff training, data collection, data processing, data analysis, interpreting data, and writing of the manuscript. PK, SG, and EB contributed to writing of the manuscript, manuscript editing, and intellectual content. CR and SR involved in data collection and editing the manuscript. BY involved in subject recruitment, staff training, data collection, data processing, data analysis, interpreting data, intellectual content, and manuscript editing. VN contributed to subject recruitment, staff training, data collection, data processing, data analysis, interpreting data, manuscript editing, and intellectual content. KC, JW, and VP are co-PIs and were involved in study conception, design, staff training, subject recruitment, manuscript editing, intellectual content, and supervised all aspects of the study. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by NIH/NINDS grants 1R01NS105 646-01A1, RC1MH090912-01, T32GM008692, UL1TR000427, K23NS086852, T32EB011434, R01EB000856-06, and R01 EB009103-01 and by the DARPA RCI Program (MTO) N660 01-12-C-4025 and HIST Program (MTO) N66001-11-1-4013. Additional funding was also provided through a Coulter Translational Research Award, the American Heart Association Grant 1T32EB011434-01A1, AHA Innovative Research Award – National (Marcus Foundation) 15IRG22760009, AHA Midwest Grant in Aid Award 15GRNT25780033, the Foundation of ASNR, UW Milwaukee-Madison Intercampus Grants, the UW Graduate School, and by Shapiro Foundation Grants.
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank all the study participants and their families and caregivers for participating in the study as well as all the hospital and support staff for their work. The authors also wish to extend thanks to Justin Sattin for referrals of stroke subjects as well as clinical documentation.
References
Ang, K. K., Guan, C., Phua, K. S., Wang, C., Zhou, L., Tang, K. Y., et al. (2014). Brain-computer interface-based robotic end effector system for wrist and hand rehabilitation: results of a three-armed randomized controlled trial for chronic stroke. Front. Neuroeng. 7:30. doi: 10.3389/fneng.2014.00030
Babaiasl, M., Mahdioun, S. H., Jaryani, P., and Yazdani, M. (2016). A review of technological and clinical aspects of robot-aided rehabilitation of upper-extremity after stroke. Disabil. Rehabil. Assist. Technol. 11, 263–280. doi: 10.3109/17483107.2014.1002539
Bai, Z., Fong, K. N. K., Zhang, J. J., Chan, J., and Ting, K. H. (2020). Immediate and long-term effects of BCI-based rehabilitation of the upper extremity after stroke: a systematic review and meta-analysis. J. Neuroeng. Rehabil. 17:57. doi: 10.1186/s12984-020-00686-2
Benjamin, E. J., Blaha, M. J., Chiuve, S. E., Cushman, M., Das, S. R., Deo, R., et al. (2017). Heart disease and stroke statistics-2017 update: a report From the American Heart Association. Circulation 135, e146–e603. doi: 10.1161/cir.0000000000000485
Benjamin, E. J., Muntner, P., Alonso, A., Bittencourt, M. S., Callaway, C. W., Carson, A. P., et al. (2019). Heart disease and stroke statistics-2019 update: a report from the american heart Association. Circulation 139, e56–e528. doi: 10.1161/cir.0000000000000659
Bergquist, A. J., Clair, J. M., Lagerquist, O., Mang, C. S., Okuma, Y., and Collins, D. F. (2011). Neuromuscular electrical stimulation: implications of the electrically evoked sensory volley. Eur. J. Appl. Physiol. 111, 2409–2426. doi: 10.1007/s00421-011-2087-9
Biasiucci, A., Leeb, R., Iturrate, I., Perdikis, S., Al-Khodairy, A., Corbet, T., et al. (2018). Brain-actuated functional electrical stimulation elicits lasting arm motor recovery after stroke. Nat. Commun. 9:2421. doi: 10.1038/s41467-018-04673-z
Birbaumer, N. (2006). Breaking the silence: brain-computer interfaces (BCI) for communication and motor control. Psychophysiology 43, 517–532. doi: 10.1111/j.1469-8986.2006.00456.x
Birbaumer, N., Weber, C., Neuper, C., Buch, E., Haapen, K., and Cohen, L. (2006). Physiological regulation of thinking: brain-computer interface (BCI) research. Prog. Brain Res. 159, 369–391. doi: 10.1016/S0079-6123(06)59024-7
Boissy, P., Bourbonnais, D., Carlotti, M. M., Gravel, D., and Arsenault, B. A. (1999). Maximal grip force in chronic stroke subjects and its relationship to global upper extremity function. Clin. Rehabil. 13, 354–362. doi: 10.1191/026921599676433080
Brodmann, K. (1909). Vergleichende Localisationslehre der Glosshirnrinde in ihren Prinzipien dargestellt auf Grund des Zeelenbaues. Leipzig: Barth, Johann Ambrosius.
Bundy, D. T., Wronkiewicz, M., Sharma, M., Moran, D. W., Corbetta, M., and Leuthardt, E. C. (2012). Using ipsilateral motor signals in the unaffected cerebral hemisphere as a signal platform for brain-computer interfaces in hemiplegic stroke survivors. J. Neural Eng. 9:036011. doi: 10.1088/1741-2560/9/3/036011
Cervera, M. A., Soekadar, S. R., Ushiba, J., Millán, J. D. R., Liu, M., Birbaumer, N., et al. (2018). Brain-computer interfaces for post-stroke motor rehabilitation: a meta-analysis. Ann. Clin. Transl. Neurol. 5, 651–663. doi: 10.1002/acn3.544
Daly, J. J., Cheng, R., Rogers, J., Litinas, K., Hrovat, K., and Dohring, M. (2009). Feasibility of a new application of noninvasive Brain Computer Interface (BCI): a case study of training for recovery of volitional motor control after stroke. J. Neurol Phys. Ther. 33, 203–211. doi: 10.1097/NPT.0b013e3181c1fc0b
Daly, J. J., and Wolpaw, J. R. (2008). Brain-computer interfaces in neurological rehabilitation. Lancet Neurol. 7, 1032–1043. doi: 10.1016/S1474-4422(08)70223-0
Dubovik, S., Pignat, J. M., Ptak, R., Aboulafia, T., Allet, L., Gillabert, N., et al. (2012). The behavioral significance of coherent resting-state oscillations after stroke. Neuroimage 61, 249–257. doi: 10.1016/j.neuroimage.2012.03.024
Edelman, B. J., Meng, J., Suma, D., Zurn, C., Nagarajan, E., Baxter, B. S., et al. (2019). Noninvasive neuroimaging enhances continuous neural tracking for robotic device control. Sci. Robot. 4:eaaw6844. doi: 10.1126/scirobotics.aaw6844
Grefkes, C., and Fink, G. R. (2014). Connectivity-based approaches in stroke and recovery of function. Lancet Neurol. 13, 206–216. doi: 10.1016/s1474-4422(13)70264-3
Hu, X. L., Tong, R. K., Ho, N. S., Xue, J. J., Rong, W., and Li, L. S. (2015). Wrist rehabilitation assisted by an electromyography-driven neuromuscular electrical stimulation robot after stroke. Neurorehabil. Neural Repair 29, 767–776. doi: 10.1177/1545968314565510
Irimia, D., Sabathiel, N., Ortner, R., Poboroniuc, M., Coon, W., Allison, B. Z., et al. (2016). recoveriX: a new BCI-based technology for persons with stroke. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2016, 1504–1507. doi: 10.1109/embc.2016.7590995
Kim, T., Kim, S., and Lee, B. (2016). Effects of action observational training plus brain-computer interface-based functional electrical stimulation on paretic arm motor recovery in patient with stroke: a randomized controlled trial. Occup. Ther. Int. 23, 39–47. doi: 10.1002/oti.1403
Kitaura, Y., Nishida, K., Yoshimura, M., Mii, H., Katsura, K., Ueda, S., et al. (2017). Functional localization and effective connectivity of cortical theta and alpha oscillatory activity during an attention task. Clin. Neurophysiol. Pract. 2, 193–200. doi: 10.1016/j.cnp.2017.09.002
Leuthardt, E. C., Schalk, G., Wolpaw, J. R., Ojemann, J. G., and Moran, D. W. (2004). A brain-computer interface using electrocorticographic signals in humans. J. Neural Eng. 1, 63–71. doi: 10.1088/1741-2560/1/2/001
Li, M., Liu, Y., Wu, Y., Liu, S., Jia, J., and Zhang, L. (2014). Neurophysiological substrates of stroke patients with motor imagery-based brain-computer interface training. Int. J. Neurosci. 124, 403–415. doi: 10.3109/00207454.2013.850082
Mazrooyisebdani, M., Nair, V. A., Loh, P. L., Remsik, A. B., Young, B. M., Moreno, B. S., et al. (2018). Evaluation of changes in the motor network following BCI therapy based on graph theory analysis. Front. Neurosci. 12:861. doi: 10.3389/fnins.2018.00861
McFarland, D. J., Miner, L. A., Vaughan, T. M., and Wolpaw, J. R. (2000). Mu and beta rhythm topographies during motor imagery and actual movements. Brain Topogr. 12, 177–186. doi: 10.1023/a:1023437823106
Mihara, M., Hattori, N., Hatakenaka, M., Yagura, H., Kawano, T., Hino, T., et al. (2013). Near-infrared spectroscopy-mediated neurofeedback enhances efficacy of motor imagery-based training in poststroke victims: a pilot study. Stroke 44, 1091–1098. doi: 10.1161/STROKEAHA.111.674507
Mitra, P. P., and Pesaran, B. (1999). Analysis of dynamic brain imaging data. Biophys. J. 76, 691–708. doi: 10.1016/S0006-3495(99)77236-X
Mohanty, R., Sinha, A. M., Remsik, A. B., Dodd, K. C., Young, B. M., Jacobson, T., et al. (2018a). Early findings on functional connectivity correlates of behavioral outcomes of brain-computer interface stroke rehabilitation using machine learning. Front. Neurosci. 12:624. doi: 10.3389/fnins.2018.00624
Mohanty, R., Sinha, A. M., Remsik, A. B., Dodd, K. C., Young, B. M., Jacobson, T., et al. (2018b). Machine learning classification to identify the stage of brain-computer interface therapy for stroke rehabilitation using functional connectivity. Front. Neurosci. 12:353. doi: 10.3389/fnins.2018.00353
Nair, V. A., Young, B. M., La, C., Reiter, P., Nadkarni, T. N., Song, J., et al. (2015). Functional connectivity changes in the language network during stroke recovery. Ann. Clin. Transl. Neurol. 2, 185–195. doi: 10.1002/acn3.165
Nam, C. S., Jeon, Y., Kim, Y. J., Lee, I., and Park, K. (2011). Movement imagery-related lateralization of event-related (de)synchronization (ERD/ERS): motor-imagery duration effects. Clin. Neurophysiol. 122, 567–577. doi: 10.1016/j.clinph.2010.08.002
Neuper, C., and Pfurtscheller, G. (2001). Event-related dynamics of cortical rhythms: frequency-specific features and functional correlates. Int. J. Psychophysiol. 43, 41–58. doi: 10.1016/s0167-8760(01)00178-7
Neuper, C., Wörtz, M., and Pfurtscheller, G. (2006). ERD/ERS patterns reflecting sensorimotor activation and deactivation. Prog. Brain Res. 159, 211–222. doi: 10.1016/s0079-6123(06)59014-4
Nicolo, P., Rizk, S., Magnin, C., Pietro, M. D., Schnider, A., and Guggisberg, A. G. (2015). Coherent neural oscillations predict future motor and language improvement after stroke. Brain 138(Pt 10), 3048–3060. doi: 10.1093/brain/awv200
Nishimoto, A., Kawakami, M., Fujiwara, T., Hiramoto, M., Honaga, K., Abe, K., et al. (2018). Feasibility of task-specific brain-machine interface training for upper-extremity paralysis in patients with chronic hemiparetic stroke. J. Rehabil. Med. 50, 52–58. doi: 10.2340/16501977-2275
Nudo, R. J. (2011). Neural bases of recovery after brain injury. J. Commun. Disord. 44, 515–520. doi: 10.1016/j.jcomdis.2011.04.004
Nudo, R. J., and McNeal, D. (2013). Plasticity of cerebral functions. Handb. Clin. Neurol. 110, 13–21. doi: 10.1016/b978-0-444-52901-5.00002-2
Ono, T., Shindo, K., Kawashima, K., Ota, N., Ito, M., Ota, T., et al. (2014). Brain-computer interface with somatosensory feedback improves functional recovery from severe hemiplegia due to chronic stroke. Front. Neuroeng. 7:19. doi: 10.3389/fneng.2014.00019
Ortner, R., Irimia, D. C., Scharinger, J., and Guger, C. (2012). A motor imagery based brain-computer interface for stroke rehabilitation. Stud. Health Technol. Inform. 181, 319–323.
Pascual-Marqui, R. D. (2007). Discrete, 3D distributed, linear imaging methods of electric neuronal activity. Part 1: exact, zero error localization. arXiv [Preprint]. Available online at: https://arxiv.org/abs/0710.3341 (accessed October 7, 2021).
Pascual-Marqui, R. D., Biscay, R. J., Bosch-Bayard, J., Lehmann, D., Kochi, K., Kinoshita, T., et al. (2014). Assessing direct paths of intracortical causal information flow of oscillatory activity with the isolated effective coherence (iCoh). Front. Hum. Neurosci. 8:448. doi: 10.3389/fnhum.2014.00448
Pascual-Marqui, R. D., Michel, C. M., and Lehmann, D. (1994). Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. Int. J. Psychophysiol. 18, 49–65. doi: 10.1016/0167-8760(84)90014-x
Pfurtscheller, G. (1999). “EEG event-related desynchronization (ERD) and event-related synchronization (ERS),” in Electroencephalography: Basic Principles, Clinical Applications and Related Fields, 4th Edn, eds E. Niedermeyer and F. H. Lopes da Silva (Baltimore: Williams and Wilkins), 958–967.
Pfurtscheller, G., and Berghold, A. (1989). Patterns of cortical activation during planning of voluntary movement. Electroencephalogr. Clin. Neurophysiol. 72, 250–258. doi: 10.1016/0013-4694(89)90250-2
Pfurtscheller, G., and Lopes da Silva, F. H. (1999). Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin. Neurophysiol. 110, 1842–1857.
Pfurtscheller, G., Neuper, C., Andrew, C., and Edlinger, G. (1997). Foot and hand area mu rhythms. Int. J. Psychophysiol. 26, 121–135. doi: 10.1016/s0167-8760(97)00760-5
Pfurtscheller, G., Neuper, C., and Birbaumer, N. (2005). “Human Brain-Computer Interface,” in Motor Cortex in Voluntary Movements: A Distributed System for Distributed Functions, eds A. Riehle and E. Vaadia (Washington, DC: CRC Press), 367–401.
Pichiorri, F., Morone, G., Petti, M., Toppi, J., Pisotta, I., Molinari, M., et al. (2015). Brain-computer interface boosts motor imagery practice during stroke recovery. Ann. Neurol. 77, 851–865. doi: 10.1002/ana.24390
Pineda, J. A. (2005). The functional significance of mu rhythms: translating “seeing” and “hearing” into “doing”. Brain Res. Brain Res. Rev. 50, 57–68. doi: 10.1016/j.brainresrev.2005.04.005
Popovic, D. B. (2014). Advances in functional electrical stimulation (FES). J. Electromyogr. Kinesiol. 24, 795–802. doi: 10.1016/j.jelekin.2014.09.008
Popovic, D. B., Sinkaer, T., and Popovic, M. B. (2009). Electrical stimulation as a means for achieving recovery of function in stroke patients. NeuroRehabilitation 25, 45–58. doi: 10.3233/NRE-2009-0498
Prasad, G., Herman, P., Coyle, D., McDonough, S., and Crosbie, J. (2010). Applying a brain-computer interface to support motor imagery practice in people with stroke for upper limb recovery: a feasibility study. J. Neuroeng. Rehabil. 7:60. doi: 10.1186/1743-0003-7-60
Pundik, S., McCabe, J. P., Hrovat, K., Fredrickson, A. E., Tatsuoka, C., Feng, I. J., et al. (2015). Recovery of post stroke proximal arm function, driven by complex neuroplastic bilateral brain activation patterns and predicted by baseline motor dysfunction severity. Front. Hum. Neurosci. 9:394. doi: 10.3389/fnhum.2015.00394
Ramos-Murguialday, A., Broetz, D., Rea, M., Laer, L., Yilmaz, O., Brasil, F. L., et al. (2013). Brain-machine interface in chronic stroke rehabilitation: a controlled study. Ann. Neurol. 74, 100–108. doi: 10.1002/ana.23879
Rehme, A. K., Eickhoff, S. B., Wang, L. E., Fink, G. R., and Grefkes, C. (2011). Dynamic causal modeling of cortical activity from the acute to the chronic stage after stroke. Neuroimage 55, 1147–1158. doi: 10.1016/j.neuroimage.2011.01.014
Remsik, A., Bjorklund, E., Williams, L., Gjini, K., Mazrooyisebdani, M., Nair, V., et al. (2020). “Abstract TP158: direct path functional connectivity changes in stroke survivors’ brain motor areas using a brain computer interface intervention for upper extremity recovery,” in Proceedings of the International Stroke Conference, New Orleans, LA.
Remsik, A., Dodd, K., Gjini, K., Van Kan, P., Nair, V., Farrar Edwards, D., et al. (2019). “Abstract WP190: ipsilesional mu changes track with behavioral changes in UE BCI intervention in stroke survivors,” in Proceedings of the International Stroke Conference, New Orleans, LA.
Remsik, A., Young, B., Vermilyea, R., Kiekhoefer, L., Abrams, J., Evander Elmore, S., et al. (2016). A review of the progression and future implications of brain-computer interface therapies for restoration of distal upper extremity motor function after stroke. Expert Rev. Med. Devices 13, 445–454. doi: 10.1080/17434440.2016.1174572
Remsik, A. B., Bjorklund, E., Gjini, K., Williams, L., Van Kan, P. L. E., Nair, V., et al. (2020). “Abstract WP191: ipsilesional mu sensorimotor rhythms track with hand function recovery in stroke survivors,” in Proceedings of the International Stroke Conference, New Orleans, LA.
Remsik, A. B., Dodd, K., Williams, L. Jr., Thoma, J., Jacobson, T., Allen, J. D., et al. (2018). Behavioral outcomes following brain-computer interface intervention for upper extremity rehabilitation in stroke: a randomized controlled trial. Front. Neurosci. 12:752. doi: 10.3389/fnins.2018.00752
Remsik, A. B., Williams, L. Jr., Gjini, K., Dodd, K., Thoma, J., Jacobson, T., et al. (2019). Ipsilesional mu rhythm desynchronization and changes in motor behavior following post stroke BCI intervention for motor rehabilitation. Front. Neurosci. 13:53. doi: 10.3389/fnins.2019.00053
Schalk, G., McFarland, D. J., Hinterberger, T., Birbaumer, N., and Wolpaw, J. R. (2004). BCI2000: a general-purpose brain-computer interface (BCI) system. IEEE Trans. Biomed. Eng. 51, 1034–1043. doi: 10.1109/TBME.2004.827072
Schalk, G., and Mellinger, J. (2010). A Practical Guide to Brain–Computer Interfacing with Bci2000: General-Purpose Software for Brain–Computer Interface Research, Data Acquisition, Stimulus Presentation, and Brain Monitoring. London: Springer.
Simon, C., Bolton, D. A. E., Kennedy, N. C., Soekadar, S. R., and Ruddy, K. L. (2021). Challenges and opportunities for the future of brain-computer interface in neurorehabilitation. Front. Neurosci. 15:699428. doi: 10.3389/fnins.2021.699428
Song, J., Nair, V. A., Young, B. M., Walton, L. M., Nigogosyan, Z., Remsik, A., et al. (2015). DTI measures track and predict motor function outcomes in stroke rehabilitation utilizing BCI technology. Front. Hum. Neurosci. 9:195. doi: 10.3389/fnhum.2015.00195
Song, J., Young, B. M., Nigogosyan, Z., Walton, L. M., Nair, V. A., Grogan, S. W., et al. (2014). Characterizing relationships of DTI, fMRI, and motor recovery in stroke rehabilitation utilizing brain-computer interface technology. Front. Neuroeng. 7:31. doi: 10.3389/fneng.2014.00031
Tabernig, C. B., Lopez, C. A., Carrere, L. C., Spaich, E. G., and Ballario, C. H. (2018). Neurorehabilitation therapy of patients with severe stroke based on functional electrical stimulation commanded by a brain computer interface. J. Rehabil. Assist. Technol. Eng. 5:2055668318789280. doi: 10.1177/2055668318789280
Takahashi, M., Takeda, K., Otaka, Y., Osu, R., Hanakawa, T., Gouko, M., et al. (2012). Event related desynchronization-modulated functional electrical stimulation system for stroke rehabilitation: a feasibility study. J. Neuroeng. Rehabil. 9:56. doi: 10.1186/1743-0003-9-56
Várkuti, B., Guan, C., Pan, Y., Phua, K. S., Ang, K. K., Kuah, C. W. K., et al. (2013). Resting state changes in functional connectivity correlate with movement recovery for BCI and robot-assisted upper-extremity training after stroke. Neurorehabil. Neural Repair 27, 53–62. doi: 10.1177/1545968312445910
Ward, N. S., Brown, M. M., Thompson, A. J., and Frackowiak, R. S. (2003). Neural correlates of motor recovery after stroke: a longitudinal fMRI study. Brain 126(Pt 11), 2476–2496.
Wilson, J. A., Schalk, G., Walton, L. M., and Williams, J. C. (2009). Using an EEG-based brain-computer interface for virtual cursor movement with BCI200. J. Vis. Exp. 29:1319. doi: 10.3791/1319
Wilson, J. A., Walton, L. M., Tyler, M., and Williams, J. (2012). Lingual electrotactile stimulation as an alternative sensory feedback pathway for brain-computer interface applications. J. Neural Eng. 9:045007. doi: 10.1088/1741-2560/9/4/045007
Wolpaw, J. R., Birbaumer, N., McFarland, D. J., Pfurtscheller, G., and Vaughan, T. M. (2002). Brain-computer interfaces for communication and control. Clin. Neurophysiol. 113, 767–791. doi: 10.1016/s1388-2457(02)00057-3
Wolpaw, J. R., McFarland, D. J., Neat, G. W., and Forneris, C. A. (1991). An EEG-based brain-computer interface for cursor control. Electroencephalogr. Clin. Neurophysiol. 78, 252–259. doi: 10.1016/0013-4694(91)90040-b
Wu, J., Quinlan, E. B., Dodakian, L., McKenzie, A., Kathuria, N., Zhou, R. J., et al. (2015). Connectivity measures are robust biomarkers of cortical function and plasticity after stroke. Brain 138(Pt 8), 2359–2369. doi: 10.1093/brain/awv156
Young, B. M., Nigogosyan, Z., Nair, V. A., Walton, L. M., Song, J., Tyler, M. E., et al. (2014a). Case report: post-stroke interventional BCI rehabilitation in an individual with preexisting sensorineural disability. Front. Neuroeng. 7:18. doi: 10.3389/fneng.2014.00018
Young, B. M., Nigogosyan, Z., Remsik, A., Walton, L. M., Song, J., Nair, V. A., et al. (2014b). Changes in functional connectivity correlate with behavioral gains in stroke patients after therapy using a brain-computer interface device. Front. Neuroeng. 7:25. doi: 10.3389/fneng.2014.00025
Young, B. M., Nigogosyan, Z., Walton, L. M., Song, J., Nair, V. A., Grogan, S. W., et al. (2014c). Changes in functional brain organization and behavioral correlations after rehabilitative therapy using a brain-computer interface. Front. Neuroeng. 7:26. doi: 10.3389/fneng.2014.00026
Young, B. M., Nigogosyan, Z., Walton, L. M., Remsik, A., Song, J., Nair, V. A., et al. (2015). Dose-response relationships using brain-computer interface technology impact stroke rehabilitation. Front. Hum. Neurosci. 9:361. doi: 10.3389/fnhum.2015.00361
Young, B. M., Stamm, J. M., Song, J., Remsik, A. B., Nair, V. A., Tyler, M. E., et al. (2016). Brain-computer interface training after stroke affects patterns of brain-behavior relationships in corticospinal motor fibers. Front. Hum. Neurosci. 10:457. doi: 10.3389/fnhum.2016.00457
Young, B. M., Williams, J., and Prabhakaran, V. (2014d). BCI-FES: could a new rehabilitation device hold fresh promise for stroke patients? Expert Rev. Med. Devices 11, 537–539. doi: 10.1586/17434440.2014.941811
Keywords: brain-computer interface, sensorimotor rhythm, Mu, Beta, stroke, motor function, iCoh
Citation: Remsik AB, Gjini K, Williams L Jr, van Kan PLE, Gloe S, Bjorklund E, Rivera CA, Romero S, Young BM, Nair VA, Caldera KE, Williams JC and Prabhakaran V (2021) Ipsilesional Mu Rhythm Desynchronization Correlates With Improvements in Affected Hand Grip Strength and Functional Connectivity in Sensorimotor Cortices Following BCI-FES Intervention for Upper Extremity in Stroke Survivors. Front. Hum. Neurosci. 15:725645. doi: 10.3389/fnhum.2021.725645
Received: 15 June 2021; Accepted: 01 October 2021;
Published: 28 October 2021.
Edited by:
Takufumi Yanagisawa, Osaka University, JapanReviewed by:
Bradley Jay Edelman, Max Planck Institute of Neurobiology (MPIN), GermanyJinung A. N., Daegu Gyeongbuk Institute of Science and Technology (DGIST), South Korea
Copyright © 2021 Remsik, Gjini, Williams, van Kan, Gloe, Bjorklund, Rivera, Romero, Young, Nair, Caldera, Williams and Prabhakaran. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) 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: Alexander B. Remsik, aremsik@wisc.edu