Reconstructing subcortical and cortical somatosensory activity via the RAMUS inverse source analysis technique using median nerve SEP data

Facilitation of Early and Middle Latency SEP after tDCS of M1: No Evidence of Primary Somatosensory Homeostatic Plasticity

Homeostatic plasticity is a mechanism that stabilizes cortical excitability within a physiological range. Most homeostatic plasticity protocols have primed and tested the homeostatic response of the primary motor cortex (M1). This study investigated if a homeostatic response could be recorded from the primary sensory cortex (S1) after inducing homeostatic plasticity in M1. In 31 healthy participants, homeostatic plasticity was induced over M1 with a priming and testing block of transcranial direct current stimulation (tDCS) in two different sessions (anodal and cathodal). S1 excitability was assessed by early (N20, P25) and middle-latency (N33-P45) somatosensory evoked potentials (SEP) extracted from 4 electrodes (CP5, CP3, P5, P3). Baseline and post-measures (post-priming, 0-min, 10-min, and 20-min after homeostatic induction) were taken. Anodal M1 homeostatic plasticity induction significantly facilitated the N20-P25, P45 peak, and N33-P45 early SEP components up to 20-min post-induction, without any indication of a homeostatic response (i.e., reduced SEP). Cathodal homeostatic induction did not induce any significant effect on early or middle latency SEPs. M1 homeostatic plasticity induction by anodal stimulation protocol to the primary motor cortex did not induce a homeostatic response in SEPs.

Further characterisation of late somatosensory evoked potentials using electroencephalogram and magnetoencephalogram source imaging

Beside the well-documented involvement of secondary somatosensory area, the cortical network underlying late somatosensory evoked potentials (P60/N60 and P100/N100) is still unknown. Electroencephalogram and magnetoencephalogram source imaging were performed to further investigate the origin of the brain cortical areas involved in late somatosensory evoked potentials, using sensory inputs of different strengths and by testing the correlation between cortical sources. Simultaneous high-density electroencephalograms and magnetoencephalograms were performed in 19 participants, and electrical stimulation was applied to the median nerve (wrist level) at intensity between 1.5 and 9 times the perceptual threshold. Source imaging was undertaken to map the stimulus-induced brain cortical activity according to each individual brain magnetic resonance imaging, during three windows of analysis covering early and late somatosensory evoked potentials. Results for P60/N60 and P100/N100 were compared with those for P20/N20 (early response). According to literature, maximal activity during P20/N20 was found in central sulcus contralateral to stimulation site. During P60/N60 and P100/N100, activity was observed in contralateral primary sensorimotor area, secondary somatosensory area (on both hemispheres) and premotor and multisensory associative cortices. Late responses exhibited similar characteristics but different from P20/N20, and no significant correlation was found between early and late generated activities. Specific clusters of cortical activities were activated with specific input/output relationships underlying early and late somatosensory evoked potentials. Cortical networks, partly common to and distinct from early somatosensory responses, contribute to late responses, all participating in the complex somatosensory brain processing.

Somatosensory evoked potentials recorded from DBS electrodes: the origin of subcortical N18

The different peaks of somatosensory-evoked potentials (SEP) originate from a variety of anatomical sites in the central nervous system. The origin of the median nerve subcortical N18 SEP has been studied under various conditions, but the exact site of its generation is still unclear. While it has been claimed to be located in the thalamic region, other studies indicated its possible origin below the pontomedullary junction. Here, we scrutinized and compared SEP recordings from median nerve stimulation through deep brain stimulation (DBS) electrodes implanted in various subcortical targets. We studied 24 patients with dystonia, Parkinson's disease, and chronic pain who underwent quadripolar electrode implantation for chronic DBS and recorded median nerve SEPs from globus pallidus internus (GPi), subthalamic nucleus (STN), thalamic ventral intermediate nucleus (Vim), and ventral posterolateral nucleus (VPL) and the centromedian-parafascicular complex (CM-Pf). The largest amplitude of the triphasic potential of the N18 complex was recorded in Vim. Bipolar recordings confirmed the origin to be close to Vim electrodes (and VPL/CM-Pf) and less close to STN electrodes. GPi recorded only far-field potentials in unipolar derivation. Recordings from DBS electrodes located in different subcortical areas allow determining the origin of certain subcortical SEP waves more precisely. The subcortical N18 of the median nerve SEP-to its largest extent-is generated ventral to the Vim in the region of the prelemniscal radiation/ zona incerta.

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

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

Modeling of blood flow in cerebral arterial circulation and its dynamic impact on electrical conductivity in a realistic multi-compartment head model

BACKGROUND AND OBJECTIVE

This study aims to assess the dynamic impact of non-Newtonian cerebral arterial circulation on electrical conductivity within a realistic multi-compartment head model. Evaluating this research question is crucial and challenging due to its relevance to electrophysiological modalities like transcranial electrical stimulation (tES), electro-/magnetoencephalography (EEG/MEG), and electrical impedance tomography (EIT). In these modalities, accurate forward modeling depends on the electrical conductivity, which is affected by complex tortuous vessel networks, limited data acquisition in Magnetic Resonance Imaging (MRI), and non-linear blood flow phenomena, including shear rate and viscosity in non-Newtonian fluid.

METHODS

To obtain an approximation for the blood concentration, we first use Navier-Stokes equations (NSEs) to solve for the pressure and velocity of the blood in the major vessels. Then Fick's law is used to solve for the blood concentration in the tissues. Finally, Archie's law is used to estimate the electrical conductivity distribution based on the blood concentration.

RESULTS

The results, obtained with an open 7 Tesla MRI dataset, suggest that a dynamic model of cerebral blood flow (CBF) for both arterial and microcirculation can be established; we find blood pressure and electrical conductivity distributions given a numerically simulated pulse sequence and approximate the blood concentration and electrical conductivity inside the brain based on those.

CONCLUSIONS

Our model provides an approximation of the dynamical blood flow and the corresponding electrical conductivity distribution in the different parts of the brain. The advantage of our approach is that it is applicable with limited a priori information about the blood flow and with an arbitrary head model distinguishing the arteries.

Pressure-Poisson equation in numerical simulation of cerebral arterial circulation and its effect on the electrical conductivity of the brain

BACKGROUND AND OBJECTIVE

This study considers dynamic modeling of the cerebral arterial circulation and reconstructing an atlas for the electrical conductivity of the brain. Electrical conductivity is a governing parameter in several electrophysiological modalities applied in neuroscience, such as electroencephalography (EEG), transcranial electrical stimulation (tES), and electrical impedance tomography (EIT). While high-resolution 7-Tesla (T) Magnetic Resonance Imaging (MRI) data allow for reconstructing the cerebral arteries with a cross-sectional diameter larger than the voxel size, electrical conductivity cannot be directly inferred from MRI data. Brain models of electrophysiology typically associate each brain tissue compartment with a constant electrical conductivity, omitting any dynamic effects of cerebral blood circulation. Incorporating those effects poses the challenge of solving a system of incompressible Navier-Stokes equations (NSEs) in a realistic multi-compartment head model. However, using a simplified circulation model is well-motivated since, on the one hand, the complete system does not always have a numerically stable solution and, on the other hand, the full set of arteries cannot be perfectly reconstructed from the MRI data, meaning that any solution will be approximative.

METHODS

We postulate that circulation in the distinguishable arteries can be estimated via the pressure-Poisson equation (PPE), which is coupled with Fick's law of diffusion for microcirculation. To establish a fluid exchange model between arteries and microarteries, a boundary condition derived from the Hagen-Poisseuille model is applied. The relationship between the estimated volumetric blood concentration and the electrical conductivity of the brain tissue is approximated through Archie's law for fluid flow in porous media.

RESULTS

Through the formulation of the PPE and a set of boundary conditions (BCs) based on the Hagen-Poisseuille model, we obtained an equivalent formulation of the incompressible Stokes equation (SE). Thus, allowing effective blood pressure estimation in cerebral arteries segmented from open 7T MRI data.

CONCLUSIONS

As a result of this research, we developed and built a useful modeling framework that accounts for the effects of dynamic blood flow on a novel MRI-based electrical conductivity atlas. The electrical conductivity perturbation obtained in numerical experiments has an appropriate overall match with previous studies on this subject. Further research to validate these results will be necessary.

Multi-compartment head modeling in EEG: Unstructured boundary-fitted tetra meshing with subcortical structures

This paper introduces an automated approach for generating a finite element (FE) discretization of a multi-compartment human head model for electroencephalographic (EEG) source localization. We aim to provide an adaptable FE mesh generation tool for EEG studies. Our technique relies on recursive solid angle labeling of a surface segmentation coupled with smoothing, refinement, inflation, and optimization procedures to enhance the mesh quality. In this study, we performed numerical meshing experiments with the three-layer Ary sphere and a magnetic resonance imaging (MRI)-based multi-compartment head segmentation which incorporates a comprehensive set of subcortical brain structures. These experiments are motivated, on one hand, by the sensitivity of non-invasive subcortical source localization to modeling errors and, on the other hand, by the present lack of open EEG software pipelines to discretize all these structures. Our approach was found to successfully produce an unstructured and boundary-fitted tetrahedral mesh with a sub-one-millimeter fitting error, providing the desired accuracy for the three-dimensional anatomical details, EEG lead field matrix, and source localization. The mesh generator applied in this study has been implemented in the open MATLAB-based Zeffiro Interface toolbox for forward and inverse processing in EEG and it allows for graphics processing unit acceleration.

An F-ratio-based method for estimating the number of active sources in MEG

Introduction

Magnetoencephalography (MEG) is a powerful technique for studying the human brain function. However, accurately estimating the number of sources that contribute to the MEG recordings remains a challenging problem due to the low signal-to-noise ratio (SNR), the presence of correlated sources, inaccuracies in head modeling, and variations in individual anatomy.

Methods

To address these issues, our study introduces a robust method for accurately estimating the number of active sources in the brain based on the F-ratio statistical approach, which allows for a comparison between a full model with a higher number of sources and a reduced model with fewer sources. Using this approach, we developed a formal statistical procedure that sequentially increases the number of sources in the multiple dipole localization problem until all sources are found.

Results

Our results revealed that the selection of thresholds plays a critical role in determining the method's overall performance, and appropriate thresholds needed to be adjusted for the number of sources and SNR levels, while they remained largely invariant to different inter-source correlations, translational modeling inaccuracies, and different cortical anatomies. By identifying optimal thresholds and validating our F-ratio-based method in simulated, real phantom, and human MEG data, we demonstrated the superiority of our F-ratio-based method over existing state-of-the-art statistical approaches, such as the Akaike Information Criterion (AIC) and Minimum Description Length (MDL).

Discussion

Overall, when tuned for optimal selection of thresholds, our method offers researchers a precise tool to estimate the true number of active brain sources and accurately model brain function.

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.

L1-norm vs. L2-norm fitting in optimizing focal multi-channel tES stimulation: linear and semidefinite programming vs. weighted least squares

BACKGROUND AND OBJECTIVE

This study focuses on Multi-Channel Transcranial Electrical Stimulation, a non-invasive brain method for stimulating neuronal activity under the influence of low-intensity currents. We introduce a mathematical formulation for finding a current pattern that optimizes an L1-norm fit between a given focal target distribution and volumetric current density inside the brain. L1-norm is well-known to favor well-localized or sparse distributions compared to L2-norm (least-squares) fitted estimates.

METHODS

We present a linear programming approach that performs L1-norm fitting and penalization of the current pattern (L1L1) to control the number of non-zero currents. The optimizer filters a large set of candidate solutions using a two-stage metaheuristic search from a pre-filtered set of candidates.

RESULTS

The numerical simulation results obtained with both 8- and 20-channel electrode montages suggest that our hypothesis on the benefits of L1-norm data fitting is valid. Compared to an L1-norm regularized L2-norm fitting (L1L2) via semidefinite programming and weighted Tikhonov least-squares method (TLS), the L1L1 results were overall preferable for maximizing the focused current density at the target position, and the ratio between focused and nuisance current magnitudes.

CONCLUSIONS

We propose the metaheuristic L1L1 optimization approach as a potential technique to obtain a well-localized stimulus with a controllable magnitude at a given target position. L1L1 finds a current pattern with a steep contrast between the anodal and cathodal electrodes while suppressing the nuisance currents in the brain, hence, providing a potential alternative to modulate the effects of the stimulation, e.g., the sensation experienced by the subject.

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.