Concurrent spinal and brain imaging with optically pumped magnetometers

BACKGROUND

The spinal cord and its interactions with the brain are fundamental for movement control and somatosensation. However, brain and spinal electrophysiology in humans have largely been treated as distinct enterprises, in part due to the relative inaccessibility of the spinal cord. Consequently, there is a dearth of knowledge on human spinal electrophysiology, including the multiple pathologies that affect the spinal cord as well as the brain.

NEW METHOD

Here we exploit recent advances in the development of wearable optically pumped magnetometers (OPMs) which can be flexibly arranged to provide coverage of both the spinal cord and the brain in relatively unconstrained environments. This system for magnetospinoencephalography (MSEG) measures both spinal and cortical signals simultaneously by employing custom-made scanning casts.

RESULTS

We evidence the utility of such a system by recording spinal and cortical evoked responses to median nerve stimulation at the wrist. MSEG revealed early (10 - 15 ms) and late (>20 ms) responses at the spinal cord, in addition to typical cortical evoked responses (i.e., N20).

COMPARISON WITH EXISTING METHODS

Early spinal evoked responses detected were in line with conventional somatosensory evoked potential recordings.

CONCLUSION

This MSEG system demonstrates the novel ability for concurrent non-invasive millisecond imaging of brain and spinal cord.

Adaptive multipole models of optically pumped magnetometer data

Multipole expansions have been used extensively in the Magnetoencephalography (MEG) literature for mitigating environmental interference and modelling brain signal. However, their application to Optically Pumped Magnetometer (OPM) data is challenging due to the wide variety of existing OPM sensor and array designs. We therefore explore how such multipole models can be adapted to provide stable models of brain signal and interference across OPM systems. Firstly, we demonstrate how prolate spheroidal (rather than spherical) harmonics can provide a compact representation of brain signal when sampling on the scalp surface with as few as 100 channels. We then introduce a type of orthogonal projection incorporating this basis set. The Adaptive Multipole Models (AMM), which provides robust interference rejection across systems, even in the presence of spatially structured nonlinearity errors (shielding factor is the reciprocal of the maximum fractional nonlinearity error). Furthermore, this projection is always stable, as it is an orthogonal projection, and will only ever decrease the white noise in the data. However, for array designs that are suboptimal for spatially separating brain signal and interference, this method can remove brain signal components. We contrast these properties with the more typically used multipole expansion, Signal Space Separation (SSS), which never reduces brain signal amplitude but is less robust to the effect of sensor nonlinearity errors on interference rejection and can increase noise in the data if the system is sub-optimally designed (as it is an oblique projection). We conclude with an empirical example utilizing AMM to maximize signal to noise ratio (SNR) for the stimulus locked neuronal response to a flickering visual checkerboard in a 128-channel OPM system and demonstrate up to 40 dB software shielding in real data.

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.

Real-time, model-based magnetic field correction for moving, wearable MEG

Most neuroimaging techniques require the participant to remain still for reliable recordings to be made. Optically pumped magnetometer (OPM) based magnetoencephalography (OP-MEG) however, is a neuroimaging technique which can be used to measure neural signals during large participant movement (approximately 1 m) within a magnetically shielded room (MSR) (Boto et al., 2018; Seymour et al., 2021). Nevertheless, environmental magnetic fields vary both spatially and temporally and OPMs can only operate within a limited magnetic field range, which constrains participant movement. Here we implement real-time updates to electromagnetic coils mounted on-board of the OPMs, to cancel out the changing background magnetic fields. The coil currents were chosen based on a continually updating harmonic model of the background magnetic field, effectively implementing homogeneous field correction (HFC) in real-time (Tierney et al., 2021). During a stationary, empty room recording, we show an improvement in very low frequency noise of 24 dB. In an auditory paradigm, during participant movement of up to 2 m within a magnetically shielded room, introduction of the real-time correction more than doubled the proportion of trials in which no sensor saturated recorded outside of a 50 cm radius from the optimally-shielded centre of the room. The main advantage of such model-based (rather than direct) feedback is that it could allow one to correct field components along unmeasured OPM axes, potentially mitigating sensor gain and calibration issues (Borna et al., 2022).

Non-invasive measurements of ictal and interictal epileptiform activity using optically pumped magnetometers

Magneto- and electroencephalography (MEG/EEG) are important techniques for the diagnosis and pre-surgical evaluation of epilepsy. Yet, in current cryogen-based MEG systems the sensors are offset from the scalp, which limits the signal-to-noise ratio (SNR) and thereby the sensitivity to activity from deep structures such as the hippocampus. This effect is amplified in children, for whom adult-sized fixed-helmet systems are typically too big. Moreover, ictal recordings with fixed-helmet systems are problematic because of limited movement tolerance and/or logistical considerations. Optically Pumped Magnetometers (OPMs) can be placed directly on the scalp, thereby improving SNR and enabling recordings during seizures. We aimed to demonstrate the performance of OPMs in a clinical population. Seven patients with challenging cases of epilepsy underwent MEG recordings using a 12-channel OPM-system and a 306-channel cryogen-based whole-head system: three adults with known deep or weak (low SNR) sources of interictal epileptiform discharges (IEDs), along with three children with focal epilepsy and one adult with frequent seizures. The consistency of the recorded IEDs across the two systems was assessed. In one patient the OPMs detected IEDs that were not found with the SQUID-system, and in two patients no IEDs were found with either system. For the other patients the OPM data were remarkably consistent with the data from the cryogenic system, noting that these were recorded in different sessions, with comparable SNRs and IED-yields overall. Importantly, the wearability of OPMs enabled the recording of seizure activity in a patient with hyperkinetic movements during the seizure. The observed ictal onset and semiology were in agreement with previous video- and stereo-EEG recordings. The relatively affordable technology, in combination with reduced running and maintenance costs, means that OPM-based MEG could be used more widely than current MEG systems, and may become an affordable alternative to scalp EEG, with the potential benefits of increased spatial accuracy, reduced sensitivity to volume conduction/field spread, and increased sensitivity to deep sources. Wearable MEG thus provides an unprecedented opportunity for epilepsy, and given its patient-friendliness, we envisage that it will not only be used for presurgical evaluation of epilepsy patients, but also for diagnosis after a first seizure.

Spherical harmonic based noise rejection and neuronal sampling with multi-axis OPMs

In this study we explore the interference rejection and spatial sampling properties of multi-axis Optically Pumped Magnetometer (OPM) data. We use both vector spherical harmonics and eigenspectra to quantify how well an array can separate neuronal signal from environmental interference while adequately sampling the entire cortex. We found that triaxial OPMs have superb noise rejection properties allowing for very high orders of interference (L=6) to be accounted for while minimally affecting the neural space (2dB attenuation for a 60-sensor triaxial system). We show that at least 11th order (143 spatial degrees of freedom) irregular solid harmonics or 95 eigenvectors of the lead field are needed to model the neural space for OPM data (regardless of number of axes measured). This can be adequately sampled with 75-100 equidistant triaxial sensors (225-300 channels) or 200 equidistant radial channels. In other words, ordering the same number of channels in triaxial (rather than purely radial) configuration may give significant advantages not only in terms of external noise rejection but also by minimizing cost, weight and cross-talk.

A lightweight magnetically shielded room with active shielding

Magnetically shielded rooms (MSRs) use multiple layers of materials such as MuMetal to screen external magnetic fields that would otherwise interfere with high precision magnetic field measurements such as magnetoencephalography (MEG). Optically pumped magnetometers (OPMs) have enabled the development of wearable MEG systems which have the potential to provide a motion tolerant functional brain imaging system with high spatiotemporal resolution. Despite significant promise, OPMs impose stringent magnetic shielding requirements, operating around a zero magnetic field resonance within a dynamic range of ± 5 nT. MSRs developed for OPM-MEG must therefore effectively shield external sources and provide a low remnant magnetic field inside the enclosure. Existing MSRs optimised for OPM-MEG are expensive, heavy, and difficult to site. Electromagnetic coils are used to further cancel the remnant field inside the MSR enabling participant movements during OPM-MEG, but present coil systems are challenging to engineer and occupy space in the MSR limiting participant movements and negatively impacting patient experience. Here we present a lightweight MSR design (30% reduction in weight and 40-60% reduction in external dimensions compared to a standard OPM-optimised MSR) which takes significant steps towards addressing these barriers. We also designed a 'window coil' active shielding system, featuring a series of simple rectangular coils placed directly onto the walls of the MSR. By mapping the remnant magnetic field inside the MSR, and the magnetic field produced by the coils, we can identify optimal coil currents and cancel the remnant magnetic field over the central cubic metre to just |B|= 670 ± 160 pT. These advances reduce the cost, installation time and siting restrictions of MSRs which will be essential for the widespread deployment of OPM-MEG.

Interference suppression techniques for OPM-based MEG: Opportunities and challenges

One of the primary technical challenges facing magnetoencephalography (MEG) is that the magnitude of neuromagnetic fields is several orders of magnitude lower than interfering signals. Recently, a new type of sensor has been developed - the optically pumped magnetometer (OPM). These sensors can be placed directly on the scalp and move with the head during participant movement, making them wearable. This opens up a range of exciting experimental and clinical opportunities for OPM-based MEG experiments, including paediatric studies, and the incorporation of naturalistic movements into neuroimaging paradigms. However, OPMs face some unique challenges in terms of interference suppression, especially in situations involving mobile participants, and when OPMs are integrated with electrical equipment required for naturalistic paradigms, such as motion capture systems. Here we briefly review various hardware solutions for OPM interference suppression. We then outline several signal processing strategies aimed at increasing the signal from neuromagnetic sources. These include regression-based strategies, temporal filtering and spatial filtering approaches. The focus is on the practical application of these signal processing algorithms to OPM data. In a similar vein, we include two worked-through experiments using OPM data collected from a whole-head sensor array. These tutorial-style examples illustrate how the steps for suppressing external interference can be implemented, including the associated data and code so that researchers can try the pipelines for themselves. With the popularity of OPM-based MEG rising, there will be an increasing need to deal with interference suppression. We hope this practical paper provides a resource for OPM-based MEG researchers to build upon.

Using OPMs to measure neural activity in standing, mobile participants

Optically pumped magnetometer-based magnetoencephalography (OP-MEG) can be used to measure neuromagnetic fields while participants move in a magnetically shielded room. Head movements in previous OP-MEG studies have been up to 20 cm translation and ∼30° rotation in a sitting position. While this represents a step-change over stationary MEG systems, naturalistic head movement is likely to exceed these limits, particularly when participants are standing up. In this proof-of-concept study, we sought to push the movement limits of OP-MEG even further. Using a 90 channel (45-sensor) whole-head OP-MEG system and concurrent motion capture, we recorded auditory evoked fields while participants were: (i) sitting still, (ii) standing up and still, and (iii) standing up and making large natural head movements continuously throughout the recording - maximum translation 120 cm, maximum rotation 198°. Following pre-processing, movement artefacts were substantially reduced but not eliminated. However, upon utilisation of a beamformer, the M100 event-related field localised to primary auditory regions. Furthermore, the event-related fields from auditory cortex were remarkably consistent across the three conditions. These results suggest that a wide range of movement is possible with current OP-MEG systems. This in turn underscores the exciting potential of OP-MEG for recording neural activity during naturalistic paradigms that involve movement (e.g. navigation), and for scanning populations who are difficult to study with stationary MEG (e.g. young children).

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

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

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.

The Role of Hippocampal-Ventromedial Prefrontal Cortex Neural Dynamics in Building Mental Representations

The hippocampus and ventromedial prefrontal cortex (vmPFC) play key roles in numerous cognitive domains including mind-wandering, episodic memory, and imagining the future. Perspectives differ on precisely how they support these diverse functions, but there is general agreement that it involves constructing representations composed of numerous elements. Visual scenes have been deployed extensively in cognitive neuroscience because they are paradigmatic multielement stimuli. However, it remains unclear whether scenes, rather than other types of multifeature stimuli, preferentially engage hippocampus and vmPFC. Here, we leveraged the high temporal resolution of magnetoencephalography to test participants as they gradually built scene imagery from three successive auditorily presented object descriptions and an imagined 3-D space. This was contrasted with constructing mental images of nonscene arrays that were composed of three objects and an imagined 2-D space. The scene and array stimuli were, therefore, highly matched, and this paradigm permitted a closer examination of step-by-step mental construction than has been undertaken previously. We observed modulation of theta power in our two regions of interest-anterior hippocampus during the initial stage and vmPFC during the first two stages, of scene relative to array construction. Moreover, the scene-specific anterior hippocampal activity during the first construction stage was driven by the vmPFC, with mutual entrainment between the two brain regions thereafter. These findings suggest that hippocampal and vmPFC neural activity is especially tuned to scene representations during the earliest stage of their formation, with implications for theories of how these brain areas enable cognitive functions such as episodic memory.

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.

The Effect of Object Type on Building Scene Imagery-an MEG Study

Previous studies have reported that some objects evoke a sense of local three-dimensional space (space-defining; SD), while others do not (space-ambiguous; SA), despite being imagined or viewed in isolation devoid of a background context. Moreover, people show a strong preference for SD objects when given a choice of objects with which to mentally construct scene imagery. When deconstructing scenes, people retain significantly more SD objects than SA objects. It, therefore, seems that SD objects might enjoy a privileged role in scene construction. In the current study, we leveraged the high temporal resolution of magnetoencephalography (MEG) to compare the neural responses to SD and SA objects while they were being used to build imagined scene representations, as this has not been examined before using neuroimaging. On each trial, participants gradually built a scene image from three successive auditorily-presented object descriptions and an imagined 3D space. We then examined the neural dynamics associated with the points during scene construction when either SD or SA objects were being imagined. We found that SD objects elicited theta changes relative to SA objects in two brain regions, the right ventromedial prefrontal cortex (vmPFC) and the right superior temporal gyrus (STG). Furthermore, using dynamic causal modeling, we observed that the vmPFC drove STG activity. These findings may indicate that SD objects serve to activate schematic and conceptual knowledge in vmPFC and STG upon which scene representations are then built.

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.

vmPFC Drives Hippocampal Processing during Autobiographical Memory Recall Regardless of Remoteness

Our ability to recall past experiences, autobiographical memories (AMs), is crucial to cognition, endowing us with a sense of self and underwriting our capacity for autonomy. Traditional views assume that the hippocampus orchestrates event recall, whereas recent accounts propose that the ventromedial prefrontal cortex (vmPFC) instigates and coordinates hippocampal-dependent processes. Here we sought to characterize the dynamic interplay between the hippocampus and vmPFC during AM recall to adjudicate between these perspectives. Leveraging the high temporal resolution of magnetoencephalography, we found that the left hippocampus and the vmPFC showed the greatest power changes during AM retrieval. Moreover, responses in the vmPFC preceded activity in the hippocampus during initiation of AM recall, except during retrieval of the most recent AMs. The vmPFC drove hippocampal activity during recall initiation and also as AMs unfolded over subsequent seconds, and this effect was evident regardless of AM age. These results recast the positions of the hippocampus and the vmPFC in the AM retrieval hierarchy, with implications for theoretical accounts of memory processing and systems-level consolidation.

Representation of probabilistic outcomes during risky decision-making

Goal-directed behaviour requires prospectively retrieving and evaluating multiple possible action outcomes. While a plethora of studies suggested sequential retrieval for deterministic choice outcomes, it remains unclear whether this is also the case when integrating multiple probabilistic outcomes of the same action. We address this question by capitalising on magnetoencephalography (MEG) in humans who made choices in a risky foraging task. We train classifiers to distinguish MEG field patterns during presentation of two probabilistic outcomes (reward, loss), and then apply these to decode such patterns during deliberation. First, decoded outcome representations have a temporal structure, suggesting alternating retrieval of the outcomes. Moreover, the probability that one or the other outcome is being represented depends on loss magnitude, but not on loss probability, and it predicts the chosen action. In summary, we demonstrate decodable outcome representations during probabilistic decision-making, which are sequentially structured, depend on task features, and predict subsequent action.

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.

Optically pumped magnetoencephalography in epilepsy

We demonstrate the first use of Optically Pumped Magnetoencephalography (OP-MEG) in an epilepsy patient with unrestricted head movement. Current clinical MEG uses a traditional SQUID system, where sensors are cryogenically cooled and housed in a helmet in which the patient's head is fixed. Here, we use a different type of sensor (OPM), which operates at room temperature and can be placed directly on the patient's scalp, permitting free head movement. We performed OP-MEG recording in a patient with refractory focal epilepsy. OP-MEG-identified analogous interictal activity to scalp EEG, and source localized this activity to an appropriate brain region.

Imaging the human hippocampus with optically-pumped magnetoencephalography

Optically-pumped (OP) magnetometers allow magnetoencephalography (MEG) to be performed while a participant's head is unconstrained. To fully leverage this new technology, and in particular its capacity for mobility, the activity of deep brain structures which facilitate explorative behaviours such as navigation, must be detectable using OP-MEG. One such crucial brain region is the hippocampus. Here we had three healthy adult participants perform a hippocampal-dependent task - the imagination of novel scene imagery - while being scanned using OP-MEG. A conjunction analysis across these three participants revealed a significant change in theta power in the medial temporal lobe. The peak of this activated cluster was located in the anterior hippocampus. We repeated the experiment with the same participants in a conventional SQUID-MEG scanner and found similar engagement of the medial temporal lobe, also with a peak in the anterior hippocampus. These OP-MEG findings indicate exciting new opportunities for investigating the neural correlates of a range of crucial cognitive functions in naturalistic contexts including spatial navigation, episodic memory and social interactions.

A tool for functional brain imaging with lifespan compliance

The human brain undergoes significant functional and structural changes in the first decades of life, as the foundations for human cognition are laid down. However, non-invasive imaging techniques to investigate brain function throughout neurodevelopment are limited due to growth in head-size with age and substantial head movement in young participants. Experimental designs to probe brain function are also limited by the unnatural environment typical brain imaging systems impose. However, developments in quantum technology allowed fabrication of a new generation of wearable magnetoencephalography (MEG) technology with the potential to revolutionise electrophysiological measures of brain activity. Here we demonstrate a lifespan-compliant MEG system, showing recordings of high fidelity data in toddlers, young children, teenagers and adults. We show how this system can support new types of experimental paradigm involving naturalistic learning. This work reveals a new approach to functional imaging, providing a robust platform for investigation of neurodevelopment in health and disease.

Magnetoencephalography (MEG) recordings are sensitive to movement and therefore are especially challenging with young participants. Here the authors develop a wearable MEG system based on a modified bicycle helmet, which enables reliable recordings in toddlers, children, teenagers and adults.

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.

Balanced, bi-planar magnetic field and field gradient coils for field compensation in wearable magnetoencephalography

To allow wearable magnetoencephalography (MEG) recordings to be made on unconstrained subjects the spatially inhomogeneous remnant magnetic field inside the magnetically shielded room (MSR) must be nulled. Previously, a large bi-planar coil system which produces uniform fields and field gradients was used for this purpose. Its construction presented a significant challenge, six distinct coils were wound on two 1.6 × 1.6 m2 planes. Here, we exploit shared coil symmetries to produce coils simultaneously optimised to generate homogenous fields and gradients. We show nulling performance comparable to that of a six-coil system is achieved with this three-coil system, decreasing the strongest field component Bx by a factor of 53, and the strongest gradient dBx/dz by a factor of 7. To allow the coils to be used in environments with temporally-varying magnetic interference a dynamic nulling system was developed with a shielding factor of 40 dB at 0.01 Hz. Reducing the number of coils required and incorporating dynamic nulling should allow for greater take-up of this technology. Interactions of the coils with the high-permeability walls of the MSR were investigated using a method of images approach. Simulations show a degrading of field uniformity which was broadly consistent with measured values. These effects should be incorporated into future designs.

Towards OPM-MEG in a virtual reality environment

Virtual reality (VR) provides an immersive environment in which a participant can experience a feeling of presence in a virtual world. Such environments generate strong emotional and physical responses and have been used for wide-ranging applications. The ability to collect functional neuroimaging data whilst a participant is immersed in VR would represent a step change for experimental paradigms; unfortunately, traditional brain imaging requires participants to remain still, limiting the scope of naturalistic interaction within VR. Recently however, a new type of magnetoencephalography (MEG) device has been developed, that employs scalp-mounted optically-pumped magnetometers (OPMs) to measure brain electrophysiology. Lightweight OPMs, coupled with precise control of the background magnetic field, enables participant movement during data acquisition. Here, we exploit this technology to acquire MEG data whilst a participant uses a virtual reality head-mounted display (VRHMD). We show that, despite increased magnetic interference from the VRHMD, we were able to measure modulation of alpha-band oscillations, and the visual evoked field. Moreover, in a VR experiment in which a participant had to move their head to look around a virtual wall and view a visual stimulus, we showed that the measured MEG signals map spatially in accordance with the known organisation of primary visual cortex. This technique could transform the type of neuroscientific experiment that can be undertaken using functional neuroimaging.

Optically pumped magnetometers: From quantum origins to multi-channel magnetoencephalography

Optically Pumped Magnetometers (OPMs) have emerged as a viable and wearable alternative to cryogenic, superconducting MEG systems. This new generation of sensors has the advantage of not requiring cryogenic cooling and as a result can be flexibly placed on any part of the body. The purpose of this review is to provide a neuroscience audience with the theoretical background needed to understand the physical basis for the signal observed by OPMs. Those already familiar with the physics of MRI and NMR should note that OPMs share much of the same theory as the operation of OPMs rely on magnetic resonance. This review establishes the physical basis for the signal equation for OPMs. We re-derive the equations defining the bounds on OPM performance and highlight the important trade-offs between quantities such as bandwidth, sensor size and sensitivity. These equations lead to a direct upper bound on the gain change due to cross-talk for a multi-channel OPM system.

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.

Using optically pumped magnetometers to measure magnetoencephalographic signals in the human cerebellum

Key points

The application of conventional cryogenic magnetoencephalography (MEG) to the study of cerebellar functions is highly limited because typical cryogenic sensor arrays are far away from the cerebellum and naturalistic movement is not allowed in the recording.

A new generation of MEG using optically pumped magnetometers (OPMs) that can be worn on the head during movement has opened up an opportunity to image the cerebellar electrophysiological activity non‐invasively.

We use OPMs to record human cerebellar MEG signals elicited by air‐puff stimulation to the eye.

We demonstrate robust responses in the cerebellum.

OPMs pave the way for studying the neurophysiology of the human cerebellum.

Abstract

We test the feasibility of an optically pumped magnetometer‐based magnetoencephalographic (OP‐MEG) system for the measurement of human cerebellar activity. This is to our knowledge the first study investigating the human cerebellar electrophysiology using optically pumped magnetometers. As a proof of principle, we use an air‐puff stimulus to the eyeball in order to elicit cerebellar activity that is well characterized in non‐human models. In three subjects, we observe an evoked component at approx. 50 ms post‐stimulus, followed by a second component at approx. 85–115 ms post‐stimulus. Source inversion localizes both components in the cerebellum, while control experiments exclude potential sources elsewhere. We also assess the induced oscillations, with time‐frequency decompositions, and identify additional sources in the occipital lobe, a region expected to be active in our paradigm, and in the neck muscles. Neither of these contributes to the stimulus‐evoked responses at 50–115 ms. We conclude that OP‐MEG technology offers a promising way to advance the understanding of the information processing mechanisms in the human cerebellum.

The application of conventional cryogenic magnetoencephalography (MEG) to the study of cerebellar functions is highly limited because typical cryogenic sensor arrays are far away from the cerebellum and naturalistic movement is not allowed in the recording.

A new generation of MEG using optically pumped magnetometers (OPMs) that can be worn on the head during movement has opened up an opportunity to image the cerebellar electrophysiological activity non‐invasively.

We use OPMs to record human cerebellar MEG signals elicited by air‐puff stimulation to the eye.

We demonstrate robust responses in the cerebellum.

OPMs pave the way for studying the neurophysiology of the human cerebellum.

A bi-planar coil system for nulling background magnetic fields in scalp mounted magnetoencephalography

Small, commercially-available Optically Pumped Magnetometers (OPMs) can be used to construct a wearable Magnetoencephalography (MEG) system that allows large head movements to be made during recording. The small dynamic range of these sensors however means that movement in the residual static magnetic field found inside typical Magnetically Shielded Rooms (MSRs) can saturate the sensor outputs, rendering the data unusable. This problem can be ameliorated by using a set of electromagnetic coils to attenuate the spatially-varying remnant field. Here, an array of bi-planar coils, which produce an open and accessible scanning environment, was designed and constructed. The coils were designed using a harmonic minimisation method previously used for gradient coil design in Magnetic Resonance Imaging (MRI). Six coils were constructed to null Bx, By and Bz as well as the three dominant field gradients dBx/dz, dBy/dz and dBz/dz. The coils produce homogeneous (within ±5%) fields or field gradients over a volume of 40 × 40 × 40 cm3. This volume is sufficient to contain an array of OPMs, mounted in a 3D-printed scanner-cast, during basic and natural movements. Automated control of the coils using reference sensor measurements allows reduction of the largest component of the static field (Bx) from 21.8 ± 0.2 nT to 0.47 ± 0.08 nT. The largest gradient (dBx/dz) was reduced from 7.4 nT/m to 0.55 nT/m. High precision optical tracking allowed experiments involving controlled and measured head movements, which revealed that a rotation of the scanner-cast by ±34° and translation of ±9.7 cm of the OPMs in this field generated only a 1 nT magnetic field variation across the OPM array, when field nulling was applied. This variation could be further reduced to 0.04 nT by linear regression of field variations that were correlated with the measured motion parameters. To demonstrate the effectiveness of the bi-planar coil field cancellation system in a real MEG experiment, a novel measurement of retinotopy was investigated, where the stimulus remains fixed and head movements made by the subject shift the visual presentation to the lower left or right quadrants of the field of view. Left and right visual field stimulation produced the expected responses in the opposing hemisphere. This simple demonstration shows that the bi-planar coil system allows accurate OPM-MEG recordings to be made on an unrestrained subject.

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.

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.

Magnetoencephalographic Correlates of Perceptual State During Auditory Bistability

Bistability occurs when two alternative percepts can be derived from the same physical stimulus. To identify the neural correlates of specific subjective experiences we used a bistable auditory stimulus and determined whether the two perceptual states could be distinguished electrophysiologically. Fourteen participants underwent magnetoencephalography while reporting their perceptual experience while listening to a continuous bistable stream of auditory tones. Participants reported bistability with a similar overall proportion of the two alternative percepts (52% vs 48%). At the individual level, sensor space electrophysiological discrimination between the percepts was possible in 9/14 participants with canonical variate analysis (CVA) or linear support vector machine (SVM) analysis over space and time dimensions. Classification was possible in 14/14 subjects with non-linear SVM. Similar effects were noted in an unconstrained source space CVA analysis (classifying 10/14 participants), linear SVM (classifying 9/14 subjects) and non-linear SVM (classifiying 13/14 participants). Source space analysis restricted to a priori ROIs showed discrimination was possible in the right and left auditory cortex with each classification approach but in the right intraparietal sulcus this was only apparent with non-linear SVM and only in a minority of particpants. Magnetoencephalography can be used to objectively classify auditory experiences from individual subjects.

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

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

Section

Analysis methods.

Classifications

Source localization: inverse problem; Source localization: other.

•

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

•

Source inversion algorithms with sparsity constraints performed best.

•

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

Reconstructing anatomy from electro-physiological data

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

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.

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.

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.

Cortical surface reconstruction based on MEG data and spherical harmonics

Estimates of coefficients of a spherical harmonic Fourier decomposition of the cortical surface can be obtained solely using MEG/EEG data and free energy as objective function. A stochastic methodology based on a Metropolis Search followed by a Bayesian Model Averaging is proposed to reconstruct cortical anatomy based functional information.

Ventromedial prefrontal cortex drives hippocampal theta oscillations induced by mismatch computations

Detecting environmental change is fundamental for adaptive behavior in an uncertain world. Previous work indicates the hippocampus supports the generation of novelty signals via implementation of a match–mismatch detector that signals when an incoming sensory input violates expectations based on past experience. While existing work has emphasized the particular contribution of the hippocampus, here we ask which other brain structures also contribute to match–mismatch detection. Furthermore, we leverage the fine-grained temporal resolution of magnetoencephalography (MEG) to investigate whether mismatch computations are spectrally confined to the theta range, based on the prominence of this range of oscillations in models of hippocampal function. By recording MEG activity while human subjects perform a task that incorporates conditions of match–mismatch novelty we show that mismatch signals are confined to the theta band and are expressed in both the hippocampus and ventromedial prefrontal cortex (vmPFC). Effective connectivity analyses (dynamic causal modeling) show that the hippocampus and vmPFC work as a functional circuit during mismatch detection. Surprisingly, our results suggest that the vmPFC drives the hippocampus during the generation and processing of mismatch signals. Our findings provide new evidence that the hippocampal–vmPFC circuit is engaged during novelty processing, which has implications for emerging theories regarding the role of vmPFC in memory.

•

Mismatch detection engages human hippocampus and ventromedial prefrontal cortex.

•

Novelty signals are spectrally confined to the theta band.

•

Ventromedial prefrontal cortex drives hippocampal theta induced by mismatches.

Resting state MEG oscillations show long-range temporal correlations of phase synchrony that break down during finger movement

The capacity of the human brain to interpret and respond to multiple temporal scales in its surroundings suggests that its internal interactions must also be able to operate over a broad temporal range. In this paper, we utilize a recently introduced method for characterizing the rate of change of the phase difference between MEG signals and use it to study the temporal structure of the phase interactions between MEG recordings from the left and right motor cortices during rest and during a finger-tapping task. We use the Hilbert transform to estimate moment-to-moment fluctuations of the phase difference between signals. After confirming the presence of scale-invariance we estimate the Hurst exponent using detrended fluctuation analysis (DFA). An exponent of >0.5 is indicative of long-range temporal correlations (LRTCs) in the signal. We find that LRTCs are present in the α/μ and β frequency bands of resting state MEG data. We demonstrate that finger movement disrupts LRTCs correlations, producing a phase relationship with a structure similar to that of Gaussian white noise. The results are validated by applying the same analysis to data with Gaussian white noise phase difference, recordings from an empty scanner and phase-shuffled time series. We interpret the findings through comparison of the results with those we obtained from an earlier study during which we adopted this method to characterize phase relationships within a Kuramoto model of oscillators in its sub-critical, critical, and super-critical synchronization states. We find that the resting state MEG from left and right motor cortices shows moment-to-moment fluctuations of phase difference with a similar temporal structure to that of a system of Kuramoto oscillators just prior to its critical level of coupling, and that finger tapping moves the system away from this pre-critical state toward a more random state.

Dynamic recruitment of resting state sub-networks

Resting state networks (RSNs) are of fundamental importance in human systems neuroscience with evidence suggesting that they are integral to healthy brain function and perturbed in pathology. Despite rapid progress in this area, the temporal dynamics governing the functional connectivities that underlie RSN structure remain poorly understood. Here, we present a framework to help further our understanding of RSN dynamics. We describe a methodology which exploits the direct nature and high temporal resolution of magnetoencephalography (MEG). This technique, which builds on previous work, extends from solving fundamental confounds in MEG (source leakage) to multivariate modelling of transient connectivity. The resulting processing pipeline facilitates direct (electrophysiological) measurement of dynamic functional networks. Our results show that, when functional connectivity is assessed in small time windows, the canonical sensorimotor network can be decomposed into a number of transiently synchronising sub-networks, recruitment of which depends on current mental state. These rapidly changing sub-networks are spatially focal with, for example, bilateral primary sensory and motor areas resolved into two separate sub-networks. The likely interpretation is that the larger canonical sensorimotor network most often seen in neuroimaging studies reflects only a temporal aggregate of these transient sub-networks. Our approach opens new frontiers to study RSN dynamics, showing that MEG is capable of revealing the spatial, temporal and spectral signature of the human connectome in health and disease.

•

The sensorimotor network consists of a series of transiently synchronising subnetworks.

•

These subnetworks are robust across multiple tasks.

•

The occurrence of these subnetworks is modulated by the current mental state.

Head-eye coordination in normals and in patients with vestibular disorders

Oculomotor response in the absence of vision was examined in 8 normal subjects, 11 unilaterally labyrinthectomised patients and 2 patients with suspected bilateral canal paresis. The experiments involved (a) voluntary oscillation of the head, (b) whole body oscillation on a turntable and (c) stimulation of neck afferents by oscillation of the body with the head fixed. In the patients with unilateral lesions there was a directional preponderance of the slow phase eye velocity towards the side of the lesion which differed significantly from that of the normal population. In the patients with bilateral paresis the oculomotor response to whole body oscillation was negligible, whereas the response to voluntary head movement had a mean gain of 0.45 and at high frequency could not be suppressed when viewing a head-fixed image. The saccadic activity during voluntary head movement was similar in all subjects and was correlated with slow phase velocity.

Structure predicts function: combining non-invasive electrophysiology with in-vivo histology

We present an approach for combining high resolution MRI-based myelin mapping with functional information from electroencephalography (EEG) or magnetoencephalography (MEG). The main contribution to the primary currents detectable with EEG and MEG comes from ionic currents in the apical dendrites of cortical pyramidal cells, aligned perpendicularly to the local cortical surface. We provide evidence from an in-vivo experiment that the variation in MRI-based myeloarchitecture measures across the cortex predicts the variation of the current density over individuals and thus is of functional relevance. Equivalent current dipole locations and moments due to pitch onset evoked response fields (ERFs) were estimated by means of a variational Bayesian algorithm. The myeloarchitecture was estimated indirectly from individual high resolution quantitative multi-parameter maps (MPMs) acquired at 800μm isotropic resolution. Myelin estimates across cortical areas correlated positively with dipole magnitude. This correlation was spatially specific: regions of interest in the auditory cortex provided significantly better models than those covering whole hemispheres. Based on the MPM data we identified the auditory cortical area TE1.2 as the most likely origin of the pitch ERFs measured by MEG. We can now proceed to exploit the higher spatial resolution of quantitative MPMs to identify the cortical origin of M/EEG signals, inform M/EEG source reconstruction and explore structure-function relationships at a fine structural level in the living human brain.

Between Thought and Expression, a Magnetoencephalography Study of the “Tip-of-the-Tongue” Phenomenon

“Tip-of-the-tongue” (TOT) is the phenomenon associated with the inaccessibility of a known word from memory. It is universally experienced, increases in frequency with age, and is most common for proper nouns. It is a good model for the symptom of anomia experienced much more frequently by some aphasic patients following brain injury. Here, we induced the TOT state in older participants while they underwent brain scanning with magnetoencephalography to investigate the changes in oscillatory brain activity associated with failed retrieval of known words. Using confrontation naming of pictures of celebrities, we successfully induced the TOT state in 29% of trials and contrasted it with two other states: “Know” where the participants both correctly recognized the celebrity's face and retrieved their name and “Don't Know” when the participants did not recognize the celebrity. We wished to test Levelt's influential model of speech output by carrying out two analyses, one epoching the data to the point in time when the picture was displayed and the other looking back in time from when the participants first articulated their responses. Our main findings supported the components of Levelt's model, but not their serial activation over time as both semantic and motor areas were identified in both analyses. We also found enduring decreases in the alpha frequency band in the left ventral temporal region during the TOT state, suggesting ongoing semantic search. Finally, we identified reduced beta power in classical peri-sylvian language areas for the TOT condition, suggesting that brain regions that encode linguistic memories are also involved in their attempted retrieval.

Modulation of alpha and gamma oscillations related to retrospectively orienting attention within working memory

Selective attention mechanisms allow us to focus on information that is relevant to the current behavior and, equally important, ignore irrelevant information. An influential model proposes that oscillatory neural activity in the alpha band serves as an active functional inhibitory mechanism. Recent studies have shown that, in the same way that attention can be selectively oriented to bias sensory processing in favor of relevant stimuli in perceptual tasks, it is also possible to retrospectively orient attention to internal representations held in working memory. However, these studies have not explored the associated oscillatory phenomena. In the current study, we analysed the patterns of neural oscillatory activity recorded with magnetoencephalography while participants performed a change detection task, in which a spatial retro-cue was presented during the maintenance period, indicating which item or items were relevant for subsequent retrieval. Participants benefited from retro-cues in terms of accuracy and reaction time. Retro-cues also modulated oscillatory activity in the alpha and gamma frequency bands. We observed greater alpha activity in a ventral visual region ipsilateral to the attended hemifield, thus supporting its suppressive role, i.e. a functional disengagement of task-irrelevant regions. Accompanying this modulation, we found an increase in gamma activity contralateral to the attended hemifield, which could reflect attentional orienting and selective processing. These findings suggest that the oscillatory mechanisms underlying attentional orienting to representations held in working memory are similar to those engaged when attention is oriented in the perceptual space.

Medial prefrontal theta phase coupling during spatial memory retrieval

Memory retrieval is believed to involve a disparate network of areas, including medial prefrontal and medial temporal cortices, but the mechanisms underlying their coordination remain elusive. One suggestion is that oscillatory coherence mediates inter-regional communication, implicating theta phase and theta-gamma phase-amplitude coupling in mnemonic function across species. To examine this hypothesis, we used non-invasive whole-head magnetoencephalography (MEG) as participants retrieved the location of objects encountered within a virtual environment. We demonstrate that, when participants are cued with the image of an object whose location they must subsequently navigate to, there is a significant increase in 4–8 Hz theta power in medial prefrontal cortex (mPFC), and the phase of this oscillation is coupled both with ongoing theta phase in the medial temporal lobe (MTL) and perceptually induced 65–85 Hz gamma amplitude in medial parietal cortex. These results suggest that theta phase coupling between mPFC and MTL and theta-gamma phase-amplitude coupling between mPFC and neocortical regions may play a role in human spatial memory retrieval. © 2014 The Authors. Hippocampus Published by Wiley Periodicals, Inc.

Algorithmic procedures for Bayesian MEG/EEG source reconstruction in SPM

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