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

The impact of EEG electrode density on the mapping of cortical activity networks in infants

OBJECTIVE

Electroencephalography (EEG) is widely used for assessing infant's brain activity, and multi-channel recordings support studies on functional cortical networks. Here, we aimed to assess how the number of recording electrodes affects the quality and level of details accessible in studying infant's cortical networks.

METHODS

Dense array EEG recordings with 124 channels from N=20 infants were used as the reference, and lower electrode numbers were subsampled to simulate recording setups with 63, 31, and 19 electrodes, respectively. Cortical activity networks were computed for each recording setup and different frequencies using amplitude and phase correlation measures. The effects of the recording setup were systematically assessed on global, nodal, and edge levels.

RESULTS

Compared to the reference 124-channel recording setup, lowering electrode density affected network measures in a modality- and frequency-specific manner. The global network features were essentially comparable with 63 or 31 channels. However, the analytic reliability of the local network measures, both at nodal and edge levels, was proportional to the electrode density. The low-frequency amplitude correlations were most robust to the number of recording electrodes, whereas higher frequency phase correlation networks were most sensitive to the density of recording electrodes.

CONCLUSIONS

Our findings suggest strong and predictable effects of recording setup on the network analyses. Higher electrode number supports studies on networks with phase correlations, higher frequency, and finer spatial details.

SIGNIFICANCE

The relationship between the recording setup and reliability of network analyses is essential for the prospective design of research data collection, as well as for guiding analytic strategies when using already collected EEG data from infants.

Improving EEG Forward Modeling Using High-Resolution Five-Layer BEM-FMM Head Models: Effect on Source Reconstruction Accuracy

Electroencephalographic (EEG) source localization is a fundamental tool for clinical diagnoses and brain-computer interfaces. We investigate the impact of model complexity on reconstruction accuracy by comparing the widely used three-layer boundary element method (BEM) as an inverse method against a five-layer BEM accelerated by the fast multipole method (BEM-FMM) and coupled with adaptive mesh refinement (AMR) as forward solver. Modern BEM-FMM with AMR can solve high-resolution multi-tissue models efficiently and accurately. We generated noiseless 256-channel EEG data from 15 subjects in the Connectome Young Adult dataset, using four anatomically relevant dipole positions, three conductivity sets, and two head segmentations; we mapped localization errors across the entire grey matter from 4000 dipole positions. The average location error among our four selected dipoles is ∼5mm (±2mm) with an orientation error of ∼12∘ (±7∘). The average source localization error across the entire grey matter is ∼9mm (±4mm), with a tendency for smaller errors on the occipital lobe. Our findings indicate that while three-layer models are robust under noiseless conditions, substantial localization errors (10-20mm) are common. Therefore, models of five or more layers may be needed for accurate source reconstruction in critical applications involving noisy EEG data.

Global sensitivity of EEG source analysis to tissue conductivity uncertainties

Introduction

To reliably solve the EEG inverse problem, accurate EEG forward solutions based on a detailed, individual volume conductor model of the head are essential. A crucial-but often neglected-aspect in generating a volume conductor model is the choice of the tissue conductivities, as these may vary from subject to subject. In this study, we investigate the sensitivity of EEG forward and inverse solutions to tissue conductivity uncertainties for sources distributed over the whole cortex surface.

Methods

We employ a detailed five-compartment head model distinguishing skin, skull, cerebrospinal fluid, gray matter, and white matter, where we consider uncertainties of skin, skull, gray matter, and white matter conductivities. We use the finite element method (FEM) to calculate EEG forward solutions and goal function scans (GFS) as inverse approach. To be able to generate the large number of EEG forward solutions, we employ generalized polynomial chaos (gPC) expansions.

Results

For sources up to a depth of 4 cm, we find the strongest influence on the signal topography of EEG forward solutions for the skull conductivity and a notable effect for the skin conductivity. For even deeper sources, e.g., located deep in the longitudinal fissure, we find an increasing influence of the white matter conductivity. The conductivity variations translate to varying source localizations particularly for quasi-tangential sources on sulcal walls, whereas source localizations of quasi-radial sources on the top of gyri are less affected. We find a strong correlation between skull conductivity and the variation of source localizations and especially the depth of the reconstructed source for quasi-tangential sources. We furthermore find a clear but weaker correlation between depth of the reconstructed source and the skin conductivity.

Discussion

Our results clearly show the influence of tissue conductivity uncertainties on EEG source analysis. We find a particularly strong influence of skull and skin conductivity uncertainties.

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.

Optimized high-definition tDCS in patients with skull defects and skull plates

Introduction

Transcranial direct current stimulation (tDCS) has been shown to benefit patients with brain lesions or traumatic brain injury (TBI). These patients usually have skull defects with different sizes and electrical conductivities. There is very little data in the literature that show how to optimally stimulate these patients with the presence of skull defects.

Methods

Here we leveraged high-resolution (1 mm) realistic head models to explore the best montages targeting right beneath the skull defects with different sizes and conductivities. Specifically, open-source software ROAST was used to solve for the lead field on the publicly available MIDA model. Four different skull defects/plates were modeled with the center above the right primary motor cortex: a larger defect (10 cm diameter) modeled as either titanium or acrylic plate, and a smaller defect (2.5 cm diameter) modeled as either acute state filled with cerebrospinal fluid (CSF) or chronic state with scar tissue. Optimized stimulation with maximal intensity was run using ROAST targeting the right primary motor cortex.

Results

We show that optimized high-definition montages can achieve an average of 0.3 V/m higher stimulation intensities at the target compared to un-optimized montages (M1-SO or 4×1). Large skull defects with titanium or acrylic plates significantly reduce the stimulation intensity by about 80%, while small defects with acute (CSF) or chronic (scar) tissues significantly increase the stimulation intensity by about 200%. Furthermore, one can use M1-SO to achieve almost the same stimulation strength as the optimized montage if the skull has a large defect with titanium plate, and there is no significant difference in stimulation intensity between 4×1 montage and the optimized montage for small skull defects with scar tissue.

Discussion

Based on this work, future modeling studies leveraging individual anatomy of skull defects may help guide tDCS practice on patients with skull defects and skull plates.

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.

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.