Transforming and comparing data between standard SQUID and OPM-MEG systems

Source localization comparison and combination of OPM-MEG and fMRI to detect sensorimotor cortex responses

BACKGROUND AND OBJECTIVES

The exploration of various neuroimaging techniques have become focal points within the field of neuroscience research. Magnetoencephalography based on optically pumped magnetometers (OPM-MEG) has shown significant potential to be the next generation of functional neuroimaging with the advantages of high signal intensity and flexible sensor arrangement. In this study, we constructed a 31-channel OPM-MEG system and performed a preliminary comparison of the temporal and spatial relationship between magnetic responses measured by OPM-MEG and blood-oxygen-level-dependent signals detected by functional magnetic resonance imaging (fMRI) during a grasping task.

METHODS

For OPM-MEG, the β-band (15-30 Hz) oscillatory activities can be reliably detected across multiple subjects and multiple session runs. To effectively localize the inhibitory oscillatory activities, a source power-spectrum ratio-based imaging method was proposed. This approach was compared with conventional source imaging methods, such as minimum norm-type and beamformer methods, and was applied in OPM-MEG source analysis. Subsequently, the spatial and temporal responses at the source-level between OPM-MEG and fMRI were analyzed.

RESULTS

The effectiveness of the proposed method was confirmed through simulations compared to benchmark methods. Our demonstration revealed an average spatial separation of 10.57 ± 4.41 mm between the localization results of OPM-MEG and fMRI across four subjects. Furthermore, the fMRI-constrained OPM-MEG localization results indicated a more focused imaging extent.

CONCLUSIONS

Taken together, the performance exhibited by OPM-MEG positions it as a potential instrument for functional surgery assessment.

Applications of OPM-MEG for translational neuroscience: a perspective

Magnetoencephalography (MEG) allows the non-invasive measurement of brain activity at millisecond precision combined with localization of the underlying generators. So far, MEG-systems consisted of superconducting quantum interference devices (SQUIDS), which suffer from several limitations. Recent technological advances, however, have enabled the development of novel MEG-systems based on optically pumped magnetometers (OPMs), offering several advantages over conventional SQUID-MEG systems. Considering potential improvements in the measurement of neuronal signals as well as reduced operating costs, the application of OPM-MEG systems for clinical neuroscience and diagnostic settings is highly promising. Here we provide an overview of the current state-of-the art of OPM-MEG and its unique potential for translational neuroscience. First, we discuss the technological features of OPMs and benchmark OPM-MEG against SQUID-MEG and electroencephalography (EEG), followed by a summary of pioneering studies of OPMs in healthy populations. Key applications of OPM-MEG for the investigation of psychiatric and neurological conditions are then reviewed. Specifically, we suggest novel applications of OPM-MEG for the identification of biomarkers and circuit deficits in schizophrenia, dementias, movement disorders, epilepsy, and neurodevelopmental syndromes (autism spectrum disorder and attention deficit hyperactivity disorder). Finally, we give an outlook of OPM-MEG for translational neuroscience with a focus on remaining methodological and technical challenges.

Expanding the clinical application of OPM-MEG using an effective automatic suppression method for the dental brace metal artifact

Optically pumped magnetometer magnetoencephalography (OPM-MEG) holds significant promise for clinical functional brain imaging due to its superior spatiotemporal resolution. However, effectively suppressing metallic artifacts, particularly from devices such as orthodontic braces and vagal nerve stimulators remains a major challenge, hindering the wider clinical application of wearable OPM-MEG devices. A comprehensive analysis of metal artifact characteristics from time, frequency, and time-frequency perspectives was conducted for the first time using an OPM-MEG device in clinical medicine. This study focused on patients with metal orthodontics, examining the modulation of metal artifacts by breath and head movement, the incomplete regular sub-Gaussian distribution, and the high absolute power ratio in the 0.5-8 Hz band. The existing metal artifact suppression algorithms applied to SQUID-MEG, such as fast independent component analysis (FastICA), information maximization (Infomax), and algorithms for multiple unknown signal extraction (AMUSE), exhibit limited efficacy. Consequently, this study introduced the second-order blind identification (SOBI) algorithm, which utilized multiple time delays for the component separation of OPM-MEG measurement signals. We modified the time delays of the SOBI method to improve its efficacy in separating artifact components, particularly those in the ultralow frequency range. This approach employs the frequency-domain absolute power ratio, root mean square (RMS) value, and mutual information methods to automate the artifact component screening process. The effectiveness of this method was validated through simulation experiments involving four subjects in both resting and evoked experiments. In addition, the proposed method was also validated by the actual OPM-MEG evoked experiments of three subjects. Comparative analyses were conducted against the FastICA, Infomax, and AMUSE algorithms. Evaluation metrics included normalized mean square error, normalized delta band power error, RMS error, and signal-to-noise ratio, demonstrating that the proposed method provides optimal suppression of metal artifacts. This advancement holds promise for enhancing data quality and expanding the clinical applications of OPM-MEG.

Decoding N400m Evoked Component: A Tutorial on Multivariate Pattern Analysis for OP-MEG Data

Multivariate pattern analysis (MVPA) has played an extensive role in interpreting brain activity, which has been applied in studies with modalities such as functional Magnetic Resonance Imaging (fMRI), Magnetoencephalography (MEG) and Electroencephalography (EEG). The advent of wearable MEG systems based on optically pumped magnetometers (OPMs), i.e., OP-MEG, has broadened the application of bio-magnetism in the realm of neuroscience. Nonetheless, it also raises challenges in temporal decoding analysis due to the unique attributes of OP-MEG itself. The efficacy of decoding performance utilizing multimodal fusion, such as MEG-EEG, also remains to be elucidated. In this regard, we investigated the impact of several factors, such as processing methods, models and modalities, on the decoding outcomes of OP-MEG. Our findings indicate that the number of averaged trials, dimensionality reduction (DR) methods, and the number of cross-validation folds significantly affect the decoding performance of OP-MEG data. Additionally, decoding results vary across modalities and fusion strategy. In contrast, decoder type, resampling frequency, and sliding window length exert marginal effects. Furthermore, we introduced mutual information (MI) to investigate how information loss due to OP-MEG data processing affect decoding accuracy. Our study offers insights for linear decoding research using OP-MEG and expand its application in the fields of cognitive neuroscience.

A simulation study: comparing independent component analysis and signal-space projection - source-informed reconstruction for rejecting muscle artifacts evoked by transcranial magnetic stimulation

Introduction

The combination of transcranial magnetic stimulation (TMS) and electroencephalography (EEG) allows researchers to explore cortico-cortical connections. To study effective connections, the first few tens of milliseconds of the TMS-evoked potentials are the most critical. Yet, TMS-evoked artifacts complicate the interpretation of early-latency data. Data-processing strategies like independent component analysis (ICA) and the combined signal-space projection-source-informed reconstruction approach (SSP-SIR) are designed to mitigate artifacts, but their objective assessment is challenging because the true neuronal EEG responses under large-amplitude artifacts are generally unknown. Through simulations, we quantified how the spatiotemporal properties of the artifacts affect the cleaning performances of ICA and SSP-SIR.

Methods

We simulated TMS-induced muscle artifacts and superposed them on pre-processed TMS-EEG data, serving as the ground truth. The simulated muscle artifacts were varied both in terms of their topography and temporal profiles. The signals were then cleaned using ICA and SSP-SIR, and subsequent comparisons were made with the ground truth data.

Results

ICA performed better when the artifact time courses were highly variable across the trials, whereas the effectiveness of SSP-SIR depended on the congruence between the artifact and neuronal topographies, with the performance of SSP-SIR being better when difference between topographies was larger. Overall, SSP-SIR performed better than ICA across the tested conditions. Based on these simulations, SSP-SIR appears to be more effective in suppressing TMS-evoked muscle artifacts. These artifacts are shown to be highly time-locked to the TMS pulse and manifest in topographies that differ substantially from the patterns of neuronal potentials.

Discussion

Selecting between ICA and SSP-SIR should be guided by the characteristics of the artifacts. SSP-SIR might be better equipped for suppressing time-locked artifacts, provided that their topographies are sufficiently different from the neuronal potential patterns of interest, and that the SSP-SIR algorithm can successfully find those artifact topographies from the high-pass-filtered data. ICA remains a powerful tool for rejecting artifacts that are not strongly time locked to the TMS pulse.

Single-trial classification of evoked responses to auditory tones using OPM- and SQUID-MEG

Objective.Optically pumped magnetometers (OPMs) are emerging as a near-room-temperature alternative to superconducting quantum interference devices (SQUIDs) for magnetoencephalography (MEG). In contrast to SQUIDs, OPMs can be placed in a close proximity to subject's scalp potentially increasing the signal-to-noise ratio and spatial resolution of MEG. However, experimental demonstrations of these suggested benefits are still scarce. Here, to compare a 24-channel OPM-MEG system to a commercial whole-head SQUID system in a data-driven way, we quantified their performance in classifying single-trial evoked responses.Approach.We measured evoked responses to three auditory tones in six participants using both OPM- and SQUID-MEG systems. We performed pairwise temporal classification of the single-trial responses with linear discriminant analysis as well as multiclass classification with both EEGNet convolutional neural network and xDAWN decoding.Main results.OPMs provided higher classification accuracies than SQUIDs having a similar coverage of the left hemisphere of the participant. However, the SQUID sensors covering the whole helmet had classification scores larger than those of OPMs for two of the tone pairs, demonstrating the benefits of a whole-head measurement.Significance.The results demonstrate that the current OPM-MEG system provides high-quality data about the brain with room for improvement for high bandwidth non-invasive brain-computer interfacing.

Enabling ambulatory movement in wearable magnetoencephalography with matrix coil active magnetic shielding

The ability to collect high-quality neuroimaging data during ambulatory participant movement would enable a wealth of neuroscientific paradigms. Wearable magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) has the potential to allow participant movement during a scan. However, the strict zero magnetic field requirement of OPMs means that systems must be operated inside a magnetically shielded room (MSR) and also require active shielding using electromagnetic coils to cancel residual fields and field changes (due to external sources and sensor movements) that would otherwise prevent accurate neuronal source reconstructions. Existing active shielding systems only compensate fields over small, fixed regions and do not allow ambulatory movement. Here we describe the matrix coil, a new type of active shielding system for OPM-MEG which is formed from 48 square unit coils arranged on two planes which can compensate magnetic fields in regions that can be flexibly placed between the planes. Through the integration of optical tracking with OPM data acquisition, field changes induced by participant movement are cancelled with low latency (25 ms). High-quality MEG source data were collected despite the presence of large (65 cm translations and 270° rotations) ambulatory participant movements.

Improved Biomagnetic Signal-To-Noise Ratio and Source Localization Using Optically Pumped Magnetometers with Synthetic Gradiometers

Optically pumped magnetometers (OPMs) can capture brain activity but are susceptible to magnetic noise. The objective of this study was to evaluate a novel methodology used to reduce magnetic noise in OPM measurements. A portable magnetoencephalography (MEG) prototype was developed with OPMs. The OPMs were divided into primary sensors and reference sensors. For each primary sensor, a synthetic gradiometer (SG) was constructed by computing a secondary sensor that simulated noise with signals from the reference sensors. MEG data from a phantom with known source signals and six human participants were used to assess the efficacy of the SGs. Magnetic noise in the OPM data appeared predominantly in a low frequency range (<4 Hz) and varied among OPMs. The SGs significantly reduced magnetic noise (p < 0.01), enhanced the signal-to-noise ratio (SNR) (p < 0.001) and improved the accuracy of source localization (p < 0.02). The SGs precisely revealed movement-evoked magnetic fields in MEG data recorded from human participants. SGs provided an effective method to enhance SNR and improve the accuracy of source localization by suppressing noise. Software-simulated SGs may provide new opportunities regarding the use of OPM measurements in various clinical and research applications, especially those in which movement is relevant.

Monitoring magnetic nanoparticle clustering and immobilization with thermal noise magnetometry using optically pumped magnetometers

Thermal noise magnetometry (TNM) is a recently developed magnetic characterization technique where thermally induced fluctuations in magnetization are measured to gain insight into nanomagnetic structures like magnetic nanoparticles (MNPs). Due to the stochastic nature of the method, its signal amplitude scales with the square of the volume of the individual fluctuators, which makes the method therefore extra attractive to study MNP clustering and aggregation processes. Until now, TNM signals have exclusively been detected by using a superconducting quantum interference device (SQUID) sensor. In contrast, we present here a tabletop setup using optically pumped magnetometers (OPMs) in a compact magnetic shield, as a flexible alternative. The agreement between results obtained with both measurement systems is shown for different commercially available MNP samples. We argue that the OPM setup with low complexity complements the SQUID setup with high sensitivity and bandwidth. Furthermore, the OPM tabletop setup is well suited to monitor aggregation processes because of its excellent sensitivity in lower frequencies. As a proof of concept, we show the changes in the noise spectrum for three different MNP immobilization and clustering processes. From our results, we conclude that the tabletop setup offers a flexible and widely adoptable measurement unit to monitor the immobilization, aggregation, and clustering of MNPs for different applications, including interactions of the particles with biological systems and the long-term stability of magnetic samples.

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.

Towards a Practical Implementation of a Single-Beam All-Optical Non-Zero-Field Magnetic Sensor for Magnetoencephalographic Complexes

We present a single-beam all-optical two-channel magnetic sensor scheme developed for biological applications such as non-zero-field magnetoencephalography and magnetocardiography. The pumping, excitation and detection of magnetic resonance in two cells are performed using a single laser beam with time-modulated linear polarization: the linear polarization of the beam switches to orthogonal every half-cycle of the Larmor frequency. Light with such characteristics can be transmitted over a single-mode polarization-maintaining fiber without any loss in the quality of the polarization characteristics. We also present an algorithm for calculating optical elements in a sensor scheme, the results of measuring the parametric dependences of magnetic resonance in cells, and the results of direct testing of a sensor in a magnetic shield. We demonstrate sensitivity at the level of 20 fT/√Hz in one sensor channel in the frequency range of 80-200 Hz.

Optimising the sensing volume of OPM sensors for MEG source reconstruction

Magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) has been hailed as the future of electrophysiological recordings from the human brain. In this work, we investigate how the dimensions of the sensing volume (the vapour cell) affect the performance of both a single OPM-MEG sensor and a multi-sensor OPM-MEG system. We consider a realistic noise model that accounts for background brain activity and residual noise. By using source reconstruction metrics such as localization accuracy and time-course reconstruction accuracy, we demonstrate that the best overall sensitivity and reconstruction accuracy are achieved with cells that are significantly longer and wider that those of the majority of current commercial OPM sensors. Our work provides useful tools to optimise the cell dimensions of OPM sensors in a wide range of environments.

Simulation Study of Different OPM-MEG Measurement Components

Magnetoencephalography (MEG) is a neuroimaging technique that measures the magnetic fields of the brain outside of the head. In the past, the most suitable magnetometer for MEG was the superconducting quantum interference device (SQUID), but in recent years, a new type has also been used, the optically pumped magnetometer (OPM). OPMs can be configured to measure multiple directions of magnetic field simultaneously. This work explored whether combining multiple directions of the magnetic field lowers the source localization error of brain sources under various conditions of noise. We simulated dipolar-like sources for multiple configurations of both SQUID- and OPM-MEG systems. To test the performance of a given layout, we calculated the average signal-to-noise ratio and the root mean square of the simulated magnetic field; furthermore, we evaluated the performance of the dipole fit. The results showed that the field direction normal to the scalp yields a higher signal-to-noise ratio and that ambient noise has a much lower impact on its localization error; therefore, this is the optimal choice for source localization when only one direction of magnetic field can be measured. For a low number of OPMs, combining multiple field directions greatly improves the source localization results. Lastly, we showed that MEG sensors that can be placed closer to the brain are more suitable for localizing deeper sources.