Effect of head shape variations among individuals on the EEG/MEG forward and inverse problems

We study the effect of the head shape variations on the EEG/magnetoencephalography (MEG) forward and inverse problems. We build a random head model such that each sample represents the head shape of a different individual and solve the forward problem assuming this random head model, using a polynomial chaos expansion. The random solution of the forward problem is then used to quantify the effect of the geometry when the inverse problem is solved with a standard head model. The results derived with this approach are valid for a continuous family of head models, rather than just for a set of cases. The random model consists of three random surfaces that define layers of different electric conductivity, and we built an example based on a set of 30 deterministic models from adults. Our results show that for a dipolar source model, the effect of the head shape variations on the EEG/MEG inverse problem due to the random head model is slightly larger than the effect of the electronic noise present in the sensors. The variations in the EEG inverse problem solutions are due to the variations in the shape of the volume conductor, while the variations in the MEG inverse problem solutions, larger than the EEG ones, are caused mainly by the variations of the absolute position of the sources in a coordinate system based on anatomical landmarks, in which the magnetometers have a fixed position.Instituto de Investigaciones en Electrónica, Control y Procesamiento de Señale

Uncertainty through polynomial chaos in the EEG problem

A sensitivity and correlation analysis of EEG sensors influenced by uncertain conductivity is conducted. We assume a three layer spherical head model with different and random layer conductivities. This randomness is modeled by Polynomial Chaos (PC). On average, we observe the least influenced electrodes along the great longitudinal fissure. Also, sensors located closer to a dipole source, are of greater influence to a change in conductivity -- this is in agreement with previous research. The highly influenced sensors were on average located temporal. This was also the case in the correlation analysis, which was made possible by our approach with PC. Sensors in the temporal parts of the brain are highly correlated. Whereas the sensors in the occipital and lower frontal region, though they are close together, are not so highly correlated as in the temporal regions

Techniques for imaging small impedance changes in the human head due to neuronal depolarisation

A new imaging modality is being developed, which may be capable of imaging small impedance changes in the human head due to neuronal depolarization. One way to do this would be by imaging the impedance changes associated with ion channels opening in neuronal membranes in the brain during activity. The results of previous modelling and experimental studies indicated that impedance changes between 0.6%and 1.7% locally in brain grey matter when recorded at DC. This reduces by a further of 10% if measured at the surface of the head, due to distance and the effect of the resistive skull. In principle, this could be measured using Electrical Impedance Tomography (ElT) but it is close to its threshold of detectability. With the inherent limitation in the use of electrodes, this work proposed two new schemes. The first is a magnetic measurement scheme based on recording the magnetic field with Superconducting Quantum Interference Devices (SQUIDs), used in Magnetoencephalography (MEG) as a result of a non-invasive injection of current into the head. This scheme assumes that the skull does not attenuate the magnetic field. The second scheme takes into consideration that the human skull is irregular in shape, with less and varying conductivity as compared to other head tissues. Therefore, a key issue is to know through which electrodes current can be injected in order to obtain high percentage changes in surface potential when there is local conductivity change in the head. This model will enable the prediction of the current density distribution at specific regions in the brain with respect to the varying skull and local conductivities. In the magnetic study, the head was modelled as concentric spheres, and realistic head shapes to mimic the scalp, skull, Cerebrospinal Auid (CSF) and brain using the Finite Element Method (FEM). An impedance change of 1 % in a 2cm-radius spherical volume depicting the physiological change in the brain was modelled as the region of depolarisation. The magnetic field, 1 cm away from the scalp, was estimated on injecting a constant current of 100 µA into the head from diametrically opposed electrodes. However, in the second scheme, only the realistic FEM of the head was used, which included a specific region of interest; the primary visual cortex (V1). The simulated physiological change was the variation in conductivity of V1 when neurons were assumed to be firing during a visual evoked response. A near DC current of 100 µA was driven through possible pairs of 31 electrodes using ElT techniques. For a fixed skull conductivity, the resulting surface potentials were calculated when the whole head remained unperturbed, or when the conductivity of V1 changed by 0.6%, 1 %, and 1.6%. The results of the magnetic measurement predicted that standing magnetic field was about 10pT and the field changed by about 3fT (0.03%) on depolarization. For the second scheme, the greatest mean current density through V1 was 0.020 ± 0.005 µAmm-2, and occurred with injection through two electrodes positioned near the occipital cortex. The corresponding maximum change in potential from baseline was 0.02%. Saline tank experiments confirmed the accuracy of the estimated standing potentials. As the noise density in a typical MEG system in the frequency band is about 7fT/√Hz, it places the change at the limit of detectability due to low signal to noise ratio. This is therefore similar to electrical recording, as in conventional ElT systems, but there may be advantages to MEG in that the magnetic field direcdy traverses the skull and instrumentation errors from the electrode-skin interface will be obviated. This has enabled the estimation of electrode positions most likely to permit recording of changes in human experiments and suggests that the changes, although tiny, may just be discernible from noise

Reconstruction of electric fields and source distributions in EEG brain imaging

In this thesis, three different approaches are developed for the estimation of focal brain activity using EEG measurements. The proposed approaches have been tested and found feasible using simulated data. First, we develop a robust solver for the recovery of focal dipole sources. The solver uses a weighted dipole strength penalty term (also called weighted L1,2 norm) as prior information in order to ensure that the sources are sparse and focal, and that both the source orientation and depth bias are reduced. The solver is based on the truncated Newton interior point method combined with a logarithmic barrier method for the approximation of the penalty term. In addition, we use a Bayesian framework to derive the depth weights in the prior that are used to reduce the tendency of the solver to favor superficial sources. In the second approach, vector field tomography (VFT) is used for the estimation of underlying electric fields inside the brain from external EEG measurements. The electric field is reconstructed using a set of line integrals. This is the first time that VFT has been used for the recovery of fields when the dipole source lies inside the domain of reconstruction. The benefit of this approach is that we do not need a mathematical model for the sources. The test cases indicated that the approach can accurately localize the source activity. In the last part of the thesis, we show that, by using the Bayesian approximation error approach (AEA), precise knowledge of the tissue conductivities and head geometry are not always needed. We deliberately use a coarse head model and we take the typical variations in the head geometry and tissue conductivities into account statistically in the inverse model. We demonstrate that the AEA results are comparable to those obtained with an accurate head model.Open Acces

Estimating neural activity from visual areas using functionally defined EEG templates

This work was supported by BBSRC grant BB/N018516/1 and Wellcome Trust ISSF 204821/Z/16/Z.Electroencephalography (EEG) is a common and inexpensive method to record neural activity in humans. However, it lacks spatial resolution making it difficult to determine which areas of the brain are responsible for the observed EEG response. Here we present a new easy-to-use method that relies on EEG topographical templates. Using MRI and fMRI scans of 50 participants, we simulated how the activity in each visual area appears on the scalp and averaged this signal to produce functionally defined EEG templates. Once created, these templates can be used to estimate how much each visual area contributes to the observed EEG activity. We tested this method on extensive simulations and on real data. The proposed procedure is as good as bespoke individual source localization methods, robust to a wide range of factors, and has several strengths. First, because it does not rely on individual brain scans, it is inexpensive and can be used on any EEG data set, past or present. Second, the results are readily interpretable in terms of functional brain regions and can be compared across neuroimaging techniques. Finally, this method is easy to understand, simple to use and expandable to other brain sources.Publisher PDFPeer reviewe

The functional neuroanatomy of auditory sensory gating and its behavioural implications

Auditory sensory gating (ASG) is the ability in individuals to suppress incoming irrelevant sensory input, indexed by evoked response to paired auditory stimuli. ASG is impaired in psychopathology such as schizophrenia, in which it has been proposed as putative endophenotype. This study aims to characterise electrophysiological properties of the phenomenon using MEG in time and frequency domains as well as to localise putative networks involved in the process at both sensor and source level. We also investigated the relationship between ASG measures and personality profiles in healthy participants in the light of its candidate endophenotype role in psychiatric disorders. Auditory evoked magnetic fields were recorded in twenty seven healthy participants by P50 ‘paired-click’ paradigm presented in pairs (conditioning stimulus S1- testing stimulus S2) at 80dB, separated by 250msec with inter trial interval of 7-10 seconds. Gating ratio in healthy adults ranged from 0.5 to 0.8 suggesting dimensional nature of P50 ASG. The brain regions active during this process were bilateral superior temporal gyrus (STG) and bilateral inferior frontal gyrus (IFG); activation was significantly stronger in IFG during S2 as compared to S1 (at p<0.05). Measures of effective connectivity between these regions using DCM modelling revealed the role of frontal cortex in modulating ASG as suggested by intracranial studies, indicating major role of inhibitory interneuron connections. Findings from this study identified a unique event-related oscillatory pattern for P50 ASG with alpha (STG)-beta (IFG) desynchronization and increase in cortical oscillatory gamma power (IFG) during S2 condition as compared to S1. These findings show that the main generator for P50 response is within temporal lobe and that inhibitory interneurons and gamma oscillations in the frontal cortex contributes substantially towards sensory gating. Our findings also show that ASG is a predictor of personality profiles (introvert vs extrovert dimension)