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

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.

Normative atlases of high-frequency oscillation and spike rates under Sevoflurane anesthesia

OBJECTIVE

We analyzed the dose-dependent effects of Sevoflurane anesthesia on high-frequency oscillations (HFOs) and spike discharges at non-epileptic sites and evaluated their effectiveness in identifying the epileptogenic zone.

METHODS

We studied 21 children with drug-resistant focal epilepsy who achieved seizure control after focal resective surgery. Open-source detectors quantified HFO and spike rates during extraoperative and intraoperative intracranial EEG recordings performed before resection. We determined under which anesthetic conditions HFO and spike rates differentiated the seizure onset zone (SOZ) within the resected area from non-epileptic sites.

RESULTS

We analyzed 925 artifact-free electrodes, including 867 at non-epileptic sites and 58 at SOZ sites. Higher Sevoflurane doses significantly increased HFO and spike rates at non-epileptic sites, exhibiting spatial variability among different detectors. These biomarkers were elevated in the SOZ more than in non-epileptic sites under 2-4 vol% Sevoflurane anesthesia, with Cohen's d effect sizes above 3.0 and Mann-Whitney U-Test r effect sizes above 0.5.

CONCLUSIONS

We provided normative atlases of HFO and spike rates under different Sevoflurane anesthesia conditions. Sevoflurane elevates HFO and spike rates preferentially in the epileptogenic zone.

SIGNIFICANCE

Assessing the relative severity of biomarker levels across sites may be relevant for localizing the epileptogenic zone under Sevoflurane anesthesia.

Technical Feasibility of Delineating the Thalamic Gustatory Tract Using Tractography

BACKGROUND AND OBJECTIVES

Magnetic resonance-guided focused ultrasound (MRgFUS) has been increasingly performed in recent years as a minimally invasive treatment of essential tremor and tremor-dominant Parkinson disease. One of the side effects after treatment is dysgeusia. Some centers use tractography to facilitate the treatment planning. However, there have been no reports of identifying gustatory tracts so far. Our aim was to investigate the technical feasibility of isolating and visualizing the gustatory tracts, as well as to explore the relationship between the gustatory tract and the MRgFUS lesion using actual patient data.

METHODS

We used 20 randomly selected individuals from the Human Connectome Project database to perform tractography of the gustatory tracts. We defined region of interest as the dorsal region of the brainstem, Brodmann area 43 associated with taste perception, and a sphere with a 3-mm radius centered around the ventral intermediate nucleus in the anterior commissure-posterior commissure plane. We also examined the position of the gustatory tract in relation with other tracts, including the medial lemniscus, the pyramidal tract, and the dentatorubrothalamic tract. In addition, using the data of real patients with essential tremor, we investigated the distance between MRgFUS lesions and the gustatory tract and its association with the development of dysgeusia.

RESULTS

We delineated a mean of 15 streamlines of the gustatory tracts per subject in each hemisphere. There was no statistical difference in the localization of the gustatory tracts between the left and right cerebral hemispheres. The gustatory tract was located anteromedial to the medial lemniscus and posteromedial to the dentatorubrothalamic tract in the anterior commissure-posterior commissure plane. The distance from the MRgFUS lesion to the gustatory tract was significantly shorter in the case where dysgeusia occurred compared with nondysgeusia cases (P-value: .0068).

CONCLUSION

The thalamic gustatory tracts can be reliably visualized using tractography.

Dynamic functional connectivity in verbal cognitive control and word reading

Cognitive control processes enable the suppression of automatic behaviors and the initiation of appropriate responses. The Stroop color naming task serves as a benchmark paradigm for understanding the neurobiological model of verbal cognitive control. Previous research indicates a predominant engagement of the prefrontal and premotor cortex during the Stroop task compared to reading. We aim to further this understanding by creating a dynamic atlas of task-preferential modulations of functional connectivity through white matter. Patients undertook word-reading and Stroop tasks during intracranial EEG recording. We quantified task-related high-gamma amplitude modulations at 547 nonepileptic electrode sites, and a mixed model analysis identified regions and timeframes where these amplitudes differed between tasks. We then visualized white matter pathways with task-preferential functional connectivity enhancements at given moments. Word reading, compared to the Stroop task, exhibited enhanced functional connectivity in inter- and intra-hemispheric white matter pathways from the left occipital-temporal region 350-600 ms before response, including the posterior callosal fibers as well as the left vertical occipital, inferior longitudinal, inferior fronto-occipital, and arcuate fasciculi. The Stroop task showed enhanced functional connectivity in the pathways from the left middle-frontal pre-central gyri, involving the left frontal u-fibers and anterior callosal fibers. Automatic word reading largely utilizes the left occipital-temporal cortices and associated white matter tracts. Verbal cognitive control predominantly involves the left middle frontal and precentral gyri and its connected pathways. Our dynamic tractography atlases may serve as a novel resource providing insights into the unique neural dynamics and pathways of automatic reading and verbal cognitive control.

Connectivity of high-frequency bursts as SOZ localization biomarker

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

How to design optimal brain stimulation to modulate phase-amplitude coupling?

Objective.Phase-amplitude coupling (PAC), the coupling of the amplitude of a faster brain rhythm to the phase of a slower brain rhythm, plays a significant role in brain activity and has been implicated in various neurological disorders. For example, in Parkinson's disease, PAC between the beta (13-30 Hz) and gamma (30-100 Hz) rhythms in the motor cortex is exaggerated, while in Alzheimer's disease, PAC between the theta (4-8 Hz) and gamma rhythms is diminished. Modulating PAC (i.e. reducing or enhancing PAC) using brain stimulation could therefore open new therapeutic avenues. However, while it has been previously reported that phase-locked stimulation can increase PAC, it is unclear what the optimal stimulation strategy to modulate PAC might be. Here, we provide a theoretical framework to narrow down the experimental optimisation of stimulation aimed at modulating PAC, which would otherwise rely on trial and error.Approach.We make analytical predictions using a Stuart-Landau model, and confirm these predictions in a more realistic model of coupled neural populations.Main results.Our framework specifies the critical Fourier coefficients of the stimulation waveform which should be tuned to optimally modulate PAC. Depending on the characteristics of the amplitude response curve of the fast population, these components may include the slow frequency, the fast frequency, combinations of these, as well as their harmonics. We also show that the optimal balance of energy between these Fourier components depends on the relative strength of the endogenous slow and fast rhythms, and that the alignment of fast components with the fast rhythm should change throughout the slow cycle. Furthermore, we identify the conditions requiring to phase-lock stimulation to the fast and/or slow rhythms.Significance.Together, our theoretical framework lays the foundation for guiding the development of innovative and more effective brain stimulation aimed at modulating PAC for therapeutic benefit.

Cortical and white matter substrates supporting visuospatial working memory

OBJECTIVE

In tasks involving new visuospatial information, we rely on working memory, supported by a distributed brain network. We investigated the dynamic interplay between brain regions, including cortical and white matter structures, to understand how neural interactions change with different memory loads and trials, and their subsequent impact on working memory performance.

METHODS

Patients undertook a task of immediate spatial recall during intracranial EEG monitoring. We charted the dynamics of cortical high-gamma activity and associated functional connectivity modulations in white matter tracts.

RESULTS

Elevated memory loads were linked to enhanced functional connectivity via occipital longitudinal tracts, yet decreased through arcuate, uncinate, and superior-longitudinal fasciculi. As task familiarity grew, there was increased high-gamma activity in the posterior inferior-frontal gyrus (pIFG) and diminished functional connectivity across a network encompassing frontal, parietal, and temporal lobes. Early pIFG high-gamma activity was predictive of successful recall. Including this metric in a logistic regression model yielded an accuracy of 0.76.

CONCLUSIONS

Optimizing visuospatial working memory through practice is tied to early pIFG activation and decreased dependence on irrelevant neural pathways.

SIGNIFICANCE

This study expands our knowledge of human adaptation for visuospatial working memory, showing the spatiotemporal dynamics of cortical network modulations through white matter tracts.

PyHFO: lightweight deep learning-powered end-to-end high-frequency oscillations analysis application

Objective. This study aims to develop and validate an end-to-end software platform, PyHFO, that streamlines the application of deep learning (DL) methodologies in detecting neurophysiological biomarkers for epileptogenic zones from EEG recordings.Approach. We introduced PyHFO, which enables time-efficient high-frequency oscillation (HFO) detection algorithms like short-term energy and Montreal Neurological Institute and Hospital detectors. It incorporates DL models for artifact and HFO with spike classification, designed to operate efficiently on standard computer hardware.Main results. The validation of PyHFO was conducted on three separate datasets: the first comprised solely of grid/strip electrodes, the second a combination of grid/strip and depth electrodes, and the third derived from rodent studies, which sampled the neocortex and hippocampus using depth electrodes. PyHFO demonstrated an ability to handle datasets efficiently, with optimization techniques enabling it to achieve speeds up to 50 times faster than traditional HFO detection applications. Users have the flexibility to employ our pre-trained DL model or use their EEG data for custom model training.Significance. PyHFO successfully bridges the computational challenge faced in applying DL techniques to EEG data analysis in epilepsy studies, presenting a feasible solution for both clinical and research settings. By offering a user-friendly and computationally efficient platform, PyHFO paves the way for broader adoption of advanced EEG data analysis tools in clinical practice and fosters potential for large-scale research collaborations.

Modulation index predicts the effect of ethosuximide on developmental and epileptic encephalopathy with spike-and-wave activation in sleep

PURPOSE

In developmental and epileptic encephalopathy with spike-and-wave activation in sleep (DEE-SWAS), the thalamocortical network is suggested to play an important role in the pathophysiology of the progression from focal epilepsy to DEE-SWAS. Ethosuximide (ESM) exerts effects by blocking T-type calcium channels in thalamic neurons. With the thalamocortical network in mind, we studied the prediction of ESM effectiveness in DEE-SWAS treatment using phase-amplitude coupling (PAC) analysis.

METHODS

We retrospectively enrolled children with DEE-SWAS who had an electroencephalogram (EEG) recorded between January 2009 and September 2022 and were prescribed ESM at Okayama University Hospital. Only patients whose EEG showed continuous spike-and-wave during sleep were included. We extracted 5-min non-rapid eye movement sleep stage N2 segments from EEG recorded before starting ESM. We calculated the modulation index (MI) as the measure of PAC in pair combination comprising one of two fast oscillation types (gamma, 40-80 Hz; ripples, 80-150 Hz) and one of five slow-wave bands (delta, 0.5-1, 1-2, 2-3, and 3-4 Hz; theta, 4-8 Hz), and compared it between ESM responders and non-responders.

RESULTS

We identified 20 children with a diagnosis of DEE-SWAS who took ESM. Fifteen were ESM responders. Regarding gamma oscillations, significant differences were seen only in MI with 0.5-1 Hz slow waves in the frontal pole and occipital regions. Regarding ripples, ESM responders had significantly higher MI in coupling with all slow waves in the frontal pole region, 0.5-1, 3-4, and 4-8 Hz slow waves in the frontal region, 3-4 Hz slow waves in the parietal region, 0.5-1, 2-3, 3-4, and 4-8 Hz slow waves in the occipital region, and 3-4 Hz slow waves in the anterior-temporal region.

SIGNIFICANCE

High MI in a wider area of the brain may represent the epileptic network mediated by the thalamus in DEE-SWAS and may be a predictor of ESM effectiveness.

Technological Vanguard: the outstanding performance of the LTY-CNN model for the early prediction of epileptic seizures

BACKGROUND

Epilepsy is a common neurological disorder that affects approximately 60 million people worldwide. Characterized by unpredictable neural electrical activity abnormalities, it results in seizures with varying intensity levels. Electroencephalography (EEG), as a crucial technology for monitoring and predicting epileptic seizures, plays an essential role in improving the quality of life for people with epilepsy.

METHOD

This study introduces an innovative deep learning model, a lightweight triscale yielding convolutional neural network" (LTY-CNN), that is specifically designed for EEG signal analysis. The model integrates a parallel convolutional structure with a multihead attention mechanism to capture complex EEG signal features across multiple scales and enhance the efficiency achieved when processing time series data. The lightweight design of the LTY-CNN enables it to maintain high performance in environments with limited computational resources while preserving the interpretability and maintainability of the model.

RESULTS

In tests conducted on the SWEC-ETHZ and CHB-MIT datasets, the LTY-CNN demonstrated outstanding performance. On the SWEC-ETHZ dataset, the LTY-CNN achieved an accuracy of 99.9%, an area under the receiver operating characteristic curve (AUROC) of 0.99, a sensitivity of 99.9%, and a specificity of 98.8%. Furthermore, on the CHB-MIT dataset, it recorded an accuracy of 99%, an AUROC of 0.932, a sensitivity of 99.1%, and a specificity of 93.2%. These results signify the remarkable ability of the LTY-CNN to distinguish between epileptic seizures and nonseizure events. Compared to other existing epilepsy detection classifiers, the LTY-CNN attained higher accuracy and sensitivity.

CONCLUSION

The high accuracy and sensitivity of the LTY-CNN model demonstrate its significant potential for epilepsy management, particularly in terms of predicting and mitigating epileptic seizures. Its value in personalized treatments and widespread clinical applications reflects the broad prospects of deep learning in the health care sector. This also highlights the crucial role of technological innovation in enhancing the quality of life experienced by patients.