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.

Consistency of electrical source imaging in presurgical evaluation of epilepsy across different vigilance states.

OBJECTIVE: The use of electrical source imaging (ESI) in assessing the source of interictal epileptic discharges (IEDs) is gaining increasing popularity in presurgical work-up of patients with drug-resistant focal epilepsy. While vigilance affects the ability to locate IEDs and identify the epileptogenic zone, we know little about its impact on ESI. METHODS: We studied overnight high-density electroencephalography recordings in focal drug-resistant epilepsy. IEDs were marked visually in each vigilance state, and examined in the sensor and source space. ESIs were calculated and compared between all vigilance states and the clinical ground truth. Two conditions were considered within each vigilance state, an unequalized and an equalized number of IEDs. RESULTS: The number, amplitude, and duration of IEDs were affected by the vigilance state, with N3 sleep presenting the highest number, amplitude, and duration for both conditions (P < 0.001), while signal-to-noise ratio only differed in the unequalized condition (P < 0.001). The vigilance state did not affect channel involvement (P > 0.05). ESI maps showed no differences in distance, quality, extent, or maxima distances compared to the clinical ground truth for both conditions (P > 0.05). Only when an absolute reference (wakefulness) was used, the channel involvement (P < 0.05) and ESI source extent (P < 0.01) were impacted during rapid-eye-movement (REM) sleep. Clustering of amplitude-sensitive and -insensitive ESI maps pointed to amplitude rather than the spatial profile as the driver (P < 0.05). INTERPRETATION: IED ESI results are stable across vigilance states, including REM sleep, if controlled for amplitude and IED number. ESI is thus stable and invariant to the vigilance state.

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.

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.