Imaging the extent and location of spatiotemporally distributed epileptiform sources from MEG measurements

Magnetoencephalography for Epilepsy Presurgical Evaluation

PURPOSE OF THE REVIEW

Magnetoencephalography (MEG) is a functional neuroimaging technique that records neurophysiology data with millisecond temporal resolution and localizes it with subcentimeter accuracy. Its capability to provide high resolution in both of these domains makes it a powerful tool both in basic neuroscience as well as clinical applications. In neurology, it has proven useful in its ability to record and localize epileptiform activity. Epilepsy workup typically begins with scalp electroencephalography (EEG), but in many situations, EEG-based localization of the epileptogenic zone is inadequate. The complementary sensitivity of MEG can be crucial in such cases, and MEG has been adopted at many centers as an important resource in building a surgical hypothesis. In this paper, we review recent work evaluating the extent of MEG influence of presurgical evaluations, novel analyses of MEG data employed in surgical workup, and new MEG instrumentation that will likely affect the field of clinical MEG.

RECENT FINDINGS

MEG consistently contributes to presurgical evaluation and these contributions often change the plan for epilepsy surgery. Extensive work has been done to develop new analytic methods for localizing the source of epileptiform activity with MEG. Systems using optically pumped magnetometry (OPM) have been successfully deployed to record and localize epileptiform activity. MEG remains an important noninvasive tool for epilepsy presurgical evaluation. Continued improvements in analytic methodology will likely increase the diagnostic yield of the test. Novel instrumentation with OPM may contribute to this as well, and may increase accessibility of MEG by decreasing cost.

Disentanglement of Resting State Brain Networks for Localizing Epileptogenic Zone in Focal Epilepsy

The objective of this study is to extract pathological brain networks from interictal period of E/MEG recordings to localize epileptic foci for presurgical evaluation. We proposed here a resting state E/MEG analysis framework, to disentangle brain functional networks represented by neural oscillations. By using an Embedded Hidden Markov Model, we constructed a state space for resting state recordings consisting of brain states with different spatiotemporal patterns. Functional connectivity analysis along with graph theory was applied on the extracted brain states to quantify the network features of the extracted brain states, based on which the source location of pathological states is determined. The method is evaluated by computer simulations and our simulation results revealed the proposed framework can extract brain states with high accuracy regarding both spatial and temporal profiles. We further evaluated the framework as compared with intracranial EEG defined seizure onset zone in 10 patients with drug-resistant focal epilepsy who underwent MEG recordings and were seizure free after surgical resection. The real patient data analysis showed very good localization results using the extracted pathological brain states in 6/10 patients, with localization error of about 15 mm as compared to the seizure onset zone. We show that the pathological brain networks can be disentangled from the resting-state electromagnetic recording and could be identified based on the connectivity features. The framework can serve as a useful tool in extracting brain functional networks from noninvasive resting state electromagnetic recordings, and promises to offer an alternative to aid presurgical evaluation guiding intracranial EEG electrodes implantation.

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

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

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.