Effects of forward model errors on EEG source localization

The Value of Normal Interictal EEGs in Epilepsy Diagnosis and Treatment Planning: A Retrospective Cohort Study using Population-level Spectral Power and Connectivity Patterns

Introduction

Scalp electroencephalography (EEG) is a cornerstone in the diagnosis and treatment of epilepsy, but routine EEG is often interpreted as normal without identification of epileptiform activity during expert visual review. The absence of interictal epileptiform activity on routine scalp EEGs can cause delays in receiving clinical treatment. These delays can be particularly problematic in the diagnosis and treatment of people with drug-resistant epilepsy (DRE) and those without structural abnormalities on MRI (i.e., MRI negative). Thus, there is a clinical need for alternative quantitative approaches that can inform diagnostic and treatment decisions when visual EEG review is inconclusive. In this study, we leverage a large population-level routine EEG database of people with and without focal epilepsy to investigate whether normal interictal EEG segments contain subtle deviations that could support the diagnosis of focal epilepsy.

Data & Methods

We identified multiple epochs representing eyes-closed wakefulness from 19-channel routine EEGs of a large and diverse neurological patient population (N=13,652 recordings, 12,134 unique patients). We then extracted the average spectral power and phase-lag-index-based connectivity within 1-45Hz of each EEG recording using these identified epochs. We decomposed the power spectral density and phase-based connectivity information of all the visually reviewed normal EEGs (N=6,242) using unsupervised tensor decompositions to extract dominant patterns of spectral power and scalp connectivity. We also identified an independent set of routine EEGs of a cohort of patients with focal epilepsy (N= 121) with various diagnostic classifications, including focal epilepsy origin (temporal, frontal), MRI (lesional, non-lesional), and response to anti-seizure medications (responsive vs. drug-resistant epilepsy). We analyzed visually normal interictal epochs from the EEGs using the power-spectral and phase-based connectivity patterns identified above and evaluated their potential in clinically relevant binary classifications.

Results

We obtained six patterns with distinct interpretable spatio-spectral signatures corresponding to putative aperiodic, oscillatory, and artifactual activity recorded on the EEG. The loadings for these patterns showed associations with patient age and expert-assigned grades of EEG abnormality. Further analysis using a physiologically relevant subset of these loadings differentiated patients with focal epilepsy from controls without history of focal epilepsy (mean AUC 0.78) but were unable to differentiate between frontal or temporal lobe epilepsy. In temporal lobe epilepsy, loadings of the power spectral patterns best differentiated drug-resistant epilepsy from drug-responsive epilepsy (mean AUC 0.73), as well as lesional epilepsy from non-lesional epilepsy (mean AUC 0.67), albeit with high variability across patients.

Significance

Our findings from a large population sample of EEGs suggest that normal interictal EEGs of patients with epilepsy contain subtle differences of predictive value that may improve the overall diagnostic yield of routine and prolonged EEGs. The presented approach for analyzing normal EEGs has the capacity to differentiate several diagnostic classifications of epilepsy, and can quantitatively characterize EEG activity in a scalable, expert-interpretable, and patient-specific fashion. Further technical development and clinical validation may yield normal EEG-derived computational biomarkers that could augment epilepsy diagnosis and assist clinical decision-making in the future.

On the concordance between electrical source imaging, anatomical and functional neuroimaging in patients with focal epilepsy

OBJECTIVE

Limited knowledge exists regarding how electrical source imaging (ESI) of interictal epileptiform discharges (IEDs) aligns with findings from other neuroimaging modalities. This study investigates the relationships of interictal ESI with MRI, 18FDG PET, SISCOM, and voxel-based morphometry (VBM) during presurgical evaluation of drug-resistant epilepsy (DRE).

METHODS

A cross-sectional study evaluated the concordance of IED locations from ESI using various inverse solutions (CLARA, LAURA, LORETA, SLORETA, SWLORETA, SSLOFO) with MRI lesions, 18FDG PET, SISCOM, and VBM grey matter abnormalities. The role of ESI in presurgical evaluation of DRE was assessed.

RESULTS

Significant relationships were identified between the localization and distribution of IEDs identified by ESI and the various sets of neuroimages. SLORETA and SWLORETA exhibited the highest concordance and interlobar associations with MRI, 18FDG PET and SISCOM. The main cluster of IEDs proved helpful in locating the epileptogenic zone (EZ).

CONCLUSIONS

The distribution of IEDs identified by the ESI technique exhibited a high degree of significant relationships with other neuroimaging sources. Its use may prove valuable in defining the epileptogenic zone.

SIGNIFICANCE

Combining ESI of IEDs with other neuroimaging techniques may be useful in the presurgical evaluation of DRE.

Localized estimation of event-related neural source activity from simultaneous MEG-EEG with a recurrent neural network

Estimating intracranial current sources underlying the electromagnetic signals observed from extracranial sensors is a perennial challenge in non-invasive neuroimaging. Established solutions to this inverse problem treat time samples independently without considering the temporal dynamics of event-related brain processes. This paper describes current source estimation from simultaneously recorded magneto- and electro-encephalography (MEEG) using a recurrent neural network (RNN) that learns sequential relationships from neural data. The RNN was trained in two phases: (1) pre-training and (2) transfer learning with L1 regularization applied to the source estimation layer. Performance of using scaled labels derived from MEEG, magnetoencephalography (MEG), or electroencephalography (EEG) were compared, as were results from volumetric source space with free dipole orientation and surface source space with fixed dipole orientation. Exact low-resolution electromagnetic tomography (eLORETA) and mixed-norm L1/L2 (MxNE) source estimation methods were also applied to these data for comparison with the RNN method. The RNN approach outperformed other methods in terms of output signal-to-noise ratio, correlation and mean-squared error metrics evaluated against reference event-related field (ERF) and event-related potential (ERP) waveforms. Using MEEG labels with fixed-orientation surface sources produced the most consistent estimates. To estimate sources of ERF and ERP waveforms, the RNN generates temporal dynamics within its internal computational units, driven by sequential structure in neural data used as training labels. It thus provides a data-driven model of computational transformations from psychophysiological events into corresponding event-related neural signals, which is unique among MEEG source reconstruction solutions.

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.

Investigating the effect of template head models on Event-Related Potential source localization: a simulation and real-data study

Introduction

Event-Related Potentials (ERPs) are valuable for studying brain activity with millisecond-level temporal resolution. While the temporal resolution of this technique is excellent, the spatial resolution is limited. Source localization aims to identify the brain regions generating the EEG data, thus increasing the spatial resolution, but its accuracy depends heavily on the head model used. This study compares the performance of subject-specific and template-based head models in both simulated and real-world ERP localization tasks.

Methods

Simulated data mimicking realistic ERPs was created to evaluate the impact of head model choice systematically, after which subject-specific and template-based head models were used for the reconstruction of the data. The different modeling approaches were also applied to a face recognition dataset.

Results

The results indicate that the template models capture the simulated activity less accurately, producing more spurious sources and identifying less true sources correctly. Furthermore, the results show that while creating more accurate and detailed head models is beneficial for the localization accuracy when using subject-specific head models, this is less the case for template head models. The main N170 source of the face recognition dataset was correctly localized to the fusiform gyrus, a known face processing area, using the subject-specific models. Apart from the fusiform gyrus, the template models also reconstructed several other sources, illustrating the localization inaccuracies.

Discussion

While template models allow researchers to investigate the neural generators of ERP components when no subject-specific MRIs are available, it could lead to misinterpretations. Therefore, it is important to consider a priori knowledge and hypotheses when interpreting results obtained with template head models, acknowledging potential localization errors.

Semi-analytic three-shell forward calculation for magnetoencephalography

In previous studies, the magnetic lead field theorem in the quasi-static approximation was derived and used for the development of a method for the forward problem of MEG. It was applied and tested on a single-shell model of the human head and the question whether one shell is adequate enough for the calculation of the magnetic field is the main reason for this study. This forward method is based on the fundamental concept that one can calculate the lead field for MEG by decomposing it into two parts: the lead field of an arbitrary volume conductor that is already known and the gradient of basis functions that have to be harmonic, here derived from spherical harmonics. The problem then is reduced to evaluating the coefficients found in the basis functions. In this research we aim to improve the accuracy of the forward model, hence improving the localization accuracy in inverse methods by introducing a more detailed realistic head model. We here generalize the algorithm developed for a single-shell volume conductor to a three-shell volume conductor representing the brain, the skull and the skin with homogenous and isotropic conductivities in realistic ratios. The expansion to three shells could be tested as the three-shell algorithm is approaching the single-shell with high accuracy in special cases where three-shell solutions can also be calculated using a single-shell solution, especially for higher levels of expansion. The deviation in the calculation of the lead field is also evaluated when using three shells with realistic conductivities. The magnetic field turned out to differ to an important measurable extend in particular for deeper sources, making the three-shell algorithm substantially more accurate for these dipole locations.

Improved Dipole Source Localization from Simultaneous MEG-EEG Data by Combining a Global Optimization Algorithm with a Local Parameter Search: A Brain Phantom Study

Dipole localization, a fundamental challenge in electromagnetic source imaging, inherently constitutes an optimization problem aimed at solving the inverse problem of electric current source estimation within the human brain. The accuracy of dipole localization algorithms is contingent upon the complexity of the forward model, often referred to as the head model, and the signal-to-noise ratio (SNR) of measurements. In scenarios characterized by low SNR, often corresponding to deep-seated sources, existing optimization techniques struggle to converge to global minima, thereby leading to the localization of dipoles at erroneous positions, far from their true locations. This study presents a novel hybrid algorithm that combines simulated annealing with the traditional quasi-Newton optimization method, tailored to address the inherent limitations of dipole localization under low-SNR conditions. Using a realistic head model for both electroencephalography (EEG) and magnetoencephalography (MEG), it is demonstrated that this novel hybrid algorithm enables significant improvements of up to 45% in dipole localization accuracy compared to the often-used dipole scanning and gradient descent techniques. Localization improvements are not only found for single dipoles but also in two-dipole-source scenarios, where sources are proximal to each other. The novel methodology presented in this work could be useful in various applications of clinical neuroimaging, particularly in cases where recordings are noisy or sources are located deep within the brain.

Mobile neuroimaging: What we have learned about the neural control of human walking, with an emphasis on EEG-based research

Our understanding of the neural control of human walking has changed significantly over the last twenty years and mobile brain imaging methods have contributed substantially to current knowledge. High-density electroencephalography (EEG) has the advantages of being lightweight and mobile while providing temporal resolution of brain changes within a gait cycle. Advances in EEG hardware and processing methods have led to a proliferation of research on the neural control of locomotion in neurologically intact adults. We provide a narrative review of the advantages and disadvantages of different mobile brain imaging methods, then summarize findings from mobile EEG studies quantifying electrocortical activity during human walking. Contrary to historical views on the neural control of locomotion, recent studies highlight the widespread involvement of many areas, such as the anterior cingulate, posterior parietal, prefrontal, premotor, sensorimotor, supplementary motor, and occipital cortices, that show active fluctuations in electrical power during walking. The electrocortical activity changes with speed, stability, perturbations, and gait adaptation. We end with a discussion on the next steps in mobile EEG research.

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

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

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

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

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

Objective

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

Methods

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

Results

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

Conclusions

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

Significance

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

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

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

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

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.

Electrical brain activity during human walking with parametric variations in terrain unevenness and walking speed

Mobile brain imaging with high-density electroencephalography (EEG) can provide insight into the cortical processes involved in complex human walking tasks. While uneven terrain is common in the natural environment and poses challenges to human balance control, there is limited understanding of the supraspinal processes involved with traversing uneven terrain. The primary objective of this study was to quantify electrocortical activity related to parametric variations in terrain unevenness for neurotypical young adults. We used high-density EEG to measure brain activity when 32 young adults walked on a novel custom-made uneven terrain treadmill surface with four levels of difficulty at a walking speed tailored to each participant. We identified multiple brain regions associated with uneven terrain walking. Alpha (8 - 13 Hz) and beta (13 - 30 Hz) spectral power decreased in the sensorimotor and posterior parietal areas with increasing terrain unevenness while theta (4 - 8 Hz) power increased in the mid/posterior cingulate area with terrain unevenness. We also found that within stride spectral power fluctuations increased with terrain unevenness. Our secondary goal was to investigate the effect of parametric changes in walking speed (0.25 m/s, 0.5 m/s, 0.75 m/s, 1.0 m/s) to differentiate the effects of walking speed from uneven terrain. Our results revealed that electrocortical activities only changed substantially with speed within the sensorimotor area but not in other brain areas. Together, these results indicate there are distinct cortical processes contributing to the control of walking over uneven terrain versus modulation of walking speed on smooth, flat terrain. Our findings increase our understanding of cortical involvement in an ecologically valid walking task and could serve as a benchmark for identifying deficits in cortical dynamics that occur in people with mobility deficits.

Finding tau rhythms in EEG: An independent component analysis approach

Tau rhythms are largely defined by sound responsive alpha band (~8-13 Hz) oscillations generated largely within auditory areas of the superior temporal gyri. Studies of tau have mostly employed magnetoencephalography or intracranial recording because of tau's elusiveness in the electroencephalogram. Here, we demonstrate that independent component analysis (ICA) decomposition can be an effective way to identify tau sources and study tau source activities in EEG recordings. Subjects (N = 18) were passively exposed to complex acoustic stimuli while the EEG was recorded from 68 electrodes across the scalp. Subjects' data were split into 60 parallel processing pipelines entailing use of five levels of high-pass filtering (passbands of 0.1, 0.5, 1, 2, and 4 Hz), three levels of low-pass filtering (25, 50, and 100 Hz), and four different ICA algorithms (fastICA, infomax, adaptive mixture ICA [AMICA], and multi-model AMICA [mAMICA]). Tau-related independent component (IC) processes were identified from this data as being localized near the superior temporal gyri with a spectral peak in the 8-13 Hz alpha band. These "tau ICs" showed alpha suppression during sound presentations that was not seen for other commonly observed IC clusters with spectral peaks in the alpha range (e.g., those associated with somatomotor mu, and parietal or occipital alpha). The choice of analysis parameters impacted the likelihood of obtaining tau ICs from an ICA decomposition. Lower cutoff frequencies for high-pass filtering resulted in significantly fewer subjects showing a tau IC than more aggressive high-pass filtering. Decomposition using the fastICA algorithm performed the poorest in this regard, while mAMICA performed best. The best combination of filters and ICA model choice was able to identify at least one tau IC in the data of ~94% of the sample. Altogether, the data reveal close similarities between tau EEG IC dynamics and tau dynamics observed in MEG and intracranial data. Use of relatively aggressive high-pass filters and mAMICA decomposition should allow researchers to identify and characterize tau rhythms in a majority of their subjects. We believe adopting the ICA decomposition approach to EEG analysis can increase the rate and range of discoveries related to auditory responsive tau rhythms.

Two-Stage Atomic Decomposition of Multichannel EEG and the Previously Undetectable Sleep Spindles

We propose a two-step procedure for atomic decomposition of multichannel EEGs, based upon multivariate matching pursuit and dipolar inverse solution, from which atoms representing relevant EEG structures are selected according to prior knowledge. We detect sleep spindles in 147 polysomnographic recordings from the Montreal Archive of Sleep Studies. Detection is compared with human scorers and two state-of-the-art algorithms, which find only about a third of the structures conforming to the definition of sleep spindles and detected by the proposed method. We provide arguments supporting the thesis that the previously undetectable sleep spindles share the same properties as those marked by human experts and previously applied methods, and were previously omitted only because of unfavorable local signal-to-noise ratios, obscuring their visibility to both human experts and algorithms replicating their markings. All detected EEG structures are automatically parametrized by their time and frequency centers, width duration, phase, and spatial location of an equivalent dipolar source within the brain. It allowed us, for the first time, to estimate the spatial gradient of sleep spindles frequencies, which not only confirmed quantitatively the well-known prevalence of higher frequencies in posterior regions, but also revealed a significant gradient in the sagittal plane. The software used in this study is freely available.

First activity and interactions in thalamus and cortex using raw single-trial EEG and MEG elicited by somatosensory stimulation

Introduction

One of the primary motivations for studying the human brain is to comprehend how external sensory input is processed and ultimately perceived by the brain. A good understanding of these processes can promote the identification of biomarkers for the diagnosis of various neurological disorders; it can also provide ways of evaluating therapeutic techniques. In this work, we seek the minimal requirements for identifying key stages of activity in the brain elicited by median nerve stimulation.

Methods

We have used a priori knowledge and applied a simple, linear, spatial filter on the electroencephalography and magnetoencephalography signals to identify the early responses in the thalamus and cortex evoked by short electrical stimulation of the median nerve at the wrist. The spatial filter is defined first from the average EEG and MEG signals and then refined using consistency selection rules across ST. The refined spatial filter is then applied to extract the timecourses of each ST in each targeted generator. These ST timecourses are studied through clustering to quantify the ST variability. The nature of ST connectivity between thalamic and cortical generators is then studied within each identified cluster using linear and non-linear algorithms with time delays to extract linked and directional activities. A novel combination of linear and non-linear methods provides in addition discrimination of influences as excitatory or inhibitory.

Results

Our method identifies two key aspects of the evoked response. Firstly, the early onset of activity in the thalamus and the somatosensory cortex, known as the P14 and P20 in EEG and the second M20 for MEG. Secondly, good estimates are obtained for the early timecourse of activity from these two areas. The results confirm the existence of variability in ST brain activations and reveal distinct and novel patterns of connectivity in different clusters.

Discussion

It has been demonstrated that we can extract new insights into stimulus processing without the use of computationally costly source reconstruction techniques which require assumptions and detailed modeling of the brain. Our methodology, thanks to its simplicity and minimal computational requirements, has the potential for real-time applications such as in neurofeedback systems and brain-computer interfaces.

EEG Source Imaging-Clinical Considerations for EEG Acquisition and Signal Processing for Improved Temporo-Spatial Resolution

SUMMARY

EEG source imaging (ESI) has gained traction in recent years as a useful clinical tool for the noninvasive surgical work-up of patients with drug-resistant focal epilepsy. Despite its proven benefits for the temporo-spatial modeling of spike and seizure sources, ESI remains widely underused in clinical practice. This partly relates to a lack of clarity around an optimal approach to the acquisition and processing of scalp EEG data for the purpose of ESI. Here, we describe some of the practical considerations for the clinical application of ESI. We focus on patient preparation, the impact of electrode number and distribution across the scalp, the benefit of averaging raw data for signal analysis, and the relevance of modeling different phases of the interictal discharge as it evolves from take-off to peak. We emphasize the importance of recording high signal-to-noise ratio data for reliable source analysis. We argue that the accuracy of modeling cortical sources can be improved using higher electrode counts that include an inferior temporal array, by averaging interictal waveforms rather than limiting ESI to single spike analysis, and by careful interrogation of earlier phase components of these waveforms. No amount of postacquisition signal processing or source modeling sophistication, however, can make up for suboptimally recorded scalp EEG data in a poorly prepared patient.

EEG electrode localization with 3D iPhone scanning using point-cloud electrode selection (PC-ES)

Objective.Electroencephalography source imaging (ESI) is a valuable tool in clinical evaluation for epilepsy patients but is underutilized in part due to sensitivity to anatomical modeling errors. Accurate localization of scalp electrodes is instrumental to ESI, but existing localization devices are expensive and not portable. As a result, electrode localization challenges further impede access to ESI, particularly in inpatient and intensive care settings.Approach.To address this challenge, we present a portable and affordable electrode digitization method using the 3D scanning feature in modern iPhone models. This technique combines iPhone scanning with semi-automated image processing using point-cloud electrode selection (PC-ES), a custom MATLAB desktop application. We compare iPhone electrode localization to state-of-the-art photogrammetry technology in a human study with over 6000 electrodes labeled using each method. We also characterize the performance of PC-ES with respect to head location and examine the relative impact of different algorithm parameters.Main Results.The median electrode position variation across reviewers was 1.50 mm for PC-ES scanning and 0.53 mm for photogrammetry, and the average median distance between PC-ES and photogrammetry electrodes was 3.4 mm. These metrics demonstrate comparable performance of iPhone/PC-ES scanning to currently available technology and sufficient accuracy for ESI.Significance.Low cost, portable electrode localization using iPhone scanning removes barriers to ESI in inpatient, outpatient, and remote care settings. While PC-ES has current limitations in user bias and processing time, we anticipate these will improve with software automation techniques as well as future developments in iPhone 3D scanning technology.

Changes in EEG alpha-band power during prehension indicates neural motor drive inhibition

Changes in alpha band activity (8-12 Hz) indicate the downregulation of brain regions during cognitive tasks, reflecting real-time cognitive load. Despite this, its feasibility to be used in a more dynamic environment with ongoing motor corrections has not been studied. This research used electroencephalography (EEG) to explore how different brain regions are engaged during a simple grasp and lift task where unexpected changes to the object's weight or surface friction are introduced. The results suggest that alpha activity changes related to motor error correction occur only in motor-related areas (i.e. central areas) but not in error processing areas (i.e., frontoparietal network) during unexpected weight changes. This suggests that oscillations over motor areas reflect the reduction of motor drive related to motor error correction, thus, being a potential cortical electrophysiological biomarker for the process and not solely as a proxy for cognitive demands. This observation is particularly relevant in scenarios where these signals are used to evaluate high cognitive demands co-occurring with high levels of motor errors and corrections, such as prosthesis use. The establishment of electrophysiological biomarkers of mental resource allocation during movement and cognition can help identify indicators of mental workload and motor drive, which may be useful for improving brain-machine interfaces.NEW & NOTEWORTHY We demonstrated that alpha suppression, an EEG phenomenon with high temporal resolution, occurs over the primary sensorimotor area during error correction during lift movements. Interpretations of alpha activity are often attributed to high cognitive demands, thus recognizing that it is also influenced by motor processes is important in situations where cognitive demands are paired with movement errors. This could further have application as a biomarker for error correction in human-machine interfaces, such as neuroprostheses.

iCanClean Removes Motion, Muscle, Eye, and Line-Noise Artifacts from Phantom EEG

The goal of this study was to test a novel approach (iCanClean) to remove non-brain sources from scalp EEG data recorded in mobile conditions. We created an electrically conductive phantom head with 10 brain sources, 10 contaminating sources, scalp, and hair. We tested the ability of iCanClean to remove artifacts while preserving brain activity under six conditions: Brain, Brain + Eyes, Brain + Neck Muscles, Brain + Facial Muscles, Brain + Walking Motion, and Brain + All Artifacts. We compared iCanClean to three other methods: Artifact Subspace Reconstruction (ASR), Auto-CCA, and Adaptive Filtering. Before and after cleaning, we calculated a Data Quality Score (0-100%), based on the average correlation between brain sources and EEG channels. iCanClean consistently outperformed the other three methods, regardless of the type or number of artifacts present. The most striking result was for the condition with all artifacts simultaneously present. Starting from a Data Quality Score of 15.7% (before cleaning), the Brain + All Artifacts condition improved to 55.9% after iCanClean. Meanwhile, it only improved to 27.6%, 27.2%, and 32.9% after ASR, Auto-CCA, and Adaptive Filtering. For context, the Brain condition scored 57.2% without cleaning (reasonable target). We conclude that iCanClean offers the ability to clear multiple artifact sources in real time and could facilitate human mobile brain-imaging studies with EEG.

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.

Effortful listening produces both enhancement and suppression of alpha in the EEG

Introduction

Adverse listening conditions can drive increased mental effort during listening. Neuromagnetic alpha oscillations (8-13 Hz) may index this listening effort, but inconsistencies regarding the direction of the relationship are abundant. We performed source analyses on high-density EEG data collected during a speech-on-speech listening task to address the possibility that opposing alpha power relationships among alpha producing brain sources drive this inconsistency.

Methods

Listeners (N=20) heard two simultaneously presented sentences of the form: Ready go to now. They either reported the color/number pair of a "Baron" call sign sentence (active: high effort), or ignored the stimuli (passive: low effort). Independent component analysis (ICA) was used to segregate temporally distinct sources in the EEG.

Results

Analysis of independent components (ICs) revealed simultaneous alpha enhancements (e.g., for somatomotor mu ICs) and suppressions (e.g., for left temporal ICs) for different brain sources. The active condition exhibited stronger enhancement for left somatomotor mu rhythm ICs, but stronger suppression for central occipital ICs.

Discussion

This study shows both alpha enhancement and suppression to be associated with increases in listening effort. Literature inconsistencies could partially relate to some source activities overwhelming others in scalp recordings.

Ultrasound measurements of interactive turn-taking in question-answer sequences: Articulatory preparation is delayed but not tied to the response

We know that speech planning in conversational turn-taking can happen in overlap with the previous turn and research suggests that it starts as early as possible, that is, as soon as the gist of the previous turn becomes clear. The present study aimed to investigate whether planning proceeds all the way up to the last stage of articulatory preparation (i.e., putting the articulators in place for the first phoneme of the response) and what the timing of this process is. Participants answered pre-recorded quiz questions (being under the illusion that they were asked live), while their tongue movements were measured using ultrasound. Planning could start early for some quiz questions (i.e., midway during the question), but late for others (i.e., only at the end of the question). The results showed no evidence for a difference between tongue movements in these two types of questions for at least two seconds after planning could start in early-planning questions, suggesting that speech planning in overlap with the current turn proceeds more slowly than in the clear. On the other hand, when time-locking to speech onset, tongue movements differed between the two conditions from up to two seconds before this point. This suggests that articulatory preparation can occur in advance and is not fully tied to the overt response itself.

Feasibility of Using Source-Level Brain Computer Interface for People with Multiple Sclerosis

This work evaluates the feasibility of using a source level Brain-Computer Interface (BCI) for people with Multiple Sclerosis (MS). Data used was previously collected EEG of eight participants (one participant with MS and seven neurotypical participants) who performed imagined movement of the right and left hand. Equivalent current dipole cluster fitting was used to assess related brain activity at the source level and assessed using dipole location and power spectrum analysis. Dipole clusters were resolved within the motor cortices with some notable spatial difference between the MS and control participants. Neural sources that generate motor imagery originated from similar motor areas in the participant with MS compared to the neurotypical participants. Power spectral analysis indicated a reduced level of alpha power in the participant with MS during imagery tasks compared to neurotypical participants. Power in the beta band may be used to distinguish between left and right imagined movement for users with MS in BCI applications.Clinical Relevance- This paper demonstrates the cortical areas activated during imagined BCI-type tasks in a participant with Multiple Sclerosis (MS), and is a proof of concept for translating BCI research to potential users with MS.

Evaluating the Influence of Anatomical Accuracy and Electrode Positions on EEG Forward Solutions

Generating realistic volume conductor models for forward calculations in electroencephalography (EEG) is not trivial and several factors contribute to the accuracy of such models, two of which are its anatomical accuracy and the accuracy with which electrode positions are known. Here, we investigate effects of anatomical accuracy by comparing forward solutions from SimNIBS, a tool which allows state-of-the-art anatomical modeling, with well-established pipelines in MNE-Python and FieldTrip. We also compare different ways of specifying electrode locations when digitized positions are not available such as transformation of measured positions from standard space and transformation of a manufacturer layout. Substantial effects of anatomical accuracy were seen throughout the entire brain both in terms of field topography and magnitude with SimNIBS generally being more accurate than the pipelines in MNE-Python and FieldTrip. Topographic and magnitude effects were particularly pronounced for MNE-Python which uses a three-layer boundary element method (BEM) model. We attribute these mainly to the coarse representation of the anatomy used in this model, in particular differences in skull and cerebrospinal fluid (CSF). Effects of electrode specification method were evident in occipital and posterior areas when using a transformed manufacturer layout whereas transforming measured positions from standard space generally resulted in smaller errors. We suggest modeling the anatomy of the volume conductor as accurately possible and we hope to facilitate this by making it easy to export simulations from SimNIBS to MNE-Python and FieldTrip for further analysis. Likewise, if digitized electrode positions are not available, a set of measured positions on a standard head template may be preferable to those specified by the manufacturer.

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

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

Automatic Quality Control of the numerical accuracy of EEG Lead fields

Precise individualized EEG source localization is predicated on having accurate subject-specific Lead Fields (LFs) obtained from their Magnetic Resonance Images (MRI). LF calculation is a complex process involving several error-prone steps that start with obtaining a realistic head model from the MRI and finalizing with computationally expensive solvers such as the Boundary Element Method (BEM) or Finite Element Method (FEM). Current Big-Data applications require the calculation of batches of hundreds or thousands of LFs. LF. Quality Control is conventionally checked subjectively by experts, a procedure not feasible in practice for larger batches. To facilitate this step, we introduce the Lead Field Automatic-Quality Control Index (LF-AQI) that flags LF with potential errors. We base our LF-AQI on the assumption that LFs obtained from simpler head models, i.e., the homogeneous head model LF (HHM-LF) or spherical head model LF (SHM-LF), deviate only moderately from a "good" realistic test LF. Since these simpler LFs are easier to compute and check for errors, they may serve as "reference LF" to detect anomalous realistic test LF. We investigated this assumption by comparing correlation-based channel ρ_min(ref,test)and source τ_min(ref,test) similarity indices (SI) between "gold standards," i.e., very accurate FEM and BEM LFs, and the proposed references (HHM-LF and SHM-LF). Surprisingly we found that the most uncomplicated possible reference, HHM-LF had high SI values with the gold standards-leading us to explore further use of the channel ρ_min(HHM-LF,test)and source τ_min(HHM-LF,test)_{SI as a basis for} our LF-AQI. Indeed, these SI successfully detected five simulated scenarios of LFs artifacts. This result encouraged us to evaluate the SI on a large dataset and thus define our LF-AQI. We thus computed the SI of 1251 LFs obtained from the Child Mind Institute (CMI) MRI dataset. When ρ_min(HHM-LF,test)and source τ_min(HHM-LF,test) were plotted for all test subjects on a 2D space, most were tightly clustered around the median of a high similarity centroid (HSC), except for a smaller proportion of outliers. We define the LF-AQI for a given LF as the log Euclidean distance between its SI and the HSC median. To automatically detect outliers, the threshold is at the 90^th percentile of the CMI LF-AQIs (-0.9755). LF-AQI greater than this threshold flag individual LF to be checked. The robustness of this LF-AQI screening was checked by repeated out-of-sample validation. Strikingly, minor corrections in re-processing the flagged cases eliminated their status as outliers. Furthermore, the "doubtful" labels assigned by LF-AQI were validated by neuroscience students using a Likert scale questionnaire designed to manually check the LF's quality. Item Response Theory (IRT) analysis was applied to the questionnaire results to compute an optimized model and a latent variable θ for that model. A linear mixed model (LMM) between the θ and LF-AQI resulted in an effect with a Cohen's d value of 1.3 and a p-value <0.001, thus validating the correspondence of LF-AQI with the visual quality control. We provide an open-source pipeline to implement both LF calculation and its quality control to allow further evaluation of our index.

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

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

Multiple Dipole Source Position and Orientation Estimation Using Non-Invasive EEG-like Signals

The problem of precisely estimating the position and orientation of multiple dipoles using synthetic EEG signals is considered in this paper. After determining a proper forward model, a nonlinear constrained optimization problem with regularization is solved, and the results are compared with a widely used research code, namely EEGLAB. A thorough sensitivity analysis of the estimation algorithm to the parameters (such as the number of samples and sensors) in the assumed signal measurement model is conducted. To confirm the efficacy of the proposed source identification algorithm on any category of data sets, three different kinds of data-synthetic model data, visually evoked clinical EEG data, and seizure clinical EEG data are used. Furthermore, the algorithm is tested on both the spherical head model and the realistic head model based on the MNI coordinates. The numerical results and comparisons with the EEGLAB show very good agreement, with little pre-processing required for the acquired data.

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

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

EEG cortical activity and connectivity correlates of early sympathetic response during cold pressor test

Previous studies have identified several brain regions involved in the sympathetic response and its integration with pain, cognition, emotions and memory processes. However, little is known about how such regions dynamically interact during a sympathetic activation task. In this study, we analyzed EEG activity and effective connectivity during a cold pressor test (CPT). A source localization analysis identified a network of common active sources including the right precuneus (r-PCu), right and left precentral gyri (r-PCG, l-PCG), left premotor cortex (l-PMC) and left anterior cingulate cortex (l-ACC). We comprehensively analyzed the network dynamics by estimating power variation and causal interactions among the network regions through the direct directed transfer function (dDTF). A connectivity pattern dominated by interactions in [Formula: see text] (8-12) Hz band was observed in the resting state, with r-PCu acting as the main hub of information flow. After the CPT onset, we observed an abrupt suppression of such [Formula: see text]-band interactions, followed by a partial recovery towards the end of the task. On the other hand, an increase of [Formula: see text]-band (1-4) Hz interactions characterized the first part of CPT task. These results provide novel information on the brain dynamics induced by sympathetic stimuli. Our findings suggest that the observed suppression of [Formula: see text] and rise of [Formula: see text] dynamical interactions could reflect non-pain-specific arousal and attention-related response linked to stimulus' salience.

Estimating neural activity from visual areas using functionally defined EEG templates

Electroencephalography (EEG) is a common and inexpensive method to record neural activity in humans. However, it lacks spatial resolution making it difficult to determine which areas of the brain are responsible for the observed EEG response. Here we present a new easy-to-use method that relies on EEG topographical templates. Using MRI and fMRI scans of 50 participants, we simulated how the activity in each visual area appears on the scalp and averaged this signal to produce functionally defined EEG templates. Once created, these templates can be used to estimate how much each visual area contributes to the observed EEG activity. We tested this method on extensive simulations and on real data. The proposed procedure is as good as bespoke individual source localization methods, robust to a wide range of factors, and has several strengths. First, because it does not rely on individual brain scans, it is inexpensive and can be used on any EEG data set, past or present. Second, the results are readily interpretable in terms of functional brain regions and can be compared across neuroimaging techniques. Finally, this method is easy to understand, simple to use and expandable to other brain sources.

Ideal current dipoles are appropriate source representations for simulating neurons for intracranial recordings

OBJECTIVE

To determine whether dipoles are an appropriate simplified representation of neural sources for stereo-EEG (sEEG).

METHODS

We compared the distributions of voltages generated by a dipole, biophysically realistic cortical neuron models, and extended regions of cortex to determine how well a dipole represented neural sources at different spatial scales and at electrode to neuron distances relevant for sEEG. We also quantified errors introduced by the dipole approximation of neural sources in sEEG source localization using standardized low-resolution electrotomography (sLORETA).

RESULTS

For pyramidal neurons, the coefficient of correlation between voltages generated by a dipole and neuron model were > 0.9 for distances > 1 mm. For small regions of cortex (∼0.1 cm2), the error in voltages between a dipole and region was < 100 µV for all distances. However, larger regions of active cortex (>5 cm2) yielded > 50 µV errors within 1.5 cm of an electrode when compared to single dipoles. Finally, source localization errors were < 5 mm when using dipoles to represent realistic neural sources.

CONCLUSIONS

Single dipoles are an appropriate source model to represent both single neurons and small regions of active cortex, while multiple dipoles are required to represent large regions of cortex.

SIGNIFICANCE

Dipoles are computationally tractable and valid source models for sEEG.

Methylphenidate, Guanfacine, and Combined Treatment Effects on Electroencephalography Correlates of Spatial Working Memory in Attention-Deficit/Hyperactivity Disorder

OBJECTIVE

The combination of d-methylphenidate and guanfacine (an α-2A adrenergic agonist) may be an effective alternative to either agent as monotherapy in children with attention-deficit/hyperactivity disorder (ADHD). This study investigated the neural mechanisms underlying medication effects using cortical source analysis of electroencephalography (EEG) data.

METHOD

A total of 172 children with ADHD (aged 7-14; 118 boys) completed an 8-week randomized, double-blind, comparative study with 3 treatment arms: d-methylphenidate, guanfacine, or their combination. EEG modulations of brain oscillations at baseline and end point were measured during a spatial working memory task from cortical sources localized within the anterior cingulate (midfrontal) and primary visual cortex (midoccipital), based on previously reported ADHD and control differences. Linear mixed models examined treatment effects on EEG and performance measures.

RESULTS

Combined treatment decreased midoccipital EEG power across most frequency bands and task phases. Several midoccipital EEG measures also showed significantly greater changes with combined treatment than with monotherapies. D-methylphenidate significantly increased midoccipital theta during retrieval, while guanfacine produced only trend-level reductions in midoccipital alpha during maintenance and retrieval. Task accuracy improved with combined treatment, was unchanged with d-methylphenidate, and worsened with guanfacine. Treatment-related changes in midoccipital power correlated with improvement in ADHD severity.

CONCLUSION

These findings show that combined treatment ameliorates midoccipital neural activity associated with treatment-related behavioral improvements and previously implicated in visuo-attentional deficits in ADHD. Both monotherapies had limited effects on EEG measures, with guanfacine further showing detrimental effects on performance. The identified midoccipital EEG profile may aid future treatment monitoring for children with ADHD.

CLINICAL TRIAL REGISTRATION INFORMATION

Single Versus Combination Medication Treatment for Children With Attention Deficit Hyperactivity Disorder (Project1); https://clinicaltrials.gov/; NCT00429273.

DIVERSITY & INCLUSION STATEMENT

We worked to ensure race, ethnic, and/or other types of diversity in the recruitment of human participants. We worked to ensure sex and gender balance in the recruitment of human participants. One or more of the authors of this paper self-identifies as a member of one or more historically underrepresented racial and/or ethnic groups in science. While citing references scientifically relevant for this work, we also actively worked to promote sex and gender balance in our reference list. We actively worked to promote inclusion of historically underrepresented racial and/or ethnic groups in science in our author group. We actively worked to promote sex and gender balance in our author group.

Closed-loop modulation of local slow oscillations in human NREM sleep

Slow-wave sleep is the deep non-rapid eye-movement (NREM) sleep stage that is most relevant for the recuperative function of sleep. Its defining property is the presence of slow oscillations (<2 Hz) in the scalp electroencephalogram (EEG). Slow oscillations are generated by a synchronous back and forth between highly active UP-states and silent DOWN-states in neocortical neurons. Growing evidence suggests that closed-loop sensory stimulation targeted at UP-states of EEG-defined slow oscillations can enhance the slow oscillatory activity, increase sleep depth, and boost sleep's recuperative functions. However, several studies failed to replicate such findings. Failed replications might be due to the use of conventional closed-loop stimulation algorithms that analyze the signal from one single electrode and thereby neglect the fact that slow oscillations vary with respect to their origins, distributions, and trajectories on the scalp. In particular, conventional algorithms nonspecifically target functionally heterogeneous UP-states of distinct origins. After all, slow oscillations at distinct sites of the scalp have been associated with distinct functions. Here we present a novel EEG-based closed-loop stimulation algorithm that allows targeting UP- and DOWN-states of distinct cerebral origins based on topographic analyses of the EEG: the topographic targeting of slow oscillations (TOPOSO) algorithm. We present evidence that the TOPOSO algorithm can detect and target local slow oscillations with specific, predefined voltage maps on the scalp in real-time. When compared to a more conventional, single-channel-based approach, TOPOSO leads to fewer but locally more specific stimulations in a simulation study. In a validation study with napping participants, TOPOSO targets auditory stimulation reliably at local UP-states over frontal, sensorimotor, and centro-parietal regions. Importantly, auditory stimulation temporarily enhanced the targeted local state. However, stimulation then elicited a standard frontal slow oscillation rather than local slow oscillations. The TOPOSO algorithm is suitable for the modulation and the study of the functions of local slow oscillations.

Lateralization of orthographic processing in fixed-gaze and natural reading conditions

Lateralized processing of orthographic information is a hallmark of proficient reading. However, how this finding obtained for fixed-gaze processing of orthographic stimuli translates to ecologically valid reading conditions remained to be clarified. To address this shortcoming, here we assessed the lateralization of early orthographic processing in fixed-gaze and natural reading conditions using concurrent eye-tracking and EEG data recorded from young adults without reading difficulties. Sensor-space analyses confirmed the well-known left-lateralized negative-going deflection of fixed-gaze EEG activity throughout the period of early orthographic processing. At the same time, fixation-related EEG activity exhibited left-lateralized followed by right-lateralized processing of text stimuli during natural reading. A strong positive relationship was found between the early leftward lateralization in fixed-gaze and natural reading conditions. Using source-space analyses, early left-lateralized brain activity was obtained in lateraloccipital and posterior ventral occipito-temporal cortices reflecting letter-level processing in both conditions. In addition, in the same time interval, left-lateralized source activity was found also in premotor and parietal brain regions during natural reading. While brain activity remained left-lateralized in later stages representing word-level processing in posterior and middle ventral temporal regions in the fixed-gaze condition, fixation-related source activity became stronger in the right hemisphere in medial and more anterior ventral temporal brain regions indicating higher-level processing of orthographic information. Although our results show a strong positive relationship between the lateralization of letter-level processing in the two reading modes and suggest lateralized brain activity as a general marker for processing of orthographic information, they also clearly indicate the need for reading research in ecologically valid conditions to identify the neural basis of visuospatial attentional, oculomotor and higher-level processes specific to natural reading.

From theory to practical fundamentals of electroencephalographic source imaging in localizing the epileptogenic zone

With continued advancement in computational technologies, the analysis of electroencephalography (EEG) has shifted from pure visual analysis to a noninvasive computational technique called EEG source imaging (ESI), which involves mathematical modeling of dipolar and distributed sources of a given scalp EEG pattern. ESI is a noninvasive phase I test for presurgical localization of the seizure onset zone in focal epilepsy. It is a relatively inexpensive modality, as it leverages scalp EEG and magnetic resonance imaging (MRI) data already collected typically during presurgical evaluation. With an adequate number of electrodes and combined with patient-specific MRI-based head models, ESI has proven to be a valuable and accurate clinical diagnostic tool for localizing the epileptogenic zone. Despite its advantages, however, ESI is routinely used at only a minority of epilepsy centers. This paper reviews the current evidence and practical fundamentals for using ESI of interictal and ictal epileptic activity during the presurgical evaluation of drug-resistant patients. We identify common errors in processing and interpreting ESI studies, describe the differences in approach needed for localizing interictal and ictal EEG discharges through practical examples, and describe best practices for optimizing the diagnostic information available from these studies.

Mapping and understanding of correlated electroencephalogram (EEG) responses to the newsvendor problem

Decision-making is one of the most critical activities of human beings. To better understand the underlying neurocognitive mechanism while making decisions under an economic context, we designed a decision-making paradigm based on the newsvendor problem (NP) with two scenarios: low-profit margins as the more challenging scenario and high-profit margins as the less difficult one. The EEG signals were acquired from healthy humans while subjects were performing the task. We adopted the Correlated Component Analysis (CorrCA) method to identify linear combinations of EEG channels that maximize the correlation across subjects ([Formula: see text]) or trials ([Formula: see text]). The inter-subject or inter-trial correlation values (ISC or ITC) of the first three components were estimated to investigate the modulation of the task difficulty on subjects' EEG signals and respective correlations. We also calculated the alpha- and beta-band power of the projection components obtained by the CorrCA to assess the brain responses across multiple task periods. Finally, the CorrCA forward models, which represent the scalp projections of the brain activities by the maximally correlated components, were further translated into source distributions of underlying cortical activity using the exact Low Resolution Electromagnetic Tomography Algorithm (eLORETA). Our results revealed strong and significant correlations in EEG signals among multiple subjects and trials during the more difficult decision-making task than the easier one. We also observed that the NP decision-making and feedback tasks desynchronized the normalized alpha and beta powers of the CorrCA components, reflecting the engagement state of subjects. Source localization results furthermore suggested several sources of neural activities during the NP decision-making process, including the dorsolateral prefrontal cortex, anterior PFC, orbitofrontal cortex, posterior cingulate cortex, and somatosensory association cortex.

Automated methodology for optimal selection of minimum electrode subsets for accurate EEG source estimation based on Genetic Algorithm optimization

High-density Electroencephalography (HD-EEG) has proven to be the EEG montage that estimates the neural activity inside the brain with highest accuracy. Multiple studies have reported the effect of electrode number on source localization for specific sources and specific electrode configurations. The electrodes for these configurations are often manually selected to uniformly cover the entire head, going from 32 to 128 electrodes, but electrode configurations are not often selected according to their contribution to estimation accuracy. In this work, an optimization-based study is proposed to determine the minimum number of electrodes that can be used and to identify the optimal combinations of electrodes that can retain the localization accuracy of HD-EEG reconstructions. This optimization approach incorporates scalp landmark positions of widely used EEG montages. In this way, a systematic search for the minimum electrode subset is performed for single- and multiple-source localization problems. The Non-dominated Sorting Genetic Algorithm II (NSGA-II) combined with source reconstruction methods is used to formulate a multi-objective optimization problem that concurrently minimizes (1) the localization error for each source and (2) the number of required EEG electrodes. The method can be used for evaluating the source localization quality of low-density EEG systems (e.g. consumer-grade wearable EEG). We performed an evaluation over synthetic and real EEG datasets with known ground-truth. The experimental results show that optimal subsets with 6 electrodes can attain an equal or better accuracy than HD-EEG (with more than 200 channels) for a single source case. This happened when reconstructing a particular brain activity in more than 88% of the cases in synthetic signals and 63% in real signals, and in more than 88% and 73% of cases when considering optimal combinations with 8 channels. For a multiple-source case of three sources (only with synthetic signals), it was found that optimized combinations of 8, 12 and 16 electrodes attained an equal or better accuracy than HD-EEG with 231 electrodes in at least 58%, 76%, and 82% of cases respectively. Additionally, for such electrode numbers, lower mean errors and standard deviations than with 231 electrodes were obtained.

Evaluation of skull conductivity using SCALE head tissue conductivity estimation using EEG

Inaccurate estimation of skull conductivity is the largest impediment to high-resolution EEG source imaging because of its strong influence and wide variability across individuals. Nonetheless, there is yet no widely applied method for noninvasively measuring individual skull conductivity. We presented a skull conductivity and source location estimation algorithm (SCALE) for simultaneously estimating skull conductivity and the cortical distributions of 18-20 effective sources derived from the EEG data by independent component analysis (ICA). SCALE combines a realistic Finite Element Method (FEM) head model built from a magnetic resonance (MR) head image with the effective source scalp maps to estimate brain-to-skull conductivity ratio (BSCR) and to map the effective sources on the cortical surface. To estimate the robustness of SCALE BSCR estimates, we applied SCALE to MR image and high-density EEG data from ten participants, five having data from 2-3 different tasks and sessions. As expected, across participants SCALE BSCR estimates differed widely (mean 32.8, range 18-78). Within-participant SCALE BSCR estimates were far more consistent than between participants. By incorporating SCALE-optimized distributed EEG source localization, stable functional imaging of cortical EEG effective sources can become routine, giving relatively low-cost EEG imaging a spatial resolution compatible with other brain imaging results and uniquely capable for studying brain dynamics supporting thought and action in laboratory, virtual, and natural environments.

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.

Continuous table tennis is associated with processing in frontal brain areas: an EEG approach

Coordinative challenging exercises in changing environments referred to as open-skill exercises seem to be beneficial on cognitive function. Although electroencephalographic research allows to investigate changes in cortical processing during movement, information about cortical dynamics during open-skill exercise is lacking. Therefore, the present study examines frontal brain activation during table tennis as an open-skill exercise compared to cycling exercise and a cognitive task. 21 healthy young adults conducted three blocks of table tennis, cycling and n-back task. Throughout the experiment, cortical activity was measured using 64-channel EEG system connected to a wireless amplifier. Cortical activity was analyzed calculating theta power (4-7.5 Hz) in frontocentral clusters revealed from independent component analysis. Repeated measures ANOVA was used to identify within subject differences between conditions (table tennis, cycling, n-back; p < .05). ANOVA revealed main-effects of condition on theta power in frontal (p < .01, ηp2 = 0.35) and frontocentral (p < .01, ηp2 = 0.39) brain areas. Post-hoc tests revealed increased theta power in table tennis compared to cycling in frontal brain areas (p < .05, d = 1.42). In frontocentral brain areas, theta power was significant higher in table tennis compared to cycling (p < .01, d = 1.03) and table tennis compared to the cognitive task (p < .01, d = 1.06). Increases in theta power during continuous table tennis may reflect the increased demands in perception and processing of environmental stimuli during open-skill exercise. This study provides important insights that support the beneficial effect of open-skill exercise on brain function and suggest that using open-skill exercise may serve as an intervention to induce activation of the frontal cortex.

Intensity-dependent modulation of cortical somatosensory processing during external, low-frequency peripheral nerve stimulation in humans

External low-frequency peripheral nerve stimulation (LFS) has been proposed as a novel method for neuropathic pain relief. Previous studies have reported that LFS elicits long-term depression-like effects on human pain perception when delivered at noxious intensities, whereas lower intensities are ineffective. To shed light on cortical regions mediating the effects of LFS, we investigated changes in somatosensory-evoked potentials (SEPs) during four LFS intensities. LFS was applied to the radial nerve (600 pulses, 1 Hz) of 24 healthy participants at perception (1 times), low (5 times), medium (10 times), and high intensities (15 times detection threshold). SEPs were recorded during LFS, and averaged SEPs in 10 consecutive 1-min epochs of LFS were analyzed using source dipole modeling. Changes in resting electroencephalography (EEG) were investigated after each LFS block. Source activity in the midcingulate cortex (MCC) decreased linearly during LFS, with greater attenuation at stronger LFS intensities, and in the ipsilateral operculo-insular cortex during the two lowest LFS stimulus intensities. Increased LFS intensities resulted in greater augmentation of contralateral primary sensorimotor cortex (SI/MI) activity. Stronger LFS intensities were followed by increased α (alpha, 9-11 Hz) band power in SI/MI and decreased θ (theta, 3-5 Hz) band power in MCC. Intensity-dependent attenuation of MCC activity with LFS is consistent with a state of long-term depression. Sustained increases in contralateral SI/MI activity suggests that effects of LFS on somatosensory processing may also be dependent on satiation of SI/MI. Further research could clarify if the activation of SI/MI during LFS competes with nociceptive processing in neuropathic pain.NEW & NOTEWORTHY Somatosensory-evoked potentials during low-frequency stimulation of peripheral nerves were examined at graded stimulus intensities. Low-frequency stimulation was associated with decreased responsiveness in the midcingulate cortex and increased responsiveness in primary sensorimotor cortex. Greater intensities were associated with increased midcingulate cortex θ band power and decreased sensorimotor cortex α band power. Results further previous evidence of an inhibition of somatosensory processing during and after low-frequency stimulation and point toward a potential augmentation of activity in somatosensory processing regions.

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

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

Electrocorticographic evidence of a common neurocognitive sequence for mentalizing about the self and others

Neuroimaging studies of mentalizing (i.e., theory of mind) consistently implicate the default mode network (DMN). Nevertheless, the social cognitive functions of individual DMN regions remain unclear, perhaps due to limited spatiotemporal resolution in neuroimaging. Here we use electrocorticography (ECoG) to directly record neuronal population activity while 16 human participants judge the psychological traits of themselves and others. Self- and other-mentalizing recruit near-identical cortical sites in a common spatiotemporal sequence. Activations begin in the visual cortex, followed by temporoparietal DMN regions, then finally in medial prefrontal regions. Moreover, regions with later activations exhibit stronger functional specificity for mentalizing, stronger associations with behavioral responses, and stronger self/other differentiation. Specifically, other-mentalizing evokes slower and longer activations than self-mentalizing across successive DMN regions, implying lengthier processing at higher levels of representation. Our results suggest a common neurocognitive pathway for self- and other-mentalizing that follows a complex spatiotemporal gradient of functional specialization across DMN and beyond.

Assessment of Biomechanical Predictors of Occurrence of Low-Amplitude N1 Potentials Evoked by Naturally Occurring Postural Instabilities

Naturally occurring postural instabilities that occur while standing and walking elicit specific cortical responses in the fronto-central regions (N1 potentials) followed by corrective balance responses to prevent falling. However, no framework could simultaneously track different biomechanical parameters preceding N1s, predict N1s, and assess their predictive power. Here, we propose a framework and show its utility by examining cortical activity (through electroencephalography [EEG]), ground reaction forces, and head acceleration in the anterior-posterior (AP) direction. Ten healthy young adults carried out a balance task of standing on a support surface with or without sway referencing in the AP direction, amplifying, or dampening natural body sway. Using independent components from the fronto-central cortical region obtained from subject-specific head models, we first robustly validated a prior approach on identifying low-amplitude N1 potentials before early signs of balance corrections. Then, a machine learning algorithm was used to evaluate different biomechanical parameters obtained before N1 potentials, to predict the occurrence of N1s. When different biomechanical parameters were directly compared, the time to boundary (TTB) was found to be the best predictor of the occurrence of upcoming low-amplitude N1 potentials during a balance task. Based on these findings, we confirm that the spatio-temporal characteristics of the center of pressure (COP) might serve as an essential parameter that can facilitate the early detection of postural instability in a balance task. Extending our framework to identify such biomarkers in dynamic situations like walking might improve the implementation of corrective balance responses through brain-machine-interfaces to reduce falls in the elderly.

Changes in brain activity of trainees during laparoscopic surgical virtual training assessed with electroencephalography

OBJECTIVE

Evaluate changes in brain activity of trainees during laparoscopic surgical training from electroencephalographic (EEG) signals in an ecological scenario with few restrictions for the user. Design Longitudinal study with two follow-up measurements in the first and last session of a 4-week training with LapSim laparoscopic surgery simulator. Variables analyzed include EEG neuronal activations in theta and alpha bands, tasks performance measures, and subjective measures such as perception of mental workload. Setting Medical School, Universidad de Antioquia, Medellin, Colombia. Participants First-year surgical residents (n = 16, age = 28.0 ± 2.6 years old, right-handed, 9 females) RESULTS: Significant improvements in tasks performance were found together with changes in neuronal activity over frontal and parietal cortex. These changes were also correlated with task performance through training sessions.

CONCLUSIONS

The use of neurophysiological measures such as electroencephalography combined with source separation techniques allows evaluating neural changes associated with motor training. The experiment proposed in this work establishes less controlled recording conditions leading to a more realistic analysis scenario to cognitive assessment in residents training.

Source imaging of high-density visual evoked potentials with multi-scale brain parcellations and connectomes

We describe the multimodal neuroimaging dataset VEPCON (OpenNeuro Dataset ds003505). It includes raw data and derivatives of high-density EEG, structural MRI, diffusion weighted images (DWI) and single-trial behavior (accuracy, reaction time). Visual evoked potentials (VEPs) were recorded while participants (n = 20) discriminated briefly presented faces from scrambled faces, or coherently moving stimuli from incoherent ones. EEG and MRI were recorded separately from the same participants. The dataset contains raw EEG and behavioral data, pre-processed EEG of single trials in each condition, structural MRIs, individual brain parcellations at 5 spatial resolutions (83 to 1015 regions), and the corresponding structural connectomes computed from fiber count, fiber density, average fractional anisotropy and mean diffusivity maps. For source imaging, VEPCON provides EEG inverse solutions based on individual anatomy, with Python and Matlab scripts to derive activity time-series in each brain region, for each parcellation level. The BIDS-compatible dataset can contribute to multimodal methods development, studying structure-function relations, and to unimodal optimization of source imaging and graph analyses, among many other possibilities.