Flexible head-casts for high spatial precision MEG

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.

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).

Diverse beta burst waveform motifs characterize movement-related cortical dynamics

Classical analyses of induced, frequency-specific neural activity typically average band-limited power over trials. More recently, it has become widely appreciated that in individual trials, beta band activity occurs as transient bursts rather than amplitude-modulated oscillations. Most studies of beta bursts treat them as unitary, and having a stereotyped waveform. However, we show there is a wide diversity of burst shapes. Using a biophysical model of burst generation, we demonstrate that waveform variability is predicted by variability in the synaptic drives that generate beta bursts. We then use a novel, adaptive burst detection algorithm to identify bursts from human MEG sensor data recorded during a joystick-based reaching task, and apply principal component analysis to burst waveforms to define a set of dimensions, or motifs, that best explain waveform variance. Finally, we show that bursts with a particular range of waveform motifs, ones not fully accounted for by the biophysical model, differentially contribute to movement-related beta dynamics. Sensorimotor beta bursts are therefore not homogeneous events and likely reflect distinct computational processes.

Spatiotemporal organisation of human sensorimotor beta burst activity

Beta oscillations in human sensorimotor cortex are hallmark signatures of healthy and pathological movement. In single trials, beta oscillations include bursts of intermittent, transient periods of high-power activity. These burst events have been linked to a range of sensory and motor processes, but their precise spatial, spectral, and temporal structure remains unclear. Specifically, a role for beta burst activity in information coding and communication suggests spatiotemporal patterns, or travelling wave activity, along specific anatomical gradients. We here show in human magnetoencephalography recordings that burst activity in sensorimotor cortex occurs in planar spatiotemporal wave-like patterns that dominate along two axes either parallel or perpendicular to the central sulcus. Moreover, we find that the two propagation directions are characterised by distinct anatomical and physiological features. Finally, our results suggest that sensorimotor beta bursts occurring before and after a movement can be distinguished by their anatomical, spectral, and spatiotemporal characteristics, indicating distinct functional roles.

The hearing hippocampus

The hippocampus has a well-established role in spatial and episodic memory but a broader function has been proposed including aspects of perception and relational processing. Neural bases of sound analysis have been described in the pathway to auditory cortex, but wider networks supporting auditory cognition are still being established. We review what is known about the role of the hippocampus in processing auditory information, and how the hippocampus itself is shaped by sound. In examining imaging, recording, and lesion studies in species from rodents to humans, we uncover a hierarchy of hippocampal responses to sound including during passive exposure, active listening, and the learning of associations between sounds and other stimuli. We describe how the hippocampus' connectivity and computational architecture allow it to track and manipulate auditory information - whether in the form of speech, music, or environmental, emotional, or phantom sounds. Functional and structural correlates of auditory experience are also identified. The extent of auditory-hippocampal interactions is consistent with the view that the hippocampus makes broad contributions to perception and cognition, beyond spatial and episodic memory. More deeply understanding these interactions may unlock applications including entraining hippocampal rhythms to support cognition, and intervening in links between hearing loss and dementia.

Modified MRI Anonymization (De-Facing) for Improved MEG Coregistration

Localising the sources of MEG/EEG signals often requires a structural MRI to create a head model, while ensuring reproducible scientific results requires sharing data and code. However, sharing structural MRI data often requires the face go be hidden to help protect the identity of the individuals concerned. While automated de-facing methods exist, they tend to remove the whole face, which can impair methods for coregistering the MRI data with the EEG/MEG data. We show that a new, automated de-facing method that retains the nose maintains good MRI-MEG/EEG coregistration. Importantly, behavioural data show that this "face-trimming" method does not increase levels of identification relative to a standard de-facing approach and has less effect on the automated segmentation and surface extraction sometimes used to create head models for MEG/EEG localisation. We suggest that this trimming approach could be employed for future sharing of structural MRI data, at least for those to be used in forward modelling (source reconstruction) of EEG/MEG data.

Brainprints: identifying individuals from magnetoencephalograms

Magnetoencephalography (MEG) is used to study a wide variety of cognitive processes. Increasingly, researchers are adopting principles of open science and releasing their MEG data. While essential for reproducibility, sharing MEG data has unforeseen privacy risks. Individual differences may make a participant identifiable from their anonymized recordings. However, our ability to identify individuals based on these individual differences has not yet been assessed. Here, we propose interpretable MEG features to characterize individual difference. We term these features brainprints (brain fingerprints). We show through several datasets that brainprints accurately identify individuals across days, tasks, and even between MEG and Electroencephalography (EEG). Furthermore, we identify consistent brainprint components that are important for identification. We study the dependence of identifiability on the amount of data available. We also relate identifiability to the level of preprocessing and the experimental task. Our findings reveal specific aspects of individual variability in MEG. They also raise concerns about unregulated sharing of brain data, even if anonymized.

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

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

A 10-hour within-participant magnetoencephalography narrative dataset to test models of language comprehension

Recently, cognitive neuroscientists have increasingly studied the brain responses to narratives. At the same time, we are witnessing exciting developments in natural language processing where large-scale neural network models can be used to instantiate cognitive hypotheses in narrative processing. Yet, they learn from text alone and we lack ways of incorporating biological constraints during training. To mitigate this gap, we provide a narrative comprehension magnetoencephalography (MEG) data resource that can be used to train neural network models directly on brain data. We recorded from 3 participants, 10 separate recording hour-long sessions each, while they listened to audiobooks in English. After story listening, participants answered short questions about their experience. To minimize head movement, the participants wore MEG-compatible head casts, which immobilized their head position during recording. We report a basic evoked-response analysis showing that the responses accurately localize to primary auditory areas. The responses are robust and conserved across 10 sessions for every participant. We also provide usage notes and briefly outline possible future uses of the resource.

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

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

MaDoPO: Magnetic Detection of Positions and Orientations of Segmented Deep-Brain Stimulation Electrodes: A Radiation-Free Method Based on Magnetoencephalography

Background: Current approaches to detect the positions and orientations of directional deep-brain stimulation (DBS) electrodes rely on radiative imaging data. In this study, we aim to present an improved version of a radiation-free method for magnetic detection of the position and the orientation (MaDoPO) of directional electrodes based on a series of magnetoencephalography (MEG) measurements and a possible future solution for optimized results using emerging on-scalp MEG systems. Methods: A directional DBS system was positioned into a realistic head–torso phantom and placed in the MEG scanner. A total of 24 measurements of 180 s each were performed with different predefined electrode configurations. Finite element modeling and model fitting were used to determine the position and orientation of the electrode in the phantom. Related measurements were fitted simultaneously, constraining solutions to the a priori known geometry of the electrode. Results were compared with the results of the high-quality CT imaging of the phantom. Results: The accuracy in electrode localization and orientation detection depended on the number of combined measurements. The localization error was minimized to 2.02 mm by considering six measurements with different non-directional bipolar electrode configurations. Another six measurements with directional bipolar stimulations minimized the orientation error to 4°. These values are mainly limited due to the spatial resolution of the MEG. Moreover, accuracies were investigated as a function of measurement time, number of sensors, and measurement direction of the sensors in order to define an optimized MEG device for this application. Conclusion: Although MEG introduces inaccuracies in the detection of the position and orientation of the electrode, these can be accepted when evaluating the benefits of a radiation-free method. Inaccuracies can be further reduced by the use of on-scalp MEG sensor arrays, which may find their way into clinics in the foreseeable future.

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).

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.

Can EEG and MEG detect signals from the human cerebellum?

The cerebellum plays a key role in the regulation of motor learning, coordination and timing, and has been implicated in sensory and cognitive processes as well. However, our current knowledge of its electrophysiological mechanisms comes primarily from direct recordings in animals, as investigations into cerebellar function in humans have instead predominantly relied on lesion, haemodynamic and metabolic imaging studies. While the latter provide fundamental insights into the contribution of the cerebellum to various cerebellar-cortical pathways mediating behaviour, they remain limited in terms of temporal and spectral resolution. In principle, this shortcoming could be overcome by monitoring the cerebellum's electrophysiological signals. Non-invasive assessment of cerebellar electrophysiology in humans, however, is hampered by the limited spatial resolution of electroencephalography (EEG) and magnetoencephalography (MEG) in subcortical structures, i.e., deep sources. Furthermore, it has been argued that the anatomical configuration of the cerebellum leads to signal cancellation in MEG and EEG. Yet, claims that MEG and EEG are unable to detect cerebellar activity have been challenged by an increasing number of studies over the last decade. Here we address this controversy and survey reports in which electrophysiological signals were successfully recorded from the human cerebellum. We argue that the detection of cerebellum activity non-invasively with MEG and EEG is indeed possible and can be enhanced with appropriate methods, in particular using connectivity analysis in source space. We provide illustrative examples of cerebellar activity detected with MEG and EEG. Furthermore, we propose practical guidelines to optimize the detection of cerebellar activity with MEG and EEG. Finally, we discuss MEG and EEG signal contamination that may lead to localizing spurious sources in the cerebellum and suggest ways of handling such artefacts. This review is to be read as a perspective review that highlights that it is indeed possible to measure cerebellum with MEG and EEG and encourages MEG and EEG researchers to do so. Its added value beyond highlighting and encouraging is that it offers useful advice for researchers aspiring to investigate the cerebellum with MEG and EEG.

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.

Prospects for Future Methodological Development and Application of Magnetoencephalography Devices in Psychiatry

Magnetoencephalography (MEG) is a functional neuroimaging tool that can record activity from the entire cortex on the order of milliseconds. MEG has been used to investigate numerous psychiatric disorders, such as schizophrenia, bipolar disorder, major depression, dementia, and autism spectrum disorder. Although several review papers on the subject have been published, perspectives and opinions regarding the use of MEG in psychiatric research have primarily been discussed from a psychiatric research point of view. Owing to a newly developed MEG sensor, the use of MEG devices will soon enter a critical period, and now is a good time to discuss the future of MEG use in psychiatric research. In this paper, we will discuss MEG devices from a methodological point of view. We will first introduce the utilization of MEG in psychiatric research and the development of its technology. Then, we will describe the principle theory of MEG and common algorithms, which are useful for applying MEG tools to psychiatric research. Next, we will consider three topics-child psychiatry, resting-state networks, and cortico-subcortical networks-and address the future use of MEG in psychiatry from a broader perspective. Finally, we will introduce the newly developed device, the optically-pumped magnetometer, and discuss its future use in MEG systems in psychiatric research from a methodological point of view. We believe that state-of-the-art electrophysiological tools, such as this new MEG system, will further contribute to our understanding of the core pathology in various psychiatric disorders and translational research.

A decade of infant neuroimaging research: What have we learned and where are we going?

The past decade has seen the emergence of neuroimaging studies of infant populations. Incorporating imaging has resulted in invaluable insights about neurodevelopment at the start of life. However, little has been enquired of the experimental specifications and study characteristics of typical findings. This review systematically screened empirical studies that used electroencephalography (EEG), magnetoencephalography (MEG), functional near-infrared spectroscopy (fNIRS), and functional magnetic resonance imaging (fMRI) on infants (max. age of 24 months). From more than 21,000 publications, a total of 710 records were included for analyses. With the exception of EEG studies, infant studies with MEG, fNIRS, and fMRI were most often conducted around birth and at 12 months. The vast majority of infant studies came from North America, with very few studies conducted in Africa, certain parts of South America, and Southeast Asia. Finally, longitudinal neuroimaging studies were inclined to adopt EEG, followed by fMRI, fNIRS, and MEG. These results show that there is compelling need for studies with larger sample sizes, studies investigating a broader range of infant developmental periods, and studies from under- and less-developed regions in the world. Addressing these shortcomings in the future will provide a more representative and accurate understanding of neurodevelopment in infancy.

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.

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.

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

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

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

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

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.

Imaging human cortical responses to intraneural microstimulation using magnetoencephalography

The sensation of touch in the glabrous skin of the human hand is conveyed by thousands of fast-conducting mechanoreceptive afferents, which can be categorised into four distinct types. The spiking properties of these afferents in the periphery in response to varied tactile stimuli are well-characterised, but relatively little is known about the spatiotemporal properties of the neural representations of these different receptor types in the human cortex. Here, we use the novel methodological combination of single-unit intraneural microstimulation (INMS) with magnetoencephalography (MEG) to localise cortical representations of individual touch afferents in humans, by measuring the extracranial magnetic fields from neural currents. We found that by assessing the modulation of the beta (13-30 Hz) rhythm during single-unit INMS, significant changes in oscillatory amplitude occur in the contralateral primary somatosensory cortex within and across a group of fast adapting type I mechanoreceptive afferents, which corresponded well to the induced response from matched vibrotactile stimulation. Combining the spatiotemporal specificity of MEG with the selective single-unit stimulation of INMS enables the interrogation of the central representations of different aspects of tactile afferent signalling within the human cortices. The fundamental finding that single-unit INMS ERD responses are robust and consistent with natural somatosensory stimuli will permit us to more dynamically probe the central nervous system responses in humans, to address questions about the processing of touch from the different classes of mechanoreceptive afferents and the effects of varying the stimulus frequency and patterning.

Characterizing hippocampal dynamics with MEG: A systematic review and evidence-based guidelines

The hippocampus, a hub of activity for a variety of important cognitive processes, is a target of increasing interest for researchers and clinicians. Magnetoencephalography (MEG) is an attractive technique for imaging spectro-temporal aspects of function, for example, neural oscillations and network timing, especially in shallow cortical structures. However, the decrease in MEG signal-to-noise ratio as a function of source depth implies that the utility of MEG for investigations of deeper brain structures, including the hippocampus, is less clear. To determine whether MEG can be used to detect and localize activity from the hippocampus, we executed a systematic review of the existing literature and found successful detection of oscillatory neural activity originating in the hippocampus with MEG. Prerequisites are the use of established experimental paradigms, adequate coregistration, forward modeling, analysis methods, optimization of signal-to-noise ratios, and protocol trial designs that maximize contrast for hippocampal activity while minimizing those from other brain regions. While localizing activity to specific sub-structures within the hippocampus has not been achieved, we provide recommendations for improving the reliability of such endeavors.

Improved EEG source localization with Bayesian uncertainty modelling of unknown skull conductivity

Electroencephalography (EEG) source imaging is an ill-posed inverse problem that requires accurate conductivity modelling of the head tissues, especially the skull. Unfortunately, the conductivity values are difficult to determine in vivo. In this paper, we show that the exact knowledge of the skull conductivity is not always necessary when the Bayesian approximation error (BAE) approach is exploited. In BAE, we first postulate a probability distribution for the skull conductivity that describes our (lack of) knowledge on its value, and model the effects of this uncertainty on EEG recordings with the help of an additive error term in the observation model. Before the Bayesian inference, the likelihood is marginalized over this error term. Thus, in the inversion we estimate only our primary unknown, the source distribution. We quantified the improvements in the source localization when the proposed Bayesian modelling was used in the presence of different skull conductivity errors and levels of measurement noise. Based on the results, BAE was able to improve the source localization accuracy, particularly when the unknown (true) skull conductivity was much lower than the expected standard conductivity value. The source locations that gained the highest improvements were shallow and originally exhibited the largest localization errors. In our case study, the benefits of BAE became negligible when the signal-to-noise ratio dropped to 20 dB.

Requirements for Coregistration Accuracy in On-Scalp MEG

Recent advances in magnetic sensing has made on-scalp magnetoencephalography (MEG) possible. In particular, optically-pumped magnetometers (OPMs) have reached sensitivity levels that enable their use in MEG. In contrast to the SQUID sensors used in current MEG systems, OPMs do not require cryogenic cooling and can thus be placed within millimetres from the head, enabling the construction of sensor arrays that conform to the shape of an individual's head. To properly estimate the location of neural sources within the brain, one must accurately know the position and orientation of sensors in relation to the head. With the adaptable on-scalp MEG sensor arrays, this coregistration becomes more challenging than in current SQUID-based MEG systems that use rigid sensor arrays. Here, we used simulations to quantify how accurately one needs to know the position and orientation of sensors in an on-scalp MEG system. The effects that different types of localisation errors have on forward modelling and source estimates obtained by minimum-norm estimation, dipole fitting, and beamforming are detailed. We found that sensor position errors generally have a larger effect than orientation errors and that these errors affect the localisation accuracy of superficial sources the most. To obtain similar or higher accuracy than with current SQUID-based MEG systems, RMS sensor position and orientation errors should be [Formula: see text] and [Formula: see text], respectively.

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

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

Quantifying the performance of MEG source reconstruction using resting state data

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

Lamina-specific cortical dynamics in human visual and sensorimotor cortices

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

Proposal for an accurate TMS-MRI co-registration process via 3D laser scanning

An important technical issue in transcranial magnetic stimulation (TMS) usage is how accurately the specific brain areas activated by TMS are assessed. However, in practice, electric field induced in TMS is dispersed and therefore actual estimation is still difficult. As a preliminary step, the projection line which is perpendicular to the TMS stimulation coil beneath the center of the coil must be accurately estimated into the brain. Therefore, we have developed a new TMS-MRI co-registration procedure that employs a 3D laser-scanner system that is very useful for general hand-manipulated TMS, and which easily estimates the TMS projection point onto the brain. The proposed system accurately captures the positional relationship between the TMS coil and anatomical images. The results of 3D image processing revealed that the registration error at each stage was kept within the submillimeter level. In addition, a motor evoked potential experiment examining the right finger motor area revealed that understandable responses were obtained when stimulation was targeted to the three different motor areas according to Penfield's map. 3D laser scanning is a technique of substantial recent interest for anatomical co-registration. The proposed method demonstrated submillimeter level accuracy of TMS-MRI co-registration.

IFCN-endorsed practical guidelines for clinical magnetoencephalography (MEG)

Magnetoencephalography (MEG) records weak magnetic fields outside the human head and thereby provides millisecond-accurate information about neuronal currents supporting human brain function. MEG and electroencephalography (EEG) are closely related complementary methods and should be interpreted together whenever possible. This manuscript covers the basic physical and physiological principles of MEG and discusses the main aspects of state-of-the-art MEG data analysis. We provide guidelines for best practices of patient preparation, stimulus presentation, MEG data collection and analysis, as well as for MEG interpretation in routine clinical examinations. In 2017, about 200 whole-scalp MEG devices were in operation worldwide, many of them located in clinical environments. Yet, the established clinical indications for MEG examinations remain few, mainly restricted to the diagnostics of epilepsy and to preoperative functional evaluation of neurosurgical patients. We are confident that the extensive ongoing basic MEG research indicates potential for the evaluation of neurological and psychiatric syndromes, developmental disorders, and the integrity of cortical brain networks after stroke. Basic and clinical research is, thus, paving way for new clinical applications to be identified by an increasing number of practitioners of MEG.

Mapping the topological organisation of beta oscillations in motor cortex using MEG

The spatial topology of the human motor cortex has been well studied, particularly using functional Magnetic Resonance Imaging (fMRI) which allows spatial separation of haemodynamic responses arising from stimulation of different body parts, individual digits and even spatially separate areas of the same digit. However, the spatial organisation of electrophysiological responses, particularly neural oscillations (rhythmic changes in electrical potential across cellular assemblies) has been less well studied. Mapping the spatial signature of neural oscillations is possible using magnetoencephalography (MEG), however spatial differentiation of responses induced by movement of separate digits is a challenge, because the brain regions involved are separated by only a few millimetres. In this paper we first show, in simulation, how to optimise experimental design and beamformer spatial filtering techniques to increase the spatial specificity of MEG derived functional images. Combining this result with experimental data, we then capture the organisation of the post-movement beta band (13–30 Hz) oscillatory response to movement of digits 2 and 5 of the dominant hand, in individual subjects. By comparing these MEG results to ultra-high field (7T) fMRI, we also show significant spatial agreement between beta modulation and the blood oxygenation level dependent (BOLD) response. Our results show that, when using an optimised inverse solution and controlling subject movement (using custom fitted foam padding) the spatial resolution of MEG can be of order 3–5 mm. The method described offers exciting potential to understand better the cortical organisation of oscillations, and to probe such organisation in patient populations where those oscillations are known to be abnormal.

•

Aim is to map the topological organisation of neural oscillations in motor cortex.

•

MEG spatial resolution optimised by temporal separation of sources.

•

Subject motion controlled using foam headcasts.

•

Cortical representation of Digit 2 and Digit 5 separated spatially.

•

Post movement beta rebound maps motortopically in agreement with BOLD responses.

Localizing on-scalp MEG sensors using an array of magnetic dipole coils

Accurate estimation of the neural activity underlying magnetoencephalography (MEG) signals requires co-registration i.e., determination of the position and orientation of the sensors with respect to the head. In modern MEG systems, an array of hundreds of low-Tc SQUID sensors is used to localize a set of small, magnetic dipole-like (head-position indicator, HPI) coils that are attached to the subject's head. With accurate prior knowledge of the positions and orientations of the sensors with respect to one another, the HPI coils can be localized with high precision, and thereby the positions of the sensors in relation to the head. With advances in magnetic field sensing technologies, e.g., high-Tc SQUIDs and optically pumped magnetometers (OPM), that require less extreme operating temperatures than low-Tc SQUID sensors, on-scalp MEG is on the horizon. To utilize the full potential of on-scalp MEG, flexible sensor arrays are preferable. Conventional co-registration is impractical for such systems as the relative positions and orientations of the sensors to each other are subject-specific and hence not known a priori. Herein, we present a method for co-registration of on-scalp MEG sensors. We propose to invert the conventional co-registration approach and localize the sensors relative to an array of HPI coils on the subject's head. We show that given accurate prior knowledge of the positions of the HPI coils with respect to one another, the sensors can be localized with high precision. We simulated our method with realistic parameters and layouts for sensor and coil arrays. Results indicate co-registration is possible with sub-millimeter accuracy, but the performance strongly depends upon a number of factors. Accurate calibration of the coils and precise determination of the positions and orientations of the coils with respect to one another are crucial. Finally, we propose methods to tackle practical challenges to further improve the method.

Non-invasive Investigation of Human Hippocampal Rhythms Using Magnetoencephalography: A Review

Hippocampal rhythms are believed to support crucial cognitive processes including memory, navigation, and language. Due to the location of the hippocampus deep in the brain, studying hippocampal rhythms using non-invasive magnetoencephalography (MEG) recordings has generally been assumed to be methodologically challenging. However, with the advent of whole-head MEG systems in the 1990s and development of advanced source localization techniques, simulation and empirical studies have provided evidence that human hippocampal signals can be sensed by MEG and reliably reconstructed by source localization algorithms. This paper systematically reviews simulation studies and empirical evidence of the current capacities and limitations of MEG “deep source imaging” of the human hippocampus. Overall, these studies confirm that MEG provides a unique avenue to investigate human hippocampal rhythms in cognition, and can bridge the gap between animal studies and human hippocampal research, as well as elucidate the functional role and the behavioral correlates of human hippocampal oscillations.

On the Physiological Modulation and Potential Mechanisms Underlying Parieto-Occipital Alpha Oscillations

The parieto-occipital alpha (8–13 Hz) rhythm is by far the strongest spectral fingerprint in the human brain. Almost 90 years later, its physiological origin is still far from clear. In this Research Topic I review human pharmacological studies using electroencephalography (EEG) and magnetoencephalography (MEG) that investigated the physiological mechanisms behind posterior alpha. Based on results from classical and recent experimental studies, I find a wide spectrum of drugs that modulate parieto-occipital alpha power. Alpha frequency is rarely affected, but this might be due to the range of drug dosages employed. Animal and human pharmacological findings suggest that both GABA enhancers and NMDA blockers systematically decrease posterior alpha power. Surprisingly, most of the theoretical frameworks do not seem to embrace these empirical findings and the debate on the functional role of alpha oscillations has been polarized between the inhibition vs. active poles hypotheses. Here, I speculate that the functional role of alpha might depend on physiological excitation as much as on physiological inhibition. This is supported by animal and human pharmacological work showing that GABAergic, glutamatergic, cholinergic, and serotonergic receptors in the thalamus and the cortex play a key role in the regulation of alpha power and frequency. This myriad of physiological modulations fit with the view that the alpha rhythm is a complex rhythm with multiple sources supported by both thalamo-cortical and cortico-cortical loops. Finally, I briefly discuss how future research combining experimental measurements derived from theoretical predictions based of biophysically realistic computational models will be crucial to the reconciliation of these disparate findings.

Quality assessment of MEG-to-MRI coregistrations

For high precision in source reconstruction of magnetoencephalography (MEG) or electroencephalography data, high accuracy of the coregistration of sources and sensors is mandatory. Usually, the source space is derived from magnetic resonance imaging (MRI). In most cases, however, no quality assessment is reported for sensor-to-MRI coregistrations. If any, typically root mean squares (RMS) of point residuals are provided. It has been shown, however, that RMS of residuals do not correlate with coregistration errors. We suggest using target registration error (TRE) as criterion for the quality of sensor-to-MRI coregistrations. TRE measures the effect of uncertainty in coregistrations at all points of interest. In total, 5544 data sets with sensor-to-head and 128 head-to-MRI coregistrations, from a single MEG laboratory, were analyzed. An adaptive Metropolis algorithm was used to estimate the optimal coregistration and to sample the coregistration parameters (rotation and translation). We found an average TRE between 1.3 and 2.3 mm at the head surface. Further, we observed a mean absolute difference in coregistration parameters between the Metropolis and iterative closest point algorithm of [Formula: see text] and [Formula: see text] m. A paired sample t-test indicated a significant improvement in goal function minimization by using the Metropolis algorithm. The sampled parameters allowed computation of TRE on the entire grid of the MRI volume. Hence, we recommend the Metropolis algorithm for head-to-MRI coregistrations.

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.

Evaluation of realistic layouts for next generation on-scalp MEG: spatial information density maps

While commercial magnetoencephalography (MEG) systems are the functional neuroimaging state-of-the-art in terms of spatio-temporal resolution, MEG sensors have not changed significantly since the 1990s. Interest in newer sensors that operate at less extreme temperatures, e.g., high critical temperature (high-T c) SQUIDs, optically-pumped magnetometers, etc., is growing because they enable significant reductions in head-to-sensor standoff (on-scalp MEG). Various metrics quantify the advantages of on-scalp MEG, but a single straightforward one is lacking. Previous works have furthermore been limited to arbitrary and/or unrealistic sensor layouts. We introduce spatial information density (SID) maps for quantitative and qualitative evaluations of sensor arrays. SID-maps present the spatial distribution of information a sensor array extracts from a source space while accounting for relevant source and sensor parameters. We use it in a systematic comparison of three practical on-scalp MEG sensor array layouts (based on high-T c SQUIDs) and the standard Elekta Neuromag TRIUX magnetometer array. Results strengthen the case for on-scalp and specifically high-T c SQUID-based MEG while providing a path for the practical design of future MEG systems. SID-maps are furthermore general to arbitrary magnetic sensor technologies and source spaces and can thus be used for design and evaluation of sensor arrays for magnetocardiography, magnetic particle imaging, etc.

An intra-neural microstimulation system for ultra-high field magnetic resonance imaging and magnetoencephalography

•

We propose an intra-neural microstimulation system for 7 T fMRI and MEG.

•

This custom-built system removes issues with existing equipment.

•

It provides efficient work-flow and improved participant comfort and safety.

•

Stimulating single mechanoreceptors evokes activity in 7 T fMRI and MEG.

•

Responses to unitary stimulation are shown for the first time in MEG.

Background

Intra-neural microstimulation (INMS) is a technique that allows the precise delivery of low-current electrical pulses into human peripheral nerves. Single unit INMS can be used to stimulate individual afferent nerve fibres during microneurography. Combining this with neuroimaging allows the unique monitoring of central nervous system activation in response to unitary, controlled tactile input, with functional magnetic resonance imaging (fMRI) providing exquisite spatial localisation of brain activity and magnetoencephalography (MEG) high temporal resolution.

New method

INMS systems suitable for use within electrophysiology laboratories have been available for many years. We describe an INMS system specifically designed to provide compatibility with both ultra-high field (7 T) fMRI and MEG. Numerous technical and safety issues are addressed. The system is fully analogue, allowing for arbitrary frequency and amplitude INMS stimulation.

Results

Unitary recordings obtained within both the MRI and MEG screened-room environments are comparable with those obtained in ‘clean’ electrophysiology recording environments. Single unit INMS (current <7 μA, 200 μs pulses) of individual mechanoreceptive afferents produces appropriate and robust responses during fMRI and MEG.

Comparison with existing method(s)

This custom-built MRI- and MEG-compatible stimulator overcomes issues with existing INMS approaches; it allows well-controlled switching between recording and stimulus mode, prevents electrical shocks because of long cable lengths, permits unlimited patterns of stimulation, and provides a system with improved work-flow and participant comfort.

Conclusions

We demonstrate that the requirements for an INMS-integrated system, which can be used with both fMRI and MEG imaging systems, have been fully met.

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.

A Review of Issues Related to Data Acquisition and Analysis in EEG/MEG Studies

Electroencephalography (EEG) and magnetoencephalography (MEG) are non-invasive electrophysiological methods, which record electric potentials and magnetic fields due to electric currents in synchronously-active neurons. With MEG being more sensitive to neural activity from tangential currents and EEG being able to detect both radial and tangential sources, the two methods are complementary. Over the years, neurophysiological studies have changed considerably: high-density recordings are becoming de rigueur; there is interest in both spontaneous and evoked activity; and sophisticated artifact detection and removal methods are available. Improved head models for source estimation have also increased the precision of the current estimates, particularly for EEG and combined EEG/MEG. Because of their complementarity, more investigators are beginning to perform simultaneous EEG/MEG studies to gain more complete information about neural activity. Given the increase in methodological complexity in EEG/MEG, it is important to gather data that are of high quality and that are as artifact free as possible. Here, we discuss some issues in data acquisition and analysis of EEG and MEG data. Practical considerations for different types of EEG and MEG studies are also discussed.

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.

Optimising experimental design for MEG resting state functional connectivity measurement

The study of functional connectivity using magnetoencephalography (MEG) is an expanding area of neuroimaging, and adds an extra dimension to the more common assessments made using fMRI. The importance of such metrics is growing, with recent demonstrations of their utility in clinical research, however previous reports suggest that whilst group level resting state connectivity is robust, single session recordings lack repeatability. Such robustness is critical if MEG measures in individual subjects are to prove clinically valuable. In the present paper, we test how practical aspects of experimental design affect the intra-subject repeatability of MEG findings; specifically we assess the effect of co-registration method and data recording duration. We show that the use of a foam head-cast, which is known to improve co-registration accuracy, increased significantly the between session repeatability of both beamformer reconstruction and connectivity estimation. We also show that recording duration is a critical parameter, with large improvements in repeatability apparent when using ten minute, compared to five minute recordings. Further analyses suggest that the origin of this latter effect is not underpinned by technical aspects of source reconstruction, but rather by a genuine effect of brain state; short recordings are simply inefficient at capturing the canonical MEG network in a single subject. Our results provide important insights on experimental design and will prove valuable for future MEG connectivity studies.