Next Article in Journal / Special Issue
Deep Brain Stimulation for Gilles de la Tourette Syndrome: Toward Limbic Targets
Previous Article in Journal
Boosting Slow Oscillations during Sleep to Improve Memory Function in Elderly People: A Review of the Literature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modulation of Fibers to Motor Cortex during Thalamic DBS in Tourette Patients Correlates with Tic Reduction

1
Department of Stereotactic and Functional Neurosurgery, University Hospital of Cologne, 50397 Cologne, Germany
2
Department of Neurosurgery, University Hospital of Cologne, 50397 Cologne, Germany
3
Department of Neuroradiology, University Hospital of Cologne, 50397 Cologne, Germany
4
Department of Psychiatry and Psychotherapy, University Hospital of Cologne, 50397 Cologne, Germany
5
Department of Neurology, University Hospital of Cologne, 50397 Cologne, Germany
*
Author to whom correspondence should be addressed.
Brain Sci. 2020, 10(5), 302; https://doi.org/10.3390/brainsci10050302
Submission received: 12 April 2020 / Revised: 6 May 2020 / Accepted: 11 May 2020 / Published: 15 May 2020
(This article belongs to the Special Issue Deep Brain Stimulation and Tourette Syndrome)

Abstract

:
Probabilistic tractography in Tourette syndrome (TS) patients have shown an alteration in the connectivity of the primary motor cortex and supplementary motor area with the striatum and thalamus, suggesting an abnormal connectivity of the cortico-striatum-thalamocortical-pathways in TS. Deep brain stimulation (DBS) of the centromedian nucleus–nucleus ventrooralis internus (CM-Voi complex) in the thalamus is an effective treatment for refractory TS patients. We investigated the connectivity of activated fibers from CM-Voi to the motor cortex and its correlation between these projections and their clinical outcome. Seven patients with TS underwent CM-Voi-DBS surgery and were clinically evaluated preoperatively and six months postoperatively. We performed diffusion tensor imaging to display the activated fibers projecting from the CM-Voi to the different motor cortex regions of interest. These analyses showed that the extent of tic reduction during DBS is associated with the degree of stimulation-dependent connectivity between CM-Voi and the motor cortex, and in particular, an increased density of projections to the presupplementary motor area (preSMA). Non-responder patients displayed the largest amount of active fibers projecting into cortical areas other than motor cortex compared to responder patients. These findings support the notion that an abnormal connectivity of thalamocortical pathways underlies TS, and that modulation of these circuits through DBS could restore the function and reduce symptoms.

1. Introduction

Tourette syndrome (TS) is defined as a chronic neuropsychiatric disorder characterized by sudden, uncontrolled, repetitive, stereotypical movements or vocalizations described as tics, typically starting during childhood [1]. The pathophysiology lying behind this disorder is still unclear. Some studies suggest a higher cortical excitability in the primary motor cortex [2], which could result from underlying alterations in the cortico-basal ganglia-thalamocortical (CBTC) loops [3]. TS patients have a higher functional disorganization and an increased global interaction in the CBTC circuits, suggesting a deficiency in brain [4]. Probabilistic tractography of TS patients have shown an alteration in the connectivity of the primary motor cortex and supplementary motor area with the striatum and the thalamus. However, these reports resulted in contradictory conclusions, with some studies showing a decreased connectivity [5,6], while others showed an enhanced connectivity [7,8].
More than half of the patients become tic-free or only present with mild symptoms after reaching early adulthood, therefore not needing further treatment [9]. In severe cases refractory to behavioral therapies and pharmacological treatment that persist during adulthood, deep brain stimulation (DBS) can be considered as a therapeutic option [10]. Up until now, there are several different targets along the CSPT circuit that have been explored in TS patients, the most commonly used being the thalamus, the globus pallidus internus (GPi), the globus pallidus externus (GPe), the anterior limb of the internal capsule/nucleus accumbens (ALIC/NA) and the subthalamic nucleus (STN) [11]. Nevertheless, among these targets, the majority of cases have been performed on the medial thalamus and the GPi. Within the medial thalamus, various possible targets have been described, the most frequently used being the centromedian nucleus–nucleus ventrooralis internus (CM-Voi) and the centromedian nucleus–parafascicular (CM-Pf) complexes [12]. The clinical outcome of TS DBS in these targets is relatively variable within the same cohorts, with none showing a clear advantage over the others [13,14].
In recent years, diffusion tensor imaging (DTI) has emerged as a useful tool for optimizing targets in DBS for different anatomical targets. DTI allows to display fibers projecting from the tissue activated at a nominated target to different cortical and subcortical regions through tractography [15]. Previous studies on different neurological and psychiatric disorders have shown direct correlations between connectivity patterns and clinical outcome in patients, for instance, in obsessive-compulsive disorder (OCD) [16], Parkinson’s disease [17] and tremor [18]. The aim of this study is firstly to examine the connectivity of activated fibers in TS patients during CM-Voi DBS and secondly, to elucidate if connectivity patterns are correlated with a better clinical outcome after long-term CM-Voi DBS.

2. Materials and Methods

2.1. Patients and Study Design

We retrospectively analyzed seven patients who underwent surgery for bilateral thalamic DBS for TS at the University Hospital of Cologne between 2016 and 2018. Patients were diagnosed based on the diagnostic and statistical manual of mental disorders fifth edition and the Tourette syndrome diagnostic confidence index [1]. The clinical evaluation was assessed using the Yale Global Tic Severity Scale (YGTSS), which was performed before, and six months after the surgery. This study was approved by the Medical Ethical Committee of the University Hospital of Cologne (Code 19-1421). All subjects granted their consent according to the Declaration of Helsinki.

2.2. Stereotactic Surgery

The CM-Voi complex was targeted for stimulation with standard coordinates (midpoint AC/PC line coordinates: x= +/− 5 mm, y= −4 mm, z= 0 mm), adapted to the patient’s anatomy and with trajectory along the longest axis of Voi, based on preoperative magnetic resonance imaging (MRI) fused with stereotactic computed tomography (CT) scan with the atlas of the human brain (Schaltenbrand Atlas) overlaid (Cranial Software, version StealthStation™ S7, Medtronic, Minneapolis, MN, USA). Surgery was performed under analgosedation with intraoperative microelectrode recordings. All patients received quadripolar electrodes (Medtronic 3389) bilaterally, stereotactically guided. The accurate localization of the electrodes was confirmed via postsurgical CT scan. For optimal stimulation, the frequency, voltage and pulse width were individually adapted for best results.

2.3. Fiber Tracking

Based on previous studies [5,8], we selected the primary motor cortex (M1), the supplementary motor area (SMA) and pre-SMA as our seed regions. We defined M1 manually in the ICBM 152 MNI 2009b asymmetrical cortical atlas and used the Human Motor Area Template [19] for SMA and preSMA. These cortical mask regions were transferred to each patient’s T2 sequence using FLIRT from the FSL (FMRIB Software Library, www.fmrib.ox.ac.uk/fsl) v6.0 program for affine linear matching, followed by deformable registration by the advanced normalization tools (ANTs, http://stnava.github.io/ANTs/) [20]. The seed regions were subsequently manually adjusted using the ITK-SNAP desktop app to each patient’s individual anatomy based on T1 sequences on the preoperative MRI [21].
Using the open-source software Lead-DBS (www.lead-dbs.org, Charité University Hospital, Berlin, Germany), we coregistered and normalized the preoperative MRI scans and the postsurgical CT scan to MNI space [22]. Afterwards, we reconstructed the electrodes and the active contacts. Using the model described by Dembek and colleagues [23], we estimated the area stimulated by the active contacts (VTA: volume of tissue activated) based on the stimulation parameters, with a general heuristic electric-field threshold of 0.2 V/mm. We transferred the VTAs, using the same transformation procedure as described above, into patient space. The DTIs were corrected for distortions and subject movement using FUGUE from FSL tools [24,25,26,27]. After this, we performed the probabilistic fiber tracking from the VTAs to the ipsilateral hemisphere as well as the ipsilateral M1, SMA and preSMA using the Oxford Centre for Functional MRI of the Brain (FMRIB), FSL probtrackx2 program [28,29].
For further analysis, we applied the seed regions as binary masks on the respective tracks using the fslmaths tool and quantified the mean and maximum intensity of voxels localized within the seed regions.

2.4. Statistical Analysis

All data were analyzed using GraphPad Prism software (version 5, GraphPad Software, San Diego, CA, USA) with p-values under 0.05 considered as significant. In order to estimate the correlation between clinical improvement in the YGTSS and the fiber tract values, data were analyzed using a linear regression analysis. A Wilcoxon signed-rank test was performed to compare pre- and post-operative clinical YGTSS score differences between the samples. A Spearman rank’s test was performed to establish correlation analyses between different parameters. A power calculation was performed in order to calculate the amount of TS patients required for a significant statistical analysis (90% power, 0.05 alpha).
The data that support the findings of this study, such as the DBS MRI datasets, are not publicly available due to data privacy regulations, but are available from the corresponding author upon reasonable request.

3. Results

3.1. Demographic and Clinical Characteristics

Demographic and individual clinical data from patients included in this study (n = 7, six men and one woman, age mean 30 years SD = 9.1) are presented in Table 1. Disease onset ranged from three to ten years of age. Comorbidities were present among the participants with the most frequent conditions including OCD and depression (n = 2, 28% each). Disability and impairment were severe in every patient with a mean total preoperative YGTSS score of 90 (SD = 8.7). The mean YGTSS score six months postoperatively was 52 (SD = 21.4), with a significant mean improvement in YGTSS score of 43.5% (SD = 18.9, p < 0.001) pre- versus post-operatively. Clinical response being defined as an improvement of at least 35% on the YGTSS, five patients were considered responders (mean improvement in YGTSS 52.7% with a range of 43.7–71.8%) and two as non-responders (mean improvement in YGTSS 20.5% with a range of 11–30%).

3.2. Neuroimaging Analysis

3.2.1. Generation of Fiber Tracts and Cortical Projections

The individual location of bilateral electrodes and the generated VTAs of all patients are shown in Figure 1A,B, respectively. A clear difference is seen among the generated VTAs between responders and non-responders, as shown in Figure 1C,D, respectively. Non-responders presented with significant (p < 0.05) higher VTA-volumes (mean volume 2597 mm3, SD = 804) compared to responders (mean volume 1313 mm3, SD = 515). Furthermore, a direct comparison between volumetric values of individual VTAs and the percentage of improvement in the YGTSS six months after surgery showed a significant negative correlation in all participants (p < 0.01, R= −0.84), as shown in Figure 2A. Thus, the patients with better clinical outcomes showed smaller magnitudes in the generated VTAs. Derived from these results, we analyzed the amount of active fibers involved in the corresponding VTAs that projected to all cortical areas. This analysis showed a significant negative correlation between the amount of cortical fiber tracts and the percentage of improvement in the YGTSS six months after surgery (p < 0.05, R = −0.58), as displayed in Figure 2B. Furthermore, the non-responder patients presented with larger amount of fibers projecting to the cortex with higher variability between individuals (mean number of fibers 45,452,500, SD = 35,585,148) when compared to the responder participants (mean number of fibers 15,038,000, SD = 4,976,307). Therefore, the patients with better clinical outcomes displayed a general reduced cortical fiber activation in the whole brain when compared to the non-responders.

3.2.2. Parcellation of Seed Regions

Parcellation of the seed regions of interest in the motor cortex (preSMA, SMA and M1) showed that non-responder patients displayed the largest amount of active fibers, but projecting to cortical areas different from the seed regions than responder patients. This means that, although cortical connectivity in non-responders was generally increased, it was primarily projecting to nonmotor areas, as shown in Figure 3. These active fibers projected in its majority to different areas of the prefrontal cortex, as shown in Figure 3A. Responders displayed generally less modulated connectivity compared with non-responders, which was highly clustered towards the preSMA, SMA and M1 cortices, as shown in Figure 3B. Moreover, the total amount of cortical fibers originated in the VTAs that projected to nonmotor seed regions showed a negative correlation with the percentage of improvement in the YGTSS six months after surgery (p < 0.05, R = −0.60), as displayed in Figure 3C.

3.2.3. Density of Fibers in the Motor Cortex

The tractography of the three seed regions analyzed independently showed that all patients displayed active fibers projecting to all cortices including preSMA, SMA and M1, as shown in Figure 4. The density of the projections to the motor cortex was analyzed by investigating the mean amount of activated fibers per voxel inside the seed regions. This analysis showed a significant positive correlation between the mean density in all seed regions together (p < 0.05, R = 0.63) as well as in preSMA (p < 0.05, R = 0.61) in association with the percentage of improvement in the YGTSS six months after surgery, demonstrating a higher density of projecting fibers within these areas in patients with better clinical outcomes, as shown in Figure 5A,B. SMA and M1 showed both a positive tendency on the mean density, however, it was not significant when analyzed independently (p > 0.05 both, R = 0.25 and R = 0.46, respectively).

4. Discussion

In the present study we investigated the connectivity between the CM-Voi complex in the thalamus and the motor cortex during DBS in TS patients. Our findings showed that an increase in the density of fiber projections to the seed regions of the motor cortex, including preSMA, SMA and M1, was accompanied with better clinical outcomes. These results support the concept that abnormal structural connectivity of thalamo-cortical pathways in TS may be involved in the generation of motor and vocal tics [8]. More specifically, activation of fibers projecting to the preSMA was significantly higher in responder patients. In line with these findings, previous studies have suggested that tic severity is correlated with activations in the premotor cortex (SMA and preSMA) [31]. Although the function of these cortical areas in TS is still to be determined, it has been proposed that its role in planning of movements, motor inhibition and therefore voluntary tic suppression may be related to symptom reduction [32,33]. Furthermore, there is evidence that the urge to tic is associated with the SMA region, thus DBS for TS may also act by reducing pathological urges to tic [34]. Although previous studies have reported apparent conflicting results about connectivity in TS patients, it is noteworthy to mention that certain differences in the design between studies can be observed [5,6,7,8]. For instance, Cheng and colleagues reflect if progression of the disease can play a role in reducing connectivity patterns in TS, since they studied a cohort of adult patients with a mean disease duration of almost three decades [5]. Nevertheless, Makki and colleagues reported comparable results with lower probability of connection between caudate nucleus and anterior-dorsolateral-frontal cortex on the left hemisphere in children with TS [6]. Another crucial difference between these two reports was the presence of comorbidities. While Cheng et al. selected a cohort of TS patients without psychiatric comorbidities, Makki et al. studied a group of individuals with accompanying psychiatric behavior like OCD. The presence of other disorders may complicate the analysis of these results, since multiple circuits may involve several cortical structures simultaneously. Although our results confirm the relevance of these cortical structures, there are still certain differences to these studies. Contrarily, Thomalla and colleagues showed enhanced connectivity between the thalamus and the SMA [7]. In this study, the cohort may be more comparable to our patients, where younger adult patients with similar duration disease duration where included. However, our group of TS patients was fairly heterogeneous, with an almost equally divided sample of patients with comorbidities and without them. Consequently, medication between participants was also quite diverse in our group. On the contrary, Thomalla et al. included TS patients with no psychiatric comorbidities and no medications at the inclusion to the study. Nevertheless, the same authors of this report discussed the very selected group of patients with these unique characteristics, which may not necessarily represent the majority of the TS population. In this regard, when we analyze the patients included in our study, two cases were also diagnosed with OCD, one of the most common comorbid conditions in TS. This finding could be interpreted by itself as a better representation of comorbidities of this disorder in our study than the ones reported in other studies. It could also be interpreted as a confounding factor, since these two patients with OCD had the worst clinical outcome. The connectivity profile of these two cases interestingly followed a similar pattern, showing a clear difference with the patterns of responders. Recapitulating, all these alterations in connectivity probably illustrate a process of adaptation between pathological input at diverse stages of the disease and the capacity of compensation according to the individual variables of each patient. In this context, it is relevant to mention that the connectivity displayed in our study is a stimulation-dependent connectivity, while the above-mentioned studies are based on structural connectivity paradigms. For this reason, it was of crucial importance that our DTIs were performed under anesthesia. This way we avoided movement-associated artifacts, since the slightest movement in TS patients could generate alterations in tractographies.
We also demonstrated that non-responders displayed larger amount of active fibers than responder patients, but projecting to cortical areas different from the seed regions. These findings suggest that diffuse stimulation of several cortical areas simultaneously may have therapeutic implications on these patients. Since non-responder patients also displayed projections into the preSMA, SMA and M1 cortices, (however, with a lower density than the good responders) these findings may suggest that a more selective connectivity to these seed regions can reproduce better clinical outcomes. These observations translated directly into the size of the generated VTAs and the amount of fibers activated by these VTAs that projected into the entire cortex (preSMA, SMA, M1 and others). This assessment showed a clear negative correlation between higher VTAs volumes and the improvement in the YGTSS scores. Therefore, the consequence of these larger VTAs was the activation of a larger quantity of fibers that projected to the general cortex. Previous studies in different psychiatric disorders have shown a similar negative effect by an extensive activation of nonrelevant fibers during stimulation of the anterior limb of the internal capsule and nucleus accumbens in OCD patients [16]. In this study, the authors suggest that the activation of nonrelated structures may counteract the positive effect reached by the targeted fibers. In order to validate this concept, after we finished the interpretation of our results, we reanalyzed our data using fixed-size VTAs of 1, 2 and 3 Volts in all patients instead of the actual individual parameters of each subject. This analysis showed the same connectivity patterns of increased fiber projections to the general cortex and not primarily to the preSMA, SMA and M1 cortices in non-responders, therefore, demonstrating that the amount of activated fibers would ultimately be the most relevant factor, and not exactly the size of the VTAs. Although larger VTAs could activate more fibers, apparently this principle is dependent on the exact anatomical area that is stimulated. Moreover, not only the amount of fibers projecting to relevant areas may be critical for a positive outcome, but also the capacity to exclusively modulate the areas of interest. In order to achieve this objective, a more precise anatomical location of the active contacts within the thalamus may prove to be an essential factor. Our target is based on the positive clinical effect on Tourette symptoms of thalamotomies performed by Hassler and Dieckmann in 1970. Hassler and Dieckmann performed more than ten lesions per hemisphere, mainly in the median and intralaminar thalamic nuclei as well as the inner part of the ventral oral thalamic nucleus [35]. The rationale for the first DBS case was to strategically define one target so that the different nuclei lesioned by Hassler could be modulated by implantation of one lead per hemisphere [36]. This target was chosen at the anterior part of the centromedian nucleus (as part of the intralaminar thalamic nuclei), at its medial border with the medial thalamic nuclei, and along the Voi. DBS of this specific target has proven to be effective in a randomized controlled trial [37].
Comparable clinical results have been reported in up to four different targets within the CM thalamic region [12], however anatomical differences between these nuclei should be taken into consideration. A recent study with a similar setting, in this case targeting the CM-Pf (parafascicular) region, showed a high variability between the clinical outcome and the connectivity [38]. In this report, the cortical projections associated with better clinical responses were the right frontal middle gyrus, the left frontal superior sulci region and the left cingulate sulci region. None of these regions were particularly involved in our study, suggesting a potential difference between the connectivity of the CM-Voi and the CM-Pf complexes.
A major limitation of our study is the limited size of our sample, in the context of such a heterogeneous population. As we mentioned before, although this factor may represent more properly the general TS population, in our case it could also fragment the data and extrapolate the outlier data. Therefore, larger cohorts with similar conditions and long-term follow-ups even after initial clinical failure are required to validate the findings of our study. However, to our knowledge, this is the first time that TS patients with ongoing stimulation of the CM-Voi have been investigated with stimulation-dependent tractographies based on high-definition dMRI data.

5. Conclusions

To summarize, we demonstrate that the extent of tic reduction in TS patients after CM-Voi DBS is associated with the degree of connectivity between the stimulated area and the motor cortex. In particular, an increased density of projections to the preSMA, SMA and M1 areas correlates with better clinical outcomes after chronic DBS. These findings reinforce the theory that an abnormal structural connectivity of thalamo-cortical pathways underlies TS, and therefore the modulation of these circuits through DBS could reduce motor and vocal tics.

Author Contributions

Conceptualization, P.A. and V.V.-V.; methodology, P.A., J.C.B., D.H. and V.V.-V.; software, P.A., P.H. and M.H.; validation, P.A., P.H., M.S. and V.V.-V; formal analysis, P.A., P.H. and V.V.-V.; investigation, P.A., J.C.B., D.H. and V.V.-V.; resources, P.A., M.H., J.C.B. and M.S.; data curation, P.A., P.H., J.C.B., D.H. and V.V.-V.; writing—original draft preparation, P.A. and P.H.; writing—review and editing, P.A., P.H., M.H., M.S., J.C.B., D.H. and V.V.-V.; visualization, P.H., M.H. and M.S.; supervision, P.A., P.H. and V.V.V.; project administration, P.A. and V.V.-V.; funding acquisition, N/A. All authors have read and agreed to the published version of the manuscript.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Acknowledgments

We would like to thank Julia Burkhardt for her expert assistance and help during the creation of the fiber tracking images.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Scatter plot depicting the correlation between postoperative improvement of YGTSS after six months expressed in percentage and the volumetric measurements of the generated VTAs ((A), p < 0.001, R = −0.74) and the numbers of fibers in the VTA that project to the cortex ((B), p < 0.01, R = −0.58). Data shown in this figure represent a segmented analysis of each hemisphere (n = 14).
Figure A1. Scatter plot depicting the correlation between postoperative improvement of YGTSS after six months expressed in percentage and the volumetric measurements of the generated VTAs ((A), p < 0.001, R = −0.74) and the numbers of fibers in the VTA that project to the cortex ((B), p < 0.01, R = −0.58). Data shown in this figure represent a segmented analysis of each hemisphere (n = 14).
Brainsci 10 00302 g0a1
Figure A2. Correlation between postoperative improvement of YGTSS after six months expressed in percentage and the mean of fibers per voxel located within the cortical seed regions. All regions of interest (ROIs, p < 0.05, R = 0.32) summed together including the preSMA, SMA and M1 are shown in (A), while the independent seed regions are shown separately in figures (B) (preSMA, p < 0.001, R = 0.71), (C) (SMA, p > 0.05, R = 0.16) and (D) (M1, p > 0.05, R = 0.07). Data shown in this figure represent a segmented analysis of each hemisphere (n = 14).
Figure A2. Correlation between postoperative improvement of YGTSS after six months expressed in percentage and the mean of fibers per voxel located within the cortical seed regions. All regions of interest (ROIs, p < 0.05, R = 0.32) summed together including the preSMA, SMA and M1 are shown in (A), while the independent seed regions are shown separately in figures (B) (preSMA, p < 0.001, R = 0.71), (C) (SMA, p > 0.05, R = 0.16) and (D) (M1, p > 0.05, R = 0.07). Data shown in this figure represent a segmented analysis of each hemisphere (n = 14).
Brainsci 10 00302 g0a2

References

  1. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders, 5th ed.; American Psychiatric Association Publishing: Arlington, VA, USA, 2013. [Google Scholar]
  2. Heise, K.-F.; Steven, B.; Liuzzi, G.; Thomalla, G.; Jonas, M.; Müller-Vahl, K.; Sauseng, P.; Münchau, A.; Gerloff, C.; Hummel, F.C. Altered modulation of intracortical excitability during movement preparation in Gilles de la Tourette syndrome. Brain 2009, 133, 580–590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ganos, C.; Roessner, V.; Münchau, A. The functional anatomy of Gilles de la Tourette syndrome. Neurosci. Biobehav. Rev. 2013, 37, 1050–1062. [Google Scholar] [CrossRef] [PubMed]
  4. Worbe, Y.; Malherbe, C.; Hartmann, A.; Pélégrini-Issac, M.; Messé, A.; Vidailhet, M.; Lehéricy, S.; Benali, H. Functional immaturity of cortico-basal ganglia networks in Gilles de la Tourette syndrome. Brain 2012, 135, 1937–1946. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cheng, B.; Braass, H.; Ganos, C.; Treszl, A.; Biermann-Ruben, K.; Hummel, F.C.; Müller-Vahl, K.; Schnitzler, A.; Gerloff, C.; Münchau, A.; et al. Altered intrahemispheric structural connectivity in Gilles de la Tourette syndrome. NeuroImage Clin. 2013, 4, 174–181. [Google Scholar] [CrossRef] [Green Version]
  6. Makki, M.I.; Govindan, R.M.; Wilson, B.J.; Behen, M.E.; Chugani, H.T. Altered Fronto-Striato-Thalamic Connectivity in Children with Tourette Syndrome Assessed with Diffusion Tensor MRI and Probabilistic Fiber Tracking. J. Child Neurol. 2009, 24, 669–678. [Google Scholar] [CrossRef]
  7. Götz, T.; Siebner, H.R.; Jonas, M.; Bäumer, T.; Biermann-Ruben, K.; Hummel, F.; Gerloff, C.; Müller-Vahl, K.; Schnitzler, A.; Orth, M.; et al. Structural changes in the somatosensory system correlate with tic severity in Gilles de la Tourette syndrome. Brain 2009, 132, 765–777. [Google Scholar] [CrossRef] [Green Version]
  8. Worbe, Y.; Marrakchi-Kacem, L.; LeComte, S.; Valabregue, R.; Poupon, F.; Guevara, P.; Tucholka, A.; Mangin, J.-F.; Vidailhet, M.; Lehericy, S.; et al. Altered structural connectivity of cortico-striato-pallido-thalamic networks in Gilles de la Tourette syndrome. Brain 2014, 138, 472–482. [Google Scholar] [CrossRef] [Green Version]
  9. Leckman, J.F.; Zhang, H.; Vitale, A.; Lahnin, F.; Lynch, K.; Bondi, C.; Kim, Y.-S.; Peterson, B.S. Course of tic severity in Tourette syndrome: The first two decades. Pediatrics 1998, 102, 14–19. [Google Scholar] [CrossRef]
  10. Jankovic, J.; Kurlan, R. Tourette syndrome: Evolving concepts. Mov. Disord. 2011, 26, 1149–1156. [Google Scholar] [CrossRef]
  11. Andrade, P.; Visser-Vandewalle, V. DBS in Tourette syndrome: Where are we standing now? J. Neural Transm. 2016, 123, 791–796. [Google Scholar] [CrossRef]
  12. A Johnson, K.; Fletcher, P.T.; Servello, D.; Bona, A.; Porta, M.; Ostrem, J.L.; Bardinet, E.; Welter, M.-L.; Lozano, A.M.; Baldermann, J.C.; 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, 1078–1090. [Google Scholar] [CrossRef] [PubMed]
  13. Baldermann, J.; Schüller, T.; Huys, D.; Becker, I.; Timmermann, L.; Jessen, F.; Visser-Vandewalle, V.; Kuhn, J.; Information, P.E.K.F.C. Deep Brain Stimulation for Tourette-Syndrome: A Systematic Review and Meta-Analysis. Brain Stimul. 2016, 9, 296–304. [Google Scholar] [CrossRef]
  14. Schrock, L.E.; Mink, J.W.; Woods, D.W.; Porta, M.; Servello, D.; Visser-Vandewalle, V.; Silburn, P.A.; Foltynie, T.; Walker, H.C.; Shahed-Jimenez, J.; et al. Tourette Syndrome Association International Deep Brain Stimulation (DBS) Database and Registry Study Group. Tourette syndrome deep brain stimulation: A review and updated recommendations. Mov. Disord. 2015, 30, 448–471. [Google Scholar] [CrossRef] [PubMed]
  15. Henderson, J.M. “Connectomic surgery”: Diffusion tensor imaging (DTI) tractography as a targeting modality for surgical modulation of neural networks. Front. Integr. Neurosci. 2012, 6, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hartmann, C.J.; Lujan, J.L.; Chaturvedi, A.; Goodman, W.K.; Okun, M.S.; McIntyre, C.C.; Haq, I.U. Tractography Activation Patterns in Dorsolateral Prefrontal Cortex Suggest Better Clinical Responses in OCD DBS. Front. Mol. Neurosci. 2016, 9, 1104. [Google Scholar] [CrossRef] [Green Version]
  17. Vanegas-Arroyave, N.; Lauro, P.M.; Huang, L.; Hallett, M.; Horovitz, S.G.; Zaghloul, K.; Lungu, C. Tractography patterns of subthalamic nucleus deep brain stimulation. Brain 2016, 139, 1200–1210. [Google Scholar] [CrossRef] [Green Version]
  18. Morishita, T.; Higuchi, M.-A.; 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–156. [Google Scholar] [CrossRef]
  19. Mayka, M.A.; Corcos, D.M.; Leurgans, S.E.; Vaillancourt, D.E. Three-dimensional locations and boundaries of motor and premotor cortices as defined by functional brain imaging: A meta-analysis. NeuroImage 2006, 31, 1453–1474. [Google Scholar] [CrossRef] [Green Version]
  20. Tustison, N.J.; Avants, B.B.; Gee, J.C. Learning image-based spatial transformations via convolutional neural networks: A review. Magn. Reson. Imaging 2019, 64, 142–153. [Google Scholar] [CrossRef]
  21. Yushkevich, P.A.; Piven, J.; Hazlett, H.C.; Smith, R.G.; Ho, S.; Gee, J.C.; Gerig, G. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. NeuroImage 2006, 31, 1116–1128. [Google Scholar] [CrossRef] [Green Version]
  22. Horn, A.; A Kühn, A. Lead-DBS: A toolbox for deep brain stimulation electrode localizations and visualizations. NeuroImage 2015, 107, 127–135. [Google Scholar] [CrossRef]
  23. Dembek, T.A.; Barbe, M.T.; Åström, M.; Hoevels, M.; Visser-Vandewalle, V.; Fink, G.R.; Timmermann, L. Probabilistic mapping of deep brain stimulation effects in essential tremor. NeuroImage Clin. 2016, 13, 164–173. [Google Scholar] [CrossRef] [PubMed]
  24. Andersson, J.L.; Skare, S.; Ashburner, J. How to correct susceptibility distortions in spin-echo echo-planar images: Application to diffusion tensor imaging. NeuroImage 2003, 20, 870–888. [Google Scholar] [CrossRef]
  25. Andersson, J.L.; Sotiropoulos, S.N. An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. NeuroImage 2016, 125, 1063–1078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Andersson, J.L.R.; Sotiropoulos, S.N. Non-parametric representation and prediction of single- and multi-shell diffusion-weighted MRI data using Gaussian processes. NeuroImage 2015, 122, 166–176. [Google Scholar] [CrossRef] [Green Version]
  27. Smith, S.M.; Jenkinson, M.; Woolrich, M.W.; Beckmann, C.F.; Behrens, T.E.; Johansen-Berg, H.; Bannister, P.R.; De Luca, M.; Drobnjak, I.; Flitney, D.E.; et al. Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage 2004, 23, S208–S219. [Google Scholar] [CrossRef] [Green Version]
  28. Behrens, T.E.; Woolrich, M.; Jenkinson, M.; Johansen-Berg, H.; Nunes, R.G.; Clare, S.; Matthews, P.M.; Brady, J.; Smith, S.M. Characterization and propagation of uncertainty in diffusion-weighted MR imaging. Magn. Reson. Med. 2003, 50, 1077–1088. [Google Scholar] [CrossRef]
  29. Behrens, T.E.; Johansen-Berg, H.; Jbabdi, S.; Rushworth, M.F.; Woolrich, M. Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? NeuroImage 2007, 34, 144–155. [Google Scholar] [CrossRef]
  30. Behrens, T.E.J.; Johansen-Berg, H.; Woolrich, M.W.; Smith, S.; Wheeler-Kingshott, C.A.G.; A Boulby, P.; Barker, G.J.; Sillery, E.L.; Sheehan, K.; Ciccarelli, O.; et al. Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat. Neurosci. 2003, 6, 750–757. [Google Scholar] [CrossRef]
  31. Polyanska, L.; Critchley, H.D.; Rae, C.L. Centrality of prefrontal and motor preparation cortices to Tourette Syndrome revealed by meta-analysis of task-based neuroimaging studies. NeuroImage Clin. 2017, 16, 257–267. [Google Scholar] [CrossRef]
  32. Zapparoli, L.; Porta, M.; Paulesu, E. The anarchic brain in action. Curr. Opin. Neurol. 2015, 28, 604–611. [Google Scholar] [CrossRef]
  33. Nachev, P.; Wydell, H.; O’Neill, K.; Husain, M.; Kennard, C. The role of the pre-supplementary motor area in the control of action. NeuroImage 2007, 36, T155–T163. [Google Scholar] [CrossRef] [Green Version]
  34. Cavanna, A.E.; Black, K.J.; Hallett, M.; Voon, V. Neurobiology of the Premonitory Urge in Tourette’s Syndrome: Pathophysiology and Treatment Implications. J. Neuropsychiatry Clin. Neurosci. 2017, 29, 95–104. [Google Scholar] [CrossRef] [Green Version]
  35. Hassler, R.; Dieckmann, G. Stereotaxic treatment of tics and inarticulate cries or coprolalia considered as motor obsessional phenomena in Gilles de la Tourette’s disease. Rev. Neurol. 1970, 123, 89–100. [Google Scholar] [PubMed]
  36. Vandewalle, V.; Van Der Linden, C.; Groenewegen, H.J.; Caemaert, J. Stereotactic treatment of Gilles de la Tourette syndrome by high frequency stimulation of thalamus. Lancet 1999, 353, 724. [Google Scholar] [CrossRef]
  37. Ackermans, L.; Duits, A.; Van Der Linden, C.; Tijssen, M.; Schruers, K.; Temel, Y.; Kleijer, M.; Nederveen, P.; Bruggeman, R.; Tromp, S.; et al. Double-blind clinical trial of thalamic stimulation in patients with Tourette syndrome. Brain 2011, 134, 832–844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Brito, M.; Teixeira, M.J.; Mendes, M.M.; França, C.; Iglesio, R.; Barbosa, E.R.; Cury, R.G. Exploring the clinical outcomes after deep brain stimulation in Tourette syndrome. J. Neurol. Sci. 2019, 402, 48–51. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Anatomical location of electrodes and volume of tissue activated (VTAs) within the thalamus. Overview of bilateral electrodes in the centromedian nucleus–nucleus ventrooralis internus (CM-Voi) complex (A) and cluster for bilateral VTAs of all patients on the MNI space (B) in sagittal, coronal and axial views. Volume of tissue activated three-dimensional representation shown in light red for the best clinical responder in the cohort, only patient 1 represented (C). Volume of tissue activated three-dimensional representation shown in light red for the patient with the worst clinical outcome in the cohort, only patient 7 represented (D). Thalamic segmentation according to cortical connectivity is shown in different colors in C and D in order to exemplify the position of the implanted electrodes [30].
Figure 1. Anatomical location of electrodes and volume of tissue activated (VTAs) within the thalamus. Overview of bilateral electrodes in the centromedian nucleus–nucleus ventrooralis internus (CM-Voi) complex (A) and cluster for bilateral VTAs of all patients on the MNI space (B) in sagittal, coronal and axial views. Volume of tissue activated three-dimensional representation shown in light red for the best clinical responder in the cohort, only patient 1 represented (C). Volume of tissue activated three-dimensional representation shown in light red for the patient with the worst clinical outcome in the cohort, only patient 7 represented (D). Thalamic segmentation according to cortical connectivity is shown in different colors in C and D in order to exemplify the position of the implanted electrodes [30].
Brainsci 10 00302 g001
Figure 2. Scatter plot depicting the correlation between postoperative improvement of the Yale Global Tic Severity Scale (YGTSS) after six months expressed in percentage and the volumetric measurements of the generated VTAs ((A), p < 0.01, R = −0.84) and the numbers of fibers in the VTA that project to the cortex ((B), p < 0.05, R = −0.58). Data shown in this figure represent individual patient analysis (n = 7). For a segmented analysis of each hemisphere (n = 14), please consult Figure A1.
Figure 2. Scatter plot depicting the correlation between postoperative improvement of the Yale Global Tic Severity Scale (YGTSS) after six months expressed in percentage and the volumetric measurements of the generated VTAs ((A), p < 0.01, R = −0.84) and the numbers of fibers in the VTA that project to the cortex ((B), p < 0.05, R = −0.58). Data shown in this figure represent individual patient analysis (n = 7). For a segmented analysis of each hemisphere (n = 14), please consult Figure A1.
Brainsci 10 00302 g002
Figure 3. Comparison of structural connectivity between a non-responder (A) and a good responder (B) patient, in sagittal and axial views. All fibers traversing the VTAs that connect into the preSMA (presupplementary motor area), SMA and M1 (primary motor cortex) cortices are shown in blue, while the rest of the fibers that project to other parts of the cortex are shown in red. Correlation between the clinical response and the total amount of cortical projections into nonmotor seed regions ((C), p < 0.05, R = −0.60).
Figure 3. Comparison of structural connectivity between a non-responder (A) and a good responder (B) patient, in sagittal and axial views. All fibers traversing the VTAs that connect into the preSMA (presupplementary motor area), SMA and M1 (primary motor cortex) cortices are shown in blue, while the rest of the fibers that project to other parts of the cortex are shown in red. Correlation between the clinical response and the total amount of cortical projections into nonmotor seed regions ((C), p < 0.05, R = −0.60).
Brainsci 10 00302 g003
Figure 4. Connectivity-based segmentation of the seed areas. Representative tractographies of each patient. Images are shown in patient-space, where the most significant density of fibers was observed. Therefore, the most symbolic MRI cut of each patient was depicted, including the most representative hemisphere in each individual. Fibers originated from the VTAs that connect to M1 are shown in green, whereas fibers connecting to SMA and preSMA are shown in orange and blue, respectively. Patients are separated in responders (green bar) and non-responders (red bar) in descending order from the best to the worst clinical outcome from left to right. Independently from their clinical outcome, all patients showed thalamic connectivity to the seed areas.
Figure 4. Connectivity-based segmentation of the seed areas. Representative tractographies of each patient. Images are shown in patient-space, where the most significant density of fibers was observed. Therefore, the most symbolic MRI cut of each patient was depicted, including the most representative hemisphere in each individual. Fibers originated from the VTAs that connect to M1 are shown in green, whereas fibers connecting to SMA and preSMA are shown in orange and blue, respectively. Patients are separated in responders (green bar) and non-responders (red bar) in descending order from the best to the worst clinical outcome from left to right. Independently from their clinical outcome, all patients showed thalamic connectivity to the seed areas.
Brainsci 10 00302 g004
Figure 5. Correlation between postoperative improvement of YGTSS after six months expressed in percentage and the mean of fibers per voxel located within the cortical seed regions. All regions of interest (ROIs, p < 0.05, R = 0.63) summed together including the preSMA, SMA and M1 are shown in (A), while the independent seed regions are shown separately in (B) (preSMA, p < 0.05, R = 0.61), (C) (SMA, p > 0.05, R = 0.25) and (D) (M1, p >0.05, R = 0.11). Data shown in this figure represent individual patient analysis (n = 7). For a segmented analysis of each hemisphere (n = 14), please consult Figure A2.
Figure 5. Correlation between postoperative improvement of YGTSS after six months expressed in percentage and the mean of fibers per voxel located within the cortical seed regions. All regions of interest (ROIs, p < 0.05, R = 0.63) summed together including the preSMA, SMA and M1 are shown in (A), while the independent seed regions are shown separately in (B) (preSMA, p < 0.05, R = 0.61), (C) (SMA, p > 0.05, R = 0.25) and (D) (M1, p >0.05, R = 0.11). Data shown in this figure represent individual patient analysis (n = 7). For a segmented analysis of each hemisphere (n = 14), please consult Figure A2.
Brainsci 10 00302 g005
Table 1. Demographic Table. Summarizes the clinical outcome of each patient six months after surgery, individual deep brain stimulation (DBS) parameters including the active contacts responsible for stimulation, comorbidities and medication. YGTSS—Yale Global Tic Severity Scale, OCD—obsessive compulsive disorder, THC—Tetrahydrocannabinol.
Table 1. Demographic Table. Summarizes the clinical outcome of each patient six months after surgery, individual deep brain stimulation (DBS) parameters including the active contacts responsible for stimulation, comorbidities and medication. YGTSS—Yale Global Tic Severity Scale, OCD—obsessive compulsive disorder, THC—Tetrahydrocannabinol.
Pat. NumberGenderAgeYGTSS PreopYGTSS PostopImprovement (%)Active Contacts
Pat. 1Male22782271.80 and 1
Pat. 2Male32894450.62 and 3
Pat. 3Female26954750.51 and 2
Pat. 4Male37814346.92 and 3
Pat. 5Male20874943.72 and 3
Pat. 6Male461007030.02 and 3
Pat. 7Male271008911.01 and 2
Pat. NumberFreq. HzPulse Width µsAmpl. VComorbidityMedication
Pat. 1130903.2NoneAripiprazole 10 mg, Citalopram 30 mg
Pat. 2130903.3NoneTiapride 800 mg, Quetiapine 150 mg
Pat. 3100903.0DepressionEscitalopram 10 mg, THC 5 mg, Risperidone 4 mg, Quetiapine 25 mg
Pat. 4110904.1Opiate dependence (in substitution)Methadone 28 mL
Pat. 51001505.2NoneNone
Pat. 6801505.0Depression, OCDAripiprazole 5 mg, Pipamperone 40 mg, Quetiapine 200 mg
Pat. 7803004.7Impulse control disorder, OCD, anxiety disorderClonazepam 1 mg, Pregabaline 300 mg, Risperidone 2 mg, THC 15 mg, Venlafaxine 225 mg

Share and Cite

MDPI and ACS Style

Andrade, P.; Heiden, P.; Hoevels, M.; Schlamann, M.; Baldermann, J.C.; Huys, D.; Visser-Vandewalle, V. Modulation of Fibers to Motor Cortex during Thalamic DBS in Tourette Patients Correlates with Tic Reduction. Brain Sci. 2020, 10, 302. https://doi.org/10.3390/brainsci10050302

AMA Style

Andrade P, Heiden P, Hoevels M, Schlamann M, Baldermann JC, Huys D, Visser-Vandewalle V. Modulation of Fibers to Motor Cortex during Thalamic DBS in Tourette Patients Correlates with Tic Reduction. Brain Sciences. 2020; 10(5):302. https://doi.org/10.3390/brainsci10050302

Chicago/Turabian Style

Andrade, Pablo, Petra Heiden, Moritz Hoevels, Marc Schlamann, Juan C. Baldermann, Daniel Huys, and Veerle Visser-Vandewalle. 2020. "Modulation of Fibers to Motor Cortex during Thalamic DBS in Tourette Patients Correlates with Tic Reduction" Brain Sciences 10, no. 5: 302. https://doi.org/10.3390/brainsci10050302

APA Style

Andrade, P., Heiden, P., Hoevels, M., Schlamann, M., Baldermann, J. C., Huys, D., & Visser-Vandewalle, V. (2020). Modulation of Fibers to Motor Cortex during Thalamic DBS in Tourette Patients Correlates with Tic Reduction. Brain Sciences, 10(5), 302. https://doi.org/10.3390/brainsci10050302

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