Is there frequency-specificity in the motor control of walking? The putative differential role of alpha and beta oscillations.

Alpha and beta oscillations have been assessed thoroughly during walking due to their potential role as proxies of the corticoreticulospinal tract (CReST) and corticospinal tract (CST), respectively. Given that damage to a descending tract after stroke can cause walking deficits, detailed knowledge of how these oscillations mechanistically contribute to walking could be utilized in strategies for post-stroke locomotor recovery. In this review, the goal was to summarize, synthesize, and discuss the existing evidence on the potential differential role of these oscillations on the motor descending drive, the effect of transcranial alternate current stimulation (tACS) on neurotypical and post-stroke walking, and to discuss remaining gaps in knowledge, future directions, and methodological considerations. Electrophysiological studies of corticomuscular, intermuscular, and intramuscular coherence during walking clearly demonstrate that beta oscillations are predominantly present in the dorsiflexors during the swing phase and may be absent post-stroke. The role of alpha oscillations, however, has not been pinpointed as clearly. We concluded that both animal and human studies should focus on the electrophysiological characterization of alpha oscillations and their potential role to the CReST. Another approach in elucidating the role of these oscillations is to modulate them and then quantify the impact on walking behavior. This is possible through tACS, whose beneficial effect on walking behavior (including boosting of beta oscillations in intramuscular coherence) has been recently demonstrated in both neurotypical adults and stroke patients. However, these studies still do not allow for specific roles of alpha and beta oscillations to be delineated because the tACS frequency used was much lower (i.e., individualized calculated gait frequency was used). Thus, we identify a main gap in the literature, which is tACS studies actually stimulating at alpha and beta frequencies during walking. Overall, we conclude that for beta oscillations there is a clear connection to descending drive in the corticospinal tract. The precise relationship between alpha oscillations and CReST remains elusive due to the gaps in the literature identified here. However, better understanding the role of alpha (and beta) oscillations in the motor control of walking can be used to progress and develop rehabilitation strategies for promoting locomotor recovery.

Tuning Neural Synchronization: The Role of Variable Oscillation Frequencies in Neural Circuits.

Brain oscillations emerge during sensory and cognitive processes and have been classified into different frequency bands. Yet, even within the same frequency band and between nearby brain locations, the exact frequencies of brain oscillations can differ. These frequency differences (detuning) have been largely ignored and play little role in current functional theories of brain oscillations. This contrasts with the crucial role that detuning plays in synchronization theory, as originally derived in physical systems. Here, we propose that detuning is equally important to understand synchronization in biological systems. Detuning is a critical control parameter in synchronization, which is not only important in shaping phase-locking, but also in establishing preferred phase relations between oscillators. We review recent evidence that frequency differences between brain locations are ubiquitous and essential in shaping temporal neural coordination. With the rise of powerful experimental techniques to probe brain oscillations, the contributions of exact frequency and detuning across neural circuits will become increasingly clear and will play a key part in developing a new understanding of the role of oscillations in brain function.

Empirically constrained network models for contrast-dependent modulation of gamma rhythm in V1.

Gamma oscillations are thought to play a key role in neuronal network function and neuronal communication, yet the underlying generating mechanisms have not been fully elucidated to date. At least partly, this may be due to the fact that even in simple network models of interconnected inhibitory (I) and excitatory (E) neurons, many parameters remain unknown and are set based on practical considerations or by convention. Here, we mitigate this problem by requiring PING (Pyramidal Interneuron Network Gamma) models to simultaneously satisfy a broad set of criteria for realistic behaviour based on empirical data spanning both the single unit (spikes) and local population (LFP) levels while unknown parameters are varied. By doing so, we were able to constrain the parameter ranges and select empirically valid models. The derived model constraints implied weak rather than strong PING as the generating mechanism for gamma, connectivity between E and I neurons within specific bounds, and variations of the external input to E but not I neurons. Constrained models showed valid behaviours, including gamma frequency increases with contrast and power saturation or decay at high contrasts. Using an empirically-validated model we studied the route to gamma instability at high contrasts. This involved increased heterogeneity of E neurons with increasing input triggering a breakdown of I neuron pacemaker function. Further, we illustrate the model's capacity to resolve disputes in the literature concerning gamma oscillation properties and GABA conductance proxies. We propose that the models derived in our study will be useful for other modelling studies, and that our approach to the empirical constraining of PING models can be expanded when richer empirical datasets become available. As local gamma networks are the building blocks of larger networks that aim to understand complex cognition through their interactions, there is considerable value in improving our models of these building blocks.

Functional brain networks of patients with epilepsy exhibit pronounced multiscale periodicities, which correlate with seizure onset.

Epileptic seizure detection and prediction by using noninvasive measurements such as scalp EEG signals or invasive, intracranial recordings, has been at the heart of epilepsy studies for at least three decades. To this end, the most common approach has been to consider short-length recordings (several seconds to a few minutes) around a seizure, aiming to identify significant changes that occur before or during seizures. An inherent assumption in this approach is the presence of a relatively constant EEG activity in the interictal period, which is interrupted by seizure occurrence. Here, we examine this assumption by using long-duration scalp EEG data (21-94 hr) in nine patients with epilepsy, based on which we construct functional brain networks. Our results reveal that these networks vary over time in a periodic fashion, exhibiting multiple peaks at periods ranging between 1 and 24 hr. The effects of seizure onset on the functional brain network properties were found to be considerably smaller in magnitude compared to the changes due to these inherent periodic cycles. Importantly, the properties of the identified network periodic components (instantaneous phase) were found to be strongly correlated to seizure onset, especially for the periodicities around 3 and 5 hr. These correlations were found to be largely absent between EEG signal periodicities and seizure onset, suggesting that higher specificity may be achieved by using network-based metrics. In turn, this implies that more robust seizure detection and prediction can be achieved if longer term underlying functional brain network periodic variations are taken into account.

Graph Theoretical Characteristics of EEG-Based Functional Brain Networks in Patients With Epilepsy: The Effect of Reference Choice and Volume Conduction.

It is well-established that both volume conduction and the choice of recording reference (montage) affect the correlation measures obtained from scalp EEG, both in the time and frequency domains. As a result, a number of correlation measures have been proposed aiming to reduce these effects. In our previous work, we have showed that scalp-EEG based functional brain networks in patients with epilepsy exhibit clear periodic patterns at different time scales and that these patterns are strongly correlated to seizure onset, particularly at shorter time scales (around 3 and 5 h), which has important clinical implications. In the present work, we use the same long-duration clinical scalp EEG data (multiple days) to investigate the extent to which the aforementioned results are affected by the choice of reference choice and correlation measure, by considering several widely used montages as well as correlation metrics that are differentially sensitive to the effects of volume conduction. Specifically, we compare two standard and commonly used linear correlation measures, cross-correlation in the time domain, and coherence in the frequency domain, with measures that account for zero-lag correlations: corrected cross-correlation, imaginary coherence, phase lag index, and weighted phase lag index. We show that the graphs constructed with corrected cross-correlation and WPLI are more stable across different choices of reference. Also, we demonstrate that all the examined correlation measures revealed similar periodic patterns in the obtained graph measures when the bipolar and common reference (Cz) montage were used. This includes circadian-related periodicities (e.g., a clear increase in connectivity during sleep periods as compared to awake periods), as well as periodicities at shorter time scales (around 3 and 5 h). On the other hand, these results were affected to a large degree when the average reference montage was used in combination with standard cross-correlation, coherence, imaginary coherence, and PLI, which is likely due to the low number of electrodes and inadequate electrode coverage of the scalp. Finally, we demonstrate that the correlation between seizure onset and the brain network periodicities is preserved when corrected cross-correlation and WPLI were used for all the examined montages. This suggests that, even in the standard clinical setting of EEG recording in epilepsy where only a limited number of scalp EEG measurements are available, graph-theoretic quantification of periodic patterns using appropriate montage, and correlation measures corrected for volume conduction provides useful insights into seizure onset.

Epileptic seizure onset correlates with long term EEG functional brain network properties.

We investigated the correlation of epileptic seizure onset times with long term EEG functional brain network properties. To do so, we constructed binary functional brain networks from long-term, multichannel electroencephalographic data recorded from nine patients with epilepsy. The corresponding network properties were quantified using the average network degree. It was found that the network degree (as well as other network properties such as the network efficiency and clustering coefficient) exhibited large fluctuations over time; however, it also exhibited specific periodic temporal structure over different time scales (1.5hr-24hr periods) that was consistent across subjects. We investigated the correlation of the phases of these network periodicities with the seizure onset by using circular statistics. The results showed that the instantaneous phases of the 3.5hr, 5.5hr, 12hr and 24hr network degree periodic components are not uniformly distributed, suggesting that functional network properties are related to seizure generation and occurrence.

Input-dependent frequency modulation of cortical gamma oscillations shapes spatial synchronization and enables phase coding.

Fine-scale temporal organization of cortical activity in the gamma range (∼25-80Hz) may play a significant role in information processing, for example by neural grouping ('binding') and phase coding. Recent experimental studies have shown that the precise frequency of gamma oscillations varies with input drive (e.g. visual contrast) and that it can differ among nearby cortical locations. This has challenged theories assuming widespread gamma synchronization at a fixed common frequency. In the present study, we investigated which principles govern gamma synchronization in the presence of input-dependent frequency modulations and whether they are detrimental for meaningful input-dependent gamma-mediated temporal organization. To this aim, we constructed a biophysically realistic excitatory-inhibitory network able to express different oscillation frequencies at nearby spatial locations. Similarly to cortical networks, the model was topographically organized with spatially local connectivity and spatially-varying input drive. We analyzed gamma synchronization with respect to phase-locking, phase-relations and frequency differences, and quantified the stimulus-related information represented by gamma phase and frequency. By stepwise simplification of our models, we found that the gamma-mediated temporal organization could be reduced to basic synchronization principles of weakly coupled oscillators, where input drive determines the intrinsic (natural) frequency of oscillators. The gamma phase-locking, the precise phase relation and the emergent (measurable) frequencies were determined by two principal factors: the detuning (intrinsic frequency difference, i.e. local input difference) and the coupling strength. In addition to frequency coding, gamma phase contained complementary stimulus information. Crucially, the phase code reflected input differences, but not the absolute input level. This property of relative input-to-phase conversion, contrasting with latency codes or slower oscillation phase codes, may resolve conflicting experimental observations on gamma phase coding. Our modeling results offer clear testable experimental predictions. We conclude that input-dependency of gamma frequencies could be essential rather than detrimental for meaningful gamma-mediated temporal organization of cortical activity.

Periodicity in functional brain networks: application to scalp EEG from epilepsy patients.

Seizure detection and prediction studies using scalp- or intracranial-EEG measurements often focus on short-length recordings around the occurrence of the seizure, normally ranging between several seconds and up to a few minutes before and after the event. The underlying assumption in these studies is the presence of a relatively constant EEG activity in the interictal period, that is presumably interrupted by the occurrence of a seizure, at the time the seizure starts or slightly earlier. In this study, we put this assumption under test, by examining long-duration scalp EEG recordings, ranging between 22 and 72 hours, of five patients with epilepsy. For each patient, we construct functional brain networks, by calculating correlations between the scalp electrodes, and examine how these networks vary in time. The results suggest not only that the network varies over time, but it does so in a periodic fashion, with periods ranging between 11 and 25 hours.

Estimation of functional connectivity from electromagnetic signals and the amount of empirical data required.

An increasing number of neuroimaging studies are concerned with the identification of interactions or statistical dependencies between brain areas. Dependencies between the activities of different brain regions can be quantified with functional connectivity measures such as the cross-correlation coefficient. An important factor limiting the accuracy of such measures is the amount of empirical data available. For event-related protocols, the amount of data also affects the temporal resolution of the analysis. We use analytical expressions to calculate the amount of empirical data needed to establish whether a certain level of dependency is significant when the time series are autocorrelated, as is the case for biological signals. These analytical results are then contrasted with estimates from simulations based on real data recorded with magnetoencephalography during a resting-state paradigm and during the presentation of visual stimuli. Results indicate that, for broadband signals, 50-100 s of data is required to detect a true underlying cross-correlations coefficient of 0.05. This corresponds to a resolution of a few hundred milliseconds for typical event-related recordings. The required time window increases for narrow band signals as frequency decreases. For instance, approximately 3 times as much data is necessary for signals in the alpha band. Important implications can be derived for the design and interpretation of experiments to characterize weak interactions, which are potentially important for brain processing.

Using brain-computer interfaces and brain-state dependent stimulation as tools in cognitive neuroscience.

Large efforts are currently being made to develop and improve online analysis of brain activity which can be used, e.g., for brain-computer interfacing (BCI). A BCI allows a subject to control a device by willfully changing his/her own brain activity. BCI therefore holds the promise as a tool for aiding the disabled and for augmenting human performance. While technical developments obviously are important, we will here argue that new insight gained from cognitive neuroscience can be used to identify signatures of neural activation which reliably can be modulated by the subject at will. This review will focus mainly on oscillatory activity in the alpha band which is strongly modulated by changes in covert attention. Besides developing BCIs for their traditional purpose, they might also be used as a research tool for cognitive neuroscience. There is currently a strong interest in how brain-state fluctuations impact cognition. These state fluctuations are partly reflected by ongoing oscillatory activity. The functional role of the brain state can be investigated by introducing stimuli in real-time to subjects depending on the actual state of the brain. This principle of brain-state dependent stimulation may also be used as a practical tool for augmenting human behavior. In conclusion, new approaches based on online analysis of ongoing brain activity are currently in rapid development. These approaches are amongst others informed by new insight gained from electroencephalography/magnetoencephalography studies in cognitive neuroscience and hold the promise of providing new ways for investigating the brain at work.

The relationship between the visual evoked potential and the gamma band investigated by blind and semi-blind methods.

Gamma Band Activity (GBA) is increasingly studied for its relation with attention, change detection, maintenance of working memory and the processing of sensory stimuli. Activity around the gamma range has also been linked with early visual processing, although the relationship between this activity and the low frequency visual evoked potential (VEP) remains unclear. This study examined the ability of blind and semi-blind source separation techniques to extract sources specifically related to the VEP and GBA in order to shed light on the relationship between them. Blind (Independent Component Analysis-ICA) and semi-Blind (Functional Source Separation-FSS) methods were applied to dense array EEG data recorded during checkerboard stimulation. FSS was performed with both temporal and spectral constraints to identify specifically the generators of the main peak of the VEP (P100) and of the GBA. Source localisation and time-frequency analyses were then used to investigate the properties and co-dependencies between VEP/P100 and GBA. Analysis of the VEP extracted using the different methods demonstrated very similar morphology and localisation of the generators. Single trial time frequency analysis showed higher GBA when a larger amplitude VEP/P100 occurred. Further examination indicated that the evoked (phase-locked) component of the GBA was more related to the P100, whilst the induced component correlated with the VEP as a whole. The results suggest that the VEP and GBA may be generated by the same neuronal populations, and implicate this relationship as a potential mediator of the correlation between the VEP and the Blood Oxygenation Level Dependent (BOLD) effect measured with fMRI.

Identifying spatially overlapping local cortical networks with MEG.

Recent modelling studies (Hadjipapas et al. [2009]: Neuroimage 44:1290-1303) have shown that it may be possible to distinguish between different neuronal populations on the basis of their macroscopically measured (EEG/MEG) mean field. We set out to test whether the different orientation columns contributing to a signal at a specific cortical location could be identified based on the measured MEG signal. We used 1.5deg square, static, obliquely oriented grating stimuli to generate sustained gamma oscillations in a focal region of primary visual cortex. We then used multivariate classifier methods to predict the orientation (left or right oblique) of the stimuli based purely on the time-series data from this one location. Both the single trial evoked response (0-300 ms) and induced post-transient power spectra (300-2,300 ms, 20-70 Hz band) due to the different stimuli were classifiable significantly above chance in 11/12 and 10/12 datasets respectively. Interestingly, stimulus-specific information is preserved in the sustained part of the gamma oscillation, long after perception has occurred and all neuronal transients have decayed. Importantly, the classification of this induced oscillation was still possible even when the power spectra were rank-transformed showing that the different underlying networks give rise to different characteristic temporal signatures.

Can we observe collective neuronal activity from macroscopic aggregate signals?

The fundamental problem faced by noninvasive neuroimaging techniques such as EEG/MEG(1) is to elucidate functionally important aspects of the microscopic neuronal network dynamics from macroscopic aggregate measurements. Due to the mixing of the activities of large neuronal populations in the observed macroscopic aggregate, recovering the underlying network that generates the signal in the absence of any additional information represents a considerable challenge. Recent MEG studies have shown that macroscopic measurements contain sufficient information to allow the differentiation between patterns of activity, which are likely to represent different stimulus-specific collective modes in the underlying network (Hadjipapas, A., Adjamian, P., Swettenham, J.B., Holliday, I.E., Barnes, G.R., 2007. Stimuli of varying spatial scale induce gamma activity with distinct temporal characteristics in human visual cortex. NeuroImage 35, 518-530). The next question arising in this context is whether aspects of collective network activity can be recovered from a macroscopic aggregate signal. We propose that this issue is most appropriately addressed if MEG/EEG signals are to be viewed as macroscopic aggregates arising from networks of coupled systems as opposed to aggregates across a mass of largely independent neural systems. We show that collective modes arising in a network of simulated coupled systems can be indeed recovered from the macroscopic aggregate. Moreover, we show that nonlinear state space methods yield a good approximation of the number of effective degrees of freedom in the network. Importantly, information about hidden variables, which do not directly contribute to the aggregate signal, can also be recovered. Finally, this theoretical framework can be applied to experimental MEG/EEG data in the future, enabling the inference of state dependent changes in the degree of local synchrony in the underlying network.

Induced Gamma activity in primary visual cortex is related to luminance and not color contrast: An MEG study.

Gamma activity in the visual cortex has been reported in numerous EEG studies of coherent and illusory figures. A dominant theme of many such findings has been that temporal synchronization in the gamma band in response to these identifiable percepts is related to perceptual binding of the common features of the stimulus. In two recent studies using magnetoencephalography (MEG) and the beamformer analysis technique, we have shown that the magnitude of induced gamma activity in visual cortex is dependent upon independent stimulus features such as spatial frequency and contrast. In particular, we showed that induced gamma activity is maximal in response to gratings of 3 cycles per degree (3 cpd) of high luminance contrast. In this work, we set out to examine stimulus contrast further by using isoluminant red/green gratings that possess color but not luminance contrast using the same cohort of subjects. We found no induced gamma activity in V1 or visual cortex in response to the isoluminant gratings in these subjects who had previously shown strong induced gamma activity in V1 for luminance contrast gratings.

Functional source separation applied to induced visual gamma activity.

Objective of this work was to explore the performance of a recently introduced source extraction method, FSS (Functional Source Separation), in recovering induced oscillatory change responses from extra-cephalic magnetoencephalographic (MEG) signals. Unlike algorithms used to solve the inverse problem, FSS does not make any assumption about the underlying biophysical source model; instead, it makes use of task-related features (functional constraints) to estimate source/s of interest. FSS was compared with blind source separation (BSS) approaches such as Principal and Independent Component Analysis, PCA and ICA, which are not subject to any explicit forward solution or functional constraint, but require source uncorrelatedness (PCA), or independence (ICA). A visual MEG experiment with signals recorded from six subjects viewing a set of static horizontal black/white square-wave grating patterns at different spatial frequencies was analyzed. The beamforming technique Synthetic Aperture Magnetometry (SAM) was applied to localize task-related sources; obtained spatial filters were used to automatically select BSS and FSS components in the spatial area of interest. Source spectral properties were investigated by using Morlet-wavelet time-frequency representations and significant task-induced changes were evaluated by means of a resampling technique; the resulting spectral behaviours in the gamma frequency band of interest (20-70 Hz), as well as the spatial frequency-dependent gamma reactivity, were quantified and compared among methods. Among the tested approaches, only FSS was able to estimate the expected sustained gamma activity enhancement in primary visual cortex, throughout the whole duration of the stimulus presentation for all subjects, and to obtain sources comparable to invasively recorded data.

Stimuli of varying spatial scale induce gamma activity with distinct temporal characteristics in human visual cortex.

Gamma activity to stationary grating stimuli was studied non-invasively using MEG recordings in humans. Using a spatial filtering technique, we localized gamma activity to primary visual cortex. We tested the hypothesis that spatial frequency properties of visual stimuli may be related to the temporal frequency characteristics of the associated cortical responses. We devised a method to assess temporal frequency differences between stimulus-related responses that typically exhibit complex spectral shapes. We applied this methodology to either single-trial (induced) or time-averaged (evoked) responses in four frequency ranges (0-40, 20-60, 40-80 and 60-100 Hz) and two time windows (either the entire duration of stimulus presentation or the first second following stimulus onset). Our results suggest that stimuli of varying spatial frequency induce responses that exhibit significantly different temporal frequency characteristics. These effects were particularly accentuated for induced responses in the classical gamma frequency band (20-60 Hz) analyzed over the entire duration of stimulus presentation. Strikingly, examining the first second of the responses following stimulus onset resulted in significant loss in stimulus specificity, suggesting that late signal components contain functionally relevant information. These findings advocate a functional role of gamma activity in sensory representation. We suggest that stimulus specific frequency characteristics of MEG signals can be mapped to processes of neuronal synchronization within the framework of coupled dynamical systems.

Imaging the dynamics of the auditory steady-state evoked response.

This study used magnetoencephalography (MEG) to examine the dynamic patterns of neural activity underlying the auditory steady-state response. We examined the continuous time-series of responses to a 32-Hz amplitude modulation. Fluctuations in the amplitude of the evoked response were found to be mediated by non-linear interactions with oscillatory processes both at the same source, in the alpha and beta frequency bands, and in the opposite hemisphere.

Assessing interactions of linear and nonlinear neuronal sources using MEG beamformers: a proof of concept.

OBJECTIVE: This study aimed to explore methods of assessing interactions between neuronal sources using MEG beamformers. However, beamformer methodology is based on the assumption of no linear long-term source interdependencies [VanVeen BD, vanDrongelen W, Yuchtman M, Suzuki A. Localization of brain electrical activity via linearly constrained minimum variance spatial filtering. IEEE Trans Biomed Eng 1997;44:867-80; Robinson SE, Vrba J. Functional neuroimaging by synthetic aperture magnetometry (SAM). In: Recent advances in Biomagnetism. Sendai: Tohoku University Press; 1999. p. 302-5]. Although such long-term correlations are not efficient and should not be anticipated in a healthy brain [Friston KJ. The labile brain. I. Neuronal transients and nonlinear coupling. Philos Trans R Soc Lond B Biol Sci 2000;355:215-36], transient correlations seem to underlie functional cortical coordination [Singer W. Neuronal synchrony: a versatile code for the definition of relations? Neuron 1999;49-65; Rodriguez E, George N, Lachaux J, Martinerie J, Renault B, Varela F. Perception's shadow: long-distance synchronization of human brain activity. Nature 1999;397:430-3; Bressler SL, Kelso J. Cortical coordination dynamics and cognition. Trends Cogn Sci 2001;5:26-36]. METHODS: Two periodic sources were simulated and the effects of transient source correlation on the spatial and temporal performance of the MEG beamformer were examined. Subsequently, the interdependencies of the reconstructed sources were investigated using coherence and phase synchronization analysis based on Mutual Information. Finally, two interacting nonlinear systems served as neuronal sources and their phase interdependencies were studied under realistic measurement conditions. RESULTS: Both the spatial and the temporal beamformer source reconstructions were accurate as long as the transient source correlation did not exceed 30-40 percent of the duration of beamformer analysis. In addition, the interdependencies of periodic sources were preserved by the beamformer and phase synchronization of interacting nonlinear sources could be detected. CONCLUSIONS: MEG beamformer methods in conjunction with analysis of source interdependencies could provide accurate spatial and temporal descriptions of interactions between linear and nonlinear neuronal sources. SIGNIFICANCE: The proposed methods can be used for the study of interactions between neuronal sources.

The missing link: analogous human and primate cortical gamma oscillations.

Recent animal studies highlighting the relationship between functional imaging signals and the underlying neuronal activity have revealed the potential capabilities of non-invasive methods. However, the valuable exchange of information between animal and human studies remains restricted by the limited evidence of direct physiological links between species. In this study we used magnetoencephalography (MEG) to investigate the occurrence of 30-70 Hz (gamma) oscillations in human visual cortex, induced by the presentation of visual stimuli of varying contrast. These oscillations, well described in the animal literature, were observed in retinotopically concordant locations of visual cortex and show striking similarity to those found in primate visual cortex using surgically implanted electrodes. The amplitude of the gamma oscillations increases linearly with stimulus contrast in strong correlation with the gamma oscillations found in the local field potential (LFP) of the macaque. We demonstrate that non-invasive magnetic field measurements of gamma oscillations in human visual cortex concur with invasive measures of activation in primate visual cortex, suggesting both a direct representation of underlying neuronal activity and a concurrence between human and primate cortical activity.