Intracranial EEG potentials estimated from MEG sources: A new approach to correlate MEG and iEEG data in epilepsy

Validating MEG estimated resting-state connectome with intracranial EEG

Magnetoencephalography (MEG) is widely used for studying resting-state brain connectivity. However, MEG source imaging is ill posed and has limited spatial resolution. This introduces source-leakage issues, making it challenging to interpret MEG-derived connectivity in resting states. To address this, we validated MEG-derived connectivity from 45 healthy participants using a normative intracranial EEG (iEEG) atlas. The MEG inverse problem was solved using the wavelet-maximum entropy on the mean method. We computed four connectivity metrics: amplitude envelope correlation (AEC), orthogonalized AEC (OAEC), phase locking value (PLV), and weighted-phase lag index (wPLI). We compared spatial correlation between MEG and iEEG connectomes across standard canonical frequency bands. We found moderate spatial correlations between MEG and iEEG connectomes for AEC and PLV. However, when considering metrics that correct/remove zero-lag connectivity (OAEC/wPLI), the spatial correlation between MEG and iEEG connectomes decreased. MEG exhibited higher zero-lag connectivity compared with iEEG. The correlations between MEG and iEEG connectomes suggest that relevant connectivity patterns can be recovered from MEG. However, since these correlations are moderate/low, MEG connectivity results should be interpreted with caution. Metrics that correct for zero-lag connectivity show decreased correlations, highlighting a trade-off; while MEG may capture more connectivity due to source-leakage, removing zero-lag connectivity can eliminate true connections.

M/EEG source localization for both subcortical and cortical sources using a convolutional neural network with a realistic head conductivity model

While electroencephalography (EEG) and magnetoencephalography (MEG) are well-established noninvasive methods in neuroscience and clinical medicine, they suffer from low spatial resolution. Electrophysiological source imaging (ESI) addresses this by noninvasively exploring the neuronal origins of M/EEG signals. Although subcortical structures are crucial to many brain functions and neuronal diseases, accurately localizing subcortical sources of M/EEG remains particularly challenging, and the feasibility is still a subject of debate. Traditional ESIs, which depend on explicitly defined regularization priors, have struggled to set optimal priors and accurately localize brain sources. To overcome this, we introduced a data-driven, deep learning-based ESI approach without the need for these priors. We proposed a four-layered convolutional neural network (4LCNN) designed to locate both subcortical and cortical sources underlying M/EEG signals. We also employed a sophisticated realistic head conductivity model using the state-of-the-art segmentation method of ten different head tissues from individual MRI data to generate realistic training data. This is the first attempt at deep learning-based ESI targeting subcortical regions. Our method showed excellent accuracy in source localization, particularly in subcortical areas compared to other methods. This was validated through M/EEG simulations, evoked responses, and invasive recordings. The potential for accurate source localization of the 4LCNNs demonstrated in this study suggests future contributions to various research endeavors such as the clinical diagnosis, understanding of the pathophysiology of various neuronal diseases, and basic brain functions.

EEG/MEG source imaging of deep brain activity within the maximum entropy on the mean framework: Simulations and validation in epilepsy

Electro/Magneto-EncephaloGraphy (EEG/MEG) source imaging (EMSI) of epileptic activity from deep generators is often challenging due to the higher sensitivity of EEG/MEG to superficial regions and to the spatial configuration of subcortical structures. We previously demonstrated the ability of the coherent Maximum Entropy on the Mean (cMEM) method to accurately localize the superficial cortical generators and their spatial extent. Here, we propose a depth-weighted adaptation of cMEM to localize deep generators more accurately. These methods were evaluated using realistic MEG/high-density EEG (HD-EEG) simulations of epileptic activity and actual MEG/HD-EEG recordings from patients with focal epilepsy. We incorporated depth-weighting within the MEM framework to compensate for its preference for superficial generators. We also included a mesh of both hippocampi, as an additional deep structure in the source model. We generated 5400 realistic simulations of interictal epileptic discharges for MEG and HD-EEG involving a wide range of spatial extents and signal-to-noise ratio (SNR) levels, before investigating EMSI on clinical HD-EEG in 16 patients and MEG in 14 patients. Clinical interictal epileptic discharges were marked by visual inspection. We applied three EMSI methods: cMEM, depth-weighted cMEM and depth-weighted minimum norm estimate (MNE). The ground truth was defined as the true simulated generator or as a drawn region based on clinical information available for patients. For deep sources, depth-weighted cMEM improved the localization when compared to cMEM and depth-weighted MNE, whereas depth-weighted cMEM did not deteriorate localization accuracy for superficial regions. For patients' data, we observed improvement in localization for deep sources, especially for the patients with mesial temporal epilepsy, for which cMEM failed to reconstruct the initial generator in the hippocampus. Depth weighting was more crucial for MEG (gradiometers) than for HD-EEG. Similar findings were found when considering depth weighting for the wavelet extension of MEM. In conclusion, depth-weighted cMEM improved the localization of deep sources without or with minimal deterioration of the localization of the superficial sources. This was demonstrated using extensive simulations with MEG and HD-EEG and clinical MEG and HD-EEG for epilepsy patients.

Spontaneous HFO Sequences Reveal Propagation Pathways for Precise Delineation of Epileptogenic Networks

Epilepsy, a neurological disorder affecting millions worldwide, poses great challenges in precisely delineating the epileptogenic zone - the brain region generating seizures - for effective treatment. High-frequency oscillations (HFOs) are emerging as promising biomarkers; however, the clinical utility is hindered by the difficulties in distinguishing pathological HFOs from non- epileptiform activities at single electrode and single patient resolution and understanding their dynamic role in epileptic networks. Here, we introduce an HFO-sequencing approach to analyze spontaneous HFOs traversing cortical regions in 40 drug-resistant epilepsy patients. This data- driven method automatically detected over 8.9 million HFOs, pinpointing pathological HFO- networks, and unveiled intricate millisecond-scale spatiotemporal dynamics, stability, and functional connectivity of HFOs in prolonged intracranial EEG recordings. These HFO sequences demonstrated a significant improvement in localization of epileptic tissue, with an 818.47% increase in concordance with seizure-onset zone (mean error: 2.92 mm), compared to conventional benchmarks. They also accurately predicted seizure outcomes for 90% AUC based on pre-surgical information using generalized linear models. Importantly, this mapping remained reliable even with short recordings (mean standard deviation: 3.23 mm for 30-minute segments). Furthermore, HFO sequences exhibited distinct yet highly repetitive spatiotemporal patterns, characterized by pronounced synchrony and predominant inward information flow from periphery towards areas involved in propagation, suggesting a crucial role for excitation-inhibition balance in HFO initiation and progression. Together, these findings shed light on the intricate organization of epileptic network and highlight the potential of HFO-sequencing as a translational tool for improved diagnosis, surgical targeting, and ultimately, better outcomes for vulnerable patients with drug-resistant epilepsy.

One Sentence Summary

Pathological fast brain oscillations travel like traffic along varied routes, outlining recurrently visited neural sites emerging as critical hotspots in epilepsy network.

The impact of ROI extraction method for MEG connectivity estimation: practical recommendations for the study of resting state data

Magnetoencephalography and electroencephalography (M/EEG) seed-based connectivity analysis requires the extraction of measures from regions of interest (ROI). M/EEG ROI-derived source activity can be treated in different ways. It is possible, for instance, to average each ROI's time series prior to calculating connectivity measures. Alternatively, one can compute connectivity maps for each element of the ROI prior to dimensionality reduction to obtain a single map. The impact of these different strategies on connectivity results is still unclear. Here, we address this question within a large MEG resting state cohort (N=113) and within simulated data. We consider 68 ROIs (Desikan-Kiliany atlas), two measures of connectivity (phase locking value-PLV, and its imaginary counterpart- ciPLV), and three frequency bands (theta 4-8 Hz, alpha 9-12 Hz, beta 15-30 Hz). We compare four extraction methods: (i) mean, or (ii) PCA of the activity within the seed or ROI before computing connectivity, map of the (iii) average, or (iv) maximum connectivity after computing connectivity for each element of the seed. Hierarchical clustering is then applied to compare connectivity outputs across multiple strategies, followed by direct contrasts across extraction methods. Finally, the results are validated by using a set of realistic simulations. We show that ROI-based connectivity maps vary remarkably across strategies in terms of connectivity magnitude and spatial distribution. Dimensionality reduction procedures conducted after computing connectivity are more similar to each-other, while PCA before approach is the most dissimilar to other approaches. Although differences across methods are consistent across frequency bands, they are influenced by the connectivity metric and ROI size. Greater differences were observed for ciPLV than PLV, and in larger ROIs. Realistic simulations confirmed that after aggregation procedures are generally more accurate but have lower specificity (higher rate of false positive connections). Though computationally demanding, after dimensionality reduction strategies should be preferred when higher sensitivity is desired. Given the remarkable differences across aggregation procedures, caution is warranted in comparing results across studies applying different methods.

Validating MEG source imaging of resting state oscillatory patterns with an intracranial EEG atlas

BACKGROUND

Magnetoencephalography (MEG) is a widely used non-invasive tool to estimate brain activity with high temporal resolution. However, due to the ill-posed nature of the MEG source imaging (MSI) problem, the ability of MSI to identify accurately underlying brain sources along the cortical surface is still uncertain and requires validation.

METHOD

We validated the ability of MSI to estimate the background resting state activity of 45 healthy participants by comparing it to the intracranial EEG (iEEG) atlas (https://mni-open-ieegatlas.

RESEARCH

mcgill.ca/). First, we applied wavelet-based Maximum Entropy on the Mean (wMEM) as an MSI technique. Next, we converted MEG source maps into intracranial space by applying a forward model to the MEG-reconstructed source maps, and estimated virtual iEEG (ViEEG) potentials on each iEEG channel location; we finally quantitatively compared those with actual iEEG signals from the atlas for 38 regions of interest in the canonical frequency bands.

RESULTS

The MEG spectra were more accurately estimated in the lateral regions compared to the medial regions. The regions with higher amplitude in the ViEEG than in the iEEG were more accurately recovered. In the deep regions, MEG-estimated amplitudes were largely underestimated and the spectra were poorly recovered. Overall, our wMEM results were similar to those obtained with minimum norm or beamformer source localization. Moreover, the MEG largely overestimated oscillatory peaks in the alpha band, especially in the anterior and deep regions. This is possibly due to higher phase synchronization of alpha oscillations over extended regions, exceeding the spatial sensitivity of iEEG but detected by MEG. Importantly, we found that MEG-estimated spectra were more comparable to spectra from the iEEG atlas after the aperiodic components were removed.

CONCLUSION

This study identifies brain regions and frequencies for which MEG source analysis is likely to be reliable, a promising step towards resolving the uncertainty in recovering intracerebral activity from non-invasive MEG studies.

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.

Comparison of beamformer and ICA for dynamic connectivity analysis: A simultaneous MEG-SEEG study

Magnetoencephalography (MEG) is a powerful tool for estimating brain connectivity with both good spatial and temporal resolution. It is particularly helpful in epilepsy to characterize non-invasively the epileptic networks. However, using MEG to map brain networks requires solving a difficult inverse problem that introduces uncertainty in the activity localization and connectivity measures. Our goal here was to compare independent component analysis (ICA) followed by dipole source localization and the linearly constrained minimum-variance beamformer (LCMV-BF) for characterizing regions with interictal epileptic activity and their dynamic connectivity. After a simulation study, we compared ICA and LCMV-BF results with intracerebral EEG (stereotaxic EEG, SEEG) recorded simultaneously in 8 epileptic patients, which provide a unique 'ground truth' to which non-invasive results can be confronted. We compared the signal time courses extracted applying ICA and LCMV-BF on MEG data to that of SEEG, both for the actual signals and the dynamic connectivity computed using cross-correlation (evolution of links in time). With our simulations, we illustrated the different effect of the temporal and spatial correlation among sources on the two methods. While ICA was more affected by the temporal correlation but robust against spatial configurations, LCMV-BF showed opposite behavior. Moreover, ICA seems more suited to retrieve the simulated networks. In case of real patient data, good MEG/SEEG correlation and good localization were obtained in 6 out of 8 patients. In 4 of them ICA had the best performance (higher correlation, lower localization distance). In terms of dynamic connectivity, the evolution in time of the cross-correlation links could be retrieved in 5 patients out of 6, however, with more variable results in terms of correlation and distance. In two patients LCMV-BF had better results than ICA. In one patient the two methods showed equally good outcomes, and in the remaining two patients ICA performed best. In conclusion, our results obtained by exploiting simultaneous MEG/SEEG recordings suggest that ICA and LCMV-BF have complementary qualities for retrieving the dynamics of interictal sources and their network interactions.

Data-driven approach for the delineation of the irritative zone in epilepsy in MEG

The reliable identification of the irritative zone (IZ) is a prerequisite for the correct clinical evaluation of medically refractory patients affected by epilepsy. Given the complexity of MEG data, visual analysis of epileptiform neurophysiological activity is highly time consuming and might leave clinically relevant information undetected. We recorded and analyzed the interictal activity from seven patients affected by epilepsy (Vectorview Neuromag), who successfully underwent epilepsy surgery (Engel > = II). We visually marked and localized characteristic epileptiform activity (VIS). We implemented a two-stage pipeline for the detection of interictal spikes and the delineation of the IZ. First, we detected candidate events from peaky ICA components, and then clustered events around spatio-temporal patterns identified by convolutional sparse coding. We used the average of clustered events to create IZ maps computed at the amplitude peak (PEAK), and at the 50% of the peak ascending slope (SLOPE). We validated our approach by computing the distance of the estimated IZ (VIS, SLOPE and PEAK) from the border of the surgically resected area (RA). We identified 25 spatiotemporal patterns mimicking the underlying interictal activity (3.6 clusters/patient). Each cluster was populated on average by 22.1 [15.0–31.0] spikes. The predicted IZ maps had an average distance from the resection margin of 8.4 ± 9.3 mm for visual analysis, 12.0 ± 16.5 mm for SLOPE and 22.7 ±. 16.4 mm for PEAK. The consideration of the source spread at the ascending slope provided an IZ closer to RA and resembled the analysis of an expert observer. We validated here the performance of a data-driven approach for the automated detection of interictal spikes and delineation of the IZ. This computational framework provides the basis for reproducible and bias-free analysis of MEG recordings in epilepsy.

Hierarchical Bayesian modeling of the relationship between task-related hemodynamic responses and cortical excitability

Investigating the relationship between task-related hemodynamic responses and cortical excitability is challenging because it requires simultaneous measurement of hemodynamic responses while applying noninvasive brain stimulation. Moreover, cortical excitability and task-related hemodynamic responses are both associated with inter-/intra-subject variability. To reliably assess such a relationship, we applied hierarchical Bayesian modeling. This study involved 16 healthy subjects who underwent simultaneous Paired Associative Stimulation (PAS10, PAS25, Sham) while monitoring brain activity using functional Near-Infrared Spectroscopy (fNIRS), targeting the primary motor cortex (M1). Cortical excitability was measured by Motor Evoked Potentials (MEPs), and the motor task-related hemodynamic responses were measured using fNIRS 3D reconstructions. We constructed three models to investigate: (1) PAS effects on the M1 excitability, (2) PAS effects on fNIRS hemodynamic responses to a finger tapping task, and (3) the correlation between PAS effects on M1 excitability and PAS effects on task-related hemodynamic responses. Significant increase in cortical excitability was found following PAS25, whereas a small reduction of the cortical excitability was shown after PAS10 and a subtle increase occurred after sham. Both HbO and HbR absolute amplitudes increased after PAS25 and decreased after PAS10. The probability of the positive correlation between modulation of cortical excitability and hemodynamic activity was 0.77 for HbO and 0.79 for HbR. We demonstrated that PAS stimulation modulates task-related cortical hemodynamic responses in addition to M1 excitability. Moreover, the positive correlation between PAS modulations of excitability and hemodynamics brought insight into understanding the fundamental properties of cortical function and cortical excitability.

The Impact of Magnetoencephalography-Directed Stereo-Electroencephalography Depth Electrode Implantation on Seizure Control Outcome in Children

Introduction The use of magnetoencephalography (MEG) in localizing epileptic foci and directing surgical treatment of medically refractory epilepsy is well established in clinical practice; however, it has not yet been incorporated into the routine planning of stereo-electroencephalography (EEG) (SEEG) depth electrode trajectories during invasive intracranial evaluation for epileptic foci localization. In this study, we assess the impact of MEG-directed SEEG on seizure outcomes in a pediatric cohort. Methods A retrospective analysis was performed on a single-institution cohort of pediatric patients with medically refractory epilepsy who underwent epilepsy surgery. The primary endpoint was the reduction in seizure burden as determined by dichotomized Engel scores (favorable outcome: Engel scores I and II; poor outcome: Engel scores III and IV). Results Thirty-seven patients met the inclusion criteria (24 males and 13 females). The median age at seizure onset was three years, the median age at surgery was 14.1 years, and the median follow-up length was 30.8 months. Concordance was noted in 7/10 (70%) patients who received MEG-directed SEEG. Good clinical outcomes were achieved in 70% of MEG-directed SEEG patients, compared to 59.4% in their counterparts; however, this difference was not statistically significant (p=0.72). We noted no statistically significant association between sex, disease laterality, or age at surgery and good clinical outcomes. Conclusions Patients who underwent MEG-directed SEEG had favorable clinical outcomes, which demonstrated the practicability of this technique for determining SEEG electrode placement. Although no significant difference in clinical outcomes was obtained between the two groups, this may have been due to low statistical power. Future prospective, multi-institutional investigations to assess the benefit of MEG-directed SEEG are warranted.

Virtual implantation using conventional scalp EEG delineates seizure onset and predicts surgical outcome in children with epilepsy

OBJECTIVE

Delineation of the seizure onset zone (SOZ) is required in children with drug resistant epilepsy (DRE) undergoing neurosurgery. Intracranial EEG (icEEG) serves as gold standard but has limitations. Here, we examine the utility of virtual implantation with electrical source imaging (ESI) on ictal scalp EEG for mapping the SOZ and predict surgical outcome.

METHODS

We retrospectively analyzed EEG data from 35 children with DRE who underwent surgery and dichotomized into seizure-free (SF) and non-seizure-free (NSF). We estimated virtual sensors (VSs) at brain locations that matched icEEG implantation and compared ictal patterns at VSs vs icEEG. We calculated the agreement between VSs SOZ and clinically defined SOZ and built receiver operating characteristic (ROC) curves to test whether it predicted outcome.

RESULTS

Twenty-one patients were SF after surgery. Moderate agreement between virtual and icEEG patterns was observed (kappa = 0.45, p < 0.001). Virtual SOZ agreement with clinically defined SOZ was higher in SF vs NSF patients (66.6% vs 41.6%, p = 0.01). Anatomical concordance of virtual SOZ with clinically defined SOZ predicted outcome (AUC = 0.73; 95% CI: 0.57-0.89; sensitivity = 66.7%; specificity = 78.6%; accuracy = 71.4%).

CONCLUSIONS

Virtual implantation on ictal scalp EEG can approximate the SOZ and predict outcome.

SIGNIFICANCE

SOZ mapping with VSs may contribute to tailoring icEEG implantation and predict outcome.

Clinical Yield of Electromagnetic Source Imaging and Hemodynamic Responses in Epilepsy: Validation With Intracerebral Data

BACKGROUND AND OBJECTIVES

Accurate delineation of the seizure-onset zone (SOZ) in focal drug-resistant epilepsy often requires stereo-EEG (SEEG) recordings. Our aims were to propose a truly objective and quantitative comparison between EEG/magnetoencephalography (MEG) source imaging (EMSI), EEG/fMRI responses for similar spikes with primary irritative zone (PIZ) and SOZ defined by SEEG and to evaluate the value of EMSI and EEG/fMRI to predict postsurgical outcome.

METHODS

We identified patients with drug-resistant epilepsy who underwent EEG/MEG, EEG/fMRI, and subsequent SEEG at the Epilepsy Service from the Montreal Neurological Institute and Hospital. We quantified multimodal concordance within the SEEG channel space as spatial overlap with PIZ/SOZ and distances to the spike-onset, spike maximum amplitude and seizure core intracerebral channels by applying a new methodology consisting of converting EMSI results into SEEG electrical potentials (EMSIe-SEEG) and projecting the most significant fMRI response on the SEEG channels (fMRIp-SEEG). Spatial overlaps with PIZ/SOZ (AUCPIZ, AUCSOZ) were assessed by using the area under the receiver operating characteristic curve (AUC). Here, AUC represents the probability that a randomly picked active contact exhibited higher amplitude when located inside the spatial reference than outside.

RESULTS

Seventeen patients were included. Mean spatial overlaps with the PIZ and SOZ were 0.71 and 0.65 for EMSIe-SEEG and 0.57 and 0.62 for fMRIp-SEEG. Good EMSIe-SEEG spatial overlap with the PIZ was associated with smaller distance from the maximum EMSIe-SEEG contact to the spike maximum amplitude channel (median distance 14 mm). Conversely, good fMRIp-SEEG spatial overlap with the SOZ was associated with smaller distances from the maximum fMRIp-SEEG contact to the spike-onset and seizure core channels (median distances 10 and 5 mm, respectively). Surgical outcomes were correctly predicted by EEG/MEG in 12 of 15 (80%) patients and EEG/fMRI in 6 of 11(54%) patients.

DISCUSSION

With the use of a unique quantitative approach estimating EMSI and fMRI results in the reference SEEG channel space, EEG/MEG and EEG/fMRI accurately localized the SOZ and the PIZ. Precisely, EEG/MEG more accurately localized the PIZ, whereas EEG/fMRI was more sensitive to the SOZ. Both neuroimaging techniques provide complementary localization that can help guide SEEG implantation and select good candidates for surgery.

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean

Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40[Formula: see text] along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.

Assessment of Effective Network Connectivity among MEG None Contaminated Epileptic Transitory Events

Characterizing epileptogenic zones EZ (sources responsible of excessive discharges) would assist a neurologist during epilepsy diagnosis. Locating efficiently these abnormal sources among magnetoencephalography (MEG) biomarker is obtained by several inverse problem techniques. These techniques present different assumptions and particular epileptic network connectivity. Here, we proposed to evaluate performances of distributed inverse problem in defining EZ. First, we applied an advanced technique based on Singular Value Decomposition (SVD) to recover only pure transitory activities (interictal epileptiform discharges). We evaluated our technique's robustness in separation between transitory and ripples versus frequency range, transitory shapes, and signal to noise ratio on simulated data (depicting both epileptic biomarkers and respecting time series and spectral properties of realistic data). We validated our technique on MEG signal using detector precision on 5 patients. Then, we applied four methods of inverse problem to define cortical areas and neural generators of excessive discharges. We computed network connectivity of each technique. Then, we confronted obtained noninvasive networks to intracerebral EEG transitory network connectivity using nodes in common, connection strength, distance metrics between concordant nodes of MEG and IEEG, and average propagation delay. Coherent Maximum Entropy on the Mean (cMEM) proved a high matching between MEG network connectivity and IEEG based on distance between active sources, followed by Exact low-resolution brain electromagnetic tomography (eLORETA), Dynamical Statistical Parametric Mapping (dSPM), and Minimum norm estimation (MNE). Clinical performance was interesting for entire methods providing in an average of 73.5% of active sources detected in depth and seen in MEG, and vice versa, about 77.15% of active sources were detected from MEG and seen in IEEG. Investigated problem techniques succeed at least in finding one part of seizure onset zone. dSPM and eLORETA depict the highest connection strength among all techniques. Propagation delay varies in this range [18, 25]ms, knowing that eLORETA ensures the lowest propagation delay (18 ms) and the closet one to IEEG propagation delay.

Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean

In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.

Presurgical Evaluation of Epilepsy Using Resting-State MEG Functional Connectivity

Epilepsy is caused by abnormal electrical discharges (clinically identified by electrophysiological recording) in a specific part of the brain [originating in only one part of the brain, namely, the epileptogenic zone (EZ)]. Epilepsy is now defined as an archetypical hyperexcited neural network disorder. It can be investigated through the network analysis of interictal discharges, ictal discharges, and resting-state functional connectivity. Currently, there is an increasing interest in embedding resting-state connectivity analysis into the preoperative evaluation of epilepsy. Among the various neuroimaging technologies employed to achieve brain functional networks, magnetoencephalography (MEG) with the excellent temporal resolution is an ideal tool for estimating the resting-state connectivity between brain regions, which can reveal network abnormalities in epilepsy. What value does MEG resting-state functional connectivity offer for epileptic presurgical evaluation? Regarding this topic, this paper introduced the origin of MEG and the workflow of constructing source-space functional connectivity based on MEG signals. Resting-state functional connectivity abnormalities correlate with epileptogenic networks, which are defined by the brain regions involved in the production and propagation of epileptic activities. This paper reviewed the evidence of altered epileptic connectivity based on low- or high-frequency oscillations (HFOs) and the evidence of the advantage of using simultaneous MEG and intracranial electroencephalography (iEEG) recordings. More importantly, this review highlighted that MEG-based resting-state functional connectivity has the potential to predict postsurgical outcomes. In conclusion, resting-state MEG functional connectivity has made a substantial progress toward serving as a candidate biomarker included in epileptic presurgical evaluations.

Dynamic analysis on simultaneous iEEG-MEG data via hidden Markov model

BACKGROUND

Intracranial electroencephalography (iEEG) recordings are used for clinical evaluation prior to surgical resection of the focus of epileptic seizures and also provide a window into normal brain function. A major difficulty with interpreting iEEG results at the group level is inconsistent placement of electrodes between subjects making it difficult to select contacts that correspond to the same functional areas. Recent work using time delay embedded hidden Markov model (HMM) applied to magnetoencephalography (MEG) resting data revealed a distinct set of brain states with each state engaging a specific set of cortical regions. Here we use a rare group dataset with simultaneously acquired resting iEEG and MEG to test whether there is correspondence between HMM states and iEEG power changes that would allow classifying iEEG contacts into functional clusters.

METHODS

Simultaneous MEG-iEEG recordings were performed at rest on 11 patients with epilepsy whose intracranial electrodes were implanted for pre-surgical evaluation. Pre-processed MEG sensor data was projected to source space. Time delay embedded HMM was then applied to MEG time series. At the same time, iEEG time series were analyzed with time-frequency decomposition to obtain spectral power changes with time. To relate MEG and iEEG results, correlations were computed between HMM probability time courses of state activation and iEEG power time course from the mid contact pair for each electrode in equally spaced frequency bins and presented as correlation spectra for the respective states and iEEG channels. Association of iEEG electrodes with HMM states based on significant correlations was compared to that based on the distance to peaks in subject-specific state topographies.

RESULTS

Five HMM states were inferred from MEG. Two of them corresponded to the left and the right temporal activations and had a spectral signature primarily in the theta/alpha frequency band. All the electrodes had significant correlations with at least one of the states (p < 0.05 uncorrected) and for 27/50 electrodes these survived within-subject FDR correction (q < 0.05). These correlations peaked in the theta/alpha band. There was a highly significant dependence between the association of states and electrodes based on functional correlations and that based on spatial proximity (p = 5.6e-6,χ2 test for independence). Despite the potentially atypical functional anatomy and physiological abnormalities related to epilepsy, HMM model estimated from the patient group was very similar to that estimated from healthy subjects.

CONCLUSION

Epilepsy does not preclude HMM analysis of interictal data. The resulting group functional states are highly similar to those reported for healthy controls. Power changes recorded with iEEG correlate with HMM state time courses in the alpha-theta band and the presence of this correlation can be related to the spatial location of electrode contacts close to the individual peaks of the corresponding state topographies. Thus, the hypothesized relation between iEEG contacts and HMM states exists and HMM could be further explored as a method for identifying comparable iEEG channels across subjects for the purposes of group analysis.

Noninvasive high-frequency oscillations riding spikes delineates epileptogenic sources

Millions of people affected by epilepsy may undergo surgical resection of the epileptic tissues to stop seizures if such epileptic foci can be accurately delineated. High-frequency oscillations (HFOs), existing in electroencephalography, are highly correlated with epileptic brain, which is promising for guiding successful neurosurgery. However, it is unclear whether and how pathological HFOs can be differentiated to localize the epileptogenic tissues given the presence of various nonepileptic high-frequency activities. Here, we show morphological and source imaging evidence that pathological HFOs can be identified by the concurrence of epileptiform spikes. We describe a framework to delineate the underlying epileptogenicity using this biomarker. Our work may offer translational tools to improve treatments by noninvasively demarking pathological activities and hence epileptic foci.

High-frequency oscillations (HFOs) are a promising biomarker for localizing epileptogenic brain and guiding successful neurosurgery. However, the utility and translation of noninvasive HFOs, although highly desirable, is impeded by the difficulty in differentiating pathological HFOs from nonepileptiform high-frequency activities and localizing the epileptic tissue using noninvasive scalp recordings, which are typically contaminated with high noise levels. Here, we show that the consistent concurrence of HFOs with epileptiform spikes (pHFOs) provides a tractable means to identify pathological HFOs automatically, and this in turn demarks an epileptiform spike subgroup with higher epileptic relevance than the other spikes in a cohort of 25 temporal epilepsy patients (including a total of 2,967 interictal spikes and 1,477 HFO events). We found significant morphological distinctions of HFOs and spikes in the presence/absence of this concurrent status. We also demonstrated that the proposed pHFO source imaging enhanced localization of epileptogenic tissue by 162% (∼5.36 mm) for concordance with surgical resection and by 186% (∼12.48 mm) with seizure-onset zone determined by invasive studies, compared to conventional spike imaging, and demonstrated superior congruence with the surgical outcomes. Strikingly, the performance of spike imaging was selectively boosted by the presence of spikes with pHFOs, especially in patients with multitype spikes. Our findings suggest that concurrent HFOs and spikes reciprocally discriminate pathological activities, providing a translational tool for noninvasive presurgical diagnosis and postsurgical evaluation in vulnerable patients.

Ant Colony System Optimization for Spatiotemporal Modelling of Combined EEG and MEG Data

Electroencephalography/Magnetoencephalography (EEG/MEG) source localization involves the estimation of neural activity inside the brain volume that underlies the EEG/MEG measures observed at the sensor array. In this paper, we consider a Bayesian finite spatial mixture model for source reconstruction and implement Ant Colony System (ACS) optimization coupled with Iterated Conditional Modes (ICM) for computing estimates of the neural source activity. Our approach is evaluated using simulation studies and a real data application in which we implement a nonparametric bootstrap for interval estimation. We demonstrate improved performance of the ACS-ICM algorithm as compared to existing methodology for the same spatiotemporal model.

Magnetoencephalography resting state connectivity patterns as indicatives of surgical outcome in epilepsy patients

OBJECTIVE

Focal epilepsy is a disorder affecting several brain networks; however, epilepsy surgery usually targets a restricted region, the so-called epileptic focus. There is a growing interest in embedding resting state (RS) connectivity analysis into pre-surgical workup.

APPROACH

In this retrospective study, we analyzed Magnetoencephalography (MEG) long-range RS functional connectivity patterns in patients with drug-resistant focal epilepsy. MEG recorded prior to surgery from seven seizure-free (Engel Ia) and five non seizure-free (Engel III or IV) patients were analyzed (minimum 2-years post-surgical follow-up). MEG segments without any detectable epileptic activity were source localized using wavelet-based Maximum Entropy on the Mean method. Amplitude envelope correlation in the theta (4-8 Hz), alpha (8-13 Hz), and beta (13-26 Hz) bands were used for assessing connectivity.

MAIN RESULTS

For seizure-free patients, we found an isolated epileptic network characterized by weaker connections between the brain region where interictal epileptic discharges (IED) are generated and the rest of the cortex, when compared to connectivity between the corresponding contralateral homologous region and the rest of the cortex. Contrarily, non seizure-free patients exhibited a widespread RS epileptic network characterized by stronger connectivity between the IED generator and the rest of the cortex, in comparison to the contralateral region and the cortex. Differences between the two seizure outcome groups concerned mainly distant long-range connections and were found in the alpha-band.

SIGNIFICANCE

Importantly, these connectivity patterns suggest specific mechanisms describing the underlying organization of the epileptic network and were detectable at the individual patient level, supporting the prospect use of MEG connectivity patterns in epilepsy to predict post-surgical seizure outcome.

Effects of Independent Component Analysis on Magnetoencephalography Source Localization in Pre-surgical Frontal Lobe Epilepsy Patients

Objective: Magnetoencephalography source imaging (MSI) of interictal epileptiform discharges (IED) is a useful presurgical tool in the evaluation of drug-resistant frontal lobe epilepsy (FLE) patients. Yet, failures in MSI can arise related to artifacts and to interference of background activity. Independent component analysis (ICA) is a popular denoising procedure but its clinical application remains challenging, as the selection of multiple independent components (IC) is controversial, operator dependent, and time consuming. We evaluated whether selecting only one IC of interest based on its similarity with the average IED field improves MSI in FLE.

Methods: MSI was performed with the equivalent current dipole (ECD) technique and two distributed magnetic source imaging (dMSI) approaches: minimum norm estimate (MNE) and coherent Maximum Entropy on the Mean (cMEM). MSI accuracy was evaluated under three conditions: (1) ICA of continuous data (Cont_ICA), (2) ICA at the time of IED (IED_ICA), and (3) without ICA (No_ICA). Localization performance was quantitatively measured as actual distance of the source maximum in relation to the focus (Dmin), and spatial dispersion (SD) for dMSI.

Results: After ICA, ECD Dmin did not change significantly (p > 0.200). For both dMSI techniques, ICA application worsened the source localization accuracy. We observed a worsening of both MNE Dmin (p < 0.05, consistently) and MNE SD (p < 0.001, consistently) for both ICA approaches. A similar behaviour was observed for cMEM, for which, however, Cont_ICA seemed less detrimental.

Conclusion: We demonstrated that a simplified ICA approach selecting one IC of interest in combination with distributed magnetic source imaging can be detrimental. More complex approaches may provide better results but would be rather difficult to apply in real-world clinical setting. In a broader perspective, caution should be taken in applying ICA for source localization of interictal activity. To ensure optimal and useful results, effort should focus on acquiring good quality data, minimizing artifacts, and determining optimal candidacy for MEG, rather than counting on data cleaning techniques.

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

Scalp ripples as prognostic biomarkers of epileptogenicity in pediatric surgery

Objective

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

Methods

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

Results

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

Interpretation

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

Magnetoencephalographic Mapping of Epileptic Spike Population Using Distributed Source Analysis: Comparison With Intracranial Electroencephalographic Spikes

INTRODUCTION

This study evaluates magnetoencephalographic (MEG) spike population as compared with intracranial electroencephalographic (IEEG) spikes using a quantitative method based on distributed source analysis.

METHODS

We retrospectively studied eight patients with medically intractable epilepsy who had an MEG and subsequent IEEG monitoring. Fifty MEG spikes were analyzed in each patient using minimum norm estimate. For individual spikes, each vertex in the source space was considered activated when its source amplitude at the peak latency was higher than a threshold, which was set at 50% of the maximum amplitude over all vertices. We mapped the total count of activation at each vertex. We also analyzed 50 IEEG spikes in the same manner over the intracranial electrodes and created the activation count map. The location of the electrodes was obtained in the MEG source space by coregistering postimplantation computed tomography to MRI. We estimated the MEG- and IEEG-active regions associated with the spike populations using the vertices/electrodes with a count over 25.

RESULTS

The activation count maps of MEG spikes demonstrated the localization associated with the spike population by variable count values at each vertex. The MEG-active region overlapped with 65 to 85% of the IEEG-active region in our patient group.

CONCLUSIONS

Mapping the MEG spike population is valid for demonstrating the trend of spikes clustering in patients with epilepsy. In addition, comparison of MEG and IEEG spikes quantitatively may be informative for understanding their relationship.

Virtual localization of the seizure onset zone: Using non-invasive MEG virtual electrodes at stereo-EEG electrode locations in refractory epilepsy patients

In some patients with medically refractory epilepsy, EEG with intracerebrally placed electrodes (stereo-electroencephalography, SEEG) is needed to locate the seizure onset zone (SOZ) for successful epilepsy surgery. SEEG has limitations and entails risk of complications because of its invasive character. Non-invasive magnetoencephalography virtual electrodes (MEG-VEs) may overcome SEEG limitations and optimize electrode placement making SOZ localization safer. Our purpose was to assess whether interictal activity measured by MEG-VEs and SEEG at identical anatomical locations were comparable, and whether MEG-VEs activity properties could determine the location of a later resected brain area (RA) as an approximation of the SOZ. We analyzed data from nine patients who underwent MEG and SEEG evaluation, and surgery for medically refractory epilepsy. MEG activity was retrospectively reconstructed using beamforming to obtain VEs at the anatomical locations corresponding to those of SEEG electrodes. Spectral, functional connectivity and functional network properties were obtained for both, MEG-VEs and SEEG time series, and their correlation and reliability were established. Based on these properties, the approximation of the SOZ was characterized by the differences between RA and non-RA (NRA). We found significant positive correlation and reliability between MEG-VEs and SEEG spectral measures (particularly in delta [0.5-4 Hz], alpha2 [10-13 Hz], and beta [13-30 Hz] bands) and broadband functional connectivity. Both modalities showed significantly slower activity and a tendency towards increased broadband functional connectivity in the RA compared to the NRA. Our findings show that spectral and functional connectivity properties of non-invasively obtained MEG-VEs match those of invasive SEEG recordings, and can characterize the SOZ. This suggests that MEG-VEs might be used for optimal SEEG planning and fewer depth electrode implantations, making the localization of the SOZ safer and more successful.

Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy

Fusion of electroencephalography (EEG) and magnetoencephalography (MEG) data using maximum entropy on the mean method (MEM-fusion) takes advantage of the complementarities between EEG and MEG to improve localization accuracy. Simulation studies demonstrated MEM-fusion to be robust especially in noisy conditions such as single spike source localizations (SSSL). Our objective was to assess the reliability of SSSL using MEM-fusion on clinical data. We proposed to cluster SSSL results to find the most reliable and consistent source map from the reconstructed sources, the so-called consensus map. Thirty-four types of interictal epileptic discharges (IEDs) were analyzed from 26 patients with well-defined epileptogenic focus. SSSLs were performed on EEG, MEG, and fusion data and consensus maps were estimated using hierarchical clustering. Qualitative (spike-to-spike reproducibility rate, SSR) and quantitative (localization error and spatial dispersion) assessments were performed using the epileptogenic focus as clinical reference. Fusion SSSL provided significantly better results than EEG or MEG alone. Fusion found at least one cluster concordant with the clinical reference in all cases. This concordant cluster was always the one involving the highest number of spikes. Fusion yielded highest reproducibility (SSR EEG = 55%, MEG = 71%, fusion = 90%) and lowest localization error. Also, using only few channels from either modality (21EEG + 272MEG or 54EEG + 25MEG) was sufficient to reach accurate fusion. MEM-fusion with consensus map approach provides an objective way of finding the most reliable and concordant generators of IEDs. We, therefore, suggest the pertinence of SSSL using MEM-fusion as a valuable clinical tool for presurgical evaluation of epilepsy.

Clinical yield of magnetoencephalography distributed source imaging in epilepsy: A comparison with equivalent current dipole method

OBJECTIVE

Source localization of interictal epileptic discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. It is usually obtained by equivalent current dipole (ECD) which localizes a point source and is the only inverse solution approved by clinical guidelines. In contrast, magnetic source imaging using distributed methods (dMSI) provides maps of the location and the extent of the generators, but its yield has not been clinically validated. We systematically compared ECD versus dMSI performed using coherent Maximum Entropy on the Mean (cMEM), a method sensitive to the spatial extent of the generators.

METHODS

340 source localizations of IEDs derived from 49 focal epilepsy patients with foci well-defined through intracranial EEG, MRI lesions, and surgery were analyzed. The comparison was based on the assessment of the sublobar concordance with the focus and of the distance between the source and the focus.

RESULTS

dMSI sublobar concordance was significantly higher than ECD (81% vs 69%, P < 0.001), especially for extratemporal lobe sources (dMSI = 84%; ECD = 67%, P < 0.001) and for seizure free patients (dMSI = 83%; ECD = 70%, P < 0.001). The median distance from the focus was 4.88 mm for ECD and 3.44 mm for dMSI (P < 0.001). ECD dipoles were often wrongly localized in deep brain regions.

CONCLUSIONS

dMSI using cMEM exhibited better accuracy. dMSI also offered the advantage of recovering more realistic maps of the generator, which could be exploited for neuronavigation aimed at targeting invasive EEG and surgical resection. Therefore, dMSI may be preferred to ECD in clinical practice. Hum Brain Mapp 39:218-231, 2018. © 2017 Wiley Periodicals, Inc.

Magnetoencephalography in Preoperative Epileptic Foci Localization: Enlightenment from Cognitive Studies

Over 30% epileptic patients are refractory to medication, who are amenable to neurosurgical treatment. Non-invasive brain imaging technologies including video-electroencephalogram (EEG), magnetic resonance imaging (MRI), and magnetoencephalography (MEG) are widely used in presurgical assessment of epileptic patients. This review mainly discussed the current development of clinical MEG imaging as a diagnose approach, and its correlations with the golden standard intracranial electroencephalogram (iEEG). More importantly, this review discussed the possible applications of functional networks in preoperative epileptic foci localization in future studies.

Automated long-term EEG analysis to localize the epileptogenic zone

Objective

We investigated the performance of automatic spike detection and subsequent electroencephalogram (EEG) source imaging to localize the epileptogenic zone (EZ) from long-term EEG recorded during video-EEG monitoring.

Methods

In 32 patients, spikes were automatically detected in the EEG and clustered according to their morphology. The two spike clusters with most single events in each patient were averaged and localized in the brain at the half-rising time and peak of the spike using EEG source imaging. On the basis of the distance from the sources to the resection and the known patient outcome after surgery, the performance of the automated EEG analysis to localize the EZ was quantified.

Results

In 28 out of the 32 patients, the automatically detected spike clusters corresponded with the reported interictal findings. The median distance to the resection in patients with Engel class I outcome was 6.5 and 15 mm for spike cluster 1 and 27 and 26 mm for cluster 2, at the peak and the half-rising time of the spike, respectively. Spike occurrence (cluster 1 vs. cluster 2) and spike timing (peak vs. half-rising) significantly influenced the distance to the resection (p < 0.05). For patients with Engel class II, III, and IV outcomes, the median distance increased to 36 and 36 mm for cluster 1. Localizing spike cluster 1 at the peak resulted in a sensitivity of 70% and specificity of 100%, positive prediction value (PPV) of 100%, and negative predictive value (NPV) of 53%. Including the results of spike cluster 2 led to an increased sensitivity of 79% NPV of 55% and diagnostic OR of 11.4, while the specificity dropped to 75% and the PPV to 90%.

Significance

We showed that automated analysis of long-term EEG recordings results in a high sensitivity and specificity to localize the epileptogenic focus.

How to record high-frequency oscillations in epilepsy: A practical guideline

OBJECTIVE

Technology for localizing epileptogenic brain regions plays a central role in surgical planning. Recent improvements in acquisition and electrode technology have revealed that high-frequency oscillations (HFOs) within the 80-500 Hz frequency range provide the neurophysiologist with new information about the extent of the epileptogenic tissue in addition to ictal and interictal lower frequency events. Nevertheless, two decades after their discovery there remain questions about HFOs as biomarkers of epileptogenic brain and there use in clinical practice.

METHODS

In this review, we provide practical, technical guidance for epileptologists and clinical researchers on recording, evaluation, and interpretation of ripples, fast ripples, and very high-frequency oscillations.

RESULTS

We emphasize the importance of low noise recording to minimize artifacts. HFO analysis, either visual or with automatic detection methods, of high fidelity recordings can still be challenging because of various artifacts including muscle, movement, and filtering. Magnetoencephalography and intracranial electroencephalography (iEEG) recordings are subject to the same artifacts.

SIGNIFICANCE

High-frequency oscillations are promising new biomarkers in epilepsy. This review provides interested researchers and clinicians with a review of current state of the art of recording and identification and potential challenges to clinical translation.

Comparison of the spatial resolution of source imaging techniques in high-density EEG and MEG

BACKGROUND

The present study aims at evaluating and comparing electrical and magnetic distributed source imaging methods applied to high-density Electroencephalography (hdEEG) and Magnetoencephalography (MEG) data. We used resolution matrices to characterize spatial resolution properties of Minimum Norm Estimate (MNE), dynamic Statistical Parametric Mapping (dSPM), standardized Low-Resolution Electromagnetic Tomography (sLORETA) and coherent Maximum Entropy on the Mean (cMEM, an entropy-based technique). The resolution matrix provides information of the Point Spread Functions (PSF) and of the Crosstalk functions (CT), this latter being also called source leakage, as it reflects the influence of a source on its neighbors.

METHODS

The spatial resolution of the inverse operators was first evaluated theoretically and then with real data acquired using electrical median nerve stimulation on five healthy participants. We evaluated the Dipole Localization Error (DLE) and the Spatial Dispersion (SD) of each PSF and CT map.

RESULTS

cMEM showed the smallest spatial spread (SD) for both PSF and CT maps, whereas localization errors (DLE) were similar for all methods. Whereas cMEM SD values were lower in MEG compared to hdEEG, the other methods slightly favored hdEEG over MEG. In real data, cMEM provided similar localization error and significantly less spatial spread than other methods for both MEG and hdEEG. Whereas both MEG and hdEEG provided very accurate localizations, all the source imaging methods actually performed better in MEG compared to hdEEG according to all evaluation metrics, probably due to the higher signal-to-noise ratio of the data in MEG.

CONCLUSION

Our overall results show that all investigated methods provide similar localization errors, suggesting very accurate localization for both MEG and hdEEG when similar number of sensors are considered for both modalities. Intrinsic properties of source imaging methods as well as their behavior for well-controlled tasks, suggest an overall better performance of cMEM in regards to spatial resolution and spatial leakage for both hdEEG and MEG. This indicates that cMEM would be a good candidate for studying source localization of focal and extended generators as well as functional connectivity studies.

Dynamic Electrical Source Imaging (DESI) of Seizures and Interictal Epileptic Discharges Without Ensemble Averaging

We propose an algorithm for electrical source imaging of epileptic discharges that takes a data-driven approach to regularizing the dynamics of solutions. The method is based on linear system identification on short time segments, combined with a classical inverse solution approach. Whereas ensemble averaging of segments or epochs discards inter-segment variations by averaging across them, our approach explicitly models them. Indeed, it may even be possible to avoid the need for the time-consuming process of marking epochs containing discharges altogether. We demonstrate that this approach can produce both stable and accurate inverse solutions in experiments using simulated data and real data from epilepsy patients. In an illustrative example, we show that we are able to image propagation using this approach. We show that when applied to imaging seizure data, our approach reproducibly localized frequent seizure activity to within the margins of surgeries that led to patients' seizure freedom. The same approach could be used in the planning of epilepsy surgeries, as a way to localize potentially epileptogenic tissue that should be resected.

Interictal High Frequency Oscillations Detected with Simultaneous Magnetoencephalography and Electroencephalography as Biomarker of Pediatric Epilepsy

Crucial to the success of epilepsy surgery is the availability of a robust biomarker that identifies the Epileptogenic Zone (EZ). High Frequency Oscillations (HFOs) have emerged as potential presurgical biomarkers for the identification of the EZ in addition to Interictal Epileptiform Discharges (IEDs) and ictal activity. Although they are promising to localize the EZ, they are not yet suited for the diagnosis or monitoring of epilepsy in clinical practice. Primary barriers remain: the lack of a formal and global definition for HFOs; the consequent heterogeneity of methodological approaches used for their study; and the practical difficulties to detect and localize them noninvasively from scalp recordings. Here, we present a methodology for the recording, detection, and localization of interictal HFOs from pediatric patients with refractory epilepsy. We report representative data of HFOs detected noninvasively from interictal scalp EEG and MEG from two children undergoing surgery.

The underlying generators of HFOs were localized by solving the inverse problem and their localization was compared to the Seizure Onset Zone (SOZ) as this was defined by the epileptologists. For both patients, Interictal Epileptogenic Discharges (IEDs) and HFOs were localized with source imaging at concordant locations. For one patient, intracranial EEG (iEEG) data were also available. For this patient, we found that the HFOs localization was concordant between noninvasive and invasive methods. The comparison of iEEG with the results from scalp recordings served to validate these findings. To our best knowledge, this is the first study that presents the source localization of scalp HFOs from simultaneous EEG and MEG recordings comparing the results with invasive recordings. These findings suggest that HFOs can be reliably detected and localized noninvasively with scalp EEG and MEG. We conclude that the noninvasive localization of interictal HFOs could significantly improve the presurgical evaluation for pediatric patients with epilepsy.

Source localization of the seizure onset zone from ictal EEG/MEG data

INTRODUCTION

Surgical treatment of drug-resistant epilepsy relies on the identification of the seizure onset zone (SOZ) and often requires intracranial EEG (iEEG). We have developed a new approach for non-invasive magnetic and electric source imaging of the SOZ (MSI-SOZ and ESI-SOZ) from ictal magnetoencephalography (MEG) and EEG recordings, using wavelet-based Maximum Entropy on the Mean (wMEM) method. We compared the performance of MSI-SOZ and ESI-SOZ with interictal spike source localization (MSI-spikes and ESI-spikes) and clinical localization of the SOZ (i.e., based on iEEG or lesion topography, denoted as clinical-SOZ).

METHODS

A total of 46 MEG or EEG seizures from 13 patients were analyzed. wMEM was applied around seizure onset, centered on the frequency band showing the strongest power change. Principal component analysis applied to spatiotemporal reconstructed wMEM sources (0.4-1 s around seizure onset) identified the main spatial pattern of ictal oscillations. Qualitative sublobar concordance and quantitative measures of distance and spatial overlaps were estimated to compare MSI/ESI-SOZ with MSI/ESI-Spikes and clinical-SOZ.

RESULTS

MSI/ESI-SOZ were concordant with clinical-SOZ in 81% of seizures (MSI 90%, ESI 64%). MSI-SOZ was more accurate and identified sources closer to the clinical-SOZ (P = 0.012) and to MSI-Spikes (P = 0.040) as compared with ESI-SOZ. MSI/ESI-SOZ and MSI/ESI-Spikes did not differ in terms of concordance and distance from the clinical-SOZ.

CONCLUSIONS

wMEM allows non-invasive localization of the SOZ from ictal MEG and EEG. MSI-SOZ performs better than ESI-SOZ. MSI/ESI-SOZ can provide important additional information to MSI/ESI-Spikes during presurgical evaluation. Hum Brain Mapp 37:2528-2546, 2016. © 2016 Wiley Periodicals, Inc.