MEG Source Localization via Deep Learning

We present a deep learning solution to the problem of localization of magnetoencephalography (MEG) brain signals. The proposed deep model architectures are tuned for single and multiple time point MEG data, and can estimate varying numbers of dipole sources. Results from simulated MEG data on the cortical surface of a real human subject demonstrated improvements against the popular RAP-MUSIC localization algorithm in specific scenarios with varying SNR levels, inter-source correlation values, and number of sources. Importantly, the deep learning models had robust performance to forward model errors and a significant reduction in computation time, to a fraction of 1 ms, paving the way to real-time MEG source localization

Measuring electrophysiological connectivity by power envelope correlation: a technical review on MEG methods

The human brain can be divided into multiple areas, each responsible for different aspects of behaviour. Healthy brain function relies upon efficient connectivity between these areas and, in recent years, neuroimaging has been revolutionised by an ability to estimate this connectivity. In this paper we discuss measurement of network connectivity using magnetoencephalography (MEG), a technique capable of imaging electrophysiological brain activity with good (~5mm) spatial resolution and excellent (~1ms) temporal resolution. The rich information content of MEG facilitates many disparate measures of connectivity between spatially separate regions and in this paper we discuss a single metric known as power envelope correlation. We review in detail the methodology required to measure power envelope correlation including i) projection of MEG data into source space, ii) removing confounds introduced by the MEG inverse problem and iii) estimation of connectivity itself. In this way, we aim to provide researchers with a description of the key steps required to assess envelope based functional networks, which are thought to represent an intrinsic mode of coupling in the human brain. We highlight the principal findings of the techniques discussed, and furthermore, we show evidence that this method can probe how the brain forms and dissolves multiple transient networks on a rapid timescale in order to support current processing demand. Overall, power envelope correlation offers a unique and verifiable means to gain novel insights into network coordination and is proving to be of significant value in elucidating the neural dynamics of the human connectome in health and disease

Testing covariance models for MEG source reconstruction of hippocampal activity

Beamforming is one of the most commonly used source reconstruction methods for magneto- and electroencephalography (M/EEG). One underlying assumption, however, is that distant sources are uncorrelated and here we tested whether this is an appropriate model for the human hippocampal data. We revised the Empirical Bayesian Beamfomer (EBB) to accommodate specific a-priori correlated source models. We showed in simulation that we could use model evidence (as approximated by Free Energy) to distinguish between different correlated and uncorrelated source scenarios. Using group MEG data in which the participants performed a hippocampal-dependent task, we explored the possibility that the hippocampus or the cortex or both were correlated in their activity across hemispheres. We found that incorporating a correlated hippocampal source model significantly improved model evidence. Our findings help to explain why, up until now, the majority of MEG-reported hippocampal activity (typically making use of beamformers) has been estimated as unilateral

The impact of MEG source reconstruction method on source-space connectivity estimation: A comparison between minimum-norm solution and beamforming.

Despite numerous important contributions, the investigation of brain connectivity with magnetoencephalography (MEG) still faces multiple challenges. One critical aspect of source-level connectivity, largely overlooked in the literature, is the putative effect of the choice of the inverse method on the subsequent cortico-cortical coupling analysis. We set out to investigate the impact of three inverse methods on source coherence detection using simulated MEG data. To this end, thousands of randomly located pairs of sources were created. Several parameters were manipulated, including inter- and intra-source correlation strength, source size and spatial configuration. The simulated pairs of sources were then used to generate sensor-level MEG measurements at varying signal-to-noise ratios (SNR). Next, the source level power and coherence maps were calculated using three methods (a) L2-Minimum-Norm Estimate (MNE), (b) Linearly Constrained Minimum Variance (LCMV) beamforming, and (c) Dynamic Imaging of Coherent Sources (DICS) beamforming. The performances of the methods were evaluated using Receiver Operating Characteristic (ROC) curves. The results indicate that beamformers perform better than MNE for coherence reconstructions if the interacting cortical sources consist of point-like sources. On the other hand, MNE provides better connectivity estimation than beamformers, if the interacting sources are simulated as extended cortical patches, where each patch consists of dipoles with identical time series (high intra-patch coherence). However, the performance of the beamformers for interacting patches improves substantially if each patch of active cortex is simulated with only partly coherent time series (partial intra-patch coherence). These results demonstrate that the choice of the inverse method impacts the results of MEG source-space coherence analysis, and that the optimal choice of the inverse solution depends on the spatial and synchronization profile of the interacting cortical sources. The insights revealed here can guide method selection and help improve data interpretation regarding MEG connectivity estimation

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

On the potential of a new generation of magnetometers for MEG: a beamformer simulation study

Magnetoencephalography (MEG) is a sophisticated tool which yields rich information on the spatial, spectral and temporal signatures of human brain function. Despite unique potential, MEG is limited by a low signal-to-noise ratio (SNR) which is caused by both the inherently small magnetic fields generated by the brain, and the scalp-to-sensor distance. The latter is limited in current systems due to a requirement for pickup coils to be cryogenically cooled. Recent work suggests that optically-pumped magnetometers (OPMs) might be a viable alternative to superconducting detectors for MEG measurement. They have the advantage that sensors can be brought to within ~4 mm of the scalp, thus offering increased sensitivity. Here, using simulations, we quantify the advantages of hypothetical OPM systems in terms of sensitivity, reconstruction accuracy and spatial resolution. Our results show that a multi-channel whole-head OPM system offers (on average) a fivefold improvement in sensitivity for an adult brain, as well as clear improvements in reconstruction accuracy and spatial resolution. However, we also show that such improvements depend critically on accurate forward models; indeed, the reconstruction accuracy of our simulated OPM system only outperformed that of a simulated superconducting system in cases where forward field error was less than 5%. Overall, our results imply that the realisation of a viable whole-head multi-channel OPM system could generate a step change in the utility of MEG as a means to assess brain electrophysiological activity in health and disease. However in practice, this will require both improved hardware and modelling algorithms

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

Measuring temporal, spectral and spatial changes in electrophysiological brain network connectivity

The topic of functional connectivity in neuroimaging is expanding rapidly and many studies now focus on coupling between spatially separate brain regions. These studies show that a relatively small number of large scale networks exist within the brain, and that healthy function of these networks is disrupted in many clinical populations. To date, the vast majority of studies probing connectivity employ techniques that compute time averaged correlation over several minutes, and between specific pre-defined brain locations. However, increasing evidence suggests that functional connectivity is non-stationary in time. Further, electrophysiological measurements show that connectivity is dependent on the frequency band of neural oscillations. It is also conceivable that networks exhibit a degree of spatial inhomogeneity, i.e. the large scale networks that we observe may result from the time average of multiple transiently synchronised sub-networks, each with their own spatial signature. This means that the next generation of neuroimaging tools to compute functional connectivity must account for spatial inhomogeneity, spectral non-uniformity and temporal non-stationarity. Here, we present a means to achieve this via application of windowed canonical correlation analysis (CCA) to source space projected MEG data. We describe the generation of time–frequency connectivity plots, showing the temporal and spectral distribution of coupling between brain regions. Moreover, CCA over voxels provides a means to assess spatial non-uniformity within short time–frequency windows. The feasibility of this technique is demonstrated in simulation and in a resting state MEG experiment where we elucidate multiple distinct spatio-temporal-spectral modes of covariation between the left and right sensorimotor areas