EEG/MEG source imaging of deep brain activity within the maximum entropy on the mean framework: Simulations and validation in epilepsy.

Electro/Magneto-EncephaloGraphy (EEG/MEG) source imaging (EMSI) of epileptic activity from deep generators is often challenging due to the higher sensitivity of EEG/MEG to superficial regions and to the spatial configuration of subcortical structures. We previously demonstrated the ability of the coherent Maximum Entropy on the Mean (cMEM) method to accurately localize the superficial cortical generators and their spatial extent. Here, we propose a depth-weighted adaptation of cMEM to localize deep generators more accurately. These methods were evaluated using realistic MEG/high-density EEG (HD-EEG) simulations of epileptic activity and actual MEG/HD-EEG recordings from patients with focal epilepsy. We incorporated depth-weighting within the MEM framework to compensate for its preference for superficial generators. We also included a mesh of both hippocampi, as an additional deep structure in the source model. We generated 5400 realistic simulations of interictal epileptic discharges for MEG and HD-EEG involving a wide range of spatial extents and signal-to-noise ratio (SNR) levels, before investigating EMSI on clinical HD-EEG in 16 patients and MEG in 14 patients. Clinical interictal epileptic discharges were marked by visual inspection. We applied three EMSI methods: cMEM, depth-weighted cMEM and depth-weighted minimum norm estimate (MNE). The ground truth was defined as the true simulated generator or as a drawn region based on clinical information available for patients. For deep sources, depth-weighted cMEM improved the localization when compared to cMEM and depth-weighted MNE, whereas depth-weighted cMEM did not deteriorate localization accuracy for superficial regions. For patients' data, we observed improvement in localization for deep sources, especially for the patients with mesial temporal epilepsy, for which cMEM failed to reconstruct the initial generator in the hippocampus. Depth weighting was more crucial for MEG (gradiometers) than for HD-EEG. Similar findings were found when considering depth weighting for the wavelet extension of MEM. In conclusion, depth-weighted cMEM improved the localization of deep sources without or with minimal deterioration of the localization of the superficial sources. This was demonstrated using extensive simulations with MEG and HD-EEG and clinical MEG and HD-EEG for 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.

Bayesian workflow for the investigation of hierarchical classification models from tau-PET and structural MRI data across the Alzheimer's disease spectrum.

Background: Alzheimer's disease (AD) diagnosis in its early stages remains difficult with current diagnostic approaches. Though tau neurofibrillary tangles (NFTs) generally follow the stereotypical pattern described by the Braak staging scheme, the network degeneration hypothesis (NDH) has suggested that NFTs spread selectively along functional networks of the brain. To evaluate this, we implemented a Bayesian workflow to develop hierarchical multinomial logistic regression models with increasing levels of complexity of the brain from tau-PET and structural MRI data to investigate whether it is beneficial to incorporate network-level information into an ROI-based predictive model for the presence/absence of AD. Methods: This study included data from the Translational Biomarkers in Aging and Dementia (TRIAD) longitudinal cohort from McGill University's Research Centre for Studies in Aging (MCSA). Baseline and 1 year follow-up structural MRI and [18F]MK-6240 tau-PET scans were acquired for 72 cognitive normal (CN), 23 mild cognitive impairment (MCI), and 18 Alzheimer's disease dementia subjects. We constructed the four following hierarchical Bayesian models in order of increasing complexity: (Model 1) a complete-pooling model with observations, (Model 2) a partial-pooling model with observations clustered within ROIs, (Model 3) a partial-pooling model with observations clustered within functional networks, and (Model 4) a partial-pooling model with observations clustered within ROIs that are also clustered within functional brain networks. We then investigated which of the models had better predictive performance given tau-PET or structural MRI data as an input, in the form of a relative annualized rate of change. Results: The Bayesian leave-one-out cross-validation (LOO-CV) estimate of the expected log pointwise predictive density (ELPD) results indicated that models 3 and 4 were substantially better than other models for both tau-PET and structural MRI inputs. For tau-PET data, model 3 was slightly better than 4 with an absolute difference in ELPD of 3.10 ± 1.30. For structural MRI data, model 4 was considerably better than other models with an absolute difference in ELPD of 29.83 ± 7.55 relative to model 3, the second-best model. Conclusion: Our results suggest that representing the data generating process in terms of a hierarchical model that encompasses both ROI-level and network-level heterogeneity leads to better predictive ability for both tau-PET and structural MRI inputs over all other model iterations.

Estimation of fMRI responses related to epileptic discharges using Bayesian hierarchical modeling.

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a unique and noninvasive method for epilepsy presurgical evaluation. When selecting voxels by null-hypothesis tests, the conventional analysis may overestimate fMRI response amplitudes related to interictal epileptic discharges (IEDs), especially when IEDs are rare. We aimed to estimate fMRI response amplitudes represented by blood oxygen level dependent (BOLD) percentage changes related to IEDs using a hierarchical model. It involves the local and distributed hemodynamic response homogeneity to regularize estimations. Bayesian inference was applied to fit the model. Eighty-two epilepsy patients who underwent EEG-fMRI and subsequent surgery were included in this study. A conventional voxel-wise general linear model was compared to the hierarchical model on estimated fMRI response amplitudes and on the concordance between the highest response cluster and the surgical cavity. The voxel-wise model overestimated fMRI responses compared to the hierarchical model, evidenced by a practically and statistically significant difference between the estimated BOLD percentage changes. Only the hierarchical model differentiated brief and long-lasting IEDs with significantly different BOLD percentage changes. Overall, the hierarchical model outperformed the voxel-wise model on presurgical evaluation, measured by higher prediction performance. When compared with a previous study, the hierarchical model showed higher performance metric values, but the same or lower sensitivity. Our results demonstrated the capability of the hierarchical model of providing more physiologically reasonable and more accurate estimations of fMRI response amplitudes induced by IEDs. To enhance the sensitivity of EEG-fMRI for presurgical evaluation, it may be necessary to incorporate more appropriate spatial priors and bespoke decision strategies.

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.

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.

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[Formula: see text] 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.

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.

The movement time analyser task investigated with functional near infrared spectroscopy: an ecologic approach for measuring hemodynamic response in the motor system.

AIMS: Movement time analyzer (MTA) is an objective instrument to evaluate the degree of motor impairment as well as to investigate the dopaminergic drug effect in Parkinson's disease patients. The aim of this study is to validate a new ecologic neuroimaging tool for quantifying MTA-related hemodynamic response of the cortical motor system by means of functional near-infrared spectroscopy (fNIRS). MATERIALS: 11 right-handed healthy volunteers (six male and five female, age range 27-64 years) were studied with fNIRS and functional magnetic resonance imaging (fMRI) while performing MTA task for each hand. RESULTS: MTA performance was better for the dominant hand and younger participants. Both fNIRS and fMRI analyses revealed MTA-related increase of haemoglobin levels in the primary motor and premotor cortices contralateral to the moving hand. This response progressively increased with aging. CONCLUSION: These findings supported the translation of fNIRS-based MTA behavioural tool in clinical practice.