Bimodal coupling of ripples and slower oscillations during sleep in patients with focal epilepsy

A universal method for seizure onset zone localization in focal epilepsy using standard deviation of spike amplitude

BACKGROUND

Precisely localizing the seizure onset zone (SOZ) is critical for focal epilepsy surgery. Existing methods mainly focus on high-frequency activities in stereo-electroencephalography, but often fail when seizures are not driven by high-frequency activities. Recognized as biomarkers of epileptic seizures, ictal spikes in SOZ induce epileptiform discharges in other brain regions. Based on this understanding, we aim to develop a universal algorithm to localize SOZ and investigate how ictal spikes within the SOZ induce seizures.

METHODS

We proposed a novel metric called standard deviation of spike amplitude (SDSA) and utilized channel-averaged SDSA to describe seizure processes and detect seizures. By integrating SDSA values in specific intervals, the score for each channel located within SOZ was calculated. Channels with high SOZ scores were clustered as SOZ. The localization accuracy was asserted using area under the receiver operating characteristic (ROC) curve. Further, we analyzed early ictal signals from SOZ channels and investigated factors influencing their duration to reveal the seizure inducing conditions.

RESULTS

We analyzed data from 15 patients with focal epilepsy. The channel-averaged SDSA successfully detected all 28 seizures without false alarms. Using SDSA integration, we achieved precise SOZ localization with an average area under ROC curve (AUC) of 0.96, significantly outperforming previous methods based on high-frequency activities. Further, we discovered that energy of ictal spikes in SOZ was concentrated at a specific frequency distributed in [6, 12 Hz]. Additionally, we found that the higher the energy per second in this frequency band, the faster ictal spikes could induce seizures.

CONCLUSION

The SDSA metric offered precise SOZ localization with robustness and low computational cost, making it suitable for clinical practice. By studying the propagation patterns of ictal spikes between the SOZ and non-SOZ, we suggest that ictal spikes from SOZ need to accumulate energy at a specific central frequency to induce epileptic spikes in non-SOZ, which may have significant implications for understanding the seizure onset pattern.

Self-Supervised Data-Driven Approach Defines Pathological High-Frequency Oscillations in Human

Objective

Interictal high-frequency oscillations (HFOs) are a promising neurophysiological biomarker of the epileptogenic zone (EZ). However, objective criteria for distinguishing pathological from physiological HFOs remain elusive, hindering clinical application. We investigated whether the distinct mechanisms underlying pathological and physiological HFOs are encapsulated in their signal morphology in intracranial EEG (iEEG) recordings and whether this mechanism-driven distinction could be simulated by a deep generative model.

Methods

In a retrospective cohort of 185 epilepsy patients who underwent iEEG monitoring, we analyzed 686,410 HFOs across 18,265 brain contacts. To learn morphological characteristics, each event was transformed into a time-frequency plot and input into a variational autoencoder. We characterized latent space clusters containing morphologically defined putative pathological HFOs (mpHFOs) using interpretability analysis, including latent space disentanglement and time-domain perturbation.

Results

mpHFOs showed strong associations with expert-defined spikes and were predominantly located within the seizure onset zone (SOZ). Discovered novel pathological features included high power in the gamma (30-80 Hz) and ripple (>80 Hz) bands centered on the event. These characteristics were consistent across multiple variables, including institution, electrode type, and patient demographics. Predicting 12-month postoperative seizure outcomes using the resection ratio of mpHFOs outperformed unclassified HFOs (F1=0.72 vs. 0.68) and matched current clinical standards using SOZ resection (F1=0.74). Combining mpHFO data with demographic and SOZ resection status further improved prediction accuracy (F1=0.83).

Interpretation

Our data-driven approach yielded a novel, explainable definition of pathological HFOs, which has the potential to further enhance the clinical use of HFOs for EZ delineation.

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

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

Clinical Utility of Sleep Recordings During Presurgical Epilepsy Evaluation With Stereo-Electroencephalography: A Systematic Review

SUMMARY

Although the role of sleep in modulating epileptic activity is well established, many epileptologists overlook the significance of considering sleep during presurgical epilepsy evaluations in cases of drug-resistant epilepsy. Here, we conducted a comprehensive literature review from January 2000 to May 2023 using the PubMed electronic database and compiled evidence to highlight the need to revise the current clinical approach. All articles were assessed for eligibility by two independent reviewers. Our aim was to shed light on the clinical value of incorporating sleep monitoring into presurgical evaluations with stereo-electroencephalography. We present the latest developments on the important bidirectional interactions between sleep and various forms of epileptic activity observed in stereo-electroencephalography recordings. Specifically, epileptic activity is modulated by different sleep stages, peaking in non-rapid eye movement sleep, while being suppressed in rapid eye movement sleep. However, this modulation can vary across different brain regions, underlining the need to account for sleep to accurately pinpoint the epileptogenic zone during presurgical assessments. Finally, we offer practical solutions, such as automated sleep scoring algorithms using stereo-electroencephalography data alone, to seamlessly integrate sleep monitoring into routine clinical practice. It is hoped that this review will provide clinicians with a readily accessible roadmap to the latest evidence concerning the clinical utility of sleep monitoring in the context of stereo-electroencephalography and aid the development of therapeutic and diagnostic strategies to improve patient surgical outcomes.

Simulated resections and RNS placement can optimize post-operative seizure outcomes when guided by fast ripple networks

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

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.

Pathological neurons generate ripples at the UP-DOWN transition disrupting information transfer

OBJECTIVE

To confirm and investigate why pathological high-frequency oscillations (pHFOs), including ripples (80-200 Hz) and fast ripples (200-600 Hz), are generated during the UP-DOWN transition of the slow wave and if information transmission mediated by ripple temporal coupling is disrupted in the seizure-onset zone (SOZ).

METHODS

We isolated 217 total units from 175.95 intracranial electroencephalography (iEEG) contact-hours of synchronized macro- and microelectrode recordings from 6 patients. Sleep slow oscillation (.1-2 Hz) epochs were identified in the iEEG recording. iEEG HFOs that occurred superimposed on the slow wave were transformed to phasors and adjusted by the phase of maximum firing in nearby units (i.e., maximum UP). We tested whether, in the SOZ, HFOs and associated action potentials (APs) occur more often at the UP-DOWN transition. We also examined ripple temporal correlations using cross-correlograms.

RESULTS

At the group level in the SOZ, HFO and HFO-associated AP probability was highest during the UP-DOWN transition of slow wave excitability (p < < .001). In the non-SOZ, HFO and HFO-associated AP was highest during the DOWN-UP transition (p < < .001). At the unit level in the SOZ, 15.6% and 20% of units exhibited more robust firing during ripples (Cohen's d = .11-.83) and fast ripples (d = .36-.90) at the UP-DOWN transition (p < .05 f.d.r. corrected), respectively. By comparison, also in the SOZ, 6.6% (d = .14-.30) and 8.5% (d = .33-.41) of units had significantly less firing during ripples and fast ripples at the UP-DOWN transition, respectively. Additional data shows that ripple and fast ripple temporal correlations, involving global slow waves, between the hippocampus, entorhinal cortex, and parahippocampal gyrus were reduced by >50% in the SOZ compared to the non-SOZ (N = 3).

SIGNIFICANCE

The UP-DOWN transition of slow wave excitability facilitates the activation of pathological neurons to generate pHFOs. Ripple temporal correlations across brain regions may be important in memory consolidation and are disrupted in the SOZ, perhaps by pHFO generation.

Widespread slow oscillations support interictal epileptiform discharge networks in focal epilepsy

Interictal epileptiform discharges (IEDs) often co-occur across spatially-separated cortical regions, forming IED networks. However, the factors prompting IED propagation remain unelucidated. We hypothesized that slow oscillations (SOs) might facilitate IED propagation. Here, the amplitude and phase synchronization of SOs preceding propagating and non-propagating IEDs were compared in 22 patients with focal epilepsy undergoing intracranial electroencephalography (EEG) evaluation. Intracranial channels were categorized into the irritative zone (IZ) and normal zone (NOZ) regarding the presence of IEDs. During wakefulness, we found that pre-IED SOs within the IZ exhibited higher amplitudes for propagating IEDs than non-propagating IEDs (delta band: p = 0.001, theta band: p < 0.001). This increase in SOs was also concurrently observed in the NOZ (delta band: p = 0.04). Similarly, the inter-channel phase synchronization of SOs prior to propagating IEDs was higher than those preceding non-propagating IEDs in the IZ (delta band: p = 0.04). Through sliding window analysis, we observed that SOs preceding propagating IEDs progressively increased in amplitude and phase synchronization, while those preceding non-propagating IEDs remained relatively stable. Significant differences in amplitude occurred approximately 1150 ms before IEDs. During non-rapid eye movement (NREM) sleep, SOs on scalp recordings also showed higher amplitudes before intracranial propagating IEDs than before non-propagating IEDs (delta band: p = 0.006). Furthermore, the analysis of IED density around sleep SOs revealed that only high-amplitude sleep SOs demonstrated correlation with IED propagation. Overall, our study highlights that transient but widely distributed SOs are associated with IED propagation as well as generation in focal epilepsy during sleep and wakefulness, providing new insight into the EEG substrate supporting IED networks.

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.

Pathological neurons generate ripples at the UP-DOWN transition disrupting information transfer

Objective

To confirm and investigate why pathological HFOs (pHFOs), including Ripples [80-200 Hz] and fast ripples [200-600 Hz], are generated during the UP-DOWN transition of the slow wave and if pHFOs interfere with information transmission.

Methods

We isolated 217 total units from 175.95 iEEG contact-hours of synchronized macro- and microelectrode recordings from 6 patients. Sleep slow oscillation (0.1-2 Hz) epochs were identified in the iEEG recording. iEEG HFOs that occurred superimposed on the slow wave were transformed to phasors and adjusted by the phase of maximum firing in nearby units (i.e., maximum UP). We tested whether, in the seizure onset zone (SOZ), HFOs and associated action potentials (AP) occur more often at the UP-DOWN transition. We also examined ripple temporal correlations using cross correlograms.

Results

At the group level in the SOZ, HFO and HFO-associated AP probability was highest during the UP-DOWN transition of slow wave excitability (p<<0.001). In the non-SOZ, HFO and HFO-associated AP was highest during the DOWN-UP transition (p<<0.001). At the unit level in the SOZ, 15.6% and 20% of units exhibited more robust firing during ripples (Cohen's d=0.11-0.83) and fast ripples (d=0.36-0.90) at the UP-DOWN transition (p<0.05 f.d.r corrected), respectively. By comparison, also in the SOZ, 6.6% (d=0.14-0.30) and 8.5% (d=0.33-0.41) of units had significantly less firing during ripples and fast ripples at the UP-DOWN transition, respectively. Additional data shows ripple temporal correlations, involving global slow waves, between the hippocampus, entorhinal cortex, and parahippocampal gyrus were reduced by ~50-80% in the SOZ compared to the non-SOZ (N=3).

Significance

The UP-DOWN transition of slow wave excitability facilitates the activation of pathological neurons to generate pHFOs. The pathological neurons and pHFOs disrupt ripple temporal correlations across brain regions that transfer information and may be important in memory consolidation.

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.

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.

Integration of multiscale entropy and BASED scale of electroencephalography after adrenocorticotropic hormone therapy predict relapse of infantile spasms

BACKGROUND

Even though adrenocorticotropic hormone (ACTH) demonstrated powerful efficacy in the initially successful treatment of infantile spasms (IS), nearly half of patients have experienced a relapse. We sought to investigate whether features of electroencephalogram (EEG) predict relapse in those IS patients without structural brain abnormalities.

METHODS

We retrospectively reviewed data from children with IS who achieved initial response after ACTH treatment, along with EEG recorded within the last two days of treatment. The recurrence of epileptic spasms following treatment was tracked for 12 months. Subjects were categorized as either non-relapse or relapse groups. General clinical and EEG recordings were collected, burden of amplitudes and epileptiform discharges (BASED) score and multiscale entropy (MSE) were carefully explored for cross-group comparisons.

RESULTS

Forty-one patients were enrolled in the study, of which 26 (63.4%) experienced a relapse. The BASED score was significantly higher in the relapse group. MSE in the non-relapse group was significantly lower than the relapse group in the γ band but higher in the lower frequency range (δ, θ, α). Sensitivity and specificity were 85.71% and 92.31%, respectively, when combining MSE in the δ/γ frequency of the occipital region, plus BASED score were used to distinguish relapse from non-relapse groups.

CONCLUSIONS

BASED score and MSE of EEG after ACTH treatment could be used to predict relapse for IS patients without brain structural abnormalities. Patients with BASED score ≥ 3, MSE increased in higher frequency, and decreased in lower frequency had a high risk of relapse.

SleepSEEG: Automatic sleep scoring using intracranial EEG recordings only

OBJECTIVE

To perform automatic sleep scoring based only on intracranial EEG, without the need for scalp electroencephalography (EEG), electrooculography (EOG) and electromyography (EMG), in order to study sleep, epilepsy, and their interaction.

APPROACH

Data from 33 adult patients was used for development and training of the automatic scoring algorithm using both oscillatory and non-oscillatory spectral features. The first step consisted in unsupervised clustering of channels based on feature variability. For each cluster the classification was done in two steps, a multiclass tree followed by binary classification trees to distinguish the more challenging stage N1. The test data consisted in 11 patients, in whom the classification was done independently for each channel and then combined to get a single stage per epoch.

MAIN RESULTS

An overall agreement of 78% was observed in the test set between the sleep scoring of the algorithm and two human experts scoring based on scalp EEG, EOG and EMG. Balanced sensitivity and specificity were obtained for the different sleep stages. The performance was excellent for stages W, N2, and N3, and good for stage R, but with high variability across patients. The performance for the challenging stage N1 was poor, but at a similar level as for published algorithms based on scalp EEG. High confidence epochs in different stages (other than N1) can be identified with median per patient specificity >80%.

SIGNIFICANCE

The automatic algorithm can perform sleep scoring of long term recordings of patients with intracranial electrodes undergoing presurgical evaluation in the absence of scalp EEG, EOG and EMG, which are normally required to define sleep stages but are difficult to use in the context of intracerebral studies. It also constitutes a valuable tool to generate hypotheses regarding local aspects of sleep, and will be significant for sleep evaluation in clinical epileptology and neuroscience research.

Graph theoretical measures of fast ripples support the epileptic network hypothesis

The epileptic network hypothesis and epileptogenic zone hypothesis are two theories of ictogenesis. The network hypothesis posits that coordinated activity among interconnected nodes produces seizures. The epileptogenic zone hypothesis posits that distinct regions are necessary and sufficient for seizure generation. High-frequency oscillations, and particularly fast ripples, are thought to be biomarkers of the epileptogenic zone. We sought to test these theories by comparing high-frequency oscillation rates and networks in surgical responders and non-responders, with no appreciable change in seizure frequency or severity, within a retrospective cohort of 48 patients implanted with stereo-EEG electrodes. We recorded inter-ictal activity during non-rapid eye movement sleep and semi-automatically detected and quantified high-frequency oscillations. Each electrode contact was localized in normalized coordinates. We found that the accuracy of seizure onset zone electrode contact classification using high-frequency oscillation rates was not significantly different in surgical responders and non-responders, suggesting that in non-responders the epileptogenic zone partially encompassed the seizure onset zone(s) (P > 0.05). We also found that in the responders, fast ripple on oscillations exhibited a higher spectral content in the seizure onset zone compared with the non-seizure onset zone (P < 1 × 10⁻⁵). By contrast, in the non-responders, fast ripple had a lower spectral content in the seizure onset zone (P < 1 × 10⁻⁵). We constructed two different networks of fast ripple with a spectral content >350 Hz. The first was a rate–distance network that multiplied the Euclidian distance between fast ripple-generating contacts by the average rate of fast ripple in the two contacts. The radius of the rate–distance network, which excluded seizure onset zone nodes, discriminated non-responders, including patients not offered resection or responsive neurostimulation due to diffuse multifocal onsets, with an accuracy of 0.77 [95% confidence interval (CI) 0.56–0.98]. The second fast ripple network was constructed using the mutual information between the timing of the events to measure functional connectivity. For most non-responders, this network had a longer characteristic path length, lower mean local efficiency in the non-seizure onset zone, and a higher nodal strength among non-seizure onset zone nodes relative to seizure onset zone nodes. The graphical theoretical measures from the rate–distance and mutual information networks of 22 non- responsive neurostimulation treated patients was used to train a support vector machine, which when tested on 13 distinct patients classified non-responders with an accuracy of 0.92 (95% CI 0.75–1). These results indicate patients who do not respond to surgery or those not selected for resection or responsive neurostimulation can be explained by the epileptic network hypothesis that is a decentralized network consisting of widely distributed, hyperexcitable fast ripple-generating nodes.

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.

Drug-resistant epilepsy and the hypothesis of intrinsic severity: What about the high-frequency oscillations?

Drug‐resistant epilepsy (DRE) affects approximately one‐third of the patients with epilepsy. Based on experimental findings from animal models and brain tissue from patients with DRE, different hypotheses have been proposed to explain the cause(s) of drug resistance. One is the intrinsic severity hypothesis that posits that drug resistance is an inherent property of epilepsy related to disease severity. Seizure frequency is one measure of epilepsy severity, but frequency alone is an incomplete measure of severity and does not fully explain basic research and clinical studies on drug resistance; thus, other measures of epilepsy severity are needed. One such measure could be pathological high‐frequency oscillations (HFOs), which are believed to reflect the neuronal disturbances responsible for the development of epilepsy and the generation of spontaneous seizures. In this manuscript, we will briefly review the intrinsic severity hypothesis, describe basic and clinical research on HFOs in the epileptic brain, and based on this evidence discuss whether HFOs could be a clinical measure of epilepsy severity. Understanding the mechanisms of DRE is critical for producing breakthroughs in the development and testing of novel strategies for treatment.

Automated preprocessing and phase-amplitude coupling analysis of scalp EEG discriminates infantile spasms from controls during wakefulness

OBJECTIVE

Delta-gamma phase-amplitude coupling in EEG is useful for localizing epileptic sources and to evaluate severity in children with infantile spasms. We (1) develop an automated EEG preprocessing pipeline to clean data using artifact subspace reconstruction (ASR) and independent component (IC) analysis (ICA) and (2) evaluate delta-gamma modulation index (MI) as a method to distinguish children with epileptic spasms (cases) from normal controls during sleep and awake.

METHODS

Using 400 scalp EEG datasets (200 sleep, 200 awake) from 100 subjects, we calculated MI after applying high-pass and line-noise filters (Clean 0), and after ASR followed by either conservative (Clean 1) or stringent (Clean 2) artifactual IC rejection. Classification of cases and controls using MI was evaluated with Receiver Operating Characteristics (ROC) to obtain area under curve (AUC).

RESULTS

The artifact rejection algorithm reduced raw signal variance by 29-45% and 38-60% for Clean 1 and Clean 2, respectively. MI derived from sleep data, with or without preprocessing, robustly classified the groups (all AUC > 0.98). In contrast, group classification using MI derived from awake data was successful only after Clean 2 (AUC = 0.85).

CONCLUSIONS

We have developed an automated EEG preprocessing pipeline to perform artifact rejection and quantify delta-gamma modulation index.

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

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

Accuracy of high-frequency oscillations recorded intraoperatively for classification of epileptogenic regions

To see whether acute intraoperative recordings using stereo EEG (SEEG) electrodes can replace prolonged interictal intracranial EEG (iEEG) recording, making the process more efficient and safer, 10 min of iEEG were recorded following electrode implantation in 16 anesthetized patients, and 1-2 days later during non-rapid eye movement (REM) sleep. Ripples on oscillations (RonO, 80-250 Hz), ripples on spikes (RonS), sharp-spikes, fast RonO (fRonO, 250-600 Hz), and fast RonS (fRonS) were semi-automatically detected. HFO power and frequency were compared between the conditions using a generalized linear mixed-effects model. HFO rates were compared using a two-way repeated measures ANOVA with anesthesia type and SOZ as factors. A receiver-operating characteristic (ROC) curve analysis quantified seizure onset zone (SOZ) classification accuracy, and the scalar product was used to assess spatial reliability. Resection of contacts with the highest rate of events was compared with outcome. During sleep, all HFOs, except fRonO, were larger in amplitude compared to intraoperatively (p < 0.01). HFO frequency was also affected (p < 0.01). Anesthesia selection affected HFO and sharp-spike rates. In both conditions combined, sharp-spikes and all HFO subtypes were increased in the SOZ (p < 0.01). However, the increases were larger during the sleep recordings (p < 0.05). The area under the ROC curves for SOZ classification were significantly smaller for intraoperative sharp-spikes, fRonO, and fRonS rates (p < 0.05). HFOs and spikes were only significantly spatially reliable for a subset of the patients (p < 0.05). A failure to resect fRonO areas in the sleep recordings trended the most sensitive and accurate for predicting failure. In summary, HFO morphology is altered by anesthesia. Intraoperative SEEG recordings exhibit increased rates of HFOs in the SOZ, but their spatial distribution can differ from sleep recordings. Recording these biomarkers during non-REM sleep offers a more accurate delineation of the SOZ and possibly the epileptogenic zone.

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

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

Sleep and Epilepsy, Clinical Spectrum and Updated Review

Electroencephalogram (EEG) recording is essential in the evaluation of complex movement and behaviors during sleep, but in particular for differentiating epileptic versus nonepileptic events. In general, epileptiform discharges occur with greater density in the first few nonerapid eye movement cycles, and approximately 12% to 20% of seizures occur exclusively at night. This review examines the epilepsy types and syndromes whose presentation is strongly influenced by the sleep state, with an appraisal about the role that sleep plays in facilitating seizures, while deleaneatign EEG findings and clinical manifestation. The review will summarize the typical semiology of sleep-related hypermotor seizures and contrasted with those occurring during none/rapid eye movement parasomnias and sleep-related movement disorders.

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

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

Scalp EEG interictal high frequency oscillations as an objective biomarker of infantile spasms

OBJECTIVE

To investigate the diagnostic utility of high frequency oscillations (HFOs) via scalp electroencephalogram (EEG) in infantile spasms.

METHODS

We retrospectively analyzed interictal slow-wave sleep EEGs sampled at 2,000 Hz recorded from 30 consecutive patients who were suspected of having infantile spasms. We measured the rate of HFOs (80-500 Hz) and the strength of the cross-frequency coupling between HFOs and slow-wave activity (SWA) at 3-4 Hz and 0.5-1 Hz as quantified with modulation indices (MIs).

RESULTS

Twenty-three patients (77%) exhibited active spasms during the overnight EEG recording. Although the HFOs were detected in all children, increased HFO rate and MIs correlated with the presence of active spasms (p < 0.001 by HFO rate; p < 0.01 by MIs at 3-4 Hz; p = 0.02 by MIs at 0.5-1 Hz). The presence of active spasms was predicted by the logistic regression models incorporating HFO-related metrics (AUC: 0.80-0.98) better than that incorporating hypsarrhythmia (AUC: 0.61). The predictive performance of the best model remained favorable (87.5% accuracy) after a cross-validation procedure.

CONCLUSIONS

Increased rate of HFOs and coupling between HFOs and SWA are associated with active epileptic spasms.

SIGNIFICANCE

Scalp-recorded HFOs may serve as an objective EEG biomarker for active epileptic spasms.

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.

Phase-amplitude coupling of delta brush unveiling neuronal modulation development in the neonatal brain

INTRODUCTION

Delta brushes are an indicator of brain maturity on a neonatal EEG. We investigated phase-amplitude coupling (PAC) between slow delta waves and superimposed alpha-beta activity in delta brushes to elucidate the spatiotemporal developments of the delta brush with post-menstrual weeks (PMW).

METHODS

The subjects were 18 neurologically intact patients (seven girls). We analyzed EEG within 42 PMW. Patients were divided into four age groups as follows: PMW ≤30w; 31-34 w; 35-38 w; and 39-42 w. We selected up to three epochs of 2-minute EEG segments including delta brushes. We calculated the modulation index (MI), direct mean vector length (dMVL), and mean of phase angle of coupling by PAC between slow waves (0.5-1.5 Hz) and fast activities (8-25 Hz) in four regions (F: Fp1 and Fp2, C: C3 and C4, T: T3 and T4, O: O1 and O2).

RESULTS

We collected data from 18 patients and 31 epochs between 29 and 42 PMW, which comprised one, four, five, and eight patients, and two, seven, eight, and 14 epochs in the ≤30w, 31-34 w, 35-38 w, and 39-42 w groups, respectively. There were significant differences in the dMVL between the four regions in age groups ≤30w (P =  0.033) and 31-34w (0.017). Both MI and dMVL showed that delta brushes became higher in the occipital region from 32 to 36 PMW. The mean phase angle of coupling concentrated around either 0° or 180° for all age groups.

CONCLUSIONS

PAC analysis revealed the spatiotemporal relations of alpha-beta activities that are modulated by slow delta waves in neonatal delta brushes. The delta brushes appeared to be at a maximum around 32-36 PMW with the predominant occipital distribution. The PAC of the delta brush might represent the cortical neuronal fast activity that is modulated by slow delta waves of subcortical regions during a particular neonatal period.

Interictal scalp fast ripple occurrence and high frequency oscillation slow wave coupling in epileptic spasms

OBJECTIVE

Intracranial high frequency oscillation (HFO) occurrence rate (OR) and slow wave activity (SWA) coupling are potential markers of epileptogenicity in epileptic spasms (ES). Scalp ripple (R) detection and SWA coupling have been described in ES; however, the feasibility of scalp fast ripple (FR) detection and measurement of scalp FR coupling to SWA is not known. We evaluated interictal scalp R and FR OR and SWA coupling in pre-treatment EEG in children with short-term treatment-refractory ES compared to short-term treatment non-refractory ES.

METHODS

We retrospectively identified children with ES and identified HFOs using a semi-automated HFO detector on pre-treatment scalp EEG during sleep. We evaluated HFO OR and event-triggered modulation index (MI) to quantify R (100-250 Hz) and FR (250-600 Hz) coupling strength with different SWA passbands (0.5-1, 1-2, 2-3, 3-4, and 4-8 Hz). We used HFO phasor transform and circular statistics to evaluate phase coupling angle distributions.

RESULTS

We identified 15 children with ES with pre-treatment EEG recorded at 2000 Hz. Thirteen out of 15 patients had HFOs and were included for analysis. There were six treatment responders and seven nonresponders three months after treatment initiation. Responders and nonresponders were similar in age (6.1 vs 7.2 mo), ES diagnosis duration (0.7 vs 2.6 mo), and HFO OR (R: 1.07 vs 2.30/min, FR: 0.43 vs 1.96/min). No differences between responders and nonresponders were seen in HFO MI at different SWA. Coupling of R and FR to 2-3 Hz SWA demonstrated increased incidence rate ratio in nonresponders relative to responders at distinct phase coupling angle distributions.

CONCLUSIONS

This study demonstrates the feasibility of interictal scalp R and FR detection and quantification of scalp R and FR coupling to SWA in ES.

SIGNIFICANCE

HFO phase coupling with SWA may be useful as a marker of potential treatment refractoriness in patients with ES.

Ripples Have Distinct Spectral Properties and Phase-Amplitude Coupling With Slow Waves, but Indistinct Unit Firing, in Human Epileptogenic Hippocampus

Ripple oscillations (80-200 Hz) in the normal hippocampus are involved in memory consolidation during rest and sleep. In the epileptic brain, increased ripple and fast ripple (200-600 Hz) rates serve as a biomarker of epileptogenic brain. We report that both ripples and fast ripples exhibit a preferred phase angle of coupling with the trough-peak (or On-Off) state transition of the sleep slow wave in the hippocampal seizure onset zone (SOZ). Ripples on slow waves in the hippocampal SOZ also had a lower power, greater spectral frequency, and shorter duration than those in the non-SOZ. Slow waves in the mesial temporal lobe modulated the baseline firing rate of excitatory neurons, but did not significantly influence the increased firing rate associated with ripples. In summary, pathological ripples and fast ripples occur preferentially during the On-Off state transition of the slow wave in the epileptogenic hippocampus, and ripples do not require the increased recruitment of excitatory neurons.

Closed-Loop Acoustic Stimulation Enhances Sleep Oscillations But Not Memory Performance

Slow oscillations and spindle activity during non-rapid eye movement sleep have been implicated in memory consolidation. Closed-loop acoustic stimulation has previously been shown to enhance slow oscillations and spindle activity during sleep and improve verbal associative memory. We assessed the effect of closed-loop acoustic stimulation during a daytime nap on a virtual reality spatial navigation task in 12 healthy human subjects in a randomized within-subject crossover design. We show robust enhancement of slow oscillation and spindle activity during sleep. However, no effects on behavioral performance were observed when comparing real versus sham stimulation. To explore whether memory enhancement effects were task specific and dependent on nocturnal sleep, in a second experiment with 19 healthy subjects, we aimed to replicate a previous study that used closed-loop acoustic stimulation to enhance memory for word pairs. The methods used were as close as possible to those used in the original study, except that we used a double-blind protocol, in which both subject and experimenter were unaware of the test condition. Again, we successfully enhanced slow oscillation and spindle power, but again did not strengthen associative memory performance with stimulation. We conclude that enhancement of sleep oscillations may be insufficient to enhance memory performance in spatial navigation or verbal association tasks, and provide possible explanations for lack of behavioral replication.

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.

Nonictal EEG biomarkers for diagnosis and treatment

There are no reliable nonictal biomarkers for epilepsy, electroencephalography (EEG) or otherwise, but efforts to identify biomarkers that would predict the development of epilepsy after a potential epileptogenic insult, diagnose the existence of epilepsy, or assess the effects of antiseizure or antiepileptogenic interventions are relying heavily on electrophysiology. The most promising EEG biomarkers to date are pathologic high‐frequency oscillations (pHFOs), brief EEG events in the range of 100 to 600 Hz, which are believed to reflect summated action potentials from synchronously bursting neurons. Studies of patients with epilepsy, and experimental animal models, have been based primarily on direct brain recording, which makes pHFOs potentially useful for localizing the epileptogenic zone for surgical resection, but application for other diagnostic and therapeutic purposes is limited. Consequently, recent efforts have involved identification of HFOs recorded with scalp electrodes, and with magnetoencephalography, which may reflect the same pathophysiologic mechanisms as pHFOs recorded directly from the brain. The search is also on for other EEG changes that might serve as epilepsy biomarkers, and candidates include arcuate rhythms, which may reflect repetitive pHFOs, reduction in theta rhythm, which correlates with epileptogenesis in several rodent models of epilepsy, and shortened sleep spindles that correlate with ictogenesis.

Phase-amplitude coupling between interictal high-frequency activity and slow waves in epilepsy surgery

OBJECTIVE

We hypothesized that the modulation index (MI), a summary measure of the strength of phase-amplitude coupling between high-frequency activity (>150 Hz) and the phase of slow waves (3-4 Hz), would serve as a useful interictal biomarker for epilepsy presurgical evaluation.

METHODS

We investigated 123 patients who underwent focal cortical resection following extraoperative electrocorticography recording and had at least 1 year of postoperative follow-up. We examined whether consideration of MI would improve the prediction of postoperative seizure outcome. MI was measured at each intracranial electrode site during interictal slow-wave sleep. We compared the accuracy of prediction of patients achieving International League Against Epilepsy class 1 outcome between the full multivariate logistic regression model incorporating MI in addition to conventional clinical, seizure onset zone (SOZ), and neuroimaging variables, and the reduced logistic regression model incorporating all variables other than MI.

RESULTS

Ninety patients had class 1 outcome at the time of most recent follow-up (mean follow-up = 5.7 years). The full model had a noteworthy outcome predictive ability, as reflected by regression model fit R² of 0.409 and area under the curve (AUC) of receiver operating characteristic plot of 0.838. Incomplete resection of SOZ (P < 0.001), larger number of antiepileptic drugs at the time of surgery (P = 0.007), and larger MI in nonresected tissues relative to that in resected tissue (P = 0.020) were independently associated with a reduced probability of class 1 outcome. The reduced model had a lower predictive ability as reflected by R² of 0.266 and AUC of 0.767. Anatomical variability in MI existed among nonepileptic electrode sites, defined as those unaffected by magnetic resonance imaging lesion, SOZ, or interictal spike discharges. With MI adjusted for anatomical variability, the full model yielded the outcome predictive ability of R² of 0.422, AUC of 0.844, and sensitivity/specificity of 0.86/0.76.

SIGNIFICANCE

MI during interictal recording may provide useful information for the prediction of postoperative seizure outcome.

Phase-amplitude coupling and epileptogenesis in an animal model of mesial temporal lobe epilepsy

Polyrhythmic coupling of oscillatory components in electrophysiological signals results from the interactions between neuronal sub-populations within and between cell assemblies. Since the mechanisms underlying epileptic disorders should affect such interactions, abnormal level of cross-frequency coupling is expected to provide a signal marker of epileptogenesis. We measured phase-amplitude coupling (PAC), a form of cross-frequency coupling between neural oscillations, in a rodent model of mesial temporal lobe epilepsy. Sprague-Dawley rats (n = 4, 250-300 g) were injected with pilocarpine (380 mg/kg, i.p) to induce a status epilepticus (SE) that was stopped after 1 h with diazepam (5 mg/kg, s.c.) and ketamine (50 mg/kg, s.c.). Control animals (n = 6) did not receive any injection or treatment. Three days after SE, all animals were implanted with bipolar electrodes in the hippocampal CA3 subfield, entorhinal cortex, dentate gyrus and subiculum. Continuous video/EEG recordings were performed 24/7 at a sampling rate of 2 kHz, over 15 consecutive days. Pilocarpine-treated animals showed interictal spikes (5.25 (±2.5) per minute) and seizures (n = 32) that appeared 7 (±0.8) days after SE. We found that CA3 was the seizure onset zone in most epileptic animals, with stronger ongoing PAC coupling between seizures than in controls (Kruskal-Wallis test: chi2 (1,36) = 46.3, Bonferroni corrected, p < 0.001). Strong PAC in CA3 occurred between the phase of slow-wave oscillations (<1 Hz) and the amplitude of faster rhythms (50-180 Hz), with the strongest bouts of high-frequency activity occurring preferentially on the ascending phase of the slow wave. We also identified that cross-frequency coupling in CA3 (rho = 0.44, p < 0.001) and subiculum (rho = 0.41, p < 0.001) was positively correlated with the daily number of seizures. Overall, our study demonstrates that cross-frequency coupling may represent a signal marker in epilepsy and suggests that this methodology could be transferred to clinical scalp MEG and EEG recordings.

Comparing spiking and slow wave activity from invasive electroencephalography in patients with and without seizures

OBJECTIVES

To develop quantitative measures for estimating seizure probability, we examine intracranial EEG data from patient groups with three qualitative seizure probabilities: patients with drug resistant focal epilepsy (high), these patients during cortical stimulation (intermediate), and patients who have no history of seizures (low).

METHODS

Patients with focal epilepsy were implanted with subdural electrodes during presurgical evaluation. Patients without seizures were implanted during treatment with motor cortex stimulation for atypical facial pain.

RESULTS

The rate and amplitude of spikes correlate with qualitative seizure probability across patient groups and with proximity to the seizure onset zone in focal epilepsy patients. Spikes occur earlier during the negative oscillation of underlying slow activity (0.5-2 Hz) when seizure probability is increased. Similarly, coupling between slow and fast activity is increased.

CONCLUSIONS

There is likely a continuum of sharply contoured activity between non-epileptiform and epileptiform. Characteristics of spiking and how spikes relate to slow activity can be combined to predict seizure onset zones.

SIGNIFICANCE

Intracranial EEG data from patients without seizures represent a unique comparison group and highlight changes seen in spiking and slow wave activity with increased seizure probability. Slow wave activity and related physiology are an important potential biomarker for estimating seizure probability.

Strong coupling between slow oscillations and wide fast ripples in children with epileptic spasms: Investigation of modulation index and occurrence rate

OBJECTIVE

Epileptic spasms (ES) often become drug-resistant. To reveal the electrophysiological difference between children with ES (ES+) and without ES (ES-), we compared the occurrence rate (OR) of high-frequency oscillations (HFOs) and the modulation index (MI) of coupling between slow and fast oscillations. In ES+, we hypothesized that (1) pathological HFOs are more widely distributed and (2) slow oscillations show stronger coupling with pathological HFOs than in ES-.

METHODS

We retrospectively reviewed 24 children with drug-resistant multilobar onset epilepsy, who underwent intracranial video electroencephalography prior to multilobar resections. We measured the OR of HFOs and determined the electrodes with a high rate of HFOs by cluster analysis. We calculated MI, which reflects the degree of coupling between HFO (ripple/fast ripple [FR]) amplitude and 5 different frequency bands of delta and theta activities (0.5-1 Hz, 1-2 Hz, 2-3 Hz, 3-4 Hz, 4-8 Hz).

RESULTS

In ES+ (n = 10), the OR_(FRs) , the number of electrodes with high-rate FRs, and the MI_{(FRs & 3-4 Hz)} in all electrodes were significantly higher than in ES- (n = 14). In both the ES+ and ES- groups, MI_{(ripples/FRs & 3-4 Hz)} was the highest among the 5 frequency bands. Within the good seizure outcome group, the OR_(FRs) and the MI_{(FRs & 3-4 Hz)} in the resected area in ES+ were significantly higher than in ES- (OR_[FRs] , P = .04; MI_{[FRs & 3-4 Hz]} , P = .04).

SIGNIFICANCE

In ES+, the larger number of high-rate FR electrodes indicates more widespread epileptogenicity than in ES-. High values of OR_(FRs) and MI_{(FRs & 3-4 Hz)} in ES+ compared to ES- are a signature of the severity of epileptogenicity. We proved that ES+ children who achieved seizure freedom following multilobar resections exhibited strong coupling between slow oscillations and FRs.