Brisk heart rate and EEG changes during execution and withholding of cue-paced foot motor imagery

Introduction

Cue-paced motor imagery (MI) is one of the most frequently used mental strategies in brain-computer interface (BCI) applications (Pfurtscheller and Neuper, 2001; Wolpaw et al., 2002; Faller et al., 2012) for severely handicapped patients (Neuper et al., 2003; Müller-Putz and Pfurtscheller, 2008; Pichiorri et al., 2011) and rehabilitation after stroke (Birbaumer et al., 2008; Kaiser et al., 2012). The major benefit of such a strategy is that the subject has to pay no attention to any externally presented stimuli during self-paced operation at free will. A prerequisite for a good performance, however, is that the subject experience some cue-paced training sessions without or with feedback. Such training consists of a cue stimulus presentation in slightly varying intervals of a few seconds (Pfurtscheller et al., 2008b; Fazli et al., 2012). The location, color or the form of cue indicates either the type of MI (e.g., right vs. left hand) or whether MI has to be performed or withheld (Solis-Escalante et al., 2012). This procedure results in a number of cortical and subcortical processing steps starting with the stimulus perception and decision making and ending with the requested execution of the MI task, whereby the user focuses attention to visually presented cue stimuli. Recently it has been reported that visually cued imagination of left/right hand revealed the highest classification accuracy ~1 s after cue-onset (Fazli et al., 2012). A similar early classification peak during discrimination between visually cued hand and foot MI was reported by Pfurtscheller et al. (2008b). Remarkable in these studies is not only that the highest classification accuracy was found so early after visual stimulation, but also that this early classification peak was observed in nearly every subject. This stability of the early MI classification within and between subjects is of interest and needs further investigation not only in the cerebral (EEG) but also in the autonomic system (heart rate, HR), because both systems are closely linked together. The first who emphasized the close interaction between brain and heart was Claude Bernard (1867), where he especially pointed to the mutual interaction between the two most important body organs (Darwin, 1872/1999; pp. 71–72). Recently, Thayer and Lane (2009) made an extensive review on the cortical control of cardiac activity and about the pathways by which the prefrontal cortex might influence the control of HR. The prefrontal cortex plays a dominant role in the temporal organization of action (Fuster et al., 2000; Haggard, 2005; Soon et al., 2008) whether physically executed or imagined.

Each externally presented stimulus results in a complex total body response, the orienting reflex (OR), first described by Pavlov in 1910 (see Pavlov, 1927) and developed in terms of psychophysiological measures by Sokolov (1960). The OR is an organismic reflex resulting in a range of changes (for details see Barry, 2006) such as EEG alpha desynchronization, cephalic vasodilatation, short-lasting HR deceleration and vasoconstriction in the periphery. The basis of this short-lasting HR deceleration over some seconds was identified in the brief prolongation of the cardiac cycle following stimulation reported by Lacey and Lacey (1978, 1980). Barry (1983) noted that the Laceys' vagally-mediated “primary bradycardia” was the beginning of the HR deceleration observed in the OR context, and suggested that this HR deceleration is an obligatory response marking the initial detection of the stimulus onset transient. Another response to external stimuli is the brisk, short-lasting (~500 ms) beta event-related desynchronization (ERD; Pfurtscheller and Lopes da Silva, 1999) after action-coded visual stimuli - a frequently reported phenomenon (Doyle et al., 2005; Pfurtscheller et al., 2008b; Zhang et al., 2008; Tzagarakis et al., 2010; Wang et al., 2010; Solis-Escalante et al., 2012) related to motor planning, response preparation, response inhibition, and response uncertainty.

The fact that only a minority of people are able to control a BCI properly, without extensive training, depends not only on factors like motivation, mental effort, and mood, but also on automatic operating processes, such as the fast interaction between brain and heart. The aim of this paper is first, to report in detail the brisk reaction of central beta oscillations originating in sensorimotor areas and the HR response after cue-paced movement imagination using a go-nogo paradigm (execution vs. withholding of motor imagery). While the MIgo task was clearly defined, the mental state during the MInogo trials is unknown. Participants were instructed not to move and to relax. Second, this paper aims to discuss how coupled phasic EEG desynchronization and HR deceleration can be interpreted, and provides an understanding of why the classification accuracy in a MI-based BCI can be highest early after the visual cue stimulus. Third, this paper underlines the importance of common EEG and HR changes in BCI research.

Materials and Methods

Subjects and Experimental Paradigm

The data of 16 healthy subjects (8 males and 8 females, age 23.6 ± 3.5 years) were recorded. Initially 20 subjects participated in this study; the data from four subjects were discarded because of EEG artifacts and recorded muscle contractions in the EMG during MI. All subjects were seated in a comfortable armchair one meter in front of a computer screen. The computer screen showed cues (duration 2 s) for a go (green circle) and nogo experiment (red circle). The interval between the cues was varied between 11 and 14 s.

The experimental paradigm consisted of two runs with cue-paced motor execution (ME) and three runs with cue-paced MI. The data were recorded in two blocks. In the first block, all participants completed two runs of motor execution. Then, after a pause of about 10 min, all participants completed three runs of MI. For the second block, the participants were instructed to perform kinesthetic MI, i.e., to imagine the sensation of moving their legs in response to the cue. Since none of the participants had previous experience with MI, the first experimental block was included to let the participants pay attention to the kinesthetic aspect of the task. Each run consisted of 40 trials with a go/nogo class probability of 50%. The participants' task was to execute or imagine a brisk movement (dorsiflexion) of both feet following a green circle (MEgo, MIgo), or to withhold the motor task (MEnogo, MInogo) if a red circle appeared.

The order of the blocks did not affect the “baseline” activity of the sensorimotor cortex, meaning that motor execution did not pre-activate the cortex, affecting the estimation of relative changes (ERD/ERS). Evidence against a possible effect of the block order comes from an analysis described in Solis-Escalante et al. (2012), in which the baseline activity preceding nogo trials in the conditions ME vs. MI was compared. No statistically significant differences were found. For further details, e.g., the frequency bands under analysis, see Solis-Escalante et al. (2012).

In this work only data from the MI sessions are reported. The protocol was approved by the Ethics committee of the Medical University of Graz and the subjects gave informed written consent before the experiments.

Signal Recording

The EEG was recorded with fifteen Ag-AgCl electrodes, arranged in three Laplacian derivations around electrode positions C3, Cz, and C4; overlaying the left hand, foot, and right hand sensorimotor cortex. Inter-electrode distance was 2.5 cm. Reference and ground electrodes were attached to the left and right mastoid, respectively. In addition to the EEG, the EMG was recorded from the tibialis anterior muscles in both legs using bipolar derivations. The ground electrode for EMG was attached to the right hip. The EMG was pre-processed before recording. Raw EMG signals were band-pass filtered between 1–1000 Hz, baseline-corrected, full-wave rectified, and integrated with a time constant of 100 ms. EEG and the integrated EMG were recorded with a biosignal amplifier (g.BSamp, Guger Technologies, Austria). Sampling rate was 250 Hz, with filters set between 0.5 and 100 Hz, and a notch filter (50 Hz), for both EEG and integrated EMG. The electrocardiogram (ECG) was recorded from a single bipolar derivation. The negative lead was attached to the chest at the left (mid) clavicular line and the 2nd intercostal space, and the positive lead was attached to the chest at the left midaxilar line and the 6th intercostal space. The ground electrode was placed on the right hip. Self-adhesive Ag-AgCl electrodes were used for these recordings. ECG was recorded with a band-pass filter between 2 and 100 Hz and also sampled at 250 Hz.

EEG Processing

The monopolar EEG data recorded with respect to the left mastoid were transformed with a Laplacian filter to improve the SNR of the signal. The inter-trial variance method (Pfurtscheller and Lopes da Silva, 1999) was used for quantification of the event-related (de)synchronization (ERD/ERS). Time-frequency representations were calculated to visualize the ERD/ERS patterns; only significant ERD/ERS values are displayed (Graimann et al., 2002). ERD/ERS quantification and visualization was managed by the BioSig toolbox (Schlögl and Brunner, 2008; available online http://biosig.sourceforge.net/).

The data were analyzed between 8 and 45 Hz in intervals of 2 Hz. Trials were filtered with a 5th order Butterworth filter, and the ensemble average was subtracted from individual trials. This operation reduces the contribution of phase-locked responses to the ERD/ERS quantification (Kalcher and Pfurtscheller, 1995). The trials were squared and averaged. A moving average window of 250 ms smoothed the ERD/ERS estimation. Significance intervals (p = 0.05) were computed with a Box-Cox transformation (Box and Cox, 1964). The reference interval for the relative power changes was 2–4 s prior to cue-onset. For the determination of the most reactive subject-specific beta bands we used the ERD/ERS maps computed from the motor execution task (MEgo) and the discriminability (Cohen's kappa) calculation (Cohen, 1960) between the beta rebound and the reference interval. The beta rebound is much larger during physical execution of movement compared to imagination of the same movement (Solis-Escalante et al., 2012). Hence, to obtain a reliable estimation of the beta band, we used the MEgo trials instead of the MIgo trials. User-specific bands were defined according to the strongest ERS response following the movement execution in MEgo trials. Details on selection of a user-specific band are given elsewhere (Solis-Escalante et al., 2012).

Calculation of the Heart Rate (HR) Changes

The first step in ECG processing is to detect the QRS complexes in the ECG signal. The interval between consecutive QRS-complexes marks the time from one heart contraction to the next (RR interval; RRI). The QRS-complex was detected automatically based on an algorithm using a filter bank to decompose the ECG signal into various sub-bands (Afonso et al., 1999; implemented in the Biosig toolbox). After detecting every QRS-complex, the RRI is estimated from every pair of successive complexes. The HR value during an RRI equals the inverse of the RRI concerned, and (in beats per minute) is linearly interpolated between two complexes; the resultant time series is averaged across trials from the same class. To obtain HR changes the individual HR signals were synchronously averaged relative to cue-onset (Pfurtscheller et al., 2007).

Artifact Rejection

Trials were rejected if the EMG from the left or right leg exceeded a threshold. For MInogo trials, this threshold was equal to the mean plus three times the standard deviation of the EMG at rest. For MIgo trials, the threshold was equal to the mean plus five times the standard deviation of the EMG at rest, since imagined movements can increase the tonus of the target muscles (Guillot et al., 2007). The “EMG at rest” period was defined as a 5 s interval before cue onset.

After computing the HR, additional trials were rejected by identification of extreme HR changes, i.e., outliers. The first derivative of the HR from each trial was computed, and the absolute value was analyzed using descriptive statistics. Trials were rejected in an iterative fashion: (i) compute the first and third quartile and the interquartile distance from all trials available (class-wise, e.g., separately for MIgo and MInogo); (ii) reject trials if any value of the HR's first derivative was: (1) less than the first quartile minus 1.5 times the interquartile distance, or (2) greater than the third quartile plus 1.5 times the interquartile distance; (iii) repeat until all outliers have been removed.

Statistical Analysis

The statistical significance of EEG band power changes (ERD, ERS) and HR changes was evaluated with a t-percentile bootstrap algorithm, using the implementation of the BioSig toolbox. Upper and lower confidence intervals were obtained from 500 resample repetitions, at confidence levels of 0.05 and 0.01.

Statistical analysis comprised a paired t-test for evaluating differences in beta ERD and HR between MIgo and MInogo and a Pearson's product moment correlation for evaluating possible relations between beta ERD and HR deceleration. Comparisons between MIgo and MInogo were analyzed for significant differences with a t-percentile bootstrap algorithm (as described above) for every sample in the time series.

Results

EEG Time-Frequency Maps During go and nogo Conditions

To obtain an overview of the dynamics of sensorimotor rhythms in the go (MI) and nogo condition (withholding of MI), time-frequency maps were calculated separately for each condition. Examples from three characteristic subjects (A–C) are shown in Figure 1. The maps in Figures 1A,B (left panels) display the classical response in the beta band during imagined dorsiflexions of both feet (MIgo) namely a long-lasting (~1–2 s), midcentrally localized, peri-imagery ERD during execution of the imagery process and a post-imagery ERS (beta rebound) after termination of this process. But subjects can also show an early beta ERD without following beta ERS during MIgo (example see Figure 1C, left panel), similar to the early ERD in the nogo condition (Figures 1A–C, right panels). Because of the variable duration of each imagery task and the alignment of the averaging process to cue-onset and not to the exact end of MI, the beta ERS (beta power increase) in the go condition is relatively long-lasting. In the nogo condition (Figure 1 right panels) the early beta ERD starts immediately after cue-onset, has a duration of ~500 ms and is terminated by a brisk beta ERS with a maximum ~1 s after cue-onset (panels A–C).

FIGURE 1

Figure 1. Time-frequency maps for MIgo and MInogo from three characteristic subjects (A–C) illustrating different reactivity patterns during execution and withholding of MI. Common for all three subjects is only the early beta ERD during MInogo (marked by stippled circles). Significant ERD values are displayed in “red” and significant ERS values in “blue.” The vertical line at second 5 indicates cue-onset. (A) Early beta ERD and beta ERS during MIgo and MInogo. (B) Early beta ERD and beta ERS during MIgo, and early ERD and weak beta ERS during MInogo. (C) Early beta ERD during MIgo and early beta ERD with beta ERS during MInogo.

The subject-specific beta band (mean ± SD) determined from the time-frequency maps and used for further analyses was between 18.8 ± 4.7 and 28.4 ± 4.8 Hz. The magnitudes and latencies of beta power minima (largest ERD magnitude) and HR minima measured within the first second after stimulation onset are summarized in Table 1 for both the go and nogo conditions.

TABLE 1

Table 1. Mean ± SD (16 subjects) for latency (ms) and magnitude (%) of early beta power minimum (beta ERD) and change of HR for go and no conditions.

The beta ERD was significantly (paired t-test) larger (39.82 vs. 32.00%; t = 5.37, p < 0.001) during MIgo compared to MInogo, but its peak latency was not significantly longer (471 vs. 397 ms; t = 1.05, ns). In contrast, HR deceleration was somewhat smaller in MIgo vs. MInogo (1.00 vs. 1.94 %, t = 2.03, p = 0.06); again, peak latencies (640 vs. 652 ms) did not differ significantly (t = 0.18, ns).

Changes of Beta Power and HR During Execution (go) and Withholding (nogo) of Motor Imagery

Grand averages of changes in percentage HR (upper panels) and EEG beta power (lower panels) for both conditions are displayed in Figure 2. Reference interval is 2–4 s before cue presentation. Beside the mean ± SE also significant changes (bootstrap p = 0.01) are indicated. Excepting the HR changes in the go condition, all other HR and EEG changes (HR deceleration and EEG beta desynchronization) were significant in the first second after cue presentation in go and nogo conditions (p < 0.01). While in the go condition the HR showed a significant (p < 0.01) acceleration after the brisk deceleration, no further significant HR changes were observed after the deceleration in the nogo condition.

FIGURE 2

Figure 2. Grand averages of percentage HR and beta power (ERD/ERS) changes for execution (MIgo) and withholding (MInogo) of foot motor imagery. In addition to the mean and the standard error (stippled lines) also the significant changes (bootstrap) are displayed (one asterisk vertically p < 0.05, and two asterisks vertically p < 0.01). Data of 16 subjects. Note the significant HR decrease in the nogo condition and the significant early beta ERD in both the go and nogo conditions.

For a direct comparison of the beta power decrease during motor planning in the go/nogo conditions after visual cue presentation with a similar 248-channel MEG study (Tzagarakis et al., 2010) the reference interval was changed to 0–1 s prior to cue onset [the standard reference interval for the processing of all EEG and HR data (Figure 2) was 2–4 s prior to cue onset in accordance with the work of Solis-Escalante et al. (2012)]. Common for both studies (Figures 3A,B) is the early beta power decrease from cue onset (Tzagarakis et al., 2010 reported a delay of 120 ms) not dependent on the type of information provided by the visual stimulus but rather determined by a constant visuomotor delay. The different behavior of the beta ERD between execution and withholding of MI starts between 200 and 300 ms after cue onset, indicating on the one hand that the cue stimulus-induced early beta ERD peaks at ~400 ms, and on the other hand that the cognitive task (execution of MI) is able to modulate this desychronization.

FIGURE 3

Figure 3. (A) Details of the grand average beta power time courses (ERD/ERS) of 16 subjects within the first second after cue onset for MIgo (dashed line) and MInogo (continuous line). Beta power was compared between MI go and MInogo (bootstrap, one asterisk p < 0.05, two asterisks p < 0.01) and revealed significant differences starting at 250 ms after cue onset. (B) Time courses of relative beta power changes after a visually instructed reaching task with one, two, or three possible target directions (with permission modified from Tzagarakis et al., 2010). Note the similar early beta ERD with a maximal magnitude at ~400 ms in EEG (A) and MEG (B) recordings during different motor tasks.

When searching for possible interpretations of the HR responses in the go and nogo condition we subtracted the nogo HR response from the go response (MIgo minus MInogo) and achieved a new response pattern termed “hypothetical” HR response (Figure 4A) with a HR acceleration starting immediately at cue onset and a maximum at ~2 s latter. Remarkable is the similarity between the “hypothetical” HR response in Figure 4A and the HR responses in a similar no/go study without any motor preparation displayed in Figure 4B (Lawrence and Barry, 2010). The brisk HR deceleration is the dominant feature in both studies.

FIGURE 4

Figure 4. Heart rate changes after go and nogo cues. (A) Grand average time courses of percentage HR changes aligned to cue onset referenced to an interval 1–0 s prior to cue onset (different from the text and analyses above). Displayed are the HR responses for MIgo, MInogo and their difference (hypothetical HR acceleration). Significant differences between conditions are marked with black asterisks (bootstrap, one asterisk p < 0.05, two asterisks p < 0.01). (B) Grand average HR response for silent counting of tones (COUNT) vs. no-count (NO COUNT) and hypothetical HR response (DIFFERENCE). Modified from Lawrence and Barry (2010).

Interdependency Between Early ERD and Brisk HR Deceleration

When testing linear inter-subject correlations between beta ERD and HR deceleration within the first second after cue-onset, no significant correlation was found for MInogo (Pearson correlation coefficient r = 0.05, p = 0.86; see Figure 5B); the correlation for MIgo was in the expected direction (greater deceleration with greater ERD) but failed to reach significance (r = −0.28, p = 0.29; see Figure 5A).

FIGURE 5

Figure 5. Scatter plot between early beta ERD (%) and brisk HR decrease (%) during motor imagery (A) and withholding of motor imagery (B).

Discussion

The following findings need special attention: (i) the slightly different reactivity patterns in the EEG during go and nogo conditions within the first second in the form of an early beta ERD, (ii) the different magnitudes of brisk HR decelerations in both conditions and (iii) that no significant correlation was found between early ERD and HR deceleration.

EEG Changes During Execution (go) and Withholding (nogo) of MI

In the EEG rhythms, various phenomena have to be differentiated during execution and withholding of MI. First the early beta ERD with a peak ~400 ms after cue-onset observed in both conditions, second the longer-lasting beta ERD associated with the conscious intention of MI (peri-imagery ERD) in the go condition, and third the beta rebound (post-imagery ERS) terminating the beta ERD. In the go condition the peri-imagery ERD becomes either superimposed on the early beta ERD within the first second (see Figures 1A,B) or the early beta ERD remains in the go condition (Figure 1C). Whether the absence of a longer-lasting beta ERD acts as indicator for a not correctly performed imagery task can only be speculated. The ERD represents a process to enable focal attention so that information processing may be optimized (Lopes da Silva, 1991) or, in other words, the ERD indicates a state of cortical activation (Gerloff et al., 1998; Pfurtscheller and Lopes da Silva, 1999; Leocani et al., 2001; Doyle et al., 2005). The beta ERS corresponds to a state of deactivation of the motor cortex circuitry and signifies active immobilization (Salmelin et al., 1995). In the nogo condition the beta ERD begins to recover 400 ms after cue onset and reaches its peak ERS at 1000 ms (see Figure 2 right lower panel and Figures 1A,C). Support for an activation of the foot representation area during MIgo came from fMRI studies. Physical execution but also mental simulation of foot/toe movement is not only accompanied by peri-imagery and mid-centrally localized beta ERD (Pfurtscheller et al., 2006a) but also by a positive fMRI BOLD signal localized to the foot representation area (Ehrsson et al., 2003).

In respect to the early ERD direct recordings from single neurons of the monkey‘s motor cortex are of interest. Such neurons start to discharge 60–80 ms after occurrence of a visual cue, indicating the direction of the upcoming movement (Georgopoulos et al., 1982). Preparatory changes in neuronal activity prior to movement execution have been documented also in monkey premotor and primary cortex 100 ms after the visual cue signal (Riehle and Requin, 1995). These studies show clearly that the motor cortex is engaged as early as about 100 ms after an informative visual cue indicates to execute or withhold a specific motor task. Similarly, in another monkey study, the beta activity displayed an early ERD in go and nogo conditions starting ~110 ms after visual stimulus presentation (Zhang et al., 2008). Therefore, we can conclude that the early phasic beta ERD observed in our data is a sign of preparatory activity in motor cortex networks triggered by the visual stimulus and modulated by the cognitive task.

Most interesting is that the early beta ERD displays a slightly but significant stronger magnitude (p < 0.01) and longer peak latency (n.s.) during MIgo than in MInogo (Table 1 and Figure 3A). The decision to execute or withhold MI is estimated to occur ~190 ms after stimulation (Zhang et al., 2008). Support for such an early stimulus induced preparatory activity in motor cortex networks came from two MEG studies (Tzagarakis et al., 2010; Wang et al., 2010). In Wang's study subjects were asked to imagine four different center-out movements using the wrist after visual target onset. The source-space multivariate test revealed that by 300–400 ms after target onset the cortical activity becomes highly modulated and can be used to decode the type of covert movement. Tzagarakis et al. reported on a 248-channel whole head MEG study while subjects performed a visually instructed reaching task with one, two, or three possible directions. The early beta ERD with a maximum at ~400 ms is heightened with decreased directional uncertainty (see Figure 3B). Remarkable is that the early beta ERD displays a peak at 500 ms in such two completely different studies, one with MEG and response uncertainty and the other with EEG and MI. This suggests that the early beta ERD can be maximal around 500 ms after cue onset not only in EEG and MEG recordings, but also can be modulated by different cognitive processes such as MI or response uncertainty.

HR Changes During Execution (go) and Withholding (nogo) of MI

The initial obligatory post-stimulus HR deceleration in the OR context mentioned in the Introduction is exemplified in response to the nogo cue, as illustrated in the right panel of Figure 2. This brisk HR response is significantly and clearly of larger magnitude when compared with the response during MI (go condition; Figure 2 left). Of interest is that this HR response is nearly identical with the “hypothetical” HR response difference in a similar go-nogo experiment where subjects were required to silently count the number of tones (go condition); in the nogo condition the subjects were told to relax during tone presentation (Lawrence and Barry, 2010; see also Figure 4B). These authors linked their hypothetical HR acceleration to the increased mental effort of the count condition. The hypothetical HR acceleration during MI is linked to the mental effort in the same fashion. A HR acceleration is characteristic for an increased mental effort during the simulation of a motor act (MI) (Decety et al., 1991; Oishi et al., 2000; Pfurtscheller et al., 2006b, 2008a).

Coupling Between Early Beta ERD and Brisk HR Deceleration

The beta ERD and concomitant HR deceleration after stimulation is not an isolated phenomenon but is also found during preparation for movement or even in the resting brain. The pre-movement beta ERD and the pre-movement HR deceleration are well documented phenomena (Damen and Brunia, 1987; Papakostopoulos et al., 1990; Florian et al., 1998; Pfurtscheller and Lopes da Silva, 1999). A slow cyclic (in intervals of ~10 s) central beta (alpha) power decrease and a nearly simultaneous HR decrease (beat-to-beat increase) in the resting state was reported recently (Pfurtscheller et al., 2012a,b). If the desynchronization of sensorimotor rhythms prior to movement can be interpreted as a correlate of preparatory activity then also the early beta ERD after cue presentation can be seen as stimulus triggered preparatory activity in the sensorimotor system. This preparatory activity, and also the HR deceleration, very likely result from routine processes operating automatically and unconsciously (Koch, 2003; Haggard, 2005). The EEG-HR coupling prior to voluntary movement, during rest and after visual stimulation can be seen as an example of the fast interaction between brain and heart mediated by pathways most likely from the prefrontal cortex to cardiovascular nuclei in the brain stem and vagally via the sino-atrial (SA) node to the heart (Thayer and Lane, 2009). The SA node responds very quickly to vagal (~150 ms latency) in contrast to sympathetic (~1–2 s) influences (Smyth et al., 1969).

Visual information is not only mediated by the lateral geniculate body, but also by the superior colliculus. The superior colliculus is important for the control of saccadic eye movements (Krauzlis et al., 2004), and is believed to be responsible for a transient increase in blood pressure and HR by drive of sympathetic activity (Iigaya et al., 2012) Electrical stimulation of the optic tract also increases blood pressure and HR, and inhibits baroreflex vagal bradycardia, mediated by the superior colliculus (Cheng et al., 2001). The participation of the superior colliculus in reflexive vagal mediated primary bradycardia (Lacey and Lacey, 1980; Barry, 1986) is therefore very unlikely.

That no significant correlation was found between early beta ERD and HR deceleration during MIgo (see Figure 5A) might be explained by different reasons. While the MIgo cardiac response can be conceptualized as an additional cardiac acceleration reflecting mental effort, superimposed on the reflex bradycardia triggered by the cue (see Figure 4A), a comparative examination of the corresponding beta power response profiles gives a more complex picture (see Figure 3A). Here the reflexive early beta ERD (peaking at ~400 ms) is superimposed on the peri-imagery beta ERD induced by the conscious mental simulation of foot movement. Although the early beta ERD with the following beta ERS in the nogo condition is an automatically induced response, only the early beta ERD can be considered as a component of the OR. The relationship between these differences in the EEG and HR response profiles during MIgo is complex and needs further exploration, possibly going beyond the linear correlation of responses in the short interval tested here. Another reason for the non-significant correlation could be that in all subjects the EEG was recorded at the vertex only (Cz: Laplacian derivation) although it is known from multichannel MEG recordings that the source of the beta band ERD can vary in the peri-Rolandic region across subjects (Tzagarakis et al., 2010). Beside the EEG also the ECG can induce some variability. The HR displays two preferential rhythms, one around 0.25 Hz (respiration) and one around 0.1 Hz (blood pressure or Mayer waves) (De Boer et al., 1985) as well as fluctuations induced by “central commands” (see review Benarroch, 1993; Thayer and Lane, 2009). Slow intrinsic HR fluctuations during rest can be as high as ~10 % of the mean HR (Pfurtscheller et al., 2012a) and mask the relatively small (~2%) brisk HR deceleration induced after stimulation.

Conclusion

The early beta ERD with the largest magnitude at ~400 ms is slightly, but significantly (p < 0.01) larger during execution (go) as compared to withholding of MI (nogo). At the same time the HR displays a significant (p < 0.01) brisk deceleration only during withholding of MI. Both the early beta ERD and the HR deceleration are the result of an automatic operating process and probably part of the OR (Barry, 2006). Of interest, however, is that both automatic reactions in brain and heart can be modulated by the increased mental effort associated with execution of MI. The calculation of the “hypothetical” HR acceleration introduced by Barry (1984) offers a new way to explore the mental effort during motor imagery early after stimulation. A high mental effort during imagined movements is accompanied by a HR increase (Decety et al., 1991; Oishi et al., 2000). Such indicator of mental effort could be used in a BCI to improve its performance (Pfurtscheller et al., 2010). For online monitoring of the mental effort in a cue-paced MI task, the averaged HR response of the MInogo condition obtained in a pre-experiment has to be subtracted from each ongoing HR response (calculation of the “hypothetical” HR). By this way the mental effort can be supervised and an intervention made if necessary.

The early beta ERD is an important brain feature that might be responsible for the early classification peak in imagery-based BCIs because of its somatotopic organization as demonstrated recently (Pfurtscheller et al., 2008b). The close relationship between early ERD and action-coded visual stimulation also was discussed elsewhere (Waldert et al., 2008; Wang et al., 2010). If the visual stimulus indicates only one specific type of MI the attention is focused on this specific imagery task and results in an early beta ERD in both conditions (go/nogo). The results suggest that the early classification peak in cue-paced BCIs is very likely the result of an automatic operating process induced by the preparatory visual cue stimulus. This hypothesis, however, needs to be further tested, possibly with an auditory cue stimulus.

The fast interaction between brain and heart is a basic factor in BCI research. HR slowing and concomitant pericentral EEG desynchronization is not only characteristic for the planning and/or preparation of a motor task (e.g., ME or MI), but also a dominant feature immediately after visual cue presentation. Still open for research and discussion is the link between slow (~0.1 Hz) spontaneous HR and EEG oscillations during awake rest (Pfurtscheller et al., 2012a,b). In this respect, the finding that the relation between early beta ERD and brisk HR deceleration is very likely the result of an automatic process, operating after visual cue presentation, is a small step to improve the success rate in BCI research and focus research attention on the importance of intrinsic HR and blood pressure oscillations at a frequency of around 0.1 Hz.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was partly supported by the European Union research project BETTER (ICT-2009.7.2-247935) and Neuro Center Styria (NCS).

References

Afonso, V., Tompkins, W., Nguyen, T., and Luo, S. (1999). ECG beat detection using filter banks. IEEE Trans. Biomed. Eng. 46, 192–202. doi: 10.1109/10.740882

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Barry, R. J. (1983). Primary bradycardia and the evoked cardiac response in the OR context. Physiol. Psychol. 11, 135–140.

Barry, R. J. (1984). The evoked cardiac response under processing load. Physiol. Psychol. 12, 35–40.

Barry, R. J. (1986). Heart rate deceleration to innocuous stimuli: an index of the orienting response or stimulus registration? Physiol. Psychol. 14, 42–48.

Barry, R. J. (2006). Promise versus reality in relation to the unitary orienting reflex: a case study examining the role of theory in psychophysiology. Int. J. Psychophysiol. 62, 353–366. doi: 10.1016/j.ijpsycho.2006.01.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Benarroch, E. E. (1993). The central autonomic network: functional organization, dysfunction, and perspective. Mayo Clin. Proc. 68, 988–1001. doi: 10.1016/S0025-6196(12)62272-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Birbaumer, N., Murguialday, A. R., and Cohen, L. (2008). Brain-computer interface in paralysis. Curr. Opin. Neurol. 21, 634–638. doi: 10.1097/WCO.0b013e328315ee2d

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Box, G. E. P., and Cox, D. R. (1964). An analysis of transformations. J. Roy. Statist. Soc. B 26, 211–252

Cheng, Z. B., Kobayashi, M., and Nosaka, S. (2001). Effects of optic tract stimulation on baroreflex vagal bradycardia in rats. Clin. Exp. Pharmacol. Physiol. 28, 721–728. doi: 10.1046/j.1440-1681.2001.03509.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cohen, J. (1960). A coefficient of agreement for nominal scales. Educ. Psychol. Meas. 20, 37–46

Damen, E., and Brunia, C. (1987). Changes in heart rate and slow brain potentials related to motor preparation and stimulus anticipation in a time estimation task. Psychophysiology 24, 700–713. doi: 10.1111/j.1469-8986.1987.tb00353.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Darwin, C. (1872/1999). The Expression of the Emotions in Man and Animals. London: Harper Collins.

De Boer, R., Karemaker, J., and Strackee, J. (1985). Relationships between short-term blood-pressure fluctuations and heart-rate variability in resting subjects. II: A simple model. Med. Biol. Eng. Comput. 23, 359–364. doi: 10.1007/BF02441590

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Decety, J., Jeannerod, M., Germain, M., and Pastene, J. (1991). Vegetative response during imagined movement is proportional to mental effort. Behav. Brain Res. 42, 1–5. doi: 10.1016/S0166-4328(05)80033-6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Doyle, L. M., Yarrow, K., and Brown, P. (2005). Lateralization of event-related beta desynchronization in the EEG during pre-cued reaction time tasks. Clin. Neurophysiol. 116, 1879–1888. doi: 10.1016/j.clinph.2005.03.017

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ehrsson, H. H., Geyer, S., and Naito, E. (2003). Imagery of voluntary movement of fingers, toes, and tongue activates corresponding body-part-specific motor representations. J. Neurophysiol. 90, 3304–3316. doi: 10.1152/jn.01113.2002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Faller, J., Vidaurre, C., Solis Escalante, T., Neuper, C., and Scherer, R. (2012). Autocalibration and recurrent adaptation: towards a plug and play online ERD-BCI. IEEE Trans. Neural Syst. Rehabil. Eng. 20, 313–319. doi: 10.1109/TNSRE.2012.2189584

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fazli, S., Mehnert, J., Steinbrink, J., Curio, G., Villringer, A., Müller, K. R., et al. (2012). Enhanced performance by a hybrid NIRS-EEG brain computer interface. Neuroimage 59, 519–529. doi: 10.1016/j.neuroimage.2011.07.084

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Florian, G., Andrew, C., and Pfurtscheller, G. (1998). Do changes in coherence always reflect changes in functional coupling. Electroencephalogr. Clin. Neurophysiol. 106, 87–91. doi: 10.1016/S0013-4694(97)00105-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fuster, J. M., Bodner, M., and Kroger, J. K. (2000). Cross-modal and cross-temporal association in neurons of frontal cortex. Nature 405, 347–351. doi: 10.1038/35012613

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Graimann, B., Huggins, J. E., Levine, S. P., and Pfurtscheller, G. (2002). Visualization of significant ERD/ERS patterns in multichannel EEG and ECoG data. Clin. Neurophysiol. 113, 43–47. doi: 10.1016/S1388-2457(01)00697-6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Georgopoulos, A. P., Kalaska, J. F., Caminiti, R., and Massey, J. T. (1982). On the relations between the direction of two-dimensional arm movements and cell discharge in primate motor cortex. J. Neurosci. 2, 1527–1537.

Pubmed Abstract | Pubmed Full Text

Gerloff, C., Uenishi, N., and Hallett, M. (1998). Cortical activation during fast repetitive finger movements in humans: dipole sources of steady-state movement-related cortical potentials. J. Clin. Neurophysiol. 15, 502–513. doi: 10.1097/00004691-199811000-00009

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Guillot, A., Lebon, F., Rouffet, D., Champely, S., Doyon, J., and Collet, C. (2007). Muscular responses during motor imagery as a function of muscle contraction types. Int. J. Psychophysiol. 66, 18–27. doi: 10.1016/j.ijpsycho.2007.05.009

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Haggard, P. (2005). Conscious intention and motor cognition. Trends Cogn. Sci. 9, 290–295. doi: 10.1016/j.tics.2005.04.012

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Iigaya, K., Mueller-Ribeiro, F. C., Horiuchi, J., McDowall, L. M., Nalivaiko, E., Fontes, M. A., et al. (2012). Synchronized activation of sympathetic vasomotor, cardiac, and respiratory outputs by neurons in the midbrain colliculi. Am. J. Physiol. Regul. Integr. Comp. Physiol. 303, 599–610. doi: 10.1152/ajpregu.00205.2012

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kaiser, V., Daly, I., Pichiorri, F., Mattia, D., Müller-Putz, G. R., and Neuper, C. (2012). Relationship between electrical brain responses to motor imagery and motor impairment in stroke. Stroke 43, 2735–2740. doi: 10.1161/STROKEAHA.112.665489

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kalcher, J., and Pfurtscheller, G. (1995). Discrimination between phase-locked and non-phase-locked event-related EEG activity. Electroencephalogr. Clin. Neurophysiol. 94, 381–384. doi: 10.1016/0013-4694(95)00040-6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Koch, C. (2003). The Quest for Consciousness: a Neurobiological Approach. Greenwood Village, CO: Roberts and Co Publishers.

Krauzlis, R. J., Liston, D., and Carello, C. D. (2004). Target selection and the superior colliculus: goals, choices and hypotheseses. Vision Res. 44, 1445–1451. doi: 10.1016/j.visres.2004.01.005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lacey, B. C., and Lacey, J. I. (1978). Two-way communication between the heart and the brain. Significance of time within the cardiac cycle. Am. Psychol. 33, 99–113. doi: 10.1037/0003-066X.33.2.99

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lacey, B. C., and Lacey, J. I. (1980). Cognitive modulation of time-dependent primary bradycardia. Psychophysiology 17, 209–222. doi: 10.1111/j.1469-8986.1980.tb00137.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lawrence, C. A., and Barry, R. J. (2010). Cognitive processing effects on auditory event-related potentials and the evoked cardiac response. Int. J. Psychophysiol. 78, 100–106. doi: 10.1016/j.ijpsycho.2010.06.027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Leocani, L., Toro, C., Zhuang, P., Gerloff, C., and Hallett, M. (2001). Event-related desynchronization in reaction time paradigms: a comparison with event-related potentials and corticospinal excitability. Clin. Neurophysiol. 112, 923–930. doi: 10.1016/S1388-2457(01)00530-2

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lopes da Silva, F. (1991). Neural mechanisms underlying brain waves: from neural membranes to networks. Electroencephalogr. Clin. Neurophysiol. 79, 81–93. doi: 10.1016/0013-4694(91)90044-5

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Müller-Putz, G., and Pfurtscheller, G. (2008). Control of an electrical prosthesis with an SSVEP-based BCI. IEEE Trans. Biomed. Eng. 55, 361–364. doi: 10.1109/TBME.2007.897815

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Neuper, C., Müller, G., Kübler, A., Birbaumer, N., and Pfurtscheller, G. (2003). Clinical application of an EEG-based brain-computer interface, a case study in a patient with severe motor impairment. Clin. Neurophysiol. 114, 399–409. doi: 10.1016/S1388-2457(02)00387-5

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Oishi, K., Kasai, T., and Maeshima, T. (2000). Autonomic response specificity during motor imagery. J. Physiol. Anthropol. Appl. Hum. Sci. 19, 255–261. doi: 10.2114/jpa.19.255

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Papakostopoulos, D., Banerji, N., and Pocock, P. (1990). Performance, EMG, brain electrical potentials and heart rate change during a self-paced skilled motor task in parkinson's disease. J. Psychophysiol. 4, 163–183.

Pavlov, I. P. (1927). Conditioned Reflexes. New York, NY: Dover.

Pfurtscheller, G., Allison, B. Z., Brunner, C., Bauernfeind, G., Solis-Escalante, T., Scherer, R., et al. (2010). The hybrid BCI. Front. Neurosci. 4:30. doi: 10.3389/fnpro.2010.00003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Bauernfeind, G., Neuper, C., and Lopes da Silva, F. H. (2012a). Does conscious intention to perform a motor act depend on slow prefrontal (de)oxyhemoglobin oscillations in the resting brain. Neurosci. Lett. 508, 89–94. doi: 10.1016/j.neulet.2011.12.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Daly, I., Bauernfeind, G., and Müller-Putz, G. (2012b). Coupling between intrinsic prefrontal HbO2 and central EEG beta power oscillations in the resting brain. PLoS ONE 7:e43640. doi: 10.1371/journal.pone.0043640

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Brunner, C., Schlögl, A., and Lopes da Silva, F. H. (2006a). Mu rhythm (de)synchronization and EEG single-trial classification of different motor imagery tasks. Neuroimage 31, 153–159. doi: 10.1016/j.neuroimage.2005.12.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Leeb, R., and Slater, M. (2006b). Cardiac responses induced during thought-based control of a virtual environment. Int. J. Psychophysiol. 62, 134–140. doi: 10.1016/j.ijpsycho.2006.03.001

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Grabner, R. H., Brunner, C., and Neuper, C. (2007). Phasic heart rate changes during word translation of different difficulties. Psychophysiology 44, 807–813. doi: 10.1111/j.1469-8986.2007.00553.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Leeb, R., Friedman, D., and Slater, M. (2008a). Centrally controlled heart rate changes during mental practice in immersive virtual environment: a case study with a tetralegic. Int. J. Psychophysiol. 68, 1–5. doi: 10.1016/j.ijpsycho.2007.11.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., Scherer, R., Müller-Putz, G. R., and Lopes da Silva, F. H. (2008b). Short-lived brain state after cued motor imagery in naive subjects. Eur. J. Neurosci. 28, 1419–1426. doi: 10.1111/j.1460-9568.2008.06441.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., and Lopes da Silva, F. H. (1999). Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin. Neurophysiol. 110, 1842–1857. doi: 10.1016/S1388-2457(99)00141-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pfurtscheller, G., and Neuper, C. (2001). Motor imagery and direct brain-computer communication (Invited paper). Proc. IEEE (Special Issue) Neural Eng. Merging Eng. Neurosci. 89, 1123–1134.

Pichiorri, F., De Vico Fallani, F., Cincotti, F., Babiloni, F., Molinari, M., Kleih, S. C., et al. (2011). Sensorimotor rhythm-based brain-computer interface training: the impact on motor cortical responsiveness. J. Neural Eng. 8:025020.

Pubmed Abstract | Pubmed Full Text

Riehle, A., and Requin, J. (1995). Neuronal correlates of the specification of movement direction and force in four cortical areas of the monkey. Behav. Brain Res. 70, 1–13. doi: 10.1016/0166-4328(94)00180-N

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Salmelin, R., Hämäläinen, M., Kajola, M., and Hari, R. (1995). Functional segregation of movement related rhythmic activity in the human brain. Neuroimage 2, 237–243. doi: 10.1006/nimg.1995.1031

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schlögl, A., and Brunner, C. (2008). BioSig: a free and open source software library for BCI research. Computer 41, 44–50. doi: 10.1109/MC.2008.407

CrossRef Full Text

Smyth, H. S., Sleight, P., and Pickering, G. W. (1969). Reflex regulation of arterial pressure during sleep in man. A quantitative method of assessing baroreflex sensitivity. Circ. Res. 24, 109–121. doi: 10.1161/01.RES.24.1.109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sokolov, E. N. (1960). “The central nervous system and behavior,” in Neuronal Models and the Orienting Reflex, ed Mary A. B. Brazier (New York, NY: Josiah Macy, Jr. Foundation), 187–276.

Solis-Escalante, T., Müller-Putz, G. R., Pfurtscheller, G., and Neuper, C. (2012). Cue-induced beta rebound during withholding of overt and covert foot movement. Clin. Neurophysiol. 123, 1182–1190. doi: 10.1016/j.clinph.2012.01.013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Soon, C. S., Brass, M., Heinze, H. J., and Haynes, J. D. (2008). Unconscious determinants of free decisions in the human brain. Nat. Neurosci. 11, 543–545. doi: 10.1038/nn.2112

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Thayer, J. F., and Lane, R. D. (2009). Claude Bernard and the heart-brain connection: further elaboration of a model of neurovisceral integration. Neurosci. Biobehav. Rev. 33, 81–88. doi: 10.1016/j.neubiorev.2008.08.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tzagarakis, C., Ince, N. F., Leuthold, A. C., and Pellizzer, G. (2010). Beta-band activity during motor planning reflects response uncertainty. J. Neurosci. 30, 11270–11277. doi: 10.1523/JNEUROSCI.6026-09.2010

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Waldert, S., Preissl, H., Demandt, E., Braun, C., Birbaumer, N., Aertsen, A., et al. (2008). Hand movement direction decoded from MEG and EEG. J. Neurosci. 28, 1000–1008. doi: 10.1523/JNEUROSCI.5171-07.2008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wang, J., Mordkoff, J. T., and Sainburg, R. L. (2010). Visuomotor learning generalizes between bilateral and unilateral conditions despite varying degrees of bilateral interference. J. Neurophysiol. 104, 2913–2921. doi: 10.1152/jn.00225.2010

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zhang, Y., Chen, Y., Bressler, S. L., and Ding, M. (2008). Response preparation and inhibition: the role of the cortical sensorimotor beta rhythm. Neuroscience 156, 238–246. doi: 10.1016/j.neuroscience.2008.06.061

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text