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

Standardized hierarchical adaptive Lp regression for noise robust focal epilepsy source reconstructions

OBJECTIVE

To investigate the ability of standardization to reduce source localization errors and measurement noise uncertainties for hierarchical Bayesian algorithms with L1- and L2-norms as priors in electroencephalography and magnetoencephalography of focal epilepsy.

METHODS

Description of the standardization methodology relying on the Hierarchical Bayesian framework, referred to as the Standardized Hierarchical Adaptive Lp-norm Regularization (SHALpR). The performance was tested using real data from two focal epilepsy patients. Simulated data that resembled the available real data was constructed for further localization and noise robustness investigation.

RESULTS

The proposed algorithms were compared to their non-standardized counterparts, Standardized low-resolution brain electromagnetic tomography, Standardized Shrinking LORETA-FOCUSS, and Dynamic statistical parametric maps. Based on the simulations, the standardized Hierarchical adaptive algorithm using L2-norm was noise robust for 10 dB signal-to-noise ratio (SNR), whereas the L1-norm prior worked robustly also with 5 dB SNR. The accuracy of the standardized L1-normed methodology to localize focal activity was under 1 cm for both patients.

CONCLUSIONS

Numerical results of the proposed methodology display improved localization and noise robustness. The proposed methodology also outperformed the compared methods when dealing with real data.

SIGNIFICANCE

The proposed standardized methodology, especially when employing the L1-norm, could serve as a valuable assessment tool in surgical decision-making.

How to assess the accuracy of volume conduction models? A validation study with stereotactic EEG data

Introduction

Volume conduction models of the human head are used in various neuroscience fields, such as for source reconstruction in EEG and MEG, and for modeling the effects of brain stimulation. Numerous studies have quantified the accuracy and sensitivity of volume conduction models by analyzing the effects of the geometrical and electrical features of the head model, the sensor model, the source model, and the numerical method. Most studies are based on simulations as it is hard to obtain sufficiently detailed measurements to compare to models. The recording of stereotactic EEG during electric stimulation mapping provides an opportunity for such empirical validation.

Methods

In the study presented here, we used the potential distribution of volume-conducted artifacts that are due to cortical stimulation to evaluate the accuracy of finite element method (FEM) volume conduction models. We adopted a widely used strategy for numerical comparison, i.e., we fixed the geometrical description of the head model and the mathematical method to perform simulations, and we gradually altered the head models, by increasing the level of detail of the conductivity profile. We compared the simulated potentials at different levels of refinement with the measured potentials in three epilepsy patients.

Results

Our results show that increasing the level of detail of the volume conduction head model only marginally improves the accuracy of the simulated potentials when compared to in-vivo sEEG measurements. The mismatch between measured and simulated potentials is, throughout all patients and models, maximally 40 microvolts (i.e., 10% relative error) in 80% of the stimulation-recording combination pairs and it is modulated by the distance between recording and stimulating electrodes.

Discussion

Our study suggests that commonly used strategies used to validate volume conduction models based solely on simulations might give an overly optimistic idea about volume conduction model accuracy. We recommend more empirical validations to be performed to identify those factors in volume conduction models that have the highest impact on the accuracy of simulated potentials. We share the dataset to allow researchers to further investigate the mismatch between measurements and FEM models and to contribute to improving volume conduction models.

Validating EEG source imaging using intracranial electrical stimulation

Electrical source imaging is used in presurgical epilepsy evaluation and in cognitive neurosciences to localize neuronal sources of brain potentials recorded on EEG. This study evaluates the spatial accuracy of electrical source imaging for known sources, using electrical stimulation potentials recorded on simultaneous stereo-EEG and 37-electrode scalp EEG, and identifies factors determining the localization error. In 11 patients undergoing simultaneous stereo-EEG and 37-electrode scalp EEG recordings, sequential series of 99-110 biphasic pulses (2 ms pulse width) were applied by bipolar electrical stimulation on adjacent contacts of implanted stereo-EEG electrodes. The scalp EEG correlates of stimulation potentials were recorded with a sampling rate of 30 kHz. Electrical source imaging of averaged stimulation potentials was calculated utilizing a dipole source model of peak stimulation potentials based on individual four-compartment finite element method head models with various skull conductivities (range from 0.0413 to 0.001 S/m). Fitted dipoles with a goodness of fit of ≥80% were included in the analysis. The localization error was calculated using the Euclidean distance between the estimated dipoles and the centre point of adjacent stimulating contacts. A total of 3619 stimulation locations, respectively, dipole localizations, were included in the evaluation. Mean localization errors ranged from 10.3 to 26 mm, depending on source depth and selected skull conductivity. The mean localization error increased with an increase in source depth (r(3617) = [0.19], P = 0.000) and decreased with an increase in skull conductivity (r(3617) = [-0.26], P = 0.000). High skull conductivities (0.0413-0.0118 S/m) yielded significantly lower localization errors for all source depths. For superficial sources (<20 mm from the inner skull), all skull conductivities yielded insignificantly different localization errors. However, for deeper sources, in particular >40 mm, high skull conductivities of 0.0413 and 0.0206 S/m yielded significantly lower localization errors. In relation to stimulation locations, the majority of estimated dipoles moved outward-forward-downward to inward-forward-downward with a decrease in source depth and an increase in skull conductivity. Multivariate analysis revealed that an increase in source depth, number of skull holes and white matter volume, while a decrease in skull conductivity independently led to higher localization error. This evaluation of electrical source imaging accuracy using artificial patterns with a high signal-to-noise ratio supports its application in presurgical epilepsy evaluation and cognitive neurosciences. In our artificial potential model, optimizing the selected skull conductivity minimized the localization error. Future studies should examine if this accounts for true neural signals.

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.

Can individually targeted and optimized multi-channel tDCS outperform standard bipolar tDCS in stimulating the primary somatosensory cortex?

BACKGROUND

Transcranial direct current stimulation (tDCS) has emerged as a non-invasive neuro-modulation technique. Most studies show that anodal tDCS increases cortical excitability, however, with variable outcomes. Previously, we have shown in computer simulations that our multi-channel tDCS (mc-tDCS) approach, the distributed constrained maximum intensity (D-CMI) method can potentially lead to better controlled tDCS results due to the improved directionality of the injected current at the target side for individually optimized D-CMI montages.

OBJECTIVE

In this study, we test the application of the D-CMI approach in an experimental study to stimulate the somatosensory P20/N20 target source in Brodmann area 3b and compare it with standard bipolar tDCS and sham conditions.

METHODS

We applied anodal D-CMI, the standard bipolar and D-CMI based Sham tDCS for 10 min to target the 20 ms post-stimulus somatosensory P20/N20 target brain source in Brodmann area 3b reconstructed using combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis in realistic head models with calibrated skull conductivity in a group-study with 13 subjects. Finger-stimulated somatosensory evoked fields (SEF) were recorded and the component at 20 ms post-stimulus (M20) was analyzed before and after the application of the three tDCS conditions in order to read out the stimulation effect on Brodmann area 3b.

RESULTS

Analysis of the finger stimulated SEF M20 peak before (baseline) and after tDCS shows a significant increase in source amplitude in Brodmann area 3b for D-CMI (6-16 min after tDCS), while no significant effects are found for standard bipolar (6-16 min after tDCS) and sham (6-16 min after tDCS) stimulation conditions. For the later time courses (16-26 and 27-37 min post-stimulation), we found a significant decrease in M20 peak source amplitude for standard bipolar and sham tDCS, while there was no effect for D-CMI.

CONCLUSION

Our results indicate that targeted and optimized, and thereby highly individualized, mc-tDCS can outperform standard bipolar stimulation and lead to better control over stimulation outcomes with, however, a considerable amount of additional work compared to standard bipolar tDCS.

Hippocampal dentate gyrus integrity revealed with ultrahigh resolution diffusion imaging predicts memory performance in older adults

Medial temporal lobe (MTL) atrophy is a core feature of age-related cognitive decline and Alzheimer's disease (AD). While regional volumes and thickness are often used as a proxy for neurodegeneration, they lack the sensitivity to serve as an accurate diagnostic test and indicate advanced neurodegeneration. Here, we used a submillimeter resolution diffusion weighted MRI sequence (ZOOMit) to quantify microstructural properties of hippocampal subfields in older adults (63-98 years old) using tensor derived measures: fractional anisotropy (FA) and mean diffusivity (MD). We demonstrate that the high-resolution sequence, and not a standard resolution sequence, identifies dissociable profiles for CA1, dentate gyrus (DG), and the collateral sulcus. Using ZOOMit, we show that advanced age is associated with increased MD of the CA1 and DG as well as decreased FA of the DG. Increased MD of the DG, reflecting decreased cellular density, mediated the relationship between age and word list recall. Further, increased MD in the DG, but not DG volume, was linked to worse spatial pattern separation. Our results demonstrate that ultrahigh-resolution diffusion imaging enables the detection of microstructural differences in hippocampal subfield integrity and will lead to novel insights into the mechanisms of age-related memory loss.

Black Line Sign in Focal Cortical Dysplasia IIB: A 7T MRI and Electroclinicopathologic Study

BACKGROUND AND OBJECTIVES

We aim to provide detailed imaging-electroclinicopathologic characterization of the black line sign, a novel MRI marker for focal cortical dysplasia (FCD) IIB.

METHODS

7T T2*-weighted gradient-echo (T2*w-GRE) images were retrospectively reviewed in a consecutive cohort of patients with medically intractable epilepsy with pathology-proven FCD II, for the occurrence of the black line sign. We examined the overlap between the black line region and the seizure-onset zone (SOZ) defined by intracranial EEG (ICEEG) and additionally assessed whether complete inclusion of the black line region in the surgical resection was associated with postoperative seizure freedom. The histopathologic specimen was aligned with the MRI to investigate the pathologic underpinning of the black line sign. Region-of-interest-based quantitative MRI (qMRI) analysis on the 7T T1 map was performed in the black line region, entire lesional gray matter (GM), and contralateral/ipsilateral normal gray and white matter (WM).

RESULTS

We included 20 patients with FCD II (14 IIB and 6 IIA). The black line sign was identified in 12/14 (85.7%) of FCD IIB and 0/6 of FCD IIA on 7T T2*w-GRE. The black line region was highly concordant with the ICEEG-defined SOZ (5/7 complete and 2/7 partial overlap). Seizure freedom was seen in 8/8 patients whose black line region was completely included in the surgical resection; in the 2 patients whose resection did not completely include the black line region, both had recurring seizures. Inclusion of the black line region in the surgical resection was significantly associated with seizure freedom (p = 0.02). QMRI analyses showed that the T1 mean value of the black line region was significantly different from the WM (p < 0.001), but similar to the GM. Well-matched histopathologic slices in one case revealed accumulated dysmorphic neurons and balloon cells in the black line region.

DISCUSSION

The black line sign may serve as a noninvasive marker for FCD IIB. Both MRI-pathology and qMRI analyses suggest that the black line region was an abnormal GM component within the FCD. Being highly concordant with ICEEG-defined SOZ and significantly associated with seizure freedom when included in resection, the black line sign may contribute to the planning of ICEEG/surgery of patients with medically intractable epilepsy with FCD IIB.

CLASSIFICATION OF EVIDENCE

This study provides Class II evidence that in individuals with intractable focal epilepsy undergoing resection who have a 7T MRI with adequate image quality, the presence of the black line sign may suggest FCD IIB, be concordant with SOZ from ICEEG, and be associated with more seizure freedom if fully included in resection.

Validating EEG, MEG and Combined MEG and EEG Beamforming for an Estimation of the Epileptogenic Zone in Focal Cortical Dysplasia

MEG and EEG source analysis is frequently used for the presurgical evaluation of pharmacoresistant epilepsy patients. The source localization of the epileptogenic zone depends, among other aspects, on the selected inverse and forward approaches and their respective parameter choices. In this validation study, we compare the standard dipole scanning method with two beamformer approaches for the inverse problem, and we investigate the influence of the covariance estimation method and the strength of regularization on the localization performance for EEG, MEG, and combined EEG and MEG. For forward modelling, we investigate the difference between calibrated six-compartment and standard three-compartment head modelling. In a retrospective study, two patients with focal epilepsy due to focal cortical dysplasia type IIb and seizure freedom following lesionectomy or radiofrequency-guided thermocoagulation (RFTC) used the distance of the localization of interictal epileptic spikes to the resection cavity resp. RFTC lesion as reference for good localization. We found that beamformer localization can be sensitive to the choice of the regularization parameter, which has to be individually optimized. Estimation of the covariance matrix with averaged spike data yielded more robust results across the modalities. MEG was the dominant modality and provided a good localization in one case, while it was EEG for the other. When combining the modalities, the good results of the dominant modality were mostly not spoiled by the weaker modality. For appropriate regularization parameter choices, the beamformer localized better than the standard dipole scan. Compared to the importance of an appropriate regularization, the sensitivity of the localization to the head modelling was smaller, due to similar skull conductivity modelling and the fixed source space without orientation constraint.

Individually optimized multi-channel tDCS for targeting somatosensory cortex

OBJECTIVE

Transcranial direct current stimulation (tDCS) is a non-invasive neuro-modulation technique that delivers current through the scalp by a pair of patch electrodes (2-Patch). This study proposes a new multi-channel tDCS (mc-tDCS) optimization method, the distributed constrained maximum intensity (D-CMI) approach. For targeting the P20/N20 somatosensory source at Brodmann area 3b, an integrated combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis is used with individualized skull conductivity calibrated realistic head modeling.

METHODS

Simulated electric fields (EF) for our new D-CMI method and the already known maximum intensity (MI), alternating direction method of multipliers (ADMM) and 2-Patch methods were produced and compared for the individualized P20/N20 somatosensory target for 10 subjects.

RESULTS

D-CMI and MI showed highest intensities parallel to the P20/N20 target compared to ADMM and 2-Patch, with ADMM achieving highest focality. D-CMI showed a slight reduction in intensity compared to MI while reducing side effects and skin level sensations by current distribution over multiple stimulation electrodes.

CONCLUSION

Individualized D-CMI montages are preferred for our follow up somatosensory experiment to provide a good balance between high current intensities at the target and reduced side effects and skin sensations.

SIGNIFICANCE

An integrated combined MEG and EEG source analysis with D-CMI montages for mc-tDCS stimulation potentially can improve control, reproducibility and reduce sensitivity differences between sham and real stimulations.

A comprehensive study on electroencephalography and magnetoencephalography sensitivity to cortical and subcortical sources

Signal-to-noise ratio (SNR) maps are a good way to visualize electroencephalography (EEG) and magnetoencephalography (MEG) sensitivity. SNR maps extend the knowledge about the modulation of EEG and MEG signals by source locations and orientations and can therefore help to better understand and interpret measured signals as well as source reconstruction results thereof. Our work has two main objectives. First, we investigated the accuracy and reliability of EEG and MEG finite element method (FEM)-based sensitivity maps for three different head models, namely an isotropic three and four-compartment and an anisotropic six-compartment head model. As a result, we found that ignoring the cerebrospinal fluid leads to an overestimation of EEG SNR values. Second, we examined and compared EEG and MEG SNR mappings for both cortical and subcortical sources and their modulation by source location and orientation. Our results for cortical sources show that EEG sensitivity is higher for radial and deep sources and MEG for tangential ones, which are the majority of sources. As to the subcortical sources, we found that deep sources with sufficient tangential source orientation are recordable by the MEG. Our work, which represents the first comprehensive study where cortical and subcortical sources are considered in highly detailed FEM-based EEG and MEG SNR mappings, sheds a new light on the sensitivity of EEG and MEG and might influence the decision of brain researchers or clinicians in their choice of the best modality for their experiment or diagnostics, respectively.

A novel method for calibrating head models to account for variability in conductivity and its evaluation in a sphere model

The accuracy in electroencephalography (EEG) and combined EEG and magnetoencephalography (MEG) source reconstructions as well as in optimized transcranial electric stimulation (TES) depends on the conductive properties assigned to the head model, and most importantly on individual skull conductivity. In this study, we present an automatic pipeline to calibrate head models with respect to skull conductivity based on the reconstruction of the P20/N20 response using somatosensory evoked potentials and fields. In order to validate in a well-controlled setup without interplay with numerical errors, we evaluate the accuracy of this algorithm in a 4-layer spherical head model using realistic noise levels as well as dipole sources at different eccentricities with strengths and orientations related to somatosensory experiments. Our results show that the reference skull conductivity can be reliably reconstructed for sources resembling the generator of the P20/N20 response. In case of erroneous assumptions on scalp conductivity, the resulting skull conductivity parameter counterbalances this effect, so that EEG source reconstructions using the fitted skull conductivity parameter result in lower errors than when using the standard value. We propose an automatized procedure to calibrate head models which only relies on non-invasive modalities that are available in a standard MEG laboratory, measures under in vivo conditions and in the low frequency range of interest. Calibrated head modeling can improve EEG and combined EEG/MEG source analysis as well as optimized TES.

Changing concepts in presurgical assessment for epilepsy surgery

Candidates for epilepsy surgery must undergo presurgical evaluation to establish whether and how surgical treatment can stop seizures without causing neurological deficits. Various techniques, including MRI, PET, single-photon emission CT, video-EEG, magnetoencephalography and invasive EEG, aim to identify the diseased brain tissue and the involved network. Recent technical and methodological developments, encompassing both advances in existing techniques and new combinations of technologies, are enhancing the ability to define the optimal resection strategy. Multimodal interpretation and predictive computer models are expected to aid surgical planning and patient counselling, and multimodal intraoperative guidance is likely to increase surgical precision. In this Review, we discuss how the knowledge derived from these new approaches is challenging our way of thinking about surgery to stop focal seizures. In particular, we highlight the importance of looking beyond the EEG seizure onset zone and considering focal epilepsy as a brain network disease in which long-range connections need to be taken into account. We also explore how new diagnostic techniques are revealing essential information in the brain that was previously hidden from view.

Localizing non-epileptiform abnormal brain function in children using high density EEG: Electric Source Imaging of focal slowing

BACKGROUND

Electric Source Imaging (ESI) of interictal epileptiform discharges (IED) is increasingly validated for localizing epileptic activity. In children, IED can be absent or multifocal even in cases of a focal epileptogenic zone and additional electrophysiological markers are needed. Here, we investigated ESI of pathological focal slowing (FS) recorded on EEG as a new localizing marker in children with drug-resistant epilepsy.

METHODS

We selected 15 children (median: 12; range: 4-18yrs), with high-density EEG (hdEEG), presurgical evaluation and surgical resection. One patient had a non-lesional MRI. ESI of patient-specific focal slow activity was performed (distributed linear inverse solution and individual head model). The maximal average power in the band of interest was considered as the source of focal slowing (ESI-FS). The Euclidian distance between ESI-FS and the resection (5 mm margin) was compared to the localization of maximal ESI of interictal epileptiform discharges (ESI-IED), interictal FDG-PET and ictal SPECT/SISCOM.

RESULTS

In 9/15 patients (60%), ESI of focal slowing (ESI-FS) was inside or ≤5 mm from resection margins. The remaining 6/15 cases had distances ≤15 mm. In 9/15 patients with interictal spikes, the ESI-IED was concordant with the resection. 6/15 patients with concordant ESI-FS showed also interictal concordant ESI of IED; in 3/15 patients, ESI-FS but not ESI-IED was concordant with the resection. In 10/15 patients, ESI-FS was concordant with MRI lesion and for ESI-IED this concordance was on 8/15 patients. Maximal hypometabolism and SISCOM were concordant with the resection for 7/15 and 7/12, respectively.

CONCLUSION

These findings suggest that "non-epileptiform" EEG activity, such as focal slowing, could be a complementary useful marker to localize the epileptogenic zone. ESI-FS may notably be applied in young patients without focal interictal spikes or multifocal spikes. This potential new marker of brain dysfunction has potential applications to other neurological disorders associated with slow EEG activity.

Evidence for the Role of Magnetic Source Imaging in the Presurgical Evaluation of Refractory Epilepsy Patients

Magnetoencephalography (MEG) in the field of epilepsy has multiple advantages; just like electroencephalography (EEG), MEG is able to measure the epilepsy specific information (i.e., the brain activity reflecting seizures and/or interictal epileptiform discharges) directly, non-invasively and with a very high temporal resolution (millisecond-range). In addition MEG has a unique sensitivity for tangential sources, resulting in a full picture of the brain activity when combined with EEG. It accurately allows to perform source imaging of focal epileptic activity and functional cortex and shows a specific high sensitivity for a source in the neocortex. In this paper the current evidence and practice for using magnetic source imaging of focal interictal and ictal epileptic activity during the presurgical evaluation of drug resistant patients is being reviewed.

The effect of stimulation type, head modeling, and combined EEG and MEG on the source reconstruction of the somatosensory P20/N20 component

Modeling and experimental parameters influence the Electro- (EEG) and Magnetoencephalography (MEG) source analysis of the somatosensory P20/N20 component. In a sensitivity group study, we compare P20/N20 source analysis due to different stimulation type (Electric-Wrist [EW], Braille-Tactile [BT], or Pneumato-Tactile [PT]), measurement modality (combined EEG/MEG - EMEG, EEG, or MEG) and head model (standard or individually skull-conductivity calibrated including brain anisotropic conductivity). Considerable differences between pairs of stimulation types occurred (EW-BT: 8.7 ± 3.3 mm/27.1° ± 16.4°, BT-PT: 9 ± 5 mm/29.9° ± 17.3°, and EW-PT: 9.8 ± 7.4 mm/15.9° ± 16.5° and 75% strength reduction of BT or PT when compared to EW) regardless of the head model used. EMEG has nearly no localization differences to MEG, but large ones to EEG (16.1 ± 4.9 mm), while source orientation differences are non-negligible to both EEG (14° ± 3.7°) and MEG (12.5° ± 10.9°). Our calibration results show a considerable inter-subject variability (3.1-14 mS/m) for skull conductivity. The comparison due to different head model show localization differences smaller for EMEG (EW: 3.4 ± 2.4 mm, BT: 3.7 ± 3.4 mm, and PT: 5.9 ± 6.8 mm) than for EEG (EW: 8.6 ± 8.3 mm, BT: 11.8 ± 6.2 mm, and PT: 10.5 ± 5.3 mm), while source orientation differences for EMEG (EW: 15.4° ± 6.3°, BT: 25.7° ± 15.2° and PT: 14° ± 11.5°) and EEG (EW: 14.6° ± 9.5°, BT: 16.3° ± 11.1° and PT: 12.9° ± 8.9°) are in the same range. Our results show that stimulation type, modality and head modeling all have a non-negligible influence on the source reconstruction of the P20/N20 component. The complementary information of both modalities in EMEG can be exploited on the basis of detailed and individualized head models.

Parallel Transmission for Ultrahigh Field MRI

Magnetic resonance imaging (MRI) has been driven toward ultrahigh magnetic fields (UHF) in order to benefit from correspondingly higher signal-to-noise ratio and spectral resolution. Technological challenges associated with UHF, such as increased radiofrequency (RF) energy deposition and RF excitation inhomogeneity, limit realization of the full potential of these benefits. Parallel RF transmission (pTx) enables decreases in the inhomogeneity of RF excitations and in RF energy deposition by using multiple-transmit RF coils driven independently and operating simultaneously. pTx plays a fundamental role in UHF MRI by bringing the potential applications of UHF into reality. In this review article, we review the recent developments in pTx pulse design and RF safety in pTx. Simultaneous multislice imaging and inner volume imaging using pTx are reviewed with a focus on UHF applications. Emerging pTx design approaches using improved pTx design frameworks and calibrations are reviewed together with calibration-free approaches that remove the necessity of time-consuming calibrations necessary for successful pTx. Lastly, we focus on the safety of pTx that is improved by using intersubject variability analysis, proactively managing pTx and temperature-based pTx approaches.

Variation in Reported Human Head Tissue Electrical Conductivity Values

Electromagnetic source characterisation requires accurate volume conductor models representing head geometry and the electrical conductivity field. Head tissue conductivity is often assumed from previous literature, however, despite extensive research, measurements are inconsistent. A meta-analysis of reported human head electrical conductivity values was therefore conducted to determine significant variation and subsequent influential factors. Of 3121 identified publications spanning three databases, 56 papers were included in data extraction. Conductivity values were categorised according to tissue type, and recorded alongside methodology, measurement condition, current frequency, tissue temperature, participant pathology and age. We found variation in electrical conductivity of the whole-skull, the spongiform layer of the skull, isotropic, perpendicularly- and parallelly-oriented white matter (WM) and the brain-to-skull-conductivity ratio (BSCR) could be significantly attributed to a combination of differences in methodology and demographics. This large variation should be acknowledged, and care should be taken when creating volume conductor models, ideally constructing them on an individual basis, rather than assuming them from the literature. When personalised models are unavailable, it is suggested weighted average means from the current meta-analysis are used. Assigning conductivity as: 0.41 S/m for the scalp, 0.02 S/m for the whole skull, or when better modelled as a three-layer skull 0.048 S/m for the spongiform layer, 0.007 S/m for the inner compact and 0.005 S/m for the outer compact, as well as 1.71 S/m for the CSF, 0.47 S/m for the grey matter, 0.22 S/m for WM and 50.4 for the BSCR.

Improved EEG source localization with Bayesian uncertainty modelling of unknown skull conductivity

Electroencephalography (EEG) source imaging is an ill-posed inverse problem that requires accurate conductivity modelling of the head tissues, especially the skull. Unfortunately, the conductivity values are difficult to determine in vivo. In this paper, we show that the exact knowledge of the skull conductivity is not always necessary when the Bayesian approximation error (BAE) approach is exploited. In BAE, we first postulate a probability distribution for the skull conductivity that describes our (lack of) knowledge on its value, and model the effects of this uncertainty on EEG recordings with the help of an additive error term in the observation model. Before the Bayesian inference, the likelihood is marginalized over this error term. Thus, in the inversion we estimate only our primary unknown, the source distribution. We quantified the improvements in the source localization when the proposed Bayesian modelling was used in the presence of different skull conductivity errors and levels of measurement noise. Based on the results, BAE was able to improve the source localization accuracy, particularly when the unknown (true) skull conductivity was much lower than the expected standard conductivity value. The source locations that gained the highest improvements were shallow and originally exhibited the largest localization errors. In our case study, the benefits of BAE became negligible when the signal-to-noise ratio dropped to 20 dB.

A realistic, accurate and fast source modeling approach for the EEG forward problem

The aim of this paper is to advance electroencephalography (EEG) source analysis using finite element method (FEM) head volume conductor models that go beyond the standard three compartment (skin, skull, brain) approach and take brain tissue inhomogeneity (gray and white matter and cerebrospinal fluid) into account. The new approach should enable accurate EEG forward modeling in the thin human cortical structures and, more specifically, in the especially thin cortices in children brain research or in pathological applications. The source model should thus be focal enough to be usable in the thin cortices, but should on the other side be more realistic than the current standard mathematical point dipole. Furthermore, it should be numerically accurate and computationally fast. We propose to achieve the best balance between these demands with a current preserving (divergence conforming) dipolar source model. We develop and investigate a varying number of current preserving source basis elements n (n=1,…,n=5). For validation, we conducted numerical experiments within a multi-layered spherical domain, where an analytical solution exists. We show that the accuracy increases along with the number of basis elements, while focality decreases. The results suggest that the best balance between accuracy and focality in thin cortices is achieved with n=4 (or in extreme cases even n=3) basis functions, while in thicker cortices n=5 is recommended to obtain the highest accuracy. We also compare the current preserving approach to two further FEM source modeling techniques, namely partial integration and St. Venant, and show that the best current preserving source model outperforms the competing methods with regard to overall balance. For all tested approaches, FEM transfer matrices enable high computational speed. We implemented the new EEG forward modeling approaches into the open source duneuro library for forward modeling in bioelectromagnetism to enable its broader use by the brain research community. This library is build upon the DUNE framework for parallel finite elements simulations and integrates with high-level toolboxes like FieldTrip. Additionally, an inversion test has been implemented using the realistic head model to demonstrate and compare the differences between the aforementioned source models.

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.

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.

The Discontinuous Galerkin Finite Element Method for Solving the MEG and the Combined MEG/EEG Forward Problem

In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches in EEG or combined MEG/EEG source analysis scenarios, which motivates a further evaluation of DG-FEM for applications in bioelectromagnetism.

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 (g2) 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.