Identification of epileptogenic foci from causal analysis of ECoG interictal spike activity

Measurable transitions during seizures in intracranial EEG: A stereoelectroencephalography and SPECT study

OBJECTIVE

Ictal Single Photon Emission Computed Tomography (SPECT) and stereo-electroencephalography (SEEG) are diagnostic techniques used for the management of patients with drug-resistant focal epilepsies. While hyperperfusion patterns in ictal SPECT studies reveal seizure onset and propagation pathways, the role of ictal hypoperfusion remains poorly understood. The goal of this study was to systematically characterize the spatio-temporal information flow dynamics between differently perfused brain regions using stereo-EEG recordings.

METHODS

We identified seizure-free patients after resective epilepsy surgery who had prior ictal SPECT and SEEG investigations. We estimated directional connectivity between the epileptogenic-zone (EZ), non-resected areas of hyperperfusion, hypoperfusion, and baseline perfusion during the interictal, preictal, ictal, and postictal periods.

RESULTS

Compared to the background, we noted significant information flow (1) during the preictal period from the EZ to the baseline and hyperperfused regions, (2) during the ictal onset from the EZ to all three regions, and (3) during the period of seizure evolution from the area of hypoperfusion to all three regions.

CONCLUSIONS

Hypoperfused brain regions were found to indirectly interact with the EZ during the ictal period.

SIGNIFICANCE

Our unique study, combining intracranial electrophysiology and perfusion imaging, presents compelling evidence of dynamic changes in directional connectivity between brain regions during the transition from interictal to ictal states.

Causal Connectivity Network Analysis of Ictal Electrocorticogram With Temporal Lobe Epilepsy Based on Dynamic Phase Transfer Entropy

Temporallobe epilepsy (TLE) has been conceptualized as a brain network disease, which generates brain connectivity dynamics within and beyond the temporal lobe structures in seizures. The hippocampus is a representative epileptogenic focus in TLE. Understanding the causal connectivity in terms of brain network during seizures is crucial in revealing the triggering mechanism of epileptic seizures originating from the hippocampus (HPC) spread to the lateral temporal cortex (LTC) by ictal electrocorticogram (ECoG), particularly in high-frequency oscillations (HFOs) bands. In this study, we proposed the unified-epoch dynamic causality analysis method to investigate the causal influence dynamics between two brain regions (HPC and LTC) at interictal and ictal phases in the frequency range of 1-500 Hz by introducing the phase transfer entropy (PTE) out/in-ratio and sliding window. We also proposed PTE-based machine learning algorithms to identify epileptogenic zone (EZ). Nine patients with a total of 26 seizures were included in this study. We hypothesized that: 1) HPC is the focus with the stronger causal connectivity than that in LTC in the ictal state at gamma and HFOs bands. 2) Causal connectivity in the ictal phase shows significant changes compared to that in the interictal phase. 3) The PTE out/in-ratio in the HFOs band can identify the EZ with the best prediction performance.

Spike propagation mapping reveals effective connectivity and predicts surgical outcome in epilepsy

Neurosurgical intervention is the best available treatment for selected patients with drug resistant epilepsy. For these patients, surgical planning requires biomarkers that delineate the epileptogenic zone, the brain area that is indispensable for the generation of seizures. Interictal spikes recorded with electrophysiological techniques are considered key biomarkers of epilepsy. Yet, they lack specificity, mostly because they propagate across brain areas forming networks. Understanding the relationship between interictal spike propagation and functional connections among the involved brain areas may help develop novel biomarkers that can delineate the epileptogenic zone with high precision. Here, we reveal the relationship between spike propagation and effective connectivity among onset and areas of spread and assess the prognostic value of resecting these areas. We analysed intracranial EEG data from 43 children with drug resistant epilepsy who underwent invasive monitoring for neurosurgical planning. Using electric source imaging, we mapped spike propagation in the source domain and identified three zones: onset, early-spread and late-spread. For each zone, we calculated the overlap and distance from surgical resection. We then estimated a virtual sensor for each zone and the direction of information flow among them via Granger causality. Finally, we compared the prognostic value of resecting these zones, the clinically-defined seizure onset zone and the spike onset on intracranial EEG channels by estimating their overlap with resection. We observed a spike propagation in source space for 37 patients with a median duration of 95 ms (interquartile range: 34-206), a spatial displacement of 14 cm (7.5-22 cm) and a velocity of 0.5 m/s (0.3-0.8 m/s). In patients with good surgical outcome (25 patients, Engel I), the onset had higher overlap with resection [96% (40-100%)] than early-spread [86% (34-100%), P = 0.01] and late-spread [59% (12-100%), P = 0.002], and it was also closer to resection than late-spread [5 mm versus 9 mm, P = 0.007]. We found an information flow from onset to early-spread in 66% of patients with good outcomes, and from early-spread to onset in 50% of patients with poor outcome. Finally, resection of spike onset, but not area of spike spread or the seizure onset zone, predicted outcome with positive predictive value of 79% and negative predictive value of 56% (P = 0.04). Spatiotemporal mapping of spike propagation reveals information flow from onset to areas of spread in epilepsy brain. Surgical resection of the spike onset disrupts the epileptogenic network and may render patients with drug resistant epilepsy seizure-free without having to wait for a seizure to occur during intracranial monitoring.

An application of dynamical directed connectivity of ictal intracranial EEG recordings in seizure onset zone localization

BACKGROUND

Identification of the seizure onset zone (SOZ) is a challenging task in epilepsy surgery. Patients with epilepsy have an altered brain network, allowing connectivity-based analyses to have a great potential in SOZ identification. We investigated a dynamical directed connectivity analysis utilizing ictal intracranial electroencephalographic (iEEG) recordings and proposed an algorithm for SOZ identification based on grouping iEEG contacts.

NEW METHODS

Granger Causality was used for directed connectivity analysis in this study. The intracranial contacts were grouped into visually detected contacts (VDCs), which were identified as SOZ by epileptologists, and non-resected contacts (NRCs). The intragroup and intergroup directed connectivity for VDCs and NRCs were calculated around seizure onset. We then proposed an algorithm for SOZ identification based on the cross-correlation of intragroup outflow and inflow of SOZ candidate contacts.

RESULTS

Our results revealed that the intragroup connectivity of VDCs (VDC→VDC) was significantly larger than the intragroup connectivity of NRCs (NRC→NRC) and the intergroup connectivity between NRCs and VDCs (NRC→VDC) around seizure onset. We found that the proposed algorithm had 90.1 % accuracy for SOZ identification in the seizure-free patients.

COMPARISON WITH EXISTING METHODS

The existing connectivity-based methods for SOZ identification often use either outflow or inflow. In this study, SOZ contacts were identified by integrating outflow and inflow based on the cross correlation between these two measures.

CONCLUSIONS

The proposed group-based dynamical connectivity analysis in this study can aid our understanding of underlying seizure network and may be used to assist in identifying the SOZ contacts before epilepsy surgery.

Interictal SEEG Resting-State Connectivity Localizes the Seizure Onset Zone and Predicts Seizure Outcome

Localization of epileptogenic zone currently requires prolonged intracranial recordings to capture seizure, which may take days to weeks. The authors developed a novel method to identify the seizure onset zone (SOZ) and predict seizure outcome using short-time resting-state stereotacticelectroencephalography (SEEG) data. In a cohort of 27 drug-resistant epilepsy patients, the authors estimated the information flow via directional connectivity and inferred the excitation-inhibition ratio from the 1/f power slope. They hypothesized that the antagonism of information flow at multiple frequencies between SOZ and non-SOZ underlying the relatively stable epilepsy resting state could be related to the disrupted excitation-inhibition balance. They found flatter 1/f power slope in non-SOZ regions compared to the SOZ, with dominant information flow from non-SOZ to SOZ regions. Greater differences in resting-state information flow between SOZ and non-SOZ regions are associated with favorable seizure outcome. By integrating a balanced random forest model with resting-state connectivity, their method localized the SOZ with an accuracy of 88% and predicted the seizure outcome with an accuracy of 92% using clinically determined SOZ. Overall, this study suggests that brief resting-state SEEG data can significantly facilitate the identification of SOZ and may eventually predict seizure outcomes without requiring long-term ictal recordings.

Microscale dynamics of electrophysiological markers of epilepsy

OBJECTIVE

Interictal discharges (IIDs) and high frequency oscillations (HFOs) are established neurophysiologic biomarkers of epilepsy, while microseizures are less well studied. We used custom poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS) microelectrodes to better understand these markers' microscale spatial dynamics.

METHODS

Electrodes with spatial resolution down to 50 µm were used to record intraoperatively in 30 subjects. IIDs' degree of spread and spatiotemporal paths were generated by peak-tracking followed by clustering. Repeating HFO patterns were delineated by clustering similar time windows. Multi-unit activity (MUA) was analyzed in relation to IID and HFO timing.

RESULTS

We detected IIDs encompassing the entire array in 93% of subjects, while localized IIDs, observed across < 50% of channels, were seen in 53%. IIDs traveled along specific paths. HFOs appeared in small, repeated spatiotemporal patterns. Finally, we identified microseizure events that spanned 50-100 µm. HFOs covaried with MUA, but not with IIDs.

CONCLUSIONS

Overall, these data suggest that irritable cortex micro-domains may form part of an underlying pathologic architecture which could contribute to the seizure network.

SIGNIFICANCE

These results, supporting the possibility that epileptogenic cortex comprises a mosaic of irritable domains, suggests that microscale approaches might be an important perspective in devising novel seizure control therapies.

Interictal Epileptiform Discharge Dynamics in Peri-sylvian Polymicrogyria Using EEG-fMRI

Polymicrogyria (PMG) is a common malformation of cortical development associated with a higher susceptibility to epileptic seizures. Seizures secondary to PMG are characterized by difficult-to-localize cerebral sources due to the complex and widespread lesion structure. Tracing the dynamics of interictal epileptiform discharges (IEDs) in patients with epilepsy has been shown to reveal the location of epileptic activity sources, crucial for successful treatment in cases of focal drug-resistant epilepsy. In this case series IED dynamics were evaluated with simultaneous EEG-fMRI recordings in four patients with unilateral peri-sylvian polymicrogyria (PSPMG) by tracking BOLD activations over time: before, during and following IED appearance on scalp EEG. In all cases, focal BOLD activations within the lesion itself preceded the activity associated with the time of IED appearance on EEG, which showed stronger and more widespread activations. We therefore propose that early hemodynamic activity corresponding to IEDs may hold important localizing information potentially leading to the cerebral sources of epileptic activity. IEDs are suggested to develop within a small area in the PSPMG lesion with structural properties obscuring the appearance of their electric field on the scalp and only later engage widespread structures which allow the production of large currents which are recognized as IEDs on EEG.

Ictal gamma-band interactions localize ictogenic nodes of the epileptic network in focal cortical dysplasia

OBJECTIVE

Epilepsy surgery fails in > 30% of patients with focal cortical dysplasia (FCD). The seizure persistence after surgery can be attributed to the inability to precisely localize the tissue with an endogenous potential to generate seizures. In this study, we aimed to identify the critical components of the epileptic network that were actively involved in seizure genesis.

METHODS

The directed transfer function was applied to intracranial EEG recordings and the effective connectivity was determined with a high temporal and frequency resolution. Pre-ictal network properties were compared with ictal epochs to identify regions actively generating ictal activity and discriminate them from the areas of propagation.

RESULTS

Analysis of 276 seizures from 30 patients revealed the existence of a seizure-related network reconfiguration in the gamma-band (25-170 Hz; p < 0.005) - ictogenic nodes. Unlike seizure onset zone, resecting the majority of ictogenic nodes correlated with favorable outcomes (p < 0.012).

CONCLUSION

The prerequisite to successful epilepsy surgery is the accurate identification of brain areas from which seizures arise. We show that in FCD-related epilepsy, gamma-band network markers can reliably identify and distinguish ictogenic areas in macroelectrode recordings, improve intracranial EEG interpretation and better delineate the epileptogenic zone.

SIGNIFICANCE

Ictogenic nodes localize the critical parts of the epileptogenic tissue and increase the diagnostic yield of intracranial evaluation.

Model-based design for seizure control by stimulation

OBJECTIVE

Current brain stimulation paradigms are largely empirical rather than theoretical. An opportunity exists to improve upon their modest effectiveness in closed-loop control strategies with the development of theoretically grounded, model-based designs.

APPROACH

Inspired by this need, here we couple experimental data and mathematical modeling with a control-theoretic strategy for seizure termination. We begin by exercising a dynamical systems approach to model seizures (n = 94) recorded using intracranial EEG (iEEG) from 21 patients with medication-resistant, localization-related epilepsy.

MAIN RESULTS

Although each patient's seizures displayed unique spatial and temporal patterns, their evolution can be parsimoniously characterized by the same model form. Idiosyncracies of the model can inform individualized intervention strategies, specifically in iEEG samples with well-localized seizure onset zones. Temporal fluctuations in the spatial profiles of the oscillatory modes show that seizure onset marks a transition into a regime in which the underlying system supports prolonged rhythmic and focal activity. Based on these observations, we propose a control-theoretic strategy that aims to stabilize ictal activity using static output feedback for linear time-invariant switching systems. Finally, we demonstrate in silico that our proposed strategy allows us to dampen the emerging focal oscillatory sources using only a small set of electrodes.

SIGNIFICANCE

Our integrative study informs the development of modulation and control algorithms for neurostimulation that could improve the effectiveness of implantable, closed-loop anti-epileptic devices.

Interictal Epileptiform Discharges and the Quality of Human Intracranial Neurophysiology Data

Intracranial electroencephalography (IEEG) involves recording from electrodes placed directly onto the cortical surface or deep brain locations. It is performed on patients with medically refractory epilepsy, undergoing pre-surgical seizure localization. IEEG recordings, combined with advancements in computational capacity and analysis tools, have accelerated cognitive neuroscience. This Perspective describes a potential pitfall latent in many of these recordings by virtue of the subject population—namely interictal epileptiform discharges (IEDs), which can cause spurious results due to the contamination of normal neurophysiological signals by pathological waveforms related to epilepsy. We first discuss the nature of IED hazards, and why they deserve the attention of neurophysiology researchers. We then describe four general strategies used when handling IEDs (manual identification, automated identification, manual-automated hybrids, and ignoring by leaving them in the data), and discuss their pros, cons, and contextual factors. Finally, we describe current practices of human neurophysiology researchers worldwide based on a cross-sectional literature review and a voluntary survey. We put these results in the context of the listed strategies and make suggestions on improving awareness and clarity of reporting to enrich both data quality and communication in the field.

Epileptic source connectivity analysis based on estimating of dynamic time series of regions of interest

We propose a new source connectivity method by focusing on estimating time courses of the regions of interest (ROIs). To this aim, it is necessary to consider the strong inherent non-stationary behavior of neural activity. We develop an iterative dynamic approach to extract a single time course for each ROI encoding the temporal non-stationary features. The proposed approach explicitly includes dynamic constraints by taking into account the evolution of the sources activities for further dynamic connectivity analysis. We simulated an epileptic network with a non-stationary structure; accordingly, EEG source reconstruction using LORETA is performed. Using the reconstructed sources, the spatially compact ROIs are selected. Then, a single time course encoding the temporal non-stationarity is extracted for each ROI. An adaptive directed transfer function (ADTF) is applied to measure the information flow of underlying brain networks. Obtained results demonstrate that the contributed approach is more efficient to estimate the ROI time series and ROI to ROI information flow in comparison with existing methods. Our work is validated in three drug-resistance epilepsy patients. The proposed ROI time series estimation directly affects the quality of connectivity analysis, leading to the best possible seizure onset zone (SOZ) localization verified by electrocorticography and post-operational results.

Effects of disparity on visual discomfort caused by short-term stereoscopic viewing based on electroencephalograph analysis

BACKGROUND

Discomfort evoked by stereoscopic depth has been widely concerned. Previous studies have proposed a comfortable disparity range and considered that disparities exceed this range would cause visual discomfort. Brain activity recordings including Electroencephalograph (EEG) monitoring enable better understanding of perceptual and cognitive processes related to stereo depth-induced visual comfort.

METHODS

EEG data was collected using a stereo-visual evoked potential (VEP) test system by providing visual stimulus to subjects aged from 21 to 25 with normal stereoscopic vision. For each type of visual stimulus, data were processed using directed transfer function (DTF) and adaptive directed transfer function (ADTF) in combination with subjective feedbacks (comfort or discomfort). The topographies of information flow were constructed to compare responses stimulated by different stereoscopic depth, and to determine the difference in comfort and discomfort situations upon stimulation with same stereoscopic depth.

RESULTS

Based on EEG analysis results, we found that the occipital P270 was moderately related to the disparity. Moreover, the ADTF of P270 showed that the information flows at frontal lobe and central-parietal lobe changed when stimulation with different stereoscopic depth applied. As to the stereo images with same stereoscopic depth, the DTF outflows at the temporal and temporal-parietal lobes in δ band, central and central-parietal lobes in α and θ bands, and the comparison of inflows in these three bands could be considered as discriminated indexes for matching the stereoscopic effect with viewers' comfort or discomfort state impacted by disparity. The subjective feedbacks indicated that the comfort judgments remained as a result of cumulative effect.

CONCLUSIONS

This study proposed a short-term stereo-VEP experiment that shorted the duration of each stimulus in the experimental scheme to minimize the interference from other factors except the disparity. The occipital P270 had a mid-relevance to the disparity and its ADTF showed the affected areas when viewers are receiving stimulations with different disparities. DTF could be considered as discriminated indexes for matching the stereoscopic effect with viewers' comfort or discomfort state induced by disparity. This study proposed a preferable experiment to observe the single effect of disparity and provided an intuitive and easy-to-read result in a more convenient manner.

Localization of Epileptogenic Zone With the Correction of Pathological Networks

Patients with focal drug-resistant epilepsy are potential candidates for surgery. Stereo-electroencephalograph (SEEG) is often considered as the “gold standard” to identify the epileptogenic zone (EZ) that accounts for the onset and propagation of epileptiform discharges. However, visual analysis of SEEG still prevails in clinical practice. In addition, epilepsy is increasingly understood to be the result of network disorder, but the specific organization of the epileptic network is still unclear. Therefore, it is necessary to quantitatively localize the EZ and investigate the nature of epileptogenic networks. In this study, intracranial recordings from 10 patients were analyzed through adaptive directed transfer function, and the out-degree of effective network was selected as the principal indicator to localize the epileptogenic area. Furthermore, a coupled neuronal population model was used to qualitatively simulate electrical activity in the brain. By removing individual populations, virtual surgery adjusting the network organization could be performed. Results suggested that the accuracy and detection rate of the EZ localization were 82.86 and 85.29%, respectively. In addition, the same stage shared a relatively stable connectivity pattern, while the patterns changed with transition to different processes. Meanwhile, eight cases of simulations indicated that networks in the ictal stage were more likely to generate rhythmic spikes. This indicated the existence of epileptogenic networks, which could enhance local excitability and facilitate synchronization. The removal of the EZ could correct these pathological networks and reduce the amount of spikes by at least 75%. This might be one reason why accurate resection could reduce or even suppress seizures. This study provides novel insights into epilepsy and surgical treatments from the network perspective.

Early alterations of social brain networks in young children with autism

Social impairments are a hallmark of Autism Spectrum Disorders (ASD), but empirical evidence for early brain network alterations in response to social stimuli is scant in ASD. We recorded the gaze patterns and brain activity of toddlers with ASD and their typically developing peers while they explored dynamic social scenes. Directed functional connectivity analyses based on electrical source imaging revealed frequency specific network atypicalities in the theta and alpha frequency bands, manifesting as alterations in both the driving and the connections from key nodes of the social brain associated with autism. Analyses of brain-behavioural relationships within the ASD group suggested that compensatory mechanisms from dorsomedial frontal, inferior temporal and insular cortical regions were associated with less atypical gaze patterns and lower clinical impairment. Our results provide strong evidence that directed functional connectivity alterations of social brain networks is a core component of atypical brain development at early stages of ASD.

Newborns are attracted to voices, faces and social gestures. Paying attention to these social cues in everyday life helps infants and young children learn how to interact with others. During this period of development, a network of connections forms between different parts of the brain that helps children to understand other people’s social behaviors.

During their first year of life, infants who later develop autism spectrum disorders (ASD) pay less attention to social cues. This early indifference to these important signals leads to social deficits in children with ASD. They are less able to understand other people’s behaviors or engage in typical social interactions. It’s not yet clear why children with ASD are less attuned to social cues. But is likely that the development of brain networks essential for understanding social behavior suffers as a result. Studying how such networks develop in typical very young children and those with ASD may help scientist learn more.

Now, Sperdin et al. confirm there are differences in the social brain-networks of very young children with ASD compared with their typical peers. In the experiment, 3-year-old children with ASD and without watched videos of other children playing, while Sperdin et al. recorded what they looked at and what happened in their brains. Eyemovements were measured with a tracker, and the brain activity was recorded using an electroencephalogram (EEG), which uses sensors placed on the scalp to measure electrical signals.

What children with ASD looked at was different than their typical peers, and these differences corresponded with alterations in the brain networks that process social information. Children with ASD who had less severe symptoms had stronger activity in these brain networks. What they looked at also was more similar to typical children. This suggests less severely affected children with ASD may be able to compensate that way.

Identifying ASD-like behaviors and brain differences early in life may help scientists to better understand what causes the condition. It may also help clinicians provide more individualized therapies early in life when the brain is most adaptable. Long-term studies of these brain-network differences in children with ASD are necessary to better understand how therapies can influence these changes.

Increased neuronal synchrony prepares mesial temporal networks for seizures of neocortical origin

OBJECTIVE

To gain understanding of the neuronal mechanisms underlying regional seizure spread, the impact of regional synchrony between seizure focus and downstream networks on neuronal activity during the transition to seizure in those downstream networks was assessed.

METHODS

Seven patients undergoing diagnostic intracranial electroencephalographic studies for surgical resection of epileptogenic regions were implanted with subdural clinical electrodes into the cortex (site of seizure initiation) and mesial temporal lobe (MTL) structures (downstream) as well as microwires into MTL. Neural activity was recorded (24/7) in parallel with the clinical intracranial electroencephalogram recordings for the duration of the patient's diagnostic stay. Changes in (1) regional synchrony (ie, coherence) between the presumptive neocortical seizure focus and MTL, (2) local synchrony between MTL neurons and their local field potential, and (3) neuronal firing rates within MTL in the time leading up to seizure were examined to study the mechanisms underlying seizure spread.

RESULTS

In seizures of neocortical origin, an increase in regional synchrony preceded the spread of seizures into MTL (predominantly hippocampal). Within frequencies similar to those of regional synchrony, MTL networks showed an increase in unit-field coherence and a decrease in neuronal firing rate, specifically for inhibitory interneuron populations but not pyramidal cell populations.

SIGNIFICANCE

These results suggest a mechanism of spreading seizures whereby the seizure focus first synchronizes local field potentials in downstream networks to the seizure activity. This change in local field coherence modifies the activity of interneuron populations in these downstream networks, which leads to the attenuation of interneuronal firing rate, effectively shutting down local interneuron populations prior to the spread of seizure. Therefore, regional synchrony may influence the failure of downstream interneurons to prevent the spread of the seizures during generalization.

Electromagnetic source imaging using simultaneous scalp EEG and intracranial EEG: An emerging tool for interacting with pathological brain networks

OBJECTIVE

The goal of this study is to investigate the performance, merits and limitations of source imaging using intracranial EEG (iEEG) recordings and to compare its accuracy to the results of EEG source imaging. Accuracy in this study, is measured both by determining the location and inter-nodal connectivity of underlying brain networks.

METHODS

Systematic computer simulation studies are conducted to evaluate iEEG-based source imaging vs. EEG-based source imaging, and source imaging using both EEG and iEEG. To test the source imaging models, networks of inter-connected nodes (in terms of activity) are simulated. The location of the network nodes is randomly selected within a realistic geometry head model and a connectivity link is created among these nodes based on a multi-variate auto-regressive (MVAR) model. Then the forward problem is solved to calculate the potentials at the electrodes and noise (white and correlated) is added to these simulated potentials to simulate realistic measurements. Subsequently, the inverse problem is solved and an algorithm based on principle component analysis is performed on the estimated source activities to determine the location of the simulated network nodes. The activity of these nodes (over time), is then extracted, and used to estimate the connectivity links among the mentioned nodes using Granger causality analysis.

RESULTS

Source imaging based on iEEG recordings may or may not improve the accuracy in localization, depending on the number and location of active nodes relative to iEEG electrodes and to other nodes within the network. However, our simulation results suggest that combining EEG and iEEG modalities (simultaneous scalp and intracranial recordings) can improve the imaging accuracy significantly.

CONCLUSIONS

While iEEG source imaging is useful in estimating the exact location of sources near the iEEG electrodes, combining EEG and iEEG recordings can achieve a more accurate imaging due to the high spatial coverage of the scalp electrodes and the added near field information provided by the iEEG electrodes.

SIGNIFICANCE

The present results suggest the feasibility of localizing brain electrical sources from iEEG recordings and improving EEG source localization using simultaneous EEG and iEEG recordings to cover the whole brain. The hybrid EEG and iEEG source imaging can assist the clinicians when unequivocal decisions about determining the epileptogenic zone cannot be reached using a single modality.

Automated diagnosis of temporal lobe epilepsy in the absence of interictal spikes

Objective

To diagnose and lateralise temporal lobe epilepsy (TLE) by building a classification system that uses directed functional connectivity patterns estimated during EEG periods without visible pathological activity.

Methods

Resting-state high-density EEG recording data from 20 left TLE patients, 20 right TLE patients and 35 healthy controls was used. Epochs without interictal spikes were selected. The cortical source activity was obtained for 82 regions of interest and whole-brain directed functional connectivity was estimated in the theta, alpha and beta frequency bands. These connectivity values were then used to build a classification system based on two two-class Random Forests classifiers: TLE vs healthy controls and left vs right TLE. Feature selection and classifier training were done in a leave-one-out procedure to compute the mean classification accuracy.

Results

The diagnosis and lateralization classifiers achieved a high accuracy (90.7% and 90.0% respectively), sensitivity (95.0% and 90.0% respectively) and specificity (85.7% and 90.0% respectively). The most important features for diagnosis were the outflows from left and right medial temporal lobe, and for lateralization the right anterior cingulate cortex. The interaction between features was important to achieve correct classification.

Significance

This is the first study to automatically diagnose and lateralise TLE based on EEG. The high accuracy achieved demonstrates the potential of directed functional connectivity estimated from EEG periods without visible pathological activity for helping in the diagnosis and lateralization of TLE.

•

Both classifiers for TLE diagnosis and lateralization achieved high accuracy (90%).

•

Outflow from left and right hippocampus was the most important for diagnosis.

•

Outflow from the right ACC was the most important for lateralization.

•

The interaction between features is important for a correct classification.

Time-Varying Networks of Inter-Ictal Discharging Reveal Epileptogenic Zone

The neuronal synchronous discharging may cause an epileptic seizure. Currently, most of the studies conducted to investigate the mechanism of epilepsy are based on EEGs or functional magnetic resonance imaging (fMRI) recorded during the ictal discharging or the resting-state, and few studies have probed into the dynamic patterns during the inter-ictal discharging that are much easier to record in clinical applications. Here, we propose a time-varying network analysis based on adaptive directed transfer function to uncover the dynamic brain network patterns during the inter-ictal discharging. In addition, an algorithm based on the time-varying outflow of information derived from the network analysis is developed to detect the epileptogenic zone. The analysis performed revealed the time-varying network patterns during different stages of inter-ictal discharging; the epileptogenic zone was activated prior to the discharge onset then worked as the source to propagate the activity to other brain regions. Consistence between the epileptogenic zones detected by our proposed approach and the actual epileptogenic zones proved that time-varying network analysis could not only reveal the underlying neural mechanism of epilepsy, but also function as a useful tool in detecting the epileptogenic zone based on the EEGs in the inter-ictal discharging.

Toward a Mechanistic Understanding of Epileptic Networks

Focal epileptic seizures have long been considered to arise from a small susceptible brain area and spread through uninvolved regions. In the past decade, the idea that focal seizures instead arise from coordinated activity across large-scale epileptic networks has become widely accepted. Understanding the network model's applicability is critical, due to its increasing influence on clinical research and surgical treatment paradigms. In this review, we examine the origins of the concept of epileptic networks as the nidus for recurring seizures. We summarize analytical and methodological elements of epileptic network studies and discuss findings from recent detailed electrophysiological investigations. Our review highlights the strengths and limitations of the epileptic network theory as a metaphor for the complex interactions that occur during seizures. We present lines of investigation that may usefully probe these interactions and thus serve to advance our understanding of the long-range effects of epileptiform activity.

Recurring Functional Interactions Predict Network Architecture of Interictal and Ictal States in Neocortical Epilepsy

Human epilepsy patients suffer from spontaneous seizures, which originate in brain regions that also subserve normal function. Prior studies demonstrate focal, neocortical epilepsy is associated with dysfunction, several hours before seizures. How does the epileptic network perpetuate dysfunction during baseline periods? To address this question, we developed an unsupervised machine learning technique to disentangle patterns of functional interactions between brain regions, or subgraphs, from dynamic functional networks constructed from approximately 100 h of intracranial recordings in each of 22 neocortical epilepsy patients. Using this approach, we found: (1) subgraphs from ictal (seizure) and interictal (baseline) epochs are topologically similar, (2) interictal subgraph topology and dynamics can predict brain regions that generate seizures, and (3) subgraphs undergo slower and more coordinated fluctuations during ictal epochs compared to interictal epochs. Our observations suggest that seizures mark a critical shift away from interictal states that is driven by changes in the dynamical expression of strongly interacting components of the epileptic network.

Spatiotemporal Mapping of Interictal Spike Propagation: A Novel Methodology Applied to Pediatric Intracranial EEG Recordings

Synchronized cortical activity is implicated in both normative cognitive functioning and many neurologic disorders. For epilepsy patients with intractable seizures, irregular synchronization within the epileptogenic zone (EZ) is believed to provide the network substrate through which seizures initiate and propagate. Mapping the EZ prior to epilepsy surgery is critical for detecting seizure networks in order to achieve postsurgical seizure control. However, automated techniques for characterizing epileptic networks have yet to gain traction in the clinical setting. Recent advances in signal processing and spike detection have made it possible to examine the spatiotemporal propagation of interictal spike discharges across the epileptic cortex. In this study, we present a novel methodology for detecting, extracting, and visualizing spike propagation and demonstrate its potential utility as a biomarker for the EZ. Eighteen presurgical intracranial EEG recordings were obtained from pediatric patients ultimately experiencing favorable (i.e., seizure-free, n = 9) or unfavorable (i.e., seizure-persistent, n = 9) surgical outcomes. Novel algorithms were applied to extract multichannel spike discharges and visualize their spatiotemporal propagation. Quantitative analysis of spike propagation was performed using trajectory clustering and spatial autocorrelation techniques. Comparison of interictal propagation patterns revealed an increase in trajectory organization (i.e., spatial autocorrelation) among Sz-Free patients compared with Sz-Persist patients. The pathophysiological basis and clinical implications of these findings are considered.

Noninvasive Electromagnetic Source Imaging and Granger Causality Analysis: An Electrophysiological Connectome (eConnectome) Approach

OBJECTIVE

Combined source-imaging techniques and directional connectivity analysis can provide useful information about the underlying brain networks in a noninvasive fashion. Source-imaging techniques have been used successfully to either determine the source of activity or to extract source time-courses for Granger causality analysis, previously. In this work, we utilize source-imaging algorithms to both find the network nodes [regions of interest (ROI)] and then extract the activation time series for further Granger causality analysis. The aim of this work is to find network nodes objectively from noninvasive electromagnetic signals, extract activation time-courses, and apply Granger analysis on the extracted series to study brain networks under realistic conditions.

METHODS

Source-imaging methods are used to identify network nodes and extract time-courses and then Granger causality analysis is applied to delineate the directional functional connectivity of underlying brain networks. Computer simulations studies where the underlying network (nodes and connectivity pattern) is known were performed; additionally, this approach has been evaluated in partial epilepsy patients to study epilepsy networks from interictal and ictal signals recorded by EEG and/or Magnetoencephalography (MEG).

RESULTS

Localization errors of network nodes are less than 5 mm and normalized connectivity errors of ∼20% in estimating underlying brain networks in simulation studies. Additionally, two focal epilepsy patients were studied and the identified nodes driving the epileptic network were concordant with clinical findings from intracranial recordings or surgical resection.

CONCLUSION

Our study indicates that combined source-imaging algorithms with Granger causality analysis can identify underlying networks precisely (both in terms of network nodes location and internodal connectivity).

SIGNIFICANCE

The combined source imaging and Granger analysis technique is an effective tool for studying normal or pathological brain conditions.

Dynamic Network Connectivity Analysis to Identify Epileptogenic Zones Based on Stereo-Electroencephalography

Objectives: Accurate localization of epileptogenic zones (EZs) is essential for successful surgical treatment of refractory focal epilepsy. The aim of the present study is to investigate whether a dynamic network connectivity analysis based on stereo-electroencephalography (SEEG) signals is effective in localizing EZs.

Methods: SEEG data were recorded from seven patients who underwent presurgical evaluation for the treatment of refractory focal epilepsy and for whom the subsequent resective surgery gave a good outcome. A time-variant multivariate autoregressive model was constructed using a Kalman filter, and the time-variant partial directed coherence was computed. This was then used to construct a dynamic directed network model of the epileptic brain. Three graph measures (in-degree, out-degree, and betweenness centrality) were used to analyze the characteristics of the dynamic network and to find the important nodes in it.

Results: In all seven patients, the indicative EZs localized by the in-degree and the betweenness centrality were highly consistent with the clinically diagnosed EZs. However, the out-degree did not indicate any significant differences between nodes in the network.

Conclusions: In this work, a method based on ictal SEEG signals and effective connectivity analysis localized EZs accurately. The results suggest that the in-degree and betweenness centrality may be better network characteristics to localize EZs than the out-degree.

Directed Functional Brain Connectivity Based on EEG Source Imaging: Methodology and Application to Temporal Lobe Epilepsy

OBJECTIVE

The importance of functional brain connectivity to study physiological and pathological brain activity has been widely recognized. Here, we aimed to 1) review a methodological pipeline to investigate directed functional connectivity between brain regions using source signals derived from high-density EEG; 2) elaborate on some methodological challenges; and 3) apply this pipeline to temporal lobe epilepsy (TLE) patients and healthy controls to investigate directed functional connectivity differences in the theta and beta frequency bands during EEG epochs without visible pathological activity.

METHODS

The methodological pipeline includes: EEG acquisition and preprocessing, electrical-source imaging (ESI) using individual head models and distributed inverse solutions, parcellation of the gray matter in regions of interest, fixation of the dipole orientation for each region, computation of the spectral power in the source space, and directed functional connectivity estimation using Granger-causal modeling. We specifically analyzed how the signal-to-noise ratio (SNR) changes using different approaches for the dipole orientation fixation. We applied this pipeline to 20 left TLE patients, 20 right TLE patients, and 20 healthy controls.

RESULTS

Projecting each dipole to the predominant dipole orientation leads to a threefold SNR increase as compared to the norm of the dipoles. By comparing connectivity in TLE versus controls, we found significant frequency-specific outflow differences in physiologically plausible regions.

CONCLUSION

The results suggest that directed functional connectivity derived from ESI can help better understand frequency-specific resting-state network alterations underlying focal epilepsy.

SIGNIFICANCE

EEG-based directed functional connectivity could contribute to the search of new biomarkers of this disorder.

Seizure-Onset Mapping Based on Time-Variant Multivariate Functional Connectivity Analysis of High-Dimensional Intracranial EEG: A Kalman Filter Approach

The visual interpretation of intracranial EEG (iEEG) is the standard method used in complex epilepsy surgery cases to map the regions of seizure onset targeted for resection. Still, visual iEEG analysis is labor-intensive and biased due to interpreter dependency. Multivariate parametric functional connectivity measures using adaptive autoregressive (AR) modeling of the iEEG signals based on the Kalman filter algorithm have been used successfully to localize the electrographic seizure onsets. Due to their high computational cost, these methods have been applied to a limited number of iEEG time-series (<60). The aim of this study was to test two Kalman filter implementations, a well-known multivariate adaptive AR model (Arnold et al. 1998) and a simplified, computationally efficient derivation of it, for their potential application to connectivity analysis of high-dimensional (up to 192 channels) iEEG data. When used on simulated seizures together with a multivariate connectivity estimator, the partial directed coherence, the two AR models were compared for their ability to reconstitute the designed seizure signal connections from noisy data. Next, focal seizures from iEEG recordings (73-113 channels) in three patients rendered seizure-free after surgery were mapped with the outdegree, a graph-theory index of outward directed connectivity. Simulation results indicated high levels of mapping accuracy for the two models in the presence of low-to-moderate noise cross-correlation. Accordingly, both AR models correctly mapped the real seizure onset to the resection volume. This study supports the possibility of conducting fully data-driven multivariate connectivity estimations on high-dimensional iEEG datasets using the Kalman filter approach.

Assessing directionality and strength of coupling through symbolic analysis: an application to epilepsy patients

Inferring strength and direction of interactions from electroencephalographic (EEG) recordings is of crucial importance to improve our understanding of dynamical interdependencies underlying various physiological and pathophysiological conditions in the human epileptic brain. We here use approaches from symbolic analysis to investigate--in a time-resolved manner--weighted and directed, short- to long-ranged interactions between various brain regions constituting the epileptic network. Our observations point to complex spatial-temporal interdependencies underlying the epileptic process and their role in the generation of epileptic seizures, despite the massive reduction of the complex information content of multi-day, multi-channel EEG recordings through symbolization. We discuss limitations and potential future improvements of this approach.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy

OBJECTIVE

There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity.

METHODS

In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments.

RESULTS

Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits.

SIGNIFICANCE

The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

Localization of epileptogenic zones in Lennox-Gastaut syndrome (LGS) using graph theoretical analysis of ictal intracranial EEG: a preliminary investigation

INTRODUCTION

Precise localization of epileptogenic zones is essential for the successful surgical treatment of refractory epilepsy including Lennox-Gastaut syndrome (LGS). The surgical resection areas are generally determined by epileptologists based on diverse neuroimaging modalities; however, exact epileptogenic zones cannot be accurately localized in many patients with LGS using the conventional methods. Therefore, new reliable algorithms are still required for enhancing the success rate of the resective epilepsy surgery. In the present study, we introduce an approach to localize epileptogenic zones in LGS based on the graph theoretical analysis of ical intracranial EEG (iEEG).

METHODS

Four patients with LGS who became seizure-free after the resective epilepsy surgery were selected. Before the surgery, their epileptogenic zones were delineated using EEG, iEEG, and several conventional imaging modalities. Phase locking value (PLV) analysis was applied to construct functional connectivity networks during ictal events, and then several graph theoretical indices including betweenness centrality (BC) were evaluated for each iEEG sensor to find the primary hubs of the ictal epileptic network. The graph theoretical index values were then overlaid on 3D individual cortical surface.

RESULTS

The iEEG channels with high BC values coincided well with the surgical resection areas. Among various graph theoretical measures such as local efficiency, participation coefficient, and eigenvector centrality, only BC showed fair correspondence with the surgical resection areas.

CONCLUSIONS

The primary hubs in the ictal epileptic networks coincided well with areas of surgical resection in LGS patients with successful surgical outcomes. This observation warrants further studies to determine if the graph theoretical network analysis of ictal iEEG recordings can serve as a new auxiliary tool to localize epileptogenic zones in LGS.

Functional brain connectivity from EEG in epilepsy: seizure prediction and epileptogenic focus localization

Today, neuroimaging techniques are frequently used to investigate the integration of functionally specialized brain regions in a network. Functional connectivity, which quantifies the statistical dependencies among the dynamics of simultaneously recorded signals, allows to infer the dynamical interactions of segregated brain regions. In this review we discuss how the functional connectivity patterns obtained from intracranial and scalp electroencephalographic (EEG) recordings reveal information about the dynamics of the epileptic brain and can be used to predict upcoming seizures and to localize the seizure onset zone. The added value of extracting information that is not visibly identifiable in the EEG data using functional connectivity analysis is stressed. Despite the fact that many studies have showed promising results, we must conclude that functional connectivity analysis has not made its way into clinical practice yet.

Ictal propagation of high frequency activity is recapitulated in interictal recordings: effective connectivity of epileptogenic networks recorded with intracranial EEG

Seizures are increasingly understood to arise from epileptogenic networks across which ictal activity is propagated and sustained. In patients undergoing invasive monitoring for epilepsy surgery, high frequency oscillations have been observed within the seizure onset zone during both ictal and interictal intervals. We hypothesized that the patterns by which high frequency activity is propagated would help elucidate epileptogenic networks and thereby identify network nodes relevant for surgical planning. Intracranial EEG recordings were analyzed with a multivariate autoregressive modeling technique (short-time direct directed transfer function--SdDTF), based on the concept of Granger causality, to estimate the directionality and intensity of propagation of high frequency activity (70-175 Hz) during ictal and interictal recordings. These analyses revealed prominent divergence and convergence of high frequency activity propagation at sites identified by epileptologists as part of the ictal onset zone. In contrast, relatively little propagation of this activity was observed among the other analyzed sites. This pattern was observed in both subdural and depth electrode recordings of patients with focal ictal onset, but not in patients with a widely distributed ictal onset. In patients with focal ictal onsets, the patterns of propagation recorded during pre-ictal (up to 5 min immediately preceding ictal onset) and interictal (more than 24h before and after seizures) intervals were very similar to those recorded during seizures. The ability to characterize epileptogenic networks from interictal recordings could have important clinical implications for epilepsy surgery planning by reducing the need for prolonged invasive monitoring to record spontaneous seizures.

Assessing dynamic spectral causality by lagged adaptive directed transfer function and instantaneous effect factor

It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The nonzero covariance of the model's residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the "causal ordering" is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In this study, we first investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in time-variant spectral causality assessment was demonstrated through the mutual causality investigation of brain activity during the VEP experiments.

Detection of interictal epileptiform discharges using signal envelope distribution modelling: application to epileptic and non-epileptic intracranial recordings

Interictal epileptiform discharges (spikes, IEDs) are electrographic markers of epileptic tissue and their quantification is utilized in planning of surgical resection. Visual analysis of long-term multi-channel intracranial recordings is extremely laborious and prone to bias. Development of new and reliable techniques of automatic spike detection represents a crucial step towards increasing the information yield of intracranial recordings and to improve surgical outcome. In this study, we designed a novel and robust detection algorithm that adaptively models statistical distributions of signal envelopes and enables discrimination of signals containing IEDs from signals with background activity. This detector demonstrates performance superior both to human readers and to an established detector. It is even capable of identifying low-amplitude IEDs which are often missed by experts and which may represent an important source of clinical information. Application of the detector to non-epileptic intracranial data from patients with intractable facial pain revealed the existence of sharp transients with waveforms reminiscent of interictal discharges that can represent biological sources of false positive detections. Identification of these transients enabled us to develop and propose secondary processing steps, which may exclude these transients, improving the detector's specificity and having important implications for future development of spike detectors in general.

Short-time windowed covariance: a metric for identifying non-stationary, event-related covariant cortical sites

BACKGROUND

Electrocorticography (ECoG) signals can provide high spatio-temporal resolution and high signal to noise ratio recordings of local neural activity from the surface of the brain. Previous studies have shown that broad-band, spatially focal, high-frequency increases in ECoG signals are highly correlated with movement and other cognitive tasks and can be volitionally modulated. However, significant additional information may be present in inter-electrode interactions, but adding additional higher order inter-electrode interactions can be impractical from a computational aspect, if not impossible.

NEW METHOD

In this paper we present a new method of calculating high frequency interactions between electrodes called Short-Time Windowed Covariance (STWC) that builds on mathematical techniques currently used in neural signal analysis, along with an implementation that accelerates the algorithm by orders of magnitude by leveraging commodity, off-the-shelf graphics processing unit (GPU) hardware.

RESULTS

Using the hardware-accelerated implementation of STWC, we identify many types of event-related inter-electrode interactions from human ECoG recordings on global and local scales that have not been identified by previous methods. Unique temporal patterns are observed for digit flexion in both low- (10mm spacing) and high-resolution (3mm spacing) electrode arrays.

COMPARISON WITH EXISTING METHODS

Covariance is a commonly used metric for identifying correlated signals, but the standard covariance calculations do not allow for temporally varying covariance. In contrast STWC allows and identifies event-driven changes in covariance without identifying spurious noise correlations.

CONCLUSIONS

STWC can be used to identify event-related neural interactions whose high computational load is well suited to GPU capabilities.

Canonical Granger causality between regions of interest

Estimating and modeling functional connectivity in the brain is a challenging problem with potential applications in the understanding of brain organization and various neurological and neuropsychological conditions. An important objective in connectivity analysis is to determine the connections between regions of interest in the brain. However, traditional functional connectivity analyses have frequently focused on modeling interactions between time series recordings at individual sensors, voxels, or vertices despite the fact that a single region of interest will often include multiple such recordings. In this paper, we present a novel measure of interaction between regions of interest rather than individual signals. The proposed measure, termed canonical Granger causality, combines ideas from canonical correlation and Granger causality analysis to yield a measure that reflects directed causality between two regions of interest. In particular, canonical Granger causality uses optimized linear combinations of signals from each region of interest to enable accurate causality measurements from substantially less data compared to alternative multivariate methods that have previously been proposed for this scenario. The optimized linear combinations are obtained using a variation of a technique developed for optimization on the Stiefel manifold. We demonstrate the advantages of canonical Granger causality in comparison to alternative causality measures for a range of different simulated datasets. We also apply the proposed measure to local field potential data recorded in a macaque brain during a visuomotor task. Results demonstrate that canonical Granger causality can be used to identify causal relationships between striate and prestriate cortexes in cases where standard Granger causality is unable to identify statistically significant interactions.

EEG source imaging in epilepsy--practicalities and pitfalls

EEG source imaging (ESI) is a model-based imaging technique that integrates temporal and spatial components of EEG to identify the generating source of electrical potentials recorded on the scalp. Recent advances in computer technologies have made the analysis of ESI data less time-consuming, and have rekindled interest in this technique as a clinical diagnostic tool. On the basis of the available body of evidence, ESI seems to be a promising tool for epilepsy evaluation; however, the precise clinical value of ESI in presurgical evaluation of epilepsy and in localization of eloquent cortex remains to be investigated. In this Review, we describe two fundamental issues in ESI; namely, the forward and inverse problems, and their solutions. The clinical application of ESI in surgical planning for patients with medically refractory focal epilepsy, and its use in source reconstruction together with invasive recordings, is also discussed. As ESI can be used to map evoked responses, we discuss the clinical utility of this technique in cortical mapping-an essential process when planning resective surgery for brain regions that are in close proximity to eloquent cortex.

Dynamic imaging of seizure activity in pediatric epilepsy patients

OBJECTIVE

To investigate the feasibility of using noninvasive EEG source imaging approach to image continuous seizure activity in pediatric epilepsy patients.

METHODS

Nine pediatric patients with medically intractable epilepsy were included in this study. Eight of the patients had extratemporal lobe epilepsy and one had temporal lobe epilepsy. All of the patients underwent resective surgery and seven of them underwent intracranial EEG (iEEG) monitoring. The ictal EEG was analyzed using a noninvasive dynamic seizure imaging (DSI) approach. The DSI approach separates scalp EEGs into independent components and extracts the spatio-temporal ictal features to achieve dynamic imaging of seizure sources. Surgical resection and intracranial recordings were used to validate the noninvasive imaging results.

RESULTS

The DSI determined seizure onset zones (SOZs) in these patients were localized within or in close vicinity to the surgically resected region. In the seven patients with intracranial monitoring, the estimated seizure onset sources were concordant with the seizure onset zones of iEEG. The DSI also localized the multiple foci involved in the later seizure propagation, which were confirmed by the iEEG recordings.

CONCLUSIONS

Dynamic seizure imaging can noninvasively image the seizure activations in pediatric patients with both temporal and extratemporal lobe epilepsy.

SIGNIFICANCE

EEG seizure imaging can potentially be used to noninvasively image the SOZs and aid the pre-surgical planning in pediatric epilepsy patients.

Localization of ictal onset zones in Lennox-Gastaut syndrome using directed transfer function method

Neuroscientists are becoming interested in the application of computational EEG analysis to the identification of ictal onset zones, however, most studies have focused on the localization of ictal onset zones in focal epilepsy. The present study aimed to estimate the ictal onset zone of Lennox-Gastaut syndrome (LGS) with bilaterally synchronous epileptiform discharges from intracranial electroencephalography (iEEG) recordings using directional connectivity analysis. We analyzed ictal iEEG data acquired from three patients with LGS. To identify the ictal onset zones, we estimated the functional directional connectivity network among the intracerebral electrodes using the directed transfer function (DTF) method. The analysis results demonstrated that areas with high average outflow values corresponded well with the surgical resection areas identified using electrophysiologic data and conventional neuroimaging modalities. Our results suggest that the DTF analysis can be a useful auxiliary tool for determining surgical resection areas prior to epilepsy surgery in LGS patients.

Partial directed coherence analysis of intracranial neural spikes in epilepsy patients

Intracranial electroencephalograms (EEG) provide a direct observation of neural activity by placing an electrode array on the cortical surface near the suspected epileptic foci. The neural spikes appeared during inter-ictal stages are mainly produced by abnormal neural discharges from epileptic foci. The topological mapping of spikes' potentials is commonly used to identify the epileptogenic zone. However, the propagations among multi-channel spikes are also important to identify the epileptic source activity. In addition, the changes of source activities in a series of consecutive spikes reveal the time-varying neural activations during discharge process, which provide alternative information for interpreting epileptic source activity. This paper proposes a spike classification based on the similarity of phase-space features to select candidate spikes from the intracranial EEGs recorded from an 8×8 electrocorticogram grid. Then, the partial directed coherence (PDC), which can provide the flow of source activity, at each spiking time point is computed. The outflow PDCs of all electrodes are therefore displayed on the grid. Our result showed that the derived source activities in the preceding spikes had high concentrated distributions but decreased in latter spikes. This implied the epileptic discharges were initially induced by a small-area cortex neurons and then spread out.

Seizure source imaging by means of FINE spatio-temporal dipole localization and directed transfer function in partial epilepsy patients

OBJECTIVE

To investigate the usage of a high-density EEG recording system and source imaging technique for localizing seizure activity in patients with medically intractable partial epilepsy.

METHODS

High-density, 76-channel scalp EEG signals were recorded in 10 patients with partial epilepsy. The patients underwent routine clinical pre-surgical evaluation and all had resective surgery with seizure free outcome. After applying a FINE (first principle vectors) spatio-temporal source localization and DTF (directed transfer function) connectivity analysis approach, ictal sources were imaged. Effects of number of scalp EEG electrodes on the seizure localization were also assessed using 76, 64, 48, 32, and 21 electrodes, respectively.

RESULTS

Surgical resections were used to assess the source imaging results. Results from the 76-channel EEG in the 10 patients showed high correlation with the surgically resected brain regions. The localization of seizure onset zone from 76-channel EEG showed improved source detection accuracy compared to other EEG configurations with fewer electrodes.

CONCLUSIONS

FINE together with DTF was able to localize seizure onset zones of partial epilepsy patients. High-density EEG recording can help achieve improved seizure source imaging.

SIGNIFICANCE

The present results suggest the promise of high-density EEG and electrical source imaging for noninvasively localizing seizure onset zones.