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

Simulated resections and responsive neurostimulator placement can optimize postoperative seizure outcomes when guided by fast ripple networks

In medication-resistant epilepsy, the goal of epilepsy surgery is to make a patient seizure free with a resection/ablation that is as small as possible to minimize morbidity. The standard of care in planning the margins of epilepsy surgery involves electroclinical delineation of the seizure-onset zone and incorporation of neuroimaging findings from MRI, PET, single-photon emission CT and magnetoencephalography modalities. Resecting cortical tissue generating high-frequency oscillations has been investigated as a more efficacious alternative to targeting the seizure-onset zone. In this study, we used a support vector machine (SVM), with four distinct fast ripple (FR: 350-600 Hz on oscillations, 200-600 Hz on spikes) metrics as factors. These metrics included the FR resection ratio, a spatial FR network measure and two temporal FR network measures. The SVM was trained by the value of these four factors with respect to the actual resection boundaries and actual seizure-free labels of 18 patients with medically refractory focal epilepsy. Leave-one-out cross-validation of the trained SVM in this training set had an accuracy of 0.78. We next used a simulated iterative virtual resection targeting the FR sites that were of highest rate and showed most temporal autonomy. The trained SVM utilized the four virtual FR metrics to predict virtual seizure freedom. In all but one of the nine patients who were seizure free after surgery, we found that the virtual resections sufficient for virtual seizure freedom were larger in volume (P < 0.05). In nine patients who were not seizure free, a larger virtual resection made five virtually seizure free. We also examined 10 medically refractory focal epilepsy patients implanted with the responsive neurostimulator system and virtually targeted the responsive neurostimulator system stimulation contacts proximal to sites generating FR at highest rates to determine if the simulated value of the stimulated seizure-onset zone and stimulated FR metrics would trend towards those patients with a better seizure outcome. Our results suggest the following: (i) FR measures can accurately predict whether a resection, defined by the standard of care, will result in seizure freedom; (ii) utilizing FR alone for planning an efficacious surgery can be associated with larger resections; (iii) when FR metrics predict the standard-of-care resection will fail, amending the boundaries of the planned resection with certain FR-generating sites may improve outcome and (iv) more work is required to determine whether targeting responsive neurostimulator system stimulation contact proximal to FR generating sites will improve seizure outcome.

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.

Information dynamics of in silico EEG Brain Waves: Insights into oscillations and functions

The relation between electroencephalography (EEG) rhythms, brain functions, and behavioral correlates is well-established. Some physiological mechanisms underlying rhythm generation are understood, enabling the replication of brain rhythms in silico. This offers a pathway to explore connections between neural oscillations and specific neuronal circuits, potentially yielding fundamental insights into the functional properties of brain waves. Information theory frameworks, such as Integrated Information Decomposition (Φ-ID), relate dynamical regimes with informational properties, providing deeper insights into neuronal dynamic functions. Here, we investigate wave emergence in an excitatory/inhibitory (E/I) balanced network of integrate and fire neurons with short-term synaptic plasticity. This model produces a diverse range of EEG-like rhythms, from low δ waves to high-frequency oscillations. Through Φ-ID, we analyze the network's information dynamics and its relation with different emergent rhythms, elucidating the system's suitability for functions such as robust information transfer, storage, and parallel operation. Furthermore, our study helps to identify regimes that may resemble pathological states due to poor informational properties and high randomness. We found, e.g., that in silico β and δ waves are associated with maximum information transfer in inhibitory and excitatory neuron populations, respectively, and that the coexistence of excitatory θ, α, and β waves is associated to information storage. Additionally, we observed that high-frequency oscillations can exhibit either high or poor informational properties, potentially shedding light on ongoing discussions regarding physiological versus pathological high-frequency oscillations. In summary, our study demonstrates that dynamical regimes with similar oscillations may exhibit vastly different information dynamics. Characterizing information dynamics within these regimes serves as a potent tool for gaining insights into the functions of complex neuronal networks. Finally, our findings suggest that the use of information dynamics in both model and experimental data analysis, could help discriminate between oscillations associated with cognitive functions and those linked to neuronal disorders.

Investigating current clinical opinions in stereoelectroencephalography-informed epilepsy surgery

OBJECTIVE

Stereoelectroencephalography (SEEG) is increasingly utilized worldwide in epilepsy surgery planning. International guidelines for SEEG terminology and interpretation are yet to be proposed. There are worldwide differences in SEEG definitions, application of features in epilepsy surgery planning, and interpretation of surgical outcomes. This hinders the clinical interpretation of SEEG findings and collaborative research. We aimed to assess the global perspectives on SEEG terminology, differences in the application of presurgical features, and variability in the interpretation of surgery outcome scores, and analyze how clinical expert demographics influenced these opinions.

METHODS

We assessed the practices and opinions of epileptologists with specialized training in SEEG using a survey. Data were qualitatively analyzed, and subgroups were examined based on geographical regions and years of experience. Primary outcomes included opinions on SEEG terminology, features used for epilepsy surgery, and interpretation of outcome scores. Additionally, we conducted a multilevel regression and poststratification analysis to characterize the nonresponders.

RESULTS

A total of 321 expert responses from 39 countries were analyzed. We observed substantial differences in terminology, practices, and use of presurgical features across geographical regions and SEEG expertise levels. The majority of experts (220, 68.5%) favored the Lüders epileptogenic zone definition. Experts were divided regarding the seizure onset zone definition, with 179 (55.8%) favoring onset alone and 135 (42.1%) supporting onset and early propagation. In terms of presurgical SEEG features, a clear preference was found for ictal features over interictal features. Seizure onset patterns were identified as the most important features by 265 experts (82.5%). We found similar trends after correcting for nonresponders using regression analysis.

SIGNIFICANCE

This study underscores the need for standardized terminology, interpretation, and outcome assessment in SEEG-informed epilepsy surgery. By highlighting the diverse perspectives and practices in SEEG, this research lays a solid foundation for developing globally accepted terminology and guidelines, advancing the field toward improved communication and standardization in epilepsy surgery.

Interictal intracranial EEG asymmetry lateralizes temporal lobe epilepsy

Patients with drug-resistant temporal lobe epilepsy often undergo intracranial EEG recording to capture multiple seizures in order to lateralize the seizure onset zone. This process is associated with morbidity and often ends in postoperative seizure recurrence. Abundant interictal (between-seizure) data are captured during this process, but these data currently play a small role in surgical planning. Our objective was to predict the laterality of the seizure onset zone using interictal intracranial EEG data in patients with temporal lobe epilepsy. We performed a retrospective cohort study (single-centre study for model development; two-centre study for model validation). We studied patients with temporal lobe epilepsy undergoing intracranial EEG at the University of Pennsylvania (internal cohort) and the Medical University of South Carolina (external cohort) between 2015 and 2022. We developed a logistic regression model to predict seizure onset zone laterality using several interictal EEG features derived from recent publications. We compared the concordance between the model-predicted seizure onset zone laterality and the side of surgery between patients with good and poor surgical outcomes. Forty-seven patients (30 female; ages 20-69; 20 left-sided, 10 right-sided and 17 bilateral seizure onsets) were analysed for model development and internal validation. Nineteen patients (10 female; ages 23-73; 5 left-sided, 10 right-sided, 4 bilateral) were analysed for external validation. The internal cohort cross-validated area under the curve for a model trained using spike rates was 0.83 for a model predicting left-sided seizure onset and 0.68 for a model predicting right-sided seizure onset. Balanced accuracies in the external cohort were 79.3% and 78.9% for the left- and right-sided predictions, respectively. The predicted concordance between the laterality of the seizure onset zone and the side of surgery was higher in patients with good surgical outcome. We replicated the finding that right temporal lobe epilepsy was harder to distinguish in a separate modality of resting-state functional MRI. In conclusion, interictal EEG signatures are distinct across seizure onset zone lateralities. Left-sided seizure onsets are easier to distinguish than right-sided onsets. A model trained on spike rates accurately identifies patients with left-sided seizure onset zones and predicts surgical outcome. A potential clinical application of these findings could be to either support or oppose a hypothesis of unilateral temporal lobe epilepsy when deciding to pursue surgical resection or ablation as opposed to device implantation.

Utility of intracranial EEG networks depends on re-referencing and connectivity choice

Studies of intracranial EEG networks have been used to reveal seizure generators in patients with drug-resistant epilepsy. Intracranial EEG is implanted to capture the epileptic network, the collection of brain tissue that forms a substrate for seizures to start and spread. Interictal intracranial EEG measures brain activity at baseline, and networks computed during this state can reveal aberrant brain tissue without requiring seizure recordings. Intracranial EEG network analyses require choosing a reference and applying statistical measures of functional connectivity. Approaches to these technical choices vary widely across studies, and the impact of these technical choices on downstream analyses is poorly understood. Our objective was to examine the effects of different re-referencing and connectivity approaches on connectivity results and on the ability to lateralize the seizure onset zone in patients with drug-resistant epilepsy. We applied 48 pre-processing pipelines to a cohort of 125 patients with drug-resistant epilepsy recorded with interictal intracranial EEG across two epilepsy centres to generate intracranial EEG functional connectivity networks. Twenty-four functional connectivity measures across time and frequency domains were applied in combination with common average re-referencing or bipolar re-referencing. We applied an unsupervised clustering algorithm to identify groups of pre-processing pipelines. We subjected each pre-processing approach to three quality tests: (i) the introduction of spurious correlations; (ii) robustness to incomplete spatial sampling; and (iii) the ability to lateralize the clinician-defined seizure onset zone. Three groups of similar pre-processing pipelines emerged: common average re-referencing pipelines, bipolar re-referencing pipelines and relative entropy-based connectivity pipelines. Relative entropy and common average re-referencing networks were more robust to incomplete electrode sampling than bipolar re-referencing and other connectivity methods (Friedman test, Dunn-Šidák test P < 0.0001). Bipolar re-referencing reduced spurious correlations at non-adjacent channels better than common average re-referencing (Δ mean from machine ref = -0.36 versus -0.22) and worse in adjacent channels (Δ mean from machine ref = -0.14 versus -0.40). Relative entropy-based network measures lateralized the seizure onset hemisphere better than other measures in patients with temporal lobe epilepsy (Benjamini-Hochberg-corrected P < 0.05, Cohen's d: 0.60-0.76). Finally, we present an interface where users can rapidly evaluate intracranial EEG pre-processing choices to select the optimal pre-processing methods tailored to specific research questions. The choice of pre-processing methods affects downstream network analyses. Choosing a single method among highly correlated approaches can reduce redundancy in processing. Relative entropy outperforms other connectivity methods in multiple quality tests. We present a method and interface for researchers to optimize their pre-processing methods for deriving intracranial EEG brain networks.

The value of linear and non-linear quantitative EEG analysis in paediatric epilepsy surgery: a machine learning approach

Epilepsy surgery is effective for patients with medication-resistant seizures, however 20-40% of them are not seizure free after surgery. Aim of this study is to evaluate the role of linear and non-linear EEG features to predict post-surgical outcome. We included 123 paediatric patients who underwent epilepsy surgery at Bambino Gesù Children Hospital (January 2009-April 2020). All patients had long term video-EEG monitoring. We analysed 1-min scalp interictal EEG (wakefulness and sleep) and extracted 13 linear and non-linear EEG features (power spectral density (PSD), Hjorth, approximate entropy, permutation entropy, Lyapunov and Hurst value). We used a logistic regression (LR) as feature selection process. To quantify the correlation between EEG features and surgical outcome we used an artificial neural network (ANN) model with 18 architectures. LR revealed a significant correlation between PSD of alpha band (sleep), Mobility index (sleep) and the Hurst value (sleep and awake) with outcome. The fifty-four ANN models gave a range of accuracy (46-65%) in predicting outcome. Within the fifty-four ANN models, we found a higher accuracy (64.8% ± 7.6%) in seizure outcome prediction, using features selected by LR. The combination of PSD of alpha band, mobility and the Hurst value positively correlate with good surgical outcome.

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.

Recent advances in clinical electroencephalography

PURPOSE OF REVIEW

Clinical electroencephalography (EEG) is a conservative medical field. This explains likely the significant gap between clinical practice and new research developments. This narrative review discusses possible causes of this discrepancy and how to circumvent them. More specifically, we summarize recent advances in three applications of clinical EEG: source imaging (ESI), high-frequency oscillations (HFOs) and EEG in critically ill patients.

RECENT FINDINGS

Recently published studies on ESI provide further evidence for the accuracy and clinical utility of this method in the multimodal presurgical evaluation of patients with drug-resistant focal epilepsy, and opened new possibilities for further improvement of the accuracy. HFOs have received much attention as a novel biomarker in epilepsy. However, recent studies questioned their clinical utility at the level of individual patients. We discuss the impediments, show up possible solutions and highlight the perspectives of future research in this field. EEG in the ICU has been one of the major driving forces in the development of clinical EEG. We review the achievements and the limitations in this field.

SUMMARY

This review will promote clinical implementation of recent advances in EEG, in the fields of ESI, HFOs and EEG in the intensive care.

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.

Interictal intracranial EEG asymmetry lateralizes temporal lobe epilepsy

Patients with drug-resistant temporal lobe epilepsy often undergo intracranial EEG recording to capture multiple seizures in order to lateralize the seizure onset zone. This process is associated with morbidity and often ends in postoperative seizure recurrence. Abundant interictal (between-seizure) data is captured during this process, but these data currently play a small role in surgical planning. Our objective was to predict the laterality of the seizure onset zone using interictal (between-seizure) intracranial EEG data in patients with temporal lobe epilepsy. We performed a retrospective cohort study (single-center study for model development; two-center study for model validation). We studied patients with temporal lobe epilepsy undergoing intracranial EEG at the University of Pennsylvania (internal cohort) and the Medical University of South Carolina (external cohort) between 2015 and 2022. We developed a logistic regression model to predict seizure onset zone laterality using interictal EEG. We compared the concordance between the model-predicted seizure onset zone laterality and the side of surgery between patients with good and poor surgical outcomes. 47 patients (30 women; ages 20-69; 20 left-sided, 10 right-sided, and 17 bilateral seizure onsets) were analyzed for model development and internal validation. 19 patients (10 women; ages 23-73; 5 left-sided, 10 right-sided, 4 bilateral) were analyzed for external validation. The internal cohort cross-validated area under the curve for a model trained using spike rates was 0.83 for a model predicting left-sided seizure onset and 0.68 for a model predicting right-sided seizure onset. Balanced accuracies in the external cohort were 79.3% and 78.9% for the left- and right-sided predictions, respectively. The predicted concordance between the laterality of the seizure onset zone and the side of surgery was higher in patients with good surgical outcome. In conclusion, interictal EEG signatures are distinct across seizure onset zone lateralities. Left-sided seizure onsets are easier to distinguish than right-sided onsets. A model trained on spike rates accurately identifies patients with left-sided seizure onset zones and predicts surgical outcome.