Neuroanatomical Considerations for Optimizing Thalamic Deep Brain Stimulation in Tourette Syndrome ================================================================================================== * Takashi Morishita * Yuki Sakai * Hitoshi Iida * Saki Yoshimura * Atsushi Ishii * Shinsuke Fujioka * Saori C. Tanaka * Tooru Inoue ## Abstract **Background** Deep brain stimulation (DBS) of the centromedian (CM) thalamic nucleus has reportedly been used to treat severe Tourette syndrome (TS) with promising outcomes; however, it remains unclear how DBS electrode position and stimulation parameters modulate the specific area and related networks. We aimed to evaluate the relationships between the anatomical location of stimulation fields and clinical responses including therapeutic and side effects. **Methods** We collected data from eight TS patients treated with DBS. We evaluated the clinical outcomes using Yale Global Tic Severity Scale (YGTSS), Yale-Brown Obsessive Compulsive Scale (Y-BOCS), and Hamilton Depression Rating Scale (HAM-D). The DBS lead location was evaluated in the normalised brain space using a 3-D atlas. The volume of tissue activated (VTA) was determined, and the associated normative connective analyses were performed to link the stimulation field with the therapeutic and side effects. **Results** The mean follow-up period was 10.9 ± 3.9 months. All clinical scale showed significant significant improvement. While the VTA associated with therapeutic effects covers the CM and ventrolateral nuclei and showed association with motor networks, those associated with paraesthesia and dizziness were associated with stimulation of the ventralis caudalis and red nucleus, respectively. Depressed mood was associated with the spread of stimulation current to the mediodorsal nucleus and showed association with limbic networks. **Conclusion** Our study addresses the importance of accurate implantation of DBS electrodes for obtaining standardised clinical outcomes and suggests that meticulous programming with careful monitoring of clinical symptoms may improve outcomes. Key Words * Tourette syndrome * deep brain stimulation * thalamus * adverse events * normative connectome ## Introduction Tourette syndrome (TS) is a neuropsychiatric disorder characterised by involuntary tic movements and psychiatric comorbidities. Although symptoms are subtle and subside spontaneously in childhood, some patients affected by this syndrome suffer from severe, debilitating involuntary movements. Common treatment options for TS include medical and behavioural therapies, and deep brain stimulation (DBS) may be indicated for the treatment of severe TS. Since the publication of the first report on the successful application of DBS therapy for TS [1], reports on its clinical outcomes have been published [2, 3]. A previous study concerning the treatment of TS using DBS has reported mixed results although overall clinical outcomes were favourable [2-4]. Additionally, a recent long-term follow-up study showed that a loss of benefit occurs after several years of continuous thalamic stimulation [5]. Several targets for DBS therapy including the centromedian (CM) thalamic nucleus, globus pallidus interna (GPi), and anterior limb of internal capsule have been proposed for the treatment of TS using DBS, and stimulation of these targets is reportedly equally effective [2, 3]. However, response to the therapy varies among patients. The heterogeneity of the therapeutic effects is potentially associated with the employed patient selection criteria, surgical technique, and DBS programme. A potential factor associated with DBS failure is lead misplacement [6], and a recent study reported that the incidence of DBS lead misplacement in Parkinson’s disease cohort was higher than 15% [7]. Even though the incidence of lead misplacement in the TS cohort has not been reported, we suspect its occurrence in some patients with TS refractory to DBS in the reported cohorts. Because the visualisation of each thalamic nucleus for stereotactic targeting has been challenging and because neurophysiological microelectrode recordings (MER) have not been established for the CM nucleus, it is difficult to determine the occurrence of lead misplacement in a patient with TS showing poor DBS outcome. Additionally, detailed subdivisions of the thalamic structures associated with favourable stimulation outcome have not been identified as in the case of the GPi or subthalamic nucleus, which are stimulated for the treatment of Parkinson’s disease although a multi-centre analysis reported preferred stimulation points [8]. Clinicians should be aware that a subtle difference in the lead position and electrically stimulated areas results in the activation of different neuronal networks; however, it remains unclear how DBS electrode position and stimulation parameters modulate the specific area and related networks. The position of DBS electrodes is conventionally reported using Cartesian coordinates relative to the midcommissural point, but recent studies have addressed the importance of considering anatomical variations while determining the optimal DBS electrode position [9]. To analyse patient-specific electrode positions, various authors have used a common population-based standard Montreal Neurological Institute (MNI) space [8, 10], and recent studies have reported the utility of three-dimensional anatomical atlases and structural connectome [11, 12]. Using these neuroimaging techniques, we aimed to evaluate the relationships between the anatomical location of implanted DBS leads and clinical responses. We also evaluated stimulation-induced side effects as they interfere with optimising stimulation settings for symptoms suppression. ## Methods ### Study Design We prospectively recorded the clinical course of the patients and reviewed the charts to obtain detailed information. To interpret the DBS lead positions associated with clinical benefits and stimulation-induced adverse events, we evaluated the DBS lead position in the postoperative images. This study was approved by our institutional review board (IRB) named Fukuoka University-Medical Ethics Review Board (IRB approval number: U02-02-001), and informed consent was obtained from the participants. This study was conducted in accordance with the Declaration of Helsinki. We included eight patients who underwent DBS surgery at our department, and all patients completed at least a six-month follow-up. DBS therapy was indicated for patients with severe, medication-refractory symptoms who were older than 12 years of age. One patient (case 7) underwent DBS surgery because of severe cervical tic with the risk of spinal cord injury despite relatively low clinical scoring owing to the absence of phonic tics and other comorbidities. All patients underwent multidisciplinary evaluation performed by a team consisting of a neurologist, psychiatrist, paediatrician, and neurosurgeon, as recommended by a recent guidelines paper [13]. The demographics of patients are summarised in Table 1. View this table: [Table 1.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/T1) Table 1. Patient Demographics. **Childhood: age 3-12; Adolescence: age 13-19, Young adult: age 20-39** **DBS = deep brain stimulation; SD = standard deviation** We evaluated the clinical outcomes using Yale Global Tic Severity Scale (YGTSS) [14], Yale-Brown Obsessive Compulsive Scale (Y-BOCS) [15], and Hamilton Depression Rating Scale (HAM-D) [16] (Table 2). These evaluations were performed prior to surgery and six months and one year after surgery. Additionally, we recorded the adverse events associated with surgery and neurostimulation. Mood changes were recorded in the medical records on the basis of self-reported complaints of the patients. View this table: [Table 2.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/T2) Table 2. Clinical outcomes. Abbreviations. F/U = follow up; HAM-D = Hamilton Depression Rating Scale; N/A = not applicable; Y-BOCS = Yale-Brown Obsessive Compulsive Scale; YGTSS = Yale Global Tic Severity Scale. ### Surgical Procedure Because we have reported our DBS procedures previously [17, 18], we describe our method briefly herein. Magnetic resonance imaging (MRI) was performed on all the patients for stereotactic planning using an MRI scanner (Ingenia 1.5T, Philips, Netherland), and the MRI sequences included volumetric T1-weighted imaging (T1WI) with contrast, volumetric fast grey matter acquisition T1 inversion recovery (FGATIR), and diffusion tensor imaging. Using commercialised software (iPlan stereotaxy, Brainlab, Germany), we identified the target structure and planned the trajectory of DBS lead implantation. Above mentioned MRI sequences were automatically fused, and the anterior commissure (AC), posterior commissure (PC), and the midline plane to anchor the Cartesian coordinate system were determined. The tentative target was initially set at 5 mm lateral and 4 mm posterior to the midcommissural point on the AC-PC plane. The trajectory was planned such that the DBS lead would not injure the blood vessels, lateral ventricle, and sulci during T1WI with contrast, and the target (presumed lead tip position at the ventral surface of the CM nucleus) was modified according to the location of the CM nucleus in each case. The CM nucleus was identified as a relatively high-intensity area on FGATIR images, as reported previously [18]. The target was meticulously checked using the mammillothalamic tract and the red nucleus (RN) as landmarks (Supplementary Figure 1). In this study, all patients underwent simultaneous bilateral implantation of DBS leads (model 3387, Medtronic, Minnesota) and implantation of an implantable pulse generator (IPG) (Activa RC, Medtronic, Minnesota) on the same day. The first four patients underwent DBS lead implantation under local anaesthesia; however, three of them could not maintain the fixed head position during the awake procedure. Thus, we used general anaesthesia throughout the surgical procedure for all system implantations, including the stereotactic frame fixation for the remaining four cases (cases 5–8). ### DBS Programming Electrical stimulation was delivered immediately following DBS system implantation in the first six patients and two weeks after DBS lead implantation in the remaining two patients (cases 7 and 8). The initial stimulation was delivered using the monopolar setting, which activates the second most ventral contacts (contacts 1 and 9) bilaterally, with the following parameters: amplitude, 2.0 volts or milliamperes; pulse width (PW), 60 microseconds; and frequency, 130 Hz. One month after surgery, we measured the threshold levels of stimulation-induced side effects in each contact at a fixed PW of 60 microseconds and frequency of 130 Hz. We carefully interviewed the patients at each programming session, and we increased the stimulation intensity within the therapeutic window such that acute side effects could be avoided. At each programming visit, new programming parameters were saved as a new group of settings so that patients can go back to the previous setting in case they experience the chronic stimulation-induced side effects such as mood changes. If the clinical response was insufficient with maximum intensity at one contact, we activated another contact to apply interleaving stimulations using multiple contacts. The frequency of programming depends on the patients’ accessibility and the clinical response to the stimulation. ### DBS Lead Localization Stereotactic computed tomography (CT) scans were performed on all patients to determine the lead location, typically on postoperative day nine at the point when the pneumocephalus was resolved. The CT image was fused to a preoperative scan to measure the DBS lead trajectory and contact positions relative to the midcommissural point in each case. We compared the stereotactic planning and DBS lead locations to measure the stereotactic errors. Stereotactic targeting error was calculated as the distance between the planned target point and trajectory of the implanted lead [19]. All imaging data were preprocessed using the Lead-DBS software ([www.lead-dbs.org](http://www.lead-dbs.org)) [20]. Briefly, all (pre- and postoperative) CT and MRIs were linearly co-registered to the preoperative T1WI using SPM ([https://www.fil.ion.ucl.ac.uk/spm/software/spm12/](https://www.fil.ion.ucl.ac.uk/spm/software/spm12/)). Registration between postoperative CT and preoperative T1WI was further refined using the linear registration within the subcortical target region of interest to minimise non-linear bias caused by the surgery [21]. All data were normalised into the standard MNI space. This procedure was performed with the whole-brain non-linear SyN registration implemented in Advanced Normalisation Tools (ANTs, [http://stnava.github.io/ANTs/](http://stnava.github.io/ANTs/)) [22] using the “effective (low variance)” setting with subcortical refinement as implemented in Lead-DBS [20]. Electrode trajectories and contacts were automatically pre-reconstructed using PaCER algorithm [23] and manually refined using Lead-DBS. We also determined the lead trajectory and contact positions using subject-independent 3D atlases of the thalamic nuclei [24]. We considered that contacts existed within those thalamic nuclei when spherical regions of 0.2 mm radius around contacts overlapped with each nucleus: CM, mediodorsal (MD), ventral lateral (VL), and ventral posterior (VP) nuclei. Detected electrodes and thalamic nuclei were visualised using Lead-DBS. ### Volume of Tissue Activated Mapping We calculated the volume of tissue activated (VTA) under the conditions with the optimised stimulation parameters (Table 3) and those with stimulation-induced side effects (Table 3, 4). All VTA calculations were conducted using Lead-DBS [25]. A volume conductor model was constructed on the basis of a tetrahedral volume mesh that included the DBS electrode and surrounding tissues. Conductivities of 0.14 S/m were assigned to grey and white matter [26]. On the basis of the volume conductor model, the electric field distribution was simulated using the FieldTrip-SimBio pipeline that was integrated into Lead-DBS ([https://www.mrt.uni-jena.de/simbio/index.php/](https://www.mrt.uni-jena.de/simbio/index.php/); [http://fieldtriptoolbox.org](http://fieldtriptoolbox.org)). The electric field distribution was thresholded for magnitudes above a commonly used value of 0.2 V/mm [27, 28] to define the extent of VTA. View this table: [Table 3.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/T3) Table 3. The stimulation parameters at last follow up and those induced mood changes. All side effects were reversal by reducing the stimulation intensity. *Interleaving stimulation settings were applied. Abbreviation. PW = pulse width View this table: [Table 4.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/T4) Table 4. Threshold levels of acute stimulation-induced side effects. Threshold levels of stimulation-induced side effects in each contact were measure at a fixed PW of 60 microseconds and frequency of 130 Hz. N/A = not applicable ### Normative Connectome To characterise the brain connectome differences from very adjacent VTAs and to elucidate the brain-wide mechanism of therapeutic stimulation and side effects, we used the population-averaged atlas of the macroscale human structural connectome derived from diffusion-weighted imaging (DWI) data (N = 842, Human Connectome Project) [29]. Considering the limited volume of data, the VTAs in the right hemisphere were non-linearly transformed into those in the left hemisphere using ANTs, and all VTAs were pooled across hemispheres. We created specific VTAs that did not overlap with the other VTAs and depicted the 300 fibres from them using DSI studio ([http://dsi-studio.labsolver.org](http://dsi-studio.labsolver.org)). All detected structural fibres were visualised using Lead DBS. We evaluated the proportion of the 300 fibres that passed through each brain region defined by Harvard-Oxford cortical/subcortical atlases [30] combined with the cerebellum from AAL atlas [31]. All small parcels of the cerebellum defined in AAL were integrated into one binarized parcel. The proportion of 300 fibres was thresholded above 0.02 and visualised using a circular plot. ### Statistical Analysis To compare the pre- and post-DBS clinical scores on YGTSS, Y-BOCS, and HAM-D, we used Wilcoxon signed rank test. Statistical analyses were performed using SPSS version 21.0 (IBM Corp., Armonk, NY, USA), and p < 0.05 indicated statistical significance. ## Results ### Clinical Outcomes The mean follow-up period was 10.9 ± 3.9 months, but one patient (case 6) missed a formal 6-month follow-up owing to the occurrence of the COVID-19 pandemic during the study period. YGTSS severity and impairment scores improved from 41.9 ± 3.9 and 47.5 ± 4.6 at the baseline to 15.8 ± 6.3 (z = −2.52, p = 0.012) and 18.8 ± 8.3 (z = −2.59, p = 0.010) at the last follow-up, respectively. Furthermore, Y-BOCS and HAM-D scores improved from 16.8 ± 11.5 and 8.6 ± 6.7 at the baseline to 6.1 ± 7.5 (z = −2.20, p = 0.028) and 4.8 ± 4.4 (z = −2.11, p = 0.035) at the last follow-up, respectively. Clinical outcomes are summarised in Table 2. It is noteworthy that one patient (Case 8) experienced an enormous microlesion effect while the tic movements had disappeared for a week after surgery. The surgical adverse events included wound dehiscence and lead misplacement. Two patients (cases 1 and 4) experienced wound dehiscence of the scalp incision site and underwent a wound revision. One patient (case 4) underwent a left DBS lead revision owing to insufficient phonic tic suppression despite the motor tic improvement. ### DBS Lead Location Coordinates of the preoperative stereotactic targeting and the measured DBS lead locations are summarised in Supplementary Table 1. Stereotactic targeting error was 1.2 ± 0.5 mm. Analysis of DBS lead positions indicated that 15 of 17 DBS leads, including a misplaced lead, penetrated the CM nucleus and that 11 leads penetrated the MD nucleus (Supplementary Table 2). The dorsal contacts of the quadripolar electrodes were located in the ventrolateral nucleus in all cases except for the initial left lead of case 4. In case 4, the tip of the misplaced DBS lead on the left was on the border between MD and CM nuclei, and three contacts were positioned in the MD nucleus. Two ventral contacts and two dorsal contacts of the revised lead in case 4 were successfully placed in the CM nucleus and ventrolateral nucleus, respectively (Supplementary Table 2 and Supplementary Movie 1). The DBS lead contact locations in the normalised brain space are shown in Fig. 1 and Supplementary Movie 1. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/30/2020.09.29.20200501/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/F1) Figure 1. Lead electrodes placement in MNI space. Lead electrodes with the CM nucleus (peach), RN (red), MD nucleus (purple), and ventral lateral dorsal/ventral nucleus (yellow). In (A), the ventral lateral dorsal/ventral nucleus was removed. The lead electrodes of each patient were displayed in different colours (case 1, orange; 2, red; 3, blue; 4, purple; 4 after the repositioning, dark purple (only in the left hemisphere); 5, light green; 6, yellow; 7, light blue: 8, green). ### Clinical Response According to the Stimulated Area Contacts 2 (10) and 3 (11) were likely to be selected as the active contacts (Table 3). Therapeutic effects were achieved by stimulation of the border between the CM and VL nuclei in our case series. Depressed mood or anxiety were reported after chronic stimulation (Table 3). Four patients commonly experienced mood changes after chronic stimulation for several days after a programming session. The programming parameters at the last follow-up and chronic side effects are summarised in Table 3. The common acute side effects were dizziness and paraesthesia in the contralateral upper extremity of the active DBS lead. Dizziness and paraesthesia tended to be observed with high-intensity stimulation of the two ventral contacts and two middle contacts, respectively. These acute stimulation-induced side effects are summarised in Table 4. ### Stimulated Areas and Normative Connectome In Fig. 2, we show the VTAs related to the therapeutic stimulation and side effects. We found the VTA related to the clinical response in the border between CM and VL nuclei (blue area in Fig. 2A). Although the paraesthesia-related region mostly overlapped the therapeutic region, it extended slightly into lateral, medial, and anterior directions (green area). The VTA related to dizziness extended into the ventral direction and overlapped with the RN (orange area). We found that chronic stimulation of the relatively medial and dorsal region entering the MD nucleus led to a depressed mood (purple area). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/30/2020.09.29.20200501/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/F2) Figure 2. VTAs related to therapeutic stimulation and side effects. Each colour of VTAs represents areas associated with the following effects: blue = therapeutic effect, orange = dizziness, green = paraesthesia, and purple = depressed mood. The peach-coloured region in the right hemisphere is the CM nucleus. (A) front, (B) medial, and (C) lateral view. Using the normative connectome, we found that each VTA related to the therapeutic stimulation and side effects showed clearly different network properties (Fig. 3). Brain regions connected with each VTA are summarised using the circular plot (Fig. 4, Supplementary Figure 2). Fibres of therapeutic stimulation were characterised by more dense connections with the precentral gyrus than were those of side effects (Fig. 4A). Dizziness- and paraesthesia-related VTAs extended fibres into relatively specific brain regions (Fig. 4B, 4C). In contrast, the stimulation-induced mood change was related to more spatially distributed brain regions (Fig. 4D). Specifically, we detected dizziness related fibres into the cerebellorubral network (Fig. 3B, 4B). The paraesthesia symptom was characterised by the fibres that connect the thalamus and insular cortex (Fig. 3C, 4C). We found a relatively dense connection with the thalamus in the amygdala and the orbitofrontal cortex, which are related to depressed mood (Fig. 3D, 4D). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/30/2020.09.29.20200501/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/F3) Figure 3. Normative connectome from VTAs related to therapeutic stimulation and side effects. The normative connectome from VTAs related to the (A) therapeutic stimulation, (B) dizziness, (C) paraesthesia, and (D) depressed mood. The peach-coloured region is the CM nucleus. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/30/2020.09.29.20200501/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/09/30/2020.09.29.20200501/F4) Figure 4. Proportion of the structural fibres projected from VTAs related to therapeutic stimulation and side effects. The circular plot represents the proportion of projected fibres that connect the thalamus and the other brain regions which have at least one connectivity (see Supplementary Figure 2 for all brain regions); (A) therapeutic stimulation, (B) dizziness, (C) paraesthesia, and (D) depressed mood. Abbreviations: ant = anterior; inf = inferior; mid = middle; post = posterior; sup = superior. ## Discussion Studies that evaluated clinical outcomes of DBS in patients with TS have reported mixed results on a case-by-case basis although overall outcomes were favourable [2-4], and our results were consistent with those of the past reports. The results of our study showed a variety of DBS contacts (Fig. 1, Supplementary Movie 1) and illustrated the relationship between the DBS lead locations and clinical effects inclusive of therapeutic and side effects (Fig. 2, 3, 4). Regarding the stereotactic planning technique, the lead trajectory is determined according to the safety issues [32]. Indirect targeting using the same template coordinates is not appropriate as this technique does not consider the anatomical variations in the subcortical structures [9]. According to the literature, DBS leads were presumed to pass through the ventral-oralis complex (Vo) nucleus to the CM nucleus; however, our series showed that there were a variety of lead locations even among the responders (Fig. 1, Supplementary Movie 1). We consider that such differences may have led to the occurrence of difference in threshold levels of stimulation-induced side effects (Table 3 and 4). Various factors associated with the surgical procedure such as brain shift and distortion of the frame are considered to deviate the lead position from the planned target [19]. Although we did not measure the electrode positions intraoperatively, an image-guided procedure using intraoperative CT or MRI may improve the stereotactic accuracy and precision [33, 34]. A potential factor predicting the incidence of a favourable outcome is the microlesion effect, and this issue also underpins the importance of the accurate lead implantation. Microlesion effect has been reported to be associated with the optimal lead position in other movement disorders such as Parkinson’s disease and essential tremor [35, 36]. Similarly, we have noted a temporary but complete resolution of tic movements immediately after surgery in a patient (case 8). These findings indicate that stimulation of a specific area is associated with tic suppression. In this context, clinical outcomes of TS DBS may be improved with accurate and precise DBS lead implantation in the specific area. Additionally, our experience that DBS lead repositioning improved the clinical outcome (case 4) addresses the importance of implanting DBS leads in the optimal position to maximise the benefit. Difficulty in DBS programming for TS exists, in that, immediate response may not be observed in the clinical setting unlike that in the case of Parkinson’s disease or essential tremor. Practitioners programming the DBS, however, should understand the relationship between the DBS lead and the surrounding structures. A subtle difference in stimulating position can lead to various effects through different properties of the brain network as the thalamus is quite a dense structure of many small nuclei that relays information between different brain regions [24]. The VTA and normative connectivity analyses linked our clinical findings with how electrical currents spreads to the surrounding neuroanatomical structures according to the situations (Fig. 2, 3, 4). The VTA associated with therapeutic effects covered the dorsal area of the CM nucleus and the VL nucleus, and our normative connectome analysis showed its association with the motor networks (Fig. 3A, 4A). This finding is consistent with that reported by a recent study showing that the relationship of tic reduction with the fibres connecting the thalamus and the motor cortex [37]. On the contrary, it is well known that high electrical intensity is usually required for tic suppression [4], and therapeutic stimulation may also modulate the limbic systems sub-clinically inducing psychiatric side effects. With regard to the immediate side effects, the ventralis caudalis (Vc) nucleus and RN are located posterolateral and ventral to the CM nucleus, respectively. Spread of stimulation current to the Vc nucleus and RN may induce paraesthesia and dizziness, respectively. These findings were partly supported by the findings of normative connectome analyses. These side effects can be observed immediately when the stimulation intensity is above the threshold levels in the clinical setting, but mood changes are likely to be detected several days after increasing the stimulation intensity. The lower threshold levels of these stimulation-induced side effects may indicate lead misplacement, and clinicians should be aware that the careful evaluation of clinical response to each stimulation parameter is important to estimate the lead location. The reported mood changes are considered to result from the use of supra-threshold intensity for the electrical stimulation of the limbic network, and this finding was supported by that of the normative connectome analysis, that is, the electrically activated area associated with the mood change is connected to limbic structures such as the amygdala (Fig. 3D, 4D). We consider that the high intensity of electrical current to the MD nucleus was associated with the side effects, and this may underpin the mechanism of loss of benefits due to side effects such as reduced levels of energy following long-term thalamic DBS reported recently [5]. In our study, the depressed moods were temporary because we decreased the stimulation intensity or changed active contacts when patients experienced the side effects. However, chronic DBS may induce irreversible effects as reported recently [38], and chronic high-intensity stimulation of the MD nucleus may result in irreversible mood changes or loss of beneficial effects. Although our study has demonstrated promising clinical outcomes of DBS therapy for severe, medication-refractory TS and provides a guide for identifying stimulation areas in the thalamus that yield desirable effects, it has several limitations. The evaluation was not blinded, and records of mood change in our study were based on patients’ reports rather than on quantified data. To confirm our findings, multi-centre studies with a higher number of patients are warranted. A registry collecting meticulous data from multiples centres may address these issues as reported by several authors [3, 8, 39]. ## Conclusion Our study addresses the importance of accurate implantation of DBS electrodes for obtaining standardised clinical outcomes and suggests that meticulous programming with careful monitoring of clinical symptoms may improve outcomes. Clinicians should attempt to detect any subtle changes in the clinical symptoms at each clinical visit for better stimulation adjustment, and in this context, findings of our study may be useful with regard to the systematic programming paradigm. Further meticulous evaluation of clinical effects of neurostimulation at specific area with a large study population is warranted. ## Data Availability Data Availability The clinical outcomes of our cohort were registered to the International Tourette Deep Brain Stimulation Database and Registry ([https://tourettedeepbrainstimulationregistry.ese.ufhealth.org/](https://tourettedeepbrainstimulationregistry.ese.ufhealth.org/)). The personal data including the imaging data obtained from the study subjects will not be distributed openly to protect the patients' privacy. We used the population-averaged atlas of the macroscale human structural connectome derived from DWI ([http://dsi-studio.labsolver.org/](http://dsi-studio.labsolver.org/)). For imaging analysis, we used free software inclusive of ANTs ([http://stnava.github.io/ANTs/](http://stnava.github.io/ANTs/)) (Avants et al., 2008), Lead-DBS ([https://www.lead-dbs.org/](https://www.lead-dbs.org/)) (Horn et al., 2019), PaCER algorithm (Husch et al., 2018), FieldTrip-SimBio ([https://www.mrt.uni-jena.de/simbio/index.php/](https://www.mrt.uni-jena.de/simbio/index.php/); [http://fieldtriptoolbox.org](http://fieldtriptoolbox.org)), and DSI studio ([http://dsi-studio.labsolver.org](http://dsi-studio.labsolver.org)). The subject-independent 3D atlases of the thalamic nuclei are also openly available online ([http://www.humanmotorthalamus.com/](http://www.humanmotorthalamus.com/)) (Ilinsky et al., 2018). ## Fundings This study was partly supported by Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Scientific Research (C) (Grant number: 18K08956), the Central Research Institute of Fukuoka University (Grant number: 201045), and JSPS KAKENHI Grant (Grant number: JP16H06396). ## Acknowledgements None ## Abbreviations AC : Anterior commissure ANTs : Advanced Normalisation Tools CM : centromedian DBS : Deep brain stimulation DWI : Diffusion-weighted imaging FGATIR : fast grey matter acquisition T1 inversion recovery GPi : globus pallidus interna IPG : Implantable pulse generator MNI : Montreal Neurological Institute MRI : Magnetic resonance imaging PC : Posterior commissure PW : Pulse width RN : Red nucleus TS : Tourette syndrome VTA : Volume of tissue activated YGTSS : Yale Global Tic Severity Scale * Received September 29, 2020. * Revision received September 29, 2020. * Accepted September 30, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].Vandewalle V, van der Linden C, Groenewegen HJ, Caemaert J. Stereotactic treatment of Gilles de la Tourette syndrome by high frequency stimulation of thalamus. Lancet 1999;353(9154):724. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(98)05964-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10073521&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078966400019&link_type=ISI) 2. [2].Baldermann JC, Schuller T, Huys D, Becker I, Timmermann L, Jessen F, et al. Deep Brain Stimulation for Tourette-Syndrome: A Systematic Review and Meta-Analysis. Brain Stimul 2015. 3. [3].Martinez-Ramirez D, Jimenez-Shahed J, Leckman JF, Porta M, Servello D, Meng FG, et al. Efficacy and Safety of Deep Brain Stimulation in Tourette Syndrome: The International Tourette Syndrome Deep Brain Stimulation Public Database and Registry. JAMA Neurol 2018;75(3):353–9. 4. [4].Dowd RS, Pourfar M, Mogilner AY. Deep brain stimulation for Tourette syndrome: a single-center series. J Neurosurg 2018;128(2):596–604. 5. [5].Smeets A, Duits AA, Leentjens AFG, Schruers K, van Kranen-Mastenbroek V, Visser-Vandewalle V, et al. Thalamic Deep Brain Stimulation for Refractory Tourette Syndrome: Clinical Evidence for Increasing Disbalance of Therapeutic Effects and Side Effects at Long-Term Follow-Up. Neuromodulation 2018;21(2):197–202. 6. [6].Morishita T, Foote KD, Burdick AP, Katayama Y, Yamamoto T, Frucht SJ, et al. Identification and management of deep brain stimulation intra- and postoperative urgencies and emergencies. Parkinsonism Relat Disord 2010;16(3):153–62. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.parkreldis.2009.10.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19896407&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 7. [7].Rolston JD, Englot DJ, Starr PA, Larson PS. An unexpectedly high rate of revisions and removals in deep brain stimulation surgery: Analysis of multiple databases. Parkinsonism Relat Disord 2016;33:72–7. 8. [8].Johnson KA, Fletcher PT, Servello D, Bona A, Porta M, Ostrem JL, et al. Image-based analysis and long-term clinical outcomes of deep brain stimulation for Tourette syndrome: a multisite study. J Neurol Neurosurg Psychiatry 2019;90(10):1078–90. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czoxMDoiOTAvMTAvMTA3OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzMwLzIwMjAuMDkuMjkuMjAyMDA1MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 9. [9].Nestor KA, Jones JD, Butson CR, Morishita T, Jacobson CEt, Peace DA, et al. Coordinate-based lead location does not predict Parkinson’s disease deep brain stimulation outcome. PLoS One 2014;9(4):e93524. 10. [10].Witt K, Granert O, Daniels C, Volkmann J, Falk D, van Eimeren T, et al. Relation of lead trajectory and electrode position to neuropsychological outcomes of subthalamic neurostimulation in Parkinson’s disease: results from a randomized trial. Brain 2013;136(Pt 7):2109–19. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awt151&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23801735&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000321060200018&link_type=ISI) 11. [11].Van Essen DC, Ugurbil K, Auerbach E, Barch D, Behrens TE, Bucholz R, et al. The Human Connectome Project: a data acquisition perspective. Neuroimage 2012;62(4):2222– 31. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.02.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22366334&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308265200006&link_type=ISI) 12. [12].Wang Q, Akram H, Muthuraman M, Gonzalez-Escamilla G, Sheth SA, Oxenford S, et al. Normative vs. patient-specific brain connectivity in Deep Brain Stimulation. Neuroimage 2020:117307. 13. [13].Pringsheim T, Okun MS, Muller-Vahl K, Martino D, Jankovic J, Cavanna AE, et al. Practice guideline recommendations summary: Treatment of tics in people with Tourette syndrome and chronic tic disorders. Neurology 2019;92(19):896–906. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 14. [14].Storch EA, Murphy TK, Geffken GR, Sajid M, Allen P, Roberti JW, et al. Reliability and validity of the Yale Global Tic Severity Scale. Psychol Assess 2005;17(4):486–91. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/1040-3590.17.4.486&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16393016&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234545500013&link_type=ISI) 15. [15].Storch EA, Shapira NA, Dimoulas E, Geffken GR, Murphy TK, Goodman WK. Yale-Brown Obsessive Compulsive Scale: the dimensional structure revisited. Depress Anxiety 2005;22(1):28–35. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/da.20088&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16130118&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231997400004&link_type=ISI) 16. [16].Hamilton M. Development of a rating scale for primary depressive illness. Br J Soc Clin Psychol 1967;6(4):278–96. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.2044-8260.1967.tb00530.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6080235&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 17. [17].Morishita T, Higuchi MA, Saita K, Tsuboi Y, Abe H, Inoue T. Changes in Motor-Related Cortical Activity Following Deep Brain Stimulation for Parkinson’s Disease Detected by Functional Near Infrared Spectroscopy: A Pilot Study. Front Hum Neurosci 2016;10:629. 18. [18].Morishita T, Higuchi MA, Kobayashi H, Abe H, Higashi T, Inoue T. A retrospective evaluation of thalamic targeting for tremor deep brain stimulation using high-resolution anatomical imaging with supplementary fiber tractography. J Neurol Sci 2019;398:148–56. 19. [19].Zrinzo L. Pitfalls in precision stereotactic surgery. Surg Neurol Int 2012;3(Suppl 1):S53–61. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22826812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 20. [20].Horn A, Li N, Dembek TA, Kappel A, Boulay C, Ewert S, et al. Lead-DBS v2: Towards a comprehensive pipeline for deep brain stimulation imaging. Neuroimage 2019;184:293–316. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 21. [21].Schonecker T, Kupsch A, Kuhn AA, Schneider GH, Hoffmann KT. Automated optimization of subcortical cerebral MR imaging-atlas coregistration for improved postoperative electrode localization in deep brain stimulation. AJNR Am J Neuroradiol 2009;30(10):1914–21. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czoxMDoiMzAvMTAvMTkxNCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzMwLzIwMjAuMDkuMjkuMjAyMDA1MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 22. [22].Avants BB, Epstein CL, Grossman M, Gee JC. Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal 2008;12(1):26–41. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.media.2007.06.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17659998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254032800004&link_type=ISI) 23. [23].Husch A, M VP, Gemmar P, Goncalves J, Hertel F. PaCER - A fully automated method for electrode trajectory and contact reconstruction in deep brain stimulation. Neuroimage Clin 2018;17:80–9. 24. [24].Ilinsky I, Horn A, Paul-Gilloteaux P, Gressens P, Verney C, Kultas-Ilinsky K. Human Motor Thalamus Reconstructed in 3D from Continuous Sagittal Sections with Identified Subcortical Afferent Territories. eNeuro 2018;5(3). 25. [25].Horn A, Reich M, Vorwerk J, Li N, Wenzel G, Fang Q, et al. Connectivity Predicts deep brain stimulation outcome in Parkinson disease. Ann Neurol 2017;82(1):67–78. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.24974&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28586141&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 26. [26].Vorwerk J, Cho JH, Rampp S, Hamer H, Knosche TR, Wolters CH. A guideline for head volume conductor modeling in EEG and MEG. Neuroimage 2014;100:590–607. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.06.040&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24971512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 27. [27].Astrom M, Diczfalusy E, Martens H, Wardell K. Relationship between neural activation and electric field distribution during deep brain stimulation. IEEE Trans Biomed Eng 2015;62(2):664–72. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TBME.2014.2363494&link_type=DOI) 28. [28].Hemm S, Mennessier G, Vayssiere N, Cif L, Coubes P. Co-registration of stereotactic MRI and isofieldlines during deep brain stimulation. Brain Res Bull 2005;68(1-2):59–61. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1016/j.brainresbull.2005.08.024&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16325005&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 29. [29].Yeh FC, Panesar S, Fernandes D, Meola A, Yoshino M, Fernandez-Miranda JC, et al. Population-averaged atlas of the macroscale human structural connectome and its network topology. Neuroimage 2018;178:57–68. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2018.05.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29758339&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) 30. [30].Makris N, Goldstein JM, Kennedy D, Hodge SM, Caviness VS, Faraone SV, et al. Decreased volume of left and total anterior insular lobule in schizophrenia. Schizophr Res 2006;83(2-3):155–71. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2005.11.020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16448806&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237012200005&link_type=ISI) 31. [31].Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 2002;15(1):273–89. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.2001.0978&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11771995&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173174900027&link_type=ISI) 32. [32].Morishita T, Okun MS, Burdick A, Jacobson CEt, Foote KD. Cerebral Venous Infarction: A Potentially Avoidable Complication of Deep Brain Stimulation Surgery. Neuromodulation 2013. 33. [33].Burchiel KJ, McCartney S, Lee A, Raslan AM. Accuracy of deep brain stimulation electrode placement using intraoperative computed tomography without microelectrode recording. J Neurosurg 2013;119(2):301–6. 34. [34].Larson PS, Starr PA, Martin AJ. Deep Brain Stimulation: Interventional and Intraoperative MRI Approaches. Prog Neurol Surg 2018;33:187–97. 35. [35].Mann JM, Foote KD, Garvan CW, Fernandez HH, Jacobson CEt, Rodriguez RL, et al. Brain penetration effects of microelectrodes and DBS leads in STN or GPi. J Neurol Neurosurg Psychiatry 2009;80(7):794–7. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI4MC83Lzc5NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzMwLzIwMjAuMDkuMjkuMjAyMDA1MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 36. [36].Morishita T, Foote KD, Wu SS, Jacobson CE, Rodriguez RL, Haq IU, et al. Brain penetration effects of microelectrodes and deep brain stimulation leads in ventral intermediate nucleus stimulation for essential tremor. J Neurosurg 2010;112(3):491–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3171/2009.7.JNS09150&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19663554&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275043900004&link_type=ISI) 37. [37].Andrade P, Heiden P, Hoevels M, Schlamann M, Baldermann JC, Huys D, et al. Modulation of Fibers to Motor Cortex during Thalamic DBS in Tourette Patients Correlates with Tic Reduction. Brain Sci 2020;10(5). 38. [38].Kimura Y, Ikegaya N, Iijima K, Takayama Y, Kaneko Y, Omori M, et al. Withdrawal of deep brain stimulation in patients with gilles de la tourette syndrome. Mov Disord 2019;34(12):1925–6. 39. [39].Schrock LE, Mink JW, Woods DW, Porta M, Servello D, Visser-Vandewalle V, et al. Tourette syndrome deep brain stimulation: a review and updated recommendations. Mov Disord 2015;30(4):448–71. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mds.26094&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25476818&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F30%2F2020.09.29.20200501.atom)