Co-registration of magnetoencephalography with magnetic resonance imaging using bite-bar-based fiducials and surface-matching

Accuracy and precision of navigated transcranial magnetic stimulation

Objective.Transcranial magnetic stimulation (TMS) induces an electric field (E-field) in the cortex. To facilitate stimulation targeting, image-guided neuronavigation systems have been introduced. Such systems track the placement of the coil with respect to the head and visualize the estimated cortical stimulation location on an anatomical brain image in real time. The accuracy and precision of the neuronavigation is affected by multiple factors. Our aim was to analyze how different factors in TMS neuronavigation affect the accuracy and precision of the coil-head coregistration and the estimated E-field.Approach.By performing simulations, we estimated navigation errors due to distortions in magnetic resonance images (MRIs), head-to-MRI registration (landmark- and surface-based registrations), localization and movement of the head tracker, and localization of the coil tracker. We analyzed the effect of these errors on coil and head coregistration and on the induced E-field as determined with simplistic and realistic head models.Main results.Average total coregistration accuracies were in the range of 2.2-3.6 mm and 1°; precision values were about half of the accuracy values. The coregistration errors were mainly due to head-to-MRI registration with average accuracies 1.5-1.9 mm/0.2-0.4° and precisions 0.5-0.8 mm/0.1-0.2° better with surface-based registration. The other major source of error was the movement of the head tracker with average accuracy of 1.5 mm and precision of 1.1 mm. When assessed within an E-field method, the average accuracies of the peak E-field location, orientation, and magnitude ranged between 1.5 and 5.0 mm, 0.9 and 4.8°, and 4.4 and 8.5% across the E-field models studied. The largest errors were obtained with the landmark-based registration. When computing another accuracy measure with the most realistic E-field model as a reference, the accuracies tended to improve from about 10 mm/15°/25% to about 2 mm/2°/5% when increasing realism of the E-field model.Significance.The results of this comprehensive analysis help TMS operators to recognize the main sources of error in TMS navigation and that the coregistration errors and their effect in the E-field estimation depend on the methods applied. To ensure reliable TMS navigation, we recommend surface-based head-to-MRI registration and realistic models for E-field computations.

A unified view on beamformers for M/EEG source reconstruction

Beamforming is a popular method for functional source reconstruction using magnetoencephalography (MEG) and electroencephalography (EEG) data. Beamformers, which were first proposed for MEG more than two decades ago, have since been applied in hundreds of studies, demonstrating that they are a versatile and robust tool for neuroscience. However, certain characteristics of beamformers remain somewhat elusive and there currently does not exist a unified documentation of the mathematical underpinnings and computational subtleties of beamformers as implemented in the most widely used academic open source software packages for MEG analysis (Brainstorm, FieldTrip, MNE, and SPM). Here, we provide such documentation that aims at providing the mathematical background of beamforming and unifying the terminology. Beamformer implementations are compared across toolboxes and pitfalls of beamforming analyses are discussed. Specifically, we provide details on handling rank deficient covariance matrices, prewhitening, the rank reduction of forward fields, and on the combination of heterogeneous sensor types, such as magnetometers and gradiometers. The overall aim of this paper is to contribute to contemporary efforts towards higher levels of computational transparency in functional neuroimaging.

Laminar dynamics of high amplitude beta bursts in human motor cortex

Motor cortical activity in the beta frequency range is one of the strongest and most studied movement-related neural signals. At the single trial level, beta band activity is often characterized by transient, high amplitude, bursting events rather than slowly modulating oscillations. The timing of these bursting events is tightly linked to behavior, suggesting a more dynamic functional role for beta activity than previously believed. However, the neural mechanisms underlying beta bursts in sensorimotor circuits are poorly understood. To address this, we here leverage and extend recent developments in high precision MEG for temporally resolved laminar analysis of burst activity, combined with a neocortical circuit model that simulates the biophysical generators of the electrical currents which drive beta bursts. This approach pinpoints the generation of beta bursts in human motor cortex to distinct excitatory synaptic inputs to deep and superficial cortical layers, which drive current flow in opposite directions. These laminar dynamics of beta bursts in motor cortex align with prior invasive animal recordings within the somatosensory cortex, and suggest a conserved mechanism for somatosensory and motor cortical beta bursts. More generally, we demonstrate the ability for uncovering the laminar dynamics of event-related neural signals in human non-invasive recordings. This provides important constraints to theories about the functional role of burst activity for movement control in health and disease, and crucial links between macro-scale phenomena measured in humans and micro-circuit activity recorded from animal models.

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.

Deep Learning-Based Localization of EEG Electrodes Within MRI Acquisitions

The simultaneous acquisition of electroencephalographic (EEG) signals and functional magnetic resonance images (fMRI) aims to measure brain activity with good spatial and temporal resolution. This bimodal neuroimaging can bring complementary and very relevant information in many cases and in particular for epilepsy. Indeed, it has been shown that it can facilitate the localization of epileptic networks. Regarding the EEG, source localization requires the resolution of a complex inverse problem that depends on several parameters, one of the most important of which is the position of the EEG electrodes on the scalp. These positions are often roughly estimated using fiducial points. In simultaneous EEG-fMRI acquisitions, specific MRI sequences can provide valuable spatial information. In this work, we propose a new fully automatic method based on neural networks to segment an ultra-short echo-time MR volume in order to retrieve the coordinates and labels of the EEG electrodes. It consists of two steps: a segmentation of the images by a neural network, followed by the registration of an EEG template on the obtained detections. We trained the neural network using 37 MR volumes and then we tested our method on 23 new volumes. The results show an average detection accuracy of 99.7% with an average position error of 2.24 mm, as well as 100% accuracy in the labeling.

MEG detection of high frequency oscillations and intracranial-EEG validation in pediatric epilepsy surgery

OBJECTIVE

To assess the feasibility of automatically detecting high frequency oscillations (HFOs) in magnetoencephalography (MEG) recordings in a group of ten paediatric epilepsy surgery patients who had undergone intracranial electroencephalography (iEEG).

METHODS

A beamforming source-analysis method was used to construct virtual sensors and an automatic algorithm was applied to detect HFOs (80-250 Hz). We evaluated the concordance of MEG findings with the sources of iEEG HFOs, the clinically defined seizure onset zone (SOZ), the location of resected brain structures, and with post-operative outcome.

RESULTS

In 8/9 patients there was good concordance between the sources of MEG HFOs and iEEG HFOs and the SOZ. Significantly more HFOs were detected in iEEG relative to MEG t(71) = 2.85, p < .05. There was good concordance between sources of MEG HFOs and the resected area in patients with good and poor outcome, however HFOs were also detected outside of the resected area in patients with poor outcome.

CONCLUSION

Our findings demonstrate the feasibility of automatically detecting HFOs non-invasively in MEG recordings in paediatric patients, and confirm compatibility of results with invasive recordings.

SIGNIFICANCE

This approach provides support for the non-invasive detection of HFOs to aid surgical planning and potentially reduce the need for invasive monitoring, which is pertinent to paediatric patients.

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.

A comparison of automated and manual co-registration for magnetoencephalography

Magnetoencephalography (MEG) is a neuroimaging technique that accurately captures the rapid (sub-millisecond) activity of neuronal populations. Interpretation of functional data from MEG relies upon registration to the participant’s anatomical MRI. The key remaining step is to transform the participant’s MRI into the MEG head coordinate space. Although both automated and manual approaches to co-registration are available, the relative accuracy of two approaches has not been systematically evaluated. The goal of the present study was to compare the accuracy of manual and automated co-registration. Resting MEG and T1-weighted MRI data were collected from 90 participants. Automated and manual co-registration were performed on the same subjects, and the inter-method reliability of the two methods assessed using the intra-class correlation. Median co-registration error for both methods was within acceptable limits. Inter-method reliability was in the “good” range for co-registration error, and the “good” to “excellent” range for translation and rotation. These results suggest that the output of the automated co-registration procedure is comparable to that achieved using manual co-registration.

Mapping language networks and their association with verbal abilities in paediatric epilepsy using MEG and graph analysis

Recent theoretical models of language have emphasised the importance of integration within distributed networks during language processing. This is particularly relevant to young patients with epilepsy, as the topology of the functional network and its dynamics may be altered by the disease, resulting in reorganisation of functional language networks. Thus, understanding connectivity within the language network in patients with epilepsy could provide valuable insights into healthy and pathological brain function, particularly when combined with clinical correlates. The objective of this study was to investigate interactions within the language network in a paediatric population of epilepsy patients using measures of MEG phase synchronisation and graph-theoretical analysis, and to examine their association with language abilities. Task dependent increases in connectivity were observed in fronto-temporal networks during verb generation across a group of 22 paediatric patients (9 males and 13 females; mean age 14 years). Differences in network connectivity were observed between patients with typical and atypical language representation and between patients with good and poor language abilities. In addition, node centrality in left frontal and temporal regions was significantly associated with language abilities, where patients with good language abilities had significantly higher node centrality within inferior frontal and superior temporal regions of the left hemisphere, compared to patients with poor language abilities. Our study is one of the first to apply task-based measures of MEG network synchronisation in paediatric epilepsy, and we propose that these measures of functional connectivity and node centrality could be used as tools to identify critical regions of the language network prior to epilepsy surgery.

Estimates of cortical column orientation improve MEG source inversion

Determining the anatomical source of brain activity non-invasively measured from EEG or MEG sensors is challenging. In order to simplify the source localization problem, many techniques introduce the assumption that current sources lie on the cortical surface. Another common assumption is that this current flow is orthogonal to the cortical surface, thereby approximating the orientation of cortical columns. However, it is not clear which cortical surface to use to define the current source locations, and normal vectors computed from a single cortical surface may not be the best approximation to the orientation of cortical columns. We compared three different surface location priors and five different approaches for estimating dipole vector orientation, both in simulations and visual and motor evoked MEG responses. We show that models with source locations on the white matter surface and using methods based on establishing correspondences between white matter and pial cortical surfaces dramatically outperform models with source locations on the pial or combined pial/white surfaces and which use methods based on the geometry of a single cortical surface in fitting evoked visual and motor responses. These methods can be easily implemented and adopted in most M/EEG analysis pipelines, with the potential to significantly improve source localization of evoked responses.

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.

Automatic Spatial Calibration of Ultra-Low-Field MRI for High-Accuracy Hybrid MEG-MRI

With a hybrid magnetoencephalography (MEG)-MRI device that uses the same sensors for both modalities, the co-registration of MRI and MEG data can be replaced by an automatic calibration step. Based on the highly accurate signal model of ultra-low-field (ULF) MRI, we introduce a calibration method that eliminates the error sources of traditional co-registration. The signal model includes complex sensitivity profiles of the superconducting pickup coils. In the ULF MRI, the profiles are independent of the sample and therefore well-defined. In the most basic form, the spatial information of the profiles, captured in parallel ULF-MR acquisitions, is used to find the exact coordinate transformation required. We assessed our calibration method by simulations assuming a helmet-shaped pickup-coil-array geometry. Using a carefully constructed objective function and sufficient approximations, even with low-SNR images, sub-voxel and sub-millimeter calibration accuracy were achieved. After the calibration, distortion-free MRI and high spatial accuracy for MEG source localization can be achieved. For an accurate sensor-array geometry, the co-registration and associated errors are eliminated, and the positional error can be reduced to a negligible level.

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

Co-registration between structural head images and functional MEG data is needed for anatomically-informed MEG data analysis. Despite the efforts to minimize the co-registration error, conventional landmark- and surface-based strategies for co-registering head and MEG device coordinates achieve an accuracy of typically 5-10 mm. Recent advances in instrumentation and technical solutions, such as the development of hybrid ultra-low-field (ULF) MRI-MEG devices or the use of 3D-printed individualized foam head-casts, promise unprecedented co-registration accuracy, i.e., 2 mm or better. In the present study, we assess through simulations the impact of such an improved co-registration on MEG connectivity analysis. We generated synthetic MEG recordings for pairs of connected cortical sources with variable locations. We then assessed the capability to reconstruct source-level connectivity from these recordings for 0-15-mm co-registration error, three levels of head modeling detail (one-, three- and four-compartment models), two source estimation techniques (linearly constrained minimum-variance beamforming and minimum-norm estimation MNE) and five separate connectivity metrics (imaginary coherency, phase-locking value, amplitude-envelope correlation, phase-slope index and frequency-domain Granger causality). We found that beamforming can better take advantage of an accurate co-registration than MNE. Specifically, when the co-registration error was smaller than 3 mm, the relative error in connectivity estimates was down to one-third of that observed with typical co-registration errors. MNE provided stable results for a wide range of co-registration errors, while the performance of beamforming rapidly degraded as the co-registration error increased. Furthermore, we found that even moderate co-registration errors (>6 mm, on average) essentially decrease the difference of four- and three- or one-compartment models. Hence, a precise co-registration is important if one wants to take full advantage of highly accurate head models for connectivity analysis. We conclude that an improved co-registration will be beneficial for reliable connectivity analysis and effective use of highly accurate head models in future MEG connectivity studies.

Optical Co-registration of MRI and On-scalp MEG

To estimate the neural generators of magnetoencephalographic (MEG) signals, MEG data have to be co-registered with an anatomical image, typically an MR image. Optically-pumped magnetometers (OPMs) enable the construction of on-scalp MEG systems providing higher sensitivity and spatial resolution than conventional SQUID-based MEG systems. We present a co-registration method that can be applied to on-scalp MEG systems, regardless of the number of sensors. We apply a structured-light scanner to create a surface mesh of the subject’s head and the sensor array, which we fit to the MR image. We quantified the reproducibility of the mesh and localised current dipoles with a phantom. Additionally, we measured somatosensory evoked fields (SEFs) to median nerve stimulation and compared the dipole positions between on-scalp and SQUID-based systems. The scanner reproduced the head surface with <1 mm error. Phantom dipoles were localised with 2.1 mm mean error. SEF dipoles corresponding to the P35m response for OPMs were well localised to the somatosensory cortex, while SQUID dipoles for two subjects were erroneously localised to the motor cortex. The developed co-registration method is inexpensive, fast and can easily be applied to on-scalp MEG. It is more convenient than traditional co-registration methods while also being more accurate.

MEG Assessment of Expressive Language in Children Evaluated for Epilepsy Surgery

Establishing language dominance is an important step in the presurgical evaluation of patients with refractory epilepsy. In the absence of a universally accepted gold-standard non-invasive method to determine language dominance in the preoperative assessment, a range of tools and methodologies have recently received attention. When applied to pediatric age, many of the proposed methods, such as functional magnetic resonance imaging (fMRI), may present some challenges due to the time-varying effects of epileptogenic lesions and of on-going seizures on maturational phenomena. Magnetoencephalography (MEG) has the advantage of being insensitive to the distortive effects of anatomical lesions on brain microvasculature and to differences in the metabolism or vascularization of the developing brain and also provides a less intimidating recording environment for younger children. In this study we investigated the reliability of lateralized synchronous cortical activation during a verb generation task in a group of 28 children (10 males and 18 females, mean age 12 years) with refractory epilepsy who were evaluated for epilepsy surgery. The verb generation task was associated with significant decreases in beta oscillatory power (13–30 Hz) in frontal and temporal lobes. The MEG data were compared with other available presurgical non-invasive data including cortical stimulation, neuropsychological and fMRI data on language lateralization where available. We found that the lateralization of MEG beta power reduction was concordant with language dominance determined by one or more different assessment methods (i.e. cortical stimulation mapping, neuropsychological, fMRI or post-operative data) in 89% of patients. Our data suggest that qualitative hemispheric differences in task-related changes of spectral power could offer a promising insight into the contribution of dominant and non-dominant hemispheres in language processing and may help to characterize the specialization and lateralization of language processes in children.

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

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

Quality assessment of MEG-to-MRI coregistrations

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

An evaluation of kurtosis beamforming in magnetoencephalography to localize the epileptogenic zone in drug resistant epilepsy patients

•

Objective localizations of interictal spikes using a kurtosis beamformer.

•

Kurtosis Beamforming can provide confidence to scattered dipoles.

•

Kurtosis beamforming can assist in localizing the epileptogenic zone.

Objective

Kurtosis beamforming is a useful technique for analysing magnetoencephalograpy (MEG) data containing epileptic spikes. However, the implementation varies and few studies measure concordance with subsequently resected areas. We evaluated kurtosis beamforming as a means of localizing spikes in drug-resistant epilepsy patients.

Methods

We retrospectively applied kurtosis beamforming to MEG recordings of 22 epilepsy patients that had previously been analysed using equivalent current dipole (ECD) fitting. Virtual electrodes were placed in the kurtosis volumetric peaks and visually inspected to select a candidate source. The candidate sources were compared to the ECD localizations and resection areas.

Results

The kurtosis beamformer produced interpretable localizations in 18/22 patients, of which the candidate source coincided with the resection lobe in 9/13 seizure-free patients and in 3/5 patients with persistent seizures. The sublobar accuracy of the kurtosis beamformer with respect to the resection zone was higher than ECD (56% and 50%, respectively), however, ECD resulted in a higher lobar accuracy (75%, 67%).

Conclusions

Kurtosis beamforming may provide additional value when spikes are not clearly discernible on the sensors and support ECD localizations when dipoles are scattered.

Significance

Kurtosis beamforming should be integrated with existing clinical protocols to assist in localizing the epileptogenic zone.

Optimized 3D co-registration of ultra-low-field and high-field magnetic resonance images

The prototypes of ultra-low-field (ULF) MRI scanners developed in recent years represent new, innovative, cost-effective and safer systems, which are suitable to be integrated in multi-modal (Magnetoencephalography and MRI) devices. Integrated ULF-MRI and MEG scanners could represent an ideal solution to obtain functional (MEG) and anatomical (ULF MRI) information in the same environment, without errors that may limit source reconstruction accuracy. However, the low resolution and signal-to-noise ratio (SNR) of ULF images, as well as their limited coverage, do not generally allow for the construction of an accurate individual volume conductor model suitable for MEG localization. Thus, for practical usage, a high-field (HF) MRI image is also acquired, and the HF-MRI images are co-registered to the ULF-MRI ones. We address here this issue through an optimized pipeline (SWIM—Sliding WIndow grouping supporting Mutual information). The co-registration is performed by an affine transformation, the parameters of which are estimated using Normalized Mutual Information as the cost function, and Adaptive Simulated Annealing as the minimization algorithm. The sub-voxel resolution of the ULF images is handled by a sliding-window approach applying multiple grouping strategies to down-sample HF MRI to the ULF-MRI resolution. The pipeline has been tested on phantom and real data from different ULF-MRI devices, and comparison with well-known toolboxes for fMRI analysis has been performed. Our pipeline always outperformed the fMRI toolboxes (FSL and SPM). The HF–ULF MRI co-registration obtained by means of our pipeline could lead to an effective integration of ULF MRI with MEG, with the aim of improving localization accuracy, but also to help exploit ULF MRI in tumor imaging.

Modulation of Neural Oscillatory Activity during Dynamic Face Processing

Various neuroimaging and neurophysiological methods have been used to examine neural activation patterns in response to faces. However, much of previous research has relied on static images of faces, which do not allow a complete description of the temporal structure of face-specific neural activities to be made. More recently, insights are emerging from fMRI studies about the neural substrates that underpin our perception of naturalistic dynamic face stimuli, but the temporal and spectral oscillatory activity associated with processing dynamic faces has yet to be fully characterized. Here, we used MEG and beamformer source localization to examine the spatiotemporal profile of neurophysiological oscillatory activity in response to dynamic faces. Source analysis revealed a number of regions showing enhanced activation in response to dynamic relative to static faces in the distributed face network, which were spatially coincident with regions that were previously identified with fMRI. Furthermore, our results demonstrate that perception of realistic dynamic facial stimuli activates a distributed neural network at varying time points facilitated by modulations in low-frequency power within alpha and beta frequency ranges (8-30 Hz). Naturalistic dynamic face stimuli may provide a better means of representing the complex nature of perceiving facial expressions in the real world, and neural oscillatory activity can provide additional insights into the associated neural processes.

Consistency of magnetoencephalographic functional connectivity and network reconstruction using a template versus native MRI for co-registration

Introduction

Studies using functional connectivity and network analyses based on magnetoencephalography (MEG) with source localization are rapidly emerging in neuroscientific literature. However, these analyses currently depend on the availability of costly and sometimes burdensome individual MR scans for co‐registration. We evaluated the consistency of these measures when using a template MRI, instead of native MRI, for the analysis of functional connectivity and network topology.

Methods

Seventeen healthy participants underwent resting‐state eyes‐closed MEG and anatomical MRI. These data were projected into source space using an atlas‐based peak voxel and a centroid beamforming approach either using (1) participants’ native MRIs or (2) the Montreal Neurological Institute's template. For both methods, time series were reconstructed from 78 cortical atlas regions. Relative power was determined in six classical frequency bands per region and globally averaged. Functional connectivity (phase lag index) between each pair of regions was calculated. The adjacency matrices were then used to reconstruct functional networks, of which regional and global metrics were determined. Intraclass correlation coefficients were calculated and Bland–Altman plots were made to quantify the consistency and potential bias of the use of template versus native MRI.

Results

Co‐registration with the template yielded largely consistent relative power, connectivity, and network estimates compared to native MRI.

Discussion

These findings indicate that there is no (systematic) bias or inconsistency between template and native MRI co‐registration of MEG. They open up possibilities for retrospective and prospective analyses to MEG datasets in the general population that have no native MRIs available. Hum Brain Mapp, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.Hum Brain Mapp 39:104–119, 2018. © 2017 Wiley Periodicals, Inc.

Photogrammetry-Based Head Digitization for Rapid and Accurate Localization of EEG Electrodes and MEG Fiducial Markers Using a Single Digital SLR Camera

The performance of EEG source reconstruction has benefited from the increasing use of advanced head modeling techniques that take advantage of MRI together with the precise positions of the recording electrodes. The prevailing technique for registering EEG electrode coordinates involves electromagnetic digitization. However, the procedure adds several minutes to experiment preparation and typical digitizers may not be accurate enough for optimal source reconstruction performance (Dalal et al., 2014). Here, we present a rapid, accurate, and cost-effective alternative method to register EEG electrode positions, using a single digital SLR camera, photogrammetry software, and computer vision techniques implemented in our open-source toolbox, janus3D. Our approach uses photogrammetry to construct 3D models from multiple photographs of the participant's head wearing the EEG electrode cap. Electrodes are detected automatically or semi-automatically using a template. The rigid facial features from these photo-based models are then surface-matched to MRI-based head reconstructions to facilitate coregistration to MRI space. This method yields a final electrode coregistration error of 0.8 mm, while a standard technique using an electromagnetic digitizer yielded an error of 6.1 mm. The technique furthermore reduces preparation time, and could be extended to a multi-camera array, which would make the procedure virtually instantaneous. In addition to EEG, the technique could likewise capture the position of the fiducial markers used in magnetoencephalography systems to register head position.

Using generative models to make probabilistic statements about hippocampal engagement in MEG

Magnetoencephalography (MEG) enables non-invasive real time characterization of brain activity. However, convincing demonstrations of signal contributions from deeper sources such as the hippocampus remain controversial and are made difficult by its depth, structural complexity and proximity to neocortex. Here, we demonstrate a method for quantifying hippocampal engagement probabilistically using simulated hippocampal activity and realistic anatomical and electromagnetic source modelling. We construct two generative models, one which supports neuronal current flow on the cortical surface, and one which supports neuronal current flow on both the cortical and hippocampal surface. Using Bayesian model comparison, we then infer which of the two models provides a more likely explanation of the dataset at hand. We also carry out a set of control experiments to rule out bias, including simulating medial temporal lobe sources to assess the risk of falsely positive results, and adding different types of displacements to the hippocampal portion of the mesh to test for anatomical specificity of the results. In addition, we test the robustness of this inference by adding co-registration error and sensor level noise. We find that the model comparison framework is sensitive to hippocampal activity when co-registration error is <3 mm and the sensor-level signal-to-noise ratio (SNR) is >-20 dB. These levels of co-registration error and SNR can now be achieved empirically using recently developed subject-specific head-casts.

The auditory evoked-gamma response and its relation with the N1m

This study explored the patterns of oscillatory activity that underpin the N1m auditory evoked response. Evoked gamma activity is a small and relatively rarely-reported component of the auditory evoked response, and the objective of this work was to determine how this component relates to the larger and more prolonged changes in lower frequency bands. An event-related beamformer analysis of MEG data from monaural click stimulation was used to reconstruct volumetric images and virtual electrode time series. Group analysis of localisations showed that activity in the gamma band originated from a source that was more medial than those for activity in the theta-to-beta band, and virtual-electrode analysis showed that the source of the gamma activity could be statistically dissociated from the lower-frequency response. These findings are in accordance with separate functional roles for the activity in each frequency band, and provide evidence that the oscillatory activity that underpins the auditory evoked response may contain important information about the physiological basis of the macroscopic signals recorded by MEG in response to auditory stimulation.

Flexible head-casts for high spatial precision MEG

BACKGROUND

In combination with magnetoencephalographic (MEG) data, accurate knowledge of the brain's structure and location provide a principled way of reconstructing neural activity with high temporal resolution. However, measuring the brain's location is compromised by head movement during scanning, and by fiducial-based co-registration with magnetic resonance imaging (MRI) data. The uncertainty from these two factors introduces errors into the forward model and limit the spatial resolution of the data.

NEW METHOD

We present a method for stabilizing and reliably repositioning the head during scanning, and for co-registering MRI and MEG data with low error.

RESULTS

Using this new flexible and comfortable subject-specific head-cast prototype, we find within-session movements of <0.25mm and between-session repositioning errors around 1mm.

COMPARISON WITH EXISTING METHOD(S)

This method is an improvement over existing methods for stabilizing the head or correcting for location shifts on- or off-line, which still introduce approximately 5mm of uncertainty at best (Adjamian et al., 2004; Stolk et al., 2013; Whalen et al., 2008). Further, the head-cast design presented here is more comfortable, safer, and easier to use than the earlier 3D printed prototype, and give slightly lower co-registration errors (Troebinger et al., 2014b).

CONCLUSIONS

We provide an empirical example of how these head-casts impact on source level reproducibility. Employment of the individual flexible head-casts for MEG recordings provide a reliable method of safely stabilizing the head during MEG recordings, and for co-registering MRI anatomical images to MEG functional data.

Minimizing noise in pediatric task-based functional MRI; Adolescents with developmental disabilities and typical development

Functional Magnetic Resonance Imaging (fMRI) represents a powerful tool with which to examine brain functioning and development in typically developing pediatric groups as well as children and adolescents with clinical disorders. However, fMRI data can be highly susceptible to misinterpretation due to the effects of excessive levels of noise, often related to head motion. Imaging children, especially with developmental disorders, requires extra considerations related to hyperactivity, anxiety and the ability to perform and maintain attention to the fMRI paradigm. We discuss a number of methods that can be employed to minimize noise, in particular movement-related noise. To this end we focus on strategies prior to, during and following the data acquisition phase employed primarily within our own laboratory. We discuss the impact of factors such as experimental design, screening of potential participants and pre-scan training on head motion in our adolescents with developmental disorders and typical development. We make some suggestions that may minimize noise during data acquisition itself and finally we briefly discuss some current processing techniques that may help to identify and remove noise in the data. Many advances have been made in the field of pediatric imaging, particularly with regard to research involving children with developmental disorders. Mindfulness of issues such as those discussed here will ensure continued progress and greater consistency across studies.

Localising the auditory N1m with event-related beamformers: localisation accuracy following bilateral and unilateral stimulation

The auditory evoked N1m-P2m response complex presents a challenging case for MEG source-modelling, because symmetrical, phase-locked activity occurs in the hemispheres both contralateral and ipsilateral to stimulation. Beamformer methods, in particular, can be susceptible to localisation bias and spurious sources under these conditions. This study explored the accuracy and efficiency of event-related beamformer source models for auditory MEG data under typical experimental conditions: monaural and diotic stimulation; and whole-head beamformer analysis compared to a half-head analysis using only sensors from the hemisphere contralateral to stimulation. Event-related beamformer localisations were also compared with more traditional single-dipole models. At the group level, the event-related beamformer performed equally well as the single-dipole models in terms of accuracy for both the N1m and the P2m, and in terms of efficiency (number of successful source models) for the N1m. The results yielded by the half-head analysis did not differ significantly from those produced by the traditional whole-head analysis. Any localisation bias caused by the presence of correlated sources is minimal in the context of the inter-individual variability in source localisations. In conclusion, event-related beamformers provide a useful alternative to equivalent-current dipole models in localisation of auditory evoked responses.

Source-Space Cross-Frequency Amplitude-Amplitude Coupling in Tinnitus

The thalamocortical dysrhythmia (TCD) model has been influential in the development of theoretical explanations for the neurological mechanisms of tinnitus. It asserts that thalamocortical oscillations lock a region in the auditory cortex into an ectopic slow-wave theta rhythm (4-8 Hz). The cortical area surrounding this region is hypothesized to generate abnormal gamma (>30 Hz) oscillations ("edge effect") giving rise to the tinnitus percept. Consequently, the model predicts enhanced cross-frequency coherence in a broad range between theta and gamma. In this magnetoencephalography study involving tinnitus and control cohorts, we investigated this prediction. Using beamforming, cross-frequency amplitude-amplitude coupling (AAC) was computed within the auditory cortices for frequencies (f1, f2) between 2 and 80 Hz. We find the AAC signal to decompose into two distinct components at low (f1, f2 < 30 Hz) and high (f1, f2 > 30 Hz) frequencies, respectively. Studying the correlation of AAC with several key covariates (age, hearing level (HL), tinnitus handicap and duration, and HL at tinnitus frequency), we observe a statistically significant association between age and low-frequency AAC. Contrary to the TCD predictions, however, we do not find any indication of statistical differences in AAC between tinnitus and controls and thus no evidence for the predicted enhancement of cross-frequency coupling in tinnitus.

Source space estimation of oscillatory power and brain connectivity in tinnitus

Tinnitus is the perception of an internally generated sound that is postulated to emerge as a result of structural and functional changes in the brain. However, the precise pathophysiology of tinnitus remains unknown. Llinas’ thalamocortical dysrhythmia model suggests that neural deafferentation due to hearing loss causes a dysregulation of coherent activity between thalamus and auditory cortex. This leads to a pathological coupling of theta and gamma oscillatory activity in the resting state, localised to the auditory cortex where normally alpha oscillations should occur. Numerous studies also suggest that tinnitus perception relies on the interplay between auditory and non-auditory brain areas. According to the Global Brain Model, a network of global fronto—parietal—cingulate areas is important in the generation and maintenance of the conscious perception of tinnitus. Thus, the distress experienced by many individuals with tinnitus is related to the top—down influence of this global network on auditory areas. In this magnetoencephalographic study, we compare resting-state oscillatory activity of tinnitus participants and normal-hearing controls to examine effects on spectral power as well as functional and effective connectivity. The analysis is based on beamformer source projection and an atlas-based region-of-interest approach. We find increased functional connectivity within the auditory cortices in the alpha band. A significant increase is also found for the effective connectivity from a global brain network to the auditory cortices in the alpha and beta bands. We do not find evidence of effects on spectral power. Overall, our results provide only limited support for the thalamocortical dysrhythmia and Global Brain models of tinnitus.

GABA-mediated changes in inter-hemispheric beta frequency activity in early-stage Parkinson's disease

•

In PD, contralateral M1 showed greater beta power than ipsilateral M1.

•

Zolpidem reduced contralateral beta power while ipsilateral power was increased.

•

This resulted in a hemispheric power ratio that approached parity.

•

Changes were reflected in pre-movement desynchronization and post-movement rebound.

•

These changes underlie the symptomatic improvements afforded by zolpidem.

In Parkinson’s disease (PD), elevated beta (15–35 Hz) power in subcortical motor networks is widely believed to promote aspects of PD symptomatology, moreover, a reduction in beta power and coherence accompanies symptomatic improvement following effective treatment with l-DOPA. Previous studies have reported symptomatic improvements that correlate with changes in cortical network activity following GABA_A receptor modulation. In this study we have used whole-head magnetoencephalography to characterize neuronal network activity, at rest and during visually cued finger abductions, in unilaterally symptomatic PD and age-matched control participants. Recordings were then repeated following administration of sub-sedative doses of the hypnotic drug zolpidem (0.05 mg/kg), which binds to the benzodiazepine site of the GABA_A receptor. A beamforming based ‘virtual electrode’ approach was used to reconstruct oscillatory power in the primary motor cortex (M1), contralateral and ipsilateral to symptom presentation in PD patients or dominant hand in control participants. In PD patients, contralateral M1 showed significantly greater beta power than ipsilateral M1. Following zolpidem administration contralateral beta power was significantly reduced while ipsilateral beta power was significantly increased resulting in a hemispheric power ratio that approached parity. Furthermore, there was highly significant correlation between hemispheric beta power ratio and Unified Parkinson’s Disease Rating Scale (UPDRS). The changes in contralateral and ipsilateral beta power were reflected in pre-movement beta desynchronization and the late post-movement beta rebound. However, the absolute level of movement-related beta desynchronization was not altered. These results show that low-dose zolpidem not only reduces contralateral beta but also increases ipsilateral beta, while rebalancing the dynamic range of M1 network oscillations between the two hemispheres. These changes appear to underlie the symptomatic improvements afforded by low-dose zolpidem.

Qualitative evaluation of fiducial markers for radiotherapy imaging

PURPOSE

To evaluate visibility, artifacts, and distortions of various commercial markers in magnetic resonance imaging (MRI), computer tomography (CT), and ultrasound imaging used for radiotherapy planning and treatment guidance.

METHODS

We compare 2 solid gold markers, 4 gold coils, and 1 polymer marker from 3 vendors. Imaging modalities used were 3-T and 1.5-T GE MRIs, Siemens Sequoia 512 Ultrasound, Phillips Big Bore CT, Varian Trilogy linear accelerator (cone-beam CT [CBCT], on-board imager kilovoltage [OBI-kV], electronic portal imaging device megavoltage [EPID-MV]), and Medtronic O-ARM CBCT. Markers were imaged in a 30 × 30 × 10 cm(3) custom bolus phantom. In one experiment, Surgilube was used around the markers to reduce air gaps. Images were saved in Digital Imaging and Communications in Medicine (DICOM) format and analyzed using an in-house software. Profiles across the markers were used for objective comparison of the markers' signals. The visibility and artifacts/distortions produced by each marker were assessed qualitatively and quantitatively.

RESULTS

All markers are visible in CT, CBCT, OBI-kV, and ultrasound. Gold markers below 0.75 mm in diameter are not visible in EPID-MV images. The larger the markers, the more CT and CBCT image artifacts there are, yet the degree of the artifact depends on scan parameters and the scanner itself. Visibility of gold coils of 0.75 mm diameter or larger is comparable across all imaging modalities studied. The polymer marker causes minimal artifacts in CT and CBCT but has poor visibility in EPID-MV. Gold coils of 0.5 mm exhibit poor visibility in MRI and EPID-MV due to their small size. Gold markers are more visible in 3-T T1 gradient-recalled echo than in 1.5-T T1 fast spin-echo, depending on the scan sequence. In this study, all markers are clearly visible on ultrasound.

CONCLUSION

All gold markers are visible in CT, CBCT, kV, and ultrasound; however, only the large diameter markers are visible in MV. When MR and EPID-MV imagers are used, the selection of fiducial markers is not straightforward. For hybrid kV/MV image-guided radiotherapy imaging, larger diameter markers are suggested. If using kV imaging alone, smaller sized markers may be used in smaller sized patients in order to reduce artifacts. Only larger diameter gold markers are visible across all imaging modalities.

Magnetoencephalography in the study of epilepsy and consciousness

The neural bases of altered consciousness in patients with epilepsy during seizures and at rest have raised significant interest in the last decade. This exponential growth has been supported by the parallel development of techniques and methods to investigate brain function noninvasively with unprecedented spatial and temporal resolution. In this article, we review the contribution of magnetoencephalography to deconvolve the bioelectrical changes associated with impaired consciousness during seizures. We use data collected from a patient with refractory absence seizures to discuss how spike-wave discharges are associated with perturbations in optimal connectivity within and between brain regions and discuss indirect evidence to suggest that this phenomenon might explain the cognitive deficits experienced during prolonged 3/s spike-wave discharges.

High precision anatomy for MEG

Precise MEG estimates of neuronal current flow are undermined by uncertain knowledge of the head location with respect to the MEG sensors. This is either due to head movements within the scanning session or systematic errors in co-registration to anatomy. Here we show how such errors can be minimized using subject-specific head-casts produced using 3D printing technology. The casts fit the scalp of the subject internally and the inside of the MEG dewar externally, reducing within session and between session head movements. Systematic errors in matching to MRI coordinate system are also reduced through the use of MRI-visible fiducial markers placed on the same cast. Bootstrap estimates of absolute co-registration error were of the order of 1 mm. Estimates of relative co-registration error were < 1.5 mm between sessions. We corroborated these scalp based estimates by looking at the MEG data recorded over a 6 month period. We found that the between session sensor variability of the subject's evoked response was of the order of the within session noise, showing no appreciable noise due to between-session movement. Simulations suggest that the between-session sensor level amplitude SNR improved by a factor of 5 over conventional strategies. We show that at this level of coregistration accuracy there is strong evidence for anatomical models based on the individual rather than canonical anatomy; but that this advantage disappears for errors of greater than 5 mm. This work paves the way for source reconstruction methods which can exploit very high SNR signals and accurate anatomical models; and also significantly increases the sensitivity of longitudinal studies with MEG.

•

We introduce MEG coregistration scheme using 3D printed subject specific head casts.

•

Using this scheme reduces relative/absolute coregistration errors to 1–2 mm levels.

•

The ability to reposition between sessions results in high SNR data sets.

•

At this level of coregistration error, we see a clear benefit for using individual cortical meshes.

Gamma oscillatory amplitude encodes stimulus intensity in primary somatosensory cortex

Gamma oscillations have previously been linked to pain perception and it has been hypothesized that they may have a potential role in encoding pain intensity. Stimulus response experiments have reported an increase in activity in the primary somatosensory cortex (SI) with increasing stimulus intensity, but the specific role of oscillatory dynamics in this change in activation remains unclear. In this study, Magnetoencephalography (MEG) was used to investigate the changes in cortical oscillations during four different intensities of a train of electrical stimuli to the right index finger, ranging from low sensation to strong pain. In those participants showing changes in evoked oscillatory gamma in SI during stimulation, the strength of the gamma power was found to increase with increasing stimulus intensity at both pain and sub-pain thresholds. These results suggest that evoked gamma oscillations in SI are not specific to pain but may have a role in encoding somatosensory stimulus intensity.

Oscillatory beta activity mediates neuroplastic effects of motor cortex stimulation in humans

Continuous theta burst stimulation (cTBS) is a repetitive transcranial magnetic stimulation protocol that can inhibit human motor cortex (M1) excitability and impair movement for ≤ 1 h. While offering valuable insights into brain function and potential therapeutic benefits, these neuroplastic effects are highly variable between individuals. The source of this variability, and the electrophysiological mechanisms underlying the inhibitory after-effects, are largely unknown. In this regard, oscillatory activity at beta frequency (15-35 Hz) is of particular interest as it is elevated in motor disorders such as Parkinson's disease and modulated during the generation of movements. Here, we used a source-level magnetoencephalography approach to investigate the hypothesis that the presence of neuroplastic effects following cTBS is associated with concurrent changes in oscillatory M1 beta activity. M1 cortices were localized with a synthetic aperture magnetometry beamforming analysis of visually cued index finger movements. Virtual electrode analysis was used to reconstruct the spontaneous and movement-related oscillatory activity in bilateral M1 cortices, before and from 10 to 45 min after cTBS. We demonstrate that 40 s of cTBS applied over left M1 reduced corticospinal excitability in the right index finger of 8/16 participants. In these responder participants only, cTBS increased the power of the spontaneous beta oscillations in stimulated M1 and delayed reaction times in the contralateral index finger. No further changes were observed in the latency or power of movement-related beta oscillations. These data provide insights into the electrophysiological mechanisms underlying cTBS-mediated impairment of motor function and demonstrate the association between spontaneous oscillatory beta activity in M1 and the inhibition of motor function.

SQUID-detected ultra-low field MRI

MRI remains the premier method for non-invasive imaging of soft-tissue. Since the first demonstration of ULF MRI the trend has been towards ever higher magnetic fields. This is because the signal, and efficiency of Faraday detectors, increases with ever higher magnetic fields and corresponding Larmor frequencies. Nevertheless, there are many compelling reasons to continue to explore MRI at much weaker magnetic fields, the so-called ultra-low field or (ULF) regime. In the past decade many excellent proof-of-concept demonstrations of ULF MRI have been made. These include combined MRI and magnetoencephalography, imaging in the presence of metal, unique tissue contrast, and implementation in situations where a high magnetic field is simply impractical. These demonstrations have routinely used pulsed pre-polarization (at magnetic fields from ~10 to 100 mT) followed by read-out in a much weaker (1-100 μT) magnetic fields using the ultra-sensitive Superconducting Quantum Interference Device (SQUID) sensor. Even with pre-polarization and SQUID detection, ULF MRI suffers from many challenges associated with lower magnetization (i.e. signal) and inherently long acquisition times compared to conventional >1 T MRI. These are fundamental limitations imposed by the low measurement and gradient fields used. In this review article we discuss some of the techniques, potential applications, and inherent challenges of ULF MRI.

SQUID-detected ultra-low field MRI

MRI remains the premier method for non-invasive imaging of soft-tissue. Since the first demonstration of ULF MRI the trend has been towards ever higher magnetic fields. This is because the signal, and efficiency of Faraday detectors, increases with ever higher magnetic fields and corresponding Larmor frequencies. Nevertheless, there are many compelling reasons to continue to explore MRI at much weaker magnetic fields, the so-called ultra-low field or (ULF) regime. In the past decade many excellent proof-of-concept demonstrations of ULF MRI have been made. These include combined MRI and magnetoencephalography, imaging in the presence of metal, unique tissue contrast, and implementation in situations where a high magnetic field is simply impractical. These demonstrations have routinely used pulsed pre-polarization (at magnetic fields from ∼10 to 100mT) followed by read-out in a much weaker (1-100μT) magnetic fields using the ultra-sensitive Superconducting Quantum Interference Device (SQUID) sensor. Even with pre-polarization and SQUID detection, ULF MRI suffers from many challenges associated with lower magnetization (i.e. signal) and inherently long acquisition times compared to conventional >1T MRI. These are fundamental limitations imposed by the low measurement and gradient fields used. In this review article we discuss some of the techniques, potential applications, and inherent challenges of ULF MRI.

Good practice for conducting and reporting MEG research

Magnetoencephalographic (MEG) recordings are a rich source of information about the neural dynamics underlying cognitive processes in the brain, with excellent temporal and good spatial resolution. In recent years there have been considerable advances in MEG hardware developments and methods. Sophisticated analysis techniques are now routinely applied and continuously improved, leading to fascinating insights into the intricate dynamics of neural processes. However, the rapidly increasing level of complexity of the different steps in a MEG study make it difficult for novices, and sometimes even for experts, to stay aware of possible limitations and caveats. Furthermore, the complexity of MEG data acquisition and data analysis requires special attention when describing MEG studies in publications, in order to facilitate interpretation and reproduction of the results. This manuscript aims at making recommendations for a number of important data acquisition and data analysis steps and suggests details that should be specified in manuscripts reporting MEG studies. These recommendations will hopefully serve as guidelines that help to strengthen the position of the MEG research community within the field of neuroscience, and may foster discussion in order to further enhance the quality and impact of MEG research.

Combining Temporal and Spectral Information with Spatial Mapping to Identify Differences between Phonological and Semantic Networks: A Magnetoencephalographic Approach

Early, lesion-based models of language processing suggested that semantic and phonological processes are associated with distinct temporal and parietal regions respectively, with frontal areas more indirectly involved. Contemporary spatial brain mapping techniques have not supported such clear-cut segregation, with strong evidence of activation in left temporal areas by both processes and disputed evidence of involvement of frontal areas in both processes. We suggest that combining spatial information with temporal and spectral data may allow a closer scrutiny of the differential involvement of closely overlapping cortical areas in language processing. Using beamforming techniques to analyze magnetoencephalography data, we localized the neuronal substrates underlying primed responses to nouns requiring either phonological or semantic processing, and examined the associated measures of time and frequency in those areas where activation was common to both tasks. Power changes in the beta (14–30 Hz) and gamma (30–50 Hz) frequency bands were analyzed in pre-selected time windows of 350–550 and 500–700 ms In left temporal regions, both tasks elicited power changes in the same time window (350–550 ms), but with different spectral characteristics, low beta (14–20 Hz) for the phonological task and high beta (20–30 Hz) for the semantic task. In frontal areas (BA10), both tasks elicited power changes in the gamma band (30–50 Hz), but in different time windows, 500–700 ms for the phonological task and 350–550 ms for the semantic task. In the left inferior parietal area (BA40), both tasks elicited changes in the 20–30 Hz beta frequency band but in different time windows, 350–550 ms for the phonological task and 500–700 ms for the semantic task. Our findings suggest that, where spatial measures may indicate overlapping areas of involvement, additional beamforming techniques can demonstrate differential activation in time and frequency domains.

Neuromagnetic indicators of tinnitus and tinnitus masking in patients with and without hearing loss

Tinnitus is an auditory phenomenon characterised by the perception of a sound in the absence of an external auditory stimulus. Chronic subjective tinnitus is almost certainly maintained via central mechanisms, and this is consistent with observed measures of altered spontaneous brain activity. A number of putative central auditory mechanisms for tinnitus have been proposed. The influential thalamocortical dysrhythmia model suggests that tinnitus can be attributed to the disruption of coherent oscillatory activity between thalamus and cortex following hearing loss. However, the extent to which this disruption specifically contributes to tinnitus or is simply a consequence of the hearing loss is unclear because the necessary matched controls have not been tested. Here, we rigorously test several predictions made by this model in four groups of participants (tinnitus with hearing loss, tinnitus with clinically normal hearing, no tinnitus with hearing loss and no tinnitus with clinically normal hearing). Magnetoencephalography was used to measure oscillatory brain activity within different frequency bands in a ‘resting’ state and during presentation of a masking noise. Results revealed that low-frequency activity in the delta band (1–4 Hz) was significantly higher in the ‘tinnitus with hearing loss’ group compared to the ‘no tinnitus with normal hearing’ group. A planned comparison indicated that this effect was unlikely to be driven by the hearing loss alone, but could possibly be a consequence of tinnitus and hearing loss. A further interpretative linkage to tinnitus was given by the result that the delta activity tended to reduce when tinnitus was masked. High-frequency activity in the gamma band (25–80 Hz) was not correlated with tinnitus (or hearing loss). The findings partly support the thalamocortical dysrhythmia model and suggest that slow-wave (delta band) activity may be a more reliable correlate of tinnitus than high-frequency activity.

A general Bayesian treatment for MEG source reconstruction incorporating lead field uncertainty

There is uncertainty introduced when a cortical surface based model derived from an anatomical MRI is used to reconstruct neural activity with MEG data. This is a specific case of a problem with uncertainty in parameters on which M/EEG lead fields depend non-linearly. Here we present a general mathematical treatment of any such problem with a particular focus on co-registration. We use a Metropolis search followed by Bayesian Model Averaging over multiple sparse prior source inversions with different headlocation/orientation parameters. Based on MEG data alone we can locate the cortex to within 4mm at empirically realistic signal to noise ratios. We also show that this process gives improved posterior distributions on the estimated current distributions, and can be extended to make inference on the locations of local maxima by providing confidence intervals for each source.

A novel binaural pitch elicited by phase-modulated noise: MEG and psychophysical observations

Binaural pitches are auditory percepts that emerge from combined inputs to the ears but that cannot be heard if the stimulus is presented to either ear alone. Here, we describe a binaural pitch that is not easily accommodated within current models of binaural processing. Convergent magnetoencephalography (MEG) and psychophysical measurements were used to characterize the pitch, heard when band-limited noise had a rapidly changing interaural phase difference. Several interesting features emerged: First, the pitch was perceptually lateralized, in agreement with the lateralization of the evoked changes in MEG spectral power, and its salience depended on dichotic binaural presentation. Second, the frequency of the pure tone that matched the binaural pitch lay within a lower spectral sideband of the phase-modulated noise and followed the frequency of that sideband when the modulation frequency or center frequency and bandwidth of the noise changed. Thus, the binaural pitch depended on the processing of binaural information in that lower sideband.