Data mining neocortical high-frequency oscillations in epilepsy and controls

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.

Automatic Detection of Scalp High-Frequency Oscillations Based on Deep Learning

Scalp high-frequency oscillations (sHFOs) are a promising non-invasive biomarker of epilepsy. However, the visual marking of sHFOs is a time-consuming and subjective process, existing automatic detectors based on single-dimensional analysis have difficulty with accurately eliminating artifacts and thus do not provide sufficient reliability to meet clinical needs. Therefore, we propose a high-performance sHFOs detector based on a deep learning algorithm. An initial detection module was designed to extract candidate high-frequency oscillations. Then, one-dimensional (1D) and two-dimensional (2D) deep learning models were designed, respectively. Finally, the weighted voting method is used to combine the outputs of the two model. In experiments, the precision, recall, specificity and F1-score were 83.44%, 83.60%, 96.61% and 83.42%, respectively, on average and the kappa coefficient was 80.02%. In addition, the proposed detector showed a stable performance on multi-centre datasets. Our sHFOs detector demonstrated high robustness and generalisation ability, which indicates its potential applicability as a clinical assistance tool. The proposed sHFOs detector achieves an accurate and robust method via deep learning algorithm.

Interictal invasive very high-frequency oscillations in resting awake state and sleep

Interictal very high-frequency oscillations (VHFOs, 500-2000 Hz) in a resting awake state seem to be, according to a precedent study of our team, a more specific predictor of a good outcome of the epilepsy surgery compared to traditional interictal high-frequency oscillations (HFOs, 80-500 Hz). In this study, we retested this hypothesis on a larger cohort of patients. In addition, we also collected patients' sleep data and hypothesized that the occurrence of VHFOs in sleep will be greater than in resting state. We recorded interictal invasive electroencephalographic (iEEG) oscillations in 104 patients with drug-resistant epilepsy in a resting state and in 35 patients during sleep. 21 patients in the rest study and 11 patients in the sleep study met the inclusion criteria (interictal HFOs and VHFOs present in iEEG recordings, a surgical intervention and a postoperative follow-up of at least 1 year) for further evaluation of iEEG data. In the rest study, patients with good postoperative outcomes had significantly higher ratio of resected contacts with VHFOs compared to HFOs. In sleep, VHFOs were more abundant than in rest and the percentage of resected contacts in patients with good and poor outcomes did not considerably differ in any type of oscillations. In conclusion, (1) our results confirm, in a larger patient cohort, our previous work about VHFOs being a specific predictor of the area which needs to be resected; and (2) that more frequent sleep VHFOs do not further improve the results.

High-resolution neural recordings improve the accuracy of speech decoding

Patients suffering from debilitating neurodegenerative diseases often lose the ability to communicate, detrimentally affecting their quality of life. One solution to restore communication is to decode signals directly from the brain to enable neural speech prostheses. However, decoding has been limited by coarse neural recordings which inadequately capture the rich spatio-temporal structure of human brain signals. To resolve this limitation, we performed high-resolution, micro-electrocorticographic (µECoG) neural recordings during intra-operative speech production. We obtained neural signals with 57× higher spatial resolution and 48% higher signal-to-noise ratio compared to macro-ECoG and SEEG. This increased signal quality improved decoding by 35% compared to standard intracranial signals. Accurate decoding was dependent on the high-spatial resolution of the neural interface. Non-linear decoding models designed to utilize enhanced spatio-temporal neural information produced better results than linear techniques. We show that high-density µECoG can enable high-quality speech decoding for future neural speech prostheses.

Flexible, high-resolution cortical arrays with large coverage capture microscale high-frequency oscillations in patients with epilepsy

OBJECTIVE

Effective surgical treatment of drug-resistant epilepsy depends on accurate localization of the epileptogenic zone (EZ). High-frequency oscillations (HFOs) are potential biomarkers of the EZ. Previous research has shown that HFOs often occur within submillimeter areas of brain tissue and that the coarse spatial sampling of clinical intracranial electrode arrays may limit the accurate capture of HFO activity. In this study, we sought to characterize microscale HFO activity captured on thin, flexible microelectrocorticographic (μECoG) arrays, which provide high spatial resolution over large cortical surface areas.

METHODS

We used novel liquid crystal polymer thin-film μECoG arrays (.76-1.72-mm intercontact spacing) to capture HFOs in eight intraoperative recordings from seven patients with epilepsy. We identified ripple (80-250 Hz) and fast ripple (250-600 Hz) HFOs using a common energy thresholding detection algorithm along with two stages of artifact rejection. We visualized microscale subregions of HFO activity using spatial maps of HFO rate, signal-to-noise ratio, and mean peak frequency. We quantified the spatial extent of HFO events by measuring covariance between detected HFOs and surrounding activity. We also compared HFO detection rates on microcontacts to simulated macrocontacts by spatially averaging data.

RESULTS

We found visually delineable subregions of elevated HFO activity within each μECoG recording. Forty-seven percent of HFOs occurred on single 200-μm-diameter recording contacts, with minimal high-frequency activity on surrounding contacts. Other HFO events occurred across multiple contacts simultaneously, with covarying activity most often limited to a .95-mm radius. Through spatial averaging, we estimated that macrocontacts with 2-3-mm diameter would only capture 44% of the HFOs detected in our μECoG recordings.

SIGNIFICANCE

These results demonstrate that thin-film microcontact surface arrays with both highresolution and large coverage accurately capture microscale HFO activity and may improve the utility of HFOs to localize the EZ for treatment of drug-resistant epilepsy.

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

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

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.

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

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

Normative intracranial EEG maps epileptogenic tissues in focal epilepsy

Planning surgery for patients with medically refractory epilepsy often requires recording seizures using intracranial EEG. Quantitative measures derived from interictal intracranial EEG yield potentially appealing biomarkers to guide these surgical procedures; however, their utility is limited by the sparsity of electrode implantation as well as the normal confounds of spatiotemporally varying neural activity and connectivity. We propose that comparing intracranial EEG recordings to a normative atlas of intracranial EEG activity and connectivity can reliably map abnormal regions, identify targets for invasive treatment and increase our understanding of human epilepsy. Merging data from the Penn Epilepsy Center and a public database from the Montreal Neurological Institute, we aggregated interictal intracranial EEG retrospectively across 166 subjects comprising >5000 channels. For each channel, we calculated the normalized spectral power and coherence in each canonical frequency band. We constructed an intracranial EEG atlas by mapping the distribution of each feature across the brain and tested the atlas against data from novel patients by generating a z-score for each channel. We demonstrate that for seizure onset zones within the mesial temporal lobe, measures of connectivity abnormality provide greater distinguishing value than univariate measures of abnormal neural activity. We also find that patients with a longer diagnosis of epilepsy have greater abnormalities in connectivity. By integrating measures of both single-channel activity and inter-regional functional connectivity, we find a better accuracy in predicting the seizure onset zones versus normal brain (area under the curve = 0.77) compared with either group of features alone. We propose that aggregating normative intracranial EEG data across epilepsy centres into a normative atlas provides a rigorous, quantitative method to map epileptic networks and guide invasive therapy. We publicly share our data, infrastructure and methods, and propose an international framework for leveraging big data in surgical planning for refractory epilepsy.

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.

Interictal high frequency background activity as a biomarker of epileptogenic tissue

High frequency oscillations (HFOs) are very brief events that are a well-established biomarker of the epileptogenic zone (EZ) but are rare and comprise only a tiny fraction of the total recorded EEG. We hypothesize that the interictal high frequency ‘background’ data, which has received little attention but represents the majority of the EEG record, also may contain additional, novel information for identifying the EZ. We analysed intracranial EEG (30–500 Hz frequency range) acquired from 24 patients who underwent resective surgery. We computed 38 quantitative features based on all usable, interictal data (63–307 h per subject), excluding all detected HFOs. We assessed association between each feature and the seizure onset zone (SOZ) and resected volume (RV) using logistic regression. A pathology score per channel was also created via principle component analysis and logistic regression, using hold-out-one-patient cross-validation to avoid in-sample training. Association of the pathology score with the SOZ and RV was quantified using an asymmetry measure. Many features were associated with the SOZ: 23/38 features had odds ratios >1.3 or <0.7 and 17/38 had odds ratios different than zero with high significance (P < 0.001/39, logistic regression with Bonferroni Correction). The pathology score, the rate of HFOs, and their channel-wise product were each strongly associated with the SOZ [median asymmetry ≥0.44, good surgery outcome patients; median asymmetry ≥0.40, patients with other outcomes; 95% confidence interval (CI) > 0.27 in both cases]. The pathology score and the channel-wise product also had higher asymmetry with respect to the SOZ than the HFO rate alone (median difference in asymmetry ≥0.18, 95% CI >0.05). These results support that the high frequency background data contains useful information for determining the EZ, distinct and complementary to information from detected HFOs. The concordance between the high frequency activity pathology score and the rate of HFOs appears to be a better biomarker of epileptic tissue than either measure alone.

Phase-Amplitude Coupling and Epileptogenic Zone Localization of Frontal Epilepsy Based on Intracranial EEG

Objective: This study aimed to explore the characteristics of phase-amplitude coupling in patients with frontal epilepsy based on their electrocorticography data, in order to identify the localization of epileptic regions and further guide clinical resection surgery.

Methods: We adopted the modulation index based on the Kullback-Leibler distance, phase-amplitude coupling co-modulogram, and time-varying phase-amplitude modulogram to explore the temporal-spatial patterns and characterization of PAC strength during the period from inter- seizure to post-seizure. Taking the resected area as the gold standard, the epileptogenic zone was located based on MI values of 7 different seizure periods, and the accuracy of localization was measured by the area under the receiver operating curve.

Results: (1) The PAC in the inter- and pre-seizure periods was weak and paroxysmal, but strong PAC channels were confined more to the seizure-onset zone and resection region. PAC during the seizure period was intense and persistent, but gradually deviated from the seizure-onset zone. (2) The characteristics of coupling strength of the inter- and pre-seizure EEG can be used to accurately locate the epileptogenic zone, which is better than that in periods after the beginning of a seizure. (3) In an epileptic seizure, the preferred phases of coupling were usually in the rising branches at the pre- and early-seizure stages, while those in the middle- and terminal-seizure were usually in the falling branch. We thus speculate that the coupling occurred in the rising branch can promote the recruitment of abnormal discharge, while the coupling occurred in the falling branch can inhibit the abnormal discharge.

Conclusion: The findings suggest that the phase-amplitude coupling during inter- and pre-seizure is a promising marker of epileptic focus location. The preferred phase of coupling changed regularly with the time of epileptic seizure, suggesting that the surge and suppression of abnormal discharges are related to different phases.

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.

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.

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.

Flexible, high-resolution thin-film electrodes for human and animal neural research

Objective.Brain functions such as perception, motor control, learning, and memory arise from the coordinated activity of neuronal assemblies distributed across multiple brain regions. While major progress has been made in understanding the function of individual neurons, circuit interactions remain poorly understood. A fundamental obstacle to deciphering circuit interactions is the limited availability of research tools to observe and manipulate the activity of large, distributed neuronal populations in humans. Here we describe the development, validation, and dissemination of flexible, high-resolution, thin-film (TF) electrodes for recording neural activity in animals and humans.Approach.We leveraged standard flexible printed-circuit manufacturing processes to build high-resolution TF electrode arrays. We used biocompatible materials to form the substrate (liquid crystal polymer; LCP), metals (Au, PtIr, and Pd), molding (medical-grade silicone), and 3D-printed housing (nylon). We designed a custom, miniaturized, digitizing headstage to reduce the number of cables required to connect to the acquisition system and reduce the distance between the electrodes and the amplifiers. A custom mechanical system enabled the electrodes and headstages to be pre-assembled prior to sterilization, minimizing the setup time required in the operating room. PtIr electrode coatings lowered impedance and enabled stimulation. High-volume, commercial manufacturing enables cost-effective production of LCP-TF electrodes in large quantities.Main Results. Our LCP-TF arrays achieve 25× higher electrode density, 20× higher channel count, and 11× reduced stiffness than conventional clinical electrodes. We validated our LCP-TF electrodes in multiple human intraoperative recording sessions and have disseminated this technology to >10 research groups. Using these arrays, we have observed high-frequency neural activity with sub-millimeter resolution.Significance.Our LCP-TF electrodes will advance human neuroscience research and improve clinical care by enabling broad access to transformative, high-resolution electrode arrays.

Region-specific Effects of Early-life Status Epilepticus on the Adult Hippocampal CA3 - Medial Entorhinal Cortex Circuitry In vitro: Focus on Interictal Spikes and Concurrent High-frequency Oscillations

Convulsive status epilepticus (SE) in immature life is often associated with lasting neurobiological changes. We provoked SE by pentylenetetrazole in postnatal day 20 rat pups and examined communication modalities between the temporal hippocampus and medial entorhinal cortex (mEC) in vitro. After a minimum of 40 days post-SE, we prepared combined temporal hippocampal - medial entorhinal cortex (mEC) slices from conditioned (SE) and naïve (N) adult rats and recorded 4-aminopyridine-induced spontaneous epileptiform interictal-like discharges (IED) simultaneously from CA3 and mEC layer V-VI. We analyzed IED frequency and high frequency oscillations (HFOs) in intact slices and after surgical separation of hippocampus from mEC, by two successive incisions (Schaffer collateral cut, Parasubiculum cut). In all slices, IED frequency was higher in CA3 vs mEC (5N, 4SE) and Raster plots indicated no temporal coincidence between them either in intact or in CA1-cut slices (4N, 4SE). IED frequency was significantly higher in SE mEC, but similar in SE and N CA3, independently of connectivity state. Ripples (R) and Fast Ripples (FR) coincided with IEDs and their power differed between SE and N intact slices (22N, 12SE), both in CA3 and mEC. CA3 FR/R ratios were higher in the absence of mEC (14N, 8SE). Moreover, SE (vs N) slices showed significantly higher FR/R ratios independently of the presence of mEC. Taken together, these findings suggest lasting effects of immature SE in network dynamics governing hippocampal-entorhinal communication which may impact adult cognitive, behavioral, and/or seizure threshold sequalae.

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.

Physiological Ripples Associated With Sleep Spindles Can Be Identified in Patients With Refractory Epilepsy Beyond Mesio-Temporal Structures

Introduction: High frequency oscillations (HFO) are promising biomarkers of epileptic tissue. While group analysis suggested a correlation between surgical removal of HFO generating tissue and seizure free outcome, HFO could not predict seizure outcome on an individual patient level. One possible explanation is the lack of differentiation between physiological and epileptic HFO. In the mesio-temporal lobe, a proportion of physiological ripples can be identified by their association with scalp sleep spindles. Spike associated ripples in contrast can be considered epileptic. This study investigated whether categorizing ripples by the co-occurrence with sleep spindles or spikes improves outcome prediction after surgery. Additionally, it aimed to investigate whether spindle-ripple association is limited to the mesio-temporal lobe structures or visible across the whole brain. Methods: We retrospectively analyzed EEG of 31 patients with chronic intracranial EEG. Sleep spindles in scalp EEG and ripples and epileptic spikes in iEEG were automatically detected. Three ripple subtypes were obtained: SpindleR, Non-SpindleR, and SpikeR. Rate ratios between removed and non-removed brain areas were calculated. We compared the distinct ripple subtypes and their rates in different brain regions, inside and outside seizure onset areas and between patients with good and poor seizure outcome. Results: SpindleR were found across all brain regions. SpikeR had significantly higher rates in the SOZ than in Non-SOZ channels. A significant positive correlation between removal of ripple-events and good outcome was found for the mixed ripple group (r_s = 0.43, p = 0.017) and for ripples not associated with spindles (r_s=0.40, p = 0.044). Also, a significantly high proportion of spikes associated with ripples were removed in seizure free patients (p = 0.036). Discussion: SpindleR are found in mesio-temporal and neocortical structures, indicating that ripple-spindle-coupling might have functional importance beyond mesio-temporal structures. Overall, the proportion of SpindleR was low and separating spindle and spike associated ripples did not improve outcome prediction in our patient group. SpindleR analysis therefore can be a tool to identify physiological events but needs to be used in combination with other methods to have clinical relevance.

Viability of Preictal High-Frequency Oscillation Rates as a Biomarker for Seizure Prediction

Motivation: There is an ongoing search for definitive and reliable biomarkers to forecast or predict imminent seizure onset, but to date most research has been limited to EEG with sampling rates <1,000 Hz. High-frequency oscillations (HFOs) have gained acceptance as an indicator of epileptic tissue, but few have investigated the temporal properties of HFOs or their potential role as a predictor in seizure prediction. Here we evaluate time-varying trends in preictal HFO rates as a potential biomarker of seizure prediction.

Methods: HFOs were identified for all interictal and preictal periods with a validated automated detector in 27 patients who underwent intracranial EEG monitoring. We used LASSO logistic regression with several features of the HFO rate to distinguish preictal from interictal periods in each individual. We then tested these models with held-out data and evaluated their performance with the area-under-the-curve (AUC) of their receiver-operating curve (ROC). Finally, we assessed the significance of these results using non-parametric statistical tests.

Results: There was variability in the ability of HFOs to discern preictal from interictal states across our cohort. We identified a subset of 10 patients in whom the presence of the preictal state could be successfully predicted better than chance. For some of these individuals, average AUC in the held-out data reached higher than 0.80, which suggests that HFO rates can significantly differentiate preictal and interictal periods for certain patients.

Significance: These findings show that temporal trends in HFO rate can predict the preictal state better than random chance in some individuals. Such promising results indicate that future prediction efforts would benefit from the inclusion of high-frequency information in their predictive models and technological architecture.

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.

Seizure Onset Zone Lateralization Using a Non-linear Analysis of Micro vs. Macro Electroencephalographic Recordings During Seizure-Free Stages of the Sleep-Wake Cycle From Epilepsy Patients

The application of non-linear signal analysis techniques to biomedical data is key to improve our knowledge about complex physiological and pathological processes. In particular, the use of non-linear techniques to study electroencephalographic (EEG) recordings can provide an advanced characterization of brain dynamics. In epilepsy these dynamics are altered at different spatial scales of neuronal organization. We therefore apply non-linear signal analysis to EEG recordings from epilepsy patients derived with intracranial hybrid electrodes, which are composed of classical macro contacts and micro wires. Thereby, these electrodes record EEG at two different spatial scales. Our aim is to test the degree to which the analysis of the EEG recorded at these different scales allows us to characterize the neuronal dynamics affected by epilepsy. For this purpose, we retrospectively analyzed long-term recordings performed during five nights in three patients during which no seizures took place. As a benchmark we used the accuracy with which this analysis allows determining the hemisphere that contains the seizure onset zone, which is the brain area where clinical seizures originate. We applied the surrogate-corrected non-linear predictability score (ψ), a non-linear signal analysis technique which was shown previously to be useful for the lateralization of the seizure onset zone from classical intracranial EEG macro contact recordings. Higher values of ψ were found predominantly for signals recorded from the hemisphere containing the seizure onset zone as compared to signals recorded from the opposite hemisphere. These differences were found not only for the EEG signals recorded with macro contacts, but also for those recorded with micro wires. In conclusion, the information obtained from the analysis of classical macro EEG contacts can be complemented by the one of micro wire EEG recordings. This combined approach may therefore help to further improve the degree to which quantitative EEG analysis can contribute to the diagnostics in epilepsy patients.

Distinguishing false and true positive detections of high frequency oscillations

Objective.

High frequency oscillations (HFOs) are a promising biomarker of tissue that instigates seizures. However, ambiguous data and random background fluctuations can cause any HFO detector (human or automated) to falsely label non-HFO data as an HFO (a false positive detection). The objective of this paper was to identify quantitative features of HFOs that distinguish between true and false positive detections.

Approach.

Feature selection was performed using background data in multi-day, interictal intracranial recordings from ten patients. We selected the feature most similar between randomly selected segments of background data and HFOs detected in surrogate background data (false positive detections by construction). We then compared these results with fuzzy clustering of detected HFOs in clinical data to verify the feature’s applicability. We validated the feature is sensitive to false versus true positive HFO detections by using an independent data set (six subjects) scored for HFOs by three human reviewers. Lastly, we compared the effect of redacting putative false positive HFO detections on the distribution of HFOs across channels and their association with seizure onset zone (SOZ) and resected volume (RV).

Main results.

Of the 15 analyzed features, the analysis selected only skewness of the curvature (skewCurve). The feature was validated in human scored data to be associated with distinguishing true and false positive HFO detections. Automated HFO detections with higher skewCurve were more focal based on entropy measures and had increased localization to both the SOZ and RV.

Significance.

We identified a quantitative feature of HFOs which helps distinguish between true and false positive detections. Redacting putative false positive HFO detections improves the specificity of HFOs as a biomarker of epileptic tissue.

Quantitative EEG biomarkers for epilepsy and their relation to chemical biomarkers

The electroencephalogram (EEG) is the most important method to diagnose epilepsy. In clinical settings, it is evaluated by experts who identify patterns visually. Quantitative EEG is the application of digital signal processing to clinical recordings in order to automatize diagnostic procedures, and to make patterns visible that are hidden to the human eye. The EEG is related to chemical biomarkers, as electrical activity is based on chemical signals. The most well-known chemical biomarkers are blood laboratory tests to identify seizures after they have happened. However, research on chemical biomarkers is much less extensive than research on quantitative EEG, and combined studies are rarely published, but highly warranted. Quantitative EEG is as old as the EEG itself, but still, the methods are not yet standard in clinical practice. The most evident application is an automation of manual work, but also a quantitative description and localization of interictal epileptiform events as well as seizures can reveal important hints for diagnosis and contribute to presurgical evaluation. In addition, the assessment of network characteristics and entropy measures were found to reveal important insights into epileptic brain activity. Application scenarios of quantitative EEG in epilepsy include seizure prediction, pharmaco-EEG, treatment monitoring, evaluation of cognition, and neurofeedback. The main challenges to quantitative EEG are poor reliability and poor generalizability of measures, as well as the need for individualization of procedures. A main hindrance for quantitative EEG to enter clinical routine is also that training is not yet part of standard curricula for clinical neurophysiologists.

Preictal variability of high-frequency oscillation rates in refractory epilepsy

OBJECTIVE

High-frequency oscillations (HFOs) have shown promising utility in the spatial localization of the seizure onset zone for patients with focal refractory epilepsy. Comparatively few studies have addressed potential temporal variations in HFOs, or their role in the preictal period. Here, we introduce a novel evaluation of the instantaneous HFO rate through interictal and peri-ictal epochs to assess their usefulness in identifying imminent seizure onset.

METHODS

Utilizing an automated HFO detector, we analyzed intracranial electroencephalographic data from 30 patients with refractory epilepsy undergoing long-term presurgical evaluation. We evaluated HFO rates both as a 30-minute average and as a continuous function of time and used nonparametric statistical methods to compare individual and population-level differences in rate during peri-ictal and interictal periods.

RESULTS

Mean HFO rate was significantly higher for all epochs in seizure onset zone channels versus other channels. Across the 30 patients of our cohort, we found no statistically significant differences in mean HFO rate during preictal and interictal epochs. For continuous HFO rates in seizure onset zone channels, however, we found significant population-wide increases in preictal trends relative to interictal periods. Using a data-driven analysis, we identified a subset of 11 patients in whom either preictal HFO rates or their continuous trends were significantly increased relative to those of interictal baseline and the rest of the population.

SIGNIFICANCE

These results corroborate existing findings that HFO rates within epileptic tissue are higher during interictal periods. We show this finding is also present in preictal, ictal, and postictal data, and identify a novel biomarker of preictal state: an upward trend in HFO rate leading into seizures in some patients. Overall, our findings provide preliminary evidence that HFOs can function as a temporal biomarker of seizure onset.

Association Between Interictal High-Frequency Oscillations and Slow Wave in Refractory Focal Epilepsy With Good Surgical Outcome

High-frequency oscillations (HFOs) have been proposed as a promising biomarker of the epileptogenic zone (EZ). But accurate delineation of EZ based on HFOs is still challenging. Our study compared HFOs from EZ and non-EZ on the basis of their associations with interictal slow waves, aiming at exploring a new way to localize EZ. Nineteen medically intractable epilepsy patients with good surgical outcome were included. Five minute interictal intracranial electroencephalography (EEG) epochs of slow-wave sleep were randomly selected; then ripples (80–200 Hz), fast ripples (FRs; 200–500 Hz), and slow waves (0.1–4 Hz) were automatically analyzed. The EZ and non-EZ were identified by resection range during the surgeries. We found that both ripples and FRs superimposed more frequently on slow waves in EZ than in non-EZ (P < 0.01). Although ripples preferred to occur on the down state of slow waves in both two groups, ripples in EZ tended to be closer to the down-state peak of slow wave than in non-EZ (-174 vs. -231 ms, P = 0.008). As for FR, no statistical difference was found between the two groups (P = 0.430). Additionally, slow wave-containing ripples in EZ had a steeper slope (1.7 vs. 1.5 μV/ms, P < 0.001) and wider distribution ratio (32.3 vs. 30.1%, P < 0.001) than those in the non-EZ. But for slow wave-containing FR, only a steeper slope (1.7 vs. 1.4 μV/ms, P < 0.001) was observed. Our study innovatively compared the different features of association between HFOs and slow wave in EZ and non-EZ from refractory focal epilepsy with good surgical outcome, proposing a new method to localize EZ and facilitating the surgical plan.

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

High-frequency oscillations (HFOs) in intracranial electroencephalography (EEG) are a promising biomarker of the epileptogenic zone and tool for surgical planning. Many studies have shown that a high rate of HFOs (number per minute) is correlated with the seizure-onset zone, and complete removal of HFO-generating brain regions has been associated with seizure-free outcome after surgery. In order to use HFOs as a biomarker, these transient events must first be detected in electrophysiological data. Because visual detection of HFOs is time-consuming and subject to low interrater reliability, many automated algorithms have been developed, and they are being used increasingly for such studies. However, there is little guidance on how to select an algorithm, implement it in a clinical setting, and validate the performance. Therefore, we aim to review automated HFO detection algorithms, focusing on conceptual similarities and differences between them. We summarize the standard steps for data pre-processing, as well as post-processing strategies for rejection of false-positive detections. We also detail four methods for algorithm testing and validation, and we describe the specific goal achieved by each one. We briefly review direct comparisons of automated algorithms applied to the same data set, emphasizing the importance of optimizing detection parameters. Then, to assess trends in the use of automated algorithms and their potential for use in clinical studies, we review evidence for the relationship between automatically detected HFOs and surgical outcome. We conclude with practical recommendations and propose standards for the selection, implementation, and validation of automated HFO-detection algorithms.

Expert consensus on clinical applications of high-frequency oscillations in epilepsy

Studies in animal models of epilepsy and pre-surgical patients have unanimously found a strong correlation between high-frequency oscillations (HFOs, > 80 Hz) and the epileptogenic zone, suggesting that HFOs can be a potential biomarker of epileptogenicity and epileptogenesis. This consensus includes the definition and standard detection techniques of HFOs, the localizing value of pathological HFOs for epileptic foci, and different ways to distinguish physiological from epileptic HFOs. The latest clinical applications of HFOs in epilepsy and the related findings are also discussed. HFOs will advance our understanding of the pathophysiology of epilepsy.

Detection of anomalous high-frequency events in human intracranial EEG

Objective

High‐frequency oscillations (HFOs) are a promising biomarker for the epileptogenic zone. However, no physiological definition of an HFO has been established, so detection relies on the empirical definition of an HFO derived from visual observation. This can bias estimates of HFO features such as amplitude and duration, thereby hindering their utility as biomarkers. Therefore, we set out to develop an algorithm that detects high‐frequency events in the intracranial EEG that are morphologically distinct from background without requiring assumptions about event amplitude or shape.

Method

We propose the anomaly detection algorithm (ADA), which uses unsupervised machine learning to identify segments of data that are distinct from the background. We apply ADA and a standard HFO detector using a root mean square amplitude threshold to intracranial EEG from 11 patients undergoing evaluation for epilepsy surgery. The rate, amplitude, and duration of the detected events and the percent overlap between the two detectors are compared.

Result

In the seizure onset zone (SOZ), ADA detected a subset of conventional HFOs. In non‐SOZ channels, ADA detected at least twice as many events as the standard approach, including some conventional HFOs; however, ADA also identified many low and intermediate amplitude events missed by the standard amplitude‐based method. The rate of ADA events was similar across all channels; however, the amplitude of ADA events was significantly higher in SOZ channels (P < .0045), and the amplitude measurement was more stable over time than the HFO rate, as indicated by a lower coefficient of variation (P < .0125).

Significance

ADA does not require human supervision, parameter optimization, or prior assumptions about event shape, amplitude, or duration. Our results suggest that the algorithm's estimate of event amplitude may differentiate SOZ and non‐SOZ channels. Further studies will examine the utility of HFO amplitude as a biomarker for epilepsy surgical outcome.

Multiband entropy-based feature-extraction method for automatic identification of epileptic focus based on high-frequency components in interictal iEEG

Presurgical investigations for categorizing focal patterns are crucial, leading to localization and surgical removal of the epileptic focus. This paper presents a machine learning approach using information theoretic features extracted from high-frequency subbands to detect the epileptic focus from interictal intracranial electroencephalogram (iEEG). It is known that high-frequency subbands (>80 Hz) include important biomarkers such as high-frequency oscillations (HFOs) for identifying epileptic focus commonly referred to as the seizure onset zone (SOZ). In this analysis, the multi-channel interictal iEEG signals were splitted into segments and each segment was decomposed into multiple high-frequency subbands. The different types of entropy were calculated for each of the subbands and the sparse linear discriminant analysis (sLDA) was applied to select the prominent entropy features. Due to the imbalance of SOZ and non-SOZ channels in iEEG data, the use of machine learning techniques is always tricky. To deal with the imbalanced learning problem, an adaptive synthetic oversampling approach (ADASYN) with radial basis function kernel-based SVM was used to detect the focal segments. Finally, the epileptic focus was identified based on detection of focal segments on SOZ and non-SOZ channels. Eight patients were examined to observe the efficiency of the automatic detector. The experimental results and statistical tests indicate that the proposed automatic detector can identify the epileptic focus accurately and efficiently.

Optically pumped magnetoencephalography in epilepsy

We demonstrate the first use of Optically Pumped Magnetoencephalography (OP-MEG) in an epilepsy patient with unrestricted head movement. Current clinical MEG uses a traditional SQUID system, where sensors are cryogenically cooled and housed in a helmet in which the patient's head is fixed. Here, we use a different type of sensor (OPM), which operates at room temperature and can be placed directly on the patient's scalp, permitting free head movement. We performed OP-MEG recording in a patient with refractory focal epilepsy. OP-MEG-identified analogous interictal activity to scalp EEG, and source localized this activity to an appropriate brain region.

Scalp ripples as prognostic biomarkers of epileptogenicity in pediatric surgery

Objective

To assess the ability of high‐density Electroencephalography (HD‐EEG) and magnetoencephalography (MEG) to localize interictal ripples, distinguish between ripples co‐occurring with spikes (ripples‐on‐spike) and independent from spikes (ripples‐alone), and evaluate their localizing value as biomarkers of epileptogenicity in children with medically refractory epilepsy.

Methods

We retrospectively studied 20 children who underwent epilepsy surgery. We identified ripples on HD‐EEG and MEG data, localized their generators, and compared them with intracranial EEG (icEEG) ripples. When ripples and spikes co‐occurred, we performed source imaging distinctly on the data above 80 Hz (to localize ripples) and below 70 Hz (to localize spikes). We assessed whether missed resection of ripple sources predicted poor outcome, separately for ripples‐on‐spikes and ripples‐alone. Similarly, predictive value of spikes was calculated.

Results

We observed scalp ripples in 16 patients (10 good outcome). Ripple sources were highly concordant to the icEEG ripples (HD‐EEG concordance: 79%; MEG: 83%). When ripples and spikes co‐occurred, their sources were spatially distinct in 83‐84% of the cases. Removing the sources of ripples‐on‐spikes predicted good outcome with 90% accuracy for HD‐EEG (P = 0.008) and 86% for MEG (P = 0.044). Conversely, removing ripples‐alone did not predict outcome. Resection of spike sources (generated at the same time as ripples) predicted good outcome for HD‐EEG (P = 0.036; accuracy = 87%), while did not reach significance for MEG (P = 0.1; accuracy = 80%).

Interpretation

HD‐EEG and MEG localize interictal ripples with high precision in children with refractory epilepsy. Scalp ripples‐on‐spikes are prognostic, noninvasive biomarkers of epileptogenicity, since removing their cortical generators predicts good outcome. Conversely, scalp ripples‐alone are most likely generated by non‐epileptogenic areas.

Temporal Pattern of Ripple Events in Temporal Lobe Epilepsy: Towards a Pattern-based Localization of the Seizure Onset Zone

Ripples (80-250 Hz) are brief high-frequency oscillations that are often detected in intracranial EEG (iEEG) and are currently investigated as a potential biomarker to facilitate the Iocalization of the seizure onset zone (SOZ) in patients with drug-resistant epilepsy. While the rate and shape of these oscillations have been positively correlated with the SOZ, the temporal pattern of these oscillations in the epileptic brain still requires investigation. In this study, we investigate the temporal pattern of ripple events in five patients with temporal lobe epilepsy (TLE), which is one of the most common forms of epilepsy. The rate of ripple events is positively correlated with the SOZ in TLE but its diagnostic utility in localizing the SOZ remains unclear, which suggests that additional ripple-related features should be investigated. By combining point process modeling and cluster analysis, we show that a recurrent, non-stationary bursting pattern characterizes the SOZ channels consistently across patients, while the non-SOZ channels have poor between-channel similarity and no consistent pattern over time nor across patients. Furthermore, the degree of separation between SOZ and non-SOZ model parameter vectors is significantly higher (ANOVA test, ${P}$-value $P\lt 0.01$) than the degree of separation between the ripple rates, which suggests that the temporal pattern more than the rate may contribute to the pre- surgical Iocalization of the SOZ.

High Frequency Oscillations in Epilepsy: Detection Methods and Considerations in Clinical Application

High frequency oscillations (HFOs) is a brain activity observed in electroencephalography (EEG) in frequency ranges between 80–500 Hz. HFOs can be classified into ripples (80–200 Hz) and fast ripples (200–500 Hz) by their distinctive characteristics. Recent studies reported that both ripples and fast fipples can be regarded as a new biomarker of epileptogenesis and ictogenesis. Previous studies verified that HFOs are clinically important both in patients with mesial temporal lobe epilepsy and neocortical epilepsy. Also, in epilepsy surgery, patients with higher resection ratio of brain regions with HFOs showed better outcome than a group with lower resection ratio. For clinical application of HFOs, it is important to delineate HFOs accurately and discriminate them from artifacts. There have been technical improvements in detecting HFOs by developing various detection algorithms. Still, there is a difficult issue on discriminating clinically important HFOs among detected HFOs, where both quantitative and subjective approaches are suggested. This paper is a review on published HFO studies focused on clinical findings and detection techniques of HFOs as well as tips for clinical applications.

Stereotyped high-frequency oscillations discriminate seizure onset zones and critical functional cortex in focal epilepsy

High-frequency oscillations in local field potentials recorded with intracranial EEG are putative biomarkers of seizure onset zones in epileptic brain. However, localized 80-500 Hz oscillations can also be recorded from normal and non-epileptic cerebral structures. When defined only by rate or frequency, physiological high-frequency oscillations are indistinguishable from pathological ones, which limit their application in epilepsy presurgical planning. We hypothesized that pathological high-frequency oscillations occur in a repetitive fashion with a similar waveform morphology that specifically indicates seizure onset zones. We investigated the waveform patterns of automatically detected high-frequency oscillations in 13 epilepsy patients and five control subjects, with an average of 73 subdural and intracerebral electrodes recorded per patient. The repetitive oscillatory waveforms were identified by using a pipeline of unsupervised machine learning techniques and were then correlated with independently clinician-defined seizure onset zones. Consistently in all patients, the stereotypical high-frequency oscillations with the highest degree of waveform similarity were localized within the seizure onset zones only, whereas the channels generating high-frequency oscillations embedded in random waveforms were found in the functional regions independent from the epileptogenic locations. The repetitive waveform pattern was more evident in fast ripples compared to ripples, suggesting a potential association between waveform repetition and the underlying pathological network. Our findings provided a new tool for the interpretation of pathological high-frequency oscillations that can be efficiently applied to distinguish seizure onset zones from functionally important sites, which is a critical step towards the translation of these signature events into valid clinical biomarkers.awx374media15721572971001.

Localizing epileptogenic regions using high-frequency oscillations and machine learning

Pathological high frequency oscillations (HFOs) are putative neurophysiological biomarkers of epileptogenic brain tissue. Utilizing HFOs for epilepsy surgery planning offers the promise of improved seizure outcomes for patients with medically refractory epilepsy. This review discusses possible machine learning strategies that can be applied to HFO biomarkers to better identify epileptogenic regions. We discuss the role of HFO rate, and utilizing features such as explicit HFO properties (spectral content, duration, and power) and phase-amplitude coupling for distinguishing pathological HFO (pHFO) events from physiological HFO events. In addition, the review highlights the importance of neuroanatomical localization in machine learning strategies.

Scalp recorded spike ripples predict seizure risk in childhood epilepsy better than spikes

In the past decade, brief bursts of fast oscillations in the ripple range have been identified in the scalp EEG as a promising non-invasive biomarker for epilepsy. However, investigation and clinical application of this biomarker have been limited because standard approaches to identify these brief, low amplitude events are difficult, time consuming, and subjective. Recent studies have demonstrated that ripples co-occurring with epileptiform discharges ('spike ripple events') are easier to detect than ripples alone and have greater pathological significance. Here, we used objective techniques to quantify spike ripples and test whether this biomarker predicts seizure risk in childhood epilepsy. We evaluated spike ripples in scalp EEG recordings from a prospective cohort of children with a self-limited epilepsy syndrome, benign epilepsy with centrotemporal spikes, and healthy control children. We compared the rate of spike ripples between children with epilepsy and healthy controls, and between children with epilepsy during periods of active disease (active, within 1 year of seizure) and after a period of sustained seizure-freedom (seizure-free, >1 year without seizure), using semi-automated and automated detection techniques. Spike ripple rate was higher in subjects with active epilepsy compared to healthy controls (P = 0.0018) or subjects with epilepsy who were seizure-free ON or OFF medication (P = 0.0018). Among epilepsy subjects with spike ripples, each month seizure-free decreased the odds of a spike ripple by a factor of 0.66 [95% confidence interval (0.47, 0.91), P = 0.021]. Comparing the diagnostic accuracy of the presence of at least one spike ripple versus a classic spike event to identify group, we found comparable sensitivity and negative predictive value, but greater specificity and positive predictive value of spike ripples compared to spikes (P = 0.016 and P = 0.006, respectively). We found qualitatively consistent results using a fully automated spike ripple detector, including comparison with an automated spike detector. We conclude that scalp spike ripple events identify disease and track with seizure risk in this epilepsy population, using both semi-automated and fully automated detection methods, and that this biomarker outperforms analysis of spikes alone in categorizing seizure risk. These data provide evidence that spike ripples are a specific non-invasive biomarker for seizure risk in benign epilepsy with centrotemporal spikes and support future work to evaluate the utility of this biomarker to guide medication trials and tapers in these children and predict seizure risk in other at-risk populations.

Redaction of false high frequency oscillations due to muscle artifact improves specificity to epileptic tissue

OBJECTIVE

High Frequency Oscillations (HFOs) are a promising biomarker of epilepsy. HFOs are typically acquired on intracranial electrodes, but contamination from muscle artifacts is still problematic in HFO analysis. This paper evaluates the effect of myogenic artifacts on intracranial HFO detection and how to remove them.

METHODS

Intracranial EEG was recorded in 31 patients. HFOs were detected for the entire recording using an automated algorithm. When available, simultaneous scalp EEG was used to identify periods of muscle artifact. Those markings were used to train an automated scalp EMG detector and an intracranial EMG detector. Specificity to epileptic tissue was evaluated by comparison with seizure onset zone and resected volume in patients with good outcome.

RESULTS

EMG artifacts are frequent and produce large numbers of false HFOs, especially in the anterior temporal lobe. The scalp and intracranial EMG detectors both had high accuracy. Removing false HFOs improved specificity to epileptic tissue.

CONCLUSIONS

Evaluation of HFOs requires accounting for the effect of muscle artifact. We present two tools that effectively mitigate the effect of muscle artifact on HFOs.

SIGNIFICANCE

Removing muscle artifacts improves the specificity of HFOs to epileptic tissue. Future HFO work should account for this effect, especially when using automated algorithms or when scalp electrodes are not present.

Localization of the Epileptogenic Zone Using High Frequency Oscillations

For patients with drug-resistant focal epilepsy, surgery is the therapy of choice in order to achieve seizure freedom. Epilepsy surgery foremost requires the identification of the epileptogenic zone (EZ), defined as the brain area indispensable for seizure generation. The current gold standard for identification of the EZ is the seizure-onset zone (SOZ). The fact, however that surgical outcomes are unfavorable in 40-50% of well-selected patients, suggests that the SOZ is a suboptimal biomarker of the EZ, and that new biomarkers resulting in better postsurgical outcomes are needed. Research of recent years suggested that high-frequency oscillations (HFOs) are a promising biomarker of the EZ, with a potential to improve surgical success in patients with drug-resistant epilepsy without the need to record seizures. Nonetheless, in order to establish HFOs as a clinical biomarker, the following issues need to be addressed. First, evidence on HFOs as a clinically relevant biomarker stems predominantly from retrospective assessments with visual marking, leading to problems of reproducibility and reliability. Prospective assessments of the use of HFOs for surgery planning using automatic detection of HFOs are needed in order to determine their clinical value. Second, disentangling physiologic from pathologic HFOs is still an unsolved issue. Considering the appearance and the topographic location of presumed physiologic HFOs could be immanent for the interpretation of HFO findings in a clinical context. Third, recording HFOs non-invasively via scalp electroencephalography (EEG) and magnetoencephalography (MEG) is highly desirable, as it would provide us with the possibility to translate the use of HFOs to the scalp in a large number of patients. This article reviews the literature regarding these three issues. The first part of the article focuses on the clinical value of invasively recorded HFOs in localizing the EZ, the detection of HFOs, as well as their separation from physiologic HFOs. The second part of the article focuses on the current state of the literature regarding non-invasively recorded HFOs with emphasis on findings and technical considerations regarding their localization.

Surgical resection of ripple onset predicts outcome in pediatric epilepsy

OBJECTIVE

In patients with medically refractory epilepsy (MRE), interictal ripples (80-250Hz) are observed in large brain areas whose resection may be unnecessary for seizure freedom. This limits their utility as epilepsy biomarkers for surgery. We assessed the spatiotemporal propagation of interictal ripples on intracranial electroencephalography (iEEG) in children with MRE, compared it with the propagation of spikes, identified ripples that initiated propagation (onset-ripples), and evaluated their clinical value as epilepsy biomarkers.

METHODS

Twenty-seven children who underwent epilepsy surgery were studied. We identified propagation sequences of ripples and spikes across multiple iEEG contacts and calculated each ripple or spike latency from the propagation onset. We classified ripples and spikes into categories (ie, onset, spread, and isolated) based on their spatiotemporal characteristics and correlated their mean rate inside and outside resection with outcome (good outcome, Engel 1 versus poor outcome, Engel≥2). We determined, as onset-zone, spread-zone, and isolated-zone, the areas generating the corresponding ripple or spike category and evaluated the predictive value of their resection.

RESULTS

We observed ripple propagation in all patients and spike propagation in 25 patients. Mean rate of onset-ripples inside resection predicted the outcome (odds ratio = 5.37; p = 0.02) and correlated with Engel class (rho = -0.55; p = 0.003). Resection of the onset-ripple-zone was associated with good outcome (p = 0.047). No association was found for the spread-ripple-zone, isolated-ripple-zone, or any spike-zone.

INTERPRETATION

Interictal ripples propagate across iEEG contacts in children with MRE. The association between the onset-ripple-zone resection and good outcome indicates that onset-ripples are promising epilepsy biomarkers, which estimate the epileptogenic tissue better than spread-ripples or onset-spikes. Ann Neurol 2018;84:331-346.

Decision Support System for Seizure Onset Zone Localization Based on Channel Ranking and High-Frequency EEG Activity

Interictal high-frequency oscillations (HFO) are a promising biomarker that can help define the seizure onset zone (SOZ) and predict the surgical outcome after the epilepsy surgery. The utility of HFO in planning the surgery, though, is unclear. Reasons include the variability of the HFO across patients and brain regions and the influence of the sleep-wake cycle, which causes large fluctuations in the ratio between the HFO observed in SOZ and non-SOZ regions. To cope with these limitations, a rank-based solution is proposed to identify the SOZ by using the HFO in multichannel intracranial EEG. A time-varying index of the epileptic susceptibility of the different brain areas is derived from the HFO rate and a support vector machine is applied on this index to identify the SOZ. The solution is trained and tested on separate groups of patients to avoid the use of patient-specific information and provides optimal SOZ prediction using as little as 30 min of recordings per channel (window). Tested on 14 patients with various combinations of seizure type, epilepsy etiology, and SOZ arrangement (172.7 ± 90.1 h/channel per patient and 75.6 ± 23.5 channels/patient, mean ± S.D.), our solution identified the SOZ with 0.92 ± 0.03 accuracy and 0.91 ± 0.03 area under the ROC curve (mean ± S.D.) across patients. For each patient, the window onset time was varied over 72 continuous hours and the prediction of the SOZ remained insensitive to the onset time, thus showing potential for surgery planning.

High-Frequency Oscillations in the Normal Human Brain

OBJECTIVE

High-frequency oscillations (HFOs) are a promising biomarker for the epileptogenic zone. It has not been possible, however, to differentiate physiological from pathological HFOs, and baseline rates of HFO occurrence vary substantially across brain regions. This project establishes region-specific normative values for physiological HFOs and high-frequency activity (HFA).

METHODS

Intracerebral stereo-encephalographic recordings with channels displaying normal physiological activity from nonlesional tissue were selected from 2 tertiary epilepsy centers. Twenty-minute sections from N2/N3 sleep were selected for automatic detection of ripples (80-250Hz), fast ripples (>250Hz), and HFA defined as long-lasting activity > 80Hz. Normative values are provided for 17 brain regions.

RESULTS

A total of 1,171 bipolar channels with normal physiological activity from 71 patients were analyzed. The highest rates of ripples were recorded in the occipital cortex, medial and basal temporal region, transverse temporal gyrus and planum temporale, pre- and postcentral gyri, and medial parietal lobe. The mean rate of fast ripples was very low (0.038/min). Only 5% of channels had a rate > 0.2/min HFA was observed in the medial occipital lobe, pre- and postcentral gyri, transverse temporal gyri and planum temporale, and lateral occipital lobe.

INTERPRETATION

This multicenter atlas is the first to provide region-specific normative values for physiological HFO rates and HFA in common stereotactic space; rates above these can now be considered pathological. Physiological ripples are frequent in eloquent cortex. In contrast, physiological fast ripples are very rare, making fast ripples a good candidate for defining the epileptogenic zone. Ann Neurol 2018;84:374-385.

Physiological and pathological high frequency oscillations in focal epilepsy

Objective

This study investigates high‐frequency oscillations (HFOs; 65–600 Hz) as a biomarker of epileptogenic brain and explores three barriers to their clinical translation: (1) Distinguishing pathological HFOs (pathHFO) from physiological HFOs (physHFO). (2) Classifying tissue under individual electrodes as epileptogenic (3) Reproducing results across laboratories.

Methods

We recorded HFOs using intracranial EEG (iEEG) in 90 patients with focal epilepsy and 11 patients without epilepsy. In nine patients with epilepsy putative physHFOs were induced by cognitive or motor tasks. HFOs were identified using validated detectors. A support vector machine (SVM) using HFO features was developed to classify tissue under individual electrodes as normal or epileptogenic.

Results

There was significant overlap in the amplitude, frequency, and duration distributions for spontaneous physHFO, task induced physHFO, and pathHFO, but the amplitudes of the pathHFO were higher (P < 0.0001). High gamma pathHFO had the strongest association with seizure onset zone (SOZ), and were elevated on SOZ electrodes in 70% of epilepsy patients (P < 0.0001). Failure to resect tissue generating high gamma pathHFO was associated with poor outcomes (P < 0.0001). A SVM classified individual electrodes as epileptogenic with 63.9% sensitivity and 73.7% specificity using SOZ as the target.

Interpretation

A broader range of interictal pathHFO (65–600 Hz) than previously recognized are biomarkers of epileptogenic brain, and are associated with SOZ and surgical outcome. Classification of HFOs into physiological or pathological remains challenging. Classification of tissue under individual electrodes was demonstrated to be feasible. The open source data and algorithms provide a resource for future studies.

Variability in the location of high frequency oscillations during prolonged intracranial EEG recordings

The rate of interictal high frequency oscillations (HFOs) is a promising biomarker of the seizure onset zone, though little is known about its consistency over hours to days. Here we test whether the highest HFO-rate channels are consistent across different 10-min segments of EEG during sleep. An automated HFO detector and blind source separation are applied to nearly 3000 total hours of data from 121 subjects, including 12 control subjects without epilepsy. Although interictal HFOs are significantly correlated with the seizure onset zone, the precise localization is consistent in only 22% of patients. The remaining patients either have one intermittent source (16%), different sources varying over time (45%), or insufficient HFOs (17%). Multiple HFO networks are found in patients with both one and multiple seizure foci. These results indicate that robust HFO interpretation requires prolonged analysis in context with other clinical data, rather than isolated review of short data segments.

Postoperative seizure outcome-guided machine learning for interictal electrocorticography in neocortical epilepsy

The objective of this study was to introduce a new machine learning guided by outcome of resective epilepsy surgery defined as the presence/absence of seizures to improve data mining for interictal pathological activities in neocortical epilepsy. Electrocorticographies for 39 patients with medically intractable neocortical epilepsy were analyzed. We separately analyzed 38 frequencies from 0.9 to 800 Hz including both high-frequency activities and low-frequency activities to select bands related to seizure outcome. An automatic detector using amplitude-duration-number thresholds was used. Interictal electrocorticography data sets of 8 min for each patient were selected. In the first training data set of 20 patients, the automatic detector was optimized to best differentiate the seizure-free group from not-seizure-free-group based on ranks of resection percentages of activities detected using a genetic algorithm. The optimization was validated in a different data set of 19 patients. There were 16 (41%) seizure-free patients. The mean follow-up duration was 21 ± 11 mo (range, 13-44 mo). After validation, frequencies significantly related to seizure outcome were 5.8, 8.4-25, 30, 36, 52, and 75 among low-frequency activities and 108 and 800 Hz among high-frequency activities. Resection for 5.8, 8.4-25, 108, and 800 Hz activities consistently improved seizure outcome. Resection effects of 17-36, 52, and 75 Hz activities on seizure outcome were variable according to thresholds. We developed and validated an automated detector for monitoring interictal pathological and inhibitory/physiological activities in neocortical epilepsy using a data-driven approach through outcome-guided machine learning. NEW & NOTEWORTHY Outcome-guided machine learning based on seizure outcome was used to improve detections for interictal electrocorticographic low- and high-frequency activities. This method resulted in better separation of seizure outcome groups than others reported in the literature. The automatic detector can be trained without human intervention and no prior information. It is based only on objective seizure outcome data without relying on an expert's manual annotations. Using the method, we could find and characterize pathological and inhibitory activities.

Learning-enhanced coupling between ripple oscillations in association cortices and hippocampus

Consolidation of declarative memories requires hippocampal-neocortical communication. Although experimental evidence supports the role of sharp-wave ripples in transferring hippocampal information to the neocortex, the exact cortical destinations and the physiological mechanisms of such transfer are not known. We used a conducting polymer-based conformable microelectrode array (NeuroGrid) to record local field potentials and neural spiking across the dorsal cortical surface of the rat brain, combined with silicon probe recordings in the hippocampus, to identify candidate physiological patterns. Parietal, midline, and prefrontal, but not primary cortical areas, displayed localized ripple (100 to 150 hertz) oscillations during sleep, concurrent with hippocampal ripples. Coupling between hippocampal and neocortical ripples was strengthened during sleep following learning. These findings suggest that ripple-ripple coupling supports hippocampal-association cortical transfer of memory traces.

Dopamine drives left-hemispheric lateralization of neural networks during human speech

Although the concept of left-hemispheric lateralization of neural processes during speech production has been known since the times of Broca, its physiological underpinnings still remain elusive. We sought to assess the modulatory influences of a major neurotransmitter, dopamine, on hemispheric lateralization during real-life speaking using a multimodal analysis of functional MRI, intracranial EEG recordings, and large-scale neural population simulations based on diffusion-weighted MRI. We demonstrate that speech-induced phasic dopamine release into the dorsal striatum and speech motor cortex exerts direct modulation of neuronal activity in these regions and drives left-hemispheric lateralization of speech production network. Dopamine-induced lateralization of functional activity and networks during speaking is not dependent on lateralization of structural nigro-striatal and nigro-motocortical pathways. Our findings provide the first mechanistic explanation for left-hemispheric lateralization of human speech that is due to left-lateralized dopaminergic modulation of brain activity and functional networks.