High-Frequency Oscillations in the Normal Human Brain

Chronobiological Spatial Clusters of Cortical Regions in the Human Brain

PURPOSE

We demonstrate that different regions of the cerebral cortex have different diurnal rhythms of spontaneously occurring high-frequency oscillations (HFOs).

METHODS

High-frequency oscillations were assessed with standard-of-care stereotactic electroencephalography in patients with drug-resistant epilepsy. To ensure generalizability of our findings beyond patients with drug-resistant epilepsy, we excluded stereotactic electroencephalography electrode contacts lying within seizure-onset zones, epileptogenic lesions, having frequent epileptiform activity, and excessive artifact. For each patient, we evaluated twenty-four 5-minute stereotactic electroencephalography epochs, sampled hourly throughout the day, and obtained the HFO rate (number of HFOs/minute) in every stereotactic electroencephalography channel. We analyzed diurnal rhythms of the HFO rates with the cosinor model and clustered neuroanatomic parcels in a standard brain space based on similarity of their cosinor parameters. Finally, we compared overlap among resting-state networks, described in the neuroimaging literature, and chronobiological spatial clusters discovered by us.

RESULTS

We found five clusters that localized predominantly or exclusively to the left perisylvian, left perirolandic and left temporal, right perisylvian and right parietal, right frontal, and right insular-opercular cortices, respectively. These clusters were characterized by similarity of the HFO rates according to the time of the day. Also, these chronobiological spatial clusters preferentially overlapped with specific resting-state networks, particularly default mode network (clusters 1 and 3), frontoparietal network (cluster 1), visual network (cluster 1), and mesial temporal network (cluster 2).

CONCLUSIONS

This is probably the first human study to report clusters of cortical regions with similar diurnal rhythms of electrographic activity. Overlap with resting-state networks attests to their functional significance and has implications for understanding cognitive functions and epilepsy-related mortality.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

SIGNIFICANCE

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

PLAIN LANGUAGE SUMMARY

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

High frequency oscillations in human memory and cognition: a neurophysiological substrate of engrams?

Despite advances in understanding the cellular and molecular processes underlying memory and cognition, and recent successful modulation of cognitive performance in brain disorders, the neurophysiological mechanisms remain underexplored. High frequency oscillations beyond the classic electroencephalogram spectrum have emerged as a potential neural correlate of fundamental cognitive processes. High frequency oscillations are detected in the human mesial temporal lobe and neocortical intracranial recordings spanning gamma/epsilon (60-150 Hz), ripple (80-250 Hz) and higher frequency ranges. Separate from other non-oscillatory activities, these brief electrophysiological oscillations of distinct duration, frequency and amplitude are thought to be generated by coordinated spiking of neuronal ensembles within volumes as small as a single cortical column. Although the exact origins, mechanisms and physiological roles in health and disease remain elusive, they have been associated with human memory consolidation and cognitive processing. Recent studies suggest their involvement in encoding and recall of episodic memory with a possible role in the formation and reactivation of memory traces. High frequency oscillations are detected during encoding, throughout maintenance, and right before recall of remembered items, meeting a basic definition for an engram activity. The temporal coordination of high frequency oscillations reactivated across cortical and subcortical neural networks is ideally suited for integrating multimodal memory representations, which can be replayed and consolidated during states of wakefulness and sleep. High frequency oscillations have been shown to reflect coordinated bursts of neuronal assembly firing and offer a promising substrate for tracking and modulation of the hypothetical electrophysiological engram.

The hemodynamic response to co-occurring interictal epileptiform discharges and high-frequency oscillations localizes the seizure-onset zone

OBJECTIVE

To use intracranial electroencephalography (EEG) to characterize functional magnetic resonance imaging (fMRI) activation maps associated with high-frequency oscillations (HFOs) (80-250 Hz) and examine their proximity to HFO- and seizure-generating tissue.

METHODS

Forty-five patients implanted with intracranial depth electrodes underwent a simultaneous EEG-fMRI study at 3 T. HFOs were detected algorithmically from cleaned EEG and visually confirmed by an experienced electroencephalographer. HFOs that co-occurred with interictal epileptiform discharges (IEDs) were subsequently identified. fMRI activation maps associated with HFOs were generated that occurred either independently of IEDs or within ±200 ms of an IED. For all significant analyses, the Maximum, Second Maximum, and Closest activation clusters were identified, and distances were measured to both the electrodes where the HFOs were observed and the electrodes involved in seizure onset.

RESULTS

We identified 108 distinct groups of HFOs from 45 patients. We found that HFOs with IEDs produced fMRI clusters that were closer to the local field potentials of the corresponding HFOs observed within the EEG than HFOs without IEDs. In addition to the fMRI clusters being closer to the location of the EEG correlate, HFOs with IEDs generated Maximum clusters with greater z-scores and larger volumes than HFOs without IEDs. We also observed that HFOs with IEDs resulted in more discrete activation maps.

SIGNIFICANCE

Intracranial EEG-fMRI can be used to probe the hemodynamic response to HFOs. The hemodynamic response associated with HFOs that co-occur with IEDs better identifies known epileptic tissue than HFOs that occur independently.

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.

Rényi entropy-complexity causality space: a novel neurocomputational tool for detecting scale-free features in EEG/iEEG data

Scale-free brain activity, linked with learning, the integration of different time scales, and the formation of mental models, is correlated with a metastable cognitive basis. The spectral slope, a key aspect of scale-free dynamics, was proposed as a potential indicator to distinguish between different sleep stages. Studies suggest that brain networks maintain a consistent scale-free structure across wakefulness, anesthesia, and recovery. Although differences in anesthetic sensitivity between the sexes are recognized, these variations are not evident in clinical electroencephalographic recordings of the cortex. Recently, changes in the slope of the power law exponent of neural activity were found to correlate with changes in Rényi entropy, an extended concept of Shannon's information entropy. These findings establish quantifiers as a promising tool for the study of scale-free dynamics in the brain. Our study presents a novel visual representation called the Rényi entropy-complexity causality space, which encapsulates complexity, permutation entropy, and the Rényi parameter q. The main goal of this study is to define this space for classical dynamical systems within theoretical bounds. In addition, the study aims to investigate how well different time series mimicking scale-free activity can be discriminated. Finally, this tool is used to detect dynamic features in intracranial electroencephalography (iEEG) signals. To achieve these goals, the study implementse the Bandt and Pompe method for ordinal patterns. In this process, each signal is associated with a probability distribution, and the causal measures of Rényi entropy and complexity are computed based on the parameter q. This method is a valuable tool for analyzing simulated time series. It effectively distinguishes elements of correlated noise and provides a straightforward means of examining differences in behaviors, characteristics, and classifications. For the iEEG experimental data, the REM state showed a greater number of significant sex-based differences, while the supramarginal gyrus region showed the most variation across different modes and analyzes. Exploring scale-free brain activity with this framework could provide valuable insights into cognition and neurological disorders. The results may have implications for understanding differences in brain function between the sexes and their possible relevance to neurological disorders.

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.

Seizure Outcome After Intraoperative Electrocorticography-Tailored Epilepsy Surgery: A Systematic Review and Meta-Analysis

BACKGROUND AND OBJECTIVES

Tailoring epilepsy surgery using intraoperative electrocorticography (ioECoG) has been debated, and modest number of epilepsy surgery centers apply this diagnostic method. We assessed the current evidence to use ioECoG-tailored epilepsy surgery for improving postsurgical outcome.

METHODS

PubMed and Embase were searched for original studies reporting on ≥10 cases who underwent ioECoG-tailored surgery for epilepsy, with a follow-up of at least 6 months. We used a random-effects model to calculate the overall rate of patients achieving favorable seizure outcome (FSO), defined as Engel class I, ILAE class 1, or seizure-free status. Meta-regression was used to investigate potential sources of heterogeneity. We calculated the odds ratio (OR) for estimating variables on FSO:ioECoG vs non-ioECoG-tailored surgery (if included studies contained patients with non-ioECoG-tailored surgery), ioECoG-tailored epilepsy surgery in children vs adults, temporal (TL) vs extratemporal lobe (eTL), MRI-positive vs MRI-negative, and complete vs incomplete resection of tissue that generated interictal epileptiform discharges (IEDs). A Bayesian network meta-analysis was conducted for underlying pathologies. We assessed the evidence certainty using the Grading of Recommendations, Assessment, Development, and Evaluation (GRADE).

RESULTS

Eighty-three studies (82 observational studies, 1 trial) comprising 3,631 patients with ioECoG-tailored surgery were included. The overall pooled rate of patients who attained FSO after ioECoG-tailored surgery was 74% (95% CI 71-77) with significant heterogeneity, which was predominantly attributed to pathologies and seizure outcome classifications. Twenty-two studies contained non-ioECoG-tailored surgeries. IoECoG-tailored surgeries reached a higher rate of FSO than non-ioECoG-tailored surgeries (OR 2.10 [95% CI 1.37-3.24]; p < 0.01; very low certainty). Complete resection of tissue that displayed IEDs in ioECoG predicted FSO better compared with incomplete resection (OR 3.04 [1.76-5.25]; p < 0.01; low certainty). We found insignificant difference in FSO after ioECoG-tailored surgery in children vs adults, TL vs eTL, or MRI-positive vs MRI-negative. The network meta-analysis showed that the odds of FSO was lower for malformations of cortical development than for tumors (OR 0.47 95% credible interval 0.25-0.87).

DISCUSSION

Although limited by low-quality evidence, our meta-analysis shows a relatively good surgical outcome (74% FSO) after epilepsy surgery with ioECoG, especially in tumors, with better outcome for ioECoG-tailored surgeries in studies describing both and better outcome after complete removal of IED areas.

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.

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.

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.

Recent advances in clinical electroencephalography

PURPOSE OF REVIEW

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

RECENT FINDINGS

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

SUMMARY

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

Scalp and hippocampal sleep correlates of memory function in drug-resistant temporal lobe epilepsy

Seminal animal studies demonstrated the role of sleep oscillations such as cortical slow waves, thalamocortical spindles, and hippocampal ripples in memory consolidation. In humans, whether ripples are involved in sleep-related memory processes is less clear. Here, we explored the interactions between sleep oscillations (measured as traits) and general episodic memory abilities in 26 adults with drug-resistant temporal lobe epilepsy who performed scalp-intracranial electroencephalographic recordings and neuropsychological testing, including two analogous hippocampal-dependent verbal and nonverbal memory tasks. We explored the relationships between hemispheric scalp (spindles, slow waves) and hippocampal physiological and pathological oscillations (spindles, slow waves, ripples, and epileptic spikes) and material-specific memory function. To differentiate physiological from pathological ripples, we used multiple unbiased data-driven clustering approaches. At the individual level, we found material-specific cerebral lateralization effects (left-verbal memory, right-nonverbal memory) for all scalp spindles (rs > 0.51, ps < 0.01) and fast spindles (rs > 0.61, ps < 0.002). Hippocampal epileptic spikes and short pathological ripples, but not physiological oscillations, were negatively (rs > -0.59, ps < 0.01) associated with verbal learning and retention scores, with left lateralizing and antero-posterior effects. However, data-driven clustering failed to separate the ripple events into defined clusters. Correlation analyses with the resulting clusters revealed no meaningful or significant associations with the memory scores. Our results corroborate the role of scalp spindles in memory processes in patients with drug-resistant temporal lobe epilepsy. Yet, physiological and pathological ripples were not separable when using data-driven clustering, and thus our findings do not provide support for a role of sleep ripples as trait-like characteristics of general memory abilities in epilepsy.

High frequency oscillation network dynamics predict outcome in non-palliative epilepsy surgery

High frequency oscillations are a promising biomarker of outcome in intractable epilepsy. Prior high frequency oscillation work focused on counting high frequency oscillations on individual channels, and it is still unclear how to translate those results into clinical care. We show that high frequency oscillations arise as network discharges that have valuable properties as predictive biomarkers. Here, we develop a tool to predict patient outcome before surgical resection is performed, based on only prospective information. In addition to determining high frequency oscillation rate on every channel, we performed a correlational analysis to evaluate the functional connectivity of high frequency oscillations in 28 patients with intracranial electrodes. We found that high frequency oscillations were often not solitary events on a single channel, but part of a local network discharge. Eigenvector and outcloseness centrality were used to rank channel importance within the connectivity network, then used to compare patient outcome by comparison with the seizure onset zone or a proportion within the proposed resected channels (critical resection percentage). Combining the knowledge of each patient's seizure onset zone resection plan along with our computed high frequency oscillation network centralities and high frequency oscillation rate, we develop a Naïve Bayes model that predicts outcome (positive predictive value: 100%) better than predicting based upon fully resecting the seizure onset zone (positive predictive value: 71%). Surgical margins had a large effect on outcomes: non-palliative patients in whom most of the seizure onset zone was resected ('definitive surgery', ≥ 80% resected) had predictable outcomes, whereas palliative surgeries (<80% resected) were not predictable. These results suggest that the addition of network properties of high frequency oscillations is more accurate in predicting patient outcome than seizure onset zone alone in patients with most of the seizure onset zone removed and offer great promise for informing clinical decisions in surgery for refractory epilepsy.

Learn how to interpret and use intracranial EEG findings

Epilepsy surgery is the therapy of choice for many patients with drug-resistant focal epilepsy. Recognizing and describing ictal and interictal patterns with intracranial electroencephalography (EEG) recordings is important in order to most efficiently leverage advantages of this technique to accurately delineate the seizure-onset zone before undergoing surgery. In this seminar in epileptology, we address learning objective "1.4.11 Recognize and describe ictal and interictal patterns with intracranial recordings" of the International League against Epilepsy curriculum for epileptologists. We will review principal considerations of the implantation planning, summarize the literature for the most relevant ictal and interictal EEG patterns within and beyond the Berger frequency spectrum, review invasive stimulation for seizure and functional mapping, discuss caveats in the interpretation of intracranial EEG findings, provide an overview on special considerations in children and in subdural grids/strips, and review available quantitative/signal analysis approaches. To be as practically oriented as possible, we will provide a mini atlas of the most frequent EEG patterns, highlight pearls for its not infrequently challenging interpretation, and conclude with two illustrative case examples. This article shall serve as a useful learning resource for trainees in clinical neurophysiology/epileptology by providing a basic understanding on the concepts of invasive intracranial EEG.

Magnetoencephalography for Epilepsy Presurgical Evaluation

PURPOSE OF THE REVIEW

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

RECENT FINDINGS

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

Mouse model of focal cortical dysplasia type II generates a wide spectrum of high-frequency activities

High-frequency oscillations (HFOs) represent an electrographic biomarker of endogenous epileptogenicity and seizure-generating tissue that proved clinically useful in presurgical planning and delineating the resection area. In the neocortex, the clinical observations on HFOs are not sufficiently supported by experimental studies stemming from a lack of realistic neocortical epilepsy models that could provide an explanation of the pathophysiological substrates of neocortical HFOs. In this study, we explored pathological epileptiform network phenomena, particularly HFOs, in a highly realistic murine model of neocortical epilepsy due to focal cortical dysplasia (FCD) type II. FCD was induced in mice by the expression of the human pathogenic mTOR gene mutation during embryonic stages of brain development. Electrographic recordings from multiple cortical regions in freely moving animals with FCD and epilepsy demonstrated that the FCD lesion generates HFOs from all frequency ranges, i.e., gamma, ripples, and fast ripples up to 800 Hz. Gamma-ripples were recorded almost exclusively in FCD animals, while fast ripples occurred in controls as well, although at a lower rate. Gamma-ripple activity is particularly valuable for localizing the FCD lesion, surpassing the utility of fast ripples that were also observed in control animals, although at significantly lower rates. Propagating HFOs occurred outside the FCD, and the contralateral cortex also generated HFOs independently of the FCD, pointing to a wider FCD network dysfunction. Optogenetic activation of neurons carrying mTOR mutation and expressing Channelrhodopsin-2 evoked fast ripple oscillations that displayed spectral and morphological profiles analogous to spontaneous oscillations. This study brings experimental evidence that FCD type II generates pathological HFOs across all frequency bands and provides information about the spatiotemporal properties of each HFO subtype in FCD. The study shows that mutated neurons represent a functionally interconnected and active component of the FCD network, as they can induce interictal epileptiform phenomena and HFOs.

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

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

Machine learning-based classification of physiological and pathological high-frequency oscillations recorded by stereoelectroencephalography

OBJECTIVE

High-frequency oscillations (HFOs) are an efficient indicator to locate the epileptogenic zone (EZ). However, physiological HFOs produced in the normal brain region may interfere with EZ localization. The present study aimed to build a machine learning-based classifier to distinguish the properties of each HFO event based on features in different domains.

METHODS

HFOs were detected in focal epilepsy patients from two different hospitals who underwent stereoelectroencephalography and subsequent resection surgery. Subsequently, 37 features in four different domains (time, frequency and time-frequency, entropy-based and nonlinear) were extracted for each HFO. After extraction, a fast correlation-based filter (FCBF) algorithm was applied for feature selection. The machine learning classifier was trained on the feature matrix with and without FCBF and then tested on the data set from patients in another hospital.

RESULTS

A dataset was compiled, consisting of 89,844 pathological HFOs and 23,613 physiological HFOs from 17 patients assigned to the training dataset. Additionally, 12,695 pathological HFOs and 5,599 physiological HFOs from 9 patients were assigned to the testing dataset. Four features (ripple band power, arithmetic mean, Petrosian fractal dimension and zero crossings) were obtained for classifier training after FCBF. The classifier showed an area under the curve (AUC) of 0.95/0.98 for FCBF/no FCBF features in the training dataset and AUC of 0.82/0.90 for FCBF/no FCBF features in the testing dataset. Our findings indicated that the classifier utilizing all features demonstrated superior performance compared to the one relying on FCBF-processed features.

CONCLUSION

Our classifier could reliably differentiate pathological HFOs from physiological ones, which could promote the development of HFOs in EZ localization.

Wake slow waves in focal human epilepsy impact network activity and cognition

Slow waves of neuronal activity are a fundamental component of sleep that are proposed to have homeostatic and restorative functions. Despite this, their interaction with pathology is unclear and there is only indirect evidence of their presence during wakefulness. Using intracortical recordings from the temporal lobe of 25 patients with epilepsy, we demonstrate the existence of local wake slow waves (LoWS) with key features of sleep slow waves, including a down-state of neuronal firing. Consistent with a reduction in neuronal activity, LoWS were associated with slowed cognitive processing. However, we also found that LoWS showed signatures of a homeostatic relationship with interictal epileptiform discharges (IEDs): exhibiting progressive adaptation during the build-up of network excitability before an IED and reducing the impact of subsequent IEDs on network excitability. We therefore propose an epilepsy homeostasis hypothesis: that slow waves in epilepsy reduce aberrant activity at the price of transient cognitive impairment.

The role of quantitative markers in surgical prognostication after stereoelectroencephalography

OBJECTIVE

Stereoelectroencephalography (SEEG) is the reference method in the presurgical exploration of drug-resistant focal epilepsy. However, prognosticating surgery on an individual level is difficult. A quantified estimation of the most epileptogenic regions by searching for relevant biomarkers can be proposed for this purpose. We investigated the performances of ictal (Epileptogenicity Index, EI; Connectivity EI, cEI), interictal (spikes, high-frequency oscillations, HFO [80-300 Hz]; Spikes × HFO), and combined (Spikes × EI; Spikes × cEI) biomarkers in predicting surgical outcome and searched for prognostic factors based on SEEG-signal quantification.

METHODS

Fifty-three patients operated on following SEEG were included. We compared, using precision-recall, the epileptogenic zone quantified using different biomarkers (EZq ) against the visual analysis (EZC ). Correlations between the EZ resection rates or the EZ extent and surgical prognosis were analyzed.

RESULTS

EI and Spikes × EI showed the best precision against EZc (0.74; 0.70), followed by Spikes × cEI and cEI, whereas interictal markers showed lower precision. The EZ resection rates were greater in seizure-free than in non-seizure-free patients for the EZ defined by ictal biomarkers and were correlated with the outcome for EI and Spikes × EI. No such correlation was found for interictal markers. The extent of the quantified EZ did not correlate with the prognosis.

INTERPRETATION

Ictal or combined ictal-interictal markers overperformed the interictal markers both for detecting the EZ and predicting seizure freedom. Combining ictal and interictal epileptogenicity markers improves detection accuracy. Resection rates of the quantified EZ using ictal markers were the only statistically significant determinants for surgical prognosis.

Rapid eye movement sleep affects interictal epileptic activity differently in mesiotemporal and neocortical areas

OBJECTIVE

Rapid eye movement (REM) sleep reduces the rate and extent of interictal epileptiform discharges (IEDs). Breakthrough epileptic activity during REM sleep is therefore thought to best localize the seizure onset zone (SOZ). We utilized polysomnography combined with direct cortical recordings to investigate the influences of anatomical locations and the time of night on the suppressive effect of REM sleep on IEDs.

METHODS

Forty consecutive patients with drug-resistant focal epilepsy underwent combined polysomnography and stereo-electroencephalography during presurgical evaluation. Ten-minute interictal epochs were selected 2 h prior to sleep onset (wakefulness), and from the first and second half of the night during non-REM (NREM) sleep and REM sleep. IEDs were detected automatically across all channels. Anatomic localization, time of night, and channel type (within or outside the SOZ) were tested as modulating factors.

RESULTS

Relative to wakefulness, there was a suppression of IEDs by REM sleep in neocortical regions (median = -27.6%), whereas mesiotemporal regions showed an increase in IEDs (19.1%, p = .01, d = .39). This effect was reversed when comparing the regional suppression of IEDs by REM sleep relative to NREM sleep (-35.1% in neocortical, -58.7% in mesiotemporal, p < .001, d = .39). Across all patients, no clinically relevant novel IED regions were observed in REM sleep versus NREM or wakefulness based on our predetermined thresholds (4 IEDs/min in REM, 0 IEDs/min in NREM and wakefulness). Finally, there was a reduction in IEDs in late (NREM: 1.08/min, REM: .61/min) compared to early sleep (NREM: 1.22/min, REM: .69/min) for both NREM (p < .001, d = .21) and REM (p = .04, d = .14).

SIGNIFICANCE

Our results demonstrate a spatiotemporal effect of IED suppression by REM sleep relative to wakefulness in neocortical but not mesiotemporal regions, and in late versus early sleep. This suggests the importance of considering sleep stage interactions and the potential influences of anatomical locations when using IEDs to define the epileptic focus.

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

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

SIGNIFICANCE

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

Has recording of seizures become obsolete?

Some patients with medically intractable epilepsy are considered for surgical treatment. In some surgical candidates, the investigation includes the placement of intracerebral electrodes and long-term monitoring to find the region of seizure onset. This region is the primary determinant of the surgical resection but about one-third of patients are not offered surgery after electrode implantation and among those operated only about 55% are seizure free after five years. This paper discusses why the primary reliance on the seizure onset maybe suboptimal and may be in part responsible for the relatively low surgical success rate. It also proposes to consider some interictal markers that may have advantages over seizure onset and may be easier to obtain.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Fast ripples reflect increased excitability that primes epileptiform spikes

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

An evolutionary conserved division-of-labor between archicortical and neocortical ripples organizes information transfer during sleep

Systems-level memory consolidation during sleep depends on the temporally precise interplay between cardinal sleep oscillations. Specifically, hippocampal ripples constitute a key substrate of the hippocampal-neocortical dialogue underlying memory formation. Recently, it became evident that ripples are not unique to archicortex, but constitute a wide-spread neocortical phenomenon. To date, little is known about the morphological similarities between archi- and neocortical ripples. Moreover, it remains undetermined if neocortical ripples fulfill distinct functional roles. Leveraging intracranial recordings from the human medial temporal lobe (MTL) and neocortex during sleep, our results reveal region-specific functional specializations, albeit a near-uniform morphology. While MTL ripples synchronize the memory network to trigger directional MTL-to-neocortical information flow, neocortical ripples reduce information flow to minimize interference. At the population level, MTL ripples confined population dynamics to a low-dimensional subspace, while neocortical ripples diversified the population response; thus, constituting an effective mechanism to functionally uncouple the MTL-neocortical network. Critically, we replicated the key findings in rodents, where the same division-of-labor between archi- and neocortical ripples was evident. In sum, these results uncover an evolutionary preserved mechanism where the precisely coordinated interplay between MTL and neocortical ripples temporally segregates MTL information transfer from subsequent neocortical processing during sleep.

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

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

SIGNIFICANCE

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

Quantitative approaches to guide epilepsy surgery from intracranial EEG

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

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

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

Fast ripples reflect increased excitability that primes epileptiform spikes

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

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

Objective

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

Methods

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

Results

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

Conclusions

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

Significance

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

HIGHLIGHTS

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

Identification and classification of pathology and artifacts for human intracranial cognitive research

Intracranial electroencephalography (iEEG) presents a unique opportunity to extend human neuroscientific understanding. However, typically iEEG is collected from patients diagnosed with focal drug-resistant epilepsy (DRE) and contains transient bursts of pathological activity. This activity disrupts performances on cognitive tasks and can distort findings from human neurophysiology studies. In addition to manual marking by a trained expert, numerous IED detectors have been developed to identify these pathological events. Even so, the versatility and usefulness of these detectors is limited by training on small datasets, incomplete performance metrics, and lack of generalizability to iEEG. Here, we employed a large annotated public iEEG dataset from two institutions to train a random forest classifier (RFC) to distinguish data segments as either 'non-cerebral artifact' (n = 73,902), 'pathological activity' (n = 67,797), or 'physiological activity' (n = 151,290). We found our model performed with an accuracy of 0.941, specificity of 0.950, sensitivity of 0.908, precision of 0.911, and F1 score of 0.910, averaged across all three event types. We extended the generalizability of our model to continuous bipolar data collected in a task-state at a different institution with a lower sampling rate and found our model performed with an accuracy of 0.789, specificity of 0.806, and sensitivity of 0.742, averaged across all three event types. Additionally, we created a custom graphical user interface to implement our classifier and enhance usability.

Chaotic dynamics of the Hénon map and neuronal input-output: A comparison with neurophysiological data

In this study, the Hénon map was analyzed using quantifiers from information theory in order to compare its dynamics to experimental data from brain regions known to exhibit chaotic behavior. The goal was to investigate the potential of the Hénon map as a model for replicating chaotic brain dynamics in the treatment of Parkinson's and epilepsy patients. The dynamic properties of the Hénon map were compared with data from the subthalamic nucleus, the medial frontal cortex, and a q-DG model of neuronal input-output with easy numerical implementation to simulate the local behavior of a population. Using information theory tools, Shannon entropy, statistical complexity, and Fisher's information were analyzed, taking into account the causality of the time series. For this purpose, different windows over the time series were considered. The findings revealed that neither the Hénon map nor the q-DG model could perfectly replicate the dynamics of the brain regions studied. However, with careful consideration of the parameters, scales, and sampling used, they were able to model some characteristics of neural activity. According to these results, normal neural dynamics in the subthalamic nucleus region may present a more complex spectrum within the complexity-entropy causality plane that cannot be represented by chaotic models alone. The dynamic behavior observed in these systems using these tools is highly dependent on the studied temporal scale. As the size of the sample studied increases, the dynamics of the Hénon map become increasingly different from those of biological and artificial neural systems.

Local neurodynamics as a signature of cortical areas: new insights from sleep

Sleep crucial for the animal survival is accompanied by huge changes in neuronal electrical activity over time, the neurodynamics. Here, drawing on intracranial stereo-electroencephalographic (sEEG) recordings from the Montreal Neurological Institute (MNI), we analyzed local neurodynamics in the waking state at rest and during the N2, N3, and rapid eye movement (REM) sleep phases. Higuchi fractal dimension (HFD)-a measure of signal complexity-was studied as a feature of the local neurodynamics of the primary motor (M1), somatosensory (S1), and auditory (A1) cortices. The key working hypothesis, that the relationships between local neurodynamics preserve in all sleep phases despite the neurodynamics complexity reduces in sleep compared with wakefulness, was supported by the results. In fact, while HFD awake > REM > N2 > N3 (P < 0.001 consistently), HFD in M1 > S1 > A1 in awake and all sleep stages (P < 0.05 consistently). Also power spectral density was studied for consistency with previous investigations. Meaningfully, we found a local specificity of neurodynamics, well quantified by the fractal dimension, expressed in wakefulness and during sleep. We reinforce the idea that neurodynamic may become a new criterion for cortical parcellation, prospectively improving the understanding and ability of compensatory interventions for behavioral disorders.

A Subpopulation of Spikes Predicts Successful Epilepsy Surgery Outcome

OBJECTIVE

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

METHODS

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

RESULTS

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

INTERPRETATION

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

SIGNIFICANCE

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

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.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Hippocampal area CA2 controls seizure dynamics, interictal EEG abnormalities and social comorbidity in mouse models of temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is characterized by spontaneous recurrent seizures, abnormal activity between seizures, and impaired behavior. CA2 pyramidal neurons (PNs) are potentially important because inhibiting them with a chemogenetic approach reduces seizure frequency in a mouse model of TLE. However, whether seizures could be stopped by timing inhibition just as a seizure begins is unclear. Furthermore, whether inhibition would reduce the cortical and motor manifestations of seizures are not clear. Finally, whether interictal EEG abnormalities and TLE comorbidities would be improved are unknown. Therefore, real-time optogenetic silencing of CA2 PNs during seizures, interictal activity and behavior were studied in 2 mouse models of TLE. CA2 silencing significantly reduced seizure duration and time spent in convulsive behavior. Interictal spikes and high frequency oscillations were significantly reduced, and social behavior was improved. Therefore, brief focal silencing of CA2 PNs reduces seizures, their propagation, and convulsive manifestations, improves interictal EEG, and ameliorates social comorbidities.

HIGHLIGHTS

Real-time CA2 silencing at the onset of seizures reduces seizure durationWhen CA2 silencing reduces seizure activity in hippocampus it also reduces cortical seizure activity and convulsive manifestations of seizuresInterictal spikes and high frequency oscillations are reduced by real-time CA2 silencingReal-time CA2 silencing of high frequency oscillations (>250Hz) rescues social memory deficits of chronic epileptic mice.

Graph theoretical measures of fast ripple networks improve the accuracy of post-operative seizure outcome prediction

Fast ripples (FR) are a biomarker of epileptogenic brain, but when larger portions of FR generating regions are resected seizure freedom is not always achieved. To evaluate and improve the diagnostic accuracy of FR resection for predicting seizure freedom we compared the FR resection ratio (RR) with FR network graph theoretical measures. In 23 patients FR were semi-automatically detected and quantified in stereo EEG recordings during sleep. MRI normalization and co-registration localized contacts and relation to resection margins. The number of FR, and graph theoretical measures, which were spatial (i.e., FR rate-distance radius) or temporal correlational (i.e., FR mutual information), were compared with the resection margins and with seizure outcome We found that the FR RR did not correlate with seizure-outcome (p > 0.05). In contrast, the FR rate-distance radius resected difference and the FR MI mean characteristic path length RR did correlate with seizure-outcome (p < 0.05). Retesting of positive FR RR patients using either FR rate-distance radius resected difference or the FR MI mean characteristic path length RR reduced seizure-free misclassifications from 44 to 22% and 17%, respectively. These results indicate that graph theoretical measures of FR networks can improve the diagnostic accuracy of the resection of FR events for predicting seizure freedom.

High-frequency oscillations (250-500 Hz) in animal models of Alzheimer's disease and two animal models of epilepsy

OBJECTIVE

To test the hypothesis that high-frequency oscillations (HFOs) between 250 and 500 Hz occur in mouse models of Alzheimer's disease (AD) and thus are not unique to epilepsy.

METHODS

Experiments were conducted in three mouse models of AD: Tg2576 mice that simulate a form of familial AD, presenilin 2 knock-out (PS2KO) mice, and the Ts65Dn model of Down's syndrome. We recorded HFOs using wideband (0.1-500 Hz, 2 kHz) intra-hippocampal and cortical surface electroencephalography (EEG) at 1 month until 24 months of age during wakefulness, slow wave sleep (SWS), and rapid eye movement (REM) sleep. In addition, interictal spikes (IISs) and seizures were analyzed for the possible presence of HFOs. Comparisons were made to the intra-hippocampal kainic acid and pilocarpine models of epilepsy.

RESULTS

We describe for the first time that hippocampal and cortical HFOs are a new EEG abnormality in AD mouse models. HFOs occurred in all transgenic mice but no controls. They were also detectable as early as 1 month of age and prior to amyloid beta plaque neuropathology. HFOs were most frequent during SWS (vs REM sleep or wakefulness). Notably, HFOs in the AD and epilepsy models were indistinguishable in both spectral frequency and duration. HFOs also occurred during IISs and seizures in the AD models, although with altered spectral properties compared to isolated HFOs.

SIGNIFICANCE

Our data demonstrate that HFOs, an epilepsy biomarker with high translational value, are not unique to epilepsy and thus not disease specific. Our findings also strengthen the idea of hyperexcitability in AD and its significant overlap with epilepsy. HFOs in AD mouse models may serve as an EEG biomarker, which is detectable from the scalp and thus amenable to noninvasive detection in people at risk for AD.

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.

Automatic Localization of Seizure Onset Zone Based on Multi-Epileptogenic Biomarkers Analysis of Single-Contact from Interictal SEEG

Successful surgery on drug-resistant epilepsy patients (DRE) needs precise localization of the seizure onset zone (SOZ). Previous studies analyzing this issue still face limitations, such as inadequate analysis of features, low sensitivity and limited generality. Our study proposed an innovative and effective SOZ localization method based on multiple epileptogenic biomarkers (spike and HFOs), and analysis of single-contact (MEBM-SC) to address the above problems. We extracted contacts epileptic features from signal distributions and signal energy based on machine learning and end-to-end deep learning. Among them, a normalized pathological ripple rate was designed to reduce the disturbance of physiological ripple and enhance the performance of SOZ localization. Then, a feature selection algorithm based on Shapley value and hypothetical testing (ShapHT+) was used to limit interference from irrelevant features. Moreover, an attention mechanism and a focal loss algorithm were used on the classifier to learn significant features and overcome the unbalance of SOZ/nSOZ contacts. Finally, we provided an SOZ prediction and visualization on magnetic resonance imaging (MRI). Ten patients with DRE were selected to verify our method. The experiment performed cross-validation and revealed that MEBM-SC obtains higher sensitivity. Additionally, the spike has better sensitivity while HFOs have better specificity, and the combination of these biomarkers can achieve the best performance. The study confirmed that MEBM-SC can increase the sensitivity and accuracy of SOZ localization and help clinicians to perform a precise and reliable preoperative evaluation based on interictal SEEG.