Hierarchical multiscale Bayesian algorithm for robust MEG/EEG source reconstruction

From bench to bedside: Overview of magnetoencephalography in basic principle, signal processing, source localization and clinical applications

Magnetoencephalography (MEG) is a non-invasive technique that can precisely capture the dynamic spatiotemporal patterns of the brain by measuring the magnetic fields arising from neuronal activity along the order of milliseconds. Observations of brain dynamics have been used in cognitive neuroscience, the diagnosis of neurological diseases, and the brain-computer interface (BCI). In this study, we outline the basic principle, signal processing, and source localization of MEG, and describe its clinical applications for cognitive assessment, the diagnoses of neurological diseases and mental disorders, preoperative evaluation, and the BCI. This review not only provides an overall perspective of MEG, ranging from practical techniques to clinical applications, but also enhances the prevalent understanding of neural mechanisms. The use of MEG is expected to lead to significant breakthroughs in neuroscience.

A tutorial on fitting joint models of M/EEG and behavior to understand cognition

We present motivation and practical steps necessary to find parameter estimates of joint models of behavior and neural electrophysiological data. This tutorial is written for researchers wishing to build joint models of human behavior and scalp and intracranial electroencephalographic (EEG) or magnetoencephalographic (MEG) data, and more specifically those researchers who seek to understand human cognition. Although these techniques could easily be applied to animal models, the focus of this tutorial is on human participants. Joint modeling of M/EEG and behavior requires some knowledge of existing computational and cognitive theories, M/EEG artifact correction, M/EEG analysis techniques, cognitive modeling, and programming for statistical modeling implementation. This paper seeks to give an introduction to these techniques as they apply to estimating parameters from neurocognitive models of M/EEG and human behavior, and to evaluate model results and compare models. Due to our research and knowledge on the subject matter, our examples in this paper will focus on testing specific hypotheses in human decision-making theory. However, most of the motivation and discussion of this paper applies across many modeling procedures and applications. We provide Python (and linked R) code examples in the tutorial and appendix. Readers are encouraged to try the exercises at the end of the document.

Electrophysiological brain imaging based on simulation-driven deep learning in the context of epilepsy

Identifying the location, the spatial extent and the electrical activity of distributed brain sources in the context of epilepsy through ElectroEncephaloGraphy (EEG) recordings is a challenging task because of the highly ill-posed nature of the underlying Electrophysiological Source Imaging (ESI) problem. To guarantee a unique solution, most existing ESI methods pay more attention to solve this inverse problem by imposing physiological constraints. This paper proposes an efficient ESI approach based on simulation-driven deep learning. Epileptic High-resolution 256-channels scalp EEG (Hr-EEG) signals are simulated in a realistic manner to train the proposed patient-specific model. More particularly, a computational neural mass model developed in our team is used to generate the temporal dynamics of the activity of each dipole while the forward problem is solved using a patient-specific three-shell realistic head model and the boundary element method. A Temporal Convolutional Network (TCN) is considered in the proposed model to capture local spatial patterns. To enable the model to observe the EEG signals from different scale levels, the multi-scale strategy is leveraged to capture the overall features and fine-grain features by adjusting the convolutional kernel size. Then, the Long Short-Term Memory (LSTM) is used to extract temporal dependencies among the computed spatial features. The performance of the proposed method is evaluated through three different scenarios of realistic synthetic interictal Hr-EEG data as well as on real interictal Hr-EEG data acquired in three patients with drug-resistant partial epilepsy, during their presurgical evaluation. A performance comparison study is also conducted with two other deep learning-based methods and four classical ESI techniques. The proposed model achieved a Dipole Localization Error (DLE) of 1.39 and Normalized Hamming Distance (NHD) of 0.28 in the case of one patch with SNR of 10 dB. In the case of two uncorrelated patches with an SNR of 10 dB, obtained DLE and NHD were respectively 1.50 and 0.28. Even in the more challenging scenario of two correlated patches with an SNR of 10 dB, the proposed approach still achieved a DLE of 3.74 and an NHD of 0.43. The results obtained on simulated data demonstrate that the proposed method outperforms the existing methods for different signal-to-noise and source configurations. The good behavior of the proposed method is also confirmed on real interictal EEG data. The robustness with respect to noise makes it a promising and alternative tool to localize epileptic brain areas and to reconstruct their electrical activities from EEG signals.

Sparse Bayesian Learning for End-to-End EEG Decoding

Decoding brain activity from non-invasive electroencephalography (EEG) is crucial for brain-computer interfaces (BCIs) and the study of brain disorders. Notably, end-to-end EEG decoding has gained widespread popularity in recent years owing to the remarkable advances in deep learning research. However, many EEG studies suffer from limited sample sizes, making it difficult for existing deep learning models to effectively generalize to highly noisy EEG data. To address this fundamental limitation, this paper proposes a novel end-to-end EEG decoding algorithm that utilizes a low-rank weight matrix to encode both spatio-temporal filters and the classifier, all optimized under a principled sparse Bayesian learning (SBL) framework. Importantly, this SBL framework also enables us to learn hyperparameters that optimally penalize the model in a Bayesian fashion. The proposed decoding algorithm is systematically benchmarked on five motor imagery BCI EEG datasets ( N=192) and an emotion recognition EEG dataset ( N=45), in comparison with several contemporary algorithms, including end-to-end deep-learning-based EEG decoding algorithms. The classification results demonstrate that our algorithm significantly outperforms the competing algorithms while yielding neurophysiologically meaningful spatio-temporal patterns. Our algorithm therefore advances the state-of-the-art by providing a novel EEG-tailored machine learning tool for decoding brain activity.

WRA-MTSI: A Robust Extended Source Imaging Algorithm Based on Multi-Trial EEG

OBJECTIVE

Reconstructing brain activities from electroencephalography (EEG) signals is crucial for studying brain functions and their abnormalities. However, since EEG signals are nonstationary and vulnerable to noise, brain activities reconstructed from single-trial EEG data are often unstable, and significant variability may occur across different EEG trials even for the same cognitive task.

METHODS

In an effort to leverage the shared information across the EEG data of multiple trials, this paper proposes a multi-trial EEG source imaging method based on Wasserstein regularization, termed WRA-MTSI. In WRA-MTSI, Wasserstein regularization is employed to perform multi-trial source distribution similarity learning, and the structured sparsity constraint is enforced to enable accurate estimation of the source extents, locations and time series. The resulting optimization problem is solved by a computationally efficient algorithm based on the alternating direction method of multipliers (ADMM).

RESULTS

Both numerical simulations and real EEG data analysis demonstrate that WRA-MTSI outperforms existing single-trial ESI methods (e.g., wMNE, LORETA, SISSY, and SBL) in mitigating the influence of artifacts in EEG data. Moreover, WRA-MTSI yields superior performance compared to other state-of-the-art multi-trial ESI methods (e.g., group lasso, the dirty model, and MTW) in estimating source extents.

CONCLUSION AND SIGNIFICANCE

WRA-MTSI may serve as an effective robust EEG source imaging method in the presence of multi-trial noisy EEG data.

One Definition to Join Them All: The N-Spherical Solution for the EEG Lead Field

Albeit its simplicity, the concentric spheres head model is widely used in EEG. The reason behind this is its simple mathematical definition, which allows for the calculation of lead fields with negligible computational cost, for example, for iterative approaches. Nevertheless, the literature shows contradictory formulations for the electrical solution of this head model. In this work, we study several different definitions for the electrical lead field of a four concentric spheres conduction model, finding that their results are contradictory. A thorough exploration of the mathematics used to build these formulations, provided in the original works, allowed for the identification of errors in some of the formulae, which proved to be the reason for the discrepancies. Moreover, this mathematical review revealed the iterative nature of some of these formulations, which allowed us to develop a formulation to solve the lead field in a head model built from an arbitrary number of concentric, homogeneous, and isotropic spheres.

Bayesian Adaptive Beamformer for Robust Electromagnetic Brain Imaging of Correlated Sources in High Spatial Resolution

Reconstructing complex brain source activity at a high spatiotemporal resolution from magnetoencephalography (MEG) or electroencephalography (EEG) remains a challenging problem. Adaptive beamformers are routinely deployed for this imaging domain using the sample data covariance. However adaptive beamformers have long been hindered by 1) high degree of correlation between multiple brain sources, and 2) interference and noise embedded in sensor measurements. This study develops a novel framework for minimum variance adaptive beamformers that uses a model data covariance learned from data using a sparse Bayesian learning algorithm (SBL-BF). The learned model data covariance effectively removes influence from correlated brain sources and is robust to noise and interference without the need for baseline measurements. A multiresolution framework for model data covariance computation and parallelization of the beamformer implementation enables efficient high-resolution reconstruction images. Results with both simulations and real datasets indicate that multiple highly correlated sources can be accurately reconstructed, and that interference and noise can be sufficiently suppressed. Reconstructions at 2-2.5mm resolution (  ∼  150K voxels) are possible with efficient run times of 1-3 minutes. This novel adaptive beamforming algorithm significantly outperforms the state-of-the-art benchmarks. Therefore, SBL-BF provides an effective framework for efficiently reconstructing multiple correlated brain sources with high resolution and robustness to interference and noise.

Bayesian Algorithms for Joint Estimation of Brain Activity and Noise in Electromagnetic Imaging

Simultaneously estimating brain source activity and noise has long been a challenging task in electromagnetic brain imaging using magneto- and electroencephalography. The problem is challenging not only in terms of solving the NP-hard inverse problem of reconstructing unknown brain activity across thousands of voxels from a limited number of sensors, but also for the need to simultaneously estimate the noise and interference. We present a generative model with an augmented leadfield matrix to simultaneously estimate brain source activity and sensor noise statistics in electromagnetic brain imaging (EBI). We then derive three Bayesian inference algorithms for this generative model (expectation-maximization (EBI-EM), convex bounding (EBI-Convex) and fixed-point (EBI-Mackay)) to simultaneously estimate the hyperparameters of the prior distribution for brain source activity and sensor noise. A comprehensive performance evaluation for these three algorithms is performed. Simulations consistently show that the performance of EBI-Convex and EBI-Mackay updates is superior to that of EBI-EM. In contrast to the EBI-EM algorithm, both EBI-Convex and EBI-Mackay updates are quite robust to initialization, and are computationally efficient with fast convergence in the presence of both Gaussian and real brain noise. We also demonstrate that EBI-Convex and EBI-Mackay update algorithms can reconstruct complex brain activity with only a few trials of sensor data, and for resting-state data, achieving significant improvement in source reconstruction and noise learning for electromagnetic brain imaging.

Empirical Bayesian localization of event-related time-frequency neural activity dynamics

Accurate reconstruction of the spatio-temporal dynamics of event-related cortical oscillations across human brain regions is an important problem in functional brain imaging and human cognitive neuroscience with magnetoencephalography (MEG) and electroencephalography (EEG). The problem is challenging not only in terms of localization of complex source configurations from sensor measurements with unknown noise and interference but also for reconstruction of transient event-related time-frequency dynamics of cortical oscillations. We recently proposed a robust empirical Bayesian algorithm for simultaneous reconstruction of complex brain source activity and noise covariance, in the context of evoked and resting-state data. In this paper, we expand upon this empirical Bayesian framework for optimal reconstruction of event-related time-frequency dynamics of regional cortical oscillations, referred to as time-frequency Champagne (TFC). This framework enables imaging of five-dimensional (space, time, and frequency) event-related brain activity from M/EEG data, and can be viewed as a time-frequency optimized adaptive Bayesian beamformer. We evaluate TFC in both simulations and several real datasets, with comparisons to benchmark standards - variants of time-frequency optimized adaptive beamformers (TFBF) as well as the sLORETA algorithm. In simulations, we demonstrate several advantages in estimating time-frequency cortical oscillatory dynamics compared to benchmarks. With real MEG data, we demonstrate across many datasets that the proposed approach is robust to highly correlated brain activity and low SNR data, and is able to accurately reconstruct cortical dynamics with data from just a few epochs.

NLGC: Network Localized Granger Causality with Application to MEG Directional Functional Connectivity Analysis

Identifying the directed connectivity that underlie networked activity between different cortical areas is critical for understanding the neural mechanisms behind sensory processing. Granger causality (GC) is widely used for this purpose in functional magnetic resonance imaging analysis, but there the temporal resolution is low, making it difficult to capture the millisecond-scale interactions underlying sensory processing. Magnetoencephalography (MEG) has millisecond resolution, but only provides low-dimensional sensor-level linear mixtures of neural sources, which makes GC inference challenging. Conventional methods proceed in two stages: First, cortical sources are estimated from MEG using a source localization technique, followed by GC inference among the estimated sources. However, the spatiotemporal biases in estimating sources propagate into the subsequent GC analysis stage, may result in both false alarms and missing true GC links. Here, we introduce the Network Localized Granger Causality (NLGC) inference paradigm, which models the source dynamics as latent sparse multivariate autoregressive processes and estimates their parameters directly from the MEG measurements, integrated with source localization, and employs the resulting parameter estimates to produce a precise statistical characterization of the detected GC links. We offer several theoretical and algorithmic innovations within NLGC and further examine its utility via comprehensive simulations and application to MEG data from an auditory task involving tone processing from both younger and older participants. Our simulation studies reveal that NLGC is markedly robust with respect to model mismatch, network size, and low signal-to-noise ratio, whereas the conventional two-stage methods result in high false alarms and mis-detections. We also demonstrate the advantages of NLGC in revealing the cortical network-level characterization of neural activity during tone processing and resting state by delineating task- and age-related connectivity changes.

Multiple Source Detection based on Spatial Clustering and Its Applications on Wearable OPM-MEG

OBJECTIVE

Magnetoencephalography (MEG) is a non-invasive technique that measures the magnetic fields of brain activity. In particular, a new type of optically pumped magnetometer (OPM)-based wearable MEG system has been developed in recent years. Source localization in MEG can provide signals and locations of brain activity. However, conventional source localization methods face the difficulty of accurately estimating multiple sources. The present study presented a new parametric method to estimate the number of sources and localize multiple sources. In addition, we applied the proposed method to a constructed wearable OPM-MEG system.

METHODS

We used spatial clustering of the dipole spatial distribution to detect sources. The spatial distribution of dipoles was obtained by segmenting the MEG data temporally into slices and then estimating the parameters of the dipoles on each data slice using the particle swarm optimization algorithm. Spatial clustering was performed using the spatial-temporal density-based spatial clustering of applications with a noise algorithm. The performance of our approach for detecting multiple sources was compared with that of four typical benchmark algorithms using the OPM-MEG sensor configuration.

RESULTS

The simulation results showed that the proposed method had the best performance for detecting multiple sources. Moreover, the effectiveness of the method was verified by a multimodel sensory stimuli experiment on a real constructed 31-channel OPM-MEG.

CONCLUSION

Our study provides an effective method for the detection of multiple sources.

SIGNIFICANCE

With the improvement of the source localization methods, MEG may have a wider range of applications in neuroscience and clinical research.

Hyperbolic trade-off: The importance of balancing trial and subject sample sizes in neuroimaging

Here we investigate the crucial role of trials in task-based neuroimaging from the perspectives of statistical efficiency and condition-level generalizability. Big data initiatives have gained popularity for leveraging a large sample of subjects to study a wide range of effect magnitudes in the brain. On the other hand, most task-based FMRI designs feature a relatively small number of subjects, so that resulting parameter estimates may be associated with compromised precision. Nevertheless, little attention has been given to another important dimension of experimental design, which can equally boost a study's statistical efficiency: the trial sample size. The common practice of condition-level modeling implicitly assumes no cross-trial variability. Here, we systematically explore the different factors that impact effect uncertainty, drawing on evidence from hierarchical modeling, simulations and an FMRI dataset of 42 subjects who completed a large number of trials of cognitive control task. We find that, due to an approximately symmetic hyperbola-relationship between trial and subject sample sizes in the presence of relatively large cross-trial variability, 1) trial sample size has nearly the same impact as subject sample size on statistical efficiency; 2) increasing both the number of trials and subjects improves statistical efficiency more effectively than focusing on subjects alone; 3) trial sample size can be leveraged alongside subject sample size to improve the cost-effectiveness of an experimental design; 4) for small trial sample sizes, trial-level modeling, rather than condition-level modeling through summary statistics, may be necessary to accurately assess the standard error of an effect estimate. We close by making practical suggestions for improving experimental designs across neuroimaging and behavioral studies.

Imaging the extent and location of spatiotemporally distributed epileptiform sources from MEG measurements

Non-invasive MEG/EEG source imaging provides valuable information about the epileptogenic brain areas which can be used to aid presurgical planning in focal epilepsy patients suffering from drug-resistant seizures. However, the source extent estimation for electrophysiological source imaging remains to be a challenge and is usually largely dependent on subjective choice. Our recently developed algorithm, fast spatiotemporal iteratively reweighted edge sparsity minimization (FAST-IRES) strategy, has been shown to objectively estimate extended sources from EEG recording, while it has not been applied to MEG recordings. In this work, through extensive numerical experiments and real data analysis in a group of focal drug-resistant epilepsy patients' interictal spikes, we demonstrated the ability of FAST-IRES algorithm to image the location and extent of underlying epilepsy sources from MEG measurements. Our results indicate the merits of FAST-IRES in imaging the location and extent of epilepsy sources for pre-surgical evaluation from MEG measurements.

Unification of sparse Bayesian learning algorithms for electromagnetic brain imaging with the majorization minimization framework

Methods for electro- or magnetoencephalography (EEG/MEG) based brain source imaging (BSI) using sparse Bayesian learning (SBL) have been demonstrated to achieve excellent performance in situations with low numbers of distinct active sources, such as event-related designs. This paper extends the theory and practice of SBL in three important ways. First, we reformulate three existing SBL algorithms under the majorization-minimization (MM) framework. This unification perspective not only provides a useful theoretical framework for comparing different algorithms in terms of their convergence behavior, but also provides a principled recipe for constructing novel algorithms with specific properties by designing appropriate bounds of the Bayesian marginal likelihood function. Second, building on the MM principle, we propose a novel method called LowSNR-BSI that achieves favorable source reconstruction performance in low signal-to-noise-ratio (SNR) settings. Third, precise knowledge of the noise level is a crucial requirement for accurate source reconstruction. Here we present a novel principled technique to accurately learn the noise variance from the data either jointly within the source reconstruction procedure or using one of two proposed cross-validation strategies. Empirically, we could show that the monotonous convergence behavior predicted from MM theory is confirmed in numerical experiments. Using simulations, we further demonstrate the advantage of LowSNR-BSI over conventional SBL in low-SNR regimes, and the advantage of learned noise levels over estimates derived from baseline data. To demonstrate the usefulness of our novel approach, we show neurophysiologically plausible source reconstructions on averaged auditory evoked potential data.

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.

Clinical Validation of the Champagne Algorithm for Epilepsy Spike Localization

Magnetoencephalography (MEG) is increasingly used for presurgical planning in people with medically refractory focal epilepsy. Localization of interictal epileptiform activity, a surrogate for the seizure onset zone whose removal may prevent seizures, is challenging and depends on the use of multiple complementary techniques. Accurate and reliable localization of epileptiform activity from spontaneous MEG data has been an elusive goal. One approach toward this goal is to use a novel Bayesian inference algorithm-the Champagne algorithm with noise learning-which has shown tremendous success in source reconstruction, especially for focal brain sources. In this study, we localized sources of manually identified MEG spikes using the Champagne algorithm in a cohort of 16 patients with medically refractory epilepsy collected in two consecutive series. To evaluate the reliability of this approach, we compared the performance to equivalent current dipole (ECD) modeling, a conventional source localization technique that is commonly used in clinical practice. Results suggest that Champagne may be a robust, automated, alternative to manual parametric dipole fitting methods for localization of interictal MEG spikes, in addition to its previously described clinical and research applications.

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.

Probabilistic Structure Learning for EEG/MEG Source Imaging With Hierarchical Graph Priors

Brain source imaging is an important method for noninvasively characterizing brain activity using Electroencephalogram (EEG) or Magnetoencephalography (MEG) recordings. Traditional EEG/MEG Source Imaging (ESI) methods usually assume the source activities at different time points are unrelated, and do not utilize the temporal structure in the source activation, making the ESI analysis sensitive to noise. Some methods may encourage very similar activation patterns across the entire time course and may be incapable of accounting the variation along the time course. To effectively deal with noise while maintaining flexibility and continuity among brain activation patterns, we propose a novel probabilistic ESI model based on a hierarchical graph prior. Under our method, a spanning tree constraint ensures that activity patterns have spatiotemporal continuity. An efficient algorithm based on an alternating convex search is presented to solve the resulting problem of the proposed model with guaranteed convergence. Comprehensive numerical studies using synthetic data on a realistic brain model are conducted under different levels of signal-to-noise ratio (SNR) from both sensor and source spaces. We also examine the EEG/MEG datasets in two real applications, in which our ESI reconstructions are neurologically plausible. All the results demonstrate significant improvements of the proposed method over benchmark methods in terms of source localization performance, especially at high noise levels.

A Novel Bayesian Approach for EEG Source Localization

We propose a new method for EEG source localization. An efficient solution to this problem requires choosing an appropriate regularization term in order to constraint the original problem. In our work, we adopt the Bayesian framework to place constraints; hence, the regularization term is closely connected to the prior distribution. More specifically, we propose a new sparse prior for the localization of EEG sources. The proposed prior distribution has sparse properties favoring focal EEG sources. In order to obtain an efficient algorithm, we use the variational Bayesian (VB) framework which provides us with a tractable iterative algorithm of closed-form equations. Additionally, we provide extensions of our method in cases where we observe group structures and spatially extended EEG sources. We have performed experiments using synthetic EEG data and real EEG data from three publicly available datasets. The real EEG data are produced due to the presentation of auditory and visual stimulus. We compare the proposed method with well-known approaches of EEG source localization and the results have shown that our method presents state-of-the-art performance, especially in cases where we expect few activated brain regions. The proposed method can effectively detect EEG sources in various circumstances. Overall, the proposed sparse prior for EEG source localization results in more accurate localization of EEG sources than state-of-the-art approaches.

Robust estimation of noise for electromagnetic brain imaging with the champagne algorithm

Robust estimation of the number, location, and activity of multiple correlated brain sources has long been a challenging task in electromagnetic brain imaging from M/EEG data, one that is significantly impacted by interference from spontaneous brain activity, sensor noise, and other sources of artifacts. Recently, we introduced the Champagne algorithm, a novel Bayesian inference algorithm that has shown tremendous success in M/EEG source reconstruction. Inherent to Champagne and most other related Bayesian reconstruction algorithms is the assumption that the noise covariance in sensor data can be estimated from “baseline” or “control” measurements. However, in many scenarios, such baseline data is not available, or is unreliable, and it is unclear how best to estimate the noise covariance. In this technical note, we propose several robust methods to estimate the contributions to sensors from noise arising from outside the brain without the need for additional baseline measurements. The incorporation of these methods for diagonal noise covariance estimation improves the robust reconstruction of complex brain source activity under high levels of noise and interference, while maintaining the performance features of Champagne. Specifically, we show that the resulting algorithm, Champagne with noise learning, is quite robust to initialization and is computationally efficient. In simulations, performance of the proposed noise learning algorithm is consistently superior to Champagne without noise learning. We also demonstrate that, even without the use of any baseline data, Champagne with noise learning is able to reconstruct complex brain activity with just a few trials or even a single trial, demonstrating significant improvements in source reconstruction for electromagnetic brain imaging.

NUTMEG: Open Source Software for M/EEG Source Reconstruction

Neurodynamic Utility Toolbox for Magnetoencephalo- and Electroencephalography (NUTMEG) is an open-source MATLAB-based toolbox for the analysis and reconstruction of magnetoencephalography/electroencephalography data in source space. NUTMEG includes a variety of options for the user in data import, preprocessing, source reconstruction, and functional connectivity. A group analysis toolbox allows the user to run a variety of inferential statistics on their data in an easy-to-use GUI-driven format. Importantly, NUTMEG features an interactive five-dimensional data visualization platform. A key feature of NUTMEG is the availability of a large menu of interference cancelation and source reconstruction algorithms. Each NUTMEG operation acts as a stand-alone MATLAB function, allowing the package to be easily adaptable and scripted for the more advanced user for interoperability with other software toolboxes. Therefore, NUTMEG enables a wide range of users access to a complete "sensor-to- source-statistics" analysis pipeline.

Robust Empirical Bayesian Reconstruction of Distributed Sources for Electromagnetic Brain Imaging

Electromagnetic brain imaging is the reconstruction of brain activity from non-invasive recordings of the magnetic fields and electric potentials. An enduring challenge in this imaging modality is estimating the number, location, and time course of sources, especially for the reconstruction of distributed brain sources with complex spatial extent. Here, we introduce a novel robust empirical Bayesian algorithm that enables better reconstruction of distributed brain source activity with two key ideas: kernel smoothing and hyperparameter tiling. Since the proposed algorithm builds upon many of the performance features of the sparse source reconstruction algorithm - Champagne and we refer to this algorithm as Smooth Champagne. Smooth Champagne is robust to the effects of high levels of noise, interference, and highly correlated brain source activity. Simulations demonstrate excellent performance of Smooth Champagne when compared to benchmark algorithms in accurately determining the spatial extent of distributed source activity. Smooth Champagne also accurately reconstructs real MEG and EEG data.

Comparison of DSSP and tSSS algorithms for removing artifacts from vagus nerve stimulators in magnetoencephalography data

OBJECTIVE

Large-amplitude artifacts from vagus nerve stimulator (VNS) implants for refractory epilepsy affect magnetoencephalography (MEG) recordings and are difficult to reject, resulting in unusable data from this important population of patients who are frequently evaluated for surgical treatment of epilepsy. Here we compare the performance of two artifact removal algorithms for MEG data: dual signal subspace projection (DSSP) and temporally extended signal space separation (tSSS).

APPROACH

Each algorithm's performance was first evaluated in simulations. We then tested the performance of each algorithm on resting-state MEG data from patients with VNS implants. We also examined how each algorithm improved source localization of somatosensory evoked fields in patients with VNS implants.

MAIN RESULTS

DSSP and tSSS algorithms have a similar ability to reject interference in both simulated and real MEG data if the origin location for tSSS is appropriately set. If the origin set for tSSS is inappropriate, the signal after tSSS can be distorted due to a mismatch between the internal region and the actual source space. Both DSSP and tSSS are able to remove large-amplitude artifacts from outside the brain. DSSP might be a better choice than tSSS when the choice of origin location for tSSS is difficult.

SIGNIFICANCE

Both DSSP and tSSS algorithms can recover distorted MEG recordings from people with intractable epilepsy and VNS implants, improving epileptic spike identification and source localization of both functional activity and epileptiform activity.

Exploring EEG Effective Connectivity Network in Estimating Influence of Color on Emotion and Memory

Color is a perceptual stimulus that has a significant impact on improving human emotion and memory. Studies have revealed that colored multimedia learning materials (MLMs) have a positive effect on learner's emotion and learning where it was assessed by subjective/objective measurements. This study aimed to quantitatively assess the influence of colored MLMs on emotion, cognitive processes during learning, and long-term memory (LTM) retention using electroencephalography (EEG). The dataset consisted of 45 healthy participants, and MLMs were designed in colored or achromatic illustrations to elicit emotion and that to assess its impact on LTM retention after 30-min and 1-month delay. The EEG signal analysis was first started to estimate the effective connectivity network (ECN) using the phase slope index and expand it to characterize the ECN pattern using graph theoretical analysis. EEG results showed that colored MLMs had influences on theta and alpha networks, including (1) an increased frontal-parietal connectivity (top-down processing), (2) a larger number of brain hubs, (3) a lower clustering coefficient, and (4) a higher local efficiency, indicating that color influences information processing in the brain, as reflected by ECN, together with a significant improvement in learner's emotion and memory performance. This is evidenced by a more positive emotional valence and higher recall accuracy for groups who learned with colored MLMs than that of achromatic MLMs. In conclusion, this paper demonstrated how the EEG ECN parameters could help quantify the influences of colored MLMs on emotion and cognitive processes during learning.

A Novel Scanning Algorithm for MEG/EEG imaging using Covariance Partitioning and Noise Learning

In this paper, we present a novel scanning algorithm, called Covariance Optimization Garnering Noise for Active Cancellation (COGNAC), for magnetoencephalography (MEG) and electroencephalography (EEG) source localization. COGNAC uses a probabilistic graphical generative model for describing sensor data. This novel generative model partitions contributions to sensor data from sources at a particular scan location and from sources outside the scan location, with corresponding multi-resolution variance parameters that are estimated from data. Maximizing a convex upper bound on the marginal likelihood of the data under this generative model results in a cost function that can be optimized efficiently. Importantly, this generative model enables learning of sensor noise without the need for additional baseline or pre-stimulus data. The resulting inference algorithm is quite robust to reconstruction of highly correlated sources and to the effect of high levels of interference and noise sources. Algorithm performance was compared to representative benchmark algorithms on both simulated and real brain activity. In simulations, performance of our novel algorithm is consistently superior to benchmarks. We also demonstrate that the new algorithm is robust to correlated brain activity present in real MEG/EEG data.