Abstract
Background Cervical dystonia is a movement disorder, characterised by involuntary head and neck muscle contractions. Although deep brain stimulation (DBS) of the globus pallidus internus (GPi) is an effective treatment option, motor outcomes can vary even when sufficient targeting accuracy is achieved. Increasing evidence supports a role of brainstem and cerebellum dysfunction in cervical dystonia pathogenesis.
Objective To determine whether morphometry of brainstem and dentate nuclei, and DBS stimulatory overlap with cerebello-thalamic tracts modelled from normative connectivity, were related to DBS clinical motor outcomes.
Methods 27 patients with idiopathic cervical dystonia underwent bilateral targeting of the GPi. and were separated into suboptimal and optimal motor outcome groups. Dentate nuclei and brainstem volumes were quantified in association with clinical outcomes. A brainstem shape analysis was conducted and used as a seed to assess connectivity from a normative structural connectome. Patient-specific electrodes were modelled to quantify stimulatory overlap with the GPi and proximity to cerebellothalamic tracts.
Results GPi implantation accuracy did not significantly differ between groups. Significantly reduced dentate nuclei and brainstem volumes were observed in patients with poorer clinical outcomes. Regional surface shape change of the brainstem was also observed in patients with poorer responses. Fibre tracking from this area intersected cerebellar, pallidal and cortical motor regions. Electrode field intersection with the non-decussating dentatorubrothalamic tract in the right, and in both hemispheres were also positively associated with clinical outcome.
Conclusions Variability in cerebellar and brainstem morphometry, and stimulation of non-decussating cerebello-thalamic pathways may contribute to the mediation of DBS motor outcomes.
Introduction
Cervical dystonia is an isolated focal movement disorder that causes involuntary muscular contractions of the neck, resulting in pain, cramping and abnormal posturing.(1) Deep brain stimulation (DBS) of the globus pallidus internus (GPi) has proven to be an effective therapeutic option for the treatment of medically refractory cervical dystonia.(2,3) Whilst the accuracy of GPi targeting in dystonia is the best predictor for a favourable outcome from DBS, in some cases, sufficient motor symptom improvements may still not be achieved.(4–6) The partial predictive ability of implantation accuracy and the poor understanding of disorder aetiology warrants the need for additional explanatory biomarkers for cervical dystonia. Identifying patients who have increased likelihoods of greater outcomes is paramount to reliably inform patient selection.
Once thought to be a sole disorder of the basal ganglia, research has elucidated dystonia to reflect a “circuitopathy”, indicating aberrant networks underpinning the pathophysiology.(7,8) Lesion-identification studies have identified spinal cord, basal ganglia, brainstem, cerebellum, and thalamic regions as most affected in patients with cervical dystonia,(9,10) with the brainstem and cerebellum as the most commonly affected regions in a recent study.(11) All lesion sites have been classified to ascertain connectivity with the cerebellum,(10) and surmounting evidence supports cerebellar involvement in dystonic sensorimotor pathophysiology.(7,8,12–14)
Brain-based imaging markers offer a potential avenue to explore variability in patients with cervical dystonia, in relation to therapeutic DBS outcomes. For example, functional MRI-derived connectivity has identified increased GPi-cerebellum and GPi-somatomotor cortex associations as indicative of greater and poorer clinical outcomes, respectively.(10) Furthermore, recent work using functional and tractography-derived connectivity identified the same network pattern, encapsulating the cerebellum and somatosensory cortex, as “optimal” pathways for therapeutic DBS outcomes in patients with cervical dystonia.(15)
Assessing the structure of the cerebellum and brainstem through imaging offers an alternative approach to explore these insights in relation to DBS outcomes. Cortico-ponto-cerebellar and cerebello-thalamo-cortical tracts form the structural connections that interlink the cerebellum, brainstem, regions of the basal ganglia and sensorimotor cortices.(16–19) Evidence of brainstem and cerebellar structural abnormalities have been observed in patients with cervical dystonia relative to healthy controls.(20–23) Pre-operative brain volumes have been seldom explored for DBS efficacy in cervical dystonia. To date, one study observed increased cerebellar vermis grey volume to be associated with improved outcomes.(24) To the best of our knowledge, no research has assessed brainstem structure in relation to DBS outcomes in cervical dystonia.
Assessment of white matter tracts offer an additional approach to assess brain structure in relation to DBS outcomes. The dentatorubrothalamic tract (DRTT) is a tract that has gained increasing popularity for importance of DBS outcomes in movement disorders including Parkinson’s disease and essential tremor,(25–28) but has yet to be explored in cervical dystonia. Originating at the dentate nucleus of the cerebellum, classically defined decussating and more recently disclosed non-decussating tracts, convey cerebellar output to the thalamus, with onward projections to sensorimotor cortex.(29–31) Evidence for structural changes in cerebello-thalamo-cortical and specific DRTT pathways have been identified in DYT1/DYT6 and cervical dystonia, respectively.(32,33) Probing this pathway in the context of DBS modulation could present as a powerful biomarker for cervical dystonia pathophysiology, aiding patient selection and treatment prediction.
The present study aimed to address these questions by quantifying brainstem volume and brainstem shape morphometry to determine relationships with clinical DBS motor improvements. We hypothesised that patients with poorer motor outcomes would have volumetric reductions in comparison to patients with improved responses. In a subsequent analysis using a structural normative connectome, we sought to determine whether connectivity from brainstem areas found to be different between outcome groups would map to subcortical and cortical sensorimotor regions. In doing so, areas of change may occupy networks involved in cervical dystonia pathophysiology and aid potential therapeutic DBS mediation. Finally, we explored the role of the decussating and non-decussating DRTT’s by first quantifying volumes of cerebellar dentate nuclei (as the DRTT origination point) and then assessing the relationship between stimulation fields and active contact proximities relative to the tracts in relation to motor improvements. We hypothesised that patients with poorer motor outcomes from DBS would have reduced dentate nuclei volumes and that electric field overlap coverage and contact proximity would be greater and closer to both DRTT pathways, respectively.
Methods
Patient Cohort
A retrospective cohort of 27 patients diagnosed with idiopathic cervical dystonia and treated with bilateral DBS was recruited from the movement disorders clinics at the Walton Centre NHS Foundation Trust (Liverpool, UK). Inclusion criteria included no structural lesions identifiable on the diagnostic MRI protocol, including pre-operative 2D T2, 2D T2-FLAIR and 3D T1-weighted magnetic resonance imaging (MRI), a confirmatory post-operative CT scan and available pre-and-post-operative clinical motor scores.
The Toronto Western Spasmodic Torticollis Rating Scale (TWSTRS) was carried out by a neurologist to measure cervical dystonia motor symptoms pre-and-post DBS implantation.(34) The TWSTRS score (Tout) was used as the outcome variable in the current study and was calculated using the following formula: The Tout for each patient was obtained at the best improvement score between three and five years post-DBS implantation to provide a long-term clinically stabilised score,(35–37) excluding two patients who due to the recency of date of implantation, best two-year post-operative scores were used. Unlike movement disorders such as Parkinson’s disease, at shorter time periods, the response of dystonia to DBS is not stable, suggested to be related to motor-plasticity effects from pallidal DBS.(38) To create groups based on motor response to DBS, patients were assigned to a sub-optimal response group given a Tout of < 66% and an optimal-response group given a Tout of ≥ 66%.
DBS Implantation
All patients included in the current study underwent DBS surgery at the Walton Centre NHS Foundation Trust (Liverpool, UK) over a period of 11 years (2009 - 2020). Patients were implanted with bilateral leads in the GPi, excluding one patient who was implanted with bilateral leads targeting the ventral intermediate nucleus/zona incerta (VIM/ZI) for cervical dystonia and additional dystonic tremor. Different electrode models were implanted across the patient cohort, namely Medtronic 3387 (n = 5), St Judes Active Tip 6142/6145 (n = 17), Boston Scientific Vercise Directed (n = 3) and St Judes Directed 6172 (n = 1).
MRI Acquisition
3 Tesla 3D T1-weighted images were acquired across three different scanners. Acquisition using a Phillips Achieva system with an 8-channel SENSE head coil (Philips Medical Systems, Best, The Netherlands), of 1 mm3 isotropic MPRAGE (gadolinium contrast enhancement (CE) in 10 patients and no CE in five patients), 256 slices (for 14 patients) and 288 slices (for one patient) flip angle of 8°, TE = 0.004, TR = 0.009. For the GE Discovery MR750 system (GE healthcare, USA) with a 16-channel head coil, isotropic 1 mm3 in 1 patient and 1.4 mm3 in 1 patient CE-SPGR with 256 slices, flip angle of 12°, TE = .003, TR = .008. For the Siemens MAGNETOM Skyra (Siemens, Healthineers, Germany) with a 20-channel head coil, 256 slices at a resolution of 1 mm × 1 mm × 0.9 mm CE-SPGR, with a flip angle of 9°, TE = .003, TR = .011.
Image Processing
An overview of the image processing pipeline is presented schematically in figure 1 and is described in detail in the following sections. As some images were acquired with contrast, all T1-weighted MRI were processed with SynthSR (https://github.com/BBillot/SynthSR), a convolutional neural network used to balance contrast and provide a T1-weighted 1mm isotropic MPRAGE image for each patient.(39) Synthesis of the images was a necessary step in the pipeline to improve grey and white matter contrast (and subsequent segmentation) of the original images and provide cross-scanner harmonisation. A post-operative computed tomography (CT) image was also obtained for each patient following DBS implantation.
Brainstem and Cerebellar Morphometry
To obtain volumetric estimates of the brainstem, the synthesised images were processed using FreeSurfer (version 7.2.0, https://surfer.nmr.mgh.harvard.edu). Bayesian segmentation was performed in FreeSurfer to obtain grey matter volumes (mm3) for pons, superior cerebellar peduncle (SCP), medulla, midbrain, and total brainstem.(40) The Spatially Unbiased Infratentorial Toolbox (SUIT; version 3.4, https://github.com/jdiedrichsen/suit/) was used to calculate grey matter volumes (mm3) of left and right dentate nuclei of the cerebellum.(41) SUIT was used to first isolate the cerebellum from the rest of the brain using the synthesised T1-weighted MRI. Diffeomorphic anatomical registration through exponentiated lie algebra (DARTEL) was used to normalise the isolated cerebellar images into Montreal Neurological Institute (MNI) space.(42)
Brainstem Shape Analysis
To measure localised points of surface change of the brainstem, we used FSL-integrated registration and segmentation toolbox (FSL-FIRST; version 5.0; https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FIRST).(43) The brainstem was automatically segmented from the synthesised T1-weighted MRI images and vertex-wise comparisons were computed along the surface of the brainstem in MNI space. The contrast was focused on assessing where vertex-wise reductions may occur with regards to sub-optimal responders relative to optimal responders.
Electrode Modelling
Lead-DBS (version 2.6, https://www.lead-dbs.org) was used to model electrode trajectories for each patient.(44–46) 25 patients were included in lead modelling, with one patient being excluded due to poor CT acquisition and one patient being excluded due to implantation of the VIM/ZI for additional dystonic tremor. First, a two-stage linear co-registration of pre-operative MRI to post-operative CT was performed using advanced normalisation tools (ANTs; https://stnava.github.io/ANTs/).(47) Following this, ANTs SyN Diffeomorphic Mapping was used to normalise volumes to ICBM 2009b nonlinear asymmetric (“MNI”) space.(48) Brain shift correction was performed using a coarse mask. Electrode trajectories were pre-reconstructed for each patient manually within Lead-DBS and were inspected and refined by a movement disorder specialist. Patient-specific volume of tissue activated (VTA) were constructed using the FieldTrip-SimBio pipeline.(46)
DBS Analyses
To visualise the structural connectivity from the area of brainstem change, the FSL FIRST shape analysis mask was used to seed connectivity in the lead connectome mapper within Lead-DBS. The normative structural connectome MGH-USC HCP-32 was used to obtain a whole-brain fibre density map at a 0.5 mm resolution that permitted the analysis of the number of fibres connecting the brainstem ROI to all other voxels.(49)
The VTA-GPi overlap was calculated to assess implantation accuracy in the cohort. The DISTAL atlas was used to calculate overlap with the whole GPi, and the sensorimotor (posteroventral) region.(50) Here, one patient was removed due to targeting of the VIM/ZI for additional dystonic tremor (n = 25).
To assess the relationship between DBS and the cerebellar tracts, electric field distribution and contact proximity approaches were used. Three patients were excluded due to implantation of directional electrodes. Tract ROIs were obtained from the DBS tractography atlas.(51) Electric field distribution overlap was correlated with Tout scores for the left and right decussating and non-decussating DRTT’s. To assess electrode proximity to these tracts, active contact positions for each electrode (left and right) were correlated in relation to the distance (mm) from these tracts with Tout scores. For each patient, the distance was averaged across active contacts for each electrode.
Statistics
Brainstem and dentate nuclei volumes were corrected for age, sex, and estimated total intracranial volume (eTIV; estimated from FreeSurfer) using the residuals from linear regression.(52) Volumes were compared between sub-optimal and optimal groups using two-sample t-tests. Statistical thresholds were Bonferroni corrected for the number of multiple comparisons. Corrections for cerebellar hemispheres (p < 0.025) and the four segmented brainstem nuclei (p < 0.0125) were performed.
Brainstem surface change was computed based on differences between the sub-optimal and optimal groups. The Threshold-Free Cluster Enhancement (TFCE) test-statistic was applied in randomise with 1000 permutations to identify clusters of voxels showing significant differences between groups. A threshold of p < .05 was used to as the level of significance. Age and sex were used as confounds in the model.
For Lead-DBS-based analysis, statistics were computed using Lead-group.(53) Linear regressions were used to assess the statistical significance of the VTA intersection with the GPi with a threshold of p < .025 to correct for laterality and the electric field overlap and electrode proximity to pre-defined tracts with a corrected threshold of p < 0.0125.
Results
Patient Characteristics
For the surface and volumetric analyses, a total of 13 subjects were identified as optimal responders (six female) with 14 subjects identified as suboptimal responders (nine female). Optimal and sub-optimal responders did not show significant differences between age, eTIV or known disease duration (p ≥ .2). The demographic information for this cohort is presented in table 1.
Brainstem and Cerebellar Morphometry
Grey matter brainstem volumetry revealed significant differences in pons (t = 2.87, p = .009) and total brainstem (t = 2.73, p = .01) volume between outcome groups. SCP volume did not reach statistical significance following the correction for multiple comparisons (t = 2.26, p = .03). Medulla and midbrain volume differences did not reach statistical significance (t = 1.88, p = .07 and t = 1.78, p = .08, respectively). Left (t = -2.73, p = .0119) and right (t = -2.85, p = .008) dentate nuclei volumes were significantly different between groups. Results are presented in figure 2.
From the vertex-wise shape analysis, an area of left and right superior ventral and dorsal brainstem at the level of midbrain and pons showed inwards deflation in patients with a sub-optimal response, relative to an optimal response (figure 3A).
Brainstem Connectivity Mapping
Normative connectivity was seeded from the brainstem area showing localised deflation in patients with sub-optimal responses (figure 3A) to create a whole-brain fibre density map. The fibre density map reflects the streamline count and outlined termination at the cerebellum with intersection of the deep cerebellar nuclei and crus lobules I and II and VIIb (SUIT atlas; figure 3B, left)(41). Termination was also observed at cortical motor regions, particularly at medial premotor cortex, delineating 100% probability of intersection using the Juelich Histological Atlas, distributed with FSL (http://www.fmrib.ox.ac.uk/fsl/; figure 3B, right). Subcortically, intersection with the left and right pallidum at 100% probability was observed using the Harvard-Oxford Subcortical Structural Atlas, also distributed with FSL (figure 3C), and primarily traversed the posteroventral GPi (figure 3D).
Electrode Modelling, VTA Overlap and Tract Analysis
The placement of electrodes for the patient cohort is presented in figure 4A, showing optimal (n = 13) and sub-optimal (n = 12) responders.
To test for lead placement between response groups, the normalised GPi-VTA overlap was computed for each patient. No significant association with clinical outcome was observed for VTA overlap with the whole GPi (R = .27, p = .09) or sensorimotor (posteroventral) GPi (R = .28, p = .09).
Normalised electric field intersections identified a significant positive relationship of overlap with the non-decussating DRTT on the right side (R = .61, p = .001) and with both hemispheres (R = .5, p = .008) with clinical outcome (figure 4B). The decussating tracts showed no significant associations with clinical outcome. When assessing the relationship between active contact proximity to the non-decussating DRTT (suboptimal: right (M: 9.4 ± 2.9 mm), left (M: 11.2 ± 2.7 mm); optimal: right (M: 8.7 ± 2.1 mm), left (M: 9.9 ± 2.4 mm)) and decussating DRTT (suboptimal: right (M: 7.4 ± 2.7 mm), left (M: 8.4 ± 2.2 mm); optimal: right (M: 6.9 ± 2.2 mm), left (M: 7.4 ± 2.3 mm)), no significant relations with Tout were observed for either tract (p > .4), nor were there any differences between groups.
Discussion
In the present study, we identified structural correlates related to the post-operative clinical motor outcome from DBS in patients with cervical dystonia. First, volumetric reduction of left and right dentate nuclei of the cerebellum and brainstem nuclei were observed in patients who experienced sub-optimal alleviation of motor symptoms from DBS treatment at long-term time-points (≥ 2-years) post-implantation. In addition, we report surface level changes, namely regional inward deflation of the brainstem, in the sub-optimal response group relative to the optimal response group. When normative structural connectivity was mapped from the surface area showing inwards deflation, probabilistic streamlines connected with cerebellar, subcortical and cortical motor regions, which corroborate a literature-defined network in cervical dystonia.(10,15) Analysis of electric field overlap with cerebellar tracts revealed that the non-decussating DRTT was significantly correlated with clinical outcome. Importantly, differences in clinical outcomes between patient groups appeared to not be driven by the accuracy of electrode placement, further enforcing the need to identify biomarkers beyond surgical accuracy. Taken together, our findings indicate potential new pathophysiological biomarkers for cervical dystonia that may aid in the identification of patients more likely to respond to DBS.
Volumetric reduction of the total brainstem, pons, and SCP (uncorrected) in patients with sub-optimal motor outcomes correspond to literature findings of affected areas in patients with cervical dystonia relative to healthy controls.(22,23) Pons and SCP connectivity, inferred using diffusion tractography, has previously been reported to show altered microstructural integrity, relative to healthy controls.(20,21) Both the pons and SCP share connectivity with the cerebellum, through the cortico-ponto- and cerebello-thalamo-cortical pathways, respectively,(55,56) and are thus heavily involved in the neuroanatomical substrate of whole-brain motor pathways. The pons has been identified as an area associated with increased regional cerebral blood flow in patients with clinical improvements of cervical dystonia symptoms from GPi-DBS.(57) Furthermore, using functional MRI, increased activation of the pons accompanied by decreased activation of the sensorimotor cortex have been observed in patients with cervical dystonia when switching from optimal to non-optimal GPi-DBS settings.(58) Damage to the cerebellopontine angle, referring to the space between the pons and cerebellum, has been shown to elicit cervical dystonia in a single case study,(59) indicating a pathophysiological role of brainstem-cerebellar coupling in symptom generation.
Structural brain changes in focal dystonia are favoured within a compensatory framework in contrast to neurodegeneration.(60) Physiological cerebellar changes in focal dystonia have been debated to arise as part of the core pathology or as compensatory changes to basal ganglia dysfunction.(61,62) Purkinje cell loss, as observed in post-mortem patients with cervical dystonia,(12) is consistent with reduced cerebellar volumes identified in patients with poorer motor response to DBS. Purkinje cell loss has been proposed to result in diminished GABAergic olivo-cerebellar output to the thalamus,(63) resulting in greater impairment of motor networks, and in the current study, reduced efficacy of therapeutic targeting. In patients with cervical dystonia, reduced levels of GABA have been identified in the thalamus relative to healthy controls, supporting a hypothesis of deficient cerebellar outflow.(64) In rodent models, pharmacological dampening of purkinje cell spiking has been shown to elicit dystonic symptoms directly.(65) Furthermore, optogenetic dentate nucleus stimulation to dampen cerebello-thalamic excitability has been shown to directly reduce dystonic motor symptoms.(66)
Forming a structural pathway between the cerebellum, subcortical regions, and areas of the cortex, the DRTT poses as a novel therapeutic target for movement disorder DBS.(67) To the best of our knowledge, the DRTT in cervical dystonia has been explored in one previous study, showing altered markers of microstructure relative to healthy controls.(33) Namely, reduced fractional anisotropy of the left DRTT and reduced axial and mean diffusivity of the right DRTT were observed. Work addressing the non-decussating pathway was not identified. The strong associations of electric field overlap with the non-decussating DRTT in the current study present interesting findings in the context of DBS response. Aside from an evolutionary mechanism of ipsilateral cerebellar hemispheric capacity for bilateral motor control, diffusion tractography has shown that the decussating and non-decussating DRTT’s project onto anatomically distinct thalamic regions. Non-decussating tracts favour targeting of ventroposterior thalamic nuclei whereas decussating tracts favour targeting of ventrolateral thalamic nuclei.(31) The ventroposterior nuclear complex has affiliation with somatosensory functioning, connecting with cortical sensorimotor regions.(68,69) It could be theorised that modulating non-decussating DRTT could elicit corrective effects on abnormal cerebellar outflow with consequent upstream normalisation of aberrant sensorimotor processing.(70– 72) The lack of significance for the left non-decussating DRTT is a novel finding and may be related to the dominant side of symptomology present in the patient cohort. However, as this data was not available in the current study, this question could not be explored further.
The observed connectivity pattern derived from the morphometric brainstem surface change supports a potential role of outcome mediation in a network-based manner. Intersection of the cerebellum and pallidum verifies lesion sites identified in patients with secondary cervical dystonia.(9–11) Posteroventral GPi intersection supports the brainstems involvement in mediating therapeutic networks, given that this area of the pallidum is widely regarded as the DBS hotspot for cervical dystonia, and occupies somatotopic cervical territory.(15,73,74) The observed cortical motor fibre intersection corroborates the somatotopic head and neck areas, as the terminus of white matter tracts associated with optimal clinical motor symptom reduction, in cervical dystonia.(15) Structural changes to primary and pre-motor regions,(22,75–77) and decreased resting-state functional MRI activity between cerebellar and primary motor cortex,(78) have been identified in patients with cervical dystonia relative to healthy controls, supporting a role of motor circuit changes in the disorder’s pathophysiology.
Limitations
We acknowledge the limitations in the present study. Firstly, the use of SynthSR as a valid contrast balancing tool has yet to be independently validated. It is unclear how much variability is introduced in the image following synthesis; however, its application has shown reliable volumetric estimations from anisotropic 2D images, which are less desirable than the 3D isotropic images used in the present study. (Iglesias et al, 2021) Also, due to scans being acquired across three separate scanners, systematic variability in volumetric estimates cannot be ruled out entirely. SynthSR was favoured for both these cases by first, providing a degree of harmonisation across scanners, and second, showing high visual consistency pre-and-post synthesis following manual inspection.
The restriction to T1-weighted images due to the routine clinical care MRI protocol was a limitation in this study. The addition of 3D T2 and 3D T2-FLAIR imaging would have been beneficial to improve segmentation accuracy for volumetry and electrode modelling.
In addition, normative connectivity is not reflective of patient-specific connectivity and may not capture disorder-related microstructural tract changes associated with cervical dystonia.(32,33) Despite this, normative and patient-specific connectivity utility has been shown to be comparable in Parkinson’s disease.(79)
Conclusions
Taken together, the results from this study support cervical dystonia as a network disorder and provide evidence for associations between the structure of key motor regions, and structural pathways with clinical motor outcomes. We show that dentate and brainstem nuclei morphometry may represent markers that determine the efficacy of DBS, outside of GPi targeting accuracy. Furthermore, we present the potential importance of the non-decussating DRTT in the pathophysiology and therapeutic targeting of cervical dystonia.
Data Availability
Data is available upon reasonable request to the corresponding author.
Ethical compliance statement
This work was sponsored by the University of Liverpool and the Walton Centre NHS Trust. Approval for this work was obtained by HRA and Health and Care Research Wales (HRCW; REC reference number: 22/PR/1326).
Informed patient consent was not necessary for this work.
Acknowledgements
Mrs Bonnie O’Sullivan, Digital Systems PACS Manager at Radiology Directorate, The Walton Centre NHS Trust
Neurologists of Movement Disorders DBS Team, The Walton Centre NHS Trust: Dr Dinesh Damodaran, Dr Sundus Alusi, Dr Jay Panicker, Dr Michael Bonello