Emergent oscillations in a realistic network: the role of inhibition and the effect of the spatiotemporal distribution of the input

Inhibitory interneuron circuits at cortical and spinal levels are associated with individual differences in corticomuscular coherence during isometric voluntary contraction

Corticomuscular coherence (CMC) is an oscillatory synchronization of 15-35 Hz (β-band) between electroencephalogram (EEG) of the sensorimotor cortex and electromyogram of contracting muscles. Although we reported that the magnitude of CMC varies among individuals, the physiological mechanisms underlying this variation are still unclear. Here, we aimed to investigate the associations between CMC and intracortical inhibition (ICI) in the primary motor cortex (M1)/recurrent inhibition (RI) in the spinal cord, which probably affect oscillatory neural activities. Firstly, we quantified ICI from changes in motor-evoked potentials induced by paired-pulse transcranial magnetic stimulation in M1 during tonic isometric voluntary contraction of the first dorsal interosseous. ICI showed a significant, negative correlation with the strength of EEG β-oscillation, but not with the magnitude of CMC across individuals. Next, we quantified RI from changes in H-reflexes induced by paired-pulse electrical nerve stimulation to the posterior tibial nerve during isometric contraction of the soleus muscle. We observed a significant, positive correlation between RI and peak CMC across individuals. These results suggest that the local inhibitory interneuron networks in cortical and spinal levels are associated with the oscillatory activity in corticospinal loop.

Intermuscular Coherence in Normal Adults: Variability and Changes with Age

We investigated beta-band intermuscular coherence (IMC) in 92 healthy adults stratified by decade of age, and analysed variability between and within subjects. In the dominant upper limb, IMC was estimated between extensor digitorum communis and first dorsal interosseous as well as between flexor digitorum superficialis and first dorsal interosseous. In the ipsilateral lower limb, IMC was measured between medial gastrocnemius and extensor digitorum brevis as well as between tibialis anterior and extensor digitorum brevis. Age-related changes in IMC were analysed with age as a continuous variable or binned by decade. Intrasession variance of IMC was examined by dividing sessions into pairs of epochs and comparing coherence estimates between these pairs. Eight volunteers returned for a further session after one year, allowing us to compare intrasession and intersession variance. We found no age-related changes in IMC amplitude across almost six decades of age, allowing us to collate data from all ages into an aggregate normative dataset. Interindividual variability ranged over two orders of magnitude. Intrasession variance was significantly greater than expected from statistical variability alone, and intersession variance was even larger. Potential contributors include fluctuations in task performance, differences in electrode montage and short-term random variation in central coupling. These factors require further exploration and, where possible, minimisation. This study provides evidence that coherence is remarkably robust to senescent changes in the nervous system and provides a large normative dataset for future applications of IMC as a biomarker in disease states.

Diversity of LFPs Activated in Different Target Regions by a Common CA3 Input

Identifying the pathways contributing to local field potential (LFP) events and oscillations is essential to determine whether synchronous interregional patterns indicate functional connectivity. Here, we studied experimentally and numerically how different target structures receiving input from a common population shape their LFPs. We focused on the bilateral CA3 that sends gamma-paced excitatory packages to the bilateral CA1, the lateral septum, and itself (recurrent input). The CA3-specific contribution was isolated from multisite LFPs in target regions using spatial discrimination techniques. We found strong modulation of LFPs by target-specific features, including the morphology and population arrangement of cells, the timing of CA3 inputs, volume conduction from nearby targets, and co-activated inhibition. Jointly they greatly affect the LFP amplitude, profile, and frequency characteristics. For instance, ipsilateral (Schaffer) LFPs occluded contralateral ones, and septal LFPs arise mostly from remote sources while local contribution from CA3 input was minor. In the CA3 itself, gamma waves have dual origin from local networks: in-phase excitatory and nearly antiphase inhibitory. Also, waves may have different duration and varying phase in different targets. These results indicate that to explore the cellular basis of LFPs and the functional connectivity between structures, besides identifying the origin population/s, target modifiers should be considered.

Inter-network interactions: impact of connections between oscillatory neuronal networks on oscillation frequency and pattern

Oscillations in electrical activity are a characteristic feature of many brain networks and display a wide variety of temporal patterns. A network may express a single oscillation frequency, alternate between two or more distinct frequencies, or continually express multiple frequencies. In addition, oscillation amplitude may fluctuate over time. The origin of this complex repertoire of activity remains unclear. Different cortical layers often produce distinct oscillation frequencies. To investigate whether interactions between different networks could contribute to the variety of oscillation patterns, we created two model networks, one generating on its own a relatively slow frequency (20 Hz; slow network) and one generating a fast frequency (32 Hz; fast network). Taking either the slow or the fast network as source network projecting connections to the other, or target, network, we systematically investigated how type and strength of inter-network connections affected target network activity. For high inter-network connection strengths, we found that the slow network was more effective at completely imposing its rhythm on the fast network than the other way around. The strongest entrainment occurred when excitatory cells of the slow network projected to excitatory or inhibitory cells of the fast network. The fast network most strongly imposed its rhythm on the slow network when its excitatory cells projected to excitatory cells of the slow network. Interestingly, for lower inter-network connection strengths, multiple frequencies coexisted in the target network. Just as observed in rat prefrontal cortex, the target network could express multiple frequencies at the same time, alternate between two distinct oscillation frequencies, or express a single frequency with alternating episodes of high and low amplitude. Together, our results suggest that input from other oscillating networks may markedly alter a network's frequency spectrum and may partly be responsible for the rich repertoire of temporal oscillation patterns observed in the brain.

Cytoarchitectonic and dynamic origins of giant positive local field potentials in the dentate gyrus

To determine why some pathways but not others produce sizable local field potentials (LFPs) and how far from the source can these be recorded, complementary experimental analyses and realistic modeling of specific brain structures are required. In the present study, we combined multiple in vivo linear recordings in rats and a tridimensional finite element model of the dentate gyrus, a curved structure displaying abnormally large positive LFPs. We demonstrate that the polarized dendritic arbour of granule cells (GCs), combined with the curved layered configuration of the population promote the spatial clustering of GC currents in the interposed hilus and project them through the open side at a distance from cell domains. LFPs grow up to 20 times larger than observed in synaptic sites. The dominant positive polarity of hilar LFPs was only produced by the synchronous activation of GCs in both blades by either somatic inhibition or dendritic excitation. Moreover, the corresponding anatomical pathways must project to both blades of the dentate gyrus as even a mild decrease in the spatial synchronization resulted in a dramatic reduction in LFP power in distant sites, yet not in the GC domains. It is concluded that the activation of layered structures may establish sharply delimited spatial domains where synaptic currents from one or another input appear to be segregated according to the topology of afferent pathways and the cytoarchitectonic features of the target population. These also determine preferred directions for volume conduction in the brain, of relevance for interpretation of surface EEG recordings.

Neural adaptation facilitates oscillatory responses to static inputs in a recurrent network of ON and OFF cells

We investigate the role of adaptation in a neural field model, composed of ON and OFF cells, with delayed all-to-all recurrent connections. As external spatially profiled inputs drive the network, ON cells receive inputs directly, while OFF cells receive an inverted image of the original signals. Via global and delayed inhibitory connections, these signals can cause the system to enter states of sustained oscillatory activity. We perform a bifurcation analysis of our model to elucidate how neural adaptation influences the ability of the network to exhibit oscillatory activity. We show that slow adaptation encourages input-induced rhythmic states by decreasing the Andronov-Hopf bifurcation threshold. We further determine how the feedback and adaptation together shape the resonant properties of the ON and OFF cell network and how this affects the response to time-periodic input. By introducing an additional frequency in the system, adaptation alters the resonance frequency by shifting the peaks where the response is maximal. We support these results with numerical experiments of the neural field model. Although developed in the context of the circuitry of the electric sense, these results are applicable to any network of spontaneously firing cells with global inhibitory feedback to themselves, in which a fraction of these cells receive external input directly, while the remaining ones receive an inverted version of this input via feedforward di-synaptic inhibition. Thus the results are relevant beyond the many sensory systems where ON and OFF cells are usually identified, and provide the backbone for understanding dynamical network effects of lateral connections and various forms of ON/OFF responses.

Developmental tuning and decay in senescence of oscillations linking the corticospinal system

There is increasing evidence of the importance of synchronous activity within the corticospinal system for motor control. We compared oscillatory activity in the primary sensorimotor cortex [EEG of sensorimotor cortex (SMC-EEG)] and a motor neuronal pool [surface electromyogram of opponens pollicis (OP-EMG)], and their coherence in children (4-12 years of age), young adults (20-35 years of age), and elderly adults (>55 years of age). The ratio between lower (2-13 Hz) and higher (14-32 Hz) frequencies in both SMC-EEG and OP-EMG decreased with age, correlating inversely with motor performance. There was evidence for larger, more distributed cortical networks in the children and elderly compared with young adults. Corticomuscular coherence (CMC) was present in all age groups and shifted between frontal and parietal cortical areas. In children, CMC was smaller and less stationary in amplitude and frequency than in adults. Young adults had single peaks of CMC clustered near the modal frequency (23 Hz); multiple peaks with a broad spread of frequencies occurred in children and the elderly; the further the frequency of the maximum peak CMC was from 23 Hz, the poorer the performance. CMC amplitude was inversely related to performance in young adults but was not modulated in relation to performance in children and the elderly. We propose that progressive fine-tuning of the frequency coding and stabilization of the dynamic properties within and between corticospinal networks occurs during adolescence, refining the capacity for efficient dynamic communication in adulthood. In old age, blurring of the tuning between networks and breakdown in their integration occurs and is likely to contribute to a decrement in motor control.

Renshaw cell recurrent inhibition improves physiological tremor by reducing corticomuscular coupling at 10 Hz

Corticomuscular coherence between the primary motor cortex (M1) and hand muscle electromyograms (EMG) occurs at approximately 20 Hz but is rarely seen at approximately 10 Hz. This is unexpected, because M1 has oscillations at both frequencies, which are effectively transmitted to the spinal cord via the corticospinal tract. We have previously speculated that a specific "neural filter" may selectively reduce coherence at approximately 10 Hz. This would have functional utility in minimizing physiological tremor, which often has a dominant component around this frequency. Recurrent inhibition via Renshaw cells in the spinal cord is a putative neural substrate for such a filter. Here we investigate this system in more detail with a biophysically based computational model. Renshaw cell recurrent inhibition reduced EMG oscillations at approximately 10 Hz, and also reduced corticomuscular coherence at this frequency (from 0.038 to 0.014). Renshaw cell inhibitory feedback also generated synchronous oscillations in the motoneuron pool at approximately 30 Hz. We show that the effects at 10 Hz and 30 Hz can both be understood from the dynamics of the inhibitory feedback loop. We conclude that recurrent inhibition certainly plays an important role in reducing 10 Hz oscillations in muscle, thereby decreasing tremor amplitude. However, our quantitative results suggest it is unlikely to be the only system for tremor reduction, and probably acts in concert with other neural circuits which remain to be elucidated.

Feedback-induced gain control in stochastic spiking networks

The joint influence of recurrent feedback and noise on gain control in a network of globally coupled spiking leaky integrate-and-fire neurons is studied theoretically and numerically. The context of our work is the origin of divisive versus subtractive gain control, as mixtures of these effects are seen in a variety of experimental systems. We focus on changes in the slope of the mean firing frequency-versus-input bias (f-I) curve when the gain control signal to the cells comes from the cells' output spikes. Feedback spikes are modeled as alpha functions that produce an additive current in the current balance equation. For generality, they occur after a fixed minimum delay. We show that purely divisive gain control, i.e. changes in the slope of the f-I curve, arises naturally with this additive negative or positive feedback, due to a linearizing actions of feedback. Negative feedback alone lowers the gain, accounting in particular for gain changes in weakly electric fish upon pharmacological opening of the feedback loop as reported by Bastian (J Neurosci 6:553-562, 1986). When negative feedback is sufficiently strong it further causes oscillatory firing patterns which produce irregularities in the f-I curve. Small positive feedback alone increases the gain, but larger amounts cause abrupt jumps to higher firing frequencies. On the other hand, noise alone in open loop linearizes the f-I curve around threshold, and produces mixtures of divisive and subtractive gain control. With both noise and feedback, the combined gain control schemes produce a primarily divisive gain control shift, indicating the robustness of feedback gain control in stochastic networks. Similar results are found when the "input" parameter is the contrast of a time-varying signal rather than the bias current. Theoretical results are derived relating the slope of the f-I curve to feedback gain and noise strength. Good agreement with simulation results are found for inhibitory and excitatory feedback. Finally, divisive feedback is also found for conductance-based feedback (shunting or excitatory) with and without noise.

Selectivity for grasp in local field potential and single neuron activity recorded simultaneously from M1 and F5 in the awake macaque monkey

The selectivity for object-specific grasp in local field potentials (LFPs) was investigated in two awake macaque monkeys trained to observe, reach out, grasp and hold one of six objects presented in a pseudorandom order. Simultaneous, multiple electrode recordings were made from the hand representations of primary motor cortex (M1) and ventral premotor cortex (area F5). LFP activity was well developed during the observation and hold periods of the task, especially in the beta-frequency range (15-30 Hz). Selectivity of LFP activity for upcoming grasp was rare in the observation period, but common during stable grasp. The majority of M1 (90 of 92) and F5 (81 of 97) sites showed selectivity for at least one frequency, which was maximal in the beta range but also present at higher frequencies (30-50 Hz). When the LFP power associated with grasp of a specific object was large in the beta-frequency range, it was usually of low power in the higher 30-50 Hz range, and vice-versa. Simple hook grips involving flexion of one or more fingers were associated with large beta power, whereas more complex grips involving the thumb (e.g., precision grip) were associated with small beta power. At many M1 sites, there was a highly significant inverse relationship between the tuning of spikes (including those of identified pyramidal tract neurons) and beta-range LFP for different grasps, whereas a positive correlation was found at higher frequencies (30-50 Hz). High levels of beta LFP and low pyramidal cell spike rate may reflect a common mechanism used to control motor set during different types of grasp.

Oscillatory interactions between sensorimotor cortex and the periphery

Field potential recordings from motor cortex show oscillations in the beta-band (∼20 Hz), which are coherent with similar oscillations in the activity of contralateral contracting muscles. Recent findings have revised concepts of how this activity might be generated in the cortex, suggesting it could achieve useful computation. Other evidence shows that these oscillations engage not just motor structures, but also return from muscle to the central nervous system via feedback afferent pathways. Somatosensory cortex has strong beta-band oscillations, which are synchronised with those in motor cortex, allowing oscillatory sensory reafference to be interpreted in the context of the oscillatory motor command which produced it.

Muscle responses to transcranial stimulation in man depend on background oscillatory activity

Muscle responses to transcranial stimulation show high sweep-to-sweep variability, which may reflect an underlying noise process in the motor system. We examined whether response amplitude correlated with the level of prestimulus background EMG, and network oscillations. Transcranial magnetic or electrical stimulation was delivered to primary motor cortex whilst human subjects performed a precision grip task known to promote beta-band ( approximately 20 Hz) cortical oscillations. Responses were recorded from two intrinsic hand muscles. Response magnitude correlated significantly with the level of background EMG (mean r(2) = 0.20). Using a novel wavelet method, we quantified the amplitude and phase of oscillations in prestimulus sensorimotor EEG. Surprisingly, response magnitude showed no significant correlation with EEG oscillations at any frequency. However, oscillations in the prestimulus EMG were significantly correlated with response size; the correlation coefficient had peaks around 20 Hz. When oscillations in one muscle were used to predict response amplitude in a different muscle, correlations were substantially smaller. Finally, for each recording, we calculated the best possible prediction of response size obtainable from up to 20 measures of prestimulus EEG and EMG oscillations. Such optimal predictions had low correlation coefficients (mean r(2) = 0.2; 76% were below 0.3). We conclude that prestimulus oscillations, mainly in the beta-band, do explain some of the variability in responses to transcranial stimulation. Oscillations may likewise increase the noise of natural motor processing, explaining why this form of network activity is usually suppressed prior to dynamic movements. However, the majority of the variation is determined by other factors, which are not accessible by noninvasive recordings.

Network oscillations and intrinsic spiking rhythmicity do not covary in monkey sensorimotor areas

We investigated the relationship between local field potential (LFP) oscillations and intrinsic spiking rhythmicity in the sensorimotor system, because intrinsic rhythmicity has the potential to enhance network oscillations. LFPs and 918 single units were recorded from primary motor cortex (M1), primary somatosensory cortex (S1, areas 3a and 2), posterior parietal cortex (area 5) and the deep cerebellar nuclei (DCN). Some cells were antidromically identified as pyramidal tract neurons (PTNs). In each area the power of approximately 20 Hz LFP oscillations was assessed during periods of steady holding, when such oscillations have previously been shown to be maximal in M1. Oscillations were strongest in area 5 and weakest in the DCN. Using a previously developed method, the postspike distance-to-threshold trajectory was determined from the interspike interval histogram for each cell. Many cells had significant peaks, suggesting an intrinsic tendency towards rhythmic firing. Surprisingly, trajectory peaks were most common for M1 PTNs (115/146 cells) and rarest for area 5 neurons (12/82 cells). The extent of intrinsic spiking rhythmicity is not therefore simply related to the strength of 20 Hz oscillations in the sensorimotor system. These results suggest that intrinsic rhythmicity is not required for the generation and maintenance of oscillatory activity.

Changes in EMG coherence between long and short thumb abductor muscles during human development

In adults, motoneurone pools of synergistic muscles that act around a common joint share a common presynaptic drive. Common drive can be revealed by both time domain and frequency domain analysis of EMG signals. Analysis in the frequency domain reveals significant coherence in the range 1-45 Hz, with maximal coherence in low (1-12 Hz) and high (16-32 Hz) ranges. The high-frequency range depends on cortical drive to motoneurones and is coherent with cortical oscillations at approximately 20 Hz frequencies. It is of interest to know whether oscillatory drive to human motoneurone pools changes with development. In the present study we examined age-related changes in coherence between rectified surface EMG signals recorded from the short and long thumb abductor muscles during steady isometric contraction obtained while subjects abducted the thumb against a manipulandum. We analysed EMG data from 36 subjects aged between 4 and 14 years, and 11 adult subjects aged between 22 and 59 years. Using the techniques of pooled coherence analysis and the chi(2) difference of coherence test we demonstrate that between the ages of 7 and 9 years, and 12 and 14 years, there are marked increases in the prevalence and magnitude of coherence at frequencies between 11 and 45 Hz. The data from subjects aged 12-14 years were similar to those obtained from adult controls. The most significant differences between younger children and the older age groups were detected at frequencies close to 20 Hz. We believe that these are the first reported results demonstrating significant late maturational changes in the approximately 20 Hz common oscillatory drive to human motoneurone pools.

Inhibitory synaptic transmission governs inspiratory motoneuron synchronization

Neurons within the intact respiratory network produce bursts of action potentials that cause inspiration or expiration. Within inspiratory bursts, activity is synchronized on a shorter timescale to generate clusters of action potentials that occur in a set frequency range and are called synchronous oscillations. We investigated how GABA and glycine modulate synchronous oscillations and respiratory rhythm during postnatal development. We recorded inspiratory activity from hypoglossal nerves using the in vitro rhythmically active mouse medullary slice preparation from P0-P11 mice. Average oscillation frequency increased with postnatal development, from 17 +/- 12 Hz in P0-P6 mice (n = 15) to 38 +/- 7 Hz in P7-P11 mice (n = 37) (P < 0.0001). Bath application of GABAA and GlyR antagonists significantly reduced oscillation power in neonates (P0-P6) and juveniles (P7-P10) and increased peak integrated activity in both age groups. To test whether elevating slice excitability is sufficient to reduce oscillation power, Substance P was bath applied alone. Substance P, although increasing peak integrated activity, had no significant effect on oscillation power. Prolonging the time course of GABAergic synaptic currents with zolpidem decreased the median oscillation frequency in P9-P10 mouse slices. These data demonstrate that oscillation frequency increases with postnatal development and that both GABAergic and glycinergic transmission contribute to synchronization of activity. Further, the time course of synaptic GABAergic currents is a determinant of oscillation frequency.

The effect of transcranial magnetic stimulation and peripheral nerve stimulation on corticomuscular coherence in humans

Cortex and muscle show coupled oscillations in the 15-35 Hz frequency band during voluntary movements. To obtain evidence of the neuronal network responsible for this rhythmicity we investigated the effect of transcranial magnetic stimulation (TMS) and peripheral nerve stimulation on the coupling between eletcroencephalographic (EEG) activity recorded from the scalp over the motor cortex and electromyographic (EMG) activity recorded from the tibialis anterior (TA) muscle in 15 healthy human subjects. TMS over the leg area at intensities between 0.95 and 1.1 x threshold for a motor evoked potential (MEP) in the TA increased corticomuscular coherence in the 15-35 Hz frequency band. This effect lasted on average for 300 ms, but could last up to 600-800 ms in some subjects. Stimulation of motor nerves from the ankle muscles suppressed corticomuscular coherence in the 15-35 Hz frequency range between leg area EEG and TA EMG for a period up to 600-800 ms. In addition, increased coherence around 10 Hz was observed for a period up to 250 ms after the stimulation. Stimulation of motor nerves in the arm and motor nerves from the ankle muscles in the other leg had no effect. The findings indicate that TMS has direct access to the neuronal circuitry in the motor cortex, which generates the corticomuscular coherence. This effect was caused either by direct activation of corticospinal cells or by activation of local neuronal circuitries in the motor cortex. The effects of peripheral nerve stimulation suggest that an alternative rhythm generator may entrain the cortical cells into a lower 10 Hz rhythm and disrupt the 15-35 Hz rhythm.

The effect of carbamazepine on human corticomuscular coherence

EEG recordings from motor cortex show oscillations at approximately 10 and 20 Hz. The 20-Hz oscillations are coherent with contralateral EMG; in most studies those at 10 Hz are not. However, significant 10-Hz coherence has recently been reported in a group of epileptic patients, all of whom were taking the anticonvulsant drug carbamazepine (CBZ). In a double blind study, we investigated the effects of CBZ on corticomuscular coherence in eight healthy human subjects (all male). Subjects performed a precision grip task against an auxotonic load, whilst left sensorimotor EEG and EMGs from five muscles in the right hand and forearm were recorded. CBZ (100 mg) or a placebo was then given orally, and 6 h later subjects were re-tested. One week separated CBZ and placebo experiments in each subject. Coherence averaged across subjects and muscles during the hold phase of the task was maximal at 21 Hz; it increased significantly (P < 0.05, Z-test) by 89% after CBZ administration. This was significantly greater than a much smaller increase following placebo, which itself may reflect an effect of the time of day when experiments were performed. There was no significant approximately 10-Hz coherence either before or after CBZ administration. CBZ did not significantly alter EEG power at either 10 or 20 Hz. Recently, we showed that diazepam markedly increases the power of approximately 20-Hz motor cortical oscillations with little effect on coherence. We show here that CBZ raises coherence without altering EEG power. This pharmacological dissociation may indicate an important role for corticomuscular coherence in motor control.

Post-spike distance-to-threshold trajectories of neurones in monkey motor cortex

A recently developed method permits calculation of the post-spike distance-to-threshold trajectory from an extracellularly recorded spontaneous spike train, using a transform of the interspike interval histogram. We applied this method to 61 single neurones recorded from the primary motor cortex of an awake behaving monkey; 39 cells were antidromically identified as pyramidal tract neurones (PTNs). The cells fell into three categories. Fifty-three trajectories (37 from PTNs) had statistically significant peaks 10-60 ms after the preceding spike. Six neurones (2 PTNs) had non-peaked trajectories which rose exponentially towards threshold. Two cells (both unidentified) had trajectories which declined monotonically away from threshold with increasing post-spike latency. The peaked trajectories were unlikely simply to be an artefact of changing firing rate, which potentially can invalidate this method. Firstly, computer simulations confirmed that the method could accurately re-create both exponential and peaked trajectories, even in the presence of the same rate modulation as seen experimentally. Secondly, the responses of eight cells to weak single pulse intracortical microstimulation (20 microA) through a nearby electrode were measured. For each cell, including representatives of all three trajectory shapes, the modulation of response probability with post-spike latency was consistent with the trajectory computed from the spontaneous discharge. We also demonstrated that cells showed a peaked trajectory during periods with either high or low spontaneous network oscillations, so that the peaks were likely to be generated in part by single cell properties rather than exclusively by network activity. We conclude that many single neurones in motor cortex have an increased probability of firing a spike around 30 ms after the previous action potential. This could act to enhance synchronized oscillatory discharge among populations of cells at functionally relevant frequencies.

EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis

We have developed a toolbox and graphic user interface, EEGLAB, running under the crossplatform MATLAB environment (The Mathworks, Inc.) for processing collections of single-trial and/or averaged EEG data of any number of channels. Available functions include EEG data, channel and event information importing, data visualization (scrolling, scalp map and dipole model plotting, plus multi-trial ERP-image plots), preprocessing (including artifact rejection, filtering, epoch selection, and averaging), independent component analysis (ICA) and time/frequency decompositions including channel and component cross-coherence supported by bootstrap statistical methods based on data resampling. EEGLAB functions are organized into three layers. Top-layer functions allow users to interact with the data through the graphic interface without needing to use MATLAB syntax. Menu options allow users to tune the behavior of EEGLAB to available memory. Middle-layer functions allow users to customize data processing using command history and interactive 'pop' functions. Experienced MATLAB users can use EEGLAB data structures and stand-alone signal processing functions to write custom and/or batch analysis scripts. Extensive function help and tutorial information are included. A 'plug-in' facility allows easy incorporation of new EEG modules into the main menu. EEGLAB is freely available (http://www.sccn.ucsd.edu/eeglab/) under the GNU public license for noncommercial use and open source development, together with sample data, user tutorial and extensive documentation.

Resonant synchronization in heterogeneous networks of inhibitory neurons

Brain rhythms arise through the synchronization of neurons and their entrainment in a regular firing pattern. In this process, networks of reciprocally connected inhibitory neurons are often involved, but what mechanism determines the oscillation frequency is not completely understood. Analytical studies predict that the emerging frequency band is primarily constrained by the decay rate of the unitary IPSC. We observed a new phenomenon of resonant synchronization in computer-simulated networks of inhibitory neurons in which the synaptic current has a delayed onset, reflecting finite spike propagation and synaptic transmission times. At the resonant level of network excitation, all neurons fire synchronously and rhythmically with a period approximately four times the mean delay of the onset of the inhibitory synaptic current. The amplitude and decay time constant of the synaptic current have relatively minor effects on the emerging frequency band. By varying the axonal delay of the inhibitory connections, networks with a realistic synaptic kinetics can be tuned to frequencies from 40 to >200 Hz. This resonance phenomenon arises in heterogeneous networks with, on average, as few as five connections per neuron. We conclude that the delay of the synaptic current is the primary parameter controlling the oscillation frequency of inhibitory networks and propose that delay-induced synchronization is a mechanism for fast brain rhythms that depend on intact inhibitory synaptic transmission.

Synchronous cortical oscillatory activity during motor action

Oscillations of the motor cortex interact with similar activity of the spinal motoneuron pool in the 15-30 Hertz frequency range. Recent observations have demonstrated how this interaction affects the firing of single corticospinal neurons. The interaction, reflected as corticomuscular coherence, occurs for both distal and proximal muscles and it constitutes one connection in a larger web of oscillatory interactions, including several other motor areas in the cortex, thalamus, and cerebellum. New results cast light on the possible functional significance of this interaction. The rhythmic interaction may reveal interesting information in several motor disorders, including essential tremor, Parkinson's disease, myoclonus epilepsy, and mirror movements.

Quantification of the factors that influence discharge correlation in model motor neurons

The purpose of this study was to quantify the influence of intrinsic properties, active dendritic conductances, and background excitation and inhibition on measures of discharge correlation in the time and frequency domains with known levels and patterns of common synaptic input. The study involved a computer simulation of a population of neurons with a range of input resistances (0.54-3.7 MOmega) and surface areas (407,000-712,000 microm(2)). The neurons were simulated with no, moderate, or high levels of active dendritic conductances and were activated with either excitatory input only or excitatory and inhibitory inputs. The patterns of common input, either branched common input or common modulation, were tested with 0, 30, 60, and 90% common input. The results confirm previous findings of an exponential relation between the level of common input and indexes of synchronization; only when the common input comprised >/=60% of the total excitatory input was there a significant effect on discharge correlation. Synchronization was greatest in models that had passive dendrites. Active dendritic conductances caused the discharge rate of the neuron to saturate and decreased motor-unit synchronization. However, the addition of 10% background inhibitory input increased synchronization in these models. In contrast, common rhythmic modulation of inputs at 24 Hz usually decreased synchronization. Significant coherence at the modulated frequency occurred in the commonly modulated neurons when >/=60% of the inputs were modulated. Furthermore, active dendritic conductances decreased coherence. Branched common input caused high levels of coherence across a broad spectrum and when combined with active dendritic conductances caused significant frequency peaks in the 30- to 50-Hz band. In conclusion, the level of inhibitory input and active dendritic conductances interact with the amount of common input to determine time- and frequency-domain discharge correlation.

Synchronization in monkey motor cortex during a precision grip task. II. effect of oscillatory activity on corticospinal output

Recordings from primary motor cortex (M1) during periods of steady contraction show oscillatory activity; these oscillations are coherent with the activity of contralateral muscles. We investigated synchronization of corticospinal output neurons with the oscillations, which could provide the pathway for their transmission to the spinal motoneurons. One hundred seventy-six antidromically identified pyramidal tract neurons (PTNs) were recorded from M1 in three macaque monkeys trained to perform a precision grip task. Local field potentials (LFP) were simultaneously recorded. All analysis was confined to the hold period of the task, where our previous work has shown that there is the strongest oscillatory activity. Coherence was calculated between LFP and PTN discharge. Significant coherence was seen in three bands, with frequencies of 10-14, 17-31, and 34-44 Hz. Coherence values were low, with the majority of PTN-LFP coherences having a peak lower than 0.05. The phase of coherence was approximately -pi/2 radians for each band (with LFP polarity defined as negative upward), although there was some dispersion of phase across the population of PTNs. Coherence was also calculated between pairs of PTNs that had been simultaneously recorded. Where there was significant coherence, it was also generally smaller than 0.05. The phase of PTN-PTN coherence clustered around zero radians. A computer model was constructed to assist the interpretation of the experimental results. It simulated an integrate-and-fire neuron responding to synaptic inputs. A fraction of the synaptic inputs was synchronized with a simulated LFP; the remainder were uncorrelated with it. The model showed that coherence between the LFP and the output spike train considerably underestimated the fraction of synchronized inputs. Additionally, for a given fraction of synchronized inputs, coherence was smaller for high- compared with low-frequency bins. Cell discharge rate also influenced the spike-LFP coherence: coherence was higher for simulations in which the cell discharged at a faster rate. Thus although levels of PTN-LFP coherence seen experimentally were low, a considerable proportion of the input to the PTN must be synchronized with the global oscillatory activity recorded by the LFP. The low LFP-PTN coherences do however indicate that cortical oscillations are transmitted with only low fidelity in the discharge of a single PTN. Using further computer simulations, it was demonstrated that a small population of PTNs could encode the cortical oscillatory signal effectively, since the action of averaging across the population improves the signal:noise ratio. The oscillations will therefore be effectively transmitted to spinal motoneurons, and this has important consequences for the possible role of oscillations in motor control of the hand.

The effect of diazepam on motor cortical oscillations and corticomuscular coherence studied in man

EEG recordings from sensorimotor cortex show oscillations around 10 and 20 Hz. These modulate with task performance, and are strongest during periods of steady contraction. The 20 Hz oscillations are coherent with contralateral EMG. Computer modelling suggests that oscillations arising within the cortex may be especially dependent on inhibitory systems. The benzodiazepine diazepam enhances the size of GABA(A) IPSPs; its effects are reversed by the antagonist flumazenil. We tested the effect of these drugs on spectral measures of EEG and EMG, whilst eight healthy human subjects performed a precision grip task containing both holding and movement phases. Either an auxotonic or isometric load was used. EEG changes following electrical stimulation of the contralateral median nerve were also assessed. The EEG power showed similar changes in all task/stimulation protocols used. Power around 20 Hz doubled at the highest dose of diazepam used (5 mg), and returned to control levels following flumazenil. EEG power at 10 Hz was by contrast little altered. The peak frequency of EEG power in both bands was not changed by diazepam. Corticomuscular coherence at ca 20 Hz was reduced following diazepam injection, but the magnitude of this effect was small (mean coherence during steady holding in the auxotonic task was 0.062 in control recordings, 0.051 after 2.5 mg and 5 mg doses of diazepam). These results imply that 20 Hz oscillations in the sensorimotor cortex are at least partially produced by local cortical circuits reliant on GABA(A)-mediated intracortical inhibition, whereas 10 Hz rhythms arise by a different mechanism. Rhythms generated during different tasks, or following nerve stimulation, are likely to arise from similar mechanisms. By examining the formulae used to calculate coherence, we show that if cortical oscillations are simply transmitted to the periphery, corticomuscular coherence should increase in parallel with the ratio of EEG to EMG power. The relative constancy of coherence even when the amplitude of cortical oscillations is perturbed suggests that corticomuscular coherence itself may have a functional role in motor control.

Inhibitory feedback required for network oscillatory responses to communication but not prey stimuli

Stimulus-induced oscillations occur in visual, olfactory and somatosensory systems. Several experimental and theoretical studies have shown how such oscillations can be generated by inhibitory connections between neurons. But the effects of realistic spatiotemporal sensory input on oscillatory network dynamics and the overall functional roles of such oscillations in sensory processing are poorly understood. Weakly electric fish must detect electric field modulations produced by both prey (spatially localized) and communication (spatially diffuse) signals. Here we show, through in vivo recordings, that sensory pyramidal neurons in these animals produce an oscillatory response to communication-like stimuli, but not to prey-like stimuli. On the basis of well-characterized circuitry, we construct a network model of pyramidal neurons that predicts that diffuse delayed inhibitory feedback is required to achieve oscillatory behaviour only in response to communication-like stimuli. This prediction is experimentally verified by reversible blockade of feedback inhibition that removes oscillatory behaviour in the presence of communication-like stimuli. Our results show that a sensory system can use inhibitory feedback as a mechanism to 'toggle' between oscillatory and non-oscillatory firing states, each associated with a naturalistic stimulus.

The effects of subthreshold 1 Hz repetitive TMS on cortico-cortical and interhemispheric coherence

OBJECTIVES

Repetitive transcranial magnetic stimulation (rTMS) shows promise as a treatment for various movement and psychiatric disorders. Just how rTMS may have persistent effects on cortical function remains unclear. We hypothesised that it may act by modulating cortico-cortical and interhemispheric connectivity. To this end we assessed cortico-cortical and interhemispheric coherence before and after low frequency, subthreshold rTMS of the left motor cortex.

METHODS

Fifteen healthy subjects received one train (1Hz, 90% of active motor threshold, 1500 stimuli) of rTMS to the left motor hand area. Spectral power and coherence estimates were calculated between different electroencephalogram (EEG) signals at rest and while muscles of the distal upper limb were tonically contracted.

RESULTS

rTMS over the left motor hand area caused a significant increase in ipsilateral EEG-EEG coherence and in the interhemispheric coherence between motor areas in the alpha band. The effects of rTMS lasted up to 25 min post-stimulation. There was no significant change in EEG-EEG coherence over the hemisphere contralateral to stimulation.

CONCLUSIONS

Low frequency, subthreshold rTMS of the motor cortex increases ipsilateral cortico-cortical and interhemispheric coherence in the alpha band. This may, in part, mediate the inhibitory effects of low frequency rTMS.

Rhythm generation in monkey motor cortex explored using pyramidal tract stimulation

We investigated whether stimulation of the pyramidal tract (PT) could reset the phase of 15-30 Hz beta oscillations observed in the macaque motor cortex. We recorded local field potentials (LFPs) and multiple single-unit activity from two conscious macaque monkeys performing a precision grip task. EMG activity was also recorded from the second animal. Single PT stimuli were delivered during the hold period of the task, when oscillations in the LFP were most prominent. Stimulus-triggered averaging of the LFP showed a phase-locked oscillatory response to PT stimulation. Frequency domain analysis revealed two components within the response: a 15-30 Hz component, which represented resetting of on-going beta rhythms, and a lower frequency 10 Hz response. Only the higher frequency could be observed in the EMG activity, at stronger stimulus intensities than were required for resetting the cortical rhythm. Stimulation of the PT during movement elicited a greatly reduced oscillatory response. Analysis of single-unit discharge confirmed that PT stimulation was capable of resetting periodic activity in motor cortex. The firing patterns of pyramidal tract neurones (PTNs) and unidentified neurones exhibited successive cycles of suppression and facilitation, time locked to the stimulus. We conclude that PTN activity directly influences the generation of the 15-30 Hz rhythm. These PTNs facilitate EMG activity in upper limb muscles, contributing to corticomuscular coherence at this same frequency. Since the earliest oscillatory effect observed following stimulation was a suppression of firing, we speculate that inhibitory feedback may be the key mechanism generating such oscillations in the motor cortex.

Rapid signaling at inhibitory synapses in a dentate gyrus interneuron network

Mutual synaptic interactions between GABAergic interneurons are thought to be of critical importance for the generation of network oscillations and for temporal encoding of information in the hippocampus. However, the functional properties of synaptic transmission between hippocampal interneurons are largely unknown. We have made paired recordings from basket cells (BCs) in the dentate gyrus of rat hippocampal slices, followed by correlated light and electron microscopical analysis. Unitary GABA(A) receptor-mediated IPSCs at BC-BC synapses recorded at the soma showed a fast rise and decay, with a mean decay time constant of 2.5 +/- 0.2 msec (32 degrees C). Synaptic transmission at BC-BC synapses showed paired-pulse depression (PPD) (32 +/- 5% for 10 msec interpulse intervals) and multiple-pulse depression during repetitive stimulation. Detailed passive cable model simulations based on somatodendritic morphology and localization of synaptic contacts further indicated that the conductance change at the postsynaptic site was even faster, decaying with a mean time constant of 1.8 +/- 0.6 msec. Sequential triple recordings revealed that the decay time course of IPSCs at BC-BC synapses was approximately twofold faster than that at BC-granule cell synapses, whereas the extent of PPD was comparable. To examine the consequences of the fast postsynaptic conductance change for the generation of oscillatory activity, we developed a computational model of an interneuron network. The model showed robust oscillations at frequencies >60 Hz if the excitatory drive was sufficiently large. Thus the fast conductance change at interneuron-interneuron synapses may promote the generation of high-frequency oscillations observed in the dentate gyrus in vivo.

Precise burst synchrony in the superior colliculus of the awake cat during moving stimulus presentation

This study aimed to characterize the synchrony that occurs between cell discharges in the superior colliculus of the awake cat. We trained cats to perform a visual fixation in the presence of a visual moving stimulus and then recorded 686 pairs of neighboring cells in the superior colliculus during task performance. A new method to assess the significance of precise discharge synchronization is described, which permits analysis of nonstationary data. Of 181 pairs with sufficient data for quantitative analysis, 125 showed a cross-correlation histogram (CCH) with features assessed as significant using this approach. CCHs frequently showed an isolated central peak (41 of 125) or a peak flanked by one or two troughs (68 of 125), and in a few cases an oscillatory pattern of approximately 65 Hz (16 of 125). This is in contrast to the oscillation frequency reported for the visual cortex and shows that oscillations in the superior colliculus probably arise from a cortex-independent mechanism. Our method also permits direct quantification of the correlation shift predictors, assessing precise time locking of spikes to the stimulus. Only 1 of 125 cross-correlation shift predictors had a significant central peak, meaning that most of the CCH features were not related to cell discharges time-locked to the stimulus presentation. Further investigation using a burst-jittering method showed that synchrony in the superior colliculus is attributable to precise synchronization of short bursts of spikes. Such synchrony could be related to the network dynamics and the common inhibitory feedback from local interneurons, which would act as temporal selectors of the cells with greatest or fastest response.

Statistical signs of common inhibitory feedback with delay

Cross-correlation histograms (CCHs) sometimes exhibit an isolated central peak flanked by two troughs. What can cause this pattern? The absence of CCH satellite peak makes an oscillatory common input doubtful. It is here shown using a simple counting model that a common inhibitory feedback with delay can account for this pattern.

Evidence for glutamatergic tectotectal neurons in the cat superior colliculus: a comparison with GABAergic tectotectal neurons

The tectotectal commissural pathway is commonly regarded as responsible for the reciprocal inhibition that takes place between the two superior colliculi (SC). Although this hypothesis has received strong support from electrophysiological studies, more recent investigations have suggested that some collicular cells, e.g. fixation neurons, may establish excitatory connections with cells in the contralateral SC through the collicular commissure. The goal of the present study was to seek immunohistochemical evidence for glutamatergic tectotectal cells in the cat SC by using a double-labelling technique. Tectotectal cells were retrogradely labelled with wheat germ agglutinin (WGA) -horseradish peroxidase (HRP) coupled to colloidal gold injected in the contralateral SC, and neurons containing glutamate or gamma-aminobutyric acid (GABA) were then identified with immunohistochemical techniques. The present study provides evidence that, in the cat SC, equal numbers of tectotectal cells are immunopositive to glutamate and GABA, suggesting that the tectotectal pathway may consist of two distinct functional components. The finding that an equal number of tectotectal cells are GABAergic and glutamatergic is somewhat surprising as electrophysiological studies have invariantly indicated that the inhibitory component of the tectotectal projection predominates. Another striking feature of the GABAergic and glutamatergic tectotectal cell populations is their identical topographic distribution in the SC. These results suggest that not only cells in the rostral fixation zone establish excitatory connections with the contralateral SC. Tectotectal projections could be potentially important to shape the spatial pattern of saccade-related activity that may occur simultaneously in the two SC during vertical and oblique orienting movements.

An accurate measure of the instantaneous discharge probability, with application to unitary joint-even analysis

We present an estimate for the instantaneous discharge probability of a neurone, based on single-trial spike-train analysis. By detecting points where the neurone abruptly changes its firing rate and treating them specially, the method is able to achieve smooth estimates yet avoid the blurring of significant changes. This estimate of instantaneous discharge probability is then applied to the method of unitary event analysis. We show that the unitary event analysis as originally conceived is highly sensitive to firing-rate nonstationarities and covariations, but that it can be considerably improved if calculations of statistical significance use an instantaneous discharge probability instead of a firing-rate estimate based on averaging across multiple trials.

Multiple single unit recording in the cortex of monkeys using independently moveable microelectrodes

Simultaneous recording from multiple single neurones presents many technical difficulties. However, obtaining such data has many advantages, which make it highly worthwhile to overcome the technical problems. This report describes methods which we have developed to permit recordings in awake behaving monkeys using the 'Eckhorn' 16 electrode microdrive. Structural magnetic resonance images are collected to guide electrode placement. Head fixation is achieved using a specially designed headpiece, modified for the multiple electrode approach, and access to the cortex is provided via a novel recording chamber. Growth of scar tissue over the exposed dura mater is reduced using an anti-mitotic compound. Control of the microdrive is achieved by a computerised system which permits several experimenters to move different electrodes simultaneously, considerably reducing the load on an individual operator. Neurones are identified as pyramidal tract neurones by antidromic stimulation through chronically implanted electrodes; stimulus control is integrated into the computerised system. Finally, analysis of multiple single unit recordings requires accurate methods to correct for non-stationarity in unit firing. A novel technique for such correction is discussed.