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

Emotion brain network topology in healthy subjects following passive listening to different auditory stimuli

A large body of research establishes the efficacy of musical intervention in many aspects of physical, cognitive, communication, social, and emotional rehabilitation. However, the underlying neural mechanisms for musical therapy remain elusive. This study aimed to investigate the potential neural correlates of musical therapy, focusing on the changes in the topology of emotion brain network. To this end, a Bayesian statistical approach and a cross-over experimental design were employed together with two resting-state magnetoencephalography (MEG) as controls. MEG recordings of 30 healthy subjects were acquired while listening to five auditory stimuli in random order. Two resting-state MEG recordings of each subject were obtained, one prior to the first stimulus (pre) and one after the final stimulus (post). Time series at the level of brain regions were estimated using depth-weighted minimum norm estimation (wMNE) source reconstruction method and the functional connectivity between these regions were computed. The resultant connectivity matrices were used to derive two topological network measures: transitivity and global efficiency which are important in gauging the functional segregation and integration of brain network respectively. The differences in these measures between pre- and post-stimuli resting MEG were set as the equivalence regions. We found that the network measures under all auditory stimuli were equivalent to the resting state network measures in all frequency bands, indicating that the topology of the functional brain network associated with emotional regulation in healthy subjects remains unchanged following these auditory stimuli. This suggests that changes in the emotion network topology may not be the underlying neural mechanism of musical therapy. Nonetheless, further studies are required to explore the neural mechanisms of musical interventions especially in the populations with neuropsychiatric disorders.

Interictal magnetoencephalography abnormalities to guide intracranial electrode implantation and predict surgical outcome

Intracranial EEG is the gold standard technique for epileptogenic zone localization but requires a preconceived hypothesis of the location of the epileptogenic tissue. This placement is guided by qualitative interpretations of seizure semiology, MRI, EEG and other imaging modalities, such as magnetoencephalography. Quantitative abnormality mapping using magnetoencephalography has recently been shown to have potential clinical value. We hypothesized that if quantifiable magnetoencephalography abnormalities were sampled by intracranial EEG, then patients' post-resection seizure outcome may be better. Thirty-two individuals with refractory neocortical epilepsy underwent magnetoencephalography and subsequent intracranial EEG recordings as part of presurgical evaluation. Eyes-closed resting-state interictal magnetoencephalography band power abnormality maps were derived from 70 healthy controls as a normative baseline. Magnetoencephalography abnormality maps were compared to intracranial EEG electrode implantation, with the spatial overlap of intracranial EEG electrode placement and cerebral magnetoencephalography abnormalities recorded. Finally, we assessed if the implantation of electrodes in abnormal tissue and subsequent resection of the strongest abnormalities determined by magnetoencephalography and intracranial EEG corresponded to surgical success. We used the area under the receiver operating characteristic curve as a measure of effect size. Intracranial electrodes were implanted in brain tissue with the most abnormal magnetoencephalography findings-in individuals that were seizure-free postoperatively (T = 3.9, P = 0.001) but not in those who did not become seizure-free. The overlap between magnetoencephalography abnormalities and electrode placement distinguished surgical outcome groups moderately well (area under the receiver operating characteristic curve = 0.68). In isolation, the resection of the strongest abnormalities as defined by magnetoencephalography and intracranial EEG separated surgical outcome groups well, area under the receiver operating characteristic curve = 0.71 and area under the receiver operating characteristic curve = 0.74, respectively. A model incorporating all three features separated surgical outcome groups best (area under the receiver operating characteristic curve = 0.80). Intracranial EEG is a key tool to delineate the epileptogenic zone and help render individuals seizure-free postoperatively. We showed that data-driven abnormality maps derived from resting-state magnetoencephalography recordings demonstrate clinical value and may help guide electrode placement in individuals with neocortical epilepsy. Additionally, our predictive model of postoperative seizure freedom, which leverages both magnetoencephalography and intracranial EEG recordings, could aid patient counselling of expected outcome.

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.

Non-invasive mapping of epileptogenic networks predicts surgical outcome

Epilepsy is increasingly considered a disorder of brain networks. Studying these networks with functional connectivity can help identify hubs that facilitate the spread of epileptiform activity. Surgical resection of these hubs may lead patients who suffer from drug-resistant epilepsy to seizure freedom. Here, we aim to map non-invasively epileptogenic networks, through the virtual implantation of sensors estimated with electric and magnetic source imaging, in patients with drug-resistant epilepsy. We hypothesize that highly connected hubs identified non-invasively with source imaging can predict the epileptogenic zone and the surgical outcome better than spikes localized with conventional source localization methods (dipoles). We retrospectively analysed simultaneous high-density electroencephalography (EEG) and magnetoencephalography data recorded from 37 children and young adults with drug-resistant epilepsy who underwent neurosurgery. Using source imaging, we estimated virtual sensors at locations where intracranial EEG contacts were placed. On data with and without spikes, we computed undirected functional connectivity between sensors/contacts using amplitude envelope correlation and phase locking value for physiologically relevant frequency bands. From each functional connectivity matrix, we generated an undirected network containing the strongest connections within sensors/contacts using the minimum spanning tree. For each sensor/contact, we computed graph centrality measures. We compared functional connectivity and their derived graph centrality of sensors/contacts inside resection for good (n = 22, ILAE I) and poor (n = 15, ILAE II-VI) outcome patients, tested their ability to predict the epileptogenic zone in good-outcome patients, examined the association between highly connected hubs removal and surgical outcome and performed leave-one-out cross-validation to support their prognostic value. We also compared the predictive values of functional connectivity with those of dipoles. Finally, we tested the reliability of virtual sensor measures via Spearman's correlation with intracranial EEG at population- and patient-level. We observed higher functional connectivity inside than outside resection (P < 0.05, Wilcoxon signed-rank test) for good-outcome patients, on data with and without spikes across different bands for intracranial EEG and electric/magnetic source imaging and few differences for poor-outcome patients. These functional connectivity measures were predictive of both the epileptogenic zone and outcome (positive and negative predictive values ≥55%, validated using leave-one-out cross-validation) outperforming dipoles on spikes. Significant correlations were found between source imaging and intracranial EEG measures (0.4 ≤ rho ≤ 0.9, P < 0.05). Our findings suggest that virtual implantation of sensors through source imaging can non-invasively identify highly connected hubs in patients with drug-resistant epilepsy, even in the absence of frank epileptiform activity. Surgical resection of these hubs predicts outcome better than dipoles.

Reduced connectivity of primary auditory and motor cortices during exposure to auditory white noise

Auditory white noise (WN) is widely used in daily life for inducing sleep, and in neuroscience to mask unwanted environmental noise and cues. However, WN was recently reported to influence corticospinal excitability and behavioral performance. Here, we expand previous preliminary findings on the influence of WN exposure on cortical functioning, and we hypothesize that it may modulate cortical connectivity. We tested our hypothesis by performing magnetoencephalography in 20 healthy subjects. WN reduces cortical connectivity of the primary auditory and motor regions with very distant cortical areas, showing a right lateralized connectivity reduction for primary motor cortex. The present results, together with previous finding concerning WN impact on corticospinal excitability and behavioral performance, further support the role of WN as a modulator of cortical function. This suggest avoiding its unrestricted use as a masking tool, while purposely designed and controlled WN application could be exploited to harness brain function and to treat neuropsychiatric conditions.

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.

MEG Node Degree for Focus Localization: Comparison with Invasive EEG

Epilepsy surgery is a viable therapy option for patients with pharmacoresistant focal epilepsies. A prerequisite for postoperative seizure freedom is the localization of the epileptogenic zone, e.g., using electro- and magnetoencephalography (EEG/MEG). Evidence shows that resting state MEG contains subtle alterations, which may add information to the workup of epilepsy surgery. Here, we investigate node degree (ND), a graph-theoretical parameter of functional connectivity, in relation to the seizure onset zone (SOZ) determined by invasive EEG (iEEG) in a consecutive series of 50 adult patients. Resting state data were subjected to whole brain, all-to-all connectivity analysis using the imaginary part of coherence. Graphs were described using parcellated ND. SOZ localization was investigated on a lobar and sublobar level. On a lobar level, all frequency bands except alpha showed significantly higher maximal ND (mND) values inside the SOZ compared to outside (ratios 1.11-1.20, alpha 1.02). Area-under-the-curve (AUC) was 0.67-0.78 for all expected alpha (0.44, ns). On a sublobar level, mND inside the SOZ was higher for all frequency bands (1.13-1.38, AUC 0.58-0.78) except gamma (1.02). MEG ND is significantly related to SOZ in delta, theta and beta bands. ND may provide new localization tools for presurgical evaluation of epilepsy surgery.

Auditory driven gamma synchrony is associated with cortical thickness in widespread cortical areas

OBJECTIVE

Gamma synchrony is a fundamental functional property of the cerebral cortex, impaired in multiple neuropsychiatric conditions (i.e. schizophrenia, Alzheimer's disease, stroke etc.). Auditory stimulation in the gamma range allows to drive gamma synchrony of the entire cortical mantle and to estimate the efficiency of the mechanisms sustaining it. As gamma synchrony depends strongly on the interplay between parvalbumin-positive interneurons and pyramidal neurons, we hypothesize an association between cortical thickness and gamma synchrony. To test this hypothesis, we employed a combined magnetoencephalography (MEG) - Magnetic Resonance Imaging (MRI) study.

METHODS

Cortical thickness was estimated from anatomical MRI scans. MEG measurements related to exposure of 40 Hz amplitude modulated tones were projected onto the cortical surface. Two measures of cortical synchrony were considered: (a) inter-trial phase consistency at 40 Hz, providing a vertex-wise estimation of gamma synchronization, and (b) phase-locking values between primary auditory cortices and whole cortical mantle, providing a measure of long-range cortical synchrony. A correlation between cortical thickness and synchronization measures was then calculated for 72 MRI-MEG scans.

RESULTS

Both inter-trial phase consistency and phase locking values showed a significant positive correlation with cortical thickness. For inter-trial phase consistency, clusters of strong associations were found in the temporal and frontal lobes, especially in the bilateral auditory and pre-motor cortices. Higher phase-locking values corresponded to higher cortical thickness in the frontal, temporal, occipital and parietal lobes.

DISCUSSION AND CONCLUSIONS

In healthy subjects, a thicker cortex corresponds to higher gamma synchrony and connectivity in the primary auditory cortex and beyond, likely reflecting underlying cell density involved in gamma circuitries. This result hints towards an involvement of gamma synchrony together with underlying brain structure in brain areas for higher order cognitive functions. This study contributes to the understanding of inherent cortical functional and structural brain properties, which might in turn constitute the basis for the definition of useful biomarkers in patients showing aberrant gamma synchronization.

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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {cm}^2$$\end{document}cm2 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.

Resting state network connectivity is attenuated by fMRI acoustic noise

INTRODUCTION

During the past decades there has been an increasing interest in tracking brain network fluctuations in health and disease by means of resting state functional magnetic resonance imaging (rs-fMRI). Rs-fMRI however does not provide the ideal environmental setting, as participants are continuously exposed to noise generated by MRI coils during acquisition of Echo Planar Imaging (EPI). We investigated the effect of EPI noise on resting state activity and connectivity using magnetoencephalography (MEG), by reproducing the acoustic characteristics of rs-fMRI environment during the recordings. As compared to fMRI, MEG has little sensitivity to brain activity generated in deep brain structures, but has the advantage to capture both the dynamic of cortical magnetic oscillations with high temporal resolution and the slow magnetic fluctuations highly correlated with BOLD signal.

METHODS

Thirty healthy subjects were enrolled in a counterbalanced design study including three conditions: a) silent resting state (Silence), b) resting state upon EPI noise (fMRI), and c) resting state upon white noise (White). White noise was employed to test the specificity of fMRI noise effect. The amplitude envelope correlation (AEC) in alpha band measured the connectivity of seven Resting State Networks (RSN) of interest (default mode network, dorsal attention network, language, left and right auditory and left and right sensory-motor). Vigilance dynamic was estimated from power spectral activity.

RESULTS

fMRI and White acoustic noise consistently reduced connectivity of cortical networks. The effects were widespread, but noise and network specificities were also present. For fMRI noise, decreased connectivity was found in the right auditory and sensory-motor networks. Progressive increase of slow theta-delta activity related to drowsiness was found in all conditions, but was significantly higher for fMRI . Theta-delta significantly and positively correlated with variations of cortical connectivity.

DISCUSSION

rs-fMRI connectivity is biased by unavoidable environmental factors during scanning, which warrant more careful control and improved experimental designs. MEG is free from acoustic noise and allows a sensitive estimation of resting state connectivity in cortical areas. Although underutilized, MEG could overcome issues related to noise during fMRI, in particular when investigation of motor and auditory networks is needed.

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.

Changes in the Functional Brain Network of Children Undergoing Repeated Epilepsy Surgery: An EEG Source Connectivity Study

About 30% of children with drug-resistant epilepsy (DRE) continue to have seizures after epilepsy surgery. Since epilepsy is increasingly conceptualized as a network disorder, understanding how brain regions interact may be critical for planning re-operation in these patients. We aimed to estimate functional brain connectivity using scalp EEG and its evolution over time in patients who had repeated surgery (RS-group, n = 9) and patients who had one successful surgery (seizure-free, SF-group, n = 12). We analyzed EEGs without epileptiform activity at varying time points (before and after each surgery). We estimated functional connectivity between cortical regions and their relative centrality within the network. We compared the pre- and post-surgical centrality of all the non-resected (untouched) regions (far or adjacent to resection) for each group (using the Wilcoxon signed rank test). In alpha, theta, and beta frequency bands, the post-surgical centrality of the untouched cortical regions increased in the SF group (p < 0.001) whereas they decreased (p < 0.05) or did not change (p > 0.05) in the RS group after failed surgeries; when re-operation was successful, the post-surgical centrality of far regions increased (p < 0.05). Our data suggest that removal of the epileptogenic focus in children with DRE leads to a gain in the network centrality of the untouched areas. In contrast, unaltered or decreased connectivity is seen when seizures persist after surgery.

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.

Evaluation of Brain Network Properties in Patients with MRI-Negative Temporal Lobe Epilepsy: An MEG Study

Abnormal functional brain networks of temporal lobe epilepsy (TLE) patients with structural abnormalities may partially reflect structural lesions rather than either TLE per se or functional compensatory processes. In this study, we sought to investigate the brain-network properties of intractable TLE patients apart from the effects of structural abnormalities. The brain network properties of 20 left and 23 right MRI-negative TLE patients and 22 healthy controls were evaluated using magnetoencephalographic recordings in six main frequency bands. A slowing of oscillatory brain activity was observed for the left or right TLE group vs. healthy controls. The TLE groups presented significantly increased functional connectivity in the delta, theta, lower alpha and beta bands, and significantly greater values in the normalized clustering coefficient and path length, and significantly smaller values in the weighted small-world measure in the theta band when compared to healthy controls. Alterations in global and regional band powers can be attributed to spectral slowing in TLE patients. The brain networks of TLE patients displayed abnormally high synchronization in multi-frequency bands and shifted toward a more regular architecture with worse network efficiency in the theta band. Without the contamination of structural lesions, these significant findings can be helpful for better understanding of the pathophysiological mechanism of TLE. The theta band can be considered as a preferred frequency band for investigating the brain-network dysfunction of MRI-negative intractable TLE patients.

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.

Efficient high resolution sLORETA in brain source localization

Objective. Estimation of the source location within the brain from electroencephalography (EEG) and magnetoencephalography measures is a challenging task. Among the existing techniques in the field, which are known as brain imaging methods, standardized low-resolution brain electromagnetic tomography (sLORETA) is the most popular method due to its simplicity and high accuracy. However, in this work we illustrate that sLORETA is still noisy and the additive noise is causing the blurry image. The existing pre-fixed/manual thresholding process after sLORETA can partially take care of denoising. However, this ad-hoc theresholding can either remove so much of the desired data or leave much of the noise in the process. Manual correction to avoid such extreme cases can be time-consuming. The objective of this paper is to automate the denoising process in the form of adaptive thresholding. Approach. The proposed method, denoted by efficient high-resolution sLORETA (EHR-sLORETA), is based on minimizing the error between the desired denoised source and the source estimates. Main results. The approach is evaluated using synthetic EEG and real EEG data. spatial dispersion (SD), and mean square error (MSE) are used as metrics to provide the quantitative performance of the method. In addition, qualitative analysis of the method is provided for real EEG data. This proposed model demonstrates advantages over the existing methods in sense of accuracy and robustness with SD and MSE comparison. Significance. EHR-sLORETA could have a significant impact on clinical studies with source estimation task, as it improves the accuracy of source estimation and eliminates the need for manual thresholding.

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.

Removal of Interictal MEG-Derived Network Hubs Is Associated With Postoperative Seizure Freedom

Objective: To investigate whether MEG network connectivity was associated with epilepsy duration, to identify functional brain network hubs in patients with refractory focal epilepsy, and assess if their surgical removal was associated with post-operative seizure freedom. Methods: We studied 31 patients with drug refractory focal epilepsy who underwent resting state magnetoencephalography (MEG), and structural magnetic resonance imaging (MRI) as part of pre-surgical evaluation. Using the structural MRI, we generated 114 cortical regions of interest, performed surface reconstruction and MEG source localization. Representative source localized signals for each region were correlated with each other to generate a functional brain network. We repeated this procedure across three randomly chosen one-minute epochs. Network hubs were defined as those with the highest intra-hemispheric mean correlations. Post-operative MRI identified regions that were surgically removed. Results: Greater mean MEG network connectivity was associated with a longer duration of epilepsy. Patients who were seizure free after surgery had more hubs surgically removed than patients who were not seizure free (AUC = 0.76, p = 0.01) consistently across three randomly chosen time segments. Conclusion: Our results support a growing literature implicating network hub involvement in focal epilepsy, the removal of which by surgery is associated with greater chance of post-operative seizure freedom.

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.