This protocol describes the method of neuronavigated electrode placement for focal, transcranial direct current stimulation (tDCS) administered during functional magnetic resonance imaging (fMRI).
Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique that allows the modulation of the excitability and plasticity of the human brain. Focalized tDCS setups use specific electrode arrangements to constrain the current flow to circumscribed brain regions. However, the effectiveness of focalized tDCS can be compromised by electrode positioning errors on the scalp, resulting in significant reductions of the current dose reaching the target brain regions for tDCS. Electrode placement guided by neuronavigation based on the individual's head and brain anatomy derived from structural magnetic resonance imaging (MRI) data may be suited to improve positioning accuracy.
This protocol describes the method of neuronavigated electrode placement for a focalized tDCS setup, which is suitable for concurrent administration during functional MRI (fMRI). We also quantify the accuracy of electrode placement and investigate electrode drift in a concurrent tDCS-fMRI experiment. Critical steps involve the optimization of electrode positions based on current modeling that considers the individual's head and brain anatomy, the implementation of neuronavigated electrode placement on the scalp, and the administration of optimized and focal tDCS during fMRI.
The regional precision of electrode placement is quantified using the Euclidean norm (L2 Norm) to determine deviations of the actual from the intended electrode positions during a concurrent tDCS-fMRI study. Any potential displacement of electrodes (drift) during the experiment is investigated by comparing actual electrode positions before and after the fMRI acquisition. In addition, we directly compare the placement accuracy of neuronavigated tDCS to that achieved by a scalp-based targeting approach (a 10-20 Electroencephalography (EEG) system). These analyses demonstrate superior placement accuracy for neuronavigation compared to scalp-based electrode placement and negligible electrode drift across a 20 min scanning period.
Transcranial direct current stimulation (tDCS) is a noninvasive brain stimulation technique that allows the modification of cognition and physiological brain functions in experimental and clinical contexts1,2,3. Acute administration of tDCS can have transient changes in neuronal excitability, with the aftereffects lasting from minutes to hours after the stimulation4,5. The applied current does not induce action potentials but rather transiently shifts the resting membrane potential of the neuron toward either de- or hyperpolarization, resulting in either increased or decreased neuronal excitability at the macroscopic level using standard protocols4,5,6. Furthermore, regarding synaptic plasticity effects of tDCS, animal and human studies have shown that tDCS induces long-term potentiation and depression (LTP and LTD)-like processes4,5.
In the motor system, the modulation of motor evoked potentials (MEPs) allows direct assessment of neurophysiological effects of tDCS on local cortical excitability7. However, this approach cannot quantify the neural effects of tDCS on higher-order cognitive functions supported by large-scale functional brain networks8. The effects on brain networks can be investigated by combining tDCS with modern functional imaging techniques9,10. Among those, functional magnetic resonance imaging (fMRI) has become the most frequently used approach because it provides excellent spatial and sufficient temporal resolution to reveal the neural mechanisms by which tDCS affects local brain activity at the stimulation site and large-scale neural networks11,12,13,14.
So far, combined fMRI-tDCS studies have mainly employed so-called conventional tDCS setups, which use relatively large rubber electrodes between 25 and 35 cm2 (5 x 5 cm2 and 5 x 7 cm2) inserted into saline-soaked sponge pockets15,16. These setups project the current between two electrodes that are typically attached over (a) a target brain region for tDCS and (b) a return electrode over non-target brain regions or extracranial areas (e.g., the shoulder). This results in widespread current flow across the brain, affecting regions other than the target region, thereby complicating causal assumptions and interpretations about the neural origin of tDCS effects17.
More precise spatial targeting can be achieved by focalized tDCS18. These setups employ arrays of smaller electrodes arranged in close proximity to each other or by using a ring-shaped cathode placed around a center anode to constrain current flow to the target region18,19. Computer simulations of electric current flow suggest that focalized tDCS can result in higher spatial precision of the current flow to the target region than conventional montages20. Moreover, behavioral studies have demonstrated regional and task-specific behavioral modulation using focalized setups19,21,22. However, only a few studies have used focalized tDCS during fMRI. These studies have been able to establish the feasibility of this approach and provided the first evidence for region-specific neural modulation19,23.
However, because of regionally precise current delivery, focalized tDCS setups may be more sensitive to electrode positioning errors on the scalp than conventional montages. For example, Seo et al. demonstrated that positioning errors of 5 mm in a focalized motor cortex setup reduced peak somatic polarization in the hand knob by up to 87%24. Moreover, a recent computational modeling study demonstrated that electrode displacement from intended positions for focal compared to conventional setups resulted in significant current dose reductions in the target regions for tDCS, ranging from 26% to 43%25. Therefore, it was concluded that future studies should routinely include appropriate methods for the improvement of electrode positioning and the verification of electrode positioning before and after fMRI5.
In the present study, we describe the method of neuronavigated electrode positioning for a novel fMRI-compatible focal 3 x 1 tDCS setup (i.e., three individual cathodes that are arranged in a circle around a single center anode), which is currently being used in a collaborative research consortium funded by the German Science Foundation (DFG Research Unit 5429, https://www.memoslap.de). The consortium investigates behavioral and neural effects of focalized tDCS on learning and memory and predictors of stimulation response across four functional domains (i.e., visual-spatial, language, motor, and executive functions). Structural T1- and T2-weighted MRI data of study participants are acquired during a baseline scan. These data are used for individualized current flow simulations26 to determine scalp positions of electrodes that maximize the current flow to the target region in individual study participants. As an example, this protocol will describe neuronavigated targeting of individually determined electrode positions centered over the right dorsolateral prefrontal cortex (rDLPFC) in one participant.
The representative results section is based on structural imaging data acquired before and after concurrent tDCS-fMRI in three subprojects of the Research Unit. These studies targeted the right occipitotemporal cortex (rOTC), left temporo-parietal cortex (lTPC), and rDLPFC. Data were acquired at the Department of Neurology of the University Medicine Greifswald. Using these data, we aimed to achieve two main objectives: (1) To quantify the spatial precision of neuronavigated electrode placement by comparing "intended" versus empirically determined "actual" electrode positions25, and (2) investigate the degree of electrode displacement over the course of the fMRI sessions (i.e., electrode drift). These factors are crucial for improving the accuracy and reliability of tDCS effects in concurrent tDCS-fMRI studies27. In addition, the targeting accuracy of neuronavigated tDCS is compared to that of a scalp-based approach using data from a previous tDCS-fMRI study of our group25.
All experimental procedures presented in this protocol have been reviewed and approved by the ethics committee of the University Medicine Greifswald. All participants provided informed consent prior to study inclusion and granted permission for their data to be published anonymously.
1. Screening of contraindications and general considerations
2. Baseline MRI scan and individualized current modeling
3. Neuronavigation
4. tDCS-fMRI
Data from 43 healthy young participants (20 men/23 women, aged 24.74 ± 5.50 years) were included. The participants completed up to four fMRI sessions. Neuronavigated placement of electrodes was conducted prior to each fMRI session. In total, 338 datasets representing the positions of the center anodes before and after fMRI were included in data analyses.
To determine the intended positions of the electrodes, individualized current modeling was performed using structural MRI data acquired in a separate MRI session (prior to the fMRI sessions). Current modeling aimed to optimize the scalp positions of the 3 x 1 setups to maximize current flow to the rDLPFC (N = 19, participants/N = 137 data points), rOTC (N = 12/104) or lTPJ (N = 12/97). Actual positions of the center anodes of the 3 x 1 setups were empirically determined using structural MRI data (PETRA images) acquired before and after each fMRI session (Note: Coordinates for targeting were provided in subject space. Hence, PETRA scans of individual participants were co-registered with their own T1-weighted image acquired during the baseline scan, on which field simulations were based). Because identification of electrode positions in 3D is prone to errors, a modified 2D "pancake" view of the reconstructed scalp with electrodes was used (see Figure 537,38) to determine six edge points of the center anodes. Subsequently, the center of mass was calculated to obtain the center coordinates of the anodes for each fMRI session and time-point (i.e., before or after tDCS). Subsequently, a reverse 2D->3D transformation was applied to obtain electrode coordinates in 3D space.
The spatial precision of the neuronavigated electrode placement was investigated by comparing "intended" versus "actual" electrode positions25 and is reported as the L2 Norm (Euclidean norm) distance39. The L2 Norm is a standard vector norm calculated as the square root of the sum of the squares of the vector components. The difference between intended and actual coordinates are calculated by the difference of all vector components, namely each dimension of the 3D space (x/y/z). Electrode displacement over the course of the experiment (i.e., drift) was investigated by comparing actual electrode positions determined before and after fMRI.
In addition, we compared targeting accuracy for neuronavigated electrode placement with data from a previous study of our group that used a scalp-based targeting approach (10-20 EEG system25). This study comprised data from 60 participants (31 men/29 women, aged 22.17 ± 2.33 years), two target regions (N = 30/30, motor cortex/inferior frontal gyrus) and two fMRI sessions for a total of 240 data points. Data and code used for all analyses are publicly available at https://github.com/LawsOfForm/VerFlu_Simulation_Niemann_2023.
Linear mixed models were used to analyze (a) deviation of actual electrode positions from intended locations and (b) drift across the fMRI session with the following formula and the L² Norm as the dependent variable:
yij = β0 + POSITIONx1ij + REGIONx2ij + POSITION * REGIONx3ij + µ0j + εij
yij is the L2 Norm distance between intended and actual electrode positions at POSITION i of subject j. β0 is the fixed effect intercept, X1 is the value for POSITION (1: intended, 2: actual pre-fMRI; 3: post-fMRI) at position i in subject j. X2 is the value for the REGION of interest (1: lTPC, 2: rDLPFC, 3:rOTC) at region i in subject j. X3 is the value for the interaction for POSITION and REGION for position and region i and subject j, u0 is the random intercept for subject j (µ0j ~ N(0, τ002)); εij are the residuals for timepoint i in subject j (εij ~ N(0,σ2)).
Comparison of targeting accuracy for neuronavigated and scalp-based electrode placement was done using a t-test based on the L2 Norm differences between the intended and actual positions in the respective groups.
Neuronavigated placement accuracy and displacement across the experiment
The reference model (intended placement of center anodes; region: lTPC) showed significant deviations of electrode positions from intended positions at both pre-fMRI and post-fMRI assessments (8.71 mm, 95%-CI: [6.57-10.84], p < .001; 8.51 mm, 95%-CI: [6.35-10.66], p < .001) Supplemental Table S1. Post-hoc analysis for estimated marginal means (EMMs) were calculated for all main and simple effects of interest (see Supplemental Table S2, Supplemental Table S3, Figure 6, Figure 7, Supplemental Figure S3, and Supplemental Figure S4). The predicted average placement error for all pre-fMRI sessions was 8.54 mm 95%-CI: [7.66-9.42] and significantly different from the intended positions (t (296) = -8.54, p < 0.0001). This was also the case for post-fMRI coordinates (8.61 mm, 95%-CI: [7.73-9.50]; t (296) = -10.97, p < 0.0001, see Figure 6 and Supplemental Table S2). Because placement errors may vary by region25,40,41, misplacement was also calculated between target regions for each position. However, no region-specific differences between rOTC, lTPC, and rDLPFC were found (Supplemental Table S3, Figure 7, and Supplemental Figure S4).
Displacement of center anodes (i.e., drift), assessed by comparison of actual electrode positions between pre- and post-fMRI was negligible (-0.08 ± 0.39 mm; t (291) =-0.20, p = 0.98). For details, see Figure 6 and Supplemental Table S2.
Comparison of electrode placement errors between neuronavigated and scalp-based placement approaches
To investigate the relevance of neuronavigated placement, L² Norm differences between intended and actual positions based on neuronavigated electrode placement were retrospectively compared to data obtained from a previous study that used a scalp-based targeting approach25. In this study, the L² Norm difference between the intended and actual positions of the anode was 12.33 ± 5.67 mm. Direct comparison between neuronavigated (8.83 ± 4.68) and scalp-based approaches revealed a significant decrease of the L² Norm distance by 3.50 mm in favor of neuronavigated electrode placement (t (397) =4.91, p < 0.001).
Figure 1: A simplified description of the sequence of steps representing the method for individualized current modeling. The method aims to determine optimized scalp positions of electrodes over the target region (right DLPFC) using a focal 3 x 1 tDCS setup. (A) Anatomical target region in gray. (B) Optimized electrode position illustrated on the reconstructed head of the participant (anode = black, cathode 1 = red, cathode 2 = green, cathode 3 = blue. (C) Simulated electric field based on individualized current modeling for this montage (in V/m). (D) Target coordinates for all intended electrode positions used during neuronavigated targeting (dimensions: x/y/z). Abbreviations: DLPFC = right dorsolateral prefrontal cortex; tDCS = transcranial direct current stimulation. Please click here to view a larger version of this figure.
Figure 2: Neuronavigation equipment. 1) Subject tracker, 2) Screwdriver, 3) Hexagonal rod, 4) Goggles, 5) Pointer. NOTE: The neuronavigation software communicates with the position sensor (the tracking camera) to obtain the location of the subject tracker (mounted on the goggles). It uses a registration matrix (obtained by identifying landmarks on the subject's head through the pointer) to align the subject's head position from the real world to the participant's image space (based on the T1-weighted image of the participant). An overview of the overall setup is provided in Supplemental Figure S1. Please click here to view a larger version of this figure.
Figure 3: Target region validation. Illustrates the procedure of validating stimulation target regions via neuronavigation. 1) The Subject tracker assembled on the left side of goggles. 2) The Pointer is used to locate the intended position of the anode (rDLPFC). 3) The green crosshair and the red dot (representing the tip of the pointer) on the computer screen. 3a) The index number is green if the pointer deviates from the crosshair by less than 3 mm; 3b) the index number is orange if the pointer deviates between 3 and 5 mm; and 3c) it is red if the pointer deviates more than 5 mm. Abbreviation: rDLPFC = right Dorsolateral Prefrontal Cortex. Please click here to view a larger version of this figure.
Figure 4: Selected components of multi-channel DC-stimulator and *custom components. (A) 1) Adapter that plugs into the DC stimulator outlet, 2) Outer box LAN cable, 3) *Outer box (filter box) that is positioned inside the waveguide of the scanner during tDCS-fMRI. The clamp is attached to the waveguide to connect the filter box with the floating mass of the shielded MRI chamber to reduce potential artifacts. 4) Inner box cable, 5) inner box with plugs for electrode cables, 6) *fuse adapter that combines three individual channels to one for the center anode, 7) *fused electrode cables and electrode of the center anode (10 mm radius) made of flexible and electrically conductive silicone, 8) individual cathode and cable, 9) example of electrode inserted into 3D printed fill aid, 10) example of electrode placed in 3D spacer, 11) measuring tape, 12) EEG cap without plastic inserts to secure electrodes on the scalp. For the complete setup and standard components not used in this protocol, readers are referred to the website of the manufacturer (link in Table of Materials). (B) 3D model of the electrode fill aid. (C) 3D model of the spacer. Abbreviations: tDCS-fMRI = transcranial direct current stimulation-functional magnetic resonance imaging; EEG = electroencephalography. Please click here to view a larger version of this figure.
Figure 5: Pancake view of scalp and electrodes based on the PETRA image. Colors represent different sizes of scalp flattened mask (red = 2 mm; green = 4 mm, blue = 6 mm of the inflated scalp for depth visualization of the head). Orientations of the transformed image are marked in white letters. White dots represent manually selected edge points of the center anode placed over the right dorsolateral prefrontal cortex. The center of mass was calculated as the midpoint coordinate of the edge points prior to retransformation into 3D space. Abbreviations: PETRA = Pointwise Encoding Time reduction with Radial Acquisition; L = left; R = right; P = posterior; A = anterior. Please click here to view a larger version of this figure.
Figure 6: Estimated Marginal Means of the L2 Norm (Euclidean norm) distance between intended and actual electrode positions of the anodal electrode across all sessions and regions (rOTC, lTPC, rDLPFC). Original data points in turquoise. EMMs 95% confidence intervals in grey. Abbreviations: EMMs = Estimated Marginal Means; rOTC = Right Occipital Temporal Cortex; lTPC = left Temporal Parietal Cortex; rDLPFC = right Dorsolateral Prefrontal Cortex; fMRI = functional magnetic resonance imaging. Please click here to view a larger version of this figure.
Figure 7: Estimated Marginal Means of L2 Norm (Euclidean norm) distance between intended and actual electrode positions of the anodal electrode over all sessions between each region (rOTC, lTPC, rDLPFC) for each position (intended, pre/post fMRI). Original data points in turquoise. EMMs 95% confidence intervals in grey. Abbreviations: EMMs = Estimated Marginal Means; rOTC = Right Occipital Temporal Cortex; lTPC = left Temporal Parietal Cortex; rDLPFC = right Dorsolateral Prefrontal Cortex; fMRI = functional magnetic resonance imaging. Please click here to view a larger version of this figure.
Supplemental Figure S1: The neuronavigation system overview. 1) Control computer of the neuronavigation system with the neuronavigation software installed; 2) Polaris camera; 3) Foot pedal; 4) Goggles with subject tracker (see Figure 2 for details); 5) Pointer. Please click here to download this File.
Supplemental Figure S2: Reconstructed skin with landmarks in yellow. 1) Nasion, 2) left nostril; 3) right nostril; 4) left preauricular pit; 5) right preauricular pit (see protocol step 3.1.7). Abbreviations: LPA = left preauricular pit; RPA = right preauricular pit. Please click here to download this File.
Supplemental Figure S3: Estimated Marginal Means of L² Norm (Euclidean norm) distance between intended and actual electrode positions of the anodal electrode across all sessions between positions (intended, pre/post fMRI) for each region (rOTC, lTPC, rDLPFC). Original data points in turquoise. EMMs 95% confidence intervals in grey. Abbreviations: EMMs = Estimated Marginal Means; rOTC = Right Occipital Temporal Cortex; lTPC = left Temporal Parietal Cortex; rDLPFC = right Dorsolateral Prefrontal Cortex; fMRI = functional magnetic resonance imaging. Please click here to download this File.
Supplemental Figure S4: Estimated Marginal Means of L² Norm (Euclidean norm) distance between intended and actual electrode positions of the anodal electrode across all sessions and positions (intended, pre/post fMRI) for each region of interest (rOTC, lTPC, rDLPFC). Original data points in turquoise. EMMs 95% confidence intervals in grey. Right Occipital Temporal Cortex (rOTC), left Temporal Parietal Cortex (lTPC), right Dorsolateral Prefrontal Cortex (rDLPFC). Please click here to download this File.
Supplemental Table S1: Results of linear mixed model analysis for L² Norm (Euclidian norm) distance between intended and actual coordinates. Regression coefficients of linear mixed models (random intercept models) and two-sided p-values are reported. Abbreviations: CI = confidence interval; rOTC = Right Occipital Temporal Cortex; lTPC = left Temporal Parietal Cortex; rDLPFC = right Dorsolateral Prefrontal Cortex; df = degrees of freedom; σ2 = pooled residual variance; τ = pooled variance explained by subjects; ICC = Intraclass correlation coefficient; N = number of subjects. Please click here to download this File.
Supplemental Table S2: Post-hoc EMMs for Linear Mixed Model analysis for L² Norm (Euclidian norm) differences). Main and simple contrasts for main effect of position (intended, pre/post fMRI, see also Supplemental Figure S2) and simple effects of position for each region (lTPC, rDLPFC, rOTC, see also Figure 6). For simple effects, only significant effects are shown. Abbreviations: EMMs = Estimated Marginal Means; SE = Standard Error; df = degrees of freedom; rOTC = right Occipital Temporal Cortex; lTPC = left Temporal Parietal Cortex; rDLPFC = right Dorsolateral Prefrontal Cortex. Please click here to download this File.
Supplemental Table S3: Post-hoc EMMs for Linear Mixed Model analysis of L² Norm (Euclidian norm) difference. Main and simple contrasts for main effect of region (lTPC, rDLPFC, rOTC, see also Supplemental Figure S3) and simple effects of region for each position (intended, pre/post fMRI, see also Figure 7). Abbreviations: EMMs = Estimated Marginal Means; SE = Standard Error; df = degrees of freedom; rOTC = right Occipital Temporal Cortex; lTPC = left Temporal Parietal Cortex; rDLPFC = right Dorsolateral Prefrontal Cortex. Please click here to download this File.
Critical steps, potential modifications, and troubleshooting of the method
Accurate positioning of electrodes is a crucial technical factor in tDCS experiments, and deviations from intended scalp positions or electrode drift can affect current flow to the intended target brain regions42,43. This is particularly relevant for focalized tDCS, as the regional specificity of the administered current makes these setups particularly susceptible to the effects of positioning errors5,25. In the current protocol, we described a method of neuronavigated, focalized tDCS administered during fMRI, which aims to enhance the placement accuracy of individually determined scalp positions. Critical steps include (1) the use of electrical current flow modeling to optimize scalp position of the 3 x 1 setup for each study participant, (2) identification of optimized scalp positions by neuronavigation, (3) positioning and securing of electrodes on the scalp using standardized procedures, and (4) investigation of placement accuracy and electrode drift before and after concurrent tDCS-fMRI.
The first step in our protocol describes the optimization of electrode positions for an individual study participant of our ongoing Research Unit (RU) by means of individualized current modeling (code to rebuild the optimization can be found here https://github.com/memoslap/SimNIBS_Optimization). This approach is highly specific to the RU projects and involved two preparatory steps not reported here: (1) Definition of anatomical target regions for eight empirical RU projects (www.memoslap.de) based on fMRI and/or tDCS studies that used the same/similar motor or cognitive tasks; (2) equalizing the average current dose across projects by adjusting the spacing between the center anode and the surrounding cathodes. This involved a computational modeling study and structural imaging data available to our group to determine the project-specific spacing of electrodes. These project-specific montages were then used to optimize scalp positions for individual participants in the RU projects to achieve maximal current intensity in the respective target regions (in this protocol: right DLPFC). Notably, there are many possible ways to determine the desired scalp positions, and those depend on the goals of the study (i.e., regional or dose control, anatomical or functional targeting).
Similarly, there are a number of different commercially available neuronavigation devices and DC stimulators that can be used for neuronavigated placement and conducting a tDCS-fMRI study. Consequently, modifications of the procedures outlined in this protocol will be required if equipment from different manufacturers is used or tDCS-fMRI is conducted with a different scanner type. Nonetheless, we believe that the protocol can be adapted in a flexible way to meet study-specific requirements. For example, for 3 x 1 setups, we recommend using a 3D printed spacer to ensure that the desired distance between the electrodes is achieved and maintained during placement and scanning. Standardization of the amount of conductive paste can be achieved using an electrode fill aid as described in this protocol. To allow for customization, the 3D models used for printing have been made publicly available (see protocol step 4.1.3.2).
The complexity of intrascanner tDCS fMRI experiments is another crucial consideration, which requires participants to move into electrodes attached to their heads. These studies carry a significant risk of subsequent electrode displacement, even with precise positioning5. Therefore, we recommend routinely implementing appropriate methods to minimize electrode displacement and drift (e.g., using standard MRI-compatible inflatable cushions to avoid displacement of electrodes inside the MRI scanner, EEG caps or contoured straps to keep electrodes in place, and other measures detailed in the protocol above).
Finally, compared to behavioral tDCS studies, administering tDCS during fMRI allows for the verification of the exact positions of electrodes on the scalp (based on 3D structural MRI images) and the quantification of potential deviations from the intended positions. However, the identification of electrode positions in 3D is prone to errors. Therefore, we suggest a 2D transformation of the scalp and electrodes to facilitate the verification of actual electrode positions. However, future studies are required to validate this novel approach and compare outcomes to data extraction from 3D images.
Significance of method relative to existing other methods
Scalp-based targeting approaches (e.g., the 10-20 EEG system44) have been used in the majority of tDCS research to date. This involves manual or automated identification of anatomical landmarks and additional measurements (e.g., head circumference, calculation of intersections between landmarks) to determine the intended scalp positions of the electrodes. While this approach accounts for the head size of individual participants to some extent, it neglects other characteristics of brain and skull morphology, resulting in a loss of accuracy5,41,45. The enhanced precision that can be achieved with MRI-guided neuronavigation is illustrated in the representative results section of this protocol. Here, we quantified the precision of neuronavigated targeting using data from an ongoing multicenter tDCS-fMRI study and retrospectively compared placement accuracy of neuronavigated targeting with a scalp-based approach. In addition, we investigated the degree of electrode displacement during fMRI. The main outcomes of the empirical part of this protocol were that (1) neuronavigated targeting is superior to scalp-based electrode placement, with approximately 40% higher precision based on the L2 Norm; and (2) there was minimal drift of electrodes across fMRI sessions, highlighting that the implemented technical and procedural steps to keep electrodes in place (i.e., cap, spacer, inflatable cushions) were successful.
Limitations and future applications of the method
There are several limiting factors for using the approach described in this protocol. For example, optimization of the intended electrode positions by computational modeling, neuronavigated targeting, and verification of electrodes require access to structural MRI data, specific technical expertise, and costly equipment that may not be available. Hence, this approach is limited to specialized research centers performing studies that require high spatial precision of current delivery (e.g., experimental studies in neurotypical populations investigating causal brain-behavior relationships or computational modeling studies investigating dose-response relationships). In this context, it is worth mentioning that the validity of MRI-based current modeling may be decreased in patients with head surgery or progressive lesions, which needs to be taken into consideration when applying the methods described in this protocol to clinical populations.
In addition, the optimization of electrode positions by computational modeling and neuronavigated placement in this protocol aims to maximize anatomical precision. As mentioned above, other ways of individualized targeting are possible (including invasive approaches limited to specific clinical populations46), and anatomical precision does not necessarily imply functional relevance of the stimulated area. The latter would require a localizer task (either on the individual participant or group level), acquired prior to tDCS-fMRI, to identify functionally relevant brain regions for neuronavigated targeting. In addition, future research needs to determine the functional relevance of electrode displacement on behavioral and neural functions, including a direct comparison of effects between focal and conventional montages.
It should be noted that computer simulations rely on the anatomical precision of the (semi-)automated tissue segmentations acquired from MRI and estimate the electric field based on assumptions about the electrical conductivity of various tissue classes5,47. Therefore, to increase the accuracy of computational models, it is essential to validate these assumptions and evaluate the segmentation accuracy. In this regard, new methods to study brain tissue's tDCS-induced current flow conductivity during simultaneous tDCS-MRI measurements include Magnetic Resonance Current Density Imaging (MRCDI) and Magnetic Resonance Electrical Impedance Tomography (MREIT)5,48,49. These approaches have the potential for validating tDCS electric field simulations and optimizing individualized current dose calculations. Combining such optimized computational modeling approaches with neuronavigated tDCS will likely improve the precision and effectiveness of tDCS in experimental and clinical contexts in the near future.
The authors have nothing to disclose.
This research was funded by the German Research Foundation (project grants: FL 379/26-1; ME 3161/3-1; CRC INST 276/741-2 and 292/155-1, Research Unit 5429/1 (467143400), FL 379/34-1, FL 379/35-1, Fl 379/37-1, Fl 379/22-1, Fl 379/26-1, ME 3161/5-1, ME 3161/6-1, AN 1103/5-1, TH 1330/6-1, TH 1330/7-1). AT was supported by the Lundbeck Foundation (grant R313-2019-622). We thank Sophie Dabelstein and Kira Hering for their help with data extraction.
Brainsight neuronavigation system | Brainsight; Rogue Research Inc., Montréal, Canada | ||
CR-5 Pro high temp 3D printer | CREALITY, Shenzhen, China | ||
DC-STIMULATOR MC | NeuroConn GmbH, Ilmenau, Germany | https://www.neurocaregroup.com/technology/dc-stimulator-mc | |
EMLA Cream 5% | Aspen, Dublin, Ireland | ||
MAGNETOM Vida 3T, syngo_MR_XA50 software | Siemens Healthineers AG, Forchheim, Germany | ||
Polaris camera | Polaris Vicra; Northern Digital Inc., Waterloo, Canada | ||
Ten20 conductive EEG paste | Weaver and Company, Aurora, USA | ||
TPU 3D printer filament | SUNLU International, Hong-Kong, China | ||
Example of alternatives | |||
Ingenia 3.0T (MR-scanner) | Phillips, Amsterdam, Netherlands | ||
Localite TMS Navigator (Neuronavigation equipment) | Localite, Bonn, Germany | ||
Neural Navigator (Neuronavigation equipment) | Soterix, New Jersey, USA | ||
PEBA 3D printer filament | Kimya, Nantes, France | ||
PLA 3D printer filament | Filamentworld, Neu-Ulm, Deutschland | ||
StarStim (Stimulator) | Neuroelectrics, Barcelona, Spain |
.