High-frequency oscillations are not better biomarkers of epileptogenic tissues than spikes

Application of HFO and scaling analysis of neuronal oscillations in the presurgical evaluation of focal epilepsy

PURPOSE

To explore the utility of high frequency oscillations (HFO) and long-range temporal correlations (LRTCs) in preoperative assessment of epilepsy.

METHODS

MEG ripples were detected in 59 drug-resistant epilepsy patients, comprising 5 with parietal lobe epilepsy (PLE), 21 with frontal lobe epilepsy (FLE), 14 with lateral temporal lobe epilepsy (LTLE), and 19 with mesial temporal lobe epilepsy (MTLE) to identify the epileptogenic zone (EZ). The results were compared with clinical MEG reports and resection area. Subsequently, LRTCs were quantified at the source-level by detrended fluctuation analysis (DFA) and life/waiting -time at 5 bands for 90 cerebral cortex regions. The brain regions with larger DFA exponents and standardized life-waiting biomarkers were compared with the resection results.

RESULTS

Compared to MEG sensor-level data, ripple sources were more frequently localized within the resection area. Moreover, source-level analysis revealed a higher proportion of DFA exponents and life-waiting biomarkers with relatively higher rankings, primarily distributed within the resection area (p<0.01). Moreover, these two LRCT indices across five distinct frequency bands correlated with EZ.

CONCLUSION

HFO and source-level LRTCs are correlated with EZ. Integrating HFO and LRTCs may be an effective approach for presurgical evaluation of epilepsy.

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.

Optimizing automated detection of high frequency oscillations using visual markings does not improve SOZ localization

OBJECTIVE

High frequency oscillations (HFOs) are a biomarker of the seizure onset zone (SOZ) and can be visually or automatically detected. In theory, one can optimize an automated algorithm's parameters to maximize SOZ localization accuracy; however, there is no consensus on whether or how this should be done. Therefore, we optimized an automated detector using visually identified HFOs and evaluated the impact on SOZ localization accuracy.

METHODS

We detected HFOs in intracranial EEG from 20 patients with refractory epilepsy from two centers using (1) unoptimized automated detection, (2) visual identification, and (3) automated detection optimized to match visually detected HFOs.

RESULTS

SOZ localization accuracy based on HFO rate was not significantly different between the three methods. Across patients, visually optimized detector settings varied, and no single set of settings produced universally accurate SOZ localization. Exploratory analysis suggests that, for many patients, detection settings exist that would improve SOZ localization.

CONCLUSIONS

SOZ localization accuracy was similar for all three methods, was not improved by visually optimizing detector settings, and may benefit from patient-specific parameter optimization.

SIGNIFICANCE

Visual HFO marking is laborious, and optimizing automated detection using visual markings does not improve localization accuracy. New patient-specific detector optimization methods are needed.

Discovering Neurophysiological Characteristics of Pathological High-Frequency Oscillations in Epilepsy with an Explainable Deep Generative Model

Interictal high-frequency oscillation (HFO) is a promising biomarker of the epileptogenic zone (EZ). However, objective definitions to distinguish between pathological and physiological HFOs have remained elusive, impeding HFOs' clinical applications. We employed self-supervised deep generative variational autoencoders to learn such discriminative HFO features directly from their morphologies in a data-driven manner. We studied a large retrospective cohort of 185 patients who underwent intracranial monitoring and analyzed 686,410 candidate HFO events collected from 18,265 brain contacts across diverse brain regions. The model automatically clustered HFOs into distinct morphological groups in the latent space. One cluster consisted of putative morphologically defined pathological HFOs (mpHFOs): HFOs in that cluster were observed to be associated with spikes and exhibited high signal intensity both in the HFO band (>80 Hz) at detection and in the sub-HFO band (10-80 Hz) surrounding the detection and were primarily localized in the seizure onset zone (SOZ). Moreover, resection of brain regions based on a higher prevalence of interictal mpHFOs better predicted postoperative seizure outcomes than current clinical standards based on SOZ removal. Our self-supervised, explainable, deep generative model distills pathological HFOs and thus potentially helps delineate the EZ purely from interictal intracranial EEG data.

Epilepsy lesion localization method based on brain function network

Objective

In the past, the localization of seizure onset zone (SOZ) primarily relied on traditional EEG signal analysis methods. However, due to their limited spatial and temporal resolution, accurately pinpointing neural activity was challenging, thereby restricting their clinical applicability. Compared with traditional EEG signals, SEEG signals have superior spatial and temporal resolution, and can more accurately record neural activity near epileptic foci, making them better suited for studying SOZ. In addition, the traditional EEG signal analysis methods still have limitations, mainly focusing on the analysis of local signal features, while ignoring the complexity and interconnection of the overall brain network. How to more accurately locate SOZ is still not well resolved. The purpose of this study is to develop an effective positioning method for more accurate positioning.

Method

To overcome these limitations, this study proposed a model integrating brain functional network analysis with nonlinear dynamics. We utilized weighted phase lag index (WPLI) to construct brain functional network, epilepic network connectivity strength (ENCS) as the feature, and introduced persistence entropy (PE) for feature fusion, subsequently employing support vector machine (SVM) classification.

Results

The proposed method was verified on the HUP-iEEG dataset, our solution identified the SOZ with 0.9440 accuracy, 0.9848 precision, 0.8974 recall rate, 0.9340 F1 score and 0.9697 area under the ROC curve across patients, which outperforms the existing approaches. It exhibits a 2.30 percentage point enhancement in localisation accuracy along with a 2.97 percentage points in AUC compared to others.

Conclusion

Our method consider the interactions between nodes in brain network connections, as well as the inherent nonlinear and non-stationary properties of neural signals, to be more robust.

Neuronal current imaging of epileptic activity: An MRI study in patients with a first unprovoked epileptic seizure

OBJECTIVE

This study evaluates the performance of the novel MRI sequence stimulus-induced rotary saturation (SIRS) to map responses to interictal epileptic activity in the human cortex. Spin-lock pulses have been applied to indirectly detect neuronal activity through magnetic field perturbations. Following initial reports about the feasibility of the method in humans and animals with epilepsy, we aimed to investigate the diagnostic yield of spin-lock MR pulses in comparison with scalp-EEG in first seizure patients.

METHODS

We employed a novel method for measurements of neuronal activity through the detection of a resonant oscillating field, stimulus-induced rotary saturation contrast (SIRS) at spin-lock frequencies of 120 and 240 Hz acquired at a single 3T MRI system. Within a prospective observational study, we conducted SIRS experiments in 55 patients within 7 days after a suspected first unprovoked epileptic seizure and 61 healthy control subjects. In this study, we report on the analysis of data from a single 3T MRI system, encompassing 35 first seizure patients and 31 controls.

RESULTS

The SIRS method was applicable in all patients and healthy controls at frequencies of 120 and 240 Hz. We did not observe any significant age- or sex-related differences. Specificity of SIRS at 120 Hz was 90.3% and 93.5% at 240 Hz. Sensitivity was 17.1% at 120 Hz and 40.0% at 240 Hz.

SIGNIFICANCE

SIRS targets neuronal oscillating magnetic fields in patients with epilepsy. The coupling of presaturated spins to epilepsy-related magnetic field perturbations may serve as a-at this stage experimental-diagnostic test in first seizure patients to complement EEG findings as a standard screening test.

PLAIN LANGUAGE SUMMARY

Routine diagnostic tests carry several limitations when applied after a suspected first seizure. SIRS is a noninvasive MRI method to enable time-sensitive diagnosis of image correlates of epileptic activity with increased sensitivity compared to routine EEG.

Spike ripples localize the epileptogenic zone best: an international intracranial study

We evaluated whether spike ripples, the combination of epileptiform spikes and ripples, provide a reliable and improved biomarker for the epileptogenic zone compared with other leading interictal biomarkers in a multicentre, international study. We first validated an automated spike ripple detector on intracranial EEG recordings. We then applied this detector to subjects from four centres who subsequently underwent surgical resection with known 1-year outcomes. We evaluated the spike ripple rate in subjects cured after resection [International League Against Epilepsy Class 1 outcome (ILAE 1)] and those with persistent seizures (ILAE 2-6) across sites and recording types. We also evaluated available interictal biomarkers: spike, spike-gamma, wideband high frequency oscillation (HFO, 80-500 Hz), ripple (80-250 Hz) and fast ripple (250-500 Hz) rates using previously validated automated detectors. The proportion of resected events was computed and compared across subject outcomes and biomarkers. Overall, 109 subjects were included. Most spike ripples were removed in subjects with ILAE 1 outcome (P < 0.001), and this was qualitatively observed across all sites and for depth and subdural electrodes (P < 0.001 and P < 0.001, respectively). Among ILAE 1 subjects, the mean spike ripple rate was higher in the resected volume (0.66/min) than in the non-removed tissue (0.08/min, P < 0.001). A higher proportion of spike ripples were removed in subjects with ILAE 1 outcomes compared with ILAE 2-6 outcomes (P = 0.06). Among ILAE 1 subjects, the proportion of spike ripples removed was higher than the proportion of spikes (P < 0.001), spike-gamma (P < 0.001), wideband HFOs (P < 0.001), ripples (P = 0.009) and fast ripples (P = 0.009) removed. At the individual level, more subjects with ILAE 1 outcomes had the majority of spike ripples removed (79%, 38/48) than spikes (69%, P = 0.12), spike-gamma (69%, P = 0.12), wideband HFOs (63%, P = 0.03), ripples (45%, P = 0.01) or fast ripples (36%, P < 0.001) removed. Thus, in this large, multicentre cohort, when surgical resection was successful, the majority of spike ripples were removed. Furthermore, automatically detected spike ripples localize the epileptogenic tissue better than spikes, spike-gamma, wideband HFOs, ripples and fast ripples.

Prognostic Value of Complete Resection of the High-Frequency Oscillation Area in Intracranial EEG: A Systematic Review and Meta-Analysis

BACKGROUND AND OBJECTIVES

High-frequency oscillations (HFOs; ripples 80-250 Hz; fast ripples [FRs] 250-500 Hz) recorded with intracranial electrodes generated excitement and debate about their potential to localize epileptogenic foci. We performed a systematic review and meta-analysis on the prognostic value of complete resection of the HFOs-area (crHFOs-area) for epilepsy surgical outcome in intracranial EEG (iEEG) accessing multiple subgroups.

METHODS

We searched PubMed, Embase, and Web of Science for original research from inception to October 27, 2022. We defined favorable surgical outcome (FSO) as Engel class I, International League Against Epilepsy class 1, or seizure-free status. The prognostic value of crHFOs-area for FSO was assessed by (1) the pooled FSO proportion after crHFOs-area; (2) FSO for crHFOs-area vs without crHFOs-area; and (3) the predictive performance. We defined high combined prognostic value as FSO proportion >80% + FSO crHFOs-area >without crHFOs-area + area under the curve (AUC) >0.75 and examined this for the clinical subgroups (study design, age, diagnostic type, HFOs-identification method, HFOs-rate thresholding, and iEEG state). Temporal lobe epilepsy (TLE) was compared with extra-TLE through dichotomous variable analysis. Individual patient analysis was performed for sex, affected hemisphere, MRI findings, surgery location, and pathology.

RESULTS

Of 1,387 studies screened, 31 studies (703 patients) met our eligibility criteria. Twenty-seven studies (602 patients) analyzed FRs and 20 studies (424 patients) ripples. Pooled FSO proportion after crHFOs-area was 81% (95% CI 76%-86%) for FRs and 82% (73%-89%) for ripples. Patients with crHFOs-area achieved more often FSO than those without crHFOs-area (FRs odds ratio [OR] 6.38, 4.03-10.09, p < 0.001; ripples 4.04, 2.32-7.04, p < 0.001). The pooled AUCs were 0.81 (0.77-0.84) for FRs and 0.76 (0.72-0.79) for ripples. Combined prognostic value was high in 10 subgroups: retrospective, children, long-term iEEG, threshold (FRs and ripples) and automated detection and interictal (FRs). FSO after complete resection of FRs-area (crFRs-area) was achieved less often in people with TLE than extra-TLE (OR 0.37, 0.15-0.89, p = 0.006). Individual patient analyses showed that crFRs-area was seen more in patients with FSO with than without MRI lesions (p = 0.02 after multiple correction).

DISCUSSION

Complete resection of the brain area with HFOs is associated with good postsurgical outcome. Its prognostic value holds, especially for FRs, for various subgroups. The use of HFOs for extra-TLE patients requires further evidence.

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.

Association between Removal of High-Frequency Oscillations and the Effect of Epilepsy Surgery: A Meta-Analysis

BACKGROUND

 High-frequency oscillations (HFOs) are spontaneous electroencephalographic (EEG) events that occur within the frequency range of 80 to 500 Hz and consist of at least four distinct oscillations that stand out from the background activity. They can be further classified into "ripples" (80-250 Hz) and "fast ripples" (FR; 250-500 Hz) based on different frequency bands. Studies have indicated that HFOs may serve as important markers for identifying epileptogenic regions and networks in patients with refractory epilepsy. Furthermore, a higher extent of removal of brain regions generating HFOs could potentially lead to improved prognosis. However, the clinical application criteria for HFOs remain controversial, and the results from different research groups exhibit inconsistencies. Given this controversy, the aim of this study was to conduct a meta-analysis to explore the utility of HFOs in predicting postoperative seizure outcomes by examining the prognosis of refractory epilepsy patients with varying ratios of HFO removal.

METHODS

 Prospective and retrospective studies that analyzed HFOs and postoperative seizure outcomes in epilepsy patients who underwent resective surgery were included in the meta-analysis. The patients in these studies were grouped based on the ratio of HFOs removed, resulting in four groups: completely removed FR (C-FR), completely removed ripples (C-Ripples), mostly removed FR (P-FR), and partial ripples removal (P-Ripples). The prognosis of patients within each group was compared to investigate the correlation between the ratio of HFO removal and patient prognosis.

RESULTS

 A total of nine studies were included in the meta-analysis. The prognosis of patients in the C-FR group was significantly better than that of patients with incomplete FR removal (odds ratio [OR] = 6.62; 95% confidence interval [CI]: 3.10-14.15; p < 0.00001). Similarly, patients in the C-Ripples group had a more favorable prognosis compared with those with incomplete ripples removal (OR = 4.45; 95% CI: 1.33-14.89; p = 0.02). Patients in the P-FR group had better prognosis than those with a majority of FR remaining untouched (OR = 6.23; 95% CI: 2.04-19.06; p = 0.001). In the P-Ripples group, the prognosis of patients with a majority of ripples removed was superior to that of patients with a majority of ripples remaining untouched (OR = 8.14; 95% CI: 2.62-25.33; p = 0.0003).

CONCLUSIONS

 There is a positive correlation between the greater removal of brain regions generating HFOs and more favorable postoperative seizure outcomes. However, further investigations, particularly through clinical trials, are necessary to justify the clinical application of HFOs in guiding epilepsy surgery.

Pathological and Physiological High-frequency Oscillations on Electroencephalography in Patients with Epilepsy

High-frequency oscillations (HFOs) encompass ripples (80 Hz-200 Hz) and fast ripples (200 Hz-600 Hz), serving as a promising biomarker for localizing the epileptogenic zone in epilepsy. Spontaneous fast ripples are always pathological, while ripples may be physiological or pathological. Distinguishing physiological from pathological ripples is important not only for designating epileptogenic brain regions, but also for investigations that study ripples in the context of memory encoding, consolidation, and recall in patients with epilepsy. Many studies have sought to identify distinguishing features between pathological and physiological ripples over the past two decades. Physiological and pathological ripples differ with respect to their spatial location, cellular mechanisms, morphology, and coupling with background electroencephalographic activity. Retrospective studies have demonstrated that differentiating between pathological and physiological ripples can improve surgical outcome prediction. In this review, we summarize the characteristics, differences, and applications of pathological and physiological HFOs and discuss strategies for their clinical translation.

Epileptic EEG patterns recognition through machine learning techniques and relevant time-frequency features

OBJECTIVES

The present study is designed to explore the process of epileptic patterns' automatic detection, specifically, epileptic spikes and high-frequency oscillations (HFOs), via a selection of machine learning (ML) techniques. The primary motivation for conducting such a research lies mainly in the need to investigate the long-term electroencephalography (EEG) recordings' visual examination process, often considered as a time-consuming and potentially error-prone procedure, requiring a great deal of mental focus and highly experimented neurologists. On attempting to resolve such a challenge, a number of state-of-the-art ML algorithms have been evaluated and compare in terms of performance, to pinpoint the most effective algorithm fit for accurately extracting epileptic EEG patterns.

CONTENT

Based on intracranial as well as simulated EEG data, the attained findings turn out to reveal that the randomforest (RF) method proved to be the most consistently effective approach, significantly outperforming the entirety of examined methods in terms of EEG recordings epileptic-pattern identification. Indeed, the RF classifier appeared to record an average balanced classification rate (BCR) of 92.38 % in regard to spikes recognition process, and 78.77 % in terms of HFOs detection.

SUMMARY

Compared to other approaches, our results provide valuable insights into the RF classifier's effectiveness as a powerful ML technique, fit for detecting EEG signals born epileptic bursts.

OUTLOOK

As a potential future work, we envisage to further validate and sustain our major reached findings through incorporating a larger EEG dataset. We also aim to explore the generative adversarial networks (GANs) application so as to generate synthetic EEG signals or combine signal generation techniques with deep learning approaches. Through this new vein of thought, we actually preconize to enhance and boost the automated detection methods' performance even more, thereby, noticeably enhancing the epileptic EEG pattern recognition area.

Generalisability of epileptiform patterns across time and patients

The complexity of localising the epileptogenic zone (EZ) contributes to surgical resection failures in achieving seizure freedom. The distinct patterns of epileptiform activity during interictal and ictal phases, varying across patients, often lead to suboptimal localisation using electroencephalography (EEG) features. We posed two key questions: whether neural signals reflecting epileptogenicity generalise from interictal to ictal time windows within each patient, and whether epileptiform patterns generalise across patients. Utilising an intracranial EEG dataset from 55 patients, we extracted a large battery of simple to complex features from stereo-EEG (SEEG) and electrocorticographic (ECoG) neural signals during interictal and ictal windows. Our features (n = 34) quantified many aspects of the signals including statistical moments, complexities, frequency-domain and cross-channel network attributes. Decision tree classifiers were then trained and tested on distinct time windows and patients to evaluate the generalisability of epileptogenic patterns across time and patients, respectively. Evidence strongly supported generalisability from interictal to ictal time windows across patients, particularly in signal power and high-frequency network-based features. Consistent patterns of epileptogenicity were observed across time windows within most patients, and signal features of epileptogenic regions generalised across patients, with higher generalisability in the ictal window. Signal complexity features were particularly contributory in cross-patient generalisation across patients. These findings offer insights into generalisable features of epileptic neural activity across time and patients, with implications for future automated approaches to supplement other EZ localisation methods.

Combining independent component analysis and source localization for improving spatial sampling of stereoelectroencephalography in epilepsy

Stereoelectroencephalography is a powerful intracerebral EEG recording method for the presurgical evaluation of epilepsy. It consists in implanting depth electrodes in the patient's brain to record electrical activity and map the epileptogenic zone, which should be resected to render the patient seizure-free. Stereoelectroencephalography has high spatial accuracy and signal-to-noise ratio but remains limited in the coverage of the explored brain regions. Thus, the implantation might provide a suboptimal sampling of epileptogenic regions. We investigate the potential of improving a suboptimal stereoelectroencephalography recording by performing source localization on stereoelectroencephalography signals. We propose combining independent component analysis, connectivity measures to identify components of interest, and distributed source modelling. This approach was tested on two patients with two implantations each, the first failing to characterize the epileptogenic zone and the second giving a better diagnosis. We demonstrate that ictal and interictal source localization performed on the first stereoelectroencephalography recordings matches the findings of the second stereo-EEG exploration. Our findings suggest that independent component analysis followed by source localization on the topographies of interest is a promising method for retrieving the epileptogenic zone in case of suboptimal implantation.

Magnetoencephalography for Epilepsy Presurgical Evaluation

PURPOSE OF THE REVIEW

Magnetoencephalography (MEG) is a functional neuroimaging technique that records neurophysiology data with millisecond temporal resolution and localizes it with subcentimeter accuracy. Its capability to provide high resolution in both of these domains makes it a powerful tool both in basic neuroscience as well as clinical applications. In neurology, it has proven useful in its ability to record and localize epileptiform activity. Epilepsy workup typically begins with scalp electroencephalography (EEG), but in many situations, EEG-based localization of the epileptogenic zone is inadequate. The complementary sensitivity of MEG can be crucial in such cases, and MEG has been adopted at many centers as an important resource in building a surgical hypothesis. In this paper, we review recent work evaluating the extent of MEG influence of presurgical evaluations, novel analyses of MEG data employed in surgical workup, and new MEG instrumentation that will likely affect the field of clinical MEG.

RECENT FINDINGS

MEG consistently contributes to presurgical evaluation and these contributions often change the plan for epilepsy surgery. Extensive work has been done to develop new analytic methods for localizing the source of epileptiform activity with MEG. Systems using optically pumped magnetometry (OPM) have been successfully deployed to record and localize epileptiform activity. MEG remains an important noninvasive tool for epilepsy presurgical evaluation. Continued improvements in analytic methodology will likely increase the diagnostic yield of the test. Novel instrumentation with OPM may contribute to this as well, and may increase accessibility of MEG by decreasing cost.

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.

High frequency oscillations in relation to interictal spikes in predicting postsurgical seizure freedom

We evaluate whether interictal spikes, epileptiform HFOs and their co-occurrence (Spike + HFO) were included in the resection area with respect to seizure outcome. We also characterise the relationship between high frequency oscillations (HFOs) and propagating spikes. We analysed intracranial EEG of 20 patients that underwent resective epilepsy surgery. The co-occurrence of ripples and fast ripples was considered an HFO event; the co-occurrence of an interictal spike and HFO was considered a Spike + HFO event. HFO distribution and spike onset were compared in cases of spike propagation. Accuracy in predicting seizure outcome was 85% for HFO, 60% for Spikes, and 79% for Spike + HFO. Sensitivity was 57% for HFO, 71% for Spikes and 67% for Spikes + HFO. Specificity was 100% for HFO, 54% for Spikes and 85% for Spikes + HFO. In 2/2 patients with spike propagation, the spike onset included the HFO area. Combining interictal spikes with HFO had comparable accuracy to HFO. In patients with propagating spikes, HFO rate was maximal at the onset of spike propagation.

Passive and active markers of cortical excitability in epilepsy

Electroencephalography (EEG) has been the primary diagnostic tool in clinical epilepsy for nearly a century. Its review is performed using qualitative clinical methods that have changed little over time. However, the intersection of higher resolution digital EEG and analytical tools developed in the past decade invites a re-exploration of relevant methodology. In addition to the established spatial and temporal markers of spikes and high-frequency oscillations, novel markers involving advanced postprocessing and active probing of the interictal EEG are gaining ground. This review provides an overview of the EEG-based passive and active markers of cortical excitability in epilepsy and of the techniques developed to facilitate their identification. Several different emerging tools are discussed in the context of specific EEG applications and the barriers we must overcome to translate these tools into clinical practice.

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.

High-rate leading spikes in propagating spike sequences predict seizure outcome in surgical patients with temporal lobe epilepsy

Inter-ictal spikes aid in the diagnosis of epilepsy and in planning surgery of medication-resistant epilepsy. However, the localizing information from spikes can be unreliable because spikes can propagate, and the burden of spikes, often assessed as a rate, does not always correlate with the seizure onset zone or seizure outcome. Recent work indicates identifying where spikes regularly emerge and spread could localize the seizure network. Thus, the current study sought to better understand where and how rates of single and coupled spikes, and especially brain regions with high-rate and leading spike of a propagating sequence, informs the extent of the seizure network. In 37 patients with medication-resistant temporal lobe seizures, who had surgery to treat their seizure disorder, an algorithm detected spikes in the pre-surgical depth inter-ictal EEG. A separate algorithm detected spike propagation sequences and identified the location of leading and downstream spikes in each sequence. We analysed the rate and power of single spikes on each electrode and coupled spikes between pairs of electrodes, and the proportion of sites with high-rate, leading spikes in relation to the seizure onset zone of patients seizure free (n = 19) and those with continuing seizures (n = 18). We found increased rates of single spikes in mesial temporal seizure onset zone (ANOVA, P < 0.001, η2 = 0.138), and increased rates of coupled spikes within, but not between, mesial-, lateral- and extra-temporal seizure onset zone of patients with continuing seizures (P < 0.001; η2 = 0.195, 0.113 and 0.102, respectively). In these same patients, there was a higher proportion of brain regions with high-rate leaders, and each sequence contained a greater number of spikes that propagated with a higher efficiency over a longer distance outside the seizure onset zone than patients seizure free (Wilcoxon, P = 0.0172). The proportion of high-rate leaders in and outside the seizure onset zone could predict seizure outcome with area under curve = 0.699, but not rates of single or coupled spikes (0.514 and 0.566). Rates of coupled spikes to a greater extent than single spikes localize the seizure onset zone and provide evidence for inter-ictal functional segregation, which could be an adaptation to avert seizures. Spike rates, however, have little value in predicting seizure outcome. High-rate spike sites leading propagation could represent sources of spikes that are important components of an efficient seizure network beyond the clinical seizure onset zone, and like the seizure onset zone these, too, need to be removed, disconnected or stimulated to increase the likelihood for seizure control.

Developmental atlas of phase-amplitude coupling between physiologic high-frequency oscillations and slow waves

We investigated the developmental changes in high-frequency oscillation (HFO) and Modulation Index (MI) - the coupling measure between HFO and slow-wave phase. We generated normative brain atlases, using subdural EEG signals from 8251 nonepileptic electrode sites in 114 patients (ages 1.0-41.5 years) who achieved seizure control following resective epilepsy surgery. We observed a higher MI in the occipital lobe across all ages, and occipital MI increased notably during early childhood. The cortical areas exhibiting MI co-growth were connected via the vertical occipital fasciculi and posterior callosal fibers. While occipital HFO rate showed no significant age-association, the temporal, frontal, and parietal lobes exhibited an age-inversed HFO rate. Assessment of 1006 seizure onset sites revealed that z-score normalized MI and HFO rate were higher at seizure onset versus nonepileptic electrode sites. We have publicly shared our intracranial EEG data to enable investigators to validate MI and HFO-centric presurgical evaluations to identify the epileptogenic zone.

HFA analysis using scalp electroencephalograms in two cases of Rasmussen's syndrome

OBJECTIVE

In recent years, wide-band EEGs have been used to assess brain activity, and their effectiveness in the pathological analysis of epilepsy has been demonstrated. This report describes two cases of Rasmussen's syndrome (RS) in which high-frequency scalp EEGs were retrospectively analyzed to assess the pathological condition of epilepsy in RS.

METHODS

The two RS cases were divided into three periods: incipient, stable, and frequent seizure periods. Using the EEG record of each period, interictal epileptiform discharges (IEDs) were visually extracted. Subsequently, a time-frequency analysis was performed to calculate the rate of high-frequency activities (HFAs) (IED-HFA rate). Finally, differences between the three periods were examined.

RESULTS

IED-HFA rates significantly increased in the frequent seizure period compared with the stable period in both cases(P < 0.05).

CONCLUSION

there was a significant increase in HFAs superimposed over IEDs during the frequent seizure period compared to the stable period. HFAs are thought to be associated with epileptogenicity. Similarly, HFAs could be a useful biomarker for the pathological condition of epilepsy in RS.

Fast ripples reflect increased excitability that primes epileptiform spikes

The neuronal circuit disturbances that drive inter-ictal and ictal epileptiform discharges remain elusive. Using a combination of extra-operative macro-electrode and micro-electrode inter-ictal recordings in six pre-surgical patients during non-rapid eye movement sleep, we found that, exclusively in the seizure onset zone, fast ripples (200-600 Hz), but not ripples (80-200 Hz), frequently occur <300 ms before an inter-ictal intra-cranial EEG spike with a probability exceeding chance (bootstrapping, P < 1e-5). Such fast ripple events are associated with higher spectral power (P < 1e-10) and correlated with more vigorous neuronal firing than solitary fast ripple (generalized linear mixed-effects model, P < 1e-9). During the intra-cranial EEG spike that follows a fast ripple, action potential firing is lower than during an intra-cranial EEG spike alone (generalized linear mixed-effects model, P < 0.05), reflecting an inhibitory restraint of intra-cranial EEG spike initiation. In contrast, ripples do not appear to prime epileptiform spikes. We next investigated the clinical significance of pre-spike fast ripple in a separate cohort of 23 patients implanted with stereo EEG electrodes, who underwent resections. In non-rapid eye movement sleep recordings, sites containing a high proportion of fast ripple preceding intra-cranial EEG spikes correlate with brain areas where seizures begin more than solitary fast ripple (P < 1e-5). Despite this correlation, removal of these sites does not guarantee seizure freedom. These results are consistent with the hypothesis that fast ripple preceding EEG spikes reflect an increase in local excitability that primes EEG spike discharges preferentially in the seizure onset zone and that epileptogenic brain regions are necessary, but not sufficient, for initiating inter-ictal epileptiform discharges.

Spontaneous fast-ultradian dynamics of polymorphic interictal events in drug-resistant focal epilepsy

OBJECTIVE

We studied the rate dynamics of interictal events occurring over fast-ultradian time scales, as commonly examined in clinics to guide surgical planning in epilepsy.

METHODS

Stereo-electroencephalography (SEEG) traces of 35 patients with good surgical outcome (Engel I) were analyzed. For this we developed a general data mining method aimed at clustering the plethora of transient waveform shapes including interictal epileptiform discharges (IEDs) and assessed the temporal fluctuations in the capability of mapping the epileptogenic zone (EZ) of each type of event.

RESULTS

We found that the fast-ultradian dynamics of the IED rate may effectively impair the precision of EZ identification, and appear to occur spontaneously, that is, not triggered by or exclusively associated with a particular cognitive task, wakefulness, sleep, seizure occurrence, post-ictal state, or antiepileptic drug withdrawal. Propagation of IEDs from the EZ to the propagation zone (PZ) could explain the observed fast-ultradian fluctuations in a reduced fraction of the analyzed patients, suggesting that other factors like the excitability of the epileptogenic tissue could play a more relevant role. A novel link was found between the fast-ultradian dynamics of the overall rate of polymorphic events and the rate of specific IEDs subtypes. We exploited this feature to estimate in each patient the 5 min interictal epoch for near-optimal EZ and resected-zone (RZ) localization. This approach produces at the population level a better EZ/RZ classification when compared to both (1) the whole time series available in each patient (p = .084 for EZ, p < .001 for RZ, Wilcoxon signed-rank test) and (2) 5 min epochs sampled randomly from the interictal recordings of each patient (p < .05 for EZ, p < .001 for RZ, 105 random samplings).

SIGNIFICANCE

Our results highlight the relevance of the fast-ultradian IED dynamics in mapping the EZ, and show how this dynamics can be estimated prospectively to inform surgical planning in epilepsy.

Seizure Onset Localization in Focal Epilepsy using intracranial-EEG data and the Schrodinger Operator's Spectrum

Epilepsy is a neurological disorder characterized by recurrent, unprovoked seizures that vary from short attention failure to convulsions. Despite its threats and limitations, existing medications target only specific types of seizures while up to 33% of epileptic conditions are drug-resistant. The best available treatment is surgical resection or neurostimulation and both require accurate localization of the Seizure Onset Zone. Its delineation is performed by analyzing neuronal activity by epileptologists, however, it is time-consuming and error-prone. Therefore, if the said zone could be located faster and more accurately, the seizure freedom of patients would be significantly enhanced. An effort within the field is aiming at developing computer-aided methods to assist medical experts and this starts with characterizing electrical neural activity. In the present paper, a new method for characterizing the epileptic intracranial EEG is proposed. The method is based on a semi-classical signal analysis (SCSA) method. Functional connectivity measures are used to compare patterns observed when feeding these measures with the raw time-series and when feeding them with SCSA features. The obtained results are undeniably promising for further investigation and improvement of the framework.Clinical relevance- The paper contributes to the design methods and algorithms to build reliable software solutions to assist medical experts in identifying Seizure Onset Zone in focal epilepsy.

The role of superficial and deep layers in the generation of high frequency oscillations and interictal epileptiform discharges in the human cortex

Describing intracortical laminar organization of interictal epileptiform discharges (IED) and high frequency oscillations (HFOs), also known as ripples. Defining the frequency limits of slow and fast ripples. We recorded potential gradients with laminar multielectrode arrays (LME) for current source density (CSD) and multi-unit activity (MUA) analysis of interictal epileptiform discharges IEDs and HFOs in the neocortex and mesial temporal lobe of focal epilepsy patients. IEDs were observed in 20/29, while ripples only in 9/29 patients. Ripples were all detected within the seizure onset zone (SOZ). Compared to hippocampal HFOs, neocortical ripples proved to be longer, lower in frequency and amplitude, and presented non-uniform cycles. A subset of ripples (≈ 50%) co-occurred with IEDs, while IEDs were shown to contain variable high-frequency activity, even below HFO detection threshold. The limit between slow and fast ripples was defined at 150 Hz, while IEDs' high frequency components form clusters separated at 185 Hz. CSD analysis of IEDs and ripples revealed an alternating sink-source pair in the supragranular cortical layers, although fast ripple CSD appeared lower and engaged a wider cortical domain than slow ripples MUA analysis suggested a possible role of infragranularly located neural populations in ripple and IED generation. Laminar distribution of peak frequencies derived from HFOs and IEDs, respectively, showed that supragranular layers were dominated by slower (< 150 Hz) components. Our findings suggest that cortical slow ripples are generated primarily in upper layers while fast ripples and associated MUA in deeper layers. The dissociation of macro- and microdomains suggests that microelectrode recordings may be more selective for SOZ-linked ripples. We found a complex interplay between neural activity in the neocortical laminae during ripple and IED formation. We observed a potential leading role of cortical neurons in deeper layers, suggesting a refined utilization of LMEs in SOZ localization.

Sevoflurane-induced high-frequency oscillations, effective connectivity and intraoperative classification of epileptic brain areas

OBJECTIVE

To determine how sevoflurane anesthesia modulates intraoperative epilepsy biomarkers on electrocorticography, including high-frequency oscillation (HFO) effective connectivity (EC), and to investigate their relation to epileptogenicity and anatomical white matter.

METHODS

We studied eight pediatric drug-resistant focal epilepsy patients who achieved seizure control after invasive monitoring and resective surgery. We visualized spatial distributions of the electrocorticography biomarkers at an oxygen baseline, three time-points while sevoflurane was increasing, and at a plateau of 2 minimum alveolar concentration (MAC) sevoflurane. HFO EC was combined with diffusion-weighted imaging, in dynamic tractography.

RESULTS

Intraoperative HFO EC diffusely increased as a function of sevoflurane concentration, although most in epileptogenic sites (defined as those included in the resection); their ability to classify epileptogenicity was optimized at sevoflurane 2 MAC. HFO EC could be visualized on major white matter tracts, as a function of sevoflurane level.

CONCLUSIONS

The results strengthened the hypothesis that sevoflurane-activated HFO biomarkers may help intraoperatively localize the epileptogenic zone.

SIGNIFICANCE

Our results help characterize how HFOs at non-epileptogenic and epileptogenic networks respond to sevoflurane. It may be warranted to establish a normative HFO atlas incorporating the modifying effects of sevoflurane and major white matter pathways, as critical reference in epilepsy presurgical evaluation.

Quantitative approaches to guide epilepsy surgery from intracranial EEG

Over the past 10 years, the drive to improve outcomes from epilepsy surgery has stimulated widespread interest in methods to quantitatively guide epilepsy surgery from intracranial EEG (iEEG). Many patients fail to achieve seizure freedom, in part due to the challenges in subjective iEEG interpretation. To address this clinical need, quantitative iEEG analytics have been developed using a variety of approaches, spanning studies of seizures, interictal periods, and their transitions, and encompass a range of techniques including electrographic signal analysis, dynamical systems modeling, machine learning and graph theory. Unfortunately, many methods fail to generalize to new data and are sensitive to differences in pathology and electrode placement. Here, we critically review selected literature on computational methods of identifying the epileptogenic zone from iEEG. We highlight shared methodological challenges common to many studies in this field and propose ways that they can be addressed. One fundamental common pitfall is a lack of open-source, high-quality data, which we specifically address by sharing a centralized high-quality, well-annotated, multicentre dataset consisting of >100 patients to support larger and more rigorous studies. Ultimately, we provide a road map to help these tools reach clinical trials and hope to improve the lives of future patients.

Neurophysiology, Neuropsychology, and Epilepsy, in 2022: Hills We Have Climbed and Hills Ahead. Neurophysiology in epilepsy

Since the discovery of the human electroencephalogram (EEG), neurophysiology techniques have become indispensable tools in our armamentarium to localize epileptic seizures. New signal analysis techniques and the prospects of artificial intelligence and big data will offer unprecedented opportunities to further advance the field in the near future, ultimately resulting in improved quality of life for many patients with drug-resistant epilepsy. This article summarizes selected presentations from Day 1 of the two-day symposium "Neurophysiology, Neuropsychology, Epilepsy, 2022: Hills We Have Climbed and the Hills Ahead". Day 1 was dedicated to highlighting and honoring the work of Dr. Jean Gotman, a pioneer in EEG, intracranial EEG, simultaneous EEG/ functional magnetic resonance imaging, and signal analysis of epilepsy. The program focused on two main research directions of Dr. Gotman, and was dedicated to "High-frequency oscillations, a new biomarker of epilepsy" and "Probing the epileptic focus from inside and outside". All talks were presented by colleagues and former trainees of Dr. Gotman. The extended summaries provide an overview of historical and current work in the neurophysiology of epilepsy with emphasis on novel EEG biomarkers of epilepsy and source imaging and concluded with an outlook on the future of epilepsy research, and what is needed to bring the field to the next level.

Fast ripples reflect increased excitability that primes epileptiform spikes

The neuronal circuit disturbances that drive interictal and ictal epileptiform discharges remains elusive. Using a combination of extraoperative macro- and micro-electrode interictal recordings in six presurgical patients during non-rapid eye movement (REM) sleep we found that, exclusively in the seizure onset zone, fast ripples (FR; 200-600Hz), but not ripples (80-200 Hz), frequently occur <300 msec before an interictal intracranial EEG (iEEG) spike with a probability exceeding chance (bootstrapping, p<1e-5). Such FR events are associated with higher spectral power (p<1e-10) and correlated with more vigorous neuronal firing than solitary FR (generalized linear mixed-effects model, GLMM, p<1e-3) irrespective of FR power. During the iEEG spike that follows a FR, action potential firing is lower than during a iEEG spike alone (GLMM, p<1e-10), reflecting an inhibitory restraint of iEEG spike initiation. In contrast, ripples do not appear to prime epileptiform spikes. We next investigated the clinical significance of pre-spike FR in a separate cohort of 23 patients implanted with stereo EEG electrodes who underwent resections. In non-REM sleep recordings, sites containing a high proportion of FR preceding iEEG spikes correlate with brain areas where seizures begin more than solitary FR (p<1e-5). Despite this correlation, removal of these sites does not guarantee seizure freedom. These results are consistent with the hypothesis that FR preceding EEG spikes reflect an increase in local excitability that primes EEG spike discharges preferentially in the seizure onset zone and that epileptogenic brain regions are necessary, but not sufficient, for initiating interictal epileptiform discharges.

Detection of pathological high-frequency oscillations in refractory epilepsy patients undergoing simultaneous stereo-electroencephalography and magnetoencephalography

BACKGROUND

Stereo-electroencephalography (SEEG) and magnetoencephalography (MEG) have generally been used independently as part of the pre-surgical evaluation of drug-resistant epilepsy (DRE) patients. However, the possibility of simultaneously employing these recording techniques to determine whether MEG has the potential of offering the same information as SEEG less invasively, or whether it could offer a greater spatial indication of the epileptogenic zone (EZ) to aid surgical planning, has not been previously evaluated.

METHODS

Data from 24 paediatric and adult DRE patients, undergoing simultaneous SEEG and MEG as part of their pre-surgical evaluation, was analysed employing manual and automated high-frequency oscillations (HFOs) detection, and spectral and source localisation analyses.

RESULTS

Twelve patients (50%) were included in the analysis (4 males; mean age=25.08 years) and showed interictal SEEG and MEG HFOs. HFOs detection was concordant between the two recording modalities, but SEEG displayed higher ability of differentiating between deep and superficial epileptogenic sources. Automated HFO detector in MEG recordings was validated against the manual MEG detection method. Spectral analysis revealed that SEEG and MEG detect distinct epileptic events. The EZ was well correlated with the simultaneously recorded data in 50% patients, while 25% patients displayed poor correlation or discordance.

CONCLUSION

MEG recordings can detect HFOs, and simultaneous use of SEEG and MEG HFO identification facilitates EZ localisation during the presurgical planning stage for DRE patients. Further studies are necessary to validate these findings and support the translation of automated HFO detectors into routine clinical practice.

Multilayer brain networks can identify the epileptogenic zone and seizure dynamics

Seizure generation, propagation, and termination occur through spatiotemporal brain networks. In this paper, we demonstrate the significance of large-scale brain interactions in high-frequency (80-200Hz) for the identification of the epileptogenic zone (EZ) and seizure evolution. To incorporate the continuity of neural dynamics, here we have modeled brain connectivity constructed from stereoelectroencephalography (SEEG) data during seizures using multilayer networks. After introducing a new measure of brain connectivity for temporal networks, named multilayer eigenvector centrality (mlEVC), we applied a consensus hierarchical clustering on the developed model to identify the EZ as a cluster of nodes with distinctive brain connectivity in the ictal period. Our algorithm could successfully predict electrodes inside the resected volume as EZ for 88% of participants, who all were seizure-free for at least 12 months after surgery. Our findings illustrated significant and unique desynchronization between EZ and the rest of the brain in the early to mid-seizure. We showed that aging and the duration of epilepsy intensify this desynchronization, which can be the outcome of abnormal neuroplasticity. Additionally, we illustrated that seizures evolve with various network topologies, confirming the existence of different epileptogenic networks in each patient. Our findings suggest not only the importance of early intervention in epilepsy but possible factors that correlate with disease severity. Moreover, by analyzing the propagation patterns of different seizures, we demonstrate the necessity of collecting sufficient data for identifying epileptogenic networks.

A Subpopulation of Spikes Predicts Successful Epilepsy Surgery Outcome

OBJECTIVE

Epileptic spikes are the traditional interictal electroencephalographic (EEG) biomarker for epilepsy. Given their low specificity for identifying the epileptogenic zone (EZ), they are given only moderate attention in presurgical evaluation. This study aims to demonstrate that it is possible to identify specific spike features in intracranial EEG that optimally define the EZ and predict surgical outcome.

METHODS

We analyzed spike features on stereo-EEG segments from 83 operated patients from 2 epilepsy centers (37 Engel IA) in wakefulness, non-rapid eye movement sleep, and rapid eye movement sleep. After automated spike detection, we investigated 135 spike features based on rate, morphology, propagation, and energy to determine the best feature or feature combination to discriminate the EZ in seizure-free and non-seizure-free patients by applying 4-fold cross-validation.

RESULTS

The rate of spikes with preceding gamma activity in wakefulness performed better for surgical outcome classification (4-fold area under receiver operating characteristics curve [AUC] = 0.755 ± 0.07) than the seizure onset zone, the current gold standard (AUC = 0.563 ± 0.05, p = 0.015) and the ripple rate, an emerging seizure-independent biomarker (AUC = 0.537 ± 0.07, p = 0.006). Channels with a spike-gamma rate exceeding 1.9/min had an 80% probability of being in the EZ. Combining features did not improve the results.

INTERPRETATION

Resection of brain regions with high spike-gamma rates in wakefulness is associated with a high probability of achieving seizure freedom. This rate could be applied to determine the minimal number of spiking channels requiring resection. In addition to quantitative analysis, this feature is easily accessible to visual analysis, which could aid clinicians during presurgical evaluation. ANN NEUROL 2023;93:522-535.

Interictal HFO and FDG-PET correlation predicts surgical outcome following SEEG

OBJECTIVE

This study aimed to investigate the quantitative relationship between interictal ¹⁸ F-fluorodeoxyglucose-positron emission tomography (FDG-PET) and interictal high-frequency oscillations (HFOs) from stereo-electroencephalography (SEEG) recordings in patients with refractory epilepsy.

METHODS

We retrospectively included 32 patients. FDG-PET data were quantified through statistical parametric mapping (SPM) t test modeling with normal controls. Interictal SEEG segments with four, 10-min segments were selected randomly. HFO detection and classification procedures were automatically performed. Channel-based HFOs separating ripple (80-250 Hz) and fast ripple (FR; 250-500 Hz) counts were correlated with the surrounding metabolism T score at the individual and group level, respectively. The association was further validated across anatomic seizure origins and sleep vs wake states. We built a joint feature FR × T reflecting the FR and hypometabolism concordance to predict surgical outcomes in 28 patients who underwent surgery.

RESULTS

We found a negative correlation between interictal FDG-PET and HFOs through the linear mixed-effects model (R²  = .346 and .457 for ripples and FRs, respectively, p < .001); these correlations were generalizable to different epileptogenic-zone lobar localizations and vigilance states. The FR × T inside the resection volume could be used as a predictor for surgical outcomes with an area under the curve of 0.81.

SIGNIFICANCE

The degree of hypometabolism is associated with HFO generation rate, especially for FRs. This relationship would be meaningful for selection of SEEG candidates and for optimizing SEEG scheme planning. The concordance between FRs and hypometabolism inside the resection volume could provide prognostic information regarding surgical outcome.

Amortized Bayesian inference on generative dynamical network models of epilepsy using deep neural density estimators

Whole-brain modeling of epilepsy combines personalized anatomical data with dynamical models of abnormal activities to generate spatio-temporal seizure patterns as observed in brain imaging data. Such a parametric simulator is equipped with a stochastic generative process, which itself provides the basis for inference and prediction of the local and global brain dynamics affected by disorders. However, the calculation of likelihood function at whole-brain scale is often intractable. Thus, likelihood-free algorithms are required to efficiently estimate the parameters pertaining to the hypothetical areas, ideally including the uncertainty. In this study, we introduce the simulation-based inference for the virtual epileptic patient model (SBI-VEP), enabling us to amortize the approximate posterior of the generative process from a low-dimensional representation of whole-brain epileptic patterns. The state-of-the-art deep learning algorithms for conditional density estimation are used to readily retrieve the statistical relationships between parameters and observations through a sequence of invertible transformations. We show that the SBI-VEP is able to efficiently estimate the posterior distribution of parameters linked to the extent of the epileptogenic and propagation zones from sparse intracranial electroencephalography recordings. The presented Bayesian methodology can deal with non-linear latent dynamics and parameter degeneracy, paving the way for fast and reliable inference on brain disorders from neuroimaging modalities.

Relative entropy is an easy-to-use invasive electroencephalographic biomarker of the epileptogenic zone

OBJECTIVE

High-frequency oscillations are considered among the most promising interictal biomarkers of the epileptogenic zone in patients suffering from pharmacoresistant focal epilepsy. However, there is no clear definition of pathological high-frequency oscillations, and the existing detectors vary in methodology, performance, and computational costs. This study proposes relative entropy as an easy-to-use novel interictal biomarker of the epileptic tissue.

METHODS

We evaluated relative entropy and high-frequency oscillation biomarkers on intracranial electroencephalographic data from 39 patients with seizure-free postoperative outcome (Engel Ia) from three institutions. We tested their capability to localize the epileptogenic zone, defined as resected contacts located in the seizure onset zone. The performance was compared using areas under the receiver operating curves (AUROCs) and precision-recall curves. Then we tested whether a universal threshold can be used to delineate the epileptogenic zone across patients from different institutions.

RESULTS

Relative entropy in the ripple band (80-250 Hz) achieved an average AUROC of .85. The normalized high-frequency oscillation rate in the ripple band showed an identical AUROC of .85. In contrast to high-frequency oscillations, relative entropy did not require any patient-level normalization and was easy and fast to calculate due to its clear and straightforward definition. One threshold could be set across different patients and institutions, because relative entropy is independent of signal amplitude and sampling frequency.

SIGNIFICANCE

Although both relative entropy and high-frequency oscillations have a similar performance, relative entropy has significant advantages such as straightforward definition, computational speed, and universal interpatient threshold, making it an easy-to-use promising biomarker of the epileptogenic zone.

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.

Clinical utility of intraoperative electrocorticography for epilepsy surgery: A systematic review and meta-analysis

Despite the widespread use of intraoperative electrocorticography (iECoG) during resective epilepsy surgery, there are conflicting data on its overall efficacy and inability to predict benefit per pathology. Given the heterogeneity of iECoG use in resective epilepsy surgery, it is important to assess the utility of interictal-based iECoG. This individual patient data (IPD) meta-analysis seeks to identify the benefit of iECoG during resective epilepsy surgery in achieving seizure freedom for various pathologies. Embase, Scopus, and PubMed were searched from inception to January 31, 2021 using the following terms: "ecog", "electrocorticography", and "epilepsy". Articles were included if they reported seizure freedom at ≥12-month follow-up in cohorts with and without iECoG for epilepsy surgery. Non-English articles, noncomparative iECoG cohorts, and studies with <10% iECoG use were excluded. This meta-analysis followed the PRISMA 2020 guidelines. The primary outcome was seizure freedom at last follow-up and time to seizure recurrence, if applicable. Forest plots with random effects modeling assessed the relationship between iECoG use and seizure freedom. Cox regression of IPD was performed to identify predictors of longer duration of seizure freedom. Kaplan-Meier curves with log-rank test were created to visualize differences in time to seizure recurrence. Of 7504 articles identified, 18 were included for study-level analysis. iECoG was not associated with higher seizure freedom at the study level (relative risk = 1.09, 95% confidence interval [CI] = 0.96-1.23, p = .19, I²  = 64%), but on IPD (n = 7 studies, 231 patients) iECoG use was independently associated with more favorable seizure outcomes (hazard ratio = 0.47, 95% CI = .23-.95, p = .037). In Kaplan-Meier analysis of specific pathologies, iECoG use was significantly associated with longer seizure freedom only for focal cortical dysplasia (FCD; p < .001) etiology. Number needed to treat for iECoG was 8.8, and for iECoG in FCD it was 4.7. We show iECoG seizure freedom is not achieved uniformly across centers. iECoG is particularly beneficial for FCD etiology in improving seizure freedom.

Intrinsic and secondary epileptogenicity in focal cortical dysplasia type II

OBJECTIVE

Favorable seizure outcome is reported following resection of bottom-of-sulcus dysplasia (BOSD). We assessed the distribution of epileptogenicity and dysplasia in and around BOSD to better understand this clinical outcome and the optimal surgical approach.

METHODS

We studied 27 children and adolescents with magnetic resonance imaging (MRI)-positive BOSD who underwent epilepsy surgery; 85% became seizure-free postresection (median = 5.0 years follow-up). All patients had resection of the dysplastic sulcus, and 11 had additional resection of the gyral crown (GC) or adjacent gyri (AG). Markers of epileptogenicity were relative cortical hypometabolism on preoperative 18 F-fluorodeoxyglucose (FDG) positron emission tomography (PET), and spiking, ripples, fast ripples, spike-high-frequency oscillation cross-rate, and phase amplitude coupling (PAC) on preresection and postresection electrocorticography (ECoG), all analyzed at the bottom-of-sulcus (BOS), top-of-sulcus (TOS), GC, and AG. Markers of dysplasia were increased cortical thickness on preoperative MRI, and dysmorphic neuron density and variant allele frequency of somatic MTOR mutations in resected tissue, analyzed at similar locations.

RESULTS

Relative cortical metabolism was significantly reduced and ECoG markers were significantly increased at the BOS compared to other regions. Apart from spiking and PAC, which were greater at the TOS compared to the GC, there were no significant differences in PET and other ECoG markers between the TOS, GC, and AG, suggesting a cutoff of epileptogenicity at the TOS rather than a tapering gradient on the cortical surface. MRI and tissue markers of dysplasia were all maximal in the BOS, reduced in the TOS, and mostly absent in the GC. Spiking and PAC reduced significantly over the GC after resection of the dysplastic sulcus.

SIGNIFICANCE

These findings support the concept that dysplasia and intrinsic epileptogenicity are mostly limited to the dysplastic sulcus in BOSD and support resection or ablation confined to the MRI-visible lesion as a first-line surgical approach. 18 F-FDG PET and ECoG abnormalities in surrounding cortex seem to be secondary phenomena.

Shadows of very high-frequency oscillations can be detected in lower frequency bands of routine stereoelectroencephalography

Very high-frequency oscillations (VHFOs, > 500 Hz) are more specific in localizing the epileptogenic zone (EZ) than high-frequency oscillations (HFOs, < 500 Hz). Unfortunately, VHFOs are not visible in standard clinical stereo-EEG (SEEG) recordings with sampling rates of 1 kHz or lower. Here we show that "shadows" of VHFOs can be found in frequencies below 500 Hz and can help us to identify SEEG channels with a higher probability of increased VHFO rates. Subsequent analysis of Logistic regression models on 141 SEEG channels from thirteen patients shows that VHFO "shadows" provide additional information to gold standard HFO analysis and can potentially help in precise EZ delineation in standard clinical recordings.

EEG biomarker candidates for the identification of epilepsy

Electroencephalography (EEG) is one of the main pillars used for the diagnosis and study of epilepsy, readily employed after a possible first seizure has occurred. The most established biomarker of epilepsy, in case seizures are not recorded, are interictal epileptiform discharges (IEDs). In clinical practice, however, IEDs are not always present and the EEG may appear completely normal despite an underlying epileptic disorder, often leading to difficulties in the diagnosis of the disease. Thus, finding other biomarkers that reliably predict whether an individual suffers from epilepsy even in the absence of evident epileptic activity would be extremely helpful, since they could allow shortening the period of diagnostic uncertainty and consequently decreasing the risk of seizure. To date only a few EEG features other than IEDs seem to be promising candidates able to distinguish between epilepsy, i.e. > 60 % risk of recurrent seizures, or other (pathological) conditions. The aim of this narrative review is to provide an overview of the EEG-based biomarker candidates for epilepsy and the techniques employed for their identification.

Characterizing physiological high-frequency oscillations using deep learning

Objective.Intracranially-recorded interictal high-frequency oscillations (HFOs) have been proposed as a promising spatial biomarker of the epileptogenic zone. However, HFOs can also be recorded in the healthy brain regions, which complicates the interpretation of HFOs. The present study aimed to characterize salient features of physiological HFOs using deep learning (DL).Approach.We studied children with neocortical epilepsy who underwent intracranial strip/grid evaluation. Time-series EEG data were transformed into DL training inputs. The eloquent cortex (EC) was defined by functional cortical mapping and used as a DL label. Morphological characteristics of HFOs obtained from EC (ecHFOs) were distilled and interpreted through a novel weakly supervised DL model.Main results.A total of 63 379 interictal intracranially-recorded HFOs from 18 children were analyzed. The ecHFOs had lower amplitude throughout the 80-500 Hz frequency band around the HFO onset and also had a lower signal amplitude in the low frequency band throughout a one-second time window than non-ecHFOs, resembling a bell-shaped template in the time-frequency map. A minority of ecHFOs were HFOs with spikes (22.9%). Such morphological characteristics were confirmed to influence DL model prediction via perturbation analyses. Using the resection ratio (removed HFOs/detected HFOs) of non-ecHFOs, the prediction of postoperative seizure outcomes improved compared to using uncorrected HFOs (area under the ROC curve of 0.82, increased from 0.76).Significance.We characterized salient features of physiological HFOs using a DL algorithm. Our results suggested that this DL-based HFO classification, once trained, might help separate physiological from pathological HFOs, and efficiently guide surgical resection using HFOs.

Local neuronal excitation and global inhibition during epileptic fast ripples in humans

Understanding the neuronal basis of epileptic activity is a major challenge in neurology. Cellular integration into larger scale networks is all the more challenging. In the local field potential, interictal epileptic discharges can be associated with fast ripples (200-600 Hz), which are a promising marker of the epileptogenic zone. Yet, how neuronal populations in the epileptogenic zone and in healthy tissue are affected by fast ripples remain unclear. Here, we used a novel 'hybrid' macro-micro depth electrode in nine drug-resistant epileptic patients, combining classic depth recording of local field potentials (macro-contacts) and two or three tetrodes (four micro-wires bundled together) enabling up to 15 neurons in local circuits to be simultaneously recorded. We characterized neuronal responses (190 single units) with the timing of fast ripples (2233 fast ripples) on the same hybrid and other electrodes that target other brain regions. Micro-wire recordings reveal signals that are not visible on macro-contacts. While fast ripples detected on the closest macro-contact to the tetrodes were always associated with fast ripples on the tetrodes, 82% of fast ripples detected on tetrodes were associated with detectable fast ripples on the nearest macro-contact. Moreover, neuronal recordings were taken in and outside the epileptogenic zone of implanted epileptic subjects and they revealed an interlay of excitation and inhibition across anatomical scales. While fast ripples were associated with increased neuronal activity in very local circuits only, they were followed by inhibition in large-scale networks (beyond the epileptogenic zone, even in healthy cortex). Neuronal responses to fast ripples were homogeneous in local networks but differed across brain areas. Similarly, post-fast ripple inhibition varied across recording locations and subjects and was shorter than typical inter-fast ripple intervals, suggesting that this inhibition is a fundamental refractory process for the networks. These findings demonstrate that fast ripples engage local and global networks, including healthy tissue, and point to network features that pave the way for new diagnostic and therapeutic strategies. They also reveal how even localized pathological brain dynamics can affect a broad range of cognitive functions.

Metabolism and Intracranial Epileptogenicity in Temporal Lobe Long-Term Epilepsy-Associated Tumor

Brain tumors are common in epilepsy surgery and frequently occur in the temporal lobe, but the optimal surgical strategies to remove the tumor and epileptogenic zone remain controversial. We aim at illustrating the positron emission tomography (PET) metabolism and the stereoelectroencephalography (SEEG) epileptogenicity of temporal lobe long-term epilepsy-associated tumors (LEAT). In this study, 70 patients and 25 healthy controls were included. Our analysis leveraged group-level analysis to reveal the whole-brain metabolic pattern of temporal lobe LEATs. The SEEG-based epileptogenicity mapping was performed to verify the PET findings in the epileptic network. Compared to controls, patients with a temporal lobe LEAT showed a more widespread epileptic network based on 18FDG-PET in patients with a mesial temporal lobe LEAT than in those with a lateral temporal lobe LEAT. The significant brain clusters mainly involved the paracentral lobule (ANOVA F = 9.731, p < 0.001), caudate nucleus (ANOVA F = 20.749, p < 0.001), putamen (Kruskal−Wallis H = 19.258, p < 0.001), and thalamus (ANOVA F = 4.754, p = 0.011). Subgroup analysis and SEEG-based epileptogenicity mapping are similar to the metabolic pattern. Our findings demonstrate the metabolic and electrophysiological organization of the temporal lobe LEAT epileptic network, which may assist in a patient-specific surgical strategy.

Focal non-invasive deep-brain stimulation with temporal interference for the suppression of epileptic biomarkers

Introduction

Neurostimulation applied from deep brain stimulation (DBS) electrodes is an effective therapeutic intervention in patients suffering from intractable drug-resistant epilepsy when resective surgery is contraindicated or failed. Inhibitory DBS to suppress seizures and associated epileptogenic biomarkers could be performed with high-frequency stimulation (HFS), typically between 100 and 165 Hz, to various deep-seated targets, such as the Mesio-temporal lobe (MTL), which leads to changes in brain rhythms, specifically in the hippocampus. The most prominent alterations concern high-frequency oscillations (HFOs), namely an increase in ripples, a reduction in pathological Fast Ripples (FRs), and a decrease in pathological interictal epileptiform discharges (IEDs).

Materials and methods

In the current study, we use Temporal Interference (TI) stimulation to provide a non-invasive DBS (130 Hz) of the MTL, specifically the hippocampus, in both mouse models of epilepsy, and scale the method using human cadavers to demonstrate the potential efficacy in human patients. Simulations for both mice and human heads were performed to calculate the best coordinates to reach the hippocampus.

Results

This non-invasive DBS increases physiological ripples, and decreases the number of FRs and IEDs in a mouse model of epilepsy. Similarly, we show the inability of 130 Hz transcranial current stimulation (TCS) to achieve similar results. We therefore further demonstrate the translatability to human subjects via measurements of the TI stimulation vs. TCS in human cadavers. Results show a better penetration of TI fields into the human hippocampus as compared with TCS.

Significance

These results constitute the first proof of the feasibility and efficiency of TI to stimulate at depth an area without impacting the surrounding tissue. The data tend to show the sufficiently focal character of the induced effects and suggest promising therapeutic applications in epilepsy.

Spikes and High Frequency Oscillations in Lateral Neocortical Temporal Lobe Epilepsy: Can They Predict the Success Chance of Hippocampus-Sparing Resections?

Purpose

We investigated the distribution of spikes and HFOs recorded during intraoperative electrocorticography (ioECoG) and tried to elaborate a predictive model for postsurgical outcomes of patients with lateral neocortical temporal lobe epilepsy (TLE) whose mesiotemporal structures are left in situ.

Methods

We selected patients with temporal lateral neocortical epilepsy focus who underwent ioECoG-tailored resections without amygdalo–hippocampectomies. We visually marked spikes, ripples (80–250 Hz), and fast ripples (FRs; 250–500 Hz) on neocortical and mesiotemporal channels before and after resections. We looked for differences in event rates and resection ratios between good (Engel 1A) and poor outcome groups and performed logistic regression analysis to identify outcome predictors.

Results

Fourteen out of 24 included patients had a good outcome. The poor-outcome patients showed higher rates of ripples on neocortical channels distant from the resection in pre- and post-ioECoG than people with good outcomes (p_pre = 0.04, p_post = 0.05). Post-ioECoG FRs were found only in poor-outcome patients (N = 3). A prediction model based on regression analysis showed low rates of mesiotemporal post-ioECoG ripples (OR_mesio = 0.13, p_mesio = 0.04) and older age at epilepsy onset (OR = 1.76, p = 0.04) to be predictors of good seizure outcome.

Conclusion

HFOs in ioECoG may help to inform the neurosurgeon of the hippocampus-sparing resection success chance in patients with lateral neocortical TLE.