Preoperative evaluation using magnetoencephalography: Experience in 382 epilepsy patients

Association Between Magnetoencephalography-Localized Epileptogenic Zone, Surgical Resection Volume, and Postsurgical Seizure Outcome

PURPOSE

Surgical resection of magnetoencephalography (MEG) dipole clusters, reconstructed from interictal epileptiform discharges, is associated with favorable seizure outcomes. However, the relation of MEG cluster resection to the surgical resection volume is not known nor is it clear whether this association is direct and causal, or it may be mediated by the resection volume or other predictive factors. This study aims to clarify these open questions and assess the diagnostic accuracy of MEG in our center.

METHODS

We performed a retrospective cohort study of 68 patients with drug-resistant epilepsy who underwent MEG followed by resective epilepsy surgery and had at least 12 months of postsurgical follow-up.

RESULTS

Good seizure outcomes were associated with monofocal localization (χ 2 = 6.94, P = 0.001; diagnostic odds ratio = 10.2) and complete resection of MEG clusters (χ 2 = 22.1, P < 0.001; diagnostic odds ratio = 42.5). Resection volumes in patients with and without removal of MEG clusters were not significantly different ( t = 0.18, P = 0.86; removed: M = 20,118 mm 3 , SD = 10,257; not removed: M = 19,566 mm 3 , SD = 10,703). Logistic regression showed that removal of MEG clusters predicts seizure-free outcome independent of the resection volume and other prognostic factors ( P < 0.001).

CONCLUSIONS

Complete resection of MEG clusters leads to favorable seizure outcomes without affecting the volume of surgical resection and independent of other prognostic factors. MEG can localize the epileptogenic zone with high accuracy. MEG interictal epileptiform discharges mapping should be used whenever feasible to improve postsurgical seizure outcomes.

From bench to bedside: Overview of magnetoencephalography in basic principle, signal processing, source localization and clinical applications

Magnetoencephalography (MEG) is a non-invasive technique that can precisely capture the dynamic spatiotemporal patterns of the brain by measuring the magnetic fields arising from neuronal activity along the order of milliseconds. Observations of brain dynamics have been used in cognitive neuroscience, the diagnosis of neurological diseases, and the brain-computer interface (BCI). In this study, we outline the basic principle, signal processing, and source localization of MEG, and describe its clinical applications for cognitive assessment, the diagnoses of neurological diseases and mental disorders, preoperative evaluation, and the BCI. This review not only provides an overall perspective of MEG, ranging from practical techniques to clinical applications, but also enhances the prevalent understanding of neural mechanisms. The use of MEG is expected to lead to significant breakthroughs in neuroscience.

Standardized hierarchical adaptive Lp regression for noise robust focal epilepsy source reconstructions

OBJECTIVE

To investigate the ability of standardization to reduce source localization errors and measurement noise uncertainties for hierarchical Bayesian algorithms with L1- and L2-norms as priors in electroencephalography and magnetoencephalography of focal epilepsy.

METHODS

Description of the standardization methodology relying on the Hierarchical Bayesian framework, referred to as the Standardized Hierarchical Adaptive Lp-norm Regularization (SHALpR). The performance was tested using real data from two focal epilepsy patients. Simulated data that resembled the available real data was constructed for further localization and noise robustness investigation.

RESULTS

The proposed algorithms were compared to their non-standardized counterparts, Standardized low-resolution brain electromagnetic tomography, Standardized Shrinking LORETA-FOCUSS, and Dynamic statistical parametric maps. Based on the simulations, the standardized Hierarchical adaptive algorithm using L2-norm was noise robust for 10 dB signal-to-noise ratio (SNR), whereas the L1-norm prior worked robustly also with 5 dB SNR. The accuracy of the standardized L1-normed methodology to localize focal activity was under 1 cm for both patients.

CONCLUSIONS

Numerical results of the proposed methodology display improved localization and noise robustness. The proposed methodology also outperformed the compared methods when dealing with real data.

SIGNIFICANCE

The proposed standardized methodology, especially when employing the L1-norm, could serve as a valuable assessment tool in surgical decision-making.

Non-invasive measurements of ictal and interictal epileptiform activity using optically pumped magnetometers

Magneto- and electroencephalography (MEG/EEG) are important techniques for the diagnosis and pre-surgical evaluation of epilepsy. Yet, in current cryogen-based MEG systems the sensors are offset from the scalp, which limits the signal-to-noise ratio (SNR) and thereby the sensitivity to activity from deep structures such as the hippocampus. This effect is amplified in children, for whom adult-sized fixed-helmet systems are typically too big. Moreover, ictal recordings with fixed-helmet systems are problematic because of limited movement tolerance and/or logistical considerations. Optically Pumped Magnetometers (OPMs) can be placed directly on the scalp, thereby improving SNR and enabling recordings during seizures. We aimed to demonstrate the performance of OPMs in a clinical population. Seven patients with challenging cases of epilepsy underwent MEG recordings using a 12-channel OPM-system and a 306-channel cryogen-based whole-head system: three adults with known deep or weak (low SNR) sources of interictal epileptiform discharges (IEDs), along with three children with focal epilepsy and one adult with frequent seizures. The consistency of the recorded IEDs across the two systems was assessed. In one patient the OPMs detected IEDs that were not found with the SQUID-system, and in two patients no IEDs were found with either system. For the other patients the OPM data were remarkably consistent with the data from the cryogenic system, noting that these were recorded in different sessions, with comparable SNRs and IED-yields overall. Importantly, the wearability of OPMs enabled the recording of seizure activity in a patient with hyperkinetic movements during the seizure. The observed ictal onset and semiology were in agreement with previous video- and stereo-EEG recordings. The relatively affordable technology, in combination with reduced running and maintenance costs, means that OPM-based MEG could be used more widely than current MEG systems, and may become an affordable alternative to scalp EEG, with the potential benefits of increased spatial accuracy, reduced sensitivity to volume conduction/field spread, and increased sensitivity to deep sources. Wearable MEG thus provides an unprecedented opportunity for epilepsy, and given its patient-friendliness, we envisage that it will not only be used for presurgical evaluation of epilepsy patients, but also for diagnosis after a first seizure.

Altered neuromagnetic activity in default mode network in childhood absence epilepsy

Purpose

The electrophysiological characterization of resting state oscillatory functional connectivity within the default mode network (DMN) during interictal periods in childhood absence epilepsy (CAE) remains unclear. Using magnetoencephalographic (MEG) recordings, this study investigated how the connectivity within the DMN was altered in CAE.

Methods

Using a cross-sectional design, we analyzed MEG data from 33 children newly diagnosed with CAE and 26 controls matched for age and sex. The spectral power and functional connectivity of the DMN were estimated using minimum norm estimation combined with the Welch technique and corrected amplitude envelope correlation.

Results

Default mode network showed stronger activation in the delta band during the ictal period, however, the relative spectral power in other bands was significantly lower than that in the interictal period (pcorrected &lt; 0.05 for DMN regions, except bilateral medial frontal cortex, left medial temporal lobe, left posterior cingulate cortex in the theta band, and the bilateral precuneus in the alpha band). It should be noted that the significant power peak in the alpha band was lost compared with the interictal data. Compared with controls, the interictal relative spectral power of DMN regions (except bilateral precuneus) in CAE patients was significantly increased in the delta band (pcorrected &lt; 0.01), whereas the values of all DMN regions in the beta-gamma 2 band were significantly decreased (pcorrected &lt; 0.01). In the higher frequency band (alpha-gamma1), especially in the beta and gamma1 band, the ictal node strength of DMN regions except the left precuneus was significantly higher than that in the interictal periods (pcorrected &lt; 0.01), and the node strength of the right inferior parietal lobe increased most significantly in the beta band (Ictal: 3.8712 vs. Interictal: 0.7503, pcorrected &lt; 0.01). Compared with the controls, the interictal node strength of DMN increased in all frequency bands, especially the right medial frontal cortex in the beta band (Controls: 0.1510 vs. Interictal: 3.527, pcorrected &lt; 0.01). Comparing relative node strength between groups, the right precuneus in CAE children decreased significantly (β: Controls: 0.1009 vs. Interictal: 0.0475; γ 1: Controls:0.1149 vs. Interictal:0.0587, pcorrected &lt; 0.01) such that it was no longer the central hub.

Conclusion

These findings indicated DMN abnormalities in CAE patients, even in interictal periods without interictal epileptic discharges. Abnormal functional connectivity in CAE may reflect abnormal anatomo-functional architectural integration in DMN, as a result of cognitive mental impairment and unconsciousness during absence seizure. Future studies are needed to examine if the altered functional connectivity can be used as a biomarker for treatment responses, cognitive dysfunction, and prognosis in CAE patients.

Pier K, Mullinger KJ, Marsh ED, Brookes MJ, Gaetz W. Mapping Interictal activity in epilepsy using a hidden Markov model: A magnetoencephalography study

Epilepsy is a highly heterogeneous neurological disorder with variable etiology, manifestation, and response to treatment. It is imperative that new models of epileptiform brain activity account for this variability, to identify individual needs and allow clinicians to curate personalized care. Here, we use a hidden Markov model (HMM) to create a unique statistical model of interictal brain activity for 10 pediatric patients. We use magnetoencephalography (MEG) data acquired as part of standard clinical care for patients at the Children's Hospital of Philadelphia. These data are routinely analyzed using excess kurtosis mapping (EKM); however, as cases become more complex (extreme multifocal and/or polymorphic activity), they become harder to interpret with EKM. We assessed the performance of the HMM against EKM for three patient groups, with increasingly complicated presentation. The difference in localization of epileptogenic foci for the two methods was 7 ± 2 mm (mean ± SD over all 10 patients); and 94% ± 13% of EKM temporal markers were matched by an HMM state visit. The HMM localizes epileptogenic areas (in agreement with EKM) and provides additional information about the relationship between those areas. A key advantage over current methods is that the HMM is a data-driven model, so the output is tuned to each individual. Finally, the model output is intuitive, allowing a user (clinician) to review the result and manually select the HMM epileptiform state, offering multiple advantages over previous methods and allowing for broader implementation of MEG epileptiform analysis in surgical decision-making for patients with intractable epilepsy.

Validating EEG, MEG and Combined MEG and EEG Beamforming for an Estimation of the Epileptogenic Zone in Focal Cortical Dysplasia

MEG and EEG source analysis is frequently used for the presurgical evaluation of pharmacoresistant epilepsy patients. The source localization of the epileptogenic zone depends, among other aspects, on the selected inverse and forward approaches and their respective parameter choices. In this validation study, we compare the standard dipole scanning method with two beamformer approaches for the inverse problem, and we investigate the influence of the covariance estimation method and the strength of regularization on the localization performance for EEG, MEG, and combined EEG and MEG. For forward modelling, we investigate the difference between calibrated six-compartment and standard three-compartment head modelling. In a retrospective study, two patients with focal epilepsy due to focal cortical dysplasia type IIb and seizure freedom following lesionectomy or radiofrequency-guided thermocoagulation (RFTC) used the distance of the localization of interictal epileptic spikes to the resection cavity resp. RFTC lesion as reference for good localization. We found that beamformer localization can be sensitive to the choice of the regularization parameter, which has to be individually optimized. Estimation of the covariance matrix with averaged spike data yielded more robust results across the modalities. MEG was the dominant modality and provided a good localization in one case, while it was EEG for the other. When combining the modalities, the good results of the dominant modality were mostly not spoiled by the weaker modality. For appropriate regularization parameter choices, the beamformer localized better than the standard dipole scan. Compared to the importance of an appropriate regularization, the sensitivity of the localization to the head modelling was smaller, due to similar skull conductivity modelling and the fixed source space without orientation constraint.

OUP accepted manuscript

Presurgical identification of the epileptogenic zone is a critical determinant of seizure control following surgical resection in epilepsy. Excitatory glutamate α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor is a major component of neurotransmission. Although elevated α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor levels are observed in surgically resected brain areas of patients with epilepsy, it remains unclear whether increased α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor-mediated currents initiate epileptic discharges. We have recently developed the first PET tracer for α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor, [¹¹C]K-2, to visualize and quantify the density of α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors in living human brains. Here, we detected elevated [¹¹C]K-2 uptake in the epileptogenic temporal lobe of patients with mesial temporal lobe epilepsy. Brain areas with high [¹¹C]K-2 uptake are closely colocalized with the location of equivalent current dipoles estimated by magnetoencephalography or with seizure onset zones detected by intracranial electroencephalogram. These results suggest that epileptic discharges initiate from brain areas with increased α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors, providing a biological basis for epileptic discharges and an additional non-invasive option to identify the epileptogenic zone in patients with mesial temporal lobe epilepsy.

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.

Identifying the epileptogenic zone by four non-invasive imaging techniques versus stereo-EEG in MRI-negative pre-surgery epilepsy patients

OBJECTIVE

We evaluated four imaging techniques, i.e. Electroencephalography (EEG)-functional Magnetic Resonance Imaging (MRI) (EEG-fMRI), High-resolution EEG (HR-EEG), Magnetoencephalography (MEG) and 2-[18F]fluoro-2-deoxy-D-glucose positron emission tomography (PET), for the identification of the epileptogenic zone (EZ) in 41 patients with negative MRI, candidate to neurosurgery.

METHODS

For each technique, results were compared to the Stereo-EEG. Diagnostic measures were calculated with respect to the post-surgical outcome, either for all the patients (39/41, two patients excluded) and for the subgroup of patients with the EZ involving more than one lobe (20/41).

RESULTS

When considered individually, each functional technique showed accuracy values ranging 54,6%-63,2%, having PET, MEG and HR-EEG higher sensitivity, and EEG-fMRI higher specificity. In patients with multilobar epileptogenic zone, functional techniques achieved the best accuracies (up to 80%) when three techniques, including EEG-fMRI, were considered together.

CONCLUSIONS

The study highlights the accuracy of a combination of functional imaging techniques in the identification of EZ in MRI negative focal epilepsy. The best diagnostic yield was obtained if the combination of PET, MEG (or HR-EEG as alternative), EEG-fMRI were considered together.

SIGNIFICANCE

The functional imaging techniques may improve the presurgical workup of MRI negative focal epilepsy, if epileptogenic zone involves more than one lobe.

Magnetoencephalography for epileptic focus localization in a series of 1000 cases

The aim of epilepsy surgery in patients with focal, pharmacoresistant epilepsies is to remove the complete epileptogenic zone to achieve long-term seizure freedom. In addition to a spectrum of diagnostic methods, magnetoencephalography focus localization is used for planning of epilepsy surgery. We present results from a retrospective observational cohort study of 1000 patients, evaluated using magnetoencephalography at the University Hospital Erlangen over the time span of 28 years. One thousand consecutive cases were included in the study, evaluated at the University Hospital Erlangen between 1990 and 2018. All patients underwent magnetoencephalography as part of clinical workup for epilepsy surgery. Of these, 405 underwent epilepsy surgery after magnetoencephalography, with postsurgical follow-ups of up to 20 years. Sensitivity for interictal epileptic activity was evaluated, in addition to concordance of localization with the consensus of presurgical workup on a lobar level. We evaluate magnetoencephalography characteristics of patients who underwent epilepsy surgery versus patients who did not proceed to surgery. In operated patients, resection of magnetoencephalography localizations were related to postsurgical seizure outcomes, including long-term results after several years. In comparison, association of lesionectomy with seizure outcomes was analysed. Measures of diagnostic accuracy were calculated for magnetoencephalography resection and lesionectomy. Sensitivity for interictal epileptic activity was 72% with significant differences between temporal and extra-temporal lobe epilepsy. Magnetoencephalography was concordant with the presurgical consensus in 51% and showed additional or more focal involvement in an additional 32%. Patients who proceeded to surgery showed a significantly higher percentage of monofocal magnetoencephalography results. Complete magnetoencephalography resection was associated with significantly higher chances to achieve seizure freedom in the short and long-term. Diagnostic accuracy was significant in temporal and extra-temporal lobe cases, but was significantly higher in extra-temporal lobe epilepsy (diagnostic odds ratios of 4.4 and 41.6). Odds ratios were also higher in non-lesional versus lesional cases (42.0 versus 6.2). The results show that magnetoencephalography provides non-redundant information, which significantly contributes to patient selection, focus localization and ultimately long-term seizure freedom after epilepsy surgery. Specifically in extra-temporal lobe epilepsy and non-lesional cases, magnetoencephalography provides excellent accuracy.

Assessing the localization accuracy and clinical utility of electric and magnetic source imaging in children with epilepsy

OBJECTIVE

To evaluate the accuracy and clinical utility of conventional 21-channel EEG (conv-EEG), 72-channel high-density EEG (HD-EEG) and 306-channel MEG in localizing interictal epileptiform discharges (IEDs).

METHODS

Twenty-four children who underwent epilepsy surgery were studied. IEDs on conv-EEG, HD-EEG, MEG and intracranial EEG (iEEG) were localized using equivalent current dipoles and dynamical statistical parametric mapping (dSPM). We compared the localization error (E_Loc) with respect to the ground-truth Irritative Zone (IZ), defined by iEEG sources, between non-invasive modalities and the distance from resection (D_res) between good- (Engel 1) and poor-outcomes. For each patient, we estimated the resection percentage of IED sources and tested whether it predicted outcome.

RESULTS

MEG presented lower E_Loc than HD-EEG and conv-EEG. For all modalities, D_res was shorter in good-outcome than poor-outcome patients, but only the resection percentage of the ground-truth IZ and MEG-IZ predicted surgical outcome.

CONCLUSIONS

MEG localizes the IZ more accurately than conv-EEG and HD-EEG. MSI may help the presurgical evaluation in terms of patient's outcome prediction. The promising clinical value of ESI for both conv-EEG and HD-EEG prompts the use of higher-density EEG-systems to possibly achieve MEG performance.

SIGNIFICANCE

Localizing the IZ non-invasively with MSI/ESI facilitates presurgical evaluation and surgical prognosis assessment.

Cognitive neuroscience using wearable magnetometer arrays: Non-invasive assessment of language function

Recent work has demonstrated that Optically Pumped Magnetometers (OPMs) can be utilised to create a wearable Magnetoencephalography (MEG) system that is motion robust. In this study, we use this system to map eloquent cortex using a clinically validated language lateralisation paradigm (covert verb generation: 120 trials, ∼10 min total duration) in healthy adults (n = 3). We show that it is possible to lateralise and localise language function on a case by case basis using this system. Specifically, we show that at a sensor and source level we can reliably detect a lateralising beta band (15-30 Hz) desynchronization in all subjects. This is the first study of human cognition using OPMs and not only highlights this technology's utility as tool for (developmental) cognitive neuroscience but also its potential to contribute to surgical planning via mapping of eloquent cortex, especially in young children.

Partial Least Square Aided Beamforming Algorithm in Magnetoencephalography Source Imaging

Beamforming techniques have played a prominent role in source imaging in neuroimaging and in locating epileptogenic zones. However, existing vector-beamformers are sensitive to noise on localization of epileptogenic zones. In this study, partial least square (PLS) was used to aid the minimum variance beamforming approach for source imaging with magnetoencephalography (MEG) arrays, and verified its effectiveness in simulated data and epilepsy data. First, PLS was employed to extract the components of the MEG arrays by maximizing the covariance between a linear combination of the predictors and the class variable. Noise was then removed by reconstructing the MEG arrays based on those components. The minimum variance beamforming method was used to estimate a source model. Simulations with a realistic head model and varying noise levels indicated that the proposed approach can provide higher spatial accuracy than other well-known beamforming methods. For real MEG recordings in 10 patients with temporal lobe epilepsy, the ratios of the number of spikes localized in the surgical excised region to the total number of spikes using the proposed method were higher than that of the dipole fitting method. These localization results using the proposed method are more consistent with the clinical evaluation. The proposed method may provide a new imaging marker for localization of epileptogenic zones.

Localization of the Epileptogenic Zone Using Interictal MEG and Machine Learning in a Large Cohort of Drug-Resistant Epilepsy Patients

Objective: Epilepsy surgery results in seizure freedom in the majority of drug-resistant patients. To improve surgery outcome we studied whether MEG metrics combined with machine learning can improve localization of the epileptogenic zone, thereby enhancing the chance of seizure freedom.

Methods: Presurgical interictal MEG recordings of 94 patients (64 seizure-free >1y post-surgery) were analyzed to extract four metrics in source space: delta power, low-to-high-frequency power ratio, functional connectivity (phase lag index), and minimum spanning tree betweenness centrality. At the group level, we estimated the overlap of the resection area with the five highest values for each metric and determined whether this overlap differed between surgery outcomes. At the individual level, those metrics were used in machine learning classifiers (linear support vector machine (SVM) and random forest) to distinguish between resection and non-resection areas and between surgery outcome groups.

Results: The highest values, for all metrics, overlapped with the resection area in more than half of the patients, but the overlap did not differ between surgery outcome groups. The classifiers distinguished the resection areas from non-resection areas with 59.94% accuracy (95% confidence interval: 59.67–60.22%) for SVM and 60.34% (59.98–60.71%) for random forest, but could not differentiate seizure-free from not seizure-free patients [43.77% accuracy (42.08–45.45%) for SVM and 49.03% (47.25–50.82%) for random forest].

Significance: All four metrics localized the resection area but did not distinguish between surgery outcome groups, demonstrating that metrics derived from interictal MEG correspond to expert consensus based on several presurgical evaluation modalities, but do not yet localize the epileptogenic zone. Metrics should be improved such that they correspond to the resection area in seizure-free patients but not in patients with persistent seizures. It is important to test such localization strategies at an individual level, for example by using machine learning or individualized models, since surgery is individually tailored.

An evaluation of kurtosis beamforming in magnetoencephalography to localize the epileptogenic zone in drug resistant epilepsy patients

•

Objective localizations of interictal spikes using a kurtosis beamformer.

•

Kurtosis Beamforming can provide confidence to scattered dipoles.

•

Kurtosis beamforming can assist in localizing the epileptogenic zone.

Objective

Kurtosis beamforming is a useful technique for analysing magnetoencephalograpy (MEG) data containing epileptic spikes. However, the implementation varies and few studies measure concordance with subsequently resected areas. We evaluated kurtosis beamforming as a means of localizing spikes in drug-resistant epilepsy patients.

Methods

We retrospectively applied kurtosis beamforming to MEG recordings of 22 epilepsy patients that had previously been analysed using equivalent current dipole (ECD) fitting. Virtual electrodes were placed in the kurtosis volumetric peaks and visually inspected to select a candidate source. The candidate sources were compared to the ECD localizations and resection areas.

Results

The kurtosis beamformer produced interpretable localizations in 18/22 patients, of which the candidate source coincided with the resection lobe in 9/13 seizure-free patients and in 3/5 patients with persistent seizures. The sublobar accuracy of the kurtosis beamformer with respect to the resection zone was higher than ECD (56% and 50%, respectively), however, ECD resulted in a higher lobar accuracy (75%, 67%).

Conclusions

Kurtosis beamforming may provide additional value when spikes are not clearly discernible on the sensors and support ECD localizations when dipoles are scattered.

Significance

Kurtosis beamforming should be integrated with existing clinical protocols to assist in localizing the epileptogenic zone.

Deep Source Localization with Magnetoencephalography Based on Sensor Array Decomposition and Beamforming

In recent years, the source localization technique of magnetoencephalography (MEG) has played a prominent role in cognitive neuroscience and in the diagnosis and treatment of neurological and psychological disorders. However, locating deep brain activities such as in the mesial temporal structures, especially in preoperative evaluation of epilepsy patients, may be more challenging. In this work we have proposed a modified beamforming approach for finding deep sources. First, an iterative spatiotemporal signal decomposition was employed for reconstructing the sensor arrays, which could characterize the intrinsic discriminant features for interpreting sensor signals. Next, a sensor covariance matrix was estimated under the new reconstructed space. Then, a well-known vector beamforming approach, which was a linearly constraint minimum variance (LCMV) approach, was applied to compute the solution for the inverse problem. It can be shown that the proposed source localization approach can give better localization accuracy than two other commonly-used beamforming methods (LCMV, MUSIC) in simulated MEG measurements generated with deep sources. Further, we applied the proposed approach to real MEG data recorded from ten patients with medically-refractory mesial temporal lobe epilepsy (mTLE) for finding epileptogenic zone(s), and there was a good agreement between those findings by the proposed approach and the clinical comprehensive results.

Automated Detection of Epileptic Biomarkers in Resting-State Interictal MEG Data

Certain differences between brain networks of healthy and epilectic subjects have been reported even during the interictal activity, in which no epileptic seizures occur. Here, magnetoencephalography (MEG) data recorded in the resting state is used to discriminate between healthy subjects and patients with either idiopathic generalized epilepsy or frontal focal epilepsy. Signal features extracted from interictal periods without any epileptiform activity are used to train a machine learning algorithm to draw a diagnosis. This is potentially relevant to patients without frequent or easily detectable spikes. To analyze the data, we use an up-to-date machine learning algorithm and explore the benefits of including different features obtained from the MEG data as inputs to the algorithm. We find that the relative power spectral density of the MEG time-series is sufficient to distinguish between healthy and epileptic subjects with a high prediction accuracy. We also find that a combination of features such as the phase-locked value and the relative power spectral density allow to discriminate generalized and focal epilepsy, when these features are calculated over a filtered version of the signals in certain frequency bands. Machine learning algorithms are currently being applied to the analysis and classification of brain signals. It is, however, less evident to identify the proper features of these signals that are prone to be used in such machine learning algorithms. Here, we evaluate the influence of the input feature selection on a clinical scenario to distinguish between healthy and epileptic subjects. Our results indicate that such distinction is possible with a high accuracy (86%), allowing the discrimination between idiopathic generalized and frontal focal epilepsy types.

Current and Emerging Potential of Magnetoencephalography in the Detection and Localization of High-Frequency Oscillations in Epilepsy

Up to one-third of patients with epilepsy are medically intractable and need resective surgery. To be successful, epilepsy surgery requires a comprehensive preoperative evaluation to define the epileptogenic zone (EZ), the brain area that should be resected to achieve seizure freedom. Due to lack of tools and methods that measure the EZ directly, this area is defined indirectly based on concordant data from a multitude of presurgical non-invasive tests and intracranial recordings. However, the results of these tests are often insufficiently concordant or inconclusive. Thus, the presurgical evaluation of surgical candidates is frequently challenging or unsuccessful. To improve the efficacy of the surgical treatment, there is an overriding need for reliable biomarkers that can delineate the EZ. High-frequency oscillations (HFOs) have emerged over the last decade as new potential biomarkers for the delineation of the EZ. Multiple studies have shown that HFOs are spatially associated with the EZ. Despite the encouraging findings, there are still significant challenges for the translation of HFOs as epileptogenic biomarkers to the clinical practice. One of the major barriers is the difficulty to detect and localize them with non-invasive techniques, such as magnetoencephalography (MEG) or scalp electroencephalography (EEG). Although most literature has studied HFOs using invasive recordings, recent studies have reported the detection and localization of HFOs using MEG or scalp EEG. MEG seems to be particularly advantageous compared to scalp EEG due to its inherent advantages of being less affected by skull conductivity and less susceptible to contamination from muscular activity. The detection and localization of HFOs with MEG would largely expand the clinical utility of these new promising biomarkers to an earlier stage in the diagnostic process and to a wider range of patients with epilepsy. Here, we conduct a thorough critical review of the recent MEG literature that investigates HFOs in patients with epilepsy, summarizing the different methodological approaches and the main findings. Our goal is to highlight the emerging potential of MEG in the non-invasive detection and localization of HFOs for the presurgical evaluation of patients with medically refractory epilepsy (MRE).

Identifying the epileptogenic zone in interictal resting-state MEG source-space networks

OBJECTIVE

In one third of patients, seizures remain after epilepsy surgery, meaning that improved preoperative evaluation methods are needed to identify the epileptogenic zone. A potential framework for such a method is network theory, as it can be applied to noninvasive recordings, even in the absence of epileptiform activity. Our aim was to identify the epileptogenic zone on the basis of hub status of local brain areas in interictal magnetoencephalography (MEG) networks.

METHODS

Preoperative eyes-closed resting-state MEG recordings were retrospectively analyzed in 22 patients with refractory epilepsy, of whom 14 were seizure-free 1 year after surgery. Beamformer-based time series were reconstructed for 90 cortical and subcortical automated anatomic labeling (AAL) regions of interest (ROIs). Broadband functional connectivity was estimated using the phase lag index in artifact-free epochs without interictal epileptiform abnormalities. A minimum spanning tree was generated to represent the network, and the hub status of each ROI was calculated using betweenness centrality, which indicates the centrality of a node in a network. The correspondence of resection cavity to hub values was evaluated on four levels: resection cavity, lobar, hemisphere, and temporal versus extratemporal areas.

RESULTS

Hubs were localized within the resection cavity in 8 of 14 seizure-free patients and in zero of 8 patients who were not seizure-free (57% sensitivity, 100% specificity, 73% accuracy). Hubs were localized in the lobe of resection in 9 of 14 seizure-free patients and in zero of 8 patients who were not seizure-free (64% sensitivity, 100% specificity, 77% accuracy). For the other two levels, the true negatives are unknown; hence, only sensitivity could be determined: hubs coincided with both the resection hemisphere and the resection location (temporal versus extratemporal) in 11 of 14 seizure-free patients (79% sensitivity).

SIGNIFICANCE

Identifying hubs noninvasively before surgery is a valuable approach with the potential of indicating the epileptogenic zone in patients without interictal abnormalities.