The impact of improved MEG-MRI co-registration on MEG connectivity analysis

High-Definition MEG Source Estimation using the Reciprocal Boundary Element Fast Multipole Method

Magnetoencephalographic (MEG) source estimation relies on the computation of the gain (lead-field) matrix, which embodies the linear relationship between the amplitudes of the sources and the recorded signals. However, with a realistic forward model, the calculation of the gain matrix in a "direct" fashion is a computationally expensive task because the number of dipolar sources in standard MEG pipelines is often limited to ~10,000. We propose a fast approach based on the reciprocal relationship between MEG and transcranial magnetic stimulation (TMS). This approach couples naturally with the charge-based boundary element fast multipole method (BEM-FMM), which allows us to efficiently generate gain matrices for high-resolution multi-layer non-nested meshes involving source spaces of up to a ~1 million dipoles. We evaluate our approach by performing MEG source reconstruction against simulated data (at varying noise levels) obtained from the direct computation of MEG readings from 2000 different dipole positions over the cortical surface of 5 healthy subjects. Additionally, we test our methods with real MEG data from evoked somatosensory fields by right-hand median nerve stimulation in these same 5 subjects. We compare our experimental source reconstruction results against the standard MNE-Python source reconstruction pipeline.

PyFaceWipe: a new defacing tool for almost any MRI contrast

RATIONALE AND OBJECTIVES

Defacing research MRI brain scans is often a mandatory step. With current defacing software, there are issues with Windows compatibility and researcher doubt regarding the adequacy of preservation of brain voxels in non-T1w scans. To address this, we developed PyFaceWipe, a multiplatform software for multiple MRI contrasts, which was evaluated based on its anonymisation ability and effect on downstream processing.

MATERIALS AND METHODS

Multiple MRI brain scan contrasts from the OASIS-3 dataset were defaced with PyFaceWipe and PyDeface and manually assessed for brain voxel preservation, remnant facial features and effect on automated face detection. Original and PyFaceWipe-defaced data from locally acquired T1w structural scans underwent volumetry with FastSurfer and brain atlas generation with ANTS.

RESULTS

214 MRI scans of several contrasts from OASIS-3 were successfully processed with both PyFaceWipe and PyDeface. PyFaceWipe maintained complete brain voxel preservation in all tested contrasts except ASL (45%) and DWI (90%), and PyDeface in all tested contrasts except ASL (95%), BOLD (25%), DWI (40%) and T2* (25%). Manual review of PyFaceWipe showed no failures of facial feature removal. Pinna removal was less successful (6% of T1 scans showed residual complete pinna). PyDeface achieved 5.1% failure rate. Automated detection found no faces in PyFaceWipe-defaced scans, 19 faces in PyDeface scans compared with 78 from the 224 original scans. Brain atlas generation showed no significant difference between atlases created from original and defaced data in both young adulthood and late elderly cohorts. Structural volumetry dice scores were ≥ 0.98 for all structures except for grey matter which had 0.93. PyFaceWipe output was identical across the tested operating systems.

CONCLUSION

PyFaceWipe is a promising multiplatform defacing tool, demonstrating excellent brain voxel preservation and competitive defacing in multiple MRI contrasts, performing favourably against PyDeface. ASL, BOLD, DWI and T2* scans did not produce recognisable 3D renders and hence should not require defacing. Structural volumetry dice scores (≥ 0.98) were higher than previously published FreeSurfer results, except for grey matter which were comparable. The effect is measurable and care should be exercised during studies. ANTS atlas creation showed no significant effect from PyFaceWipe defacing.

Investigating the impact of the regularization parameter on EEG resting-state source reconstruction and functional connectivity using real and simulated data

Accurate EEG source localization is crucial for mapping resting-state network dynamics and it plays a key role in estimating source-level functional connectivity. However, EEG source estimation techniques encounter numerous methodological challenges, with a key one being the selection of the regularization parameter in minimum norm estimation. This choice is particularly intricate because the optimal amount of regularization for EEG source estimation may not align with the requirements of EEG connectivity analysis, highlighting a nuanced trade-off. In this study, we employed a methodological approach to determine the optimal regularization coefficient that yields the most effective reconstruction outcomes across all simulations involving varying signal-to-noise ratios for synthetic EEG signals. To this aim, we considered three resting state networks: the Motor Network, the Visual Network, and the Dorsal Attention Network. The performance was assessed using three metrics, at different regularization parameters: the Region Localization Error, source extension, and source fragmentation. The results were validated using real functional connectivity data. We show that the best estimate of functional connectivity is obtained using 10-2, while 10-1 has to be preferred when source localization only is at target.

A bicoherence approach to analyze multi-dimensional cross-frequency coupling in EEG/MEG data

We introduce a blockwise generalisation of the Antisymmetric Cross-Bicoherence (ACB), a statistical method based on bispectral analysis. The Multi-dimensional ACB (MACB) is an approach that aims at detecting quadratic lagged phase-interactions between vector time series in the frequency domain. Such a coupling can be empirically observed in functional neuroimaging data, e.g., in electro/magnetoencephalographic signals. MACB is invariant under orthogonal trasformations of the data, which makes it independent, e.g., on the choice of the physical coordinate system in the neuro-electromagnetic inverse procedure. In extensive synthetic experiments, we prove that MACB performance is significantly better than that obtained by ACB. Specifically, the shorter the data length, or the higher the dimension of the single data space, the larger the difference between the two methods.

An advance in novel intelligent sensory technologies: From an implicit-tracking perspective of food perception

Food sensory evaluation mainly includes explicit and implicit measurement methods. Implicit measures of consumer perception are gaining significant attention in food sensory and consumer science as they provide effective, subconscious, objective analysis. A wide range of advanced technologies are now available for analyzing physiological and psychological responses, including facial analysis technology, neuroimaging technology, autonomic nervous system technology, and behavioral pattern measurement. However, researchers in the food field often lack systematic knowledge of these multidisciplinary technologies and struggle with interpreting their results. In order to bridge this gap, this review systematically describes the principles and highlights the applications in food sensory and consumer science of facial analysis technologies such as eye tracking, facial electromyography, and automatic facial expression analysis, as well as neuroimaging technologies like electroencephalography, magnetoencephalography, functional magnetic resonance imaging, and functional near-infrared spectroscopy. Furthermore, we critically compare and discuss these advanced implicit techniques in the context of food sensory research and then accordingly propose prospects. Ultimately, we conclude that implicit measures should be complemented by traditional explicit measures to capture responses beyond preference. Facial analysis technologies offer a more objective reflection of sensory perception and attitudes toward food, whereas neuroimaging techniques provide valuable insight into the implicit physiological responses during food consumption. To enhance the interpretability and generalizability of implicit measurement results, further sensory studies are needed. Looking ahead, the combination of different methodological techniques in real-life situations holds promise for consumer sensory science in the field of food research.

The impact of ROI extraction method for MEG connectivity estimation: practical recommendations for the study of resting state data

Magnetoencephalography and electroencephalography (M/EEG) seed-based connectivity analysis requires the extraction of measures from regions of interest (ROI). M/EEG ROI-derived source activity can be treated in different ways. It is possible, for instance, to average each ROI's time series prior to calculating connectivity measures. Alternatively, one can compute connectivity maps for each element of the ROI prior to dimensionality reduction to obtain a single map. The impact of these different strategies on connectivity results is still unclear. Here, we address this question within a large MEG resting state cohort (N=113) and within simulated data. We consider 68 ROIs (Desikan-Kiliany atlas), two measures of connectivity (phase locking value-PLV, and its imaginary counterpart- ciPLV), and three frequency bands (theta 4-8 Hz, alpha 9-12 Hz, beta 15-30 Hz). We compare four extraction methods: (i) mean, or (ii) PCA of the activity within the seed or ROI before computing connectivity, map of the (iii) average, or (iv) maximum connectivity after computing connectivity for each element of the seed. Hierarchical clustering is then applied to compare connectivity outputs across multiple strategies, followed by direct contrasts across extraction methods. Finally, the results are validated by using a set of realistic simulations. We show that ROI-based connectivity maps vary remarkably across strategies in terms of connectivity magnitude and spatial distribution. Dimensionality reduction procedures conducted after computing connectivity are more similar to each-other, while PCA before approach is the most dissimilar to other approaches. Although differences across methods are consistent across frequency bands, they are influenced by the connectivity metric and ROI size. Greater differences were observed for ciPLV than PLV, and in larger ROIs. Realistic simulations confirmed that after aggregation procedures are generally more accurate but have lower specificity (higher rate of false positive connections). Though computationally demanding, after dimensionality reduction strategies should be preferred when higher sensitivity is desired. Given the remarkable differences across aggregation procedures, caution is warranted in comparing results across studies applying different methods.

Optical Co-Registration Method of Triaxial OPM-MEG and MRI

The advent of optically pumped magnetometers (OPMs) facilitates the development of on-scalp magnetoencephalography (MEG). In particular, the triaxial OPM emerged recently, making simultaneous measurements of all three orthogonal components of vector fields possible. The detection of triaxial magnetic fields improves the interference suppression capability and achieves higher source localization accuracy using fewer sensors. The source localization accuracy of MEG is based on the accurate co-registration of MEG and MRI. In this study, we proposed a triaxial co-registration method according to combined principal component analysis and iterative closest point algorithms for use of a flexible cap. A reference phantom with known sensor positions and orientations was designed and constructed to evaluate the accuracy of the proposed method. Experiments showed that the average co-registered position errors of all sensors were approximately 1 mm and average orientation errors were less than 2.5° in the X -and Y orientations and less than 1.6° in the Z orientation. Furthermore, we assessed the influence of co-registration errors on the source localization using simulations. The average source localization error of approximately 1 mm reflects the effectiveness of the co-registration method. The proposed co-registration method facilitates future applications of triaxial sensors on flexible caps.

On electromagnetic head digitization in MEG and EEG

In MEG and EEG studies, the accuracy of the head digitization impacts the co-registration between functional and structural data. The co-registration is one of the major factors that affect the spatial accuracy in MEG/EEG source imaging. Precisely digitized head-surface (scalp) points do not only improve the co-registration but can also deform a template MRI. Such an individualized-template MRI can be used for conductivity modeling in MEG/EEG source imaging if the individual's structural MRI is unavailable. Electromagnetic tracking (EMT) systems (particularly Fastrak, Polhemus Inc., Colchester, VT, USA) have been the most common solution for digitization in MEG and EEG. However, they may occasionally suffer from ambient electromagnetic interference which makes it challenging to achieve (sub-)millimeter digitization accuracy. The current study-(i) evaluated the performance of the Fastrak EMT system under different conditions in MEG/EEG digitization, and (ii) explores the usability of two alternative EMT systems (Aurora, NDI, Waterloo, ON, Canada; Fastrak with a short-range transmitter) for digitization. Tracking fluctuation, digitization accuracy, and robustness of the systems were evaluated in several test cases using test frames and human head models. The performance of the two alternative systems was compared against the Fastrak system. The results showed that the Fastrak system is accurate and robust for MEG/EEG digitization if the recommended operating conditions are met. The Fastrak with the short-range transmitter shows comparatively higher digitization error if digitization is not carried out very close to the transmitter. The study also evinces that the Aurora system can be used for MEG/EEG digitization within a constrained range; however, some modifications would be required to make the system a practical and easy-to-use digitizer. Its real-time error estimation feature can potentially improve digitization accuracy.

Mutual beta power modulation in dyadic entrainment

Across a broad spectrum of interactions, humans exhibit a prominent tendency to synchronize their movements with one another. Traditionally, this phenomenon has been explained from the perspectives of predictive coding or dynamical systems theory. While these theories diverge with respect to whether individuals hold internal models of each other, they both assume a predictive or anticipatory mechanism enabling rhythmic interactions. However, the neural bases underpinning interpersonal synchronization are still a subject under active investigation. Here we provide evidence that the brain relies on a common oscillatory mechanism to pace self-generated rhythmic movements and to track the movements produced by a partner. By performing dual-electroencephalography recordings during a joint finger-tapping task, we identified an oscillatory component in the beta range (∼ 20 Hz), which was significantly modulated by both self-generated and other-generated movement. In conditions where the partners perceived each other, we observed periodic fluctuations of beta power as a function of the reciprocal movement cycles. Crucially, this modulation occurred both in visually and in auditorily coupled conditions, and was accompanied by recurrent periods of dyadic synchronized behavior. Our results show that periodic beta power modulations may be a critical mechanism underlying interpersonal synchronization, possibly enabling mutual predictions between coupled individuals, leading to co-regulation of timing and overt mutual adaptation. Our findings thus provide a potential bridge between influential theories attempting to explain interpersonal coordination, and a concrete connection to its neurophysiological bases.

Infant brain imaging using magnetoencephalography: Challenges, solutions, and best practices

The excellent temporal resolution and advanced spatial resolution of magnetoencephalography (MEG) makes it an excellent tool to study the neural dynamics underlying cognitive processes in the developing brain. Nonetheless, a number of challenges exist when using MEG to image infant populations. There is a persistent belief that collecting MEG data with infants presents a number of limitations and challenges that are difficult to overcome. Due to this notion, many researchers either avoid conducting infant MEG research or believe that, in order to collect high-quality data, they must impose limiting restrictions on the infant or the experimental paradigm. In this article, we discuss the various challenges unique to imaging awake infants and young children with MEG, and share general best-practice guidelines and recommendations for data collection, acquisition, preprocessing, and analysis. The current article is focused on methodology that allows investigators to test the sensory, perceptual, and cognitive capacities of awake and moving infants. We believe that such methodology opens the pathway for using MEG to provide mechanistic explanations for the complex behavior observed in awake, sentient, and dynamically interacting infants, thus addressing core topics in developmental cognitive neuroscience.

Looking through the windows: a study about the dependency of phase-coupling estimates on the data length

OBJECTIVE

Being able to characterize functional connectivity (FC) state dynamics in a real-time setting, such as in brain-computer interface, neurofeedback or closed-loop neurostimulation frameworks, requires the rapid detection of the statistical dependencies that quantify FC in short windows of data. The aim of this study is to characterize, through extensive realistic simulations, the reliability of FC estimation as a function of the data length. In particular, we focused on FC as measured by phase-coupling (PC) of neuronal oscillations, one of the most functionally relevant neural coupling modes.

APPROACH

We generated synthetic data corresponding to different scenarios by varying the data length, the signal-to-noise ratio, the phase difference value, the spectral analysis approach (Hilbert or Fourier) and the fractional bandwidth. We compared seven PC metrics, i.e. imaginary part of phase locking value (PLV), PLV of orthogonalized signals, phase lag index (PLI), debiased weighted PLI, imaginary part of coherency, coherence of orthogonalized signals and lagged coherence.

MAIN RESULTS

Our findings show that, for a signal-to-noise-ratio of at least 10 dB, a data window that contains 5 to 8 cycles of the oscillation of interest (e.g. a 500-800ms window at 10Hz) is generally required to achieve reliable PC estimates. In general, Hilbert-based approaches were associated with higher performance than Fourier-based approaches. Furthermore, the results suggest that, when the analysis is performed in a narrow frequency range, a larger window is required.

SIGNIFICANCE

The achieved results pave the way to the introduction of best-practice guidelines to be followed when a real-time frequency-specific PC assessment is at target.

Co-registration Comparison of On-Scalp Magnetoencephalography and Magnetic Resonance Imaging

Magnetoencephalography (MEG) can non-invasively measure the electromagnetic activity of the brain. A new type of MEG, on-scalp MEG, has attracted the attention of researchers recently. Compared to the conventional SQUID-MEG, on-scalp MEG constructed with optically pumped magnetometers is wearable and has a high signal-to-noise ratio. While the co-registration between MEG and magnetic resonance imaging (MRI) significantly influences the source localization accuracy, co-registration error requires assessment, and quantification. Recent studies have evaluated the co-registration error of on-scalp MEG mainly based on the surface fit error or the repeatability error of different measurements, which do not reflect the true co-registration error. In this study, a three-dimensional-printed reference phantom was constructed to provide the ground truth of MEG sensor locations and orientations relative to MRI. The co-registration performances of commonly used three devices—electromagnetic digitization system, structured-light scanner, and laser scanner—were compared and quantified by the indices of final co-registration errors in the reference phantom and human experiments. Furthermore, the influence of the co-registration error on the performance of source localization was analyzed via simulations. The laser scanner had the best co-registration accuracy (rotation error of 0.23° and translation error of 0.76 mm based on the phantom experiment), whereas the structured-light scanner had the best cost performance. The results of this study provide recommendations and precautions for researchers regarding selecting and using an appropriate device for the co-registration of on-scalp MEG and MRI.

Longitudinal consistency of source-space spectral power and functional connectivity using different magnetoencephalography recording systems

Longitudinal analyses of magnetoencephalography (MEG) data are essential for a full understanding of the pathophysiology of brain diseases and the development of brain activity over time. However, time-dependent factors, such as the recording environment and the type of MEG recording system may affect such longitudinal analyses. We hypothesized that, using source-space analysis, hardware and software differences between two recordings systems may be overcome, with the aim of finding consistent neurophysiological results. We studied eight healthy subjects who underwent three consecutive MEG recordings over 7 years, using two different MEG recordings systems; a 151-channel VSM-CTF system for the first two time points and a 306-channel Elekta Vectorview system for the third time point. We assessed the within (longitudinal) and between-subject (cross-sectional) consistency of power spectra and functional connectivity matrices. Consistency of within-subject spectral power and functional connectivity matrices was good and was not significantly different when using different MEG recording systems as compared to using the same system. Importantly, we confirmed that within-subject consistency values were higher than between-subject values. We demonstrated consistent neurophysiological findings in healthy subjects over a time span of seven years, despite using data recorded on different MEG systems and different implementations of the analysis pipeline.

The Role of Spectral Complexity in Connectivity Estimation

The study of functional connectivity from magnetoecenphalographic (MEG) data consists of quantifying the statistical dependencies among time series describing the activity of different neural sources from the magnetic field recorded outside the scalp. This problem can be addressed by utilizing connectivity measures whose computation in the frequency domain often relies on the evaluation of the cross-power spectrum of the neural time series estimated by solving the MEG inverse problem. Recent studies have focused on the optimal determination of the cross-power spectrum in the framework of regularization theory for ill-posed inverse problems, providing indications that, rather surprisingly, the regularization process that leads to the optimal estimate of the neural activity does not lead to the optimal estimate of the corresponding functional connectivity. Along these lines, the present paper utilizes synthetic time series simulating the neural activity recorded by an MEG device to show that the regularization of the cross-power spectrum is significantly correlated with the signal-to-noise ratio of the measurements and that, as a consequence, this regularization correspondingly depends on the spectral complexity of the neural activity.

Temporal Dynamics on Decoding Target Stimuli in Rapid Serial Visual Presentation using Magnetoencephalography

Rapid serial visual presentation (RSVP) is a high efficient paradigm in brain-computer interface (BCI). Target detection accuracy is the first consideration of RSVP-BCI. But the influence of different frequency bands and time ranges on decoding accuracy are still an open questions. Moreover, the underlying neural dynamic of the rapid target detecting process is still unclear. Methods: This work focused the temporal dynamic of the responses triggered by target stimuli in a static RSVP paradigm using paired structural Magnetic Resonance Imaging (MRI) and magnetoencephalography (MEG) signals with different frequency bands. Multivariate pattern analysis (MVPA) was applied on the MEG signal with different frequency bands and time points after stimuli onset. Cortical neuronal activation estimation technology was also applied to present the temporal-spatial dynamic on cortex surface. Results: The MVPA results showed that the low frequency signals (0.1 - 7 Hz) yield highest decoding accuracy, and the decoding power reached its peak at 0.4 second after target stimuli onset. The cortical neuronal activation method identified the target stimuli triggered regions, like bilateral parahippocampal cortex, precentral gyrus and insula cortex, and the averaged time series were presented.

Multi-channel whole-head OPM-MEG: Helmet design and a comparison with a conventional system

Magnetoencephalography (MEG) is a powerful technique for functional neuroimaging, offering a non-invasive window on brain electrophysiology. MEG systems have traditionally been based on cryogenic sensors which detect the small extracranial magnetic fields generated by synchronised current in neuronal assemblies, however, such systems have fundamental limitations. In recent years, non-cryogenic quantum-enabled sensors, called optically-pumped magnetometers (OPMs), in combination with novel techniques for accurate background magnetic field control, have promised to lift those restrictions offering an adaptable, motion-robust MEG system, with improved data quality, at reduced cost. However, OPM-MEG remains a nascent technology, and whilst viable systems exist, most employ small numbers of sensors sited above targeted brain regions. Here, building on previous work, we construct a wearable OPM-MEG system with ‘whole-head’ coverage based upon commercially available OPMs, and test its capabilities to measure alpha, beta and gamma oscillations. We design two methods for OPM mounting; a flexible (EEG-like) cap and rigid (additively-manufactured) helmet. Whilst both designs allow for high quality data to be collected, we argue that the rigid helmet offers a more robust option with significant advantages for reconstruction of field data into 3D images of changes in neuronal current. Using repeat measurements in two participants, we show signal detection for our device to be highly robust. Moreover, via application of source-space modelling, we show that, despite having 5 times fewer sensors, our system exhibits comparable performance to an established cryogenic MEG device. While significant challenges still remain, these developments provide further evidence that OPM-MEG is likely to facilitate a step change for functional neuroimaging.

Comparison of beamformer implementations for MEG source localization

Beamformers are applied for estimating spatiotemporal characteristics of neuronal sources underlying measured MEG/EEG signals. Several MEG analysis toolboxes include an implementation of a linearly constrained minimum-variance (LCMV) beamformer. However, differences in implementations and in their results complicate the selection and application of beamformers and may hinder their wider adoption in research and clinical use. Additionally, combinations of different MEG sensor types (such as magnetometers and planar gradiometers) and application of preprocessing methods for interference suppression, such as signal space separation (SSS), can affect the results in different ways for different implementations. So far, a systematic evaluation of the different implementations has not been performed. Here, we compared the localization performance of the LCMV beamformer pipelines in four widely used open-source toolboxes (MNE-Python, FieldTrip, DAiSS (SPM12), and Brainstorm) using datasets both with and without SSS interference suppression.

We analyzed MEG data that were i) simulated, ii) recorded from a static and moving phantom, and iii) recorded from a healthy volunteer receiving auditory, visual, and somatosensory stimulation. We also investigated the effects of SSS and the combination of the magnetometer and gradiometer signals. We quantified how localization error and point-spread volume vary with the signal-to-noise ratio (SNR) in all four toolboxes.

When applied carefully to MEG data with a typical SNR (3–15 ​dB), all four toolboxes localized the sources reliably; however, they differed in their sensitivity to preprocessing parameters. As expected, localizations were highly unreliable at very low SNR, but we found high localization error also at very high SNRs for the first three toolboxes while Brainstorm showed greater robustness but with lower spatial resolution. We also found that the SNR improvement offered by SSS led to more accurate localization.

•

Different beamformer implementations are reported to sometimes yield differing source estimates for the same MEG data.

•

We compared beamformers in four major open-source MEG analysis toolboxes.

•

All toolboxes provide consistent and accurate results with 3–15-dB input SNR.

•

However, localization errors are high at very high input SNR for the tested scalar beamformers.

•

We discuss the critical differences between the implementations.

Noninvasive muscle activity imaging using magnetography

Historically, the activity of muscle groups has been evaluated indirectly, through strength and range of movement measurements. While direct evaluation of muscle activity has been achieved using surface or invasive electrodes, the latter methodologies are limited to individual muscles. By contrast, magnetic recording of muscle activity, as illustrated here, offers a noninvasive approach to evaluate the activity of large muscle groups. This includes deep structures such as the heart and intraabdominal musculature such as the iliacus muscle. An unexpected finding was the ability to image activity associated with muscle pain. Thus, this method offers a noninvasive biomarker for muscular pain and, significantly, of neck and back pain, which now relies mainly on the patient’s report of signs and symptoms.

A spectroscopic paradigm has been developed that allows the magnetic field emissions generated by the electrical activity in the human body to be imaged in real time. The growing significance of imaging modalities in biology is evident by the almost exponential increase of their use in research, from the molecular to the ecological level. The method of analysis described here allows totally noninvasive imaging of muscular activity (heart, somatic musculature). Such imaging can be obtained without additional methodological steps such as the use of contrast media.

Non-Invasive Functional-Brain-Imaging with an OPM-based Magnetoencephalography System

A non-invasive functional-brain-imaging system based on optically-pumped-magnetometers (OPM) is presented. The OPM-based magnetoencephalography (MEG) system features 20 OPM channels conforming to the subject's scalp. We have conducted two MEG experiments on three subjects: assessment of somatosensory evoked magnetic field (SEF) and auditory evoked magnetic field (AEF) using our OPM-based MEG system and a commercial MEG system based on superconducting quantum interference devices (SQUIDs). We cross validated the robustness of our system by calculating the distance between the location of the equivalent current dipole (ECD) yielded by our OPM-based MEG system and the ECD location calculated by the commercial SQUID-based MEG system. We achieved sub-centimeter accuracy for both SEF and AEF responses in all three subjects. Due to the proximity (12 mm) of the OPM channels to the scalp, it is anticipated that future OPM-based MEG systems will offer enhanced spatial resolution as they will capture finer spatial features compared to traditional MEG systems employing SQUIDs.

Brain Functional Connectivity Through Phase Coupling of Neuronal Oscillations: A Perspective From Magnetoencephalography

Magnetoencephalography has gained an increasing importance in systems neuroscience thanks to the possibility it offers of unraveling brain networks at time-scales relevant to behavior, i.e., frequencies in the 1-100 Hz range, with sufficient spatial resolution. In the first part of this review, we describe, in a unified mathematical framework, a large set of metrics used to estimate MEG functional connectivity at the same or at different frequencies. The different metrics are presented according to their characteristics: same-frequency or cross-frequency, univariate or multivariate, directed or undirected. We focus on phase coupling metrics given that phase coupling of neuronal oscillations is a putative mechanism for inter-areal communication, and that MEG is an ideal tool to non-invasively detect such coupling. In the second part of this review, we present examples of the use of specific phase methods on real MEG data in the context of resting state, visuospatial attention and working memory. Overall, the results of the studies provide evidence for frequency specific and/or cross-frequency brain circuits which partially overlap with brain networks as identified by hemodynamic-based imaging techniques, such as functional Magnetic Resonance (fMRI). Additionally, the relation of these functional brain circuits to anatomy and to behavior highlights the usefulness of MEG phase coupling in systems neuroscience studies. In conclusion, we believe that the field of MEG functional connectivity has made substantial steps forward in the recent years and is now ready for bringing the study of brain networks to a more mechanistic understanding.

How Sensitive Are Conventional MEG Functional Connectivity Metrics With Sliding Windows to Detect Genuine Fluctuations in Dynamic Functional Connectivity?

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a priori defined ground truths to systematically analyze the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data.