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.

How cerebral cortex protects itself from interictal spikes: The alpha/beta inhibition mechanism.

Interactions between interictal epileptiform discharges (IEDs) and distant cortical regions subserve potential effects on cognition of patients with focal epilepsy. We hypothesize that "healthy" brain areas at a distance from the epileptic focus may respond to the interference of IEDs by generating inhibitory alpha and beta oscillations. We predict that more prominent alpha-beta oscillations can be found in patients with less impaired neurocognitive profile. We performed a source imaging magnetoencephalography study, including 41 focal epilepsy patients: 21 with frontal lobe epilepsy (FLE) and 20 with mesial temporal lobe epilepsy. We investigated the effect of anterior (i.e., frontal and temporal) IEDs on the oscillatory pattern over posterior head regions. We compared cortical oscillations (5-80 Hz) temporally linked to 3,749 IEDs (1,945 frontal and 1,803 temporal) versus an equal number of IED-free segments. We correlated results from IED triggered oscillations to global neurocognitive performance. Only frontal IEDs triggered alpha-beta oscillations over posterior head regions. IEDs with higher amplitude triggered alpha-beta oscillations of higher magnitude. The intensity of posterior head region alpha-beta oscillations significantly correlated with a better neuropsychological profile. Our study demonstrated that cerebral cortex protects itself from IEDs with generation of inhibitory alpha-beta oscillations at distant cortical regions. The association of more prominent oscillations with a better cognitive status suggests that this mechanism might play a role in determining the cognitive resilience in patients with FLE.

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.

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.

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.

Detection and Magnetic Source Imaging of Fast Oscillations (40-160 Hz) Recorded with Magnetoencephalography in Focal Epilepsy Patients.

We present a framework to detect fast oscillations (FOs) in magnetoencephalography (MEG) and to perform magnetic source imaging (MSI) to determine the location and extent of their generators in the cortex. FOs can be of physiologic origin associated to sensory processing and memory consolidation. In epilepsy, FOs are of pathologic origin and biomarkers of the epileptogenic zone. Seventeen patients with focal epilepsy previously confirmed with identified FOs in scalp electroencephalography (EEG) were evaluated. To handle data deriving from large number of sensors (275 axial gradiometers) we used an automatic detector with high sensitivity. False positives were discarded by two human experts. MSI of the FOs was performed with the wavelet based maximum entropy on the mean method. We found FOs in 11/17 patients, in only one patient the channel with highest FO rate was not concordant with the epileptogenic region and might correspond to physiologic oscillations. MEG FOs rates were very low: 0.02-4.55 per minute. Compared to scalp EEG, detection sensitivity was lower, but the specificity higher in MEG. MSI of FOs showed concordance or partial concordance with proven generators of seizures and epileptiform activity in 10/11 patients. We have validated the proposed framework for the non-invasive study of FOs with MEG. The excellent overall concordance with other clinical gold standard evaluation tools indicates that MEG FOs can provide relevant information to guide implantation for intracranial EEG pre-surgical evaluation and for surgical treatment, and demonstrates the important added value of choosing appropriate FOs detection and source localization methods.

Localization Accuracy of Distributed Inverse Solutions for Electric and Magnetic Source Imaging of Interictal Epileptic Discharges in Patients with Focal Epilepsy.

Distributed inverse solutions aim to realistically reconstruct the origin of interictal epileptic discharges (IEDs) from noninvasively recorded electroencephalography (EEG) and magnetoencephalography (MEG) signals. Our aim was to compare the performance of different distributed inverse solutions in localizing IEDs: coherent maximum entropy on the mean (cMEM), hierarchical Bayesian implementations of independent identically distributed sources (IID, minimum norm prior) and spatially coherent sources (COH, spatial smoothness prior). Source maxima (i.e., the vertex with the maximum source amplitude) of IEDs in 14 EEG and 19 MEG studies from 15 patients with focal epilepsy were analyzed. We visually compared their concordance with intracranial EEG (iEEG) based on 17 cortical regions of interest and their spatial dispersion around source maxima. Magnetic source imaging (MSI) maxima from cMEM were most often confirmed by iEEG (cMEM: 14/19, COH: 9/19, IID: 8/19 studies). COH electric source imaging (ESI) maxima co-localized best with iEEG (cMEM: 8/14, COH: 11/14, IID: 10/14 studies). In addition, cMEM was less spatially spread than COH and IID for ESI and MSI (p < 0.001 Bonferroni-corrected post hoc t test). Highest positive predictive values for cortical regions with IEDs in iEEG could be obtained with cMEM for MSI and with COH for ESI. Additional realistic EEG/MEG simulations confirmed our findings. Accurate spatially extended sources, as found in cMEM (ESI and MSI) and COH (ESI) are desirable for source imaging of IEDs because this might influence surgical decision. Our simulations suggest that COH and IID overestimate the spatial extent of the generators compared to cMEM.

MEG-EEG Information Fusion and Electromagnetic Source Imaging: From Theory to Clinical Application in Epilepsy.

The purpose of this study is to develop and quantitatively assess whether fusion of EEG and MEG (MEEG) data within the maximum entropy on the mean (MEM) framework increases the spatial accuracy of source localization, by yielding better recovery of the spatial extent and propagation pathway of the underlying generators of inter-ictal epileptic discharges (IEDs). The key element in this study is the integration of the complementary information from EEG and MEG data within the MEM framework. MEEG was compared with EEG and MEG when localizing single transient IEDs. The fusion approach was evaluated using realistic simulation models involving one or two spatially extended sources mimicking propagation patterns of IEDs. We also assessed the impact of the number of EEG electrodes required for an efficient EEG-MEG fusion. MEM was compared with minimum norm estimate, dynamic statistical parametric mapping, and standardized low-resolution electromagnetic tomography. The fusion approach was finally assessed on real epileptic data recorded from two patients showing IEDs simultaneously in EEG and MEG. Overall the localization of MEEG data using MEM provided better recovery of the source spatial extent, more sensitivity to the source depth and more accurate detection of the onset and propagation of IEDs than EEG or MEG alone. MEM was more accurate than the other methods. MEEG proved more robust than EEG and MEG for single IED localization in low signal-to-noise ratio conditions. We also showed that only few EEG electrodes are required to bring additional relevant information to MEG during MEM fusion.

Spatial correlation of hemodynamic changes related to interictal epileptic discharges with electric and magnetic source imaging.

INTRODUCTION: Blood oxygenation level-dependent (BOLD) signal changes at the time of interictal epileptic discharges (IEDs) identify their associated vascular/hemodynamic responses. BOLD activations and deactivations can be found within the epileptogenic zone but also at a distance. Source imaging identifies electric (ESI) and magnetic (MSI) sources of IEDs, with the advantage of a higher temporal resolution. Therefore, the objective of our study was to evaluate the spatial concordance between ESI/MSI and BOLD responses for similar IEDs. METHODS: Twenty-one patients with similar IEDs in simultaneous electroencephalogram/functional magnetic resonance imaging (EEG/fMRI) and in simultaneous EEG/magnetoencephalogram (MEG) recordings were studied. IEDs in EEG/fMRI acquisition were analyzed in an event-related paradigm within a general linear model (GLM). ESI/MSI of averaged IEDs was performed using the Maximum Entropy on the Mean. We assessed the spatial concordance between ESI/MSI and clusters of BOLD activations/deactivations with surface-based metrics. RESULTS: ESI/MSI were concordant with one BOLD cluster for 20/21 patients (concordance with activation: 14/21 patients, deactivation: 6/21 patients, no concordance: 1/21 patients; concordance with MSI only: 3/21, ESI only: 2/21). These BOLD clusters exhibited in 19/20 cases the most significant voxel. BOLD clusters that were spatially concordant with ESI/MSI were concordant with IEDs from invasive recordings in 8/11 patients (activations: 5/8, deactivations: 3/8). CONCLUSION: As the results of BOLD, ESI and MSI are often concordant, they reinforce our confidence in all of them. ESI and MSI confirm the most significant BOLD cluster within BOLD maps, emphasizing the importance of these clusters for the definition of the epileptic focus.

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

OBJECTIVES: Accurate delineation of the seizure-onset zone (SOZ) in focal drug-resistant epilepsy often requires stereo-electroencephalography (SEEG) recordings. We aimed at: (1) proposing a truly objective and quantitative comparison between electro-encephalography/magnetoencephalography (EEG/MEG) source-imaging (EMSI), EEG/functional MRI (EEG/fMRI) responses for similar spikes with primary-irritative zone (PIZ) and SOZ defined by SEEG and (2) evaluating 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 primary-irritative zone and seizure-onset zone 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 primary-irritative zone 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 seizure-onset zone was associated with smaller distances from the maximum  fMRIp-SEEG contact to the Spike-onset and Seizure-core channels (median distances 10 mm and 5mm respectively). Surgical outcomes were correctly predicted by EEG/MEG in 12/15 (80%) patients and EEG/fMRI in 6/11(54%) patients. CONCLUSIONS: Using a unique quantitative approach estimating EMSI and fMRI results in the reference SEEG channel-space, EEG/MEG and EEG/fMRI accurately localized the seizure-onset zone as well as the primary-irritative zone. Precisely, EEG/MEG more accurately localized the primary-irritative zone, whereas EEG/fMRI was more sensitive to the seizure-onset zone. Both neuro-imaging techniques provide complementary localization that can help guiding SEEG implantation and selecting good candidates for surgery.

Fast oscillations >40 Hz localize the epileptogenic zone: An electrical source imaging study using high-density electroencephalography.

OBJECTIVE: Fast Oscillations (FO) >40 Hz are a promising biomarker of the epileptogenic zone (EZ). Evidence using scalp electroencephalography (EEG) remains scarce. We assessed if electrical source imaging of FO using 256-channel high-density EEG (HD-EEG) is useful for EZ identification. METHODS: We analyzed HD-EEG recordings of 10 focal drug-resistant epilepsy patients with seizure-free postsurgical outcome. We marked FO candidate events at the time of epileptic spikes and verified them by screening for an isolated peak in the time-frequency plot. We performed electrical source imaging of spikes and FO within the Maximum Entropy of the Mean framework. Source localization maps were validated against the surgical cavity. RESULTS: We identified FO in five out of 10 patients who had a superficial or intermediate deep generator. The maximum of the FO maps was localized inside the cavity in all patients (100%). Analysis with a reduced electrode coverage using the 10-10 and 10-20 system showed a decreased localization accuracy of 60% and 40% respectively. CONCLUSIONS: FO recorded with HD-EEG localize the EZ. HD-EEG is better suited to detect and localize FO than conventional EEG approaches. SIGNIFICANCE: This study acts as proof-of-concept that FO localization using 256-channel HD-EEG is a viable marker of the EZ.

Spectral bandwidth of interictal fast epileptic activity characterizes the seizure onset zone.

The foremost aim of presurgical epilepsy evaluation is the delineation of the seizure onset zone (SOZ). There is increasing evidence that fast epileptic activity (FEA, 14-250 Hz) occurring interictally, i.e. between seizures, is predominantly localized within the SOZ. Currently it is unknown, which frequency band of FEA performs best in identifying the SOZ, although prior studies suggest highest concordance of spectral changes with the SOZ for high frequency changes. We suspected that FEA reflects dampened oscillations in local cortical excitatory-inhibitory neural networks, and that interictal FEA in the SOZ is a consequence of reduced oscillatory damping. We therefore predict a narrowing of the spectral bandwidth alongside increased amplitudes of spectral peaks during interictal FEA events. To test this hypothesis, we evaluated spectral changes during interictal FEA in invasive EEG (iEEG) recordings of 13 patients with focal epilepsy. In relative spectra of beta and gamma band changes (14-250 Hz) during FEA, we found that spectral peaks within the SOZ indeed were significantly more narrow-banded and their power changes were significantly higher than outside the SOZ. In contrast, the peak frequency did not differ within and outside the SOZ. Our results show that bandwidth and power changes of spectral modulations during FEA both help localizing the SOZ. We propose the spectral bandwidth as new source of information for the evaluation of EEG data.