High-resolution EEG source localization in personalized segmentation-free head model with multi-dipole fitting

Objective. Electroencephalograms (EEGs) are often used to monitor brain activity. Several source localization methods have been proposed to estimate the location of brain activity corresponding to EEG readings. However, only a few studies evaluated source localization accuracy from measured EEG using personalized head models in a millimeter resolution. In this study, based on a volume conductor analysis of a high-resolution personalized human head model constructed from magnetic resonance images, a finite difference method was used to solve the forward problem and to reconstruct the field distribution.Approach. We used a personalized segmentation-free head model developed using machine learning techniques, in which the abrupt change of electrical conductivity occurred at the tissue interface is suppressed. Using this model, a smooth field distribution was obtained to address the forward problem. Next, multi-dipole fitting was conducted using EEG measurements for each subject (N= 10 male subjects, age: 22.5 ± 0.5), and the source location and electric field distribution were estimated.Main results.For measured somatosensory evoked potential for electrostimulation to the wrist, a multi-dipole model with lead field matrix computed with the volume conductor model was found to be superior than a single dipole model when using personalized segmentation-free models (6/10). The correlation coefficient between measured and estimated scalp potentials was 0.89 for segmentation-free head models and 0.71 for conventional segmented models. The proposed method is straightforward model development and comparable localization difference of the maximum electric field from the target wrist reported using fMR (i.e. 16.4 ± 5.2 mm) in previous study. For comparison, DUNEuro based on sLORETA was (EEG: 17.0 ± 4.0 mm). In addition, somatosensory evoked magnetic fields obtained by Magnetoencephalography was 25.3 ± 8.5 mm using three-layer sphere and sLORETA.Significance. For measured EEG signals, our procedures using personalized head models demonstrated that effective localization of the somatosensory cortex, which is located in a non-shallower cortex region. This method may be potentially applied for imaging brain activity located in other non-shallow regions.

Assessing the effects of head modelling errors and measurement noise on EEG source localization accuracy in preterm newborns: A single-subject study

The accuracy of electroencephalogram (EEG) source localization is compromised because of head modelling errors. In this study, we investigated the effect of inaccuracy in the conductivity of head tissues and head model structural deficiencies on the accuracy of EEG source analysis in premature neonates. A series of EEG forward and inverse simulations was performed by introducing structural deficiencies into the reference head models to generate test models, which were then used to investigate head modelling errors caused by cerebrospinal fluid (CSF) exclusion, lack of grey matter (GM)-white matter (WM) distinction, fontanel exclusion and inaccuracy in skull conductivity. The modelling errors were computed between forward and inverse solutions obtained using the reference and test models generated for each deficiency. Our results showed that the exclusion of CSF from the head model had a strong widespread effect on the accuracy of the EEG source localization with position errors lower than 4.17 mm. The GM and WM distinction also caused strong localization errors (up to 3.5 mm). The exclusion of fontanels from the head model also strongly affected the accuracy of the EEG source localization for sources located beneath the fontanels with a maximum localization error of 4.37 mm. Similarly, inaccuracies in the skull conductivity caused errors in EEG forward and inverse modelling in sources beneath cranial bones. Our results indicate that the accuracy of EEG source imaging in premature neonates can be largely improved by using head models, which include not only the brain, skull and scalp but also the CSF, GM, WM and fontanels.

Evaluating the Influence of Anatomical Accuracy and Electrode Positions on EEG Forward Solutions

Generating realistic volume conductor models for forward calculations in electroencephalography (EEG) is not trivial and several factors contribute to the accuracy of such models, two of which are its anatomical accuracy and the accuracy with which electrode positions are known. Here, we investigate effects of anatomical accuracy by comparing forward solutions from SimNIBS, a tool which allows state-of-the-art anatomical modeling, with well-established pipelines in MNE-Python and FieldTrip. We also compare different ways of specifying electrode locations when digitized positions are not available such as transformation of measured positions from standard space and transformation of a manufacturer layout. Substantial effects of anatomical accuracy were seen throughout the entire brain both in terms of field topography and magnitude with SimNIBS generally being more accurate than the pipelines in MNE-Python and FieldTrip. Topographic and magnitude effects were particularly pronounced for MNE-Python which uses a three-layer boundary element method (BEM) model. We attribute these mainly to the coarse representation of the anatomy used in this model, in particular differences in skull and cerebrospinal fluid (CSF). Effects of electrode specification method were evident in occipital and posterior areas when using a transformed manufacturer layout whereas transforming measured positions from standard space generally resulted in smaller errors. We suggest modeling the anatomy of the volume conductor as accurately possible and we hope to facilitate this by making it easy to export simulations from SimNIBS to MNE-Python and FieldTrip for further analysis. Likewise, if digitized electrode positions are not available, a set of measured positions on a standard head template may be preferable to those specified by the manufacturer.

Bioelectromagnetism in Human Brain Research: New Applications, New Questions

Bioelectromagnetism has contributed some of the most commonly used techniques to human neuroscience such as magnetoencephalography (MEG), electroencephalography (EEG), transcranial magnetic stimulation (TMS), and transcranial electric stimulation (TES). The considerable differences in their technical design and practical use give rise to the impression that these are quite different techniques altogether. Here, we review, discuss and illustrate the fundamental principle of Helmholtz reciprocity that provides a common ground for all four techniques. We show that, more than 150 years after its discovery by Helmholtz in 1853, reciprocity is important to appreciate the strengths and limitations of these four classical tools in neuroscience. We build this case by explaining the concept of Helmholtz reciprocity, presenting a methodological account of this principle for all four methods and, finally, by illustrating its application in practical clinical studies.

An investigation of the effect of brain atrophy on brain injury in multiple sclerosis

Multiple sclerosis (MS) is a disease of the central nervous system (CNS) that affects the brain and spinal cord. It is estimated that the average prevalence of MS is 35.9 cases per 100,000 and a total of 2.8 million people worldwide have MS. Brain atrophy is usually seen in the early stages of MS, and its progress is faster than healthy people. The present study was a numerical study that uses the Fluid-structure interaction (FSI) model to investigate the effect of brain atrophy on brain injury in MS. Firstly, a healthy model was constructed from MRI images and validated by experimental data. Then three models with different degrees of brain atrophy, which showed the rate of brain atrophy in different years in MS patients, were developed to model the brain atrophy in MS. The models were subjected to two different types of impact conditions. Type I, which only produced a translational motion and the HIC value of 744, was applied to each model. Type II produced both translational and rotational motion. In this type of impact, the experimental kinematics, with peaks of 450 g for the translational acceleration and 26.2 krad/s2 for the rotational acceleration, were applied to the nodes that located in the center of gravity of the head models and the results were extracted from each one. According to the results of impact type I, the pressure of the frontal lobe of the brain is 149,647 Pa in the health model and 137,690 Pa in the model with severe atrophy.

Patient-specific solution of the electrocorticography forward problem in deforming brain

Invasive intracranial electroencephalography (iEEG), or electrocorticography (ECoG), measures electric potential directly on the surface of the brain and can be used to inform treatment planning for epilepsy surgery. Combined with numerical modeling it can further improve accuracy of epilepsy surgery planning. Accurate solution of the iEEG forward problem, which is a crucial prerequisite for solving the iEEG inverse problemin epilepsy seizure onset zone localization, requires accurate representation of the patient's brain geometry and tissue electrical conductivity after implantation of electrodes. However, implantation of subdural grid electrodes causes the brain to deform, which invalidates preoperatively acquired image data. Moreover, postoperative magnetic resonance imaging (MRI) is incompatible with implanted electrodes and computed tomography (CT) has insufficient range of soft tissue contrast, which precludes both MRI and CT from being used to obtain the deformed postoperative geometry. In this paper, we present a biomechanics-based image warping procedure using preoperative MRI for tissue classification and postoperative CT for locating implanted electrodes to perform non-rigid registration of the preoperative image data to the postoperative configuration. We solve the iEEG forward problem on the predicted postoperative geometry using the finite element method (FEM) which accounts for patient-specific inhomogeneity and anisotropy of tissue conductivity. Results for the simulation of a current source in the brain show large differences in electric potential predicted by the models based on the original images and the deformed images corresponding to the brain geometry deformed by placement of invasive electrodes. Computation of the lead field matrix (useful for solution of the iEEG inverse problem) also showed significant differences between the different models. The results suggest that rapid and accurate solution of the forward problem in a deformed brain for a given patient is achievable.

Cortical Source Analysis of Event-Related Potentials: A Developmental Approach

Cortical source analysis of electroencephalographic (EEG) signals has become an important tool in the analysis of brain activity. The aim of source analysis is to reconstruct the cortical generators (sources) of the EEG signal recorded on the scalp. The quality of the source reconstruction relies on the accuracy of the forward problem, and consequently the inverse problem. An accurate forward solution is obtained when an appropriate imaging modality (i.e., structural magnetic resonance imaging - MRI) is used to describe the head geometry, precise electrode locations are identified with 3D maps of the sensor positions on the scalp, and realistic conductivity values are determined for each tissue type of the head model. Together these parameters contribute to the definition of realistic head models. Here, we describe the steps necessary to reconstruct the cortical generators of the EEG signal recorded on the scalp. We provide an example of source reconstruction of event-related potentials (ERPs) during a face-processing task performed by a 6-month-old infant. We discuss the adjustments necessary to perform source analysis with measures different from the ERPs. The proposed pipeline can be applied to the investigation of different cognitive tasks in both younger and older participants.

Impact of skull sutures, spongiform bone distribution, and aging skull conductivities on the EEG forward and inverse problems

Source imaging is a principal objective for electroencephalography (EEG), the solutions of which require forward problem (FP) computations characterising the electric potential distribution on the scalp due to known sources. Additionally, the EEG-FP is dependent upon realistic, anatomically correct volume conductors and accurate tissue conductivities, where the skull is particularly important. Skull conductivity, however, deviates according to bone composition and the presence of adult sutures. The presented study therefore analyses the effect the presence of adult sutures and differing bone composition have on the EEG-FP and inverse problem (IP) solutions. Utilising a well-established head atlas, detailed head models were generated including compact and spongiform bone and adult sutures. The true skull conductivity was considered as inhomogeneous according to spongiform bone proportion and sutures. The EEG-FP and EEG-IP were solved and compared to results employing homogeneous skull models, with varying conductivities and omitting sutures, as well as using a hypothesised aging skull conductivity model. Significant localised FP errors, with relative error up to 85%, were revealed, particularly evident along suture lines and directly related to the proportion of spongiform bone. This remained evident at various ages. Similar EEG-IP inaccuracies were found, with the largest (maximum 4.14 cm) across suture lines. It is concluded that modelling the skull as an inhomogeneous layer that varies according to spongiform bone proportion and includes differing suture conductivity is imperative for accurate EEG-FP and source localisation calculations. Their omission can result in significant errors, relevant for EEG research and clinical diagnosis.

Investigating the neural signature of microsleeps using EEG

A microsleep (MS) is a complete lapse of responsiveness due to an episode of brief sleep (≲ 15 s) with eyes partially or completely closed. MSs are highly correlated with the risk of car accidents, severe injuries, and death. To investigate EEG changes during MSs, we used a 2D continuous visuomotor tracking (CVT) task and eye-video to identify MSs in 20 subjects performing the 50-min task. Following pre-processing, FFT spectral analysis was used to calculate the activity in the EEG delta, theta, alpha, beta, and gamma bands, followed by eLORETA for source reconstruction. A group statistical analysis was performed to compare the change in activity over EEG bands of an MS to its baseline. After correction for multiple comparisons, we found maximum increases in delta, theta, and alpha activities over the frontal lobe, and beta over the parietal and occipital lobes. There were no significant changes in the gamma band, and no significant decreases in any band. Our results are in agreement with previous studies which reported increased alpha activity in MSs. However, this is the first study to have reported increased beta activity during MSs, which, due to the usual association of beta activity with wakefulness, was unexpected.

Age- and gender-specific characteristics of the resting-state brain activity: a magnetoencephalography study

Aging and gender influence regional brain activities. Although these biases should be considered during the clinical examinations using magnetoencephalography, they have yet to be standardized. In the present study, resting-state magnetoencephalography data were recorded from 54 healthy females and 48 males aged 22 to 75 years, who were controlled for cognitive performance. The regional oscillatory power was estimated for each frequency band (delta, theta, alpha, beta, low-gamma, and high-gamma) using the sLORETA-like algorithm and the biases of age and gender were evaluated, respectively. The results showed that faster oscillatory powers increased with age in the rostral regions and decreased in the caudal regions, while few slower oscillatory powers changed with age. Gender differences in oscillatory powers were found in a broad frequency range, mostly in the caudal brain regions. The present study characterized the effects of healthy aging and gender asymmetricity on the regional resting-state brain activity, with the aim to facilitate the accurate and efficient use of magnetoencephalography in clinical practice.

Effect of structural complexities in head modeling on the accuracy of EEG source localization in neonates

OBJECTIVE

Neonatal electroencephalography (EEG) source localization is highly prone to errors due to head modeling deficiencies. In this study, we investigated the effect of head model complexities on the accuracy of EEG source localization in full term neonates using a realistic volume conductor head model.

APPROACH

We performed numerical simulations to investigate source localization errors caused by cerebrospinal fluid (CSF) and fontanel exclusion and gray matter (GM)/white matter (WM) distinction using the finite element method.

MAIN RESULTS

Our results showed that the exclusion of CSF from the head model could cause significant localization errors mostly for sources closer to the inner surface of the skull. With a less pronounced effect compared to the CSF exclusion, the discrimination between GM and WM also widely affected all sources, especially those located in deeper structures. The exclusion of the fontanels from the head model led to source localization errors for sources located in areas beneath the fontanels. Our finding clearly shows that the CSF inclusion and GM/WM distinction in EEG inverse modeling can substantially reduce EEG source localization errors. Moreover, fontanels should be included in neonatal head models, particularly in source localization applications, in which sources of interest are located beneath or in vicinity of fontanels.

SIGNIFICANCE

Our findings have practical implications for a better understanding of the impact of head model complexities on the accuracy of EEG source localization in neonates.

Mechanism of Coup and Contrecoup Injuries Induced by a Knock-Out Punch

Primary Objective: The interaction of cerebrospinal fluid with the brain parenchyma in an impact scenario is studied. Research Design: A computational fluid-structure interaction model is used to simulate the interaction of cerebrospinal fluid with a comprehensive brain model. Methods and Procedures: The method of smoothed particle hydrodynamics is used to simulate the fluid flow, induced by the impact, simultaneously with finite element analysis to solve the large deformations in the brain model. Main Outcomes and Results: Mechanism of injury resulting in concussion is demonstrated. The locations with the highest stress values on the brain parenchyma are shown. Conclusions: Our simulations found that the damage to the brain resulting from the contrecoup injury is more severe than that resulting from the coup injury. Additionally, we show that the contrecoup injury does not always appear on the side opposite from where impact occurs.

A hybrid analytical–numerical algorithm for determining the neuronal current via electroencephalography

Specific mental processes are associated with brain activation of a unique form, which are, in turn, expressed via the generation of specific neuronal electric currents. Electroencephalography (EEG) is based on measurements on the scalp of the electric potential generated by the neuronal current flowing in the cortex. This specific form of EEG data has been employed for a plethora of medical applications, from sleep studies to diagnosing focal epilepsy. In recent years, there have been efforts to use EEG data for a more ambitious purpose, namely to determine the underlying neuronal current. Although it has been known since 1853, from the studies by Helmholtz, that the knowledge of the electric potential of the external surface of a conductor is insufficient for the determination of the electric current that gave rise to this potential, the important question of which part of the current can actually be determined from the knowledge of this potential remained open until work published in 1997, when it was shown that EEG provides information only about the irrotational part of the current, which will be denoted by Ψ ; moreover, an explicit formula was derived in the above work relating this part of the current, the measured electric potential, and a certain auxiliary function, v s , that depends on the geometry of the various compartments of the brain–head system and their conductivities. In the present paper: (i) Motivated by recent results which show that, in the case of ellipsoidal geometry, the assumption of the L 2 minimization of the current yields a unique solution, we derive an analogous analytic formula characterizing this minimization for arbitrary geometry. (ii) We show that the above auxiliary function can be computed numerically via a line integral from the values of a related function v s computed via OpenMEEG; moreover, we propose an alternative approach to computing the auxiliary function v s based on the construction of a certain surrogate model. (iii) By expanding Ψ in terms of an inverse multiquadric radial basis we implement the relevant formulae numerically. The above algorithm performs well for synthetic data; its implementation with real data only requires the knowledge of the coordinates of the positions where the given EEG data are obtained.

Benchmarking transcranial electrical stimulation finite element models: a comparison study

OBJECTIVE

To compare field measure differences in simulations of transcranial electrical stimulation (tES) generated by variations in finite element (FE) models due to boundary condition specification, use of tissue compartment smoothing filters, and use of free or structured tetrahedral meshes based on magnetic resonance imaging (MRI) data.

APPROACH

A structural MRI head volume was acquired at 1 mm³ resolution and segmented into ten tissue compartments. Predicted current densities and electric fields were computed in segmented models using modeling pipelines involving either an in-house (block) or a commercial platform commonly used in previous FE tES studies involving smoothed compartments and free meshing procedures (smooth). The same boundary conditions were used for both block and smooth pipelines. Differences caused by varying boundary conditions were examined using a simple geometry. Percentage differences of median current density values in five cortical structures were compared between the two pipelines for three electrode montages (F3-right supraorbital, T7-T8 and Cz-Oz).

MAIN RESULTS

Use of boundary conditions commonly used in previous tES FE studies produced asymmetric current density profiles in the simple geometry. In head models, median current density differences produced by the two pipelines, using the same boundary conditions, were up to 6% (isotropic) and 18% (anisotropic) in structures targeted by each montage. Tangential electric field measures calculated via either pipeline were within the range of values reported in the literature, when averaged over cortical surface patches.

SIGNIFICANCE

Apparently equivalent boundary settings may affect predicted current density outcomes and care must be taken in their specification. Smoothing FE model compartments may not be necessary, and directly translated, voxellated tissue boundaries at 1 mm³ resolution may be sufficient for use in tES FE studies, greatly reducing processing times. The findings here may be used to inform future current density modeling studies.

A Finite Element Solution of the Forward Problem in EEG for Multipolar Sources

Multipolar source models have been presented in the context of electro/magnetoencephalography (E/MEG) to compensate for the limitations of the classical equivalent current dipole to represent realistic generators of brain activity. Although there exist several reports accounting for the advantages of multipolar components over single dipoles, there is still no available numerical implementation in fully personalized scenarios. In this paper, we present, for the first time, a finite element framework for simulating EEG signals generated by multipolar current sources in individualized, heterogeneous, and anisotropic head models. This formulation is based on the subtraction approach, guaranteeing the existence and uniqueness of the solution. In particular, we analyze the cases of monopolar, dipolar, and quadrupolar source components, for which we study their performance in idealized and realistic head models. Numerical solutions are compared with analytical formulas in multi-layered spherical models. Such formulas are available in the case of monopolar and dipolar sources, and here derived for the quadrupolar components. We finally illustrate their advantages in the description of extended current generators using a realistic head model. The framework presented here enables further analysis towards the estimation of biophysically principled source parameters from standard E/MEG experiments.

The effects of asymmetric volume conductor modeling on non-invasive fetal ECG extraction

OBJECTIVE

Non-invasive fetal electrocardiography (NI-FECG) shows promise for capturing novel physiological information that may indicate signs of fetal distress. However, significant deterioration in NI-FECG signal quality occurs during the presence of a highly non-conductive layer known as vernix caseosa which forms on the fetal body surface beginning in approximately the 28th week of gestation. This work investigates asymmetric modeling of vernix caseosa and other maternal-fetal tissues in accordance with clinical observations and assesses their impacts for NI-FECG signal processing.

APPROACH

We develop a process for simulating dynamic maternal-fetal abdominal ECG mixtures using a synthetic cardiac source model embedded in a finite element volume conductor. Using this process, changes in NI-FECG signal morphology are assessed in an extensive set of finite element models including spatially variable distributions of vernix caseosa.

MAIN RESULTS

Our simulations show that volume conductor asymmetry can result in over 70% error in the observed T/QRS ratio and significant changes to signal morphology compared to a homogeneous volume conductor model. Volume conductor effects must be considered when analyzing T/QRS ratios obtained via NI-FECG and should be considered in future algorithm benchmarks using simulated data.

SIGNIFICANCE

This work shows that without knowledge of the influence of volume conductor effects, clinical evaluation of the T/QRS ratio derived via NI-FECG should be avoided.

New Strategy for Finite Element Mesh Generation for Accurate Solutions of Electroencephalography Forward Problems

The finite element method (FEM) is a numerical method that is often used for solving electroencephalography (EEG) forward problems involving realistic head models. In this study, FEM solutions obtained using three different mesh structures, namely coarse, densely refined, and adaptively refined meshes, are compared. The simulation results showed that the accuracy of FEM solutions could be significantly enhanced by adding a small number of elements around regions with large estimated errors. Moreover, it was demonstrated that the adaptively refined regions were always near the current dipole sources, suggesting that selectively generating additional elements around the cortical surface might be a new promising strategy for more efficient FEM-based EEG forward analysis.

Sensitivity of the Projected Subtraction Approach to Mesh Degeneracies and Its Impact on the Forward Problem in EEG

OBJECTIVE

Subtraction-based techniques are known for being theoretically rigorous and accurate methods for solving the forward problem in electroencephalography (EEG-FP) by means of the finite-element method. Within them, the projected subtraction (PS) approach is generally adopted because of its computational efficiency. Although this technique received the attention of the community, its sensitivity to degenerated elements is still poorly understood. In this paper, we investigate the impact of low-quality tetrahedra on the results computed with the PS approach.

METHODS

We derived upper bounds on the relative error of the element source vector as a function of geometrical features describing the tetrahedral discretization of the domain. These error bounds were then utilized for showing the instability of the PS method with regards to the mesh quality. To overcome this issue, we proposed an alternative technique, coined projected gradient subtraction (PGS) approach, that exploits the stability of the corresponding bounds.

RESULTS

Computer simulations showed that the PS method is extremely sensitive to the mesh shape and size, leading to unacceptable solutions of the EEG-FP in case of using suboptimal tessellations. This was not the case of the PGS approach, which led to stable and accurate results in a comparable amount of time.

CONCLUSION

Solutions of the EEG-FP computed with the PS method are highly sensitive to degenerated elements. Such errors can be mitigated by the PGS approach, which showed better performance than the PS technique.

SIGNIFICANCE

The PGS is an efficient method for computing high-quality lead field matrices even in the presence of degenerated elements.

The FieldTrip-SimBio pipeline for EEG forward solutions

Background

Accurately solving the electroencephalography (EEG) forward problem is crucial for precise EEG source analysis. Previous studies have shown that the use of multicompartment head models in combination with the finite element method (FEM) can yield high accuracies both numerically and with regard to the geometrical approximation of the human head. However, the workload for the generation of multicompartment head models has often been too high and the use of publicly available FEM implementations too complicated for a wider application of FEM in research studies. In this paper, we present a MATLAB-based pipeline that aims to resolve this lack of easy-to-use integrated software solutions. The presented pipeline allows for the easy application of five-compartment head models with the FEM within the FieldTrip toolbox for EEG source analysis.

Methods

The FEM from the SimBio toolbox, more specifically the St. Venant approach, was integrated into the FieldTrip toolbox. We give a short sketch of the implementation and its application, and we perform a source localization of somatosensory evoked potentials (SEPs) using this pipeline. We then evaluate the accuracy that can be achieved using the automatically generated five-compartment hexahedral head model [skin, skull, cerebrospinal fluid (CSF), gray matter, white matter] in comparison to a highly accurate tetrahedral head model that was generated on the basis of a semiautomatic segmentation with very careful and time-consuming manual corrections.

Results

The source analysis of the SEP data correctly localizes the P20 component and achieves a high goodness of fit. The subsequent comparison to the highly detailed tetrahedral head model shows that the automatically generated five-compartment head model performs about as well as a highly detailed four-compartment head model (skin, skull, CSF, brain). This is a significant improvement in comparison to a three-compartment head model, which is frequently used in praxis, since the importance of modeling the CSF compartment has been shown in a variety of studies.

Conclusion

The presented pipeline facilitates the use of five-compartment head models with the FEM for EEG source analysis. The accuracy with which the EEG forward problem can thereby be solved is increased compared to the commonly used three-compartment head models, and more reliable EEG source reconstruction results can be obtained.

Electronic supplementary material

The online version of this article (10.1186/s12938-018-0463-y) contains supplementary material, which is available to authorized users.

Robust EEG/MEG Based Functional Connectivity with the Envelope of the Imaginary Coherence: Sensor Space Analysis

The brain's functional connectivity (FC) estimated at sensor level from electromagnetic (EEG/MEG) signals can provide quick and useful information towards understanding cognition and brain disorders. Volume conduction (VC) is a fundamental issue in FC analysis due to the effects of instantaneous correlations. FC methods based on the imaginary part of the coherence (iCOH) of any two signals are readily robust to VC effects, but neglecting the real part of the coherence leads to negligible FC when the processes are truly connected but with zero or π-phase (modulus 2π) interaction. We ameliorate this issue by proposing a novel method that implements an envelope of the imaginary coherence (EIC) to approximate the coherence estimate of supposedly active underlying sources. We compare EIC with state-of-the-art FC measures that included lagged coherence, iCOH, phase lag index (PLI) and weighted PLI (wPLI), using bivariate autoregressive and stochastic neural mass models. Additionally, we create realistic simulations where three and five regions were mapped on a template cortical surface and synthetic MEG signals were obtained after computing the electromagnetic leadfield. With this simulation and comparison study, we also demonstrate the feasibility of sensor FC analysis using receiver operating curve analysis whilst varying the signal's noise level. However, these results should be interpreted with caution given the known limitations of the sensor-based FC approach. Overall, we found that EIC and iCOH demonstrate superior results with most accurate FC maps. As they complement each other in different scenarios, that will be important to study normal and diseased brain activity.

Changes in network connectivity during motor imagery and execution

Background

Recent studies of functional or effective connectivity in the brain have reported that motor-related brain regions were activated during motor execution and motor imagery, but the relationship between motor and cognitive areas has not yet been completely understood. The objectives of our study were to analyze the effective connectivity between motor and cognitive networks in order to define network dynamics during motor execution and motor imagery in healthy individuals. Second, we analyzed the differences in effective connectivity between correct and incorrect responses during motor execution and imagery using dynamic causal modeling (DCM) of electroencephalography (EEG) data.

Method

Twenty healthy subjects performed a sequence of finger tapping trials using either motor execution or motor imagery, and the performances were recorded. Changes in effective connectivity between the primary motor cortex (M1), supplementary motor area (SMA), premotor cortex (PMC), and dorsolateral prefrontal cortex (DLPFC) were estimated using dynamic causal modeling. Bayesian model averaging with family-level inference and fixed-effects analysis was applied to determine the most likely connectivity model for these regions.

Results

Motor execution and imagery showed inputs to distinct brain regions, the premotor cortex and the supplementary motor area, respectively. During motor execution, the coupling strength of a feedforward network from the DLPFC to the PMC was greater than that during motor imagery. During motor imagery, the coupling strengths of a feedforward network from the PMC to the SMA and of a feedback network from M1 to the PMC were higher than that during motor execution. In imagined movement, although there were connectivity differences between correct and incorrect task responses, each motor imagery task that included correct and incorrect responses showed similar network connectivity characteristics. Correct motor imagery responses showed connectivity from the PMC to the DLPFC, while the incorrect responses had characteristic connectivity from the SMA to the DLPFC.

Conclusions

These findings provide an understanding of effective connectivity between motor and cognitive areas during motor execution and imagery as well as the basis for future connectivity studies for patients with stroke.

The Effect of Head Model Simplification on Beamformer Source Localization

Beamformers are a widely-used tool in brain analysis with magnetoencephalography (MEG) and electroencephalography (EEG). For the construction of the beamformer filters realistic head volume conductor modeling is necessary for accurately computing the EEG and MEG leadfields, i.e., for solving the EEG and MEG forward problem. In this work, we investigate the influence of including realistic head tissue compartments into a finite element method (FEM) model on the beamformer's localization ability. Specifically, we investigate the effect of including cerebrospinal fluid, gray matter, and white matter distinction, as well as segmenting the skull bone into compacta and spongiosa, and modeling white matter anisotropy. We simulate an interictal epileptic measurement with white sensor noise. Beamformer filters are constructed with unit gain, unit array gain, and unit noise gain constraint. Beamformer source positions are determined by evaluating power and excess sample kurtosis (g₂) of the source-waveforms at all source space nodes. For both modalities, we see a strong effect of modeling the cerebrospinal fluid and white and gray matter. Depending on the source position, both effects can each be in the magnitude of centimeters, rendering their modeling necessary for successful localization. Precise skull modeling mainly effected the EEG up to a few millimeters, while both modalities could profit from modeling white matter anisotropy to a smaller extent of 5-10 mm. The unit noise gain or neural activity index beamformer behaves similarly to the array gain beamformer when noise strength is sufficiently high. Variance localization seems more robust against modeling errors than kurtosis.

Zoomed MRI Guided by Combined EEG/MEG Source Analysis: A Multimodal Approach for Optimizing Presurgical Epilepsy Work-up and its Application in a Multi-focal Epilepsy Patient Case Study

In recent years, the use of source analysis based on electroencephalography (EEG) and magnetoencephalography (MEG) has gained considerable attention in presurgical epilepsy diagnosis. However, in many cases the source analysis alone is not used to tailor surgery unless the findings are confirmed by lesions, such as, e.g., cortical malformations in MRI. For many patients, the histology of tissue resected from MRI negative epilepsy shows small lesions, which indicates the need for more sensitive MR sequences. In this paper, we describe a technique to maximize the synergy between combined EEG/MEG (EMEG) source analysis and high resolution MRI. The procedure has three main steps: (1) construction of a detailed and calibrated finite element head model that considers the variation of individual skull conductivities and white matter anisotropy, (2) EMEG source analysis performed on averaged interictal epileptic discharges (IED), (3) high resolution (0.5 mm) zoomed MR imaging, limited to small areas centered at the EMEG source locations. The proposed new diagnosis procedure was then applied in a particularly challenging case of an epilepsy patient: EMEG analysis at the peak of the IED coincided with a right frontal focal cortical dysplasia (FCD), which had been detected at standard 1 mm resolution MRI. Of higher interest, zoomed MR imaging (applying parallel transmission, 'ZOOMit') guided by EMEG at the spike onset revealed a second, fairly subtle, FCD in the left fronto-central region. The evaluation revealed that this second FCD, which had not been detectable with standard 1 mm resolution, was the trigger of the seizures.

Combined EMD-sLORETA Analysis of EEG Data Collected during a Contour Integration Task

Lately, Ensemble Empirical Mode Decomposition (EEMD) techniques receive growing interest in biomedical data analysis. Event-Related Modes (ERMs) represent features extracted by an EEMD from electroencephalographic (EEG) recordings. We present a new approach for source localization of EEG data based on combining ERMs with inverse models. As the first step, 64 channel EEG recordings are pooled according to six brain areas and decomposed, by applying an EEMD, into their underlying ERMs. Then, based upon the problem at hand, the most closely related ERM, in terms of frequency and amplitude, is combined with inverse modeling techniques for source localization. More specifically, the standardized low resolution brain electromagnetic tomography (sLORETA) procedure is employed in this work. Accuracy and robustness of the results indicate that this approach deems highly promising in source localization techniques for EEG data.

Electrical Stimulation of the Human Cerebral Cortex by Extracranial Muscle Activity: Effect Quantification With Intracranial EEG and FEM Simulations

OBJECTIVE

Electric fields (EF) of approx. 0.2 V/m have been shown to be sufficiently strong to both modulate neuronal activity in the cerebral cortex and have measurable effects on cognitive performance. We hypothesized that the EF caused by the electrical activity of extracranial muscles during natural chewing may reach similar strength in the cerebral cortex and hence might act as an endogenous modality of brain stimulation. Here, we present first steps toward validating this hypothesis.

METHODS

Using a realistic volume conductor head model of an epilepsy patient having undergone intracranial electrode placement and utilizing simultaneous intracranial and extracranial electrical recordings during chewing, we derive predictions about the chewing-related cortical EF strength to be expected in healthy individuals.

RESULTS

We find that in the region of the temporal poles, the expected EF strength may reach amplitudes in the order of 0.1-1 V/m.

CONCLUSION

The cortical EF caused by natural chewing could be large enough to modulate ongoing neural activity in the cerebral cortex and influence cognitive performance.

SIGNIFICANCE

The present study lends first support for the assumption that extracranial muscle activity might represent an endogenous source of electrical brain stimulation. This offers a new potential explanation for the puzzling effects of gum chewing on cognition, which have been repeatedly reported in the literature.

The Unfitted Discontinuous Galerkin Method for Solving the EEG Forward Problem

OBJECTIVE

The purpose of this study is to introduce and evaluate the unfitted discontinuous Galerkin finite element method (UDG-FEM) for solving the electroencephalography (EEG) forward problem.

METHODS

This new approach for source analysis does not use a geometry conforming volume triangulation, but instead uses a structured mesh that does not resolve the geometry. The geometry is described using level set functions and is incorporated implicitly in its mathematical formulation. As no triangulation is necessary, the complexity of a simulation pipeline and the need for manual interaction for patient-specific simulations can be reduced and is comparable with that of the FEM for hexahedral meshes. In addition, it maintains conservation laws on a discrete level. Here, we present the theory for UDG-FEM forward modeling, its verification using quasi-analytical solutions in multilayer sphere models and an evaluation in a comparison with a discontinuous Galerkin (DG-FEM) method on hexahedral and on conforming tetrahedral meshes. We furthermore apply the UDG-FEM forward approach in a realistic head model simulation study.

RESULTS

The results show convergence to the quasi-analytical solution and indicate a good accuracy of UDG-FEM. UDG-FEM performs comparable or even better than DG-FEM on a conforming tetrahedral mesh while providing a less complex simulation pipeline. When compared to DG-FEM on hexahedral meshes, an overall better accuracy is achieved.

CONCLUSION

The UDG-FEM approach is an accurate, flexible, and promising method to solve the EEG forward problem.

SIGNIFICANCE

This study shows the first application of the UDG-FEM approach to the EEG forward problem.

Combined EEG/MEG can outperform single modality EEG or MEG source reconstruction in presurgical epilepsy diagnosis

We investigated two important means for improving source reconstruction in presurgical epilepsy diagnosis. The first investigation is about the optimal choice of the number of epileptic spikes in averaging to (1) sufficiently reduce the noise bias for an accurate determination of the center of gravity of the epileptic activity and (2) still get an estimation of the extent of the irritative zone. The second study focuses on the differences in single modality EEG (80-electrodes) or MEG (275-gradiometers) and especially on the benefits of combined EEG/MEG (EMEG) source analysis. Both investigations were validated with simultaneous stereo-EEG (sEEG) (167-contacts) and low-density EEG (ldEEG) (21-electrodes). To account for the different sensitivity profiles of EEG and MEG, we constructed a six-compartment finite element head model with anisotropic white matter conductivity, and calibrated the skull conductivity via somatosensory evoked responses. Our results show that, unlike single modality EEG or MEG, combined EMEG uses the complementary information of both modalities and thereby allows accurate source reconstructions also at early instants in time (epileptic spike onset), i.e., time points with low SNR, which are not yet subject to propagation and thus supposed to be closer to the origin of the epileptic activity. EMEG is furthermore able to reveal the propagation pathway at later time points in agreement with sEEG, while EEG or MEG alone reconstructed only parts of it. Subaveraging provides important and accurate information about both the center of gravity and the extent of the epileptogenic tissue that neither single nor grand-averaged spike localizations can supply.

Influence of the head model on EEG and MEG source connectivity analyses

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

Combining EEG and MEG for the reconstruction of epileptic activity using a calibrated realistic volume conductor model

To increase the reliability for the non-invasive determination of the irritative zone in presurgical epilepsy diagnosis, we introduce here a new experimental and methodological source analysis pipeline that combines the complementary information in EEG and MEG, and apply it to data from a patient, suffering from refractory focal epilepsy. Skull conductivity parameters in a six compartment finite element head model with brain anisotropy, constructed from individual MRI data, are estimated in a calibration procedure using somatosensory evoked potential (SEP) and field (SEF) data. These data are measured in a single run before acquisition of further runs of spontaneous epileptic activity. Our results show that even for single interictal spikes, volume conduction effects dominate over noise and need to be taken into account for accurate source analysis. While cerebrospinal fluid and brain anisotropy influence both modalities, only EEG is sensitive to skull conductivity and conductivity calibration significantly reduces the difference in especially depth localization of both modalities, emphasizing its importance for combining EEG and MEG source analysis. On the other hand, localization differences which are due to the distinct sensitivity profiles of EEG and MEG persist. In case of a moderate error in skull conductivity, combined source analysis results can still profit from the different sensitivity profiles of EEG and MEG to accurately determine location, orientation and strength of the underlying sources. On the other side, significant errors in skull modeling are reflected in EEG reconstruction errors and could reduce the goodness of fit to combined datasets. For combined EEG and MEG source analysis, we therefore recommend calibrating skull conductivity using additionally acquired SEP/SEF data.