Discrimination of cortical laminae using MEG

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.

The promise of layer-specific neuroimaging for testing predictive coding theories of psychosis

Predictive coding potentially provides an explanatory model for understanding the neurocognitive mechanisms of psychosis. It proposes that cognitive processes, such as perception and inference, are implemented by a hierarchical system, with the influence of each level being a function of the estimated precision of beliefs at that level. However, predictive coding models of psychosis are insufficiently constrained-any phenomenon can be explained in multiple ways by postulating different changes to precision at different levels of processing. One reason for the lack of constraint in these models is that the core processes are thought to be implemented by the function of specific cortical layers, and the technology to measure layer specific neural activity in humans has until recently been lacking. As a result, our ability to constrain the models with empirical data has been limited. In this review we provide a brief overview of predictive processing models of psychosis and then describe the potential for newly developed, layer specific neuroimaging techniques to test and thus constrain these models. We conclude by discussing the most promising avenues for this research as well as the technical and conceptual challenges which may limit its application.

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.

An Impending Paradigm Shift in Motor Imagery Based Brain-Computer Interfaces

The development of reliable assistive devices for patients that suffer from motor impairments following central nervous system lesions remains a major challenge in the field of non-invasive Brain-Computer Interfaces (BCIs). These approaches are predominated by electroencephalography and rely on advanced signal processing and machine learning methods to extract neural correlates of motor activity. However, despite tremendous and still ongoing efforts, their value as effective clinical tools remains limited. We advocate that a rather overlooked research avenue lies in efforts to question neurophysiological markers traditionally targeted in non-invasive motor BCIs. We propose an alternative approach grounded by recent fundamental advances in non-invasive neurophysiology, specifically subject-specific feature extraction of sensorimotor bursts of activity recorded via (possibly magnetoencephalography-optimized) electroencephalography. This path holds promise in overcoming a significant proportion of existing limitations, and could foster the wider adoption of online BCIs in rehabilitation protocols.

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.

Non-invasively measured brain activity and radiological progression in diffuse glioma

Non-invasively measured brain activity is related to progression-free survival in glioma patients, suggesting its potential as a marker of glioma progression. We therefore assessed the relationship between brain activity and increasing tumor volumes on routine clinical magnetic resonance imaging (MRI) in glioma patients. Postoperative magnetoencephalography (MEG) was recorded in 45 diffuse glioma patients. Brain activity was estimated using three measures (absolute broadband power, offset and slope) calculated at three spatial levels: global average, averaged across the peritumoral areas, and averaged across the homologues of these peritumoral areas in the contralateral hemisphere. Tumors were segmented on MRI. Changes in tumor volume between the two scans surrounding the MEG were calculated and correlated with brain activity. Brain activity was compared between patient groups classified into having increasing or stable tumor volume. Results show that brain activity was significantly increased in the tumor hemisphere in general, and in peritumoral regions specifically. However, none of the measures and spatial levels of brain activity correlated with changes in tumor volume, nor did they differ between patients with increasing versus stable tumor volumes. Longitudinal studies in more homogeneous subgroups of glioma patients are necessary to further explore the clinical potential of non-invasively measured brain activity.

Co-registration Comparison of On-Scalp Magnetoencephalography and Magnetic Resonance Imaging

Magnetoencephalography (MEG) can non-invasively measure the electromagnetic activity of the brain. A new type of MEG, on-scalp MEG, has attracted the attention of researchers recently. Compared to the conventional SQUID-MEG, on-scalp MEG constructed with optically pumped magnetometers is wearable and has a high signal-to-noise ratio. While the co-registration between MEG and magnetic resonance imaging (MRI) significantly influences the source localization accuracy, co-registration error requires assessment, and quantification. Recent studies have evaluated the co-registration error of on-scalp MEG mainly based on the surface fit error or the repeatability error of different measurements, which do not reflect the true co-registration error. In this study, a three-dimensional-printed reference phantom was constructed to provide the ground truth of MEG sensor locations and orientations relative to MRI. The co-registration performances of commonly used three devices—electromagnetic digitization system, structured-light scanner, and laser scanner—were compared and quantified by the indices of final co-registration errors in the reference phantom and human experiments. Furthermore, the influence of the co-registration error on the performance of source localization was analyzed via simulations. The laser scanner had the best co-registration accuracy (rotation error of 0.23° and translation error of 0.76 mm based on the phantom experiment), whereas the structured-light scanner had the best cost performance. The results of this study provide recommendations and precautions for researchers regarding selecting and using an appropriate device for the co-registration of on-scalp MEG and MRI.

Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function

BACKGROUND

The specificity of gradient echo (GE)-BOLD laminar fMRI activation profiles is degraded by intracortical veins that drain blood from lower to upper cortical layers, propagating activation signal in the same direction. This work describes an approach to obtain layer specific profiles by deconvolving the measured profiles with a physiological Point Spread Function (PSF).

NEW METHOD

It is shown that the PSF can be characterised by a TE-dependent peak to tail (p2t) value that is independent of cortical depth and can be estimated by simulation. An experimental estimation of individual p2t values and the sensitivity of the deconvolved profiles to variations in p2t is obtained using laminar data measured with a multi-echo 3D-FLASH sequence. These profiles are echo time dependent, but the underlying neuronal response is the same, allowing a data-based estimation of the PSF.

RESULTS

The deconvolved profiles are highly similar to the gold-standard obtained from extremely high resolution 3D-EPI data, for a range of p2t values of 5-9, which covers both the empirically determined value (6.8) and the value obtained by simulation (6.3). -Comparison with Existing Method(s) Corrected profiles show a flatter shape across the cortex and a high level of similarity with the gold-standard, defined as a subset of profiles that are unaffected by intracortical veins.

CONCLUSIONS

We conclude that deconvolution is a robust approach for removing the effect of signal propagation through intracortical veins. This makes it possible to obtain profiles with high laminar specificity while benefitting from the higher efficiency of GE-BOLD sequences.

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.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

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.

BigBrain 3D atlas of cortical layers: Cortical and laminar thickness gradients diverge in sensory and motor cortices

Histological atlases of the cerebral cortex, such as those made famous by Brodmann and von Economo, are invaluable for understanding human brain microstructure and its relationship with functional organization in the brain. However, these existing atlases are limited to small numbers of manually annotated samples from a single cerebral hemisphere, measured from 2D histological sections. We present the first whole-brain quantitative 3D laminar atlas of the human cerebral cortex. It was derived from a 3D histological atlas of the human brain at 20-micrometer isotropic resolution (BigBrain), using a convolutional neural network to segment, automatically, the cortical layers in both hemispheres. Our approach overcomes many of the historical challenges with measurement of histological thickness in 2D, and the resultant laminar atlas provides an unprecedented level of precision and detail. We utilized this BigBrain cortical atlas to test whether previously reported thickness gradients, as measured by MRI in sensory and motor processing cortices, were present in a histological atlas of cortical thickness and which cortical layers were contributing to these gradients. Cortical thickness increased across sensory processing hierarchies, primarily driven by layers III, V, and VI. In contrast, motor-frontal cortices showed the opposite pattern, with decreases in total and pyramidal layer thickness from motor to frontal association cortices. These findings illustrate how this laminar atlas will provide a link between single-neuron morphology, mesoscale cortical layering, macroscopic cortical thickness, and, ultimately, functional neuroanatomy.

Using deep learning to segment the layers of the cerebral cortex, this study presents the first whole brain quantitative atlas of cortical and laminar structure. This laminar atlas provides a novel framework for bridging between the scales of neuroscience.

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.

Increased segregation of functional networks in developing brains

A growing literature conceptualises typical brain development from a network perspective. However, largely due to technical and methodological challenges inherent in paediatric functional neuroimaging, there remains an important gap in our knowledge regarding the typical development of functional brain networks in "preschool" childhood (i.e., children younger than 6 years of age). In this study, we recorded brain oscillatory activity using age-appropriate magnetoencephalography in 24 children, including 14 preschool children aged from 4 to 6 years and 10 school children aged from 7 to 12 years. We compared the topology of the resting-state brain networks in these children, estimated using minimum spanning tree (MST) constructed from phase synchrony between beamformer-reconstructed time-series, with that of 24 adults. Our results show that during childhood the MST topology shifts from a star-like (centralised) toward a more line-like (de-centralised) configuration, indicating the functional brain networks become increasingly segregated. In addition, the increasing global network segregation is frequency-independent and accompanied by decreases in centrality (or connectedness) of cortical regions with age, especially in areas of the default mode network. We propose a heuristic MST model of "network space", which posits a clear developmental trajectory for the emergence of complex brain networks. Our results not only revealed topological reorganisation of functional networks across multiple temporal and spatial scales in childhood, but also fill a gap in the literature regarding neurophysiological mechanisms of functional brain maturation during the preschool years of childhood.

Bayesian Modelling of Induced Responses and Neuronal Rhythms

Neural rhythms or oscillations are ubiquitous in neuroimaging data. These spectral responses have been linked to several cognitive processes; including working memory, attention, perceptual binding and neuronal coordination. In this paper, we show how Bayesian methods can be used to finesse the ill-posed problem of reconstructing-and explaining-oscillatory responses. We offer an overview of recent developments in this field, focusing on (i) the use of MEG data and Empirical Bayes to build hierarchical models for group analyses-and the identification of important sources of inter-subject variability and (ii) the construction of novel dynamic causal models of intralaminar recordings to explain layer-specific activity. We hope to show that electrophysiological measurements contain much more spatial information than is often thought: on the one hand, the dynamic causal modelling of non-invasive (low spatial resolution) electrophysiology can afford sub-millimetre (hyper-acute) resolution that is limited only by the (spatial) complexity of the underlying (dynamic causal) forward model. On the other hand, invasive microelectrode recordings (that penetrate different cortical layers) can reveal laminar-specific responses and elucidate hierarchical message passing and information processing within and between cortical regions at a macroscopic scale. In short, the careful and biophysically grounded modelling of sparse data enables one to characterise the neuronal architectures generating oscillations in a remarkable detail.

A future for neuronal oscillation research

Neuronal oscillations represent the most obvious feature of electrical activity in the brain. They are linked in general with global brain state (awake, asleep, etc.) and specifically with organisation of neuronal outputs during sensory perception and cognitive processing. Oscillations can be generated by individual neurons on the basis of interaction between inputs and intrinsic conductances but are far more commonly seen at the local network level in populations of interconnected neurons with diverse arrays of functional properties. It is at this level that the brain’s rich and diverse library of oscillatory time constants serve to temporally organise large-scale neural activity patterns. The discipline is relatively mature at the microscopic (cell, local network) level – although novel discoveries are still commonplace – but requires a far greater understanding of mesoscopic and macroscopic brain dynamics than we currently hold. Without this, extrapolation from the temporal properties of neurons and their communication strategies up to whole brain function will remain largely theoretical. However, recent advances in large-scale neuronal population recordings and more direct, higher fidelity, non-invasive measurement of whole brain function suggest much progress is just around the corner.

Spatial and spectral trajectories in typical neurodevelopment from childhood to middle age

Detailed characterization of typical human neurodevelopment is key if we are to understand the nature of mental and neurological pathology. While research on the cellular processes of neurodevelopment has made great advances, in vivo human imaging is crucial to understand our uniquely human capabilities, as well as the pathologies that affect them. Using magnetoencephalography data in the largest normative sample currently available (324 participants aged 6-45 years), we assess the developmental trajectory of resting-state oscillatory power and functional connectivity from childhood to middle age. The maturational course of power, indicative of local processing, was found to both increase and decrease in a spectrally dependent fashion. Using the strength of phase-synchrony between parcellated regions, we found significant linear and nonlinear (quadratic and logarithmic) trajectories to be characterized in a spatially heterogeneous frequency-specific manner, such as a superior frontal region with linear and nonlinear trajectories in theta and gamma band respectively. Assessment of global efficiency revealed similar significant nonlinear trajectories across all frequency bands. Our results link with the development of human cognitive abilities; they also highlight the complexity of neurodevelopment and provide quantitative parameters for replication and a robust footing from which clinical research may map pathological deviations from these typical trajectories.

Cognitive neuroscience using wearable magnetometer arrays: Non-invasive assessment of language function

Recent work has demonstrated that Optically Pumped Magnetometers (OPMs) can be utilised to create a wearable Magnetoencephalography (MEG) system that is motion robust. In this study, we use this system to map eloquent cortex using a clinically validated language lateralisation paradigm (covert verb generation: 120 trials, ∼10 min total duration) in healthy adults (n = 3). We show that it is possible to lateralise and localise language function on a case by case basis using this system. Specifically, we show that at a sensor and source level we can reliably detect a lateralising beta band (15-30 Hz) desynchronization in all subjects. This is the first study of human cognition using OPMs and not only highlights this technology's utility as tool for (developmental) cognitive neuroscience but also its potential to contribute to surgical planning via mapping of eloquent cortex, especially in young children.

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.

Short timescale abnormalities in the states of spontaneous synchrony in the functional neural networks in Alzheimer's disease

Alzheimer's disease (AD) is a prevalent neurodegenerative condition that can lead to severe cognitive and functional deterioration. Functional magnetic resonance imaging (fMRI) revealed abnormalities in AD in intrinsic synchronization between spatially separate regions in the so-called default mode network (DMN) of the brain. To understand the relationship between this disruption in large-scale synchrony and the cognitive impairment in AD, it is critical to determine whether and how the deficit in the low frequency hemodynamic fluctuations recorded by fMRI translates to much faster timescales of memory and other cognitive processes. The present study employed magnetoencephalography (MEG) and a Hidden Markov Model (HMM) approach to estimate spontaneous synchrony variations in the functional neural networks with high temporal resolution. In a group of cognitively healthy (CH) older adults, we found transient (mean duration of 150-250 ms) network activity states, which were visited in a rapid succession, and were characterized by spatially coordinated changes in the amplitude of source-localized electrophysiological oscillations. The inferred states were similar to those previously observed in younger participants using MEG, and the estimated cortical source distributions of the state-specific activity resembled the classic functional neural networks, such as the DMN. In patients with AD, inferred network states were different from those of the CH group in short-scale timing and oscillatory features. The state of increased oscillatory amplitudes in the regions overlapping the DMN was visited less often in AD and for shorter periods of time, suggesting that spontaneous synchronization in this network was less likely and less stable in the patients. During the visits to this state, in some DMN nodes, the amplitude change in the higher-frequency (8-30 Hz) oscillations was less robust in the AD than CH group. These findings highlight relevance of studying short-scale temporal evolution of spontaneous activity in functional neural networks to understanding the AD pathophysiology. Capacity of flexible intrinsic synchronization in the DMN may be crucial for memory and other higher cognitive functions. Our analysis yielded metrics that quantify distinct features of the neural synchrony disorder in AD and may offer sensitive indicators of the neural network health for future investigations.

Phase reorganization leads to transient β-LFP spatial wave patterns in motor cortex during steady-state movement preparation

Previous studies on the origin and properties of spatial patterns in motor cortex β-local field potential (β-LFP) oscillations have focused on planar traveling waves. However, it is unclear 1) whether β-LFP waves are limited to plane waves, or even 2) whether they are propagating waves of excito-excitatory activity, i.e., primarily traveling waves in excitable media; they could reflect, instead, reorganization in the relative phases of transient oscillations at different spatial sites. We addressed these two problems in β-LFPs recorded via microelectrode arrays implanted in three adjacent motor cortex areas of nonhuman primates during steady-state movement preparation. Our findings are fourfold: 1) β-LFP wave patterns emerged as transient events, despite stable firing rates of single neurons concurrently recorded during the same periods. 2) β-LFP waves showed a richer variety of spatial dynamics, including rotating and complex waves. 3) β-LFP wave patterns showed no characteristic wavelength, presenting instead a range of scales with global zero-lag phase synchrony as a limiting case, features surprising for purely excito-excitatory waves but consistent with waves in coupled oscillator systems. 4) Furthermore, excito-excitatory traveling waves induced by optogenetic stimulation in motor cortex showed, in contrast, a characteristic wavelength and reduced phase synchrony. Overall, β-LFP wave statistics differed from those of induced traveling waves in excitable media recorded under the same microelectrode array setup. Our findings suggest phase reorganization in neural coupled oscillators contribute significantly to the origin of transient β-LFP spatial dynamics during preparatory steady states and outline important constraints for spatially extended models of β-LFP dynamics in motor cortex. NEW & NOTEWORTHY We show that a rich variety of transient β-local field potential (β-LFP) wave patterns emerge in motor cortex during preparatory steady states, despite stable neuronal firing rates. Furthermore, unlike optogenetically induced traveling waves, β-LFP waves showed no characteristic wavelength, presenting instead a range of scales with global phase synchrony as a limiting case. Overall, our statistical analyses suggest that transient phase reorganization in neural coupled oscillators, beyond purely excito-excitatory traveling waves, contribute significantly to the origin of motor cortex β-LFP wave patterns.

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.

Cortical layers, rhythms and BOLD signals

This review investigates how laminar fMRI can complement insights into brain function derived from the study of rhythmic neuronal synchronization. Neuronal synchronization in various frequency bands plays an important role in neuronal communication between brain areas, and it does so on the backbone of layer-specific interareal anatomical projections. Feedforward projections originate predominantly in supragranular cortical layers and terminate in layer 4, and this pattern is reflected in inter-laminar and interareal directed gamma-band influences. Thus, gamma-band synchronization likely subserves feedforward signaling. By contrast, anatomical feedback projections originate predominantly in infragranular layers and terminate outside layer 4, and this pattern is reflected in inter-laminar and interareal directed alpha- and/or beta-band influences. Thus, alpha-beta band synchronization likely subserves feedback signaling. Furthermore, these rhythms explain part of the BOLD signal, with independent contributions of alpha-beta and gamma. These findings suggest that laminar fMRI can provide us with a potentially useful method to test some of the predictions derived from the study of neuronal synchronization. We review central findings regarding the role of layer-specific neuronal synchronization for brain function, and regarding the link between neuronal synchronization and the BOLD signal. We discuss the role that laminar fMRI could play by comparing it to invasive and non-invasive electrophysiological recordings. Compared to direct electrophysiological recordings, this method provides a metric of neuronal activity that is slow and indirect, but that is uniquely non-invasive and layer-specific with potentially whole brain coverage.

Dissecting the Function of Hippocampal Oscillations in a Human Anxiety Model

Neural oscillations in hippocampus and medial prefrontal cortex (mPFC) are a hallmark of rodent anxiety models that build on conflict between approach and avoidance. Yet, the function of these oscillations, and their expression in humans, remain elusive. Here, we used magnetoencephalography (MEG) to investigate neural oscillations in a task that simulated approach–avoidance conflict, wherein 23 male and female human participants collected monetary tokens under a threat of virtual predation. Probability of threat was signaled by color and learned beforehand by direct experience. Magnitude of threat corresponded to a possible monetary loss, signaled as a quantity. We focused our analyses on an a priori defined region-of-interest, the bilateral hippocampus. Oscillatory power under conflict was linearly predicted by threat probability in a location consistent with right mid-hippocampus. This pattern was specific to the hippocampus, most pronounced in the gamma band, and not explained by spatial movement or anxiety-like behavior. Gamma power was modulated by slower theta rhythms, and this theta modulation increased with threat probability. Furthermore, theta oscillations in the same location showed greater synchrony with mPFC theta with increased threat probability. Strikingly, these findings were not seen in relation to an increase in threat magnitude, which was explicitly signaled as a quantity and induced similar behavioral responses as learned threat probability. Thus, our findings suggest that the expression of hippocampal and mPFC oscillatory activity in the context of anxiety is specifically linked to threat memory. These findings resonate with neurocomputational accounts of the role played by hippocampal oscillations in memory.

SIGNIFICANCE STATEMENT We use a biologically relevant approach–avoidance conflict test in humans while recording neural oscillations with magnetoencephalography to investigate the expression and function of hippocampal oscillations in human anxiety. Extending nonhuman studies, we can assign a possible function to hippocampal oscillations in this task, namely threat memory communication. This blends into recent attempts to elucidate the role of brain synchronization in defensive responses to threat.

Magnetoencephalography as a Tool in Psychiatric Research: Current Status and Perspective

The application of neuroimaging to provide mechanistic insights into circuit dysfunctions in major psychiatric conditions and the development of biomarkers are core challenges in current psychiatric research. We propose that recent technological and analytic advances in magnetoencephalography (MEG), a technique that allows measurement of neuronal events directly and noninvasively with millisecond resolution, provides novel opportunities to address these fundamental questions. Because of its potential in delineating normal and abnormal brain dynamics, we propose that MEG provides a crucial tool to advance our understanding of pathophysiological mechanisms of major neuropsychiatric conditions, such as schizophrenia, autism spectrum disorders, and the dementias. We summarize the mechanisms underlying the generation of MEG signals and the tools available to reconstruct generators and underlying networks using advanced source-reconstruction techniques. We then surveyed recent studies that have used MEG to examine aberrant rhythmic activity in neuropsychiatric disorders. This was followed by links with preclinical research that has highlighted possible neurobiological mechanisms, such as disturbances in excitation/inhibition parameters, that could account for measured changes in neural oscillations. Finally, we discuss challenges as well as novel methodological developments that could pave the way for widespread application of MEG in translational research with the aim of developing biomarkers for early detection and diagnosis.

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.

Dynamic causal modelling revisited

This paper revisits the dynamic causal modelling of fMRI timeseries by replacing the usual (Taylor) approximation to neuronal dynamics with a neural mass model of the canonical microcircuit. This provides a generative or dynamic causal model of laminar specific responses that can generate haemodynamic and electrophysiological measurements. In principle, this allows the fusion of haemodynamic and (event related or induced) electrophysiological responses. Furthermore, it enables Bayesian model comparison of competing hypotheses about physiologically plausible synaptic effects; for example, does attentional modulation act on superficial or deep pyramidal cells – or both? In this technical note, we describe the resulting dynamic causal model and provide an illustrative application to the attention to visual motion dataset used in previous papers. Our focus here is on how to answer long-standing questions in fMRI; for example, do haemodynamic responses reflect extrinsic (afferent) input from distant cortical regions, or do they reflect intrinsic (recurrent) neuronal activity? To what extent do inhibitory interneurons contribute to neurovascular coupling? What is the relationship between haemodynamic responses and the frequency of induced neuronal activity? This paper does not pretend to answer these questions; rather it shows how they can be addressed using neural mass models of fMRI timeseries.

•

This paper describes a DCM for fMRI based on neural mass models and canonical microcircuits.

•

This enables the (Bayesian) fusion of EEG and fMRI data.

•

That encompasses the formal modelling of neurovascular coupling.

•

Offers a surprising insight into the relationship between haemodynamic and electrophysiological responses.

Flexible head-casts for high spatial precision MEG

BACKGROUND

In combination with magnetoencephalographic (MEG) data, accurate knowledge of the brain's structure and location provide a principled way of reconstructing neural activity with high temporal resolution. However, measuring the brain's location is compromised by head movement during scanning, and by fiducial-based co-registration with magnetic resonance imaging (MRI) data. The uncertainty from these two factors introduces errors into the forward model and limit the spatial resolution of the data.

NEW METHOD

We present a method for stabilizing and reliably repositioning the head during scanning, and for co-registering MRI and MEG data with low error.

RESULTS

Using this new flexible and comfortable subject-specific head-cast prototype, we find within-session movements of <0.25mm and between-session repositioning errors around 1mm.

COMPARISON WITH EXISTING METHOD(S)

This method is an improvement over existing methods for stabilizing the head or correcting for location shifts on- or off-line, which still introduce approximately 5mm of uncertainty at best (Adjamian et al., 2004; Stolk et al., 2013; Whalen et al., 2008). Further, the head-cast design presented here is more comfortable, safer, and easier to use than the earlier 3D printed prototype, and give slightly lower co-registration errors (Troebinger et al., 2014b).

CONCLUSIONS

We provide an empirical example of how these head-casts impact on source level reproducibility. Employment of the individual flexible head-casts for MEG recordings provide a reliable method of safely stabilizing the head during MEG recordings, and for co-registering MRI anatomical images to MEG functional data.

A new generation of magnetoencephalography: Room temperature measurements using optically-pumped magnetometers

Advances in the field of quantum sensing mean that magnetic field sensors, operating at room temperature, are now able to achieve sensitivity similar to that of cryogenically cooled devices (SQUIDs). This means that room temperature magnetoencephalography (MEG), with a greatly increased flexibility of sensor placement can now be considered. Further, these new sensors can be placed directly on the scalp surface giving, theoretically, a large increase in the magnitude of the measured signal. Here, we present recordings made using a single optically-pumped magnetometer (OPM) in combination with a 3D-printed head-cast designed to accurately locate and orient the sensor relative to brain anatomy. Since our OPM is configured as a magnetometer it is highly sensitive to environmental interference. However, we show that this problem can be ameliorated via the use of simultaneous reference sensor recordings. Using median nerve stimulation, we show that the OPM can detect both evoked (phase-locked) and induced (non-phase-locked oscillatory) changes when placed over sensory cortex, with signals ~4 times larger than equivalent SQUID measurements. Using source modelling, we show that our system allows localisation of the evoked response to somatosensory cortex. Further, source-space modelling shows that, with 13 sequential OPM measurements, source-space signal-to-noise ratio (SNR) is comparable to that from a 271-channel SQUID system. Our results highlight the opportunity presented by OPMs to generate uncooled, potentially low-cost, high SNR MEG systems.

The laminar integration of sensory inputs with feedback signals in human cortex

•

Understanding how the cortex integrates feedback and feedforward signals is central to understanding brain function.

•

The data-driven framework of apical amplification which is hypothesized to have a central role in cognition is highlighted.

•

Human neuroimaging data provides evidence for layer-specific cortical feedback relevant for theories of predictive feedback.

The cortex constitutes the largest area of the human brain. Yet we have only a basic understanding of how the cortex performs one vital function: the integration of sensory signals (carried by feedforward pathways) with internal representations (carried by feedback pathways). A multi-scale, multi-species approach is essential for understanding the site of integration, computational mechanism and functional role of this processing. To improve our knowledge we must rely on brain imaging with improved spatial and temporal resolution and paradigms which can measure internal processes in the human brain, and on the bridging of disciplines in order to characterize this processing at cellular and circuit levels. We highlight apical amplification as one potential mechanism for integrating feedforward and feedback inputs within pyramidal neurons in the rodent brain. We reflect on the challenges and progress in applying this model neuronal process to the study of human cognition. We conclude that cortical-layer specific measures in humans will be an essential contribution for better understanding the landscape of information in cortical feedback, helping to bridge the explanatory gap.

On the Potential of a New Generation of Magnetometers for MEG: A Beamformer Simulation Study

Magnetoencephalography (MEG) is a sophisticated tool which yields rich information on the spatial, spectral and temporal signatures of human brain function. Despite unique potential, MEG is limited by a low signal-to-noise ratio (SNR) which is caused by both the inherently small magnetic fields generated by the brain, and the scalp-to-sensor distance. The latter is limited in current systems due to a requirement for pickup coils to be cryogenically cooled. Recent work suggests that optically-pumped magnetometers (OPMs) might be a viable alternative to superconducting detectors for MEG measurement. They have the advantage that sensors can be brought to within ~4 mm of the scalp, thus offering increased sensitivity. Here, using simulations, we quantify the advantages of hypothetical OPM systems in terms of sensitivity, reconstruction accuracy and spatial resolution. Our results show that a multi-channel whole-head OPM system offers (on average) a fivefold improvement in sensitivity for an adult brain, as well as clear improvements in reconstruction accuracy and spatial resolution. However, we also show that such improvements depend critically on accurate forward models; indeed, the reconstruction accuracy of our simulated OPM system only outperformed that of a simulated superconducting system in cases where forward field error was less than 5%. Overall, our results imply that the realisation of a viable whole-head multi-channel OPM system could generate a step change in the utility of MEG as a means to assess brain electrophysiological activity in health and disease. However in practice, this will require both improved hardware and modelling algorithms.

Let the Rhythm Guide You: Non-invasive Tracking of Cortical Communication Channels

In this issue of Neuron, a study by Michalareas et al. (2016) uses magnetoencephalography (MEG) to characterize the hierarchical organization of human visual areas based on their causal connectivity profiles.

The relationship between oscillatory EEG activity and the laminar-specific BOLD signal

Electrophysiological recordings in animals have indicated that visual cortex γ-band oscillatory activity is predominantly observed in superficial cortical layers, whereas α- and β-band activity is stronger in deep layers. These rhythms, as well as the different cortical layers, have also been closely related to feedforward and feedback streams of information. Recently, it has become possible to measure laminar activity in humans with high-resolution functional MRI (fMRI). In this study, we investigated whether these different frequency bands show a differential relation with the laminar-resolved blood-oxygen level-dependent (BOLD) signal by combining data from simultaneously recorded EEG and fMRI from the early visual cortex. Our visual attention paradigm allowed us to investigate how variations in strength over trials and variations in the attention effect over subjects relate to each other in both modalities. We demonstrate that γ-band EEG power correlates positively with the superficial layers' BOLD signal and that β-power is negatively correlated to deep layer BOLD and α-power to both deep and superficial layer BOLD. These results provide a neurophysiological basis for human laminar fMRI and link human EEG and high-resolution fMRI to systems-level neuroscience in animals.

A multi-layer network approach to MEG connectivity analysis

Recent years have shown the critical importance of inter-regional neural network connectivity in supporting healthy brain function. Such connectivity is measurable using neuroimaging techniques such as MEG, however the richness of the electrophysiological signal makes gaining a complete picture challenging. Specifically, connectivity can be calculated as statistical interdependencies between neural oscillations within a large range of different frequency bands. Further, connectivity can be computed between frequency bands. This pan-spectral network hierarchy likely helps to mediate simultaneous formation of multiple brain networks, which support ongoing task demand. However, to date it has been largely overlooked, with many electrophysiological functional connectivity studies treating individual frequency bands in isolation. Here, we combine oscillatory envelope based functional connectivity metrics with a multi-layer network framework in order to derive a more complete picture of connectivity within and between frequencies. We test this methodology using MEG data recorded during a visuomotor task, highlighting simultaneous and transient formation of motor networks in the beta band, visual networks in the gamma band and a beta to gamma interaction. Having tested our method, we use it to demonstrate differences in occipital alpha band connectivity in patients with schizophrenia compared to healthy controls. We further show that these connectivity differences are predictive of the severity of persistent symptoms of the disease, highlighting their clinical relevance. Our findings demonstrate the unique potential of MEG to characterise neural network formation and dissolution. Further, we add weight to the argument that dysconnectivity is a core feature of the neuropathology underlying schizophrenia.

Measuring electrophysiological connectivity by power envelope correlation: a technical review on MEG methods

The human brain can be divided into multiple areas, each responsible for different aspects of behaviour. Healthy brain function relies upon efficient connectivity between these areas and, in recent years, neuroimaging has been revolutionised by an ability to estimate this connectivity. In this paper we discuss measurement of network connectivity using magnetoencephalography (MEG), a technique capable of imaging electrophysiological brain activity with good (~5 mm) spatial resolution and excellent (~1 ms) temporal resolution. The rich information content of MEG facilitates many disparate measures of connectivity between spatially separate regions and in this paper we discuss a single metric known as power envelope correlation. We review in detail the methodology required to measure power envelope correlation including (i) projection of MEG data into source space, (ii) removing confounds introduced by the MEG inverse problem and (iii) estimation of connectivity itself. In this way, we aim to provide researchers with a description of the key steps required to assess envelope based functional networks, which are thought to represent an intrinsic mode of coupling in the human brain. We highlight the principal findings of the techniques discussed, and furthermore, we show evidence that this method can probe how the brain forms and dissolves multiple transient networks on a rapid timescale in order to support current processing demand. Overall, power envelope correlation offers a unique and verifiable means to gain novel insights into network coordination and is proving to be of significant value in elucidating the neural dynamics of the human connectome in health and disease.

Comparison of different spatial transformations applied to EEG data: A case study of error processing

The purpose of this paper is to compare the effects of different spatial transformations applied to the same scalp-recorded EEG data. The spatial transformations applied are two referencing schemes (average and linked earlobes), the surface Laplacian, and beamforming (a distributed source localization procedure). EEG data were collected during a speeded reaction time task that provided a comparison of activity between error vs. correct responses. Analyses focused on time-frequency power, frequency band-specific inter-electrode connectivity, and within-subject cross-trial correlations between EEG activity and reaction time. Time-frequency power analyses showed similar patterns of midfrontal delta-theta power for errors compared to correct responses across all spatial transformations. Beamforming additionally revealed error-related anterior and lateral prefrontal beta-band activity. Within-subject brain-behavior correlations showed similar patterns of results across the spatial transformations, with the correlations being the weakest after beamforming. The most striking difference among the spatial transformations was seen in connectivity analyses: linked earlobe reference produced weak inter-site connectivity that was attributable to volume conduction (zero phase lag), while the average reference and Laplacian produced more interpretable connectivity results. Beamforming did not reveal any significant condition modulations of connectivity. Overall, these analyses show that some findings are robust to spatial transformations, while other findings, particularly those involving cross-trial analyses or connectivity, are more sensitive and may depend on the use of appropriate spatial transformations.