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.

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

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

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

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.

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

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

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

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

Using OPMs to measure neural activity in standing, mobile participants.

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

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

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

Mouth magnetoencephalography: A unique perspective on the human hippocampus.

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

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.

Flexible proton density (PD) mapping using multi-contrast variable flip angle (VFA) data.

Quantitative proton density (PD) maps measure the amount of free water, which is important for non-invasive tissue characterization in pathology and across lifespan. PD mapping requires the estimation and subsequent removal of factors influencing the signal intensity other than PD. These factors include the T1, T2* relaxation effects, transmit field inhomogeneities, receiver coil sensitivity profile (RP) and the spatially invariant factor that is required to scale the data. While the transmit field can be reliably measured, the RP estimation is usually based on image post-processing techniques due to limitations of its measurement at magnetic fields higher than 1.5 T. The post-processing methods are based on unified bias-field/tissue segmentation, fitting the sensitivity profile from images obtained with different coils, or on the linear relationship between T1 and PD. The scaling factor is derived from the signal within a specific tissue compartment or reference object. However, these approaches for calculating the RP and scaling factor have limitations particularly in severe pathology or over a wide age range, restricting their application. We propose a new approach for PD mapping based on a multi-contrast variable flip angle acquisition protocol and a data-driven estimation method for the RP correction and map scaling. By combining all the multi-contrast data acquired at different echo times, we are able to fully correct the MRI signal for T2* relaxation effects and to decrease the variance and the entropy of PD values within tissue class of the final map. The RP is determined from the corrected data applying a non-parametric bias estimation, and the scaling factor is based on the median intensity of an external calibration object. Finally, we compare the signal intensity and homogeneity of the multi-contrast PD map with the well-established effective PD (PD*) mapping, for which the RP is based on concurrent bias field estimation and tissue classification, and the scaling factor is estimated from the mean white matter signal. The multi-contrast PD values homogeneity and accuracy within the cerebrospinal fluid (CSF) and deep brain structures are increased beyond that obtained using PD* maps. We demonstrate that the multi-contrast RP approach is insensitive to anatomical or a priori tissue information by applying it in a patient with extensive brain abnormalities and for whole body PD mapping in post-mortem foetal imaging.

Imaging the human hippocampus with optically-pumped magnetoencephalography.

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

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

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

Wearable neuroimaging: Combining and contrasting magnetoencephalography and electroencephalography.

One of the most severe limitations of functional neuroimaging techniques, such as magnetoencephalography (MEG), is that participants must maintain a fixed head position during data acquisition. This imposes restrictions on the characteristics of the experimental cohorts that can be scanned and the experimental questions that can be addressed. For these reasons, the use of 'wearable' neuroimaging, in which participants can move freely during scanning, is attractive. The most successful example of wearable neuroimaging is electroencephalography (EEG), which employs lightweight and flexible instrumentation that makes it useable in almost any experimental setting. However, EEG has major technical limitations compared to MEG, and therefore the development of wearable MEG, or hybrid MEG/EEG systems, is a compelling prospect. In this paper, we combine and compare EEG and MEG measurements, the latter made using a new generation of optically-pumped magnetometers (OPMs). We show that these new second generation commercial OPMs, can be mounted on the scalp in an 'EEG-like' cap, enabling the acquisition of high fidelity electrophysiological measurements. We show that these sensors can be used in conjunction with conventional EEG electrodes, offering the potential for the development of hybrid MEG/EEG systems. We compare concurrently measured signals, showing that, whilst both modalities offer high quality data in stationary subjects, OPM-MEG measurements are less sensitive to artefacts produced when subjects move. Finally, we show using simulations that OPM-MEG offers a fundamentally better spatial specificity than EEG. The demonstrated technology holds the potential to revolutionise the utility of functional brain imaging, exploiting the flexibility of wearable systems to facilitate hitherto impractical experimental paradigms.

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.

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

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

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

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

Combined electroencephalography-functional magnetic resonance imaging and electrical source imaging improves localization of pediatric focal epilepsy.

OBJECTIVE: Surgical treatment in epilepsy is effective if the epileptogenic zone (EZ) can be correctly localized and characterized. Here we use simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) data to derive EEG-fMRI and electrical source imaging (ESI) maps. Their yield and their individual and combined ability to (1) localize the EZ and (2) predict seizure outcome were then evaluated. METHODS: Fifty-three children with drug-resistant epilepsy underwent EEG-fMRI. Interictal discharges were mapped using both EEG-fMRI hemodynamic responses and ESI. A single localization was derived from each individual test (EEG-fMRI global maxima [GM]/ESI maximum) and from the combination of both maps (EEG-fMRI/ESI spatial intersection). To determine the localization accuracy and its predictive performance, the individual and combined test localizations were compared to the presumed EZ and to the postsurgical outcome. RESULTS: Fifty-two of 53 patients had significant maps: 47 of 53 for EEG-fMRI, 44 of 53 for ESI, and 34 of 53 for both. The EZ was well characterized in 29 patients; 26 had an EEG-fMRI GM localization that was correct in 11, 22 patients had ESI localization that was correct in 17, and 12 patients had combined EEG-fMRI and ESI that was correct in 11. Seizure outcome following resection was correctly predicted by EEG-fMRI GM in 8 of 20 patients, and by the ESI maximum in 13 of 16. The combined EEG-fMRI/ESI region entirely predicted outcome in 9 of 9 patients, including 3 with no lesion visible on MRI. INTERPRETATION: EEG-fMRI combined with ESI provides a simple unbiased localization that may predict surgery better than each individual test, including in MRI-negative patients. Ann Neurol 2017;82:278-287.

Thalamic volume reduction in drug-naive patients with new-onset genetic generalized epilepsy.

OBJECTIVE: Patients with genetic generalized epilepsy (GGE) have subtle morphologic abnormalities of the brain revealed with magnetic resonance imaging (MRI), particularly in the thalamus. However, it is unclear whether morphologic abnormalities of the brain in GGE are a consequence of repeated seizures over the duration of the disease, or are a consequence of treatment with antiepileptic drugs (AEDs), or are independent of these factors. Therefore, we measured brain morphometry in a cohort of AED-naive patients with GGE at disease onset. We hypothesize that drug-naive patients at disease onset have gray matter changes compared to age-matched healthy controls. METHODS: We performed quantitative measures of gray matter volume in the thalamus, putamen, caudate, pallidum, hippocampus, precuneus, prefrontal cortex, precentral cortex, and cingulate in 29 AED-naive patients with new-onset GGE and compared them to 32 age-matched healthy controls. We subsequently compared the shape of any brain structures found to differ in gray matter volume between the groups. RESULTS: The thalamus was the only structure to show reduced gray matter volume in AED-naive patients with new-onset GGE compared to healthy controls. Shape analysis revealed that the thalamus showed deflation, which was not uniformly distributed, but particularly affected a circumferential strip involving anterior, superior, posterior, and inferior regions with sparing of medial and lateral regions. SIGNIFICANCE: Structural abnormalities in the thalamus are present at the initial onset of GGE in AED-naive patients, suggesting that thalamic structural abnormality is an intrinsic feature of GGE and not a consequence of AEDs or disease duration.

Phase-amplitude coupling and the BOLD signal: A simultaneous intracranial EEG (icEEG) - fMRI study in humans performing a finger-tapping task.

Although it has been consistently found that local blood-oxygen-level-dependent (BOLD) changes are better modelled by a combination of the power of multiple EEG frequency bands rather than by the power of a unique band alone, the local electro-haemodynamic coupling function is not yet fully characterised. Electrophysiological studies have revealed that the strength of the coupling between the phase of low- and the amplitude of high- frequency EEG activities (phase-amplitude coupling - PAC) has an important role in brain function in general, and in preparation and execution of movement in particular. Using electrocorticographic (ECoG) and functional magnetic resonance imaging (fMRI) data recorded simultaneously in humans performing a finger-tapping task, we investigated the single-trial relationship between the amplitude of the BOLD signal and the strength of PAC and the power of α, ß, and γ bands, at a local level. In line with previous studies, we found a positive correlation for the γ band, and negative correlations for the PACßγ strength, and the α and ß bands. More importantly, we found that the PACßγ strength explained variance of the amplitude of the BOLD signal that was not explained by a combination of the α, ß, and γ band powers. Our main finding sheds further light on the distinct nature of PAC as a functionally relevant mechanism and suggests that the sensitivity of EEG-informed fMRI studies may increase by including the PAC strength in the BOLD signal model, in addition to the power of the low- and high- frequency EEG bands.