Next Article in Journal
Advances in Personalized Medicines along with Functional Genomics and Pharmacogenomics
Next Article in Special Issue
Clinical Neurophysiology, Neuroimaging, and Neuromodulation of Neuropsychiatric Disorders
Previous Article in Journal
Predictors of the CD24/CD11b Biomarker among Healthy Subjects
Previous Article in Special Issue
Changes in EEG Alpha Activity during Attention Control in Patients: Association with Sleep Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Focality-Oriented Selection of Current Dose for Transcranial Direct Current Stimulation

by
Rajan Kashyap
1,*,
Sagarika Bhattacharjee
2,
Ramaswamy Arumugam
1,3,
Rose Dawn Bharath
4,
Kaviraja Udupa
5,
Kenichi Oishi
6,
John E. Desmond
6,†,
S. H. Annabel Chen
2,3,7,8,‡ and
Cuntai Guan
1,*,‡
1
School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore
2
School of Social Sciences (SSS), Nanyang Technological University, Singapore 639818, Singapore
3
Centre for Research and Development in Learning (CRADLE), Nanyang Technological University, Singapore 637460, Singapore
4
Department of Neuroimaging and Interventional Radiology, National Institute of Mental Health and Neurosciences, Hosur Road, Bangalore 560029, India
5
Department of Neurophysiology, National Institute of Mental Health and Neurosciences, Hosur Road, Bangalore 560029, India
6
The Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
7
Lee Kong Chian School of Medicine (LKC Medicine), Nanyang Technological University, Singapore 637553, Singapore
8
National Institute of Education, Nanyang Technological University, Singapore 637553, Singapore
*
Authors to whom correspondence should be addressed.
Senior author.
Equal contribution.
J. Pers. Med. 2021, 11(9), 940; https://doi.org/10.3390/jpm11090940
Submission received: 22 June 2021 / Revised: 18 September 2021 / Accepted: 18 September 2021 / Published: 21 September 2021

Abstract

:
Background: In transcranial direct current stimulation (tDCS), the injected current becomes distributed across the brain areas. The objective is to stimulate the target region of interest (ROI) while minimizing the current in non-target ROIs (the ‘focality’ of tDCS). For this purpose, determining the appropriate current dose for an individual is difficult. Aim: To introduce a dose–target determination index (DTDI) to quantify the focality of tDCS and examine the dose–focality relationship in three different populations. Method: Here, we extended our previous toolbox i-SATA to the MNI reference space. After a tDCS montage is simulated for a current dose, the i-SATA(MNI) computes the average (over voxels) current density for every region in the brain. DTDI is the ratio of the average current density at the target ROI to the ROI with a maximum value (the peak region). Ideally, target ROI should be the peak region, so DTDI shall range from 0 to 1. The higher the value, the better the dose. We estimated the variation of DTDI within and across individuals using T1-weighted brain images of 45 males and females distributed equally across three age groups: (a) young adults (20 ≤ x ˂ 40 years), (b) mid adults (40 ≤ x ˂ 60 years), and (c) older adults (60 ≤ x ˂ 80 years). DTDI’s were evaluated for the frontal montage with electrodes at F3 and the right supraorbital for three current doses of 1 mA, 2 mA, and 3 mA, with the target ROI at the left middle frontal gyrus. Result: As the dose is incremented, DTDI may show (a) increase, (b) decrease, and (c) no change across the individuals depending on the relationship (nonlinear or linear) between the injected tDCS current and the distribution of current density in the target ROI. The nonlinearity is predominant in older adults with a decrease in focality. The decline is stronger in males. Higher current dose at older age can enhance the focality of stimulation. Conclusion: DTDI provides information on which tDCS current dose will optimize the focality of stimulation. The recommended DTDI dose should be prioritized based on the age (>40 years) and sex (especially for males) of an individual. The toolbox i-SATA(MNI) is freely available.

1. Introduction

Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique that could alleviate symptoms of several neurological and psychiatric brain disorders [1,2,3]. A conventional tDCS setup consists of an anode and cathode placed over the scalp (referred to as a ‘montage’) with a low intensity of current (~1–3 mA) being injected to stimulate the target region of interest (ROI) [4,5]. However, the injected current becomes diffused in the intermediary regions of the brain and might not effectively stimulate the target ROI with the desired intensity [6,7]. Computational models that predict the pattern of current flow across the brain of an individual are used to optimize the tDCS stimulation parameters [8,9,10,11,12,13,14]. The amount of injected current (referred as the ‘current dose’) plays an important role in the dispersal of the stimulation’s intensity across the brain regions [15,16]. The distribution may vary from person to person and within a person based on the quantity of the dose [17,18,19]. Therefore, the selection of the optimal current dose for an individual’s brain that could sufficiently stimulate the target ROI while minimizing the current in non-target ROIs is important.
In recent years, there has been a growing interest in the individualization of the current dose [15,16,20]. It has been reported that varying the current intensity on the scalp for each individual can reduce the interindividual variability in the electric field intensity (or current density) at the target ROI [20]. The current dose calculated through inverse modelling of the tDCS-induced electric field at the target ROI correlates with the motor thresholds generated by transcranial magnetic stimulation [15]. In a recent tDCS experiment using a frontal montage and a 2 mA (fixed) current dose, individuals with a high current density at the target ROI (left dorsolateral prefrontal cortex) were found to have stronger improvements in working memory compared to those with a low current density [21]. They also showed that individualizing the current dose by fixing the desired current density at the target region can maximize the benefits of tDCS [21]. Though the models are a step towards individualizing the current dose, they do not consider the spread of the field to intermediary (non-target) regions. The current flow in the intermediary regions have a vital role to play in determining the outcome of tDCS [6,12,22,23,24,25]. It has been found that some brain regions may act as conduits, clustering most of the current to a specific location that can deter the intensity of the stimulation expected at the target ROI [6,26]. At this point, it is important to mention that other stimulation techniques (like peripheral nerve stimulation) are also intended to increase the stimulation intensity at target region while minimizing the stimulation received at non-target regions [27,28,29,30]. With tDCS, poor focality in stimulating the target region has constrained its efficacy. Therefore, the approaches to individualize the current dose [31,32,33,34] should consider the focality of stimulation in order to recommend the optimal intensity of input current.
In our previous work, we developed an individual-Systematic-Approach-for-tDCS-Analysis (i-SATA) toolbox [35] that estimates the average current density received by target ROIs and intermediary regions of an individual’s brain after a montage has been simulated in a realistic volumetric approach-based simulator for transcranial electric stimulation (ROAST) toolbox [10]. The ROAST-simulated current density in the ROIs has been found to be strongly correlated with electrophysiological measurements performed in vivo [9]. Integrated with ROAST, the i-SATA toolbox can be applied on an individual brain to reverse-calculate the current dose that can be used to stimulate the target ROI with the desired intensity [35]. This was performed based on the assumption that electric field intensity at the target ROI increases linearly with increase in the current dose by following the procedure laid down by Evans et al. [20]. Since we will be using it throughout the study, it will be helpful to familiarize our readers with an example. Suppose the calculated current density at the target ROI is 0.25 mA/m2 when 1 mA of current is applied on the scalp. To achieve a desired density of 0.5 mA/m2 at the target ROI, the required dosage (individualized) can be reverse-calculated as I n d i v i d u a l i s e d   d o s e = ( D e s i r e d   C u r r e n t   D e n s i t y A c t u a l   C u r r e n t   D e n s i t y ) × F i x e d   d o s e   [ i . e . , ( 0.5 0.25 ) × 1 = 2 mA].
In i-SATA, we used the Talairach client toolbox [36] to map an individual brain to the Talairach atlas space [37]. Another widely used brain template that provides detailed stereotaxic information on the location and variability of cortical areas is provided by the Montreal Neurological Institute (MNI) reference space [38,39,40,41,42]. Simon et al. [43,44,45] had developed the SPM anatomy toolbox that integrates the cytoarchitectonic maps in the MNI space. Here, we leveraged the potential of the SPM anatomy toolbox to extend i-SATA to the MNI space. The extended i-SATA(MNI) toolbox, which integrates the SPM anatomy toolbox with i-SATA, will enable researchers to visualize the comprehensive overview of the current density distribution across the cortex (target and intermediary regions) in the MNI space.
With i-SATA(MNI), we introduce the dose–target determination index (DTDI), a simple estimate that will quantify the focality of stimulation and facilitate the selection of optimal current dose required to stimulate the target ROI in an individual’s brain. A similar metric defined as the ‘selectivity index’ that measures the recruitment of the targeted region compared to other non-targeted regions is used to quantify the effectiveness of peripheral nerve stimulation [27,28,29,30]. For tDCS, the DTDI will aim to provide a comprehensive overview of the intensity of the stimulation received by the target ROI and intermediary regions after a montage has been postprocessed in i-SATA(MNI). To explain DTDI, we will use the montage with an anode positioned at F3 and a cathode at the right supraorbital (RSO) (referred to as F3-RSO). The montage has been shown to stimulate the left middle frontal gyrus [22,25] and is effective for depression [3,22,46] and working memory [47]. To make it easy for our readers to interpret how DTDI facilitates selection of the current dose, we will show the interindividual as well as the intraindividual variation in the index by uniformly increasing the current dose. Finally, we will evaluate the variation in DTDI by the age and sex of individuals. The purpose will be to explore if dose selection should be prioritised for any category (age and sex) of individuals.

2. Methods

2.1. Data

We obtained the T1-weighted (T1WI) magnetic resonance image (MRI) of the brain of 90 age–sex matched healthy individuals (45 male) from Cambridge Centre for Ageing and Neuroscience (Cam-CAN) study (available at http://www.mrc-cbu.cam.ac.uk/datasets/camcan/, accessed on 21 October 2020 [48,49]). This study was approved by the local ethics committee, Cambridgeshire 2, Research Ethics Committee (reference: 10/H0308/50). In this study, the Cam-CAN team recruited adult participants (aged 18–87 years old) in three stages that comprised of a home-based interview (stage one), followed by an evaluation of their health status (stage two). Subjects that were cognitively healthy (determined by a mini-mental state exam (MMSE) score ≥ 27), who met hearing, vision, and English language ability criteria, and who were free of MRI contraindications and neurologic or psychiatric conditions were recruited for stage three. In stage three, multimodal data (functional and structural MRI, magnetoencephalography, and behavioural) were collected from each participant. The T1WIs were collected from a 3T Siemens TIM Trio scanner with a 32-channel head coil using an MPRAGE sequence, TR = 2250 milliseconds (ms), TE = 2.99 ms, flip angle = 9°, Voxel size = 1 × 1 × 1 mm3, FOV = 256 × 240 × 192 mm3, GRAPPA: 2; TI: 900 ms. We selected 90 T1WIs from the following three age groups with 30 individuals (15 right-handed males and females) in each group:(a) young adults (20 ≤ x ˂ 40 years), (b) mid adults (40 ≤ x ˂ 60 years), and (c) older adults (60 ≤ x ˂ 80 years) were selected. The equal grouping across the three groups would allow an evaluation of the relationship of tDCS current dosage with sex and age.

2.2. Preprocessing with ROAST

We simulated the montage F3-RSO with the electrode size 5 × 5 cm2 (Figure 1A). For each individual MRI, the montages were simulated for three current doses of 1 mA, 2 mA, and 3 mA. In total, 270 simulations were performed in ROAST (Total = 90 MRI × 3 current doses = 270) [10]. Default conductivity values of the tissues (white matter (default 0.126 S/m); grey matter (default 0.276 S/m); cerebrospinal fluid (default 1.65 S/m); bone (default 0.01 S/m); skin (default 0.465 S/m); air (default 2.5 × 10−14 S/m); gel (default 0.3 S/m); electrode (default 5.9 × 107 S/m)) were used for each MRI that we simulated in ROAST. The ROAST simulation outputs the locations (x, y, and z coordinates) of the brain regions and the current density (mA/m2) value at each location in the native space.

2.3. i-SATA(MNI)

The i-SATA(MNI) is similar to i-SATA (available for download at https://doi.org/10.21979/N9/5W3RIM, [35]) except for the atlas space. In short, for each montage simulated in ROAST, i-SATA extracts the location (x, y, and z coordinates) of all the points in the cortex to detect the location of three anatomical landmarks (anterior commissure, posterior commissure, and mid-sagittal) using the acpcdetect toolbox [35,50]. With these landmarks, the individual’s native space was mapped to the reference space (i.e., the Talairach atlas space) using the fieldtrip toolbox [51] followed by the Talairach client toolbox [36]. Details on the methodology and application can be obtained from previous works [11,35,52]. For i-SATA(MNI), instead of the Talairach atlas space, we mapped the outputs (x, y, and z coordinates) to the MNI reference space using the SPM anatomy toolbox [43,44,45]. The SPM anatomy toolbox has an option for using the gyri/sulci-based labelling system wherein the automated anatomical labeling atlas with 116 regions outlined on the Colin27 brain template is implemented (for details, see [53]). The i-SATA(MNI) extracts and uses the labels provided by this atlas for assigning the cortical and subcortical region corresponding to each location. A detailed explanation on the nomenclature of the delineated regions can be found at [53]. We developed i-SATA(MNI) using SPM12 (revision 6470, available at https://www.fil.ion.ucl.ac.uk/spm/software/spm12/) that has the SPM Anatomy toolbox (version 2.2b) inbuilt into the framework. The magnitude of current density corresponding to each location (voxel) was then used to obtain the average magnitude of the current density received by each cortical region of the brain. This will provide an estimate of the current density induced in the target and intermediary region due to tDCS. As an example, we will postprocess the standard MNI 152 averaged head in i-SATA(MNI) for the three current doses (1 mA, 2 mA, and 3 mA) using the F3-RSO montage to show the distribution of the average current density across the cortical regions (Figure 1B–D).

2.4. Dose–Target Determination Index (DTDI)

The output of i-SATA(MNI) (i.e., the average current density in the target ROI and the non-target regions) is used to calculate the DTDI for a montage simulated at a current dose. For this, we will find the ROI that has the maximum value of average current density (peak region) amongst all the ROIs. DTDI is then calculated as
D T D I = A v e r a g e   C u r r e n t   d e n s i t y   a t   t h e   T a r g e t   R O I M a x i m u m   v a l u e   o f   a v e r a g e   c u r r e n t   d e n s i t y   f o r m e d   a t   a n y   R O I
DTDI will lie within the range of 0 to 1. An ideal tDCS setup will deliver the maximum intensity of stimulation (average current density) to the target ROI, thereby generating a DTDI value equal to 1. However, the peak intensity may be received at non-targeted ROI. A DTDI value = 0 will indicate no stimulation of the target ROI. For an individual, the current dose for which DTDI is higher should be preferred over other doses. To make this clear, we will estimate the DTDI of three individuals across three current doses. Hypothetically, the value of the DTDI should remain constant across doses, since it is assumed that the current flow in the brain increases linearly with an increase in current intensity [15,16,20,23,54].

2.5. Statistical Analysis of Variation in DTDI

All individual MRIs were post-processed in i-SATA(MNI) for the three current doses using the F3-RSO montage to estimate the DTDI’s (Total = 90 MRI × 3 current doses = 270). We show the interindividual and intraindividual variation in the DTDI for both sexes across the three age groups. We performed a three-way mixed ANOVA with age and sex as between subject and dose as within subject factor. Post hoc analyses were performed to further characterize the nature of the main effects and interactions.

2.6. Code Availability

The i-SATA(MNI) is a Linux-based MATLAB toolbox integrating acpcdetect v2.0, fieldtrip, and SPM12 (version 6470) with the integrated SPM Anatomy toolbox (version 2.2b). The package can be downloaded at (https://doi.org/10.21979/N9/KWTCWK). A reference manual is also provided to help users run each step with ease.

3. Results

3.1. Output of i-SATA(MNI) on the Standard Head Model

The montage F3-RSO that we applied to the MNI 152 averaged head model and simulated in ROAST is shown in Figure 1A. The output of i-SATA(MNI), i.e., the distribution of the average current density across the cortical regions, are shown in Figure 1B–D for the three current doses (1 mA, 2 mA, and 3 mA). The average current density in the target ROI (the left middle frontal gyrus) varies linearly with the current dose. Therefore, the DTDI remains constant (at approximately ~ 0.85) across the doses. Of note, similarly to i-SATA [35] and SATA (the standard head model with a graphical user interface is available for download at https://doi.org/10.21979/N9/DMWPZK [11,52,55]), users can visualize the i-SATA(MNI) outputs on the brain surface as well (Figure not shown).

3.2. Interpretation of DTDI for Appropriate Selection of Current Dose

For any individual, DTDI can guide the selection of the appropriate current dose that will sufficiently stimulate the target ROI with minimal spread of current to other regions. For interpretation, we have shown the variation of DTDI for three individuals across the three current doses (Figure 2). For the first individual, the current intensity at the target region increases with the increase in dose, and the DTDI remains fairly constant (Figure 2A). This implies that the target ROI will be sufficiently stimulated by any current dose, and that the user can tune it according to the extent of stimulation desired. For the second individual, a low DTDI (0.43) is seen at lower dose (1 mA) suggesting that target ROI is receiving minimal current and non-target regions are receiving most of the current. With the increase in dose, it can be seen that the current intensity at the target ROI is increasing, but fewer regions are receiving a current higher than the target ROI. As a result, the DTDI is increasing with increase in the dose, suggesting that higher current dose should be beneficial (Figure 2B). Finally for the third individual, a decrease in DTDI is seen with the increase in dose (Figure 2C). The drop in DTDI from 1 mA to 2 and 3 mA seems to be due to an increase in current in the right superior parietal lobule at 2 mA and 3 mA only. Although the current intensity at the target ROI increases with the increases in the dose, the maximal amount of current also becomes dissipated to other brain regions. Thus, the conventional way of increasing the current dose to attain the desired stimulation intensity at the target ROI might result in the stimulation of unwanted brain regions (as seen for the superior parietal lobule). For this individual, a lower dose showing higher DTDI can maximize the advantages of stimulation.

3.3. Statistical Analysis of Variance in DTDI

The change in DTDI as a function of the dose for males and females across three age groups is shown in Figure 3. The mixed ANOVA revealed a significant main effect of age (F(2, 84) = 43.98, p = 8.51 × 10−14, effect size = 0.405 (generalised eta squared)), with DTDI significantly decreasing in older adults compared to young adults (p = 0.0008). The main effect of sex (F(1, 84) = 12.14, p < 7.85 × 10−4, effect size = 0.086) and its interaction with age (F(2, 84) = 3.78, p= 2.70 × 10−2, effect size = 0.05) was also found to be significant. The post hoc analysis shows that females had higher DTDI values than males for both mid (p = 2.89 × 10−6) and older adults (p = 3.18 × 10−3). The interaction effect of age and dose was also found to be significant (F(3.34, 140.48) = 7.269, p = 7.65 × 10−5, effect size = 0.05). In older adults only, the post hoc analysis revealed that there is a significant increase in the DTDI values at 3 mA compared to 1 mA for both males (p = 0.01) and females (p = 0.02). All the post hoc pairwise comparisons were Bonferroni-corrected. This shows that the focality of stimulation could be enhanced in older adults by increasing the dose.

4. Discussion

In this paper, we extended the toolbox i-SATA [35] to the MNI reference space to enable users to obtain the average current density induced at each cortical region of an individual’s brain during tDCS. We then used these values to estimate the DTDI as an objective measure to quantify the focality of stimulation and aid the selection of appropriate current dose for an individual. We demonstrated the utility of DTDI across three subjects and doses (Figure 2). The optimal stimulation of the target ROI was experienced differently across the three subjects: (a) subject 1 was neutral to changes in current dose, and the electric field intensity remained proportional to the injected current; (b) subject 2 had better focality from a dose of 2 mA or more (but not from 1 mA); and (c) subject 3 gained better stimulation from 1 mA compared to 2 or 3 mA of current dose. Such interindividual inconsistency in tDCS due to the current dose has been widely reported in previous studies [17,18,19,56,57,58,59,60,61,62]. With this i-SATA(MNI) framework post-processing the structural scans simulated in ROAST, tDCS users can configure personalized protocols for montage selection (refer to [11,35]) and identify the optimal current dose for cortical targeting (guided by DTDI). We applied the framework on a wide age range (20 to 80 years) of individuals from both sexes to highlight the importance of DTDI and the need for a focality-based selection of the current dose.
Previous tDCS based studies have combined electroencephalography, functional MRI, or transcranial magnetic stimulation to determine the current dose for optimal targeting [63,64,65,66]. Recently, a computational study has put forward the model to reverse calculate the current dose from the simulated electric field [16] based on the assumption that the intensity of current flow increases linearly with current dose [54]. A similar prototype was also put forward by Evans et al. [20]. All these studies used young, healthy subjects to delineate the model. On one hand, we see this linearity (constant value of DTDI) being followed in the MNI standard head model (Figure 1) and to an extent in young and middle-aged individuals (Figure 2A and Figure 3A). However, on the other hand, linearity appears to diminish with advancement in age (Figure 2B,C and Figure 3A), suggesting a potential nonlinear relationship.
The different values of DTDI as a function of current dose across different subjects could be because the injected current might become clustered in brain areas (referred to as hotspots), a phenomenon that has been widely reported in tDCS studies [6,25,26,67,68,69]. Hotspots cause shunting of the current towards the surrounding brain tissue and a surge in the electric field strength at localised areas [70]. Areas that form hotspots can be away from the electrode site as well [70]. In the two cases presented in Figure 2 (Subject 2 and 3), the superior parietal lobule appears to be a hotspot. Here, it is difficult to comprehend the neuroanatomical factors that contribute to the formation of such hotspots. It has been found that tissue heterogeneities and pathological alterations (like neurodegeneration and cerebral infarcts) are the primary contributors [26,70]. As we age, atrophy in the neural configuration escalates the nonlinearities in the spatial distribution of induced electric fields [71,72]. Care must be taken to limit the possibilities of such hotspots forming for clinical application of tDCS [73]. DTDI that considers the current density in target and non-target areas incorporates an understanding of the effect of hotspots to provide a rigorous estimate of the optimal current dose.
Since maximum stimulation might not be received at the target ROI, nor in a consistent location [6,17,18,74], the interindividual and intraindividual variation in DTDI can provide insights for the appropriate determination of current dose based on the age and sex of a healthy individual. In young adults, the focality of stimulation remains intact (approximately) across the doses, indicating that there is flexibility in choosing (individualizing) a dose depending on the extent of current density desired at the target ROI. However, the focality declines with advancing age (middle age onwards, see Figure 3). This decline is higher for males compared to females. Such sexual dimorphism in tDCS-related effects have been reported in previous studies [75] and several factors related to cortical anatomy like volume [76], bone density [77], hormonal levels [77], and electrode location have been postulated to account for it. We also found that higher current doses can enhance the focality in older adults. This is in support of a recent study [78] which reported that cerebral atrophy in older adults was causing the reduction in the amount of current reaching the target ROI. With that being said, we feel that it is important to guide our readers on how to use DTDI in their tDCS study (at an individual or group level), wherein their interest may lie in stimulating single or multiple target ROIs. We will discuss this in the next section.

How to Use DTDI

  • (A). DTDI for Group-level Study
We have shown how DTDI can be used to titrate the current dose at the individual level. This can be performed at the group level also. Evans et al. [20] have suggested that the input current should be varied across individuals to maintain a constant current density at the target ROI. While we agree with them, we also suggest that the focality of the stimulation needs to be considered, especially when older individuals are recruited for the study. For primary clinical/therapeutic applications of tDCS, the focality, as revealed by DTDI, could be especially useful for setting tDCS dosage. Although the compatibility of the patient with the computationally recommended dose is always important [79], recent studies have indicated that participants readily tolerate tDCS current up to 4 mA [80,81].
In group studies in which researchers do not want to vary the current from subject to subject, DTDI values may still be used in two different ways to improve the efficacy of the study. The first would be to include a threshold for DTDI (e.g., DTDI ≥ 0.75) as a precautionary measure when individualizing the current dose. We recommend this value because studies pertaining to peripheral nerve stimulation have also found that such a threshold value (≥ 0.75) for the ‘selectivity index’ could ensure that the targeted region of the nerve is well activated [28,30]. While this may narrow down the suitability of subjects, such inclusion criteria could reduce the variability of tDCS. The second would be to use DTDI analyses for the populations under study to determine, at the start of the study, the optimal value of the tDCS current dose to be used on all subjects that will produce the greatest focality and the least amount of subject–subject variability in DTDI. For example, the current study suggests that for the F3-RSO montage, if you are including older and younger subjects, a higher current value (for the study overall) might produce the least variability in terms of focality of tDCS.
  • (B). DTDI for Multiple Target ROIs
We have shown how DTDI can be used for a target ROI. However, a user may be interested in stimulating multiple target ROIs, or a network of ROIs specific to a particular functionality [82,83,84,85,86,87,88]. This may be because studies have shown that more than one brain area can be involved in neuropsychiatric disorders like depression [89,90]. Alternatively (in an unexpected scenario), it may also happen that two or more ROIs receive the same current density after a montage is simulated for an individual’s brain. In such scenarios, the DTDI can be averaged across all the target ROIs ( D T D I   f o r   f i r s t   T a r g e t R O I + D T D I   f o r   s e c o n d   T a r g e t R O I + .. + D T D I   f o r   n t h   t a r g e t R O I N u m b e r   o f   a r e a s   ( n )   ). The averaged DTDI will also lie in between 0 and 1. The optimal dose should focus on determining the current at which the value of the averaged DTDI exhibits an increasing trend.

5. Limitations

Here, we show how DTDI aids selection of focality-based current dose in tDCS. Although the findings on the relationship between tDCS dose, sex, and age could be limited by the small sample size, they reveal that, while brains which have a linear relationship between the electric field intensity and the tDCS current are flexible to the selection of current dose, brains that have a nonlinear relationship are dose sensitive. Therefore, biomarkers that could identify brains that can behave nonlinearly prior to the onset of an experiment could improve the precision of tDCS. However, this will require large datasets for validation of the biomarker. We leave it for future studies to investigate the exact nature of such nonlinearities and to determine the factors (structural, functional, and behavioural) that contribute to them.
Finally, we would like to highlight that DTDI can be estimated from i-SATA as well. The Talairach atlas space has several benefits [36], and similarly the i-SATA(MNI) that was developed for users who prefer the MNI space has its own advantages [43]. However, simulation in i-SATA(MNI) is considerably faster than in i-SATA. This is because both i-SATA(MNI) and the integrated SPM anatomy toolbox for cortical labelling are MATLAB-based and automated. This makes i-SATA(MNI) efficient at post-processing large data sets, a trend that is emerging in neuroscientific research.

6. Conclusions

The study extends the i-SATA framework to the MNI atlas space. With i-SATA(MNI), it will be easier to calculate the individualized dose as suggested in previous studies [15,16,20,21]. Here, we introduce the DTDI as a measure to titrate the individualized current doses and select the optimum dose that has high focality and could appropriately stimulate the target ROI in an individual. This will facilitate the personalized application of tDCS so that the desired stimulation benefits are achieved. Using a montage that has been found to be optimal for a DLPFC stimulation, DTDI analysis across a broad spectrum of men and women of different age groups has revealed that focality decreases with advancing age, especially in males with more than 40 years of age. Finally, the study reveals that the selection of a current dose that increases the focality is strictly necessary for older (> 60 years) individuals irrespective of sex.

Author Contributions

Conceptualization, R.K. and S.B.; methodology, R.K. and R.A.; software, R.A. and R.K.; validation, R.A., S.B. and R.K.; formal analysis, R.K. and S.B.; investigation, R.K.; resources, S.H.A.C. and C.G.; data curation, S.B.; writing—original draft preparation, R.K.; writing—review and editing, J.E.D., S.H.A.C., K.O., R.D.B., K.U. and C.G.; visualization, K.O. and J.E.D.; supervision, S.H.A.C. and C.G.; project administration, R.K. and C.G.; funding acquisition, C.G. All authors have read and agreed to the published version of the manuscript.

Funding

The work is funded by the grant RIE2020 AME Programmatic Fund, Singapore (No. A20G8b0102).

Institutional Review Board Statement

Not Applicable to the study.

Informed Consent Statement

Not Applicable to the study.

Data Availability Statement

We would like to convey our gratitude to the CAM-CAN team (https://camcan-archive.mrc-cbu.cam.ac.uk/dataaccess/) for providing access to the dataset.

Acknowledgments

We acknowledge the contribution from the RIE2020 AME Programmatic Fund, Singapore (No. A20G8b0102) to our work. The work was also supported by the NTU-JHU grant from Nanyang Technological University, Singapore. J.D. received additional support from the NIH/NICHD grant P50 HD103538. R.K. and S.B. have support from the DBT Ramalingaswami Re-entry fellowship scheme (2021), sponsored by the Government of India.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Antal, A.; Alekseichuk, I.; Bikson, M.; Brockmöller, J.; Brunoni, A.R.; Chen, R.; Cohen, L.G.; Dowthwaite, G.; Ellrich, J.; Flöel, A. Low intensity transcranial electric stimulation: Safety, ethical, legal regulatory and application guidelines. Clin. Neurophysiol. 2017, 128, 1774–1809. [Google Scholar] [CrossRef] [Green Version]
  2. Filmer, H.L.; Dux, P.E.; Mattingley, J.B. Applications of transcranial direct current stimulation for understanding brain function. Trends Neurosci. 2014, 37, 742–753. [Google Scholar] [CrossRef] [Green Version]
  3. Razza, L.B.; Palumbo, P.; Moffa, A.H.; Carvalho, A.F.; Solmi, M.; Loo, C.K.; Brunoni, A.R. A systematic review and meta-analysis on the effects of transcranial direct current stimulation in depressive episodes. Depress. Anxiety 2020, 37, 594–608. [Google Scholar] [CrossRef] [PubMed]
  4. Nitsche, M.A.; Doemkes, S.; Karakose, T.; Antal, A.; Liebetanz, D.; Lang, N.; Tergau, F.; Paulus, W. Shaping the effects of transcranial direct current stimulation of the human motor cortex. J. Neurophysiol. 2007, 97, 3109–3117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Nitsche, M.A.; Paulus, W. Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation. J. Physiol. 2000, 527, 633–639. [Google Scholar] [CrossRef] [PubMed]
  6. Datta, A.; Bansal, V.; Diaz, J.; Patel, J.; Reato, D.; Bikson, M. Gyri-precise head model of transcranial direct current stimulation: Improved spatial focality using a ring electrode versus conventional rectangular pad. Brain Stimul. 2009, 2, 201–207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Batsikadze, G.; Moliadze, V.; Paulus, W.; Kuo, M.-F.; Nitsche, M.A. Partially non-linear stimulation intensity-dependent effects of direct current stimulation on motor cortex excitability in humans. J. Physiol. 2013, 591, 1987–2000. [Google Scholar] [CrossRef] [PubMed]
  8. Indahlastari, A.; Chauhan, M.; Sadleir, R. Benchmarking transcranial electrical stimulation finite element simulations: A comparison study. J. Neural Eng. 2019, 16, 026019. [Google Scholar] [CrossRef] [PubMed]
  9. Huang, Y.; Liu, A.A.; Lafon, B.; Friedman, D.; Dayan, M.; Wang, X.; Bikson, M.; Doyle, W.K.; Devinsky, O.; Parra, L.C. Measurements and models of electric fields in the in vivo human brain during transcranial electric stimulation. eLife 2017, 6, e18834. [Google Scholar] [CrossRef] [PubMed]
  10. Huang, Y.; Datta, A.; Bikson, M.; Parra, L.C. Realistic volumetric-approach to simulate transcranial electric stimulation—ROAST—a fully automated open-source pipeline. J. Neural Eng. 2019, 16, 056006. [Google Scholar] [CrossRef] [PubMed]
  11. Bhattacharjee, S.; Kashyap, R.; Rapp, B.; Oishi, K.; Desmond, J.E.; Chen, S.A. Simulation Analyses of tDcS Montages for the investigation of Dorsal and Ventral pathways. Sci. Rep. 2019, 9, 1–17. [Google Scholar] [CrossRef]
  12. Bikson, M.; Datta, A.; Rahman, A.; Scaturro, J. Electrode montages for tDCS and weak transcranial electrical stimulation: Role of “return” electrode’s position and size. Clin. Neurophysiol. Off. J. Int. Fed. Clin. Neurophysiol. 2010, 121, 1976. [Google Scholar] [CrossRef] [Green Version]
  13. Opitz, A.; Falchier, A.; Yan, C.-G.; Yeagle, E.M.; Linn, G.S.; Megevand, P.; Thielscher, A.; Milham, M.P.; Mehta, A.D.; Schroeder, C.E. Spatiotemporal structure of intracranial electric fields induced by transcranial electric stimulation in humans and nonhuman primates. Sci. Rep. 2016, 6, 1–11. [Google Scholar]
  14. Laakso, I.; Mikkonen, M.; Koyama, S.; Hirata, A.; Tanaka, S. Can electric fields explain inter-individual variability in transcranial direct current stimulation of the motor cortex? Sci. Rep. 2019, 9, 1–10. [Google Scholar]
  15. Caulfield, K.A.; Badran, B.W.; Li, X.; Bikson, M.; George, M.S. Can transcranial electrical stimulation motor threshold estimate individualized tDCS doses over the prefrontal cortex? Evidence from reverse-calculation electric field modeling. Brain Stimul. 2020, 13, 1150–1152. [Google Scholar] [CrossRef]
  16. Caulfield, K.A.; Badran, B.W.; DeVries, W.H.; Summers, P.M.; Kofmehl, E.; Li, X.; Borckardt, J.J.; Bikson, M.; George, M.S. Transcranial electrical stimulation motor threshold can estimate individualized tDCS dosage from reverse-calculation electric-field modeling. Brain Stimul. 2020, 13, 961–969. [Google Scholar] [CrossRef] [PubMed]
  17. Chew, T.; Ho, K.-A.; Loo, C.K. Inter-and intra-individual variability in response to transcranial direct current stimulation (tDCS) at varying current intensities. Brain Stimul. 2015, 8, 1130–1137. [Google Scholar] [CrossRef]
  18. López-Alonso, V.; Fernández-del-Olmo, M.; Costantini, A.; Gonzalez-Henriquez, J.J.; Cheeran, B. Intra-individual variability in the response to anodal transcranial direct current stimulation. Clin. Neurophysiol. 2015, 126, 2342–2347. [Google Scholar] [CrossRef] [PubMed]
  19. López-Alonso, V.; Cheeran, B.; Río-Rodríguez, D.; Fernández-del-Olmo, M. Inter-individual variability in response to non-invasive brain stimulation paradigms. Brain Stimul. 2014, 7, 372–380. [Google Scholar] [CrossRef] [PubMed]
  20. Evans, C.; Bachmann, C.; Lee, J.S.A.; Gregoriou, E.; Ward, N.; Bestmann, S. Dose-controlled tDCS reduces electric field intensity variability at a cortical target site. Brain Stimul. 2020, 13, 125–136. [Google Scholar] [CrossRef] [Green Version]
  21. Caulfield, K.A.; Indahlastari, A.; Nissim, N.R.; Lopez, J.W.; Fleischmann, H.H.; Woods, A.J.; George, M.S. Electric Field Strength From Prefrontal Transcranial Direct Current Stimulation Determines Degree of Working Memory Response: A Potential Application of Reverse-Calculation Modeling? Neuromodulation 2020. [Google Scholar] [CrossRef]
  22. Ammann, C.; Lindquist, M.A.; Celnik, P.A. Response variability of different anodal transcranial direct current stimulation intensities across multiple sessions. Brain Stimul. 2017, 10, 757–763. [Google Scholar] [CrossRef]
  23. Esmaeilpour, Z.; Marangolo, P.; Hampstead, B.M.; Bestmann, S.; Galletta, E.; Knotkova, H.; Bikson, M. Incomplete evidence that increasing current intensity of tDCS boosts outcomes. Brain Stimul. 2018, 11, 310–321. [Google Scholar] [CrossRef]
  24. Antal, A.; Polania, R.; Schmidt-Samoa, C.; Dechent, P.; Paulus, W. Transcranial direct current stimulation over the primary motor cortex during fMRI. Neuroimage 2011, 55, 590–596. [Google Scholar] [CrossRef] [PubMed]
  25. Bai, S.; Dokos, S.; Ho, K.-A.; Loo, C. A computational modelling study of transcranial direct current stimulation montages used in depression. Neuroimage 2014, 87, 332–344. [Google Scholar] [CrossRef]
  26. Wagner, T.; Fregni, F.; Fecteau, S.; Grodzinsky, A.; Zahn, M.; Pascual-Leone, A. Transcranial direct current stimulation: A computer-based human model study. NeuroImage 2007, 35, 1113–1124. [Google Scholar] [CrossRef] [PubMed]
  27. Deurloo, K.E.I.; Holsheimer, J.; Bergveld, P. Nerve stimulation with a multi-contact cuff electrode: Validation of model predictions. Arch. Physiol. Biochem. 2000, 108, 349–359. [Google Scholar] [CrossRef] [PubMed]
  28. Kosta, P.; Warren, D.J.; Lazzi, G. Selective stimulation of rat sciatic nerve using an array of mm-size magnetic coils: A simulation study. Healthc. Technol. Lett. 2019, 6, 70–75. [Google Scholar] [CrossRef]
  29. Tigra, W.; Guiraud, D.; Andreu, D.; Coulet, B.; Gelis, A.; Fattal, C.; Maciejasz, P.; Picq, C.; Rossel, O.; Teissier, J. Exploring selective neural electrical stimulation for upper limb function restoration. Eur. J. Transl. Myol. 2016, 26, 6035. [Google Scholar] [CrossRef] [Green Version]
  30. Kosta, P.; Mize, J.; Warren, D.J.; Lazzi, G. Simulation-Based Optimization of Figure-of-Eight Coil Designs and Orientations for Magnetic Stimulation of Peripheral Nerve. IEEE Trans. Neural Syst. Rehabil. Eng. 2020, 28, 2901–2913. [Google Scholar] [CrossRef]
  31. Mikkonen, M.; Laakso, I.; Tanaka, S.; Hirata, A. Cost of focality in TDCS: Interindividual variability in electric fields. Brain Stimul. 2020, 13, 117–124. [Google Scholar] [CrossRef] [Green Version]
  32. Opitz, A.; Yeagle, E.; Thielscher, A.; Schroeder, C.; Mehta, A.D.; Milham, M.P. On the importance of precise electrode placement for targeted transcranial electric stimulation. Neuroimage 2018, 181, 560–567. [Google Scholar] [CrossRef] [PubMed]
  33. Fischer, D.B.; Fried, P.J.; Ruffini, G.; Ripolles, O.; Salvador, R.; Banus, J.; Ketchabaw, W.T.; Santarnecchi, E.; Pascual-Leone, A.; Fox, M.D. Multifocal tDCS targeting the resting state motor network increases cortical excitability beyond traditional tDCS targeting unilateral motor cortex. Neuroimage 2017, 157, 34–44. [Google Scholar] [CrossRef]
  34. Bortoletto, M.; Rodella, C.; Salvador, R.; Miranda, P.C.; Miniussi, C. Reduced current spread by concentric electrodes in transcranial electrical stimulation (tES). Brain Stimul. 2016, 9, 525–528. [Google Scholar] [CrossRef]
  35. Kashyap, R.; Bhattacharjee, S.; Arumugam, R.; Oishi, K.; Desmond, J.E.; Chen, S.A. i-SATA: A MATLAB based toolbox to estimate current density generated by transcranial direct current stimulation in an individual brain. J. Neural Eng. 2020, 17, 056034. [Google Scholar] [CrossRef]
  36. Lancaster, J.L.; Woldorff, M.G.; Parsons, L.M.; Liotti, M.; Freitas, C.S.; Rainey, L.; Kochunov, P.V.; Nickerson, D.; Mikiten, S.A.; Fox, P.T. Automated Talairach Atlas labels for functional brain mapping. Hum. Brain Mapp. 2000, 10, 120–131. [Google Scholar] [CrossRef]
  37. Talairach, J.; Tournoux, P. Co-Planar Stereotaxic Atlas of the Human Brain 3-Dimensional Proportional System: An Approach to Cerebral Imaging; Georg Thieme Verlag: Stuttgart, Germany; New York, NY, USA, 1988. [Google Scholar]
  38. Collins, D.L.; Neelin, P.; Peters, T.M.; Evans, A.C. Automatic 3D intersubject registration of MR volumetric data in standardized Talairach space. J. Comput. Assist. Tomogr. 1994, 18, 192–205. [Google Scholar] [CrossRef]
  39. Evans, A.C.; Marrett, S.; Neelin, P.; Collins, L.; Worsley, K.; Dai, W.; Milot, S.; Meyer, E.; Bub, D. Anatomical mapping of functional activation in stereotactic coordinate space. Neuroimage 1992, 1, 43–53. [Google Scholar] [CrossRef]
  40. Holmes, C.J.; Hoge, R.; Collins, L.; Woods, R.; Toga, A.W.; Evans, A.C. Enhancement of MR images using registration for signal averaging. J. Comput. Assist. Tomogr. 1998, 22, 324–333. [Google Scholar] [CrossRef] [PubMed]
  41. Amunts, K.; Zilles, K. Advances in cytoarchitectonic mapping of the human cerebral cortex. Neuroimaging Clin. North Am. 2001, 11, 151–169. [Google Scholar]
  42. Zilles, K.; Schleicher, A.; Palomero-Gallagher, N.; Amunts, K. Quantitative analysis of cyto-and receptor architecture of the human brain. In Brain Mapping: The Methods; Elsevier: Amsterdam, The Netherlands, 2002; pp. 573–602. [Google Scholar]
  43. Eickhoff, S.B.; Stephan, K.E.; Mohlberg, H.; Grefkes, C.; Fink, G.R.; Amunts, K.; Zilles, K. A new SPM toolbox for combining probabilistic cytoarchitectonic maps and functional imaging data. Neuroimage 2005, 25, 1325–1335. [Google Scholar] [CrossRef] [PubMed]
  44. Eickhoff, S.B.; Heim, S.; Zilles, K.; Amunts, K. Testing anatomically specified hypotheses in functional imaging using cytoarchitectonic maps. Neuroimage 2006, 32, 570–582. [Google Scholar] [CrossRef] [PubMed]
  45. Eickhoff, S.B.; Paus, T.; Caspers, S.; Grosbras, M.-H.; Evans, A.C.; Zilles, K.; Amunts, K. Assignment of functional activations to probabilistic cytoarchitectonic areas revisited. Neuroimage 2007, 36, 511–521. [Google Scholar] [CrossRef]
  46. Moffa, A.H.; Martin, D.; Alonzo, A.; Bennabi, D.; Blumberger, D.M.; Benseñor, I.M.; Daskalakis, Z.; Fregni, F.; Haffen, E.; Lisanby, S.H. Efficacy and acceptability of transcranial direct current stimulation (tDCS) for major depressive disorder: An individual patient data meta-analysis. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 2020, 99, 109836. [Google Scholar] [CrossRef] [PubMed]
  47. Nikolin, S.; Martin, D.; Loo, C.K.; Boonstra, T.W. Effects of TDCS dosage on working memory in healthy participants. Brain Stimul. 2018, 11, 518–527. [Google Scholar] [CrossRef] [PubMed]
  48. Taylor, J.R.; Williams, N.; Cusack, R.; Auer, T.; Shafto, M.A.; Dixon, M.; Tyler, L.K.; Henson, R.N. The Cambridge Centre for Ageing and Neuroscience (Cam-CAN) data repository: Structural and functional MRI, MEG, and cognitive data from a cross-sectional adult lifespan sample. Neuroimage 2017, 144, 262–269. [Google Scholar] [CrossRef] [PubMed]
  49. Shafto, M.A.; Tyler, L.K.; Dixon, M.; Taylor, J.R.; Rowe, J.B.; Cusack, R.; Calder, A.J.; Marslen-Wilson, W.D.; Duncan, J.; Dalgleish, T. The Cambridge Centre for Ageing and Neuroscience (Cam-CAN) study protocol: A cross-sectional, lifespan, multidisciplinary examination of healthy cognitive ageing. BMC Neurol. 2014, 14, 204. [Google Scholar] [CrossRef] [Green Version]
  50. Ardekani, B.A.; Bachman, A.H. Model-based automatic detection of the anterior and posterior commissures on MRI scans. Neuroimage 2009, 46, 677–682. [Google Scholar] [CrossRef] [Green Version]
  51. Oostenveld, R.; Fries, P.; Maris, E.; Schoffelen, J.-M. FieldTrip: Open Source Software for Advanced Analysis of MEG, EEG, and Invasive Electrophysiological Data. Comput. Intell. Neurosci. 2011, 2011, 156869. [Google Scholar] [CrossRef]
  52. Bhattacharjee, S.; Kashyap, R.; O’Brien, B.A.; McCloskey, M.; Oishi, K.; Desmond, J.E.; Rapp, B.; Chen, S.H.A. Reading proficiency influences the effects of transcranial direct current stimulation: Evidence from selective modulation of dorsal and ventral pathways of reading in bilinguals. Brain Lang. 2020, 210, 104850. [Google Scholar] [CrossRef]
  53. Tzourio-Mazoyer, N.; Landeau, B.; Papathanassiou, D.; Crivello, F.; Etard, O.; Delcroix, N.; Mazoyer, B.; Joliot, M. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 2002, 15, 273–289. [Google Scholar] [CrossRef]
  54. Bikson, M.; Truong, D.Q.; Mourdoukoutas, A.P.; Aboseria, M.; Khadka, N.; Adair, D.; Rahman, A. Modeling sequence and quasi-uniform assumption in computational neurostimulation. Prog. Brain Res. 2015, 222, 1–23. [Google Scholar] [PubMed]
  55. Bhattacharjee, S.; Chew, A.; Kashyap, R.; Wu, C.; Yeo, M.; O’Brien, B.; Rapp, B.; McCloskey, M.; Oishi, K.; Desmond, J.; et al. Could tDCS Modulate Bilingual Reading? Brain Stimul. 2019, 12, 569. [Google Scholar] [CrossRef]
  56. Horvath, J.C.; Carter, O.; Forte, J.D. Transcranial direct current stimulation: Five important issues we aren’t discussing (but probably should be). Front. Syst. Neurosci. 2014, 8, 2. [Google Scholar] [CrossRef] [PubMed]
  57. Horvath, J.C.; Forte, J.D.; Carter, O. Quantitative review finds no evidence of cognitive effects in healthy populations from single-session transcranial direct current stimulation (tDCS). Brain Stimul. 2015, 8, 535–550. [Google Scholar] [CrossRef]
  58. Horvath, J.C.; Carter, O.; Forte, J.D. No significant effect of transcranial direct current stimulation (tDCS) found on simple motor reaction time comparing 15 different simulation protocols. Neuropsychologia 2016, 91, 544–552. [Google Scholar] [CrossRef]
  59. Goldsworthy, M.R.; Hordacre, B. Dose dependency of transcranial direct current stimulation: Implications for neuroplasticity induction in health and disease. J. Physiol. 2017, 595, 3265–3266. [Google Scholar] [CrossRef] [PubMed]
  60. Jamil, A.; Batsikadze, G.; Kuo, H.-I.; Labruna, L.; Hasan, A.; Paulus, W.; Nitsche, M.A. Systematic evaluation of the impact of stimulation intensity on neuroplastic after-effects induced by transcranial direct current stimulation. J. Physiol. 2017, 595, 1273–1288. [Google Scholar] [CrossRef]
  61. Labruna, L.; Jamil, A.; Fresnoza, S.; Batsikadze, G.; Kuo, M.-F.; Vanderschelden, B.; Ivry, R.B.; Nitsche, M.A. Efficacy of anodal transcranial direct current stimulation is related to sensitivity to transcranial magnetic stimulation. Brain Stimul. 2016, 9, 8–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Kidgell, D.J.; Daly, R.M.; Young, K.; Lum, J.; Tooley, G.; Jaberzadeh, S.; Zoghi, M.; Pearce, A.J. Different current intensities of anodal transcranial direct current stimulation do not differentially modulate motor cortex plasticity. Neural Plast. 2013, 2013, 603502. [Google Scholar] [CrossRef] [Green Version]
  63. Saleem, G.T.; Ewen, J.B.; Crasta, J.E.; Slomine, B.S.; Cantarero, G.L.; Suskauer, S.J. Single-arm, open-label, dose escalation phase I study to evaluate the safety and feasibility of transcranial direct current stimulation with electroencephalography biomarkers in paediatric disorders of consciousness: A study protocol. BMJ Open 2019, 9, e029967. [Google Scholar] [CrossRef] [Green Version]
  64. Dmochowski, J.P.; Datta, A.; Huang, Y.; Richardson, J.D.; Bikson, M.; Fridriksson, J.; Parra, L.C. Targeted transcranial direct current stimulation for rehabilitation after stroke. Neuroimage 2013, 75, 12–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Dmochowski, J.P.; Koessler, L.; Norcia, A.M.; Bikson, M.; Parra, L.C. Optimal use of EEG recordings to target active brain areas with transcranial electrical stimulation. Neuroimage 2017, 157, 69–80. [Google Scholar] [CrossRef] [PubMed]
  66. Cancelli, A.; Cottone, C.; Tecchio, F.; Truong, D.Q.; Dmochowski, J.; Bikson, M. A simple method for EEG guided transcranial electrical stimulation without models. J. Neural Eng. 2016, 13, 036022. [Google Scholar] [CrossRef] [PubMed]
  67. Shahid, S.S.; Bikson, M.; Salman, H.; Wen, P.; Ahfock, T. The value and cost of complexity in predictive modelling: Role of tissue anisotropic conductivity and fibre tracts in neuromodulation. J. Neural Eng. 2014, 11, 036002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Saturnino, G.B.; Antunes, A.; Thielscher, A. On the importance of electrode parameters for shaping electric field patterns generated by tDCS. NeuroImage 2015, 120, 25–35. [Google Scholar] [CrossRef]
  69. Gomez-Tames, J.; Asai, A.; Hirata, A. Significant group-level hotspots found in deep brain regions during transcranial direct current stimulation (tDCS): A computational analysis of electric fields. Clin. Neurophysiol. 2020, 131, 755–765. [Google Scholar] [CrossRef]
  70. Minjoli, S.; Saturnino, G.B.; Blicher, J.U.; Stagg, C.J.; Siebner, H.R.; Antunes, A.; Thielscher, A. The impact of large structural brain changes in chronic stroke patients on the electric field caused by transcranial brain stimulation. NeuroImage Clin. 2017, 15, 106–117. [Google Scholar] [CrossRef] [PubMed]
  71. Habich, A.; Fehér, K.D.; Antonenko, D.; Boraxbekk, C.-J.; Flöel, A.; Nissen, C.; Siebner, H.R.; Thielscher, A.; Klöppel, S. Stimulating aged brains with transcranial direct current stimulation: Opportunities and challenges. Psychiatry Res. Neuroimaging 2020, 306, 111179. [Google Scholar] [CrossRef] [PubMed]
  72. Hsu, W.-Y.; Ku, Y.; Zanto, T.P.; Gazzaley, A. Effects of noninvasive brain stimulation on cognitive function in healthy aging and Alzheimer’s disease: A systematic review and meta-analysis. Neurobiol. Aging 2015, 36, 2348–2359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Brunoni, A.R.; Nitsche, M.A.; Bolognini, N.; Bikson, M.; Wagner, T.; Merabet, L.; Edwards, D.J.; Valero-Cabre, A.; Rotenberg, A.; Pascual-Leone, A.; et al. Clinical Research with Transcranial Direct Current Stimulation (tDCS): Challenges and Future Directions. Brain Stimul. 2012, 5, 175–195. [Google Scholar] [CrossRef] [Green Version]
  74. Datta, A. Inter-individual variation during transcranial direct current stimulation and normalization of dose using MRI-derived computational models. Front. Psychiatry 2012, 3, 91. [Google Scholar] [CrossRef] [Green Version]
  75. Rudroff, T.; Workman, C.D.; Fietsam, A.C.; Kamholz, J. Response variability in transcranial direct current stimulation: Why sex matters. Front. Psychiatry 2020, 11, 585. [Google Scholar] [CrossRef]
  76. Thomas, C.; Ghodratitoostani, I.; Delbem, A.C.B.; Ali, A.; Datta, A. Influence of gender-related differences in transcranial direct current stimulation: A Computational Study. In Proceedings of the 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Berlin, Germany, 23–27 July 2019; pp. 5196–5199. [Google Scholar]
  77. Russell, M.; Goodman, T.; Wang, Q.; Groshong, B.; Lyeth, B.G. Gender differences in current received during transcranial electrical stimulation. Front. Psychiatry 2014, 5, 104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Indahlastari, A.; Albizu, A.; O’Shea, A.; Forbes, M.A.; Nissim, N.R.; Kraft, J.N.; Evangelista, N.D.; Hausman, H.K.; Woods, A.J.; Initiative, A.D.N. Modeling Transcranial Electrical Stimulation in the Aging Brain. Brain Stimul. 2020, 13, 664–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Wallace, D.; Cooper, N.R.; Paulmann, S.; Fitzgerald, P.B.; Russo, R. Perceived comfort and blinding efficacy in randomised sham-controlled transcranial direct current stimulation (tDCS) trials at 2 mA in young and older healthy adults. PLoS ONE 2016, 11, e0149703. [Google Scholar]
  80. Chhatbar, P.Y.; Chen, R.; Deardorff, R.; Dellenbach, B.; Kautz, S.A.; George, M.S.; Feng, W. Safety and tolerability of transcranial direct current stimulation to stroke patients–A phase I current escalation study. Brain Stimul. 2017, 10, 553–559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Workman, C.D.; Kamholz, J.; Rudroff, T. The tolerability and efficacy of 4 mA transcranial direct current stimulation on leg muscle fatigability. Brain Sci. 2020, 10, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Uddin, L.Q.; Yeo, B.T.; Spreng, R.N. Towards a universal taxonomy of macro-scale functional human brain networks. Brain Topogr. 2019, 32, 926–942. [Google Scholar] [CrossRef] [PubMed]
  83. Thomas Yeo, B.T.; Krienen, F.M.; Sepulcre, J.; Sabuncu, M.R.; Lashkari, D.; Hollinshead, M.; Roffman, J.L.; Smoller, J.W.; Zöllei, L.; Polimeni, J.R.; et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. 2011, 106, 1125–1165. [Google Scholar] [CrossRef]
  84. Kashyap, R.; Ouyang, G.; Sommer, W.; Zhou, C. Neuroanatomic localization of priming effects for famous faces with latency-corrected event-related potentials. Brain Res. 2016, 1632, 58–72. [Google Scholar] [CrossRef] [PubMed]
  85. Kashyap, R.; Kong, R.; Bhattacharjee, S.; Li, J.; Zhou, J.; Yeo, B.T. Individual-specific fMRI-Subspaces improve functional connectivity prediction of behavior. NeuroImage 2019, 189, 804–812. [Google Scholar] [CrossRef] [PubMed]
  86. Kashyap, R.; Bhattacharjee, S.; Sommer, W.; Zhou, C. Repetition priming effects for famous faces through dynamic causal modelling of latency-corrected event-related brain potentials. Eur. J. Neurosci. 2019, 49, 1330–1347. [Google Scholar] [CrossRef] [PubMed]
  87. Kashyap, R.; Eng, G.K.; Bhattacharjee, S.; Gupta, B.; Ho, R.; Ho, C.S.; Zhang, M.; Mahendran, R.; Sim, K.; Chen, S.A. Individual-fMRI-approaches reveal cerebellum and visual communities to be functionally connected in obsessive compulsive disorder. Sci. Rep. 2021, 11, 1–15. [Google Scholar] [CrossRef]
  88. Kashyap, R.; Bhattacharjee, S.; Yeo, B.T.; Chen, S.A. Maximizing dissimilarity in resting state detects heterogeneous subtypes in healthy population associated with high substance use and problems in antisocial personality. Hum. Brain Mapp. 2020, 41, 1261–1273. [Google Scholar] [CrossRef] [PubMed]
  89. Csifcsák, G.; Boayue, N.M.; Puonti, O.; Thielscher, A.; Mittner, M. Effects of transcranial direct current stimulation for treating depression: A modeling study. J. Affect. Disord. 2018, 234, 164–173. [Google Scholar] [CrossRef] [Green Version]
  90. Siddiqi, S.H.; Schaper, F.L.W.V.J.; Horn, A.; Hsu, J.; Padmanabhan, J.L.; Brodtmann, A.; Cash, R.F.H.; Corbetta, M.; Choi, K.S.; Dougherty, D.D.; et al. Brain stimulation and brain lesions converge on common causal circuits in neuropsychiatric disease. Nat. Hum. Behav. 2021. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Illustration of the applied montage F3-RSO (shown in (A)) and output of i-SATA(MNI) for the MNI152 standard head image across the three current doses: (B) 1 mA; (C) 2 mA; and (D) 3 mA. The average current density at the target ROI (left middle frontal gyrus) is shown in the dark-gray colored bar. The DTDI index (~ 0.85) remains fairly constant across the three current doses, indicating a linear relationship between the injected current and induced electric field. For the standard head, tDCS users can choose any current dose, and their choice depends on the intensity of the desired stimulation in the target ROI. Note that the i-SATA(MNI) outputs shown in the figure highlight only the ROIs that received the highest current density (top 10%) amongst all other areas. A user can use the i-SATA(MNI) toolbox to obtain the average current density across all the brain regions.
Figure 1. Illustration of the applied montage F3-RSO (shown in (A)) and output of i-SATA(MNI) for the MNI152 standard head image across the three current doses: (B) 1 mA; (C) 2 mA; and (D) 3 mA. The average current density at the target ROI (left middle frontal gyrus) is shown in the dark-gray colored bar. The DTDI index (~ 0.85) remains fairly constant across the three current doses, indicating a linear relationship between the injected current and induced electric field. For the standard head, tDCS users can choose any current dose, and their choice depends on the intensity of the desired stimulation in the target ROI. Note that the i-SATA(MNI) outputs shown in the figure highlight only the ROIs that received the highest current density (top 10%) amongst all other areas. A user can use the i-SATA(MNI) toolbox to obtain the average current density across all the brain regions.
Jpm 11 00940 g001
Figure 2. Illustration of the variation in DTDI for three current doses (1 mA, 2 mA, and 3 mA) with F3-RSO montage applied over three individuals showing: (A) no change in DTDI with an increase in current dose; (B) increase in DTDI with an increase in current dose; and (C) decrease in DTDI with an increase in current dose. Subject 1, with no change in DTDI, is neutral to variation in current dose, and the dose can be adjusted based on the intensity of the stimulation desired at the target ROI. Subject 2, showing increase in DTDI, would receive adequate stimulation from a higher dose (above 1 mA), whereas subject 3, showing a decrease, will most likely benefit from the lower dose (1 mA).
Figure 2. Illustration of the variation in DTDI for three current doses (1 mA, 2 mA, and 3 mA) with F3-RSO montage applied over three individuals showing: (A) no change in DTDI with an increase in current dose; (B) increase in DTDI with an increase in current dose; and (C) decrease in DTDI with an increase in current dose. Subject 1, with no change in DTDI, is neutral to variation in current dose, and the dose can be adjusted based on the intensity of the stimulation desired at the target ROI. Subject 2, showing increase in DTDI, would receive adequate stimulation from a higher dose (above 1 mA), whereas subject 3, showing a decrease, will most likely benefit from the lower dose (1 mA).
Jpm 11 00940 g002
Figure 3. Illustration of the variation in DTDI at individual and group level with: (A) The individual variation of DTDI values (0 to 1) in the females (font in red ‘F’) and males (font in blue ‘M’) distributed equally across the three age groups ((i) young adults (20 ≤ x ˂ 40 years), (ii) mid adults (40 ≤ x ˂ 60 years), and (iii) older adults (60 ≤ x ˂ 80 years)) for the three current doses (1 mA, 2 mA, and 3 mA). The interindividual and intraindividual variation in DTDI clearly shows the current dose that could be appropriate for an individual to stimulate the target ROI after a montage has been fixed; (B) The variation of DTDI for both sexes across the three age groups using a three-way mixed ANOVA. The DTDI decreases with an increase in age. In mid and older adults, females show higher focality compared to males for the three current doses (1 mA, 2 mA, and 3 mA). In older adults only, the significant (p < 0.05, Bonferroni-corrected) difference between DTDI at 1 mA and 3 mA for both sexes conveys that higher current doses are required to appropriately stimulate the target ROI.
Figure 3. Illustration of the variation in DTDI at individual and group level with: (A) The individual variation of DTDI values (0 to 1) in the females (font in red ‘F’) and males (font in blue ‘M’) distributed equally across the three age groups ((i) young adults (20 ≤ x ˂ 40 years), (ii) mid adults (40 ≤ x ˂ 60 years), and (iii) older adults (60 ≤ x ˂ 80 years)) for the three current doses (1 mA, 2 mA, and 3 mA). The interindividual and intraindividual variation in DTDI clearly shows the current dose that could be appropriate for an individual to stimulate the target ROI after a montage has been fixed; (B) The variation of DTDI for both sexes across the three age groups using a three-way mixed ANOVA. The DTDI decreases with an increase in age. In mid and older adults, females show higher focality compared to males for the three current doses (1 mA, 2 mA, and 3 mA). In older adults only, the significant (p < 0.05, Bonferroni-corrected) difference between DTDI at 1 mA and 3 mA for both sexes conveys that higher current doses are required to appropriately stimulate the target ROI.
Jpm 11 00940 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kashyap, R.; Bhattacharjee, S.; Arumugam, R.; Bharath, R.D.; Udupa, K.; Oishi, K.; Desmond, J.E.; Chen, S.H.A.; Guan, C. Focality-Oriented Selection of Current Dose for Transcranial Direct Current Stimulation. J. Pers. Med. 2021, 11, 940. https://doi.org/10.3390/jpm11090940

AMA Style

Kashyap R, Bhattacharjee S, Arumugam R, Bharath RD, Udupa K, Oishi K, Desmond JE, Chen SHA, Guan C. Focality-Oriented Selection of Current Dose for Transcranial Direct Current Stimulation. Journal of Personalized Medicine. 2021; 11(9):940. https://doi.org/10.3390/jpm11090940

Chicago/Turabian Style

Kashyap, Rajan, Sagarika Bhattacharjee, Ramaswamy Arumugam, Rose Dawn Bharath, Kaviraja Udupa, Kenichi Oishi, John E. Desmond, S. H. Annabel Chen, and Cuntai Guan. 2021. "Focality-Oriented Selection of Current Dose for Transcranial Direct Current Stimulation" Journal of Personalized Medicine 11, no. 9: 940. https://doi.org/10.3390/jpm11090940

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop