Influence of tissue resistivities on neuromagnetic fields and electric potentials studied with a finite element model of the head

Causal Contributions and Interhemispheric Interactions of the Left and Right Supramarginal Gyri in Vocal Feedback Control: Insights From Dual-Site Transcranial Magnetic Stimulation

The bilateral supramarginal gyri (SMGs) have been implicated in sensorimotor control of speech production, yet their precise roles and interhemispheric interactions are poorly understood. This event-related potential study employed dual-site continuous theta burst stimulation (c-TBS) over the bilateral SMGs simultaneously to investigate their functional dynamics in vocal motor control. Following unilateral and bilateral c-TBS over the SMG as well as sham stimulation, participants vocalized the vowel sounds while exposed to unexpected pitch perturbations in auditory feedback. Unilateral real c-TBS paired with contralateral sham stimulation led to reduced vocal compensation magnitudes and latencies and decreased P2 responses compared to bilateral sham stimulation, with no differences between left and right SMG stimulation. Source localization revealed that decreased P2 responses following left SMG stimulation localized to left-lateralized dorsolateral prefrontal cortex, supplementary motor area, SMG, middle temporal gyrus, and temporo-parietal junction, whereas such decreases following right SMG stimulation involved left-lateralized primary motor cortex, premotor cortex, and middle temporal gyrus. These findings suggest that both SMGs are causally involved in vocal feedback control through distinct but interconnected networks. Surprisingly, dual-site c-TBS over the bilateral SMG did not alter vocal compensation or cortical activity, suggesting an interhemispheric balancing mechanism for fine-tuning vocal production. Our results offer novel insights into the bihemispheric coordination of auditory-vocal integration, highlighting potential treatment for speech disorders by modulating interhemispheric interactions.

Are neurasthenia and depression the same disease entity? An electroencephalography study

BACKGROUND

The neurasthenia-depression controversy has lasted for several decades. It is challenging to solve the argument by symptoms alone for syndrome-based disease classification. Our aim was to identify objective electroencephalography (EEG) measures that can differentiate neurasthenia from major depressive disorder (MDD).

METHODS

Both electronic medical information records and EEG records from patients with neurasthenia and MDD were gathered. The demographic and clinical characteristics, EEG power spectral density, and functional connectivity were compared between the neurasthenia and MDD groups. Machine Learning methods such as random forest, logistic regression, support vector machines, and k nearest neighbors were also used for classification between groups to extend the identification that there is a significant different pattern between neurasthenia and MDD.

RESULTS

We analyzed 305 patients with neurasthenia and 45 patients with MDD. Compared with the MDD group, patients with neurasthenia reported more somatic symptoms and less emotional symptoms (p < 0.05). Moreover, lower theta connectivity was observed in patients with neurasthenia compared to those with MDD (p < 0.01). Among the classification models, random forest performed best with an accuracy of 0.93, area under the receiver operating characteristic curve of 0.97, and area under the precision-recall curve of 0.96. The essential feature contributing to the model was the theta connectivity.

LIMITATIONS

This is a retrospective study, and medical records may not include all the details of a patient's syndrome. The sample size of the MDD group was smaller than that of the neurasthenia group.

CONCLUSION

Neurasthenia and MDD are different not only in symptoms but also in brain activities.

M/EEG source localization for both subcortical and cortical sources using a convolutional neural network with a realistic head conductivity model

While electroencephalography (EEG) and magnetoencephalography (MEG) are well-established noninvasive methods in neuroscience and clinical medicine, they suffer from low spatial resolution. Electrophysiological source imaging (ESI) addresses this by noninvasively exploring the neuronal origins of M/EEG signals. Although subcortical structures are crucial to many brain functions and neuronal diseases, accurately localizing subcortical sources of M/EEG remains particularly challenging, and the feasibility is still a subject of debate. Traditional ESIs, which depend on explicitly defined regularization priors, have struggled to set optimal priors and accurately localize brain sources. To overcome this, we introduced a data-driven, deep learning-based ESI approach without the need for these priors. We proposed a four-layered convolutional neural network (4LCNN) designed to locate both subcortical and cortical sources underlying M/EEG signals. We also employed a sophisticated realistic head conductivity model using the state-of-the-art segmentation method of ten different head tissues from individual MRI data to generate realistic training data. This is the first attempt at deep learning-based ESI targeting subcortical regions. Our method showed excellent accuracy in source localization, particularly in subcortical areas compared to other methods. This was validated through M/EEG simulations, evoked responses, and invasive recordings. The potential for accurate source localization of the 4LCNNs demonstrated in this study suggests future contributions to various research endeavors such as the clinical diagnosis, understanding of the pathophysiology of various neuronal diseases, and basic brain functions.

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.

Human brain imaging with high-density electroencephalography: Techniques and applications

Electroencephalography (EEG) is a technique for non-invasively measuring neuronal activity in the human brain using electrodes placed on the participant's scalp. With the advancement of digital technologies, EEG analysis has evolved over time from the qualitative analysis of amplitude and frequency modulations to a comprehensive analysis of the complex spatiotemporal characteristics of the recorded signals. EEG is now considered a powerful tool for measuring neural processes in the same time frame in which they happen (i.e. the subsecond range). However, it is commonly argued that EEG suffers from low spatial resolution, which makes it difficult to localize the generators of EEG activity accurately and reliably. Today, the availability of high-density EEG (hdEEG) systems, combined with methods for incorporating information on head anatomy and sophisticated source-localization algorithms, has transformed EEG into an important neuroimaging tool. hdEEG offers researchers and clinicians a rich and varied range of applications. It can be used not only for investigating neural correlates in motor and cognitive neuroscience experiments, but also for clinical diagnosis, particularly in the detection of epilepsy and the characterization of neural impairments in a wide range of neurological disorders. Notably, the integration of hdEEG systems with other physiological recordings, such as kinematic and/or electromyography data, might be especially beneficial to better understand the neuromuscular mechanisms associated with deconditioning in ageing and neuromotor disorders, by mapping the neurokinematic and neuromuscular connectivity patterns directly in the brain.

Large-scale analysis of neural activity and connectivity from high-density electroencephalographic data

INTRODUCTION

High-density electroencephalography (hdEEG) is a technique used for the characterization of the neural activity and connectivity in the human brain. The analysis of EEG data involves several steps, including signal pre-processing, head modelling, source localization and activity/connectivity quantification. Visual check of the analysis steps is often necessary, making the process time- and resource-consuming and, therefore, not feasible for large datasets.

FINDINGS

Here we present the Noninvasive Electrophysiology Toolbox (NET), an open-source software for large-scale analysis of hdEEG data, running on the cross-platform MATLAB environment. NET combines all the tools required for a complete hdEEG analysis workflow, from raw signals to final measured values. By relying on reconstructed neural signals in the brain, NET can perform traditional analyses of time-locked neural responses, as well as more advanced functional connectivity and brain mapping analyses. The extracted quantitative neural data can be exported to provide broad compatibility with other software.

CONCLUSIONS

NET is freely available (https://github.com/bind-group-kul/net) under the GNU public license for non-commercial use and open-source development, together with a graphical user interface (GUI) and a user tutorial. While NET can be used interactively with the GUI, it is primarily aimed at unsupervised automation to process large hdEEG datasets efficiently. Its implementation creates indeed a highly customizable program suitable for analysis automation and tight integration into existing workflows.

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.

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.

Role of the volume conductor on simulations of local field potential recordings from deep brain stimulation electrodes

OBJECTIVE

Local field potential (LFP) recordings from deep brain stimulation (DBS) electrodes are commonly used in research analyses, and are beginning to be used in clinical practice. Computational models of DBS LFPs provide tools for investigating the biophysics and neural synchronization that underlie LFP signals. However, technical standards for DBS LFP model parameterization remain to be established. Therefore, the goal of this study was to evaluate the role of the volume conductor (VC) model complexity on simulated LFP signals in the subthalamic nucleus (STN).

APPROACH

We created a detailed human head VC model that explicitly represented the inhomogeneity and anisotropy associated with 12 different tissue structures. This VC model represented our "gold standard" for technical detail and electrical realism. We then incrementally decreased the complexity of the VC model and quantified the impact on the simulated LFP recordings. Identical STN neural source activity was used when comparing the different VC model variants. Results Ignoring tissue anisotropy reduced the simulated LFP amplitude by ~12%, while eliminating soft tissue heterogeneity had a negligible effect on the recordings. Simplification of the VC model to consist of a single homogenous isotropic tissue medium with a conductivity of 0.215 S/m contributed an additional ~3% to the error.

SIGNIFICANCE

Highly detailed VC models do generate different results than simplified VC models. However, with errors in the range of ~15%, the use of a well-parameterized simple VC model is likely to be acceptable in most contexts for DBS LFP modeling.

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.

Lumbar trans-spinal direct current stimulation: A modeling-experimental approach to dorsal root ganglia stimulation

Introduction

The excitability of spinal motor neurons (MN) can be altered through subthreshold currents, such as transcutaneous spinal direct-current stimulation (tsDCS). Current evidence shows that tsDCS can interfere with ascending somatosensory pathways and lower motor neurons' (LMN) excitability, which points to its therapeutic potential for repairing altered spinal responses. We aim to define the best tsDCS montage for maximizing the electric field (E-field) in the lumbar spinal cord (L-SC) by computer modeling; and to apply this montage to measure the effect on LMN excitability and somatosensory evoked potentials (SSEP).

Methods

A human volume conductor model was obtained from an available database. The E-field distribution was calculated considering three different electrode settings aiming at maximizing the field at L-SC and right dorsal root ganglia (DRG). The best electrode setting was then selected and applied in a blind crossover pseudo-randomized study including 14 subjects. tsDCS was delivered for 15 min (cathodal vs. sham) over L2 vertebra level (4 mA, 144 mC/cm²), and its effect on F-waves, H-reflex (including homosynaptic depression, HD) and SSEPs was investigated in the lower limbs.

Results

All simulated montages showed higher current density and E-field magnitudes between the electrodes (>0.15 V/m), with a major longitudinal component and with rostral-caudal direction. The induced E-field involved the sensory ganglia and was maximum in the right T8-left L2 montage, which was the one selected for the experimental protocol. We disclosed a statistically significant increase of the H-reflex amplitude at 0.1 Hz, after cathodal tsDCS (c-tsDCS) on both sides. No other significant change was observed.

Discussion

Our results can suggest the c-tsDCS applied to the L-SC and DRG can modulate synaptic efficiency increasing lower motor neurons response to Ia fibers excitation. The possible implications of our findings for treating clinical conditions will be addressed in future studies.

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.

Variation of cerebrospinal fluid in specific regions regulates focality in transcranial direct current stimulation

Background

Conventionally, transcranial direct current stimulation (tDCS) aims to focalize the current reaching the target region-of-interest (ROI). The focality can be quantified by the dose-target-determination-index (DTDI). Despite having a uniform tDCS setup, some individuals receive focal stimulation (high DTDI) while others show reduced focality ("non-focal"). The volume of cerebrospinal fluid (CSF), gray matter (GM), and white matter (WM) underlying each ROI govern the tDCS current distribution inside the brain, thereby regulating focality.

Aim

To determine the regional volume parameters that differentiate the focal and non-focal groups.

Methods

T1-weighted images of the brain from 300 age-sex matched adults were divided into three equal groups- (a) Young (20 ≤ × < 40 years), (b) Middle (40 ≤ × < 60 years), and (c) Older (60 ≤ × < 80 years). For each group, inter and intra-hemispheric montages with electrodes at (1) F3 and right supraorbital region (F3-RSO), and (2) CP5 and Cz (CP5-Cz) were simulated, targeting the left- Dorsolateral Prefrontal Cortex (DLPFC) and -Inferior Parietal Lobule (IPL), respectively. Both montages were simulated for two current doses (1 and 2 mA). For each individual head simulated for a tDCS configuration (montage and dose), the current density at each region-of-interest (ROI) and their DTDI were calculated. The individuals were categorized into two groups- (1) Focal (DTDI ≥ 0.75), and (2) Non-focal (DTDI < 0.75). The regional volume of CSF, GM, and WM of all the ROIs was determined. For each tDCS configuration and ROI, three 3-way analysis of variance was performed considering- (i) GM, (ii) WM, and (iii) CSF as the dependent variable (DV). The age group, sex, and focality group were the between-subject factors. For a given ROI, if any of the 3 DV's showed a significant main effect or interaction involving the focality group, then that ROI was classified as a "focal ROI."

Results

Regional CSF was the principal determinant of focality. For interhemispheric F3-RSO montage, interaction effect (p < 0.05) of age and focality was observed at Left Caudate Nucleus, with the focal group exhibiting higher CSF volume. The CSF volume of focal ROI correlated positively (r ∼ 0.16, p < 0.05) with the current density at the target ROI (DLPFC). For intrahemispheric CP5-Cz montage, a significant (p < 0.05) main effect was observed at the left pre- and post-central gyrus, with the focal group showing lower CSF volume. The CSF volume correlated negatively (r ∼ -0.16, p < 0.05) with current density at left IPL. The results were consistent for both current doses.

Conclusion

The CSF channels the flow of tDCS current between electrodes with focal ROIs acting like reservoirs of current. The position of focal ROI in the channel determines the stimulation intensity at the target ROI. For focal stimulation in interhemispheric F3-RSO, the proximity of focal ROI reserves the current density at the target ROI (DLPFC). In contrast, for intrahemispheric montage (CP5-Cz), the far-end location of focal ROI reduces the current density at the target (IPL).

Anatomical harmonics basis based brain source localization with application to epilepsy

Brain Source Localization (BSL) using Electroencephalogram (EEG) has been a useful noninvasive modality for the diagnosis of epileptogenic zones, study of evoked related potentials, and brain disorders. The inverse solution of BSL is limited by high computational cost and localization error. The performance is additionally limited by head shape assumption and the corresponding harmonics basis function. In this work, an anatomical harmonics basis (Spherical Harmonics (SH), and more particularly Head Harmonics (H²)) based BSL is presented. The spatio-temporal four shell head model is formulated in SH and H² domain. The anatomical harmonics domain formulation leads to dimensionality reduction and increased contribution of source eigenvalues, resulting in decreased computation and increased accuracy respectively. The performance of spatial subspace based Multiple Signal Classification (MUSIC) and Recursively Applied and Projected (RAP)-MUSIC method is compared with the proposed SH and H² counterparts on simulated data. SH and H² domain processing effectively resolves the problem of high computational cost without sacrificing the inverse source localization accuracy. The proposed H² MUSIC was additionally validated for epileptogenic zone localization on clinical EEG data. The proposed framework offers an effective solution to clinicians in automated and time efficient seizure localization.

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.

Estimating the influence of stroke lesions on MEG source reconstruction

Source reconstruction of magnetoencephalography (MEG) has been used to assess brain reorganization after brain damage, such as stroke. Lesions result in parts of the brain having an electrical conductivity that differs from the normal values. The effect this has on the forward solutions (i.e., the propagation of electric currents and magnetic fields generated by cortical activity) is well predictable. However, their influence on source localization results is not well characterized and understood. This is specifically a concern for patient studies with asymmetric (i.e., within one hemisphere) lesions focusing on asymmetric and lateralized brain activity, such as language. In particular, it is good practice to consider the level of geometrical detail that is necessary to compute and interpret reliable source reconstruction results. To understand the effect of lesions on source estimates and propose recommendations to researchers working with clinical data, in this study we consider the trade off between improved accuracy and the additional effort to compute more realistic head models, with the aim to answer the question whether the additional effort is worth it. We simulated and analyzed the effects of a stroke lesion (i.e., an asymmetrically distributed CSF-filled cavity) in the head model with three different sizes and locations when performing MEG source reconstruction using a finite element method (FEM). We compared the effect of the lesion with a homogeneous head model that neglects the lesion. We computed displacement and attenuation/amplification maps to quantify the localization errors and signal magnitude modulation. We conclude that brain lesions leading to asymmetrically distributed CSF-filled cavities should be modeled when performing MEG source reconstruction, especially when investigating deep sources or post-stroke hemispheric lateralization of functions. The strongest effects are not only visible in perilesional areas, but can extend up to 20 mm from the lesion. Bigger lesions lead to stronger effects impacting larger areas, independently from the lesion location. Lastly, we conclude that more priority should be given to usability and accessibility of the required computational tools, to allow researchers with less technical expertise to use the improved methods that are available but currently not widely adopted yet.

Assessing Neurokinematic and Neuromuscular Connectivity During Walking Using Mobile Brain-Body Imaging

Gait is a common but rather complex activity that supports mobility in daily life. It requires indeed sophisticated coordination of lower and upper limbs, controlled by the nervous system. The relationship between limb kinematics and muscular activity with neural activity, referred to as neurokinematic and neuromuscular connectivity (NKC/NMC) respectively, still needs to be elucidated. Recently developed analysis techniques for mobile high-density electroencephalography (hdEEG) recordings have enabled investigations of gait-related neural modulations at the brain level. To shed light on gait-related neurokinematic and neuromuscular connectivity patterns in the brain, we performed a mobile brain/body imaging (MoBI) study in young healthy participants. In each participant, we collected hdEEG signals and limb velocity/electromyography signals during treadmill walking. We reconstructed neural signals in the alpha (8–13 Hz), beta (13–30 Hz), and gamma (30–50 Hz) frequency bands, and assessed the co-modulations of their power envelopes with myogenic/velocity envelopes. Our results showed that the myogenic signals have larger discriminative power in evaluating gait-related brain-body connectivity with respect to kinematic signals. A detailed analysis of neuromuscular connectivity patterns in the brain revealed robust responses in the alpha and beta bands over the lower limb representation in the primary sensorimotor cortex. There responses were largely contralateral with respect to the body sensor used for the analysis. By using a voxel-wise analysis of variance on the NMC images, we revealed clear modulations across body sensors; the variability across frequency bands was relatively lower, and below significance. Overall, our study demonstrates that a MoBI platform based on hdEEG can be used for the investigation of gait-related brain-body connectivity. Future studies might involve more complex walking conditions to gain a better understanding of fundamental neural processes associated with gait control, or might be conducted in individuals with neuromotor disorders to identify neural markers of abnormal gait.

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.

Connectivity in Large-Scale Resting-State Brain Networks Is Related to Motor Learning: A High-Density EEG Study

Previous research has shown that resting-state functional connectivity (rsFC) between different brain regions (seeds) is related to motor learning and motor memory consolidation. Using high-density electroencephalography (hdEEG), we addressed this question from a brain network perspective. Specifically, we examined frequency-dependent functional connectivity in resting-state networks from twenty-nine young healthy participants before and after they were trained on a motor sequence learning task. Consolidation was assessed with an overnight retest on the motor task. Our results showed training-related decreases in gamma-band connectivity within the motor network, and between the motor and functionally distinct resting-state networks including the attentional network. Brain-behavior correlation analyses revealed that baseline beta, delta, and theta rsFC were related to subsequent motor learning and memory consolidation such that lower connectivity within the motor network and between the motor and several distinct resting-state networks was correlated with better learning and overnight consolidation. Lastly, training-related increases in beta-band connectivity between the motor and the visual networks were related to greater consolidation. Altogether, our results indicate that connectivity in large-scale resting-state brain networks is related to—and modulated by—motor learning and memory consolidation processes. These finding corroborate previous seed-based connectivity research and provide evidence that frequency-dependent functional connectivity in resting-state networks is critically linked to motor learning and memory consolidation.

Frequency-dependent modulation of neural oscillations across the gait cycle

Balance and walking are fundamental to support common daily activities. Relatively accurate characterizations of normal and impaired gait features were attained at the kinematic and muscular levels. Conversely, the neural processes underlying gait dynamics still need to be elucidated. To shed light on gait‐related modulations of neural activity, we collected high‐density electroencephalography (hdEEG) signals and ankle acceleration data in young healthy participants during treadmill walking. We used the ankle acceleration data to segment each gait cycle in four phases: initial double support, right leg swing, final double support, left leg swing. Then, we processed hdEEG signals to extract neural oscillations in alpha, beta, and gamma bands, and examined event‐related desynchronization/synchronization (ERD/ERS) across gait phases. Our results showed that ERD/ERS modulations for alpha, beta, and gamma bands were strongest in the primary sensorimotor cortex (M1), but were also found in premotor cortex, thalamus and cerebellum. We observed a modulation of neural oscillations across gait phases in M1 and cerebellum, and an interaction between frequency band and gait phase in premotor cortex and thalamus. Furthermore, an ERD/ERS lateralization effect was present in M1 for the alpha and beta bands, and in the cerebellum for the beta and gamma bands. Overall, our findings demonstrate that an electrophysiological source imaging approach based on hdEEG can be used to investigate dynamic neural processes of gait control. Future work on the development of mobile hdEEG‐based brain–body imaging platforms may enable overground walking investigations, with potential applications in the study of gait disorders.

Neural oscillations during motor imagery of complex gait: an HdEEG study

The aim of this study was to investigate differences between usual and complex gait motor imagery (MI) task in healthy subjects using high-density electroencephalography (hdEEG) with a MI protocol. We characterized the spatial distribution of α- and β-bands oscillations extracted from hdEEG signals recorded during MI of usual walking (UW) and walking by avoiding an obstacle (Dual-Task, DT). We applied a source localization algorithm to brain regions selected from a large cortical-subcortical network, and then we analyzed α and β bands Event-Related Desynchronizations (ERDs). Nineteen healthy subjects visually imagined walking on a path with (DT) and without (UW) obstacles. Results showed in both gait MI tasks, α- and β-band ERDs in a large cortical-subcortical network encompassing mostly frontal and parietal regions. In most of the regions, we found α- and β-band ERDs in the DT compared with the UW condition. Finally, in the β band, significant correlations emerged between ERDs and scores in imagery ability tests. Overall we detected MI gait-related α- and β-band oscillations in cortical and subcortical areas and significant differences between UW and DT MI conditions. A better understanding of gait neural correlates may lead to a better knowledge of pathophysiology of gait disturbances in neurological diseases.

In-vivo estimation of tissue electrical conductivities of a rabbit eye for precise simulation of electric field distributions during ocular iontophoresis

Precise estimation of electrical conductivity of the eyes is important for the accurate analysis of electric field distributions in the eyes during ocular iontophoresis. In this study, we estimated the tissue electrical conductivities of a rabbit eye, which has been widely employed for neuro-ophthalmological experiments, through an in vivo experiment for the first time. Electrical potentials were measured at multiple locations on the skin, while weak currents were transmitted into the skin via two surface electrodes attached to the skin around the eye. A finite element model was constructed to calculate the electric potentials at the measurement locations. The conductivity values of different tissues were then estimated using an optimization procedure to minimize the difference between the measured and calculated electric potentials. The accuracy of the estimated tissue conductivity values of the rabbit eye was validated by comparing the measured and calculated electric potential values for different electrode montages. Further multi-physical analyses of iontophoretic drug delivery to the rabbit eye showed a significant influence of the conductivity profile on the resultant particle distribution. Overall, our results provide an important reference for the tissue electrical conductivity values of the rabbit eye, which could be further utilized for designing new medical devices for delivering electric fields to the eyes, such as transorbital and transscleral electrical stimulations.

A multi-step blind source separation approach for the attenuation of artifacts in mobile high-density electroencephalography data

Objective.Electroencephalography (EEG) is a widely used technique to address research questions about brain functioning, from controlled laboratorial conditions to naturalistic environments. However, EEG data are affected by biological (e.g. ocular, myogenic) and non-biological (e.g. movement-related) artifacts, which-depending on their extent-may limit the interpretability of the study results. Blind source separation (BSS) approaches have demonstrated to be particularly promising for the attenuation of artifacts in high-density EEG (hdEEG) data. Previous EEG artifact removal studies suggested that it may not be optimal to use the same BSS method for different kinds of artifacts.Approach.In this study, we developed a novel multi-step BSS approach to optimize the attenuation of ocular, movement-related and myogenic artifacts from hdEEG data. For validation purposes, we used hdEEG data collected in a group of healthy participants in standing, slow-walking and fast-walking conditions. During part of the experiment, a series of tone bursts were used to evoke auditory responses. We quantified event-related potentials (ERPs) using hdEEG signals collected during an auditory stimulation, as well as the event-related desynchronization (ERD) by contrasting hdEEG signals collected in walking and standing conditions, without auditory stimulation. We compared the results obtained in terms of auditory ERP and motor-related ERD using the proposed multi-step BSS approach, with respect to two classically used single-step BSS approaches.Main results. The use of our approach yielded the lowest residual noise in the hdEEG data, and permitted to retrieve stronger and more reliable modulations of neural activity than alternative solutions. Overall, our study confirmed that the performance of BSS-based artifact removal can be improved by using specific BSS methods and parameters for different kinds of artifacts.Significance.Our technological solution supports a wider use of hdEEG-based source imaging in movement and rehabilitation studies, and contributes to the further development of mobile brain/body imaging applications.

Probabilistic comparison of gray and white matter coverage between depth and surface intracranial electrodes in epilepsy

In this study, we quantified the coverage of gray and white matter during intracranial electroencephalography in a cohort of epilepsy patients with surface and depth electrodes. We included 65 patients with strip electrodes (n = 12), strip and grid electrodes (n = 24), strip, grid, and depth electrodes (n = 7), or depth electrodes only (n = 22). Patient-specific imaging was used to generate probabilistic gray and white matter maps and atlas segmentations. Gray and white matter coverage was quantified using spherical volumes centered on electrode centroids, with radii ranging from 1 to 15 mm, along with detailed finite element models of local electric fields. Gray matter coverage was highly dependent on the chosen radius of influence (RoI). Using a 2.5 mm RoI, depth electrodes covered more gray matter than surface electrodes; however, surface electrodes covered more gray matter at RoI larger than 4 mm. White matter coverage and amygdala and hippocampal coverage was greatest for depth electrodes at all RoIs. This study provides the first probabilistic analysis to quantify coverage for different intracranial recording configurations. Depth electrodes offer increased coverage of gray matter over other recording strategies if the desired signals are local, while subdural grids and strips sample more gray matter if the desired signals are diffuse.

Influence of layered skull modeling on the frequency sensitivity and target accuracy in simulations of transcranial current stimulation

With the development of electrical stimulation technology, especially the emergence of temporally interfering (TI) stimulation, it is necessary to discuss the influence of current frequency on stimulation intensity. Accurate skull modeling is important for transcranial current stimulation (tCS) simulation prediction because of its large role in dispersing current. In this study, we simulated different frequencies of transcranial alternating current stimulation (tACS) and TI stimulation in single-layer and layered skull model, compared the electric field via error parameters such as the relative difference measure and relative magnification factor. Pearson correlation analysis and t-test were used to measure the differences in envelope amplitude. The results showed that the intensity of electric field in the brain generated by per unit of stimulation current will increase with current frequency, and the layered skull model had a better response to frequency. An obvious pattern difference was found between the electric fields of the layered and single-layer skull individualized models. For TI stimulation, the Pearson correlation coefficient between the envelope distribution of the layered skull model and the single-layer skull was only 0.746 in the individualized model, which is clearly lower than the correlation coefficient of 0.999 determined from the spherical model. Higher carrier frequencies seemed to be easier to generate a large enough brain electric field envelope in TI stimulation. In conclusion, we recommend using layered skull models instead of single-layer skull models in tCS (particularly TI stimulation) simulation studies in order to improve the accuracy of the prediction of stimulus intensity and stimulus target.

Influence of Anisotropic White Matter on Electroosmotic Flow Induced by Direct Current

Treatment of cerebral edema remains a major challenge in clinical practice and new innovative therapies are needed. This study presents a novel approach for mitigating cerebral edema by inducing bulk fluid transport utilizing the brain’s electroosmotic property using an anatomically detailed finite element head model incorporating anisotropy in the white matter (WM). Three representative anisotropic conductivity algorithms are employed for the WM and compared with isotropic WM. The key results are (1) the electroosmotic flow (EOF) is driven from the edema region to the subarachnoid space under an applied electric field with its magnitude linearly correlated to the electric field and direction following current flow pathways; (2) the extent of EOF distribution variation correlates highly with the degree of the anisotropic ratio of the WM regions; (3) the directions of the induced EOF in the anisotropic models deviate from its isotropically defined pathways and tend to move along the principal fiber direction. The results suggest WM anisotropy should be incorporated in head models for more reliable EOF evaluations for cerebral edema mitigation and demonstrate the promise of the electroosmosis based approach to be developed as a new therapy for edema treatment as evaluated with enhanced head models incorporating WM anisotropy.

Validating Patient-Specific Finite Element Models of Direct Electrocortical Stimulation

Direct electrocortical stimulation (DECS) with electrocorticography electrodes is an established therapy for epilepsy and an emerging application for stroke rehabilitation and brain-computer interfaces. However, the electrophysiological mechanisms that result in a therapeutic effect remain unclear. Patient-specific computational models are promising tools to predict the voltages in the brain and better understand the neural and clinical response to DECS, but the accuracy of such models has not been directly validated in humans. A key hurdle to modeling DECS is accurately locating the electrodes on the cortical surface due to brain shift after electrode implantation. Despite the inherent uncertainty introduced by brain shift, the effects of electrode localization parameters have not been investigated. The goal of this study was to validate patient-specific computational models of DECS against in vivo voltage recordings obtained during DECS and quantify the effects of electrode localization parameters on simulated voltages on the cortical surface. We measured intracranial voltages in six epilepsy patients during DECS and investigated the following electrode localization parameters: principal axis, Hermes, and Dykstra electrode projection methods combined with 0, 1, and 2 mm of cerebral spinal fluid (CSF) below the electrodes. Greater CSF depth between the electrode and cortical surface increased model errors and decreased predicted voltage accuracy. The electrode localization parameters that best estimated the recorded voltages across six patients with varying amounts of brain shift were the Hermes projection method and a CSF depth of 0 mm (r = 0.92 and linear regression slope = 1.21). These results are the first to quantify the effects of electrode localization parameters with in vivo intracranial recordings and may serve as the basis for future studies investigating the neuronal and clinical effects of DECS for epilepsy, stroke, and other emerging closed-loop applications.

Electroosmosis Based Novel Treatment Approach for Cerebral Edema

OBJECTIVE

Cerebral edema characterized as an abnormal accumulation of interstitial fluid has not been treated effectively. We propose a novel edema treatment approach to drive edematous fluid out of the brain by direct current utilizing brain tissue's electroosmotic property.

METHODS

A finite element (FE) head model is developed and employed to assess the feasibility of the approach. First, the capacity of the model for electric field prediction is validated against human experiments. Second, two electrode configurations (S and D-montage) are designed to evaluate the distribution of the electric field, electroosmotic flow (EOF), current density, and temperature across the brain under an applied direct current.

RESULTS

The S-montage is shown to induce an average EOF velocity of 7e-4 mm/s underneath the anode by a voltage of 15 V, and the D-montage induces a velocity of 9e-4 mm/s by a voltage of 5 V. Meanwhile, the brain temperature in both configurations is below 38 °C, which is within the safety range. Further, the magnitude of EOF is proportional to the electric field, and the EOF direction follows the current flow from anode to cathode. The EOF velocity in the white matter is significantly higher than that in the gray matter under the anode where the fluid is to be drawn out.

CONCLUSION

The proposed electroosmosis based approach allows alleviating brain edema within the critical time window by direct current.

SIGNIFICANCE

The approach may be further developed as a new treatment solely or as a complement to existing conventional treatments of edema.

RT-NET: real-time reconstruction of neural activity using high-density electroencephalography

High-density electroencephalography (hdEEG) has been successfully used for large-scale investigations of neural activity in the healthy and diseased human brain. Because of their high computational demand, analyses of source-projected hdEEG data are typically performed offline. Here, we present a real-time noninvasive electrophysiology toolbox, RT-NET, which has been specifically developed for online reconstruction of neural activity using hdEEG. RT-NET relies on the Lab Streaming Layer for acquiring raw data from a large number of EEG amplifiers and for streaming the processed data to external applications. RT-NET estimates a spatial filter for artifact removal and source activity reconstruction using a calibration dataset. This spatial filter is then applied to the hdEEG data as they are acquired, thereby ensuring low latencies and computation times. Overall, our analyses show that RT-NET can estimate real-time neural activity with performance comparable to offline analysis methods. It may therefore enable the development of novel brain–computer interface applications such as source-based neurofeedback.

Detection of Resting-State Functional Connectivity from High-Density Electroencephalography Data: Impact of Head Modeling Strategies

Recent technological advances have been permitted to use high-density electroencephalography (hdEEG) for the estimation of functional connectivity and the mapping of resting-state networks (RSNs). The reliable estimate of activity and connectivity from hdEEG data relies on the creation of an accurate head model, defining how neural currents propagate from the cortex to the sensors placed over the scalp. To the best of our knowledge, no study has been conducted yet to systematically test to what extent head modeling accuracy impacts on EEG-RSN reconstruction. To address this question, we used 256-channel hdEEG data collected in a group of young healthy participants at rest. We first estimated functional connectivity in EEG-RSNs by means of band-limited power envelope correlations, using neural activity estimated with an optimized analysis workflow. Then, we defined a series of head models with different levels of complexity, specifically testing the effect of different electrode positioning techniques and head tissue segmentation methods. We observed that robust EEG-RSNs can be obtained using a realistic head model, and that inaccuracies due to head tissue segmentation impact on RSN reconstruction more than those due to electrode positioning. Additionally, we found that EEG-RSN robustness to head model variations had space and frequency specificity. Overall, our results may contribute to defining a benchmark for assessing the reliability of hdEEG functional connectivity measures.

Automated Head Tissue Modelling Based on Structural Magnetic Resonance Images for Electroencephalographic Source Reconstruction

In the last years, technological advancements for the analysis of electroencephalography (EEG) recordings have permitted to investigate neural activity and connectivity in the human brain with unprecedented precision and reliability. A crucial element for accurate EEG source reconstruction is the construction of a realistic head model, incorporating information on electrode positions and head tissue distribution. In this paper, we introduce MR-TIM, a toolbox for head tissue modelling from structural magnetic resonance (MR) images. The toolbox consists of three modules: 1) image pre-processing – the raw MR image is denoised and prepared for further analyses; 2) tissue probability mapping – template tissue probability maps (TPMs) in individual space are generated from the MR image; 3) tissue segmentation – information from all the TPMs is integrated such that each voxel in the MR image is assigned to a specific tissue. MR-TIM generates highly realistic 3D masks, five of which are associated with brain structures (brain and cerebellar grey matter, brain and cerebellar white matter, and brainstem) and the remaining seven with other head tissues (cerebrospinal fluid, spongy and compact bones, eyes, muscle, fat and skin). Our validation, conducted on MR images collected in healthy volunteers and patients as well as an MR template image from an open-source repository, demonstrates that MR-TIM is more accurate than alternative approaches for whole-head tissue segmentation. We hope that MR-TIM, by yielding an increased precision in head modelling, will contribute to a more widespread use of EEG as a brain imaging technique.

Effect of skin conductivity on the electric field induced by transcranial stimulation techniques in different head models

This study aims at quantifying the effect that using different skin conductivity values has on the estimation of the electric (E)-field distribution induced by transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS) in the brain of two anatomical models. The induced E-field was calculated with numerical simulations inside MIDA and Duke models, assigning to the skin a conductivity value estimated from a multi-layered skin model and three values taken from literature. The effect of skin conductivity variations on the local E-field induced by tDCS in the brain was up to 70%. In TMS, minor local differences, in the order of 20%, were obtained in regions of interest for the onset of possible side effects. Results suggested that an accurate model of the skin is necessary in all numerical studies that aim at precisely estimating the E-field induced during TMS and tDCS applications. This also highlights the importance of further experimental studies on human skin characterization, especially at low frequencies.

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.

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.

Frequency-dependent functional connectivity in resting state networks

Functional magnetic resonance imaging studies have documented the resting human brain to be functionally organized in multiple large‐scale networks, called resting‐state networks (RSNs). Other brain imaging techniques, such as electroencephalography (EEG) and magnetoencephalography (MEG), have been used for investigating the electrophysiological basis of RSNs. To date, it is largely unclear how neural oscillations measured with EEG and MEG are related to functional connectivity in the resting state. In addition, it remains to be elucidated whether and how the observed neural oscillations are related to the spatial distribution of the network nodes over the cortex. To address these questions, we examined frequency‐dependent functional connectivity between the main nodes of several RSNs, spanning large part of the cortex. We estimated connectivity using band‐limited power correlations from high‐density EEG data collected in healthy participants. We observed that functional interactions within RSNs are characterized by a specific combination of neuronal oscillations in the alpha (8–13 Hz), beta (13–30 Hz), and gamma (30–80 Hz) bands, which highly depend on the position of the network nodes. This finding may contribute to a better understanding of the mechanisms through which neural oscillations support functional connectivity in the brain.

A Lead Field Two-Domain Model for Longitudinal Neural Tracts-Analytical Framework and Implications for Signal Bandwidth

Somatosensory evoked potentials are a well-established tool for assessing volley conduction in afferent neural pathways. However, from a clinical perspective, recording of spinal signals is still a demanding task due to the low amplitudes compared to relevant noise sources. Computer modeling is a powerful tool for gaining insight into signal genesis and, thus, for promoting future innovations in signal extraction. However, due to the complex structure of neural pathways, modeling is computationally demanding. We present a theoretical framework which allows computing the electric potential generated by a single axon in a body surface lead by the convolution of the neural lead field function with a propagating action potential term. The signal generated by a large cohort of axons was obtained by convoluting a single axonal signal with the statistical distribution of temporal dispersion of individual axonal signals. For establishing the framework, analysis was based on an analytical model. Our approach was further adopted for a numerical computation of body surface neuropotentials employing the lead field theory. Double convolution allowed straightforward analysis in the frequency domain. The highest frequency components occurred at the cellular membrane. A bandpass type spectral shape and a peak frequency of 1800 Hz was observed. The volume conductor transmitting the signal to the recording lead acted as an additional bandpass reducing the axonal peak frequency from 200 Hz to 500 Hz. The superposition of temporally dispersed axonal signals acted as an additional low-pass filter further reducing the compound action potential peak frequency from 90 Hz to 170 Hz. Our results suggest that the bandwidth of spinal evoked potentials might be narrower than the bandwidth requested by current clinical guidelines. The present findings will allow the optimization of noise suppression. Furthermore, our theoretical framework allows the adaptation in numerical methods and application in anatomically realistic geometries in future studies.

Encoding strategy mediates the effect of electrical stimulation over posterior parietal cortex on visual short-term memory

Over past decades, converging neuroimaging and electrophysiological findings have suggested a crucial role of posterior parietal cortex (PPC) in supporting the storage capacity of visual short-term memory (VSTM). Moreover, a few recent studies have shown that electrical stimulation over PPC can enhance VSTM capacity, making it a promising method for improving VSTM function. However, the reliability of these results is still in question because null findings have also been observed. Among studies that reported significant effects, some found increased VSTM capacity only in people with low capacity. Here, we hypothesized that subjects' encoding strategy might be a key source of these variable results. To directly test this hypothesis, we stimulated PPC using transcranial direct-current stimulation (tDCS) in male and female human subjects instructed to employ different encoding strategies during a VSTM recall task. We found that VSTM capacity was higher in subjects who were instructed to remember all items in the supra-capacity array of visual stimuli (i.e., the remember-all group), compared to subjects who were told to focus on a subset of these stimuli (i.e., the remember-subset group). As predicted, anodal tDCS over PPC significantly enhanced VSTM capacity only in the remember-subset group, but not in the remember-all group. Additionally, no effect of encoding strategy or its interaction with electrical stimulation was found on VSTM precision. Together, these results suggest that encoding strategy has a selective influence on VSTM capacity and this influence of encoding strategy mediates the effect of electrical stimulation over PPC on VSTM function.

Hand, foot and lip representations in primary sensorimotor cortex: a high-density electroencephalography study

The primary sensorimotor cortex plays a major role in the execution of movements of the contralateral side of the body. The topographic representation of different body parts within this brain region is commonly investigated through functional magnetic resonance imaging (fMRI). However, fMRI does not provide direct information about neuronal activity. In this study, we used high-density electroencephalography (hdEEG) to map the representations of hand, foot, and lip movements in the primary sensorimotor cortex, and to study their neural signatures. Specifically, we assessed the event-related desynchronization (ERD) in the cortical space. We found that the performance of hand, foot, and lip movements elicited an ERD in beta and gamma frequency bands. The primary regions showing significant beta- and gamma-band ERD for hand and foot movements, respectively, were consistent with previously reported using fMRI. We observed relatively weaker ERD for lip movements, which may be explained by the fact that less fine movement control was required. Overall, our study demonstrated that ERD based on hdEEG data can support the study of motor-related neural processes, with relatively high spatial resolution. An interesting avenue may be the use of hdEEG for deeper investigations into the pathophysiology of neuromotor disorders.

A retrospective evaluation of automated optimization of deep brain stimulation parameters

OBJECTIVE

We performed a retrospective analysis of an optimization algorithm for the computation of patient-specific multipolar stimulation configurations employing multiple independent current/voltage sources. We evaluated whether the obtained stimulation configurations align with clinical data and whether the optimized stimulation configurations have the potential to lead to an equal or better stimulation of the target region as manual programming, while reducing the time required for programming sessions.

APPROACH

For three patients (five electrodes) diagnosed with essential tremor, we derived optimized multipolar stimulation configurations using an approach that is suitable for the application in clinical practice. To evaluate the automatically derived stimulation settings, we compared them to the results of the monopolar review.

MAIN RESULTS

We observe a good agreement between the findings of the monopolar review and the optimized stimulation configurations, with the algorithm assigning the maximal voltage in the optimized multipolar pattern to the contact that was found to lead to the best therapeutic effect in the clinical monopolar review in all cases. Additionally, our simulation results predict that the optimized stimulation settings lead to the activation of an equal or larger volume fraction of the target compared to the manually determined settings in all cases.

SIGNIFICANCE

Our results demonstrate the feasibility of an automatic determination of optimal DBS configurations and motivate a further evaluation of the applied optimization algorithm.

Cervical trans-spinal direct current stimulation: a modelling-experimental approach

Background

Trans-spinal direct current stimulation (tsDCS) is a non-invasive technique with promising neuromodulatory effects on spinal cord (SC) circuitry. Computational studies are essential to guide effective tsDCS protocols for specific clinical applications. This study aims to combine modelling and experimental studies to determine the electrode montage that maximizes electric field (E-field) delivery during cervical tsDCS.

Methods

Current and E-field distributions in the cervical SC were predicted for four electrode montages in a human realistic model using computational methods. A double-blind crossover and randomized exploratory study was conducted using the montage that maximized E-field delivery. tsDCS was applied for 15 min in 10 healthy subjects (anodal, cathodal, sham, with polarity assigned to the cervical electrode), with a current intensity of 2.5 mA, resulting in a total current charge density delivery of 90 mC/cm². Upper limb motor (transcranial magnetic stimulation) and sensory evoked potentials (MEP, SEP), M-waves, H-reflex and F-wave responses were analysed. Central and peripheral conduction times were determined using MEP. Repeated measures ANOVA and Friedman test were used for statistical analysis (significance level α = 0.05).

Results

All montages presented higher current density and E-field magnitudes in the cervical SC region between the electrodes. However, electrodes at C3 and T3 spinous processes (C3-T3) originated the highest E-field magnitude (0.50 V/m). Using C3-T3 montage we observed significant changes in N9 SEP latency (p = 0.006), but significance did not persist in pairwise comparisons (sham-anodal: p = 0.022; sham-cathodal: p = 0.619; anodal-cathodal: p = 0.018; α = 0.017, Bonferroni corrected). MEP latency and central motor conduction time (CMCT) modified significantly on stimulation (p = 0.007 and p = 0.015, respectively). In addition, pairwise comparisons confirmed significant differences between sham and cathodal conditions after Bonferroni correction for MEP latency (sham-anodal: p = 0.868; sham-cathodal: p = 0.011; anodal-cathodal: p = 0.023) and CMCT (sham-anodal: p = 0.929; sham-cathodal: p = 0.010; anodal-cathodal: p = 0.034).

Conclusions

Computational models predicted higher E-field delivery in the cervical SC for the C3-T3 montage. Polarity-dependent effects in motor responses were reported using this montage consistent with spinal motor modulation. tsDCS experimental protocol designs should be guided by modelling studies to improve effectiveness.

Hemodynamic Correlates of Electrophysiological Activity in the Default Mode Network

Hemodynamic fluctuations in the default mode network (DMN), observed through functional magnetic resonance imaging (fMRI), have been linked to electrophysiological oscillations detected by electroencephalography (EEG). It has been reported that, among the electrophysiological oscillations, those in the alpha frequency range (8–13 Hz) are the most dominant during resting state. We hypothesized that DMN spatial configuration closely depends on the specific neuronal oscillations considered, and that alpha oscillations would mainly correlate with increased blood oxygen-level dependent (BOLD) signal in the DMN. To test this hypothesis, we used high-density EEG (hdEEG) data simultaneously collected with fMRI scanning in 20 healthy volunteers at rest. We first detected the DMN from source reconstructed hdEEG data for multiple frequency bands, and we then mapped the correlation between temporal profile of hdEEG-derived DMN activity and fMRI–BOLD signals on a voxel-by-voxel basis. In line with our hypothesis, we found that the correlation map associated with alpha oscillations, more than with any other frequency bands, displayed a larger overlap with DMN regions. Overall, our study provided further evidence for a primary role of alpha oscillations in supporting DMN functioning. We suggest that simultaneous EEG–fMRI may represent a powerful tool to investigate the neurophysiological basis of human brain networks.

Human motor fatigability as evoked by repetitive movements results from a gradual breakdown of surround inhibition

Motor fatigability emerges when demanding tasks are executed over an extended period of time. Here, we used repetitive low-force movements that cause a gradual reduction in movement speed (or 'motor slowing') to study the central component of fatigability in healthy adults. We show that motor slowing is associated with a gradual increase of net excitability in the motor network and, specifically, in primary motor cortex (M1), which results from overall disinhibition. Importantly, we link performance decrements to a breakdown of surround inhibition in M1, which is associated with high coactivation of antagonistic muscle groups. This is consistent with the model that a loss of inhibitory control might broaden the tuning of population vectors such that movement patterns become more variable, ill-timed and effortful. We propose that the release of inhibition in M1 is an important mechanism underpinning motor fatigability and, potentially, also pathological fatigue as frequently observed in patients with brain disorders.

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.

Influence of Head Tissue Conductivity Uncertainties on EEG Dipole Reconstruction

Reliable EEG source analysis depends on sufficiently detailed and accurate head models. In this study, we investigate how uncertainties inherent to the experimentally determined conductivity values of the different conductive compartments influence the results of EEG source analysis. In a single source scenario, the superficial and focal somatosensory P20/N20 component, we analyze the influence of varying conductivities on dipole reconstructions using a generalized polynomial chaos (gPC) approach. We find that in particular the conductivity uncertainties for skin and skull have a significant influence on the EEG inverse solution, leading to variations in source localization by several centimeters. The conductivity uncertainties for gray and white matter were found to have little influence on the source localization, but a strong influence on the strength and orientation of the reconstructed source, respectively. As the CSF conductivity is most accurately determined of all conductivities in a realistic head model, CSF conductivity uncertainties had a negligible influence on the source reconstruction. This small uncertainty is a further benefit of distinguishing the CSF in realistic volume conductor models.