Moving magnetoencephalography towards real-world applications with a wearable system.

Imaging human brain function with techniques such as magnetoencephalography typically requires a subject to perform tasks while their head remains still within a restrictive scanner. This artificial environment makes the technique inaccessible to many people, and limits the experimental questions that can be addressed. For example, it has been difficult to apply neuroimaging to investigation of the neural substrates of cognitive development in babies and children, or to study processes in adults that require unconstrained head movement (such as spatial navigation). Here we describe a magnetoencephalography system that can be worn like a helmet, allowing free and natural movement during scanning. This is possible owing to the integration of quantum sensors, which do not rely on superconducting technology, with a system for nulling background magnetic fields. We demonstrate human electrophysiological measurement at millisecond resolution while subjects make natural movements, including head nodding, stretching, drinking and playing a ball game. Our results compare well to those of the current state-of-the-art, even when subjects make large head movements. The system opens up new possibilities for scanning any subject or patient group, with myriad applications such as characterization of the neurodevelopmental connectome, imaging subjects moving naturally in a virtual environment and investigating the pathophysiology of movement disorders.

Mouth magnetoencephalography: A unique perspective on the human hippocampus.

Traditional magnetoencephalographic (MEG) brain imaging scanners consist of a rigid sensor array surrounding the head; this means that they are maximally sensitive to superficial brain structures. New technology based on optical pumping means that we can now consider more flexible and creative sensor placement. Here we explored the magnetic fields generated by a model of the human hippocampus not only across scalp but also at the roof of the mouth. We found that simulated hippocampal sources gave rise to dipolar field patterns with one scalp surface field extremum at the temporal lobe and a corresponding maximum or minimum at the roof of the mouth. We then constructed a fitted dental mould to accommodate an Optically Pumped Magnetometer (OPM). We collected data using a previously validated hippocampal-dependant task to test the empirical utility of a mouth-based sensor, with an accompanying array of left and right temporal lobe OPMs. We found that the mouth sensor showed the greatest task-related theta power change. We found that this sensor had a mild effect on the reconstructed power in the hippocampus (~10% change) but that coherence images between the mouth sensor and reconstructed source images showed a global maximum in the right hippocampus. We conclude that augmenting a scalp-based MEG array with sensors in the mouth shows unique promise for both basic scientists and clinicians interested in interrogating the hippocampus.

Impaired theta phase coupling underlies frontotemporal dysconnectivity in schizophrenia.

Frontotemporal dysconnectivity is a key pathology in schizophrenia. The specific nature of this dysconnectivity is unknown, but animal models imply dysfunctional theta phase coupling between hippocampus and medial prefrontal cortex (mPFC). We tested this hypothesis by examining neural dynamics in 18 participants with a schizophrenia diagnosis, both medicated and unmedicated; and 26 age, sex and IQ matched control subjects. All participants completed two tasks known to elicit hippocampal-prefrontal theta coupling: a spatial memory task (during magnetoencephalography) and a memory integration task. In addition, an overlapping group of 33 schizophrenia and 29 control subjects underwent PET to measure the availability of GABAARs expressing the α5 subunit (concentrated on hippocampal somatostatin interneurons). We demonstrate-in the spatial memory task, during memory recall-that theta power increases in left medial temporal lobe (mTL) are impaired in schizophrenia, as is theta phase coupling between mPFC and mTL. Importantly, the latter cannot be explained by theta power changes, head movement, antipsychotics, cannabis use, or IQ, and is not found in other frequency bands. Moreover, mPFC-mTL theta coupling correlated strongly with performance in controls, but not in subjects with schizophrenia, who were mildly impaired at the spatial memory task and no better than chance on the memory integration task. Finally, mTL regions showing reduced phase coupling in schizophrenia magnetoencephalography participants overlapped substantially with areas of diminished α5-GABAAR availability in the wider schizophrenia PET sample. These results indicate that mPFC-mTL dysconnectivity in schizophrenia is due to a loss of theta phase coupling, and imply α5-GABAARs (and the cells that express them) have a role in this process.

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.

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.

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.

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.

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.

Medial Prefrontal-Medial Temporal Theta Phase Coupling in Dynamic Spatial Imagery.

Hippocampal-medial prefrontal interactions are thought to play a crucial role in mental simulation. Notably, the frontal midline/medial pFC (mPFC) theta rhythm in humans has been linked to introspective thought and working memory. In parallel, theta rhythms have been proposed to coordinate processing in the medial temporal cortex, retrosplenial cortex (RSc), and parietal cortex during the movement of viewpoint in imagery, extending their association with physical movement in rodent models. Here, we used noninvasive whole-head MEG to investigate theta oscillatory power and phase-locking during the 18-sec postencoding delay period of a spatial working memory task, in which participants imagined previously learned object sequences either on a blank background (object maintenance), from a first-person viewpoint in a scene (static imagery), or moving along a path past the objects (dynamic imagery). We found increases in 4- to 7-Hz theta power in mPFC when comparing the delay period with a preencoding baseline. We then examined whether the mPFC theta rhythm was phase-coupled with ongoing theta oscillations elsewhere in the brain. The same mPFC region showed significantly higher theta phase coupling with the posterior medial temporal lobe/RSc for dynamic imagery versus either object maintenance or static imagery. mPFC theta phase coupling was not observed with any other brain region. These results implicate oscillatory coupling between mPFC and medial temporal lobe/RSc theta rhythms in the dynamic mental exploration of imagined scenes.

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.

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.

Neural Competitive Queuing of Ordinal Structure Underlies Skilled Sequential Action.

Fluent retrieval and execution of movement sequences is essential for daily activities, but the neural mechanisms underlying sequence planning remain elusive. Here participants learned finger press sequences with different orders and timings and reproduced them in a magneto-encephalography (MEG) scanner. We classified the MEG patterns for each press in the sequence and examined pattern dynamics during preparation and production. Our results demonstrate the "competitive queuing" (CQ) of upcoming action representations, extending previous computational and non-human primate recording studies to non-invasive measures in humans. In addition, we show that CQ reflects an ordinal template that generalizes across specific motor actions at each position. Finally, we demonstrate that CQ predicts participants' production accuracy and originates from parahippocampal and cerebellar sources. These results suggest that the brain learns and controls multiple sequences by flexibly combining representations of specific actions and interval timing with high-level, parallel representations of sequence position.

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.

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.

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.