Source reconstruction accuracy of MEG and EEG Bayesian inversion approaches

Feedforward connectivity patterns from visual areas to the front of the brain contain information about sensory stimuli regardless of awareness or report

Current theories of consciousness can be categorized to some extent by their predictions about the putative role of the prefrontal cortex (PFC) in conscious perception. One family of the theories proposes that the PFC is necessary for conscious perception. The other postulates that the PFC is not necessary and that other areas (e.g., posterior cortical areas) are more important for conscious perception. No-report paradigms could potentially arbitrate the debate as they disentangle task reporting from conscious perception. While previous no-report paradigms tend to point to a reduction in PFC activity, they have not examined the critical role of the PFC in "monitoring" or "reading out" the patterns of activity in the sensory cortex to generate conscious perception. To address this, we reanalysed electroencephalography (EEG) data from a no-report inattentional blindness paradigm (Shafto & Pitts, 2015). We examined the role of feedforward input patterns to the PFC from sensory cortices. We employed nonparametric spectral Granger causality and quantified the amount of information that reflected the contents of consciousness using multivariate classifiers. Unexpectedly, regardless of whether the stimulus was consciously seen or not, we found that information relating to the current sensory stimulus was present in the pattern of inputs from visual areas to the PFC. In light of these findings, we suggest various theories of consciousness need to be revised to accommodate the fact that the contents of consciousness are decodable from the input patterns from posterior sensory regions to the PFC, regardless of awareness (or report).

Combining OPM and lesion mapping data for epilepsy surgery planning: a simulation study

When planning for epilepsy surgery, multiple potential sites for resection may be identified through anatomical imaging. Magnetoencephalography (MEG) using optically pumped sensors (OP-MEG) is a non-invasive functional neuroimaging technique which could be used to help identify the epileptogenic zone from these candidate regions. Here we test the utility of a-priori information from anatomical imaging for differentiating potential lesion sites with OP-MEG. We investigate a number of scenarios: whether to use rigid or flexible sensor arrays, with or without a-priori source information and with or without source modelling errors. We simulated OP-MEG recordings for 1309 potential lesion sites identified from anatomical images in the Multi-centre Epilepsy Lesion Detection (MELD) project. To localise the simulated data, we used three source inversion schemes: unconstrained, prior source locations at centre of the candidate sites, and prior source locations within a volume around the lesion location. We found that prior knowledge of the candidate lesion zones made the inversion robust to errors in sensor gain, orientation and even location. When the reconstruction was too highly restricted and the source assumptions were inaccurate, the utility of this a-priori information was undermined. Overall, we found that constraining the reconstruction to the region including and around the participant's potential lesion sites provided the best compromise of robustness against modelling or measurement error.

Distinctive modes of cortical communications in tactile temporal order judgment

Temporal order judgment of two successive tactile stimuli delivered to our hands is often inverted when we cross our hands. The present study aimed to identify time-frequency profiles of the interactions across the cortical network associated with the crossed-hand tactile temporal order judgment task using magnetoencephalography. We found that the interactions across the cortical network were channeled to a low-frequency band (5-10 Hz) when the hands were uncrossed. However, the interactions became activated in a higher band (12-18 Hz) when the hands were crossed. The participants with fewer inverted judgments relied mainly on the higher band, whereas those with more frequent inverted judgments (reversers) utilized both. Moreover, reversers showed greater cortical interactions in the higher band when their judgment was correct compared to when it was inverted. Overall, the results show that the cortical network communicates in two distinctive frequency modes during the crossed-hand tactile temporal order judgment task. A default mode of communications in the low-frequency band encourages inverted judgments, and correct judgment is robustly achieved by recruiting the high-frequency mode.

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.

A Crucial Role of the Frontal Operculum in Task-Set Dependent Visuomotor Performance Monitoring

For adaptive goal-directed action, the brain needs to monitor action performance and detect errors. The corresponding information may be conveyed via different sensory modalities; for instance, visual and proprioceptive body position cues may inform about current manual action performance. Thereby, contextual factors such as the current task set may also determine the relative importance of each sensory modality for action guidance. Here, we analyzed human behavioral, functional magnetic resonance imaging (fMRI), and magnetoencephalography (MEG) data from two virtual reality-based hand-target phase-matching studies to identify the neuronal correlates of performance monitoring and error processing under instructed visual or proprioceptive task sets. Our main result was a general, modality-independent response of the bilateral frontal operculum (FO) to poor phase-matching accuracy, as evident from increased BOLD signal and increased source-localized gamma power. Furthermore, functional connectivity of the bilateral FO to the right posterior parietal cortex (PPC) increased under a visual versus proprioceptive task set. These findings suggest that the bilateral FO generally monitors manual action performance; and, moreover, that when visual action feedback is used to guide action, the FO may signal an increased need for control to visuomotor regions in the right PPC following errors.

Optimizing EEG Source Reconstruction with Concurrent fMRI-Derived Spatial Priors

Reconstructing EEG sources involves a complex pipeline, with the inverse problem being the most challenging. Multiple inversion algorithms are being continuously developed, aiming to tackle the non-uniqueness of this problem, which has been shown to be partially circumvented by including prior information in the inverse models. Despite a few efforts, there are still current and persistent controversies regarding the inversion algorithm of choice and the optimal set of spatial priors to be included in the inversion models. The use of simultaneous EEG-fMRI data is one approach to tackle this problem. The spatial resolution of fMRI makes fMRI derived spatial priors very convenient for EEG reconstruction, however, only task activation maps and resting-state networks (RSNs) have been explored so far, overlooking the recent, but already accepted, notion that brain networks exhibit dynamic functional connectivity fluctuations. The lack of a systematic comparison between different source reconstruction algorithms, considering potentially more brain-informative priors such as fMRI, motivates the search for better reconstruction models. Using simultaneous EEG-fMRI data, here we compared four different inversion algorithms (minimum norm, MN; low resolution electromagnetic tomography, LORETA; empirical Bayes beamformer, EBB; and multiple sparse priors, MSP) under a Bayesian framework (as implemented in SPM), each with three different sets of priors consisting of: (1) those specific to the algorithm; (2) those specific to the algorithm plus fMRI task activation maps and RSNs; and (3) those specific to the algorithm plus fMRI task activation maps and RSNs and network modules of task-related dFC states estimated from the dFC fluctuations. The quality of the reconstructed EEG sources was quantified in terms of model-based metrics, namely the expectation of the posterior probability P(model|data) and variance explained of the inversion models, and the overlap/proportion of brain regions known to be involved in the visual perception tasks that the participants were submitted to, and RSN templates, with/within EEG source components. Model-based metrics suggested that model parsimony is preferred, with the combination MSP and priors specific to this algorithm exhibiting the best performance. However, optimal overlap/proportion values were found using EBB and priors specific to this algorithm and fMRI task activation maps and RSNs or MSP and considering all the priors (algorithm priors, fMRI task activation maps and RSNs and dFC state modules), respectively, indicating that fMRI spatial priors, including dFC state modules, might contain useful information to recover EEG source components reflecting neuronal activity of interest. Our main results show that providing fMRI spatial derived priors that reflect the dynamics of the brain might be useful to map neuronal activity more accurately from EEG-fMRI. Furthermore, this work paves the way towards a more informative selection of the optimal EEG source reconstruction approach, which may be critical in future studies.

Modelling optically pumped magnetometer interference in MEG as a spatially homogeneous magnetic field

Here we propose that much of the magnetic interference observed when using optically pumped magnetometers for MEG experiments can be modeled as a spatially homogeneous magnetic field. We show that this approximation reduces sensor level variance and substantially improves statistical power. This model does not require knowledge of the underlying neuroanatomy nor the sensor positions. It only needs information about the sensor orientation. Due to the model's low rank there is little risk of removing substantial neural signal. However, we provide a framework to assess this risk for any sensor number, design or subject neuroanatomy. We find that the risk of unintentionally removing neural signal is reduced when multi-axis recordings are performed. We validated the method using a binaural auditory evoked response paradigm and demonstrated that removing the homogeneous magnetic field increases sensor level SNR by a factor of 3. Considering the model's simplicity and efficacy, we suggest that this homogeneous field correction can be a powerful preprocessing step for arrays of optically pumped magnetometers.

Information redundancy across spatial scales modulates early visual cortical processing

Visual images contain redundant information across spatial scales where low spatial frequency contrast is informative towards the location and likely content of high spatial frequency detail. Previous research suggests that the visual system makes use of those redundancies to facilitate efficient processing. In this framework, a fast, initial analysis of low-spatial frequency (LSF) information guides the slower and later processing of high spatial frequency (HSF) detail. Here, we used multivariate classification as well as time-frequency analysis of MEG responses to the viewing of intact and phase scrambled images of human faces to demonstrate that the availability of redundant LSF information, as found in broadband intact images, correlates with a reduction in HSF representational dominance in both early and higher-level visual areas as well as a reduction of gamma-band power in early visual cortex. Our results indicate that the cross spatial frequency information redundancy that can be found in all natural images might be a driving factor in the efficient integration of fine image details.

Laminar dynamics of high amplitude beta bursts in human motor cortex

Motor cortical activity in the beta frequency range is one of the strongest and most studied movement-related neural signals. At the single trial level, beta band activity is often characterized by transient, high amplitude, bursting events rather than slowly modulating oscillations. The timing of these bursting events is tightly linked to behavior, suggesting a more dynamic functional role for beta activity than previously believed. However, the neural mechanisms underlying beta bursts in sensorimotor circuits are poorly understood. To address this, we here leverage and extend recent developments in high precision MEG for temporally resolved laminar analysis of burst activity, combined with a neocortical circuit model that simulates the biophysical generators of the electrical currents which drive beta bursts. This approach pinpoints the generation of beta bursts in human motor cortex to distinct excitatory synaptic inputs to deep and superficial cortical layers, which drive current flow in opposite directions. These laminar dynamics of beta bursts in motor cortex align with prior invasive animal recordings within the somatosensory cortex, and suggest a conserved mechanism for somatosensory and motor cortical beta bursts. More generally, we demonstrate the ability for uncovering the laminar dynamics of event-related neural signals in human non-invasive recordings. This provides important constraints to theories about the functional role of burst activity for movement control in health and disease, and crucial links between macro-scale phenomena measured in humans and micro-circuit activity recorded from animal models.

A systematic evaluation of source reconstruction of resting MEG of the human brain with a new high-resolution atlas: Performance, precision, and parcellation

Noninvasive functional neuroimaging of the human brain can give crucial insight into the mechanisms that underpin healthy cognition and neurological disorders. Magnetoencephalography (MEG) measures extracranial magnetic fields originating from neuronal activity with high temporal resolution, but requires source reconstruction to make neuroanatomical inferences from these signals. Many source reconstruction algorithms are available, and have been widely evaluated in the context of localizing task-evoked activities. However, no consensus yet exists on the optimum algorithm for resting-state data. Here, we evaluated the performance of six commonly-used source reconstruction algorithms based on minimum-norm and beamforming estimates. Using human resting-state MEG, we compared the algorithms using quantitative metrics, including resolution properties of inverse solutions and explained variance in sensor-level data. Next, we proposed a data-driven approach to reduce the atlas from the Human Connectome Project's multi-modal parcellation of the human cortex based on metrics such as MEG signal-to-noise-ratio and resting-state functional connectivity gradients. This procedure produced a reduced cortical atlas with 230 regions, optimized to match the spatial resolution and the rank of MEG data from the current generation of MEG scanners. Our results show that there is no "one size fits all" algorithm, and make recommendations on the appropriate algorithms depending on the data and aimed analyses. Our comprehensive comparisons and recommendations can serve as a guide for choosing appropriate methodologies in future studies of resting-state MEG.

Testing covariance models for MEG source reconstruction of hippocampal activity

Beamforming is one of the most commonly used source reconstruction methods for magneto- and electroencephalography (M/EEG). One underlying assumption, however, is that distant sources are uncorrelated and here we tested whether this is an appropriate model for the human hippocampal data. We revised the Empirical Bayesian Beamfomer (EBB) to accommodate specific a-priori correlated source models. We showed in simulation that we could use model evidence (as approximated by Free Energy) to distinguish between different correlated and uncorrelated source scenarios. Using group MEG data in which the participants performed a hippocampal-dependent task, we explored the possibility that the hippocampus or the cortex or both were correlated in their activity across hemispheres. We found that incorporating a correlated hippocampal source model significantly improved model evidence. Our findings help to explain why, up until now, the majority of MEG-reported hippocampal activity (typically making use of beamformers) has been estimated as unilateral.

A review of magnetoencephalography use in pediatric epilepsy: an update on best practice

Introduction: Magnetoencephalography (MEG) is a noninvasive technique that is used for presurgical evaluation of children with drug-resistant epilepsy (DRE).Areas covered: The contributions of MEG for localizing the epileptogenic zone are discussed, in particular in extra-temporal lobe epilepsy and focal cortical dysplasia, which are common in children, as well as in difficult to localize epilepsy such as operculo-insular epilepsy. Further, the authors review current evidence on MEG for mapping eloquent cortex, its performance, application in clinical practice, and potential challenges.Expert opinion: MEG could change the clinical management of children with DRE by directing placement of intracranial electrodes thereby enhancing their yield. With improved identification of a circumscribed epileptogenic zone, MEG could render more patients as suitable candidates for epilepsy surgery and increase utilization of surgery.

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.

Pragmatic spatial sampling for wearable MEG arrays

Several new technologies have emerged promising new Magnetoencephalography (MEG) systems in which the sensors can be placed close to the scalp. One such technology, Optically Pumped MEG (OP-MEG) allows for a scalp mounted system that provides measurements within millimetres of the scalp surface. A question that arises in developing on-scalp systems is: how many sensors are necessary to achieve adequate performance/spatial discrimination? There are many factors to consider in answering this question such as the signal to noise ratio (SNR), the locations and depths of the sources, density of spatial sampling, sensor gain errors (due to interference, subject movement, cross-talk, etc.) and, of course, the desired spatial discrimination. In this paper, we provide simulations which show the impact these factors have on designing sensor arrays for wearable MEG. While OP-MEG has the potential to provide high information content at dense spatial samplings, we find that adequate spatial discrimination of sources (< 1 cm) can be achieved with relatively few sensors (< 100) at coarse spatial samplings (~ 30 mm) at high SNR. After this point approximately 50 more sensors are required for every 1 mm improvement in spatial discrimination. Comparable discrimination for traditional cryogenic systems require more channels by these same metrics. We also show that sensor gain errors have the greatest impact on discrimination between deep sources at high SNR. Finally, we also examine the limitation that aliasing due to undersampling has on the effective SNR of on-scalp sensors.

Estimates of cortical column orientation improve MEG source inversion

Determining the anatomical source of brain activity non-invasively measured from EEG or MEG sensors is challenging. In order to simplify the source localization problem, many techniques introduce the assumption that current sources lie on the cortical surface. Another common assumption is that this current flow is orthogonal to the cortical surface, thereby approximating the orientation of cortical columns. However, it is not clear which cortical surface to use to define the current source locations, and normal vectors computed from a single cortical surface may not be the best approximation to the orientation of cortical columns. We compared three different surface location priors and five different approaches for estimating dipole vector orientation, both in simulations and visual and motor evoked MEG responses. We show that models with source locations on the white matter surface and using methods based on establishing correspondences between white matter and pial cortical surfaces dramatically outperform models with source locations on the pial or combined pial/white surfaces and which use methods based on the geometry of a single cortical surface in fitting evoked visual and motor responses. These methods can be easily implemented and adopted in most M/EEG analysis pipelines, with the potential to significantly improve source localization of evoked responses.

Data-driven model optimization for optically pumped magnetometer sensor arrays

Optically pumped magnetometers (OPMs) have reached sensitivity levels that make them viable portable alternatives to traditional superconducting technology for magnetoencephalography (MEG). OPMs do not require cryogenic cooling and can therefore be placed directly on the scalp surface. Unlike cryogenic systems, based on a well-characterised fixed arrays essentially linear in applied flux, OPM devices, based on different physical principles, present new modelling challenges. Here, we outline an empirical Bayesian framework that can be used to compare between and optimise sensor arrays. We perturb the sensor geometry (via simulation) and with analytic model comparison methods estimate the true sensor geometry. The width of these perturbation curves allows us to compare different MEG systems. We test this technique using simulated and real data from SQUID and OPM recordings using head-casts and scanner-casts. Finally, we show that given knowledge of underlying brain anatomy, it is possible to estimate the true sensor geometry from the OPM data themselves using a model comparison framework. This implies that the requirement for accurate knowledge of the sensor positions and orientations a priori may be relaxed. As this procedure uses the cortical manifold as spatial support there is no co-registration procedure or reliance on scalp landmarks.

Human motor cortical beta bursts relate to movement planning and response errors

Motor cortical beta activity (13-30 Hz) is a hallmark signature of healthy and pathological movement, but its behavioural relevance remains unclear. Using high-precision magnetoencephalography (MEG), we show that during the classical event-related desynchronisation (ERD) and event-related synchronisation (ERS) periods, motor cortical beta activity in individual trials (n > 12,000) is dominated by high amplitude, transient, and infrequent bursts. Beta burst probability closely matched the trial-averaged beta amplitude in both the pre- and post-movement periods, but individual bursts were spatially more focal than the classical ERS peak. Furthermore, prior to movement (ERD period), beta burst timing was related to the degree of motor preparation, with later bursts resulting in delayed response times. Following movement (ERS period), the first beta burst was delayed by approximately 100 milliseconds when an incorrect response was made. Overall, beta burst timing was a stronger predictor of single trial behaviour than beta burst rate or single trial beta amplitude. This transient nature of motor cortical beta provides new constraints for theories of its role in information processing within and across cortical circuits, and its functional relevance for behaviour in both healthy and pathological movement.

High-precision magnetoencephalography for reconstructing amygdalar and hippocampal oscillations during prediction of safety and threat

Learning to associate neutral with aversive events in rodents is thought to depend on hippocampal and amygdala oscillations. In humans, oscillations underlying aversive learning are not well characterised, largely due to the technical difficulty of recording from these two structures. Here, we used high‐precision magnetoencephalography (MEG) during human discriminant delay threat conditioning. We constructed generative anatomical models relating neural activity with recorded magnetic fields at the single‐participant level, including the neocortex with or without the possibility of sources originating in the hippocampal and amygdalar structures. Models including neural activity in amygdala and hippocampus explained MEG data during threat conditioning better than exclusively neocortical models. We found that in both amygdala and hippocampus, theta oscillations during anticipation of an aversive event had lower power compared to safety, both during retrieval and extinction of aversive memories. At the same time, theta synchronisation between hippocampus and amygdala increased over repeated retrieval of aversive predictions, but not during safety. Our results suggest that high‐precision MEG is sensitive to neural activity of the human amygdala and hippocampus during threat conditioning and shed light on the oscillation‐mediated mechanisms underpinning retrieval and extinction of fear memories in humans.

Quantitative Evaluation in Estimating Sources Underlying Brain Oscillations Using Current Source Density Methods and Beamformer Approaches

Brain oscillations from EEG and MEG shed light on neurophysiological mechanisms of human behavior. However, to extract information on cortical processing, researchers have to rely on source localization methods that can be very broadly classified into current density estimates such as exact low-resolution brain electromagnetic tomography (eLORETA), minimum norm estimates (MNE), and beamformers such as dynamic imaging of coherent sources (DICS) and linearly constrained minimum variance (LCMV). These algorithms produce a distributed map of brain activity underlying sustained and transient responses during neuroimaging studies of behavior. On the other hand, there are very few comparative analyses that evaluates the “ground truth detection” capabilities of these methods. The current article evaluates the reliability in estimation of sources of spectral event generators in the cortex using a two-pronged approach. First, simulated EEG data with point dipoles and distributed dipoles are used to validate the accuracy and sensitivity of each one of these methods of source localization. The abilities of the techniques were tested by comparing the localization error, focal width, false positive (FP) ratios while detecting already known location of neural activity generators under varying signal-to-noise ratios (SNRs). Second, empirical EEG data during auditory steady state responses (ASSRs) in human participants were used to compare the distributed nature of source localization. All methods were successful in recovery of point sources in favorable signal to noise scenarios and could achieve high hit rates if FPs are ignored. Interestingly, focal activation map is generated by LCMV and DICS when compared to eLORETA while control of FPs is much superior in eLORETA. Subsequently drawbacks and strengths of each method are highlighted with a detailed discussion on how to choose a technique based on empirical requirements.

Rapid Extraction of Emotion Regularities from Complex Scenes in the Human Brain

Adaptive behavior requires the rapid extraction of behaviorally relevant information in the environment, with particular emphasis on emotional cues. However, the speed of emotional feature extraction from complex visual environments is largely undetermined. Here we use objective electrophysiological recordings in combination with frequency tagging to demonstrate that the extraction of emotional information from neutral, pleasant, or unpleasant naturalistic scenes can be completed at a presentation speed of 167 ms (i.e., 6 Hz) under high perceptual load. Emotional compared to neutral pictures evoked enhanced electrophysiological responses with distinct topographical activation patterns originating from different neural sources. Cortical facilitation in early visual cortex was also more pronounced for scenes with pleasant compared to unpleasant or neutral content, suggesting a positivity offset mechanism dominating under conditions of rapid scene processing. These results significantly advance our knowledge of complex scene processing in demonstrating rapid integrative content identification, particularly for emotional cues relevant for adaptive behavior in complex environments.

Multimodal Integration of M/EEG and f/MRI Data in SPM12

We describe the steps involved in analysis of multi-modal, multi-subject human neuroimaging data using the SPM12 free and open source software (https://www.fil.ion.ucl.ac.uk/spm/) and a publically-available dataset organized according to the Brain Imaging Data Structure (BIDS) format (https://openneuro.org/datasets/ds000117/). The dataset contains electroencephalographic (EEG), magnetoencephalographic (MEG), and functional and structural magnetic resonance imaging (MRI) data from 16 subjects who undertook multiple runs of a simple task performed on a large number of famous, unfamiliar and scrambled faces. We demonstrate: (1) batching and scripting of preprocessing of multiple runs/subjects of combined MEG and EEG data, (2) creation of trial-averaged evoked responses, (3) source-reconstruction of the power (induced and evoked) across trials within a time-frequency window around the "N/M170" evoked component, using structural MRI for forward modeling and simultaneous inversion (fusion) of MEG and EEG data, (4) group-based optimisation of spatial priors during M/EEG source reconstruction using fMRI data on the same paradigm, and (5) statistical mapping across subjects of cortical source power increases for faces vs. scrambled faces.

Quantifying the performance of MEG source reconstruction using resting state data

In magnetoencephalography (MEG) research there are a variety of inversion methods to transform sensor data into estimates of brain activity. Each new inversion scheme is generally justified against a specific simulated or task scenario. The choice of this scenario will however have a large impact on how well the scheme performs. We describe a method with minimal selection bias to quantify algorithm performance using human resting state data. These recordings provide a generic, heterogeneous, and plentiful functional substrate against which to test different MEG recording and reconstruction approaches. We used a Hidden Markov model to spatio-temporally partition data into self-similar dynamic states. To test the anatomical precision that could be achieved, we then inverted these data onto libraries of systematically distorted subject-specific cortical meshes and compared the quality of the fit using cross validation and a Free energy metric. This revealed which inversion scheme was able to identify the least distorted (most accurate) anatomical models, and allowed us to quantify an upper bound on the mean anatomical distortion accordingly. We used two resting state datasets, one recorded with head-casts and one without. In the head-cast data, the Empirical Bayesian Beamformer (EBB) algorithm showed the best mean anatomical discrimination (3.7 mm) compared with Minimum Norm/LORETA (6.0 mm) and Multiple Sparse Priors (9.4 mm). This pattern was replicated in the second (conventional dataset) although with a marginally poorer (non-significant) prediction of the missing (cross-validated) data. Our findings suggest that the abundant resting state data now commonly available could be used to refine and validate MEG source reconstruction methods and/or recording paradigms.

Lamina-specific cortical dynamics in human visual and sensorimotor cortices

Distinct anatomical and spectral channels are thought to play specialized roles in the communication within cortical networks. While activity in the alpha and beta frequency range (7 - 40 Hz) is thought to predominantly originate from infragranular cortical layers conveying feedback-related information, activity in the gamma range (>40 Hz) dominates in supragranular layers communicating feedforward signals. We leveraged high precision MEG to test this proposal, directly and non-invasively, in human participants performing visually cued actions. We found that visual alpha mapped onto deep cortical laminae, whereas visual gamma predominantly occurred more superficially. This lamina-specificity was echoed in movement-related sensorimotor beta and gamma activity. These lamina-specific pre- and post- movement changes in sensorimotor beta and gamma activity suggest a more complex functional role than the proposed feedback and feedforward communication in sensory cortex. Distinct frequency channels thus operate in a lamina-specific manner across cortex, but may fulfill distinct functional roles in sensory and motor processes.

Frequency-based brain networks: From a multiplex framework to a full multilayer description

We explore how to study dynamical interactions between brain regions by using functional multilayer networks whose layers represent different frequency bands at which a brain operates. Specifically, we investigate the consequences of considering the brain as (i) a multilayer network, in which all brain regions can interact with each other at different frequency bands; and as (ii) a multiplex network, in which interactions between different frequency bands are allowed only within each brain region and not between them. We study the second-smallest eigenvalue λ 2 of the combinatorial supra-Laplacian matrix of both the multiplex and multilayer networks, as λ 2 has been used previously as an indicator of network synchronizability and as a biomarker for several brain diseases. We show that the heterogeneity of interlayer edge weights and, especially, the fraction of missing edges crucially modify the value of λ 2, and we illustrate our results with both synthetic network models and real data obtained from resting-state magnetoencephalography. Our work highlights the differences between using a multiplex approach and a full multilayer approach when studying frequency-based multilayer brain networks.

NMDA-receptor antibodies alter cortical microcircuit dynamics

Recently, autoantibodies against NMDA receptors (NMDARs) were identified as a major cause of autoimmune encephalitis. They cause abnormalities in brain function often associated with significant changes in patients’ brain dynamics. Here we use computational modeling to identify how NMDAR dysfunction causes abnormalities in brain dynamics using patient EEGs and local field potential recordings in a mouse model of NMDAR-Ab encephalitis. NMDAR autoantibodies cause a specific shift in excitatory coupling within cortical circuits that places the circuits closer to pathological transitions between dynamic brain states. Because of the proximity to these phase transitions, otherwise benign fluctuations in neuronal coupling cause abnormal EEG responses in the presence of the antibodies. Our modeling results thus explain fluctuating abnormalities in brain dynamics observed in patients.

NMDA-receptor antibodies (NMDAR-Abs) cause an autoimmune encephalitis with a diverse range of EEG abnormalities. NMDAR-Abs are believed to disrupt receptor function, but how blocking this excitatory synaptic receptor can lead to paroxysmal EEG abnormalities—or even seizures—is poorly understood. Here we show that NMDAR-Abs change intrinsic cortical connections and neuronal population dynamics to alter the spectral composition of spontaneous EEG activity and predispose brain dynamics to paroxysmal abnormalities. Based on local field potential recordings in a mouse model, we first validate a dynamic causal model of NMDAR-Ab effects on cortical microcircuitry. Using this model, we then identify the key synaptic parameters that best explain EEG paroxysms in pediatric patients with NMDAR-Ab encephalitis. Finally, we use the mouse model to show that NMDAR-Ab–related changes render microcircuitry critically susceptible to overt EEG paroxysms when these key parameters are changed, even though the same parameter fluctuations are tolerated in the in silico model of the control condition. These findings offer mechanistic insights into circuit-level dysfunction induced by NMDAR-Ab.

Inhibition in the face of emotion: Characterization of the spatial-temporal dynamics that facilitate automatic emotion regulation

Emotion regulation mediates socio-cognitive functions and is essential for interactions with others. The capacity to automatically inhibit responses to emotional stimuli is an important aspect of emotion regulation; the underlying neural mechanisms of this ability have been rarely investigated. Forty adults completed a Go/No-go task during magnetoencephalographic (MEG) recordings, where they responded rapidly to either a blue or purple frame which contained angry or happy faces. Subjects responded to the target color in an inhibition (75% Go trials) and a vigilance condition (25% Go trials). As expected, inhibition processes showed early, sustained activation (200-450 ms) in the right inferior frontal gyrus (IFG). Emotion-related inhibition processes showed greater activity with angry faces bilaterally in the orbital-frontal gyri (OFG) starting at 225 ms and temporal poles from 250 ms, with right hemisphere dominance. The presence of happy faces elicited earlier activity in the right OFG. This study demonstrates that the timing of inhibition processes varies with the emotional context and that there is much greater activation in the presence of angry faces. It underscores the importance of the right IFG for inhibition processes, but the OFG in automatic emotion regulation.

Non-invasive laminar inference with MEG: Comparison of methods and source inversion algorithms

Magnetoencephalography (MEG) is a direct measure of neuronal current flow; its anatomical resolution is therefore not constrained by physiology but rather by data quality and the models used to explain these data. Recent simulation work has shown that it is possible to distinguish between signals arising in the deep and superficial cortical laminae given accurate knowledge of these surfaces with respect to the MEG sensors. This previous work has focused around a single inversion scheme (multiple sparse priors) and a single global parametric fit metric (free energy). In this paper we use several different source inversion algorithms and both local and global, as well as parametric and non-parametric fit metrics in order to demonstrate the robustness of the discrimination between layers. We find that only algorithms with some sparsity constraint can successfully be used to make laminar discrimination. Importantly, local t-statistics, global cross-validation and free energy all provide robust and mutually corroborating metrics of fit. We show that discrimination accuracy is affected by patch size estimates, cortical surface features, and lead field strength, which suggests several possible future improvements to this technique. This study demonstrates the possibility of determining the laminar origin of MEG sensor activity, and thus directly testing theories of human cognition that involve laminar- and frequency-specific mechanisms. This possibility can now be achieved using recent developments in high precision MEG, most notably the use of subject-specific head-casts, which allow for significant increases in data quality and therefore anatomically precise MEG recordings.

Section

Analysis methods.

Classifications

Source localization: inverse problem; Source localization: other.

•

Laminar inferences can be made with MEG using both local and global fit metrics.

•

Source inversion algorithms with sparsity constraints performed best.

•

Classification is affected by patch size estimates, anatomy, and lead field strength.

Reconstructing anatomy from electro-physiological data

Here we show how it is possible to make estimates of brain structure based on MEG data. We do this by reconstructing functional estimates onto distorted cortical manifolds parameterised in terms of their spherical harmonics. We demonstrate that both empirical and simulated MEG data give rise to consistent and plausible anatomical estimates. Importantly, the estimation of structure from MEG data can be quantified in terms of millimetres from the true brain structure. We show, for simulated data, that the functional assumptions which are closer to the functional ground-truth give rise to anatomical estimates that are closer to the true anatomy.

Low-frequency oscillations employ a general coding of the spatio-temporal similarity of dynamic faces

Brain networks use neural oscillations as information transfer mechanisms. Although the face perception network in occipitotemporal cortex is well-studied, contributions of oscillations to face representation remain an open question. We tested for links between oscillatory responses that encode facial dimensions and the theoretical proposal that faces are encoded in similarity-based "face spaces". We quantified similarity-based encoding of dynamic faces in magnetoencephalographic sensor-level oscillatory power for identity, expression, physical and perceptual similarity of facial form and motion. Our data show that evoked responses manifest physical and perceptual form similarity that distinguishes facial identities. Low-frequency induced oscillations (< 20Hz) manifested more general similarity structure, which was not limited to identity, and spanned physical and perceived form and motion. A supplementary fMRI-constrained source reconstruction implicated fusiform gyrus and V5 in this similarity-based representation. These findings introduce a potential link between "face space" encoding and oscillatory network communication, which generates new hypotheses about the potential oscillation-mediated mechanisms that might encode facial dimensions.

The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming

Despite numerous important contributions, the investigation of brain connectivity with magnetoencephalography (MEG) still faces multiple challenges. One critical aspect of source-level connectivity, largely overlooked in the literature, is the putative effect of the choice of the inverse method on the subsequent cortico-cortical coupling analysis. We set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, thousands of randomly located pairs of sources were created. Several parameters were manipulated, including inter- and intra-source correlation strength, source size and spatial configuration. The simulated pairs of sources were then used to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, the source level power and coherence maps were calculated using three methods (a) L2-Minimum-Norm Estimate (MNE), (b) Linearly Constrained Minimum Variance (LCMV) beamforming, and (c) Dynamic Imaging of Coherent Sources (DICS) beamforming. The performances of the methods were evaluated using Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions if the interacting cortical sources consist of point-like sources. On the other hand, MNE provides better connectivity estimation than beamformers, if the interacting sources are simulated as extended cortical patches, where each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each patch of active cortex is simulated with only partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. The insights revealed here can guide method selection and help improve data interpretation regarding MEG connectivity estimation.

Using generative models to make probabilistic statements about hippocampal engagement in MEG

Magnetoencephalography (MEG) enables non-invasive real time characterization of brain activity. However, convincing demonstrations of signal contributions from deeper sources such as the hippocampus remain controversial and are made difficult by its depth, structural complexity and proximity to neocortex. Here, we demonstrate a method for quantifying hippocampal engagement probabilistically using simulated hippocampal activity and realistic anatomical and electromagnetic source modelling. We construct two generative models, one which supports neuronal current flow on the cortical surface, and one which supports neuronal current flow on both the cortical and hippocampal surface. Using Bayesian model comparison, we then infer which of the two models provides a more likely explanation of the dataset at hand. We also carry out a set of control experiments to rule out bias, including simulating medial temporal lobe sources to assess the risk of falsely positive results, and adding different types of displacements to the hippocampal portion of the mesh to test for anatomical specificity of the results. In addition, we test the robustness of this inference by adding co-registration error and sensor level noise. We find that the model comparison framework is sensitive to hippocampal activity when co-registration error is <3 mm and the sensor-level signal-to-noise ratio (SNR) is >-20 dB. These levels of co-registration error and SNR can now be achieved empirically using recently developed subject-specific head-casts.

Kernel temporal enhancement approach for LORETA source reconstruction using EEG data

Reconstruction of brain sources from magnetoencephalography and electroencephalography (M/EEG) data is a well known problem in the neuroengineering field. A inverse problem should be solved and several methods have been proposed. Low Resolution Electromagnetic Tomography (LORETA) and the different variations proposed as standardized LORETA (sLORETA) and the standardized weighted LORETA (swLORETA) have solved the inverse problem following a non-parametric approach, that is by setting dipoles in the whole brain domain in order to estimate the dipole positions from the M/EEG data and assuming some spatial priors. Errors in the reconstruction of sources are presented due the low spatial resolution of the LORETA framework and the influence of noise in the observable data. In this work a kernel temporal enhancement (kTE) is proposed in order to build a preprocessing stage of the data that allows in combination with the swLORETA method a improvement in the source reconstruction. The results are quantified in terms of three dipole error localization metrics and the strategy of swLORETA + kTE obtained the best results across different signal to noise ratio (SNR) in random dipoles simulation from synthetic EEG data.

Delta and gamma oscillations in operculo-insular cortex underlie innocuous cold thermosensation

Using magnetoencephalography, we identified spatiotemporal features of central cold processing, with respect to the time course, oscillatory profile, and neural generators of cold-evoked responses in healthy human volunteers. Cold thermosensation was associated with low- and high-frequency oscillatory rhythms, both originating in operculo-insular regions. These results support further investigations of central cold processing using magnetoencephalography or EEG and the clinical utility of cold-evoked potentials for neurophysiological assessment of cold-related small-fiber function and damage.

Cold-sensitive and nociceptive neural pathways interact to shape the quality and intensity of thermal and pain perception. Yet the central processing of cold thermosensation in the human brain has not been extensively studied. Here, we used magnetoencephalography and EEG in healthy volunteers to investigate the time course (evoked fields and potentials) and oscillatory activity associated with the perception of cold temperature changes. Nonnoxious cold stimuli consisting of Δ3°C and Δ5°C decrements from an adapting temperature of 35°C were delivered on the dorsum of the left hand via a contact thermode. Cold-evoked fields peaked at around 240 and 500 ms, at peak latencies similar to the N1 and P2 cold-evoked potentials. Importantly, cold-related changes in oscillatory power indicated that innocuous thermosensation is mediated by oscillatory activity in the range of delta (1–4 Hz) and gamma (55–90 Hz) rhythms, originating in operculo-insular cortical regions. We suggest that delta rhythms coordinate functional integration between operculo-insular and frontoparietal regions, while gamma rhythms reflect local sensory processing in operculo-insular areas.

NEW & NOTEWORTHY Using magnetoencephalography, we identified spatiotemporal features of central cold processing, with respect to the time course, oscillatory profile, and neural generators of cold-evoked responses in healthy human volunteers. Cold thermosensation was associated with low- and high-frequency oscillatory rhythms, both originating in operculo-insular regions. These results support further investigations of central cold processing using magnetoencephalography or EEG and the clinical utility of cold-evoked potentials for neurophysiological assessment of cold-related small-fiber function and damage.

Comparison of group-level, source localized activity for simultaneous functional near-infrared spectroscopy-magnetoencephalography and simultaneous fNIRS-fMRI during parametric median nerve stimulation

Functional near-infrared spectroscopy (fNIRS) is a noninvasive neuroimaging technique, which uses light to measure changes in cerebral blood oxygenation through sensors placed on the surface of the scalp. We recorded concurrent fNIRS with magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI) in order to investigate the group-level correspondence of these measures with source-localized fNIRS estimates. Healthy participants took part in both a concurrent fNIRS-MEG and fNIRS-fMRI neuroimaging session during two somatosensory stimulation tasks, a blocked design median nerve localizer and parametric pulsed-pair median nerve stimulation using interpulse intervals from 100 to 500 ms. We found the spatial correlation for estimated activation patterns from the somatosensory task was [Formula: see text], 0.57, and [Formula: see text] and the amplitude correlation was [Formula: see text], 0.52, and [Formula: see text] for fMRI-MEG, fMRI-fNIRS oxy-hemoglobin, and fMRI-fNIRS deoxy-hemoglobin signals, respectively. Taken together, these results show good correspondence among the fMRI, fNIRS, and MEG with the great majority of the difference across modalities being driven by lower sensitivity for deeper brain sources in MEG and fNIRS. These results provide an important validation of source-localized fNIRS in the context of concurrent multimodal imaging for future studies of the relationship between physiological effects in the human brain.

Combined EMD-sLORETA Analysis of EEG Data Collected during a Contour Integration Task

Lately, Ensemble Empirical Mode Decomposition (EEMD) techniques receive growing interest in biomedical data analysis. Event-Related Modes (ERMs) represent features extracted by an EEMD from electroencephalographic (EEG) recordings. We present a new approach for source localization of EEG data based on combining ERMs with inverse models. As the first step, 64 channel EEG recordings are pooled according to six brain areas and decomposed, by applying an EEMD, into their underlying ERMs. Then, based upon the problem at hand, the most closely related ERM, in terms of frequency and amplitude, is combined with inverse modeling techniques for source localization. More specifically, the standardized low resolution brain electromagnetic tomography (sLORETA) procedure is employed in this work. Accuracy and robustness of the results indicate that this approach deems highly promising in source localization techniques for EEG data.

Non-linear Parameter Estimates from Non-stationary MEG Data

We demonstrate a method to estimate key electrophysiological parameters from resting state data. In this paper, we focus on the estimation of head-position parameters. The recovery of these parameters is especially challenging as they are non-linearly related to the measured field. In order to do this we use an empirical Bayesian scheme to estimate the cortical current distribution due to a range of laterally shifted head-models. We compare different methods of approaching this problem from the division of M/EEG data into stationary sections and performing separate source inversions, to explaining all of the M/EEG data with a single inversion. We demonstrate this through estimation of head position in both simulated and empirical resting state MEG data collected using a head-cast.

Simplified EEG inverse solution for BCI real-time implementation

EEG brain imaging has become a promising approach in Brain-computer interface applications. However, accurate reconstruction of active regions and computational burden are still open issues. In this paper, we propose to use a simplified forward model that includes the reduction of the cortical dipoles based on Brodmann areas together with state-of-the-art EEG brain imaging techniques. With this approach the well known Beamformers and Greedy Search inverse solutions become feasible for real-time implementation, while guaranteeing lower localization error than previous approaches used in BCI. This methodology was tested with synthetic and real EEG data from a visual attention study. Results show zero localization error in terms of active cortical regions estimation in single 1 s trial datasets, with a computation time of 1.1 s in a non-specialized personal computer. These results open the possibility to obtain in real-time information of active cortical regions in Brain-computer interfaces.

Emotion Discrimination Using Spatially Compact Regions of Interest Extracted from Imaging EEG Activity

Lately, research on computational models of emotion had been getting much attention due to their potential for understanding the mechanisms of emotions and their promising broad range of applications that potentially bridge the gap between human and machine interactions. We propose a new method for emotion classification that relies on features extracted from those active brain areas that are most likely related to emotions. To this end, we carry out the selection of spatially compact regions of interest that are computed using the brain neural activity reconstructed from Electroencephalography data. Throughout this study, we consider three representative feature extraction methods widely applied to emotion detection tasks, including Power spectral density, Wavelet, and Hjorth parameters. Further feature selection is carried out using principal component analysis. For validation purpose, these features are used to feed a support vector machine classifier that is trained under the leave-one-out cross-validation strategy. Obtained results on real affective data show that incorporation of the proposed training method in combination with the enhanced spatial resolution provided by the source estimation allows improving the performed accuracy of discrimination in most of the considered emotions, namely: dominance, valence, and liking.

Reconstruction of Neural Activity from EEG Data Using Dynamic Spatiotemporal Constraints

We present a novel iterative regularized algorithm (IRA) for neural activity reconstruction that explicitly includes spatiotemporal constraints, performing a trade-off between space and time resolutions. For improving the spatial accuracy provided by electroencephalography (EEG) signals, we explore a basis set that describes the smooth, localized areas of potentially active brain regions. In turn, we enhance the time resolution by adding the Markovian assumption for brain activity estimation at each time period. Moreover, to deal with applications that have either distributed or localized neural activity, the spatiotemporal constraints are expressed through [Formula: see text] and [Formula: see text] norms, respectively. For the purpose of validation, we estimate the neural reconstruction performance in time and space separately. Experimental testing is carried out on artificial data, simulating stationary and non-stationary EEG signals. Also, validation is accomplished on two real-world databases, one holding Evoked Potentials and another with EEG data of focal epilepsy. Moreover, responses of functional magnetic resonance imaging for the former EEG data have been measured in advance, allowing to contrast our findings. Obtained results show that the [Formula: see text]-based IRA produces a spatial resolution that is comparable to the one achieved by some widely used sparse-based estimators of brain activity. At the same time, the [Formula: see text]-based IRA outperforms other similar smooth solutions, providing a spatial resolution that is lower than the sparse [Formula: see text]-based solution. As a result, the proposed IRA is a promising method for improving the accuracy of brain activity reconstruction.

Dynamic causal modelling of electrographic seizure activity using Bayesian belief updating

Seizure activity in EEG recordings can persist for hours with seizure dynamics changing rapidly over time and space. To characterise the spatiotemporal evolution of seizure activity, large data sets often need to be analysed. Dynamic causal modelling (DCM) can be used to estimate the synaptic drivers of cortical dynamics during a seizure; however, the requisite (Bayesian) inversion procedure is computationally expensive. In this note, we describe a straightforward procedure, within the DCM framework, that provides efficient inversion of seizure activity measured with non-invasive and invasive physiological recordings; namely, EEG/ECoG. We describe the theoretical background behind a Bayesian belief updating scheme for DCM. The scheme is tested on simulated and empirical seizure activity (recorded both invasively and non-invasively) and compared with standard Bayesian inversion. We show that the Bayesian belief updating scheme provides similar estimates of time-varying synaptic parameters, compared to standard schemes, indicating no significant qualitative change in accuracy. The difference in variance explained was small (less than 5%). The updating method was substantially more efficient, taking approximately 5–10 min compared to approximately 1–2 h. Moreover, the setup of the model under the updating scheme allows for a clear specification of how neuronal variables fluctuate over separable timescales. This method now allows us to investigate the effect of fast (neuronal) activity on slow fluctuations in (synaptic) parameters, paving a way forward to understand how seizure activity is generated.

•

We describe a DCM procedure that provides efficient inversion of seizure activity.

•

Similar accuracy but substantially more efficient compared to standard DCM methods.

•

Physiological fluctuations over different timescales can be specified.

•

This scheme should contribute to understanding seizure activity using DCM.

Characterising seizures in anti-NMDA-receptor encephalitis with dynamic causal modelling

We characterised the pathophysiology of seizure onset in terms of slow fluctuations in synaptic efficacy using EEG in patients with anti-N-methyl-d-aspartate receptor (NMDA-R) encephalitis. EEG recordings were obtained from two female patients with anti-NMDA-R encephalitis with recurrent partial seizures (ages 19 and 31). Focal electrographic seizure activity was localised using an empirical Bayes beamformer. The spectral density of reconstructed source activity was then characterised with dynamic causal modelling (DCM). Eight models were compared for each patient, to evaluate the relative contribution of changes in intrinsic (excitatory and inhibitory) connectivity and endogenous afferent input. Bayesian model comparison established a role for changes in both excitatory and inhibitory connectivity during seizure activity (in addition to changes in the exogenous input). Seizures in both patients were associated with a sequence of changes in inhibitory and excitatory connectivity; a transient increase in inhibitory connectivity followed by a transient increase in excitatory connectivity and a final peak of excitatory–inhibitory balance at seizure offset. These systematic fluctuations in excitatory and inhibitory gain may be characteristic of (anti NMDA-R encephalitis) seizures. We present these results as a case study and replication to motivate analyses of larger patient cohorts, to see whether our findings generalise and further characterise the mechanisms of seizure activity in anti-NMDA-R encephalitis.

•

We characterised seizures in patient with anti-NMDA-R encephalitis using EEG.

•

Dynamic causal modelling was used to estimate causes of seizure activity.

•

Characteristic variation of excitatory–inhibitory balance during seizure activity.

•

This variation was seen for seizures within and between patients.

Discrimination of cortical laminae using MEG

Typically MEG source reconstruction is used to estimate the distribution of current flow on a single anatomically derived cortical surface model. In this study we use two such models representing superficial and deep cortical laminae. We establish how well we can discriminate between these two different cortical layer models based on the same MEG data in the presence of different levels of co-registration noise, Signal-to-Noise Ratio (SNR) and cortical patch size. We demonstrate that it is possible to make a distinction between superficial and deep cortical laminae for levels of co-registration noise of less than 2 mm translation and 2° rotation at SNR > 11 dB. We also show that an incorrect estimate of cortical patch size will tend to bias layer estimates. We then use a 3D printed head-cast (Troebinger et al., 2014) to achieve comparable levels of co-registration noise, in an auditory evoked response paradigm, and show that it is possible to discriminate between these cortical layer models in real data.

•

We evaluate necessary recording precision to distinguish superficial/deep laminae.

•

For coregistration error of < 2 mm/2° we can distinguish between these laminar models.

•

Incorrect assumptions about cortical patch size bias these layer estimates.

•

Initial results suggest that the auditory M100 derives from deep cortical layers.

Algorithmic procedures for Bayesian MEG/EEG source reconstruction in SPM

The MEG/EEG inverse problem is ill-posed, giving different source reconstructions depending on the initial assumption sets. Parametric Empirical Bayes allows one to implement most popular MEG/EEG inversion schemes (Minimum Norm, LORETA, etc.) within the same generic Bayesian framework. It also provides a cost-function in terms of the variational Free energy-an approximation to the marginal likelihood or evidence of the solution. In this manuscript, we revisit the algorithm for MEG/EEG source reconstruction with a view to providing a didactic and practical guide. The aim is to promote and help standardise the development and consolidation of other schemes within the same framework. We describe the implementation in the Statistical Parametric Mapping (SPM) software package, carefully explaining each of its stages with the help of a simple simulated data example. We focus on the Multiple Sparse Priors (MSP) model, which we compare with the well-known Minimum Norm and LORETA models, using the negative variational Free energy for model comparison. The manuscript is accompanied by Matlab scripts to allow the reader to test and explore the underlying algorithm.

Optimal spatial filtering for brain oscillatory activity using the Relevance Vector Machine

Over the past decade, various techniques have been proposed for localization of cerebral sources of oscillatory activity on the basis of magnetoencephalography (MEG) or electroencephalography recordings. Beamformers in the frequency domain, in particular, have proved useful in this endeavor. However, the localization accuracy and efficacy of such spatial filters can be markedly limited by bias from correlation between cerebral sources and short duration of source activity, both essential issues in the localization of brain data. Here, we evaluate a method for frequency-domain localization of oscillatory neural activity based on the relevance vector machine (RVM). RVM is a Bayesian algorithm for learning sparse models from possibly overcomplete data sets. The performance of our frequency-domain RVM method (fdRVM) was compared with that of dynamic imaging of coherent sources (DICS), a frequency-domain spatial filter that employs a minimum variance adaptive beamformer (MVAB) approach. The methods were tested both on simulated and real data. Two types of simulated MEG data sets were generated, one with continuous source activity and the other with transiently active sources. The real data sets were from slow finger movements and resting state. Results from simulations show comparable performance for DICS and fdRVM at high signal-to-noise ratios and low correlation. At low SNR or in conditions of high correlation between sources, fdRVM performs markedly better. fdRVM was successful on real data as well, indicating salient focal activations in the sensorimotor area. The resulting high spatial resolution of fdRVM and its sensitivity to low-SNR transient signals could be particularly beneficial when mapping event-related changes of oscillatory activity.