Sparse current source estimation for MEG using loose orientation constraints

A comparison of automated and manual co-registration for magnetoencephalography

Magnetoencephalography (MEG) is a neuroimaging technique that accurately captures the rapid (sub-millisecond) activity of neuronal populations. Interpretation of functional data from MEG relies upon registration to the participant’s anatomical MRI. The key remaining step is to transform the participant’s MRI into the MEG head coordinate space. Although both automated and manual approaches to co-registration are available, the relative accuracy of two approaches has not been systematically evaluated. The goal of the present study was to compare the accuracy of manual and automated co-registration. Resting MEG and T1-weighted MRI data were collected from 90 participants. Automated and manual co-registration were performed on the same subjects, and the inter-method reliability of the two methods assessed using the intra-class correlation. Median co-registration error for both methods was within acceptable limits. Inter-method reliability was in the “good” range for co-registration error, and the “good” to “excellent” range for translation and rotation. These results suggest that the output of the automated co-registration procedure is comparable to that achieved using manual co-registration.

MEG Source Imaging and Group Analysis Using VBMEG

Variational Bayesian Multimodal EncephaloGraphy (VBMEG) is a MATLAB toolbox that estimates distributed source currents from magnetoencephalography (MEG)/electroencephalography (EEG) data by integrating functional MRI (fMRI) (https://vbmeg.atr.jp/). VBMEG also estimates whole-brain connectome dynamics using anatomical connectivity derived from a diffusion MRI (dMRI). In this paper, we introduce the VBMEG toolbox and demonstrate its usefulness. By collaborating with VBMEG's tutorial page (https://vbmeg.atr.jp/docs/v2/static/vbmeg2_tutorial_neuromag.html), we show its full pipeline using an open dataset recorded by Wakeman and Henson (2015). We import the MEG data and preprocess them to estimate the source currents. From the estimated source currents, we perform a group analysis and examine the differences of current amplitudes between conditions by controlling the false discovery rate (FDR), which yields results consistent with previous studies. We highlight VBMEG's characteristics by comparing these results with those obtained by other source imaging methods: weighted minimum norm estimate (wMNE), dynamic statistical parametric mapping (dSPM), and linearly constrained minimum variance (LCMV) beamformer. We also estimate source currents from the EEG data and the whole-brain connectome dynamics from the MEG data and dMRI. The observed results indicate the reliability, characteristics, and usefulness of VBMEG.

Sparse EEG/MEG source estimation via a group lasso

Non-invasive recordings of human brain activity through electroencephalography (EEG) or magnetoencelphalography (MEG) are of value for both basic science and clinical applications in sensory, cognitive, and affective neuroscience. Here we introduce a new approach to estimating the intra-cranial sources of EEG/MEG activity measured from extra-cranial sensors. The approach is based on the group lasso, a sparse-prior inverse that has been adapted to take advantage of functionally-defined regions of interest for the definition of physiologically meaningful groups within a functionally-based common space. Detailed simulations using realistic source-geometries and data from a human Visual Evoked Potential experiment demonstrate that the group-lasso method has improved performance over traditional ℓ2 minimum-norm methods. In addition, we show that pooling source estimates across subjects over functionally defined regions of interest results in improvements in the accuracy of source estimates for both the group-lasso and minimum-norm approaches.

The Iterative Reweighted Mixed-Norm Estimate for Spatio-Temporal MEG/EEG Source Reconstruction

Source imaging based on magnetoencephalography (MEG) and electroencephalography (EEG) allows for the non-invasive analysis of brain activity with high temporal and good spatial resolution. As the bioelectromagnetic inverse problem is ill-posed, constraints are required. For the analysis of evoked brain activity, spatial sparsity of the neuronal activation is a common assumption. It is often taken into account using convex constraints based on the l1-norm. The resulting source estimates are however biased in amplitude and often suboptimal in terms of source selection due to high correlations in the forward model. In this work, we demonstrate that an inverse solver based on a block-separable penalty with a Frobenius norm per block and a l0.5-quasinorm over blocks addresses both of these issues. For solving the resulting non-convex optimization problem, we propose the iterative reweighted Mixed Norm Estimate (irMxNE), an optimization scheme based on iterative reweighted convex surrogate optimization problems, which are solved efficiently using a block coordinate descent scheme and an active set strategy. We compare the proposed sparse imaging method to the dSPM and the RAP-MUSIC approach based on two MEG data sets. We provide empirical evidence based on simulations and analysis of MEG data that the proposed method improves on the standard Mixed Norm Estimate (MxNE) in terms of amplitude bias, support recovery, and stability.

Functional characterisation of letter-specific responses in time, space and current polarity using magnetoencephalography

Recent neurophysiological evidence suggests that a hierarchical neural network of low-to-high level processing subserves written language comprehension. While a considerable amount of research has identified distinct regions and stages of processing, the relations between them and to this hierarchical model remain unclear. Magnetoencephalography (MEG) is a technique frequently employed in such investigations; however, no studies have sought to test whether the conventional method of reconstructing currents at the source of the magnetic field is best suited for such across-subject designs. The present study details the results of three MEG experiments addressing these issues. Neuronal populations supporting responses to low-level orthographic properties were housed posteriorly near the primary visual cortex. More anterior regions along the fusiform gyrus encoded higher-level processes and became active ~80ms later. A functional localiser of these early letter-specific responses was developed for the production of functional regions of interest in future studies. Previously established response components were successfully grouped based on proximity to the localiser, which characterised location, latency and functional sensitivity. Unconventional anatomically constrained signed minimum norm estimates of MEG data were most sensitive to the primary experimental manipulation, suggesting that the conventional unsigned unconstrained method is sub-optimal for studying written word processing.

Frequency specific interactions of MEG resting state activity within and across brain networks as revealed by the multivariate interaction measure

Resting state networks (RSNs) are sets of brain regions exhibiting temporally coherent activity fluctuations in the absence of imposed task structure. RSNs have been extensively studied with fMRI in the infra-slow frequency range (nominally <10(-1)Hz). The topography of fMRI RSNs reflects stationary temporal correlation over minutes. However, neuronal communication occurs on a much faster time scale, at frequencies nominally in the range of 10(0)-10(2)Hz. We examined phase-shifted interactions in the delta (2-3.5 Hz), theta (4-7 Hz), alpha (8-12 Hz) and beta (13-30 Hz) frequency bands of resting-state source space MEG signals. These analyses were conducted between nodes of the dorsal attention network (DAN), one of the most robust RSNs, and between the DAN and other networks. Phase shifted interactions were mapped by the multivariate interaction measure (MIM), a measure of true interaction constructed from the maximization of imaginary coherency in the virtual channels comprised of voxel signals in source space. Non-zero-phase interactions occurred between homologous left and right hemisphere regions of the DAN in the delta and alpha frequency bands. Even stronger non-zero-phase interactions were detected between networks. Visual regions bilaterally showed phase-shifted interactions in the alpha band with regions of the DAN. Bilateral somatomotor regions interacted with DAN nodes in the beta band. These results demonstrate the existence of consistent, frequency specific phase-shifted interactions on a millisecond time scale between cortical regions within RSN as well as across RSNs.