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

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.

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

OBJECTIVE

The use of electrical source imaging (ESI) in assessing the source of interictal epileptic discharges (IEDs) is gaining increasing popularity in presurgical work-up of patients with drug-resistant focal epilepsy. While vigilance affects the ability to locate IEDs and identify the epileptogenic zone, we know little about its impact on ESI.

METHODS

We studied overnight high-density electroencephalography recordings in focal drug-resistant epilepsy. IEDs were marked visually in each vigilance state, and examined in the sensor and source space. ESIs were calculated and compared between all vigilance states and the clinical ground truth. Two conditions were considered within each vigilance state, an unequalized and an equalized number of IEDs.

RESULTS

The number, amplitude, and duration of IEDs were affected by the vigilance state, with N3 sleep presenting the highest number, amplitude, and duration for both conditions (P < 0.001), while signal-to-noise ratio only differed in the unequalized condition (P < 0.001). The vigilance state did not affect channel involvement (P > 0.05). ESI maps showed no differences in distance, quality, extent, or maxima distances compared to the clinical ground truth for both conditions (P > 0.05). Only when an absolute reference (wakefulness) was used, the channel involvement (P < 0.05) and ESI source extent (P < 0.01) were impacted during rapid-eye-movement (REM) sleep. Clustering of amplitude-sensitive and -insensitive ESI maps pointed to amplitude rather than the spatial profile as the driver (P < 0.05).

INTERPRETATION

IED ESI results are stable across vigilance states, including REM sleep, if controlled for amplitude and IED number. ESI is thus stable and invariant to the vigilance state.

Targeted density electrode placement achieves high concordance with traditional high-density EEG for electrical source imaging in epilepsy

OBJECTIVE

High-density (HD) electroencephalography (EEG) is increasingly used in presurgical epilepsy evaluation, but it is demanding in time and resources. To overcome these issues, we compared EEG source imaging (ESI) solutions with a targeted density and HD-EEG montage.

METHODS

HD-EEGs from patients undergoing presurgical evaluation were analyzed. A low-density recording was created by selecting the 25 electrodes of a standard montage from the 83 electrodes of the HD-EEG and adding 8-11 electrodes around the electrode with the highest amplitude interictal epileptiform discharges. The ESI solution from this "targeted" montage was compared to that from the HD-EEG using the distance between peak vertices, sublobar concordance and a qualitative similarity measure.

RESULTS

Fifty-eight foci of forty-three patients were included. The median distance between the peak vertices of the two montages was 13.2 mm, irrespective of focus' location. Tangential generators (n = 5/58) showed a higher distance than radial generators (p = 0.04). We found sublobar concordance in 54/58 of the foci (93%). Map similarity, assessed by an epileptologist, had a median score of 4/5.

CONCLUSIONS

ESI solutions obtained from a targeted density montage show high concordance with those calculated from HD-EEG.

SIGNIFICANCE

Requiring significantly fewer electrodes, targeted density EEG allows obtaining similar ESI solutions as traditional HD-EEG montage.

Feasibility of high-density electric source imaging in the presurgical workflow: Effect of number of spikes and automated spike detection

OBJECTIVE

Presurgical high-density electric source imaging (hdESI) of interictal epileptic discharges (IEDs) is only used by few epilepsy centers. One obstacle is the time-consuming workflow both for recording as well as for visual review. Therefore, we analyzed the effect of (a) an automated IED detection and (b) the number of IEDs on the accuracy of hdESI and time-effectiveness.

METHODS

In 22 patients with pharmacoresistant focal epilepsy receiving epilepsy surgery (Engel 1) we retrospectively detected IEDs both visually and semi-automatically using the EEG analysis software Persyst in 256-channel EEGs. The amount of IEDs, the Euclidean distance between hdESI maximum and resection zone, and the operator time were compared. Additionally, we evaluated the intra-individual effect of IED quantity on the distance between hdESI maximum of all IEDs and hdESI maximum when only a reduced amount of IEDs were included.

RESULTS

There was no significant difference in the number of IEDs between visually versus semi-automatically marked IEDs (74 ± 56 IEDs/patient vs 116 ± 115 IEDs/patient). The detection method of the IEDs had no significant effect on the mean distances between resection zone and hdESI maximum (visual: 26.07 ± 31.12 mm vs semi-automated: 33.6 ± 34.75 mm). However, the mean time needed to review the full datasets semi-automatically was shorter by 275 ± 46 min (305 ± 72 min vs 30 ± 26 min, P < 0.001). The distance between hdESI of the full versus reduced amount of IEDs of the same patient was smaller than 1 cm when at least a mean of 33 IEDs were analyzed. There was a significantly shorter intraindividual distance between resection zone and hdESI maximum when 30 IEDs were analyzed as compared to the analysis of only 10 IEDs (P < 0.001).

SIGNIFICANCE

Semi-automatized processing and limiting the amount of IEDs analyzed (~30-40 IEDs per cluster) appear to be time-saving clinical tools to increase the practicability of hdESI in the presurgical work-up.

Integration of interictal EEG source localization in presurgical epilepsy evaluation - A single-center prospective study

OBJECTIVE

To investigate cost in working hours for initial integration of interictal EEG source localization (ESL) into clinical practice of a tertiary epilepsy center, and to examine concordance of results obtained with three different ESL pipelines.

METHODS

This prospective study covered the first year of using ESL in the Epilepsy-Center Berlin-Brandenburg. Patients aged ≥14 years with drug-resistant focal epilepsy referred for noninvasive presurgical evaluation were included. Interictal ESL was based on low-density EEG and individual head models. Source maxima were obtained from two freely available software packages and one commercial provider. One physician and computer scientist documented their working hours for setting up and processing ESL. Additionally, a survey was conducted among epilepsy centers in Germany to assess the current role of ESL in presurgical evaluation.

RESULTS

Of 40 patients included, 22 (55%) had enough interictal spikes for ESL. The physician's working times decreased from median 4.7 hours [interquartile range 3.9-6.4] in the first third of cases to 2.0 hours [1.9-2.4] in the remaining two thirds; P < 0.01. In addition, computer scientist and physician spent a total of 35.5 and 33.0 working hours on setting up the digital infrastructure, and on training and testing. Sublobar agreement between all three pipelines was 20%, mean measurement of agreement (kappa) 0.13. Finally, the survey revealed that 53% of epilepsy centers in Germany currently use ESL for presurgical evaluation.

SIGNIFICANCE

This study provides information regarding expected effort and costs for integration of ESL into an epilepsy surgery program. Low result agreement across different ESL pipelines calls for further standardization.

Auditory white noise exposure results in intrinsic cortical excitability changes

Cortical excitability is commonly measured by applying magnetic stimulation in combination with measuring behavioral response. This measure has, however, some shortcomings including spatial limitation to the primary motor cortex and not accounting for intrinsic excitability fluctuations. Here, we use a measure for intrinsic excitability based on phase synchronization previously validated for epilepsy. We apply this measure in 30 healthy participants' magnetoencephalography (MEG) recordings during the exposure of auditory white noise, a stimulus that has been suggested to modify cortical excitability. Using cortical parcellation of the MEG source data, we could find a specific pattern of increased and decreased excitability while participants are exposed to white noise vs. silence. Specifically, excitability during white noise exposure decreases in the frontal lobe and increases in the temporal lobe. This study thus adds to the understanding of cortical excitability changes due to specific environmental stimuli as well as the spatial extent of these effects.

An in-vivo validation of ESI methods with focal sources

Electrophysiological source imaging (ESI) aims at reconstructing the precise origin of brain activity from measurements of the electric field on the scalp. Across laboratories/research centers/hospitals, ESI is performed with different methods, partly due to the ill-posedness of the underlying mathematical problem. However, it is difficult to find systematic comparisons involving a wide variety of methods. Further, existing comparisons rarely take into account the variability of the results with respect to the input parameters. Finally, comparisons are typically performed using either synthetic data, or in-vivo data where the ground-truth is only roughly known. We use an in-vivo high-density EEG dataset recorded during intracranial single pulse electrical stimulation, in which the true sources are substantially dipolar and their locations are precisely known. We compare ten different ESI methods, using their implementation in the MNE-Python package: MNE, dSPM, LORETA, sLORETA, eLORETA, LCMV beamformers, irMxNE, Gamma Map, SESAME and dipole fitting. We perform comparisons under multiple choices of input parameters, to assess the accuracy of the best reconstruction, as well as the impact of such parameters on the localization performance. Best reconstructions often fall within 1 cm from the true source, with most accurate methods hitting an average localization error of 1.2 cm and outperforming least accurate ones erring by 2.5 cm. As expected, dipolar and sparsity-promoting methods tend to outperform distributed methods. For several distributed methods, the best regularization parameter turned out to be the one in principle associated with low SNR, despite the high SNR of the available dataset. Depth weighting played no role for two out of the six methods implementing it. Sensitivity to input parameters varied widely between methods. While one would expect high variability being associated with low localization error at the best solution, this is not always the case, with some methods producing highly variable results and high localization error, and other methods producing stable results with low localization error. In particular, recent dipolar and sparsity-promoting methods provide significantly better results than older distributed methods. As we repeated the tests with "conventional" (32 channels) and dense (64, 128, 256 channels) EEG recordings, we observed little impact of the number of channels on localization accuracy; however, for distributed methods denser montages provide smaller spatial dispersion. Overall findings confirm that EEG is a reliable technique for localization of point sources and therefore reinforce the importance that ESI may have in the clinical context, especially when applied to identify the surgical target in potential candidates for epilepsy surgery.

Validating MEG source imaging of resting state oscillatory patterns with an intracranial EEG atlas

BACKGROUND

Magnetoencephalography (MEG) is a widely used non-invasive tool to estimate brain activity with high temporal resolution. However, due to the ill-posed nature of the MEG source imaging (MSI) problem, the ability of MSI to identify accurately underlying brain sources along the cortical surface is still uncertain and requires validation.

METHOD

We validated the ability of MSI to estimate the background resting state activity of 45 healthy participants by comparing it to the intracranial EEG (iEEG) atlas (https://mni-open-ieegatlas.

RESEARCH

mcgill.ca/). First, we applied wavelet-based Maximum Entropy on the Mean (wMEM) as an MSI technique. Next, we converted MEG source maps into intracranial space by applying a forward model to the MEG-reconstructed source maps, and estimated virtual iEEG (ViEEG) potentials on each iEEG channel location; we finally quantitatively compared those with actual iEEG signals from the atlas for 38 regions of interest in the canonical frequency bands.

RESULTS

The MEG spectra were more accurately estimated in the lateral regions compared to the medial regions. The regions with higher amplitude in the ViEEG than in the iEEG were more accurately recovered. In the deep regions, MEG-estimated amplitudes were largely underestimated and the spectra were poorly recovered. Overall, our wMEM results were similar to those obtained with minimum norm or beamformer source localization. Moreover, the MEG largely overestimated oscillatory peaks in the alpha band, especially in the anterior and deep regions. This is possibly due to higher phase synchronization of alpha oscillations over extended regions, exceeding the spatial sensitivity of iEEG but detected by MEG. Importantly, we found that MEG-estimated spectra were more comparable to spectra from the iEEG atlas after the aperiodic components were removed.

CONCLUSION

This study identifies brain regions and frequencies for which MEG source analysis is likely to be reliable, a promising step towards resolving the uncertainty in recovering intracerebral activity from non-invasive MEG studies.

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

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

Neurophysiology, Neuropsychology, and Epilepsy, in 2022: Hills We Have Climbed and Hills Ahead. Neurophysiology in epilepsy

Since the discovery of the human electroencephalogram (EEG), neurophysiology techniques have become indispensable tools in our armamentarium to localize epileptic seizures. New signal analysis techniques and the prospects of artificial intelligence and big data will offer unprecedented opportunities to further advance the field in the near future, ultimately resulting in improved quality of life for many patients with drug-resistant epilepsy. This article summarizes selected presentations from Day 1 of the two-day symposium "Neurophysiology, Neuropsychology, Epilepsy, 2022: Hills We Have Climbed and the Hills Ahead". Day 1 was dedicated to highlighting and honoring the work of Dr. Jean Gotman, a pioneer in EEG, intracranial EEG, simultaneous EEG/ functional magnetic resonance imaging, and signal analysis of epilepsy. The program focused on two main research directions of Dr. Gotman, and was dedicated to "High-frequency oscillations, a new biomarker of epilepsy" and "Probing the epileptic focus from inside and outside". All talks were presented by colleagues and former trainees of Dr. Gotman. The extended summaries provide an overview of historical and current work in the neurophysiology of epilepsy with emphasis on novel EEG biomarkers of epilepsy and source imaging and concluded with an outlook on the future of epilepsy research, and what is needed to bring the field to the next level.

Hierarchical Bayesian modeling of the relationship between task-related hemodynamic responses and cortical excitability

Investigating the relationship between task-related hemodynamic responses and cortical excitability is challenging because it requires simultaneous measurement of hemodynamic responses while applying noninvasive brain stimulation. Moreover, cortical excitability and task-related hemodynamic responses are both associated with inter-/intra-subject variability. To reliably assess such a relationship, we applied hierarchical Bayesian modeling. This study involved 16 healthy subjects who underwent simultaneous Paired Associative Stimulation (PAS10, PAS25, Sham) while monitoring brain activity using functional Near-Infrared Spectroscopy (fNIRS), targeting the primary motor cortex (M1). Cortical excitability was measured by Motor Evoked Potentials (MEPs), and the motor task-related hemodynamic responses were measured using fNIRS 3D reconstructions. We constructed three models to investigate: (1) PAS effects on the M1 excitability, (2) PAS effects on fNIRS hemodynamic responses to a finger tapping task, and (3) the correlation between PAS effects on M1 excitability and PAS effects on task-related hemodynamic responses. Significant increase in cortical excitability was found following PAS25, whereas a small reduction of the cortical excitability was shown after PAS10 and a subtle increase occurred after sham. Both HbO and HbR absolute amplitudes increased after PAS25 and decreased after PAS10. The probability of the positive correlation between modulation of cortical excitability and hemodynamic activity was 0.77 for HbO and 0.79 for HbR. We demonstrated that PAS stimulation modulates task-related cortical hemodynamic responses in addition to M1 excitability. Moreover, the positive correlation between PAS modulations of excitability and hemodynamics brought insight into understanding the fundamental properties of cortical function and cortical excitability.

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.

Contributions of Magnetoencephalography to Understanding Mechanisms of Generalized Epilepsies: Blurring the Boundary Between Focal and Generalized Epilepsies?

According to the latest operational 2017 ILAE classification of epileptic seizures, the generalized epileptic seizure is still conceptualized as "originating at some point within and rapidly engaging, bilaterally distributed networks." In contrast, the focal epileptic seizure is defined as "originating within networks limited to one hemisphere." Hence, one of the main concepts of "generalized" and "focal" epilepsy comes from EEG descriptions before the era of source localization, and a presumed simultaneous bilateral onset and bi-synchrony of epileptiform discharges remains a hallmark for generalized seizures. Current literature on the pathophysiology of generalized epilepsy supports the concept of a cortical epileptogenic focus triggering rapidly generalized epileptic discharges involving intact corticothalamic and corticocortical networks, known as the cortical focus theory. Likewise, focal epilepsy with rich connectivity can give rise to generalized spike and wave discharges resulting from widespread bilateral synchronization. Therefore, making this key distinction between generalized and focal epilepsy may be challenging in some cases, and for the first time, a combined generalized and focal epilepsy is categorized in the 2017 ILAE classification. Nevertheless, treatment options, such as the choice of antiseizure medications or surgical treatment, are the reason behind the importance of accurate epilepsy classification. Over the past several decades, plentiful scientific research on the pathophysiology of generalized epilepsy has been conducted using non-invasive neuroimaging and postprocessing of the electromagnetic neural signal by measuring the spatiotemporal and interhemispheric latency of bi-synchronous or generalized epileptiform discharges as well as network analysis to identify diagnostic and prognostic biomarkers for accurate diagnosis of the two major types of epilepsy. Among all the advanced techniques, magnetoencephalography (MEG) and multiple other methods provide excellent temporal and spatial resolution, inherently suited to analyzing and visualizing the propagation of generalized EEG activities. This article aims to provide a comprehensive literature review of recent innovations in MEG methodology using source localization and network analysis techniques that contributed to the literature of idiopathic generalized epilepsy in terms of pathophysiology and clinical prognosis, thus further blurring the boundary between focal and generalized epilepsy.

Detection of interictal epileptiform discharges in an extended scalp EEG array and high-density EEG - A prospective multicenter study

OBJECTIVES

High counts of averaged interictal epileptiform discharges (IEDs) are key components of accurate interictal electric source imaging (ESI) in patients with focal epilepsy. Automated detections may be time-efficient, but they need to identify the correct IED types. Thus, we compared semiautomated and automated detection of IED types in long-term video EEG monitoring (LTM) using an extended scalp EEG array and short-term high-density EEG (hdEEG) with visual detection of IED types and the seizure onset zone (SOZ).

METHODS

We prospectively recruited consecutive patients from four epilepsy centers who underwent both LTM with 40 electrodes scalp EEG and short-term hdEEG with 256 electrodes. Only patients with a single circumscribed SOZ in LTM were included. In LTM and hdEEG, IED types were identified visually, semiautomatically and automatically. Concordances of semiautomated and automated detections in LTM and hdEEG as well as visual detections in hdEEG were compared against visually detected IED types and the SOZ in LTM.

RESULTS

Fifty-two out of 62 patients with LTM and hdEEG were included. The most frequent IED types per patient, detected semiautomatically and automatically in LTM and visually in hdEEG, were significantly concordant with the most frequently visually identified IED type in LTM and the SOZ. Semiautomated and automated detections of IED types in hdEEG were significantly concordant with visually identified IED types in LTM only when IED types with more than 50 detected single IEDs were selected. The threshold of 50 detected IED in hdEEG was reached in half of the patients. For all IED types per patient, agreement between visual and semiautomated detections in LTM was high.

SIGNIFICANCE

Semiautomated and automated detections of IED types in LTM show significant agreement with visually detected IED types and the SOZ. In short-term hdEEG, semiautomated detections of IED types are concordant with visually detected IED types and the SOZ in LTM if high IED counts were detected.

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.

Case Report: Laser Ablation Guided by State of the Art Source Imaging Ends an Adolescent's 16-Year Quest for Seizure Freedom

Epilepsy surgery is the most effective therapeutic approach for children with drug resistant epilepsy (DRE). Recent advances in neurosurgery, such as the Laser Interstitial Thermal Therapy (LITT), improved the safety and non-invasiveness of this method. Electric and magnetic source imaging (ESI/MSI) plays critical role in the delineation of the epileptogenic focus during the presurgical evaluation of children with DRE. Yet, they are currently underutilized even in tertiary epilepsy centers. Here, we present a case of an adolescent who suffered from DRE for 16 years and underwent surgery at Cook Children's Medical Center (CCMC). The patient was previously evaluated in a level 4 epilepsy center and treated with multiple antiseizure medications for several years. Presurgical evaluation at CCMC included long-term video electroencephalography (EEG), magnetoencephalography (MEG) with simultaneous conventional EEG (19 channels) and high-density EEG (256 channels) in two consecutive sessions, MRI, and fluorodeoxyglucose - positron emission tomography (FDG-PET). Video long-term EEG captured nine focal-onset clinical seizures with a maximal evolution over the right frontal/frontal midline areas. MRI was initially interpreted as non-lesional. FDG-PET revealed a small region of hypometabolism at the anterior right superior temporal gyrus. ESI and MSI performed with dipole clustering showed a tight cluster of dipoles in the right anterior insula. The patient underwent intracranial EEG which indicated the right anterior insular as seizure onset zone. Eventually LITT rendered the patient seizure free (Engel 1; 12 months after surgery). Retrospective analysis of ESI and MSI clustered dipoles found a mean distance of dipoles from the ablated volume ranging from 10 to 25 mm. Our findings highlight the importance of recent technological advances in the presurgical evaluation and surgical treatment of children with DRE, and the underutilization of epilepsy surgery in children with DRE.

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.

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.

Linear distributed inverse solutions for interictal EEG source localisation

OBJECTIVE

To compare the spatial accuracy of 6 linear distributed inverse solutions for EEG source localisation of interictal epileptic discharges: Minimum Norm, Weighted Minimum Norm, Low-Resolution Electromagnetic Tomography (LORETA), Local Autoregressive Average (LAURA), Standardised LORETA, and Exact LORETA.

METHODS

Spatial accuracy was assessed clinically by retrospectively comparing the maximum source of averaged interictal discharges to the resected brain area in 30 patients with successful epilepsy surgery, based on 204-channel EEG. Additionally, localisation errors of the inverse solutions were assessed in computer simulations, with different levels of noise added to the signal in both sensor space and source space.

RESULTS

In the clinical evaluations, the source maximum was located inside the resected brain area in 50-57% of patients when using LORETA or LAURA, while all other inverse solutions performed significantly worse (17-30%; corrected p < 0.01). In the simulation studies, when noise levels exceeded 10%, LORETA and LAURA had substantially smaller localisation errors than the other inverse solutions.

CONCLUSIONS

LORETA and LAURA provided the highest spatial accuracy both in clinical and simulated data, alongside with a comparably high robustness towards noise.

SIGNIFICANCE

Among the different linear inverse solution algorithms tested, LORETA and LAURA might be preferred for interictal EEG source localisation.

Electric Source Imaging on Intracranial EEG Localizes Spatiotemporal Propagation of Interictal Spikes in Children with Epilepsy

Interictal epileptiform discharges (IEDs) serve as sensitive but not specific biomarkers of epilepsy that can delineate the epileptogenic zone (EZ) in patients with drug resistant epilepsy (DRE) undergoing surgery. Intracranial EEG (icEEG) studies have shown that IEDs propagate in time across large areas of the brain. The onset of this propagation is regarded as a more specific biomarker of epilepsy than areas of spread. Yet, the limited spatial resolution of icEEG does not allow to identify the onset of this activity with high precision. Here, we propose a new method of mapping the spatiotemporal propagation of IEDs (and identify its onset) by using Electrical Source Imaging (ESI) on icEEG bypassing the spatial limitations of icEEG. We validated our method on icEEG recordings from 8 children with DRE who underwent surgery with good outcome (Engel score =1). On each icEEG channel, we detected IEDs and identified the propagation onset using an automated algorithm. We localized the propagation of IEDs with dynamic Statistical Parametric Mapping (dSPM) using a time-sliding window approach. We defined two brain regions: the ESI-onset and ESI-spread zone. We estimated the overlap of these regions with resection volume (in percentage), which served as the gold-standard of the EZ. We also estimated the mean distance of these regions from resection and clinically defined seizure onset zone (SOZ). We observed spatio-temporal propagation of IEDs in all patients across several channels (98 [85-102]) with a mean duration of 155 ms [96-186 ms]. A higher overlap with resection was seen for the ESI-onset zone compared to spread (73.3 % [ 47.4-100 %], 36.5 % [20.3-59.9 %], p = 0.008). The distance of the ESI-onset from resection was shorter compared to the ESI-spread zone (4.3 mm [3.4-5.5 mm], 7.4 mm [6.0-20.6 mm], p = 0.008) and the same trend was observed for the distance from the SOZ (11.9 mm [7.2-15.1 mm], 20.6 mm [15.4-27.2 mm], p = 0.02). These findings show that our method can map the spatiotemporal propagation of IEDs and de-lineate its onset, which is a reliable and focal biomarker of the EZ in children with DRE.Clinical Relevance - ESI on icEEG recordings of children with DRE can localize the spikes propagation phenomenon and help in the delineation of the EZ.

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.

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.

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.

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.

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.