iEEG Recordings Research Articles

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

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. 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. 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). 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.

Maintenance and transformation of representational formats during working memory prioritization

Visual working memory depends on both material-specific brain areas in the ventral visual stream (VVS) that support the maintenance of stimulus representations and on regions in the prefrontal cortex (PFC) that control these representations. How executive control prioritizes working memory contents and whether this affects their representational formats remains an open question, however. Here, we analyzed intracranial EEG (iEEG) recordings in epilepsy patients with electrodes in VVS and PFC who performed a multi-item working memory task involving a retro-cue. We employed Representational Similarity Analysis (RSA) with various Deep Neural Network (DNN) architectures to investigate the representational format of prioritized VWM content. While recurrent DNN representations matched PFC representations in the beta band (15–29 Hz) following the retro-cue, they corresponded to VVS representations in a lower frequency range (3–14 Hz) towards the end of the maintenance period. Our findings highlight the distinct coding schemes and representational formats of prioritized content in VVS and PFC.

Speech-induced suppression and vocal feedback sensitivity in human cortex

Across the animal kingdom, neural responses in the auditory cortex are suppressed during vocalization, and humans are no exception. A common hypothesis is that suppression increases sensitivity to auditory feedback, enabling the detection of vocalization errors. This hypothesis has been previously confirmed in non-human primates, however a direct link between auditory suppression and sensitivity in human speech monitoring remains elusive. To address this issue, we obtained intracranial electroencephalography (iEEG) recordings from 35 neurosurgical participants during speech production. We first characterized the detailed topography of auditory suppression, which varied across superior temporal gyrus (STG). Next, we performed a delayed auditory feedback (DAF) task to determine whether the suppressed sites were also sensitive to auditory feedback alterations. Indeed, overlapping sites showed enhanced responses to feedback, indicating sensitivity. Importantly, there was a strong correlation between the degree of auditory suppression and feedback sensitivity, suggesting suppression might be a key mechanism that underlies speech monitoring. Further, we found that when participants produced speech with simultaneous auditory feedback, posterior STG was selectively activated if participants were engaged in a DAF paradigm, suggesting that increased attentional load can modulate auditory feedback sensitivity.

Speech-induced suppression and vocal feedback sensitivity in human cortex.

Across the animal kingdom, neural responses in the auditory cortex are suppressed during vocalization, and humans are no exception. A common hypothesis is that suppression increases sensitivity to auditory feedback, enabling the detection of vocalization errors. This hypothesis has been previously confirmed in non-human primates, however a direct link between auditory suppression and sensitivity in human speech monitoring remains elusive. To address this issue, we obtained intracranial electroencephalography (iEEG) recordings from 35 neurosurgical participants during speech production. We first characterized the detailed topography of auditory suppression, which varied across superior temporal gyrus (STG). Next, we performed a delayed auditory feedback (DAF) task to determine whether the suppressed sites were also sensitive to auditory feedback alterations. Indeed, overlapping sites showed enhanced responses to feedback, indicating sensitivity. Importantly, there was a strong correlation between the degree of auditory suppression and feedback sensitivity, suggesting suppression might be a key mechanism that underlies speech monitoring. Further, we found that when participants produced speech with simultaneous auditory feedback, posterior STG was selectively activated if participants were engaged in a DAF paradigm, suggesting that increased attentional load can modulate auditory feedback sensitivity.

Intracranial neurophysiological correlates of rumination: A cross-sectional observational study

BackgroundRumination is a transdiagnostic psychological process that plays a prominent role in many common psychiatric disorders, albeit its neurophysiological basis remains elusive. Existing neuroimaging studies have highlighted the precuneus and hippocampus as two pivotal brain regions in rumination’s neural underpinnings. MethodsHere, we examined the intracranial electroencephalogram (iEEG) recordings from 21 patients with epilepsy during a naturalistic, continuous, active rumination state and measured the low frequency (1–8 Hz) and high gamma (70–150 Hz) band oscillation powers. ResultsWe observed enhanced low-frequency power in the precuneus and reduced high gamma power in the hippocampus during the rumination condition compared to the control condition. The hippocampal high gamma power reduction was associated with the self-reported reflection tendency. ConclusionsOur findings provided empirical evidence of the neurophysiological underpinnings of rumination and implicated a precuneus-hippocampus coupling across neural oscillation bands during an active rumination state.

Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study.

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

High-gamma and beta bursts in the Left Supramarginal Gyrus can accurately differentiate verbal memory states and performance.

The left supramarginal gyrus (LSMG) may mediate attention to memory, and gauge memory state and performance. We performed a secondary analysis of 142 verbal delayed free recall experiments, in patients with medically-refractory epilepsy with electrode contacts implanted in the LSMG. In 14 of 142 experiments (in 14 of 113 patients), the cross-validated convolutional neural networks (CNNs) that used 1-dimensional(1-D) pairs of convolved high-gamma and beta tensors, derived from the LSMG recordings, could label recalled words with an area under the receiver operating curve (AUROC) of greater than 60% [range: 60-90%]. These 14 patients were distinguished by: 1) higher amplitudes of high-gamma bursts; 2) distinct electrode placement within the LSMG; and 3) superior performance compared with a CNN that used a 1-D tensor of the broadband recordings in the LSMG. In a pilot study of 7 of these patients, we also cross-validated CNNs using paired 1-D convolved high-gamma and beta tensors, from the LSMG, to: a) distinguish word encoding epochs from free recall epochs [AUC 0.6-1]; and distinguish better performance from poor performance during delayed free recall [AUC 0.5-0.86]. These experiments show that bursts of high-gamma and beta generated in the LSMG are biomarkers of verbal memory state and performance.

Speech-induced suppression and vocal feedback sensitivity in human cortex.

Across the animal kingdom, neural responses in the auditory cortex are suppressed during vocalization, and humans are no exception. A common hypothesis is that suppression increases sensitivity to auditory feedback, enabling the detection of vocalization errors. This hypothesis has been previously confirmed in non-human primates, however a direct link between auditory suppression and sensitivity in human speech monitoring remains elusive. To address this issue, we obtained intracranial electroencephalography (iEEG) recordings from 35 neurosurgical participants during speech production. We first characterized the detailed topography of auditory suppression, which varied across superior temporal gyrus (STG). Next, we performed a delayed auditory feedback (DAF) task to determine whether the suppressed sites were also sensitive to auditory feedback alterations. Indeed, overlapping sites showed enhanced responses to feedback, indicating sensitivity. Importantly, there was a strong correlation between the degree of auditory suppression and feedback sensitivity, suggesting suppression might be a key mechanism that underlies speech monitoring. Further, we found that when participants produced speech with simultaneous auditory feedback, posterior STG was selectively activated if participants were engaged in a DAF paradigm, suggesting that increased attentional load can modulate auditory feedback sensitivity.

An electroencephalogram microdisplay to visualize neuronal activity on the brain surface.

Functional mapping during brain surgery is applied to define brain areas that control critical functions and cannot be removed. Currently, these procedures rely on verbal interactions between the neurosurgeon and electrophysiologist, which can be time-consuming. In addition, the electrode grids that are used to measure brain activity and to identify the boundaries of pathological versus functional brain regions have low resolution and limited conformity to the brain surface. Here, we present the development of an intracranial electroencephalogram (iEEG)-microdisplay that consists of freestanding arrays of 2048 GaN light-emitting diodes laminated on the back of micro-electrocorticography electrode grids. With a series of proof-of-concept experiments in rats and pigs, we demonstrate that these iEEG-microdisplays allowed us to perform real-time iEEG recordings and display cortical activities by spatially corresponding light patterns on the surface of the brain in the surgical field. Furthermore, iEEG-microdisplays allowed us to identify and display cortical landmarks and pathological activities from rat and pig models. Using a dual-color iEEG-microdisplay, we demonstrated coregistration of the functional cortical boundaries with one color and displayed the evolution of electrical potentials associated with epileptiform activity with another color. The iEEG-microdisplay holds promise to facilitate monitoring of pathological brain activity in clinical settings.

Comparative effectiveness of stereotactic, subdural, or hybrid intracranial EEG monitoring in epilepsy surgery.

Surgical intervention can be curative or palliative for drug-resistant focal epilepsy. However, if the seizure onset zone (SOZ) cannot be adequately localized via noninvasive tests, intracranial EEG (iEEG) recordings are often carried out to develop surgical plans in appropriate candidates. Stereotactic EEG (SEEG), subdural EEG (SDE), and SDE with depth electrodes (hybrid) are major tools used for investigation, but there is no class 1 or 2 evidence comparing the effectiveness of these modalities. The authors identified an institutional cohort of patients who underwent iEEG monitoring between 2001 and 2022. Demographic data, preoperative clinical features, iEEG intervention, and follow-up data were identified. Primary study endpoints included the following: 1) likelihood of SOZ localization; 2) likelihood of surgical treatment after iEEG; 3) seizure outcomes; and 4) complications. A total of 329 patients were identified (176 in the SEEG, 60 in the SDE, and 93 in the hybrid cohort) who were followed for a median of 5.4 (IQR 6.8) years. Baseline characteristics, including demographics, mean age at epilepsy diagnosis, mean age at iEEG investigation, number of preoperative antiseizure medications, and preoperative seizure frequency, were not statistically different across the 3 cohorts. Patients in the SEEG cohort were more likely to have their SOZ localized than were the patients in the SDE group (OR 2.3) and were less likely to undergo subsequent resection (OR 0.3) or to have complications (OR 0.4), although there was no statistical difference with respect to likelihood of undergoing any subsequent neurosurgical treatment, or with respect to favorable seizure outcomes. Patients in the hybrid cohort were more likely to have SOZ localized than were patients in the SDE group (OR 3.1), but were more likely to undergo resection (OR 4.9) or any neurosurgical treatment (OR 2.5) compared to patients in the SEEG group. Patients in the hybrid cohort had better seizure outcomes compared to the SDE (OR 2.3) but not to the SEEG group. Patients in the SEEG group were more likely to have their SOZ localized and patients in the SDE group were more likely to undergo resection, but they did not differ with respect to seizure outcomes.

Timing matters for accurate identification of the epileptogenic zone

ObjectiveInterictal biomarkers of the epileptogenic zone (EZ) and their use in machine learning models open promising avenues for improvement of epilepsy surgery evaluation. Currently, most studies restrict their analysis to short segments of intracranial EEG (iEEG). MethodsWe used 2381 hours of iEEG data from 25 patients to systematically select 5-minute segments across various interictal conditions. Then, we tested machine learning models for EZ localization using iEEG features calculated within these individual segments or across them and evaluated the performance by the area under the precision-recall curve (PRAUC). ResultsOn average, models achieved a score of 0.421 (the result of the chance classifier was 0.062). However, the PRAUC varied significantly across the segments (0.323–0.493). Overall, NREM sleep achieved the highest scores, with the best results of 0.493 in N2. When using data from all segments, the model performed significantly better than single segments, except NREM sleep segments. ConclusionsThe model based on a short segment of iEEG recording can achieve similar results as a model based on prolonged recordings. The analyzed segment should, however, be carefully and systematically selected, preferably from NREM sleep. SignificanceRandom selection of short iEEG segments may give rise to inaccurate localization of the EZ.

Benchmarking signal quality and spatiotemporal distribution of interictal spikes in prolonged human iEEG recordings using CorTec wireless brain interchange

Neuromodulation through implantable pulse generators (IPGs) represents an important treatment approach for neurological disorders. While the field has observed the success of state-of-the-art interventions, such as deep brain stimulation (DBS) or responsive neurostimulation (RNS), implantable systems face various technical challenges, including the restriction of recording from a limited number of brain sites, power management, and limited external access to the assessed neural data in a continuous fashion. To the best of our knowledge, for the first time in this study, we investigated the feasibility of recording human intracranial EEG (iEEG) using a benchtop version of the Brain Interchange (BIC) unit of CorTec, which is a portable, wireless, and externally powered implant with sensing and stimulation capabilities. We developed a MATLAB/SIMULINK-based rapid prototyping environment and a graphical user interface (GUI) to acquire and visualize the iEEG captured from all 32 channels of the BIC unit. We recorded prolonged iEEG (~ 24 h) from three human subjects with externalized depth leads using the BIC and commercially available clinical amplifiers simultaneously in the epilepsy monitoring unit (EMU). The iEEG signal quality of both streams was compared, and the results demonstrated a comparable power spectral density (PSD) in all the systems in the low-frequency band (< 80 Hz). However, notable differences were primarily observed above 100 Hz, where the clinical amplifiers were associated with lower noise floor (BIC-17 dB vs. clinical amplifiers < − 25 dB). We employed an established spike detector to assess and compare the spike rates in each iEEG stream. We observed over 90% conformity between the spikes rates and their spatial distribution captured with BIC and clinical systems. Additionally, we quantified the packet loss characteristic in the iEEG signal during the wireless data transfer and conducted a series of simulations to compare the performance of different interpolation methods for recovering the missing packets in signals at different frequency bands. We noted that simple linear interpolation has the potential to recover the signal and reduce the noise floor with modest packet loss levels reaching up to 10%. Overall, our results indicate that while tethered clinical amplifiers exhibited noticeably better noise floor above 80 Hz, epileptic spikes can still be detected successfully in the iEEG recorded with the externally powered wireless BIC unit opening the road for future closed-loop neuromodulation applications with continuous access to brain activity.

Epileptic tissue localization using graph-based networks in the high frequency oscillation range of intracranial electroencephalography

PurposeHigh frequency oscillations (HFOs) are an emerging biomarker of epilepsy. However, very few studies have investigated the functional connectivity of interictal iEEG signals in the frequency range of HFOs. Here, we study the corresponding functional networks using graph theory, and we assess their predictive value for automatic electrode classification in a cohort of 20 drug resistant patients. MethodsCoherence-based connectivity analysis was performed on the iEEG recordings, and six different local graph measures were computed in both sub-bands of the HFO frequency range (80–250 Hz and 250–500 Hz). Correlation analysis was implemented between the local graph measures and the ripple and fast ripple rates. Finally, the WEKA software was employed for training and testing different predictive models on the aforementioned local graph measures. ResultsThe ripple rate was significantly correlated with five out of six local graph measures in the functional network. For fast ripples, their rate was also significantly (but negatively) correlated with most of the local metrics. The results from WEKA showed that the Logistic Regression algorithm was able to classify highly HFO-contaminated electrodes with an accuracy of 82.5 % for ripples and 75.4 % for fast ripples. ConclusionFunctional connectivity networks in the HFO band could represent an alternative to the direct use of distinct HFO events, while also providing important insights about hub epileptic areas that can represent possible surgical targets. Automatic electrode classification through FC-based classifiers can help bypass the burden of manual HFO annotation, providing at the same time similar amount of information about the epileptic tissue.

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

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). 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. 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). 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.

Physiological and pathological neuronal connectivity in the living human brain based on intracranial EEG signals: the current state of research.

Over the past decades, studies of human brain networks have received growing attention as the assessment and modelling of connectivity in the brain is a topic of high impact with potential application in the understanding of human brain organization under both physiological as well as various pathological conditions. Under specific diagnostic settings, human neuronal signal can be obtained from intracranial EEG (iEEG) recording in epilepsy patients that allows gaining insight into the functional organisation of living human brain. There are two approaches to assess brain connectivity in the iEEG-based signal: evaluation of spontaneous neuronal oscillations during ongoing physiological and pathological brain activity, and analysis of the electrophysiological cortico-cortical neuronal responses, evoked by single pulse electrical stimulation (SPES). Both methods have their own advantages and limitations. The paper outlines available methodological approaches and provides an overview of current findings in studies of physiological and pathological human brain networks, based on intracranial EEG recordings.

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

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

Interictal magnetoencephalography abnormalities to guide intracranial electrode implantation and predict surgical outcome.

Intracranial EEG is the gold standard technique for epileptogenic zone localization but requires a preconceived hypothesis of the location of the epileptogenic tissue. This placement is guided by qualitative interpretations of seizure semiology, MRI, EEG and other imaging modalities, such as magnetoencephalography. Quantitative abnormality mapping using magnetoencephalography has recently been shown to have potential clinical value. We hypothesized that if quantifiable magnetoencephalography abnormalities were sampled by intracranial EEG, then patients' post-resection seizure outcome may be better. Thirty-two individuals with refractory neocortical epilepsy underwent magnetoencephalography and subsequent intracranial EEG recordings as part of presurgical evaluation. Eyes-closed resting-state interictal magnetoencephalography band power abnormality maps were derived from 70 healthy controls as a normative baseline. Magnetoencephalography abnormality maps were compared to intracranial EEG electrode implantation, with the spatial overlap of intracranial EEG electrode placement and cerebral magnetoencephalography abnormalities recorded. Finally, we assessed if the implantation of electrodes in abnormal tissue and subsequent resection of the strongest abnormalities determined by magnetoencephalography and intracranial EEG corresponded to surgical success. We used the area under the receiver operating characteristic curve as a measure of effect size. Intracranial electrodes were implanted in brain tissue with the most abnormal magnetoencephalography findings-in individuals that were seizure-free postoperatively (T = 3.9, P = 0.001) but not in those who did not become seizure-free. The overlap between magnetoencephalography abnormalities and electrode placement distinguished surgical outcome groups moderately well (area under the receiver operating characteristic curve = 0.68). In isolation, the resection of the strongest abnormalities as defined by magnetoencephalography and intracranial EEG separated surgical outcome groups well, area under the receiver operating characteristic curve = 0.71 and area under the receiver operating characteristic curve = 0.74, respectively. A model incorporating all three features separated surgical outcome groups best (area under the receiver operating characteristic curve = 0.80). Intracranial EEG is a key tool to delineate the epileptogenic zone and help render individuals seizure-free postoperatively. We showed that data-driven abnormality maps derived from resting-state magnetoencephalography recordings demonstrate clinical value and may help guide electrode placement in individuals with neocortical epilepsy. Additionally, our predictive model of postoperative seizure freedom, which leverages both magnetoencephalography and intracranial EEG recordings, could aid patient counselling of expected outcome.

Intracranial EEG Structure-Function Coupling and Seizure Outcomes After Epilepsy Surgery.

Surgery is an effective treatment for drug-resistant epilepsy, which modifies the brain's structure and networks to regulate seizure activity. Our objective was to examine the relationship between brain structure and function to determine the extent to which this relationship affects the success of the surgery in controlling seizures. We hypothesized that a stronger association between brain structure and function would lead to improved seizure control after surgery. We constructed functional and structural brain networks in patients with drug-resistant focal epilepsy by using presurgery functional data from intracranial EEG (iEEG) recordings, presurgery and postsurgery structural data from T1-weighted MRI, and presurgery diffusion-weighted MRI. We quantified the relationship (coupling) between structural and functional connectivity by using the Spearman rank correlation and analyzed this structure-function coupling at 2 spatial scales: (1) global iEEG network level and (2) individual iEEG electrode contacts using virtual surgeries. We retrospectively predicted postoperative seizure freedom by incorporating the structure-function connectivity coupling metrics and routine clinical variables into a cross-validated predictive model. We conducted a retrospective analysis on data from 39 patients who met our inclusion criteria. Brain areas implanted with iEEG electrodes had stronger structure-function coupling in seizure-free patients compared with those with seizure recurrence (p = 0.002, d = 0.76, area under the receiver operating characteristic curve [AUC] = 0.78 [95% CI 0.62-0.93]). Virtual surgeries on brain areas that resulted in stronger structure-function coupling of the remaining network were associated with seizure-free outcomes (p = 0.007, d = 0.96, AUC = 0.73 [95% CI 0.58-0.89]). The combination of global and local structure-function coupling measures accurately predicted seizure outcomes with a cross-validated AUC of 0.81 (95% CI 0.67-0.94). These measures were complementary to other clinical variables and, when included for prediction, resulted in a cross-validated AUC of 0.91 (95% CI 0.82-1.0), accuracy of 92%, sensitivity of 93%, and specificity of 91%. Our study showed that the strength of structure-function connectivity coupling may play a crucial role in determining the success of epilepsy surgery. By quantitatively incorporating structure-function coupling measures and standard-of-care clinical variables into presurgical evaluations, we may be able to better localize epileptogenic tissue and select patients for epilepsy surgery. This is a Class IV retrospective case series showing that structure-function mapping may help determine the outcome from surgical resection for treatment-resistant focal epilepsy.

Deep learning based source imaging provides strong sublobar localization of epileptogenic zone from MEG interictal spikes

Electromagnetic source imaging (ESI) offers unique capability of imaging brain dynamics for studying brain functions and aiding the clinical management of brain disorders. Challenges exist in ESI due to the ill-posedness of the inverse problem and thus the need of modeling the underlying brain dynamics for regularizations. Advances in generative models provide opportunities for more accurate and realistic source modeling that could offer an alternative approach to ESI for modeling the underlying brain dynamics beyond equivalent physical source models. However, it is not straightforward to explicitly formulate the knowledge arising from these generative models within the conventional ESI framework. Here we investigate a novel source imaging framework based on mesoscale neuronal modeling and deep learning (DL) that can learn the sensor-source mapping relationship directly from MEG data for ESI. Two DL-based ESI models were trained based on data generated by neural mass models and either generic or personalized head models. The robustness of the two DL models was evaluated by systematic computer simulations and clinical validation in a cohort of 29 drug-resistant focal epilepsy patients who underwent intracranial EEG (iEEG) evaluation or surgical resection. Results estimated from pre-operative MEG interictal spikes were quantified using the overlap with resection regions and the distance to the seizure-onset zone (SOZ) defined by iEEG recordings. The DL-based ESI provided robust results when no personalized head geometry is considered, reaching a spatial dispersion of 21.90 ± 19.03 mm, sublobar concordance of 83 %, and sublobar sensitivity and specificity of 66 and 97 % respectively. When using personalized head geometry derived from individual patients’ MRI in the training data, personalized DL-based ESI model can further improve the performance and reached a spatial dispersion of 8.19 ± 8.14 mm, sublobar concordance of 93 %, and sublobar sensitivity and specificity of 77 and 99 % respectively. When compared to the SOZ, the localization error of the personalized approach is 15.78 ± 5.54 mm, outperforming the conventional benchmarks. This work demonstrates that combining generative models and deep learning enables an accurate and robust imaging of epileptogenic zone from MEG recordings with strong sublobar precision, suggesting its added value to enhancing MEG source localization and imaging, and to epilepsy source localization and other clinical applications.

Unsupervised Multitaper Spectral Method for Identifying REM Sleep in Intracranial EEG Recordings Lacking EOG/EMG Data.

A large number of human intracranial EEG (iEEG) recordings have been collected for clinical purposes, in institutions all over the world, but the vast majority of these are unaccompanied by EOG and EMG recordings which are required to separate Wake episodes from REM sleep using accepted methods. In order to make full use of this extremely valuable data, an accurate method of classifying sleep from iEEG recordings alone is required. Existing methods of sleep scoring using only iEEG recordings accurately classify all stages of sleep, with the exception that wake (W) and rapid-eye movement (REM) sleep are not well distinguished. A novel multitaper (Wake vs. REM) alpha-rhythm classifier is developed by generalizing K-means clustering for use with multitaper spectral eigencoefficients. The performance of this unsupervised method is assessed on eight subjects exhibiting normal sleep architecture in a hold-out analysis and is compared against a classical power detector. The proposed multitaper classifier correctly identifies 36±6 min of REM in one night of recorded sleep, while incorrectly labeling less than 10% of all labeled 30 s epochs for all but one subject (human rater reliability is estimated to be near 80%), and outperforms the equivalent statistical-power classical test. Hold-out analysis indicates that when using one night's worth of data, an accurate generalization of the method on new data is likely. For the purpose of studying sleep, the introduced multitaper alpha-rhythm classifier further paves the way to making available a large quantity of otherwise unusable IEEG data.