Physiological and pathological high frequency oscillations in focal epilepsy

Classifying High-Frequency Oscillations by Morphologic Contrast to Background, With Surgical Outcome Correlates

PURPOSE

Ictal high-frequency oscillations (HFOs) are a reliable indicator of a seizure onset zone for intracranial EEG recordings. Interictal HFOs often are also observed and may be a useful biomarker to supplement ictal data, but distinguishing pathologic from physiologic HFOs continues to be a challenging task. We present a method of classifying HFOs based on morphologic contrast to the background.

METHODS

We retrospectively screened 31 consecutive patients who underwent intracranial recordings for epilepsy at Stanford Medical Center during a 2-year period, and 13 patients met the criteria for inclusion. Interictal EEG data were analyzed using an automated event detector followed by morphologic feature extraction and k-means clustering. Instead of only using event features, the algorithm also incorporated features of the background adjacent to the events. High-frequency oscillations with higher morphologic contrast to the background were labeled as pathologic, and "hotspots" with the most active pathologic HFOs were identified and compared with clinically determined seizure onset zones.

RESULTS

Clustering with contrast features produced groups with better separation and more consistent boundaries. Eleven of the 13 patients proceeded to surgery, and patients whose hotspots matched seizure onset zones had better outcomes, with 4 out of 5 "match" patients having no disabling seizures at 1+ year postoperatively (Engel I or International League Against Epilepsy Class 1-2), while all "mismatch" patients continued to have disabling seizures (Fisher exact test P-value = 0.015).

CONCLUSIONS

High-frequency oscillations with higher contrast to background more likely represent paroxysmal bursts of pathologic activity. Patients with HFO hotspots outside of identified seizure onset zones may not respond as well to surgery.

Connectivity of high-frequency bursts as SOZ localization biomarker

For patients with refractory epilepsy, the seizure onset zone (SOZ) plays an essential role in determining the specific regions of the brain that will be surgically resected. High-frequency oscillations (HFOs) and connectivity-based approaches have been identified among the potential biomarkers to localize the SOZ. However, there is no consensus on how connectivity between HFO events should be estimated, nor on its subject-specific short-term reliability. Therefore, we propose the channel-level connectivity dispersion (CLCD) as a metric to quantify the variability in synchronization between individual electrodes and to identify clusters of electrodes with abnormal synchronization, which we hypothesize to be associated with the SOZ. In addition, we developed a specialized filtering method that reduces oscillatory components caused by filtering broadband artifacts, such as sharp transients, spikes, or direct current shifts. Our connectivity estimates are therefore robust to the presence of these waveforms. To calculate our metric, we start by creating binary signals indicating the presence of high-frequency bursts in each channel, from which we calculate the pairwise connectivity between channels. Then, the CLCD is calculated by combining the connectivity matrices and measuring the variability in each electrode's combined connectivity values. We test our method using two independent open-access datasets comprising intracranial electroencephalography signals from 89 to 15 patients with refractory epilepsy, respectively. Recordings in these datasets were sampled at approximately 1000 Hz, and our proposed CLCDs were estimated in the ripple band (80-200 Hz). Across all patients in the first dataset, the average ROC-AUC was 0.73, and the average Cohen's d was 1.05, while in the second dataset, the average ROC-AUC was 0.78 and Cohen's d was 1.07. On average, SOZ channels had lower CLCD values than non-SOZ channels. Furthermore, based on the second dataset, which includes surgical outcomes (Engel I-IV), our analysis suggested that higher CLCD interquartile (as a measure of CLCD distribution spread) is associated with favorable outcomes (Engel I). This suggests that CLCD could significantly assist in identifying SOZ clusters and, therefore, provide an additional tool in surgical planning for epilepsy patients.

A data augmentation procedure to improve detection of spike ripples in brain voltage recordings

Epilepsy is a major neurological disorder characterized by recurrent, spontaneous seizures. For patients with drug-resistant epilepsy, treatments include neurostimulation or surgical removal of the epileptogenic zone (EZ), the brain region responsible for seizure generation. Precise targeting of the EZ requires reliable biomarkers. Spike ripples - high-frequency oscillations that co-occur with large amplitude epileptic discharges - have gained prominence as a candidate biomarker. However, spike ripple detection remains a challenge. The gold-standard approach requires an expert manually visualize and interpret brain voltage recordings, which limits reproducibility and high-throughput analysis. Addressing these limitations requires more objective, efficient, and automated methods for spike ripple detection, including approaches that utilize deep neural networks. Despite advancements, dataset heterogeneity and scarcity severely limit machine learning performance. Our study explores long-short term memory (LSTM) neural network architectures for spike ripple detection, leveraging data augmentation to improve classifier performance. We highlight the potential of combining training on augmented and in vivo data for enhanced spike ripple detection and ultimately improving diagnostic accuracy in epilepsy treatment.

Practical measurements distinguishing physiological and pathological stereoelectroencephalography channels based on high-frequency oscillations in the human brain

OBJECTIVE

The present study aimed to identify various distinguishing features for use in the accurate classification of stereoelectroencephalography (SEEG) channels based on high-frequency oscillations (HFOs) inside and outside the epileptogenic zone (EZ).

METHODS

HFOs were detected in patients with focal epilepsy who underwent SEEG. Subsequently, HFOs within the seizure-onset and early spread zones were defined as pathological HFOs, whereas others were defined as physiological. Three features of HFOs were identified at the channel level, namely, morphological repetition, rhythmicity, and phase-amplitude coupling (PAC). A machine-learning (ML) classifier was then built to distinguish two HFO types at the channel level by application of the above-mentioned features, and the contributions were quantified. Further verification of the characteristics and classifier performance was performed in relation to various conscious states, imaging results, EZ location, and surgical outcomes.

RESULTS

Thirty-five patients were included in this study, from whom 166  104 pathological HFOs in 255 channels and 53 374 physiological HFOs in 282 channels were entered into the analysis pipeline. The results revealed that the morphological repetitions of pathological HFOs were markedly higher than those of the physiological HFOs; this was also observed for rhythmicity and PAC. The classifier exhibited high accuracy in differentiating between the two forms of HFOs, as indicated by an area under the curve (AUC) of 0.89. Both PAC and rhythmicity contributed significantly to this distinction. The subgroup analyses supported these findings.

SIGNIFICANCE

The suggested HFO features can accurately distinguish between pathological and physiological channels substantially improving its usefulness in clinical localization.

PLAIN LANGUAGE SUMMARY

In this study, we computed three quantitative features associated with HFOs in each SEEG channel and then constructed a machine learning-based classifier for the classification of pathological and physiological channels. The classifier performed well in distinguishing the two channel types under different levels of consciousness as well as in terms of imaging results, EZ location, and patient surgical outcomes.

WONOEP 2022: Neurotechnology for the diagnosis of epilepsy

Epilepsy's myriad causes and clinical presentations ensure that accurate diagnoses and targeted treatments remain a challenge. Advanced neurotechnologies are needed to better characterize individual patients across multiple modalities and analytical techniques. At the XVIth Workshop on Neurobiology of Epilepsy: Early Onset Epilepsies: Neurobiology and Novel Therapeutic Strategies (WONOEP 2022), the session on "advanced tools" highlighted a range of approaches, from molecular phenotyping of genetic epilepsy models and resected tissue samples to imaging-guided localization of epileptogenic tissue for surgical resection of focal malformations. These tools integrate cutting edge research, clinical data acquisition, and advanced computational methods to leverage the rich information contained within increasingly large datasets. A number of common challenges and opportunities emerged, including the need for multidisciplinary collaboration, multimodal integration, potential ethical challenges, and the multistage path to clinical translation. Despite these challenges, advanced epilepsy neurotechnologies offer the potential to improve our understanding of the underlying causes of epilepsy and our capacity to provide patient-specific treatment.

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.

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.

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.

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.

Weak coupling of neurons enables very high-frequency and ultra-fast oscillations through the interplay of synchronized phase shifts

Recently, in the past decade, high-frequency oscillations (HFOs), very high-frequency oscillations (VHFOs), and ultra-fast oscillations (UFOs) were reported in epileptic patients with drug-resistant epilepsy. However, to this day, the physiological origin of these events has yet to be understood. Our study establishes a mathematical framework based on bifurcation theory for investigating the occurrence of VHFOs and UFOs in depth EEG signals of patients with focal epilepsy, focusing on the potential role of reduced connection strength between neurons in an epileptic focus. We demonstrate that synchronization of a weakly coupled network can generate very and ultra high-frequency signals detectable by nearby microelectrodes. In particular, we show that a bistability region enables the persistence of phase-shift synchronized clusters of neurons. This phenomenon is observed for different hippocampal neuron models, including Morris-Lecar, Destexhe-Paré, and an interneuron model. The mechanism seems to be robust for small coupling, and it also persists with random noise affecting the external current. Our findings suggest that weakened neuronal connections could contribute to the production of oscillations with frequencies above 1000 Hz, which could advance our understanding of epilepsy pathology and potentially improve treatment strategies. However, further exploration of various coupling types and complex network models is needed.

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.

Differential patterns of very high-frequency oscillations in two seizure types of the pilocarpine-induced TLE model

AIMS

Very high-frequency oscillations (VHFOs, >500 Hz) are considered a highly sensitive biomarker of seizures. We hypothesized that VHFOs may exhibit specificity towards hypersynchronous (HYP) seizures and low-voltage fast (LVF) seizures in temporal lobe epilepsy (TLE).

METHODS

Local field potentials were recorded from the hippocampal network in TLE mice induced by pilocarpine. Subsequently, we analyzed the VHFO features, including their temporal-frequency characteristics and VHFO/theta coupling, during three states: baseline, preictal, and postictal for both HYP- and LVF-seizure groups.

RESULTS

Significant changes in most of the VHFO features were observed during the preictal state in both seizure groups. In the postictal state, VHFO features in the HYP-seizure group exhibited inverse alterations and appeared to align with those observed during baseline conditions. However, such phenomena were not observed after TLE seizures in the LVF-seizure group.

CONCLUSION

Our findings highlight distinct patterns of VHFO feature changes across different states of HYP seizures and LVF seizures. These results suggest that VHFOs could serve as indicative biomarkers for seizure alterations specifically associated with HYP-seizure states.

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.

Evolution of interictal activity in models of mesial temporal lobe epilepsy

Interictal activity and seizures are the hallmarks of focal epileptic disorders (which include mesial temporal lobe epilepsy, MTLE) in humans and in animal models. Interictal activity, which is recorded with cortical and intracerebral EEG recordings, comprises spikes, sharp waves and high-frequency oscillations, and has been used in clinical practice to identify the epileptic zone. However, its relation with seizures remains debated. Moreover, it is unclear whether specific EEG changes in interictal activity occur during the time preceding the appearance of spontaneous seizures. This period, which is termed "latent", has been studied in rodent models of MTLE in which spontaneous seizures start to occur following an initial insult (most often a status epilepticus induced by convulsive drugs such as kainic acid or pilocarpine) and may mirror epileptogenesis, i.e., the process leading the brain to develop an enduring predisposition to seizure generation. Here, we will address this topic by reviewing experimental studies performed in MTLE models. Specifically, we will review data highlighting the dynamic changes in interictal spiking activity and high-frequency oscillations occurring during the latent period, and how optogenetic stimulation of specific cell populations can modulate them in the pilocarpine model. These findings indicate that interictal activity: (i) is heterogeneous in its EEG patterns and thus, presumably, in its underlying neuronal mechanisms; and (ii) can pinpoint to the epileptogenic processes occurring in focal epileptic disorders in animal models and, perhaps, in epileptic patients.

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.

Intracranial electrophysiological recordings on a swine model of mesial temporal lobe epilepsy

Objective

To test the feasibility and reliability of intracranial electrophysiological recordings in an acute status epilepticus model on laboratory swine.

Method

Intrahippocampal injection of kainic acid (KA) was performed on 17 male Bama pigs (Sus scrofa domestica) weighing between 25 and 35 kg. Two stereoelectroencephalography (SEEG) electrodes with a total of 16 channels were implanted bilaterally along the sensorimotor cortex to the hippocampus. Brain electrical activity was recorded 2 h daily for 9-28 days. Three KA dosages were tested to evaluate the quantities capable of evoking status epilepticus. Local field potentials (LFPs) were recorded and compared before and after the KA injection. We quantified the epileptic patterns, including the interictal spikes, seizures, and high-frequency oscillations (HFOs), up to 4 weeks after the KA injection. Test-retest reliability using intraclass correlation coefficients (ICCs) were performed on interictal HFO rates to evaluate the recording stability of this model.

Results

The KA dosage test suggested that a 10 μl (1.0 μg/μl) intrahippocampal injection could successfully evoke status epilepticus lasting from 4 to 12 h. At this dosage, eight pigs (50% of total) had prolonged epileptic events (tonic-chronic seizures + interictal spikes n = 5, interictal spikes alone n = 3) in the later 4 weeks of the video-SEEG recording period. Four pigs (25% of total) had no epileptic activities, and another four (25%) had lost the cap or did not complete the experiments. Animals that showed epileptiform events were grouped as E + (n = 8) and the four animals showing no signs of epileptic events were grouped as E- (n = 4). A total of 46 electrophysiological seizures were captured in the 4-week post-KA period from 4 E + animals, with the earliest onset on day 9. The seizure durations ranged from 12 to 45 s. A significant increase of hippocampal HFOs rate (num/min) was observed in the E+ group during the post-KA period (weeks 1, 2,4, p < 0.05) compared to the baseline. But the E-showed no change or a decrease (in week 2, p = 0.43) compared to their baseline rate. The between-group comparison showed much higher HFO rates in E + vs. E - (F = 35, p < 0.01). The high ICC value [ICC (1, k) = 0.81, p < 0.05] quantified from the HFO rate suggested that this model had a stable measurement of HFOs during the four-week post-KA periods.

Significance

This study measured intracranial electrophysiological activity in a swine model of KA-induced mesial temporal lobe epilepsy (mTLE). Using the clinical SEEG electrode, we distinguished abnormal EEG patterns in the swine brain. The high test-retest reliability of HFO rates in the post-KA period suggests the utility of this model for studying mechanisms of epileptogenesis. The use of swine may provide satisfactory translational value for clinical epilepsy research.

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

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

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.

A Systematic Review of Time Series Classification Techniques Used in Biomedical Applications

Background: Digital clinical measures collected via various digital sensing technologies such as smartphones, smartwatches, wearables, and ingestible and implantable sensors are increasingly used by individuals and clinicians to capture the health outcomes or behavioral and physiological characteristics of individuals. Time series classification (TSC) is very commonly used for modeling digital clinical measures. While deep learning models for TSC are very common and powerful, there exist some fundamental challenges. This review presents the non-deep learning models that are commonly used for time series classification in biomedical applications that can achieve high performance. Objective: We performed a systematic review to characterize the techniques that are used in time series classification of digital clinical measures throughout all the stages of data processing and model building. Methods: We conducted a literature search on PubMed, as well as the Institute of Electrical and Electronics Engineers (IEEE), Web of Science, and SCOPUS databases using a range of search terms to retrieve peer-reviewed articles that report on the academic research about digital clinical measures from a five-year period between June 2016 and June 2021. We identified and categorized the research studies based on the types of classification algorithms and sensor input types. Results: We found 452 papers in total from four different databases: PubMed, IEEE, Web of Science Database, and SCOPUS. After removing duplicates and irrelevant papers, 135 articles remained for detailed review and data extraction. Among these, engineered features using time series methods that were subsequently fed into widely used machine learning classifiers were the most commonly used technique, and also most frequently achieved the best performance metrics (77 out of 135 articles). Statistical modeling (24 out of 135 articles) algorithms were the second most common and also the second-best classification technique. Conclusions: In this review paper, summaries of the time series classification models and interpretation methods for biomedical applications are summarized and categorized. While high time series classification performance has been achieved in digital clinical, physiological, or biomedical measures, no standard benchmark datasets, modeling methods, or reporting methodology exist. There is no single widely used method for time series model development or feature interpretation, however many different methods have proven successful.

Epileptogenic high-frequency oscillations present larger amplitude both in mesial temporal and neocortical regions

High-frequency oscillations (HFO) are a promising biomarker for the identification of epileptogenic tissue. While HFO rates have been shown to predict seizure outcome, it is not yet clear whether their morphological features might improve this prediction. We validated HFO rates against seizure outcome and delineated the distribution of HFO morphological features. We collected stereo-EEG recordings from 20 patients (231 electrodes; 1,943 contacts). We computed HFO rates (the co-occurrence of ripples and fast ripples) through a validated automated detector during non-rapid eye movement sleep. Applying machine learning, we delineated HFO morphological features within and outside epileptogenic tissue across mesial temporal lobe (MTL) and Neocortex. HFO rates predicted seizure outcome with 85% accuracy, 79% specificity, 100% sensitivity, 100% negative predictive value, and 67% positive predictive value. The analysis of HFO features showed larger amplitude in the epileptogenic tissue, similar morphology for epileptogenic HFO in MTL and Neocortex, and larger amplitude for physiological HFO in MTL. We confirmed HFO rates as a reliable biomarker for epilepsy surgery and characterized the potential clinical relevance of HFO morphological features. Our results support the prospective use of HFO in epilepsy surgery and contribute to the anatomical mapping of HFO morphology.

A sparse representation strategy to eliminate pseudo-HFO events from intracranial EEG for seizure onset zone localization

Objective.High-frequency oscillations (HFOs) are considered a biomarker of the epileptogenic zone in intracranial EEG recordings. However, automated HFO detectors confound true oscillations with spurious events caused by the presence of artifacts.Approach.We hypothesized that, unlike pseudo-HFOs with sharp transients or arbitrary shapes, real HFOs have a signal characteristic that can be represented using a small number of oscillatory bases. Based on this hypothesis using a sparse representation framework, this study introduces a new classification approach to distinguish true HFOs from the pseudo-events that mislead seizure onset zone (SOZ) localization. Moreover, we further classified the HFOs into ripples and fast ripples by introducing an adaptive reconstruction scheme using sparse representation. By visualizing the raw waveforms and time-frequency representation of events recorded from 16 patients, three experts labeled 6400 candidate events that passed an initial amplitude-threshold-based HFO detector. We formed a redundant analytical multiscale dictionary built from smooth oscillatory Gabor atoms and represented each event with orthogonal matching pursuit by using a small number of dictionary elements. We used the approximation error and residual signal at each iteration to extract features that can distinguish the HFOs from any type of artifact regardless of their corresponding source. We validated our model on sixteen subjects with thirty minutes of continuous interictal intracranial EEG recording from each.Main results.We showed that the accuracy of SOZ detection after applying our method was significantly improved. In particular, we achieved a 96.65% classification accuracy in labeled events and a 17.57% improvement in SOZ detection on continuous data. Our sparse representation framework can also distinguish between ripples and fast ripples.Significance.We show that by using a sparse representation approach we can remove the pseudo-HFOs from the pool of events and improve the reliability of detected HFOs in large data sets and minimize manual artifact elimination.

Phase-amplitude coupling between low-frequency scalp EEG and high-frequency intracranial EEG during working memory task

Objective. Revealing the relationship between simultaneous scalp electroencephalography (EEG) and intracranial electroencephalography (iEEG) is of great importance for both neuroscientific research and translational applications. However, whether prominent iEEG features in the high-gamma band can be reflected by scalp EEG is largely unknown. To address this, we investigated the phase-amplitude coupling (PAC) phenomenon between the low-frequency band of scalp EEG and the high-gamma band of iEEG.Approach. We analyzed a simultaneous iEEG and scalp EEG dataset acquired under a verbal working memory paradigm from nine epilepsy subjects. The PAC values between pairs of scalp EEG channel and identified iEEG channel were explored. After identifying the frequency combinations and electrode locations that generated the most significant PAC values, we compared the PAC values of different task periods (encoding, maintenance, and retrieval) and memory loads.Main results. We demonstrated that the amplitude of high-gamma activities in the entorhinal cortex, hippocampus, and amygdala was correlated to the delta or theta phase at scalp locations such as Cz and Pz. In particular, the frequency bin that generated the maximum PAC value centered at 3.16-3.84 Hz for the phase and 50-85 Hz for the amplitude. Moreover, our results showed that PAC values for the retrieval period were significantly higher than those of the encoding and maintenance periods, and the PAC was also influenced by the memory load.Significance. This is the first human simultaneous iEEG and scalp EEG study demonstrating that the amplitude of iEEG high-gamma components is associated with the phase of low-frequency components in scalp EEG. These findings enhance our understanding of multiscale neural interactions during working memory, and meanwhile, provide a new perspective to estimate intracranial high-frequency features with non-invasive neural recordings.

Are High Frequency Oscillations in Scalp EEG Related to Age?

BACKGROUND

High-frequency oscillations (HFOs) have received much attention in recent years, particularly in the clinical context. In addition to their application as a marker for pathological changes in patients with epilepsy, HFOs have also been brought into context with several physiological mechanisms. Furthermore, recent studies reported a relation between an increase of HFO rate and age in invasive EEG recordings. The present study aimed to investigate whether this relation can be replicated in scalp-EEG.

METHODS

We recorded high-density EEG from 11 epilepsy patients at rest as well as during motor performance. Manual detection of HFOs was performed by two independent raters following a standardized protocol. Patients were grouped by age into younger (<25 years) and older (>50 years) participants.

RESULTS

No significant difference of HFO-rates was found between groups [U = 10.5, p = 0.429, r = 0.3].

CONCLUSIONS

Lack of replicability of the age effect of HFOs may be due to the local propagation patterns of age-related HFOs occurring in deep structures. However, limitations such as small sample size, decreased signal-to-noise ratio as compared to invasive recordings, as well as HFO-mimicking artifacts must be considered.

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.

Deep learning for epileptogenic zone delineation from the invasive EEG: challenges and lookouts

This scientific commentary refers to 'Refining epileptogenic high-frequency oscillations using deep learning: a reverse engineering approach' by Zhang et al. (https://doi.org/10.1093/braincomms/fcab267).

Correcting for physiological ripples improves epileptic focus identification and outcome prediction

Objective

The integration of high‐frequency oscillations (HFOs; ripples [80–250 Hz], fast ripples [250–500 Hz]) in epilepsy evaluation is hampered by physiological HFOs, which cannot be reliably differentiated from pathological HFOs. We evaluated whether defining abnormal HFO rates by statistical comparison to region‐specific physiological HFO rates observed in the healthy brain improves identification of the epileptic focus and surgical outcome prediction.

Methods

We detected HFOs in 151 consecutive patients who underwent stereo‐electroencephalography and subsequent resective epilepsy surgery at two tertiary epilepsy centers. We compared how HFOs identified the resection cavity and predicted seizure‐free outcome using two thresholds from the literature (HFO rate > 1/min; 50% of the total number of a patient's HFOs) and three thresholds based on normative rates from the Montreal Neurological Institute Open iEEG Atlas (https://mni‐open‐ieegatlas.research.mcgill.ca/): global Atlas threshold, regional Atlas threshold, and regional + 10% threshold after regional Atlas correction.

Results

Using ripples, the regional + 10% threshold performed best for focus identification (77.3% accuracy, 27% sensitivity, 97.1% specificity, 80.6% positive predictive value [PPV], 78.2% negative predictive value [NPV]) and outcome prediction (69.5% accuracy, 58.6% sensitivity, 76.3% specificity, 60.7% PPV, 74.7% NPV). This was an improvement for focus identification (+1.1% accuracy, +17.0% PPV; p < .001) and outcome prediction (+12.0% sensitivity, +1.0% PPV; p = .05) compared to the 50% threshold. The improvement was particularly marked for foci in cortex, where physiological ripples are frequent (outcome: +35.3% sensitivity, +5.3% PPV; p = .014). In these cases, the regional + 10% threshold outperformed fast ripple rate > 1/min (+3.6% accuracy, +26.5% sensitivity, +21.6% PPV; p < .001) and seizure onset zone (+13.5% accuracy, +29.4% sensitivity, +17.0% PPV; p < .05–.01) for outcome prediction. Normalization did not improve the performance of fast ripples.

Significance

Defining abnormal HFO rates by statistical comparison to rates in healthy tissue overcomes an important weakness in the clinical use of ripples. It improves focus identification and outcome prediction compared to standard HFO measures, increasing their clinical applicability.

Refining epileptogenic high-frequency oscillations using deep learning: a reverse engineering approach

Intracranially recorded interictal high-frequency oscillations have been proposed as a promising spatial biomarker of the epileptogenic zone. However, its visual verification is time-consuming and exhibits poor inter-rater reliability. Furthermore, no method is currently available to distinguish high-frequency oscillations generated from the epileptogenic zone (epileptogenic high-frequency oscillations) from those generated from other areas (non-epileptogenic high-frequency oscillations). To address these issues, we constructed a deep learning-based algorithm using chronic intracranial EEG data via subdural grids from 19 children with medication-resistant neocortical epilepsy to: (i) replicate human expert annotation of artefacts and high-frequency oscillations with or without spikes, and (ii) discover epileptogenic high-frequency oscillations by designing a novel weakly supervised model. The ‘purification power’ of deep learning is then used to automatically relabel the high-frequency oscillations to distill epileptogenic high-frequency oscillations. Using 12 958 annotated high-frequency oscillation events from 19 patients, the model achieved 96.3% accuracy on artefact detection (F1 score = 96.8%) and 86.5% accuracy on classifying high-frequency oscillations with or without spikes (F1 score = 80.8%) using patient-wise cross-validation. Based on the algorithm trained from 84 602 high-frequency oscillation events from nine patients who achieved seizure-freedom after resection, the majority of such discovered epileptogenic high-frequency oscillations were found to be ones with spikes (78.6%, P < 0.001). While the resection ratio of detected high-frequency oscillations (number of resected events/number of detected events) did not correlate significantly with post-operative seizure freedom (the area under the curve = 0.76, P = 0.06), the resection ratio of epileptogenic high-frequency oscillations positively correlated with post-operative seizure freedom (the area under the curve = 0.87, P = 0.01). We discovered that epileptogenic high-frequency oscillations had a higher signal intensity associated with ripple (80–250 Hz) and fast ripple (250–500 Hz) bands at the high-frequency oscillation onset and with a lower frequency band throughout the event time window (the inverted T-shaped), compared to non-epileptogenic high-frequency oscillations. We then designed perturbations on the input of the trained model for non-epileptogenic high-frequency oscillations to determine the model’s decision-making logic. The model confidence significantly increased towards epileptogenic high-frequency oscillations by the artificial introduction of the inverted T-shaped signal template (mean probability increase: 0.285, P < 0.001), and by the artificial insertion of spike-like signals into the time domain (mean probability increase: 0.452, P < 0.001). With this deep learning-based framework, we reliably replicated high-frequency oscillation classification tasks by human experts. Using a reverse engineering technique, we distinguished epileptogenic high-frequency oscillations from others and identified its salient features that aligned with current knowledge.

Identification of Epileptogenic and Non-epileptogenic High-Frequency Oscillations Using a Multi-Feature Convolutional Neural Network Model

Accurately identifying epileptogenic zone (EZ) using high-frequency oscillations (HFOs) is a challenge that must be mastered to transfer HFOs into clinical use. We analyzed the ability of a convolutional neural network (CNN) model to distinguish EZ and non-EZ HFOs. Nineteen medically intractable epilepsy patients with good surgical outcomes 2 years after surgery were studied. Five-minute interictal intracranial electroencephalogram epochs of slow-wave sleep were selected randomly. Then 5 s segments of ripples (80–200 Hz) and fast ripples (FRs, 200–500 Hz) were detected automatically. The EZs and non-EZs were identified using the surgery resection range. We innovatively converted all epochs into four types of images using two scales: original waveforms, filtered waveforms, wavelet spectrum images, and smoothed pseudo Wigner–Ville distribution (SPWVD) spectrum images. Two scales were fixed and fitted scales. We then used a CNN model to classify the HFOs into EZ and non-EZ categories. As a result, 7,000 epochs of ripples and 2,000 epochs of FRs were randomly selected from the EZ and non-EZ data for analysis. Our CNN model can distinguish EZ and non-EZ HFOs successfully. Except for original ripple waveforms, the results from CNN models that are trained using fixed-scale images are significantly better than those from models trained using fitted-scale images (p < 0.05). Of the four fixed-scale transformations, the CNN based on the adjusted SPWVD (ASPWVD) produced the best accuracies (80.89 ± 1.43% and 77.85 ± 1.61% for ripples and FRs, respectively, p < 0.05). The CNN using ASPWVD transformation images is an effective deep learning method that can be used to classify EZ and non-EZ HFOs.

Low frequency novel interictal EEG biomarker for localizing seizures and predicting outcomes

Localizing hyperexcitable brain tissue to treat focal seizures remains challenging. We want to identify the seizure onset zone from interictal EEG biomarkers. We hypothesize that a combination of interictal EEG biomarkers, including a novel low frequency marker, can predict mesial temporal involvement and can assist in prognosis related to surgical resections. Interictal direct current wide bandwidth invasive EEG recordings from 83 patients implanted with 5111 electrodes were retrospectively studied. Logistic regression was used to classify electrodes and patient outcomes. A feed-forward neural network was implemented to understand putative mechanisms. Interictal infraslow frequency EEG activity was decreased for seizure onset zone electrodes while faster frequencies such as delta (2-4 Hz) and beta-gamma (20-50 Hz) activity were increased. These spectral changes comprised a novel interictal EEG biomarker that was significantly increased for mesial temporal seizure onset zone electrodes compared to non-seizure onset zone electrodes. Interictal EEG biomarkers correctly classified mesial temporal seizure onset zone electrodes with a specificity of 87% and positive predictive value of 80%. These interictal EEG biomarkers also correctly classified patient outcomes after surgical resection with a specificity of 91% and positive predictive value of 87%. Interictal infraslow EEG activity is decreased near the seizure onset zone while higher frequency power is increased, which may suggest distinct underlying physiologic mechanisms. Narrowband interictal EEG power bands provide information about the seizure onset zone and can help predict mesial temporal involvement in seizure onset. Narrowband interictal EEG power bands may be less useful for predictions related to non-mesial temporal electrodes. Together with interictal epileptiform discharges and high-frequency oscillations, these interictal biomarkers may provide prognostic information prior to surgical resection. Computational modelling suggests changes in neural adaptation may be related to the observed low frequency power changes.

A study on the relationship between non-epileptic fast (40 - 200 Hz) oscillations in scalp EEG and development in children

OBJECTIVE

Physiological gamma and ripple activities may be linked to neurocognitive functions. This study investigated the relationship between development and non-epileptic, probably physiological, fast (40-200 Hz) oscillations (FOs) including gamma (40 - 80 Hz) and ripple (80 - 200 Hz) oscillations in scalp EEG in children with neurodevelopmental disorders.

METHODS

Participants were 124 children with autism spectrum disorder (ASD) and/or attention deficit/hyperactivity disorder (ADHD). Gamma and ripple oscillations were explored from 60-second-long sleep EEG data in each subject using a semi-automatic detection tool supplemented with visual confirmation and time-frequency analysis.

RESULTS

Gamma and ripple oscillations were detected in 25 (20.2%) and 22 (17.7%) children, respectively. The observation of one or more occurrence(s) of ripple oscillations, but not gamma oscillations, was significantly related to lower age at EEG recording (odds ratio, OR: 0.727 [95% confidence interval, CI: 0.568-0.929]), higher intelligence/developmental quotient (OR: 1.041, 95% CI: 1.002-1.082), and lack of a diagnosis with ADHD (OR: 0.191, 95% CI: 0.039 - 0.937) according to a binominal logistic regression analysis that included diagnosis with ASD, sex, history of perinatal complications, history of febrile seizures, and use of a sedative agent for the EEG recording as the other non-significant parameters. Diagnostic group was not related to frequency or power of spectral peaks of FOs.

CONCLUSION

The production of non-epileptic scalp ripples was confirmed to be associated with brain development and function/dysfunction in childhood. Further investigation is necessary to interpret all of the information on higher brain functions that may be embedded in scalp FOs.

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.

Leveraging electrophysiologic correlates of word encoding to map seizure onset zone in focal epilepsy: Task-dependent changes in epileptiform activity, spectral features, and functional connectivity

OBJECTIVE

Verbal memory dysfunction is common in focal, drug-resistant epilepsy (DRE). Unfortunately, surgical removal of seizure-generating brain tissue can be associated with further memory decline. Therefore, localization of both the circuits generating seizures and those underlying cognitive functions is critical in presurgical evaluations for patients who may be candidates for resective surgery. We used intracranial electroencephalographic (iEEG) recordings during a verbal memory task to investigate word encoding in focal epilepsy. We hypothesized that engagement in a memory task would exaggerate local iEEG feature differences between the seizure onset zone (SOZ) and neighboring tissue as compared to wakeful rest ("nontask").

METHODS

Ten participants undergoing presurgical iEEG evaluation for DRE performed a free recall verbal memory task. We evaluated three iEEG features in SOZ and non-SOZ electrodes during successful word encoding and compared them with nontask recordings: interictal epileptiform spike (IES) rates, power in band (PIB), and relative entropy (REN; a functional connectivity measure).

RESULTS

We found a complex pattern of PIB and REN changes in SOZ and non-SOZ electrodes during successful word encoding compared to nontask. Successful word encoding was associated with a reduction in local electrographic functional connectivity (increased REN), which was most exaggerated in temporal lobe SOZ. The IES rates were reduced during task, but only in the non-SOZ electrodes. Compared with nontask, REN features during task yielded marginal improvements in SOZ classification.

SIGNIFICANCE

Previous studies have supported REN as a biomarker for epileptic brain. We show that REN differences between SOZ and non-SOZ are enhanced during a verbal memory task. We also show that IESs are reduced during task in non-SOZ, but not in SOZ. These findings support the hypothesis that SOZ and non-SOZ respond differently to task and warrant further exploration into the use of cognitive tasks to identify functioning memory circuits and localize SOZ.

HFO Detection in Epilepsy: A Stacked Denoising Autoencoder and Sample Weight Adjusting Factors-Based Method

High-frequency oscillations (HFOs) recorded by the intracranial electroencephalography (iEEG) are the promising biomarkers of epileptogenic zones. Accurate detection of HFOs is the key to pre-operative assessment for epilepsy. Due to the subjective bias caused by manual features and the class imbalance between HFOs and false HFOs, it is difficult to obtain satisfactory detection performance by the existing methods. To solve these problems, we put forward a novel method to accurately detect HFOs based on the stacked denoising autoencoder (SDAE) and the ensemble classifier with sample weight adjusting factors. First, the adjustable threshold of Hilbert envelopes is proposed to isolate the events of interest (EoIs) from background activities. Then, the SDAE network is utilized to automatically extract features of EoIs in the time-frequency domain. Finally, the AdaBoost-based support vector machine ensemble classifier with sample weight adjusting factors is devised to separate HFOs from EoIs by using the extracted features. These adjusting factors are used to solve the class imbalance problem by adjusting sample weights when learning the base classifiers. Our HFO detection method is evaluated by using clinical iEEG data recorded from 20 patients with medically refractory epilepsy. The experimental results show that our detection method outperforms some existing methods in terms of sensitivity and false discovery rate. In addition, the HFOs detected by our method are effective for localizing seizure onset zones.

Phase-amplitude coupling between infraslow and high-frequency activities well discriminates between the preictal and interictal states

Infraslow activity (ISA) and high-frequency activity (HFA) are key biomarkers for studying epileptic seizures. We aimed to elucidate the relationship between ISA and HFA around seizure onset. We enrolled seven patients with drug-resistant focal epilepsy who underwent intracranial electrode placement. We comparatively analyzed the ISA, HFA, and ISA-HFA phase-amplitude coupling (PAC) in the seizure onset zone (SOZ) or non-SOZ (nSOZ) in the interictal, preictal, and ictal states. We recorded 15 seizures. HFA and ISA were larger in the ictal states than in the interictal or preictal state. During seizures, the HFA and ISA of the SOZ were larger and occurred earlier than those of nSOZ. In the preictal state, the ISA-HFA PAC of the SOZ was larger than that of the interictal state, and it began increasing at approximately 87 s before the seizure onset. The receiver-operating characteristic curve revealed that the ISA-HFA PAC of the SOZ showed the highest discrimination performance in the preictal and interictal states, with an area under the curve of 0.926. This study demonstrated the novel insight that ISA-HFA PAC increases before the onset of seizures. Our findings indicate that ISA-HFA PAC could be a useful biomarker for discriminating between the preictal and interictal states.

Spontaneous modulations of high-frequency cortical activity

OBJECTIVE

We clarified the clinical and mechanistic significance of physiological modulations of high-frequency broadband cortical activity associated with spontaneous saccadic eye movements during a resting state.

METHODS

We studied 30 patients who underwent epilepsy surgery following extraoperative electrocorticography and electrooculography recordings. We determined whether high-gamma activity at 70-110 Hz preceding saccade onset would predict upcoming ocular behaviors. We assessed how accurately the model incorporating saccade-related high-gamma modulations would localize the primary visual cortex defined by electrical stimulation.

RESULTS

The dynamic atlas demonstrated transient high-gamma suppression in the striatal cortex before saccade onset and high-gamma augmentation subsequently involving the widespread posterior brain regions. More intense striatal high-gamma suppression predicted the upcoming saccade directed to the ipsilateral side and lasting longer in duration. The bagged-tree-ensemble model demonstrated that intense saccade-related high-gamma modulations localized the visual cortex with an accuracy of 95%.

CONCLUSIONS

We successfully animated the neural dynamics supporting saccadic suppression, a principal mechanism minimizing the perception of blurred vision during rapid eye movements. The primary visual cortex per se may prepare actively in advance for massive image motion expected during upcoming prolonged saccades.

SIGNIFICANCE

Measuring saccade-related electrocorticographic signals may help localize the visual cortex and avoid misperceiving physiological high-frequency activity as epileptogenic.

Microscale dynamics of electrophysiological markers of epilepsy

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

SIGNIFICANCE

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

Presurgical Evaluation of Epilepsy Using Resting-State MEG Functional Connectivity

Epilepsy is caused by abnormal electrical discharges (clinically identified by electrophysiological recording) in a specific part of the brain [originating in only one part of the brain, namely, the epileptogenic zone (EZ)]. Epilepsy is now defined as an archetypical hyperexcited neural network disorder. It can be investigated through the network analysis of interictal discharges, ictal discharges, and resting-state functional connectivity. Currently, there is an increasing interest in embedding resting-state connectivity analysis into the preoperative evaluation of epilepsy. Among the various neuroimaging technologies employed to achieve brain functional networks, magnetoencephalography (MEG) with the excellent temporal resolution is an ideal tool for estimating the resting-state connectivity between brain regions, which can reveal network abnormalities in epilepsy. What value does MEG resting-state functional connectivity offer for epileptic presurgical evaluation? Regarding this topic, this paper introduced the origin of MEG and the workflow of constructing source-space functional connectivity based on MEG signals. Resting-state functional connectivity abnormalities correlate with epileptogenic networks, which are defined by the brain regions involved in the production and propagation of epileptic activities. This paper reviewed the evidence of altered epileptic connectivity based on low- or high-frequency oscillations (HFOs) and the evidence of the advantage of using simultaneous MEG and intracranial electroencephalography (iEEG) recordings. More importantly, this review highlighted that MEG-based resting-state functional connectivity has the potential to predict postsurgical outcomes. In conclusion, resting-state MEG functional connectivity has made a substantial progress toward serving as a candidate biomarker included in epileptic presurgical evaluations.

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.

Neuromagnetic high frequency spikes are a new and noninvasive biomarker for localization of epileptogenic zones

OBJECTIVE

One barrier hindering high frequency brain signals (HFBS, >80 Hz) from wide clinical applications is that the brain generates both pathological and physiological HFBS. This study was to find specific biomarkers for localizing epileptogenic zones (EZs).

METHODS

Twenty three children with drug-resistant epilepsy and age/sex matched healthy controls were studied with magnetoencephalography (MEG). High frequency oscillations (HFOs, > 4 oscillatory waveforms) and high frequency spikes (HFSs, > 1 spiky or sharp waveforms) in 80-250 Hz and 250-600 Hz bands were blindly detected with an artificial intelligence method and validated with visual inspection. The magnitude of HFOs and HFSs were quantified with spectral analyses. Sources of HFSs and HFOs were localized and compared with clinical EZs determined by invasive recordings and surgical outcomes.

RESULTS

HFOs in 80-250 Hz and 250-600 Hz were identified in both epilepsy patients (18/23, 12/23, respectively) and healthy controls (6/23, 4/23, respectively). HFSs in 80-250 Hz and 250-600 Hz were detected in patients (16/23, 11/23, respectively) but not in healthy controls. A combination of HFOs and HFSs localized EZs for 22 (22/23, 96%) patients.

CONCLUSIONS

The results indicate, for the first time, that HFSs are a newer and more specific biomarker than HFOs for localizing EZs because HFOs appeared in both epilepsy patients and healthy controls while HFSs appeared only in epilepsy patients.

Noninvasive Mapping of Ripple Onset Predicts Outcome in Epilepsy Surgery

OBJECTIVE

Intracranial electroencephalographic (icEEG) studies show that interictal ripples propagate across the brain of children with medically refractory epilepsy (MRE), and the onset of this propagation (ripple onset zone [ROZ]) estimates the epileptogenic zone. It is still unknown whether we can map this propagation noninvasively. The goal of this study is to map ripples (ripple zone [RZ]) and their propagation onset (ROZ) using high-density EEG (HD-EEG) and magnetoencephalography (MEG), and to estimate their prognostic value in pediatric epilepsy surgery.

METHODS

We retrospectively analyzed simultaneous HD-EEG and MEG data from 28 children with MRE who underwent icEEG and epilepsy surgery. Using electric and magnetic source imaging, we estimated virtual sensors (VSs) at brain locations that matched the icEEG implantation. We detected ripples on VSs, defined the virtual RZ and virtual ROZ, and estimated their distance from icEEG. We assessed the predictive value of resecting virtual RZ and virtual ROZ for postsurgical outcome. Interictal spike localization on HD-EEG and MEG was also performed and compared with ripples.

RESULTS

We mapped ripple propagation in all patients with HD-EEG and in 27 (96%) patients with MEG. The distance from icEEG did not differ between HD-EEG and MEG when mapping the RZ (26-27mm, p = 0.6) or ROZ (22-24mm, p = 0.4). Resecting the virtual ROZ, but not virtual RZ or the sources of spikes, was associated with good outcome for HD-EEG (p = 0.016) and MEG (p = 0.047).

INTERPRETATION

HD-EEG and MEG can map interictal ripples and their propagation onset (virtual ROZ). Noninvasively mapping the ripple onset may augment epilepsy surgery planning and improve surgical outcome of children with MRE. ANN NEUROL 2021;89:911-925.

Distinction of Physiologic and Epileptic Ripples: An Electrical Stimulation Study

Ripple oscillations (80-250 Hz) are a promising biomarker of epileptic activity, but are also involved in memory consolidation, which impairs their value as a diagnostic tool. Distinguishing physiologic from epileptic ripples has been particularly challenging because usually, invasive recordings are only performed in patients with refractory epilepsy. Here, we identified 'healthy' brain areas based on electrical stimulation and hypothesized that these regions specifically generate 'pure' ripples not coupled to spikes. Intracranial electroencephalography (EEG) recorded with subdural grid electrodes was retrospectively analyzed in 19 patients with drug-resistant focal epilepsy. Interictal spikes and ripples were automatically detected in slow-wave sleep using the publicly available Delphos software. We found that rates of spikes, ripples and ripples coupled to spikes ('spike-ripples') were higher inside the seizure-onset zone (p < 0.001). A comparison of receiver operating characteristic curves revealed that spike-ripples slightly delineated the seizure-onset zone channels, but did this significantly better than spikes (p < 0.001). Ripples were more frequent in the eloquent neocortex than in the remaining non-seizure onset zone areas (p < 0.001). This was due to the higher rates of 'pure' ripples (p < 0.001; median rates 3.3/min vs. 1.4/min), whereas spike-ripple rates were not significantly different (p = 0.87). 'Pure' ripples identified 'healthy' channels significantly better than chance (p < 0.001). Our findings suggest that, in contrast to epileptic spike-ripples, 'pure' ripples are mainly physiological. They may be considered, in addition to electrical stimulation, to delineate eloquent cortex in pre-surgical patients. Since we applied open source software for detection, our approach may be generally suited to tackle a variety of research questions in epilepsy and cognitive science.

Delineation of epileptogenic zones with high frequency magnetic source imaging based on kurtosis and skewness

BACKGROUND

Neuromagnetic high frequency brain signals (HFBS, > 80 Hz) are a new biomarker for localization of epileptogenic zones (EZs) for pediatric epilepsy.

METHODS

Twenty three children with drug-resistant epilepsy and age/sex matched healthy controls were studied with magnetoencephalography (MEG). Epileptic HFBS in 80-250 Hz and 250-600 Hz were quantitatively determined by comparing with normative controls in terms of kurtosis and skewness. Magnetic sources of epileptic HFBS were localized and then compared to clinical EZs determined by invasive recordings and surgical outcomes.

RESULTS

Kurtosis and skewness of HFBS were significantly elevated in epilepsy patients compared to healthy controls (p < 0,001 and p < 0.0001, respectively). Sources of elevated MEG signals in comparison to normative data were co-localized to EZs for 22 (22/23, 96 %) patients.

CONCLUSIONS

The results indicate, for the first time, that epileptic HFBS can be noninvasively quantified by measuring kurtosis and skewness in MEG data. Magnetic source imaging based on kurtosis and skewness can accurately localize EZs.

SIGNIFICANCE

Source imaging of kurtosis and skewness of MEG HFBS provides a novel way for preoperative localization of EZs for epilepsy surgery.

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.

High-Frequency Oscillations in Epilepsy: What Have We Learned and What Needs to be Addressed

For the past 2 decades, high-frequency oscillations (HFOs) have been enthusiastically studied by the epilepsy community. Emerging evidence shows that HFOs harbor great promise to delineate epileptogenic brain areas and possibly predict the likelihood of seizures. Investigations into HFOs in clinical epilepsy have advanced from small retrospective studies relying on visual identification and correlation analysis to larger prospective assessments using automatic detection and prediction strategies. Although most studies have yielded promising results, some have revealed significant obstacles to clinical application of HFOs, thus raising debate about the reliability and practicality of HFOs as clinical biomarkers. In this review, we give an overview of the current state of HFO research and pinpoint the conceptual and methodological issues that have hampered HFO translation. We highlight recent insights gained from long-term data, high-density recordings, and multicenter collaborations and discuss the open questions that need to be addressed in future research.

Scalp Ripples Can Predict Development of Epilepsy After First Unprovoked Seizure in Childhood

OBJECTIVE

Identification of children at risk of developing epilepsy after a first unprovoked seizure can be challenging. Interictal epileptiform discharges are associated with higher risk but have limited sensitivity and specificity. High frequency oscillations (HFOs) are newer biomarkers for epileptogenesis. We prospectively evaluated the predictive value of HFOs for developing epilepsy in scalp electroencephalogram (EEG) of children after a first unprovoked seizure.

METHODS

After their first seizure, 56 children were followed prospectively over 12 months and then grouped in "epilepsy" or "no epilepsy." Initial EEGs were visually analyzed for spikes, spike ripples, and ripples. Inter-group comparisons of spike-rates and HFO-rates were done by Mann-Whitney U test. Predictive values and optimal thresholds were calculated by receiver operating characteristic (ROC) curves.

RESULTS

In the epilepsy group (n = 26, 46%), mean rates of ripples (0.3 vs 0.09 / minute, p < 0.0001) and spike ripples (0.6 vs 0.06 / minute, p < 0.05) were significantly higher, with no difference in spike rates (1.7 vs 3.0 / minute, p = 0.38). Of those 3 markers, ripples showed the best predictive value (area under the curve [AUC]_ripples = 0.88). The optimal threshold for ripples was calculated to be ≥ 0.125 / minute with a sensitivity of 87% and specificity of 85%. Ripple rates were negatively correlated to days passing before epilepsy-diagnosis (R = -0.59, p < 0.0001) and time to a second seizure (R = -0.64, 95% confidence interval [CI] = -0.77 to 0.43, p < 0.0001).

INTERPRETATION

We could show that in a cohort of children with a first unprovoked seizure, ripples predict the development of epilepsy better than spikes or spike ripples and might be useful biomarkers in the estimation of prognosis and question of treatment. ANN NEUROL 2021;89:134-142.