2417

Magnetic flux density comparisons between in-vivo TACS human Magnetic Resonance Electrical Impedance Tomography measurements and MRI-derived human computational models
Aprinda Indahlastari1, Aditya Kumar Kasinadhuni2, Kevin Castellano2, Christopher Saar1, Casey Weigel1, Bakir Mousa1, Munish Chauhan1, Thomas H Mareci2, and Rosalind J Sadleir1

1Arizona State University, Tempe, AZ, United States, 2University of Florida

Synopsis

actual current densities applied to the brain in neuromodulation therapies. MREIT can be used to determine actual current densities delivered to the brain by measuring Bz. Here we present Bz distributions resulting from trans-temporal current injection obtained from MREIT imaging during 10 Hz TACS in a healthy human subject. Comparisons of MREIT results to MRI-derived computational models suggested that actual contact areas between electrodes and scalp may be smaller than electrode surface areas, and in-vivo tissue conductivity values, particularly skin and skull, may be different than assumed.

Purpose

Effects of changing electrode-scalp contact areas and tissue conductivity values in TACS finite element models were compared with in-vivo human MREIT results to determine which parameters produced the best qualitative agreement in magnetic flux density data. While MREIT can be used to reconstruct conductivity distributions as well as current density patterns, MREIT phase data quality in skin and bone is low because of low signal and the narrow structure of these compartments. If saline soaked sponges are used to make electrical contact with the scalp it is difficult to determine contact quality a priori, and there is some evidence that only portions of the electrode held down with fastening straps may conduct current [1].

Methods

Imaging was performed using the 3T Achieva Phillips MRI system at the McKnight Brain Institute, University of Florida. A pair of MR safe TACS electrodes (25 cm2 carbon rubber inside saline soaked sponges) was placed in a trans-temporal (T7-T8) configuration over the subject’s scalp and hair. Structural (T1) data were acquired as 160 sagittal slices, with 1 mm thickness per slice. Each slice had 256 x 256 pixels and a resolution of 0.9375x0.9375mm. T1 data was used to construct computational models of the subject. Subsequently, 1.5 mA amplitude current pulses at approximately 10 Hz were passed between the electrodes during Magnetic Resonance Electrical Impedance Tomography (MREIT) acquisitions. The MREIT dataset was acquired using a modified Philips mFFE-ssfp sequence (224 x 224, pixel size 1mm2) in a total of three axial slices (5 mm thickness). Phase images from MREIT acquisitions were then processed and rescaled to Bz [2]. The T1 data were segmented into ten tissue types and two electrodes for simulations and assigned literature referenced tissue conductivities3. Simulated current densities were converted into Bz using an FFT implementation of the Biot-Savart law, and resampled to the experimental resolution for comparisons. Electrode areas, bone and skin conductivity values in the models were then modified separately, and in combination to determine if large differences between synthetic and experimental Bz measurements may be caused by overestimation of electrode area, underestimation of bone conductivity or overestimation of skin conductivity or all three. Electrode contact areas were simulated to be their full original apparent size or reduced to a quarter of their original size. Bone tissue conductivity was chosen to be 0.0109 S/m or 0.1 S/m. Skin tissue conductivity was simulated to be either 0.43 S/m or 0.01 S/m as reported by Parazzini et al. [4]

Results

Figure 1 shows brain-masked MREIT Bz and predicted Bz distributions in the MREIT central axial slice for a T7-T8 electrode montage, and line profile plots. All Bz patterns were in agreement with the Biot-Savart law for current direction from left (T7) to right (T8) with high Bz values observed in the anterior and low Bz values observed in the posterior. Reducing electrode contact area by four as shown in Figure 1C did not increase synthetic Bz values, but caused slight change in Bz profiles nearby the electrodes. Increasing bone and decreasing skin conductivity values individually, as shown in Figures 1D and 1E, respectively, did not produce clear changes in Bz patterns, but increased Bz range, as shown in Figure 1G. The closest agreement between MREIT and synthetic Bz were observed when incorporating changes in electrode model, bone and skin conductivity into simulations altogether as shown in Figure 1F-G.

Discussion

Smaller electrode contact areas should cause larger current densities to enter the head, but this modification alone did not increase Bz range inside the brain. This may indicate that the current density was mainly shunted by the high conductivity skin tissue. Decreasing skin and increasing bone conductivity values individually altered the Bz range as expected, allowing more current to enter the brain. The combination of decreased electrode contact area, high bone and low skin conductivity parameters produced synthetic Bz maps that best matched the scale of experimental Bz distributions. While not conclusive, these findings indicated that our previous computational parameters used in Bz simulations of TACS may be inaccurate. Thus, our study could contribute to improvements of TACS finite element modeling. Our method can also be used to reconstruct conductivity distributions, though reconstructions in skin and bone are likely to be poor. Further processing of these data may yield useful findings of interest to workers involved in validation of computational models of transcranial electrical stimulation or in inverse source imaging applications.

Acknowledgements

Research reported in this publication was supported by the National Institute Of Neurological Disorders And Stroke of the National Institutes of Health under Award Number R21NS081646 to RJS. In addition, a portion of this work was performed in the Advanced MRI/S (AMRIS) Facility at the McKnight Brain Institute of the University of Florida, which is part of the National High Magnetic Field Laboratory (supported by National Science Foundation Cooperative Agreement DMR-1157490, the State of Florida, and the U.S. Department of Energy).

References

[1] Woods, A.J., et al., "A technical guide to tDCS, and related non-invasive brain stimulation tools." Clinical Neurophysiology vol. 127, pp. 1031–1048, 2016.

[2] E. J. Woo and J. K. Seo, "Magnetic resonance electrical impedance tomography (MREIT) for high-resolution conductivity imaging," Physiological measurement, vol. 29, no. 10, pp. 1-26, 2008.

[3] A. Indahlastari, et al., “Changing head model extents affects finite element predictions of transcranial direct current stimulation distributions.” Journal of Neural Engineering vol 13 no. 6, pp. 066006, 2016.

[4] M. Parazzini, et al., "Transcranial direct current stimulation: Estimation of the electric field and of the current density in an anatomical human head model," IEEE Trans Biomed Eng, vol. 58, pp. 1773-1780, 2011.

Figures

Figure 1. Magnetic flux density (Bz) distributions for trans-temporal current injection in experimental and simulation results. Brain-masked Bz distributions are shown for results obtained from (A) MREIT and simulations with (B) high skin, low bone conductivity with full electrode, (C) high skin, low bone conductivity and quarter electrode area (D) high bone and high skin conductivity with full electrode (E) low skin and low bone conductivity with full electrode (F) low skin, high bone and quarter electrode area. Electrode locations are indicated in panel A. Bz profiles values along the profile lines overlaid on panels A-F are shown in G.

Proc. Intl. Soc. Mag. Reson. Med. 25 (2017)
2417