High Precision Anatomy for MEG

Magnetoencephalography (MEG) is a non-invasive brain imaging method with high temporal resolution but relatively poor spatial resolution as compared to some other non-invasive techniques. This thesis examines how the spatial resolution of MEG can be improved using new recording paradigms in which the subject’s head position is fixed and known in advance. This is accomplished by using subject-specific head casts made using a combination of structural MRI and 3D printing technology. The resulting high-precision spatial models allow one to make inference at spatial scales of the order of cortical laminae. This thesis outlines the design of the head casts and examines the potential theoretical and empirical advances they offer. Specifically I outline simulation and then empirical investigations showing it is possible to non-invasively distinguish between electrophysiological signals in different layers of the cortex

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 1mm. Estimates of relative co-registration error were <1.5mm between sessions. We corroborated these scalp based estimates by looking at the MEG data recorded over a 6month 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 5mm. 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

Magnetoencephalography as a tool in psychiatric research: current status and perspective

The application of neuroimaging to provide mechanistic insights into circuit dysfunctions in major psychiatric conditions and the development of biomarkers are core challenges in current psychiatric research. In this review, we propose that recent technological and analytic advances in Magnetoencephalography (MEG), a technique which allows the measurement of neuronal events directly and non-invasively with millisecond resolution, provides novel opportunities to address these fundamental questions. Because of its potential in delineating normal and abnormal brain dynamics, we propose that MEG provides a crucial tool to advance our understanding of pathophysiological mechanisms of major neuropsychiatric conditions, such as Schizophrenia, Autism Spectrum Disorders, and the dementias. In our paper, we summarize the mechanisms underlying the generation of MEG signals and the tools available to reconstruct generators and underlying networks using advanced source-reconstruction techniques. We then survey recent studies that have utilized MEG to examine aberrant rhythmic activity in neuropsychiatric disorders. This is followed by links with preclinical research, which have highlighted possible neurobiological mechanisms, such as disturbances in excitation/inhibition parameters, which could account for measured changes in neural oscillations. In the final section of the paper, challenges as well as novel methodological developments are discussed which could pave the way for a widespread application of MEG in translational research with the aim of developing biomarkers for early detection and diagnosis

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

With a hybrid 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 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 was 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.Comment: 11 pages, 8 figures. This work is part of the BREAKBEN project and has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 68686

Neuromagnetic effects of pico-Tesla stimulation

We used a double-blind experimental design to look for an effect of pico-Tesla magnetic stimulation in healthy subjects. Pico-Tesla stimulation is thought to increase the dominant frequency of 2–7 Hz oscillations in the human brain. We used magnetoencephalography to measure resting state brain activity. Each subject had two separate recording sessions consisting of three runs in between which they were given real or sham pT stimulation. We then tried to predict the real and sham stimulation sessions based on changes in the mean peak frequency in the 2–7 Hz band. Our predictions for these individual runs were 8 out of 14 at chance level (p = 0.39). After unblinding, we found no significant effect (p = 0.11) of an increase in the frequency range (2–7 Hz) across the subject group. Finally, we performed a Bayesian model comparison between the effect size predicted from previous clinical studies and a null model. Even though this study had a sensitivity advantage of at least one order of magnitude over previous work, we found the null model to be significantly (2000 times) more likely

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

Medial prefrontal-medial temporal theta phase coupling in dynamic spatial imagery

Hippocampal–medial prefrontal interactions are thought to play a crucial role in mental simulation. Notably, the frontal midline/medial pFC (mPFC) theta rhythm in humans has been linked to introspective thought and working memory. In parallel, theta rhythms have been proposed to coordinate processing in the medial temporal cortex, retrosplenial cortex (RSc), and parietal cortex during the movement of viewpoint in imagery, extending their association with physical movement in rodent models. Here, we used noninvasive whole-head MEG to investigate theta oscillatory power and phase-locking during the 18-sec postencoding delay period of a spatial working memory task, in which participants imagined previously learned object sequences either on a blank background (object maintenance), from a first-person viewpoint in a scene (static imagery), or moving along a path past the objects (dynamic imagery). We found increases in 4- to 7-Hz theta power in mPFC when comparing the delay period with a preencoding baseline. We then examined whether the mPFC theta rhythm was phase-coupled with ongoing theta oscillations elsewhere in the brain. The same mPFC region showed significantly higher theta phase coupling with the posterior medial temporal lobe/RSc for dynamic imagery versus either object maintenance or static imagery. mPFC theta phase coupling was not observed with any other brain region. These results implicate oscillatory coupling between mPFC and medial temporal lobe/RSc theta rhythms in the dynamic mental exploration of imagined scenes

Data?driven model optimization for optically pumped magnetometer sensor arrays

© 2019 The Authors. Human Brain Mapping published by Wiley Periodicals, Inc. Optically pumped magnetometers (OPMs) have reached sensitivity levels that make them viable portable alternatives to traditional superconducting technology for magnetoencephalography (MEG). OPMs do not require cryogenic cooling and can therefore be placed directly on the scalp surface. Unlike cryogenic systems, based on a well-characterised fixed arrays essentially linear in applied flux, OPM devices, based on different physical principles, present new modelling challenges. Here, we outline an empirical Bayesian framework that can be used to compare between and optimise sensor arrays. We perturb the sensor geometry (via simulation) and with analytic model comparison methods estimate the true sensor geometry. The width of these perturbation curves allows us to compare different MEG systems. We test this technique using simulated and real data from SQUID and OPM recordings using head-casts and scanner-casts. Finally, we show that given knowledge of underlying brain anatomy, it is possible to estimate the true sensor geometry from the OPM data themselves using a model comparison framework. This implies that the requirement for accurate knowledge of the sensor positions and orientations a priori may be relaxed. As this procedure uses the cortical manifold as spatial support there is no co-registration procedure or reliance on scalp landmarks

Optimising experimental design for MEG resting state functional connectivity measurement

The study of functional connectivity using magnetoencephalography (MEG) is an expanding area of neuroimaging, and adds an extra dimension to the more common assessments made using fMRI. The importance of such metrics is growing, with recent demonstrations of their utility in clinical research, however previous reports suggest that whilst group level resting state connectivity is robust, single session recordings lack repeatability. Such robustness is critical if MEG measures in individual subjects are to prove clinically valuable. In the present paper, we test how practical aspects of experimental design affect the intra-subject repeatability of MEG findings; specifically we assess the effect of co-registration method and data recording duration. We show that the use of a foam head-cast, which is known to improve co-registration accuracy, increased significantly the between session repeatability of both beamformer reconstruction and connectivity estimation. We also show that recording duration is a critical parameter, with large improvements in repeatability apparent when using ten minute, compared to five minute recordings. Further analyses suggest that the origin of this latter effect is not underpinned by technical aspects of source reconstruction, but rather by a genuine effect of brain state; short recordings are simply inefficient at capturing the canonical MEG network in a single subject. Our results provide important insights on experimental design and will prove valuable for future MEG connectivity studies

Algorithmic procedures for Bayesian MEG/EEG source reconstruction in SPM

AbstractThe MEG/EEG inverse problem is ill-posed, giving different source reconstructions depending on the initial assumption sets. Parametric Empirical Bayes allows one to implement most popular MEG/EEG inversion schemes (Minimum Norm, LORETA, etc.) within the same generic Bayesian framework. It also provides a cost-function in terms of the variational Free energy—an approximation to the marginal likelihood or evidence of the solution. In this manuscript, we revisit the algorithm for MEG/EEG source reconstruction with a view to providing a didactic and practical guide. The aim is to promote and help standardise the development and consolidation of other schemes within the same framework. We describe the implementation in the Statistical Parametric Mapping (SPM) software package, carefully explaining each of its stages with the help of a simple simulated data example. We focus on the Multiple Sparse Priors (MSP) model, which we compare with the well-known Minimum Norm and LORETA models, using the negative variational Free energy for model comparison. The manuscript is accompanied by Matlab scripts to allow the reader to test and explore the underlying algorithm