Over the last decade, the search for a method able to reliably predict seizures hours in advance has been largely replaced by a more realistic goal of very early detection of seizure onset which would allow therapeutic or warning devices to be triggered prior to the onset of disabling clinical symptoms. We explore in this article the steps along the pathway from data acquisition to closed loop applications that can and should be considered to design the most efficient early seizure detection. Microelectrodes, high-frequency oscillations, high sampling rate, high-density arrays, and modern analysis techniques are all elements of the recording and detection process that in combination with modeling studies can provide new insights into the dynamics of seizure onsets. Each of these step needs to be considered if one wants to implement improved detection devices that will favorably impact the quality of life of patients.

Intracranial EEG studies in humans have shown that functional brain activation in a variety of functional-anatomic domains of human cortex is associated with an increase in power at a broad range of high gamma (> 60 Hz) frequencies. Although these electrophysiological responses are highly specific for the location and timing of cortical processing and in animal recordings are highly correlated with increased population firing rates, there has been little direct empirical evidence for causal interactions between different recording sites at high gamma frequencies. Such causal interactions are hypothesized to occur during cognitive tasks that activate multiple brain regions. To determine whether such causal interactions occur at high gamma frequencies and to investigate their functional significance, we used event-related causality (ERC) analysis to estimate the dynamics, directionality, and magnitude of event-related causal interactions using subdural electrocorticography (ECoG) recorded during two word production tasks: picture naming and auditory word repetition. A clinical subject who had normal hearing but was skilled in American Signed Language (ASL) provided a unique opportunity to test our hypothesis with reference to a predictable pattern of causal interactions, i.e. that language cortex interacts with different areas of sensorimotor cortex during spoken vs. signed responses. Our ERC analyses confirmed this prediction. During word production with spoken responses, perisylvian language sites had prominent causal interactions with mouth/tongue areas of motor cortex, and when responses were gestured in sign language, the most prominent interactions involved hand and arm areas of motor cortex. Furthermore, we found that the sites from which the most numerous and prominent causal interactions originated, i.e. sites with a pattern of ERC “divergence”, were also sites where high gamma power increases were most prominent and where electrocortical stimulation mapping interfered with word production. These findings suggest that the number, strength and directionality of event-related causal interactions may help identify network nodes that are not only activated by a task but are critical to its performance.

Objectives

Partial seizures are often believed to be associated with EEG signals of low complexity because seizures are associated with increased neural network synchrony. The investigations reported here provide an assessment of the signal complexity of epileptic seizure onsets using newly developed quantitative measures.

Methods

Using the Gabor atom density (GAD) measure of signal complexity, 339 partial seizures in 45 patients with intracranial electrode arrays were analyzed. Segmentation procedures were applied to determine the timing and amplitude of GAD changes relative to the electrographic onset of the seizure.

Results

330 out of 339 seizures have significant complexity level changes, with 319 (97%) having an increase in complexity. GAD increases occur within seconds of the onset of the partial seizure but are not observed in channels remote from the focus. The complexity increase is similar for seizures from mesial temporal origin, neocortical temporal and extra-temporal origin.

Conclusions

Partial onset seizures are associated with early increases in signal complexity as measured by GAD. This increase is independent of the location of the seizure focus.

Significance

Despite the often predominant rhythmic activity that characterizes onset and early evolution of epileptic seizures, partial seizure onset is associated with an early increase in complexity. These changes are common to partial seizures originating from different brain regions, indicating a similar seizure dynamic.

Multilingual patients pose a unique challenge when planning epilepsy surgery near language cortex because the cortical representations of each language may be distinct. These distinctions may not be evident with routine electrocortical stimulation mapping (ESM). Electrocorticography (ECoG) has recently been used to detect task-related spectral perturbations associated with functional brain activation. We hypothesized that using broadband high gamma augmentation (HGA, 60–150 Hz) as an index of cortical activation, ECoG would complement ESM in discriminating the cortical representations of first (L1) and second (L2) languages. We studied four adult patients for whom English was a second language, in whom subdural electrodes (a total of 358) were implanted to guide epilepsy surgery. Patients underwent ECoG recordings and ESM while performing the same visual object naming task in L1 and L2. In three of four patients, ECoG found sites activated during naming in one language but not the other. These language-specific sites were not identified using ESM. In addition, ECoG HGA was observed at more sites during L2 versus L1 naming in two patients, suggesting that L2 processing required additional cortical resources compared to L1 processing in these individuals. Post-operative language deficits were identified in three patients (one in L2 only). These deficits were predicted by ECoG spectral mapping but not by ESM. These results suggest that pre-surgical mapping should include evaluation of all utilized languages to avoid post-operative functional deficits. Finally, this study suggests that ECoG spectral mapping may potentially complement the results of ESM of language.

Multichannel intracranial recordings are used increasingly to study the functional organization of human cortex. Intracranial recordings of event-related activity, or electrocorticography (ECoG), are based on high density electrode arrays implanted directly over cortex, combining good temporal and spatial resolution. Developing appropriate statistical methods for analyzing event-related responses in these high dimensional ECoG datasets remains a major challenge for clinical and systems neuroscience. We present a novel methodological framework that combines complementary, existing methods adapted for statistical analysis of auditory event-related responses in multichannel ECoG recordings. This analytic framework integrates single-channel (time-domain, time–frequency) and multichannel analyses of event-related ECoG activity to determine statistically significant evoked responses, induced spectral responses, and effective (causal) connectivity. Implementation of this quantitative approach is illustrated using multichannel ECoG data from recent studies of auditory processing in patients with epilepsy. Methods described include a time–frequency matching pursuit algorithm adapted for modeling brief, transient cortical spectral responses to sound, and a recently developed method for estimating effective connectivity using multivariate autoregressive modeling to measure brief event-related changes in multichannel functional interactions. A semi-automated spatial normalization method for comparing intracranial electrode locations across patients is also described. The individual methods presented are published and readily accessible. We discuss the benefits of integrating multiple complementary methods in a unified and comprehensive quantitative approach. Methodological considerations in the analysis of multichannel ECoG data, including corrections for multiple comparisons are discussed, as well as remaining challenges in the development of new statistical approaches.

Summary

Purpose

A neural network simulation with realistic cortical architecture has been used to study synchronized bursting as a seizure representation. This model has the property that bursting epochs arise and cease spontaneously, and bursting epochs can be induced by external stimulation. We have used this simulation to study the time-frequency properties of the evolving bursting activity, as well as effects due to network stimulation.

Methods

The model represents a cortical region of 1.6 mm × 1.6 mm, and includes seven neuron classes organized by cortical layer, inhibitory or excitatory properties, and electrophysiological characteristics. There are a total of 65, 536 modeled single compartment neurons that operate according to a version of Hodgkin-Huxley dynamics. The intercellular wiring is based on histological studies and our previous modeling efforts.

Results

The bursting phase is characterized by a flat frequency spectrum. Stimulation pulses are applied to this modeled network, with an electric field provided by a 1 mm radius circular electrode represented mathematically in the simulation. A phase dependence to the post-stimulation quiescence is demonstrated, with local relative maxima in efficacy occurring before or during the network depolarization phase in the underlying activity. Brief periods of network insensitivity to stimulation are also demonstrated. The phase dependence was irregular and did not reach statistical significance when averaged over the full 2.5 seconds of simulated bursting investigated. This result provides comparison with previous in vivo studies which have also demonstrated increased efficacy of stimulation when pulses are applied at the peak of the local field potential during cortical afterdischarges. The network bursting is synchronous when comparing the different neuron classes represented up to an uncertainty of 10 msec. Studies performed with an excitatory chandelier cell component demonstrated increased synchronous bursting in the model, as predicted from experimental work.

Conclusions

This large scale multi-neuron neural network simulation reproduces many aspects of evolving cortical bursting behaviour as well as the timing-dependent effects of electrical stimulation on that bursting.

Recent studies using electrocorticographic (ECoG) recordings in humans have shown that functional activation of cortex is associated with an increase in power in the high-gamma frequency range (∼60–200 Hz). Here we investigate the neural correlates of this high-gamma activity in local field potential (LFP). Single units and LFP were recorded with microelectrodes from the hand region of macaque SII cortex while vibrotactile stimuli of varying intensities were presented to the hand. We found that high-gamma power in the LFP was strongly correlated with the average firing rate recorded by the microelectrodes, both temporally and on a trial-by-trial basis. In comparison, the correlation between firing rate and low-gamma power (40–80 Hz) was much smaller. In order to explore the potential effects of neuronal firing on ECoG, we developed a model to estimate ECoG power generated by different firing patterns of the underlying cortical population and studied how ECoG power varies with changes in firing rate versus the degree of synchronous firing between neurons in the population. Both an increase in firing rate and neuronal synchrony increased high-gamma power in the simulated ECoG data. However, ECoG high-gamma activity was much more sensitive to increases in neuronal synchrony than firing rate.

Neuronal oscillations in the gamma frequency range have been reported in many cortical areas, but the role they play in cortical processing remains unclear. We tested a recently proposed hypothesis that the intensity of sensory input is coded in the timing of action potentials relative to the phase of gamma oscillations, thus converting amplitude information to a temporal code. We recorded spikes and local field potential (LFP) from secondary somatosensory (SII) cortex in awake monkeys while presenting a vibratory stimulus at different amplitudes. We developed a novel technique based on matching pursuit to study the interaction between the highly transient gamma oscillations and spikes with high time-frequency resolution. We found that spikes were weakly coupled to LFP oscillations in the gamma frequency range (40−80 Hz), and strongly coupled to oscillations in higher gamma frequencies. However, the phase relationship of neither low-gamma nor high-gamma oscillations changed with stimulus intensity, even with a ten-fold increase. We conclude that, in SII, gamma oscillations are synchronized with spikes, but their phase does not vary with stimulus intensity. Furthermore, high-gamma oscillations (>60 Hz) appear to be closely linked to the occurrence of action potentials, suggesting that LFP high-gamma power could be a sensitive index of the population firing rate near the microelectrode.