Original Research ARTICLE
Gender differences in current received during transcranial electrical stimulation
- 1Aaken Laboratories, Davis, CA, USA
- 2Sutter Center for Psychiatry, Sacramento, CA, USA
- 3Center for Neuroscience, University of California Davis, Davis, CA, USA
Low current transcranial electrical stimulation (tCS) is an effective but somewhat inconsistent tool for augmenting neuromodulation. In this study, we used 3D MRI guided electrical transcranial stimulation modeling to estimate the range of current intensities received at cortical brain tissues. Combined T1, T2, and proton density MRIs from 24 adult subjects (12 male and 12 female) were modeled with virtual electrodes placed at F3, F4, C3, and C4. Two sizes of electrodes 20 mm round and 50 mm × 45 mm were examined at 0.5, 1, and 2 mA input currents. The intensity of current received was sampled in a 1-cm sphere placed at the cortex directly under each scalp electrode. There was a 10-fold difference in the amount of current received by individuals. A large gender difference was observed with female subjects receiving significantly less current at targeted parietal cortex than male subjects when stimulated at identical current levels (P < 0.05). Larger electrodes delivered somewhat larger amounts of current than the smaller ones (P < 0.01). Electrodes in the frontal regions delivered less current than those in the parietal region (P < 0.05). There were large individual differences in current levels that the subjects received. Analysis of the cranial bone showed that the gender difference and the frontal parietal differences are due to differences in cranial bone. Males have more cancelous parietal bone and females more dense parietal bone (P < 0.01). These differences should be considered when planning tCS studies and call into question earlier reports of gender differences due to hormonal influences.
Low intensity transcranial electrical stimulation (tCS) is a non-invasive procedure for modifying neural networks within the brain. There is an extensive literature reporting effective treatment for a number of conditions (1–7) (and included references). However, the findings are often inconsistent with a large degree of variability between individual subjects. Earlier reports from our laboratory (8) and others (9–11) have suggested that much of the inconsistency can be accounted for by individual differences in head anatomy and differences in the amount of current received by the brain. Due to the need for more precise targeting of stimulation recent efforts have turned to modeling of current levels and pathways. Modeling current within the head has progressed steadily. Early investigators modeled current using a single-slice of MRI within the head (12), a single subject (9, 13–18), or simplified anatomic models of cranial components (13, 19, 20). Studies by Ruffini et al. (18), Miranda et al. (21), and Sadleir et al. (9) have addressed some of these issues by using MRI modeling from representative heads. Our laboratory has attempted to take modeling somewhat further by using multiple subjects and verifying the model with physical current intensity measurements to establish validity (8). We have demonstrated that it is possible to determine tissue resistance and cranial conductivity from a combination of MRI measurements (8). This MRI based technique assesses the amount of hydrogen within a tissue as an index of water and makes the assumption that the conductivity of a tissue is largely determined by the amount of salt water within the tissue. It then calculates the MRI values of hydrogen within voxels to estimate resistivity. The present report uses our previously validated MRI modeling (8) to estimate the current received at commonly used transcranial stimulation locations (22–25). We chose to constrain the electrode locations, electrode types, and stimulus intensities to produce a representative but manageable data set. Four locations were selected, one from each side of the head above the left and right parietal bones, and each side of the frontal bone. We estimated the current received at the cortex with two commonly used sizes of electrodes at three current levels. The brain region targeted by investigators is generally envisioned to be directly under the electrode site, so we modeled the current levels at the cortex directly under the electrodes.
Materials and Methods
The subjects were 24 adults (12 male and 12 female) without known anatomical anomalies. Males had a mean age of 53 years ± 11.5 (SD) and an age range of 34–68 years. The females had a mean age of 50.5 years ± 14.3 (SD) and an age range of 21–75 years. The human subject protocol was approved by the Institutional Review Board of the Sutter Institute for Medical Research. Written informed consent was obtained from all participants.
MRI scans were performed on each subject in order to model current delivered to cortical targets. Virtual electrodes were placed at four scalp location C3, C4, F3, and F4 as defined by the 10–10 system (Figure 1) and the 10–20 system (26). The modeling estimated current density delivered to a region of interest (ROI) consisting of a sphere of cortical tissue directly under the virtual electrodes (Figure 2). To encompass the range of current stimulation intensities typically used in tCS electrical current levels of 0.5, 1, and 2 mA were modeled. The current traveled in a single direction for each electrode pair C3-anode to C4-cathode and F3-anode F4-cathode. Two sizes of commonly used electrodes were modeled as 20 mm round and 50 mm × 45 mm (Figure 1). A 1 cm sphere of cortical tissue within the cranium and directly under virtual electrodes (Figure 2) was placed at four scalp locations C3, C4, F3, and F4 as defined by the 10–10 system for scalp placement.
Figure 1. Above are illustrations of the large and small electrodes placed on the head by the 10–10 software. The large electrode is 50 mm × 45 mm the small electrode is 20 mm diameter.
Figure 2. The above left image is a T2 MRI slice with the areas of interested outlined in green directly under each parietal electrode. The image on the right is the same slice after the simulation.
To map the conductivity of brain tissue, a three-dimensional (3D) measurement of the hydrogen distribution in head and brain is needed rather than the typical individual slice MRI record. T1, T2, and proton density (PD) imaging each capture different aspects of the water (hydrogen) present in tissues, but each alone lacks the power to provide a precise brain conductivity model. Consequently for each subject, the data from the three types of MRI were combined to model conductivity. Diffusion weighted recordings were not used in this instance because they are considered to be unreliable at the cortex. The recordings were obtained on a three Tesla General Electric MRI (Discovery M750) machine with a 1 mm × 1 mm × 1 mm slice spacing. This combination of MRIs constitutes the Aaken Insite Protocol (27). Image data were combined into a single 3D representation and a conversion equation was applied to achieve an index of resistivity, yielding the subject’s resistivity model. The resistivity to the MRI intensity is expressed by the formula:
Where v ε[0,1] is the normalized intensity of the combined image at the given voxel; R is the resistivity; K = 16,000; E = 4; and D = 65 are the adjustable parameters.
Finite Element Analysis
Finite element and boundary conditions are given in detail at Russell et al. (8). Briefly, the subject’s resistivity distribution was translated to a rectangular prismatic linear finite element model. The model matrix equation and boundary conditions were formulated from the Galerkin equations (28). Solutions to the system matrix equations were obtained by using the conjugate gradient method (29). The finite element models solved the Laplace equation and current densities within the models were determined from the finite element model solution by calculating current densities for each voxel.
Modeling Current Density
Current density maps were achieved by placing virtual electrodes [red (+) and blue (−)] according to the 10–10 system for electrode placement at four locations. Two were placed above the frontal pole F3–F4 (above frontal bone) and two above the region of the motor cortex at C3–C4 (above parietal bone). Two sizes of virtual electrodes were modeled as either 20 mm diameter round or a 50 mm × 45 mm.
The area sampled for current density was a 10 mm virtual sphere under the cranium (green circle in Figure 2) directly under each electrode. The sphere defined the voxels in the cortical region assessed for the analysis.
Data analysis was performed using IBM SPSS software (Version 22, Chicago, IL, USA), which adheres to a general linear model. Alpha level for Type I error was set at 0.05 for rejecting null hypotheses. All data were expressed as mean ± standard error of the mean (SEM).
Differences in current density between stimulation intensity, independent of gender, location, and electrode size, was analyzed using one-way analysis of variance (ANOVA). Differences in current density between electrode size and electrode location were analyzed using repeated measures ANOVA. Differences in cancelous bone thickness were analyzed with independent t-test.
Table 1 presents the minimums (min), maximums (max), and means of the modeled current densities across all conditions (gender electrode placement, electrode size, and stimulation intensity in microampere per square centimeter.
Table 1. The table below shows the minimum, maximum and mean values by gender and stimulus intensity for each size of electrode 50 mm by 45 square and 20 mm round.
The range of minimum and maximum currents predicted at the ROIs was between 0.5 and 44.8 μA/cm2. There was as much as a 10-fold difference among voxels within the spherical ROIs that reflects the complex resistivity within the cortical structures. The range of mean current density values was between 1.7 μA/cm2 for the females with 0.5 mA stimulation current at C3 to a high of 17.9 μA/cm2 for the males with a 2 mA stimulation current at C3.
The relationship between stimulation current and current density in the ROIs was analyzed in a stepwise manner in order to determine contributions of stimulation intensity, electrode size, and gender to the modeled current density in brain.
First, modeled current density was averaged over the four ROIs per subject for both electrode sizes and for combined males and females. Thus, independent of gender, electrode location, and electrode size, modeled current density increased correspondingly to the stimulation intensity (Figure 3). Thus, with each doubling of stimulation intensity, our model predicted significant doubling of current density in the ROIs [F(2,69) = 52.3, P < 0.001].
Figure 3. Above is a graph of the increase in mean current densities for all subjects, electrodes and locations at 0.5, 1, and 2 mA input currents. Error bars represent ±standard error of the mean.
Modeled current density was next analyzed for differences in electrode size by averaging values over the four ROIs per subject with male and female subjects combined. The relationship between the amount of input stimulation current and the mean current density was nearly linear for both the 20 and 50 mm electrodes (Figure 4). The 50 mm electrodes produced significantly higher current densities than the 20 mm electrodes [F(1,46) = 7.138, P = 0.010] across all stimulation intensities.
Figure 4. Above are the mean current densities for the 50 mm × 45 mm electrodes and the 20 mm round electrodes at the three current intensities 0.5, 1, and 2 mA. Error bars represent ±standard error of the mean.
Modeled current density was next analyzed for differences in gender in parietal (Figure 5A) and frontal stimulation (Figure 5B) sites averaged across electrode size and hemisphere. Males received significantly higher modeled current density than females after parietal stimulation [F(1,22) = 4.853, P = 0.04]. There was no difference between males and females with frontal stimulation [F(1,22) = 0.333, P = 0.57].
Figure 5. Above are graphs of the mean current densities for the males (blue) and females (red) at the parietal locations [(A) top] and the frontal locations [(B) bottom]. Error bars represent ±standard error of the mean.
We measured the distance between the electrodes and the cortical surface, the thickness of the bone under the electrode, and the thickness of the cancelous bone using the MRI scans. No gender differences were found in the distance between the scalp surface and the cortex at the parietal bone (female mean = 15.03 mm; male mean = 15.08 mm) nor in the frontal bone thickness (female mean = 10.5; male mean = 10.1). No gender differences were found in the total parietal bone thickness (female mean = 7.34 mm; male mean = 7.62 mm). A highly significant difference was found when cancelous parietal bone thickness was compared by gender (Figure 7). The women had denser cortical bone with very little cancelous bone (Figures 6 and 7). Males had significantly thicker cancelous bone in the parietal region [t(22) = 4.245, P < 0.001] and in the frontal region [t(22) = 2.875, P = 0.009].
Figure 6. Above are two parietal coronal MRI sections through male (right) and female (left) heads. The arrows point to the cancelous bone region.
Figure 7. The graphs show the mean measurements in millimeters of the cancelous bone in the parietal and frontal regions by gender (blue male red female). The upper graph is for cancelous bone thickness and the lower graph is total bone thickness. Error bars represent ±standard error of the mean.
Knowledge of the levels of current received at a target during electrical stimulation is a basic element for the objective evaluation of outcomes. The skull is a high resistance and anatomically complex organ with multiple bone thicknesses and densities that impact current flow. This study suggests that when planning tCS treatments consideration of the regional differences in bone and tissue resistivity should be considered. The cortical targets directly under the stimulating electrodes are receiving a wide range of current in the normal subject population. The large disparity in current densities that we have observed is likely to account for a significant portion of the variation seen in reported patient outcomes. Because of its high resistivity, the bone under the scalp is the most important factor for determining the amount of current received at the targeted cortex. The gender differences in this study were unanticipated. The large differences in parietal current density values were particularly significant. The male subjects received approximately 45% greater current to the cortex than female subjects. Given the large number of investigators using tCS, a more a detailed normative study of bone density, gender, and development should be undertaken. Such a study should examine a larger number of electrode locations and skull regions to provide a detailed map of skull resistivity. There may be other factors that influence skull density and resistivity such as chronic alcohol or drug abuse, smoking, handedness, etc.
Studies that have previously attributed sex differences in tCS results to differences in neural processing, neural plasticity, or hormonal conditions (30, 31) may need to be reconsidered as those differences may have been due to differences in bone density. Studies that have used subjects of both genders may find results more meaningful if the results are analyzed by gender.
The data in this study raise some additional questions. An individual three-dimensional MRI or other analysis of cranial bone is clearly desirable to plan for optimal stimulus intensity, but usually such data is unavailable, and when available it may add significant cost to tCS studies. Should adjustments in stimulation levels be made to compensate for gender and electrode locations on the head? If such a strategy is adopted, exposing women to higher current densities may pose problems with scalp sensations and controlling for placebo effects. Should men receive lower dosages of current than they have been receiving past? In studies where women had a successful outcome, but males did not were the male current levels too high? Recently, there has been some recognition of the need to model stimulation in children (32, 33). These studies may provide us an idea of when gender differences appear.
Limitations of the Study
We did not analyze the effect of distance between electrodes, but decreasing the inter electrode distance reduces the resistance and increases the shunting through the skin. Smaller heads move electrodes closer together and lead to increased shunting. This can be compensated for to some degree by adjusting the stimulation intensity to compensate for shunting through the skin.
The MRI based modeling used in this study was previously validated with 25 surface measurements from the heads of each of three individuals (8), but validation of the results using intracranial measurements have not yet been completed. It is possible that absolute values may need to be adjusted when internal tissue validation is available, but relative values are not likely to change.
The current density values modeled for each spherical 1 cm diameter ROI was the mean of 524 mm3 voxels recorded by the MRI. There was considerable variation of current densities between the voxels within that sphere and thus the predicted stimulation by the modeling of individual voxel-sized micro environments are somewhat heterogeneous. Similarly, the resistance modeled for each voxel is a composite value formed from the multiple nervous tissues (neurons, glia, interstitial space, etc.) within each voxel. Thus, the modeling is limited to the precision of the MRI as improvements are made in the quality of MRI recording the precision of the modeling will also improve.
This report is focused on current density, but voltage is a more appropriate measure when the concern is a stimulus locked neuronal discharge. The data in this study are best applied to direct current and slow wave alternating currents and less applicable to transcranial square wave pulse stimulation or studies of stimulation rates above 150 Hz because impedance rather than simple tissue resistance becomes an issue. The optimal current dosage therapeutic window needed to elicit a response likely varies for AC, DC, and random stimulation and depend on the targets selected, but knowing the range of currents that subjects are receiving should help in beginning to establishing optimal levels.
These simulations may be useful in helping investigators evaluate the results published in the literature and in guiding their selection of stimulation levels for future studies.
This study only sampled four locations using the 10–10 system. Many investigators use the 10–20 system. The locations we have chosen C3, C4, F3, and F4 that are identical in the two systems (26). A replication in a larger sample size is warranted to cover the additional locations are specific to the two systems.
Given the linear nature of the relationship of resistivity to current density, it is reasonable to assume that input currents that differ from what we have modeled can be inferred. However, caution should be applied when stimuli are substantially above or below the sampled range.
There are likely larger differences between individuals than what is represented here. Our 24 subjects do not represent the broad range of cranial differences that is seen in the normal population. Although race was not assessed, all of the subjects in this population appeared to be of European decent. The study does not represent the range of body types seen in the human population. For a recent analysis of the diversity, see Ref. (34–36). Clinical factors like bone density, age, cranial anomalies, tumors, inflammation, etc. are likely to further influence current received.
Current density dosing levels that are too high or too low may produce poor tCS outcomes for patients. The range of current intensities received at the cortex with tCS varies with stimulus intensity, electrode size, scalp location, and subject gender. There can be a significant gender difference in the current received when the brain is stimulated transcranially. In general, men receive more current than women when stimulated above parietal bone. Larger electrodes deliver more current than smaller electrodes. Study designs should consider the influence of electrode size, subject gender, scalp location, and other issues of skull density that may influence stimulus intensity. These variables should be factored in when designing or interpreting study outcomes, and may have important consequences in clinical practice, and for patients who receive tCS treatments for neurological or psychiatric disorders.
Conflict of Interest Statement
Dr. Michael Russell has two patents that are related to this study. The other authors declare that the research was conducted in absence of any additional commercial or financial conflict of interest.
This research was supported by the Sutter Institute for Medical Research and the Aaken Research Institute. We thank Radiological Associates of Sacramento and David Seidenwurm for help with the MRI data collection. We thank Ken Van for his help in consulting and making the graphs and figures used in this paper. We also thank Joseph Visse for his help in editing and reviewing the manuscript.
1. Utz KS, Dimova V, Oppenlander K, Kerkhoff G. Electrified minds: transcranial direct current stimulation (tDCS) and galvanic vestibular stimulation (GVS) as methods of non-invasive brain stimulation in neuropsychology – a review of current data and future implications. Neuropsychologia (2010) 48:2789–810. doi:10.1016/j.neuropsychologia.2010.06.002
3. Brunoni AR, Fregni F, Pagano RL. Translational research in transcranial direct current stimulation (tDCS): a systematic review of studies in animals. Rev Neurosci (2011) 22:471–81. doi:10.1515/RNS.2011.042
5. Heimrath K, Sandmann P, Becke A, Muller NG, Zaehle T. Behavioral and electrophysiological effects of transcranial direct current stimulation of the parietal cortex in a visuo-spatial working memory task. Front Psychiatry (2012) 3:56. doi:10.3389/fpsyt.2012.00056
6. Medeiros LF, de Souza IC, Vidor LP, de Souza A, Deitos A, Volz MS, et al. Neurobiological effects of transcranial direct current stimulation: a review. Front Psychiatry (2012) 3:110. doi:10.3389/fpsyt.2012.00110
7. Brunoni AR, Shiozawa P, Truong D, Javitt DC, Elkis H, Fregni F, et al. Understanding tDCS effects in schizophrenia: a systematic review of clinical data and an integrated computation modeling analysis. Expert Rev Med Devices (2014) 11(4):383–94. doi:10.1586/17434440.2014.911082
8. Russell MJ, Goodman T, Pierson R, Shepherd S, Wang Q, Groshong B, et al. Individual differences in transcranial electrical stimulation current density. J Biomed Res (2013) 27:495–508. doi:10.7555/JBR.27.20130074
10. Horvath JC, Carter O, Forte JD. Transcranial direct current stimulation: five important issues we aren’t discussing (but probably should be). Front Syst Neurosci (2014) 8:2. doi:10.3389/fnsys.2014.00002
11. Kim JH, Kim DW, Chang WH, Kim YH, Kim K, Im CH. Inconsistent outcomes of transcranial direct current stimulation may originate from anatomical differences among individuals: electric field simulation using individual MRI data. Neurosci Lett (2014) 564:6–10. doi:10.1016/j.neulet.2014.01.054
14. Lee WH, Deng ZD, Kim TS, Laine AF, Lisanby SH, Peterchev AV. Regional electric field induced by electroconvulsive therapy: a finite element simulation study. Conf Proc IEEE Eng Med Biol Soc (2010) 2010:2045–8. doi:10.1109/IEMBS.2010.5626553
15. Lee WH, Deng ZD, Laine AF, Lisanby SH, Peterchev AV. Influence of white matter conductivity anisotropy on electric field strength induced by electroconvulsive therapy. Conf Proc IEEE Eng Med Biol Soc (2011) 2011:5473–6. doi:10.1109/IEMBS.2011.6091396
17. Parazzini M, Rossi E, Ferrucci R, Liorni I, Priori A, Ravazzani P. Modelling the electric field and the current density generated by cerebellar transcranial DC stimulation in humans. Clin Neurophysiol (2014) 125:577–84. doi:10.1016/j.clinph.2013.09.039
18. Ruffini G, Fox MD, Ripolles O, Miranda PC, Pascual-Leone A. Optimization of multifocal transcranial current stimulation for weighted cortical pattern targeting from realistic modeling of electric fields. Neuroimage (2014) 89:216–25. doi:10.1016/j.neuroimage.2013.12.002
19. Datta A, Elwassif M, Battaglia F, Bikson M. Transcranial current stimulation focality using disc and ring electrode configurations: FEM analysis. J Neural Eng (2008) 5:163–74. doi:10.1088/1741-2560/5/2/007
20. Im CH, Jung HH, Choi JD, Lee SY, Jung KY. Determination of optimal electrode positions for transcranial direct current stimulation (tDCS). Phys Med Biol (2008) 53:N219–25. doi:10.1088/0031-9155/53/11/N03
22. Cerruti C, Schlaug G. Anodal transcranial direct current stimulation of the prefrontal cortex enhances complex verbal associative thought. J Cogn Neurosci (2009) 21:1980–7. doi:10.1162/jocn.2008.21143
23. Bradnam LV, Stinear CM, Byblow WD. Cathodal transcranial direct current stimulation suppresses ipsilateral projections to presumed propriospinal neurons of the proximal upper limb. J Neurophysiol (2011) 105:2582–9. doi:10.1152/jn.01084.2010
24. Teo F, Hoy KE, Daskalakis ZJ, Fitzgerald PB. Investigating the role of current strength in tDCS modulation of working memory performance in healthy controls. Front Psychiatry (2011) 2:45. doi:10.3389/fpsyt.2011.00045
25. Park SH, Koh EJ, Choi HY, Ko MH. A double-blind, sham-controlled, pilot study to assess the effects of the concomitant use of transcranial direct current stimulation with the computer assisted cognitive rehabilitation to the prefrontal cortex on cognitive functions in patients with stroke. J Korean Neurosurg Soc (2013) 54:484–8. doi:10.3340/jkns.2013.54.6.484
26. Jurcak V, Tsuzuki D, Dan I. 10/20, 10/10, and 10/5 systems revisited: their validity as relative head-surface-based positioning systems. Neuroimage (2007) 34:1600–11. doi:10.1016/j.neuroimage.2006.09.024
30. Chaieb L, Antal A, Paulus W. Gender-specific modulation of short-term neuroplasticity in the visual cortex induced by transcranial direct current stimulation. Vis Neurosci (2008) 25:77–81. doi:10.1017/S0952523808080097
31. Lapenta OM, Fregni F, Oberman LM, Boggio PS. Bilateral temporal cortex transcranial direct current stimulation worsens male performance in a multisensory integration task. Neurosci Lett (2012) 527:105–9. doi:10.1016/j.neulet.2012.08.076
32. Minhas P, Bikson M, Woods AJ, Rosen AR, Kessler SK. Transcranial direct current stimulation in pediatric brain: a computational modeling study. Conf Proc IEEE Eng Med Biol Soc (2012) 2012:859–62. doi:10.1109/EMBC.2012.6346067
33. Kessler SK, Minhas P, Woods AJ, Rosen A, Gorman C, Bikson M. Dosage considerations for transcranial direct current stimulation in children: a computational modeling study. PLoS One (2013) 8:e76112. doi:10.1371/journal.pone.0076112
Keywords: transcranial electrical stimulation, MRI, modeling, DC Stimulation, AC Stimulation
Citation: Russell M, Goodman T, Wang Q, Groshong B and Lyeth BG (2014) Gender differences in current received during transcranial electrical stimulation. Front. Psychiatry 5:104. doi: 10.3389/fpsyt.2014.00104
Received: 26 May 2014; Accepted: 02 August 2014;
Published online: 15 August 2014.
Edited by:Nora S. Vyas, National Institutes of Health, USA
Reviewed by:Uner Tan, Cukurova University, Turkey
René Martin Müri, University of Bern, Switzerland
Jerome Brunelin, Centre Hospitalier le Vinatier, France
Copyright: © 2014 Russell, Goodman, Wang, Groshong and Lyeth. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael Russell, Aaken Laboratories, 216 F Street #76, Davis, CA, USA e-mail: firstname.lastname@example.org