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

User-defined virtual sensors: A new solution to the problem of temporal plus epilepsy sources

OBJECTIVE

The most common medically resistant epilepsy (MRE) involves the temporal lobe (TLE), and children designated as temporal plus epilepsy (TLE+) have a five-times increased risk of postoperative surgical failure. This retrospective, blinded, cross-sectional study aimed to correlate visual and computational analyses of magnetoencephalography (MEG) virtual sensor waveforms with surgical outcome and epilepsy classification (TLE and TLE+).

METHODS

Patients with MRE who underwent MEG and iEEG monitoring and had at least 1 year of postsurgical follow-up were included in this retrospective analysis. User-defined virtual sensor (UDvs) beamforming was completed with virtual sensors placed manually and symmetrically in the bilateral amygdalohippocampi, inferior/middle/superior temporal gyri, insula, suprasylvian operculum, orbitofrontal cortex, and temporoparieto-occipital junction. Additionally, MEG effective connectivity was computed and quantified using eigenvector centrality (EC) to identify hub regions. More conventional MEG methods (equivalent current dipole [ECD], standardized low-resolution brain electromagnetic tomography, synthetic aperture magnetometry beamformer), UDvs beamformer, and EC hubs were compared to iEEG.

RESULTS

Eighty patients (38 female, 42 male) with MRE (mean age = 11.3 ± 6.2 years, range = 1.0-31.5) were identified and included. Twenty-five patients (31.3%) were classified as TLE, whereas 55 (68.8%) were TLE+. When modeling the association between MEG method, iEEG, and postoperative surgical outcome (odds of a worse [International League Against Epilepsy (ILAE) class > 2] outcome), a significant result was seen only for UDvs beamformer (odds ratio [OR] = 1.22, 95% confidence interval [CI] = 1.01-1.48). Likewise, when the relationship between MEG method, iEEG, and classification (TLE and TLE+) was modeled, only UDvs beamformer had a significant association (OR = 1.47, 95% CI = 1.13-1.92). When modeling the association between EC hub location and resection/ablation to postoperative surgical outcome (odds of a good [ILAE 1-2] outcome), a significant association was seen (OR = 1.22, 95% CI = 1.05-1.43).

SIGNIFICANCE

This study demonstrates a concordance between UDvs beamforming and iEEG that is related to both postsurgical seizure outcome and presurgical classification of epilepsy (TLE and TLE+). UDvs beamforming could be a complementary approach to the well-established ECD, improving invasive electrode and surgical resection planning for patients undergoing epilepsy surgery evaluations and treatments.

Genetic fingerprinting with heritable phenotypes of the resting-state brain network topology

Cognitive, behavioral, and disease traits are influenced by both genetic and environmental factors. Individual differences in these traits have been associated with graph theoretical properties of resting-state networks, indicating that variations in connectome topology may be driven by genetics. In this study, we establish the heritability of global and local graph properties of resting-state networks derived from functional MRI (fMRI) and magnetoencephalography (MEG) using a large sample of twins and non-twin siblings from the Human Connectome Project. We examine the heritability of MEG in the source space, providing a more accurate estimate of genetic influences on electrophysiological networks. Our findings show that most graph measures are more heritable for MEG compared to fMRI and the heritability for MEG is greater for amplitude compared to phase synchrony in the delta, high beta, and gamma frequency bands. This suggests that the fast neuronal dynamics in MEG offer unique insights into the genetic basis of brain network organization. Furthermore, we demonstrate that brain network features can serve as genetic fingerprints to accurately identify pairs of identical twins within a cohort. These results highlight novel opportunities to relate individual connectome signatures to genetic mechanisms underlying brain function.

Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

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.

Hub overload and failure as a final common pathway in neurological brain network disorders

Understanding the concept of network hubs and their role in brain disease is now rapidly becoming important for clinical neurology. Hub nodes in brain networks are areas highly connected to the rest of the brain, which handle a large part of all the network traffic. They also show high levels of neural activity and metabolism, which makes them vulnerable to many different types of pathology. The present review examines recent evidence for the prevalence and nature of hub involvement in a variety of neurological disorders, emphasizing common themes across different types of pathology. In focal epilepsy, pathological hubs may play a role in spreading of seizure activity, and removal of such hub nodes is associated with improved outcome. In stroke, damage to hubs is associated with impaired cognitive recovery. Breakdown of optimal brain network organization in multiple sclerosis is accompanied by cognitive dysfunction. In Alzheimer's disease, hyperactive hub nodes are directly associated with amyloid-beta and tau pathology. Early and reliable detection of hub pathology and disturbed connectivity in Alzheimer's disease with imaging and neurophysiological techniques opens up opportunities to detect patients with a network hyperexcitability profile, who could benefit from treatment with anti-epileptic drugs.

Towards the automated detection of interictal epileptiform discharges with magnetoencephalography

BACKGROUND

The analysis of clinical magnetoencephalography (MEG) in patients with epilepsy traditionally relies on visual identification of interictal epileptiform discharges (IEDs), which is time consuming and dependent on subjective criteria.

NEW METHOD

Here, we explore the ability of Independent Components Analysis (ICA) and Hidden Markov Modeling (HMM) to automatically detect and localize IEDs. We tested our pipelines on resting-state MEG recordings from 10 school-aged children with (multi)focal epilepsy.

RESULTS

In focal epilepsy patients, both pipelines successfully detected visually identified IEDs, but also revealed unidentified low-amplitude IEDs. Success was more mitigated in patients with multifocal epilepsy, as our automated pipeline missed IED activity associated with some foci-an issue that could be alleviated by post-hoc manual selection of epileptiform ICs or HMM states.

COMPARISON WITH EXISTING METHODS

We compared our results with visual IED detection by an experienced clinical magnetoencephalographer, getting heightened sensitivity and requiring minimal input from clinical practitioners.

CONCLUSIONS

IED detection based on ICA or HMM represents an efficient way to identify IED localization and timing. The development of these automatic IED detection algorithms provide a step forward in clinical MEG practice by decreasing the duration of MEG analysis and enhancing its sensitivity.

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.

Connectome-based approaches in pediatric epilepsy surgery: "State-of-the art" and future perspectives

Modern epilepsy science has overcome the traditional interpretation of a strict region-specific origin of epilepsy, highlighting the involvement of wider patterns of altered neuronal circuits. In selected cases, surgery may constitute a valuable option to achieve both seizure freedom and neurocognitive improvement. Although epilepsy is now considered as a brain network disease, the most relevant literature concerning the "connectome-based" epilepsy surgery mainly refers to adults, with a limited number of studies dedicated to the pediatric population. In this review, the Authors summarized the main current available knowledge on the relevance of WM surgical anatomy in epilepsy surgery, the post-surgical modifications of brain structural connectivity and the related clinical impact of such modifications within the pediatric context. In the last part, possible implications and future perspectives of this approach have been discussed, especially concerning the optimization of surgical strategies and the predictive value of the epilepsy network analysis for planning tailored approaches, with the final aim of improving case selection, presurgical planning, intraoperative management, and postoperative results.

Profiles of resting state functional connectivity in temporal lobe epilepsy associated with post-laser interstitial thermal therapy seizure outcomes and semiologies

Introduction

It is now understood that in focal epilepsy, impacted neural regions are not limited to the epileptogenic zone. As such, further investigation into the underlying functional connectivity (FC) patterns in those enduring Temporal Lobe Epilepsy (TLE) with Mesial Temporal Sclerosis (MTS) is imperative to understanding the intricacies of the disease.

Methods

The rsfMRIs of 17 healthy participants, 10 left-sided TLE-MTS patients with a pre-operative history of focal impaired awareness seizures (FIA), and 13 left-sided TLE-MTS patients with a pre-operative history of focal aware seizures (FA) were compared to determine the existence of distinct FC patterns with respect to seizure types. Similarly, the rsfMRIs of the above-mentioned healthy participants, 16 left-sided TLE-MTS individuals who were seizure-free (SF) 12 months postoperatively, and 16 left-sided TLE-MTS persons without seizure freedom (nSF) were interrogated. The ROI-to-ROI connectivity analysis included a total of 175 regions of interest (ROIs) and accounted for both age and duration of epileptic activity. Significant correlations were determined via two-sample t-tests and Bonferroni correction (α = 0.05).

Results

Comparisons of FA and FIA groups depicted significant correlations between the contralateral anterior cingulate gyrus, subgenual region, and the contralateral cerebellum, lobule III (p-value = 2.26e-4, mean z-score = -0.05 ± 0.28, T = -4.23). Comparisons of SF with nSF depicted two significantly paired-ROIs; the contralateral amygdala and the contralateral precuneus (p-value = 2.9e-5, mean z-score = -0.12 ± 0.19, T = 4.98), as well as the contralateral locus coeruleus and the ipsilateral intralaminar nucleus (p-value= 1.37e-4, mean z-score = 0.06 ± 0.17, T = -4.41).

Significance

FC analysis proves to be a lucrative modality for exploring unique signatures with respect to seizure types and postoperative outcomes. By furthering our understanding of the differences between epileptic phenotypes, we can achieve improvement in future treatment modalities not limited to targeting advancements.

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.

Identifying epileptogenic abnormalities through spatial clustering of MEG interictal band power

Successful epilepsy surgery depends on localizing and resecting cerebral abnormalities and networks that generate seizures. Abnormalities, however, may be widely distributed across multiple discontiguous areas. We propose spatially constrained clusters as candidate areas for further investigation and potential resection. We quantified the spatial overlap between the abnormality cluster and subsequent resection, hypothesizing a greater overlap in seizure-free patients. Thirty-four individuals with refractory focal epilepsy underwent pre-surgical resting-state interictal magnetoencephalography (MEG) recording. Fourteen individuals were totally seizure-free (ILAE 1) after surgery and 20 continued to have some seizures post-operatively (ILAE 2+). Band power abnormality maps were derived using controls as a baseline. Patient abnormalities were spatially clustered using the k-means algorithm. The tissue within the cluster containing the most abnormal region was compared with the resection volume using the dice score. The proposed abnormality cluster overlapped with the resection in 71% of ILAE 1 patients. Conversely, an overlap only occurred in 15% of ILAE 2+ patients. This effect discriminated outcome groups well (AUC = 0.82). Our novel approach identifies clusters of spatially similar tissue with high abnormality. This is clinically valuable, providing (a) a data-driven framework to validate current hypotheses of the epileptogenic zone localization or (b) to guide further investigation.

Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy

Normal brain functioning emerges from a complex interplay among regions forming networks. In epilepsy, these networks are disrupted causing seizures. Highly connected nodes in these networks are epilepsy surgery targets. Here, we assess whether functional connectivity (FC) using intracranial electroencephalography can quantify brain regions epileptogenicity and predict surgical outcome in children with drug resistant epilepsy (DRE). We computed FC between electrodes on different states (i.e. interictal without spikes, interictal with spikes, pre-ictal, ictal, and post-ictal) and frequency bands. We then estimated the electrodes' nodal strength. We compared nodal strength between states, inside and outside resection for good- (n = 22, Engel I) and poor-outcome (n = 9, Engel II-IV) patients, respectively, and tested their utility to predict the epileptogenic zone and outcome. We observed a hierarchical epileptogenic organization among states for nodal strength: lower FC during interictal and pre-ictal states followed by higher FC during ictal and post-ictal states (p < 0.05). We further observed higher FC inside resection (p < 0.05) for good-outcome patients on different states and bands, and no differences for poor-outcome patients. Resection of nodes with high FC was predictive of outcome (positive and negative predictive values: 47-100%). Our findings suggest that FC can discriminate epileptogenic states and predict outcome in patients with DRE.

An individual data-driven virtual resection model based on epileptic network dynamics in children with intractable epilepsy: a magnetoencephalography interictal activity application

Epilepsy surgery continues to be a recommended treatment for intractable (medication-resistant) epilepsy; however, 30-70% of epilepsy surgery patients can continue to have seizures. Surgical failures are often associated with incomplete resection or inaccurate localization of the epileptogenic zone. This retrospective study aims to improve surgical outcome through in silico testing of surgical hypotheses through a personalized computational neurosurgery model created from individualized patient's magnetoencephalography recording and MRI. The framework assesses the extent of the epileptic network and evaluates underlying spike dynamics, resulting in identification of one single brain volume as a candidate for resection. Dynamic-locked networks were utilized for virtual cortical resection. This in silico protocol was tested in a cohort of 24 paediatric patients with focal drug-resistant epilepsy who underwent epilepsy surgery. Of 24 patients who were included in the analysis, 79% (19 of 24) of the models agreed with the patient's clinical surgery outcome and 21% (5 of 24) were considered as model failures (accuracy 0.79, sensitivity 0.77, specificity 0.82). Patients with unsuccessful surgery outcome typically showed a model cluster outside of the resected cavity, while those with successful surgery showed the cluster model within the cavity. Two of the model failures showed the cluster in the vicinity of the resected tissue and either a functional disconnection or lack of precision of the magnetoencephalography-MRI overlapping could explain the results. Two other cases were seizure free for 1 year but developed late recurrence. This is the first study that provides in silico personalized protocol for epilepsy surgery planning using magnetoencephalography spike network analysis. This model could provide complementary information to the traditional pre-surgical assessment methods and increase the proportion of patients achieving seizure-free outcome from surgery.

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.

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.

MEG abnormalities and mechanisms of surgical failure in neocortical epilepsy

OBJECTIVE

Epilepsy surgery fails to achieve seizure freedom in 30%-40% of cases. It is not fully understood why some surgeries are unsuccessful. By comparing interictal magnetoencephalography (MEG) band power from patient data to normative maps, which describe healthy spatial and population variability, we identify patient-specific abnormalities relating to surgical failure. We propose three mechanisms contributing to poor surgical outcome: (1) not resecting the epileptogenic abnormalities (mislocalization), (2) failing to remove all epileptogenic abnormalities (partial resection), and (3) insufficiently impacting the overall cortical abnormality. Herein we develop markers of these mechanisms, validating them against patient outcomes.

METHODS

Resting-state MEG recordings were acquired for 70 healthy controls and 32 patients with refractory neocortical epilepsy. Relative band-power spatial maps were computed using source-localized recordings. Patient and region-specific band-power abnormalities were estimated as the maximum absolute z-score across five frequency bands using healthy data as a baseline. Resected regions were identified using postoperative magnetic resonance imaging (MRI). We hypothesized that our mechanistically interpretable markers would discriminate patients with and without postoperative seizure freedom.

RESULTS

Our markers discriminated surgical outcome groups (abnormalities not targeted: area under the curve [AUC] = 0.80, p = .003; partial resection of epileptogenic zone: AUC = 0.68, p = .053; and insufficient cortical abnormality impact: AUC = 0.64, p = .096). Furthermore, 95% of those patients who were not seizure-free had markers of surgical failure for at least one of the three proposed mechanisms. In contrast, of those patients without markers for any mechanism, 80% were ultimately seizure-free.

SIGNIFICANCE

The mapping of abnormalities across the brain is important for a wide range of neurological conditions. Here we have demonstrated that interictal MEG band-power mapping has merit for the localization of pathology and improving our mechanistic understanding of epilepsy. Our markers for mechanisms of surgical failure could be used in the future to construct predictive models of surgical outcome, aiding clinical teams during patient pre-surgical evaluations.

The association of magnetoencephalography high-frequency oscillations with epilepsy types and a ripple-based method with source-level connectivity for mapping epilepsy sources

OBJECTIVE

To explore the association between high-frequency oscillations (HFOs) and epilepsy types and to improve the accuracy of source localization.

METHODS

Magnetoencephalography (MEG) ripples of 63 drug-resistant epilepsy patients were detected. Ripple rates, distribution, spatial complexity, and the clustering coefficient of ripple channels were used for the preliminary classification of lateral temporal lobe epilepsy (LTLE), mesial temporal lobe epilepsy (MTLE), and nontemporal lobe epilepsy (NTLE), mainly frontal lobe epilepsy (FLE). Furthermore, the seizure site identification was improved using the Tucker LCMV method and source-level betweenness centrality.

RESULTS

Ripple rates were significantly higher in MTLE than in LTLE and NTLE (p < 0.05). The LTLE and MTLE were mainly distributed in the temporal lobe, followed by the parietal lobe, occipital lobe, and frontal lobe, whereas MTLE ripples were mainly distributed in the frontal lobe, then parietal lobe and occipital lobe. Nevertheless, the NTLE ripples were primarily in the frontal lobe and partially in the occipital lobe (p < 0.05). Meanwhile, the spatial complexity of NTLE was significantly higher than that of LTLE and MTLE and was lowest in MTLE (p < 0.01). However, an opposite trend was observed for the standardized clustering coefficient compared with spatial complexity (p < 0.01). Finally, the tucker algorithm showed a higher percentage of ripples at the surgical site when the betweenness centrality was added (p < 0.01).

CONCLUSION

This study demonstrated that HFO rates, distribution, spatial complexity, and clustering coefficient of ripple channels varied considerably among the three epilepsy types. Additionally, tucker MEG estimation combined with ripple rates based on the source-level functional connectivity is a promising approach for presurgical epilepsy evaluation.

Interictal Functional Connectivity in Focal Refractory Epilepsies Investigated by Intracranial EEG

Introduction: Focal epilepsies are diseases of neuronal excitability affecting macroscopic networks of cortical and subcortical neural structures. These networks ("epileptogenic networks") can generate pathological electrophysiological activities during seizures, and also between seizures (interictal period). Many works attempt to describe these networks by using quantification methods, particularly based on the estimation of statistical relationships between signals produced by brain regions, namely functional connectivity (FC). Results: FC has been shown to be greatly altered during seizures and in the immediate peri-ictal period. An increasing number of studies have shown that FC is also altered during the interictal period depending on the degree of epileptogenicity of the structures. Furthermore, connectivity values could be correlated with other clinical variables including surgical outcome. Significance: This leads to a conceptual change and to consider epileptic areas as both hyperexcitable and abnormally connected. These data open the door to the use of interictal FC as a marker of epileptogenicity and as a complementary tool for predicting the effect of surgery. Aim: In this article, we review the available data concerning interictal FC estimated from intracranial electroencephalograhy (EEG) in focal epilepsies and discuss it in the light of data obtained from other modalities (EEG imaging) and modeling studies.

Virtual MEG sensors based on beamformer and independent component analysis can reconstruct epileptic activity as measured on simultaneous intracerebral recordings

The prevailing gold standard for presurgical determination of epileptogenic brain networks is intracerebral EEG, a potent yet invasive approach. Magnetoencephalography (MEG) is a state-of-the art non-invasive method for investigating epileptiform discharges. However, it is not clear at what level the precision offered by MEG can reach that of SEEG. Here, we present a strategy for non-invasively retrieving the constituents of the interictal network, with high spatial and temporal precision. Our method is based on MEG and a combination of spatial filtering and independent component analysis (ICA). We validated this approach in twelve patients with drug-resistant focal epilepsy, thanks to the unprecedented ground truth provided by simultaneous recordings of MEG and SEEG. A minimum variance adaptive beamformer estimated the source time series and ICA was used to further decompose these time series into network constituents (MEG-ICs), each having a time series (virtual electrode) and a topography (spatial distribution of amplitudes in the brain). We show that MEG has a considerable sensitivity of 0.80 and 0.84 and a specificity of 0.93 and 0.91 for reconstructing deep and superficial sources, respectively, when compared to the ground truth (SEEG). For each epileptic MEG-IC (n = 131), we found at least one significantly correlating SEEG contact close to zero lag after correcting for multiple comparisons. All the patients except one had at least one epileptic component that was highly correlated (Spearman rho>0.3) with that of SEEG traces. MEG-ICs correlated well with SEEG traces. The strength of correlation coefficients did not depend on the depth of the SEEG contacts or the clinical outcome of the patient. A significant proportion of the MEG-ICs (n = 83/131) were localized in proximity with their maximally correlating SEEG, within a mean distance of 20±12.18mm. Our research is the first to validate the MEG-retrieved beamformer IC sources against SEEG-derived ground truth in a simultaneous MEG-SEEG framework. Observations from the present study suggest that non-invasive MEG source components may potentially provide additional information, comparable to SEEG in a number of instances.

Alterations in the default mode network in rolandic epilepsy with mild spike-wave index in non-rapid eye movement sleep

Purpose

Rolandic epilepsy (RE) is one of the most common epilepsy syndromes during childhood. The aim of this study was to investigate the alterations in the default mode network (DMN) of RE patients whose spike-wave index (SWI) was within the 50–85% range during non-rapid eye movement (NREM) during sleep, as well as to detect early neuroimaging markers.

Methods

Resting-state data was recorded for each subject using magnetoencephalography (MEG). DMN-related brain regions were chosen as regions of interest. The spectral power and functional connectivity (FC) of the DMN were estimated through the use of minimum norm estimation (MNE) combined with Welch technique and corrected amplitude envelope correlation (AEC-c).

Results

The patient group included 20 patients with NREM phase 50% ≤ SWI < 85% (mild SWI group), and 18 typical RE patients (SWI < 50% group). At the regional level, the mild SWI group exhibited enhanced spectral power in the delta band of the bilateral posterial cingulate cortex and attenuated the spectral power in the alpha band of the bilateral posterial cingulate cortex. Enhanced spectral power in the bilateral precuneus (PCu) in the delta band and attenuated spectral power in the right lateral temporal cortex (LTC) in the alpha band were common across all RE patients. At the FC level, patients in the mild SWI group indicated increased AEC-c values between the bilateral posterial cingulate cortex in the delta band and between the left medial frontal cortex (MFC) and bilateral posterial cingulate cortex in the alpha band. Increased AEC-c values between the right PCu and left MFC in the delta band, and between the left PCu and right MFC in the theta band, were common across all RE patients. Moreover, the spectral power in the bilateral posterial cingulate cortex in the alpha band and the AEC-c value between the bilateral posterial cingulate cortex in the delta band demonstrated good discrimination ability.

Conclusion

The spectral power of the bilateral posterior cingulate cortex (PCC) in the alpha band and the AEC-c value between the bilateral PCC in the delta band may be promising indicators of early differentiation between mild SWI and typical RE.

Virtual implantation using conventional scalp EEG delineates seizure onset and predicts surgical outcome in children with epilepsy

OBJECTIVE

Delineation of the seizure onset zone (SOZ) is required in children with drug resistant epilepsy (DRE) undergoing neurosurgery. Intracranial EEG (icEEG) serves as gold standard but has limitations. Here, we examine the utility of virtual implantation with electrical source imaging (ESI) on ictal scalp EEG for mapping the SOZ and predict surgical outcome.

METHODS

We retrospectively analyzed EEG data from 35 children with DRE who underwent surgery and dichotomized into seizure-free (SF) and non-seizure-free (NSF). We estimated virtual sensors (VSs) at brain locations that matched icEEG implantation and compared ictal patterns at VSs vs icEEG. We calculated the agreement between VSs SOZ and clinically defined SOZ and built receiver operating characteristic (ROC) curves to test whether it predicted outcome.

RESULTS

Twenty-one patients were SF after surgery. Moderate agreement between virtual and icEEG patterns was observed (kappa = 0.45, p < 0.001). Virtual SOZ agreement with clinically defined SOZ was higher in SF vs NSF patients (66.6% vs 41.6%, p = 0.01). Anatomical concordance of virtual SOZ with clinically defined SOZ predicted outcome (AUC = 0.73; 95% CI: 0.57-0.89; sensitivity = 66.7%; specificity = 78.6%; accuracy = 71.4%).

CONCLUSIONS

Virtual implantation on ictal scalp EEG can approximate the SOZ and predict outcome.

SIGNIFICANCE

SOZ mapping with VSs may contribute to tailoring icEEG implantation and predict outcome.

Epidemic models characterize seizure propagation and the effects of epilepsy surgery in individualized brain networks based on MEG and invasive EEG recordings

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients. However, seizure-freedom is currently achieved in only 2/3 of the patients after surgery. In this study we have developed an individualized computational model based on MEG brain networks to explore seizure propagation and the efficacy of different virtual resections. Eventually, the goal is to obtain individualized models to optimize resection strategy and outcome. We have modelled seizure propagation as an epidemic process using the susceptible-infected (SI) model on individual brain networks derived from presurgical MEG. We included 10 patients who had received epilepsy surgery and for whom the surgery outcome at least one year after surgery was known. The model parameters were tuned in in order to reproduce the patient-specific seizure propagation patterns as recorded with invasive EEG. We defined a personalized search algorithm that combined structural and dynamical information to find resections that maximally decreased seizure propagation for a given resection size. The optimal resection for each patient was defined as the smallest resection leading to at least a 90% reduction in seizure propagation. The individualized model reproduced the basic aspects of seizure propagation for 9 out of 10 patients when using the resection area as the origin of epidemic spreading, and for 10 out of 10 patients with an alternative definition of the seed region. We found that, for 7 patients, the optimal resection was smaller than the resection area, and for 4 patients we also found that a resection smaller than the resection area could lead to a 100% decrease in propagation. Moreover, for two cases these alternative resections included nodes outside the resection area. Epidemic spreading models fitted with patient specific data can capture the fundamental aspects of clinically observed seizure propagation, and can be used to test virtual resections in silico. Combined with optimization algorithms, smaller or alternative resection strategies, that are individually targeted for each patient, can be determined with the ultimate goal to improve surgery outcome. MEG-based networks can provide a good approximation of structural connectivity for computational models of seizure propagation, and facilitate their clinical use.

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.

Reproducibility of graph measures derived from resting-state MEG functional connectivity metrics in sensor and source spaces

Prior studies have used graph analysis of resting‐state magnetoencephalography (MEG) to characterize abnormal brain networks in neurological disorders. However, a present challenge for researchers is the lack of guidance on which network construction strategies to employ. The reproducibility of graph measures is important for their use as clinical biomarkers. Furthermore, global graph measures should ideally not depend on whether the analysis was performed in the sensor or source space. Therefore, MEG data of the 89 healthy subjects of the Human Connectome Project were used to investigate test–retest reliability and sensor versus source association of global graph measures. Atlas‐based beamforming was used for source reconstruction, and functional connectivity (FC) was estimated for both sensor and source signals in six frequency bands using the debiased weighted phase lag index (dwPLI), amplitude envelope correlation (AEC), and leakage‐corrected AEC. Reliability was examined over multiple network density levels achieved with proportional weight and orthogonal minimum spanning tree thresholding. At a 100% density, graph measures for most FC metrics and frequency bands had fair to excellent reliability and significant sensor versus source association. The greatest reliability and sensor versus source association was obtained when using amplitude metrics. Reliability was similar between sensor and source spaces when using amplitude metrics but greater for the source than the sensor space in higher frequency bands when using the dwPLI. These results suggest that graph measures are useful biomarkers, particularly for investigating functional networks based on amplitude synchrony.

Mapping Functional Connectivity of Epileptogenic Networks through Virtual Implantation

Children with medically refractory epilepsy (MRE) require resective neurosurgery to achieve seizure freedom, whose success depends on accurate delineation of the epileptogenic zone (EZ). Functional connectivity (FC) can assess the extent of epileptic brain networks since intracranial EEG (icEEG) studies have shown its link to the EZ and predictive value for surgical outcome in these patients. Here, we propose a new noninvasive method based on magnetoencephalography (MEG) and high-density (HD-EEG) data that estimates FC metrics at the source level through an "implantation" of virtual sensors (VSs). We analyzed MEG, HD-EEG, and icEEG data from eight children with MRE who underwent surgery having good outcome and performed source localization (beamformer) on noninvasive data to build VSs at the icEEG electrode locations. We analyzed data with and without Interictal Epileptiform Discharges (IEDs) in different frequency bands, and computed the following FC matrices: Amplitude Envelope Correlation (AEC), Correlation (CORR), and Phase Locking Value (PLV). Each matrix was used to generate a graph using Minimum Spanning Tree (MST), and for each node (i.e., each sensor) we computed four centrality measures: betweenness, closeness, degree, and eigenvector. We tested the reliability of VSs measures with respect to icEEG (regarded as benchmark) via linear correlation, and compared FC values inside vs. outside resection. We observed higher FC inside than outside resection (p<0.05) for AEC [alpha (8-12 Hz), beta (12-30 Hz), and broadband (1-50 Hz)] on data with IEDs and AEC theta (4-8 Hz) on data without IEDs for icEEG, AEC broadband (1-50 Hz) on data without IEDs for MEG-VSs, as well as for all centrality measures of icEEG and MEG/HD-EEG-VSs. Additionally, icEEG and VSs metrics presented high correlation (0.6-0.9, p<0.05). Our data support the notion that the proposed method can potentially replicate the icEEG ability to map the epileptogenic network in children with MRE.Clinical Relevance - The estimation of FC with noninvasive techniques, such as MEG and HD-EEG, via VSs is a promising tool that would help the presurgical evaluation by delineating the EZ without waiting for a seizure to occur, and potentially improve the surgical outcome of patients with MRE undergoing surgery.

Localizing the Epileptogenic Zone with Novel Biomarkers

Several noninvasive methods, such as high-density EEG or magnetoencephalography, are currently used to delineate the epileptogenic zone (EZ) during the presurgical evaluation of patients with drug resistant epilepsy (DRE). Yet, none of these methods can reliably identify the EZ by their own. In most cases a multimodal approach is needed. Challenging cases often require the implantation of intracranial electrodes, either through stereo-taxic EEG or electro-corticography. Recently, a growing body of literature introduces novel biomarkers of epilepsy that can be used for analyzing both invasive as well as noninvasive electrophysiological data. Some of these biomarkers are able to delineate the EZ with high precision, augment the presurgical evaluation, and predict the surgical outcome of patients with DRE undergoing surgery. However, the use of these epilepsy biomarkers in clinical practice is limited. Here, we summarize and discuss the latest technological advances in the presurgical neurophysiological evaluation of children with DRE with emphasis on electric and magnetic source imaging, high frequency oscillations, and functional connectivity.

Cellular Substrates of Functional Network Integration and Memory in Temporal Lobe Epilepsy

Temporal lobe epilepsy (TLE) patients are at risk of memory deficits, which have been linked to functional network disturbances, particularly of integration of the default mode network (DMN). However, the cellular substrates of functional network integration are unknown. We leverage a unique cross-scale dataset of drug-resistant TLE patients (n = 31), who underwent pseudo resting-state functional magnetic resonance imaging (fMRI), resting-state magnetoencephalography (MEG) and/or neuropsychological testing before neurosurgery. fMRI and MEG underwent atlas-based connectivity analyses. Functional network centrality of the lateral middle temporal gyrus, part of the DMN, was used as a measure of local network integration. Subsequently, non-pathological cortical tissue from this region was used for single cell morphological and electrophysiological patch-clamp analysis, assessing integration in terms of total dendritic length and action potential rise speed. As could be hypothesized, greater network centrality related to better memory performance. Moreover, greater network centrality correlated with more integrative properties at the cellular level across patients. We conclude that individual differences in cognitively relevant functional network integration of a DMN region are mirrored by differences in cellular integrative properties of this region in TLE patients. These findings connect previously separate scales of investigation, increasing translational insight into focal pathology and large-scale network disturbances in TLE.

Optimization of epilepsy surgery through virtual resections on individual structural brain networks

The success of epilepsy surgery in patients with refractory epilepsy depends upon correct identification of the epileptogenic zone (EZ) and an optimal choice of the resection area. In this study we developed individualized computational models based upon structural brain networks to explore the impact of different virtual resections on the propagation of seizures. The propagation of seizures was modelled as an epidemic process [susceptible-infected-recovered (SIR) model] on individual structural networks derived from presurgical diffusion tensor imaging in 19 patients. The candidate connections for the virtual resection were all connections from the clinically hypothesized EZ, from which the seizures were modelled to start, to other brain areas. As a computationally feasible surrogate for the SIR model, we also removed the connections that maximally reduced the eigenvector centrality (EC) (large values indicate network hubs) of the hypothesized EZ, with a large reduction meaning a large effect. The optimal combination of connections to be removed for a maximal effect were found using simulated annealing. For comparison, the same number of connections were removed randomly, or based on measures that quantify the importance of a node or connection within the network. We found that 90% of the effect (defined as reduction of EC of the hypothesized EZ) could already be obtained by removing substantially less than 90% of the connections. Thus, a smaller, optimized, virtual resection achieved almost the same effect as the actual surgery yet at a considerably smaller cost, sparing on average 27.49% (standard deviation: 4.65%) of the connections. Furthermore, the maximally effective connections linked the hypothesized EZ to hubs. Finally, the optimized resection was equally or more effective than removal based on structural network characteristics both regarding reducing the EC of the hypothesized EZ and seizure spreading. The approach of using reduced EC as a surrogate for simulating seizure propagation can suggest more restrictive resection strategies, whilst obtaining an almost optimal effect on reducing seizure propagation, by taking into account the unique topology of individual structural brain networks of patients.

Investigation of dynamic functional connectivity of the source reconstructed epileptiform discharges in focal epilepsy: A graph theory approach

OBJECTIVE

The aim of the present study is to investigate with noninvasive methods the modulation of dynamic functional connectivity during interictal epileptiform discharge (IED).

METHOD

We reconstructed the cortical source of the EEG recorded IED of 17 patients with focal epilepsy. We then computed dynamic connectivity using the time resolved phase locking value (PLV). We derived graph theory indices (i.e. degree, strength, local efficiency, clustering coefficient and global efficiency). Finally, we selected the atlas node with the maximum activation as the IED cortical source investigating the graph indices dynamics in theta, alpha, beta and gamma frequency bands.

RESULTS

We observed IED-locked modulations of the graph indexes depending on the frequency bands. We detected a modulation of the strength, clustering coefficient, local and global efficiency both in theta and in alpha bands, which also displayed modulations of the degree index. In the beta band only the global efficiency was modulated by the IED, while no effects were detected in the gamma band. Finally, we found a correlation between alpha and theta local efficiency, as well as alpha global efficiency, and the epilepsy duration.

SIGNIFICANCE

Our findings suggest that the neural synchronization is not limited to the IED cortical source, but implies a phase synchronization across multiple brain areas. We hypothesize that the aberrant electrical activity originating from the IED locus is spread amongst the other network nodes throughout the low frequency bands (i.e. theta and alpha). Moreover, IED-dependent increase in the global efficiency indicates that the IED interfere with the whole network functioning. We finally discussed possible application of this methodology for future investigation.

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.

Parieto-premotor functional connectivity changes during parietal lobe seizures are associated with motor semiology

OBJECTIVES

Parietal lobe seizures (PLS) are characterized by multiple clinical manifestations including motor signs. The mechanisms underlying the occurrence of motor signs are poorly understood. The main objective of this work was to estimate the functional coupling of brain regions associated with this clinical presentation.

METHODS

We retrospectively selected patients affected by drug-resistant epilepsy who underwent Stereoelectroencephalography (SEEG) for pre-surgical evaluation and in whom the seizure onset zone (SOZ) was located in the parietal cortex. The SOZ was defined visually and quantitatively by the epileptogenicity index (EI) method. Two groups of seizures were defined according to the presence ("motor seizures") or the absence ("non-motor seizures") of motor signs. Functional connectivity (FC) estimation was based on pairwise nonlinear regression analysis (h² coefficient). To study FC changes between parietal, frontal and temporal regions, for each patient, z-score values of 16 cortico-cortical interactions were obtained comparing h² coefficients of pre-ictal, seizure onset and seizure propagation periods.

RESULTS

We included 22 patients, 13 with "motor seizures" and 9 with "non-motor seizures". Resective surgery was performed in 14 patients, 8 patients had a positive surgical outcome (Engel's class I and II). During seizure onset period, a decrease of FC was observed and was significantly more important (in comparison with background period) in "motor" seizures. This was particularly observed between parietal operculum/post-central gyrus (OP/PoCg) and mesial temporal areas. During seizure propagation, a FC increase was significantly more important (in comparison with seizure onset) in "motor seizures", in particular between lateral pre-motor (pmL) area and precuneus, pmL and superior parietal lobule (SPL) and between inferior parietal lobule (IPL) and supplementary motor area (SMA).

CONCLUSIONS

Our study shows that motor semiology in PLS is accompanied by an increase of FC between parietal and premotor cortices, significantly different than what is observed in PLS without motor semiology.

SIGNIFICANCE

Our results indicate that preferential routes of coupling between parietal and premotor cortices are responsible for the prominent motor presentation during PLS.

Noninvasive Mapping of Ripple Onset Predicts Outcome in Epilepsy Surgery

OBJECTIVE

Intracranial electroencephalographic (icEEG) studies show that interictal ripples propagate across the brain of children with medically refractory epilepsy (MRE), and the onset of this propagation (ripple onset zone [ROZ]) estimates the epileptogenic zone. It is still unknown whether we can map this propagation noninvasively. The goal of this study is to map ripples (ripple zone [RZ]) and their propagation onset (ROZ) using high-density EEG (HD-EEG) and magnetoencephalography (MEG), and to estimate their prognostic value in pediatric epilepsy surgery.

METHODS

We retrospectively analyzed simultaneous HD-EEG and MEG data from 28 children with MRE who underwent icEEG and epilepsy surgery. Using electric and magnetic source imaging, we estimated virtual sensors (VSs) at brain locations that matched the icEEG implantation. We detected ripples on VSs, defined the virtual RZ and virtual ROZ, and estimated their distance from icEEG. We assessed the predictive value of resecting virtual RZ and virtual ROZ for postsurgical outcome. Interictal spike localization on HD-EEG and MEG was also performed and compared with ripples.

RESULTS

We mapped ripple propagation in all patients with HD-EEG and in 27 (96%) patients with MEG. The distance from icEEG did not differ between HD-EEG and MEG when mapping the RZ (26-27mm, p = 0.6) or ROZ (22-24mm, p = 0.4). Resecting the virtual ROZ, but not virtual RZ or the sources of spikes, was associated with good outcome for HD-EEG (p = 0.016) and MEG (p = 0.047).

INTERPRETATION

HD-EEG and MEG can map interictal ripples and their propagation onset (virtual ROZ). Noninvasively mapping the ripple onset may augment epilepsy surgery planning and improve surgical outcome of children with MRE. ANN NEUROL 2021;89:911-925.

Astroglial Glutamine Synthetase and the Pathogenesis of Mesial Temporal Lobe Epilepsy

The enzyme glutamine synthetase (GS), also referred to as glutamate ammonia ligase, is abundant in astrocytes and catalyzes the conversion of ammonia and glutamate to glutamine. Deficiency or dysfunction of astrocytic GS in discrete brain regions have been associated with several types of epilepsy, including medically-intractable mesial temporal lobe epilepsy (MTLE), neocortical epilepsies, and glioblastoma-associated epilepsy. Moreover, experimental inhibition or deletion of GS in the entorhinal-hippocampal territory of laboratory animals causes an MTLE-like syndrome characterized by spontaneous, recurrent hippocampal-onset seizures, loss of hippocampal neurons, and in some cases comorbid depressive-like features. The goal of this review is to summarize and discuss the possible roles of astroglial GS in the pathogenesis of epilepsy.

A review of magnetoencephalography use in pediatric epilepsy: an update on best practice

Introduction: Magnetoencephalography (MEG) is a noninvasive technique that is used for presurgical evaluation of children with drug-resistant epilepsy (DRE).Areas covered: The contributions of MEG for localizing the epileptogenic zone are discussed, in particular in extra-temporal lobe epilepsy and focal cortical dysplasia, which are common in children, as well as in difficult to localize epilepsy such as operculo-insular epilepsy. Further, the authors review current evidence on MEG for mapping eloquent cortex, its performance, application in clinical practice, and potential challenges.Expert opinion: MEG could change the clinical management of children with DRE by directing placement of intracranial electrodes thereby enhancing their yield. With improved identification of a circumscribed epileptogenic zone, MEG could render more patients as suitable candidates for epilepsy surgery and increase utilization of surgery.

Interictal and Ictal MEG in presurgical evaluation for epilepsy surgery

Although presurgical evaluation of patients with pharamacoresistent focal epilepsies provides essential information for successful epilepsy surgery, there is still a need for further improvement. Developments of noninvasive electrophysiological recording and analysis techniques offer additional information based on interictal and ictal epileptic activities. In this review, we provide an overview on the application of ictal magnetoencephalography (MEG). The results of a literature research for published interictal/ictal MEG findings and experiences with own cases are demonstrated and discussed. Ictal MEG may provide added value in comparison to interictal recordings. The results may be more focal and closer to the invasively determined seizure onset zone. In some patients without clear interictal findings, ictal MEG could provide correct localization. Novel recording and analysis techniques facilitate ictal recordings. However, extended recording durations, movement and artifacts still represent practical limitations. Ictal MEG may provide added value regarding the localization of the seizure onset zone but depends on the selection of patients and the application of optimal analysis techniques.

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.

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

OBJECTIVE

Focal epilepsy is a disorder affecting several brain networks; however, epilepsy surgery usually targets a restricted region, the so-called epileptic focus. There is a growing interest in embedding resting state (RS) connectivity analysis into pre-surgical workup.

APPROACH

In this retrospective study, we analyzed Magnetoencephalography (MEG) long-range RS functional connectivity patterns in patients with drug-resistant focal epilepsy. MEG recorded prior to surgery from seven seizure-free (Engel Ia) and five non seizure-free (Engel III or IV) patients were analyzed (minimum 2-years post-surgical follow-up). MEG segments without any detectable epileptic activity were source localized using wavelet-based Maximum Entropy on the Mean method. Amplitude envelope correlation in the theta (4-8 Hz), alpha (8-13 Hz), and beta (13-26 Hz) bands were used for assessing connectivity.

MAIN RESULTS

For seizure-free patients, we found an isolated epileptic network characterized by weaker connections between the brain region where interictal epileptic discharges (IED) are generated and the rest of the cortex, when compared to connectivity between the corresponding contralateral homologous region and the rest of the cortex. Contrarily, non seizure-free patients exhibited a widespread RS epileptic network characterized by stronger connectivity between the IED generator and the rest of the cortex, in comparison to the contralateral region and the cortex. Differences between the two seizure outcome groups concerned mainly distant long-range connections and were found in the alpha-band.

SIGNIFICANCE

Importantly, these connectivity patterns suggest specific mechanisms describing the underlying organization of the epileptic network and were detectable at the individual patient level, supporting the prospect use of MEG connectivity patterns in epilepsy to predict post-surgical seizure outcome.

Lateralization of epilepsy using intra-hemispheric brain networks based on resting-state MEG data

Focal epilepsy originates within networks in one hemisphere. However, previous studies have investigated network topologies for the entire brain. In this study, magnetoencephalography (MEG) was used to investigate functional intra‐hemispheric networks of healthy controls (HCs) and patients with left‐ or right‐hemispheric temporal lobe or temporal plus extra‐temporal lobe epilepsy. 22 HCs, 25 left patients (LPs), and 16 right patients (RPs) were enrolled. The debiased weighted phase lag index was used to calculate functional connectivity between 246 brain regions in six frequency bands. Global efficiency, characteristic path length, and transitivity were computed for left and right intra‐hemispheric networks. The right global graph measures (GGMs) in the theta band were significantly different (p < .005) between RPs and both LPs and HCs. Right and left GGMs in higher frequency bands were significantly different (p < .05) between HCs and the patients. Right GGMs were used as input features of a Naïve‐Bayes classifier to classify LPs and RPs (78.0% accuracy) and all three groups (75.5% accuracy). The complete theta band brain networks were compared between LPs and RPs with network‐based statistics (NBS) and with the clustering coefficient (CC), nodal efficiency (NE), betweenness centrality (BC), and eigenvector centrality (EVC). NBS identified a subnetwork primarily composed of right intra‐hemispheric connections. Significantly different (p < .05) nodes were primarily in the right hemisphere for the CC and NE and primarily in the left hemisphere for the BC and EVC. These results indicate that intra‐hemispheric MEG networks may be incorporated in the diagnosis and lateralization of focal epilepsy.

Networks in Temporal Lobe Epilepsy

Temporal lobe epilepsy (TLE) is the most common type of drug-resistant focal epilepsy. Epilepsy can be conceptualized as a network disorder with the epileptogenic zone a critical node of the network. Temporal lobe networks can be identified on the microscale and macroscale, both during the interictal and ictal periods. This review summarizes the current understanding of TLE networks as studied by neurophysiological and imaging techniques discussing both functional and structural connectivity.

Astrocytes and Glutamine Synthetase in Epileptogenesis

The cellular, molecular, and metabolic mechanisms that underlie the development of mesial temporal lobe epilepsy are incompletely understood. Here we review the role of astrocytes in epilepsy development (a.k.a. epileptogenesis), particularly astrocyte pathologies related to: aquaporin 4, the inwardly rectifying potassium channel Kir4.1, monocarboxylate transporters MCT1 and MCT2, excitatory amino acid transporters EAAT1 and EAAT2, and glutamine synthetase. We propose that inhibition, dysfunction or loss of astrocytic glutamine synthetase is an important causative factor for some epilepsies, particularly mesial temporal lobe epilepsy and glioblastoma-associated epilepsy. We postulate that the regulatory mechanisms of glutamine synthetase as well as the downstream effects of glutamine synthetase dysfunction, represent attractive, new targets for antiepileptogenic interventions. Currently, no antiepileptogenic therapies are available for human use. The discovery of such interventions is important as it will fundamentally change the way we approach epilepsy by preventing the disease from ever becoming manifest after an epileptogenic insult to the brain.

Seizure prediction - ready for a new era

Epilepsy is a common disorder characterized by recurrent seizures. An overwhelming majority of people with epilepsy regard the unpredictability of seizures as a major issue. More than 30 years of international effort have been devoted to the prediction of seizures, aiming to remove the burden of unpredictability and to couple novel, time-specific treatment to seizure prediction technology. A highly influential review published in 2007 concluded that insufficient evidence indicated that seizures could be predicted. Since then, several advances have been made, including successful prospective seizure prediction using intracranial EEG in a small number of people in a trial of a real-time seizure prediction device. In this Review, we examine advances in the field, including EEG databases, seizure prediction competitions, the prospective trial mentioned and advances in our understanding of the mechanisms of seizures. We argue that these advances, together with statistical evaluations, set the stage for a resurgence in efforts towards the development of seizure prediction methodologies. We propose new avenues of investigation involving a synergy between mechanisms, models, data, devices and algorithms and refine the existing guidelines for the development of seizure prediction technology to instigate development of a solution that removes the burden of the unpredictability of seizures.

Background EEG Connectivity Captures the Time-Course of Epileptogenesis in a Mouse Model of Epilepsy

Large-scale brain networks are increasingly recognized as important for the generation of seizures in epilepsy. However, how a network evolves from a healthy state through the process of epileptogenesis remains unclear. To address this question, here, we study longitudinal epicranial background EEG recordings (30 electrodes, EEG free from epileptiform activity) of a mouse model of mesial temporal lobe epilepsy. We analyze functional connectivity networks and observe that over the time course of epileptogenesis the networks become increasingly asymmetric. Furthermore, computational modelling reveals that a set of nodes, located outside of the region of initial insult, emerges as particularly important for the network dynamics. These findings are consistent with experimental observations, thus demonstrating that ictogenic mechanisms can be revealed on the EEG, that computational models can be used to monitor unfolding epileptogenesis and that both the primary focus and epileptic network play a role in epileptogenesis.

Electrophysiological Brain Connectivity: Theory and Implementation

We review the theory and algorithms of electrophysiological brain connectivity analysis. This tutorial is aimed at providing an introduction to brain functional connectivity from electrophysiological signals, including electroencephalography (EEG), magnetoencephalography (MEG), electrocorticography (ECoG), stereoelectroencephalography (SEEG). Various connectivity estimators are discussed, and algorithms introduced. Important issues for estimating and mapping brain functional connectivity with electrophysiology are discussed.

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.

A comparison between scalp- and source-reconstructed EEG networks

EEG can be used to characterise functional networks using a variety of connectivity (FC) metrics. Unlike EEG source reconstruction, scalp analysis does not allow to make inferences about interacting regions, yet this latter approach has not been abandoned. Although the two approaches use different assumptions, conclusions drawn regarding the topology of the underlying networks should, ideally, not depend on the approach. The aim of the present work was to find an answer to the following questions: does scalp analysis provide a correct estimate of the network topology? how big are the distortions when using various pipelines in different experimental conditions? EEG recordings were analysed with amplitude- and phase-based metrics, founding a strong correlation for the global connectivity between scalp- and source-level. In contrast, network topology was only weakly correlated. The strongest correlations were obtained for MST leaf fraction, but only for FC metrics that limit the effects of volume conduction/signal leakage. These findings suggest that these effects alter the estimated EEG network organization, limiting the interpretation of results of scalp analysis. Finally, this study also suggests that the use of metrics that address the problem of zero lag correlations may give more reliable estimates of the underlying network topology.

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.