Quantifying the performance of MEG source reconstruction using resting state data

"What" and "when" predictions modulate auditory processing in a mutually congruent manner

Introduction

Extracting regularities from ongoing stimulus streams to form predictions is crucial for adaptive behavior. Such regularities exist in terms of the content of the stimuli and their timing, both of which are known to interactively modulate sensory processing. In real-world stimulus streams such as music, regularities can occur at multiple levels, both in terms of contents (e.g., predictions relating to individual notes vs. their more complex groups) and timing (e.g., pertaining to timing between intervals vs. the overall beat of a musical phrase). However, it is unknown whether the brain integrates predictions in a manner that is mutually congruent (e.g., if "beat" timing predictions selectively interact with "what" predictions falling on pulses which define the beat), and whether integrating predictions in different timing conditions relies on dissociable neural correlates.

Methods

To address these questions, our study manipulated "what" and "when" predictions at different levels - (local) interval-defining and (global) beat-defining - within the same stimulus stream, while neural activity was recorded using electroencephalogram (EEG) in participants (N = 20) performing a repetition detection task.

Results

Our results reveal that temporal predictions based on beat or interval timing modulated mismatch responses to violations of "what" predictions happening at the predicted time points, and that these modulations were shared between types of temporal predictions in terms of the spatiotemporal distribution of EEG signals. Effective connectivity analysis using dynamic causal modeling showed that the integration of "what" and "when" predictions selectively increased connectivity at relatively late cortical processing stages, between the superior temporal gyrus and the fronto-parietal network.

Discussion

Taken together, these results suggest that the brain integrates different predictions with a high degree of mutual congruence, but in a shared and distributed cortical network. This finding contrasts with recent studies indicating separable mechanisms for beat-based and memory-based predictive processing.

High precision magnetoencephalography reveals increased right-inferior frontal gyrus beta power during response conflict

Flexibility of behavior and the ability to rapidly switch actions is critical for adaptive living in humans. It is well established that the right-inferior frontal gyrus (R-IFG) is recruited during outright action-stopping, relating to increased beta (12-30 Hz) power. It has also been posited that inhibiting incorrect response tendencies and switching is central to motor flexibility. However, it is not known if the commonly reported R-IFG beta signature of response inhibition in action-stopping is also recruited during response conflict, which would suggest overlapping networks for stopping and switching. In the current study, we analyzed high precision magnetoencephalography (hpMEG) data recorded with multiple within subject recording sessions (trials n > 10,000) from 8 subjects during different levels of response conflict. We hypothesized that a R-IFG-triggered network for response inhibition is domain general and therefore also involved in mediating response conflict. We tested whether R-IFG showed increased beta power dependent on the level of response conflict. Using event-related spectral perturbations and linear mixed modeling, we found that R-IFG beta power increased for response conflict trials. The R-IFG beta increase was specific to trials with strong response conflict, and increased R-IFG beta power related to less error. This supports a more generalized role for R-IFG beta, beyond simple stopping behavior towards response switching.

Beyond dimension reduction: Stable electric fields emerge from and allow representational drift

It is known that the exact neurons maintaining a given memory (the neural ensemble) change from trial to trial. This raises the question of how the brain achieves stability in the face of this representational drift. Here, we demonstrate that this stability emerges at the level of the electric fields that arise from neural activity. We show that electric fields carry information about working memory content. The electric fields, in turn, can act as "guard rails" that funnel higher dimensional variable neural activity along stable lower dimensional routes. We obtained the latent space associated with each memory. We then confirmed the stability of the electric field by mapping the latent space to different cortical patches (that comprise a neural ensemble) and reconstructing information flow between patches. Stable electric fields can allow latent states to be transferred between brain areas, in accord with modern engram theory.

Virtual intracranial EEG signals reconstructed from MEG with potential for epilepsy surgery

Modelling the interactions that arise from neural dynamics in seizure genesis is challenging but important in the effort to improve the success of epilepsy surgery. Dynamical network models developed from physiological evidence offer insights into rapidly evolving brain networks in the epileptic seizure. A limitation of previous studies in this field is the dependence on invasive cortical recordings with constrained spatial sampling of brain regions that might be involved in seizure dynamics. Here, we propose virtual intracranial electroencephalography (ViEEG), which combines non-invasive ictal magnetoencephalographic imaging (MEG), dynamical network models and a virtual resection technique. In this proof-of-concept study, we show that ViEEG signals reconstructed from MEG alone preserve critical temporospatial characteristics for dynamical approaches to identify brain areas involved in seizure generation. We show the non-invasive ViEEG approach may have some advantage over intracranial electroencephalography (iEEG). Future work may be designed to test the potential of the virtual iEEG approach for use in surgical management of epilepsy.

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.

Age- and gender-specific characteristics of the resting-state brain activity: a magnetoencephalography study

Aging and gender influence regional brain activities. Although these biases should be considered during the clinical examinations using magnetoencephalography, they have yet to be standardized. In the present study, resting-state magnetoencephalography data were recorded from 54 healthy females and 48 males aged 22 to 75 years, who were controlled for cognitive performance. The regional oscillatory power was estimated for each frequency band (delta, theta, alpha, beta, low-gamma, and high-gamma) using the sLORETA-like algorithm and the biases of age and gender were evaluated, respectively. The results showed that faster oscillatory powers increased with age in the rostral regions and decreased in the caudal regions, while few slower oscillatory powers changed with age. Gender differences in oscillatory powers were found in a broad frequency range, mostly in the caudal brain regions. The present study characterized the effects of healthy aging and gender asymmetricity on the regional resting-state brain activity, with the aim to facilitate the accurate and efficient use of magnetoencephalography in clinical practice.

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.