The influence of brain tissue anisotropy on human EEG and MEG

Rain may improve survival from direct lightning strikes to the human head

There is evidence that humans can survive a direct lightning strike to the head. Our question is: could water (rain) on the skin contribute to an increase in the survival rate? We measure the influence of rain during high-energy direct lightning strikes on a realistic three-compartment human head phantom. We find a lower number of perforations and eroded areas near the lightning strike impact points on the head phantom when rain was applied compared to no rain. Current amplitudes in the brain were lower with rain compared to no rain before a fully formed flashover. We conclude that rain on the scalp potentially contributes to the survival rate of 70-90% due to: (1) lower current exposition in the brain before a fully formed flashover, and (2) reduced mechanical and thermal damage.

Thalamic Segmentation and Neural Activation Modeling Based on Individual Tissue Microstructure in Deep Brain Stimulation for Essential Tremor

OBJECTIVE

Thalamic deep brain stimulation (DBS) is the primary surgical therapy for essential tremor (ET). Thalamic DBS traditionally uses an atlas-based targeting approach, which, although nominally accurate, may obscure individual anatomic differences from population norms. The objective of this study was to compare this traditional atlas-based approach with a novel quantitative modeling methodology grounded in individual tissue microstructure (N-of-1 approach).

MATERIALS AND METHODS

The N-of-1 approach uses individual patient diffusion tensor imaging (DTI) data to perform thalamic segmentation and volume of tissue activation (VTA) modeling. For each patient, the thalamus was individually segmented into 13 nuclei using DTI-based k-means clustering. DBS-induced VTAs associated with tremor suppression and side effects were then computed for each patient with finite-element electric-field models incorporating DTI microstructural data. Results from N-of-1 and traditional atlas-based modeling were compared for a large cohort of patients with ET treated with thalamic DBS.

RESULTS

The size and shape of individual N-of-1 thalamic nuclei and VTAs varied considerably across patients (N = 22). For both methods, tremor-improving therapeutic VTAs showed similar overlap with motor thalamic nuclei and greater motor than sensory nucleus overlap. For VTAs producing undesirable sustained paresthesia, 94% of VTAs overlapped with N-of-1 sensory thalamus estimates, whereas 74% of atlas-based segmentations overlapped. For VTAs producing dysarthria/motor contraction, the N-of-1 approach predicted greater spread beyond the thalamus into the internal capsule and adjacent structures than the atlas-based method.

CONCLUSIONS

Thalamic segmentation and VTA modeling based on individual tissue microstructure explain therapeutic stimulation equally well and side effects better than a traditional atlas-based method in DBS for ET. The N-of-1 approach may be useful in DBS targeting and programming, particularly when patient neuroanatomy deviates from population norms.

Sex-related patterns in the electroencephalogram and their relevance in machine learning classifiers

Deep learning is increasingly being proposed for detecting neurological and psychiatric diseases from electroencephalogram (EEG) data but the method is prone to inadvertently incorporate biases from training data and exploit illegitimate patterns. The recent demonstration that deep learning can detect the sex from EEG implies potential sex-related biases in deep learning-based disease detectors for the many diseases with unequal prevalence between males and females. In this work, we present the male- and female-typical patterns used by a convolutional neural network that detects the sex from clinical EEG (81% accuracy in a separate test set with 142 patients). We considered neural sources, anatomical differences, and non-neural artifacts as sources of differences in the EEG curves. Using EEGs from 1140 patients, we found electrocardiac artifacts to be leaking into the supposedly brain activity-based classifiers. Nevertheless, the sex remained detectable after rejecting heart-related and other artifacts. In the cleaned data, EEG topographies were critical to detect the sex, but waveforms and frequencies were not. None of the traditional frequency bands was particularly important for sex detection. We were able to determine the sex even from EEGs with shuffled time points and therewith completely destroyed waveforms. Researchers should consider neural and non-neural sources as potential origins of sex differences in their data, they should maintain best practices of artifact rejection, even when datasets are large, and they should test their classifiers for sex biases.

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.

Reconstruction and localization of auditory sources from intracerebral SEEG using independent component analysis

Stereo-electroencephalography (SEEG) is the surgical implantation of electrodes in the brain to better localize the epileptic network in pharmaco-resistant epileptic patients. This technique has exquisite spatial and temporal resolution. Still, the number and the position of the electrodes in the brain is limited and determined by the semiology and/or preliminary non-invasive examinations, leading to a large number of unexplored brain structures in each patient. Here, we propose a new approach to reconstruct the activity of non-sampled structures in SEEG, based on independent component analysis (ICA) and dipole source localization. We have tested this approach with an auditory stimulation dataset in ten patients. The activity directly recorded from the auditory cortex served as ground truth and was compared to the ICA applied on all non-auditory electrodes. Our results show that the activity from the auditory cortex can be reconstructed at the single trial level from contacts as far as ∼40 mm from the source. Importantly, this reconstructed activity is localized via dipole fitting in the proximity of the original source. In addition, we show that the size of the confidence interval of the dipole fitting is a good indicator of the reliability of the result, which depends on the geometry of the SEEG implantation. Overall, our approach allows reconstructing the activity of structures far from the electrode locations, partially overcoming the spatial sampling limitation of intracerebral recordings.

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

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

Bioelectromagnetism in Human Brain Research: New Applications, New Questions

Bioelectromagnetism has contributed some of the most commonly used techniques to human neuroscience such as magnetoencephalography (MEG), electroencephalography (EEG), transcranial magnetic stimulation (TMS), and transcranial electric stimulation (TES). The considerable differences in their technical design and practical use give rise to the impression that these are quite different techniques altogether. Here, we review, discuss and illustrate the fundamental principle of Helmholtz reciprocity that provides a common ground for all four techniques. We show that, more than 150 years after its discovery by Helmholtz in 1853, reciprocity is important to appreciate the strengths and limitations of these four classical tools in neuroscience. We build this case by explaining the concept of Helmholtz reciprocity, presenting a methodological account of this principle for all four methods and, finally, by illustrating its application in practical clinical studies.

Artificial neural network-based rapid predictor of biological nerve fiber activation for DBS applications

Objective.Computational models are powerful tools that can enable the optimization of deep brain stimulation (DBS). To enhance the clinical practicality of these models, their computational expense and required technical expertise must be minimized. An important aspect of DBS models is the prediction of neural activation in response to electrical stimulation. Existing rapid predictors of activation simplify implementation and reduce prediction runtime, but at the expense of accuracy. We sought to address this issue by leveraging the speed and generalization abilities of artificial neural networks (ANNs) to create a novel predictor of neural fiber activation in response to DBS.Approach.We developed six variations of an ANN-based predictor to predict the response of individual, myelinated axons to extracellular electrical stimulation. ANNs were trained using datasets generated from a finite-element model of an implanted DBS system together with multi-compartment cable models of axons. We evaluated the ANN-based predictors using three white matter pathways derived from group-averaged connectome data within a patient-specific tissue conductivity field, comparing both predicted stimulus activation thresholds and pathway recruitment across a clinically relevant range of stimulus amplitudes and pulse widths.Main results.The top-performing ANN could predict the thresholds of axons with a mean absolute error (MAE) of 0.037 V, and pathway recruitment with an MAE of 0.079%, across all parameters. The ANNs reduced the time required to predict the thresholds of 288 axons by four to five orders of magnitude when compared to multi-compartment cable models.Significance.We demonstrated that ANNs can be fast, accurate, and robust predictors of neural activation in response to DBS.

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

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

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

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

The Influence of the Head Model Conductor on the Source Localization of Auditory Evoked Potentials

The accuracy of EEG source analysis reconstruction improves when a realistic head volume conductor is modeled. In this study we investigated how the progressively more complex head representations influence the spatial localization of auditory-evoked potentials (AEPs). Fourteen young-adult participants with normal hearing performed the AEP task. Individualized head models were obtained from structural MRI and diffusion-weighted imaging scans collected in a separate session. AEPs were elicited by 1 k Hz and 4 k Hz tone bursts during a passive-listening tetanizing paradigm. We compared the amplitude of the N1 and P2 components before and after 4 min of tetanic-stimulation with 1 k Hz sounds. Current density reconstruction values of both components were investigated in the primary auditory cortex and adjacent areas. Furthermore, we compared the signal topography and magnitude obtained with 10 different head models on the EEG forward solution. Starting from the simplest model (scalp, skull, brain), we investigated the influence of modeling the CSF, distinguishing between GM and WM conductors, and including anisotropic WM values. We localized the activity of AEPs within the primary auditory cortex, but not in adjacent areas. The inclusion of the CSF compartment had the strongest influence on the source reconstruction, whereas white matter anisotropy led to a smaller improvement. We conclude that individualized realistic head models provide the best solution for the forward solution when modeling the CSF conductor.

Corticospinal excitability and conductivity are related to the anatomy of the corticospinal tract

Probing the brain structure-function relationship is at the heart of modern neuroscientific explorations, enabled by recent advances in brain mapping techniques. This study aimed to explore the anatomical blueprint of corticospinal excitability and shed light on the structure-function relationship within the human motor system. Using diffusion magnetic resonance imaging tractography, based on the spherical deconvolution approach, and transcranial magnetic stimulation (TMS), we show that anatomical inter-individual variability of the corticospinal tract (CST) modulates the corticospinal excitability and conductivity. Our findings show for the first time the relationship between increased corticospinal excitability and conductivity in individuals with a bigger CST (i.e., number of streamlines), as well as increased corticospinal microstructural organization (i.e., fractional anisotropy). These findings can have important implications for the understanding of the neuroanatomical basis of TMS as well as the study of the human motor system in both health and disease.

Evaluating Current Density Modeling of Non-Invasive Eye and Brain Electrical Stimulation Using Phosphene Thresholds

Because current flow cannot be measured directly in the intact retina or brain, current density distribution models were developed to estimate it during magnetic or electrical stimulation. A paradigm is now needed to evaluate if current flow modeling can be related to physiologically meaningful signs of true current distribution in the human brain. We used phosphene threshold measurements (PTs) as surrogate markers of current-flow to determine if PTs, evoked by transcranial alternating current stimulation (tACS), can be matched with current density estimates generated by head model-based computer simulations. Healthy, male subjects (n=15) were subjected to three-staged PT measurements comparing six unilateral and one bilateral stimulation electrode montages according to the 10/20 system: Fp2-Suborbital right (So), Fp2-right shoulder (rS), Fp2-Cz, Fp2- O2, So-rS, Cz-F8 and F7-F8. The stimulation frequency was set at 16 Hz. Subjects were asked to report the appearance and localization of phosphenes in their visual field for every montage. Current density models were built using multi-modal imaging data of a standard brain, meshed with isotropic conductivities of different tissues of the head using the SimBio and SCIRun software packages. We observed that lower PTs were associated with higher simulated current levels in the unilateral montages of the model head, and shorter electrode distances to the eye had lower PTs. The lowest mean PT and the lowest variability were found in the F7-F8 montage (95±33 μA). Our results confirm the hypothesis that phosphenes are primarily of retinal origin, and they provide the first in vivo evidence that computer models of current flow using head models are a valid tool to estimate real current flow in the human eye and brain.

An insight into electrical resistivity of white matter and brain tumors

BACKGROUND

There is a lack of information regarding electrical properties of white matter and brain tumors.

OBJECTIVE

To investigate the feasibility of in-vivo measurement of electrical resistivity during brain surgery and establish a better understanding of the resistivity patterns of brain tumors in correlation to the white matter.

METHODS

A bipolar probe was used to measure electrical resistivity during surgery in a prospective cohort of patients with brain tumors. For impedance measurement, the probe applied a constant current of 0.7 μA with a frequency of 140 Hz. The measurement was performed in the white matter within and outside peritumoral edema as well as in non-enhancing, enhancing and necrotic tumor areas. Resistivity values expressed in ohmmeter (Ω∗m) were compared between different intracranial tissues and brain tumors.

RESULTS

Ninety-two patients (gliomas WHO II:16, WHO III:10, WHO IV:33, metastasis:33) were included. White matter outside peritumoral edema had higher resistivity values (13.3 ± 1.7 Ω∗m) than within peritumoral edema (8.5 ± 1.6 Ω∗m), and both had higher values than brain tumors including non-enhancing (WHO II:6.4 ± 1.3 Ω∗m, WHO III:6.3 ± 0.9 Ω∗m), enhancing (WHO IV:5 ± 1 Ω∗m, metastasis:5.4 ± 1.3 Ω∗m) and necrotic tumor areas (WHO IV:3.9 ± 1.1 Ω∗m, metastasis:4.3 ± 1.3 Ω∗m), p=<0.001. No difference was found between low-grade and anaplastic gliomas, p = 0.808, while resistivity values in both were higher than the highest values found in glioblastomas, p = 0.003 and p = 0.004, respectively.

CONCLUSIONS

The technique we applied enabled us to measure electrical resistivity of white matter and brain tumors in-vivo presumably with a significant effect with regard to dielectric polarization. Our results suggest that there are significant differences within different areas and subtypes of brain tumors and that white matter exhibits higher electrical resistivity than brain tumors.

Magnetic-resonance-based measurement of electromagnetic fields and conductivity in vivo using single current administration-A machine learning approach

Diffusion tensor magnetic resonance electrical impedance tomography (DT-MREIT) is a newly developed technique that combines MR-based measurements of magnetic flux density with diffusion tensor MRI (DT-MRI) data to reconstruct electrical conductivity tensor distributions. DT-MREIT techniques normally require injection of two independent current patterns for unique reconstruction of conductivity characteristics. In this paper, we demonstrate an algorithm that can be used to reconstruct the position dependent scale factor relating conductivity and diffusion tensors, using flux density data measured from only one current injection. We demonstrate how these images can also be used to reconstruct electric field and current density distributions. Reconstructions were performed using a mimetic algorithm and simulations of magnetic flux density from complementary electrode montages, combined with a small-scale machine learning approach. In a biological tissue phantom, we found that the method reduced relative errors between single-current and two-current DT-MREIT results to around 10%. For in vivo human experimental data the error was about 15%. These results suggest that incorporation of machine learning may make it easier to recover electrical conductivity tensors and electric field images during neuromodulation therapy without the need for multiple current administrations.

High-frequency electrical properties tomography at 9.4T as a novel contrast mechanism for brain tumors

PURPOSE

To establish high-frequency magnetic resonance electrical properties tomography (MREPT) as a novel contrast mechanism for the assessment of glioblastomas using a rat brain tumor model.

METHODS

Six F98 intracranial tumor bearing rats were imaged longitudinally 8, 11 and 14 days after tumor cell inoculation. Conductivity and mean diffusivity maps were generated using MREPT and Diffusion Tensor Imaging. These maps were co-registered with T2 -weighted images and volumes of interests (VOIs) were segmented from the normal brain, ventricles, edema, viable tumor, tumor rim, and tumor core regions. Longitudinal changes in conductivity and mean diffusivity (MD) values were compared in these regions. A correlation analysis was also performed between conductivity and mean diffusivity values.

RESULTS

The conductivity of ventricles, edematous area and tumor regions (tumor rim, viable tumor, tumor core) was significantly higher (P < .01) compared to the contralateral cortex. The conductivity of the tumor increased over time while MD from the tumor did not change. A marginal positive correlation was noted between conductivity and MD values for tumor rim and viable tumor, whereas this correlation was negative for the tumor core.

CONCLUSION

We demonstrate a novel contrast mechanism based on ionic concentration and mobility, which may aid in providing complementary information to water diffusion in probing the microenvironment of brain tumors.

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

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

Brain Waves: Emergence of Localized, Persistent, Weakly Evanescent Cortical Loops

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of nonevanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full-brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the commonly agreed paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation because of the cortex wave modes in the highly folded areas will exhibit no apparent correlation with the fiber directions. Nonlinear coupling of those linear weakly evanescent wave modes then provides a universal mechanism for the emergence of synchronized brain wave field activity. The resonant and nonresonant terms of nonlinear coupling between multiple modes produce both synchronous spiking-like high-frequency wave activity as well as low-frequency wave rhythms. Numerical simulation of forced multiple-mode dynamics shows that, as forcing increases, there is a transition from damped to oscillatory regime that can then transition quickly to a nonoscillatory state when a critical excitation threshold is reached. The resonant nonlinear coupling results in the emergence of low-frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these weakly evanescent cortical wave modes have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and magnetoencephalography, and for the understanding of brain activity in general, including mechanisms of memory.

Universal theory of brain waves: from linear loops to nonlinear synchronized spiking and collective brain rhythms

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of non-evanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and to the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the standard paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation due to the cortex wave modes in highly folded areas will exhibit no apparent correlation with the fiber directions. The results are consistent with numerous recent experimental animal and human brain studies demonstrating the existence of electrostatic field activity in the form of traveling waves (including studies where neuronal connections were severed) and with wave loop induced peaks observed in EEG spectra. In addition, we demonstrate that the resonant and non-resonant terms of the nonlinear coupling between multiple modes produce both synchronous spiking-like high frequency wave activity as well as low frequency wave rhythms as a result of their unique dispersion properties. Numerical simulation of forced multiple mode dynamics shows that as forcing increases there is a transition from damped to oscillatory regime that subsequently decays away as over-excitation is reached. The resonant nonlinear coupling results in the emergence of low frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these cortical wave modes, and this new mechanism for understanding the nature of spiking behavior, have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and MEG, and in general for the understanding of brain activity, including mechanisms of memory.

The Study of Influence of Sound on Visual ERP-Based Brain Computer Interface

The performance of the event-related potential (ERP)-based brain–computer interface (BCI) declines when applying it into the real environment, which limits the generality of the BCI. The sound is a common noise in daily life, and whether it has influence on this decline is unknown. This study designs a visual-auditory BCI task that requires the subject to focus on the visual interface to output commands and simultaneously count number according to an auditory story. The story is played at three speeds to cause different workloads. Data collected under the same or different workloads are used to train and test classifiers. The results show that when the speed of playing the story increases, the amplitudes of P300 and N200 potentials decrease by 0.86 μV (p = 0.0239) and 0.69 μV (p = 0.0158) in occipital-parietal area, leading to a 5.95% decline (p = 0.0101) of accuracy and 9.53 bits/min decline (p = 0.0416) of information transfer rate. The classifier that is trained by the high workload data achieves higher accuracy than the one trained by the low workload if using the high workload data to test the performance. The result indicates that the sound could affect the visual ERP-BCI by increasing the workload. The large similarity of the training data and testing data is as important as the amplitudes of the ERP on obtaining high performance, which gives us an insight on how make to the ERP-BCI generalized.

Measurement and analysis of partial lightning currents in a head phantom

Direct lightning strikes to the human head can lead to various effects, ranging from burnings to death. The biological and physical mechanisms of a direct lightning strike in the human head are not well understood. The aim of this paper is to design an experimental setup to measure the spatial and temporal current distribution during a direct lightning strike to physical head phantoms to establish normative values for personal lightning protection equipment design and testing. We created head phantoms made of agarose, replicating the geometric and dielectric properties of scalp, skull, and intracranial volume. The bases of the three compartments were galvanically contacted via copper electrodes to measure the current per compartment. We used pulse generators to apply aperiodic voltage and current signals that modelled lightning components. Our experiments indicated that the scalp compartment was exposed to the current with a fraction of 80–90%. The brain and skull compartments were exposed between 6–13% and 3–6% of the total measured current respectively. In case of a flashover, most of the current (98–99%) flowed through the discharge channel. Unlike previous theoretical estimates and measurements in technical setups, we observed considerably longer times for the flashover to build up. In our experiments, the time to build up a fully formed flashover varied from approximately 30–700 μs. The observed current patterns in cases without and with flashover provided information on regions of possible damage in the human head. Consequently, we identified the phenomenon of a flashover as a potential mechanism for humans to survive a lightning strike. Our measured current distributions and amplitudes formed the base for normative values, which can be used in later experimental investigations regarding the possibilities of individual lightning protection equipment for humans.

Variation in Reported Human Head Tissue Electrical Conductivity Values

Electromagnetic source characterisation requires accurate volume conductor models representing head geometry and the electrical conductivity field. Head tissue conductivity is often assumed from previous literature, however, despite extensive research, measurements are inconsistent. A meta-analysis of reported human head electrical conductivity values was therefore conducted to determine significant variation and subsequent influential factors. Of 3121 identified publications spanning three databases, 56 papers were included in data extraction. Conductivity values were categorised according to tissue type, and recorded alongside methodology, measurement condition, current frequency, tissue temperature, participant pathology and age. We found variation in electrical conductivity of the whole-skull, the spongiform layer of the skull, isotropic, perpendicularly- and parallelly-oriented white matter (WM) and the brain-to-skull-conductivity ratio (BSCR) could be significantly attributed to a combination of differences in methodology and demographics. This large variation should be acknowledged, and care should be taken when creating volume conductor models, ideally constructing them on an individual basis, rather than assuming them from the literature. When personalised models are unavailable, it is suggested weighted average means from the current meta-analysis are used. Assigning conductivity as: 0.41 S/m for the scalp, 0.02 S/m for the whole skull, or when better modelled as a three-layer skull 0.048 S/m for the spongiform layer, 0.007 S/m for the inner compact and 0.005 S/m for the outer compact, as well as 1.71 S/m for the CSF, 0.47 S/m for the grey matter, 0.22 S/m for WM and 50.4 for the BSCR.

Handling anisotropic conductivities in the EEG forward problem with a symmetric formulation

The electroencephalography (EEG) forward problem, the computation of the electric potential generated by a known electric current source configuration in the brain, is a key step of EEG source analysis. In this problem, it is often desired to model the anisotropic conductivity profiles of the skull and of the white matter. These profiles, however, cannot be handled by standard surface integral formulations and the use of volume finite elements is required. Leveraging on the representation theorem using an anisotropic fundamental solution, this paper proposes a modified symmetric formulation for solving the EEG forward problem by a surface integral equation which can take into account anisotropic conductivity profiles. A set of numerical results is presented to corroborate theoretical treatments and to show the impact of the proposed approach on both canonical and real case scenarios.

A realistic, accurate and fast source modeling approach for the EEG forward problem

The aim of this paper is to advance electroencephalography (EEG) source analysis using finite element method (FEM) head volume conductor models that go beyond the standard three compartment (skin, skull, brain) approach and take brain tissue inhomogeneity (gray and white matter and cerebrospinal fluid) into account. The new approach should enable accurate EEG forward modeling in the thin human cortical structures and, more specifically, in the especially thin cortices in children brain research or in pathological applications. The source model should thus be focal enough to be usable in the thin cortices, but should on the other side be more realistic than the current standard mathematical point dipole. Furthermore, it should be numerically accurate and computationally fast. We propose to achieve the best balance between these demands with a current preserving (divergence conforming) dipolar source model. We develop and investigate a varying number of current preserving source basis elements n (n=1,…,n=5). For validation, we conducted numerical experiments within a multi-layered spherical domain, where an analytical solution exists. We show that the accuracy increases along with the number of basis elements, while focality decreases. The results suggest that the best balance between accuracy and focality in thin cortices is achieved with n=4 (or in extreme cases even n=3) basis functions, while in thicker cortices n=5 is recommended to obtain the highest accuracy. We also compare the current preserving approach to two further FEM source modeling techniques, namely partial integration and St. Venant, and show that the best current preserving source model outperforms the competing methods with regard to overall balance. For all tested approaches, FEM transfer matrices enable high computational speed. We implemented the new EEG forward modeling approaches into the open source duneuro library for forward modeling in bioelectromagnetism to enable its broader use by the brain research community. This library is build upon the DUNE framework for parallel finite elements simulations and integrates with high-level toolboxes like FieldTrip. Additionally, an inversion test has been implemented using the realistic head model to demonstrate and compare the differences between the aforementioned source models.

New Strategy for Finite Element Mesh Generation for Accurate Solutions of Electroencephalography Forward Problems

The finite element method (FEM) is a numerical method that is often used for solving electroencephalography (EEG) forward problems involving realistic head models. In this study, FEM solutions obtained using three different mesh structures, namely coarse, densely refined, and adaptively refined meshes, are compared. The simulation results showed that the accuracy of FEM solutions could be significantly enhanced by adding a small number of elements around regions with large estimated errors. Moreover, it was demonstrated that the adaptively refined regions were always near the current dipole sources, suggesting that selectively generating additional elements around the cortical surface might be a new promising strategy for more efficient FEM-based EEG forward analysis.

Spatial Localization of Sources in the Rat Subthalamic Motor Region Using an Inverse Current Source Density Method

Objective: In this study we introduce the use of the current source density (CSD) method as a way to visualize the spatial organization of evoked responses in the rat subthalamic nucleus (STN) at fixed time stamps resulting from motor cortex stimulation. This method offers opportunities to visualize neuronal input and study the relation between the synaptic input and the neural output of neural populations.

Approach: Motor cortex evoked local field potentials and unit activity were measured in the subthalamic region, with a 3D measurement grid consisting of 320 measurement points and high spatial resolution. This allowed us to visualize the evoked synaptic input by estimating the current source density (CSD) from the measured local field potentials, using the inverse CSD method. At the same time, the neuronal output of the cells within the grid is assessed by calculating post stimulus time histograms.

Main results: The CSD method resulted in clear and distinguishable sources and sinks of the neuronal input activity in the STN after motor cortex stimulation. We showed that the center of the synaptic input of the STN from the motor cortex is located dorsal to the input from globus pallidus.

Significance: For the first time we have performed CSD analysis on motor cortex stimulation evoked LFP responses in the rat STN as a proof of principle. Our results suggest that the CSD method can be used to gain new insights into the spatial extent of synaptic pathways in brain structures.

Overview of MEG

Magnetoencephalography (MEG) is a method to study electrical activity in the human brain by recording the neuromagnetic field outside the head. MEG, like electroencephalography (EEG), provides an excellent, millisecond-scale time resolution, and allows the estimation of the spatial distribution of the underlying activity, in favorable cases with a localization accuracy of a few millimeters. To detect the weak neuromagnetic signals, superconducting sensors, magnetically shielded rooms, and advanced signal processing techniques are used. The analysis and interpretation of MEG data typically involves comparisons between subject groups and experimental conditions using various spatial, temporal, and spectral measures of cortical activity and connectivity. The application of MEG to cognitive neuroscience studies is illustrated with studies of spoken language processing in subjects with normal and impaired reading ability. The mapping of spatiotemporal patterns of activity within networks of cortical areas can provide useful information about the functional architecture of the brain related to sensory and cognitive processing, including language, memory, attention, and perception.

In-vivo measurements of human brain tissue conductivity using focal electrical current injection through intracerebral multicontact electrodes

In-vivo measurements of human brain tissue conductivity at body temperature were conducted using focal electrical currents injected through intracerebral multicontact electrodes. A total of 1,421 measurements in 15 epileptic patients (age: 28 ± 10) using a radiofrequency generator (50 kHz current injection) were analyzed. Each contact pair was classified as being from healthy (gray matter, n = 696; white matter, n = 530) or pathological (epileptogenic zone, n = 195) tissue using neuroimaging analysis of the local tissue environment and intracerebral EEG recordings. Brain tissue conductivities were obtained using numerical simulations based on conductivity estimates that accounted for the current flow in the local brain volume around the contact pairs (a cube with a side length of 13 mm). Conductivity values were 0.26 S/m for gray matter and 0.17 S/m for white matter. Healthy gray and white matter had statistically different median impedances (P < 0.0001). White matter conductivity was found to be homogeneous as normality tests did not find evidence of multiple subgroups. Gray matter had lower conductivity in healthy tissue than in the epileptogenic zone (0.26 vs. 0.29 S/m; P = 0.012), even when the epileptogenic zone was not visible in the magnetic resonance image (MRI) (P = 0.005). The present in-vivo conductivity values could serve to create more accurate volume conduction models and could help to refine the identification of relevant intracerebral contacts, especially when located within the epileptogenic zone of an MRI-invisible lesion. Hum Brain Mapp 38:974-986, 2017. © 2016 Wiley Periodicals, Inc.

Semi-Automatic Segmentation of Optic Radiations and LGN, and Their Relationship to EEG Alpha Waves

At rest, healthy human brain activity is characterized by large electroencephalography (EEG) fluctuations in the 8-13 Hz range, commonly referred to as the alpha band. Although it is well known that EEG alpha activity varies across individuals, few studies have investigated how this may be related to underlying morphological variations in brain structure. Specifically, it is generally believed that the lateral geniculate nucleus (LGN) and its efferent fibres (optic radiation, OR) play a key role in alpha activity, yet it is unclear whether their shape or size variations contribute to its inter-subject variability. Given the widespread use of EEG alpha in basic and clinical research, addressing this is important, though difficult given the problems associated with reliably segmenting the LGN and OR. For this, we employed a multi-modal approach and combined diffusion magnetic resonance imaging (dMRI), functional magnetic resonance imaging (fMRI) and EEG in 20 healthy subjects to measure structure and function, respectively. For the former, we developed a new, semi-automated approach for segmenting the OR and LGN, from which we extracted several structural metrics such as volume, position and diffusivity. Although these measures corresponded well with known morphology based on previous post-mortem studies, we nonetheless found that their inter-subject variability was not significantly correlated to alpha power or peak frequency (p >0.05). Our results therefore suggest that alpha variability may be mediated by an alternative structural source and our proposed methodology may in general help in better understanding the influence of anatomy on function such as measured by EEG or fMRI.

Alternating Current Stimulation for Vision Restoration after Optic Nerve Damage: A Randomized Clinical Trial

BACKGROUND

Vision loss after optic neuropathy is considered irreversible. Here, repetitive transorbital alternating current stimulation (rtACS) was applied in partially blind patients with the goal of activating their residual vision.

METHODS

We conducted a multicenter, prospective, randomized, double-blind, sham-controlled trial in an ambulatory setting with daily application of rtACS (n = 45) or sham-stimulation (n = 37) for 50 min for a duration of 10 week days. A volunteer sample of patients with optic nerve damage (mean age 59.1 yrs) was recruited. The primary outcome measure for efficacy was super-threshold visual fields with 48 hrs after the last treatment day and at 2-months follow-up. Secondary outcome measures were near-threshold visual fields, reaction time, visual acuity, and resting-state EEGs to assess changes in brain physiology.

RESULTS

The rtACS-treated group had a mean improvement in visual field of 24.0% which was significantly greater than after sham-stimulation (2.5%). This improvement persisted for at least 2 months in terms of both within- and between-group comparisons. Secondary analyses revealed improvements of near-threshold visual fields in the central 5° and increased thresholds in static perimetry after rtACS and improved reaction times, but visual acuity did not change compared to shams. Visual field improvement induced by rtACS was associated with EEG power-spectra and coherence alterations in visual cortical networks which are interpreted as signs of neuromodulation. Current flow simulation indicates current in the frontal cortex, eye, and optic nerve and in the subcortical but not in the cortical regions.

CONCLUSION

rtACS treatment is a safe and effective means to partially restore vision after optic nerve damage probably by modulating brain plasticity. This class 1 evidence suggests that visual fields can be improved in a clinically meaningful way.

TRIAL REGISTRATION

ClinicalTrials.gov NCT01280877.

The role of blood vessels in high-resolution volume conductor head modeling of EEG

Reconstruction of the electrical sources of human EEG activity at high spatio-temporal accuracy is an important aim in neuroscience and neurological diagnostics. Over the last decades, numerous studies have demonstrated that realistic modeling of head anatomy improves the accuracy of source reconstruction of EEG signals. For example, including a cerebro-spinal fluid compartment and the anisotropy of white matter electrical conductivity were both shown to significantly reduce modeling errors. Here, we for the first time quantify the role of detailed reconstructions of the cerebral blood vessels in volume conductor head modeling for EEG. To study the role of the highly arborized cerebral blood vessels, we created a submillimeter head model based on ultra-high-field-strength (7T) structural MRI datasets. Blood vessels (arteries and emissary/intraosseous veins) were segmented using Frangi multi-scale vesselness filtering. The final head model consisted of a geometry-adapted cubic mesh with over 17×10(6) nodes. We solved the forward model using a finite-element-method (FEM) transfer matrix approach, which allowed reducing computation times substantially and quantified the importance of the blood vessel compartment by computing forward and inverse errors resulting from ignoring the blood vessels. Our results show that ignoring emissary veins piercing the skull leads to focal localization errors of approx. 5 to 15mm. Large errors (>2cm) were observed due to the carotid arteries and the dense arterial vasculature in areas such as in the insula or in the medial temporal lobe. Thus, in such predisposed areas, errors caused by neglecting blood vessels can reach similar magnitudes as those previously reported for neglecting white matter anisotropy, the CSF or the dura - structures which are generally considered important components of realistic EEG head models. Our findings thus imply that including a realistic blood vessel compartment in EEG head models will be helpful to improve the accuracy of EEG source analyses particularly when high accuracies in brain areas with dense vasculature are required.

Tractography Activation Patterns in Dorsolateral Prefrontal Cortex Suggest Better Clinical Responses in OCD DBS

Background: Medication resistant obsessive-compulsive disorder (OCD) patients can be successfully treated with Deep Brain Stimulation (DBS) which targets the anterior limb of the internal capsule (ALIC) and the nucleus accumbens (NA). Growing evidence suggests that in patients who respond to DBS, axonal fiber bundles surrounding the electrode are activated, but it is currently unknown which discrete pathways are critical for optimal benefit. Our aim was to identify axonal pathways mediating clinical effects of ALIC-NA DBS.

Methods: We created computational models of ALIC-NA DBS to simulate the activation of fiber tracts and to identify connected cerebral regions. The pattern of activated axons and their cortical targets was investigated in six OCD patients who underwent ALIC-NA DBS.

Results: Modulation of the right anterior middle frontal gyrus (dorsolateral prefrontal cortex) was associated with an excellent response. In contrast, non-responders showed high activation in the orbital part of the right inferior frontal gyrus (lateral orbitofrontal cortex/anterior ventrolateral prefrontal cortex). Factor analysis followed by step-wise linear regression indicated that YBOCS improvement was inversely associated with factors that were predominantly determined by gray matter activation results.

Discussion: Our findings support the hypothesis that optimal therapeutic results are associated with the activation of distinct fiber pathways. This suggests that in DBS for OCD, focused stimulation of specific fiber pathways, which would allow for stimulation with lower amplitudes, may be superior to activation of a wide array of pathways, typically associated with higher stimulation amplitudes.

Data-driven model comparing the effects of glial scarring and interface interactions on chronic neural recordings in non-human primates

OBJECTIVE

We characterized electrode stability over twelve weeks of impedance and neural recording data from four chronically-implanted Utah arrays in two rhesus macaques, and investigated the effects of glial scarring and interface interactions at the electrode recording site on signal quality using a computational model.

APPROACH

A finite-element model of a Utah array microelectrode in neural tissue was coupled with a multi-compartmental model of a neuron to quantify the effects of encapsulation thickness, encapsulation resistivity, and interface resistivity on electrode impedance and waveform amplitude. The coupled model was then reconciled with the in vivo data. Histology was obtained seventeen weeks post-implantation to measure gliosis.

MAIN RESULTS

From week 1-3, mean impedance and amplitude increased at rates of 115.8 kΩ/week and 23.1 μV/week, respectively. This initial ramp up in impedance and amplitude was observed across all arrays, and is consistent with biofouling (increasing interface resistivity) and edema clearing (increasing tissue resistivity), respectively, in the model. Beyond week 3, the trends leveled out. Histology showed that thin scars formed around the electrodes. In the model, scarring could not match the in vivo data. However, a thin interface layer at the electrode tip could. Despite having a large effect on impedance, interface resistivity did not have a noticeable effect on amplitude.

SIGNIFICANCE

This study suggests that scarring does not cause an electrical problem with regard to signal quality since it does not appear to be the main contributor to increasing impedance or significantly affect amplitude unless it displaces neurons. This, in turn, suggests that neural signals can be obtained reliably despite scarring as long as the recording site has sufficiently low impedance after accumulating a thin layer of biofouling. Therefore, advancements in microelectrode technology may be expedited by focusing on improvements to the recording site-tissue interface rather than elimination of the glial scar.

Effects of head models and dipole source parameters on EEG fields

Head model and an efficient method for computing the forward EEG (electroencephalography)problem are essential to dipole source localization(DSL). In this paper, we use less expensive ovoid geometry to approximate human head, aiming at investigating the effects of head shape and dipole source parameters on EEG fields. The application of point least squares (PLS) based on meshless method was introduced for solving EEG forward problem and numerical simulation is implemented in three kinds of ovoid head models. We present the performances of the surface potential in the face of varying dipole source parameters in detail. The results show that the potential patterns are similar for different dipole position in different head shapes, but the peak value of potential is significantly influenced by the head shape. Dipole position induces a great effect on the peak value of potential and shift of peak potential. The degree of variation between sphere head model and non-sphere head models is seen at the same time. We also show that PLS method with the trigonometric basis is superior to the constant basis, linear basis, and quadratic basis functions in accuracy and efficiency.

Influence of the head model on EEG and MEG source connectivity analyses

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problems. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only the skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC.

Visualizing simulated electrical fields from electroencephalography and transcranial electric brain stimulation: a comparative evaluation

Electrical activity of neuronal populations is a crucial aspect of brain activity. This activity is not measured directly but recorded as electrical potential changes using head surface electrodes (electroencephalogram - EEG). Head surface electrodes can also be deployed to inject electrical currents in order to modulate brain activity (transcranial electric stimulation techniques) for therapeutic and neuroscientific purposes. In electroencephalography and noninvasive electric brain stimulation, electrical fields mediate between electrical signal sources and regions of interest (ROI). These fields can be very complicated in structure, and are influenced in a complex way by the conductivity profile of the human head. Visualization techniques play a central role to grasp the nature of those fields because such techniques allow for an effective conveyance of complex data and enable quick qualitative and quantitative assessments. The examination of volume conduction effects of particular head model parameterizations (e.g., skull thickness and layering), of brain anomalies (e.g., holes in the skull, tumors), location and extent of active brain areas (e.g., high concentrations of current densities) and around current injecting electrodes can be investigated using visualization. Here, we evaluate a number of widely used visualization techniques, based on either the potential distribution or on the current-flow. In particular, we focus on the extractability of quantitative and qualitative information from the obtained images, their effective integration of anatomical context information, and their interaction. We present illustrative examples from clinically and neuroscientifically relevant cases and discuss the pros and cons of the various visualization techniques.

A guideline for head volume conductor modeling in EEG and MEG

For accurate EEG/MEG source analysis it is necessary to model the head volume conductor as realistic as possible. This includes the distinction of the different conductive compartments in the human head. In this study, we investigated the influence of modeling/not modeling the conductive compartments skull spongiosa, skull compacta, cerebrospinal fluid (CSF), gray matter, and white matter and of the inclusion of white matter anisotropy on the EEG/MEG forward solution. Therefore, we created a highly realistic 6-compartment head model with white matter anisotropy and used a state-of-the-art finite element approach. Starting from a 3-compartment scenario (skin, skull, and brain), we subsequently refined our head model by distinguishing one further of the above-mentioned compartments. For each of the generated five head models, we measured the effect on the signal topography and signal magnitude both in relation to a highly resolved reference model and to the model generated in the previous refinement step. We evaluated the results of these simulations using a variety of visualization methods, allowing us to gain a general overview of effect strength, of the most important source parameters triggering these effects, and of the most affected brain regions. Thereby, starting from the 3-compartment approach, we identified the most important additional refinement steps in head volume conductor modeling. We were able to show that the inclusion of the highly conductive CSF compartment, whose conductivity value is well known, has the strongest influence on both signal topography and magnitude in both modalities. We found the effect of gray/white matter distinction to be nearly as big as that of the CSF inclusion, and for both of these steps we identified a clear pattern in the spatial distribution of effects. In comparison to these two steps, the introduction of white matter anisotropy led to a clearly weaker, but still strong, effect. Finally, the distinction between skull spongiosa and compacta caused the weakest effects in both modalities when using an optimized conductivity value for the homogenized compartment. We conclude that it is highly recommendable to include the CSF and distinguish between gray and white matter in head volume conductor modeling. Especially for the MEG, the modeling of skull spongiosa and compacta might be neglected due to the weak effects; the simplification of not modeling white matter anisotropy is admissible considering the complexity and current limitations of the underlying modeling approach.

On conductivity, permittivity, apparent diffusion coefficient, and their usefulness as cancer markers at MRI frequencies

PURPOSE

To investigate the permittivity and conductivity of cancerous and normal tissues, their correlation to the apparent diffusion coefficient (ADC), and the specificity that they could add to cancer detection.

THEORY

Breast and prostate carcinomas were induced in rats. Conductivity and permittivity measurements were performed in the anesthetized animals using a dielectric probe and an impedance analyzer between 50 and 270 MHz. The correlations between ADCs (measured at 128 MHz) and conductivity values were investigated. Frequency-dependent discriminant functions were computed to assess the value that each parameter adds to cancer detection.

METHODS

Tumors exhibited higher permittivity than muscle tissue by 27%/12%/5% at 64/128/270MHz. Frequency independent, 15-20% higher conductivity was also noted in tumors compared to muscle tissue over the same frequency range. Strong negative correlation was observed between tissue conductivity and ADC. Whereas permittivity had the strongest discriminatory power at 64 MHz, it became comparable to ADC at 128 MHz and less important than ADC at 270 MHz.

CONCLUSION

Conductivity measurements offered limited advantages in separating cancer from normal tissue beyond what ADC already provided; conversely, permittivity added separation power when added to the discriminant function. The moderately high cancerous tissue permittivity and conductivity impose strong constraints on the capability of MRI-based tissue electrical property measurements.

Anisotropic conductivity tensor imaging in MREIT using directional diffusion rate of water molecules

Magnetic resonance electrical impedance tomography (MREIT) is an emerging method to visualize electrical conductivity and/or current density images at low frequencies (below 1 KHz). Injecting currents into an imaging object, one component of the induced magnetic flux density is acquired using an MRI scanner for isotropic conductivity image reconstructions. Diffusion tensor MRI (DT-MRI) measures the intrinsic three-dimensional diffusion property of water molecules within a tissue. It characterizes the anisotropic water transport by the effective diffusion tensor. Combining the DT-MRI and MREIT techniques, we propose a novel direct method for absolute conductivity tensor image reconstructions based on a linear relationship between the water diffusion tensor and the electrical conductivity tensor. We first recover the projected current density, which is the best approximation of the internal current density one can obtain from the measured single component of the induced magnetic flux density. This enables us to estimate a scale factor between the diffusion tensor and the conductivity tensor. Combining these values at all pixels with the acquired diffusion tensor map, we can quantitatively recover the anisotropic conductivity tensor map. From numerical simulations and experimental verifications using a biological tissue phantom, we found that the new method overcomes the limitations of each method and successfully reconstructs both the direction and magnitude of the conductivity tensor for both the anisotropic and isotropic regions.

Linking a neural mass model with a 3D model of the human brain to reproduce EEG signals

Electroencephalography (EEG) is often employed to measure electrical activity in the living human brain. Simulation studies can help unravel how the brain electrical activity pattern generates the EEG signal, still a widely unresolved question. This article describes a method to simulate brain electrical activity by using neuronal populations of a neural mass model. Implementing these populations in a finite element model of the head offers the opportunity to investigate the influence of each group of neurons to the scalp potential. This model is based on structural magnetic resonance imaging data to specify tissue composition, and diffusion tensor imaging data to model local anisotropy. We simulated the EEG signals of five neuronal populations generating α waves in the visual cortex. Our results indicate that radially oriented sources dominate over tangential sources in the generation of the scalp signal. Investigating the influence of anisotropic conductivity, we found small differences in topography and phase and larger ones for the potential amplitude compared with an isotropic conductivity distribution. The outcome of this article is a fast method based on superposition of sources for simulating time-dependent EEG signals, which can be used for further studies of neurodegenerative diseases.

Transcranial current brain stimulation (tCS): models and technologies

In this paper, we provide a broad overview of models and technologies pertaining to transcranial current brain stimulation (tCS), a family of related noninvasive techniques including direct current (tDCS), alternating current (tACS), and random noise current stimulation (tRNS). These techniques are based on the delivery of weak currents through the scalp (with electrode current intensity to area ratios of about 0.3-5 A/m2) at low frequencies (typically < 1 kHz) resulting in weak electric fields in the brain (with amplitudes of about 0.2-2 V/m). Here we review the biophysics and simulation of noninvasive, current-controlled generation of electric fields in the human brain and the models for the interaction of these electric fields with neurons, including a survey of in vitro and in vivo related studies. Finally, we outline directions for future fundamental and technological research.

Targeting of white matter tracts with transcranial magnetic stimulation

BACKGROUND

TMS activations of white matter depend not only on the distance from the coil, but also on the orientation of the axons relative to the TMS-induced electric field, and especially on axonal bends that create strong local field gradient maxima. Therefore, tractography contains potentially useful information for TMS targeting.

OBJECTIVE/METHODS

Here, we utilized 1-mm resolution diffusion and structural T1-weighted MRI to construct large-scale tractography models, and localized TMS white matter activations in motor cortex using electromagnetic forward modeling in a boundary element model (BEM).

RESULTS

As expected, in sulcal walls, pyramidal cell axonal bends created preferred sites of activation that were not found in gyral crowns. The model agreed with the well-known coil orientation sensitivity of motor cortex, and also suggested unexpected activation distributions emerging from the E-field and tract configurations. We further propose a novel method for computing the optimal coil location and orientation to maximally stimulate a pre-determined axonal bundle.

CONCLUSIONS

Diffusion MRI tractography with electromagnetic modeling may improve spatial specificity and efficacy of TMS.

Subspace electrode selection methodology for EEG multiple source localization error reduction due to uncertain conductivity values

This paper proposes a modification of the subspace correlation cost function and the Recursively Applied and Projected Multiple Signal Classification (RAP-MUSIC) method for electroencephalography (EEG) source analysis in epilepsy. This enables to reconstruct neural source locations and orientations that are less degraded due to the uncertain knowledge of the head conductivity values. An extended linear forward model is used in the subspace correlation cost function that incorporates the sensitivity of the EEG potentials to the uncertain conductivity value parameter. More specifically, the principal vector of the subspace correlation function is used to provide relevant information for solving the EEG inverse problems. A simulation study is carried out on a simplified spherical head model with uncertain skull to soft tissue conductivity ratio. Results show an improvement in the reconstruction accuracy of source parameters compared to traditional methodology, when using conductivity ratio values that are different from the actual conductivity ratio.

Particle Swarm Optimization Applied to EEG Source Localization of Somatosensory Evoked Potentials

One of the most important steps in presurgical diagnosis of medically intractable epilepsy is to find the precise location of the epileptogenic foci. Electroencephalography (EEG) is a noninvasive tool commonly used at epilepsy surgery centers for presurgical diagnosis. In this paper, a modified particle swarm optimization (MPSO) method is used to solve the EEG source localization problem. The method is applied to noninvasive EEG recording of somatosensory evoked potentials (SEPs) for a healthy subject. A 1 mm hexahedra finite element volume conductor model of the subject's head was generated using T1-weighted magnetic resonance imaging data. Special consideration was made to accurately model the skull and cerebrospinal fluid. An exhaustive search pattern and the MPSO method were then applied to the peak of the averaged SEP data and both identified the same region of the somatosensory cortex as the location of the SEP source. A clinical expert independently identified the expected source location, further corroborating the source analysis methods. The MPSO converged to the global minima with significantly lower computational complexity compared to the exhaustive search method that required almost 3700 times more evaluations.

Electroencephalographic inverse localization of brain activity in acute traumatic brain injury as a guide to surgery, monitoring and treatment

OBJECTIVE

To inverse-localize epileptiform cortical electrical activity recorded from severe traumatic brain injury (TBI) patients using electroencephalography (EEG).

METHODS

Three acute TBI cases were imaged using computed tomography (CT) and multimodal magnetic resonance imaging (MRI). Semi-automatic segmentation was performed to partition the complete TBI head into 25 distinct tissue types, including 6 tissue types accounting for pathology. Segmentations were employed to generate a finite element method model of the head, and EEG activity generators were modeled as dipolar currents distributed over the cortical surface.

RESULTS

We demonstrate anatomically faithful localization of EEG generators responsible for epileptiform discharges in severe TBI. By accounting for injury-related tissue conductivity changes, our work offers the most realistic implementation currently available for the inverse estimation of cortical activity in TBI.

CONCLUSION

Whereas standard localization techniques are available for electrical activity mapping in uninjured brains, they are rarely applied to acute TBI. Modern models of TBI-induced pathology can inform the localization of epileptogenic foci, improve surgical efficacy, contribute to the improvement of critical care monitoring and provide guidance for patient-tailored treatment. With approaches such as this, neurosurgeons and neurologists can study brain activity in acute TBI and obtain insights regarding injury effects upon brain metabolism and clinical outcome.

A coaxial optrode as multifunction write-read probe for optogenetic studies in non-human primates

BACKGROUND

Advances in optogenetics have led to first reports of expression of light-gated ion-channels in non-human primates (NHPs). However, a major obstacle preventing effective application of optogenetics in NHPs and translation to optogenetic therapeutics is the absence of compatible multifunction optoelectronic probes for (1) precision light delivery, (2) low-interference electrophysiology, (3) protein fluorescence detection, and (4) repeated insertion with minimal brain trauma.

NEW METHOD

Here we describe a novel brain probe device, a "coaxial optrode", designed to minimize brain tissue damage while microfabricated to perform simultaneous electrophysiology, light delivery and fluorescence measurements in the NHP brain. The device consists of a tapered, gold-coated optical fiber inserted in a polyamide tube. A portion of the gold coating is exposed at the fiber tip to allow electrophysiological recordings in addition to light delivery/collection at the tip.

RESULTS

Coaxial optrode performance was demonstrated by experiments in rodents and NHPs, and characterized by computational models. The device mapped opsin expression in the brain and achieved precisely targeted optical stimulation and electrophysiology with minimal cortical damage.

COMPARISON WITH EXISTING METHODS

Overall, combined electrical, optical and mechanical features of the coaxial optrode allowed a performance for NHP studies which was not possible with previously existing devices.

CONCLUSIONS

Coaxial optrode is currently being used in two NHP laboratories as a major tool to study brain function by inducing light modulated neural activity and behavior. By virtue of its design, the coaxial optrode can be extended for use as a chronic implant and multisite neural stimulation/recording.