Neuroelectromagnetic forward head modeling toolbox

Patient-specific solution of the electrocorticography forward problem in deforming brain

Invasive intracranial electroencephalography (iEEG), or electrocorticography (ECoG), measures electric potential directly on the surface of the brain and can be used to inform treatment planning for epilepsy surgery. Combined with numerical modeling it can further improve accuracy of epilepsy surgery planning. Accurate solution of the iEEG forward problem, which is a crucial prerequisite for solving the iEEG inverse problemin epilepsy seizure onset zone localization, requires accurate representation of the patient's brain geometry and tissue electrical conductivity after implantation of electrodes. However, implantation of subdural grid electrodes causes the brain to deform, which invalidates preoperatively acquired image data. Moreover, postoperative magnetic resonance imaging (MRI) is incompatible with implanted electrodes and computed tomography (CT) has insufficient range of soft tissue contrast, which precludes both MRI and CT from being used to obtain the deformed postoperative geometry. In this paper, we present a biomechanics-based image warping procedure using preoperative MRI for tissue classification and postoperative CT for locating implanted electrodes to perform non-rigid registration of the preoperative image data to the postoperative configuration. We solve the iEEG forward problem on the predicted postoperative geometry using the finite element method (FEM) which accounts for patient-specific inhomogeneity and anisotropy of tissue conductivity. Results for the simulation of a current source in the brain show large differences in electric potential predicted by the models based on the original images and the deformed images corresponding to the brain geometry deformed by placement of invasive electrodes. Computation of the lead field matrix (useful for solution of the iEEG inverse problem) also showed significant differences between the different models. The results suggest that rapid and accurate solution of the forward problem in a deformed brain for a given patient is achievable.

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.

A Graph Fourier Transform Based Bidirectional Long Short-Term Memory Neural Network for Electrophysiological Source Imaging

Electrophysiological source imaging (ESI) refers to the process of reconstructing underlying activated sources on the cortex given the brain signal measured by Electroencephalography (EEG) or Magnetoencephalography (MEG). Due to the ill-posed nature of ESI, solving ESI requires the design of neurophysiologically plausible regularization or priors to guarantee a unique solution. Recovering focally extended sources is more challenging, and traditionally uses a total variation regularization to promote spatial continuity of the activated sources. In this paper, we propose to use graph Fourier transform (GFT) based bidirectional long-short term memory (BiLSTM) neural network to solve the ESI problem. The GFT delineates the 3D source space into spatially high, medium and low frequency subspaces spanned by corresponding eigenvectors. The low frequency components can naturally serve as a spatially low-band pass filter to reconstruct extended areas of source activation. The BiLSTM is adopted to learn the mapping relationship between the projection of low-frequency graph space and the recorded EEG. Numerical results show the proposed GFT-BiLSTM outperforms other benchmark algorithms in synthetic data under varied signal-to-noise ratios (SNRs). Real data experiments also demonstrate its capability of localizing the epileptogenic zone of epilepsy patients with good accuracy.

An EEG-Based Investigation of the Effect of Perceived Observation on Visual Memory in Virtual Environments

The presence of external observers has been shown to affect performance on cognitive tasks, but the parameters of this impact for different types of tasks and the underlying neural dynamics are less understood. The current study examined the behavioral and brain activity effects of perceived observation on participants’ visual working memory (VWM) in a virtual reality (VR) classroom setting, using the task format as a moderating variable. Participants (n = 21) were equipped with a 57-channel EEG cap, and neural data were collected as they completed two VWM tasks under two observation conditions (observed and not observed) in a within-subjects experimental design. The “observation” condition was operationalized through the addition of a static human avatar in the VR classroom. The avatar’s presence was associated with a significant effect on extending the task response time, but no effect was found on task accuracy. This outcome may have been due to a ceiling effect, as the mean participant task scores were quite high. EEG data analysis supported the behavioral findings by showing consistent differences between the no-observation and observation conditions for one of the VWM tasks only. These neural differences were identified in the dorsolateral prefrontal cortex (dlPFC) and the occipital cortex (OC) regions, with higher theta-band activity occurring in the dlPFC during stimulus encoding and in the OC during response selection when the “observing” avatar was present. These findings provide evidence that perceived observation can inhibit performance during visual tasks by altering attentional focus, even in virtual contexts.

Unsupervised Learning of Brain State Dynamics during Emotion Imagination using High-Density EEG

This study applies adaptive mixture independent component analysis (AMICA) to learn a set of ICA models, each optimized by fitting a distributional model for each identified component process while maximizing component process independence within some subsets of time points of a multi-channel EEG dataset. Here, we applied 20-model AMICA decomposition to long-duration (1-2 hr), high-density (128-channel) EEG data recorded while participants used guided imagination to imagine situations stimulating the experience of 15 specified emotions. These decompositions tended to return models identifying spatiotemporal EEG patterns or states within single emotion imagination periods. Model probability transitions reflected time-courses of EEG dynamics during emotion imagination, which varied across emotions. Transitions between models accounting for imagined "grief" and "happiness" were more abrupt and better aligned with participant reports, while transitions for imagined "contentment" extended into adjoining "relaxation" periods. The spatial distributions of brain-localizable independent component processes (ICs) were more similar within participants (across emotions) than emotions (across participants). Across participants, brain regions with differences in IC spatial distributions (i.e., dipole density) between emotion imagination versus relaxation were identified in or near the left rostrolateral prefrontal, posterior cingulate cortex, right insula, bilateral sensorimotor, premotor, and associative visual cortex. No difference in dipole density was found between positive versus negative emotions. AMICA models of changes in high-density EEG dynamics may allow data-driven insights into brain dynamics during emotional experience, possibly enabling the improved performance of EEG-based emotion decoding and advancing our understanding of emotion.

Large-scale EEG neural network changes in response to therapeutic TMS

BACKGROUND

Transcranial magnetic stimulation (TMS) is an effective therapy for patients with treatment-resistant depression. TMS likely induces functional connectivity changes in aberrant circuits implicated in depression. Electroencephalography (EEG) "microstates" are topographies hypothesized to represent large-scale resting networks. Canonical microstates have recently been proposed as markers for major depressive disorder (MDD), but it is not known if or how they change following TMS.

METHODS

Resting EEG was obtained from 49 MDD patients at baseline and following six weeks of daily TMS. Polarity-insensitive modified k-means clustering was used to segment EEGs into constituent microstates. Microstates were localized via sLORETA. Repeated-measures mixed models tested for within-subject differences over time and t-tests compared microstate features between TMS responder and non-responder groups.

RESULTS

Six microstates (MS-1 - MS-6) were identified from all available EEG data. Clinical response to TMS was associated with increases in features of MS-2, along with decreased metrics of MS-3. Nonresponders showed no significant changes in any microstate. Change in occurrence and coverage of both MS-2 (increased) and MS-3 (decreased) correlated with symptom change magnitude over the course of TMS treatment.

CONCLUSIONS

We identified EEG microstates associated with clinical improvement following a course of TMS therapy. Results suggest selective modulation of resting networks observable by EEG, which is inexpensive and easily acquired in the clinic setting.

Novel Electrode Placement in Electrical Bioimpedance-Based Stroke Detection: Effects on Current Penetration and Injury Characterization in a Finite Element Model

OBJECTIVE

Reducing time-to-treatment and providing acute management in stroke are essential for patient recovery. Electrical bioimpedance (EBI) is an inexpensive and non-invasive tissue measurement approach that has the potential to provide novel continuous intracranial monitoring-something not possible in current standard-of-care. While extensive previous work has evaluated the feasibility of EBI in diagnosing stroke, high-impedance anatomical features in the head have limited clinical translation.

METHODS

The present study introduces novel electrode placements near highly-conductive cerebral spinal fluid (CSF) pathways to enhance electrical current penetration through the skull and increase detection accuracy of neurologic damage. Simulations were conducted on a realistic finite element model (FEM). Novel electrode placements at the tear ducts, soft palate and base of neck were evaluated. Classification accuracy was assessed in the presence of signal noise, patient variability, and electrode positioning.

RESULTS

Algorithms were developed to successfully determine stroke etiology, location, and size relative to impedance measurements from a baseline scan. Novel electrode placements significantly increased stroke classification accuracy at various levels of signal noise (e.g. p < 0.001 at 40 dB). Novel electrodes also amplified current penetration, with up to 30% increase in current density and 57% increased sensitivity in central intracranial regions (p<0.001).

CONCLUSION

These findings support the use of novel electrode placements in EBI to overcome prior limitations, indicating a potential approach to increasing the technology's clinical utility in stroke identification.

SIGNIFICANCE

A non-invasive EBI monitor for stroke could provide essential timely intervention and care of stroke patients.

Mobile Electroencephalography for Studying Neural Control of Human Locomotion

Walking or running in real-world environments requires dynamic multisensory processing within the brain. Studying supraspinal neural pathways during human locomotion provides opportunities to better understand complex neural circuity that may become compromised due to aging, neurological disorder, or disease. Knowledge gained from studies examining human electrical brain dynamics during gait can also lay foundations for developing locomotor neurotechnologies for rehabilitation or human performance. Technical barriers have largely prohibited neuroimaging during gait, but the portability and precise temporal resolution of non-invasive electroencephalography (EEG) have expanded human neuromotor research into increasingly dynamic tasks. In this narrative mini-review, we provide a (1) brief introduction and overview of modern neuroimaging technologies and then identify considerations for (2) mobile EEG hardware, (3) and data processing, (4) including technical challenges and possible solutions. Finally, we summarize (5) knowledge gained from human locomotor control studies that have used mobile EEG, and (6) discuss future directions for real-world neuroimaging research.

Comparing physiological responses during cognitive tests in virtual environments vs. in identical real-world environments

Immersive virtual environments (VEs) are increasingly used to evaluate human responses to design variables. VEs provide a tremendous capacity to isolate and readily adjust specific features of an architectural or product design. They also allow researchers to safely and effectively measure performance factors and physiological responses. However, the success of this form of design-testing depends on the generalizability of response measurements between VEs and real-world contexts. At the current time, there is very limited research evaluating the consistency of human response data across identical real and virtual environments. Rendering tools were used to precisely replicate a real-world classroom in virtual space. Participants were recruited and asked to complete a series of cognitive tests in the real classroom and in the virtual classroom. Physiological data were collected during these tests, including electroencephalography (EEG), electrocardiography (ECG), electrooculography (EOG), galvanic skin response (GSR), and head acceleration. Participants' accuracy on the cognitive tests did not significantly differ between the real classroom and the identical VE. However, the participants answered the tests more rapidly in the VE. No significant differences were found in eye blink rate and heart rate between the real and VR settings. Head acceleration and GSR variance were lower in the VE setting. Overall, EEG frequency band-power was not significantly altered between the real-world classroom and the VE. Analysis of EEG event-related potentials likewise indicated strong similarity between the real-world classroom and the VE, with a single exception related to executive functioning in a color-mismatch task.

Developing a Multi-channel Beamformer by Enhancing Spatially Constrained ICA for Recovery of Correlated EEG Sources

Background:

Brain source imaging based on electroencephalogram (EEG) data aims to recover the neuron populations’ activity producing the scalp potentials. This procedure is known as the EEG inverse problem. Recently, beamformers have gained a lot of consideration in the EEG inverse problem.

Objective:

Beamformers lack acceptable performance in the case of correlated brain sources. These sources happen when some regions of the brain have simultaneous or correlated activities such as auditory stimulation or moving left and right extremities of the body at the same time. In this paper, we have developed a multichannel beamformer robust to correlated sources.

Material and Methods:

In this simulation study, we have looked at the problem of brain source imaging and beamforming from a blind source separation point of view. We focused on the spatially constraint independent component analysis (scICA) algorithm, which generally benefits from the pre-known partial information of mixing matrix, and modified the steps of the algorithm in a way that makes it more robust to correlated sources. We called the modified scICA algorithm Multichannel ICA based EEG Beamformer (MIEB).

Results:

We evaluated the proposed algorithm on simulated EEG data and compared its performance quantitatively with three algorithms scICA, linearly-constrained minimum-variance (LCMV) and Dual-Core beamformers; it is considered that the latter is specially designed to reconstruct correlated sources.

Conclusion:

The MIEB algorithm has much better performance in terms of normalized mean squared error in recovering the correlated/uncorrelated sources both in noise free and noisy synthetic EEG signals. Therefore, it could be used as a robust beamformer in recovering correlated brain sources.

The open EEGLAB portal Interface: High-Performance computing with EEGLAB

EEGLAB signal processing environment is currently the leading open-source software for processing electroencephalographic (EEG) data. The Neuroscience Gateway (NSG, nsgportal.org) is a web and API-based portal allowing users to easily run a variety of neuroscience-related software on high-performance computing (HPC) resources in the U.S. XSEDE network. We have reported recently (Delorme et al., 2019) on the Open EEGLAB Portal expansion of the free NSG services to allow the neuroscience community to build and run MATLAB pipelines using the EEGLAB tool environment. We are now releasing an EEGLAB plug-in, nsgportal, that interfaces EEGLAB with NSG directly from within EEGLAB running on MATLAB on any personal lab computer. The plug-in features a flexible MATLAB graphical user interface (GUI) that allows users to easily submit, interact with, and manage NSG jobs, and to retrieve and examine their results. Command line nsgportal tools supporting these GUI functionalities allow EEGLAB users and plug-in tool developers to build largely automated functions and workflows that include optional NSG job submission and processing. Here we present details on nsgportal implementation and documentation, provide user tutorials on example applications, and show sample test results comparing computation times using HPC versus laptop processing.

Electroencephalographic Study on Sensory Integration in Visually Induced Postural Sway

A periodically reversing optic flow animation, experienced while standing, induces an involuntary sway termed visually induced postural sway (VIPS). Interestingly, VIPS is suppressed during light finger touch to a stationary object. Here, we explored whether VIPS is mediated by parietal field activity in the dorsal visual stream as well as by activity in early visual areas, as has been suggested. We performed a mobile brain/body imaging study using high-density electroencephalographic recording from human participants (11 men and 3 women) standing during exposure to periodically reversing optic flow with and without light finger touch to a stable surface. We also performed recording their visuo-postural tracking movements as a typical visually guided movement to explore differences of cortical process of VIPS from the voluntary visuomotor process involving the dorsal stream. In the visuo-postural tracking condition, the participants moved their center of pressure in time with a slowly oscillating (expanding, shrinking) target rectangle. Source-resolved results showed that alpha band (8-13 Hz) activity in the medial and right occipital cortex during VIPS was modulated by the direction and velocity of optic flow and increased significantly during light finger touch. However, source-resolved potentials from the parietal association cortex showed no such modulation. During voluntary postural sway with feedback (but no visual flow) in which the dorsal stream is involved, sensorimotor areas produced more theta band (4-7 Hz) and less beta band (14-35 Hz) activity than during involuntary VIPS. These results suggest that VIPS involves cortical field dynamic changes in the early visual cortex rather than in the posterior parietal cortex of the visual dorsal stream.

EEG Effective Source Projections Are More Bilaterally Symmetric in Infants Than in Adults

Although anatomical brain hemispheric asymmetries have been clearly documented in the infant brain, findings concerning functional hemispheric specialization have been inconsistent. The present report aims to assess whether bilaterally symmetric synchronous activity between the two hemispheres is a characteristic of the infant brain. To asses cortical bilateral synchronicity, we used decomposition by independent component analysis (ICA) of high-density electroencephalographic (EEG) data collected in an auditory passive oddball paradigm. Decompositions of concatenated 64-channel EEG data epochs from each of 34 typically developing 6-month-old infants and from 18 healthy young adults participating in the same passive auditory oddball protocol were compared to characterize differences in functional brain organization between early life and adulthood. Our results show that infant EEG decompositions comprised a larger number of independent component (IC) effective source processes compatible with a cortical origin and having bilaterally near-symmetric scalp projections (13.8% of the infant data ICs presented a bilateral pattern vs. 4.3% of the adult data ICs). These IC projections could be modeled as the sum of potentials volume-conducted to the scalp from synchronous locally coherent field activities in corresponding left and right cortical source areas. To conclude, in this paradigm, source-resolved infant brain EEG exhibited more bilateral synchronicity than EEG produced by the adult brain, supporting the hypothesis that more strongly unilateral and likely more functionally specialized unihemispheric cortical field activities are concomitants of brain maturation.

More Reliable EEG Electrode Digitizing Methods Can Reduce Source Estimation Uncertainty, but Current Methods Already Accurately Identify Brodmann Areas

Electroencephalography (EEG) and source estimation can be used to identify brain areas activated during a task, which could offer greater insight on cortical dynamics. Source estimation requires knowledge of the locations of the EEG electrodes. This could be provided with a template or obtained by digitizing the EEG electrode locations. Operator skill and inherent uncertainties of a digitizing system likely produce a range of digitization reliabilities, which could affect source estimation and the interpretation of the estimated source locations. Here, we compared the reliabilities of five digitizing methods (ultrasound, structured-light 3D scan, infrared 3D scan, motion capture probe, and motion capture) and determined the relationship between digitization reliability and source estimation uncertainty, assuming other contributors to source estimation uncertainty were constant. We digitized a mannequin head using each method five times and quantified the reliability and validity of each method. We created five hundred sets of electrode locations based on our reliability results and applied a dipole fitting algorithm (DIPFIT) to perform source estimation. The motion capture method, which recorded the locations of markers placed directly on the electrodes had the best reliability with an average electrode variability of 0.001 cm. Then, in order of decreasing reliability were the method using a digitizing probe in the motion capture system, an infrared 3D scanner, a structured-light 3D scanner, and an ultrasound digitization system. Unsurprisingly, uncertainty of the estimated source locations increased with greater variability of EEG electrode locations and less reliable digitizing systems. If EEG electrode location variability was ∽1 cm, a single source could shift by as much as 2 cm. To help translate these distances into practical terms, we quantified Brodmann area accuracy for each digitizing method and found that the average Brodmann area accuracy for all digitizing methods was >80%. Using a template of electrode locations reduced the Brodmann area accuracy to ∽50%. Overall, more reliable digitizing methods can reduce source estimation uncertainty, but the significance of the source estimation uncertainty depends on the desired spatial resolution. For accurate Brodmann area identification, any of the digitizing methods tested can be used confidently.

ROAST: An Open-Source, Fully-Automated, Realistic Volumetric-Approach-Based Simulator For TES

Research in the area of transcranial electrical stimulation (TES) often relies on computational models of current flow in the brain. Models are built on magnetic resonance images (MRI) of the human head to capture detailed individual anatomy. To simulate current flow, MRIs have to be segmented, virtual electrodes have to be placed on the scalp, the volume is tessellated into a mesh, and the finite element model is solved numerically to estimate the current flow. Various software tools are available for each step, as well as processing pipelines that connect these tools for automated or semi-automated processing. The goal of the present tool - ROAST - is to provide an end-to-end pipeline that can automatically process individual heads with realistic volumetric anatomy leveraging open-source software (SPM8, iso2mesh and getDP) and custom scripts to improve segmentation and execute electrode placement. When we compare the results on a standard head with other major commercial software for finite element modeling (ScanIP, Abaqus), ROAST only leads to a small difference of 9% in the estimated electric field in the brain. We obtain a larger difference of 47% when comparing results with SimNIBS, an automated pipeline that is based on surface segmentation of the head. We release ROAST as an open-source, fully-automated pipeline at https://www.parralab.org/roast/.

Ear-EEG Forward Models: Improved Head-Models for Ear-EEG

Computational models for mapping electrical sources in the brain to potentials on the scalp have been widely explored. However, current models do not describe the external ear anatomy well, and is therefore not suitable for ear-EEG recordings. Here we present an extension to existing computational models, by incorporating an improved description of the external ear anatomy based on 3D scanned impressions of the ears. The result is a method to compute an ear-EEG forward model, which enables mapping of sources in the brain to potentials in the ear. To validate the method, individualized ear-EEG forward models were computed for four subjects, and ear-EEG and scalp EEG were recorded concurrently from the subjects in a study comprising both auditory and visual stimuli. The EEG recordings were analyzed with independent component analysis (ICA) and using the individualized ear-EEG forward models, single dipole fitting was performed for each independent component (IC). A subset of ICs were selected, based on how well they were modeled by a single dipole in the brain volume. The correlation between the topographic IC map and the topographic map predicted by the forward model, was computed for each IC. Generally, the correlation was high in the ear closest to the dipole location, showing that the ear-EEG forward models provided a good model to predict ear potentials. In addition, we demonstrated that the developed forward models can be used to explore the sensitivity to brain sources for different ear-EEG electrode configurations. We consider the proposed method to be an important step forward in the characterization and utilization of ear-EEG.

Reduced premovement positivity during the stimulus-response interval precedes errors: Using single-trial and regression ERPs to understand performance deficits in ADHD

Brain mechanisms linked to incorrect response selections made under time pressure during cognitive task performance are poorly understood, particularly in adolescents with attention-deficit hyperactivity disorder (ADHD). Using subject-specific multimodal imaging (electroencephalogram, magnetic resonance imaging, behavior) during flanker task performance by a sample of 94 human adolescents (mean age = 15.5 years, 50% female) with varying degrees of ADHD symptomatology, we examined the degree to which amplitude features of source-resolved event-related potentials (ERPs) from brain-independent component processes within a critical (but often ignored) period in the action selection process, the stimulus-response interval, were associated with motor response errors (across trials) and error rates (across individuals). Response errors were typically preceded by two smaller peaks in both trial-level and trial-averaged ERP projections from posterior medial frontal cortex (pMFC): a frontocentral P3 peaking about 390 ms after stimulus onset, and a premovement positivity (PMP) peaking about 110 ms before the motor response. Separating overlapping stimulus-locked and response-locked ERP contributions using a "regression ERP" approach showed that trial errors and participant error rates were primarily associated with smaller PMP, and not with frontocentral P3. Moreover, smaller PMP mediated the association between larger numbers of errors and ADHD symptoms, suggesting the possible value of using PMP as an intervention target to remediate performance deficits in ADHD.

Realistic volumetric-approach to simulate transcranial electric stimulation-ROAST-a fully automated open-source pipeline

OBJECTIVE

Research in the area of transcranial electrical stimulation (TES) often relies on computational models of current flow in the brain. Models are built based on magnetic resonance images (MRI) of the human head to capture detailed individual anatomy. To simulate current flow on an individual, the subject's MRI is segmented, virtual electrodes are placed on this anatomical model, the volume is tessellated into a mesh, and a finite element model (FEM) is solved numerically to estimate the current flow. Various software tools are available for each of these steps, as well as processing pipelines that connect these tools for automated or semi-automated processing. The goal of the present tool-realistic volumetric-approach to simulate transcranial electric simulation (ROAST)-is to provide an end-to-end pipeline that can automatically process individual heads with realistic volumetric anatomy leveraging open-source software and custom scripts to improve segmentation and execute electrode placement.

APPROACH

ROAST combines the segmentation algorithm of SPM12, a Matlab script for touch-up and automatic electrode placement, the finite element mesher iso2mesh and the solver getDP. We compared its performance with commercial FEM software, and SimNIBS, a well-established open-source modeling pipeline.

MAIN RESULTS

The electric fields estimated with ROAST differ little from the results obtained with commercial meshing and FEM solving software. We also do not find large differences between the various automated segmentation methods used by ROAST and SimNIBS. We do find bigger differences when volumetric segmentation are converted into surfaces in SimNIBS. However, evaluation on intracranial recordings from human subjects suggests that ROAST and SimNIBS are not significantly different in predicting field distribution, provided that users have detailed knowledge of SimNIBS.

SIGNIFICANCE

We hope that the detailed comparisons presented here of various choices in this modeling pipeline can provide guidance for future tool development. We released ROAST as an open-source, easy-to-install and fully-automated pipeline for individualized TES modeling.

Cortical Functional Domains Show Distinctive Oscillatory Dynamic in Bimanual and Mirror Visual Feedback Tasks

It is believed that Mirror Visual Feedback (MVF) increases the interlimb transfer but the exact mechanism is still a matter of debate. The aim of this study was to compare between a bimanual task (BM) and a MVF task, within functionally rather than geometrically defined cortical domains. Measure Projection Analysis (MPA) approach was applied to compare the dynamic oscillatory activity (event-related synchronization/desynchronization ERS/ERD) between and within domains. EEG was recorded in 14 healthy participants performing a BM and an MVF task with the right hand. The MPA was applied on fitted equivalent current dipoles based on independent components to define domains containing functionally similar areas. The measure of intradomain similarity was a "signed mutual information," a parameter based on the coherence. Domain analysis was performed for joint tasks (BM and MVF) and for each task separately. MVF created 9 functional domains while MB task had only 4 functionally distinctive domains, two over the left hemispheres and two bilateraly. For all domains identified for BM task alone, similar domains could be identified in MVF and joint tasks analysis. In addition MVF had domains related to motor planning on the right hemisphere and to self-recognition of action. For joint tasks analysis, seven domains were identified, with similar functions for the left and the right hand with exception of a domain covering BA32 (self-recognition of action) of the left hand only. In joint task domain analysis, the ERD/ERS showed a larger difference between domains than between tasks. All domains which involved the sensory cortex had a visible beta ERS at the onset of movement, and post movement beta ERS. The frequency of ERD varied between domains. Largest difference between tasks existed in domains responsible for the awareness of action. In conclusion, functionally distinctive domains have different ERD/ERS patterns, similar for both tasks. MVF activates contralateral hemisphere in similar manner to BM movements, while at the same time also activating the ipsilateral hemisphere. Significance: Following stroke cortical activation and interhemispheric inhibition from the contralesional side is reduced. MVF creates stronger ipsilateral activity than BM, which is highly relevant of neurorehabilitation of movements.

Lingering Sound: Event-Related Phase-Amplitude Coupling and Phase-Locking in Fronto-Temporo-Parietal Functional Networks During Memory Retrieval of Music Melodies

Brain oscillations and connectivity have emerged as promising measures of evaluating memory processes, including encoding, maintenance, and retrieval, as well as the related executive function. Although many studies have addressed the neural mechanisms underlying working memory, most of these studies have focused on the visual modality. Neurodynamics and functional connectivity related to auditory working memory are yet to be established. In this study, we explored the dynamic of high density (128-channel) electroencephalography (EEG) in a musical delayed match-to-sample task (DMST), in which 36 participants were recruited and were instructed to recognize and distinguish the target melodies from similar distractors. Event-related spectral perturbations (ERSPs), event-related phase-amplitude couplings (ERPACs), and phase-locking values (PLVs) were used to determine the corresponding brain oscillations and connectivity. First, we observed that low-frequency oscillations in the frontal, temporal, and parietal regions were increased during the processing of both target and distracting melodies. Second, the cross-frequency coupling between low-frequency phases and high-frequency amplitudes was elevated in the frontal and parietal regions when the participants were distinguishing between the target from distractor, suggesting that the phase-amplitude coupling could be an indicator of neural mechanisms underlying memory retrieval. Finally, phase-locking, an index evaluating brain functional connectivity, revealed that there was fronto-temporal phase-locking in the theta band and fronto-parietal phase-locking in the alpha band during the recognition of the two stimuli. These findings suggest the existence of functional connectivity and the phase-amplitude coupling in the neocortex during musical memory retrieval, and provide a highly resolved timeline to evaluate brain dynamics. Furthermore, the inter-regional phase-locking and phase-amplitude coupling among the frontal, temporal and parietal regions occurred at the very beginning of musical memory retrieval, which might reflect the precise timing when cognitive resources were involved in the retrieval of targets and the rejection of similar distractors. To the best of our knowledge, this is the first EEG study employing a naturalistic task to study auditory memory processes and functional connectivity during memory retrieval, results of which can shed light on the use of natural stimuli in studies that are closer to the real-life applications of cognitive evaluations, mental treatments, and brain-computer interface.

The effects of asymmetric volume conductor modeling on non-invasive fetal ECG extraction

OBJECTIVE

Non-invasive fetal electrocardiography (NI-FECG) shows promise for capturing novel physiological information that may indicate signs of fetal distress. However, significant deterioration in NI-FECG signal quality occurs during the presence of a highly non-conductive layer known as vernix caseosa which forms on the fetal body surface beginning in approximately the 28th week of gestation. This work investigates asymmetric modeling of vernix caseosa and other maternal-fetal tissues in accordance with clinical observations and assesses their impacts for NI-FECG signal processing.

APPROACH

We develop a process for simulating dynamic maternal-fetal abdominal ECG mixtures using a synthetic cardiac source model embedded in a finite element volume conductor. Using this process, changes in NI-FECG signal morphology are assessed in an extensive set of finite element models including spatially variable distributions of vernix caseosa.

MAIN RESULTS

Our simulations show that volume conductor asymmetry can result in over 70% error in the observed T/QRS ratio and significant changes to signal morphology compared to a homogeneous volume conductor model. Volume conductor effects must be considered when analyzing T/QRS ratios obtained via NI-FECG and should be considered in future algorithm benchmarks using simulated data.

SIGNIFICANCE

This work shows that without knowledge of the influence of volume conductor effects, clinical evaluation of the T/QRS ratio derived via NI-FECG should be avoided.

Modeling brain dynamic state changes with adaptive mixture independent component analysis

There is a growing interest in neuroscience in assessing the continuous, endogenous, and nonstationary dynamics of brain network activity supporting the fluidity of human cognition and behavior. This non-stationarity may involve ever-changing formation and dissolution of active cortical sources and brain networks. However, unsupervised approaches to identify and model these changes in brain dynamics as continuous transitions between quasi-stable brain states using unlabeled, noninvasive recordings of brain activity have been limited. This study explores the use of adaptive mixture independent component analysis (AMICA) to model multichannel electroencephalographic (EEG) data with a set of ICA models, each of which decomposes an adaptively learned portion of the data into statistically independent sources. We first show that AMICA can segment simulated quasi-stationary EEG data and accurately identify ground-truth sources and source model transitions. Next, we demonstrate that AMICA decomposition, applied to 6-13 channel scalp recordings from the CAP Sleep Database, can characterize sleep stage dynamics, allowing 75% accuracy in identifying transitions between six sleep stages without use of EEG power spectra. Finally, applied to 30-channel data from subjects in a driving simulator, AMICA identifies models that account for EEG during faster and slower response to driving challenges, respectively. We show changes in relative probabilities of these models allow effective prediction of subject response speed and moment-by-moment characterization of state changes within single trials. AMICA thus provides a generic unsupervised approach to identifying and modeling changes in EEG dynamics. Applied to continuous, unlabeled multichannel data, AMICA may likely be used to detect and study any changes in cognitive states.

The Discontinuous Galerkin Finite Element Method for Solving the MEG and the Combined MEG/EEG Forward Problem

In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches in EEG or combined MEG/EEG source analysis scenarios, which motivates a further evaluation of DG-FEM for applications in bioelectromagnetism.

EEG source reconstruction evidence for the noun-verb neural dissociation along semantic dimensions

One of the long-standing issues in neurolinguistic research is about the neural basis of word representation, concerning whether grammatical classification or semantic difference causes the neural dissociation of brain activity patterns when processing different word categories, especially nouns and verbs. To disentangle this puzzle, four orthogonalized word categories in Chinese: unambiguous nouns (UN), unambiguous verbs (UV), ambiguous words with noun-biased semantics (AN), and ambiguous words with verb-biased semantics (AV) were adopted in an auditory task for recording electroencephalographic (EEG) signals from 128 electrodes on the scalps of twenty-two subjects. With the advanced current density reconstruction (CDR) algorithm and the constraint of standardized low-resolution electromagnetic tomography, the spatiotemporal brain dynamics of word processing were explored with the results that in multiple time periods including P1 (60-90ms), N1 (100-140ms), P200 (150-250ms) and N400 (350-450ms), noun-verb dissociation over the parietal-occipital and frontal-central cortices appeared not only between the UN-UV grammatical classes but also between the grammatically identical but semantically different AN-AV pairs. The apparent semantic dissociation within one grammatical class strongly suggests that the semantic difference rather than grammatical classification could be interpreted as the origin of the noun-verb neural dissociation. Our results also revealed that semantic dissociation occurs from an early stage and repeats in multiple phases, thus supporting a functionally hierarchical word processing mechanism.

A low-dimensional representation for individual head geometries

In the estimation of individual head geometries for source localization and electrical stimulation in neuroelectric investigations and application, mostly complex geometrical models are directly extracted from anatomical images. We present a novel method that uses a dimensionality reduction from thousands down to the range of tens of parameters to successfully represent an individual 4-shell Boundary Element Method (BEM) head model, which can successively be fitted to any kind of data from an individual head (e.g. headshape, impedances) and then used for individual head model creation. The results show, that around 15 - 20 components can lead to satisfactory results.

High-resolution EEG source imaging of one-year-old children

Recently we described an iterative skull conductivity and source location estimation (SCALE) algorithm for simultaneously estimating head tissue conductivities and brain source locations. SCALE uses a realistic FEM forward problem head model and scalp maps of 10 or more near-dipolar sources identified by independent component analysis (ICA) decomposition of sufficient high-density EEG data. In this study, we applied SCALE to 20 minutes of 64-channel EEG data and magnetic resonance (MR) head images from four twelve-months-of-age infants. For each child, we selected 15-16 near-dipolar independent components from multiple-model adaptive mixture ICA (AMICA) decomposition of their EEG data. SCALE converged to brain-to-skull conductivity ratio (BSCR) estimates in the 10-12 range and mostly compact gyral or sulcal cortical distributions for the IC sources.

Measurements and models of electric fields in the in vivo human brain during transcranial electric stimulation

Transcranial electric stimulation aims to stimulate the brain by applying weak electrical currents at the scalp. However, the magnitude and spatial distribution of electric fields in the human brain are unknown. We measured electric potentials intracranially in ten epilepsy patients and estimated electric fields across the entire brain by leveraging calibrated current-flow models. When stimulating at 2 mA, cortical electric fields reach 0.8 V/m, the lower limit of effectiveness in animal studies. When individual whole-head anatomy is considered, the predicted electric field magnitudes correlate with the recorded values in cortical (r = 0.86) and depth (r = 0.88) electrodes. Accurate models require adjustment of tissue conductivity values reported in the literature, but accuracy is not improved when incorporating white matter anisotropy or different skull compartments. This is the first study to validate and calibrate current-flow models with in vivo intracranial recordings in humans, providing a solid foundation to target stimulation and interpret clinical trials.

DOI:http://dx.doi.org/10.7554/eLife.18834.001

Forward and inverse effects of the complete electrode model in neonatal EEG

The effect of the complete electrode model on electroencephalography forward and inverse computations is explored. A realistic neonatal head model, including a skull structure with fontanels and sutures, is used. The electrode and skull modeling differences are analyzed and compared with each other. The results suggest that the complete electrode model can be considered as an integral part of the outer head model. To achieve optimal source localization results, accurate electrode modeling might be necessary.

This paper investigates finite element method-based modeling in the context of neonatal electroencephalography (EEG). In particular, the focus lies on electrode boundary conditions. We compare the complete electrode model (CEM) with the point electrode model (PEM), which is the current standard in EEG. In the CEM, the voltage experienced by an electrode is modeled more realistically as the integral average of the potential distribution over its contact surface, whereas the PEM relies on a point value. Consequently, the CEM takes into account the subelectrode shunting currents, which are absent in the PEM. In this study, we aim to find out how the electrode voltage predicted by these two models differ, if standard size electrodes are attached to a head of a neonate. Additionally, we study voltages and voltage variation on electrode surfaces with two source locations: 1) next to the C6 electrode and 2) directly under the Fz electrode and the frontal fontanel. A realistic model of a neonatal head, including a skull with fontanels and sutures, is used. Based on the results, the forward simulation differences between CEM and PEM are in general small, but significant outliers can occur in the vicinity of the electrodes. The CEM can be considered as an integral part of the outer head model. The outcome of this study helps understanding volume conduction of neonatal EEG, since it enlightens the role of advanced skull and electrode modeling in forward and inverse computations.

NEW & NOTEWORTHY The effect of the complete electrode model on electroencephalography forward and inverse computations is explored. A realistic neonatal head model, including a skull structure with fontanels and sutures, is used. The electrode and skull modeling differences are analyzed and compared with each other. The results suggest that the complete electrode model can be considered as an integral part of the outer head model. To achieve optimal source localization results, accurate electrode modeling might be necessary.

Electroencephalography (EEG) forward modeling via H(div) finite element sources with focal interpolation

The goal of this study is to develop focal, accurate and robust finite element method (FEM) based approaches which can predict the electric potential on the surface of the computational domain given its structure and internal primary source current distribution. While conducting an EEG evaluation, the placement of source currents to the geometrically complex grey matter compartment is a challenging but necessary task to avoid forward errors attributable to tissue conductivity jumps. Here, this task is approached via a mathematically rigorous formulation, in which the current field is modeled via divergence conforming H(div) basis functions. Both linear and quadratic functions are used while the potential field is discretized via the standard linear Lagrangian (nodal) basis. The resulting model includes dipolar sources which are interpolated into a random set of positions and orientations utilizing two alternative approaches: the position based optimization (PBO) and the mean position/orientation (MPO) method. These results demonstrate that the present dipolar approach can reach or even surpass, at least in some respects, the accuracy of two classical reference methods, the partial integration (PI) and St. Venant (SV) approach which utilize monopolar loads instead of dipolar currents.

The New York Head-A precise standardized volume conductor model for EEG source localization and tES targeting

In source localization of electroencephalograpic (EEG) signals, as well as in targeted transcranial electric current stimulation (tES), a volume conductor model is required to describe the flow of electric currents in the head. Boundary element models (BEM) can be readily computed to represent major tissue compartments, but cannot encode detailed anatomical information within compartments. Finite element models (FEM) can capture more tissue types and intricate anatomical structures, but with the higher precision also comes the need for semi-automated segmentation, and a higher computational cost. In either case, adjusting to the individual human anatomy requires costly magnetic resonance imaging (MRI), and thus head modeling is often based on the anatomy of an 'arbitrary' individual (e.g. Colin27). Additionally, existing reference models for the human head often do not include the cerebro-spinal fluid (CSF), and their field of view excludes portions of the head and neck-two factors that demonstrably affect current-flow patterns. Here we present a highly detailed FEM, which we call ICBM-NY, or "New York Head". It is based on the ICBM152 anatomical template (a non-linear average of the MRI of 152 adult human brains) defined in MNI coordinates, for which we extended the field of view to the neck and performed a detailed segmentation of six tissue types (scalp, skull, CSF, gray matter, white matter, air cavities) at 0.5mm(3) resolution. The model was solved for 231 electrode locations. To evaluate its performance, additional FEMs and BEMs were constructed for four individual subjects. Each of the four individual FEMs (regarded as the 'ground truth') is compared to its BEM counterpart, the ICBM-NY, a BEM of the ICBM anatomy, an 'individualized' BEM of the ICBM anatomy warped to the individual head surface, and FEMs of the other individuals. Performance is measured in terms of EEG source localization and tES targeting errors. Results show that the ICBM-NY outperforms FEMs of mismatched individual anatomies as well as the BEM of the ICBM anatomy according to both criteria. We therefore propose the New York Head as a new standard head model to be used in future EEG and tES studies whenever an individual MRI is not available. We release all model data online at neuralengr.com/nyhead/ to facilitate broad adoption.

Conscious and Non-conscious Representations of Emotional Faces in Asperger's Syndrome

Several neuroimaging studies have suggested that the low spatial frequency content in an emotional face mainly activates the amygdala, pulvinar, and superior colliculus especially with fearful faces(1-3). These regions constitute the limbic structure in non-conscious perception of emotions and modulate cortical activity either directly or indirectly(2). In contrast, the conscious representation of emotions is more pronounced in the anterior cingulate, prefrontal cortex, and somatosensory cortex for directing voluntary attention to details in faces(3,4). Asperger's syndrome (AS)(5,6) represents an atypical mental disturbance that affects sensory, affective and communicative abilities, without interfering with normal linguistic skills and intellectual ability. Several studies have found that functional deficits in the neural circuitry important for facial emotion recognition can partly explain social communication failure in patients with AS(7-9). In order to clarify the interplay between conscious and non-conscious representations of emotional faces in AS, an EEG experimental protocol is designed with two tasks involving emotionality evaluation of either photograph or line-drawing faces. A pilot study is introduced for selecting face stimuli that minimize the differences in reaction times and scores assigned to facial emotions between the pretested patients with AS and IQ/gender-matched healthy controls. Information from the pretested patients was used to develop the scoring system used for the emotionality evaluation. Research into facial emotions and visual stimuli with different spatial frequency contents has reached discrepant findings depending on the demographic characteristics of participants and task demands(2). The experimental protocol is intended to clarify deficits in patients with AS in processing emotional faces when compared with healthy controls by controlling for factors unrelated to recognition of facial emotions, such as task difficulty, IQ and gender.

Incorporating and Compensating Cerebrospinal Fluid in Surface-Based Forward Models of Magneto- and Electroencephalography

MEG/EEG source imaging is usually done using a three-shell (3-S) or a simpler head model. Such models omit cerebrospinal fluid (CSF) that strongly affects the volume currents. We present a four-compartment (4-C) boundary-element (BEM) model that incorporates the CSF and is computationally efficient and straightforward to build using freely available software. We propose a way for compensating the omission of CSF by decreasing the skull conductivity of the 3-S model, and study the robustness of the 4-C and 3-S models to errors in skull conductivity. We generated dense boundary meshes using MRI datasets and automated SimNIBS pipeline. Then, we built a dense 4-C reference model using Galerkin BEM, and 4-C and 3-S test models using coarser meshes and both Galerkin and collocation BEMs. We compared field topographies of cortical sources, applying various skull conductivities and fitting conductivities that minimized the relative error in 4-C and 3-S models. When the CSF was left out from the EEG model, our compensated, unbiased approach improved the accuracy of the 3-S model considerably compared to the conventional approach, where CSF is neglected without any compensation (mean relative error < 20% vs. > 40%). The error due to the omission of CSF was of the same order in MEG and compensated EEG. EEG has, however, large overall error due to uncertain skull conductivity. Our results show that a realistic 4-C MEG/EEG model can be implemented using standard tools and basic BEM, without excessive workload or computational burden. If the CSF is omitted, compensated skull conductivity should be used in EEG.

The Unfitted Discontinuous Galerkin Method for Solving the EEG Forward Problem

OBJECTIVE

The purpose of this study is to introduce and evaluate the unfitted discontinuous Galerkin finite element method (UDG-FEM) for solving the electroencephalography (EEG) forward problem.

METHODS

This new approach for source analysis does not use a geometry conforming volume triangulation, but instead uses a structured mesh that does not resolve the geometry. The geometry is described using level set functions and is incorporated implicitly in its mathematical formulation. As no triangulation is necessary, the complexity of a simulation pipeline and the need for manual interaction for patient-specific simulations can be reduced and is comparable with that of the FEM for hexahedral meshes. In addition, it maintains conservation laws on a discrete level. Here, we present the theory for UDG-FEM forward modeling, its verification using quasi-analytical solutions in multilayer sphere models and an evaluation in a comparison with a discontinuous Galerkin (DG-FEM) method on hexahedral and on conforming tetrahedral meshes. We furthermore apply the UDG-FEM forward approach in a realistic head model simulation study.

RESULTS

The results show convergence to the quasi-analytical solution and indicate a good accuracy of UDG-FEM. UDG-FEM performs comparable or even better than DG-FEM on a conforming tetrahedral mesh while providing a less complex simulation pipeline. When compared to DG-FEM on hexahedral meshes, an overall better accuracy is achieved.

CONCLUSION

The UDG-FEM approach is an accurate, flexible, and promising method to solve the EEG forward problem.

SIGNIFICANCE

This study shows the first application of the UDG-FEM approach to the EEG forward problem.

Quantifying the Test-Retest Reliability of Magnetoencephalography Resting-State Functional Connectivity

The coordinated activity of the resting-state brain can be evaluated with magnetoencephalography (MEG) for distinct brain rhythms by performing source reconstruction to estimate the activities of target brain regions and employing one of the many existent functional connectivity (FC) algorithms. Although this procedure has been applied in a great amount of studies both with healthy and pathological populations, the reliability of such FC estimates is unknown, and this impairs the use of resting-state MEG FC at the individual level. In this study, the test-retest reliability of MEG resting FC was evaluated by exploring both within- and between-subject variability in FC in 16 healthy subjects who underwent three resting-state MEG scans. FC was computed after beamforming source reconstruction with four popular FC metrics: phase-locking value (PLV), phase lag index (PLI), direct envelope correlation (d-ecor), and envelope correlation with leakage correction (lc-ecor). Then, test-restest reliability and within- and between-subject agreement were evaluated with the intraclass correlation coefficient (ICC) and Kendall's W, respectively. Reliability was found to depend on the FC metric, the frequency band, and the specific link. As a general trend, greater test-retest reliability was found for PLV in theta to gamma, and for lc-ecor and d-ecor in beta. Further inspection of the ICC distribution revealed that volume conduction effects could be contributing to high ICC in PLV and d-ecor. In addition, stronger links were found to be more reliable. Overall, this encourages the further use of resting-state MEG FC for individual-level studies, especially with PLV or envelope correlation metrics.

Unsupervised Segmentation of Head Tissues from Multi-modal MR Images for EEG Source Localization

In this paper, we present and evaluate an automatic unsupervised segmentation method, hierarchical segmentation approach (HSA)-Bayesian-based adaptive mean shift (BAMS), for use in the construction of a patient-specific head conductivity model for electroencephalography (EEG) source localization. It is based on a HSA and BAMS for segmenting the tissues from multi-modal magnetic resonance (MR) head images. The evaluation of the proposed method was done both directly in terms of segmentation accuracy and indirectly in terms of source localization accuracy. The direct evaluation was performed relative to a commonly used reference method brain extraction tool (BET)-FMRIB's automated segmentation tool (FAST) and four variants of the HSA using both synthetic data and real data from ten subjects. The synthetic data includes multiple realizations of four different noise levels and several realizations of typical noise with a 20% bias field level. The Dice index and Hausdorff distance were used to measure the segmentation accuracy. The indirect evaluation was performed relative to the reference method BET-FAST using synthetic two-dimensional (2D) multimodal magnetic resonance (MR) data with 3% noise and synthetic EEG (generated for a prescribed source). The source localization accuracy was determined in terms of localization error and relative error of potential. The experimental results demonstrate the efficacy of HSA-BAMS, its robustness to noise and the bias field, and that it provides better segmentation accuracy than the reference method and variants of the HSA. They also show that it leads to a more accurate localization accuracy than the commonly used reference method and suggest that it has potential as a surrogate for expert manual segmentation for the EEG source localization problem.

Intraneural stimulation elicits discrimination of textural features by artificial fingertip in intact and amputee humans

Restoration of touch after hand amputation is a desirable feature of ideal prostheses. Here, we show that texture discrimination can be artificially provided in human subjects by implementing a neuromorphic real-time mechano-neuro-transduction (MNT), which emulates to some extent the firing dynamics of SA1 cutaneous afferents. The MNT process was used to modulate the temporal pattern of electrical spikes delivered to the human median nerve via percutaneous microstimulation in four intact subjects and via implanted intrafascicular stimulation in one transradial amputee. Both approaches allowed the subjects to reliably discriminate spatial coarseness of surfaces as confirmed also by a hybrid neural model of the median nerve. Moreover, MNT-evoked EEG activity showed physiologically plausible responses that were superimposable in time and topography to the ones elicited by a natural mechanical tactile stimulation. These findings can open up novel opportunities for sensory restoration in the next generation of neuro-prosthetic hands.

DOI:http://dx.doi.org/10.7554/eLife.09148.001

Our hands provide us with a wide variety of information about our surroundings, enabling us to detect pain, temperature and pressure. Our sense of touch also allows us to interact with objects by feeling their texture and solidity. However, completely reproducing a sense of touch in artificial or prosthetic hands has proven challenging. While commercial prostheses can mimic the range of movements of natural limbs, even the latest experimental prostheses have only a limited ability to ‘feel’ the objects being manipulated. Oddo, Raspopovic et al. have now brought this ability a step closer by exploiting an artificial fingertip and appropriate neural interfaces through which different textures can be identified.

The initial experiments were performed in four healthy volunteers with intact limbs. Oddo, Raspopovic et al. connected the artificial fingertip to the volunteers via an electrode inserted into a nerve in the arm. When moved over a rough surface, sensors in the fingertip produced patterns of electrical pulses that stimulated the nerve, causing the volunteers to feel like they were touching the surface. The volunteers were even able to tell the difference between the different surface textures the artificial fingertip moved across.

The temporary electrodes used in this group of volunteers are unsuitable for use with prosthetic limbs because they can easily be knocked out of position. Therefore, in a further experiment involving a volunteer who had undergone an arm amputation a number of years previously, Oddo, Raspopovic et al. tested an implanted electrode array that could, in principle, remain in place long-term. This volunteer could also identify the different textures the artificial fingertip touched, with a slightly higher degree of accuracy than the previous group of intact volunteers. Further studies are now required to explore the potential of this approach in larger groups of volunteers.

DOI:http://dx.doi.org/10.7554/eLife.09148.002

ICA-derived cortical responses indexing rapid multi-feature auditory processing in six-month-old infants

The abilities of infants to perceive basic acoustic differences, essential for language development, can be studied using auditory event-related potentials (ERPs). However, scalp-channel averaged ERPs sum volume-conducted contributions from many cortical areas, reducing the functional specificity and interpretability of channel-based ERP measures. This study represents the first attempt to investigate rapid auditory processing in infancy using independent component analysis (ICA), allowing exploration of source-resolved ERP dynamics and identification of ERP cortical generators. Here, we recorded 60-channel EEG data in 34 typically developing 6-month-old infants during a passive acoustic oddball paradigm presenting 'standard' tones interspersed with frequency- or duration-deviant tones. ICA decomposition was applied to single-subject EEG data. The best-fitting equivalent dipole or bilaterally symmetric dipole pair was then estimated for each resulting independent component (IC) process using a four-layer infant head model. Similar brain-source ICs were clustered across subjects. Results showed ERP contributions from auditory cortex and multiple extra-auditory cortical areas (often, bilaterally paired). Different cortical source combinations contributed to the frequency- and duration-deviant ERP peak sequences. For ICs in an ERP-dominant source cluster located in or near the mid-cingulate cortex, source-resolved frequency-deviant response N2 latency and P3 amplitude at 6 months-of-age predicted vocabulary size at 20 months-of-age. The same measures for scalp channel F6 (though not for other frontal channels) showed similar but weaker correlations. These results demonstrate the significant potential of ICA analyses to facilitate a deeper understanding of the neural substrates of infant sensory processing.

Simultaneous head tissue conductivity and EEG source location estimation

Accurate electroencephalographic (EEG) source localization requires an electrical head model incorporating accurate geometries and conductivity values for the major head tissues. While consistent conductivity values have been reported for scalp, brain, and cerebrospinal fluid, measured brain-to-skull conductivity ratio (BSCR) estimates have varied between 8 and 80, likely reflecting both inter-subject and measurement method differences. In simulations, mis-estimation of skull conductivity can produce source localization errors as large as 3cm. Here, we describe an iterative gradient-based approach to Simultaneous tissue Conductivity And source Location Estimation (SCALE). The scalp projection maps used by SCALE are obtained from near-dipolar effective EEG sources found by adequate independent component analysis (ICA) decomposition of sufficient high-density EEG data. We applied SCALE to simulated scalp projections of 15cm(2)-scale cortical patch sources in an MR image-based electrical head model with simulated BSCR of 30. Initialized either with a BSCR of 80 or 20, SCALE estimated BSCR as 32.6. In Adaptive Mixture ICA (AMICA) decompositions of (45-min, 128-channel) EEG data from two young adults we identified sets of 13 independent components having near-dipolar scalp maps compatible with a single cortical source patch. Again initialized with either BSCR 80 or 25, SCALE gave BSCR estimates of 34 and 54 for the two subjects respectively. The ability to accurately estimate skull conductivity non-invasively from any well-recorded EEG data in combination with a stable and non-invasively acquired MR imaging-derived electrical head model could remove a critical barrier to using EEG as a sub-cm(2)-scale accurate 3-D functional cortical imaging modality.

Research Review: Use of EEG biomarkers in child psychiatry research - current state and future directions

BACKGROUND

Electroencephalography (EEG) and related measures have a long and productive history in child psychopathology research and are currently experiencing a renaissance in interest, particularly for use as putative biomarkers.

METHOD AND SCOPE

First, the recent history leading to the use of EEG measures as endophenotypes and biomarkers for disease and treatment response are reviewed. Two key controversies within the area of noninvasive human electrophysiology research are discussed, and problems that currently either function as barriers or provide gateways to progress. First, the differences between the main types of EEG measurements (event-related potentials, quantitative EEG, and time-frequency measures) and how they can contribute collectively to better understanding of cortical dynamics underlying cognition and behavior are highlighted. Second, we focus on the ongoing shift in analytic focus to specific cortical sources and source networks whose dynamics are relevant to the clinical and experimental focus of the study, and the effective increase in source signal-to-noise ratio (SNR) that may be obtained in the process.

CONCLUSIONS

Understanding of these issues informs any discussion of current trends in EEG research. We highlight possible ways to evolve our understanding of brain dynamics beyond the apparent contradictions in understanding and modeling EEG activity highlighted by these controversies. Finally, we summarize some promising future directions of EEG biomarker research in child psychopathology.

Test-retest reliability of resting-state magnetoencephalography power in sensor and source space

Several studies have reported changes in spontaneous brain rhythms that could be used as clinical biomarkers or in the evaluation of neuropsychological and drug treatments in longitudinal studies using magnetoencephalography (MEG). There is an increasing necessity to use these measures in early diagnosis and pathology progression; however, there is a lack of studies addressing how reliable they are. Here, we provide the first test-retest reliability estimate of MEG power in resting-state at sensor and source space. In this study, we recorded 3 sessions of resting-state MEG activity from 24 healthy subjects with an interval of a week between each session. Power values were estimated at sensor and source space with beamforming for classical frequency bands: delta (2-4 Hz), theta (4-8 Hz), alpha (8-13 Hz), low beta (13-20 Hz), high beta (20-30 Hz), and gamma (30-45 Hz). Then, test-retest reliability was evaluated using the intraclass correlation coefficient (ICC). We also evaluated the relation between source power and the within-subject variability. In general, ICC of theta, alpha, and low beta power was fairly high (ICC > 0.6) while in delta and gamma power was lower. In source space, fronto-posterior alpha, frontal beta, and medial temporal theta showed the most reliable profiles. Signal-to-noise ratio could be partially responsible for reliability as low signal intensity resulted in high within-subject variability, but also the inherent nature of some brain rhythms in resting-state might be driving these reliability patterns. In conclusion, our results described the reliability of MEG power estimates in each frequency band, which could be considered in disease characterization or clinical trials.

Evaluation of adaptive PARAFAC alogorithms for tracking of simulated moving brain sources

In this paper, we proposed an online 2D localization method for tracking of dynamic moving brain sources. For this purpose, we used an adaptive version of PARAllel FACtor (PARAFAC) analysis for factorization of electroencephalographic (EEG) signals. We utilized Boundary Element Method (BEM) with four layers to solve the forward problem for the simulated EEG signals caused by two moving dipoles within the brain. Then, we created an appropriate tensor built by second order statistics of EEG signals. We adopted an online method to brain source localization called the Recursive Least Squares Tracking (RLST) as an adaptive PARAFAC algorithm with two windowing schemes. Finally, we evaluated the performance of the method applied to EEG signals.

EMDLAB: A toolbox for analysis of single-trial EEG dynamics using empirical mode decomposition

BACKGROUND

Empirical mode decomposition (EMD) is an empirical data decomposition technique. Recently there is growing interest in applying EMD in the biomedical field.

NEW METHOD

EMDLAB is an extensible plug-in for the EEGLAB toolbox, which is an open software environment for electrophysiological data analysis.

RESULTS

EMDLAB can be used to perform, easily and effectively, four common types of EMD: plain EMD, ensemble EMD (EEMD), weighted sliding EMD (wSEMD) and multivariate EMD (MEMD) on EEG data. In addition, EMDLAB is a user-friendly toolbox and closely implemented in the EEGLAB toolbox.

COMPARISON WITH EXISTING METHODS

EMDLAB gains an advantage over other open-source toolboxes by exploiting the advantageous visualization capabilities of EEGLAB for extracted intrinsic mode functions (IMFs) and Event-Related Modes (ERMs) of the signal.

CONCLUSIONS

EMDLAB is a reliable, efficient, and automated solution for extracting and visualizing the extracted IMFs and ERMs by EMD algorithms in EEG study.

On the fully automatic construction of a realistic head model for EEG source localization

Accurate multi-tissue segmentation of magnetic resonance (MR) images is an essential first step in the construction of a realistic finite element head conductivity model (FEHCM) for electroencephalography (EEG) source localization. All of the segmentation approaches proposed to date for this purpose require manual intervention or correction and are thus laborious, time-consuming, and subjective. In this paper we propose and evaluate a fully automatic method based on a hierarchical segmentation approach (HSA) incorporating Bayesian-based adaptive mean-shift segmentation (BAMS). An evaluation of HSA-BAMS, as well as two reference methods, in terms of both segmentation accuracy and the source localization accuracy of the resulting FEHCM is also presented. The evaluation was performed using (i) synthetic 2D multi-modal MRI head data and synthetic EEG (generated for a prescribed source), and (ii) real 3D T1-weighted MRI head data and real EEG data (with expert determined source localization). Expert manual segmentation served as segmentation ground truth. The results show that HSA-BAMS outperforms the two reference methods and that it can be used as a surrogate for manual segmentation for the construction of a realistic FEHCM for EEG source localization.

Fully automated whole-head segmentation with improved smoothness and continuity, with theory reviewed

Individualized current-flow models are needed for precise targeting of brain structures using transcranial electrical or magnetic stimulation (TES/TMS). The same is true for current-source reconstruction in electroencephalography and magnetoencephalography (EEG/MEG). The first step in generating such models is to obtain an accurate segmentation of individual head anatomy, including not only brain but also cerebrospinal fluid (CSF), skull and soft tissues, with a field of view (FOV) that covers the whole head. Currently available automated segmentation tools only provide results for brain tissues, have a limited FOV, and do not guarantee continuity and smoothness of tissues, which is crucially important for accurate current-flow estimates. Here we present a tool that addresses these needs. It is based on a rigorous Bayesian inference framework that combines image intensity model, anatomical prior (atlas) and morphological constraints using Markov random fields (MRF). The method is evaluated on 20 simulated and 8 real head volumes acquired with magnetic resonance imaging (MRI) at 1 mm3 resolution. We find improved surface smoothness and continuity as compared to the segmentation algorithms currently implemented in Statistical Parametric Mapping (SPM). With this tool, accurate and morphologically correct modeling of the whole-head anatomy for individual subjects may now be feasible on a routine basis. Code and data are fully integrated into SPM software tool and are made publicly available. In addition, a review on the MRI segmentation using atlas and the MRF over the last 20 years is also provided, with the general mathematical framework clearly derived.

EEG imaging of toddlers during dyadic turn-taking: Mu-rhythm modulation while producing or observing social actions

Contemporary active-EEG and EEG-imaging methods show particular promise for studying the development of action planning and social-action representation in infancy and early childhood. Action-related mu suppression was measured in eleven 3-year-old children and their mothers during a 'live,' largely unscripted social interaction. High-density EEG was recorded from children and synchronized with motion-captured records of children's and mothers' hand actions, and with video recordings. Independent Component Analysis (ICA) was used to separate brain and non-brain source signals in toddlers' EEG records. EEG source dynamics were compared across three kinds of epochs: toddlers' own actions (execution), mothers' actions (observation), and between-turn intervals (no action). Mu (6-9Hz) power was suppressed in left and right somatomotor cortex during both action execution and observation, as reflected by independent components of individual children's EEG data. These mu rhythm components were accompanied by beta-harmonic (~16Hz) suppression, similar to findings from adults. The toddlers' power spectrum and scalp density projections provide converging evidence of adult-like mu-suppression features. Mu-suppression components' source locations were modeled using an age-specific 4-layer forward head model. Putative sources clustered around somatosensory cortex, near the hand/arm region. The results demonstrate that action-locked, event-related EEG dynamics can be measured, and source-resolved, from toddlers during social interactions with relatively unrestricted social behaviors.