ORIGINAL RESEARCH article

Front. Hum. Neurosci., 06 July 2017

Sec. Brain Imaging and Stimulation

Volume 11 - 2017 | https://doi.org/10.3389/fnhum.2017.00356

Resting-State Fluctuations of EEG Sensorimotor Rhythm Reflect BOLD Activities in the Pericentral Areas: A Simultaneous EEG-fMRI Study

  • 1. School of Fundamental Science and Technology, Graduate School of Keio University Kanagawa, Japan

  • 2. The Japan Society for the Promotion of Science Tokyo, Japan

  • 3. Department of Biosciences and Informatics, Faculty of Science and Technology, Keio University Kanagawa, Japan

  • 4. Department of Rehabilitation Medicine, Keio University School of Medicine Tokyo, Japan

  • 5. Saiseikai Kanagawa-ken Hospital Kanagawa, Japan

  • 6. Integrative Brain Imaging Center, National Center of Neurology and Psychiatry Tokyo, Japan

  • 7. Japan Science and Technology Agency, Precursory Research for Embryonic Science and Technology Saitama, Japan

  • 8. Keio Institute of Pure and Applied Sciences Kanagawa, Japan

Abstract

Blockade of the scalp electroencephalographic (EEG) sensorimotor rhythm (SMR) is a well-known phenomenon following attempted or executed motor functions. Such a frequency-specific power attenuation of the SMR occurs in the alpha and beta frequency bands and is spatially registered at primary somatosensory and motor cortices. Here, we hypothesized that resting-state fluctuations of the SMR in the alpha and beta frequency bands also covary with resting-state sensorimotor cortical activity, without involving task-related neural dynamics. The present study employed functional magnetic resonance imaging (fMRI) to investigate the neural regions whose activities were correlated with the simultaneously recorded SMR power fluctuations. The SMR power fluctuations were convolved with a canonical hemodynamic response function and correlated with blood-oxygen-level dependent (BOLD) signals obtained from the entire brain. Our findings show that the alpha and beta power components of the SMR correlate with activities of the pericentral area. Furthermore, brain regions with correlations between BOLD signals and the alpha-band SMR fluctuations were located posterior to those with correlations between BOLD signals and the beta-band SMR. These results are consistent with those of event-related studies of SMR modulation induced by sensory input or motor output. Our findings may help to understand the role of the sensorimotor cortex activity in contributing to the amplitude modulation of SMR during the resting state. This knowledge may be applied to the diagnosis of pathological conditions in the pericentral areas or the refinement of brain–computer interfaces using SMR in the future.

Introduction

Since Berger’s first electroencephalogram (EEG) recordings from the human scalp in the late 1920s, a number of studies have led to new insights into the function and mechanisms of intrinsic oscillations underlying brain activities (Niedermeyer and da Silva, 2005). One of the characteristic EEG rhythms is the sensorimotor rhythm (SMR), a spontaneous EEG (i.e., in the absence of somatosensory input or motor output) with an arch-shaped waveform in alpha (e.g., 7–11 Hz) and beta (e.g., 12–30 Hz) frequency bands (Kuhlman, 1978; Kozelka and Pedley, 1990; Arroyo et al., 1993; Kuhlman and Allison, 1977; Laufs et al., 2003; Hasan et al., 2016). The SMR signals are recorded from the C3 channel that is located closest to the hand area in sensorimotor area and its amplitude is known to be reduced by desynchronized neural activities in association with motor-related events, such as kinesthetic motor imagery or actual muscle contraction, suggesting a possible electrophysiological sign of sensorimotor excitability (Pfurtscheller et al., 2006). Actually, such SMR event-related desynchronizations (ERDs) are often localized over the pericentral gyri, and the degree of SMR-ERD is associated with cortico-spinal tract excitability (Takemi et al., 2013) as well as with intracortical disinhibition of the primary motor cortex.

Modulations of the SMR can be segregated into two physiologically different components. The alpha frequency band is located slightly posterior to the beta (Pfurtscheller and Lopes da Silva, 1999; Ritter et al., 2009) and shows prominent reaction to somatosensory events (Babiloni et al., 2008; van Ede et al., 2014). On the other hand, the amplitude of the beta frequency band EEG is attenuated, but the phase synchronization is maintained between the beta frequency band EEG and the spinal motoneuronal pool activities during tonic motor contraction tasks (Conway et al., 1995; Leocani et al., 1997; Bayraktaroglu et al., 2013). Furthermore, in the post-movement period, the amplitude of the beta-band consistently returns to and exceeds pre-movement levels via synchronization (Pfurtscheller et al., 1996a; Bauer et al., 2006; Jurkiewicz et al., 2006; Parkes et al., 2006). These studies suggest that the beta component of the SMR is more motor-related.

Does SMR also show some spontaneous fluctuations during the resting state, and, if so, is the resting-state SMR fluctuation correlated with blood-oxygen-level dependent (BOLD) signal changes in the sensorimotor areas? Spatial localization of alpha and beta components of the resting-state EEG-SMR has yet to be directly confirmed, but several lines of collateral evidence have supported such correlations. For instance, studies on sensory stimulation have shown that the prestimulus amplitude of the alpha-SMR has a significant impact on sensory stimulus detection (Linkenkaer-Hansen et al., 2004; Zhang and Ding, 2009). A few studies using single-pulse transcranial magnetic stimulation (TMS) over the primary motor cortex have shown an association between motor cortex activity or corticospinal tract excitability and the amplitude of the spontaneous beta-SMR during the resting state (Mäki and Ilmoniemi, 2010; Berger et al., 2014). Spontaneous fluctuations of the resting-state EEG may be much subtler than the motor task-related variations. However, similar to the fact that the level of the SMR synchronization reflects sensorimotor cortex activity during the course of a movement task (Pfurtscheller et al., 2000), the characteristics of the SMR may also be indicative of sensorimotor cortex activity during the resting state.

Here, based upon previous knowledge suggesting an association between spontaneous SMR fluctuation during the resting state and activity of the pericentral brain regions, we addressed the following hypotheses: (1) the spontaneous SMR power modulations were correlated with a surrogate marker of brain activity that covaries with resting-state sensorimotor cortical activity as measured by BOLD-fMRI in the pericentral area, and (2) the area correlated with alpha-band SMR was located posterior to the area correlated with beta-band SMR. To test such hypotheses, we employed EEG-fMRI simultaneous recording to identify the relationship between SMR modulations and whole-brain activity during the resting state.

Materials and Methods

Subjects

Nineteen healthy subjects (13 men and 6 women; aged 21–25 years) participated in this study. None had any sleep, medical, or psychiatric disorders. The purpose and experimental procedure were explained to the subjects, and all subjects gave informed written consent. The study was conducted in accordance with the Declaration of Helsinki and was approved by Keio University Faculty of Science and Technology Bioethics Committee and Saiseikai Kanagawa-ken Hospital Ethics Committee.

Data Acquisition

The subjects were asked to lie still on a scanner bed in the dark for 10 min with their eyes opened and fixed on a small black cross. EEG signals were recorded with an MR-compatible amplifier (BrainAmp MR plus, Brain Products GmbH, Germany) and an EEG electrode cap (BrainCap MR, Brain Products GmbH, Germany) providing 63 EEG channels and 1 electrocardiogram (ECG) channel. Their electrodes were placed according to the modified International 10–10 system (Jurcak et al., 2007). The ground electrode was placed at the AFz, the reference electrode at the FCz, and the impedance of all electrodes was kept lower than 5 kΩ throughout the experiment. The ECG electrode was placed on the back of subjects to obtain the electrocardiographic data and to subsequently correct ballistocardiographic artifacts. Raw record was sampled at 5 kHz with a bandpass filter between 0.1 and 2500 Hz using the Brain Vision Recorder, (Brain Products GmbH, Germany). The amplifier system was set beside the subject’s head within the scanner during fMRI scanning.

Functional images were acquired on a 1.5-T MR scanner (Signa Excite, GE Medical Systems, United Kingdom) and whole-brain, T2-weighted BOLD data were acquired using an echo planar image sequence (repetition time TR = 3000 ms; echo time TE = 40 ms; flip angle = 70 degrees; voxel size = 3.75 mm × 3.75 mm × 4.0 mm; 30 axial slices with no gap, 200 scans). To achieve phase synchronization clocks for digital sampling between the MRI data and the EEG system, the EEG system clock was synchronized with a SyncBox device (Brain Products GmbH, Germany) and the MRI scanner’s 10 MHz master synthesizer. The scanner also delivered each TR trigger signal that marked the onset time of every fMRI volume acquisition. These markers were used for fMRI scanning artifact correction of the EEG data.

Data Analyses

An outline of the analyses is shown in Figure 1. The several-seconds scale of modulation in the SMR resting-state amplitude was first extracted from EEGs of healthy volunteers, and the SMR power time-series was convolved with the canonical hemodynamic response function (HRF) to model the BOLD signals co-varying with the SMR modulation. A voxel-wise correlation analysis was then conducted in the whole-brain to identify the cortical and subcortical activities with which the modulation of SMR amplitude was associated.

FIGURE 1

Electroencephalogram data were processed by Brain Vision Analyzer 2.0 (Brain Products GmbH, Germany) and the MR gradient-artifacts and ballistocardiogram in the EEG signals were corrected using the average template subtraction method (Allen et al., 2000). First, to create the gradient-artifact template, EEG data were averaged over three volumes (9 s) with the reference to each TR trigger signal. The template was computed for each TR interval using the sliding window. The gradient-artifact template was subsequently subtracted from the original EEG data. Second, the ballistocardiographic artifacts were removed in a similar way. The template in this case was constructed by averaging of the EEG data with the reference to the peaks of the R-waves in the ECG signal. After the gradient-artifacts and ballistocardiographic corrections, the EEG data was bandpass filtered (3–70 Hz with 48 dB/octave phase shift-free Butterworth) with a notch (50 Hz for avoiding power line contamination), and the sampling rate of the data was decreased to 250 Hz.

It is known that the EEG signal obtained around the pericentral area contains cognitive components, such as activity related to the “mirror neuron systems,” in addition to the sensorimotor idling activity component (Yin et al., 2016). Because we focused on the sensorimotor processing component, we obtained SMR signals from the C3 channel that is located closest to the pericentral area. Furthermore, the re-reference to the common average reference should remove global-signal to emphasize the relevant localized activities around the pericentral area (McFarland et al., 1997). The power time courses of EEG data in the C3 channel were calculated by a Fast Fourier Transform (FFT) with a 1-s moving time window and no overlap. We first selected four different frequency bands [7–11 Hz (alpha), 12–16 Hz (beta1), 17–21 Hz (beta2), 22–30 Hz (beta3)] as candidates for the sensorimotor-related features (Arroyo et al., 1993; Rangaswamy et al., 2002, 2004; Ritter et al., 2009), and the power time courses of SMR to represent the SMR amplitude fluctuations were obtained by calculating the band power in each time window by averaging the magnitudes of the Fourier coefficients in the related frequency bins. To confirm that these frequency components were localized over the sensorimotor cortex, simple linear regression coefficients at a zero time lag were calculated between the C3 and the other channels. The coefficients were assigned to each channel position where the correlation with C3 was calculated, and a topographic map was obtained using the nearest neighbor interpolation with the MATLAB 4 griddata method. After calculating the regression coefficients for each individual subject, a one-sample t-test was performed for group analysis.

The fMRI data were pre-processed with Statistical Parametric Mapping (SPM12, Wellcome Department of Imaging Neuroscience1) on MATLAB (MathWorks, Natick, MA, United States). First, the slice-timing correction was applied with the reference slice based on the middle of each TR, and each volume was spatially realigned with the first volume. After the timing and spatial corrections, each volume was normalized according to the Montreal Neurological Institute (MNI) template, and was resampled to 2-mm cubic voxels. The normalized volumes were then smoothed with an isotropic Gaussian kernel [Full width at half maximum (FWHM) = 8 mm], and were 1/128 Hz high-pass filtered.

For the correlation analysis of BOLD signals with SMR power modulation in each frequency band, we followed a previous approach in the simultaneous EEG-fMRI study (Huster et al., 2012). First of all, the SMR power time series for bandwidth was convolved with HRF provided by SPM12 to take into account the hemodynamic delay of the BOLD signal. After the convolution, to conform the temporal resolution of the fMRI data, the SMR power time courses were down-sampled with the data point sampled at the reference of the fMRI slice timing correction; namely, one data point was sampled in the middle of every TR. To identify the brain regions with correlations between BOLD signals and SMR fluctuations, each alpha and beta frequency band of SMR power was used as a predictor for the BOLD signal in the General Linear Model (GLM). Because the alpha- and beta-band SMR power modulations were significantly correlated, a GLM model was created for each frequency of SMR power to avoid multicollinearity. The six realignment parameters (three translation and three rotations) were used as nuisance variables to remove head motion artifacts, and the cerebrospinal fluid (CSF) signal was also used as a nuisance variable to exclude non-neural noise originating from cardiac effects and respiration. The CSF time series were obtained by averaging the signal over all voxels within the Lateral Ventricles Mask, which were anatomically defined by the SPM12 atlas.

In the first-level analysis, we generated regression coefficient contrast images of alpha- and beta-band SMR modulations for each subject. Then, for the second-level, group data analysis, the contrast images were fed into a one-sample t-test. For all data, the fMRI results were subjected to peak- and cluster-extent-based thresholding, the most popular thresholding method for dealing with the issue of multiple comparisons over many voxels (Woo et al., 2014). The initial voxel threshold was set to p < 0.001 uncorrected for multiple comparisons (this level was considered reasonably conservative against a false positive finding). Thus, clusters were considered significant if they survived an extent threshold of p = 0.05 family-wise error (FWE) corrected for multiple comparisons for statistical inference. In each region of sensorimotor cortex correlated with alpha- or beta-power fluctuations of the SMR, we performed statistical evaluation tests to validate that these correlated regions were spatially different using a two-sample Kolmogorov–Smirnov test.

Results

EEG Results

A typical SMR time series of a single subject is shown in Figure 2A, and the power spectrums of the SMRs in all subjects are depicted in Figure 2B. Although EEG signal power varied among the subjects, the averaged power spectrum density showed that the overall frequency structures were shared across the subjects, with the peak-frequencies at 10–12 Hz and 22–26 Hz. The averaged powers in the alpha and beta1 frequency band were 6.8 × 10-5 V2 and 2.9 × 10-5 V2, respectively. The powers of beta2 and beta3 were 6.0 × 10-6 V2 and 7.5 × 10-6 V2, respectively, and were relatively small in comparison with the former two frequency bands. For each subject, the power-modulation-depth ratios were calculated as coefficients of variation; namely, the ratio of the standard deviation (SD) to the mean of time-fluctuation. The inter-subject average and SD of the power-modulation-depth ratios were 85 ± 16% in alpha; 112 ± 20% in beta1; 88 ± 39% in beta2 and 68 ± 19% in beta3 during the resting state. Figure 3A shows the spatial distribution of the selected frequency modulation components from C3. All of the distribution maps showed positive correlations locally around C3 and its contralateral (C4) electrode, with notable localization in the alpha and beta2 components.

FIGURE 2

FIGURE 3

fMRI Results

Figure 3B depicts the statistically significant t-maps of fMRI correlates (BOLD signals) of SMR modulation, which were assessed separately for the alpha and each beta component. The alpha-band SMR modulation was negatively correlated with the activities of the bilateral pericentral cortices (p < 0.001 peak level uncorrected; p < 0.05 voxel level FWE corrected). The t-values and MNI coordinates are shown in Table 1. Similar to the correlation of BOLD signals with the alpha-band SMR modulation, the beta2-band SMR modulation was negatively correlated with the activities of the bilateral pericentral cortices (Table 2). The beta3 SMR modulation was not correlated with any region. The beta1 SMR modulation was negatively correlated with the BOLD signal in the right middle frontal gyrus; this correlation was also seen in the analysis of the alpha- and beta2-band SMR modulations.

Table 1

ModelCorrelation typeBrain regionSideMNI coordinatest-scorep-valueCluster size

xyz
AlphaNegativeTriangular part of the inferior frontal gyrusR5038-65.956.30 × 10-5316
Postcentral gyrusL-52-32525.887.23 × 10-62002
R46-28625.857.61 × 10-61250
Middle frontal gyrusL-4428285.571.36 × 10-5430
Superior parietal lobuleR22-68505.531.50 × 10-5791
L-16-58645.252.72 × 10-5642
Superior frontal gyrusR262625.252.73 × 10-5241
Opercular part of the inferior frontal gyrusR5622165.034.34 × 10-5838
Middle occipital gyrusL-34-84304.541.27 × 10-4202

Brain regions whose activity correlated with the power of the alpha-band SMR modulation.

Brain regions are named based on the SPM12 atlas (Statistical Parametric Mapping, Wellcome Department of Imaging Neuroscience, http://www.fil.ion.ucl.ac.uk/spm) p-value, cluster-level family-wise error rate of 0.05. MNI, Montreal Neurological Institute.

Table 2

ModelCorrelation typeBrain regionSideMNI coordinatest-scorep-valueCluster size

xyz
Beta1NegativeMiddle frontal gyrusR3820306.382.61 × 10-5314
Beta2PositiveThalamusR2-1244.421.65 × 10-4270
L-6-863.718.06 × 10-4270
NegativeSuperior frontal gyrusL-1024586.512.00 × 10-6326
R2254304.728.47 × 10-5202
Postcentral gyrusL-44-28506.332.86 × 10-61262
R52-18525.946.41 × 10-61024
Opercular part of the inferior frontal gyrusR4214246.094.64 × 10-61758
Posterior orbital gyrusR3224-185.887.22 × 10-6313
Superior temporal gyrusL-60-2025.401.98 × 10-5293
Lingual gyrusL-20-68-105.014.50 × 10-5219
R20-66-44.689.34 × 10-5376
Middle occipital gyrusR32-82324.728.52 × 10-5224
Superior frontal gyrus medial partR650184.581.17 × 10-4340
Middle temporal gyrusR60-3004.401.73 × 10-4228

Brain regions whose activity correlated with the power of the SMR beta1- and beta2-bands.

Brain regions are named based on the SPM12 atlas. p-value, cluster-level family-wise error rate of 0.05. MNI, Montreal Neurological Institute.

We further examined a spatial topography of fMRI signal correlated with both the alpha- and beta2-band SMR modulations in the left pericentral area, where the C3 electrode was located. Figure 4 shows the pericentral activity correlated with alpha-band SMR, beta2-SMR, and overlap of both. The peak coordinate of the alpha-band SMR-correlated activity was MNI (x, y, z) = (-52, -32, 52) and that of the beta2-band SMR was MNI (x, y, z) = (-44, -28, 50). Furthermore, Figure 5 displays the distribution of these correlated areas along the anterior–posterior (y) axis. The spatial distributions of alpha-SMR and beta2-SMR correlated activities were different in both the postcentral (Figure 5A) and precentral (Figure 5B) areas (p < 0.001, two-sample Kolmogorov–Smirnov test). More precisely, the alpha-band SMR-correlated BOLD signals were distributed posterior to those correlated with the beta2-SMR.

FIGURE 4

FIGURE 5

Discussion

In the present study, we employed simultaneous EEG-fMRI recordings in healthy human subjects to investigate the whole-brain hemodynamic responses associated with spontaneous SMR fluctuations. We found significant correlations between BOLD signals of activities in the sensorimotor regions and modulation of EEG-SMRs in the resting state. In addition, the regions with negative correlations between BOLD signals and the alpha-band SMR were distributed posterior to those correlated with the beta-band SMR.

SMR Frequency Components in Alpha- and Beta-Bands

Two of the four frequency components, alpha- and beta2-band, showed distributions of significant correlations locally around C3 and its periphery. According to Ritter et al. (2009) study using a blind source separation algorithm, the SMR power modulations of the alpha (∼10 Hz) and beta (∼20 Hz) frequency bands were associated with the motor task. Thus, the amplitude modulation of the alpha- and beta2-bands of the SMR may reflect fluctuations in activity in the sensorimotor regions. Moreover, both the alpha- and beta2-band SMR components were correlated with BOLD signals around the bilateral pericentral areas, that periphery area of around C3 as well as the contralateral C4 area. The homolog C4 of bilateral pericentral areas, which is consistent with the fact that the two cerebral hemispheres have anatomical connections via the corpus callosum and each hemisphere communicates with the contralateral hemisphere by inhibitory and excitatory inputs (Bloom and Hynd, 2005; Palmer et al., 2012).

Correlation between Intrinsic SMR Fluctuations and BOLD Signals

Both the alpha- and beta2-band SMR components were correlated with BOLD signals around the bilateral pericentral areas. Previous fMRI studies also found that in the resting state, BOLD signals in the bilateral sensorimotor cortices were correlated with one another, a phenomenon referred to as the sensorimotor network (Biswal et al., 1995; Beckmann et al., 2005). Furthermore, the alpha or beta synchronizations are interpreted as correlates of “idling” motor cortex neurons (Pfurtscheller et al., 1996b; Ritter et al., 2009). Therefore, the intrinsic alpha- and beta-band SMR modulations might reflect the sensorimotor network in the resting state. In addition, another simultaneous EEG-fMRI study showed that alpha- and beta-band SMR modulations were negatively correlated with the sensorimotor network (Yin et al., 2016). Our results are compatible with these previous findings, suggesting that it is valid to draw correlations between the intrinsic SMR modulations and activity in the sensorimotor cortex.

Functional Properties of Alpha- and Beta-Band SMR

Although the topographic map of the EEG alpha-band SMR correlation was similar to that of the beta2 correlation, it may reflect somatosensory-dominant activities of the cortex. Our fMRI results may indicate functional differences between the alpha- and beta-band components of the SMR. Although both SMR components were correlated with BOLD signals around the pericentral area, the area correlated with the alpha-band SMR was in the parietal region and that correlated with the beta-band SMR was the frontal region. Recent findings suggest that the functional role of alpha oscillations is closely related to the activity of intrinsic cortical networks, indicating that the alpha-band signal occurs in different cortical layers (Klimesch et al., 2007; Palva and Palva, 2007; Ben-Simon et al., 2008; Bazanova and Vernon, 2014; Sigala et al., 2014). Activities in the beta-band, however, are known to be connected with motor functions (Brovelli et al., 2004; Jurkiewicz et al., 2006; Tzagarakis et al., 2015). In addition, in a combined TMS-electromyography (EMG) and concurrent EEG recording study, the Farzan group reported that EEG components correlated with the TMS-induced silence period in EMG, and that the cortical inhibitory processes in humans were related to the power of beta oscillations more locally than to that of the alpha oscillations (Farzan et al., 2013).

Spatial Distribution of Alpha- and Beta-Bands SMR

Previous studies reported the SMR to be most prominent in central scalp regions in the area of the sensorimotor cortex (Pfurtscheller and Neuper, 1997; Birbaumer et al., 2000; Blankertz et al., 2010), while our results depict the spatial differences of correlation areas in intrinsic alpha- and beta-bands of the SMR. In studies of SMR power modulation associated with motor execution or motor imaginary tasks, the attenuated alpha- and beta-band SMR amplitudes were localized around the contralateral sensorimotor area (Crone et al., 1998; Sauseng et al., 2013). Furthermore, the center of ERD of the alpha-band SMR was located more posterior than that of the beta-band SMR in EEG/MEG studies (Pfurtscheller and Lopes da Silva, 1999; Jurkiewicz et al., 2006). The task-related reactive sources of the alpha-band SMR were localized in the postcentral area and those of the beta-band SMR were in the precentral area during sensory input (Cheyne et al., 2003) or motor output (Salmelin et al., 1995). Although it is not obvious that the neural processing of intrinsic SMR modulation during the resting state is equivalent to that of event-related task modulation, our results indicate that the spontaneous modulations of the SMR may be based upon a similar mechanism with that of the SMR modulated by sensory inputs or movements which were already shown by event-related studies to represent pericentral activities; thus the SMR modulation in the resting state is also likely to reflect pericentral activities. Therefore, insight into the relationship between SMR modulation and the pericentral area, which has been examined in event-related studies, has now been extended to the activities of the intrinsic SMR modulation during the resting state. In the present study, neural activity correlated with SMR modulation at rest showed similar spatial localization with the neural activity correlated with movement-related SMR modulation in the literature. Thus, we have provided a testable hypothesis that a similar mechanism may underlie SMR modulation during a motor task and during a resting state. A future study directly testing this hypothesis is warranted.

Study Limitation

No EMG/visual monitoring was employed to rule out potential bodily movements during EEG-fMRI.

Statements

Author contributions

The work presented here was carried out as a collaboration among all authors. ST and JU conceived and designed the research; ST and SS performed experiments; ST analyzed data; ST, TH, and JU interpreted results of experiments; ST prepared figures and drafted the manuscript; all authors edited and revised the manuscript; all authors have read and approved the final version of the manuscript.

Funding

This study was partially supported by the Strategic Research Program for Brain Sciences from the Japan Agency for Medical Research and Development (16dm0107033h005) and Keio Institute of Pure and Applied Sciences (KiPAS) research program.

Acknowledgments

We thank Sayoko Ishii, Kumi Nanjo, and Sawako Ohtaki for their technical support.

Conflict of interest

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

References

  • 1

    AllenP. J.JosephsO.TurnerR. (2000). A method for removing imaging artifact from continuous EEG recorded during functional MRI.Neuroimage12230239. 10.1006/nimg.2000.0599

  • 2

    ArroyoS.LesserR. P.GordonB.UematsuS.JacksonD.WebberR. (1993). Functional significance of the mu rhythm of human cortex: an electrophysiologic study with subdural electrodes.Electroencephalogr. Clin. Neurophysiol.877687. 10.1016/0013-4694(93)90114-B

  • 3

    BabiloniC.CapotostoP.BrancucciA.Del PercioC.PetriniL.ButtiglioneM.et al (2008). Cortical alpha rhythms are related to the anticipation of sensorimotor interaction between painful stimuli and movements: a high-resolution EEG study.J. Pain9902911. 10.1016/j.jpain.2008.05.007

  • 4

    BauerM.OostenveldR.PeetersM.FriesP. (2006). Tactile spatial attention enhances gamma-band activity in somatosensory cortex and reduces low-frequency activity in parieto-occipital areas.J. Neurosci.26490501. 10.1523/JNEUROSCI.5228-04.2006

  • 5

    BayraktarogluZ.von Carlowitz-GhoriK.CurioG.NikulinV. V. (2013). It is not all about phase: amplitude dynamics in corticomuscular interactions.Neuroimage64496504. 10.1016/j.neuroimage.2012.08.069

  • 6

    BazanovaO. M.VernonD. (2014). Interpreting EEG alpha activity.Neurosci. Biobehav. Rev.4494110. 10.1016/j.neubiorev.2013.05.007

  • 7

    BeckmannC. F.DeLucaM.DevlinJ. T.SmithS. M. (2005). Investigations into resting-state connectivity using independent component analysis.Philos. Trans. R. Soc. Lond. B Biol. Sci.36010011013. 10.1098/rstb.2005.1634

  • 8

    Ben-SimonE.PodlipskyI.ArieliA.ZhdanovA.HendlerT. (2008). Never resting brain: simultaneous representation of two alpha related processes in humans.PLoS ONE3:e3984. 10.1371/journal.pone.0003984

  • 9

    BergerB.MinarikT.LiuzziG.HummelF. C.SausengP. (2014). EEG oscillatory phase-dependent markers of corticospinal excitability in the resting brain.BioMed Res. Int.2014:936096. 10.1155/2014/936096

  • 10

    BirbaumerN.KublerA.GhanayimN.HinterbergerT.PerelmouterJ.KaiserJ.et al (2000). The thought translation device (TTD) for completely paralyzed patients.IEEE Trans. Rehabil. Eng.8190193. 10.1109/86.847812

  • 11

    BiswalB.YetkinF. Z.HaughtonV. M.HydeJ. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI.Magn. Reson. Med.34537541. 10.1002/mrm.1910340409

  • 12

    BlankertzB.TangermannM.VidaurreC.FazliS.SannelliC.HaufeS.et al (2010). The Berlin brain–computer interface: non-medical uses of BCI technology.Front. Neurosci.4:198. 10.3389/fnins.2010.00198

  • 13

    BloomJ. S.HyndG. W. (2005). The role of the corpus callosum in interhemispheric transfer of information: excitation or inhibition?Neuropsychol. Rev.155971. 10.1007/s11065-005-6252-y

  • 14

    BrovelliA.DingM.LedbergA.ChenY.NakamuraR.BresslerS. L. (2004). Beta oscillations in a large-scale sensorimotor cortical network: directional influences revealed by Granger causality.Proc. Natl. Acad. Sci. U.S.A.10198499854. 10.1073/pnas.0308538101

  • 15

    CheyneD.GaetzW.GarneroL.LachauxJ.-P.DucorpsA.SchwartzD.et al (2003). Neuromagnetic imaging of cortical oscillations accompanying tactile stimulation.Cogn. Brain Res.17599611. 10.1016/S0926-6410(03)00173-3

  • 16

    ConwayB. A.HallidayD. M.FarmerS. F.ShahaniU.MaasP.WeirA. I.et al (1995). Synchronization between motor cortex and spinal motoneuronal pool during the performance of a maintained motor task in man.J. Physiol.489917924. 10.1113/jphysiol.1995.sp021104

  • 17

    CroneN. E.MigliorettiD. L.GordonB.SierackiJ. M.WilsonM. T.UematsuS.et al (1998). Functional mapping of human sensorimotor cortex with electrocorticographic spectral analysis. I. Alpha and beta event-related desynchronization.Brain121(Pt 12), 22712299. 10.1093/brain/121.12.2271

  • 18

    FarzanF.BarrM. S.HoppenbrouwersS. S.FitzgeraldP. B.ChenR.Pascual-LeoneA.et al (2013). The EEG correlates of the TMS induced EMG silent period in humans.Neuroimage83120134. 10.1016/j.neuroimage.2013.06.059

  • 19

    HasanM. A.FraserM.ConwayB. A.AllanD. B.VučkovićA. (2016). Reversed cortical over-activity during movement imagination following neurofeedback treatment for central neuropathic pain.Clin. Neurophysiol.12731183127. 10.1016/j.clinph.2016.06.012

  • 20

    HusterR. J.DebenerS.EicheleT.HerrmannC. S. (2012). Methods for simultaneous EEG-fMRI: an introductory review.J. Neurosci.3260536060. 10.1523/JNEUROSCI.0447-12.2012

  • 21

    JurcakV.TsuzukiD.DanI. (2007). 10/20, 10/10, and 10/5 systems revisited: their validity as relative head-surface-based positioning systems.Neuroimage3416001611. 10.1016/j.neuroimage.2006.09.024

  • 22

    JurkiewiczM. T.GaetzW. C.BostanA. C.CheyneD. (2006). Post-movement beta rebound is generated in motor cortex: evidence from neuromagnetic recordings.Neuroimage3212811289. 10.1016/j.neuroimage.2006.06.005

  • 23

    KlimeschW.SausengP.HanslmayrS. (2007). EEG alpha oscillations: the inhibition–timing hypothesis.Brain Res. Rev.536388. 10.1016/j.brainresrev.2006.06.003

  • 24

    KozelkaJ. W.PedleyT. A. (1990). Beta and mu rhythms.J. Clin. Neurophysiol.7191208. 10.1097/00004691-199004000-00004

  • 25

    KuhlmanW. N. (1978). Functional topography of the human mu rhythm.Electroencephalogr. Clin. Neurophysiol.448393. 10.1016/0013-4694(78)90107-4

  • 26

    KuhlmanW. N.AllisonT. (1977). EEG feedback training in the treatment of epilepsy: some questions and some answers.Pavlov. J. Biol. Sci. Off. J. Pavlov.12112122. 10.1007/BF03004498

  • 27

    LaufsH.KrakowK.SterzerP.EgerE.BeyerleA.Salek-HaddadiA.et al (2003). Electroencephalographic signatures of attentional and cognitive default modes in spontaneous brain activity fluctuations at rest.Proc. Natl. Acad. Sci. U.S.A.1001105311058. 10.1073/pnas.1831638100

  • 28

    LeocaniL.ToroC.ManganottiP.ZhuangP.HallettM. (1997). Event-related coherence and event-related desynchronization/synchronization in the 10 Hz and 20 Hz EEG during self-paced movements.Electroencephalogr. Clin. Neurophysiol.104199206. 10.1016/S0168-5597(96)96051-7

  • 29

    Linkenkaer-HansenK.NikulinV. V.PalvaS.IlmoniemiR. J.PalvaJ. M. (2004). Prestimulus oscillations enhance psychophysical performance in humans.J. Neurosci.241018610190. 10.1523/JNEUROSCI.2584-04.2004

  • 30

    MäkiH.IlmoniemiR. J. (2010). EEG oscillations and magnetically evoked motor potentials reflect motor system excitability in overlapping neuronal populations.Clin. Neurophysiol.121492501. 10.1016/j.clinph.2009.11.078

  • 31

    McFarlandD. J.McCaneL. M.DavidS. V.WolpawJ. R. (1997). Spatial filter selection for EEG-based communication.Electroencephalogr. Clin. Neurophysiol.103386394. 10.1016/S0013-4694(97)00022-2

  • 32

    NiedermeyerE.da SilvaF. H. L. (2005). Electroencephalography: Basic Principles, Clinical Applications, and Related Fields.Philadelphia, PA: Lippincott Williams & Wilkins.

  • 33

    PalmerL. M.SchulzJ. M.MurphyS. C.LedergerberD.MurayamaM.LarkumM. E. (2012). The cellular basis of GABAB-mediated interhemispheric inhibition.Science335989993. 10.1126/science.1217276

  • 34

    PalvaS.PalvaJ. M. (2007). New vistas for α-frequency band oscillations.Trends Neurosci.30150158. 10.1016/j.tins.2007.02.001

  • 35

    ParkesL. M.BastiaansenM. C. M.NorrisD. G. (2006). Combining EEG and fMRI to investigate the post-movement beta rebound.Neuroimage29685696. 10.1016/j.neuroimage.2005.08.018

  • 36

    PfurtschellerG.BrunnerC.SchlöglA.Lopes da SilvaF. H. (2006). Mu rhythm (de)synchronization and EEG single-trial classification of different motor imagery tasks.Neuroimage31153159. 10.1016/j.neuroimage.2005.12.003

  • 37

    PfurtschellerG.Lopes da SilvaF. H. (1999). Event-related EEG/MEG synchronization and desynchronization: basic principles.Clin. Neurophysiol.11018421857. 10.1016/S1388-2457(99)00141-8

  • 38

    PfurtschellerG.NeuperC. (1997). Motor imagery activates primary sensorimotor area in humans.Neurosci. Lett.2396568. 10.1016/S0304-3940(97)00889-6

  • 39

    PfurtschellerG.NeuperC.KrauszG. (2000). Functional dissociation of lower and upper frequency mu rhythms in relation to voluntary limb movement.Clin. Neurophysiol.11118731879. 10.1016/S1388-2457(00)00428-4

  • 40

    PfurtschellerG.StancákA.Jr.NeuperC. (1996a). Post-movement beta synchronization. A correlate of an idling motor area?Electroencephalogr. Clin. Neurophysiol.98281293. 10.1016/0013-4694(95)00258-8

  • 41

    PfurtschellerG.StancákA.Jr.NeuperC. (1996b). Event-related synchronization (ERS) in the alpha band — an electrophysiological correlate of cortical idling: a review.Int. J. Psychophysiol.243946. 10.1016/S0167-8760(96)00066-9

  • 42

    RangaswamyM.PorjeszB.ChorlianD. B.WangK.JonesK. A.BauerL. O.et al (2002). Beta power in the EEG of alcoholics.Biol. Psychiatry52831842. 10.1016/S0006-3223(02)01362-8

  • 43

    RangaswamyM.PorjeszB.ChorlianD. B.WangK.JonesK. A.KupermanS.et al (2004). Resting EEG in offspring of male alcoholics: beta frequencies.Int. J. Psychophysiol.51239251. 10.1016/j.ijpsycho.2003.09.003

  • 44

    RitterP.MoosmannM.VillringerA. (2009). Rolandic alpha and beta EEG rhythms’ strengths are inversely related to fMRI-BOLD signal in primary somatosensory and motor cortex.Hum. Brain Mapp.3011681187. 10.1002/hbm.20585

  • 45

    SalmelinR.HámáaláinenM.KajolaM.HariR. (1995). Functional segregation of movement-related rhythmic activity in the human brain.Neuroimage2237243. 10.1006/nimg.1995.1031

  • 46

    SausengP.GerloffC.HummelF. C. (2013). Two brakes are better than one: the neural bases of inhibitory control of motor memory traces.Neuroimage655258. 10.1016/j.neuroimage.2012.09.048

  • 47

    SigalaR.HaufeS.RoyD.DinseH. R.RitterP. (2014). The role of alpha-rhythm states in perceptual learning: insights from experiments and computational models.Front. Comput. Neurosci.8:36. 10.3389/fncom.2014.00036

  • 48

    TakemiM.MasakadoY.LiuM.UshibaJ. (2013). Event-related desynchronization reflects downregulation of intracortical inhibition in human primary motor cortex.J. Neurophysiol.11011581166. 10.1152/jn.01092.2012

  • 49

    TzagarakisC.WestS.PellizzerG. (2015). Brain oscillatory activity during motor preparation: effect of directional uncertainty on beta, but not alpha, frequency band.Front. Neurosci.9:246. 10.3389/fnins.2015.00246

  • 50

    van EdeF.SzebényiS.MarisE. (2014). Attentional modulations of somatosensory alpha, beta and gamma oscillations dissociate between anticipation and stimulus processing.Neuroimage97134141. 10.1016/j.neuroimage.2014.04.047

  • 51

    WooC.-W.KrishnanA.WagerT. D. (2014). Cluster-extent based thresholding in fMRI analyses: pitfalls and recommendations.Neuroimage91412419. 10.1016/j.neuroimage.2013.12.058

  • 52

    YinS.LiuY.DingM. (2016). Amplitude of sensorimotor mu rhythm is correlated with BOLD from multiple brain regions: a simultaneous EEG-fMRI study.Front. Hum. Neurosci.10:364. 10.3389/fnhum.2016.00364

  • 53

    ZhangY.DingM. (2009). Detection of a weak somatosensory stimulus: role of the prestimulus mu rhythm and its top–down modulation.J. Cogn. Neurosci.22307322. 10.1162/jocn.2009.21247

Summary

Keywords

resting state, sensorimotor rhythm, alpha-band, beta-band, spontaneous fluctuation, pericentral area, simultaneous EEG-fMRI

Citation

Tsuchimoto S, Shibusawa S, Mizuguchi N, Kato K, Ebata H, Liu M, Hanakawa T and Ushiba J (2017) Resting-State Fluctuations of EEG Sensorimotor Rhythm Reflect BOLD Activities in the Pericentral Areas: A Simultaneous EEG-fMRI Study. Front. Hum. Neurosci. 11:356. doi: 10.3389/fnhum.2017.00356

Received

23 January 2017

Accepted

21 June 2017

Published

06 July 2017

Volume

11 - 2017

Edited by

Tetsuo Kida, National Institute for Physiological Sciences (NIPS), Japan

Reviewed by

Tamer Demiralp, Istanbul University, Turkey; Anna Barnes, University College London, United Kingdom

Updates

Copyright

*Correspondence: Junichi Ushiba,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics