Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies

The interictal suppression hypothesis is the dominant differentiator of seizure onset zones in focal epilepsy

Successful surgical treatment of drug-resistant epilepsy traditionally relies on the identification of seizure onset zones (SOZs). Connectome-based analyses of electrographic data from stereo electroencephalography (SEEG) may empower improved detection of SOZs. Specifically, connectome-based analyses based on the interictal suppression hypothesis posit that when the patient is not having a seizure, SOZs are inhibited by non-SOZs through high inward connectivity and low outward connectivity. However, it is not clear whether there are other motifs that can better identify potential SOZs. Thus, we sought to use unsupervised machine learning to identify network motifs that elucidate SOZs and investigate if there is another motif that outperforms the ISH. Resting-state SEEG data from 81 patients with drug-resistant epilepsy undergoing a pre-surgical evaluation at Vanderbilt University Medical Center were collected. Directed connectivity matrices were computed using the alpha band (8-13 Hz). Principal component analysis (PCA) was performed on each patient's connectivity matrix. Each patient's components were analysed qualitatively to identify common patterns across patients. A quantitative definition was then used to identify the component that most closely matched the observed pattern in each patient. A motif characteristic of the interictal suppression hypothesis (high-inward and low-outward connectivity) was present in all individuals and found to be the most robust motif for identification of SOZs in 64/81 (79%) patients. This principal component demonstrated significant differences in SOZs compared to non-SOZs. While other motifs for identifying SOZs were present in other patients, they differed for each patient, suggesting that seizure networks are patient specific, but the ISH is present in nearly all networks. We discovered that a potentially suppressive motif based on the interictal suppression hypothesis was present in all patients, and it was the most robust motif for SOZs in 79% of patients. Each patient had additional motifs that further characterized SOZs, but these motifs were not common across all patients. This work has the potential to augment clinical identification of SOZs to improve epilepsy treatment.

Interictal stereo-electroencephalography features below 45 Hz are sufficient for correct localization of the epileptogenic zone and postsurgical outcome prediction

OBJECTIVE

Evidence suggests that the most promising results in interictal localization of the epileptogenic zone (EZ) are achieved by a combination of multiple stereo-electroencephalography (SEEG) biomarkers in machine learning models. These biomarkers usually include SEEG features calculated in standard frequency bands, but also high-frequency (HF) bands. Unfortunately, HF features require extra effort to record, store, and process. Here we investigate the added value of these HF features for EZ localization and postsurgical outcome prediction.

METHODS

In 50 patients we analyzed 30 min of SEEG recorded during non-rapid eye movement sleep and tested a logistic regression model with three different sets of features. The first model used broadband features (1-500 Hz); the second model used low-frequency features up to 45 Hz; and the third model used HF features above 65 Hz. The EZ localization by each model was evaluated by various metrics including the area under the precision-recall curve (AUPRC) and the positive predictive value (PPV). The differences between the models were tested by the Wilcoxon signed-rank tests and Cliff's Delta effect size. The differences in outcome predictions based on PPV values were further tested by the McNemar test.

RESULTS

The AUPRC score of the random chance classifier was .098. The models (broad-band, low-frequency, high-frequency) achieved median AUPRCs of .608, .582, and .522, respectively, and correctly predicted outcomes in 38, 38, and 33 patients. There were no statistically significant differences in AUPRC or any other metric between the three models. Adding HF features to the model did not have any additional contribution.

SIGNIFICANCE

Low-frequency features are sufficient for correct localization of the EZ and outcome prediction with no additional value when considering HF features. This finding allows significant simplification of the feature calculation process and opens the possibility of using these models in SEEG recordings with lower sampling rates, as commonly performed in clinical routines.

Identifying epileptogenic abnormality by decomposing intracranial EEG and MEG power spectra

BACKGROUND

Accurate identification of abnormal electroencephalographic (EEG) activity is pivotal for diagnosing and treating epilepsy. Recent studies indicate that decomposing brain activity into periodic (oscillatory) and aperiodic (trend across all frequencies) components can illuminate the drivers of spectral activity changes.

NEW METHODS

We analysed intracranial EEG (iEEG) data from 234 subjects, creating a normative map. This map was compared to a cohort of 63 patients with refractory focal epilepsy under consideration for neurosurgery. The normative map was computed using three approaches: (i) relative complete band power, (ii) relative band power with the aperiodic component removed, and (iii) the aperiodic exponent. Abnormalities were calculated for each approach in the patient cohort. We evaluated the spatial profiles, assessed their ability to localize abnormalities, and replicated the findings using magnetoencephalography (MEG).

RESULTS

Normative maps of relative complete band power and relative periodic band power exhibited similar spatial profiles, while the aperiodic normative map revealed higher exponent values in the temporal lobe. Abnormalities estimated through complete band power effectively distinguished between good and bad outcome patients. Combining periodic and aperiodic abnormalities enhanced performance, like the complete band power approach.

COMPARISON WITH EXISTING METHODS AND CONCLUSIONS

Sparing cerebral tissue with abnormalities in both periodic and aperiodic activity may result in poor surgical outcomes. Both periodic and aperiodic components do not carry sufficient information in isolation. The relative complete band power solution proved to be the most reliable method for this purpose. Future studies could investigate how cerebral location or pathology influences periodic or aperiodic abnormalities.

Gamma amplitude-envelope correlations are strongly elevated within hyperexcitable networks in focal epilepsy

Methods to quantify cortical hyperexcitability are of enormous interest for mapping epileptic networks in patients with focal epilepsy. We hypothesize that, in the resting state, cortical hyperexcitability increases firing-rate correlations between neuronal populations within seizure onset zones (SOZs). This hypothesis predicts that in the gamma frequency band (40-200 Hz), amplitude envelope correlations (AECs), a relatively straightforward measure of functional connectivity, should be elevated within SOZs compared to other areas. To test this prediction, we analyzed archived samples of interictal electrocorticographic (ECoG) signals recorded from patients who became seizure-free after surgery targeting SOZs identified by multiday intracranial recordings. We show that in the gamma band, AECs between nodes within SOZs are markedly elevated relative to those elsewhere. AEC-based node strength, eigencentrality, and clustering coefficient are also robustly increased within the SOZ with maxima in the low-gamma band (permutation test Z-scores > 8) and yield moderate discriminability of the SOZ using ROC analysis (maximal mean AUC ~ 0.73). By contrast to AECs, phase locking values (PLVs), a measure of narrow-band phase coupling across sites, and PLV-based graph metrics discriminate the seizure onset nodes weakly. Our results suggest that gamma band AECs may provide a clinically useful marker of cortical hyperexcitability in focal epilepsy.

Understanding epileptogenesis from molecules to network alteration

Epilepsy is characterized by recurrent seizures. Following an initial insult, a latent period precedes the onset of spontaneous seizures, a process referred to as epileptogenesis. This period plays a critical role in halting the progression toward epilepsy before the onset of abnormal molecular and network alterations. In this study, the fundamental concepts of epileptogenesis as well as the associated molecular and cellular targets are reviewed.

Spontaneous HFO Sequences Reveal Propagation Pathways for Precise Delineation of Epileptogenic Networks

Epilepsy, a neurological disorder affecting millions worldwide, poses great challenges in precisely delineating the epileptogenic zone - the brain region generating seizures - for effective treatment. High-frequency oscillations (HFOs) are emerging as promising biomarkers; however, the clinical utility is hindered by the difficulties in distinguishing pathological HFOs from non- epileptiform activities at single electrode and single patient resolution and understanding their dynamic role in epileptic networks. Here, we introduce an HFO-sequencing approach to analyze spontaneous HFOs traversing cortical regions in 40 drug-resistant epilepsy patients. This data- driven method automatically detected over 8.9 million HFOs, pinpointing pathological HFO- networks, and unveiled intricate millisecond-scale spatiotemporal dynamics, stability, and functional connectivity of HFOs in prolonged intracranial EEG recordings. These HFO sequences demonstrated a significant improvement in localization of epileptic tissue, with an 818.47% increase in concordance with seizure-onset zone (mean error: 2.92 mm), compared to conventional benchmarks. They also accurately predicted seizure outcomes for 90% AUC based on pre-surgical information using generalized linear models. Importantly, this mapping remained reliable even with short recordings (mean standard deviation: 3.23 mm for 30-minute segments). Furthermore, HFO sequences exhibited distinct yet highly repetitive spatiotemporal patterns, characterized by pronounced synchrony and predominant inward information flow from periphery towards areas involved in propagation, suggesting a crucial role for excitation-inhibition balance in HFO initiation and progression. Together, these findings shed light on the intricate organization of epileptic network and highlight the potential of HFO-sequencing as a translational tool for improved diagnosis, surgical targeting, and ultimately, better outcomes for vulnerable patients with drug-resistant epilepsy.

One Sentence Summary

Pathological fast brain oscillations travel like traffic along varied routes, outlining recurrently visited neural sites emerging as critical hotspots in epilepsy network.

Timing matters for accurate identification of the epileptogenic zone

OBJECTIVE

Interictal biomarkers of the epileptogenic zone (EZ) and their use in machine learning models open promising avenues for improvement of epilepsy surgery evaluation. Currently, most studies restrict their analysis to short segments of intracranial EEG (iEEG).

METHODS

We used 2381 hours of iEEG data from 25 patients to systematically select 5-minute segments across various interictal conditions. Then, we tested machine learning models for EZ localization using iEEG features calculated within these individual segments or across them and evaluated the performance by the area under the precision-recall curve (PRAUC).

RESULTS

On average, models achieved a score of 0.421 (the result of the chance classifier was 0.062). However, the PRAUC varied significantly across the segments (0.323-0.493). Overall, NREM sleep achieved the highest scores, with the best results of 0.493 in N2. When using data from all segments, the model performed significantly better than single segments, except NREM sleep segments.

CONCLUSIONS

The model based on a short segment of iEEG recording can achieve similar results as a model based on prolonged recordings. The analyzed segment should, however, be carefully and systematically selected, preferably from NREM sleep.

SIGNIFICANCE

Random selection of short iEEG segments may give rise to inaccurate localization of the EZ.

Correspondence between scalp-EEG and stereoelectroencephalography seizure-onset patterns in patients with MRI-negative drug-resistant focal epilepsy

OBJECTIVE

Our objective was to evaluate the relationship between scalp-EEG and stereoelectroencephalography (SEEG) seizure-onset patterns (SOP) in patients with MRI-negative drug-resistant focal epilepsy.

METHODS

We analyzed retrospectively 41 patients without visible lesion on brain MRI who underwent video-EEG followed by SEEG. We defined five types of SOPs on scalp-EEG and eight types on SEEG. We examined how various clinical variables affected scalp-EEG SOPs.

RESULTS

The most prevalent scalp SOPs were rhythmic sinusoidal activity (56.8%), repetitive epileptiform discharges (22.7%), and paroxysmal fast activity (15.9%). The presence of paroxysmal fast activity on scalp-EEG was always seen without delay from clinical onset and correlated with the presence of low-voltage fast activity in SEEG (sensitivity = 22.6%, specificity = 100%). The main factor explaining the discrepancy between the scalp and SEEG SOPs was the delay between clinical and scalp-EEG onset. There was a correlation between the scalp and SEEG SOPs when the scalp onset was simultaneous with the clinical onset (p = 0.026). A significant delay between clinical and scalp discharge onset was observed in 25% of patients and featured always with a rhythmic sinusoidal activity on scalp, corresponding to similar morphology of the discharge on SEEG. The presence of repetitive epileptiform discharges on scalp was associated with an underlying focal cortical dysplasia (sensitivity = 30%, specificity = 90%). There was no significant association between the scalp SOP and the epileptogenic zone location (deep or superficial), or surgical outcome.

SIGNIFICANCE

In patients with MRI-negative focal epilepsy, scalp SOP could suggest the SEEG SOP and some etiology (focal cortical dysplasia) but has no correlation with surgical prognosis. Scalp SOP correlates with the SEEG SOP in cases of simultaneous EEG and clinical onset; otherwise, scalp SOP reflects the propagation of the SEEG discharge.

PLAIN LANGUAGE SUMMARY

We looked at the correspondence between the electrical activity recorded during the start of focal seizure using scalp and intracerebral electrodes in patients with no visible lesion on MRI. If there is a fast activity on scalp, it reflects similar activity inside the brain. We found a good correspondence between scalp and intracerebral electrical activity for cases without significant delay between clinical and scalp electrical onset (seen in 75% of the cases we studied). Visualizing repetitive epileptic activity on scalp could suggest a particular cause of the epilepsy: a subtype of brain malformation called focal cortical dysplasia.

High-density electroencephalographic functional networks in genetic generalized epilepsy: Preserved whole-brain topology hides local reorganization

OBJECTIVE

Genetic generalized epilepsy (GGE) accounts for approximately 20% of adult epilepsy cases and is considered a disorder of large brain networks, involving both hemispheres. Most studies have not shown any difference in functional whole-brain network topology when compared to healthy controls. Our objective was to examine whether this preserved global network topology could hide local reorganizations that balance out at the global network level.

METHODS

We recorded high-density electroencephalograms from 20 patients and 20 controls, and reconstructed the activity of 118 regions. We computed functional connectivity in windows free of interictal epileptiform discharges in broad, delta, theta, alpha, and beta frequency bands, characterized the network topology, and used the Hub Disruption Index (HDI) to quantify the topological reorganization. We examined the generalizability of our results by reproducing a 25-electrode clinical system.

RESULTS

Our study did not reveal any significant change in whole-brain network topology among GGE patients. However, the HDI was significantly different between patients and controls in all frequency bands except alpha (p < .01, false discovery rate [FDR] corrected, d < -1), and accompanied by an increase in connectivity in the prefrontal regions and default mode network. This reorganization suggests that regions that are important in transferring the information in controls were less so in patients. Inversely, the crucial regions in patients are less so in controls. These findings were also found in delta and theta frequency bands when using 25 electrodes (p < .001, FDR corrected, d < -1).

SIGNIFICANCE

In GGE patients, the overall network topology is similar to that of healthy controls but presents a balanced local topological reorganization. This reorganization causes the prefrontal areas and default mode network to be more integrated and segregated, which may explain executive impairment associated with GGE. Additionally, the reorganization distinguishes patients from controls even when using 25 electrodes, suggesting its potential use as a diagnostic tool.

Network alterations in temporal lobe epilepsy during non-rapid eye movement sleep and wakefulness

OBJECTIVE

Investigate sleep and temporal lobe epilepsy (TLE) effects on brain networks derived from electroencephalography (EEG).

METHODS

High-density EEG was recorded during non-rapid eye movement (NREM) sleep stage 2 (N2) and wakefulness in 23 patients and healthy controls (HC). Epochs without epileptic discharges were source-reconstructed in 72 brain regions and connectivity was estimated. We calculated network integration and segregation at global (global efficiency, GE; average clustering coefficient, avgCC) and hemispheric level. These were compared between groups across frequency bands and correlated with the individual proportion of wakefulness- or sleep-related seizures.

RESULTS

At the global level, patients had higher delta GE, delta avgCC and theta avgCC than controls, irrespective of the vigilance state. During wakefulness, theta GE of patients was higher than controls and, for patients, theta GE during wakefulness was higher than during N2. Wake-to-sleep differences in TLE were notable only in the ipsilateral hemisphere. Only measures from wakefulness recordings correlated with the proportion of wakefulness- or sleep-related seizures.

CONCLUSIONS

TLE network alterations are more prominent during wakefulness and at lower frequencies. Increased integration and segregation suggest a pathological 'small world' configuration with a possible inhibitory role.

SIGNIFICANCE

Network alterations in TLE occur and are easier to detect during wakefulness.

Closed-loop electrical stimulation to prevent focal epilepsy progression and long-term memory impairment

Interictal epileptiform discharges (IEDs) are ubiquitously expressed in epileptic networks and disrupt cognitive functions. It is unclear whether addressing IED-induced dysfunction could improve epilepsy outcomes as most therapeutics target seizures. We show in a model of progressive hippocampal epilepsy that IEDs produce pathological oscillatory coupling which is associated with prolonged, hypersynchronous neural spiking in synaptically connected cortex and expands the brain territory capable of generating IEDs. A similar relationship between IED-mediated oscillatory coupling and temporal organization of IEDs across brain regions was identified in human subjects with refractory focal epilepsy. Spatiotemporally targeted closed-loop electrical stimulation triggered on hippocampal IED occurrence eliminated the abnormal cortical activity patterns, preventing spread of the epileptic network and ameliorating long-term spatial memory deficits in rodents. These findings suggest that stimulation-based network interventions that normalize interictal dynamics may be an effective treatment of epilepsy and its comorbidities, with a low barrier to clinical translation.

One-Sentence Summary

Targeted closed-loop electrical stimulation prevents spread of the epileptic network and ameliorates long-term spatial memory deficits.

Localizing seizure onset zone by a cortico-cortical evoked potentials-based machine learning approach in focal epilepsy

OBJECTIVE

We aimed to develop a new approach for identifying the localization of the seizure onset zone (SOZ) based on corticocortical evoked potentials (CCEPs) and to compare the connectivity patterns in patients with different clinical phenotypes.

METHODS

Fifty patients who underwent stereoelectroencephalography and CCEP procedures were included. Logistic regression was used in the model, and six CCEP metrics were input as features: root mean square of the first peak (N1RMS) and second peak (N2RMS), peak latency, onset latency, width duration, and area.

RESULTS

The area under the curve (AUC) for localizing the SOZ ranged from 0.88 to 0.93. The N1RMS values in the hippocampus sclerosis (HS) group were greater than that of the focal cortical dysplasia (FCD) IIa group (p < 0.001), independent of the distance between the recorded and stimulated sites. The sensitivity of localization was higher in the seizure-free group than in the non-seizure-free group (p = 0.036).

CONCLUSIONS

This new method can be used to predict the SOZ localization in various focal epilepsy phenotypes.

SIGNIFICANCE

This study proposed a machine-learning approach for localizing the SOZ. Moreover, we examined how clinical phenotypes impact large-scale abnormality of the epileptogenic networks.

Widespread slow oscillations support interictal epileptiform discharge networks in focal epilepsy

Interictal epileptiform discharges (IEDs) often co-occur across spatially-separated cortical regions, forming IED networks. However, the factors prompting IED propagation remain unelucidated. We hypothesized that slow oscillations (SOs) might facilitate IED propagation. Here, the amplitude and phase synchronization of SOs preceding propagating and non-propagating IEDs were compared in 22 patients with focal epilepsy undergoing intracranial electroencephalography (EEG) evaluation. Intracranial channels were categorized into the irritative zone (IZ) and normal zone (NOZ) regarding the presence of IEDs. During wakefulness, we found that pre-IED SOs within the IZ exhibited higher amplitudes for propagating IEDs than non-propagating IEDs (delta band: p = 0.001, theta band: p < 0.001). This increase in SOs was also concurrently observed in the NOZ (delta band: p = 0.04). Similarly, the inter-channel phase synchronization of SOs prior to propagating IEDs was higher than those preceding non-propagating IEDs in the IZ (delta band: p = 0.04). Through sliding window analysis, we observed that SOs preceding propagating IEDs progressively increased in amplitude and phase synchronization, while those preceding non-propagating IEDs remained relatively stable. Significant differences in amplitude occurred approximately 1150 ms before IEDs. During non-rapid eye movement (NREM) sleep, SOs on scalp recordings also showed higher amplitudes before intracranial propagating IEDs than before non-propagating IEDs (delta band: p = 0.006). Furthermore, the analysis of IED density around sleep SOs revealed that only high-amplitude sleep SOs demonstrated correlation with IED propagation. Overall, our study highlights that transient but widely distributed SOs are associated with IED propagation as well as generation in focal epilepsy during sleep and wakefulness, providing new insight into the EEG substrate supporting IED networks.

Delayed effective connectivity characterizes the epileptogenic zone during stereo-EEG

OBJECTIVE

Single-pulse electrical stimulations (SPES) can elicit normal and abnormal responses that might characterize the epileptogenic zone, including spikes, high-frequency oscillations and cortico-cortical evoked potentials (CCEPs). In this study, we investigate their association with the epileptogenic zone during stereoelectroencephalography (SEEG) in 28 patients with refractory focal epilepsy.

METHODS

Characteristics of CCEPs (distance-corrected or -uncorrected latency, amplitude and the connectivity index) and the occurrence of spikes and ripples were assessed. Responses within the epileptogenic zone and within the non-involved zone were compared using receiver operating characteristics curves and analysis of variance (ANOVA) either in all patients, patients with well-delineated epileptogenic zone, and patients older than 15 years old.

RESULTS

We found an increase in distance-corrected CCEPs latency after stimulation within the epileptogenic zone (area under the curve = 0.71, 0.72, 0.70, ANOVA significant after false discovery rate correction).

CONCLUSIONS

The increased distance-corrected CCEPs latency suggests that neuronal propagation velocity is altered within the epileptogenic network. This association might reflect effective connectivity changes at cortico-cortical or cortico-subcortico-cortical levels. Other responses were not associated with the epileptogenic zone, including the CCEPs amplitude, the connectivity index, the occurrences of induced ripples and spikes. The discrepancy with previous descriptions may be explained by different spatial brain sampling between subdural and depth electrodes.

SIGNIFICANCE

Increased distance-corrected CCEPs latency, indicating delayed effective connectivity, characterizes the epileptogenic zone. This marker could be used to help tailor surgical resection limits after SEEG.

Acute Effect of Vagus Nerve Stimulation in Patients with Drug-Resistant Epilepsy: A Preliminary Exploration via Stereoelectroencephalogram

As the pathophysiological mechanisms of vagus nerve stimulation (VNS) causing individual differences in the vagal ascending network remains unclear, stereoelectroencephalography (SEEG) provides a unique platform to explore the brain networks affected by VNS and helps to understand the anti-seizure mechanism of VNS more comprehensively. This study presents a preliminary exploration of the acute effect of VNS. SEEG signals were collected to assess the acute effect of VNS on neural synchronization in patients with drug-resistant epilepsy, especially in epileptogenic networks. The results show that the better the efficacy of VNS, the wider the spread of desynchronization assessed by weighted phase lag index at a high frequency band caused by VNS. Future studies should focus on the association between the change in synchronization and the efficacy of VNS, exploring the possibility of synchronization as a biomarker for patient screening and parameter programming.

Imaging and Stereotactic Electroencephalography Functional Networks to Guide Epilepsy Surgery

Epilepsy surgery is a potentially curative treatment of drug-resistant epilepsy that has remained underutilized both due to inadequate referrals and incomplete localization hypotheses. The complexity of patients evaluated for epilepsy surgery has increased, thus new approaches are necessary to treat these patients. The paradigm of epilepsy surgery has evolved to match this challenge, now considering the entire seizure network with the goal of disrupting it through resection, ablation, neuromodulation, or a combination. The network paradigm has the potential to aid in identification of the seizure network as well as treatment selection.

Interaction of interictal epileptiform activity with sleep spindles is associated with cognitive deficits and adverse surgical outcome in pediatric focal epilepsy

OBJECTIVE

Temporal coordination between oscillations enables intercortical communication and is implicated in cognition. Focal epileptic activity can affect distributed neural networks and interfere with these interactions. Refractory pediatric epilepsies are often accompanied by substantial cognitive comorbidity, but mechanisms and predictors remain mostly unknown. Here, we investigate oscillatory coupling across large-scale networks in the developing brain.

METHODS

We analyzed large-scale intracranial electroencephalographic recordings in children with medically refractory epilepsy undergoing presurgical workup (n = 25, aged 3-21 years). Interictal epileptiform discharges (IEDs), pathologic high-frequency oscillations (HFOs), and sleep spindles were detected. Spatiotemporal metrics of oscillatory coupling were determined and correlated with age, cognitive function, and postsurgical outcome.

RESULTS

Children with epilepsy demonstrated significant temporal coupling of both IEDs and HFOs to sleep spindles in discrete brain regions. HFOs were associated with stronger coupling patterns than IEDs. These interactions involved tissue beyond the clinically identified epileptogenic zone and were ubiquitous across cortical regions. Increased spatial extent of coupling was most prominent in older children. Poor neurocognitive function was significantly correlated with high IED-spindle coupling strength and spatial extent; children with strong pathologic interactions additionally had decreased likelihood of postoperative seizure freedom.

SIGNIFICANCE

Our findings identify pathologic large-scale oscillatory coupling patterns in the immature brain. These results suggest that such intercortical interactions could predict risk for adverse neurocognitive and surgical outcomes, with the potential to serve as novel therapeutic targets to restore physiologic development.

Identification of epileptic networks with graph convolutional network incorporating oscillatory activities and evoked synaptic responses

Stereoelectroencephalography (SEEG) offers unique neural data from in-depth brain structures with fine temporal resolutions to better investigate the origin of epileptic brain activities. Although oscillatory patterns from different frequency bands and functional connectivity computed from the SEEG datasets are employed to study the epileptic zones, direct electrical stimulation-evoked electrophysiological recordings of synaptic responses, namely cortical-cortical evoked potentials (CCEPs), from the same SEEG electrodes are not explored for the localization of epileptic zones. Here we proposed a two-stream model with unsupervised learning and graph convolutional network tailored to the SEEG and CCEP datasets in individual patients to perform localization of epileptic zones. We compared our localization results with the clinically marked electrode sites determined for surgical resections. Our model had good classification capability when compared to other state-of-the-art methods. Furthermore, based on our prediction results we performed group-level brain-area mapping analysis for temporal, frontal and parietal epilepsy patients and found that epileptic and non-epileptic brain networks were distinct in patients with different types of focal epilepsy. Our unsupervised data-driven model provides personalized localization analysis for the epileptic zones. The epileptic and non-epileptic brain areas disclosed by the prediction model provide novel insights into the network-level pathological characteristics of epilepsy.

Physiological and pathological neuronal connectivity in the living human brain based on intracranial EEG signals: the current state of research

Over the past decades, studies of human brain networks have received growing attention as the assessment and modelling of connectivity in the brain is a topic of high impact with potential application in the understanding of human brain organization under both physiological as well as various pathological conditions. Under specific diagnostic settings, human neuronal signal can be obtained from intracranial EEG (iEEG) recording in epilepsy patients that allows gaining insight into the functional organisation of living human brain. There are two approaches to assess brain connectivity in the iEEG-based signal: evaluation of spontaneous neuronal oscillations during ongoing physiological and pathological brain activity, and analysis of the electrophysiological cortico-cortical neuronal responses, evoked by single pulse electrical stimulation (SPES). Both methods have their own advantages and limitations. The paper outlines available methodological approaches and provides an overview of current findings in studies of physiological and pathological human brain networks, based on intracranial EEG recordings.

The role of quantitative markers in surgical prognostication after stereoelectroencephalography

OBJECTIVE

Stereoelectroencephalography (SEEG) is the reference method in the presurgical exploration of drug-resistant focal epilepsy. However, prognosticating surgery on an individual level is difficult. A quantified estimation of the most epileptogenic regions by searching for relevant biomarkers can be proposed for this purpose. We investigated the performances of ictal (Epileptogenicity Index, EI; Connectivity EI, cEI), interictal (spikes, high-frequency oscillations, HFO [80-300 Hz]; Spikes × HFO), and combined (Spikes × EI; Spikes × cEI) biomarkers in predicting surgical outcome and searched for prognostic factors based on SEEG-signal quantification.

METHODS

Fifty-three patients operated on following SEEG were included. We compared, using precision-recall, the epileptogenic zone quantified using different biomarkers (EZq ) against the visual analysis (EZC ). Correlations between the EZ resection rates or the EZ extent and surgical prognosis were analyzed.

RESULTS

EI and Spikes × EI showed the best precision against EZc (0.74; 0.70), followed by Spikes × cEI and cEI, whereas interictal markers showed lower precision. The EZ resection rates were greater in seizure-free than in non-seizure-free patients for the EZ defined by ictal biomarkers and were correlated with the outcome for EI and Spikes × EI. No such correlation was found for interictal markers. The extent of the quantified EZ did not correlate with the prognosis.

INTERPRETATION

Ictal or combined ictal-interictal markers overperformed the interictal markers both for detecting the EZ and predicting seizure freedom. Combining ictal and interictal epileptogenicity markers improves detection accuracy. Resection rates of the quantified EZ using ictal markers were the only statistically significant determinants for surgical prognosis.

Intracranial EEG Structure-Function Coupling and Seizure Outcomes After Epilepsy Surgery

BACKGROUND AND OBJECTIVES

Surgery is an effective treatment for drug-resistant epilepsy, which modifies the brain's structure and networks to regulate seizure activity. Our objective was to examine the relationship between brain structure and function to determine the extent to which this relationship affects the success of the surgery in controlling seizures. We hypothesized that a stronger association between brain structure and function would lead to improved seizure control after surgery.

METHODS

We constructed functional and structural brain networks in patients with drug-resistant focal epilepsy by using presurgery functional data from intracranial EEG (iEEG) recordings, presurgery and postsurgery structural data from T1-weighted MRI, and presurgery diffusion-weighted MRI. We quantified the relationship (coupling) between structural and functional connectivity by using the Spearman rank correlation and analyzed this structure-function coupling at 2 spatial scales: (1) global iEEG network level and (2) individual iEEG electrode contacts using virtual surgeries. We retrospectively predicted postoperative seizure freedom by incorporating the structure-function connectivity coupling metrics and routine clinical variables into a cross-validated predictive model.

RESULTS

We conducted a retrospective analysis on data from 39 patients who met our inclusion criteria. Brain areas implanted with iEEG electrodes had stronger structure-function coupling in seizure-free patients compared with those with seizure recurrence (p = 0.002, d = 0.76, area under the receiver operating characteristic curve [AUC] = 0.78 [95% CI 0.62-0.93]). Virtual surgeries on brain areas that resulted in stronger structure-function coupling of the remaining network were associated with seizure-free outcomes (p = 0.007, d = 0.96, AUC = 0.73 [95% CI 0.58-0.89]). The combination of global and local structure-function coupling measures accurately predicted seizure outcomes with a cross-validated AUC of 0.81 (95% CI 0.67-0.94). These measures were complementary to other clinical variables and, when included for prediction, resulted in a cross-validated AUC of 0.91 (95% CI 0.82-1.0), accuracy of 92%, sensitivity of 93%, and specificity of 91%.

DISCUSSION

Our study showed that the strength of structure-function connectivity coupling may play a crucial role in determining the success of epilepsy surgery. By quantitatively incorporating structure-function coupling measures and standard-of-care clinical variables into presurgical evaluations, we may be able to better localize epileptogenic tissue and select patients for epilepsy surgery.

CLASSIFICATION OF EVIDENCE

This is a Class IV retrospective case series showing that structure-function mapping may help determine the outcome from surgical resection for treatment-resistant focal epilepsy.

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.

Structure-function coupling increases during interictal spikes in temporal lobe epilepsy: A graph signal processing study

OBJECTIVE

Structure-function coupling remains largely unknown in brain disorders. We studied this coupling during interictal epileptic discharges (IEDs), using graph signal processing in temporal lobe epilepsy (TLE).

METHODS

We decomposed IEDs of 17 patients on spatial maps, i.e. network harmonics, extracted from a structural connectome. Harmonics were split in smooth maps (long-range interactions reflecting integration) and coarse maps (short-range interactions reflecting segregation) and were used to reconstruct the part of the signal coupled (Xc) and decoupled (Xd) from the structure, respectively. We analysed how Xc and Xd embed the IED energy over time, at global and regional level.

RESULTS

For Xc, the energy was smaller than for Xd before the IED onset (p < .001), but became larger around the first IED peak (p < .05, cluster 2, C2). Locally, the ipsilateral mesial regions were significantly coupled to the structure over the whole epoch. The ipsilateral hippocampus increased its coupling during C2 (p < .01).

CONCLUSIONS

At whole-brain level, segregation gives way to integrative processes during the IED. Locally, brain regions commonly involved in the TLE epileptogenic network increase their reliance on long-range couplings during IED (C2).

SIGNIFICANCE

In TLE, integration mechanisms prevail during the IED and are localized in the ipsilateral mesial temporal regions.

Does seizure propagate within or across intrinsic brain networks? An intracranial EEG study

BACKGROUND

Understanding the spatiotemporal propagation profiles of seizures is crucial for the preoperative assessment of epilepsy patients. The present study aimed to investigate whether seizures exhibit propagation patterns that align with intrinsic networks (INs).

METHODS

A quantitative analysis was conducted to examine ictal fast activity (IFA). The Epileptogenicity Index (EI) was employed to assess the epileptogenicity, spectral features, and temporal characteristics of IFA. Intra-network and inter-network comparisons were made regarding the IFA-related metrics. Additionally, the metrics were correlated with Euclidean distance. Network connection maps were generated to visualize seizures originating from different INs, allowing for comparisons between distinct groups.

RESULTS

Data for 81 seizures in 43 subjects were captured using stereoelectroencephalography implantation. Three metrics were compared: EI, time involvement (TI), and energy ratio index (ERI). Intra-network channels exhibited higher EI, earlier involvement of IFA, and stronger high-frequency energy. These findings were further validated through subgroup analyses stratified by neuropathology, seizure type, and seizure origination lobe. Correlation analyses revealed a negative association between distance and both EI and ERI, while distance exhibited a positive correlation with TI. Seizures originating from different INs exhibited varying propagation characteristics.

CONCLUSIONS

The study findings highlight the dominant role of intra-network dynamics over inter-network during seizure propagation. These results contribute to our understanding of seizure dynamics and their relationship with INs.

Pre-ictal fluctuation of EEG functional connectivity discriminates seizure phenotypes in mesial temporal lobe epilepsy

OBJECTIVE

We explored whether quantifiable differences between clinical seizures (CSs) and subclinical seizures (SCSs) occur in the pre-ictal state.

METHODS

We analyzed pre-ictal stereo-electroencephalography (SEEG) retrospectively across mesial temporal lobe epilepsy patients with recorded CSs and SCSs. Power spectral density and functional connectivity (FC) were quantified within and between the seizure onset zone (SOZ) and the early propagation zone (PZ), respectively. To evaluate the fluctuation of neural connectivity, FC variability was computed. Measures were further verified by a logistic regression model to evaluate their classification potentiality through the area under the receiver-operating-characteristics curve (AUC).

RESULTS

Fifty-four pre-ictal SEEG epochs (27 CSs and 27 SCSs) were selected among 14 patients. Within the SOZ, pre-ictal FC variability of CSs was larger than SCSs in 1-45 Hz during 30 seconds before seizure onset. Pre-ictal FC variability between the SOZ and PZ was larger in SCSs than CSs in 55-80 Hz within 1 minute before onset. Using these two variables, the logistic regression model achieved an AUC of 0.79 when classifying CSs and SCSs.

CONCLUSIONS

Pre-ictal FC variability within/between epileptic zones, not signal power or FC value, distinguished SCSs from CSs.

SIGNIFICANCE

Pre-ictal epileptic network stability possibly marks seizure phenotypes, contributing insights into ictogenesis and potentially helping seizure prediction.

Alteration of Excitation/Inhibition Imbalance in the Hippocampus and Amygdala of Drug-Resistant Epilepsy Patients Treated with Acute Vagus Nerve Stimulation

An imbalance between excitation (E) and inhibition (I) in the brain has been identified as a key pathophysiology of epilepsy over the years. The hippocampus and amygdala in the limbic system play a crucial role in the initiation and conduction of epileptic seizures and are often referred to as the transfer station and amplifier of seizure activities. Existing animal and imaging studies reveal that the hippocampus and amygdala, which are significant parts of the vagal afferent network, can be modulated in order to generate an antiepileptic effect. Using stereo-electroencephalography (SEEG) data, we examined the E/I imbalance in the hippocampus and amygdala of ten drug-resistant epilepsy children treated with acute vagus nerve stimulation (VNS) by estimating the 1/f power slope of hippocampal and amygdala signals in the range of 1-80 Hz. While the change in the 1/f power slope from VNS-BASE varied between different stimulation amplitudes and brain regions, it was more prominent in the hippocampal region. In the hippocampal region, we found a flatter 1/f power slope during VNS-ON in patients with good responsiveness to VNS under the optimal stimulation amplitude, indicating that the E/I imbalance in the region was improved. There was no obvious change in 1/f power slope for VNS poor responders. For VNS non-responders, the 1/f power slope slightly increased when the stimulation was applied. Overall, this study implies that the regulation of E/I imbalance in the epileptic brain, especially in the hippocampal region, may be an acute intracranial effect of VNS.

Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy

Normal brain functioning emerges from a complex interplay among regions forming networks. In epilepsy, these networks are disrupted causing seizures. Highly connected nodes in these networks are epilepsy surgery targets. Here, we assess whether functional connectivity (FC) using intracranial electroencephalography can quantify brain regions epileptogenicity and predict surgical outcome in children with drug resistant epilepsy (DRE). We computed FC between electrodes on different states (i.e. interictal without spikes, interictal with spikes, pre-ictal, ictal, and post-ictal) and frequency bands. We then estimated the electrodes' nodal strength. We compared nodal strength between states, inside and outside resection for good- (n = 22, Engel I) and poor-outcome (n = 9, Engel II-IV) patients, respectively, and tested their utility to predict the epileptogenic zone and outcome. We observed a hierarchical epileptogenic organization among states for nodal strength: lower FC during interictal and pre-ictal states followed by higher FC during ictal and post-ictal states (p < 0.05). We further observed higher FC inside resection (p < 0.05) for good-outcome patients on different states and bands, and no differences for poor-outcome patients. Resection of nodes with high FC was predictive of outcome (positive and negative predictive values: 47-100%). Our findings suggest that FC can discriminate epileptogenic states and predict outcome in patients with DRE.

Non-invasive mapping of epileptogenic networks predicts surgical outcome

Epilepsy is increasingly considered a disorder of brain networks. Studying these networks with functional connectivity can help identify hubs that facilitate the spread of epileptiform activity. Surgical resection of these hubs may lead patients who suffer from drug-resistant epilepsy to seizure freedom. Here, we aim to map non-invasively epileptogenic networks, through the virtual implantation of sensors estimated with electric and magnetic source imaging, in patients with drug-resistant epilepsy. We hypothesize that highly connected hubs identified non-invasively with source imaging can predict the epileptogenic zone and the surgical outcome better than spikes localized with conventional source localization methods (dipoles). We retrospectively analysed simultaneous high-density electroencephalography (EEG) and magnetoencephalography data recorded from 37 children and young adults with drug-resistant epilepsy who underwent neurosurgery. Using source imaging, we estimated virtual sensors at locations where intracranial EEG contacts were placed. On data with and without spikes, we computed undirected functional connectivity between sensors/contacts using amplitude envelope correlation and phase locking value for physiologically relevant frequency bands. From each functional connectivity matrix, we generated an undirected network containing the strongest connections within sensors/contacts using the minimum spanning tree. For each sensor/contact, we computed graph centrality measures. We compared functional connectivity and their derived graph centrality of sensors/contacts inside resection for good (n = 22, ILAE I) and poor (n = 15, ILAE II-VI) outcome patients, tested their ability to predict the epileptogenic zone in good-outcome patients, examined the association between highly connected hubs removal and surgical outcome and performed leave-one-out cross-validation to support their prognostic value. We also compared the predictive values of functional connectivity with those of dipoles. Finally, we tested the reliability of virtual sensor measures via Spearman's correlation with intracranial EEG at population- and patient-level. We observed higher functional connectivity inside than outside resection (P < 0.05, Wilcoxon signed-rank test) for good-outcome patients, on data with and without spikes across different bands for intracranial EEG and electric/magnetic source imaging and few differences for poor-outcome patients. These functional connectivity measures were predictive of both the epileptogenic zone and outcome (positive and negative predictive values ≥55%, validated using leave-one-out cross-validation) outperforming dipoles on spikes. Significant correlations were found between source imaging and intracranial EEG measures (0.4 ≤ rho ≤ 0.9, P < 0.05). Our findings suggest that virtual implantation of sensors through source imaging can non-invasively identify highly connected hubs in patients with drug-resistant epilepsy, even in the absence of frank epileptiform activity. Surgical resection of these hubs predicts outcome better than dipoles.

Epileptic Tissue Localization through Skewness-Based Functional Connectivity in the High-Frequency Band of Intracranial EEG

Functional connectivity analysis of intracranial electroencephalography (iEEG) plays an important role in understanding the mechanism of epilepsy and seizure dynamics. However, existing connectivity analysis is only suitable for low-frequency bands below 80 Hz. High-frequency oscillations (HFOs) and high-frequency activity (HFA) in the high-frequency band (80-500 Hz) are thought to be specific biomarkers in epileptic tissue localization. However, the transience in duration and variability of occurrence time and amplitudes of these events pose a challenge for conducting effective connectivity analysis. To deal with this problem, we proposed skewness-based functional connectivity (SFC) in the high-frequency band and explored its utility in epileptic tissue localization and surgical outcome evaluation. SFC comprises three main steps. The first step is the quantitative measurement of amplitude distribution asymmetry between HFOs/HFA and baseline activity. The second step is functional network construction on the basis of rank correlation of asymmetry across time. The third step is connectivity strength extraction from the functional network. Experiments were conducted in two separate datasets which consist of iEEG recordings from 59 patients with drug-resistant epilepsy. Significant difference (p<0.001) in connectivity strength was found between epileptic and non-epileptic tissue. Results were quantified via the receiver operating characteristic curve and the area under the curve (AUC). Compared with low-frequency bands, SFC demonstrated superior performance. With respect to pooled and individual epileptic tissue localization for seizure-free patients, AUCs were 0.66 (95% confidence interval (CI): 0.63-0.69) and (0.63 95% CI 0.56-0.71), respectively. For surgical outcome classification, the AUC was 0.75 (95% CI 0.59-0.85). Therefore, SFC can act as a promising assessment tool in characterizing the epileptic network and potentially provide better treatment options for patients with drug-resistant epilepsy.

Personalised virtual brain models in epilepsy

Individuals with drug-resistant focal epilepsy are candidates for surgical treatment as a curative option. Before surgery can take place, the patient must have a presurgical evaluation to establish whether and how surgical treatment might stop their seizures without causing neurological deficits. Virtual brains are a new digital modelling technology that map the brain network of a person with epilepsy, using data derived from MRI. This technique produces a computer simulation of seizures and brain imaging signals, such as those that would be recorded with intracranial EEG. When combined with machine learning, virtual brains can be used to estimate the extent and organisation of the epileptogenic zone (ie, the brain regions related to seizure generation and the spatiotemporal dynamics during seizure onset). Virtual brains could, in the future, be used for clinical decision making, to improve precision in localisation of seizure activity, and for surgical planning, but at the moment these models have some limitations, such as low spatial resolution. As evidence accumulates in support of the predictive power of personalised virtual brain models, and as methods are tested in clinical trials, virtual brains might inform clinical practice in the near future.

Resting-state SEEG-based brain network analysis for the detection of epileptic area

BACKGROUND

Most epilepsy research is based on interictal or ictal functional connectivity. However, prolonged electrode implantation may affect patients' health and the accuracy of epileptic zone identification. Brief resting-state SEEG recordings reduce the observation of epileptic discharges by reducing electrode implantation and other seizure-inducing interventions.

NEW METHOD

The location coordinates of SEEG in the brain were identified using CT and MRI. Based on undirected brain network connectivity, five functional connectivity measures and data feature vector centrality were calculated. Network connectivity was calculated from multiple perspectives of linear correlation, information theory, phase, and frequency, and the relative influence of nodes on network connectivity was considered. We investigated the potential value of resting-state SEEG for epileptic zone identification by comparing the differences between epileptic and non-epileptic zones, as well as the differences between patients with different surgical outcomes.

RESULTS

By comparing the centrality of brain network connectivity between epileptic and non-epileptic zones, we found significant differences in the distribution of brain networks between the two zones. There was a significant difference in brain network between patients with good surgical outcomes and those with poor surgical outcomes (p < 0.01). By combining support vector machines with static node importance, we predicted an AUC of 0.94 ± 0.08 for the epilepsy zone.

CONCLUSIONS AND SIGNIFICANCE

The results illustrated that nodes in epileptic zones are distinct from those in non-epileptic zones. Analysis of resting-state SEEG data and the importance of nodes in the brain network may contribute to identifying the epileptic zone and predicting the outcome.

Novel noninvasive identification of patient-specific epileptic networks in focal epilepsies: Linking single-photon emission computed tomography perfusion during seizures with resting-state magnetoencephalography dynamics

Single-photon emission computed tomography (SPECT) during seizures and magnetoencephalography (MEG) during the interictal state are noninvasive modalities employed in the localization of the epileptogenic zone in patients with drug-resistant focal epilepsy (DRFE). The present study aims to investigate whether there exists a preferentially high MEG functional connectivity (FC) among those regions of the brain that exhibit hyperperfusion or hypoperfusion during seizures. We studied MEG and SPECT data in 30 consecutive DRFE patients who had resective epilepsy surgery. We parcellated each ictal perfusion map into 200 regions of interest (ROIs) and generated ROI time series using source modeling of MEG data. FC between ROIs was quantified using coherence and phase-locking value. We defined a generalized linear model to relate the connectivity of each ROI, ictal perfusion z score, and distance between ROIs. We compared the coefficients relating perfusion z score to FC of each ROI and estimated the connectivity within and between resected and unresected ROIs. We found that perfusion z scores were strongly correlated with the FC of hyper-, and separately, hypoperfused ROIs across patients. High interictal connectivity was observed between hyperperfused brain regions inside and outside the resected area. High connectivity was also observed between regions of ictal hypoperfusion. Importantly, the ictally hypoperfused regions had a low interictal connectivity to regions that became hyperperfused during seizures. We conclude that brain regions exhibiting hyperperfusion during seizures highlight a preferentially connected interictal network, whereas regions of ictal hypoperfusion highlight a separate, discrete and interconnected, interictal network.

Basal temporal lobe epilepsy: SEEG electroclinical characteristics

OBJECTIVE

Temporal lobe epilepsy is the most common type of focal drug-resistant epilepsy. Seizures with predominant involvement of basal temporal regions (BTR) are not well characterized. In this stereo electroencephalography (SEEG) study, we aimed at describing the ictal networks involving BTR and the associated clinical features.

METHODS

We studied 24 patients explored with SEEG in our center with BTR sampling. We analyzed their seizures using a quantitative method: the "epileptogenicity index". Then we reported the features of the patients with maximal epileptogenicity within BTR, especially ictal network involved, ictal semiology and post-surgical outcome.

RESULTS

We found that rhinal cortex, parahippocampal cortex and posterior fusiform gyrus were the most epileptogenic structures within the BTR (mean EI: 0.57, 0.55, 0.54 respectively). Three main groups of epileptogenic zone organization were found: anterior (23% of total seizures) posterior (30%) and global (47%, both anterior and posterior). Contralateral spread was found in 35% of left seizures and 20% of right seizures. Naming deficit was more prevalent in left BTR (71% vs 29% in right seizures; p = 0.01) whereas automatic speech production was preferentially represented in right seizures (11% vs 54%; p = 0.001). Surgery was proposed for 11 patients (45.8%), leading to seizure freedom in 72% (Engel Class I). One patient presented post-operative permanent functional deficit.

CONCLUSION

Basal-temporal lobe epilepsy seems to be a specific entity among the temporal epilepsy spectrum with specific clinical characteristics. Resective surgery can be proposed with good outcomes in a significant proportion of patients and is safe provided that adequate language assessment has been preoperatively made.

Diverse nature of interictal oscillations: EEG-based biomarkers in epilepsy

Interictal electroencephalogram (EEG) patterns, including high-frequency oscillations (HFOs), interictal spikes (ISs), and slow wave activities (SWAs), are defined as specific oscillations between seizure events. These interictal oscillations reflect specific dynamic changes in network excitability and play various roles in epilepsy. In this review, we briefly describe the electrographic characteristics of HFOs, ISs, and SWAs in the interictal state, and discuss the underlying cellular and network mechanisms. We also summarize representative evidence from experimental and clinical epilepsy to address their critical roles in ictogenesis and epileptogenesis, indicating their potential as electrophysiological biomarkers of epilepsy. Importantly, we put forwards some perspectives for further research in the field.

Comparison of beamformer and ICA for dynamic connectivity analysis: A simultaneous MEG-SEEG study

Magnetoencephalography (MEG) is a powerful tool for estimating brain connectivity with both good spatial and temporal resolution. It is particularly helpful in epilepsy to characterize non-invasively the epileptic networks. However, using MEG to map brain networks requires solving a difficult inverse problem that introduces uncertainty in the activity localization and connectivity measures. Our goal here was to compare independent component analysis (ICA) followed by dipole source localization and the linearly constrained minimum-variance beamformer (LCMV-BF) for characterizing regions with interictal epileptic activity and their dynamic connectivity. After a simulation study, we compared ICA and LCMV-BF results with intracerebral EEG (stereotaxic EEG, SEEG) recorded simultaneously in 8 epileptic patients, which provide a unique 'ground truth' to which non-invasive results can be confronted. We compared the signal time courses extracted applying ICA and LCMV-BF on MEG data to that of SEEG, both for the actual signals and the dynamic connectivity computed using cross-correlation (evolution of links in time). With our simulations, we illustrated the different effect of the temporal and spatial correlation among sources on the two methods. While ICA was more affected by the temporal correlation but robust against spatial configurations, LCMV-BF showed opposite behavior. Moreover, ICA seems more suited to retrieve the simulated networks. In case of real patient data, good MEG/SEEG correlation and good localization were obtained in 6 out of 8 patients. In 4 of them ICA had the best performance (higher correlation, lower localization distance). In terms of dynamic connectivity, the evolution in time of the cross-correlation links could be retrieved in 5 patients out of 6, however, with more variable results in terms of correlation and distance. In two patients LCMV-BF had better results than ICA. In one patient the two methods showed equally good outcomes, and in the remaining two patients ICA performed best. In conclusion, our results obtained by exploiting simultaneous MEG/SEEG recordings suggest that ICA and LCMV-BF have complementary qualities for retrieving the dynamics of interictal sources and their network interactions.

Consciousness alteration in focal epilepsy is related to loss of signal complexity and information processing

Alteration of awareness is a main feature of focal epileptic seizures. In this work, we studied how the information contained in EEG signals was modified during temporal lobe seizures with altered awareness by using permutation entropy (PE) as a measure of the complexity of the signal. PE estimation was performed in thirty-six seizures of sixteen patients with temporal lobe epilepsy who underwent SEEG recordings. We tested whether altered awareness (based on the Consciousness Seizure Score) was correlated with a loss of signal complexity. We estimated global changes in PE as well as regional changes to gain insight into the mechanisms associated with awareness impairment. Our results reveal a positive correlation between the decrease of entropy and the consciousness score as well as the existence of a threshold on entropy that could discriminate seizures with no alteration of awareness from seizures with profound alteration of awareness. The loss of signal complexity was diffuse, extending bilaterally and to the associative cortices, in patients with profound alteration of awareness and limited to the temporal mesial structures in patients with no alteration of awareness. Thus PE is a promising tool to discriminate between the different subgroups of awareness alteration in TLE.

Ictal EEG desynchronization during low-voltage fast activity for prediction of surgical outcomes in focal epilepsy

OBJECTIVE

The authors investigated alterations in functional connectivity (FC) and EEG power during ictal onset patterns of low-voltage fast activity (LVFA) in drug-resistant focal epilepsy. They hypothesized that such changes would be useful to classify epilepsy surgical outcomes.

METHODS

In a cohort of 79 patients with drug-resistant focal epilepsy who underwent stereoelectroencephalography (SEEG) evaluation as well as resective surgery, FC changes during the peri-LVFA period were measured using nonlinear regression (h2) and power spectral properties within/between three regions: the seizure onset zone (SOZ), early propagation zone (PZ), and noninvolved zone (NIZ). Desynchronization and power desynchronization h2 indices were calculated to assess the degree of EEG desynchronization during LVFA. Multivariate logistic regression was employed to control for confounding factors. Finally, receiver operating characteristic curves were generated to evaluate the performance of desynchronization indices in predicting surgical outcome.

RESULTS

Fifty-three patients showed ictal LVFA and distinct zones of the SOZ, PZ, and NIZ. Among them, 39 patients (73.6%) achieved seizure freedom by the final follow-up. EEG desynchronization, measured by h2 analysis, was found in the seizure-free group during LVFA: FC decreased within the SOZ and between regions compared with the pre-LVFA and post-LVFA periods. In contrast, the non-seizure-free group showed no prominent EEG desynchronization. The h2 desynchronization index, but not the power desynchronization index, enabled classification of seizure-free versus non-seizure-free patients after resective surgery.

CONCLUSIONS

EEG desynchronization during the peri-LVFA period, measured by within-zone and between-zone h2 analysis, may be helpful for identifying patients with favorable postsurgical outcomes and also may potentially improve epileptogenic zone identification in the future.

Interictal Gamma Event Connectivity Differentiates the Seizure Network and Outcome in Patients after Temporal Lobe Epilepsy Surgery

Studies of interictal EEG functional connectivity in the epileptic brain seek to identify abnormal interactions between brain regions involved in generating seizures, which clinically often is defined by the seizure onset zone (SOZ). However, there is evidence for abnormal connectivity outside the SOZ (NSOZ), and removal of the SOZ does not always result in seizure control, suggesting, in some cases, that the extent of abnormal connectivity indicates a larger seizure network than the SOZ. To better understand the potential differences in interictal functional connectivity in relation to the seizure network and outcome, we computed event connectivity in the theta (4-8 Hz, ThEC), low-gamma (30-55 Hz, LGEC), and high-gamma (65-95 Hz, HGEC) bands from interictal depth EEG recorded in surgical patients with medication-resistant seizures suspected to begin in the temporal lobe. Analysis finds stronger LGEC and HGEC in SOZ than NSOZ of seizure-free (SF) patients (p = 1.10e-9, 0.0217), but no difference in not seizure-free (NSF) patients. There were stronger LGEC and HGEC between mesial and lateral temporal SOZ of SF than NSF patients (p = 0.00114, 0.00205), and stronger LGEC and ThEC in NSOZ of NSF than SF patients (p = 0.0089, 0.0111). These results show that event connectivity is sensitive to differences in the interactions between regions in SOZ and NSOZ and SF and NSF patients. Patients with differential strengths in event connectivity could represent a well-localized seizure network, whereas an absence of differences could indicate a larger seizure network than the one localized by the SOZ and higher likelihood for seizure recurrence.

Interictal Functional Connectivity in Focal Refractory Epilepsies Investigated by Intracranial EEG

Introduction: Focal epilepsies are diseases of neuronal excitability affecting macroscopic networks of cortical and subcortical neural structures. These networks ("epileptogenic networks") can generate pathological electrophysiological activities during seizures, and also between seizures (interictal period). Many works attempt to describe these networks by using quantification methods, particularly based on the estimation of statistical relationships between signals produced by brain regions, namely functional connectivity (FC). Results: FC has been shown to be greatly altered during seizures and in the immediate peri-ictal period. An increasing number of studies have shown that FC is also altered during the interictal period depending on the degree of epileptogenicity of the structures. Furthermore, connectivity values could be correlated with other clinical variables including surgical outcome. Significance: This leads to a conceptual change and to consider epileptic areas as both hyperexcitable and abnormally connected. These data open the door to the use of interictal FC as a marker of epileptogenicity and as a complementary tool for predicting the effect of surgery. Aim: In this article, we review the available data concerning interictal FC estimated from intracranial electroencephalograhy (EEG) in focal epilepsies and discuss it in the light of data obtained from other modalities (EEG imaging) and modeling studies.

Virtual MEG sensors based on beamformer and independent component analysis can reconstruct epileptic activity as measured on simultaneous intracerebral recordings

The prevailing gold standard for presurgical determination of epileptogenic brain networks is intracerebral EEG, a potent yet invasive approach. Magnetoencephalography (MEG) is a state-of-the art non-invasive method for investigating epileptiform discharges. However, it is not clear at what level the precision offered by MEG can reach that of SEEG. Here, we present a strategy for non-invasively retrieving the constituents of the interictal network, with high spatial and temporal precision. Our method is based on MEG and a combination of spatial filtering and independent component analysis (ICA). We validated this approach in twelve patients with drug-resistant focal epilepsy, thanks to the unprecedented ground truth provided by simultaneous recordings of MEG and SEEG. A minimum variance adaptive beamformer estimated the source time series and ICA was used to further decompose these time series into network constituents (MEG-ICs), each having a time series (virtual electrode) and a topography (spatial distribution of amplitudes in the brain). We show that MEG has a considerable sensitivity of 0.80 and 0.84 and a specificity of 0.93 and 0.91 for reconstructing deep and superficial sources, respectively, when compared to the ground truth (SEEG). For each epileptic MEG-IC (n = 131), we found at least one significantly correlating SEEG contact close to zero lag after correcting for multiple comparisons. All the patients except one had at least one epileptic component that was highly correlated (Spearman rho>0.3) with that of SEEG traces. MEG-ICs correlated well with SEEG traces. The strength of correlation coefficients did not depend on the depth of the SEEG contacts or the clinical outcome of the patient. A significant proportion of the MEG-ICs (n = 83/131) were localized in proximity with their maximally correlating SEEG, within a mean distance of 20±12.18mm. Our research is the first to validate the MEG-retrieved beamformer IC sources against SEEG-derived ground truth in a simultaneous MEG-SEEG framework. Observations from the present study suggest that non-invasive MEG source components may potentially provide additional information, comparable to SEEG in a number of instances.

EEG Markers of Treatment Resistance in Idiopathic Generalized Epilepsy: From Standard EEG Findings to Advanced Signal Analysis

Idiopathic generalized epilepsy (IGE) represents a common form of epilepsy in both adult and pediatric epilepsy units. Although IGE has been long considered a relatively benign epilepsy syndrome, a remarkable proportion of patients could be refractory to treatment. While some clinical prognostic factors have been largely validated among IGE patients, the impact of routine electroencephalography (EEG) findings in predicting drug resistance is still controversial and a growing number of authors highlighted the potential importance of capturing the sleep state in this setting. In addition, the development of advanced computational techniques to analyze EEG data has opened new opportunities in the identification of reliable and reproducible biomarkers of drug resistance in IGE patients. In this manuscript, we summarize the EEG findings associated with treatment resistance in IGE by reviewing the results of studies considering standard EEGs, 24-h EEG recordings, and resting-state protocols. We discuss the role of 24-h EEG recordings in assessing seizure recurrence in light of the potential prognostic relevance of generalized fast discharges occurring during sleep. In addition, we highlight new and promising biomarkers as identified by advanced EEG analysis, including hypothesis-driven functional connectivity measures of background activity and data-driven quantitative findings revealed by machine learning approaches. Finally, we thoroughly discuss the methodological limitations observed in existing studies and briefly outline future directions to identify reliable and replicable EEG biomarkers in IGE patients.

Metabolic connectivity is associated with seizure outcome in surgically treated temporal lobe epilepsies: A 18F-FDG PET seed correlation analysis

¹⁸F-FDG PET provides high sensitivity for the pre-surgical assessment of drug-resistant temporal lobe epilepsy (TLE). However, little is known about the metabolic connectivity of epileptogenic networks involved. This study therefore aimed to evaluate the association between metabolic connectivity and seizure outcome in surgically treated TLE.

METHODS

The study included 107 right-handed patients that had undergone a presurgical interictal ¹⁸F-FDG PET assessment followed by an anterior temporal lobectomy and were classified according to seizure outcome 2 years after surgery. Metabolic connectivity was evaluated by seed correlation analysis in left and right epilepsy patients with a Class Engel IA or > IA outcome and compared to age-, sex- and handedness-matched healthy controls.

RESULTS

Increased metabolic connectivity was observed in the >IA compared to the IA group within the operated temporal lobe (respective clusters of 7.5 vs 3.3 cm³ and 2.6 cm³ vs 2.2 cm³ in left and right TLE), and to a lower extent with the contralateral temporal lobe (1.2 vs 0.7 cm³ and 1.7 cm³ vs 0.7 cm³ in left and right TLE). Seed correlations provided added value for the estimated individual performance of seizure outcome over the group comparisons in left TLE (AUC of 0.74 vs 0.67).

CONCLUSION

Metabolic connectivity is associated with outcome in surgically treated TLE with a strengthened epileptogenic connectome in patients with non-free-seizure outcomes. The added value of seed correlation analysis in left TLE underlines the importance of evaluating metabolic connectivity in network related diseases.

Addressing spatial bias in intracranial EEG functional connectivity analyses for epilepsy surgical planning

Objective.To determine the effect of epilepsy on intracranial electroencephalography (EEG) functional connectivity, and the ability of functional connectivity to localize the seizure onset zone (SOZ), controlling for spatial biases.Approach.We analyzed intracranial EEG data from patients with drug-resistant epilepsy admitted for pre-surgical planning. We calculated intracranial EEG functional networks and determined whether changes in functional connectivity lateralized the SOZ using a spatial subsampling method to control for spatial bias. We developed a 'spatial null model' to localize the SOZ electrode using only spatial sampling information, ignoring EEG data. We compared the performance of this spatial null model against models incorporating EEG functional connectivity and interictal spike rates.Main results.About 110 patients were included in the study, although the number of patients differed across analyses. Controlling for spatial sampling, the average connectivity was lower in the SOZ region relative to the same anatomic region in the contralateral hemisphere. A model using intra-hemispheric connectivity accurately lateralized the SOZ (average accuracy 75.5%). A spatial null model incorporating spatial sampling information alone achieved moderate accuracy in classifying SOZ electrodes (mean AUC = 0.70, 95% CI 0.63-0.77). A model incorporating intracranial EEG functional connectivity and spike rate data further outperformed this spatial null model (AUC 0.78,p= 0.002 compared to spatial null model). However, a model incorporating functional connectivity without spike rate data did not significantly outperform the null model (AUC 0.72,p= 0.38).Significance.Intracranial EEG functional connectivity is reduced in the SOZ region, and interictal data predict SOZ electrode localization and laterality, however a predictive model incorporating functional connectivity without interictal spike rates did not significantly outperform a spatial null model. We propose constructing a spatial null model to provide an estimate of the pre-implant hypothesis of the SOZ, and to serve as a benchmark for further machine learning algorithms in order to avoid overestimating model performance because of electrode sampling alone.

Cortico-cortical and thalamo-cortical connectivity during non-REM and REM sleep: Insights from intracranial recordings in humans

OBJECTIVE

To study changes of thalamo-cortical and cortico-cortical connectivity during wakefulness, non-Rapid Eye Movement (non-REM) sleep, including N2 and N3 stages, and REM sleep, using stereoelectroencephalography (SEEG) recording in humans.

METHODS

We studied SEEG recordings of ten patients during wakefulness, non-REM sleep and REM sleep, in seven brain regions of interest including the thalamus. We calculated directed and undirected functional connectivity using a measure of non-linear correlation coefficient h².

RESULTS

The thalamus was more connected to other brain regions during N2 stage and REM sleep than during N3 stage during which cortex was more connected than the thalamus. We found two significant directed links: the first from the prefrontal region to the lateral parietal region in the delta band during N3 sleep and the second from the thalamus to the insula during REM sleep.

CONCLUSIONS

These results showed that cortico-cortical connectivity is more prominent in N3 stage than in N2 and REM sleep. During REM sleep we found significant thalamo-insular connectivity, with a driving role of the thalamus.

SIGNIFICANCE

We found a pattern of cortical connectivity during N3 sleep concordant with antero-posterior traveling slow waves. The thalamus seemed particularly involved as a hub of connectivity during REM sleep.

Epileptic spasms are associated with increased stereo-electroencephalography derived functional connectivity in tuberous sclerosis complex

OBJECTIVE

Epileptic spasms (ES) are common in tuberous sclerosis complex (TSC). However, the underlying network alterations and relationship with epileptogenic tubers are poorly understood. We examined interictal functional connectivity (FC) using stereo-electroencephalography (SEEG) in patients with TSC to investigate the relationship between tubers, epileptogenicity, and ES.

METHODS

We analyzed 18 patients with TSC who underwent SEEG (mean age = 11.5 years). The dominant tuber (DT) was defined as the most epileptogenic tuber using the epileptogenicity index. Epileptogenic zone (EZ) organization was quantitatively separated into focal (isolated DT) and complex (all other patterns). Using a 20-min interictal recording, FC was estimated with nonlinear regression, h² . We calculated (1) intrazone FC within all sampled tubers and normal-appearing cortical zones, respectively; and (2) interzone FC involving connections between DT, other tubers, and normal cortex. The relationship between FC and (1) presence of ES as a current seizure type at the time of SEEG, (2) EZ organization, and (3) epileptogenicity was analyzed using a mixed generalized linear model. Spike rate and distance between zones were considered in the model as covariates.

RESULTS

Six patients had ES as a current seizure type at time of SEEG. ES patients had a greater number of tubers with a fluid-attenuated inversion recovery hypointense center (p < .001), and none had TSC1 mutations. The presence of ES was independently associated with increased FC within both intrazone (p = .033) and interzone (p = .011) networks. Post hoc analyses identified that increased FC was associated with ES across tuber and nontuber networks. EZ organization and epileptogenicity biomarkers were not associated with FC.

SIGNIFICANCE

Increased cortical synchrony among both tuber and nontuber networks is characteristic of patients with ES and independent of both EZ organization and tuber epileptogenicity. This further supports the prospect of FC biomarkers aiding treatment paradigms in TSC.

Latest Views on the Mechanisms of Action of Surgically Implanted Cervical Vagal Nerve Stimulation in Epilepsy

BACKGROUND

Vagus nerve stimulation (VNS) is approved as an adjunctive treatment for drug-resistant epilepsy. Although there is a substantial amount of literature aiming at unraveling the mechanisms of action of VNS in epilepsy, it is still unclear how the cascade of events triggered by VNS leads to its antiepileptic effect.

OBJECTIVE

In this review, we integrated available peer-reviewed data on the effects of VNS in clinical and experimental research to identify those that are putatively responsible for its therapeutic effect. The topic of transcutaneous VNS will not be covered owing to the current lack of data supporting the differences and commonalities of its mechanisms of action in relation to invasive VNS.

SUMMARY OF THE MAIN FINDINGS

There is compelling evidence that the effect is obtained through the stimulation of large-diameter afferent myelinated fibers that project to the solitary tract nucleus, then to the parabrachial nucleus, which in turn alters the activity of the limbic system, thalamus, and cortex. VNS-induced catecholamine release from the locus coeruleus in the brainstem plays a pivotal role. Functional imaging studies tend to point toward a common vagal network that comes into play, made up of the amygdalo-hippocampal regions, left thalamus, and insular cortex.

CONCLUSIONS

Even though some crucial pieces are missing, neurochemical, molecular, cellular, and electrophysiological changes occur within the vagal afferent network at three main levels (the brainstem, the limbic system [amygdala and hippocampus], and the cortex). At this final level, VNS notably alters functional connectivity, which is known to be abnormally high within the epileptic zone and was shown to be significantly decreased by VNS in responders. The effect of crucial VNS parameters such as frequency or current amplitude on functional connectivity metrics is of utmost importance and requires further investigation.

Transcranial current stimulation in epilepsy: A systematic review of the fundamental and clinical aspects

Purpose

Transcranial electrical current stimulation (tES or tCS, as it is sometimes referred to) has been proposed as non-invasive therapy for pharmacoresistant epilepsy. This technique, which includes direct current (tDCS) and alternating current (tACS) stimulation involves the application of weak currents across the cortex to change cortical excitability. Although clinical trials have demonstrated the therapeutic efficacy of tES, its specific effects on epileptic brain activity are poorly understood. We sought to summarize the clinical and fundamental effects underlying the application of tES in epilepsy.

Methods

A systematic review was performed in accordance with the PRISMA guidelines. A database search was performed in PUBMED, MEDLINE, Web of Science and Cochrane CENTRAL for articles corresponding to the keywords “epilepsy AND (transcranial current stimulation OR transcranial electrical stimulation)”.

Results

A total of 56 studies were included in this review. Through these records, we show that tDCS and tACS epileptic patients are safe and clinically relevant techniques for epilepsy. Recent articles reported changes of functional connectivity in epileptic patients after tDCS. We argue that tDCS may act by affecting brain networks, rather than simply modifying local activity in the targeted area. To explain the mechanisms of tES, various cellular effects have been identified. Among them, reduced cell loss, mossy fiber sprouting, and hippocampal BDNF protein levels. Brain modeling and human studies highlight the influence of individual brain anatomy and physiology on the electric field distribution. Computational models may optimize the stimulation parameters and bring new therapeutic perspectives.

Conclusion

Both tDCS and tACS are promising techniques for epilepsy patients. Although the clinical effects of tDCS have been repeatedly assessed, only one clinical trial has involved a consistent number of epileptic patients and little knowledge is present about the clinical outcome of tACS. To fill this gap, multicenter studies on tES in epileptic patients are needed involving novel methods such as personalized stimulation protocols based on computational modeling. Furthermore, there is a need for more in vivo studies replicating the tES parameters applied in patients. Finally, there is a lack of clinical studies investigating changes in intracranial epileptiform discharges during tES application, which could clarify the nature of tES-related local and network dynamics in epilepsy.

Interictal epileptiform discharges show distinct spatiotemporal and morphological patterns across wake and sleep

Presurgical evaluation of mesial temporal and neocortical focal pharmacoresistant epilepsy patients using intracranial EEG recordings has led to the generation of extensive data on interictal epileptiform discharges, located within or remotely from seizure onset zones. In this study, we used this data to investigate how interictal epileptiform discharges are modulated and how their spatial distribution changes during wake and sleep and analysed the relationship between these discharge events and seizure onset zones. Preoperative evaluation data from 11 adult patients with focal pharmacoresistant epilepsy were extracted from the Epilepsiae database. Interictal epileptiform discharges were automatically detected during wakefulness and over several hours of continuous seizure-free sleep (total duration of EEG recordings:106.7 h; mean per patient: 9.7 h), and analysed across four brain areas (mesial temporal, lateral neocortical, basal cortical and the temporal pole). Sleep stages were classified manually from scalp EEG. Discharge events were characterized according to their rate and morphology (amplitude, sharpness and duration). Eight patients had a seizure onset zone over mesial areas and three patients over lateral neocortical areas. Overall, discharge rates varied across brain areas during wakefulness and sleep [wake/sleep stages × brain areas interaction; Wald χ 2(df = 6) = 31.1, P < 0.0001]. N2-N3 non-rapid eye movement sleep increased interictal epileptiform discharges in mesial areas compared with wakefulness and rapid eye movement sleep (P < 0.0001), and to other areas (P < 0.0001 for all comparisons). This mesial pattern was observed both within and outside of seizure onset zones. During wakefulness, the rate of interictal epileptiform discharges was significantly higher than during N2-N3 non-rapid eye movement sleep (P = 0.04), and rapid eye movement sleep (P = 0.01) in lateral neocortical areas (referred to as lateral neocortical pattern), a finding that was more pronounced in seizures onset zones (P = 0.004). The morphological characteristics of the discharge events were modulated during wakefulness and sleep stages across brain areas. The effect of seizure onset zones on discharge morphology was conditioned by brain area and was particularly marked in temporal pole areas. Our analysis of discharge patterns in relation to cerebral localization, vigilance state and the anatomical affiliation of seizure onset zones revealed the global and local aspects of the complex relationship between interictal discharges, sleep and seizure onset zones. This novel approach may lead to a better understanding of cognitive decline and responses to therapy, as well as to adaptation of surgical interventions for epileptic patients.