Influence of Head Tissue Conductivity Uncertainties on EEG Dipole Reconstruction

Exploring Electrocortical Signatures of Gait Adaptation: Differential Neural Dynamics in Slow and Fast Gait Adapters

Individuals exhibit significant variability in their ability to adapt locomotor skills, with some adapting quickly and others more slowly. Differences in brain activity likely contribute to this variability, but direct neural evidence is lacking. We investigated individual differences in electrocortical activity that led to faster locomotor adaptation rates. We recorded high-density electroencephalography while young, neurotypical adults adapted their walking on a split-belt treadmill and grouped them based on how quickly they restored their gait symmetry. Results revealed unique spectral signatures within the posterior parietal, bilateral sensorimotor, and right visual cortices that differ between fast and slow adapters. Specifically, fast adapters exhibited lower alpha power in the posterior parietal and right visual cortices during early adaptation, associated with quicker attainment of steady-state step length symmetry. Decreased posterior parietal alpha may reflect enhanced spatial attention, sensory integration, and movement planning to facilitate faster locomotor adaptation. Conversely, slow adapters displayed greater alpha and beta power in the right visual cortex during late adaptation, suggesting potential differences in visuospatial processing. Additionally, fast adapters demonstrated reduced spectral power in the bilateral sensorimotor cortices compared with slow adapters, particularly in the theta band, which may suggest variations in perception of the split-belt perturbation. These findings suggest that alpha and beta oscillations in the posterior parietal and visual cortices and theta oscillations in the sensorimotor cortex are related to the rate of gait adaptation.

Electrocortical theta activity may reflect sensory prediction errors during adaptation to a gradual gait perturbation

Locomotor adaptation to abrupt and gradual perturbations are likely driven by fundamentally different neural processes. The aim of this study was to quantify brain dynamics associated with gait adaptation to a gradually introduced gait perturbation, which typically results in smaller behavioral errors relative to an abrupt perturbation. Loss of balance during standing and walking elicits transient increases in midfrontal theta oscillations that have been shown to scale with perturbation intensity. We hypothesized there would be no significant change in anterior cingulate theta power (4-7 Hz) with respect to pre-adaptation when a gait perturbation is introduced gradually because the gradual perturbation acceleration and stepping kinematic errors are small relative to an abrupt perturbation. Using mobile electroencephalography (EEG), we measured gait-related spectral changes near the anterior cingulate, posterior cingulate, sensorimotor, and posterior parietal cortices as young, neurotypical adults (n = 30) adapted their gait to an incremental split-belt treadmill perturbation. Most cortical clusters we examined (>70%) did not exhibit changes in electrocortical activity between 2-50 Hz. However, we did observe gait-related theta synchronization near the left anterior cingulate cortex during strides with the largest errors, as measured by step length asymmetry. These results suggest gradual adaptation with small gait asymmetry and perturbation magnitude may not require significant cortical resources beyond normal treadmill walking. Nevertheless, the anterior cingulate may remain actively engaged in error monitoring, transmitting sensory prediction error information via theta oscillations.

Accuracy of dipole source reconstruction in the 3-layer BEM model against the 5-layer BEM-FMM model

Objective

To compare cortical dipole fitting spatial accuracy between the widely used yet highly simplified 3-layer and modern more realistic 5-layer BEM-FMM models with and without adaptive mesh refinement (AMR) methods.

Methods

We generate simulated noiseless 256-channel EEG data from 5-layer (7-compartment) meshes of 15 subjects from the Connectome Young Adult dataset. For each subject, we test four dipole positions, three sets of conductivity values, and two types of head segmentation. We use the boundary element method (BEM) with fast multipole method (FMM) acceleration, with or without (AMR), for forward modeling. Dipole fitting is carried out with the FieldTrip MATLAB toolbox.

Results

The average position error (across all tested dipoles, subjects, and models) is ~4 mm, with a standard deviation of ~2 mm. The orientation error is ~20° on average, with a standard deviation of ~15°. Without AMR, the numerical inaccuracies produce a larger disagreement between the 3- and 5-layer models, with an average position error of ~8 mm (6 mm standard deviation), and an orientation error of 28° (28° standard deviation).

Conclusions

The low-resolution 3-layer models provide excellent accuracy in dipole localization. On the other hand, dipole orientation is retrieved less accurately. Therefore, certain applications may require more realistic models for practical source reconstruction. AMR is a critical component for improving the accuracy of forward EEG computations using a high-resolution 5-layer volume conduction model.

Significance

Improving EEG source reconstruction accuracy is important for several clinical applications, including epilepsy and other seizure-inducing conditions.

A hybrid boundary element-finite element approach for solving the EEG forward problem in brain modeling

This article introduces a hybrid BE-FE method for solving the EEG forward problem, leveraging the strengths of both the Boundary Element Method (BEM) and Finite Element Method (FEM). FEM accurately models complex and anisotropic tissue properties for realistic head geometries, while BEM excels in handling isotropic tissue regions and dipolar sources efficiently. The proposed hybrid method divides regions into homogeneous boundary element (BE) regions that include sources and heterogeneous anisotropic finite element (FE) regions. So, BEM models the brain, including dipole sources, and FEM models other head layers. Validation includes inhomogeneous isotropic/anisotropic three- and four-layer spherical head models, and a four-layer MRI-based realistic head model. Results for six dipole eccentricities and two orientations are computed using BEM, FEM, and hybrid BE-FE method. Statistical analysis, comparing error criteria of RDM and MAG, reveals notable improvements using the hybrid FE-BE method. In the spherical head model, the hybrid BE-FE method compared with FEM demonstrates enhancements of at least 1.05 and 38.31% in RDM and MAG criteria, respectively. Notably, in the anisotropic four-layer head model, improvements reach a maximum of 88.3% for RDM and 93.27% for MAG over FEM. Moreover, in the anisotropic four-layer realistic head model, the proposed hybrid method exhibits 55.4% improvement in RDM and 89.3% improvement in MAG compared to FEM. These findings underscore the proposed method is a promising approach for solving the realistic EEG forward problems, advancing neuroimaging techniques and enhancing understanding of brain function.

Comparing the performance of beamformer algorithms in estimating orientations of neural sources

The efficacy of transcranial electric stimulation (tES) to effectively modulate neuronal activity depends critically on the spatial orientation of the targeted neuronal population. Therefore, precise estimation of target orientation is of utmost importance. Different beamforming algorithms provide orientation estimates; however, a systematic analysis of their performance is still lacking. For fixed brain locations, EEG and MEG data from sources with randomized orientations were simulated. The orientation was then estimated (1) with an EEG and (2) with a combined EEG-MEG approach. Three commonly used beamformer algorithms were evaluated with respect to their abilities to estimate the correct orientation: Unit-Gain (UG), Unit-Noise-Gain (UNG), and Array-Gain (AG) beamformer. Performance depends on the signal-to-noise ratios for the modalities and on the chosen beamformer. Overall, the UNG and AG beamformers appear as the most reliable. With increasing noise, the UG estimate converges to a vector determined by the leadfield, thus leading to insufficient orientation estimates.

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.

Neurobiologically realistic neural network enables cross-scale modeling of neural dynamics

Fundamental principles underlying computation in multi-scale brain networks illustrate how multiple brain areas and their coordinated activity give rise to complex cognitive functions. Whereas brain activity has been studied at the micro- to meso-scale to reveal the connections between the dynamical patterns and the behaviors, investigations of neural population dynamics are mainly limited to single-scale analysis. Our goal is to develop a cross-scale dynamical model for the collective activity of neuronal populations. Here we introduce a bio-inspired deep learning approach, termed NeuroBondGraph Network (NBGNet), to capture cross-scale dynamics that can infer and map the neural data from multiple scales. Our model not only exhibits more than an 11-fold improvement in reconstruction accuracy, but also predicts synchronous neural activity and preserves correlated low-dimensional latent dynamics. We also show that the NBGNet robustly predicts held-out data across a long time scale (2 weeks) without retraining. We further validate the effective connectivity defined from our model by demonstrating that neural connectivity during motor behaviour agrees with the established neuroanatomical hierarchy of motor control in the literature. The NBGNet approach opens the door to revealing a comprehensive understanding of brain computation, where network mechanisms of multi-scale activity are critical.

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.

Dual-layer electroencephalography data during real-world table tennis

Real-world settings are necessary to improve the ecological validity of neuroscience research, and electroencephalography (EEG) facilitates mobile electrocortical recordings because of its easy portability and high temporal resolution. Table tennis is a whole-body, goal-directed sport that requires constant visuomotor feedback, anticipation, strategic decision-making, object interception, and performance monitoring - making it an interesting testbed for a variety of neuroscience studies. Although traditionally plagued by artifact contamination, recent advances in signal processing and hardware approaches, such as the dual-layer approach, have allowed high fidelity EEG recordings during whole-body maneuvers. Here, we present a dual-layer EEG dataset from 25 healthy human participants playing table tennis with a human opponent and a ball machine. Our dataset includes synchronized, multivariate time series recordings from 120 scalp electrodes, 120 noise electrodes, 8 neck electromyography electrodes, and inertial measurement units on the participant, paddles, and ball machine to record hit events. We also include de-identified T1 anatomical MR images and digitized electrode locations to create subject-specific head models for source localization. In addition, we provide anonymized video recordings and Adobe Premiere project files with hit events labeled (originally used to mark successful/missed hits). Researchers could use the videos to mark their own events of interest. We formatted our dataset in the Brain Imaging Data Structure (BIDS) format to facilitate data reuse and to adhere to the scientific community's new organization standard.

Diagnostic deep learning algorithms that use resting EEG to distinguish major depressive disorder, bipolar disorder, and schizophrenia from each other and from healthy volunteers

BACKGROUND

Mood disorders and schizophrenia affect millions worldwide. Currently, diagnosis is primarily determined by reported symptomatology. As symptoms may overlap, misdiagnosis is common, potentially leading to ineffective or destabilizing treatment. Diagnostic biomarkers could significantly improve clinical care by reducing dependence on symptomatic presentation.

METHODS

We used deep learning analysis (DLA) of resting electroencephalograph (EEG) to differentiate healthy control (HC) subjects (N = 239), from those with major depressive disorder (MDD) (N = 105), MDD-atypical (MDD-A) (N = 27), MDD-psychotic (MDD-P) (N = 35), bipolar disorder-depressed episode (BD-DE) (N = 71), BD-manic episode (BD-ME) (N = 49), and schizophrenia (SCZ) (N = 122) and also differentiate subjects with mental disorders on a pair-wise basis. DSM-III-R diagnoses were determined and supplemented by computerized Quick Diagnostic Interview Schedule. After EEG preprocessing, robust exact low-resolution electromagnetic tomography (ReLORETA) computed EEG sources for 82 brain regions. 20 % of all subjects were then set aside for independent testing. Feature selection methods were then used for the remaining subjects to identify brain source regions that are discriminating between diagnostic categories.

RESULTS

Pair-wise classification accuracies between 90 % and 100 % were obtained using independent test subjects whose data were not used for training purposes. The most frequently selected features across various pairs are in the postcentral, supramarginal, and fusiform gyri, the hypothalamus, and the left cuneus. Brain sites discriminating SCZ from HC were mainly in the left hemisphere while those separating BD-ME from HC were on the right.

LIMITATIONS

The use of superseded DSM-III-R diagnostic system and relatively small sample size in some disorder categories that may increase the risk of overestimation.

CONCLUSIONS

DLA of EEG could be trained to autonomously classify psychiatric disorders with over 90 % accuracy compared to an expert clinical team using standardized operational methods.

A table tennis serve versus rally hit elicits differential hemispheric electrocortical power fluctuations

Human visuomotor control requires coordinated interhemispheric interactions to exploit the brain's functional lateralization. In right-handed individuals, the left hemisphere (right arm) is better for dynamic control and the right hemisphere (left arm) is better for impedance control. Table tennis is a game that requires precise movements of the paddle, whole body coordination, and cognitive engagement, providing an ecologically valid way to study visuomotor integration. The sport has many different types of strokes (e.g., serve, return, and rally shots), which should provide unique cortical dynamics given differences in the sensorimotor demands. The goal of this study was to determine the hemispheric specialization of table tennis serving - a sequential, self-paced, bimanual maneuver. We used time-frequency analysis, event-related potentials, and functional connectivity measures of source-localized electrocortical clusters and compared serves with other types of shots, which varied in the types of movement required, attentional focus, and other task demands. We found greater alpha (8-12 Hz) and beta (13-30 Hz) power in the right sensorimotor cortex than in the left sensorimotor cortex, and we found a greater magnitude of spectral power fluctuations in the right sensorimotor cortex for serve hits than return or rally hits, in all right-handed participants. Surprisingly, we did not find a difference in interhemispheric functional connectivity between a table tennis serve and return or rally hits, even though a serve could arguably be a more complex maneuver. Studying real-world brain dynamics of table tennis provides insight into bilateral sensorimotor integration.NEW & NOTEWORTHY We found different spectral power fluctuations in the left and right sensorimotor cortices during table tennis serves, returns, and rallies. Our findings contribute to the basic science understanding of hemispheric specialization in a real-world context.

Electrocortical activity correlated with locomotor adaptation during split-belt treadmill walking

Locomotor adaptation is crucial for daily gait adjustments to changing environmental demands and obstacle avoidance. Mobile brain imaging with high-density electroencephalography (EEG) now permits quantification of electrocortical dynamics during human locomotion. To determine the brain areas involved in human locomotor adaptation, we recorded high-density EEG from healthy, young adults during split-belt treadmill walking. We incorporated a dual-electrode EEG system and neck electromyography to decrease motion and muscle artefacts. Voluntary movement preparation and execution have been linked to alpha (8-13 Hz) and beta band (13-30 Hz) desynchronizations in the sensorimotor and posterior parietal cortices, whereas theta band (4-7 Hz) modulations in the anterior cingulate have been correlated with movement error monitoring. We hypothesized that relative to normal walking, split-belt walking would elicit: (1) decreases in alpha and beta band power in sensorimotor and posterior parietal cortices, reflecting enhanced motor flexibility; and (2) increases in theta band power in anterior cingulate cortex, reflecting instability and balance errors that will diminish with practice. We found electrocortical activity in multiple regions that was associated with stages of gait adaptation. Data indicated that sensorimotor and posterior parietal cortices had decreased alpha and beta band spectral power during early adaptation to split-belt treadmill walking that gradually returned to pre-adaptation levels by the end of the adaptation period. Our findings emphasize that multiple brain areas are involved in adjusting gait under changing environmental demands during human walking. Future studies could use these findings on healthy, young participants to identify dysfunctional supraspinal mechanisms that may be impairing gait adaptation. KEY POINTS: Identifying the location and time course of electrical changes in the brain correlating with gait adaptation increases our understanding of brain function and provides targets for brain stimulation interventions. Using high-density EEG in combination with 3D biomechanics, we found changes in neural oscillations localized near the sensorimotor, posterior parietal and cingulate cortices during split-belt treadmill adaptation. These findings suggest that multiple cortical mechanisms may be associated with locomotor adaptation, and their temporal dynamics can be quantified using mobile EEG. Results from this study can serve as a reference model to examine brain dynamics in individuals with movement disorders that cause gait asymmetry and reduced gait adaptation.

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.

Exploring the Neural Processes behind Narrative Engagement: An EEG Study

Past cognitive neuroscience studies using naturalistic stimuli have considered narratives holistically and focused on cognitive processes. In this study, we incorporated the narrative structure, the dramatic arc, as an object of investigation, to examine how engagement levels fluctuate across a narrative-aligned dramatic arc. We explored the possibility of predicting self-reported engagement ratings from neural activity and investigated the idiosyncratic effects of each phase of the dramatic arc on brain responses as well as the relationship between engagement and brain responses. We presented a movie excerpt following the six-phase narrative arc structure to female and male participants while collecting EEG signals. We then asked this group of participants to recall the excerpt, another group to segment the video based on the dramatic arc model, and a third to rate their engagement levels while watching the movie. The results showed that the self-reported engagement ratings followed the pattern of the narrative dramatic arc. Moreover, while EEG amplitude could not predict group-averaged engagement ratings, other features comprising dynamic intersubject correlation (dISC), including certain frequency bands, dynamic functional connectivity patterns and graph features were able to achieve this. Furthermore, neural activity in the last two phases of the dramatic arc significantly predicted engagement patterns. This study is the first to explore the cognitive processes behind the dramatic arc and its phases. By demonstrating how neural activity predicts self-reported engagement, which itself aligns with the narrative structure, this study provides insights on the interrelationships between narrative structure, neural responses, and viewer engagement.

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.

Comparison of EEG Source Localization Using Simplified and Anatomically Accurate Head Models in Younger and Older Adults

Accuracy of electroencephalography (EEG) source localization relies on the volume conduction head model. A previous analysis of young adults has shown that simplified head models have larger source localization errors when compared with head models based on magnetic resonance images (MRIs). As obtaining individual MRIs may not always be feasible, researchers often use generic head models based on template MRIs. It is unclear how much error would be introduced using template MRI head models in older adults that likely have differences in brain structure compared to young adults. The primary goal of this study was to determine the error caused by using simplified head models without individual-specific MRIs in both younger and older adults. We collected high-density EEG during uneven terrain walking and motor imagery for 15 younger (22±3 years) and 21 older adults (74±5 years) and obtained [Formula: see text]-weighted MRI for each individual. We performed equivalent dipole fitting after independent component analysis to obtain brain source locations using four forward modeling pipelines with increasing complexity. These pipelines included: 1) a generic head model with template electrode positions or 2) digitized electrode positions, 3) individual-specific head models with digitized electrode positions using simplified tissue segmentation, or 4) anatomically accurate segmentation. We found that when compared to the anatomically accurate individual-specific head models, performing dipole fitting with generic head models led to similar source localization discrepancies (up to 2 cm) for younger and older adults. Co-registering digitized electrode locations to the generic head models reduced source localization discrepancies by  ∼  6 mm. Additionally, we found that source depths generally increased with skull conductivity for the representative young adult but not as much for the older adult. Our results can help inform a more accurate interpretation of brain areas in EEG studies when individual MRIs are unavailable.

Optimized Method of Unilateral Spinal Anesthesia: A Double-blind, Randomized Clinical Study

Background

Unilateral spinal anesthesia is often accompanied by technical difficulties in implementation, multiple puncture attempts, unsuccessful punctures, and, as a result, insufficient anesthesia, along with various complications associated with a dural puncture.

Objectives

This work compares the efficacy and safety of conventional Unilateral Spinal Anesthesia (USpA) and unilateral spinal anesthesia with electrical nerve stimulation (USpA+ENS).

Methods

A total of 134 patients with an upcoming vascular surgery on one lower extremity were randomly assigned to two groups. All the patients were positioned with the operated limb below and used 7.5 mg of Bupivacaine-Spinal®. In the UsPA group, anesthesia was performed according to the standard technique. In the USpA+ENS group, electrical nerve stimulation was additionally used. Primary outcomes were the presence or absence of post-dural puncture headache (PDPH), number of puncture attempts, lateralization, and anesthesia adequacy. Secondary outcomes were intraoperative pain scores, the presence or absence of nausea and vomiting, and the need for hemodynamics correction.

Results

The frequency of puncture complications was sufficiently lower in the USpA+ENS group than in the UsPA group. The local anesthetic solution distribution, pain score indicators, and secondary outcomes were comparable in both groups with a slight difference.

Conclusions

We showed that USpA+ENS reduces the incidence of puncture complications and improves the quality of anesthesia and adherence of both patients and anesthesiologists to the unilateral spinal anesthesia technique.

In silicostudy of the effects of cerebral circulation on source localization using a dynamical anatomical atlas of the human head

Objective.This study focuses on the effects of dynamical vascular modeling on source localization errors in electroencephalography (EEG). Our aim of thisin silicostudy is to (a) find out the effects of cerebral circulation on the accuracy of EEG source localization estimates, and (b) evaluate its relevance with respect to measurement noise and interpatient variation.Approach.We employ a four-dimensional (3D + T) statistical atlas of the electrical properties of the human head with a cerebral circulation model to generate virtual patients with different cerebral circulatory conditions for EEG source localization analysis. As source reconstruction techniques, we use the linearly constraint minimum variance (LCMV) beamformer, standardized low-resolution brain electromagnetic tomography (sLORETA), and the dipole scan (DS).Main results.Results indicate that arterial blood flow affects source localization at different depths and with varying significance. The average flow rate plays an important role in source localization performance, while the pulsatility effects are very small. In cases where a personalized model of the head is available, blood circulation mismodeling causes localization errors, especially in the deep structures of the brain where the main cerebral arteries are located. When interpatient variations are considered, the results show differences up to 15 mm for sLORETA and LCMV beamformer and 10 mm for DS in the brainstem and entorhinal cortices regions. In regions far from the main arteries vessels, the discrepancies are smaller than 3 mm. When measurement noise is added and interpatient differences are considered in a deep dipolar source, the results indicate that the effects of conductivity mismatch are detectable even for moderate measurement noise. The signal-to-noise ratio limit for sLORETA and LCMV beamformer is 15 dB, while the limit is under 30 dB for DS.Significance.Localization of the brain activity via EEG constitutes an ill-posed inverse problem, where any modeling uncertainty, e.g. a slight amount of noise in the data or material parameter discrepancies, can lead to a significant deviation of the estimated activity, especially in the deep structures of the brain. Proper modeling of the conductivity distribution is necessary in order to obtain an appropriate source localization. In this study, we show that the conductivity of the deep brain structures is particularly impacted by blood flow-induced changes in conductivity because large arteries and veins access the brain through that region.

Neural stochastic differential equations network as uncertainty quantification method for EEG source localization

EEG source localization remains a challenging problem given the uncertain conductivity values of the volume conductor models (VCMs). As uncertain conductivities vary across people, they may considerably impact the forward and inverse solutions of the EEG, leading to an increase in localization mistakes and misdiagnoses of brain disorders. Calibration of conductivity values using uncertainty quantification (UQ) techniques is a promising approach to reduce localization errors. The widely-known UQ methods involve Bayesian approaches, which utilize prior conductivity values to derive their posterior inference and estimate their optimal calibration. However, these approaches have two significant drawbacks: solving for posterior inference is intractable, and choosing inappropriate priors may lead to increased localization mistakes. This study used the Neural Stochastic Differential equations Network (SDE-Net), a combination of dynamical systems and deep learning techniques that utilizes the Wiener process to minimize conductivity uncertainties in the VCM and improve the inverse problem. Results revealed that SDE-Net generated a lower localization error rate in the inverse problem compared to Bayesian techniques. Future studies may employ new stochastic dynamical systems-based techniques as a UQ technique to address further uncertainties in the EEG Source Localization problem. Our code can be found here:https://github.com/rrwabina/SDENet-UQ-ESL.

Brainstorm-DUNEuro: An integrated and user-friendly Finite Element Method for modeling electromagnetic brain activity

Human brain activity generates scalp potentials (electroencephalography - EEG), intracranial potentials (iEEG), and external magnetic fields (magnetoencephalography - MEG). These electrophysiology (e-phys) signals can often be measured simultaneously for research and clinical applications. The forward problem involves modeling these signals at their sensors for a given equivalent current dipole configuration within the brain. While earlier researchers modeled the head as a simple set of isotropic spheres, today's magnetic resonance imaging (MRI) data allow for a detailed anatomic description of brain structures and anisotropic characterization of tissue conductivities. We present a complete pipeline, integrated into the Brainstorm software, that allows users to automatically generate an individual and accurate head model based on the subject's MRI and calculate the electromagnetic forward solution using the finite element method (FEM). The head model generation is performed by integrating the latest tools for MRI segmentation and FEM mesh generation. The final head model comprises the five main compartments: white-matter, gray-matter, CSF, skull, and scalp. The anisotropic brain conductivity model is based on the effective medium approach (EMA), which estimates anisotropic conductivity tensors from diffusion-weighted imaging (DWI) data. The FEM electromagnetic forward solution is obtained through the DUNEuro library, integrated into Brainstorm, and accessible with either a user-friendly graphical interface or scripting. With tutorials and example data sets available in an open-source format on the Brainstorm website, this integrated pipeline provides access to advanced FEM tools for electromagnetic modeling to a broader neuroscience community.

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.

Source-based artifact-rejection techniques for TMS-EEG

Neuronal electroencephalography (EEG) signals arise from the cortical postsynaptic currents. Due to the conductive properties of the head, these neuronal sources produce relatively smeared spatial patterns in EEG. We can model these topographies to deduce which signals reflect genuine TMS-evoked cortical activity and which data components are merely noise and artifacts. This review will concentrate on two source-based artifact-rejection techniques developed for TMS-EEG data analysis, signal-space-projection-source-informed reconstruction (SSP-SIR), and the source-estimate-utilizing noise-discarding algorithm (SOUND). The former method was designed for rejecting TMS-evoked muscle artifacts, while the latter was developed to suppress noise signals from EEG and magnetoencephalography (MEG) in general. We shall cover the theoretical background for both methods, but most importantly, we will describe some essential practical perspectives for using these techniques effectively. We demonstrate and explain what approaches produce the most reliable inverse estimates after cleaning the data or how to perform non-biased comparisons between cleaned datasets. All noise-cleaning algorithms compromise the signals of interest to a degree. We elaborate on how the source-based methods allow objective quantification of the overcorrection. Finally, we consider possible future directions. While this article concentrates on TMS-EEG data analysis, many theoretical and practical aspects, presented here, can be readily applied in other EEG/MEG applications. Overall, the source-based cleaning methods provide a valuable set of TMS-EEG preprocessing tools. We can objectively evaluate their performance regarding possible overcorrection. Furthermore, the overcorrection can always be taken into account to compare cleaned datasets reliably. The described methods are based on current electrophysiological and anatomical understanding of the head and the EEG generators; strong assumptions of the statistical properties of the noise and artifact signals, such as independence, are not needed.

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.

Shamo: A Tool for Electromagnetic Modeling, Simulation and Sensitivity Analysis of the Head

Accurate electromagnetic modeling of the head of a subject is of main interest in the fields of source reconstruction and brain stimulation. Those processes rely heavily on the quality of the model and, even though the geometry of the tissues can be extracted from magnetic resonance images (MRI) or computed tomography (CT), their physical properties such as the electrical conductivity are difficult to measure with non intrusive techniques. In this paper, we propose a tool to assess the uncertainty in the model parameters, the tissue conductivity, as well as compute a parametric forward models for electroencephalography (EEG) and transcranial direct current stimulation (tDCS) current distribution.

Personalized tDCS for Focal Epilepsy—A Narrative Review: A Data-Driven Workflow Based on Imaging and EEG Data

Conventional transcranial electric stimulation(tES) using standard anatomical positions for the electrodes and standard stimulation currents is frequently not sufficiently selective in targeting and reaching specific brain locations, leading to suboptimal application of electric fields. Recent advancements in in vivo electric field characterization may enable clinical researchers to derive better relationships between the electric field strength and the clinical results. Subject-specific electric field simulations could lead to improved electrode placement and more efficient treatments. Through this narrative review, we present a processing workflow to personalize tES for focal epilepsy, for which there is a clear cortical target to stimulate. The workflow utilizes clinical imaging and electroencephalography data and enables us to relate the simulated fields to clinical outcomes. We review and analyze the relevant literature for the processing steps in the workflow, which are the following: tissue segmentation, source localization, and stimulation optimization. In addition, we identify shortcomings and ongoing trends with regard to, for example, segmentation quality and tissue conductivity measurements. The presented processing steps result in personalized tES based on metrics like focality and field strength, which allow for correlation with clinical outcomes.

The interaction between endogenous GABA, functional connectivity, and behavioral flexibility is critically altered with advanced age

The flexible adjustment of ongoing behavior challenges the nervous system’s dynamic control mechanisms and has shown to be specifically susceptible to age-related decline. Previous work links endogenous gamma-aminobutyric acid (GABA) with behavioral efficiency across perceptual and cognitive domains, with potentially the strongest impact on those behaviors that require a high level of dynamic control. Our analysis integrated behavior and modulation of interhemispheric phase-based connectivity during dynamic motor-state transitions with endogenous GABA concentration in adult human volunteers. We provide converging evidence for age-related differences in the preferred state of endogenous GABA concentration for more flexible behavior. We suggest that the increased interhemispheric connectivity observed in the older participants represents a compensatory neural mechanism caused by phase-entrainment in homotopic motor cortices. This mechanism appears to be most relevant in the presence of a less optimal tuning of the inhibitory tone as observed during healthy aging to uphold the required flexibility of behavioral action. Future work needs to validate the relevance of this interplay between neural connectivity and GABAergic inhibition for other domains of flexible human behavior.

An integrated analysis of GABA concentration, interhemispheric phase-based connectivity behavior and in healthy human participants revealed a task-specific association between GABA concentration and motor-cortical connectivity. These results shed light on how GABA is involved in neural oscillations and behavioral performance.

Directionality of the injected current targeting the P20/N20 source determines the efficacy of 140 Hz transcranial alternating current stimulation (tACS)-induced aftereffects in the somatosensory cortex

Interindividual anatomical differences in the human cortex can lead to suboptimal current directions and may result in response variability of transcranial electrical stimulation methods. These differences in brain anatomy require individualized electrode stimulation montages to induce an optimal current density in the targeted area of each individual subject. We aimed to explore the possible modulatory effects of 140 Hz transcranial alternating current stimulation (tACS) on the somatosensory cortex using personalized multi-electrode stimulation montages. In two randomized experiments using either tactile finger or median nerve stimulation, we measured by evoked potentials the plasticity aftereffects and oscillatory power changes after 140 Hz tACS at 1.0 mA as compared to sham stimulation (n = 17, male = 9). We found a decrease in the power of oscillatory mu-rhythms during and immediately after tactile discrimination tasks, indicating an engagement of the somatosensory system during stimulus encoding. On a group level both the oscillatory power and the evoked potential amplitudes were not modulated by tACS neither after tactile finger stimulation nor after median nerve stimulation as compared to sham stimulation. On an individual level we could however demonstrate that lower angular difference (i.e., differences between the injected current vector in the target region and the source orientation vector) is associated with significantly higher changes in both P20/N20 and N30/P30 source activities. Our findings suggest that the higher the directionality of the injected current correlates to the dipole orientation the greater the tACS-induced aftereffects are.

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.

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.

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.

Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI

Imaging of the electrical conductivity distribution inside the human body has been investigated for numerous clinical applications. The conductivity tensors of biological tissue have been obtained from water diffusion tensors by applying several models, which may not cover the entire phenomenon. Recently, a new conductivity tensor imaging (CTI) method was developed through a combination of B1 mapping, and multi-b diffusion weighted imaging. In this study, we compared the most recent CTI method with the four existing models of conductivity tensors reconstruction. Two conductivity phantoms were designed to evaluate the accuracy of the models. Applied to five human brains, the conductivity tensors using the four existing models and CTI were imaged and compared with the values from the literature. The conductivity image of the phantoms by the CTI method showed relative errors between 1.10% and 5.26%. The images by the four models using DTI could not measure the effects of different ion concentrations subsequently due to prior information of the mean conductivity values. The conductivity tensor images obtained from five human brains through the CTI method were comparable to previously reported literature values. The images by the four methods using DTI were highly correlated with the diffusion tensor images, showing a coefficient of determination (R2) value of 0.65 to 1.00. However, the images by the CTI method were less correlated with the diffusion tensor images and exhibited an averaged R2 value of 0.51. The CTI method could handle the effects of different ion concentrations as well as mobilities and extracellular volume fractions by collecting and processing additional B1 map data. It is necessary to select an application-specific model taking into account the pros and cons of each model. Future studies are essential to confirm the usefulness of these conductivity tensor imaging methods in clinical applications, such as tumor characterization, EEG source imaging, and treatment planning for electrical stimulation.

Impact of skull-to-brain conductivity ratio for high resolution EEG source localization

Objective. To measure the impact of skull-to-brain conductivity ratios on interictal spikes source localizations, using high resolution EEG (HR EEG). In previous studies, two ratios were mainly employed: 1/80 and 1/40. Consequences of the employed ratios on source localization results are poorly studied.Methods. Twenty patients with drug-resistant epilepsy were studied using HR EEG (sixty-four scalp electrodes). For each patient, three-layers realistic head models based on individual MRI were elaborated using boundary element model. For each interictal spike, source localization was performed six times, using six skull-to-brain conductivity ratios (1/80, 1/50, 1/40, 1/30, 1/20 and 1/10), exploring all the spectrum of values reported in the literature. We then measured distances between the different sources obtained and between the sources and the anterior commissure (in order to estimate sources depth).Results. We measured a mean distance of 5.3 mm (sd: 3 mm) between the sources obtained with 1/40 versus 1/80 ratio. This distance increased when the discrepancy between the two evaluated ratios increased. We measured a mean distance of 14.2 mm (sd: 4.9 mm) between sources obtained with 1/10 ratio versus 1/80 ratio. Sources localized using 1/40 ratio were 4.3 mm closer to the anterior commissure than sources localized using 1/80 ratio.Significance. Skull-to-brain conductivity ratio is an often-neglected parameter in source localization studies. The different ratios mainly used in the litterature (1/80 and 1/40) lead to significant differences in source localizations. These variations mainly occur in source depth. A more accurate estimation of skull-to-brain conductivity is needed to increase source localization accuracy.Abbreviations. ECD: equivalent current dipole; EIT: electric impedance tomography, HR EEG: High resolution Electroencephalography, IIS: Inter ictal spikes, MEG: Magnetoencephalography, MRI: Magnetic resonance imaging, mS/m: milli-Siemens/m, S/m: Siemens/m, SD: Standard deviation.

Theta but not beta power is positively associated with better explicit motor task learning

Neurophysiologic correlates of motor learning that can be monitored during neurorehabilitation interventions can facilitate the development of more effective learning methods. Previous studies have focused on the role of the beta band (14-30 Hz) because of its clear response during motor activity. However, it is difficult to discriminate between beta activity related to learning a movement and performing the movement. In this study, we analysed differences in the electroencephalography (EEG) power spectra of complex and simple explicit sequential motor tasks in healthy young subjects. The complex motor task (CMT) allowed EEG measurement related to motor learning. In contrast, the simple motor task (SMT) made it possible to control for EEG activity associated with performing the movement without significant motor learning. Source reconstruction of the EEG revealed task-related activity from 5 clusters covering both primary motor cortices (M1) and 3 clusters localised to different parts of the cingulate cortex (CC). We found no association between M1 beta power and learning, but the CMT produced stronger bilateral beta suppression compared to the SMT. However, there was a positive association between contralateral M1 theta (5-8 Hz) and alpha (8-12 Hz) power and motor learning, and theta and alpha power in the posterior mid-CC and posterior CC were positively associated with greater motor learning. These findings suggest that the theta and alpha bands are more related to motor learning than the beta band, which might merely relate to the level of perceived difficulty during learning.

Does participant's age impact on tDCS induced fields? Insights from computational simulations

Objective: Understanding the induced current flow from transcranial direct current stimulation (tDCS) is essential for determining the optimal dose and treatment. Head tissue conductivities play a key role in the resulting electromagnetic fields. However, there exists a complicated relationship between skull conductivity and participant age, that remains unclear. We explored how variations in skull electrical conductivities, particularly as a suggested function of age, affected tDCS induced electric fields.Approach: Simulations were employed to compare tDCS outcomes for different intensities across head atlases of varying age. Three databases were chosen to demonstrate differing variability in skull conductivity with age and how this may affect induced fields. Differences in tDCS electric fields due to proposed age-dependent skull conductivity variation, as well as deviations in grey matter, white matter and scalp, were compared and the most influential tissues determined.Main results: tDCS induced peak electric fields significantly negatively correlated with age, exacerbated by employing proposed age-appropriate skull conductivity (according to all three datasets). Uncertainty in skull conductivity was the most sensitive to changes in peak fields with increasing age. These results were revealed to be directly due to changing skull conductivity, rather than head geometry alone. There was no correlation between tDCS focality and age.Significance: Accurate and individualised head anatomy andin vivoskull conductivity measurements are essential for modelling tDCS induced fields. In particular, age should be taken into account when considering stimulation dose to precisely predict outcomes.

Age-Related EEG Power Reductions Cannot Be Explained by Changes of the Conductivity Distribution in the Head Due to Brain Atrophy

Electroencephalogram (EEG) power reductions in the aging brain have been described by numerous previous studies. However, the underlying mechanism for the observed brain signal power reduction remains unclear. One possible cause for reduced EEG signals in elderly subjects might be the increased distance from the primary neural electrical currents on the cortex to the scalp electrodes as the result of cortical atrophies. While brain shrinkage itself reflects age-related neurological changes, the effects of changes in the distribution of electrical conductivity are often not distinguished from altered neural activity when interpreting EEG power reductions. To address this ambiguity, we employed EEG forward models to investigate whether brain shrinkage is a major factor for the signal attenuation in the aging brain. We simulated brain shrinkage in spherical and realistic brain models and found that changes in the conductor geometry cannot fully account for the EEG power reductions even when the brain was shrunk to unrealistic sizes. Our results quantify the extent of power reductions from brain shrinkage and pave the way for more accurate inferences about deficient neural activity and circuit integrity based on EEG power reductions in the aging population.

Preclinical Neuropathic Pain Assessment; the Importance of Translatability and Bidirectional Research

For patients suffering with chronic neuropathic pain the need for suitable novel therapies is imperative. Over recent years a contributing factor for the lack of development of new analgesics for neuropathic pain has been the mismatch of primary neuropathic pain assessment endpoints in preclinical vs. clinical trials. Despite continuous forward translation failures across diverse mechanisms, reflexive quantitative sensory testing remains the primary assessment endpoint for neuropathic pain and analgesia in animals. Restricting preclinical evaluation of pain and analgesia to exclusively reflexive outcomes is over simplified and can be argued not clinically relevant due to the continued lack of forward translation and failures in the clinic. The key to developing new analgesic treatments for neuropathic pain therefore lies in the development of clinically relevant endpoints that can translate preclinical animal results to human clinical trials. In this review we discuss this mismatch of primary neuropathic pain assessment endpoints, together with clinical and preclinical evidence that supports how bidirectional research is helping to validate new clinically relevant neuropathic pain assessment endpoints. Ethological behavioral endpoints such as burrowing and facial grimacing and objective measures such as electroencephalography provide improved translatability potential together with currently used quantitative sensory testing endpoints. By tailoring objective and subjective measures of neuropathic pain the translatability of new medicines for patients suffering with neuropathic pain will hopefully be improved.

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.

Parametrizing the Conditionally Gaussian Prior Model for Source Localization with Reference to the P20/N20 Component of Median Nerve SEP/SEF

In this article, we focused on developing the conditionally Gaussian hierarchical Bayesian model (CG-HBM), which forms a superclass of several inversion methods for source localization of brain activity using somatosensory evoked potential (SEP) and field (SEF) measurements. The goal of this proof-of-concept study was to improve the applicability of the CG-HBM as a superclass by proposing a robust approach for the parametrization of focal source scenarios. We aimed at a parametrization that is invariant with respect to altering the noise level and the source space size. The posterior difference between the gamma and inverse gamma hyperprior was minimized by optimizing the shape parameter, while a suitable range for the scale parameter can be obtained via the prior-over-measurement signal-to-noise ratio, which we introduce as a new concept in this study. In the source localization experiments, the primary generator of the P20/N20 component was detected in the Brodmann area 3b using the CG-HBM approach and a parameter range derived from the existing knowledge of the Tikhonov-regularized minimum norm estimate, i.e., the classical Gaussian prior model. Moreover, it seems that the detection of deep thalamic activity simultaneously with the P20/N20 component with the gamma hyperprior can be enhanced while using a close-to-optimal shape parameter value.

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.

Neural Representation Enhanced for Speech and Reduced for Background Noise With a Hearing Aid Noise Reduction Scheme During a Selective Attention Task

Objectives

Selectively attending to a target talker while ignoring multiple interferers (competing talkers and background noise) is more difficult for hearing-impaired (HI) individuals compared to normal-hearing (NH) listeners. Such tasks also become more difficult as background noise levels increase. To overcome these difficulties, hearing aids (HAs) offer noise reduction (NR) schemes. The objective of this study was to investigate the effect of NR processing (inactive, where the NR feature was switched off, vs. active, where the NR feature was switched on) on the neural representation of speech envelopes across two different background noise levels [+3 dB signal-to-noise ratio (SNR) and +8 dB SNR] by using a stimulus reconstruction (SR) method.

Design

To explore how NR processing supports the listeners’ selective auditory attention, we recruited 22 HI participants fitted with HAs. To investigate the interplay between NR schemes, background noise, and neural representation of the speech envelopes, we used electroencephalography (EEG). The participants were instructed to listen to a target talker in front while ignoring a competing talker in front in the presence of multi-talker background babble noise.

Results

The results show that the neural representation of the attended speech envelope was enhanced by the active NR scheme for both background noise levels. The neural representation of the attended speech envelope at lower (+3 dB) SNR was shifted, approximately by 5 dB, toward the higher (+8 dB) SNR when the NR scheme was turned on. The neural representation of the ignored speech envelope was modulated by the NR scheme and was mostly enhanced in the conditions with more background noise. The neural representation of the background noise was modulated (i.e., reduced) by the NR scheme and was significantly reduced in the conditions with more background noise. The neural representation of the net sum of the ignored acoustic scene (ignored talker and background babble) was not modulated by the NR scheme but was significantly reduced in the conditions with a reduced level of background noise. Taken together, we showed that the active NR scheme enhanced the neural representation of both the attended and the ignored speakers and reduced the neural representation of background noise, while the net sum of the ignored acoustic scene was not enhanced.

Conclusion

Altogether our results support the hypothesis that the NR schemes in HAs serve to enhance the neural representation of speech and reduce the neural representation of background noise during a selective attention task. We contend that these results provide a neural index that could be useful for assessing the effects of HAs on auditory and cognitive processing in HI populations.

Influence of Patient-Specific Head Modeling on EEG Source Imaging

Electromagnetic source imaging (ESI) techniques have become one of the most common alternatives for understanding cognitive processes in the human brain and for guiding possible therapies for neurological diseases. However, ESI accuracy strongly depends on the forward model capabilities to accurately describe the subject's head anatomy from the available structural data. Attempting to improve the ESI performance, we enhance the brain structure model within the individual-defined forward problem formulation, combining the head geometry complexity of the modeled tissue compartments and the prior knowledge of the brain tissue morphology. We validate the proposed methodology using 25 subjects, from which a set of magnetic-resonance imaging scans is acquired, extracting the anatomical priors and an electroencephalography signal set needed for validating the ESI scenarios. Obtained results confirm that incorporating patient-specific head models enhances the performed accuracy and improves the localization of focal and deep sources.

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.