Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy

Assessing cortical excitability with electroencephalography: A pilot study with EEG-iTBS

BACKGROUND

Cortical excitability measures neural reactivity to stimuli, usually delivered via Transcranial Magnetic Stimulation (TMS). Excitation/inhibition balance (E/I) is the ongoing equilibrium between excitatory and inhibitory activity of neural circuits. According to some studies, E/I could be estimated in-vivo and non-invasively through the modeling of electroencephalography (EEG) signals and termed 'intrinsic excitability' measures. Several measures have been proposed (phase consistency in the gamma band, sample entropy, exponent of the power spectral density 1/f curve, E/I index extracted from detrend fluctuation analysis, and alpha power). Intermittent theta burst stimulation (iTBS) of the primary motor cortex (M1) is a non-invasive neuromodulation technique allowing controlled and focal enhancement of TMS cortical excitability and E/I of the stimulated hemisphere.

OBJECTIVE

Investigating to what extent E/I estimates scale with TMS excitability and how they relate to each other.

METHODS

M1 excitability (TMS) and several E/I estimates extracted from resting state EEG recordings were assessed before and after iTBS in a cohort of healthy subjects.

RESULTS

Enhancement of TMS M1 excitability, as measured through motor-evoked potentials (MEPs), and phase consistency of the cortex in high gamma band correlated with each other. Other measures of E/I showed some expected results, but no correlation with TMS excitability measures or strong consistency with each other.

CONCLUSIONS

EEG E/I estimates offer an intriguing opportunity to map cortical excitability non-invasively, with high spatio-temporal resolution and with a stimulus independent approach. While different EEG E/I estimates may reflect the activity of diverse excitatory-inhibitory circuits, spatial phase synchrony in the gamma band is the measure that best captures excitability changes in the primary motor cortex.

An artificial intelligence-based pipeline for automated detection and localisation of epileptic sources from magnetoencephalography

Objective.Magnetoencephalography (MEG) is a powerful non-invasive diagnostic modality for presurgical epilepsy evaluation. However, the clinical utility of MEG mapping for localising epileptic foci is limited by its low efficiency, high labour requirements, and considerable interoperator variability. To address these obstacles, we proposed a novel artificial intelligence-based automated magnetic source imaging (AMSI) pipeline for automated detection and localisation of epileptic sources from MEG data.Approach.To expedite the analysis of clinical MEG data from patients with epilepsy and reduce human bias, we developed an autolabelling method, a deep-learning model based on convolutional neural networks and a hierarchical clustering method based on a perceptual hash algorithm, to enable the coregistration of MEG and magnetic resonance imaging, the detection and clustering of epileptic activity, and the localisation of epileptic sources in a highly automated manner. We tested the capability of the AMSI pipeline by assessing MEG data from 48 epilepsy patients.Main results.The AMSI pipeline was able to rapidly detect interictal epileptiform discharges with 93.31% ± 3.87% precision based on a 35-patient dataset (with sevenfold patientwise cross-validation) and robustly rendered accurate localisation of epileptic activity with a lobar concordance of 87.18% against interictal and ictal stereo-electroencephalography findings in a 13-patient dataset. We also showed that the AMSI pipeline accomplishes the necessary processes and delivers objective results within a much shorter time frame (∼12 min) than traditional manual processes (∼4 h).Significance.The AMSI pipeline promises to facilitate increased utilisation of MEG data in the clinical analysis of patients with epilepsy.

Incorporation of Fast-Elimination Chemicals in Hair Is Governed by Pharmacokinetics-Implications for Exposure Assessment

Mechanisms governing chemicals' incorporation in hair are incompletely understood, and gaps remain to link the concentration of chemicals in hair to level of exposure and internal dose present in the body. This study assesses the relevance of hair analysis for the biomonitoring of exposure to fast-elimination compounds and investigates the role of pharmacokinetics (PK) in their incorporation in hair. Rats were administered with pesticides, bisphenols, phthalates, and DINCH over 2 months. Hairs were analyzed for 28 chemicals/metabolites to investigate correlations between their concentration in hair and the dose administered to the animals. Urine collected over 24 h after gavage was used to determine chemicals' PK and to investigate their influence on incorporation into hair by means of linear mixed models (LMMs). Eighteen chemicals presented a significant correlation between concentration in hair and level of exposure. In models combining all chemicals, agreement between concentration in hair predicted by LMM and experimental values was moderate (R² = 0.19) but significantly increased when PK were included in the models (R² = 0.37), and even more when chemical families were considered separately (e.g., R² = 0.98 for pesticides). This study shows that pharmacokinetics mediate incorporation of chemicals in hair and suggests the relevance of hair for assessing exposure to fast-elimination chemicals.

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.

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

BACKGROUND AND OBJECTIVES

Accurate delineation of the seizure-onset zone (SOZ) in focal drug-resistant epilepsy often requires stereo-EEG (SEEG) recordings. Our aims were to propose a truly objective and quantitative comparison between EEG/magnetoencephalography (MEG) source imaging (EMSI), EEG/fMRI responses for similar spikes with primary irritative zone (PIZ) and SOZ defined by SEEG and to evaluate the value of EMSI and EEG/fMRI to predict postsurgical outcome.

METHODS

We identified patients with drug-resistant epilepsy who underwent EEG/MEG, EEG/fMRI, and subsequent SEEG at the Epilepsy Service from the Montreal Neurological Institute and Hospital. We quantified multimodal concordance within the SEEG channel space as spatial overlap with PIZ/SOZ and distances to the spike-onset, spike maximum amplitude and seizure core intracerebral channels by applying a new methodology consisting of converting EMSI results into SEEG electrical potentials (EMSIe-SEEG) and projecting the most significant fMRI response on the SEEG channels (fMRIp-SEEG). Spatial overlaps with PIZ/SOZ (AUCPIZ, AUCSOZ) were assessed by using the area under the receiver operating characteristic curve (AUC). Here, AUC represents the probability that a randomly picked active contact exhibited higher amplitude when located inside the spatial reference than outside.

RESULTS

Seventeen patients were included. Mean spatial overlaps with the PIZ and SOZ were 0.71 and 0.65 for EMSIe-SEEG and 0.57 and 0.62 for fMRIp-SEEG. Good EMSIe-SEEG spatial overlap with the PIZ was associated with smaller distance from the maximum EMSIe-SEEG contact to the spike maximum amplitude channel (median distance 14 mm). Conversely, good fMRIp-SEEG spatial overlap with the SOZ was associated with smaller distances from the maximum fMRIp-SEEG contact to the spike-onset and seizure core channels (median distances 10 and 5 mm, respectively). Surgical outcomes were correctly predicted by EEG/MEG in 12 of 15 (80%) patients and EEG/fMRI in 6 of 11(54%) patients.

DISCUSSION

With the use of a unique quantitative approach estimating EMSI and fMRI results in the reference SEEG channel space, EEG/MEG and EEG/fMRI accurately localized the SOZ and the PIZ. Precisely, EEG/MEG more accurately localized the PIZ, whereas EEG/fMRI was more sensitive to the SOZ. Both neuroimaging techniques provide complementary localization that can help guide SEEG implantation and select good candidates for surgery.

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean

Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {cm}^2$$\end{document}cm2 along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.

Resting state network connectivity is attenuated by fMRI acoustic noise

INTRODUCTION

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

METHODS

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

RESULTS

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

DISCUSSION

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

Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean

In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.

Functional Neurophysiological Biomarkers of Early-Stage Alzheimer's Disease: A Perspective of Network Hyperexcitability in Disease Progression

Network hyperexcitability (NH) has recently been suggested as a potential neurophysiological indicator of Alzheimer’s disease (AD), as new, more accurate biomarkers of AD are sought. NH has generated interest as a potential indicator of certain stages in the disease trajectory and even as a disease mechanism by which network dysfunction could be modulated. NH has been demonstrated in several animal models of AD pathology and multiple lines of evidence point to the existence of NH in patients with AD, strongly supporting the physiological and clinical relevance of this readout. Several hypotheses have been put forward to explain the prevalence of NH in animal models through neurophysiological, biochemical, and imaging techniques. However, some of these hypotheses have been built on animal models with limitations and caveats that may have derived NH through other mechanisms or mechanisms without translational validity to sporadic AD patients, potentially leading to an erroneous conclusion of the underlying cause of NH occurring in patients with AD. In this review, we discuss the substantiation for NH in animal models of AD pathology and in human patients, as well as some of the hypotheses considering recently developed animal models that challenge existing hypotheses and mechanisms of NH. In addition, we provide a preclinical perspective on how the development of animal models incorporating AD-specific NH could provide physiologically relevant translational experimental data that may potentially aid the discovery and development of novel therapies for AD.

Exploring the extent of source imaging: Recent advances in noninvasive electromagnetic brain imaging

Electrophysiological source imaging (ESI) has been successfully employed in many brain imaging applications during the last 20 years. ESI estimates of underlying brain networks provide millisecond resolution of dynamic brain processes; yet, it remains to be a challenge to further improve the spatial resolution of ESI modality, in particular on its capability of imaging the extent of underlying brain sources. In this review, we discuss the recent developments in signal processing and machine learning that have made it possible to image the extent, i.e. size, of underlying brain sources noninvasively, using scalp electromagnetic measurements from electroencephalogram (EEG) and magnetoencephalogram (MEG) recordings.

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

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

Localization of Epileptic Foci Based on Simultaneous EEG-fMRI Data

Combining functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) enables a non-invasive investigation of the human brain function and evaluation of the correlation of these two important modalities of brain activity. This paper explores recent reports on using advanced simultaneous EEG–fMRI methods proposed to map the regions and networks involved in focal epileptic seizure generation. One of the applications of EEG and fMRI combination as a valuable clinical approach is the pre-surgical evaluation of patients with epilepsy to map and localize the precise brain regions associated with epileptiform activity. In the process of conventional analysis using EEG–fMRI data, the interictal epileptiform discharges (IEDs) are visually extracted from the EEG data to be convolved as binary events with a predefined hemodynamic response function (HRF) to provide a model of epileptiform BOLD activity and use as a regressor for general linear model (GLM) analysis of the fMRI data. This review examines the methodologies involved in performing such studies, including techniques used for the recording of EEG inside the scanner, artifact removal, and statistical analysis of the fMRI signal. It then discusses the results reported for patients with primary generalized epilepsy and patients with different types of focal epileptic disorders. An important matter that these results have brought to light is that the brain regions affected by interictal epileptic discharges might not be limited to the ones where they have been generated. The developed methods can help reveal the regions involved in or affected by a seizure onset zone (SOZ). As confirmed by the reviewed literature, EEG–fMRI provides information that comes particularly useful when evaluating patients with refractory epilepsy for surgery.

Ant Colony System Optimization for Spatiotemporal Modelling of Combined EEG and MEG Data

Electroencephalography/Magnetoencephalography (EEG/MEG) source localization involves the estimation of neural activity inside the brain volume that underlies the EEG/MEG measures observed at the sensor array. In this paper, we consider a Bayesian finite spatial mixture model for source reconstruction and implement Ant Colony System (ACS) optimization coupled with Iterated Conditional Modes (ICM) for computing estimates of the neural source activity. Our approach is evaluated using simulation studies and a real data application in which we implement a nonparametric bootstrap for interval estimation. We demonstrate improved performance of the ACS-ICM algorithm as compared to existing methodology for the same spatiotemporal model.

Magnetoencephalography: physics, techniques, and applications in the basic and clinical neurosciences

Magnetoencephalography (MEG) is a technique used to measure the magnetic fields generated from neuronal activity in the brain. MEG has a high temporal resolution on the order of milliseconds and provides a more direct measure of brain activity when compared with hemodynamic-based neuroimaging methods such as magnetic resonance imaging and positron emission tomography. The current review focuses on basic features of MEG such as the instrumentation and the physics that are integral to the signals that can be measured, and the principles of source localization techniques, particularly the physics of beamforming and the techniques that are used to localize the signal of interest. In addition, we review several metrics that can be used to assess functional coupling in MEG and describe the advantages and disadvantages of each approach. Lastly, we discuss the current and future applications of MEG.

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

OBJECTIVE

Fast Oscillations (FO) >40 Hz are a promising biomarker of the epileptogenic zone (EZ). Evidence using scalp electroencephalography (EEG) remains scarce. We assessed if electrical source imaging of FO using 256-channel high-density EEG (HD-EEG) is useful for EZ identification.

METHODS

We analyzed HD-EEG recordings of 10 focal drug-resistant epilepsy patients with seizure-free postsurgical outcome. We marked FO candidate events at the time of epileptic spikes and verified them by screening for an isolated peak in the time-frequency plot. We performed electrical source imaging of spikes and FO within the Maximum Entropy of the Mean framework. Source localization maps were validated against the surgical cavity.

RESULTS

We identified FO in five out of 10 patients who had a superficial or intermediate deep generator. The maximum of the FO maps was localized inside the cavity in all patients (100%). Analysis with a reduced electrode coverage using the 10-10 and 10-20 system showed a decreased localization accuracy of 60% and 40% respectively.

CONCLUSIONS

FO recorded with HD-EEG localize the EZ. HD-EEG is better suited to detect and localize FO than conventional EEG approaches.

SIGNIFICANCE

This study acts as proof-of-concept that FO localization using 256-channel HD-EEG is a viable marker of the EZ.

Quantitative EEG biomarkers for epilepsy and their relation to chemical biomarkers

The electroencephalogram (EEG) is the most important method to diagnose epilepsy. In clinical settings, it is evaluated by experts who identify patterns visually. Quantitative EEG is the application of digital signal processing to clinical recordings in order to automatize diagnostic procedures, and to make patterns visible that are hidden to the human eye. The EEG is related to chemical biomarkers, as electrical activity is based on chemical signals. The most well-known chemical biomarkers are blood laboratory tests to identify seizures after they have happened. However, research on chemical biomarkers is much less extensive than research on quantitative EEG, and combined studies are rarely published, but highly warranted. Quantitative EEG is as old as the EEG itself, but still, the methods are not yet standard in clinical practice. The most evident application is an automation of manual work, but also a quantitative description and localization of interictal epileptiform events as well as seizures can reveal important hints for diagnosis and contribute to presurgical evaluation. In addition, the assessment of network characteristics and entropy measures were found to reveal important insights into epileptic brain activity. Application scenarios of quantitative EEG in epilepsy include seizure prediction, pharmaco-EEG, treatment monitoring, evaluation of cognition, and neurofeedback. The main challenges to quantitative EEG are poor reliability and poor generalizability of measures, as well as the need for individualization of procedures. A main hindrance for quantitative EEG to enter clinical routine is also that training is not yet part of standard curricula for clinical neurophysiologists.

Enhanced Fast-VESTAL for Magnetoencephalography Source Imaging: From Theory to Clinical Application in Epilepsy

A novel magnetoencephalography source imaging approach called Fast Vector-based Spatio-Temporal Analysis (Fast-VESTAL) has been successfully applied in creating source images from evoked and resting-state data from both healthy subjects and individuals with neurological and/or psychiatric disorders, but its reconstructed source images may show false-positive activations, especially under low signal-to-noise ratio conditions. Here, to effectively reduce false-positive artifacts, we introduced an enhanced Fast-VESTAL (eFast-VESTAL) approach that adopts generalized second-order cone programming. We compared the spatiotemporal characteristics of the eFast-VESTAL approach to those of the popular distributed source approaches (e.g., the minimum L2-norm/ mixed-norm methods) using computer simulations and auditory experiments. More importantly, we applied eFast-VESTAL to the presurgical evaluation of epilepsy. Our results demonstrated that eFast-VESTAL exhibited a lower dipole localization error and/or a higher correlation coefficient (CC) between the estimated source time series and ground truth under various conditions of source waveforms. Experimentally, eFast-VESTAL displayed more focal activation maps and a higher CC between the raw and predicted sensor data in response to auditory stimulation. Notably, eFast-VESTAL was the most accurate method for noninvasively detecting the epileptic zones determined using more invasive stereo-electroencephalography in the comparison.

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.

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).

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

Cortical gamma-synchrony measured with magnetoencephalography is a marker of clinical status and predicts clinical outcome in stroke survivors

•

The outcome of stroke survivors is difficult to anticipate.

•

Gamma synchrony is a reliable measure of brain function and reserve.

•

Gamma synchrony is measured with MEG in stroke survivors undergoing rehab.

•

Auditory-entrained gamma synchrony correlates with clinical status and outcome.

Background

The outcome of stroke survivors is difficult to anticipate. While the extent of the anatomical brain lesion is only poorly correlated with the prognosis, functional measures of cortical synchrony, brain networks and cortical plasticity seem to be good predictors of clinical recovery. In this field, gamma (>30 Hz) cortical synchrony is an ideal marker of brain function, as it plays a crucial role for the integration of information, it is an indirect marker of Glutamate/GABA balance and it directly estimates the reserve of parvalbulin-positive neurons, key players in synaptic plasticity. In this study we measured gamma synchronization driven by external auditory stimulation with magnetoencephalography and tested whether it was predictive of the clinical outcome in stroke survivors undergoing intensive rehabilitation in a tertiary rehabilitation center.

Material and methods

Eleven stroke survivors undergoing intensive rehabilitation were prospectively recruited. Gamma synchrony was measured non-invasively within one month from stroke onset with magnetoencephalography, both at rest and during entrainment with external 40 Hz amplitude modulated binaural sounds. Lesion location and volume were quantitatively assessed through a high-resolution anatomical MRI. Barthel index (BI) and Functional Independence Measure (FIM) scales were measured at the beginning and at the end of the admission to the rehabilitation unit.

Results

The spatial distribution of cortical gamma synchrony was altered, and the physiological right hemispheric dominance observed in healthy controls was attenuated or lost. Entrained gamma synchronization (but not resting state gamma synchrony) showed a very high correlation with the clinical status at both admission and discharge (both BI and FIM). Neither clinical status nor gamma synchrony showed a correlation with lesion volume.

Conclusions

Cortical gamma synchrony related to auditory entrainment can be reliably measured in stroke patients. Gamma synchrony is strongly associated with the clinical outcome of stroke survivors undergoing rehabilitation.

Transcranial direct current stimulation over the sensory-motor regions inhibits gamma synchrony

Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique able to induce plasticity phenomena. Although tDCS application has been spreading over a variety of neuroscience domains, the mechanisms by which the stimulation acts are largely unknown. We investigated tDCS effects on cortical gamma synchrony, which is a crucial player in cortical function. We performed a randomized, sham-controlled, double-blind study on healthy subjects, combining tDCS and magnetoencephalography. By driving brain activity via 40 Hz auditory stimulation during magnetoencephalography, we experimentally tuned cortical gamma synchrony and measured it before and after bilateral tDCS of the primary sensory-motor hand regions (anode left, cathode right). We demonstrated that the stimulation induces a remarkable decrease of gamma synchrony (13 out of 15 subjects), as measured by gamma phase at 40 Hz. tDCS has strong remote effects, as the cortical region mostly affected was located far away from the stimulation site and covered a large area of the right centro-temporal cortex. No significant differences between stimulations were found for baseline gamma synchrony, as well as early transient auditory responses. This suggests a specific tDCS effect on externally driven gamma synchronization. This study sheds new light on the effect of tDCS on cortical function showing that the net effect of the stimulation on cortical gamma synchronization is an inhibition.

Diagnostic accuracy of interictal source imaging in presurgical epilepsy evaluation: A systematic review from the E-PILEPSY consortium

OBJECTIVE

Interictal high resolution (HR-) electric source imaging (ESI) and magnetic source imaging (MSI) are non-invasive tools to aid epileptogenic zone localization in epilepsy surgery candidates. We carried out a systematic review on the diagnostic accuracy and quality of evidence of these modalities.

METHODS

Embase, Pubmed and the Cochrane database were searched on 13 February 2017. Diagnostic accuracy studies taking post-surgical seizure outcome as reference standard were selected. Quality appraisal was based on the QUADAS-2 framework.

RESULTS

Eleven studies were included: eight MSI (n = 267), three HR-ESI (n = 127) studies. None was free from bias. This mostly involved: selection of operated patients only, interference of source imaging with surgical decision, and exclusion of indeterminate results. Summary sensitivity and specificity estimates were 82% (95% CI: 75-88%) and 53% (95% CI: 37-68%) for overall source imaging, with no statistical difference between MSI and HR-ESI. Specificity is higher when partially concordant results were included as non-concordant (p < 0.05). Inclusion of indeterminate test results as non-concordant lowered sensitivity (p < 0.05).

CONCLUSIONS

Source imaging has a relatively high sensitivity but low specificity for identification of the epileptogenic zone.

SIGNIFICANCE

We need higher quality studies allowing unbiased test evaluation to determine the added value and diagnostic accuracy of source imaging in the presurgical workup of refractory focal epilepsy.

Widespread interictal epileptic discharge more likely than focal discharges to unveil the seizure onset zone in EEG-fMRI

OBJECTIVE

We hypothesized that the number of interictal epileptic discharges (IEDs) during scan and their spatial extent are contributing factors in obtaining appropriate activations that reveal the seizure onset zone (SOZ) in EEG-fMRI.

METHODS

157 IED types, each corresponding to one EEG scalp distribution, in 64 consecutive EEG-fMRI studies from 64 patients with refractory localization-related epilepsy were reviewed. To determine reliable activation, we used the threshold corresponding to corrected whole-brain topological false discovery rate (FDR). The location with maximum activation was compared to the presumed SOZ as defined by a comprehensive evaluation for each patient.

RESULTS

The number of IEDs was significantly higher in the types with t-value above FDR than with t-value below FDR. The presumed SOZ could be delineated in 30 of the 64 patients. Among these patients, the types of IED concordant with the SOZ had significantly larger extent on scalp EEG than the IED types discordant with the SOZ.

CONCLUSIONS

The number of IEDs is important factor in obtaining reliable activations in EEG-fMRI. IEDs with larger spatial extent are more likely to reveal, on maximum BOLD, accurate location of the SOZ.

SIGNIFICANCE

Widespread discharges are more likely to yield a reliable activation for SOZ in EEG-fMRI.

Theta Activity in the Left Dorsal Premotor Cortex During Action Re-Evaluation and Motor Reprogramming

The ability to rapidly adjust our actions to changes in the environment is a key function of human motor control. Previous work implicated the dorsal premotor cortex (dPMC) in the up-dating of action plans based on environmental cues. Here we used electroencephalography (EEG) to identify neural signatures of up-dating cue-action relationships in the dPMC and connected frontoparietal areas. Ten healthy subjects performed a pre-cued alternate choice task. Simple geometric shapes cued button presses with the right or left index finger. The shapes of the pre-cue and go-cue differed in two third of trials. In these incongruent trials, the go-cue prompted a re-evaluation of the pre-cued action plan, slowing response time relative to trials with identical cues. This re-evaluation selectively increased theta band activity without modifying activity in alpha and beta band. Source-based analysis revealed a widespread theta increase in dorsal and mesial frontoparietal areas, including dPMC, supplementary motor area (SMA), primary motor and posterior parietal cortices (PPC). Theta activity scaled positively with response slowing and increased more strongly when the pre-cue was invalid and required subjects to select the alternate response. Together, the results indicate that theta activity in dPMC and connected frontoparietal areas is involved in the re-adjustment of cue-induced action tendencies.

Bilateral Transcranial Direct Current Stimulation Reshapes Resting-State Brain Networks: A Magnetoencephalography Assessment

Transcranial direct current stimulation (tDCS) can noninvasively induce brain plasticity, and it is potentially useful to treat patients affected by neurological conditions. However, little is known about tDCS effects on resting-state brain networks, which are largely involved in brain physiological functions and in diseases. In this randomized, sham-controlled, double-blind study on healthy subjects, we have assessed the effect of bilateral tDCS applied over the sensorimotor cortices on brain and network activity using a whole-head magnetoencephalography system. Bilateral tDCS, with the cathode (−) centered over C4 and the anode (+) centered over C3, reshapes brain networks in a nonfocal fashion. Compared to sham stimulation, tDCS reduces left frontal alpha, beta, and gamma power and increases global connectivity, especially in delta, alpha, beta, and gamma frequencies. The increase of connectivity is consistent across bands and widespread. These results shed new light on the effects of tDCS and may be of help in personalizing treatments in neurological disorders.

Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy

Fusion of electroencephalography (EEG) and magnetoencephalography (MEG) data using maximum entropy on the mean method (MEM-fusion) takes advantage of the complementarities between EEG and MEG to improve localization accuracy. Simulation studies demonstrated MEM-fusion to be robust especially in noisy conditions such as single spike source localizations (SSSL). Our objective was to assess the reliability of SSSL using MEM-fusion on clinical data. We proposed to cluster SSSL results to find the most reliable and consistent source map from the reconstructed sources, the so-called consensus map. Thirty-four types of interictal epileptic discharges (IEDs) were analyzed from 26 patients with well-defined epileptogenic focus. SSSLs were performed on EEG, MEG, and fusion data and consensus maps were estimated using hierarchical clustering. Qualitative (spike-to-spike reproducibility rate, SSR) and quantitative (localization error and spatial dispersion) assessments were performed using the epileptogenic focus as clinical reference. Fusion SSSL provided significantly better results than EEG or MEG alone. Fusion found at least one cluster concordant with the clinical reference in all cases. This concordant cluster was always the one involving the highest number of spikes. Fusion yielded highest reproducibility (SSR EEG = 55%, MEG = 71%, fusion = 90%) and lowest localization error. Also, using only few channels from either modality (21EEG + 272MEG or 54EEG + 25MEG) was sufficient to reach accurate fusion. MEM-fusion with consensus map approach provides an objective way of finding the most reliable and concordant generators of IEDs. We, therefore, suggest the pertinence of SSSL using MEM-fusion as a valuable clinical tool for presurgical evaluation of epilepsy.