- 1Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 2Department of Biomedical Engineering, The City College of New York, City University of New York, New York, NY, United States
We use auditory verbal hallucinations (AVH) to illustrate the challenges in defining and assessing target engagement in the context of transcranial direct current stimulation (tDCS) for psychiatric disorders. We defined the target network as the cluster of regions of interest (ROIs) that are consistently implicated in AVH based on the conjunction of multimodal meta-analytic neuroimaging data. These were prescribed in the New York Head (a population derived model) and head models of four single individuals. We appraised two potential measures of target engagement, tDCS-induced peak electric field strength and tDCS-modulated volume defined as the percentage of the volume of the AVH network exposed to electric field magnitude stronger than the postulated threshold for neuronal excitability. We examined a left unilateral (LUL) montage targeting the prefrontal cortex (PFC) and temporoparietal junction (TPJ), a bilateral (BL) prefrontal montage, and a 2 × 1 montage targeting the left PFC and the TPJ bilaterally. Using computational modeling, we estimated the peak electric field strength and modulated volume induced by each montage for current amplitudes ranging 1–4 mA. We found that the LUL montage was inferior to both other montages in terms of peak electric field strength in right-sided AVH-ROIs. The BL montage was inferior to both other montages in terms of modulated volume of the left-sided AVH-ROIs. As the modulated volume is non-linear, its variability between montages reduced for current amplitudes above 3 mA. These findings illustrate how computational target engagement for tDCS can be tailored to specific networks and provide a principled approach for future study design.
Introduction
Transcranial direct current stimulation (tDCS) is a non-invasive neuromodulation technique currently under evaluation for the treatment of neuropsychiatric conditions. tDCS involves the application of a weak electric current that flows through the brain from anodal to cathodal scalp electrodes. Proposed mechanisms of action for tDCS involve polarity-dependent changes in the resting membrane potential (depolarization at the anode and hyperpolarization at the cathode) and changes in synaptic plasticity (1–5).
As the application of tDCS to neuropsychiatric disorders is expanding, the field has begun to recognize significant conceptual gaps that hamper the design and evaluation of tDCS treatment protocols. In particular, there are currently no assays for assessing target engagement in tDCS studies, which makes it very difficult to discern the basis for therapeutic efficacy or lack thereof. This contrasts with drug development, where measures of target engagement are an essential and integral part of study design. In human studies, the in vivo assessment of drugs is critically dependent on proximal markers (e.g., receptor occupancy) as these enable a direct correlation between target engagement and measurements of drug efficacy or toxicity (6). In tDCS, target engagement depends on patient-specific factors related to head anatomy and functional state and on operator-controlled factors related to stimulation parameters (montage, current amplitude) and administration protocol (frequency and duration of treatment sessions) (7). In psychiatry, an additional challenge involves uncertainties in the definition of the target brain regions as a direct and consistent correspondence between symptoms and brain networks has yet to be established (8).
Here, we present an approach for the assessment of target engagement measures in tDCS studies of psychiatric symptoms. We focus on auditory verbal hallucinations (AVH) in patients with schizophrenia because the neural correlates of AVH are arguably better defined than those of other psychiatric symptoms. Current models of AVH propose a dual pathology involving reduced cognitive control due to hypofunction of prefrontal cortex (PFC) and abnormal activation of speech-related regions, in the superior temporal gyrus and temporoparietal junction (TPJ) (9, 10). This model has informed the tDCS montages that have been used to treat AVH that conventionally target these two regions. Despite initial positive results in case series and open label trials (Table S1 in Supplementary Material), randomized, double-blind, sham-controlled clinical trials (RCTs) have not been consistently differentiated between the active tDCS and the sham condition (11–15). Without methods to assess tDCS target engagement, it is difficult to interpret these results and improve study design.
In response, we outline an approach for improving the definition of the target network for AVH and for the assessment of two measures that could putatively provide a quantitative assessment of target engagement in tDCS. First, we identified the brain regions that comprise the AVH network using data from meta-analyses of the relevant structural and functional imaging studies. Second, to accommodate concerns about anatomical variation, we prescribed the brain regions in the AVH network as regions of interest (ROIs) in a population-based standardized volume conductor head model and in four-head models from single individuals. We then used computational modeling to quantify two putative measures of target engagement, peak electric field strength, and percentage-modulated volume derived from three different tDCS protocols. Peak electric field strength in AVH-network regions was chosen as it tracks cortical excitability in tDCS studies of motor regions (16, 17). We introduced modulated volume as another potential index of target engagement. As AVH involve a network of brain regions (as opposed to a single region like the motor cortex), the efficacy of tDCS may depend on the percentage of the volume of the network being modulated. Our primary aim was to quantify the degree of engagement of the AVH network by different tDCS montages based on the peak electric field strength and percentage-modulated brain volume. Our secondary aim was to evaluate the degree of variation in these as a function of the anatomical variability of the head models used. For modulated volume, which is a non-linear measure, we also examined interindividual variability for each montage as a function of increasing current amplitude within the tolerable current range of 1–4 mA (18). Taken together, these steps outline a novel approach to the design of tDCS interventions for target identification and engagement that is image-guided, multi-target, and supports rational testing of therapeutic hypotheses.
Materials and Methods
Identification of the AVH Network in the Brain
We interrogated databases available through the National Center for Biotechnology using relevant expanded subject headings and free text searches to identify meta-analyses of structural and functional neuroimaging studies of AVH published by February 1, 2017. Coordinates of AVH correlates reported in Talairach space were transformed to Montreal Neurological Institute space, using the “tal2icbm_fsl” transform.1 All coordinates were mapped to the Automated Anatomical Labeling digital atlas (19) to identify anatomical regions of cross-modal convergence. The regions thus identified comprise the target regions of interest (AVH-ROIs) used in subsequent computational models.
Prescription of AVH-ROIs in Population-Level and Individual Head Models
Based on the meta-analytic evidence described earlier, we prescribed the AVH-ROIs identified in five previously developed three-dimensional (3D) realistic models of the human head (20–22). These comprised the New York Head (S1) and four head models belonging to single individuals (S2–S5). The New York Head (S1) is a publicly available standardized volume conductor head model (0.5 mm3 isotropic resolution) that was constructed using T1-weighted MRI scans of 152 adult human brains acquired at 0.5 mm3 isotropic resolution and segmented into 6 tissue compartments comprising the skin, skull, cerebrospinal fluid (CSF), gray matter, white matter, and air cavities (20). The individual head models belonged to three men (S2, S3, and S4; aged 34–41 years) (21–23) and one woman (S5; aged 34 years) (22). They were generated from structural T1-weighted MRI scan with 1 mm3 isotropic resolution and segmented into the same six tissue compartments as the New York Head. These models are available upon request. For S2, individual tissue probability maps corresponding to gray matter, white matter, and CSF were automatically created using FAST (FMRIB Analysis Group, Oxford, UK) (24). The non-brain regions were semiautomatically segmented into three tissue compartments including skin, skull, and air using an in-house segmentation algorithm (21, 25–27), followed by correction of segmentation errors using tools from ITK-SNAP (28). For S3–S5, automated segmentation was performed using SPM (29), followed by correction of segmentation errors using an in-house MATLAB script (30) and ScanIP tools (Simpleware Ltd., Exeter, UK).
tDCS Electrode Montages
The search space for optimal electrode placement in the treatment of AVH is too vast for systematic empirical evaluation. Taking a pragmatic approach, we focused on three montages: (A) a left unilateral (LUL) montage with anode over the left PFC (F3-FP1) and cathode over the left TPJ (P3), (B) a bilateral (BL) prefrontal 1 × 1 montage with the anode over the left PFC (F3-FP1) and cathode over the right PFC (FP2), and (C) a third montage (2 × 1) with one anode over the left dorsolateral PFC (F3-FP1) and two cathodes, one on the left and the other on the right TPJ (T3-P3 and T4-P4). The first two montages have been most widely used in tDCS studies (Tables S1 and S2 in Supplementary Material), while the third was included because of its theoretical potential to provide better target engagement of right-sided target AVH-ROIs.
Electric Field Computational Modeling
Computational modeling was performed separately for each of the three montages in each of the five head models. In addition, we assessed the effect of varying current amplitude across all three montages and head models within the range of 1–4 mA, which is known to be tolerable in humans (31). We used the same computational modeling approach based on our previous work (21, 26). The tDCS sponge electrodes were modeled with 5 cm × 7 cm surface intersecting the scalp. For each montage, the complete 3D head models incorporating the tDCS electrodes were adaptively tessellated to produce finite element models using the restricted Delaunay triangulation algorithm (32). The electric conductivities used for each tissue type across all five head models were as follows: gray matter = 0.33 S/m, white matter = 0.14 S/m, CSF = 1.79 S/m, skin = 0.43 S/m, skull = 0.0132 S/m, and air = 0 S/m (21, 23, 27, 33, 34). The electrodes were assumed to have the conductivity of saline (1.4 S/m) (22). Constant electric current was applied to the electrode surfaces away from the head (26). The LUL and BL montages correspond to those that have been used in the current literature, and accordingly current intensity was set at 2 mA, unless otherwise stated (Table S1 in Supplementary Material). For the 2 × 1 montage, 2 mA was applied to the frontal anode and −1 mA to each of the posterior cathodes, unless otherwise stated. The electric field distribution was computed by solving the Laplace equation using the preconditioned conjugate solver within ANSYS (ANSYS Inc., Canonsburg, PA, USA) (21, 26). In addition, we modeled radial cortical electric field that represents the inward and outward component of the electric field relative to the cortical surface (22).
Data Analysis
We computed the electric field magnitude (E) and the percentage modulated volume for each AVH-ROI. The modulated brain volume was defined as the percentage of the brain volume (within the left or right AVH networks) exposed to electric field magnitude (E) stronger than the modulation threshold (Eth), i.e., the volume where E/Eth ≥ 1 (21, 34). We acknowledge that the modulation threshold in living humans is yet to be conclusively determined. For these analyses, we used a published estimate of the modulation threshold (Eth = 0.2 V/m) based on previous evidence regarding the minimum electric field strength likely to change the firing rate of neurons in model in vitro systems (35). Several studies have tried to define the minimum applied electric field sufficient to modulate transmembrane potential using different methodologies involving in vitro recordings from single neurons and brain slices and in vivo recordings in animals (36–38); the values reported in the case of direct current stimulation ranged from 0.18 to 0.5 V/m, with 0.2 V/m being most commonly reported. We therefore chose this threshold which theoretically should yield the maximum modulated brain volume for each montage. Subsequently, interindividual variation in peak electric field strength and modulated brain volume across the five head models were determined using coefficient of variation (CV). To examine the effect on the modulated brain volume of altering the current amplitude, the electric field simulations for all montages were scaled linearly to span the range of 1–4 mA and the volume where E/Eth ≥ 1 was estimated for these ranges. Finally, we used analyses of variance, followed by Bonferroni-corrected post hoc pairwise comparisons, to determine if the peak electric field strength and modulated brain volume in the AVH-ROIs differed across the three tDCS montages. The effect size for the peak electric field strength and modulated brain volume between each pair of electrode montages was computed using the Cohen’s d (39).
Results
Identification and Prescription of the Target Brain Network for AVH
We identified six meta-analyses that examined structural (40, 41) and functional differences in patients with AVH compared to healthy individuals (42–45). Details of the primary studies included in each meta-analysis are provided in Tables S3–S10 in Supplementary Material. Across all meta-analyses, brain regions consistently associated with AVH comprised the primary and secondary auditory cortex located in the superior temporal gyrus and in the Heschl’s gyrus; Broca’s and Wernicke’s area, respectively, located in the inferior frontal gyrus (pars opercularis and pars triangularis) and in the posterior superior temporal gyrus; the anterior cingulate gyrus; the somatosensory and motor cortices, respectively, located in the postcentral and precentral gyrus; the insula; and the hippocampus and the thalamus (Figure 1; Table S11 in Supplementary Material). Of note, most of the regions identified were BL (Table S11 in Supplementary Material). Based on these results, AVH-ROIs were prescribed in each of the five head models (Figure 1B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-g001.jpg)
Figure 1. (A) Spatial distribution of the coordinates from meta-analyses of studies of patients with auditory verbal hallucinations (AVH) derived from Table S11 in Supplementary Material. Red = during AVH, green = during auditory or language tasks, and blue = morphometric studies. (B) Three-dimensional rendering of the brain regions of interest associated with AVH in a representative model. L, left.
Effect of Anatomic Variation and Current Amplitude on Peak Electric Field Strength
Figure 2 shows the simulated LUL, BL, and 2 × 1 tDCS configurations. The spatial distributions of the electric field on the cortical surface and directional electric field normal to the cortical surface (inward and outward) are plotted for each modality in each individual. For each montage and each AVH-ROI, we estimated the mean peak electric field strength across the five head models. Figures 3A,B present the peak electric field strengths for each montage averaged across the five head models at the conventional 2 mA current amplitude. Figures 3C,D present the CV of the peak electric field strength in each AVH-ROI across the five head models for each montage. The LUL 1 × 1 montage predominantly induced electric fields in the left hemisphere; the inward current flow occurred mostly through the left dorsolateral prefrontal cortex while the outward current flow occurred through the left temporoparietal and occipital regions. The BL prefrontal 1 × 1 montage generated comparable electric field magnitude in both hemispheres; the inward and outward flow occurred mostly over the prefrontal and frontal areas. The 2 × 1 montage induced electric fields in the left- and in the right-sided AVH-ROIs; the inward and outflow flow occurred over the frontal cortex and the temporoparietal regions, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-g002.jpg)
Figure 2. Electric field modeling of transcranial direct current stimulation (tDCS) for three montages: (A) left unilateral 1 × 1 montage with F3-FP1 as the anode and T3-P3 as the cathode, (B) bilateral prefrontal 1 × 1 montage with F3-FP1 as the anode and FP2 as the cathode, and (C) 2 × 1 montage with F3-FP1 as the anode and cathodes at T3-P3 and T4-P4. First column shows tDCS electrode placements representing anode (red) and cathode (blue) electrodes in the New York head model (Subject 1; S1). Rows from top to bottom show electric field magnitude (E) and electric field normal to the cortical surface (En) for the models of Subjects 1–5 (S1–S5) at current amplitude of 2 mA.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-g003.jpg)
Figure 3. (A) Regional electric field strength and spatial distribution at a current of 2 mA generated by each transcranial direct current stimulation montage in a representative head model. (B) Peak electric field magnitude at a current of 2 mA in each auditory verbal hallucinations (AVH) region of interest (ROI). Bars show mean values, and error bars show SD across the five head models for each target AVH-ROI. (C) Spatial distribution of the coefficient of variation (CV) of peak electric field magnitude in each target AVH-ROI across the five head models for each montage. (D) CV of the peak electric field magnitude in target AVH-ROI across the five head models for each montage. LUL, left unilateral 1 × 1 montage with F3-FP1 as the anode and T3-P3 as the cathode; BL, bilateral prefrontal 1 × 1 montage with F3-FP1 as the anode and FP2 as the cathode; 2 × 1, montage with one anode at F3-FP1 and cathodes at T3-P3 and T4-P4. IFGo, inferior frontal gyrus (pars opercularis); IFGt, inferior frontal gyrus (pars triangularis); ACG, anterior cingulate gyrus; PreCG, precentral gyrus; PoCG, postcentral gyrus; HES, Heschl’s gyrus; STG, superior temporal gyrus; HIP, hippocampus; INS, insula; THA, thalamus.
We found significant differences in peak electric field strength across the three electrode montages (LUL, BL, and 2 × 1) as shown in Table 1. Bonferroni-corrected pairwise comparisons showed the montages induced significantly different electric field strengths in several right-sided AVH-ROIs. Specifically, the BL montage generated higher peak electric field strengths in the right inferior frontal gyrus and insula than the other two montages (p < 0.001). The 2 × 1 montage produced higher peak electric field in the right postcentral gyrus than the others (p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-t001.jpg)
Table 1. Differences in peak electric field strength for each region of interest in the auditory verbal hallucinations network across electrode montages.
Effect of Anatomic Variation and Current Amplitude on Percentage Modulated Brain Volume
Figure 4 shows the electric field maps relative to the threshold electric field (E/Eth) for current amplitudes in the range of 1–4 mA for each montage. For each montage and each AVH-ROI in each head model, we estimated the percentage of modulated volume [i.e., percentage of the AVH network volume in which the electric field strength (E) was above 0.2 V/m]. We used the threshold 0.2 V/m as it is the postulated threshold for neuronal excitability (35). Figures 5A,B show the percentage modulated volume (E ≥ Eth) for each of the three montages for the left- and right-sided AVH-ROIs for current amplitudes ranging from 1 to 4 mA averaged across the five head models. Figures 5C,D show the corresponding CV in the left- and right-sided AVH-ROIs for each montage for current amplitudes ranging from 1 to 4 mA. The corresponding data for each AVH-ROI are shown in Figures S1–S4 in Supplementary Material. Increasing the current amplitude increased modulated volumes and reduced interindividual variability, particularly for the LUL and 2 × 1 montages (Figures 4 and 5). Variability between montages reduced in general for current amplitudes above 3 mA. Nevertheless, we found significant differences in modulated volume across the three electrode montages (LUL, BL, and 2 × 1) as shown in Table 2. The modulated brain volume in the left-side AVH-ROIs was significantly different across the montages for the range of 2–4 mA at which lower modulated brain volume was produced by the BL compared to the other two montages (p < 0.05). In contrast, the BL generated statistically higher modulated volume in the right-side AVH-ROIs than the other two montages at low current of 1 mA (p < 0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-g004.jpg)
Figure 4. Brain regions with electric field strength above the modulation threshold (0.2 V/m) in a representative model as a function of current amplitude. Spatial distribution of the electric field shown for the cortical surface and in a representative coronal slice for current amplitudes ranging from 1 to 4 mA for each montage. LUL, left unilateral 1 × 1 montage with F3-FP1 as the anode and T3-P3 as the cathode; BL, bilateral prefrontal 1 × 1 montage with F3-FP1 as the anode and FP2 as the cathode; 2 × 1, montage with one anode at F3-FP1 and cathodes at T3-P3 and T4-P4. R, right.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-g005.jpg)
Figure 5. Percentage modulated volume in (A) left-sided and (B) right-sided auditory verbal hallucinations (AVH) regions of interest (ROIs) across all head models for each of the three montages as a function of current amplitude. Coefficient of variation in percentage-modulated volume in (C) left-sided and (D) right-sided AVH-ROIs across all head models for each of the three montages. LUL, left unilateral 1 × 1 montage with F3-FP1 as the anode and T3-P3 as the cathode; BL, bilateral prefrontal 1 × 1 montage with F3-FP1 as the anode and FP2 as the cathode; 2 × 1, montage with one anode at F3-FP1 and cathodes at T3-P3 and T4-P4.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/342600/fpsyt-09-00048-HTML/image_m/fpsyt-09-00048-t002.jpg)
Table 2. Differences in modulated brain volume for current amplitudes in the range of 1–4 mA for the left- and right-sided regions of interest (ROIs) in the auditory verbal hallucinations network across electrode montages.
Discussion
Here, we use AVH to illustrate the challenges in defining and assessing target engagement in the context of tDCS for psychiatric disorders and to propose a strategy for the evaluation of tDCS protocols.
There are three main challenges. First, AVH arise from dysfunction in a brain network (Figure 1) rather than a single brain region which leads to uncertainty about which regions within the AVH-related network may constitute ideal targets for tDCS (46). Moreover, commonly used montages often target left-sided AVH-ROIs while evidence from the meta-analyses (Figure 1; Table S11 in Supplementary Material) and from individual studies point to additional involvement of right-sided brain regions to the pathophysiology of AVH (47–51). Second, there is inter-individual variability in brain structure and function in general (52) and in the AVH-related network in particular (50, 51). Individual differences in the spatial distribution of tDCS-induced electric fields may impact on the engagement of AVH-related regions and hence clinical response. Third, there is no behavioral, cognitive or neurophysiological measure that can be used to assess AVH-related neural changes during or following a single tDCS session, in contrast to electroconvulsive therapy (53) and repetitive transcranial magnetic stimulation (54) where seizure and motor excitability threshold can be, respectively, used as a threshold for titrating treatment.
In response to these uncertainties, we tested whether computational modeling of electric field measures induced by tDCS can be meaningfully used to assess the engagement of the AVH network. We mined the available neuroimaging literature to define target ROIs implicated in AVH (Figure 1) and used computational modeling approach to assess the engagement of these ROIs in five head models using the three different tDCS montages (Figure 2). The search space for optimal electrode placement in the treatment of AVH is too vast for empirical evaluation of brute force. Taking a pragmatic approach, we considered the LUL and BL prefrontal montages as these are most commonly used in the literature (Tables S1 and S2 in Supplementary Material). We also evaluated an experimental 2 × 1 montage, because of its theoretical potential to provide better target engagement of right-sided target AVH-ROIs. We used two measures of target engagement, specifically peak electric field magnitude in each AVH-ROI and modulated brain volume of the AVH-network.
In terms of peak electrical field strength, the three montages differed only in a few right-sided AVH-ROIs. Specifically, the BL montage generated higher peak electric field strengths in the right inferior frontal gyrus and insula, while the 2 × 1 montage produced higher peak electric field in the right postcentral gyrus. Both of these regions have been associated with increased AVH frequency in patients who experience persistent AVH (44). The strength of the tDCS-induced electric field is a conventional measure of assessing the spatial distribution of the electric fields generated by different electrode configurations. However, its functional significance in terms of clinical efficacy has not been tested and is currently unknown. If we assume that peak electric field strength is associated with clinical efficacy, the current results suggest that the commonly used LUL montage is perhaps the least advantageous if one is interested in electrode montages that have the potential to modulate right-sided AVH-ROIs.
The three montages also differed in terms of modulated volume. In this case, however, differences were noted for the left-sided AVH network where both the LUL and 2 × 1 montages produced higher modulated volumes than the BL for currents ranging from 2 to 4 mA. As expected, the BL montage produced the highest focality within frontal regions due to the proximity of the electrodes (21). Consistent with previous findings (21), increasing the current amplitude increased modulated volume within AVH network. For the threshold-based modulated volume measure, variability between montages reduced in general for current amplitudes above 3 mA. Reduced interindividual variation at higher currents resulted from increased modulated volume within AVH network across multiple realistic head models. This is important for clinical applications, since it implies that increasing the current amplitude may be one possible option to overcoming anatomical variability and uncertainty about optimal electrode placement and potentially improve the efficacy of tDCS for AVH.
This study provides a computational evidence for two different measures of target engagement that behave differently across different montages. There are other parameters in tDCS protocols involving duration of tDCS session, interval of administration, and length of trial (55) that may contribute to efficacy but could not be examined here. We did not address the issue of tolerability for higher current amplitudes although the excellent tolerability of tDCS allows optimism (18, 31, 56). Nevertheless, we have provided a road map for the development of target engagement measures for tDCS that is based on neuroimaging evidence regarding the target networks and accommodates uncertainty about patient-specific abnormalities. This translational approach not only combines advances in computational modeling with knowledge gained from behavioral and functional neuroimaging studies regarding target definition and target engagement but also highlights knowledge gaps and points to avenues for future research.
Generally, the principled computational approach developed here could inform target engagement across tDCS application. Prior computational models typically considered electric field magnitude (or normal direction) in a single ROI. We show that consideration of multiple nodes in a functional network (based on imaging and trial meta-analysis) forces decisions on region stimulation and sparring that then suggest divergent optimal montages. The approach outlined here is particularly relevant for study designs that adopt the Research Domain Criteria (RDoC) framework.2 The RDoC initiative specifies neural circuits that may be theoretically and empirically linked to clinical symptoms and cognitive constructs. This framework therefore can be used to aid in the specification of target networks in future tDCS studies in neuropsychiatry. Moreover, considering a simple linear electric field (magnitude or normal) vs non-linear measures of ROI influence (such as our threshold-based volume measures) leads to different conclusion on the value of current intensity and best montage. As such, tDCS interventions rationalized based on computational modeling should be explicit about the underlying assumptions regarding network (dys)function and side effects (i.e., selection of target and avoid ROIs), biophysics of stimulation (i.e., measure of ROI engagement), and goals (e.g., reduce interindividual variation).
Our approach also addresses issues relevant to personalized psychiatry because, as we have shown, individual variations in electrical field distribution in tDCS can be overcome by protocol modifications such as increasing the current amplitude to 4 mA.
Ethics Statement
This study makes a secondary use of experimental data already published (20–22). All subjects gave written informed consent in accordance with the Declaration of Helsinki.
Author Contributions
WL, NK, and SF conceived and designed the study. WL performed the analysis and computation. NK and SF performed the literature survey and meta-analyses. MB provided three individual head models and helped the data interpretation. All the authors wrote the manuscript, analyzed the data, and reviewed the manuscript.
Conflict of Interest Statement
WL, NK, and SF declare no competing financial interest. The City University of New York has patents on brain stimulation with MB as an inventor. MB has equity in Soterix Medical Inc.
Funding
This work was supported by the National Institutes of Health under grants R01MH104284, CA137788, MH109289, and MH111896.
Supplementary Material
The Supplementary Material for this article can be found online at http://www.frontiersin.org/articles/10.3389/fpsyt.2018.00048/full#supplementary-material.
Footnotes
References
1. Brunoni AR, Nitsche MA, Bolognini N, Bikson M, Wagner T, Merabet L, et al. Clinical research with transcranial direct current stimulation (tDCS): challenges and future directions. Brain Stimul (2012) 5(3):175–95. doi:10.1016/j.brs.2011.03.002
2. Fritsch B, Reis J, Martinowich K, Schambra HM, Ji Y, Cohen LG, et al. Direct current stimulation promotes BDNF-dependent synaptic plasticity: potential implications for motor learning. Neuron (2010) 66(2):198–204. doi:10.1016/j.neuron.2010.03.035
3. Nitsche MA, Fricke K, Henschke U, Schlitterlau A, Liebetanz D, Lang N, et al. Pharmacological modulation of cortical excitability shifts induced by transcranial direct current stimulation in humans. J Physiol (2003) 553(Pt 1):293–301. doi:10.1113/jphysiol.2003.049916
4. Nitsche MA, Jaussi W, Liebetanz D, Lang N, Tergau F, Paulus W. Consolidation of human motor cortical neuroplasticity by d-cycloserine. Neuropsychopharmacology (2004) 29(8):1573–8. doi:10.1038/sj.npp.1300517
5. Nitsche MA, Paulus W. Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation. J Physiol (2000) 527(Pt 3):633–9. doi:10.1111/j.1469-7793.2000.t01-1-00633.x
6. Wagner JA. Strategic approach to fit-for-purpose biomarkers in drug development. Annu Rev Pharmacol Toxicol (2008) 48:631–51. doi:10.1146/annurev.pharmtox.48.113006.094611
7. Woods AJ, Antal A, Bikson M, Boggio PS, Brunoni AR, Celnik P, et al. A technical guide to tDCS, and related non-invasive brain stimulation tools. Clin Neurophysiol (2016) 127(2):1031–48. doi:10.1016/j.clinph.2015.11.012
8. Frangou S. A systems neuroscience perspective of schizophrenia and bipolar disorder. Schizophr Bull (2014) 40(3):523–31. doi:10.1093/schbul/sbu017
9. Hugdahl K. "Hearing voices": auditory hallucinations as failure of top-down control of bottom-up perceptual processes. Scand J Psychol (2009) 50(6):553–60. doi:10.1111/j.1467-9450.2009.00775.x
10. Ford JM. Studying auditory verbal hallucinations using the RDoC framework. Psychophysiology (2016) 53(3):298–304. doi:10.1111/psyp.12457
11. Brunelin J, Mondino M, Gassab L, Haesebaert F, Gaha L, Suaud-Chagny MF, et al. Examining transcranial direct-current stimulation (tDCS) as a treatment for hallucinations in schizophrenia. Am J Psychiatry (2012) 169(7):719–24. doi:10.1176/appi.ajp.2012.11071091
12. Fitzgerald PB, McQueen S, Daskalakis ZJ, Hoy KE. A negative pilot study of daily bimodal transcranial direct current stimulation in schizophrenia. Brain Stimul (2014) 7(6):813–6. doi:10.1016/j.brs.2014.08.002
13. Smith RC, Boules S, Mattiuz S, Youssef M, Tobe RH, Sershen H, et al. Effects of transcranial direct current stimulation (tDCS) on cognition, symptoms, and smoking in schizophrenia: a randomized controlled study. Schizophr Res (2015) 168(1–2):260–6. doi:10.1016/j.schres.2015.06.011
14. Frohlich F, Burrello TN, Mellin JM, Cordle AL, Lustenberger CM, Gilmore JH, et al. Exploratory study of once-daily transcranial direct current stimulation (tDCS) as a treatment for auditory hallucinations in schizophrenia. Eur Psychiatry (2016) 33:54–60. doi:10.1016/j.eurpsy.2015.11.005
15. Mondino M, Jardri R, Suaud-Chagny MF, Saoud M, Poulet E, Brunelin J. Effects of fronto-temporal transcranial direct current stimulation on auditory verbal hallucinations and resting-state functional connectivity of the left temporo-parietal junction in patients with schizophrenia. Schizophr Bull (2016) 42(2):318–26. doi:10.1093/schbul/sbv114
16. Bastani A, Jaberzadeh S. Does anodal transcranial direct current stimulation enhance excitability of the motor cortex and motor function in healthy individuals and subjects with stroke: a systematic review and meta-analysis. Neurophysiol Clin (2012) 123(4):644–57. doi:10.1016/j.clinph.2011.08.029
17. Edwards D, Cortes M, Datta A, Minhas P, Wassermann EM, Bikson M. Physiological and modeling evidence for focal transcranial electrical brain stimulation in humans: a basis for high-definition tDCS. Neuroimage (2013) 74:266–75. doi:10.1016/j.neuroimage.2013.01.042
18. Bikson M, Grossman P, Thomas C, Zannou AL, Jiang J, Adnan T, et al. Safety of transcranial direct current stimulation: evidence based update 2016. Brain Stimul (2016) 9(5):641–61. doi:10.1016/j.brs.2016.06.004
19. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage (2002) 15(1):273–89. doi:10.1006/nimg.2001.0978
20. Huang Y, Parra LC, Haufe S. The New York Head-A precise standardized volume conductor model for EEG source localization and tES targeting. Neuroimage (2016) 140:150–62. doi:10.1016/j.neuroimage.2015.12.019
21. Lee WH, Lisanby SH, Laine AF, Peterchev AV. Comparison of electric field strength and spatial distribution of electroconvulsive therapy and magnetic seizure therapy in a realistic human head model. Eur Psychiatry (2016) 36:55–64. doi:10.1016/j.eurpsy.2016.03.003
22. Datta A, Truong D, Minhas P, Parra LC, Bikson M. Inter-individual variation during transcranial direct current stimulation and normalization of dose using MRI-derived computational models. Front Psychiatry (2012) 3:91. doi:10.3389/fpsyt.2012.00091
23. Lee WH, Kim TS, Cho MH, Ahn YB, Lee SY. Methods and evaluations of MRI content-adaptive finite element mesh generation for bioelectromagnetic problems. Phys Med Biol (2006) 51(23):6173–86. doi:10.1088/0031-9155/51/23/016
24. Zhang YY, Brady M, Smith S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging (2001) 20(1):45–57. doi:10.1109/42.906424
25. Lee WH, Lisanby SH, Laine AF, Peterchev AV. Stimulation strength and focality of electroconvulsive therapy and magnetic seizure therapy in a realistic head model. Conf Proc IEEE Eng Med Biol Soc (2014) 2014:410–3. doi:10.1109/EMBC.2014.6943615
26. Lee WH, Deng ZD, Kim TS, Laine AF, Lisanby SH, Peterchev AV. Regional electric field induced by electroconvulsive therapy in a realistic finite element head model: influence of white matter anisotropic conductivity. Neuroimage (2012) 59(3):2110–23. doi:10.1016/j.neuroimage.2011.10.029
27. Lee WH, Kim TS. Methods for high-resolution anisotropic finite element modeling of the human head: automatic MR white matter anisotropy-adaptive mesh generation. Med Eng Phys (2012) 34(1):85–98. doi:10.1016/j.medengphy.2011.07.002
28. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage (2006) 31(3):1116–28. doi:10.1016/j.neuroimage.2006.01.015
29. Ashburner J. Computational anatomy with the SPM software. Magn Reson Imaging (2009) 27(8):1163–74. doi:10.1016/j.mri.2009.01.006
30. Huang Y, Dmochowski JP, Su Y, Datta A, Rorden C, Parra LC. Automated MRI segmentation for individualized modeling of current flow in the human head. J Neural Eng (2013) 10(6):066004. doi:10.1088/1741-2560/10/6/066004
31. Chhatbar PY, Chen R, Deardorff R, Dellenbach B, Kautz SA, George MS, et al. Safety and tolerability of transcranial direct current stimulation to stroke patients – a phase I current escalation study. Brain Stimul (2017) 10(3):553–9. doi:10.1016/j.brs.2017.02.007
32. Pons JP, Segonne E, Boissonnat JD, Rineau L, Yvinec M, Keriven R. High-quality consistent meshing of multi-label datasets. Inf Process Med Imaging (2007) 20:198–210. doi:10.1007/978-3-540-73273-0_17
33. Lee WH, Lisanby SH, Laine AF, Peterchev AV. Electric field model of transcranial electric stimulation in nonhuman primates: correspondence to individual motor threshold. IEEE Trans Biomed Eng (2015) 62(9):2095–105. doi:10.1109/TBME.2015.2425406
34. Lee WH, Lisanby SH, Laine AF, Peterchev AV. Minimum electric field exposure for seizure induction with electroconvulsive therapy and magnetic seizure therapy. Neuropsychopharmacology (2017) 42(6):1192–200. doi:10.1038/npp.2016.276
35. Frohlich F, McCormick DA. Endogenous electric fields may guide neocortical network activity. Neuron (2010) 67(1):129–43. doi:10.1016/j.neuron.2010.06.005
36. Radman T, Ramos RL, Brumberg JC, Bikson M. Role of cortical cell type and morphology in subthreshold and suprathreshold uniform electric field stimulation in vitro. Brain Stimul (2009) 2(4):e1–3. doi:10.1016/j.brs.2009.03.007
37. Francis JT, Gluckman BJ, Schiff SJ. Sensitivity of neurons to weak electric fields. J Neurosci (2003) 23(19):7255–61.
38. Deans JK, Powell AD, Jefferys JG. Sensitivity of coherent oscillations in rat hippocampus to AC electric fields. J Physiol (2007) 583(Pt 2):555–65. doi:10.1113/jphysiol.2007.137711
39. Nakagawa S, Cuthill IC. Effect size, confidence interval and statistical significance: a practical guide for biologists. Biol Rev Camb Philos Soc (2007) 82(4):591–605. doi:10.1111/j.1469-185X.2007.00027.x
40. Palaniyappan L, Balain V, Radua J, Liddle PF. Structural correlates of auditory hallucinations in schizophrenia: a meta-analysis. Schizophr Res (2012) 137(1–3):169–73. doi:10.1016/j.schres.2012.01.038
41. Modinos G, Costafreda SG, van Tol MJ, McGuire PK, Aleman A, Allen P. Neuroanatomy of auditory verbal hallucinations in schizophrenia: a quantitative meta-analysis of voxel-based morphometry studies. Cortex (2013) 49(4):1046–55. doi:10.1016/j.cortex.2012.01.009
42. Jardri R, Pouchet A, Pins D, Thomas P. Cortical activations during auditory verbal hallucinations in schizophrenia: a coordinate-based meta-analysis. Am J Psychiatry (2011) 168(1):73–81. doi:10.1176/appi.ajp.2010.09101522
43. Kompus K, Westerhausen R, Hugdahl K. The “paradoxical” engagement of the primary auditory cortex in patients with auditory verbal hallucinations: a meta-analysis of functional neuroimaging studies. Neuropsychologia (2011) 49(12):3361–9. doi:10.1016/j.neuropsychologia.2011.08.010
44. Kuhn S, Gallinat J. Quantitative meta-analysis on state and trait aspects of auditory verbal hallucinations in schizophrenia. Schizophr Bull (2012) 38(4):779–86. doi:10.1093/schbul/sbq152
45. Zmigrod L, Garrison JR, Carr J, Simons JS. The neural mechanisms of hallucinations: a quantitative meta-analysis of neuroimaging studies. Neurosci Biobehav Rev (2016) 69:113–23. doi:10.1016/j.neubiorev.2016.05.037
46. Moseley P, Alderson-Day B, Ellison A, Jardri R, Fernyhough C. Non-invasive brain stimulation and auditory verbal hallucinations: new techniques and future directions. Front Neurosci (2015) 9:515. doi:10.3389/fnins.2015.00515
47. Dierks T, Linden DE, Jandl M, Formisano E, Goebel R, Lanfermann H, et al. Activation of Heschl’s gyrus during auditory hallucinations. Neuron (1999) 22(3):615–21. doi:10.1016/S0896-6273(00)80715-1
48. Lennox BR, Park SB, Medley I, Morris PG, Jones PB. The functional anatomy of auditory hallucinations in schizophrenia. Psychiatry Res (2000) 100(1):13–20. doi:10.1016/S0925-4927(00)00068-8
49. Sommer IE, Diederen KM, Blom JD, Willems A, Kushan L, Slotema K, et al. Auditory verbal hallucinations predominantly activate the right inferior frontal area. Brain (2008) 131(Pt 12):3169–77. doi:10.1093/brain/awn251
50. van de Ven VG, Formisano E, Roder CH, Prvulovic D, Bittner RA, Dietz MG, et al. The spatiotemporal pattern of auditory cortical responses during verbal hallucinations. Neuroimage (2005) 27(3):644–55. doi:10.1016/j.neuroimage.2005.04.041
51. Diederen KM, Charbonnier L, Neggers SF, van Lutterveld R, Daalman K, Slotema CW, et al. Reproducibility of brain activation during auditory verbal hallucinations. Schizophr Res (2013) 146(1–3):320–5. doi:10.1016/j.schres.2013.01.025
52. Kanai R, Rees G. The structural basis of inter-individual differences in human behaviour and cognition. Nat Rev Neurosci (2011) 12(4):231–42. doi:10.1038/nrn3000
53. Peterchev AV, Rosa MA, Deng ZD, Prudic J, Lisanby SH. Electroconvulsive therapy stimulus parameters: rethinking dosage. J ECT (2010) 26(3):159–74. doi:10.1097/YCT.0b013e3181e48165
54. Groppa S, Oliviero A, Eisen A, Quartarone A, Cohen LG, Mall V, et al. A practical guide to diagnostic transcranial magnetic stimulation: report of an IFCN committee. Clin Neurophysiol (2012) 123(5):858–82. doi:10.1016/j.clinph.2012.01.010
55. Dedoncker J, Brunoni AR, Baeken C, Vanderhasselt MA. The effect of the interval-between-sessions on prefrontal transcranial direct current stimulation (tDCS) on cognitive outcomes: a systematic review and meta-analysis. J Neural Transm (Vienna) (2016) 123(10):1159–72. doi:10.1007/s00702-016-1558-x
Keywords: auditory hallucinations, neuroimaging, computational modeling, transcranial direct current stimulation, schizophrenia
Citation: Lee WH, Kennedy NI, Bikson M and Frangou S (2018) A Computational Assessment of Target Engagement in the Treatment of Auditory Hallucinations with Transcranial Direct Current Stimulation. Front. Psychiatry 9:48. doi: 10.3389/fpsyt.2018.00048
Received: 17 December 2017; Accepted: 05 February 2018;
Published: 22 February 2018
Edited by:
Christoph Mulert, University Medical Center Hamburg-Eppendorf, GermanyReviewed by:
Till R. Schneider, University of Hamburg, GermanyDaniel Keeser, Ludwig-Maximilians-Universität München, Germany
Copyright: © 2018 Lee, Kennedy, Bikson and Frangou. 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) and the copyright owner 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: Sophia Frangou, sophia.frangou@mssm.edu