Interictal networks in magnetoencephalography

Beyond Stereo-EEG: Is It Worth Combining Stereo-EEG With Other Diagnostic Methods?

SUMMARY

Stereo-EEG is a widely used method to improve the diagnostic precision of presurgical workup in patients with refractory epilepsy. Its ability to detect epileptic activity and identify epileptic networks largely depends on the chosen implantation strategy. Even in an ideal situation, electrodes record activity generated in <10% of the brain and contacts only record from brain tissue in their immediate proximity. In this article, the authors discuss how recording stereo-EEG simultaneously with other diagnostic methods can improve its diagnostic value in clinical and research settings. It can help overcome the limited spatial coverage of intracranial recording and better understand the sources of epileptic activity. Simultaneous scalp EEG is the most widely available method, often used to understand large epileptic networks, seizure propagation, and EEG activity occurring on the contralateral hemisphere. Simultaneous magnetoencephalography allows for more precise source localization and identification of deep sources outside the stereo-EEG coverage. Finally, simultaneous functional MRI can highlight metabolic changes following epileptic activity and help understand the widespread network changes associated with interictal activity. This overview highlights advantages and methodological challenges for all these methods. Clinical use and research applications are presented for each approach.

Pathological and Physiological High-frequency Oscillations on Electroencephalography in Patients with Epilepsy

High-frequency oscillations (HFOs) encompass ripples (80 Hz-200 Hz) and fast ripples (200 Hz-600 Hz), serving as a promising biomarker for localizing the epileptogenic zone in epilepsy. Spontaneous fast ripples are always pathological, while ripples may be physiological or pathological. Distinguishing physiological from pathological ripples is important not only for designating epileptogenic brain regions, but also for investigations that study ripples in the context of memory encoding, consolidation, and recall in patients with epilepsy. Many studies have sought to identify distinguishing features between pathological and physiological ripples over the past two decades. Physiological and pathological ripples differ with respect to their spatial location, cellular mechanisms, morphology, and coupling with background electroencephalographic activity. Retrospective studies have demonstrated that differentiating between pathological and physiological ripples can improve surgical outcome prediction. In this review, we summarize the characteristics, differences, and applications of pathological and physiological HFOs and discuss strategies for their clinical translation.

Identification of Early Hippocampal Dynamics during Recognition Memory with Independent Component Analysis

The hippocampus is generally considered to have relatively late involvement in recognition memory, its main electrophysiological signature being between 400 and 800 ms after stimulus onset. However, most electrophysiological studies have analyzed the hippocampus as a single responsive area, selecting only a single-site signal exhibiting the strongest effect in terms of amplitude. These classical approaches may not capture all the dynamics of this structure, hindering the contribution of other hippocampal sources that are not located in the vicinity of the selected site. We combined intracerebral electroencephalogram recordings from epileptic patients with independent component analysis during a recognition memory task involving the recognition of old and new images. We identified two sources with different responses emerging from the hippocampus: a fast one (maximal amplitude at ∼250 ms) that could not be directly identified from raw recordings and a latter one, peaking at ∼400 ms. The former component presented different amplitudes between old and new items in 6 out of 10 patients. The latter component had different delays for each condition, with a faster activation (∼290 ms after stimulus onset) for recognized items. We hypothesize that both sources represent two steps of hippocampal recognition memory, the faster reflecting the input from other structures and the latter the hippocampal internal processing. Recognized images evoking early activations would facilitate neural computation in the hippocampus, accelerating memory retrieval of complementary information. Overall, our results suggest that the hippocampal activity is composed of several sources with an early activation related to recognition memory.

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.

From cryogenic to on-scalp magnetoencephalography for the evaluation of paediatric epilepsy

Magnetoencephalography (MEG) is a neurophysiological technique based on the detection of brain magnetic fields. Whole-head MEG systems typically house a few hundred sensors requiring cryogenic cooling in a rigid one-size-fits-all (commonly adult-sized) helmet to keep a thermal insulation space. This leads to an increased brain-to-sensor distance in children, because of their smaller head circumference, and decreased signal-to-noise ratio. MEG allows detection and localization of interictal and ictal epileptiform discharges, and pathological high frequency oscillations, as a part of the presurgical assessment of children with refractory focal epilepsy, where electroencephalography is not contributive. MEG can also map the eloquent cortex before surgical resection. MEG also provides insights into the physiopathology of both generalized and focal epilepsy. On-scalp recordings based on cryogenic-free sensors have demonstrated their use in the field of childhood focal epilepsy and should become a reference technique for diagnosing epilepsy in the paediatric population. WHAT THIS PAPER ADDS: Magnetoencephalography (MEG) contributes to the diagnosis and understanding of paediatric epilepsy. On-scalp MEG recordings demonstrate some advantages over cryogenic MEG.

Combining independent component analysis and source localization for improving spatial sampling of stereoelectroencephalography in epilepsy

Stereoelectroencephalography is a powerful intracerebral EEG recording method for the presurgical evaluation of epilepsy. It consists in implanting depth electrodes in the patient's brain to record electrical activity and map the epileptogenic zone, which should be resected to render the patient seizure-free. Stereoelectroencephalography has high spatial accuracy and signal-to-noise ratio but remains limited in the coverage of the explored brain regions. Thus, the implantation might provide a suboptimal sampling of epileptogenic regions. We investigate the potential of improving a suboptimal stereoelectroencephalography recording by performing source localization on stereoelectroencephalography signals. We propose combining independent component analysis, connectivity measures to identify components of interest, and distributed source modelling. This approach was tested on two patients with two implantations each, the first failing to characterize the epileptogenic zone and the second giving a better diagnosis. We demonstrate that ictal and interictal source localization performed on the first stereoelectroencephalography recordings matches the findings of the second stereo-EEG exploration. Our findings suggest that independent component analysis followed by source localization on the topographies of interest is a promising method for retrieving the epileptogenic zone in case of suboptimal implantation.

Widespread slow oscillations support interictal epileptiform discharge networks in focal epilepsy

Interictal epileptiform discharges (IEDs) often co-occur across spatially-separated cortical regions, forming IED networks. However, the factors prompting IED propagation remain unelucidated. We hypothesized that slow oscillations (SOs) might facilitate IED propagation. Here, the amplitude and phase synchronization of SOs preceding propagating and non-propagating IEDs were compared in 22 patients with focal epilepsy undergoing intracranial electroencephalography (EEG) evaluation. Intracranial channels were categorized into the irritative zone (IZ) and normal zone (NOZ) regarding the presence of IEDs. During wakefulness, we found that pre-IED SOs within the IZ exhibited higher amplitudes for propagating IEDs than non-propagating IEDs (delta band: p = 0.001, theta band: p < 0.001). This increase in SOs was also concurrently observed in the NOZ (delta band: p = 0.04). Similarly, the inter-channel phase synchronization of SOs prior to propagating IEDs was higher than those preceding non-propagating IEDs in the IZ (delta band: p = 0.04). Through sliding window analysis, we observed that SOs preceding propagating IEDs progressively increased in amplitude and phase synchronization, while those preceding non-propagating IEDs remained relatively stable. Significant differences in amplitude occurred approximately 1150 ms before IEDs. During non-rapid eye movement (NREM) sleep, SOs on scalp recordings also showed higher amplitudes before intracranial propagating IEDs than before non-propagating IEDs (delta band: p = 0.006). Furthermore, the analysis of IED density around sleep SOs revealed that only high-amplitude sleep SOs demonstrated correlation with IED propagation. Overall, our study highlights that transient but widely distributed SOs are associated with IED propagation as well as generation in focal epilepsy during sleep and wakefulness, providing new insight into the EEG substrate supporting IED networks.

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.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

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

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

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.

Interictal discharges in the human brain are travelling waves arising from an epileptogenic source

While seizure activity may be electrographically widespread, increasing evidence has suggested that ictal discharges may in fact represent travelling waves propagated from a focal seizure source. Interictal epileptiform discharges (IEDs) are an electrographic manifestation of excessive hypersynchronization of cortical activity that occur between seizures and are considered a marker of potentially epileptogenic tissue. The precise relationship between brain regions demonstrating IEDs and those involved in seizure onset, however, remains poorly understood. Here, we hypothesize that IEDs likewise reflect the receipt of travelling waves propagated from the same regions which give rise to seizures. Forty patients from our institution who underwent invasive monitoring for epilepsy, proceeded to surgery and had at least one year of follow-up were included in our study. Interictal epileptiform discharges were detected using custom software, validated by a clinical epileptologist. We show that IEDs reach electrodes in sequences with a consistent temporal ordering, and this ordering matches the timing of receipt of ictal discharges, suggesting that both types of discharges spread as travelling waves. We use a novel approach for localization of ictal discharges, in which time differences of discharge receipt at nearby electrodes are used to compute source location; similar algorithms have been used in acoustics and geophysics. We find that interictal discharges co-localize with ictal discharges. Moreover, interictal discharges tend to localize to the resection territory in patients with good surgical outcome and outside of the resection territory in patients with poor outcome. The seizure source may originate at, and also travel to, spatially distinct IED foci. Our data provide evidence that interictal discharges may represent travelling waves of pathological activity that are similar to their ictal counterparts, and that both ictal and interictal discharges emerge from common epileptogenic brain regions. Our findings have important clinical implications, as they suggest that seizure source localizations may be derived from interictal discharges, which are much more frequent than seizures.

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.

Reconstruction and localization of auditory sources from intracerebral SEEG using independent component analysis

Stereo-electroencephalography (SEEG) is the surgical implantation of electrodes in the brain to better localize the epileptic network in pharmaco-resistant epileptic patients. This technique has exquisite spatial and temporal resolution. Still, the number and the position of the electrodes in the brain is limited and determined by the semiology and/or preliminary non-invasive examinations, leading to a large number of unexplored brain structures in each patient. Here, we propose a new approach to reconstruct the activity of non-sampled structures in SEEG, based on independent component analysis (ICA) and dipole source localization. We have tested this approach with an auditory stimulation dataset in ten patients. The activity directly recorded from the auditory cortex served as ground truth and was compared to the ICA applied on all non-auditory electrodes. Our results show that the activity from the auditory cortex can be reconstructed at the single trial level from contacts as far as ∼40 mm from the source. Importantly, this reconstructed activity is localized via dipole fitting in the proximity of the original source. In addition, we show that the size of the confidence interval of the dipole fitting is a good indicator of the reliability of the result, which depends on the geometry of the SEEG implantation. Overall, our approach allows reconstructing the activity of structures far from the electrode locations, partially overcoming the spatial sampling limitation of intracerebral recordings.

Novel noninvasive identification of patient-specific epileptic networks in focal epilepsies: Linking single-photon emission computed tomography perfusion during seizures with resting-state magnetoencephalography dynamics

Single-photon emission computed tomography (SPECT) during seizures and magnetoencephalography (MEG) during the interictal state are noninvasive modalities employed in the localization of the epileptogenic zone in patients with drug-resistant focal epilepsy (DRFE). The present study aims to investigate whether there exists a preferentially high MEG functional connectivity (FC) among those regions of the brain that exhibit hyperperfusion or hypoperfusion during seizures. We studied MEG and SPECT data in 30 consecutive DRFE patients who had resective epilepsy surgery. We parcellated each ictal perfusion map into 200 regions of interest (ROIs) and generated ROI time series using source modeling of MEG data. FC between ROIs was quantified using coherence and phase-locking value. We defined a generalized linear model to relate the connectivity of each ROI, ictal perfusion z score, and distance between ROIs. We compared the coefficients relating perfusion z score to FC of each ROI and estimated the connectivity within and between resected and unresected ROIs. We found that perfusion z scores were strongly correlated with the FC of hyper-, and separately, hypoperfused ROIs across patients. High interictal connectivity was observed between hyperperfused brain regions inside and outside the resected area. High connectivity was also observed between regions of ictal hypoperfusion. Importantly, the ictally hypoperfused regions had a low interictal connectivity to regions that became hyperperfused during seizures. We conclude that brain regions exhibiting hyperperfusion during seizures highlight a preferentially connected interictal network, whereas regions of ictal hypoperfusion highlight a separate, discrete and interconnected, interictal network.

Integration of white matter architecture to stereo-EEG better describes epileptic spike propagation

OBJECTIVE

Stereo-electroencephalography (SEEG)-derived epilepsy networks are used to better understand a patient's epilepsy; however, a unimodal approach provides an incomplete picture. We combine tractography and SEEG to determine the relationship between spike propagation and the white matter architecture and to improve our understanding of spike propagation mechanisms.

METHODS

Probablistic tractography from diffusion imaging (dMRI) of matched subjects from the Human Connectome Project (HCP) was combined with patient-specific SEEG-derived spike propagation networks. Two regions-of-interest (ROIs) with a significant spike propagation relationship constituted a Propagation Pair.

RESULTS

In 56 of 59 patients, Propagation Pairs were more often tract-connected as compared to all ROI pairs (p < 0.01; d = -1.91). The degree of spike propagation between tract-connected ROIs was greater (39 ± 21%) compared to tract-unconnected ROIs (31 ± 18%; p < 0.0001). Within the same network, ROIs receiving propagation earlier were more often tract-connected to the source (59.7%) as compared to late receivers (25.4%; p < 0.0001).

CONCLUSIONS

Brain regions involved in spike propagation are more likely to be connected by white matter tracts. Between nodes, presence of tracts suggests a direct course of propagation, whereas the absence of tracts suggests an indirect course of propagation.

SIGNIFICANCE

We demonstrate a logical and consistent relationship between spike propagation and the white matter architecture.

Comparison of beamformer and ICA for dynamic connectivity analysis: A simultaneous MEG-SEEG study

Magnetoencephalography (MEG) is a powerful tool for estimating brain connectivity with both good spatial and temporal resolution. It is particularly helpful in epilepsy to characterize non-invasively the epileptic networks. However, using MEG to map brain networks requires solving a difficult inverse problem that introduces uncertainty in the activity localization and connectivity measures. Our goal here was to compare independent component analysis (ICA) followed by dipole source localization and the linearly constrained minimum-variance beamformer (LCMV-BF) for characterizing regions with interictal epileptic activity and their dynamic connectivity. After a simulation study, we compared ICA and LCMV-BF results with intracerebral EEG (stereotaxic EEG, SEEG) recorded simultaneously in 8 epileptic patients, which provide a unique 'ground truth' to which non-invasive results can be confronted. We compared the signal time courses extracted applying ICA and LCMV-BF on MEG data to that of SEEG, both for the actual signals and the dynamic connectivity computed using cross-correlation (evolution of links in time). With our simulations, we illustrated the different effect of the temporal and spatial correlation among sources on the two methods. While ICA was more affected by the temporal correlation but robust against spatial configurations, LCMV-BF showed opposite behavior. Moreover, ICA seems more suited to retrieve the simulated networks. In case of real patient data, good MEG/SEEG correlation and good localization were obtained in 6 out of 8 patients. In 4 of them ICA had the best performance (higher correlation, lower localization distance). In terms of dynamic connectivity, the evolution in time of the cross-correlation links could be retrieved in 5 patients out of 6, however, with more variable results in terms of correlation and distance. In two patients LCMV-BF had better results than ICA. In one patient the two methods showed equally good outcomes, and in the remaining two patients ICA performed best. In conclusion, our results obtained by exploiting simultaneous MEG/SEEG recordings suggest that ICA and LCMV-BF have complementary qualities for retrieving the dynamics of interictal sources and their network interactions.

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.

Multilevel Feature Learning Method for Accurate Interictal Epileptiform Spike Detection

Interictal epileptiform spike (referred to as spike) detected from electroencephalograms lasting only 20- to 200-ms can provide a reliable evidence-based indicator for clinical seizure type diagnosis. Recent feature representation approaches focus either on the concrete-level or on abstract-level information mining of the spike, thus demonstrating suboptimal detection performance. Additionally, existing abstract-level information mining methods of the spike based deep learning networks have not realized the effective feature representation of long-term dependent distinguished information within similar waveform cycles caused by morphological heterogeneity, which affects detection performance. Thus, a multilevel feature learning method for accurate spike detection was proposed in this study. Specifically, the spatio-temporal-frequency multidomain information in concrete-level first are inferred the common mimetic properties of the spike using the multidomain feature extractors. Then, the effective feature representation of long-term dependent distinguished information within similar waveform cycles caused by morphological heterogeneity is suitably captured using the temporal convolutional network. Finally, the spatio-temporal-frequency multidomain long-term dependent feature representation of spike is calculated using the element-wise manner to fuse the feature representation in concrete- and abstract-levels. The experimental results indicate that the proposed method can achieve an accuracy of 90.62±1.38%, sensitivity of 90.38±1.52%, specificity of 91.00±1.60%, precision of 90.33±4.71%, and the false detection rate per minute is 0.148±0.020m^-1, which are higher than when using the feature representation in the concrete- or abstract-level alone. Additionally, the detection results indicate that the proposed method avoids the subjectivity and inefficiency of visual inspection, and it enables a highly accurate detection of the spike.

Magnetoencephalography can reveal deep brain network activities linked to memory processes

Recording from deep neural structures such as hippocampus noninvasively and yet with high temporal resolution remains a major challenge for human neuroscience. Although it has been proposed that deep neuronal activity might be recordable during cognitive tasks using magnetoencephalography (MEG), this remains to be demonstrated as the contribution of deep structures to MEG recordings may be too small to be detected or might be eclipsed by the activity of large‐scale neocortical networks. In the present study, we disentangled mesial activity and large‐scale networks from the MEG signals thanks to blind source separation (BSS). We then validated the MEG BSS components using intracerebral EEG signals recorded simultaneously in patients during their presurgical evaluation of epilepsy. In the MEG signals obtained during a memory task involving the recognition of old and new images, we identified with BSS a putative mesial component, which was present in all patients and all control subjects. The time course of the component selectively correlated with stereo‐electroencephalography signals recorded from hippocampus and rhinal cortex, thus confirming its mesial origin. This finding complements previous studies with epileptic activity and opens new possibilities for using MEG to study deep brain structures in cognition and in brain disorders.

Assessment of Effective Network Connectivity among MEG None Contaminated Epileptic Transitory Events

Characterizing epileptogenic zones EZ (sources responsible of excessive discharges) would assist a neurologist during epilepsy diagnosis. Locating efficiently these abnormal sources among magnetoencephalography (MEG) biomarker is obtained by several inverse problem techniques. These techniques present different assumptions and particular epileptic network connectivity. Here, we proposed to evaluate performances of distributed inverse problem in defining EZ. First, we applied an advanced technique based on Singular Value Decomposition (SVD) to recover only pure transitory activities (interictal epileptiform discharges). We evaluated our technique's robustness in separation between transitory and ripples versus frequency range, transitory shapes, and signal to noise ratio on simulated data (depicting both epileptic biomarkers and respecting time series and spectral properties of realistic data). We validated our technique on MEG signal using detector precision on 5 patients. Then, we applied four methods of inverse problem to define cortical areas and neural generators of excessive discharges. We computed network connectivity of each technique. Then, we confronted obtained noninvasive networks to intracerebral EEG transitory network connectivity using nodes in common, connection strength, distance metrics between concordant nodes of MEG and IEEG, and average propagation delay. Coherent Maximum Entropy on the Mean (cMEM) proved a high matching between MEG network connectivity and IEEG based on distance between active sources, followed by Exact low-resolution brain electromagnetic tomography (eLORETA), Dynamical Statistical Parametric Mapping (dSPM), and Minimum norm estimation (MNE). Clinical performance was interesting for entire methods providing in an average of 73.5% of active sources detected in depth and seen in MEG, and vice versa, about 77.15% of active sources were detected from MEG and seen in IEEG. Investigated problem techniques succeed at least in finding one part of seizure onset zone. dSPM and eLORETA depict the highest connection strength among all techniques. Propagation delay varies in this range [18, 25]ms, knowing that eLORETA ensures the lowest propagation delay (18 ms) and the closet one to IEEG propagation delay.

Individual localization value of resting-state fMRI in epilepsy presurgical evaluation: A combined study with stereo-EEG

OBJECTIVE

To examine the individual-patient-level localization value of resting-state functional MRI (rsfMRI) metrics for the seizure onset zone (SOZ) defined by stereo-electroencephalography (SEEG) in patients with medically intractable focal epilepsies.

METHODS

We retrospectively included 19 patients who underwent SEEG implantation for epilepsy presurgical evaluation. Voxel-wise whole-brain analysis was performed on 3.0 T rsfMRI to generate clusters for amplitude of low-frequency fluctuations (ALFF), regional homogeneity (ReHo) and degree centrality (DC), which were co-registered with the SEEG-defined SOZ to evaluate their spatial overlap. Subgroup and correlation analyses were conducted for various clinical characteristics.

RESULTS

ALFF demonstrated concordant clusters with SEEG-defined SOZ in 73.7% of patients, with 93.3% sensitivity and 77.8% PPV. The concordance rate showed no significant difference when subgrouped by lesional/non-lesional MRI, SOZ location, interictal epileptiform discharges on scalp EEG, pathology or seizure outcomes. No significant correlation was seen between ALFF concordance rate and epilepsy duration, seizure-onset age, seizure frequency or number of antiseizure medications. ReHo and DC did not achieve favorable concordance results (10.5% and 15.8%, respectively). All concordant clusters showed regional activation, representing increased neural activities.

CONCLUSION

ALFF had high concordance rate with SEEG-defined SOZ at individual-patient level.

SIGNIFICANCE

ALFF activation on rsfMRI can add localizing information for the noninvasive presurgical workup of intractable focal epilepsies.

Functional Interactions between Entorhinal Cortical Pathways Modulate Theta Activity in the Hippocampus

Simple Summary

The activity in the hippocampus is characterized by a strong oscillation at theta frequency that organizes the neuronal firing. We have recently shown that different theta oscillations are present in the hippocampus, opening the possibility to multiple interactions between theta rhythms. In this work, we analyzed the functional connectivity between theta generators during the exploration of a known environment with or without a novel stimulus. The directionality of the interactions was determined using tools based on Granger causality and transfer entropy. We found significant interactions between activity components originated in CA3 and in layers II and III of the entorhinal cortex. During exploration with a novel stimulus, the connectivity from the entorhinal cortex layer II increased, while the influence of CA3 decreased. These results suggest that the entorhinal cortex layer II may drive theta interactions and synchronization in the hippocampus during novelty exploration.

Abstract

Theta oscillations organize neuronal firing in the hippocampus during context exploration and memory formation. Recently, we have shown that multiple theta rhythms coexist in the hippocampus, reflecting the activity in their afferent regions in CA3 (Schaffer collaterals) and the entorhinal cortex layers II (EC-II, perforant pathway) and III (EC-III, temporoammonic pathway). Frequency and phase coupling between theta rhythms were modulated by the behavioral state, with synchronized theta rhythmicity preferentially occurring in tasks involving memory updating. However, information transmission between theta generators was not investigated. Here, we used source separation techniques to disentangle the current generators recorded in the hippocampus of rats exploring a known environment with or without a novel stimulus. We applied analytical tools based on Granger causality and transfer entropy to investigate linear and non-linear directed interactions, respectively, between the theta activities. Exploration in the novelty condition was associated with increased theta power in the generators with EC origin. We found a significant directed interaction from the Schaffer input over the EC-III input in CA1, and a bidirectional interaction between the inputs in the hippocampus originating in the EC, likely reflecting the connection between layers II and III. During novelty exploration, the influence of the EC-II over the EC-III generator increased, while the Schaffer influence decreased. These results associate the increase in hippocampal theta activity and synchrony during novelty exploration with an increase in the directed functional connectivity from EC-II to EC-III.

Diagnostic added value of interictal magnetic source imaging in presurgical evaluation of persons with epilepsy: A prospective blinded study

BACKGROUND AND PURPOSE

In presurgical evaluation for epilepsy surgery, information is sourced from various imaging modalities to accurately localize the epileptogenic zone. Magnetoencephalography (MEG) is a newer noninvasive technique for localization. However, there is limited literature to evaluate if MEG provides additional advantage over the conventional imaging modalities in clinical decision making. The objective of this study was to assess the diagnostic added value of MEG in decision making before epilepsy surgery.

METHOD

This was a prospective observational study. Patients underwent 3 h of recording in a MEG scanner, and the resulting localizations were compared with other complimentary investigations. Added value of MEG (considered separately from high-density electroencephalography) was defined as the frequency of cases in which (i) the information provided by magnetic source imaging (MSI) avoided implantation of intracranial electrodes and the patient was directly cleared for surgery, and (ii) MSI indicated additional substrates for implantation of intracranial electrodes. Postoperative seizure freedom was used as the diagnostic reference by which to measure the localizing accuracy of MSI.

RESULTS

A total of 102 patients underwent epilepsy surgery. MEG provided nonredundant information, which contributed to deciding the course of surgery in 33% of the patients, and prevented intracranial recordings in 19%. A total of 76% of the patients underwent surgical resection in sublobes concordant with MSI localization, and the diagnostic odds ratio for good (Engel I) outcome in these patients was 2.3 (95% confidence interval 0.68, 7.86; p = 0.183) after long-term follow-up of 36 months.

CONCLUSION

Magnetic source imaging yields additional useful information which can significantly alter as well as improve the surgical strategy for persons with epilepsy.

Dynamic analysis on simultaneous iEEG-MEG data via hidden Markov model

BACKGROUND

Intracranial electroencephalography (iEEG) recordings are used for clinical evaluation prior to surgical resection of the focus of epileptic seizures and also provide a window into normal brain function. A major difficulty with interpreting iEEG results at the group level is inconsistent placement of electrodes between subjects making it difficult to select contacts that correspond to the same functional areas. Recent work using time delay embedded hidden Markov model (HMM) applied to magnetoencephalography (MEG) resting data revealed a distinct set of brain states with each state engaging a specific set of cortical regions. Here we use a rare group dataset with simultaneously acquired resting iEEG and MEG to test whether there is correspondence between HMM states and iEEG power changes that would allow classifying iEEG contacts into functional clusters.

METHODS

Simultaneous MEG-iEEG recordings were performed at rest on 11 patients with epilepsy whose intracranial electrodes were implanted for pre-surgical evaluation. Pre-processed MEG sensor data was projected to source space. Time delay embedded HMM was then applied to MEG time series. At the same time, iEEG time series were analyzed with time-frequency decomposition to obtain spectral power changes with time. To relate MEG and iEEG results, correlations were computed between HMM probability time courses of state activation and iEEG power time course from the mid contact pair for each electrode in equally spaced frequency bins and presented as correlation spectra for the respective states and iEEG channels. Association of iEEG electrodes with HMM states based on significant correlations was compared to that based on the distance to peaks in subject-specific state topographies.

RESULTS

Five HMM states were inferred from MEG. Two of them corresponded to the left and the right temporal activations and had a spectral signature primarily in the theta/alpha frequency band. All the electrodes had significant correlations with at least one of the states (p < 0.05 uncorrected) and for 27/50 electrodes these survived within-subject FDR correction (q < 0.05). These correlations peaked in the theta/alpha band. There was a highly significant dependence between the association of states and electrodes based on functional correlations and that based on spatial proximity (p = 5.6e-6,χ2 test for independence). Despite the potentially atypical functional anatomy and physiological abnormalities related to epilepsy, HMM model estimated from the patient group was very similar to that estimated from healthy subjects.

CONCLUSION

Epilepsy does not preclude HMM analysis of interictal data. The resulting group functional states are highly similar to those reported for healthy controls. Power changes recorded with iEEG correlate with HMM state time courses in the alpha-theta band and the presence of this correlation can be related to the spatial location of electrode contacts close to the individual peaks of the corresponding state topographies. Thus, the hypothesized relation between iEEG contacts and HMM states exists and HMM could be further explored as a method for identifying comparable iEEG channels across subjects for the purposes of group analysis.

Interictal spike networks predict surgical outcome in patients with drug-resistant focal epilepsy

OBJECTIVE

To determine if properties of epileptic networks could be delineated using interictal spike propagation seen on stereo-electroencephalography (SEEG) and if these properties could predict surgical outcome in patients with drug-resistant epilepsy.

METHODS

We studied the SEEG of 45 consecutive drug-resistant epilepsy patients who underwent subsequent epilepsy surgery: 18 patients with good post-surgical outcome (Engel I) and 27 with poor outcome (Engel II-IV). Epileptic networks were derived from interictal spike propagation; these networks described the generation and propagation of interictal epileptic activity. We compared the regions in which spikes were frequent and the regions responsible for generating spikes to the area of resection and post-surgical outcome. We developed a measure termed source spike concordance, which integrates information about both spike rate and region of spike generation.

RESULTS

Inclusion in the resection of regions with high spike rate is associated with good post-surgical outcome (sensitivity = 0.82, specificity = 0.73). Inclusion in the resection of the regions responsible for generating interictal epileptic activity independently of rate is also associated with good post-surgical outcome (sensitivity = 0.88, specificity = 0.82). Finally, when integrating the spike rate and the generators, we find that the source spike concordance measure has strong predictability (sensitivity = 0.91, specificity = 0.94).

INTERPRETATIONS

Epileptic networks derived from interictal spikes can determine the generators of epileptic activity. Inclusion of the most active generators in the resection is strongly associated with good post-surgical outcome. These epileptic networks may aid clinicians in determining the area of resection during pre-surgical evaluation.

Extracting seizure onset from surface EEG with independent component analysis: Insights from simultaneous scalp and intracerebral EEG

•

Independent component analysis (ICA) is able to identify seizure generators.

•

Simultaneous long-term scalp-SEEG allows validation of the ICA results.

•

Ability to record seizure onset patterns on scalp depends on generator depth.

The success of stereoelectroencephalographic (SEEG) investigations depends crucially on the hypotheses on the putative location of the seizure onset zone. This information is derived from non-invasive data either based on visual analysis or advanced source localization algorithms. While source localization applied to interictal spikes recorded on scalp is the classical method, it does not provide unequivocal information regarding the seizure onset zone. Raw ictal activity contains a mixture of signals originating from several regions of the brain as well as EMG artifacts, hampering direct input to the source localization algorithms. We therefore introduce a methodology that disentangles the various sources contributing to the scalp ictal activity using independent component analysis and uses equivalent current dipole localization as putative locus of ictal sources. We validated the results of our analysis pipeline by performing long-term simultaneous scalp – intracerebral (SEEG) recordings in 14 patients and analyzing the wavelet coherence between the independent component encoding the ictal discharge and the SEEG signals in 8 patients passing the inclusion criteria. Our results show that invasively recorded ictal onset patterns, including low-voltage fast activity, can be captured by the independent component analysis of scalp EEG. The visibility of the ictal activity strongly depends on the depth of the sources. The equivalent current dipole localization can point to the seizure onset zone (SOZ) with an accuracy that can be as high as 10 mm for superficially located sources, that gradually decreases for deeper seizure generators, averaging at 47 mm in the 8 analyzed patients. Independent component analysis is therefore shown to have a promising SOZ localizing value, indicating whether the seizure onset zone is neocortical, and its approximate location, or located in mesial structures. That may contribute to a better crafting of the hypotheses used as basis of the stereo-EEG implantations.

Large-Scale Desynchronization During Interictal Epileptic Discharges Recorded With Intracranial EEG

It is increasingly recognized that deep understanding of epileptic seizures requires both localizing and characterizing the functional network of the region where they are initiated, i. e., the epileptic focus. Previous investigations of the epileptogenic focus' functional connectivity have yielded contrasting results, reporting both pathological increases and decreases during resting periods and seizures. In this study, we shifted paradigm to investigate the time course of connectivity in relation to interictal epileptiform discharges. We recruited 35 epileptic patients undergoing intracranial EEG (iEEG) investigation as part of their presurgical evaluation. For each patient, 50 interictal epileptic discharges (IEDs) were marked and iEEG signals were epoched around those markers. Signals were narrow-band filtered and time resolved phase-locking values were computed to track the dynamics of functional connectivity during IEDs. Results show that IEDs are associated with a transient decrease in global functional connectivity, time-locked to the peak of the discharge and specific to the high range of the gamma frequency band. Disruption of the long-range connectivity between the epileptic focus and other brain areas might be an important process for the generation of epileptic activity. Transient desynchronization could be a potential biomarker of the epileptogenic focus since 1) the functional connectivity involving the focus decreases significantly more than the connectivity outside the focus and 2) patients with good surgical outcome appear to have a significantly more disconnected focus than patients with bad outcomes.

The Ictal Signature of Thalamus and Basal Ganglia in Focal Epilepsy: A SEEG Study

OBJECTIVE

To determine the involvement of subcortical regions in human epilepsy by analyzing direct recordings from these regions during epileptic seizures using stereo-EEG (SEEG).

METHODS

We studied the SEEG recordings of a large series of patients (74 patients, 157 seizures) with an electrode sampling the thalamus and in some cases also the basal ganglia (caudate nucleus, 22 patients; and putamen, 4 patients). We applied visual analysis and signal quantification methods (Epileptogenicity Index [EI]) to their ictal recordings and compared electrophysiologic with clinical data.

RESULTS

We found that in 86% of patients, thalamus was involved during seizures (visual analysis) and 20% showed high values of epileptogenicity (EI >0.3). Basal ganglia may also disclose high values of epileptogenicity (9% in caudate nucleus) but to a lesser degree than thalamus (p < 0.01). We observed different seizure onset patterns including low voltage high frequency activities. We found high values of thalamic epileptogenicity in different epilepsy localizations, including opercular and motor epilepsies. We found no difference between epilepsy etiologies (cryptogenic vs malformation of cortical development, p = 0.77). Thalamic epileptogenicity was correlated with the extension of epileptogenic networks (p = 0.02, ρ 0.32). We found a significant effect (p < 0.05) of thalamic epileptogenicity regarding the postsurgical outcome (higher thalamic EI corresponding to higher probability of surgical failure).

CONCLUSIONS

Thalamic involvement during seizures is common in different seizure types. The degree of thalamic epileptogenicity is a possible marker of the epileptogenic network extension and of postsurgical prognosis.

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.

EMS-Net: A Deep Learning Method for Autodetecting Epileptic Magnetoencephalography Spikes

Epilepsy is a neurological disorder characterized by sudden and unpredictable epileptic seizures, which incurs significant negative impacts on patients' physical, psychological and social health. A practical approach to assist with the clinical assessment and treatment planning for patients is to process magnetoencephalography (MEG) data to identify epileptogenic zones. As a widely accepted biomarker of epileptic foci, epileptic MEG spikes need to be precisely detected. Given that the visual inspection of spikes is time consuming, an automatic and efficient system with adequate accuracy for spike detection is valuable in clinical practice. However, current approaches for MEG spike autodetection are dependent on hand-engineered features. Here, we propose a novel multiview Epileptic MEG Spikes detection algorithm based on a deep learning Network (EMS-Net) to accurately and efficiently recognize the spike events from MEG raw data. The results of the leave-k-subject-out validation tests for multiple datasets (i.e., balanced and realistic datasets) showed that EMS-Net achieved state-of-the-art classification performance (i.e., accuracy: 91.82% - 99.89%; precision: 91.90% - 99.45%; sensitivity: 91.61% - 99.53%; specificity: 91.60% - 99.96%; f1 score: 91.70% - 99.48%; and area under the curve: 0.9688 - 0.9998).

A novel method for extracting interictal epileptiform discharges in multi-channel MEG: Use of fractional type of blind source separation

OBJECTIVE

Visual inspection of interictal epileptiform discharges (IEDs) in multi-channel MEG requires a time-consuming evaluation process and often leads to inconsistent results due to variability of IED waveforms. Here, we propose a novel extraction method for IEDs using a T/k type of blind source separation (BSS_T/k).

METHODS

We applied BSS_T/k with seven patients with focal epilepsy to test the accuracy of identification of IEDs. We conducted comparisons of the results of BSS components with those obtained by visual inspection in sensor-space analysis.

RESULTS

BSS_T/k provided better signal estimation of IEDs compared with sensor-space analysis. Importantly, BSS_T/k was able to uncover IEDs that could not be detected by visual inspection. Furthermore, IED components were clearly extracted while preserving spike and wave morphology. Variable IED waveforms were decomposed into one dominant component.

CONCLUSIONS

BSS_T/k was able to visualize the spreading signals over multiple channels into a single component from a single epileptogenic zone. BSS_T/k can be applied to focal epilepsy with a simple parameter setting.

SIGNIFICANCE

Our novel method was able to highlight IEDs with increased accuracy for identification of IEDs from multi-channel MEG data.

Evidence for the Role of Magnetic Source Imaging in the Presurgical Evaluation of Refractory Epilepsy Patients

Magnetoencephalography (MEG) in the field of epilepsy has multiple advantages; just like electroencephalography (EEG), MEG is able to measure the epilepsy specific information (i.e., the brain activity reflecting seizures and/or interictal epileptiform discharges) directly, non-invasively and with a very high temporal resolution (millisecond-range). In addition MEG has a unique sensitivity for tangential sources, resulting in a full picture of the brain activity when combined with EEG. It accurately allows to perform source imaging of focal epileptic activity and functional cortex and shows a specific high sensitivity for a source in the neocortex. In this paper the current evidence and practice for using magnetic source imaging of focal interictal and ictal epileptic activity during the presurgical evaluation of drug resistant patients is being reviewed.

Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies

Drug-refractory focal epilepsies are network diseases associated with functional connectivity alterations both during ictal and interictal periods. A large majority of studies on the interictal/resting state have focused on functional MRI-based functional connectivity. Few studies have used electrophysiology, despite its high temporal capacities. In particular, stereotactic-EEG is highly suitable to study functional connectivity because it permits direct intracranial electrophysiological recordings with relative large-scale sampling. Most previous studies in stereotactic-EEG have been directed towards temporal lobe epilepsy, which does not represent the whole spectrum of drug-refractory epilepsies. The present study aims at filling this gap, investigating interictal functional connectivity alterations behind cortical epileptic organization and its association with post-surgical prognosis. To this purpose, we studied a large cohort of 59 patients with malformation of cortical development explored by stereotactic-EEG with a wide spatial sampling (76 distinct brain areas were recorded, median of 13.2 per patient). We computed functional connectivity using non-linear correlation. We focused on three zones defined by stereotactic-EEG ictal activity: the epileptogenic zone, the propagation zone and the non-involved zone. First, we compared within-zone and between-zones functional connectivity. Second, we analysed the directionality of functional connectivity between these zones. Third, we measured the associations between functional connectivity measures and clinical variables, especially post-surgical prognosis. Our study confirms that functional connectivity differs according to the zone under investigation. We found: (i) a gradual decrease of the within-zone functional connectivity with higher values for epileptogenic zone and propagation zone, and lower for non-involved zones; (ii) preferential coupling between structures of the epileptogenic zone; (iii) preferential coupling between epileptogenic zone and propagation zone; and (iv) poorer post-surgical outcome in patients with higher functional connectivity of non-involved zone (within- non-involved zone, between non-involved zone and propagation zone functional connectivity). Our work suggests that, even during the interictal state, functional connectivity is reinforced within epileptic cortices (epileptogenic zone and propagation zone) with a gradual organization. Moreover, larger functional connectivity alterations, suggesting more diffuse disease, are associated with poorer post-surgical prognosis. This is consistent with computational studies suggesting that connectivity is crucial in order to model the spatiotemporal dynamics of seizures.10.1093/brain/awy214_video1awy214media15833456182001.

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.

Reproducibility of interictal spike propagation in children with refractory epilepsy

OBJECTIVE

Interictal spikes are a characteristic feature of invasive electroencephalography (EEG) recordings in children with refractory epilepsy. Spikes frequently co-occur across multiple brain regions with discernable latencies, suggesting that spikes can propagate through distributed neural networks. The purpose of this study was to examine the long-term reproducibility of spike propagation patterns over hours to days of interictal recording.

METHODS

Twelve children (mean age 13.1 years) were retrospectively studied. A mean ± standard deviation (SD) of 47.2 ± 40.1 hours of interictal EEG recordings were examined per patient (range 17.5-166.5 hours). Interictal recordings were divided into 30-minute segments. Networks were extracted based on the frequency of spike coactivation between pairs of electrodes. For each 30-minute segment, electrodes were assigned a "Degree Preference (DP)" based on the tendency to appear upstream or downstream within propagation sequences. The consistency of DPs across segments ("DP-Stability") was quantified using the Spearman rank correlation.

RESULTS

Regions exhibited highly stable preferences to appear upstream, intermediate, or downstream in spike propagation sequences. Across networks, the mean ± SD DP-Stability was 0.88 ± 0.07, indicating that propagation patterns observed in 30-minute segments were representative of the patterns observed in the full interictal window. At the group level, regions involved in seizure generation appeared more upstream in spike propagation sequences.

SIGNIFICANCE

Interictal spike propagation is a highly reproducible output of epileptic networks. These findings shed new light on the spatiotemporal dynamics that may constrain the network mechanisms of refractory epilepsy.

High-Frequency Oscillation Networks and Surgical Outcome in Adult Focal Epilepsy

OBJECTIVE

To investigate whether high-frequency oscillations (HFOs) show spatiotemporal propagation and assess the relevance of the earliest oscillations in relation to the seizure onset zone (SOZ) and postsurgical outcome.

METHODS

We retrospectively investigated the intracerebral electroencephalography (EEG) of patients who became seizure free after subsequent surgery. We marked HFOs during 1 hour of recordings. We calculated the time delay between pairs of channels as the median delay between their HFOs and constructed a time line of the delay of each channel with respect to the earliest channel (first source channel). A network was defined when a temporal order could be established among the channels based on the existence of statistically significant delays.

RESULTS

Fifteen patients with good surgical outcome were included. We found ripple networks in all patients, and fast ripple networks in 9. For ripples, first source channels were found in a higher proportion in the SOZ than the rest of the network channels (15 of 27 [56%] versus 93 of 262 [35%]; p = 0.04). For both ripples and fast ripples, first source channels were resected more often that the rest of the network channels (ripples: 13 of 27 [48%] versus 65 of 262 [25%]; p = 0.01; fast ripples: 8 of 9 [89%] versus 17 of 40 [43%]; p = 0.002); channels with the highest rates of ripples and fast ripples were resected in a similar proportion.

INTERPRETATION

These results demonstrate that interictal HFOs are organized in networks and indicate a possible need for the resection of first source channels. However, resecting them is not superior to resecting channels with highest rates of HFOs. Ann Neurol 2019;85:485-494.

Deep brain activities can be detected with magnetoencephalography

The hippocampus and amygdala are key brain structures of the medial temporal lobe, involved in cognitive and emotional processes as well as pathological states such as epilepsy. Despite their importance, it is still unclear whether their  neural activity can be recorded non-invasively. Here, using simultaneous intracerebral and magnetoencephalography (MEG) recordings in patients with focal drug-resistant epilepsy, we demonstrate a direct contribution of amygdala and hippocampal activity to surface MEG recordings. In particular, a method of blind source separation, independent component analysis, enabled activity arising from large neocortical networks to be disentangled from that of deeper structures, whose amplitude at the surface was small but significant. This finding is highly relevant for our understanding of hippocampal and amygdala brain activity as it implies that their activity could potentially be measured non-invasively.

Presurgical intracranial investigations in epilepsy surgery

Identification and localization of the "epileptogenic process" in the brain of patients with drug-resistant epilepsy for surgical cure is the goal of presurgical investigations. Intracranial recordings are required when conflicting data between seizure clinical semiology and EEG prevent precise localization within one hemisphere or lateralization, when a visible lesion on MRI seems unrelated to the electroclinical data, or in MRI-negative cases. Two methods are currently used. The objective of the subdural grid electrocorticography with or without depth electrodes (SDG/DE) is the best possible identification of the area of onset of spontaneous seizures and localization of the eloquent cortex. The objective of stereoelectroencephalography (SEEG) is to define the epileptogenic zone (configured as a network) and its relation to an unmasked lesion. Two-dimensional (SDG) and three-dimensional (SEEG) brain sampling dictate different strategies for noninvasive presurgical phase I goals as well as for data analysis. SEEG must resolve several potential localization hypotheses in a manner that cannot be achieved with SDG. SDG operates through brain surface coverage, unlike SEEG, which samples networks. SDG estimates the extent of cortical resection through a lobar or sublobar localization of ictal onset and constraints from functional mapping. SEEG defines a tailored resection according to the results of anatomo-electro-clinical correlations in stereotaxic space that will guide the ablation of the epileptogenic zone. SEEG is currently expanding faster than SDG. The prerequisites (especially in the preimplantation hypothetical strategy) and technical tools (especially stimulation and functional mapping) in the two methods are very different. This chapter presents a comparative review of the rationale, indications, electrode implantation strategies, interpretation, and surgical decision making of these two approaches of presurgical evaluation for epilepsy surgery.

Brain regions and epileptogenicity influence epileptic interictal spike production and propagation during NREM sleep in comparison with wakefulness

OBJECTIVE

Non-rapid eye movement (NREM) sleep is known to be a brain state associated with an activation of interictal epileptic activity. The goal of this work was to quantify topographic changes occurring during NREM sleep in comparison with wakefulness.

METHOD

We studied intracerebral recordings of 20 patients who underwent stereo-electroencephalography (SEEG) during presurgical evaluation for pharmacoresistant focal epilepsy. We measured the number of interictal spikes (IS) and quantified the co-occurrence of IS between brain regions during 1 hour of NREM sleep and 1 hour of wakefulness. Co-occurrence is a method to estimate IS networks based on a temporal concordance between IS of different brain regions. Each studied region was labeled as "seizure-onset zone" (SOZ), "propagation zone" (PZ), or "not involved region" (NIR).

RESULTS

During NREM sleep, the number of interictal spikes significantly increased in all regions (mean of 68%). This increase was higher in medial temporal regions than in other regions, whether involved in the SOZ. Spike co-occurrence increased significantly in all regions during NREM sleep in comparison with wakefulness but was greater in neocortical regions. Spike co-occurrence in medial temporal regions was not higher than in other regions, suggesting that the increase of the number of spikes in this region was in great part a local effect.

SIGNIFICANCE

This study demonstrated that medial temporal regions show a greater propensity to spike production or propagation during NREM sleep compared to other brain regions, even when the medial temporal lobe is not involved in the SOZ.

Technical solutions for simultaneous MEG and SEEG recordings: towards routine clinical use

OBJECTIVE

The simultaneous recording of intracerebral EEG (stereotaxic EEG, SEEG) and magnetoencephalography (MEG) is a promising strategy that provides both local and global views on brain pathological activity. Yet, acquiring simultaneous signals poses difficult technical issues that hamper their use in clinical routine. Our objective was thus to develop a set of solutions for recording a high number of SEEG channels while preserving signal quality.

APPROACH

We recorded data in a patient with drug resistant epilepsy during presurgical evaluation. We used dedicated insertion screws and optically insulated amplifiers. We recorded 137 SEEG contacts on 10 depth electrodes (5-15 contacts each) and 248 MEG channels (magnetometers). Signal quality was assessed by comparing the distribution of RMS values in different frequency bands to a reference set of MEG acquisitions.

MAIN RESULTS

The quality of signals was excellent for both MEG and SEEG; for MEG, it was comparable to that of MEG signals without concurrent SEEG. Discharges involving several structures on SEEG were visible on MEG, whereas discharges limited in space were not seen at the surface.

SIGNIFICANCE

SEEG can now be recorded simultaneously with whole-head MEG in routine. This opens new avenues, both methodologically for understanding signals and improving signal processing methods, and clinically for future combined analyses.

Magnetoencephalography in Preoperative Epileptic Foci Localization: Enlightenment from Cognitive Studies

Over 30% epileptic patients are refractory to medication, who are amenable to neurosurgical treatment. Non-invasive brain imaging technologies including video-electroencephalogram (EEG), magnetic resonance imaging (MRI), and magnetoencephalography (MEG) are widely used in presurgical assessment of epileptic patients. This review mainly discussed the current development of clinical MEG imaging as a diagnose approach, and its correlations with the golden standard intracranial electroencephalogram (iEEG). More importantly, this review discussed the possible applications of functional networks in preoperative epileptic foci localization in future studies.

Multimodal investigation of epileptic networks: The case of insular cortex epilepsy

The insula is a deep cortical structure sharing extensive synaptic connections with a variety of brain regions, including several frontal, temporal, and parietal structures. The identification of the insular connectivity network is obviously valuable for understanding a number of cognitive processes, but also for understanding epilepsy since insular seizures involve a number of remote brain regions. Ultimately, knowledge of the structure and causal relationships within the epileptic networks associated with insular cortex epilepsy can offer deeper insights into this relatively neglected type of epilepsy enabling the refining of the clinical approach in managing patients affected by it. In the present chapter, we first review the multimodal noninvasive tests performed during the presurgical evaluation of epileptic patients with drug refractory focal epilepsy, with particular emphasis on their value for the detection of insular cortex epilepsy. Second, we review the emerging multimodal investigation techniques in the field of epilepsy, that aim to (1) enhance the detection of insular cortex epilepsy and (2) unveil the architecture and causal relationships within epileptic networks. We summarize the results of these approaches with emphasis on the specific case of insular cortex epilepsy.