MEG/EEG are non-invasive imaging techniques that record brain activity with high temporal resolution. However, estimation of brain source currents from surface recordings requires solving an ill-conditioned inverse problem. Converging lines of evidence in neuroscience, from neuronal network models to resting-state imaging and neurophysiology, suggest that cortical activation is a distributed spatiotemporal dynamic process, supported by both local and long-distance neuroanatomic connections. Because spatiotemporal dynamics of this kind are central to brain physiology, inverse solutions could be improved by incorporating models of these dynamics. In this article, we present a model for cortical activity based on nearest-neighbor autoregression that incorporates local spatiotemporal interactions between distributed sources in a manner consistent with neurophysiology and neuroanatomy. We develop a dynamic Maximum a Posteriori Expectation-Maximization (dMAP-EM) source localization algorithm for estimation of cortical sources and model parameters based on the Kalman Filter, the Fixed Interval Smoother, and the EM algorithms. We apply the dMAP-EM algorithm to simulated experiments as well as to human experimental data. Furthermore, we derive expressions to relate our dynamic estimation formulas to those of standard static models, and show how dynamic methods optimally assimilate past and future data. Our results establish the feasibility of spatiotemporal dynamic estimation in large-scale distributed source spaces with several thousand source locations and hundreds of sensors, with resulting inverse solutions that provide substantial performance improvements over static methods.

It has been long appreciated that anesthetic drugs induce stereotyped changes in electroencephalogram (EEG), but the relationships between EEG and underlying brain function remain poorly understood. Functional imaging methods including positron emission tomography (PET) and functional magnetic resonance imaging (fMRI), have become important tools for studying how anesthetic drugs act in the human brain to induce the state of general anesthesia. To date, no investigation has combined functional MRI with EEG to study general anesthesia. We report here a paradigm for conducting combined fMRI and EEG studies of human subjects under general anesthesia. We discuss the several technical and safety problems that must be solved to undertake this type of multimodal functional imaging and show combined recordings from a human subject. Combined fMRI and EEG exploits simultaneously the high spatial resolution of fMRI and the high temporal resolution of EEG. In addition, combined fMRI and EEG offers a direct way to relate established EEG patterns induced by general anesthesia to changes neural activity in specific brain regions as measured by changes in fMRI blood oxygen level dependent (BOLD) signals.

Magnetoencephalography (MEG) is an important non-invasive method for studying activity within the human brain. Source localization methods can be used to estimate spatiotemporal activity from MEG measurements with high temporal resolution, but the spatial resolution of these estimates is poor due to the ill-posed nature of the MEG inverse problem. Recent developments in source localization methodology have emphasized temporal as well as spatial constraints to improve source localization accuracy, but these methods can be computationally intense. Solutions emphasizing spatial sparsity hold tremendous promise, since the underlying neurophysiological processes generating MEG signals are often sparse in nature, whether in the form of focal sources, or distributed sources representing large-scale functional networks. Recent developments in the theory of compressed sensing (CS) provide a rigorous framework to estimate signals with sparse structure. In particular, a class of CS algorithms referred to as greedy pursuit algorithms can provide both high recovery accuracy and low computational complexity. Greedy pursuit algorithms are difficult to apply directly to the MEG inverse problem because of the high-dimensional structure of the MEG source space and the high spatial correlation in MEG measurements. In this paper, we develop a novel greedy pursuit algorithm for sparse MEG source localization that overcomes these fundamental problems. This algorithm, which we refer to as the Subspace Pursuit-based Iterative Greedy Hierarchical (SPIGH) inverse solution, exhibits very low computational complexity while achieving very high localization accuracy. We evaluate the performance of the proposed algorithm using comprehensive simulations, as well as the analysis of human MEG data during spontaneous brain activity and somatosensory stimuli. These studies reveal substantial performance gains provided by the SPIGH algorithm in terms of computational complexity, localization accuracy, and robustness.

Our approach provides substantial performance improvements over existing methods using computationally efficient algorithms that will facilitate practical applications in both neuroscience and medicine.

Dynamic estimation methods based on linear state-space models have been applied to the inverse problem of magnetoencephalography (MEG), and can improve source localization compared with static methods by incorporating temporal continuity as a constraint. The efficacy of these methods is influenced by how well the state-space model approximates the dynamics of the underlying brain current sources. While some components of the state-space model can be inferred from brain anatomy and knowledge of the MEG instrument noise structure, parameters governing the temporal evolution of underlying current sources are unknown and must be selected on an ad-hoc basis or estimated from data. In this work, we apply the Expectation-Maximization (EM) algorithm to estimate parameters and sources in an MEG state-space model, and demonstrate in simulation studies that the resulting source estimates are superior to those provided by static methods or dynamic methods employing ad hoc parameter selection.

We develop a state space approach to multimodal integration of simultaneously recorded EEG and fMRI. The EEG is represented with a distributed current source model using realistic MRI-based forward models, whose temporal evolution is governed by a linear state space model. The fMRI signal is similarly modeled by a linear state space model describing the hemodynamic response to underlying EEG current activity. We explore the feasibility of high dimensional dynamic estimation of simultaneous EEG/fMRI using simulation studies of the alpha wave.

Recent dynamic source localization algorithms for the Magnetoencephalographic inverse problem use cortical spatio-temporal dynamics to enhance the quality of the estimation. However, these methods suffer from high computational complexity due to the large number of sources that must be estimated. In this work, we introduce a fast iterative greedy algorithm incorporating the class of subspace pursuit algorithms for sparse source localization. The algorithm employs a reduced order state-space model resulting in significant computational savings. Simulation studies on MEG source localization reveal substantial gains provided by the proposed method over the widely used minimum-norm estimate, in terms of localization accuracy, with a negligible increase in computational complexity.

Cortical activity can be estimated from electroencephalogram (EEG) or magnetoencephalogram (MEG) data by solving an ill-conditioned inverse problem that is regularized using neuroanatomical, computational, and dynamic constraints. Recent methods have incorporated spatio-temporal dynamics into the inverse problem framework. In this approach, spatio-temporal interactions between neighboring sources enforce a form of spatial smoothing that enhances source localization quality. However, spatial smoothing could also occur by way of correlations within the state noise process that drives the underlying dynamic model. Estimating the spatial covariance structure of this state noise is challenging, particularly in EEG and MEG data where the number of underlying sources is far greater than the number of sensors. However, the EEG/MEG data are sparse compared to the large number of sources, and thus sparse constraints could be used to simplify the form of the state noise spatial covariance. In this work, we introduce an empirically tailored basis to represent the spatial covariance structure within the state noise processes of a cortical dynamic model for EEG source localization. We augment the method presented in Lamus, et al. (2011) to allow for sparsity enforcing priors on the covariance parameters. Simulation studies as well as analysis of real data reveal significant gains in the source localization performance over existing algorithms.