Imaging of current flow in the human head during transcranial electrical therapy

Machine learning-optimized non-invasive brain stimulation and treatment response classification for major depression

BACKGROUND/OBJECTIVES

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation intervention that shows promise as a potential treatment for depression. However, the clinical efficacy of tDCS varies, possibly due to individual differences in head anatomy affecting tDCS dosage. While functional changes in brain activity are more commonly reported in major depressive disorder (MDD), some studies suggest that subtle macroscopic structural differences, such as cortical thickness or brain volume reductions, may occur in MDD and could influence tDCS electric field (E-field) distributions. Therefore, accounting for individual anatomical differences may provide a pathway to optimize functional gains in MDD by formulating personalized tDCS dosage.

METHODS

To address the dosing variability of tDCS, we examined a subsample of sixteen active-tDCS participants' data from the larger ELECT clinical trial (NCT01894815). With this dataset, individualized neuroimaging-derived computational models of tDCS current were generated for (1) classifying treatment response, (2) elucidating essential stimulation features associated with treatment response, and (3) computing a personalized dose of tDCS to maximize the likelihood of treatment response in MDD.

RESULTS

In the ELECT trial, tDCS was superior to placebo (3.2 points [95% CI, 0.7 to 5.5; P = 0.01]). Our algorithm achieved over 90% overall accuracy in classifying treatment responders from the active-tDCS group (AUC = 0.90, F1 = 0.92, MCC = 0.79). Computed precision doses also achieved an average response likelihood of 99.981% and decreased dosing variability by 91.9%.

CONCLUSION

These findings support our previously developed precision-dosing method for a new application in psychiatry by optimizing the statistical likelihood of tDCS treatment response in MDD.

Precise and Rapid Whole-Head Segmentation from Magnetic Resonance Images of Older Adults using Deep Learning

Whole-head segmentation from Magnetic Resonance Images (MRI) establishes the foundation for individualized computational models using finite element method (FEM). This foundation paves the path for computer-aided solutions in fields, particularly in non-invasive brain stimulation. Most current automatic head segmentation tools are developed using healthy young adults. Thus, they may neglect the older population that is more prone to age-related structural decline such as brain atrophy. In this work, we present a new deep learning method called GRACE, which stands for General, Rapid, And Comprehensive whole-hEad tissue segmentation. GRACE is trained and validated on a novel dataset that consists of 177 manually corrected MR-derived reference segmentations that have undergone meticulous manual review. Each T1-weighted MRI volume is segmented into 11 tissue types, including white matter, grey matter, eyes, cerebrospinal fluid, air, blood vessel, cancellous bone, cortical bone, skin, fat, and muscle. To the best of our knowledge, this work contains the largest manually corrected dataset to date in terms of number of MRIs and segmented tissues. GRACE outperforms five freely available software tools and a traditional 3D U-Net on a five-tissue segmentation task. On this task, GRACE achieves an average Hausdorff Distance of 0.21, which exceeds the runner-up at an average Hausdorff Distance of 0.36. GRACE can segment a whole-head MRI in about 3 seconds, while the fastest software tool takes about 3 minutes. In summary, GRACE segments a spectrum of tissue types from older adults T1-MRI scans at favorable accuracy and speed. The trained GRACE model is optimized on older adult heads to enable high-precision modeling in age-related brain disorders. To support open science, the GRACE code and trained weights are made available online and open to the research community at https://github.com/lab-smile/GRACE.

Volumetric measurements of weak current-induced magnetic fields in the human brain at high resolution

PURPOSE

Clinical use of transcranial electrical stimulation (TES) requires accurate knowledge of the injected current distribution in the brain. MR current density imaging (MRCDI) uses measurements of the TES-induced magnetic fields to provide this information. However, sufficient sensitivity and image quality in humans in vivo has only been documented for single-slice imaging.

METHODS

A recently developed, optimally spoiled, acquisition-weighted, gradient echo-based 2D-MRCDI method has now been advanced for volume coverage with densely or sparsely distributed slices: The 3D rectilinear sampling (3D-DENSE) and simultaneous multislice acquisition (SMS-SPARSE) were optimized and verified by cable-loop experiments and tested with 1-mA TES experiments for two common electrode montages.

RESULTS

Comparisons between the volumetric methods against the 2D-MRCDI showed that relatively long acquisition times of 3D-DENSE using a single slab with six slices hindered the expected sensitivity improvement in the current-induced field measurements but improved sensitivity by 61% in the Laplacian of the field, on which some MRCDI reconstruction methods rely. Also, SMS-SPARSE acquisition of three slices, with a factor 2 CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) acceleration, performed best against the 2D-MRCDI with sensitivity improvements for the∆ B z , c $$ \Delta {B}_{z,c} $$ and Laplacian noise floors of 56% and 78% (baseline without current flow) as well as 43% and 55% (current injection into head). SMS-SPARSE reached a sensitivity of 67 pT for three distant slices at 2 × 2 × 3 mm3 resolution in 10 min of total scan time, and consistently improved image quality.

CONCLUSION

Volumetric MRCDI measurements with high sensitivity and image quality are well suited to characterize the TES field distribution in the human brain.

Transcranial MEP threshold voltages and current densities simulated with finite element modelling

OBJECTIVE

The aim of this study was to compare stimulation thresholds and current densities in the brain for transcranial motor evoked potentials (tcMEPs) from the hands and feet with linked quadripolar (LQP), M3-M4 and C1-C2 electrode montages.

METHODS

Twenty-five patients underwent cerebral vascular surgery with tcMEP monitoring. tcMEP voltage thresholds were compared between LQP (C1, M3, C2, M4), C1-C2, and M3-M4 montages. In a finite element model (FEM), hand, arm, and leg regions of interest (ROIs) on the cortical motor homunculus were segmented. Current densities in these ROIs at tcMEP thresholds were compared across tcMEP electrode montages.

RESULTS

LQP tcMEP thresholds were 61.5 volts for hands and 95.2 volts for feet. Thresholds were higher for M3-M4 (hands, 89.4 V; feet, 141.3 V) and C1-C2 (hands: 137.3 V; feet: 194.7 V). Total current at threshold voltage was greater for LQP (hands, 210.9 mA; feet, 311.3 mA) compared to M3-M4 (hands, 166.8 mA; feet, 256.6 mA), but similar to C1-C2 (hands, 246.7 mA; feet, 341.1 mA). In FEM simulations, current density and local current density topography in the hand ROI at threshold were very similar for LQP, M3-M4 and C1-C2.

CONCLUSIONS

TcMEP voltage thresholds were least for LQP, and lesser for M3-M4 compared to C1-C2. In FEM simulations, resistance to current to hand ROI was ordered the same (LQP < M3-M4 < C1-C2). The local distribution of current density in motor cortex with tcMEP was mainly determined by cortical geometry.

SIGNIFICANCE

Current densities and resistance to current simulated with FEM may explain threshold requirements for tcMEP electrode montages.

Estimation of brain tissue response by electrical stimulation in a subject-specific model implemented by conductivity tensor imaging

Electrical stimulation such as transcranial direct current stimulation (tDCS) is widely used to treat neuropsychiatric diseases and neurological disorders. Computational modeling is an important approach to understand the mechanisms underlying tDCS and optimize treatment planning. When applying computational modeling to treatment planning, uncertainties exist due to insufficient conductivity information inside the brain. In this feasibility study, we performed in vivo MR-based conductivity tensor imaging (CTI) experiments on the entire brain to precisely estimate the tissue response to the electrical stimulation. A recent CTI method was applied to obtain low-frequency conductivity tensor images. Subject-specific three-dimensional finite element models (FEMs) of the head were implemented by segmenting anatomical MR images and integrating a conductivity tensor distribution. The electric field and current density of brain tissues following electrical stimulation were calculated using a conductivity tensor-based model and compared to results using an isotropic conductivity model from literature values. The current density by the conductivity tensor was different from the isotropic conductivity model, with an average relative difference |rD| of 52 to 73%, respectively, across two normal volunteers. When applied to two tDCS electrode montages of C3-FP2 and F4-F3, the current density showed a focused distribution with high signal intensity which is consistent with the current flowing from the anode to the cathode electrodes through the white matter. The gray matter tended to carry larger amounts of current densities regardless of directional information. We suggest this CTI-based subject-specific model can provide detailed information on tissue responses for personalized tDCS treatment planning.

Neurophysiological mechanisms of transcranial alternating current stimulation

Neuronal oscillations are the primary basis for precise temporal coordination of neuronal processing and are linked to different brain functions. Transcranial alternating current stimulation (tACS) has demonstrated promising potential in improving cognition by entraining neural oscillations. Despite positive findings in recent decades, the results obtained are sometimes rife with variance and replicability problems, and the findings translation to humans is quite challenging. A thorough understanding of the mechanisms underlying tACS is necessitated for accurate interpretation of experimental results. Animal models are useful for understanding tACS mechanisms, optimizing parameter administration, and improving rational design for broad horizons of tACS. Here, we review recent electrophysiological advances in tACS from animal models, as well as discuss some critical issues for results coordination and translation. We hope to provide an overview of neurophysiological mechanisms and recommendations for future consideration to improve its validity, specificity, and reproducibility.

Inconsistencies in mapping current distribution in transcranial direct current stimulation

INTRODUCTION

tDCS is a non-invasive neuromodulation technique that has been widely studied both as a therapy for neuropsychiatric diseases and for cognitive enhancement. However, recent meta-analyses have reported significant inconsistencies amongst tDCS studies. Enhancing empirical understanding of current flow in the brain may help elucidate some of these inconsistencies.

METHODS

We investigated tDCS-induced current distribution by injecting a low frequency current waveform in a phantom and in vivo. MR phase images were collected during the stimulation and a time-series analysis was used to reconstruct the magnetic field. A current distribution map was derived from the field map using Ampere's law.

RESULTS

The current distribution map in the phantom showed a clear path of current flow between the two electrodes, with more than 75% of the injected current accounted for. However, in brain, the results did evidence a current path between the two target electrodes but only some portion ( 25%) of injected current reached the cortex demonstrating that a significant fraction of the current is bypassing the brain and traveling from one electrode to the other external to the brain, probably due to conductivity differences in brain tissue types. Substantial inter-subject and intra-subject (across consecutive scans) variability in current distribution maps were also observed in human but not in phantom scans.

DISCUSSIONS

An in-vivo current mapping technique proposed in this study demonstrated that much of the injected current in tDCS was not accounted for in human brain and deviated to the edge of the brain. These findings would have ramifications in the use of tDCS as a neuromodulator and may help explain some of the inconsistencies reported in other studies.

Phantom Construction and Equipment Configurations for Characterizing Electrical Properties Using MRI

Phantom objects are commonly employed in MRI systems as stable substitutes for biological tissues to ensure systems for measuring images are operating correctly and safely. For magnetic resonance electrical impedance tomography (MREIT) and magnetic resonance electrical property tomography (MREPT), conductivity or permittivity phantoms play an important role in checking MRI pulse sequences, MREIT equipment performance, and algorithm validation. The construction of these phantoms is explained in this chapter. In the first part, materials used for phantom construction are introduced. Ingredients for modifying the electromagnetic properties and relaxation times are presented, and the advantages and disadvantages of aqueous, gel, and hybrid conductivity phantoms are explained. The devices and methods used to confirm phantom electromagnetic properties are explained. Next, different types of MREIT electrode materials and the constant current sources used for MREIT studies are discussed. In the last section, we present the results of previous MREIT and MREPT studies.

On the reconstruction of magnetic resonance current density images of the human brain: Pitfalls and perspectives

Magnetic resonance current density imaging (MRCDI) of the human brain aims to reconstruct the current density distribution caused by transcranial electric stimulation from MR-based measurements of the current-induced magnetic fields. So far, the MRCDI data acquisition achieves only a low signal-to-noise ratio, does not provide a full volume coverage and lacks data from the scalp and skull regions. In addition, it is only sensitive to the component of the current-induced magnetic field parallel to the scanner field. The reconstruction problem thus involves coping with noisy and incomplete data, which makes it mathematically challenging. Most existing reconstruction methods have been validated using simulation studies and measurements in phantoms with simplified geometries. Only one reconstruction method, the projected current density algorithm, has been applied to human in-vivo data so far, however resulting in blurred current density estimates even when applied to noise-free simulated data. We analyze the underlying causes for the limited performance of the projected current density algorithm when applied to human brain data. In addition, we compare it with an approach that relies on the optimization of the conductivities of a small number of tissue compartments of anatomically detailed head models reconstructed from structural MR data. Both for simulated ground truth data and human in-vivo MRCDI data, our results indicate that the estimation of current densities benefits more from using a personalized volume conductor model than from applying the projected current density algorithm. In particular, we introduce a hierarchical statistical testing approach as a principled way to test and compare the quality of reconstructed current density images that accounts for the limited signal-to-noise ratio of the human in-vivo MRCDI data and the fact that the ground truth of the current density is unknown for measured data. Our results indicate that the statistical testing approach constitutes a valuable framework for the further development of accurate volume conductor models of the head. Our findings also highlight the importance of tailoring the reconstruction approaches to the quality and specific properties of the available data.

Sensitivity and resolution improvement for in vivo magnetic resonance current-density imaging of the human brain

PURPOSE

Magnetic resonance current-density imaging (MRCDI) combines MRI with low-intensity transcranial electrical stimulation (TES; 1-2 mA) to map current flow in the brain. However, usage of MRCDI is still hampered by low measurement sensitivity and image quality.

METHODS

Recently, a multigradient-echo-based MRCDI approach has been introduced that presently has the best-documented efficiency. This MRCDI approach has now been advanced in three directions and has been validated by phantom and in vivo experiments. First, the importance of optimum spoiling for brain imaging was verified. Second, the sensitivity and spatial resolution were improved by using acquisition weighting. Third, navigators were added as a quality control measure for tracking physiological noise. Combining these advancements, the optimized MRCDI method was tested by using 1 mA TES for two different injection profiles.

RESULTS

For a session duration of 4:20 min, the new MRCDI method was able to detect TES-induced magnetic fields at a sensitivity level of 84 picotesla, representing a twofold efficiency increase against our original method. A comparison between measurements and simulations based on personalized head models showed a consistent increase in the coefficient of determination of ΔR² = 0.12 for the current-induced magnetic fields and ΔR² = 0.22 for the current flow reconstructions. Interestingly, some of the simulations still clearly deviated from the measurements despite the strongly improved measurement quality. This highlights the utility of MRCDI to improve head models for TES simulations.

CONCLUSION

The achieved sensitivity improvement is an important step from proof-of-concept studies toward a broader application of MRCDI in clinical and basic neuroscience research.

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.

In-vivo imaging of targeting and modulation of depression-relevant circuitry by transcranial direct current stimulation: a randomized clinical trial

Recent clinical trials of transcranial direct current stimulation (tDCS) in depression have shown contrasting results. Consequently, we used in-vivo neuroimaging to confirm targeting and modulation of depression-relevant neural circuitry by tDCS. Depressed participants (N = 66, Baseline Hamilton Depression Rating Scale (HDRS) 17-item scores ≥14 and <24) were randomized into Active/Sham and High-definition (HD)/Conventional (Conv) tDCS groups using a double-blind, parallel design, and received tDCS individually targeted at the left dorsolateral prefrontal cortex (DLPFC). In accordance with Ampere's Law, tDCS currents were hypothesized to induce magnetic fields at the stimulation-target, measured in real-time using dual-echo echo-planar-imaging (DE-EPI) MRI. Additionally, the tDCS treatment trial (consisting of 12 daily 20-min sessions) was hypothesized to induce cerebral blood flow (CBF) changes post-treatment at the DLPFC target and in the reciprocally connected anterior cingulate cortex (ACC), measured using pseudo-continuous arterial spin labeling (pCASL) MRI. Significant tDCS current-induced magnetic fields were observed at the left DLPFC target for both active stimulation montages (Brodmann's area (BA) 46: p_HD = 0.048, Cohen's d_HD = 0.73; p_Conv = 0.018, d_Conv = 0.86; BA 9: p_HD = 0.011, d_HD = 0.92; p_Conv = 0.022, d_Conv = 0.83). Significant longitudinal CBF increases were observed (a) at the left DLPFC stimulation-target for both active montages (p_HD = 3.5E-3, d_HD = 0.98; p_Conv = 2.8E-3, d_Conv = 1.08), and (b) at ACC for the HD-montage only (p_HD = 2.4E-3, d_HD = 1.06; p_Conv = 0.075, d_Conv = 0.64). These results confirm that tDCS-treatment (a) engages the stimulation-target, and (b) modulates depression-relevant neural circuitry in depressed participants, with stronger network-modulations induced by the HD-montage. Although not primary outcomes, active HD-tDCS showed significant improvements of anhedonia relative to sham, though HDRS scores did not differ significantly between montages post-treatment.

Stimulating aged brains with transcranial direct current stimulation: Opportunities and challenges

Ageing involves significant neurophysiological changes that are both systematic while at the same time exhibiting divergent trajectories across individuals. These changes underlie cognitive impairments in elderly while also affecting the response of aged brains to interventions like transcranial direct current stimulation (tDCS). While the cognitive benefits of tDCS are more variable in elderly, older adults also respond differently to stimulation protocols compared to young adults. The age-related neurophysiological changes influencing the responsiveness to tDCS remain to be addressed in-depth. We review and discuss the premise that, in comparison to the better calibrated brain networks present in young adults, aged systems perform further away from a homoeostatic set-point. We argue that this age-related neurophysiological deviation from the homoeostatic optimum extends the leeway for tDCS to modulate the aged brain. This promotes the potency of immediate tDCS effects to induce directional plastic changes towards the homoeostatic equilibrium despite the impaired plasticity induction in elderly. We also consider how age-related neurophysiological changes pose specific challenges for tDCS that necessitate proper adaptations of stimulation protocols. Appreciating the distinctive properties of aged brains and the accompanying adjustment of stimulation parameters can increase the potency and reliability of tDCS as a treatment avenue in older adults.

Low frequency conductivity reconstruction based on a single current injection via MREIT

Conventional magnetic resonance electrical impedance tomography (MREIT) reconstruction methods require administration of two linearly independent currents via at least two electrode pairs. This requires long scanning times and inhibits coordination of MREIT measurements with electrical neuromodulation strategies. We sought to develop an isotropic conductivity reconstruction algorithm in MREIT based on a single current injection, both to decrease scanning time by a factor of two and enable MREIT measurements to be conveniently adapted to general transcranial- or implanted-electrode neurostimulation protocols. In this work, we propose and demonstrate an iterative algorithm that extends previously published MREIT work using two-current administration approaches. The proposed algorithm is a single-current adaptation of the harmonic B _z algorithm. Forward modeling of electric potentials is used to capture changes of conductivity along current directions that would normally be invisible using data from a single-current administration. Computational and experimental results show that the reconstruction algorithm is capable of reconstructing isotropic conductivity images that agree well in terms of L ² error and structural similarity with exact conductivity distributions or two-current-based MREIT reconstructions. We conclude that it is possible to reconstruct high quality electrical conductivity images using MREIT techniques and one current injection only.

Machine learning and individual variability in electric field characteristics predict tDCS treatment response

BACKGROUND

Transcranial direct current stimulation (tDCS) is widely investigated as a therapeutic tool to enhance cognitive function in older adults with and without neurodegenerative disease. Prior research demonstrates that electric current delivery to the brain can vary significantly across individuals. Quantification of this variability could enable person-specific optimization of tDCS outcomes. This pilot study used machine learning and MRI-derived electric field models to predict working memory improvements as a proof of concept for precision cognitive intervention.

METHODS

Fourteen healthy older adults received 20 minutes of 2 mA tDCS stimulation (F3/F4) during a two-week cognitive training intervention. Participants performed an N-back working memory task pre-/post-intervention. MRI-derived current models were passed through a linear Support Vector Machine (SVM) learning algorithm to characterize crucial tDCS current components (intensity and direction) that induced working memory improvements in tDCS responders versus non-responders.

MAIN RESULTS

SVM models of tDCS current components had 86% overall accuracy in classifying treatment responders vs. non-responders, with current intensity producing the best overall model differentiating changes in working memory performance. Median current intensity and direction in brain regions near the electrodes were positively related to intervention responses (r=0.811,p<0.001 and r=0.774,p=0.001).

CONCLUSIONS

This study provides the first evidence that pattern recognition analyses of MRI-derived tDCS current models can provide individual prognostic classification of tDCS treatment response with 86% accuracy. Individual differences in current intensity and direction play important roles in determining treatment response to tDCS. These findings provide important insights into mechanisms of tDCS response as well as proof of concept for future precision dosing models of tDCS intervention.

No Effect of Anodal tDCS on Verbal Episodic Memory Performance and Neurotransmitter Levels in Young and Elderly Participants

Healthy ageing is accompanied by cognitive decline that affects episodic memory processes in particular. Studies showed that anodal transcranial direct current stimulation (tDCS) to the left dorsolateral prefrontal cortex (DLPFC) may counteract this cognitive deterioration by increasing excitability and inducing neuroplasticity in the targeted cortical region. While stimulation gains are more consistent in initial low performers, relying solely on behavioural measures to predict treatment benefits does not suffice for a reliable implementation of this method as a therapeutic option. Hence, an exploration of the underlying neurophysiological mechanisms regarding the differential stimulation effect is warranted. Glutamatergic metabolites (Glx) and γ-aminobutyric acid (GABA) are involved in learning and memory processes and can be influenced with tDCS; wherefore, they present themselves as potential biomarkers for tDCS-induced behavioural gains, which are affiliated with neuroplasticity processes. In the present randomized, double-blind, sham-controlled, crossover study, 33 healthy young and 22 elderly participants received anodal tDCS to their left DLPFC during the encoding phase of a verbal episodic memory task. Using MEGA-PRESS edited magnetic resonance spectroscopy (MRS), Glx and GABA levels were measured in the left DLPFC before and after the stimulation period. Further, we tested whether baseline performance and neurotransmitter levels predicted subsequent gains. No beneficial group effects of tDCS emerged in either verbal retrieval performances or neurotransmitter concentrations. Moreover, baseline performance levels did not predict stimulation-induced cognitive gains, nor did Glx or GABA levels. Nevertheless, exploratory analyses suggested a predictive value of the Glx : GABA ratio, with lower ratios at baseline indicating greater tDCS-related gains in delayed recall performance. This highlights the importance of further studies investigating neurophysiological mechanisms underlying previously observed stimulation-induced cognitive benefits and their respective interindividual heterogeneity.

Modeling regional changes in dynamic stability during sleep and wakefulness

Global brain states are frequently placed within a unidimensional continuum by correlational studies, ranging from states of deep unconsciousness to ordinary wakefulness. An alternative is their multidimensional and mechanistic characterization in terms of different cognitive capacities, using computational models to reproduce the underlying neural dynamics. We explore this alternative by introducing a semi-empirical model linking regional activation and long-range functional connectivity in the different brain states visited during the natural wake-sleep cycle. Our model combines functional magnetic resonance imaging (fMRI) data, in vivo estimates of structural connectivity, and anatomically-informed priors to constrain the independent variation of regional activation. The best fit to empirical data was achieved using priors based on functionally coherent networks, with the resulting model parameters dividing the cortex into regions presenting opposite dynamical behavior. Frontoparietal regions approached a bifurcation from dynamics at a fixed point governed by noise, while sensorimotor regions approached a bifurcation from oscillatory dynamics. In agreement with human electrophysiological experiments, sleep onset induced subcortical deactivation with low correlation, which was subsequently reversed for deeper stages. Finally, we introduced periodic forcing of variable intensity to simulate external perturbations, and identified the key regions relevant for the recovery of wakefulness from deep sleep. Our model represents sleep as a state with diminished perceptual gating and the latent capacity for global accessibility that is required for rapid arousals. To the extent that the qualitative characterization of local dynamics is exhausted by the dichotomy between unstable and stable behavior, our work highlights how expanding the model parameter space can describe states of consciousness in terms of multiple dimensions with interpretations given by the choice of anatomically-informed priors.

A flexible workflow for simulating transcranial electric stimulation in healthy and lesioned brains

Simulating transcranial electric stimulation is actively researched as knowledge about the distribution of the electrical field is decisive for understanding the variability in the elicited stimulation effect. Several software pipelines comprehensively solve this task in an automated manner for standard use-cases. However, simulations for non-standard applications such as uncommon electrode shapes or the creation of head models from non-optimized T1-weighted imaging data and the inclusion of irregular structures are more difficult to accomplish. We address these limitations and suggest a comprehensive workflow to simulate transcranial electric stimulation based on open-source tools. The workflow covers the head model creation from MRI data, the electrode modeling, the modeling of anisotropic conductivity behavior of the white matter, the numerical simulation and visualization. Skin, skull, air cavities, cerebrospinal fluid, white matter, and gray matter are segmented semi-automatically from T1-weighted MR images. Electrodes of arbitrary number and shape can be modeled. The meshing of the head model is implemented in a way to preserve the feature edges of the electrodes and is free of topological restrictions of the considered structures of the head model. White matter anisotropy can be computed from diffusion-tensor imaging data. Our solver application was verified analytically and by contrasting the tDCS simulation results with that of other simulation pipelines (SimNIBS 3.0, ROAST 3.0). An agreement in both cases underlines the validity of our workflow. Our suggested solutions facilitate investigations of irregular structures in patients (e.g. lesions, implants) or new electrode types. For a coupled use of the described workflow, we provide documentation and disclose the full source code of the developed tools.

Modeling transcranial electrical stimulation in the aging brain

BACKGROUND

Varying treatment outcomes in transcranial electrical stimulation (tES) recipients may depend on the amount of current reaching the brain. Brain atrophy associated with normal aging may affect tES current delivery to the brain. Computational models have been employed to compute predicted tES current inside the brain. This study is the largest study that uses computational models to investigate tES field distribution in healthy older adults.

METHODS

Individualized head models from 587 healthy older adults (mean = 73.9years, 51-95 years) were constructed to create field maps. Two electrode montages (F3-F4, M1-SO) with 2 mA input current were modeled using ROAST with modified codes. A customized template of healthy older adults, the UFAB-587, was created from the same dataset and used to warp individual brains into the same space. Warped models were analyzed to determine the relationship between computed field measures, brain atrophy and age.

MAIN RESULTS

Computed field measures were inversely correlated with brain atrophy (R2 = 0.0829, p = 1.14e-12). Field pattern showed negative correlation with age in brain sub-regions including part of DLPFC and precentral gyrus. Mediation analysis revealed that the negative correlation between age and current density is partially mediated by brain-to-CSF ratio.

CONCLUSIONS

Computed field measures showed decreasing amount of tES current reaching the brain with increasing atrophy. Therefore, adjusting current dose by modifying tES stimulation parameters in older adults based on degree of atrophy may be necessary to achieve desired stimulation benefits. Results from this study may inform future tES application in healthy older adults.

Concurrent Imaging of Markers of Current Flow and Neurophysiological Changes During tDCS

Despite being a popular neuromodulation technique, clinical translation of transcranial direct current stimulation (tDCS) is hampered by variable responses observed within treatment cohorts. Addressing this challenge has been difficult due to the lack of an effective means of mapping the neuromodulatory electromagnetic fields together with the brain's response. In this study, we present a novel imaging technique that provides the capability of concurrently mapping markers of tDCS currents, as well as the brain's response to tDCS. A dual-echo echo-planar imaging (DE-EPI) sequence is used, wherein the phase of the acquired MRI-signal encodes the tDCS current induced magnetic field, while the magnitude encodes the blood oxygenation level dependent (BOLD) contrast. The proposed technique was first validated in a custom designed phantom. Subsequent test-retest experiments in human participants showed that tDCS-induced magnetic fields can be detected reliably in vivo. The concurrently acquired BOLD data revealed large-scale networks characteristic of a brain in resting-state as well as a 'cathodal' and an 'anodal' resting-state component under each electrode. Moreover, 'cathodal's BOLD-signal was observed to significantly decrease with the applied current at the group level in all datasets. With its ability to image markers of electromagnetic cause as well as neurophysiological changes, the proposed technique may provide an effective means to visualize neural engagement in tDCS at the group level. Our technique also contributes to addressing confounding factors in applying BOLD fMRI concurrently with tDCS.

Development and testing of implanted carbon electrodes for electromagnetic field mapping during neuromodulation

PURPOSE

Deep brain stimulation electrodes composed of carbon fibers were tested as a means of administering and imaging magnetic resonance electrical impedance tomography (MREIT) currents. Artifacts and heating properties of custom carbon-fiber deep brain stimulation (DBS) electrodes were compared with those produced with standard DBS electrodes.

METHODS

Electrodes were constructed from multiple strands of 7-μm carbon-fiber stock. The insulated carbon electrodes were matched to DBS electrode diameter and contact areas. Images of DBS and carbon electrodes were collected with and without current flow and were compared in terms of artifact and thermal effects in phantoms or tissue samples in 7T imaging conditions. Effects on magnetic flux density and current density distributions were also assessed.

RESULTS

Carbon electrodes produced magnitude artifacts with smaller FWHM values compared to the magnitude artifacts around DBS electrodes in spin echo and gradient echo imaging protocols. DBS electrodes appeared 269% larger than actual size in gradient echo images, in sharp contrast to the negligible artifact observed in diameter-matched carbon electrodes. As expected, larger temperature changes were observed near DBS electrodes during extended RF excitations compared with carbon electrodes in the same phantom. Magnitudes and distribution of magnetic flux density and current density reconstructions were comparable for carbon and DBS electrodes.

CONCLUSION

Carbon electrodes may offer a safer, MR-compatible method for administering neuromodulation currents. Use of carbon-fiber electrodes should allow imaging of structures close to electrodes, potentially allowing better targeting, electrode position revision, and the facilitation of functional imaging near electrodes during neuromodulation.

Montages for Noninvasive EEG Recording

Identifying the localization, distribution, and polarity of waveforms are the prime goals of clinical scalp EEG analysis. Appropriate choices of bipolar and referential montages are keys to emphasizing the diagnostic features of interest, and demand some understanding of the spatiotemporal physical behavior of the underlying neuronal generators. Several examples drawn from canonical epilepsy syndromes are used to illustrate this general message.

Methods to monitor accurate and consistent electrode placements in conventional transcranial electrical stimulation

BACKGROUND

Inaccurate electrode placement and electrode drift during a transcranial electrical stimulation (tES) session have been shown to alter predicted field distributions in the brain and thus may contribute to a large variation in tES study outcomes. Currently, there is no objective and independent measure to quantify electrode placement accuracy/drift in tES clinical studies.

OBJECTIVE/HYPOTHESIS

We proposed and tested novel methods to quantify accurate and consistent electrode placements in tES using models generated from a 3D scanner.

METHODS

Accurate electrode placements were quantified as Discrepancy in eight tES participants by comparing landmark distances of physical electrode locations F3/F4 to their model counterparts. Distances in models were computed using curve and linear based methods. Variability of landmark locations in a single subject was computed for multiple stimulation sessions to determine consistent electrode placements across four experimenters.

MAIN RESULTS

We obtained an average of 0.4 cm in Discrepancy, which was within the placement accuracy/drift threshold (1 cm) for conventional tES electrodes (∼35 cm2) to achieve reliable tES sessions suggested in the literature. Averaged Variability was 5.2%, with F4 electrode location as the least consistent placement.

CONCLUSIONS

These methods provide objective feedback for experimenters on their performance in placing tES electrodes. Applications of these methods can be used to monitor electrode locations in tES studies of a larger cohort using F3/F4 montage and other conventional electrode arrangements. Future studies may include co-registering the landmark locations with imaging-derived head models to quantify the effects of electrode accuracy/drift on predicted field distributions in the brain.

Benchmarking transcranial electrical stimulation finite element models: a comparison study

OBJECTIVE

To compare field measure differences in simulations of transcranial electrical stimulation (tES) generated by variations in finite element (FE) models due to boundary condition specification, use of tissue compartment smoothing filters, and use of free or structured tetrahedral meshes based on magnetic resonance imaging (MRI) data.

APPROACH

A structural MRI head volume was acquired at 1 mm³ resolution and segmented into ten tissue compartments. Predicted current densities and electric fields were computed in segmented models using modeling pipelines involving either an in-house (block) or a commercial platform commonly used in previous FE tES studies involving smoothed compartments and free meshing procedures (smooth). The same boundary conditions were used for both block and smooth pipelines. Differences caused by varying boundary conditions were examined using a simple geometry. Percentage differences of median current density values in five cortical structures were compared between the two pipelines for three electrode montages (F3-right supraorbital, T7-T8 and Cz-Oz).

MAIN RESULTS

Use of boundary conditions commonly used in previous tES FE studies produced asymmetric current density profiles in the simple geometry. In head models, median current density differences produced by the two pipelines, using the same boundary conditions, were up to 6% (isotropic) and 18% (anisotropic) in structures targeted by each montage. Tangential electric field measures calculated via either pipeline were within the range of values reported in the literature, when averaged over cortical surface patches.

SIGNIFICANCE

Apparently equivalent boundary settings may affect predicted current density outcomes and care must be taken in their specification. Smoothing FE model compartments may not be necessary, and directly translated, voxellated tissue boundaries at 1 mm³ resolution may be sufficient for use in tES FE studies, greatly reducing processing times. The findings here may be used to inform future current density modeling studies.

Methods to Compare Predicted and Observed Phosphene Experience in tACS Subjects

Background

Phosphene generation is an objective physical measure of potential transcranial alternating current stimulation (tACS) biological side effects. Interpretations from phosphene analysis can serve as a first step in understanding underlying mechanisms of tACS in healthy human subjects and assist validation of computational models.

Objective/Hypothesis

This preliminary study introduces and tests methods to analyze predicted phosphene occurrence using computational head models constructed from tACS recipients against verbal testimonies of phosphene sensations. Predicted current densities in the eyes and the occipital lobe were also verified against previously published threshold values for phosphenes.

Methods

Six healthy subjects underwent 10 Hz tACS while being imaged in an MRI scanner. Two different electrode montages, T7-T8 and Fpz-Oz, were used. Subject ratings of phosphene experience were collected during tACS and compared against current density distributions predicted in eye and occipital lobe regions of interest (ROIs) determined for each subject. Calculated median current densities in each ROI were compared to minimum thresholds for phosphene generation.

Main Results

All subjects reported phosphenes, and predicted median current densities in ROIs exceeded minimum thresholds for phosphenes found in the literature. Higher current densities in the eyes were consistently associated with decreased phosphene generation for the Fpz-Oz montage. There was an overall positive association between phosphene perceptions and current densities in the occipital lobe.

Conclusions

These methods may have promise for predicting phosphene generation using data collected during in-scanner tACS sessions and may enable better understanding of phosphene origin. Additional empirical data in a larger cohort is required to fully test the robustness of the proposed methods. Future studies should include additional montages that could dissociate retinal and occipital stimulation.

Evaluation of magnetohydrodynamic effects in magnetic resonance electrical impedance tomography at ultra-high magnetic fields

PURPOSE

Artifacts observed in experimental magnetic resonance electrical impedance tomography images were hypothesized to be because of magnetohydrodynamic (MHD) effects.

THEORY AND METHODS

Simulations of MREIT acquisition in the presence of MHD and electrical current flow were performed to confirm findings. Laminar flow and (electrostatic) electrical conduction equations were bidirectionally coupled via Lorentz force equations, and finite element simulations were performed to predict flow velocity as a function of time. Gradient sequences used in spin-echo and gradient echo acquisitions were used to calculate overall effects on MR phase images for different electrical current application or phase-encoding directions.

RESULTS

Calculated and experimental phase images agreed relatively well, both qualitatively and quantitatively, with some exceptions. Refocusing pulses in spin echo sequences did not appear to affect experimental phase images.

CONCLUSION

MHD effects were confirmed as the cause of observed experimental phase changes in MREIT images obtained at high fields. These findings may have implications for quantitative measurement of viscosity using MRI techniques. Methods developed here may be also important in studies of safety and in vivo artifacts observed in high field MRI systems.

Non-invasive Brain Stimulation: Probing Intracortical Circuits and Improving Cognition in the Aging Brain

The impact of cognitive aging on brain function and structure is complex, and the relationship between aging-related structural changes and cognitive function are not fully understood. Physiological and pathological changes to the aging brain are highly variable, making it difficult to estimate a cognitive trajectory with which to monitor the conversion to cognitive decline. Beyond the information on the structural and functional consequences of cognitive aging gained from brain imaging and neuropsychological studies, non-invasive brain stimulation techniques such as transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS) can enable stimulation of the human brain in vivo, offering useful insights into the functional integrity of intracortical circuits using electrophysiology and neuromodulation. TMS measurements can be used to identify and monitor changes in cortical reactivity, the integrity of inhibitory and excitatory intracortical circuits, the mechanisms of long-term potentiation (LTP)/depression-like plasticity and central cholinergic function. Repetitive TMS and tDCS can be used to modulate neuronal excitability and enhance cortical function, and thus offer a potential means to slow or reverse cognitive decline. This review will summarize and critically appraise relevant literature regarding the use of TMS and tDCS to probe cortical areas affected by the aging brain, and as potential therapeutic tools to improve cognitive function in the aging population. Challenges arising from intra-individual differences, limited reproducibility, and methodological differences will be discussed.

The effect of potassium chloride on Aplysia Californica abdominal ganglion activity

Objective

Spontaneous activity in the abdominal ganglion of Aplysia can be used as a convenient bioelectricity source in tests of novel MRI-based functional imaging methods, such as functional Magnetic Resonance Electrical Impedance Tomography (fMREIT). In these tests, it is necessary to find a consistent treatment that modulates neural activity, so that these results can be compared with control data. Effects of MREIT imaging currents combined with treatment were also of interest.

Approach

Potassium chloride (KCl) was employed as a rhythm modulator. In a series of experiments, effects of adding different volumes of KCl solution were tested and compared with experiments on control groups that had artificial sea water administered. In all cases, neuronal activity was measured with micro electrode arrays.

Main results

It was possible to reversibly stop spontaneous activity in ganglia by increasing the extracellular potassium chloride concentration to 89 mM. There was no effect on experimental outcomes when current was administered to the sample chamber between recordings.

Significance

KCl can be used as a reversible neural modulator for testing neural detection methods.

Low-Frequency Conductivity Tensor Imaging of the Human Head In Vivo Using DT-MREIT: First Study

We present the first in vivo images of anisotropic conductivity distribution in the human head, measured at a frequency of approximately 10 Hz. We used magnetic resonance electrical impedance tomography techniques to encode phase changes caused by current flow within the head via two independent electrode pairs. These results were then combined with diffusion tensor imaging data to reconstruct full anisotropic conductivity distributions in 5-mm-thick slices of the brains of two participants. Conductivity values recovered in this paper were broadly consistent with literature values. We anticipate that this technique will be of use in many areas of neuroscience, most importantly in functional imaging via inverse electroencephalogram. Future studies will involve pulse sequence acceleration to maximize brain coverage and resolution.