Trends in the use of automated algorithms for the detection of high-frequency oscillations associated with human epilepsy

Detection of High-Frequency Oscillations from Intracranial EEG Data with Switching State Space Model

High Frequency Oscillations (HFOs) is an important biomarker that can potentially pinpoint the epileptogenic zones (EZs). However, the duration of HFO is short with around 4 cycles, which might be hard to recognize when embedded within signals of lower frequency oscillatory background. In addition, annotating HFOs manually can be time-consuming given long-time recordings and up to hundreds of intracranial electrodes. We propose to leverage a Switching State Space Model (SSSM) to identify the HFOs events automatically and instantaneously without relying on extracting features from sliding windows. The effectiveness of the SSSM for HFOs detection is fully validated in the intracranial EEG recording from human subjects undergoing the presurgical evaluations and showed improved accuracy when capturing the HFOs occurrence and their duration.

Robust compression and detection of epileptiform patterns in ECoG using a real-time spiking neural network hardware framework

Interictal Epileptiform Discharges (IED) and High Frequency Oscillations (HFO) in intraoperative electrocorticography (ECoG) may guide the surgeon by delineating the epileptogenic zone. We designed a modular spiking neural network (SNN) in a mixed-signal neuromorphic device to process the ECoG in real-time. We exploit the variability of the inhomogeneous silicon neurons to achieve efficient sparse and decorrelated temporal signal encoding. We interface the full-custom SNN device to the BCI2000 real-time framework and configure the setup to detect HFO and IED co-occurring with HFO (IED-HFO). We validate the setup on pre-recorded data and obtain HFO rates that are concordant with a previously validated offline algorithm (Spearman's ρ = 0.75, p = 1e-4), achieving the same postsurgical seizure freedom predictions for all patients. In a remote on-line analysis, intraoperative ECoG recorded in Utrecht was compressed and transferred to Zurich for SNN processing and successful IED-HFO detection in real-time. These results further demonstrate how automated remote real-time detection may enable the use of HFO in clinical practice.

Epileptic tissue localization using graph-based networks in the high frequency oscillation range of intracranial electroencephalography

PURPOSE

High frequency oscillations (HFOs) are an emerging biomarker of epilepsy. However, very few studies have investigated the functional connectivity of interictal iEEG signals in the frequency range of HFOs. Here, we study the corresponding functional networks using graph theory, and we assess their predictive value for automatic electrode classification in a cohort of 20 drug resistant patients.

METHODS

Coherence-based connectivity analysis was performed on the iEEG recordings, and six different local graph measures were computed in both sub-bands of the HFO frequency range (80-250 Hz and 250-500 Hz). Correlation analysis was implemented between the local graph measures and the ripple and fast ripple rates. Finally, the WEKA software was employed for training and testing different predictive models on the aforementioned local graph measures.

RESULTS

The ripple rate was significantly correlated with five out of six local graph measures in the functional network. For fast ripples, their rate was also significantly (but negatively) correlated with most of the local metrics. The results from WEKA showed that the Logistic Regression algorithm was able to classify highly HFO-contaminated electrodes with an accuracy of 82.5 % for ripples and 75.4 % for fast ripples.

CONCLUSION

Functional connectivity networks in the HFO band could represent an alternative to the direct use of distinct HFO events, while also providing important insights about hub epileptic areas that can represent possible surgical targets. Automatic electrode classification through FC-based classifiers can help bypass the burden of manual HFO annotation, providing at the same time similar amount of information about the epileptic tissue.

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

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

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.

Optimizing detection and deep learning-based classification of pathological high-frequency oscillations in epilepsy

OBJECTIVE

This study aimed to explore sensitive detection methods for pathological high-frequency oscillations (HFOs) to improve seizure outcomes in epilepsy surgery.

METHODS

We analyzed interictal HFOs (80-500 Hz) in 15 children with medication-resistant focal epilepsy who underwent chronic intracranial electroencephalogram via subdural grids. The HFOs were assessed using the short-term energy (STE) and Montreal Neurological Institute (MNI) detectors and examined for spike association and time-frequency plot characteristics. A deep learning (DL)-based classification was applied to purify pathological HFOs. Postoperative seizure outcomes were correlated with HFO-resection ratios to determine the optimal HFO detection method.

RESULTS

The MNI detector identified a higher percentage of pathological HFOs than the STE detector, but some pathological HFOs were detected only by the STE detector. HFOs detected by both detectors had the highest spike association rate. The Union detector, which detects HFOs identified by either the MNI or STE detector, outperformed other detectors in predicting postoperative seizure outcomes using HFO-resection ratios before and after DL-based purification.

CONCLUSIONS

HFOs detected by standard automated detectors displayed different signal and morphological characteristics. DL-based classification effectively purified pathological HFOs.

SIGNIFICANCE

Enhancing the detection and classification methods of HFOs will improve their utility in predicting postoperative seizure outcomes.

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.

Optimizing Detection and Deep Learning-based Classification of Pathological High-Frequency Oscillations in Epilepsy

Objective

This study aimed to explore sensitive detection methods and deep learning (DL)-based classification for pathological high-frequency oscillations (HFOs).

Methods

We analyzed interictal HFOs (80-500 Hz) in 15 children with medication-resistant focal epilepsy who underwent resection after chronic intracranial electroencephalogram via subdural grids. The HFOs were assessed using the short-term energy (STE) and Montreal Neurological Institute (MNI) detectors and examined for pathological features based on spike association and time-frequency plot characteristics. A DL-based classification was applied to purify pathological HFOs. Postoperative seizure outcomes were correlated with HFO-resection ratios to determine the optimal HFO detection method.

Results

The MNI detector identified a higher percentage of pathological HFOs than the STE detector, but some pathological HFOs were detected only by the STE detector. HFOs detected by both detectors exhibited the most pathological features. The Union detector, which detects HFOs identified by either the MNI or STE detector, outperformed other detectors in predicting postoperative seizure outcomes using HFO-resection ratios before and after DL-based purification.

Conclusions

HFOs detected by standard automated detectors displayed different signal and morphological characteristics. DL-based classification effectively purified pathological HFOs.

Significance

Enhancing the detection and classification methods of HFOs will improve their utility in predicting postoperative seizure outcomes.

HIGHLIGHTS

HFOs detected by the MNI detector showed different traits and higher pathological bias than those detected by the STE detectorHFOs detected by both MNI and STE detectors (the Intersection HFOs) were deemed the most pathologicalA deep learning-based classification was able to distill pathological HFOs, regard-less of the initial HFO detection methods.

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

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

Source-sink connectivity: a novel interictal EEG marker for seizure localization

Over 15 million epilepsy patients worldwide have drug-resistant epilepsy. Successful surgery is a standard of care treatment but can only be achieved through complete resection or disconnection of the epileptogenic zone, the brain region(s) where seizures originate. Surgical success rates vary between 20% and 80%, because no clinically validated biological markers of the epileptogenic zone exist. Localizing the epileptogenic zone is a costly and time-consuming process, which often requires days to weeks of intracranial EEG (iEEG) monitoring. Clinicians visually inspect iEEG data to identify abnormal activity on individual channels occurring immediately before seizures or spikes that occur interictally (i.e. between seizures). In the end, the clinical standard mainly relies on a small proportion of the iEEG data captured to assist in epileptogenic zone localization (minutes of seizure data versus days of recordings), missing opportunities to leverage these largely ignored interictal data to better diagnose and treat patients. IEEG offers a unique opportunity to observe epileptic cortical network dynamics but waiting for seizures increases patient risks associated with invasive monitoring. In this study, we aimed to leverage interictal iEEG data by developing a new network-based interictal iEEG marker of the epileptogenic zone. We hypothesized that when a patient is not clinically seizing, it is because the epileptogenic zone is inhibited by other regions. We developed an algorithm that identifies two groups of nodes from the interictal iEEG network: those that are continuously inhibiting a set of neighbouring nodes ('sources') and the inhibited nodes themselves ('sinks'). Specifically, patient-specific dynamical network models were estimated from minutes of iEEG and their connectivity properties revealed top sources and sinks in the network, with each node being quantified by source-sink metrics. We validated the algorithm in a retrospective analysis of 65 patients. The source-sink metrics identified epileptogenic regions with 73% accuracy and clinicians agreed with the algorithm in 93% of seizure-free patients. The algorithm was further validated by using the metrics of the annotated epileptogenic zone to predict surgical outcomes. The source-sink metrics predicted outcomes with an accuracy of 79% compared to an accuracy of 43% for clinicians' predictions (surgical success rate of this dataset). In failed outcomes, we identified brain regions with high metrics that were untreated. When compared with high frequency oscillations, the most commonly proposed interictal iEEG feature for epileptogenic zone localization, source-sink metrics outperformed in predictive power (by a factor of 1.2), suggesting they may be an interictal iEEG fingerprint of the epileptogenic zone.

An approach for reliably identifying high-frequency oscillations and reducing false-positive detections

OBJECTIVE

Aiming to improve the feasibility and reliability of using high-frequency oscillations (HFOs) for translational studies of epilepsy, we present a pipeline with features specifically designed to reject false positives for HFOs to improve the automatic HFO detector.

METHODS

We presented an integrated, multi-layered procedure capable of automatically rejecting HFOs from a variety of common false positives, such as motion, background signals, and sharp transients. This method utilizes a time-frequency contour approach that embeds three different layers including peak constraints, power thresholds, and morphological identification to discard false positives. Four experts were involved in rating detected HFO events that were randomly selected from different posttraumatic epilepsy (PTE) animals for a comprehensive evaluation.

RESULTS

The algorithm was run on 768-h recordings of intracranial electrodes in 48 PTE animals. A total of 453 917 HFOs were identified by initial HFO detection, of which 450 917 were implemented for HFO refinement and 203 531 events were retained. Random sampling was used to evaluate the performance of the detector. The HFO detection yielded an overall accuracy of 0.95 ± 0.03 , with precision, recall, and F1 scores of 0.92 ± 0.05 , 0.99 ± 0.01 , and 0.94 ± 0.03 , respectively. For the HFO classification, our algorithm obtained an accuracy of 0.97 ± 0.02 . For the inter-rater reliability of algorithm evaluation, the agreement among four experts was 0.94 ± 0.03 for HFO detection and 0.85 ± 0.04 for HFO classification.

SIGNIFICANCE

Our approach shows that a segregated pipeline design with a focus on false-positive rejection can improve the detection efficiency and provide reliable results. This pipeline does not require customization and uses fixed parameters, making it highly feasible and translatable for basic and clinical applications of epilepsy.

EEG biomarkers for the diagnosis and treatment of infantile spasms

Early diagnosis and treatment are critical for young children with infantile spasms (IS), as this maximizes the possibility of the best possible child-specific outcome. However, there are major barriers to achieving this, including high rates of misdiagnosis or failure to recognize the seizures, medication failure, and relapse. There are currently no validated tools to aid clinicians in assessing objective diagnostic criteria, predicting or measuring medication response, or predicting the likelihood of relapse. However, the pivotal role of EEG in the clinical management of IS has prompted many recent studies of potential EEG biomarkers of the disease. These include both visual EEG biomarkers based on human visual interpretation of the EEG and computational EEG biomarkers in which computers calculate quantitative features of the EEG. Here, we review the literature on both types of biomarkers, organized based on the application (diagnosis, treatment response, prediction, etc.). Visual biomarkers include the assessment of hypsarrhythmia, epileptiform discharges, fast oscillations, and the Burden of AmplitudeS and Epileptiform Discharges (BASED) score. Computational markers include EEG amplitude and power spectrum, entropy, functional connectivity, high frequency oscillations (HFOs), long-range temporal correlations, and phase-amplitude coupling. We also introduce each of the computational measures and provide representative examples. Finally, we highlight remaining gaps in the literature, describe practical guidelines for future biomarker discovery and validation studies, and discuss remaining roadblocks to clinical implementation, with the goal of facilitating future work in this critical area.

Elimination of pseudo-HFOs in iEEG using sparse representation and Random Forest classifier

High-Frequency Oscillation (HFO) is a promising biomarker of the epileptogenic zone. However, sharp artifacts might easily pass the conventional HFO detectors as real HFOs and reduce the seizure onset zone (SOZ) localization. We hypothesize that, unlike pseudo-HFOs, which originates from artifacts with sharp changes or arbitrary waveform characteristic, real HFOs could be represented by a limited number of oscillatory waveforms. Accordingly, to distinguish true ones from pseudo-HFOs, we established a new classification method based on sparse representation of candidate events that passed an initial detector with high sensitivity but low specificity. Specifically, using the Orthogonal Matching Pursuit (OMP) and a redundant Gabor dictionary, each event was represented sparsely in an iterative fashion. The approximation error was estimated over 30 iterations which were concatenated to form a 30-dimensional feature vector and fed to a random forest classifier. Based on the selected dictionary elements, our method can further classify HFOs into Ripples (R) and Fast Ripples (FR). In this scheme, two experts visually inspected 2075 events captured in iEEG recordings from 5 different subjects and labeled them as true-HFO or Pseudo-HFO. We reached 90.22% classification accuracy in labeled events and a 21.16% SOZ localization improvement compared to the conventional amplitude-threshold-based detector. Our sparse representation framework also classified the detected HFOs into R and FR subcategories. We reached 91.24% SOZ accuracy with the detected [Formula: see text] events. Clinical Relevance---This sparse representation framework establishes a new approach to distinguish real from pseudo-HFOs in prolonged iEEG recordings. It also provides reliable SOZ identification without the selection of artifact-free segments.

Protocol for multicentre comparison of interictal high-frequency oscillations as a predictor of seizure freedom

In drug-resistant focal epilepsy, interictal high-frequency oscillations (HFOs) recorded from intracranial EEG (iEEG) may provide clinical information for delineating epileptogenic brain tissue. The iEEG electrode contacts that contain HFO are hypothesized to delineate the epileptogenic zone; their resection should then lead to postsurgical seizure freedom. We test whether our prospective definition of clinically relevant HFO is in agreement with postsurgical seizure outcome. The algorithm is fully automated and is equally applied to all data sets. The aim is to assess the reliability of the proposed detector and analysis approach. We use an automated data-independent prospective definition of clinically relevant HFO that has been validated in data from two independent epilepsy centres. In this study, we combine retrospectively collected data sets from nine independent epilepsy centres. The analysis is blinded to clinical outcome. We use iEEG recordings during NREM sleep with a minimum of 12 epochs of 5 min of NREM sleep. We automatically detect HFO in the ripple (80-250 Hz) and in the fast ripple (250-500 Hz) band. There is no manual rejection of events in this fully automated algorithm. The type of HFO that we consider clinically relevant is defined as the simultaneous occurrence of a fast ripple and a ripple. We calculate the temporal consistency of each patient's HFO rates over several data epochs within and between nights. Patients with temporal consistency <50% are excluded from further analysis. We determine whether all electrode contacts with high HFO rate are included in the resection volume and whether seizure freedom (ILAE 1) was achieved at ≥2 years follow-up. Applying a previously validated algorithm to a large cohort from several independent epilepsy centres may advance the clinical relevance and the generalizability of HFO analysis as essential next step for use of HFO in clinical practice.

Automatic Detection of High-Frequency Oscillations With Neuromorphic Spiking Neural Networks

Interictal high-frequency oscillations (HFO) detected in electroencephalography recordings have been proposed as biomarkers of epileptogenesis, seizure propensity, disease severity, and treatment response. Automatic HFO detectors typically analyze the data offline using complex time-consuming algorithms, which limits their clinical application. Neuromorphic circuits offer the possibility of building compact and low-power processing systems that can analyze data on-line and in real time. In this review, we describe a fully automated detection pipeline for HFO that uses, for the first time, spiking neural networks and neuromorphic technology. We demonstrated that our HFO detection pipeline can be applied to recordings from different modalities (intracranial electroencephalography, electrocorticography, and scalp electroencephalography) and validated its operation in a custom-designed neuromorphic processor. Our HFO detection approach resulted in high accuracy and specificity in the prediction of seizure outcome in patients implanted with intracranial electroencephalography and electrocorticography, and in the prediction of epilepsy severity in patients recorded with scalp electroencephalography. Our research provides a further step toward the real-time detection of HFO using compact and low-power neuromorphic devices. The real-time detection of HFO in the operation room may improve the seizure outcome of epilepsy surgery, while the use of our neuromorphic processor for non-invasive therapy monitoring might allow for more effective medication strategies to achieve seizure control. Therefore, this work has the potential to improve the quality of life in patients with epilepsy by improving epilepsy diagnostics and treatment.

Beyond rates: time-varying dynamics of high frequency oscillations as a biomarker of the seizure onset zone

Objective. High frequency oscillations (HFOs) recorded by intracranial electrodes have generated excitement for their potential to help localize epileptic tissue for surgical resection. However, the number of HFOs per minute (i.e. the HFO 'rate') is not stable over the duration of intracranial recordings; for example, the rate of HFOs increases during periods of slow-wave sleep. Moreover, HFOs that are predictive of epileptic tissue may occur in oscillatory patterns due to phase coupling with lower frequencies. Therefore, we sought to further characterize between-seizure (i.e. 'interictal') HFO dynamics both within and outside the seizure onset zone (SOZ).Approach. Using long-term intracranial EEG (mean duration 10.3 h) from 16 patients, we automatically detected HFOs using a new algorithm. We then fit a hierarchical negative binomial model to the HFO counts. To account for differences in HFO dynamics and rates between sleep and wakefulness, we also fit a mixture model to the same data that included the ability to switch between two discrete brain states that were automatically determined during the fitting process. The ability to predict the SOZ by model parameters describing HFO dynamics (i.e. clumping coefficients and coefficients of variation) was assessed using receiver operating characteristic curves.Main results. Parameters that described HFO dynamics were predictive of SOZ. In fact, these parameters were found to be more consistently predictive than HFO rate. Using concurrent scalp EEG in two patients, we show that the model-found brain states corresponded to (1) non-REM sleep and (2) awake and rapid eye movement sleep. However the brain state most likely corresponding to slow-wave sleep in the second model improved SOZ prediction compared to the first model for only some patients.Significance. This work suggests that delineation of SOZ with interictal data can be improved by the inclusion of time-varying HFO dynamics.

Analysis of power spectrum and phase lag index changes following deep brain stimulation of the anterior nucleus of the thalamus in patients with drug-resistant epilepsy: A retrospective study

OBJECTIVES

The mechanisms underlying the anterior nucleus of the thalamus (ANT) deep brain stimulation (DBS) for the treatment of drug-resistant epilepsy (DRE) have not been fully explored. The present study aimed to measure the changes in whole-brain activity generated by ANT DBS using interictal electroencephalography (EEG).

MATERIALS AND METHODS

Interictal EEG signals were retrospectively collected in 20 DRE patients who underwent ANT DBS surgery. Patients were classified as responders or non-responders depending on their response to ANT DBS treatment. The power spectrum (PS) and Phase Lag Index (PLI) were determined and data analyzed using a paired sample t-test to evaluate activity differences between pre-and-post-treatment on different frequency categories. Student's t-test, Mann-Whitney test (non-parametric test) and Fisher exact test were used to compare groups in terms of clinical variables and EEG metrics. P values < 0.05 were considered statistically significant, and FDR-corrected values were used for multiple testing.

RESULTS

PS analysis revealed that whole-brain spectral power had a significant decrease in the beta (p = 0.005) and gamma (p = 0.037) bands following ANT DBS treatment in responders. The analysis of scalp topographic images of all patients showed that ANT DBS decreases PS in the beta band at the F3, F7 and Cz electrode sites. The findings indicated a decrease in PS in the gamma band at the Fp2, F3, Cz, T3, T5 and Oz electrode sites. After ANT DBS treatment, PLI analysis showed a significant decrease in PLI between Fp1 and T3 in the gamma band in responders.

CONCLUSION

The findings showed that ANT DBS induces a decrease in power in the left frontal lobe, left temporal lobe and midline areas in the beta and gamma bands. Lower whole-brain power in the beta and gamma bands can be used as biomarkers for a favorable therapeutic response to ANT DBS, and decreased synchronization between the left frontal pole and temporal lobe in the gamma band can also be used as a biomarker for effective clinical stimulation to guide postoperative programming.

BrainQuake: An Open-Source Python Toolbox for the Stereoelectroencephalography Spatiotemporal Analysis

Intracranial stereoelectroencephalography (SEEG) is broadly used in the presurgical evaluation of intractable epilepsy, due to its high temporal resolution in neural activity recording and high spatial resolution within suspected epileptogenic zones. Neurosurgeons or technicians face the challenge of conducting a workflow of post-processing operations with the multimodal data (e.g., MRI, CT, and EEG) after the implantation surgery, such as brain surface reconstruction, electrode contact localization, and SEEG data analysis. Several software or toolboxes have been developed to take one or more steps in the workflow but without an end-to-end solution. In this study, we introduced BrainQuake, an open-source Python software for the SEEG spatiotemporal analysis, integrating modules and pipelines in surface reconstruction, electrode localization, seizure onset zone (SOZ) prediction based on ictal and interictal SEEG analysis, and final visualizations, each of which is highly automated with a user-friendly graphical user interface (GUI). BrainQuake also supports remote communications with a public server, which is facilitated with automated and standardized preprocessing pipelines, high-performance computing power, and data curation management to provide a time-saving and compatible platform for neurosurgeons and researchers.

Automatic Detection of High-Frequency Oscillations Based on an End-to-End Bi-Branch Neural Network and Clinical Cross-Validation

Accurate identification of high-frequency oscillation (HFO) is an important prerequisite for precise localization of epileptic foci and good prognosis of drug-refractory epilepsy. Exploring a high-performance automatic detection method for HFOs can effectively help clinicians reduce the error rate and reduce manpower. Due to the limited analysis perspective and simple model design, it is difficult to meet the requirements of clinical application by the existing methods. Therefore, an end-to-end bi-branch fusion model is proposed to automatically detect HFOs. With the filtered band-pass signal (signal branch) and time-frequency image (TFpic branch) as the input of the model, two backbone networks for deep feature extraction are established, respectively. Specifically, a hybrid model based on ResNet1d and long short-term memory (LSTM) is designed for signal branch, which can focus on both the features in time and space dimension, while a ResNet2d with a Convolutional Block Attention Module (CBAM) is constructed for TFpic branch, by which more attention is paid to useful information of TF images. Then the outputs of two branches are fused to realize end-to-end automatic identification of HFOs. Our method is verified on 5 patients with intractable epilepsy. In intravalidation, the proposed method obtained high sensitivity of 94.62%, specificity of 92.7%, and F1-score of 93.33%, and in cross-validation, our method achieved high sensitivity of 92.00%, specificity of 88.26%, and F1-score of 89.11% on average. The results show that the proposed method outperforms the existing detection paradigms of either single signal or single time-frequency diagram strategy. In addition, the average kappa coefficient of visual analysis and automatic detection results is 0.795. The method shows strong generalization ability and high degree of consistency with the gold standard meanwhile. Therefore, it has great potential to be a clinical assistant tool.

High-frequency oscillations detected by electroencephalography as biomarkers to evaluate treatment outcome, mirror pathological severity and predict susceptibility to epilepsy

High-frequency oscillations (HFOs) in the electroencephalography (EEG) have been extensively investigated as a potential biomarker of epileptogenic zones. The understanding of the role of HFOs in epilepsy has been advanced considerably over the past decade, and the use of scalp EEG facilitates recordings of HFOs. HFOs were initially applied in large scale in epilepsy surgery and are now being utilized in other applications. In this review, we summarize applications of HFOs in 3 subtopics: (1) HFOs as biomarkers to evaluate epilepsy treatment outcome; (2) HFOs as biomarkers to measure seizure propensity; (3) HFOs as biomarkers to reflect the pathological severity of epilepsy. Nevertheless, knowledge regarding the above clinical applications of HFOs remains limited at present. Further validation through prospective studies is required for its reliable application in the clinical management of individual epileptic patients.

HFOApp: A MATLAB Graphical User Interface for High-Frequency Oscillation Marking

Epilepsy affects 3.4 million people in the United States, and, despite the availability of numerous antiepileptic drugs, 36% of patients have uncontrollable seizures, which severely impact quality of life. High-frequency oscillations (HFOs) are a potential biomarker of epileptogenic tissue that could be useful in surgical planning. As a result, research into the efficacy of HFOs as a clinical tool has increased over the last 2 decades. However, detection and identification of these transient rhythms in intracranial electroencephalographic recordings remain time-consuming and challenging. Although automated detection algorithms have been developed, their results are widely inconsistent, reducing reliability. Thus, manual marking of HFOs remains the gold standard, and manual review of automated results is required. However, manual marking and review are time consuming and can still produce variable results because of their subjective nature and the limitations in functionality of existing open-source software. Our goal was to develop a new software with broad application that improves on existing open-source HFO detection applications in usability, speed, and accuracy. Here, we present HFOApp: a free, open-source, easy-to-use MATLAB-based graphical user interface for HFO marking. This toolbox offers a high degree of intuitive and ergonomic usability and integrates interactive automation-assist options with manual marking, significantly reducing the time needed for review and manual marking of recordings, while increasing inter-rater reliability. The toolbox also features simultaneous multichannel detection and marking. HFOApp was designed as an easy-to-use toolbox for clinicians and researchers to quickly and accurately mark, quantify, and characterize HFOs within electrophysiological datasets.

Blinded study: prospectively defined high-frequency oscillations predict seizure outcome in individual patients

Interictal high-frequency oscillations are discussed as biomarkers for epileptogenic brain tissue that should be resected in epilepsy surgery to achieve seizure freedom. The prospective classification of tissue sampled by individual electrode contacts remains a challenge. We have developed an automated, prospective definition of clinically relevant high-frequency oscillations in intracranial EEG from Montreal and tested it in recordings from Zurich. We here validated the algorithm on intracranial EEG that was recorded in an independent epilepsy centre so that the analysis was blinded to seizure outcome. We selected consecutive patients who underwent resective epilepsy surgery in Geneva with post-surgical follow-up > 12 months. We analysed long-term recordings during sleep that we segmented into intervals of 5 min. High-frequency oscillations were defined in the ripple (80–250 Hz) and the fast ripple (250–500 Hz) frequency bands. Contacts with the highest rate of ripples co-occurring with fast ripples designated the relevant area. As a validity criterion, we calculated the test–retest reliability of the high-frequency oscillations area between the 5 min intervals (dwell time ≥50%). If the area was not fully resected and the patient suffered from recurrent seizures, this was classified as a true positive prediction. We included recordings from 16 patients (median age 32 years, range 18–53 years) with stereotactic depth electrodes and/or with subdural electrode grids (median follow-up 27 months, range 12–55 months). For each patient, we included several 5 min intervals (median 17 intervals). The relevant area had high test–retest reliability across intervals (median dwell time 95%). In two patients, the test–retest reliability was too low (dwell time < 50%) so that outcome prediction was not possible. The area was fully included in the resected volume in 2/4 patients who achieved post-operative seizure freedom (specificity 50%) and was not fully included in 9/10 patients with recurrent seizures (sensitivity 90%), leading to an accuracy of 79%. An additional exploratory analysis suggested that high-frequency oscillations were associated with interictal epileptic discharges only in channels within the relevant area and not associated in channels outside the area. We thereby validated the automated procedure to delineate the clinically relevant area in each individual patient of an independently recorded dataset and achieved the same good accuracy as in our previous studies. The reproducibility of our results across datasets is promising for a multicentre study to test the clinical application of high-frequency oscillations to guide epilepsy surgery.

Improving the ripple classification in focal pediatric epilepsy: identifying pathological high-frequency oscillations by Gaussian mixture model clustering

Objective. High-frequency oscillations (HFOs) have emerged as a promising clinical biomarker for presurgical evaluation in childhood epilepsy. HFOs are commonly classified in stereo-encephalography as ripples (80-200 Hz) and fast ripples (200-500 Hz). Ripples are less specific and not so directly associated with epileptogenic activity because of their physiological and pathological origin. The aim of this paper is to distinguish HFOs in the ripple band and to improve the evaluation of the epileptogenic zone (EZ).Approach. This study constitutes a novel modeling approach evaluated in ten patients from Sant Joan de Deu Pediatric Hospital (Barcelona, Spain), with clearly-defined seizure onset zones (SOZ) during presurgical evaluation. A subject-by-subject basis analysis is proposed: a probabilistic Gaussian mixture model (GMM) based on the combination of specific ripple features is applied for estimating physiological and pathological ripple subpopulations.Main Results. Clear pathological and physiological ripples are identified. Features differ considerably among patients showing within-subject variability, suggesting that individual models are more appropriate than a traditional whole-population approach. The difference in rates inside and outside the SOZ for pathological ripples is significantly higher than when considering all the ripples. These significant differences also appear in signal segments without epileptiform activity. Pathological ripple rates show a sharp decline from SOZ to non-SOZ contacts and a gradual decrease with distance.Significance. This novel individual GMM approach improves ripple classification and helps to refine the delineation of the EZ, as well as being appropriate to investigate the interaction of epileptogenic and propagation networks.

Application of a convolutional neural network for fully-automated detection of spike ripples in the scalp electroencephalogram

BACKGROUND

A reliable biomarker to identify cortical tissue responsible for generating epileptic seizures is required to guide prognosis and treatment in epilepsy. Combined spike ripple events are a promising biomarker for epileptogenic tissue that currently require expert review for accurate identification. This expert review is time consuming and subjective, limiting reproducibility and high-throughput applications.

NEW METHOD

To address this limitation, we develop a fully-automated method for spike ripple detection. The method consists of a convolutional neural network trained to compute the probability that a spectrogram image contains a spike ripple.

RESULTS

We validate the proposed spike ripple detector on expert-labeled data and show that this detector accurately separates subjects with low and high seizure risks.

COMPARISON WITH EXISTING METHOD

The proposed method performs as well as existing methods that require manual validation of candidate spike ripple events. The introduction of a fully automated method reduces subjectivity and increases rigor and reproducibility of this epilepsy biomarker.

CONCLUSION

We introduce and validate a fully-automated spike ripple detector to support utilization of this epilepsy biomarker in clinical and translational work.

An electronic neuromorphic system for real-time detection of high frequency oscillations (HFO) in intracranial EEG

The analysis of biomedical signals for clinical studies and therapeutic applications can benefit from embedded devices that can process these signals locally and in real-time. An example is the analysis of intracranial EEG (iEEG) from epilepsy patients for the detection of High Frequency Oscillations (HFO), which are a biomarker for epileptogenic brain tissue. Mixed-signal neuromorphic circuits offer the possibility of building compact and low-power neural network processing systems that can analyze data on-line in real-time. Here we present a neuromorphic system that combines a neural recording headstage with a spiking neural network (SNN) processing core on the same die for processing iEEG, and show how it can reliably detect HFO, thereby achieving state-of-the-art accuracy, sensitivity, and specificity. This is a first feasibility study towards identifying relevant features in iEEG in real-time using mixed-signal neuromorphic computing technologies.

Objective interictal electrophysiology biomarkers optimize prediction of epilepsy surgery outcome

Researchers have looked for rapidly- and objectively-measurable electrophysiology biomarkers that accurately localize the epileptogenic zone. Promising candidates include interictal high-frequency oscillation and phase-amplitude coupling. Investigators have independently created the toolboxes that compute the high-frequency oscillation rate and the severity of phase-amplitude coupling. This study of 135 patients determined what toolboxes and analytic approaches would optimally classify patients achieving post-operative seizure control. Four different detector toolboxes computed the rate of high-frequency oscillation at ≥80 Hz at intracranial EEG channels. Another toolbox calculated the modulation index reflecting the strength of phase-amplitude coupling between high-frequency oscillation and slow-wave at 3–4 Hz. We defined the completeness of resection of interictally-abnormal regions as the subtraction of high-frequency oscillation rate (or modulation index) averaged across all preserved sites from that averaged across all resected sites. We computed the outcome classification accuracy of the logistic regression-based standard model considering clinical, ictal intracranial EEG and neuroimaging variables alone. We then determined how well the incorporation of high-frequency oscillation/modulation index would improve the standard model mentioned above. To assess the anatomical variability across non-epileptic sites, we generated the normative atlas of detector-specific high-frequency oscillation and modulation index. Each atlas allowed us to compute the statistical deviation of high-frequency oscillation/modulation index from the non-epileptic mean. We determined whether the model accuracy would be improved by incorporating absolute or normalized high-frequency oscillation/modulation index as a biomarker assessing interictally-abnormal regions. We finally determined whether the model accuracy would be improved by selectively incorporating high-frequency oscillation verified to have high-frequency oscillatory components unattributable to a high-pass filtering effect. Ninety-five patients achieved successful seizure control, defined as International League against Epilepsy class 1 outcome. Multivariate logistic regression analysis demonstrated that complete resection of interictally-abnormal regions additively increased the chance of success. The model accuracy was further improved by incorporating z-score normalized high-frequency oscillation/modulation index or selective incorporation of verified high-frequency oscillation. The standard model had a classification accuracy of 0.75. Incorporation of normalized high-frequency oscillation/modulation index or verified high-frequency oscillation improved the classification accuracy up to 0.82. These outcome prediction models survived the cross-validation process and demonstrated an agreement between the model-based likelihood of success and the observed success on an individual basis. Interictal high-frequency oscillation and modulation index had a comparably additive utility in epilepsy presurgical evaluation. Our empirical data support the theoretical notion that the prediction of post-operative seizure outcomes can be optimized with the consideration of both interictal and ictal abnormalities.

Epileptic High-Frequency Oscillations in Intracranial EEG Are Not Confounded by Cognitive Tasks

Rationale: High-frequency oscillations (HFOs) in intracranial EEG (iEEG) are used to delineate the epileptogenic zone during presurgical diagnostic assessment in patients with epilepsy. HFOs are historically divided into ripples (80-250 Hz), fast ripples (FR, >250 Hz), and their co-occurrence (FRandR). In a previous study, we had validated the rate of FRandRs during deep sleep to predict seizure outcome. Here, we ask whether epileptic FRandRs might be confounded by physiological FRandRs that are unrelated to epilepsy. Methods: We recorded iEEG in the medial temporal lobe MTL (hippocampus, entorhinal cortex, and amygdala) in 17 patients while they performed cognitive tasks. The three cognitive tasks addressed verbal working memory, visual working memory, and emotional processing. In our previous studies, these tasks activated the MTL. We re-analyzed the data of these studies with the automated detector that focuses on the co-occurrence of ripples and FRs (FRandR). Results: For each task, we identified those channels in which the HFO rate was modulated during the task condition compared to the control condition. However, the number of these channels did not exceed the chance level. Interestingly, even during wakefulness, the HFO rate was higher for channels within the seizure onset zone (SOZ) than for channels outside the SOZ. Conclusion: Our prospective definition of an epileptic HFO, the FRandR, is not confounded by physiological HFOs that might be elicited by our cognitive tasks. This is reassuring for the clinical use of FRandR as a biomarker of the EZ.