The role of virtual reality in improving motor performance as revealed by EEG: a randomized clinical trial
Journal of NeuroEngineering and Rehabilitation volume 14, Article number: 53 (2017)
Many studies have demonstrated the usefulness of repetitive task practice by using robotic-assisted gait training (RAGT) devices, including Lokomat, for the treatment of lower limb paresis. Virtual reality (VR) has proved to be a valuable tool to improve neurorehabilitation training. The aim of our pilot randomized clinical trial was to understand the neurophysiological basis of motor function recovery induced by the association between RAGT (by using Lokomat device) and VR (an animated avatar in a 2D VR) by studying electroencephalographic (EEG) oscillations.
Twenty-four patients suffering from a first unilateral ischemic stroke in the chronic phase were randomized into two groups. One group performed 40 sessions of Lokomat with VR (RAGT + VR), whereas the other group underwent Lokomat without VR (RAGT-VR). The outcomes (clinical, kinematic, and EEG) were measured before and after the robotic intervention.
As compared to the RAGT-VR group, all the patients of the RAGT + VR group improved in the Rivermead Mobility Index and Tinetti Performance Oriented Mobility Assessment. Moreover, they showed stronger event-related spectral perturbations in the high-γ and β bands and larger fronto-central cortical activations in the affected hemisphere.
The robotic-based rehabilitation combined with VR in patients with chronic hemiparesis induced an improvement in gait and balance. EEG data suggest that the use of VR may entrain several brain areas (probably encompassing the mirror neuron system) involved in motor planning and learning, thus leading to an enhanced motor performance.
Retrospectively registered in Clinical Trials on 21-11-2016, n.NCT02971371.
Virtual reality (VR) is the simulation of a real environment generated by a computer software and experienced by the user through a human–machine interface . This interface enables the patient to perceive the environment as real and 3D (i.e., the sense of presence), thus increasing patient’s engagement (i.e., embodiment) . Hence, VR can be used to provide the patient with repetitive, task-specific training (as opposed to simply using a limb by chance) that are effective for motor learning functions [3,4,5,6]. In fact, VR provides the patient with multisensory feedbacks that can potentiate the use-dependent plasticity processes within the sensory-motor cortex, thus promoting/enhancing functional motor recovery [7,8,9,10,11,12,13,14]. Furthermore, VR can increase patients’ motivation during rehabilitation by decreasing the perception of exertion , thus allowing patients to exercise more effortlessly and regularly .
It is possible to magnify the sense of presence by manipulating the characteristics of the VR, including screen size, duration of exposure, the realism of the presentation, and the use of animated avatar, i.e., a third-person view of the user that appears as a player in the VR . About that, the use of an avatar may strengthen the use-dependent plastic changes within higher sensory-motor areas belonging to the mirror neuron system (MNS) [16,17,18]. In fact, the observation of an action, even simulated (on a screen, as in the case of VR), allows the recruitment of stored motor programs that would promote, in turn, movement execution recovery [19, 20]. These processes are expressed by wide changes in α and β oscillation magnitude at the electroencephalography (EEG) (including an α activity decrease and a β activity increase) across the brain areas putatively belonging to the MNS (including the inferior frontal gyrus, the lower part of the precentral gyrus, the rostral part of the inferior parietal lobule, and the temporal, occipital and parietal visual areas) [8, 9, 21, 22].
In the last years, motor function recovery has benefited from the use of robotic devices. In particular, robot-assisted gait training (RAGT) provides the patient with highly repeated movement execution, whose feedback, in turn, permits to boost the abovementioned use-dependent plasticity processes . RAGT has been combined with VR to further improve gait in patients suffering from different neurologic diseases . Nonetheless, the knowledge of the neurophysiologic substrate underpinning neurorobotic and VR interaction is still poor [25, 26]. Indeed, a better understanding of this interaction would allow physician to design more personalized rehabilitative approaches concerning the individual brain plasticity potential to be harnessed to gain functional recovery .
The relative suppression of the μ rhythm is considered as the main index of MNS activity . Nonetheless, conjugating VR and neurorobotic could make brain dynamics more complex, because of many factors related to motor control and psychological aspects come into play, including intrinsic motivation, selective attention, goal setting, working memory, decision making, positive self-concept, and self-control. Altogether, these aspects may modify and extend the range of brain rhythms deriving from different cortical areas related to MNS activation by locomotion, including theta and gamma oscillations [29,30,31]. Specifically, theta activity has been related to the retrieval of stored motor memory traces, whereas the gamma may be linked to the conscious access to visual target representations [30, 31]. Such broadband involvement may be due to the recruitment of multiple brain pathways expressing both bottom-up (automatic recruitment of movement simulation) and top-down (task-driven) neural processes within the MNS implicated in locomotion recognition . A recent work has shown that observed, executed, and imagined action representations are decoded from putative mirror neuron areas, including Broca’s area and ventral premotor cortex, which have a complex interplay with the traditional MNS areas generating the μ rhythm .
Therefore, we hypothesized that the combined use of VR and RAGT may induce a stronger and wider modification of the brain oscillations deriving from the putative MNS areas, thus augmenting locomotor function gain [34, 35]. The aim of our pilot randomized clinical trial was to understand the neurophysiological basis underpinning gait recovery induced by the observation of an animated avatar in a 2D VR while performing RAGT by studying the temporal patterns of broadband cortical activations.
The present randomized clinical trial was conducted according to the CONSORT guidelines . The trial was designed as a pilot, prospective, assessor blinded, parallel group study, and was performed at the IRCCS Centro Neurolesi “Bonino-Pulejo” (Messina, Italy). Eligible patients were enrolled between October 2015 and February 2016, according to the following criteria: (i) age ≥ 55 years; (ii) a first-ever ischemic supra-tentorial stroke (confirmed by magnetic resonance imaging -MRI) at least 6 months before their enrollment; (iii) an unilateral hemiparesis, with a Muscle Research Council -MRC- score ≤ 3 (a score of 3 indicates that muscle strength is reduced so that the joint can be moved only against gravity with the examiner’s resistance completely removed; 2 = muscle can move only if the resistance of gravity is removed; 1 = only a trace of movement; 0 = no movement observed) ; (iv) ability to follow verbal instructions, with a Mini-Mental State Examination >24; (v) a mild to moderate spasticity of muscles of hip, knee, and ankle (according to a Modified Ashworth Scale, MAS, ≤2) ; (vi) ability to perform manual gait training with or without external devices (Functional Ambulatory Categories 0-4). (vii) no severe bone or joint disease; and (viii) no history of concomitant neurodegenerative diseases or brain surgery.
We preferred to select patients with a first-ever ischemic supra-tentorial stroke, as this represents a better model of stroke lesion to perform EEG analysis. In fact, multiple vascular lesions could represent a significant limit for data interpretation, as they can generate variable signals that can interfere with signal recording. For instance, it has been reported that hemispheric powers differ clearly in single acute ischemic episodes, but correlate less well with the subtle, multifocal, or more gradual changes . Moreover, single lesion model is more suitable to study interhemispheric balance. Last, patients with multiple vascular lesions may have different functional recovery. For the same reason, we limited the age of inclusion to >55, because beneath this age it is necessary to take into account other several, additional risk factors (including hemostatic, inflammatory and autoimmune factors, and cardioembolic sources, such as patent foramen ovale) that altogether increase the risk of multiple vascular lesions and may account for heterogeneity of sample .
Clinic-demographic characteristics are reported in Table 1. All participants gave informed consent before study participation. Additionally, written informed consent for publication of clinical images was obtained from the participants. Approval was obtained from our local Ethics Committee before beginning the study (study number registration 43/2013). The study was retrospectively registered in Clinical Trials on 21-11-2016, n.NCT02971371.
The sample size estimate was based on extrapolations from previous studies examining the effects of VR on gait in patients with stroke [41,42,43,44]. Accordingly, we used the effect size (0.9) of the primary composite endpoint for calculations. Power was set at 80%, alpha at 5%; we accounted for a dropout rate of 10%. Using a relatively conservative estimation, a total of 50 subjects (25 in each group) would be required to detect a difference of at least 20% in the primary outcome between the two treatment groups assuming non-inferiority with moderate correlations among covariates (R-squared = 0.5).
Thirty-five patients were assessed for eligibility according to the inclusion/exclusion criteria. Then, 24 patients were equally randomized into the RAGT + VR or RAGT-VR group with a 1:1 allocation ratio. For randomization, sealed envelopes were prepared in advance and marked on the inside with a + VR or -VR. The intervention period of both groups was identical, five sessions per week for eight consecutive weeks, 45-min for each session. We registered the EEG in the morning of the second day of the first week (TPRE) and the last day of the last week of treatment (TPOST). Clinic and kinematic assessment was performed immediately before and after the Lokomat training. In addition, a physiotherapist supervised patients’ cooperation and participation to Lokomat training. The experimenters who analyzed the data and the therapists who performed the clinical tests were blind on patient allocation. Moreover, patients were not informed on the content of the VR adopted in the two groups.
Prior to the Lokomat training, all participants became familiarized with the Lokomat-Pro (RAGT + VR) or Lokomat-Nanos (RAGT-VR) (Hocoma AG; Volketswil, Switzerland), and had individually adapted their body weight support (BWS), leg guidance force (GF), and foot-lifting straps (that assisted ankle dorsiflexion for adequate toe clearance during the swing phase). Gait speed was set at a maximum of 1.8 m/s. Patients were asked to walk with their maximal effort. Then, patients were randomly assigned to RAGT + VR or RAGT-VR group. Both groups performed 40 Lokomat sessions (each lasting 40-45 min), five times a week, between 9am and 11am. The RAGT + VR group received a visual feedback showing a VR run game where the patient had to collect or avoid objects, to motivate him/her to walk actively (Fig. 1). Each avatar’s leg movement corresponded to that performed by the patient. To this end, we employed an impedance-based control strategy [45, 46], which allows variable deviations from a given leg trajectory, thus making the orthosis feel compliant (instead of making the orthosis stiff, i.e., allowing only negligible deviations from the set trajectory). The Lokomat is equipped with potentiometers and force transducers, capable of providing feedback about joint movement and joint moment production, respectively. In the impedance-based control strategy, the joint moment acting on the subject’s leg (which is generated by the force transducers of the orthoses) is proportional to the angular deviation of the performed movement from the preprogrammed movement. Therefore, the patient perceives the stiffness of the orthosis based on the setting of the mechanical impedance value (i.e., zero impedance makes no feedback experienceable, whereas maximal impedance results in maximal stiffness, equal to the position controller). The stiffness of the orthoses was set to the lowest point where subjects reported feeling the inertia of the Lokomat segments (about 30% of the maximal Lokomat impedance).
On the other hand, the RAGT-VR group was not provided an avatar but a smile, indicating the goodness of the leg movements. The biofeedback of the Lokomat gait orthosis is based on the interaction torques between the participant and the orthosis. To this end, hip and knee linear drives are equipped with force sensors that measure the human-machine interaction forces at the hip and knee joints for stance and swing phase, which are required to keep the participant on a predefined gait trajectory. The weighting functions of the gait cycle were defined for each part so that the resulting biofeedback values increased for therapeutically desirable movements, i.e. hip and knee flexion-extension force. Thus, biofeedback force values (representing the physical activity of the participants) were positive when the patient was actively participating and negative when passively participating in Lokomat training (or when muscles were inappropriately or involuntary activated) . We obtained an average biofeedback value for each joint in the hemiparetic leg.
The 2D VR set-up consisted of a 42-in. flat-screen placed in front of the Lokomat and a 7.1 Dolby Surround system (Fig. 1). The Lokomat device served as a multimodal feedback system: the human-machine interaction forces measured from the Lokomat device were used as an input device for the patient’s movements into the VR (i.e., to animate the motion of the human figure in VR at a 60 Hz refresh rate in real time on the screen – virtual mirror. There was no lag between motions of the subject and virtual figure) or the smiling face. The orthosis guided subject’s leg movements in the sagittal plane within individually adapted hip and knee joint trajectories. Lokomat potentiometers provided real-time information of the subject’s hip and knee angles. The measured joint angles were used to animate the subject’s human figure in the virtual mirror. The distance between the subject and projection screen was 1.5 m. Furthermore, since the run lane was not always straight, the VR running game used an asymmetrical physical activity of the legs to induce turning in the virtual environment. In particular, turning right and left was induced by increasing the activity of the contralateral leg of the desired direction and, respectively, decreasing the activity of the ipsilateral one. Lokomat device provided visual and acoustic feedback that reflected the interactions with objects represented in the virtual environment (e.g., the boundaries of the run lane, or objects to be avoided or collected). Further, Lokomat provided haptic feedback by the gait orthosis, so that the subjects using the device were provided with a haptic experience from proprioceptive (joint angles) feedback about their movements [47, 48].
The use of 2D displays, which are not as realistic as the true stereo 3D ones (full-3D VR), are akin to looking at a scene through a window and offer a limited sense of presence, potentially limiting the significance of our findings. We potentially provided a higher sense of presence by using depth cues, such as perspective, relative motion, occlusion, and aerial perspective, despite the use of a 2D VR .
EEG recording and preprocessing
EEG was recorded using a Brain-Quick System (Micromed; Mogliano Veneto, Italy), from standard 19 electrodes headset according to the International 10-20 system (Fp1, Fp2, F7, F3, Fz, F4, F8, T3, C3, Cz, C4, T4, T5, P3, Pz, P4, T6, O1, O2, ground on the forehead), for 10 min while performing Lokomat training. To monitor eye movements, an electro-oculogram (EOG) with a bipolar montage (one pair of electrodes traced horizontal eye movements, a second pair the vertical ones) was also collected. EEG end EOG were sampled at 500 Hz, high pass filtered at 1 Hz using a zerophase FIR filter (order 7500) to minimize drifts, low pass filtered at 200 Hz (zerophase FIR filter order 36), and referenced to Cz . An adaptive filter was applied to allow real-time filtering of signals recorded from EOG [50, 51].
Electrode impedance was kept below 5kΩ. During the entire EEG recording (as well as during the entire gait training), an experimenter checked for possible signs of drowsiness (e.g., abrupt worsening in gait performance, closed eyes, increase of proportion of theta and alpha activity in the eyes-open condition) , which were counted (given that monotonous gait pattern provided by RAGT may tend to induce sleepiness, thus decreasing arousal that negatively affects gait training progress). Patients were prohibited from drinking coffee, smoking, and change their bedtime during the three days prior EEG recording.
Infomax independent component analysis (ICA) was computed on the preprocessed EEG signal to decompose neural and artefactual sources [53,54,55,56,57]. In detail, ICA was computed two times. First, 500 ms-segmented EEG signals were removed if its probability distribution exceeded the average distribution by 5 ± SD. Then, ICA was computed to reject epochs based on the probability distribution of the IC projections. Thus, EEG segments were re-filtered (8-40 Hz) and a second ICA was computed a second time. The so-obtained IC were grouped into clusters using a k-means algorithm (based on the feature vector of dipole location, power spectra, and scalp map). The IC closest to the cluster centroid was remained for each subject, so to have equal contribution of each subject to the cluster-wise analysis.
EEG analysis consisted of the computation of the power spectral density (PSD) (using Welch’s Method) and the time-frequency analysis to evaluate Event-related-spectral-perturbations (ERSPs) for each IC .
EEG was segmented into 1.4 s epochs (−700;700)ms with regard to the heel strike (HS) (i.e., the first moment the foot comes into contact with the floor) , thus obtaining 428 epochs. About that, the force-sensing resistor of the Lokomat device detected the movement onset of both lower limbs, which was synchronized with the EEG data. Epochs were rejected by using an automatic artifact rejection method (epochs with values of [−100;100]μV, ≥5SD of the mean kurtosis value, ≥5SD of the mean probability distribution, drifts of ≥50 μV/epoch and with a R2 limit ≤0.3, spectra deviating from the mean by ±50 dB in the 0-2 Hz frequency window and by [−100;50]dB in the 20-100 Hz frequency window), visually inspection for artifacts, and if the power perturbation in the 20–40 Hz band deviated by +25 or -100 dB from the baseline at least for one IC [57, 60]. Rejection rate was 5%. This low rate is not surprising, given that it has been shown that scalp EEG recording during low-speed treadmill walking is not invalidated by excessive artefacts . After this, the segmented data were time-warped and averaged together for all strides, so that initial affected-side heel strike, unaffected-side toe off, unaffected-side heel strike, affected-side toe off, and the subsequent affected-side heel strike occurred at the same times . Spectrum analysis was carried using a standard fast Fourier transform (FFT) algorithm (Hanning-window, frequency resolution 0.7 Hz) within ϑ (4-7 Hz), μ (8–12 Hz), β (12–30 Hz), low-γ (Lγ) (31-45 Hz), and high-γ (Hγ) (46-70 Hz) bands [62, 63], and related to the phases of the gait cycle . We opted to analyze these rhythms as it has been reported a different, specific role of each oscillation in sensory-motor pattern [27,28,29,30,31,32,33,34,35, 64]. For instance, there is evidence for a difference between the low and high α oscillations, which express action execution and observation, respectively [65, 66].
Single trial spectograms were computed and time-warped (thus aligning the time-points for right and left heel strike) over trials using a linear interpolation function to generate gait cycle ERSPs (i.e., epochs were based on the heel strike events, being the unaffected-side, the affected-side, and the next unaffected-side heel strike time-warped to 0, 50%, and 100% of the gait cycle, respectively). Relative changes in spectral power were obtained by averaging the difference between each single-trial log spectogram and baseline (the mean IC log spectrum over all gait cycles per training) . To visualize significant ERSP changes, deviations from the average gait cycle log spectrum were computed with a bootstrap method [56, 57]. For statistical concern, bandwidth ERSP of each IC were averaged within each 10% of the gait cycle (10-point ERPS curve, frequency resolution 0.7 Hz) . The average log spectrum for all movement cycles was subtracted from the log spectrogram for each movement cycle. We thus calculated the resulting PSD changes from this baseline (defined as the percentage decrement, event-related desynchronization –ERD- and increment, event-related synchronization –ERS- as a function of the percentage of the normalized gait cycle)  for each band and electrode-group of interest (with regard to the areas of activation of MNS reported in the literature , i.e., ipsi and contralesional frontal -Fp1/F7/F3, Fp2/F8/F4-, central -T3/C3, T4/C4- and parieto-occipital -T5/P3/O1, T6/P4/O2) [69,70,71,72,73].
The source localization approach allows examining brain activities in various sources at different temporal phases of motor control. Because of high temporal resolution of the EEG signals, brain activities before and after the movement onset can be localized in order to distinguish cortical activities related to both motor planning (movement preparation) and motor execution (corticospinal pathway activation). The Estimation of Current Densities was carried by using Low Resolution Brain Electromagnetic Tomography (LORETA; free release of LORETA-KEY alpha-software) [63, 74,75,76]. The main components detected with the ICA (signal-to-noise ratio > 1) were chosen for the source reconstruction. The distributed current density model (LORETA) with L1 norm method (based on the Montréal Neurological Institute (MNI) brain MRI) was then applied to the ICA data [77, 78]. The sources were constrained to the reconstructed layer of the folded cortex [79, 80].
The primary endpoint, with respect to VR efficacy in post-stroke condition, was the proportion of patients achieving a 20% improvement in lower limb gait and balance at the end of the training, as measured by the Rivermead Mobility Index (RMI), the Tinetti Performance Oriented Mobility Assessment (POMA), and the gait cycle-related ERSPs. Indeed, a 20% improvement correspond to a significant minimal detectable change in RMI  and POMA . According to previous work on assisted gait training in post-stroke patients, these changes are paralleled by EEG signal modification of at least 20% to be significant [83, 84].
As secondary outcomes, we considered the global MAS score derived from the muscles of hip, knee, and ankle, the Hamilton Rating Scale for Depression (HRS), the hip and knee flexion/extension force measured by the RAGT device, the extent to which a patient felt him/herself entrained in the VR training (reported on a visual analogue scale -VAS- ranging from zero -not at all- to ten -very much), and the mean of the episodes of drowsiness.
The normal distribution of the data was evaluated with the Kolmogorov-Smirnov test. Baseline data were compared between the two groups using a Student t-test for continuous variables if data were normally distributed, whereas a Mann-Whitney U test was used for non-normally distributed ordinal scale. Likewise, Wilcoxon test, Mann-Whitney U test, or t-test were used for within-group and between-group comparisons, depending on the types of data measurements.
The ERD/ERS changes for each frequency band were assessed by means of three-way ANOVA for repeated measures, employing the factor time (two levels: TPRE and TPOST) and electrode-set (three levels: frontal, central, and parieto-occipital) as within-subject factors, and group (two levels: RAGT + VR and RAGT-VR) as between-subject factor. Based on the significance of the F-value, post-hoc paired-sample t-tests were carried out to assess the significance of interactions (Bonferroni correction). A p-value <0.05 was considered significant.
BWS and GF were included as covariates in the ANOVA analysis. In fact, it has been reported that RAGT training usually implies a steady progression of BWS and GF across the training program, so it is necessary to update and analyze these parameters in relation to the assessment of the outcomes . In fact, both BWS and GF can influence spatiotemporal movement characteristics, thus affecting functional gait pattern. Indeed, finely tuning BWS and GF may somehow improve possible spatiotemporal gait asymmetries. On the other hand, missing the correction of these parameters augments inter-limb gait asymmetry for an extended duration in people with stroke . Besides the factor electrode-set (which was employed in the ANOVA analysis to carry the spatio-temporal analysis of EEG signals at scalp level), we also added the factors lesion localization as covariate in the ANOVA analysis (according to the localization within left or right frontal, parietal, occipital, and temporal lobe). Such factor was added to augment inter-subject evaluation, as the sample was non-homogeneous for stroke localization, which can affect EEG signals beyond the overhead electrodes , also influencing both motor deficit degree and recovery .
ERSPs were computed in each frequency range for RAGT + VR and RAGT-VR using the average gait cycle log spectrum computed from the RAGT-VR as common baseline. The gait cycle was divided into in two stationary (S1, 10–30%, and S2, 60–80%) and two transition phases (T1, 30–60%, and T2, 80–10%). The stationary phases correspond to the midstance (10–30%), initial swing (60–70%), and midswing phases (70–90%), whereas the transition phases correspond to the terminal stance (30–50%), preswing (50–60%), terminal swing (90–100%), and loading response (0–10%) . An ANOVA for repeated measures with the factors time (in relation with the gait cycle phases) (eight levels: two PRE and POST stationary and two PRE and POST transition phases), electrode-set (three levels: frontal, central, and parieto-occipital), and group (two levels: RAGT + VR and RAGT-VR) was computed for each frequency band. Multiple comparisons were corrected controlling for false discovery rate (p < 0.05) . Sphericity assumption violations were Greenhouse-Geisser corrected.
We summarized in the CONSORT flow diagram (Fig. 2): the numbers of participants who were randomly assigned, who received the intended treatment, and who were analyzed for the primary outcome; the losses and exclusions during periods of recruitment, randomization, and follow-up. All treated patients completed the robotic training without reporting any side effect (Fig. 2). The main analysis focused on the consequences of the avatar observation on the RMI, POMA, and gait-related ERS/ERD, and involved all the patients who were randomly assigned to the two groups. Twelve patients were enrolled in each group (RAGT + VR: mean age 60 ± 4 years, 7 males and 5 females, disease duration 8 ± 2 months; RAGT-VR: mean age, 63 ± 6 years, 7 males and 5 females, disease duration 8 ± 2 months) (Table 1). There were no significant differences concerning any parameter between the two groups. Mean and standard deviation values of all outcome measures are reported in Table 1.
Clinical and kinematic data
RMI improved more in the RAGT + VR than the RAGT-VR group, whereas POMA improved only in the RAGT + VR group (Table 2). HRS score equally decreased in both groups. MAS did not change significantly in both groups. VAS improved only in RAGT + VR group (Table 2). Knee force showed a greater improvement in RAGT + VR than RAGT-VR group, whereas hip force improved only in RAGT + VR group (Table 2).
We found three main areas of brain activation that were more evident in the RAGT + VR group as compared to the RAGT-VR group across the gait cycle (Fig. 3a): (1) BA6 (Tailarach coordinates -x,y,z- 5, −1, 60); (2) BA7 (−14, −56, 53); and (3) BA17 (−20, −88, 3).
Concerning frontal activation (Fig. 3b; Table 3), RAGT + VR induced significant μ/β ERSP changes from baseline to S1, T1, and S2 phase of the gait cycle, as reflected by the significant interaction between gait phases and training for gait cycle related modulations concerning RAGT + VR group (Table 4). In particular, we observed a frontal ERD before the heel strike, followed by a frontal ERS during the heel strike (Fig. 4; Table 4). This changes also presented a significant difference in scalp projections between the trainings (Table 5; Fig. 5), with a greater magnitude in RAGT + VR group.
μ/β changes were paralleled by a significant Lγ-ERS around the heel strike in both the groups and an Hγ-ERD from S1 to S2 (i.e., ending soon after the heel strike) only in RAGT + VR. Detailed statistical data of the main interactions and effects, and the temporal ERSP changes across the gait cycle are summarized in Tables 4 and 5 and Figs. 4 and 5.
Concerning central areas activation (Fig. 3b; Table 3), we found a significant μ/β ERD/ERS at S1, T1, and S2, a brief ϑ-ERS in the S1 phase, as revealed by the ANOVA and the post-hoc tests (Tables 4 and 5; Figs. 4 and 5) showing that spectral power in μ/β band was significantly reduced in T1 and increased in S1, and this was more evident in RAGT + VR than RAGT-VR.
Finally, parieto-occipital activation (Fig. 3b; Table 3) was characterized by a significant difference for μ and Hγ (ERD) between the groups during the entire gait cycle (Tables 4 and 5; Figs. 4 and 5). Post-hoc tests showed that scalp projections were significantly strengthened in RAGT + VR as compared to RAGT-VR.
We observed a significant correlation between RMI and POMA score improvement and central β-ERS (r = 0.895, p = 0.001 and r = 0.570, p = 0.04, respectively) and frontal Hγ-ERD magnitude (r = −0.831, p = 0.003 and r = −0.615, p = 0.04, respectively). Moreover, the improvement in hip force was significantly correlated with frontal Hγ-ERD magnitude (r = −0.802, p = 0.004). Clinical-demographic characteristics (age, gender, disease duration, comorbidity) did not influence the primary outcomes.
The main finding of our pilot study consists in the more evident activation of premotor, precuneus, and associative visual areas in the RAGT + VR group as compared to RAGT-VR group.
All the patients belonging to RAGT + VR showed a significant decrease of central μ/β power during the phase preceding the heel strike, followed by a power increase (as shown by the gait cycle phase dependent ERSP modulation), thus indicating higher neuronal activation . Importantly, we observed that the stronger the μ/β ERSPs were, the higher the clinical amelioration. Given that these ERSPs are a marker of activation and deactivation/inhibition of sensorimotor areas concerning motor planning, postural stabilization, and the prediction of potential actions [92,93,94,95,96,97,98,99,100,101,102], our findings suggest the importance of enhancing μ/β ERSPs to foster locomotor training. In addition, to monitor these brain activations would allow a better patient-tailored walking training.
The novelty of our study is the significant fronto-parietooccipital Hγ-ERD and parietooccipital α-ERD only in the RAGT + VR group. The premotor-parietooccipital desynchronization of γ-oscillations is thought to be a marker of activation of sensorimotor and visuo-spatial associative areas concerning motor planning and selective muscle activation [91, 100, 102,103,104,105,106,107,108,109,110,111,112] even during active and passive RAGT .
We also found that the magnitude of γ-band modulation was significantly correlated with the clinical amelioration and the improvement in muscle strength, and it was paralleled by a more selective μ/β-band modulation concerning either the temporal patterns of activation across the gait cycle or the hemispheric distribution of ERSPs. We may argue that VR may induce a functional fronto-parietooccipital α/γ-band activation that, in turn, allows a more efficient motor planning and execution, as shown by a stronger and selective modulation of μ rhythm across the phases of gait cycle. Such a selective modulation allows the patient to complete better the gait training (e.g., to better steer, avoid objects, and keep the line during walking) . These data are in keeping with the role of the premotor areas in planning limb movements  and initiating and adapting gait [115,116,117,118], and of the parieto-occipital cortex in spatial attention, decision making, sensorimotor integration, and movement planning in visually guided movements under both feedforward and feedback control [119,120,121,122,123,124]. The specific entrainment of γ rhythms when observing a human avatar may depend on a different entrainment of visuomotor networks as compared to the control condition (RAGT-VR). According to the canonical microcircuit model , the superficial pyramidal neurons generating γ-responses act as a dynamic filter on the visual inputs, thus affecting both the configuration of γ-oscillations (depending on the stimulus properties, including movement, contrast, localization, and size of visual cues) and the μ and β band output dynamics (which are generated by deeper pyramidal neurons) [126,127,128]. The use of an avatar may have thus specifically increased the frontal-posterior γ oscillations. The parieto-occipital α-ERD may be instead linked to basic visual processing. In addition, it has been reported that μ/γ ERSP provided by VR feedback is related to the participants’ monitoring of their own movements [129,130,131]. We can therefore hypothesize that these ERSPs in fronto-parietooccipital regions during the observation of performed movements and during visually-guided gait adaptation task potentially express the activation of the MNS.
One could argue that Hγ ERSPs may purely reflect motor activation and not specific cognitive processes related to VR, given that γ-band ERSPs express also a higher cortico-muscular connection during ambulation [132,133,134]. Nonetheless, this concern sounds unlikely, since BWS and DGF, which both change muscle activity [135, 136], were individually adapted in all patients. Consequently, Lγ band (which is instead strongly related to motor activity level) [132, 133, 137] was similar in the two groups, despite BWS and DGF individual adaptation, whereas Hγ ERSPs reflected the presence of VR rather than to motor practice.
A brief ϑ-ERS (at the beginning of the gait cycle) was present in both groups. It is hypothesizable that ϑ-ERS in a non-specific event during gait and it is probably related to the sensorimotor area demand, the basic locomotor control, and the timing of muscular activation patterns [138, 139].
ERSPs were lateralized in the affected hemisphere in the RAGT + VR group but not in the group RAGT-VR, despite the lesion localization was similar in both the groups, with the exception of parieto-occipital ERSPs, which were bilateral in both groups, as formerly reported . In fact, visuomotor information processing is distributed symmetrically during walking , except for some specialized areas located in the right hemisphere, which are crucial for the closed-loop aspects of the movements depending on the sensory feedback . We may argue that the bi-hemispheric distribution of ERSPs in the RAGT-VR may depend on a dysfunctional reshape of interhemispheric connectivity, which was instead recovered, at least partially, in the RAGT + VR group [142,143,144,145,146,147].
As limiting factors in our work, we have to acknowledge that patients were provided with objects appearing in different corners of the screen during RAGT + VR. This fact may force eye-movements planning, which is expressed by a decreased α/β power within parietooccipital regions . Nonetheless, the extent of difference in brain activation between RAGT + VR and RAGT-VR is sufficiently high to exclude a biasing effect of the activity related to saccades on our data. Moreover, it has been shown that first-person perspective is superior to third-person perspective VR [149, 150], owing to an enhanced feeling of agency . We could speculate that the increased performances might have induced a greater feeling of agency in the third person perspective. However, studies comparing first- and person perspective are needed to confirm this issue.
Further, the clinical improvement we reported may also depend on factors not directly related to the ERSPs, including a stronger motivation for active participation in the movement provided by the VR , as suggested by the few episodes of drowsiness and the high sense of entrainment in the VR setting. A stronger motivation is, however, of notable importance, given that it allows the patient to exercise more regularly, precisely, and intensely [153,154,155,156] and, at least indirectly by enhancing the voluntary drive, to improve motor planning, learning and execution [157, 158]. Thus, our results show anyway the possible benefit of goal directed walking tasks that recruit brain areas involved in motor planning, learning and execution by using VR.
Finally, we employed a relatively low-speed RAGT, which could have affected the timing of muscle activation and amplitude, thus potentially reducing the level of sensorimotor cortex activation. Nonetheless, we preferred to adopt a low-speed RAGT to avoid excessive EEG contamination due to movement artefacts.
In conclusion, our findings suggest that VR feedback during RAGT elicits stronger cortical activations within the fronto-parietooccipital areas potentially belonging to the MNS, and involved in motor intention and planning. These activations were paralleled by an evident improvement in walking ability. We may thus argue the use more demanding and interactive task during RAGT by using VR may be of benefit to the patients with stroke. Moreover, monitoring the EEG in this context allows clinicians to realize better patient-tailored rehabilitative approaches.
Body weight support
Fast Fourier transform
Leg guidance force
Hamilton Rating Scale for Depression
Independent component analysis
Low Resolution Brain Electromagnetic Tomography
Modified Ashworth Scale
Mirror neuron system
Muscle Research Council
Magnetic resonance imaging
Tinetti Performance Oriented Mobility Assessment
Power spectral density
Robot-assisted gait training
Visual analogue scale
Holden MK. Virtual environments for motor rehabilitation: review. Cyber Psychol Behav. 2005;8:187–211.
Riva G, Anguera MT, Wiederhold BK, Mantovani F. From communication to presence: cognition, emotions and culture towards the ultimate communicative experience. Amsterdam: IOS Press; 2006.
Nudo R, Milliken G, Jenkins W, Merzenich M. Use-dependent alterations of movement representations in primary motor cortex of adult squirrel monkeys. J Neurosci. 1996;16(2):785–807.
Boyd L, Vidoni E, Wessel B. Motor learning after stroke: is skill acquisition a prerequisite for contralesional neuroplastic change. Neurosci Lett. 2010;482:21–5.
Nudo R, Milliken G. Reorganization of movement representations in primary motor cortex following focal ischemic infarcts in adult squirrel monkeys. J Neurophysiol. 1996;75(5):2144–9.
Plautz E, Milliken G, Nudo R. Effects of repetitive motor training on movement representations in adult squirrel monkeys: role of use versus learning. Neurobiol Learn Mem. 2000;74:27–55.
Sveistrup H. Motor rehabilitation using virtual reality. J Neuroeng Rehabil. 2004;11:10.
Fu MJ, Knutson JS, Chae J. Stroke rehabilitation using virtual environments. Phys Med Rehabil Clin N Am. 2015;26:747–57.
Calabrò RS, Russo M, Naro A, De Luca R, Leo A, Tomasello P, Molonia F, Dattola V, Bramanti A, Bramanti P. Robotic gait training in multiple sclerosis rehabilitation: Can virtual reality make the difference? Findings from a randomized controlled trial. J Neurol Sci. 2017;377:25–30.
Kwakkel G, Kollen B, Wagenaar R. Therapy impact on functional recovery in stroke rehabilitation: a critical review of the literature. Physiotherapy. 1999;85(7):377–91.
Rossini P, Forno G. Neuronal post-stroke plasticity in the adult. Restorative Neurology and Neurosci. 2004;22(3-5):193–206.
Maier I, Schwab M. Sprouting, regeneration and circuit formation in the injured spinal cord: factors and activity. Philos Trans of the R Soc B: Biol Sci. 2006;361(1473):1611–34.
Nudo R, Plautz E, Frost S. Role of adaptive plasticity in recovery of function after damage to motor cortex. Muscle Nerve. 2001;24(8):1000–19.
Stroemer R, Kent T, Hulsebosch C. Neocortical neural sprouting, synaptogenesis, and behavioral recovery after neocortical infarction in rats. Stroke. 1995;26(11):2135–44.
Pavone EF, Tieri G, Rizza G, Tidoni E, Grisoni L, Aglioti SM. Embodying others in immersive virtual reality: electro-cortical signatures of monitoring the errors in the actions of an avatar seen from a first-person perspective. J Neurosci. 2016;36:268–79.
Fadiga L, Craighero L. Electrophysiology of action representation. Clin Neurophysiol. 2004;21:157–69.
Buccino G, Binkofski F, Fink GR, Fadiga L, Fogassi L, Gallese V, et al. Action observation activates premotor and parietal areas in a somatotopic manner: an fMRI study. Eur J Neurosci. 2001;13:400–4.
Ramachandran VS, Hirstein W. The perception of phantom limbs. Brain. 1998;121:1603–30.
Pavlova M, Lutzenberger W, Sokolov A, Birbaumer N. Dissociable cortical processing of recognizable and non-recognizable biological movement: analyzing gamma MEG activity. Cereb Cortex. 2004;14:181–8.
Modroño C, Navarrete G, Rodríguez-Hernández AF, González-Mora JL. Activation of the human mirror neuron system during the observation of the manipulation of virtual tools in the absence of a visible effector limb. Neurosci Lett. 2013;555:220–4.
Buccino G, Solodkin A, Small S. Functions of the mirror neuron system: implications for neurorehabilitation. Cogn Behav Neurol. 2006;19:55–63.
Laver K, George S, Thomas S, Deutsch JE, Crotty M. Virtual reality for stroke rehabilitation: an abridged version of a Cochrane review. Eur J Phys Rehabil Med. 2015;51:497–506.
Krakauer J. Motor learning: its relevance to stroke recovery and neurorehabilitation. Curr Opin Neurol. 2006;19:84–90.
Jäncke L, Cheetham M, Baumgartner T. Virtual reality and the role of the prefrontal cortex in adults and children. Front Neurosci. 2009;3:52–9.
Kravitz DJ, Saleem KS, Baker CI, Mishkin M. A new neural framework for visuospatial processing. Nat Rev Neurosci. 2011;12:217–30.
Miller R. Theory of the normal waking EEG: from single neurones to waveforms in the alpha, beta and gamma frequency ranges. Int J Psychophysiol. 2007;64:18–23.
Sharma N, Classen J, Cohen LG. Neural plasticity and its contribution to functional recovery. Handb Clin Neurol. 2013;110:3–12.
Oberman LM, Pineda JA, Ramachandran VS. The human mirror neuron system: a link between action observation and social skills. Soc Cogn Affect Neurosci. 2007;2(1):62–6.
Zarka D, Cevallos C, Petieau M, Hoellinger T, Dan B, Cheron G. Neural rhythmic symphony of human walking observation: upside-down and uncoordinated condition on cortical theta, alpha, beta and gamma oscillations. Front Syst Neurosci. 2014;8:169.
Cheron G, Petit G, Cheron J, et al. Brain oscillations in sport: toward EEG biomarkers of performance. Front Psychol. 2016;7:246.
Kazuhiro U, Masataka Y, Hongbo L, Chi Z. Effect of motor intensity on motion imagery with electroencephalogram signal analysis in mirror neuron system. J Neurosci Neuroeng. 2017;4(1):38–43.
McGarry LMJ. The role of the mirror neuron system in bottom-up and top-down perception of human action. Canada: Toronto University Press; 2015.
Filimon F, Rieth CA, Sereno MI, Cottrell GW. Observed, executed, and imagined action representations can be decoded from ventral and dorsal areas. Cereb Cortex. 2015;25(9):3144–58.
Ferreira Dos Santos L, Christ O, Mate K, Schmidt H, Krüger J, Dohle C. Movement visualization in virtual reality rehabilitation of the lower limb: a systematic review. Biomed Eng Online. 2016;15(Suppl 3):144.
de Rooij IJ, van de Port IG, Meijer JG. Effect of virtual reality training on balance and gait ability in patients with stroke: systematic review and meta-analysis. Phys Ther. 2016;96(12):1905–18.
Chan L, Heinemann AW, Roberts J. Elevating the quality of disability and rehabilitation research: mandatory use of the reporting guidelines. Ann Phys Rehabil Med. 2014;57:558–60.
Medical Research Council. Aids to the examination of the peripheral nervous system, memorandum no. 45. London: Pengragon House; 1978.
Bohannon RW, Smith MB. Interrater reliability of a modified Ashworth scale of muscle spasticity. Phys Ther. 1987;67(2):206–7.
Foreman B, Claassen J. Quantitative EEG for the detection of brain ischemia. Crit Care. 2012;16(2):216.
You RX, McNeil JJ, O'Malley HM, Davis SM, Thrift AG, Donnan GA. Risk factors for stroke due to cerebral infarction in young adults. Stroke. 1997;28(10):1913–8.
Song GB, Park EC. Effect of virtual reality games on stroke patients' balance, gait, depression, and interpersonal relationships. J Phys Ther Sci. 2015;27(7):2057–60.
de Rooij IJ, van de Port IG, Meijer JG. Effect of virtual reality training on balance and gait ability in patients with stroke: systematic review and meta-analysis. J Neuroeng Rehabil. 2015;12:46.
Darekar A, McFadyen BJ, Lamontagne A, Fung J. Efficacy of virtual reality-based intervention on balance and mobility disorders post-stroke: a scoping review. Biomed Res Int. 2016;2016:7309272.
Ling C, Wai LAL, Yu RM, Ming HD, Qiang L, Hai L, et al. Effect of virtual reality on postural and balance control in patients with stroke: a systematic literature review. Biomed Res Int. 2016;2016:7309272.
Riener R, Lunenburger L, Jezernik S, Anderschitz M, Colombo G, Dietz V. Patient cooperative strategies for robot-aided treadmill training: first experimental results. IEEE Trans Neur Syst Reh. 2005;13(3):380–94.
Riener R, Lunenburger L, Colombo G. Human-centered robotics applied to gait training and assessment. J Rehabil Res Dev. 2006;43(5):679–94.
Wirz M, Zemon DH, Rupp R, Scheel A, Colombo G, Dietz V, et al. Effectiveness of automated locomotor training in patients with chronic incomplete spinal cord injury: a multicenter trial. Arch Phys Med Rehabil. 2005;86(4):672–80.
Hornby TG, Zemon DH, Campbell D. Robotic-assisted, body-weight-supported treadmill training in individuals following motor incomplete spinal cord injury. Phys Ther. 2005;85(1):52–66.
Light GA, Williams LE, Minow F, Sprock J, Rissling A, Sharp R, Swerdlow NR, Braff DL. Electroencephalography (EEG) and Event-Related Potentials (ERP’s) with Human Participants. Curr Protoc Neurosci. 2010;6:2524.
Lisi G, Morimoto J. EEG single-trial detection of gait speed changes during treadmill walk. PLoS One. 2015;10(5):e0125479.
Kilicarslan A, Grossman RG. Contreras-Vidal JL. A robust adaptive denoising framework for real-time artifact removal in scalp EEG measurements. J Neural Eng. 2016;13:026013.
Fruhstorfer H, Langanke P, Meinzer K, Peter JH, Pfaff U. Neurophysiological vigilance indicators and operational analysis of a train vigilance monitoring device: a laboratory and field study. In: Mackie RR, editor. Vigilance. New York: Plenum Press; 1977. p. 147–62.
Makeig S, Debener S, Onton J, Delorme A. Mining event-related brain dynamics. Trends Cogn Sci. 2004;8(5):204–10.
Wagner J, Solis-Escalante T, Grieshofer P, Neuper C, Müller-Putz G, Scherer R. Level of participation in robotic-assisted treadmill walking modulates midline sensorimotor EEG rhythms in able-bodied subjects. NeuroImage. 2012;63:1203–11.
Onton J, Westerfield M, Townsend J, Makeig S. Imaging human EEG dynamics using independent component analysis. Neurosci Biobehav Rev. 2006;30:808–22.
Gwin JT, Gramann K, Makeig S, Ferris DP. Electrocortical activity is coupled to gait cycle phase during treadmill walking. NeuroImage. 2010;54:1289–96.
Delorme A, Makeig S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Methods. 2004;134:9–21.
Makeig S. Auditory event-related dynamics of the EEG spectrum and effects of exposure to tones. Electroencephalogr Clin Neurophysiol. 1993;86:283–93.
Perry J, Burnfield JM. Gait analysis, normal and pathological function. J Sports Sci Med. 2010;9(2):353.
Delorme A, Sejnowski T, Makeig S. Enhanced detection of artifacts in EEG data using higher-order statistics and independent component analysis. NeuroImage. 2007;34(4):1443–9.
Nathan K, Contreras-Vidal JL. Negligible motion artifacts in scalp electroencephalography (EEG) during treadmill walking. Front Hum Neurosci. 2016;9:708.
Castermans T, Duvinage M, Cheron G, Dutoit T. Towards effective non-invasive brain-computer interfaces dedicated to gait rehabilitation systems. Brain Sci. 2014;4:1–48.
Pfurtscheller G, Neuper C. Event-related synchronization of mu rhythm in the EEG over the cortical hand area in man. Neurosci Lett. 1994;174:93–6.
Andrew C, Pfurtscheller G. On the existence of different alpha band rhythms in the hand area of man. Neurosci Lett. 1997;222(2):103–6.
Frenkel-Toledo S, Bentin S, Perry A, Liebermann DG, Soroker N. Dynamics of the EEG power in the frequency and spatial domains during observation and execution of manual movements. Brain Res. 2013;1509:43–57.
Frenkel-Toledo S, Bentin S, Perry A, Liebermann DG, Soroker N. Mirror-neuron system recruitment by action observation: effects of focal brain damage on mu suppression. NeuroImage. 2014;87:127–37.
Taborri J, Palermo E, Rossi S, Cappa P. Gait Partitioning Methods: A Systematic Review. Sensors (Basel, Switzerland). 2016;16:66.
Rajmohan V, Mohandas E. Mirror neuron system. Indian J Psychiatry. 2007;49:66–69.
Kilner JM, Lemon RN. What we know currently about mirror neurons. Curr Biol. 2013;23(23):R1057–62.
Schippers MB, Keysers C. Mapping the flow of information within the putative mirror neuron system during gesture observation. NeuroImage. 2011;57(1):37–44.
Braadbaart L, Williams JH, Waiter GD. Do mirror neuron areas mediate mu rhythm suppression during imitation and action observation? Int J Psychophysiol. 2013;89(1):99–105.
Demiralp T, Bayraktaroglu Z, Lenz D, Junge S, Busch NA, Maess B, et al. Gamma amplitudes are coupled to theta phase in human EEG during visual perception. Int J Psychophysiol. 2007;64(1):24–30.
de Lange FP, Jensen O, Bauer M, Toni I. Interactions between posterior gamma and frontal alpha/beta oscillations during imagined actions. Front Hum Neurosci. 2008;2:7.
Pfurtscheller G, Neuper C. Event-related synchronization of mu rhythm in the EEG over the cortical hand area in man. Neurosci Lett. 2011;174:93–6.
Pascual-Marqui RD. Review of methods for solving the EEG inverse problem. Int J Bioelectromagn. 1999;1:75–86.
Pascual-Marqui RD, Michel CM, Lehmann D. Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. Int J Psychophysiol. 1994;18:49–65.
Fuchs M, Drenckhahn R, Wischmann HA, Wagner M. An improved boundary element method for realistic volume-conductor modeling. IEEE Trans Biomed Eng. 1998;45:980–97.
Fuchs M, Kastner J, Wagner M, et al. A standardized boundary element method volume conductor model. Clin Neurophysiol. 2002;113:702–12.
Basile LF, Yacubian J, Castro CC, Grattaz WF. Widespread electrical cortical dysfunction in schizophrenia. Schizophr Res. 2004;69:255–66.
Yao J, Dewald JP. Evaluation of different cortical source localization methods using simulated and experimental EEG data. NeuroImage. 2005;25:369–82.
Chen HM, Hsieh CL. Sing Kai lo, Liaw LJ, Chen SM, Lin JH. The test-retest reliability of 2 mobility performance tests in patients with chronic stroke. Neurorehabil Neural Repair. 2007;21(4):347–52.
Canbek J, Fulk G, Nof L, Echternach J. Test-retest reliability and construct validity of the tinetti performance-oriented mobility assessment in people with stroke. J Neurol Phys Ther. 2013;37(1):14–9.
Monge-Pereira E, Molina-Rueda F, Rivas-Montero FM, Ibáñez J, Serrano JI, Alguacil-Diego IM, et al. Electroencephalography as a post-stroke assessment method: an updated review. Neurologia. 2017;32(1):40–9.
Knaepen K, Mierau A, Swinnen E, Fernandez Tellez H, Michielsen M, Kerckhofs E, et al. Human-robot interaction: does robotic guidance force affect gait-related brain dynamics during robot-assisted treadmill walking? PLoS One. 2015;10(10):e0140626.
Dobkin BH, Duncan PW. Should body weight-supported treadmill training and robotic-assistive steppers for locomotor training trot back to the starting gate? Neurorehabil Neural Repair. 2012;26:308–17.
Wutzke CJ, Mercer VS, Lewek MD. Influence of lower extremity sensory function on locomotor adaptation following stroke: a review. Top Stroke Rehabil. 2013;20(3):233–40.
Park W, Kwon GH, Kim YH, Lee JH, Kim L. EEG response varies with lesion location in patients with chronic stroke. J Neuroeng Rehabil. 2016;13:21.
Wu J, Quinlan EB, Dodakian L, et al. Connectivity measures are robust biomarkers of cortical function and plasticity after stroke. Brain. 2015;138(8):2359–69. doi:10.1093/brain/awv156.
Perry J. Gait analysis: normal and pathological function. Thorofare: Slack Incorporated; 1992.
Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001;29:1165–88.
Pfurtscheller G, Lopes da Silva F. Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin Neurophysiol. 1999;110:1842–57.
Crone NE, Miglioretti DL, Gordon B, Sieracki JM, Wilson MT, Uematsu S, et al. Functional mapping of human sensorimotor cortex with electrocorticographic spectral analysis. I. Alpha and beta event-related desynchronization. Brain. 1998;121:2271–99.
Gilbertson T, Lalo E, Doyle L, Di Lazzaro V, Cioni B, Brown P. Existing motor state is favored at the expense of new movement during 13-35 Hz oscillatory synchrony in the human corticospinal system. J Neurosci. 2005;25:7771–9.
Jenkinson N, Brown P. New insights into the relationship between dopamine, beta oscillations and motor function. Trends Neurosci. 2011;34:611–8.
Joundi RA, Jenkinson N, Brittain JS, Aziz TZ, Brown P. Driving oscillatory activity in the human cortex enhances motor performance. Curr Biol. 2012;22:403–7.
Klimesch W, Sauseng P, Hanslmayr S. EEG alpha oscillations: the inhibition-timing hypothesis. J Neurosci. 2007;53:63–88.
Miller KJ, Leuthardt EC, Schalk G, Rao RP, Anderson NR, Moran DW, et al. Spectral changes in cortical surface potentials during motor movement. J Neurosci. 2007;27:2424–32.
Neuper C, Pfurtscheller G. Event-related dynamics of cortical rhythms: frequency-specific features and functional correlates. Int J Psychophysiol. 2001;43:41–58.
Neuper C, Wörtz M, Pfurtscheller G. ERD/ERS patterns reflecting sensorimotor activation and deactivation. J Neurosci. 2007;27:2424–32.
Knill DC, Maloney LT, Trommershäuser J. Sensorimotor Processing and Goal-Directed Movement. J Vis. 2007;7:1–2.
Pfurtscheller G, Graimann B, Huggins JE, Levine SP, Schuh LA. Spatiotemporal patterns of beta desynchronization and gamma synchronization in corticographic data during self-paced movement. Clin Neurophysiol. 2003;114:1226–36.
Pfurtscheller G, Stancak J, Neuper C. Post-movement beta synchronization. A correlate of an idling area. Electroencephalogr. Clin. Neurophysiology. 1996;98:281–93.
Alegre M, Labarga A, Gurtubay IG, Iriarte J, Malanda A, Artieda J. Movement-related changes in cortical oscillatory activity in ballistic, sustained and negative movements. Exp Brain Res. 2003;148:17–25.
Brown P. Cortical drives to human muscle: the piper and related rhythms. Prog Neurobiol. 2000;60:97–108.
Cavanna AE, Trimble MR. The precuneus: a review of its functional anatomy and behavioural correlates. Brain. 2006;129:564–83.
Mima T, Steger J, Schulman AE, Gerloff C, Hallett M. Electroencephalographic measurement of motor cortex control of muscle activity in humans. Clin Neurophysiol. 2000;111:326–37.
Müller-Putz P, Neuper C, Rupp R, Keinrath C, Gerner H, Pfurtscheller G. Event-related beta electroencephalographic changes during wrist movements induced by functional electrical stimulation of forearm muscles in man. Neurosci Lett. 2003;340:143–7.
Müller-Putz P, Zimmermann D, Graimann B, Nestinger K, Korisek G, Pfurtscheller G. Event-related beta EEG-changes during passive and attempted foot movements in paraplegic patients. Brain Res. 2007;1137:84–91.
Pfurtscheller G, Neuper C. Simultaneous EEG 10 Hz desynchronization and 40Hz synchronization during finger movements. Neuroreport. 1992;3:1057–60.
Pfurtscheller G, Neuper C, Kalcher J. 40-hz oscillations during motor behavior in man. Neurosci Lett. 1993;164:179–82.
Raethjen J, Govindan RB, Binder S, Zeuner KE, Deuschl G, Stolze H. Cortical representation of rhythmic foot movements. Brain Res. 2008;1236:79–84.
Sahyoun C, Floyer-Lea A, Johansen-Berg H, Matthews PM. Towards an understanding of gait control: brain activation during the anticipation, preparation and execution of foot movements. NeuroImage. 2004;21(2):568–75.
Gwin JT, Gramann K, Makeig S, Ferris DP. Electrocortical activity is coupled to gait cycle phase during treadmill walking. NeuroImage. 2011;54(2):1289–96.
Haggard P. Human volition: towards a neuroscience of will. Nat Rev Neurosci. 2008;9:934–46.
Haefeli J, Vögeli S, Michel J, Dietz V. Preparation and performance of obstacle steps: interaction between brain and spinal neuronal activity. Eur J Neurosci. 2011;33:338–48.
Koenraadt KL, Roelofsen EG, Duysens J, Keijsers NL. Cortical control of normal gait and precision stepping: an fNIRS study. NeuroImage. 2013;85:415–22.
Suzuki M, Miyai I, Ono T, Kubota K. Activities in the frontal cortex and gait performance are modulated by preparation. An fNIRS study. NeuroImage. 2008;39:600–7.
Suzuki M, Miyai I, Ono T, Oda I, Konishi I, Kochiyama T, et al. Prefrontal and premotor cortices are involved in adapting walking and running speed on the treadmill: an optical imaging study. NeuroImage. 2004;23:1020–6.
Capotosto P, Babiloni C, Romani GL, Corbetta M. Frontoparietal cortex controls spatial attention through modulation of anticipatory alpha rhythms. J Neurosci. 2009;29:5863–72.
Capotosto P, Babiloni C, Romani GL, Corbetta M. Differential contribution of right and left parietal cortex to the control of spatial attention: a simultaneous EEG-rTMS study. Cereb Cortex. 2012;22:446–54.
Donner TH, Siegel M. A framework for local cortical oscillation patterns. Trends Cogn Sci. 2011;15:191–9.
Hipp JF, Engel AK, Siegel M. Oscillatory synchronization in large-scale cortical networks predicts perception. Neuron. 2011;69:387–96.
Perfetti B, Moisello C, Landsness EC, Kvint S, Pruski A, Onofrj M, et al. Temporal evolution of oscillatory activity predicts performance in a choice-reaction time reaching task. J Neurophysiol. 2011;105:18–27.
Tombini M, Zappasodi F, Zollo L, Pellegrino G, Cavallo G, Tecchio F, et al. Brain activity preceding a 2D manual catching task. NeuroImage. 2009;47:1735–46.
Bastos AM, Usrey WM, Adams RA, Mangun GR, Fries P, Friston KJ. Canonical microcircuits for predictive coding. Neuron. 2012;76:695–711.
Ray S, Maunsell JHR. Differences in gamma frequencies across visual cortex restrict their possible use in computation. Neuron. 2010;67:885–96.
Brunel N, Wang XJ. What determines the frequency of fast network oscillations with irregular neural discharges? I. Synaptic dynamics and excitation-inhibition balance. J Neurophysiol. 2003;90:415–30.
Roberts MJ, Lowet E, Brunet NM, Ter Wal M, Tiesinga P, Fries P, et al. Robust gamma coherence between macaque V1 and V2 by dynamic frequency matching. Neuron. 2013;78:523–36.
Desmurget M, Reilly KT, Richard N, Szathmari A, Mottolese C, Sirigu A. Movement intention after parietal cortex stimulation in humans. Science. 2009;324:811–3.
Drew T, Andujar JE, Lajoie K, Yakovenko S. Cortical mechanisms involved in visuomotor coordination during precision walking. Brain Res Rev. 2008;57:199–211.
Fogassi L, Ferrari PF, Gesierich B, Rozzi S, Chersi F, Rizzolatti G. Parietal lobe: from action organization to intention understanding. Science. 2005;308:662–7.
Wagner J, Solis-Escalante T, Scherer R, Neuper C, Muller-Putz G. It's how you get there: walking down a virtual alley activates premotor and parietal areas. Front Hum Neurosci. 2014;8:93.
Seeber M, Scherer R, Wagner J, Solis-Escalante T, Muller-Putz GR. EEG beta suppression and low gamma modulation are different elements of human upright walking. Front Hum Neurosci. 2014;8:485.
Fries P. Neuronal gamma-band synchronization as a fundamental process in cortical computation. Annu Rev Neurosci. 2009;32:209–24.
Presacco A, Goodman R, Forrester L, Contreras-Vidal JL. Neural decoding of treadmill walking from noninvasive electroencephalographic signals. J Neurophysiol. 2011;106:1875–87.
Moreno JC, Barroso F, Farina D, Gizzi L, Santos C, Molinari M, et al. Effects of robotic guidance on the coordination of locomotion. J Neuroeng Rehabil. 2013;10:79.
Aru J, Axmacher N, Do Lam ATA, Fell J, Elger CE, Singer W, et al. Local category-specific gamma band responses in the visual cortex do not reflect conscious perception. J Neurosci. 2012;32:14909–14.
Petersen TH, Willerslev-Olsen M, Conway BA, Nielsen JB. The motor cortex drives the muscles during walking in human subjects. J Physiol. 2012;590:2443–52.
Vidal JR, Perrone-Bertolotti M, Levy J, De Palma L, Minotti L, Kahane P, et al. Neural repetition suppression in ventral occipito-temporal cortex occurs during conscious and unconscious. NeuroImage. 2012;95:129–35.
Nobre AC, Sebestyen GN, Gitelman DR, Mesulam MM, Frackowiak RS, Frith CD. Functional localization of the system for visuospatial attention using positron emission tomography. Brain. 1997;120:515–33.
Serrien DJ, Ivry RB, Swinnen SP. Dynamics of hemispheric specialization and integration in the context of motor control. Nat Rev Neurosci. 2006;7:160–6.
Cao Y, D’Olhaberriague L, Vikingstad EM, et al. Pilot study of functional MRI to assess cerebral activation of motor function after poststroke hemiparesis. Stroke. 1998;29:112–22.
Chen R, Gerloff C, Hallett M, Cohen LG. Involvement of the ipsilateral motor cortex in finger movements of different complexities. Ann Neurol. 1997;41:247–54.
Mima T, Toma K, Koshy B, Hallett M. Coherence between cortical and muscular activities after subcortical stroke. Stroke. 2001;32:2597–601.
Gerloff C, Bushara K, Sailer A, et al. Multimodal imaging of brain reorganization in motor areas of the contralesional hemisphere of well recovered patients after capsular stroke. Brain. 2006;129:791–808.
Roland PE, Skinhoj E, Lassen NA, Larsen B. Different cortical areas in man in organization of voluntary movements in extrapersonal space. J Neurophysiol. 1980;43:137–50.
Kim SG, Ashe J, Hendrich K, et al. Functional magnetic resonance imaging of motor cortex: hemispheric asymmetry and handedness. Science. 1993;261:615–7.
Snyder LH, Batista AP, Andersen RA. Coding of intention in the posterior parietal cortex. Nature. 1997;386:167–70.
Petkova VI, Khoshnevis M, Ehrsson HH. The perspective matters! Multisensory integration in ego-centric reference frames determines full-body ownership. Front Psychol. 2011;2:35.
Slater M, Spanlang B, Sanchez-Vives MV, Blanke O. First per- son experience of body transfer in virtual reality. PLoS One. 2010;5:e10564.
Tsakiris M, Longo MR, Haggard P. Having a body versus moving your body: neural signatures of agency and body-ownership. Neuropsychologia. 2010;48:2740–9.
Meadows CC, Gable PA, Lohse KR, Miller MW. Motivation and motor cortical activity can independently affect motor performance. Neuroscience. 2016;339:174–9.
Brütsch K, Schuler T, Koenig A, Zimmerli L, Merillat-Koeneke S, Lünenburger L, et al. Influence of virtual reality soccer game on walking performance in robotic assisted gait training for children. J Neuroeng Rehabil. 2010;7:15.
Brütsch K, Koenig A, Zimmerli L, Merillat-Koeneke S, Riener R, Jäncke L, et al. Virtual reality for enhancement of robot-assisted gait training in children with neurological gait disorders. J Rehabil Med. 2011;43:493–9.
Schuler T, Brütsch K, Müller R, van Hedel H, Meyer-Heim A. Virtual realities as motivational tools for robotic assisted gait training in children: a surface electromyography study. NeuroRehabilitation. 2011;28:401–11.
Zimmerli L, Jacky M, Lünenburger L, Riener R, Bolliger M. Increasing patient engagement during virtual reality based motor rehabilitation. Arch Phys Med. 2013;94:1737–46.
Kaelin-Lang A, Sawaki L, Cohen LG. Role of voluntary drive in encoding an elementary motor memory. J Neurophysiol. 2005;93:1099–103.
Lotze M, Braun C, Birbaumer N, Anders S, Cohen LG. Motor learning elicited by voluntary drive. Brain. 2003;126:866–72.
Availability of data and materials
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
No funding to be reported.
RSC, AN, MR, AB, PB: 1a, 2b, 3, and 4; AL, RDL, TB, AB, GLR: 1b, 2a, 3, and 4 (legend: 1a Substantial contributions to the conception and design of the work, and interpretation of data; 1b acquisition and analysis of data. 2a Drafting the work; 2b revising the work critically for important intellectual content. 3 Final approval of the version to be published. 4 Agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved). All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Informed consent was obtained from all individual participants included in the study.
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Calabrò, R.S., Naro, A., Russo, M. et al. The role of virtual reality in improving motor performance as revealed by EEG: a randomized clinical trial. J NeuroEngineering Rehabil 14, 53 (2017). https://doi.org/10.1186/s12984-017-0268-4