Abstract
Moving in synchrony with an auditory rhythm requires predictive action based on neurodynamic representation of temporal information. Although it is known that a regular auditory rhythm can facilitate rhythmic movement, the neural mechanisms underlying this phenomenon remain poorly understood. In this experiment using human magnetoencephalography, 12 young healthy adults listened passively to an isochronous auditory rhythm without producing rhythmic movement. We hypothesized that the dynamics of neuromagnetic beta-band oscillations (∼20 Hz)—which are known to reflect changes in an active status of sensorimotor functions—would show modulations in both power and phase-coherence related to the rate of the auditory rhythm across both auditory and motor systems. Despite the absence of an intention to move, modulation of beta amplitude as well as changes in cortico-cortical coherence followed the tempo of sound stimulation in auditory cortices and motor-related areas including the sensorimotor cortex, inferior-frontal gyrus, supplementary motor area, and the cerebellum. The time course of beta decrease after stimulus onset was consistent regardless of the rate or regularity of the stimulus, but the time course of the following beta rebound depended on the stimulus rate only in the regular stimulus conditions such that the beta amplitude reached its maximum just before the occurrence of the next sound. Our results suggest that the time course of beta modulation provides a mechanism for maintaining predictive timing, that beta oscillations reflect functional coordination between auditory and motor systems, and that coherence in beta oscillations dynamically configure the sensorimotor networks for auditory-motor coupling.
Introduction
Neural representation of temporal information is essential for performing rhythmic movements. A privileged connection between auditory and motor systems has been shown in entraining temporal information and using it for motor planning when subjects consciously perform a repetitive movement task in synchrony with external stimuli (for review, see Repp, 2005). Movement timing with a metronome beat becomes less variable and is increasingly based on prediction especially around the rate of 1–3 Hz (Toma et al., 2002), showing clear facilitation effects of the sound. The 1–3 Hz frequency range coincides with the optimal range of musical tempo in auditory perception (Large, 2008). Neural processing for timing in this range also appears to involve both basal ganglia and cerebellum (Ivry and Spencer, 2004; Mauk and Buonomano, 2004). However, the neural mechanisms underlying auditory-motor facilitation are poorly understood. In functional magnetic resonance imaging (fMRI) studies, auditory-guided or self-paced repetitive motions (Rao et al., 1997; Jäncke et al., 2000; Chen et al., 2008) and listening to and encoding auditory rhythms internally (Grahn and Brett, 2007; Chen et al., 2008) similarly activated primary sensorimotor cortex (SM1), premotor cortex (PMC), supplementary motor area (SMA), striatum, and cerebellum. Thus, auditory processing for timing may rely on a distributed sensorimotor network, which would enable auditory-motor facilitation.
Here, we hypothesized that beta-band activity (∼20 Hz) may play an important role in functional auditory-motor communication even when no movement is performed. Beta oscillations in central and peripheral sensorimotor systems (Murthy and Fetz, 1992; Sanes and Donoghue, 1993; Baker et al., 1997; Courtemanche et al., 2002, 2003) show dynamic modulation representing active states of the system (Engel and Fries, 2010). In human electroencephalography (EEG) and magnetoencephalography (MEG), beta amplitude decreases before and during a movement, followed by a recovery after task completion (Pfurtscheller, 1981; Salmelin et al., 1995). In our previous work investigating auditory processing of isochronous metronome beats without movement tasks, the signal power of the induced beta oscillations in bilateral auditory cortices was periodically modulated at the rate of the metronome (Snyder and Large, 2005; Fujioka et al., 2009). In a similar manner, beta activity in SM1 and SMA was modulated during repetitive finger motion (Gerloff et al., 1998; Toma et al., 2002). Beta phase coherence during auditory-paced motion has been found as another index of functional connectivity between primary auditory cortex (A1), SM1, SMA, cerebellum, thalamus, PMC, and posterior parietal cortex (Pollok et al., 2005). However, the contribution of auditory input to this network and the dynamics of linkages across regions have not been investigated.
Thus, the current study investigated the dynamics of beta activity to examine auditory-motor communication. Specifically, we tested the hypotheses that (1) periodic modulation of induced beta activities would anticipate acoustic events, providing a predictive temporal representation of stimulus rate, and that both the beta (2) power comodulation and (3) phase-coherence would span auditory- and motor-related brain areas, providing evidence for long-range communication (Kopell et al., 2000).
Materials and Methods
Participants.
Twelve healthy right-handed young adults (age, 22–33; seven females) participated in the study. None reported any history of neurological or psychological disorders. Normal hearing between 250 and 4000 Hz was confirmed with clinical audiometry. The Research Ethics Board at Baycrest Centre approved the study, and participants signed consent forms after written information about the study had been provided.
Stimuli.
Auditory stimuli were 262 Hz pure tones with 10 ms rise and fall times and 40 ms steady-state duration. In the first three conditions, sounds were rhythmically presented with onset-to-onset intervals of 390, 585, and 780 ms, equivalent to stimulus rates of ∼2.5, 1.7, and 1.3 Hz. In the fourth condition, the stimuli were presented at random onset-to-onset intervals chosen uniformly between 390 and 780 ms. Sequences of tones were presented in short trains of 7 s duration, containing 18, 12, and 9 tones for the regular conditions, respectively, and ∼12 tones for the irregular condition. Trains were separated by 3 s of silence. Stimulus trains of the same type were presented in blocks of ∼8 min duration. The total number of stimuli per block was 900, 600, 450, and 600, respectively. All conditions were repeated twice with the order of the blocks randomized across participants. The stimulus intensity was set to 60 dB above individual sensation threshold, which was measured immediately before each MEG recording. All stimulus sequences were presented under control of STIM software (NeuroScan). The sound was delivered binaurally through Etymotic ER3A transducers connected with 1.5 m of matched plastic tubing and foam insert earphones to the participant's ears.
MEG recording.
Neuromagnetic activity was recorded with a 151 channel whole-head axial gradiometer type MEG system (OMEGA; VSM MedTech) in a quiet magnetically shielded room. The MEG was recorded continuously with a sampling rate of 1250 Hz after 400 Hz low-pass filtering. Participants were comfortably seated in an upright position with their head resting in the helmet-shaped MEG device and their arms on the arm rests of the chair to avoid any body movements. They were instructed to pay no particular attention to the sounds while watching a subtitled movie of their choice. This is a common practice to maintain participants in an alert state over the time course of a long EEG/MEG experiment investigating automatic sensory functions with passive stimulus experience. Since the participants did not perform any tasks throughout the recording session, there was no explicit association between the auditory stimuli and actions. The compliance of the participants was verified by video monitoring with the visual field covering the arms, hands, and legs to ensure they did not move any body parts in a visible manner. The position of the head relative to the MEG sensor was measured before and after each recording block. Only blocks with head movement <5 mm were accepted for the further data analysis.
Data analysis.
We used two types of source analysis. First, to examine the behavior of auditory induced beta oscillations in bilateral auditory cortices, we obtained a dipole source model and examined the time series of auditory evoked responses and its time–frequency decomposition representing event-related changes in oscillatory activity. Second, to examine functional communication between auditory and sensorimotor systems, we applied a model-free source analysis using a spatial filter based on a MEG beamformer to investigate beta activity across the whole brain. For these data, we used two measures of functional communication based on beta oscillations: first, the comodulation of the magnitudes of oscillation across the whole brain; and second, the temporal changes in coherence specifically linked to the auditory and sensorimotor cortices, and supplementary motor areas.
Beta oscillations in auditory cortices based on dipole source analysis.
For all analyses, we segmented the MEG data into 1900 ms stimulus-related epochs beginning 800 ms before stimulus onsets. Eye blink artifacts were eliminated using principal component analysis (PCA) in which components exceeding the threshold of 1.5 pT at any MEG sensor were subtracted. From time domain averaging of the data in each experimental block, we obtained the auditory evoked responses. Visual inspection of the auditory evoked response revealed that there was no baseline shift before the stimulus onset. This means that there was no sensorimotor evoked response prior to the sound onset, thus indicating the absence of any synchronized movements (see Fig. 1E). The magnetic field maps showed exclusively dipolar patterns of the auditory evoked response, which were modeled with equivalent current dipoles in left and right auditory cortices. Specifically, the locations and orientations of the dipoles were estimated for the P1m peak around 80 ms after stimulus onset because at the fast stimulus rates in our experiment (>1/s) the P1m was the most dominant peak whereas the N1m response was strongly reduced (Näätänen and Picton, 1987) (see Fig. 1E). Based on individual dipole models, the source activities of single trials in auditory cortices were extracted from all epochs using signal space projection.
Time–frequency decompositions were calculated for each subject on each single-trial source waveform in left and right auditory cortices and for each stimulus condition using a Morlet wavelet transform (Bertrand et al., 1994). The wavelet was designed such that the half-maximum width was equal to 2 periods of the center frequency (e.g., 100 ms at 20 Hz). The wavelet transform was calculated for each time point between −800 and 1100 ms with 128 frequency bins between 6 and 100 Hz. The further analysis eliminated both ends of the data segment by 200 ms to avoid aliasing caused by the time–frequency analysis. The mean signal power was calculated at each frequency bin as the mean variance of wavelet coefficients across trials. Normalizing this to its mean value across the whole epoch for each frequency resulted in relative signal power changes expressed in percentage and was visualized as function of time and frequency in color-coded maps of event-related synchronization (ERS) and desynchronization (ERD) (Pfurtscheller and Lopes da Silva, 1999). Time courses of beta power change were calculated as the mean across the frequencies from 13 to 25 Hz. To describe the onset and offset of beta modulation, we measured the latency of the maximum beta decrement and the latencies of the midpoints of the falling and rising slopes of beta power change for each condition. Bootstrap resampling was used to estimate the 95% confidence intervals of the group mean latencies. For this procedure, we treated the datasets as independent samples because the noise was independent from trial to trial.
Beta oscillations across brain based on beamformer source analysis.
For calculating temporally correlated beta power modulation across the whole brain, we used a two-step procedure. First, a spatial filter algorithm, a beamformer called synthetic aperture magnetometry (SAM), was applied to the magnetic field data to calculate the time series of local source activity at 8 × 8 × 8 mm volume elements covering the whole brain. The SAM approach uses a linearly constrained minimum variance beamformer algorithm (Van Veen et al., 1997; Robinson and Vrba, 1999), normalizes source power across the whole cortical volume (Robinson, 2004), and is capable of revealing deep brain sources (Vrba and Robinson, 2001; Vrba, 2002). SAM source analysis has been successfully applied for identifying activities in auditory (Ross et al., 2009) and sensorimotor cortices (Jurkiewicz et al., 2006), and deep sources such as hippocampus (Riggs et al., 2009) and fusiform gyrus and amygdala (Cornwell et al., 2008). In examining beta oscillation here, the beamformer was constructed specifically based on the covariances within the magnetic field data, for all epochs of 15–25 Hz bandpass filtered data from −200 to 1000 ms relative to stimulus onset. Also to construct the SAM spatial filter, we used a template brain magnetic resonance image (MRI) (Holmes et al., 1998), which was transformed to Talairach standard coordinates (positive axes toward anterior, right, and superior direction) using the AFNI software package (Cox, 1996). It has been shown recently with simulations and actual data that MEG source analysis based on individual coregistration with a spherical head model and group analysis based on a template brain is reasonably accurate (Steinstraeter et al., 2009). Specifically, their study showed that the localization errors were not larger than the typical spatial uncertainty in group analysis based on Talairach normalization of an order of ∼5 mm in the visual cortex (Hasnain et al., 1998). Accordingly, MEG source localization based on a template brain MRI has been used as the best approximate when the individual MRI is not available (Jensen et al., 2005; Ross et al., 2009; Fujioka et al., 2010).
Time series of ERD/ERS were calculated from the mean signal power in the beta band at each volume element using the SAM-filtered single-trial data, similar to the described analysis of beta oscillations in the auditory cortex source activity. Second, we applied a spatiotemporal PCA to extract the source locations of temporally correlated beta amplitude modulation across the whole brain volume. The PCA decomposed the four-dimensional source data into principal components (PCs), each of which consisted of a spatial map of sources (factor loading) sharing the same time course (factor score) of modulation of beta power. The normalized intensity at each location in the map indicated how strongly the identified time course is observed at that location. PCA results in PCs that are ordered according to the amount of variance explained in the data. The first PC, which accounted for the most of the signal variance, was used for further analysis as the most representative time course of the beta power modulation across brain volume. The variance accounted for by the first PC was 45% or more across conditions in each individual case. For demonstrating consistency across the group, we applied the described PCA to individual data and compared the time courses of factor scores. Local maxima in the group mean factor-loading map were identified, and a t test was applied to the corresponding individual factor loads at the identified source location to examine whether it was significantly different from zero across the group. Locations with values of p < 0.001 were accepted as significant.
Event-related changes in beta coherence referenced to auditory and motor sources.
To characterize possible functional communication between auditory and motor areas, we were interested in cortico-cortical coherence at beta frequencies across the brain volume specifically from the five seed locations in bilateral A1 and SM1, and SMA, as shown in previous findings of beta-band coherence in auditory-guided repetitive movements (Pollok et al., 2005). Instead of taking a list of a priori-defined locations as the seed locations, we used the five areas, which were identified in our analysis of the beta amplitude comodulation (390 ms condition) as having the largest PCA factor loading (for coordinates, see Table 1). All source waveforms across the brain volume obtained from the beamformer source analysis in each individual were used to compute beta coherence between a given location and one of the seed locations, defined as follows: The complex waveforms x(t) and y(t) were obtained with bandpass filtering between 18 and 25 Hz, and the Hilbert transform was applied to all single trials (i = 1 … N) of source waveforms. The coherence is a complex measure of phase synchrony between the waveforms x(t) and y(t) with its absolute value between 0 and 1 and is calculated for each point in time. However, because the MEG source analysis is limited in spatial resolution, a spread of neural activity across adjacent areas may cause an apparent strong coherence specifically in the vicinity of the seed location. To avoid this, we confined our result to event-related changes in coherence but not the absolute level of the coherence itself. This approach was based on the assumption that functionally relevant coherence fluctuates over time according to the brain process, while apparent coherence due to the spatial filtering should remain constant over the course of time interval of interest rather than dynamically changing with stimulation. Therefore, we subtracted the temporal mean of coherence and analyzed temporal fluctuations only. Moreover, we applied a discrete Fourier transform to the coherence changes and considered only changes that followed the rhythm of the auditory stimulus and termed the resulting measure the event-related coherence. For group statistics, we applied a Rayleigh test (Mardia, 1972; Fisher, 1993), which tests for consistency in the phase of the time course of coherence change across the group.
To reveal the global time courses of event-related coherence, we again applied PCA to decompose the time series of coherence data into principal time series and related volumetric maps. Although the absolute value of coherence measures ranges from 0 to 1, it decreases with increasing number of trials (Maris et al., 2007), and numerically small coherence values can be highly significant in human EEG and MEG data (Srinivasan et al., 1999; Fries et al., 2001; Tallon-Baudry et al., 2001). According to phase statistics (Fisher, 1993), the null distribution of the phase coherence statistics is the circular normal distribution, and for large number of trials (e.g., n > 50), p values can be obtained as p = exp(−nR), R = n · coh2. Local maxima of the squared event-related coherence, which were significant according to the Rayleigh test, were accepted as nodes of a network of synchronous beta oscillations. Time series of event-related coherence were shown for the identified locations.
Results
Beta oscillation in the auditory cortex
The auditory evoked responses were modeled successfully with a pair of equivalent current dipoles in bilateral temporal lobes. The group mean Talairach coordinates of dipole locations were (−42, −22, 2 mm) (left) and (40, −24, 2 mm) (right) and correspond to the medial parts of bilateral primary auditory cortices. Time–frequency analysis of activity from these locations revealed periodic patterns of signal power increase and decrease according to the periodicity of the stimuli (Fig. 1A–D). Power modulation, or ERS/ERD, occurred predominantly in the beta frequency band. The time courses of signal changes in the narrow 20–22 Hz band shown in Figure 1F provide a closer look at its temporal dynamics. Beta activity decreased immediately after stimulus onset, reached a minimum around 200 ms latency. Notably, the time course of this initial beta-ERD did not change with the stimulus conditions. More importantly, only in the periodic stimulus conditions, the following beta-ERS reached the maximum again around the time of next stimulus presentation (Fig. 1A–C), whereas in the case of the randomized stimulation the beta-ERD/ERS was aperiodic and transient in its nature (Fig. 1D). Thus, overall, periodic stimulation resulted in a periodic pattern of beta modulation, and that was not case for the irregular stimulus.
We further examined whether the latency of beta-ERD and the following beta-ERS systematically changed across conditions. As shown in Figure 2, the latencies of the midpoints of the falling slope (t1) and the troughs of ERD (t2) did not change as stimulus timing changed. By contrast, the midpoints of the rising slopes (t3) differed significantly across stimulus timing conditions, such that the rising slope was increasingly shallow with increasing stimulus periodicity. Moreover, despite the fact that, in the irregular condition, the intervals were uniformly randomized between 390 and 780 ms, the beta rebound in that condition was significantly earlier compared with the mean of t3 across the three periodic conditions. This indicates that the periodic pattern of beta modulation was only observed with periodic stimulation and it suggests further that, when the time of the next auditory beat is uncertain, the system prepares for the possibility of an early beat rather than simply expecting the next sound to occur at the mean time of previous beats.
The amount of amplitude change was largest for the slowest stimulus rates; group mean ERD maxima and its SD were 1.63 ± 0.25% at 390 ms, 2.68 ± 0.45% at 585 ms, and 3.52 ± 0.52% at 780 ms. Pairwise comparisons between the different rates revealed significant differences between the fastest condition (390 ms) and each of the slower conditions (390 vs 585 ms, t(17) = 2.033, p = 0.029; 390 vs 780 ms, t(16) = 3.298, p = 0.023; for this comparison, the sample size of N = 12 was corrected based on the Welch–Satterthwaite equation to account for unequal number of trials in each condition). There was no significant difference between the two slower conditions, partly reflecting larger residual noise because of the smaller number of trials.
As for the baseline level of beta amplitude, we compared the mean beta amplitude during the 2.0 s interval preceding a stimulus train to that during the stimulus presentation, and found no significant differences between conditions.
Beta power comodulation across the whole brain
From spatiotemporal PCAs for each stimulus condition, we identified the first PC with its factor score, which reflected common temporal patterns of source activity in the beta power, and its corresponding factor-loading map of source locations, which indicated the strength of the identified temporal pattern at each volume element. Figure 3 illustrates the time series of the first PC, the most common activity across the brain, for the group averaged and individual data. The first PC explained the majority of the variance in beta-power changes for each individual and all conditions (390 ms, 54–80%; 585 ms, 47–75%; 780 ms, 47–77%; random, 52–78%). Furthermore, the time courses of the first PC were consistent across individuals (Fig. 3, bottom panel of each plot), especially for the faster tempi. For the slowest tempo and the randomized tones, the interindividual variability in the time course increased particularly for the upward slope after the initial decrease around 100–300 ms. This apparent larger variability likely reflects the smaller number of trials in these conditions compared the fastest rate condition. However, importantly, this did not affect the timing of the beta rising slope as indicated by similar confidence intervals for the midpoints of the rising slope (Fig. 2). The corresponding map of brain areas in which the beta decrease was statistically significant (p < 0.001) is illustrated in Figure 4. The complete list of Talairach coordinates of local maxima in the factor-loading map and the p values is given in Table 1. The beta network identified by the first PC involved the bilateral auditory cortices [superior temporal gyrus (STG) and/or Heschl's gyrus (HG)] and cuneus/precuneus for the three stimulus rates consistently. Inclusion of motor-related areas became evident mainly for the two faster stimulation rates. The spatial maps for the two fastest rates included the SM1 (precentral and/or postcentral gyrus) and inferior frontal gyrus (IFG). In the fastest tempo, additional motor regions in SMA and cerebellum were also significant. In contrast, parahippocampal gyrus and middle temporal gyrus participated in the beta modulation in the slowest condition only. In the randomized condition, the posterior cingulate gyrus and superior frontal gyrus (SFG) contributed to the beta network. It is important to note that these locations are limited to the areas mostly contributing to the first PC of the primary beta-power time course described here, while beta oscillations across the whole brain involve other brain areas such as those described below in the beta coherence.
Cortico-cortical beta phase synchrony
To examine functional communication between the auditory and the motor systems, the time courses of event-related changes in phase coherence across the whole brain were calculated with five seed areas in bilateral A1, bilateral SM1, and SMA for all conditions. In this analysis, we again applied PCA to extract the first PC as the temporal factor score and the corresponding factor-loading map for each of the event-related coherence data with each of the five different seeds in each stimulus condition. Figure 5 indicates the time course of the factor scores showing periodic patterns of beta coherence modulated with the stimulus rhythm in all conditions. Interestingly, coherence with right and left primary auditory cortex areas was characterized by a decrease following stimulus onset. Thus, it appears that the time series of change in coherence resembled closely the modulation of the beta amplitude. In contrast, coherence with sensorimotor areas increased after stimulus onset. It is important to note that event-related coherence showed maxima/minima at almost the same latency consistently across all periodic stimulation (Fig. 5, dashed lines).
A snapshot of the spatiotemporal pattern of these coherence data is shown for the seed in the right HG in Figure 6, in which the brain areas were identified as significantly involved in stimulus-locked fluctuation in coherence. Strongest coherence changes were found in auditory areas close to the seed (Fig. 6B) and the left cingulate gyrus, likely related to SMA and/or pre-SMA (Fig. 6A). Further local maxima were found in the right postcentral gyrus, left anterior cingulated gyrus, left inferior temporal gyrus, and right parahippocampal gyrus (Fig. 6C–F).
The complete list of brain areas connected through event-related coherence changes is indicated in Table 2. As shown in the list, ipsilateral connections to STG, middle/inferior frontal gyrus, and inferior parietal lobules and contralateral connections to medial frontal gyri were found with seed locations in primary auditory cortices. However, no significant coherence between left and right auditory cortices was found. More importantly, if the seed location was in the sensorimotor areas, coherence was mainly significant with motor-related areas, but not in the auditory areas.
Discussion
The current study revealed three main findings regarding beta oscillations induced through passive listening to isochronous sound stimuli. First, auditory cortices showed a periodic beta modulation consisting of an initial rapid beta decrease, immediately following the stimulus onset and the subsequent rebound. Importantly, only the latter depended on the stimulus rate in the regular conditions, suggesting representation of the internalized interval. Second, the brain areas showing temporally correlated beta modulation with that in the auditory cortices included motor-related cortical and subcortical areas even in the absence of movement. Third, neural synchrony measured as cortico-cortical phase coherence at beta frequencies, was modulated with the sound rhythm in a consistent manner, also spanning across auditory and motor systems.
Rebound of beta activity as representation of internalized timing
The rising slope of the beta rebound but not the initial falling slope varied according to the stimulus rate in the auditory cortex as well as across an extensive brain network. This suggests an endogenous nature of the beta rebound, in line with our previous MEG study (Fujioka et al., 2009) in which a beta decrement in the auditory cortex was absent after an occasional stimulus omission whereas the subsequent rising slope was not affected until the timing of the omission. After the omission, the rising slope continued toward the next stimulus onset. Thus, in both studies, we observed that the amplitude of beta oscillation peaked just before the next stimulus during isochronous stimulation. The faster the tempo of the stimulus, the steeper the rising slope became. Thus, the beta rebound may reflect a neural mechanism for predictive timing information about when the next stimulus is expected. This is a sharp contrast to the case in the irregular condition in which the slope was significantly shorter than in any of the regular conditions and much shorter than expected from the mean of the onset-to-onset intervals. Thus, the predictable periodic characteristic of the auditory input is critical to the beta dynamics described here. Nevertheless, the beta rebound after the initial decrease in the irregular condition may also indicate some anticipatory process because the next stimulus can be still expected within the narrow time window.
Recent studies also associated induced beta increase with endogenous processes: internally driven choices in a search task (Buschman and Miller, 2009), anticipation to the predictive timing of the task-relevant cue (Saleh et al., 2010), and top-down attentional processing in ambiguous visual perception (Okazaki et al., 2008). Our results particularly resemble those of human intracranial recordings in which the beta increase peaked just before an informative cue presented with predictable timing (Saleh et al., 2010). Thus, extending our previous findings (Snyder and Large, 2005; Fujioka et al., 2009), here we propose that the observed time course of beta rebound likely reflects endogenous processing of internalized timing.
Dissociation between beta suppression and rebound
In contrast to the beta rebound, the initial beta decrease after the stimulus onset was almost identical regardless of the regularity or rate of the stimulus sequence. This suggests that this beta decrease reflects an exogenous process driven by the stimulus input, as opposed to the endogenous nature of the following beta increase described above. The functional difference between beta suppression and rebound found here is in line with the literature. It is widely accepted that beta desynchronization indicates an “active change” in the status of sensorimotor processing (Pfurtscheller and Lopes da Silva, 1999), whereas the beta increase is associated with “maintaining an existing steady state in motor control” (Gilbertson et al., 2005; Pogosyan et al., 2009), “immobilization” (Salmelin et al., 1995), or “reset of an existing cortical network” (Pfurtscheller et al., 2005), reflecting a variety of task and attentional demands. For sensorimotor tasks, slightly different cortical areas seem to be involved in beta suppression and rebound (Salmelin and Hari, 1994; Salmelin et al., 1995; Jurkiewicz et al., 2006; Caetano et al., 2007; Bardouille et al., 2010). Our previous and current findings extend these studies by demonstrating that the difference between suppression and rebound is also manifested in their temporal dynamics.
Using a passive listening paradigm, our study discovered that periodic beta power modulation consisting of these two different processes does not require an explicit motor task. This suggests the possible leading role of the auditory system for propagation of beta oscillation in auditory-guided movement, of which the beneficial effect is found in Parkinson's disease (PD) patients (Freeman et al., 1993; McIntosh et al., 1997). A regular pattern of beta modulation caused by the regular sound might help temporarily alleviate an excessive beta power associated with movement initiation problems in PD (Brown, 2003; Hutchison et al., 2004) by possibly modulating membrane hyperpolarization in thalamic nuclei (Llinás et al., 2005). However, as only a few studies with auditory stimuli have reported beta suppression (Kaiser et al., 2002; Pesonen et al., 2006; Kim and Chung, 2008) and beta increase (Caetano et al., 2007), future research is necessary to examine the nature of beta activities in auditory tasks more in detail. Also, it remains unclear how the observed beta activity is differently influenced by the stimulus interval and its perceptual correlates as the latter is influenced by auditory sequence contexts or individual preference with musical expertise (Jones and McAuley, 2005; Repp and Doggett, 2007). Another remaining question is how the two processes might use different pathways to carry out auditory-sensorimotor communication. A different form of auditory-motor transformation such as voice perceptions involving posterior auditory cortices (Warren et al., 2005) could also activate beta modulation, in line with available evidence for speech perception based on articulatory characteristics (Eulitz and Obleser, 2007) and audiovisual multimodal integration (Fingelkurts et al., 2007).
Rate effects on beta modulation in the motor system
We found temporally correlated beta amplitude modulations in motor-related brain areas and in auditory cortices. There was extended involvement of IFG, SMA, and cerebellum in the beta modulation at faster stimulus rates. To our knowledge, this is the first demonstration of rate-specific dynamic configuration of oscillatory beta networks in passive listening to auditory rhythms. The results are generally in line with previous fMRI/PET studies showing that movement rate affects the spatial extent and magnitude of contribution to the activity in different motor regions (Blinkenberg et al., 1996; Rao et al., 1996; Sadato et al., 1996; Schlaug et al., 1996; Jäncke et al., 1998; Wildgruber et al., 2001; Riecker et al., 2003), with decreased activity in the basal ganglia and increased activity in cerebellum and primary motor cortex with increased rate of movements.
Coherence measured by beta phase synchrony
Our measure of cortico-cortical coherence was limited to event-related changes over time, to strictly avoid spatial leakage effects caused by MEG source analysis. This entirely data-driven analysis successfully identified a shared time course of beta coherence that followed the sound rhythm in auditory and sensorimotor networks. Our results indicate that the sensorimotor cortices clearly contributes the modulation of beta coherence with numerous brain areas during passive listening, and that common modulation of synchrony was mainly expressed “within” a modality-specific neural network (e.g., both for auditory and sensorimotor systems) but less expressed between the two systems and almost absent between left and right auditory areas. Thus, phase synchrony in beta oscillations may play an important role for connectivity within specific functional networks, whereas coupling between networks may rely on different mechanisms. Nevertheless, the most important finding here is that the synchrony was modulated directly by the auditory input across the different brain areas as a function of the stimulus rate. This finding is in agreement with the current model of a functional network that modulation of synchrony is a mechanism of dynamic network reconfiguration (Womelsdorf et al., 2007). Mechanisms underlying the beta dynamics could be a frequency transition of oscillatory circuits through synaptic coupling alterations (Kopell et al., 2000; Jensen et al., 2005), changes in nonsynaptic coupling through gap junction conductance in cortical neurons (Roopun et al., 2006), and synaptic plasticity in cortico-striatum connections (Fino et al., 2005). The latter is particularly interesting in that the striato-thalamo-cortical circuit is involved in coincidence detection of oscillatory activities for temporal coding (Matell and Meck, 2004).
Conclusions
In conclusion, our findings support the concept of auditory-motor facilitation in that the periodic modulation of beta activity following fast-paced regular auditory stimuli could aid the initiation of movement. More importantly, the beta rebound before the onset of the next beat appears to encode an internal representation of the predictable time interval, which would likely help predictive movement planning in a time-sensitive manner. Phase synchrony of beta oscillations plays a role in configuring motor-related networks, whereas functional coupling between auditory and motor systems is mainly expressed as correlated changes in beta power. Future studies are necessary to verify the exact nature of sensorimotor interaction expressed in beta processes and to define the precise pathways involved in auditory-motor communication in temporal processing.
Footnotes
This work was supported by Canadian Institutes of Health Research Grant 81135 (B.R.), the Centre for Stroke Rehabilitation of the Heart and Stroke Foundation of Ontario (T.F.), and the Canadian Foundation for Innovation. We sincerely thank Drs. William D. Hutchison and Norihiro Sadato for the insightful comments on an earlier version of this manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Takako Fujioka, Rotman Research Institute, 3560 Bathurst Street, Toronto, ON M6A 2E1, Canada. tfujioka{at}rotman-baycrest.on.ca