Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Target Identification for Stereotactic Thalamotomy Using Diffusion Tractography

  • Zsigmond Tamás Kincses ,

    kincsesz@nepsy.szote.u-szeged.hu

    Affiliations Department of Neurology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary, International Clinical Research Center, St. Anne's University Hospital Brno, Brno, Czech Republic

  • Nikoletta Szabó,

    Affiliation Department of Neurology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • István Valálik,

    Affiliation Department of Neurosurgery, St. John's Hospital, Budapest, Hungary

  • Zsolt Kopniczky,

    Affiliation Department of Neurosurgery, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • Lívia Dézsi,

    Affiliation Department of Neurology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • Péter Klivényi,

    Affiliation Department of Neurology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • Mark Jenkinson,

    Affiliation FMRIB Centre, Department of Clinical Neurology, John Radcliffe Hospital, University of Oxford, Oxford, United Kingdom

  • András Király,

    Affiliation Department of Neurology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • Magor Babos,

    Affiliation Euromedic Diagnostic Hungary Ltd, Szeged, Hungary

  • Erika Vörös,

    Affiliation Department of Radiology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • Pál Barzó,

    Affiliation Department of Neurosurgery, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

  • László Vécsei

    Affiliation Department of Neurology, Albert Szent-György Clinical Center, University of Szeged, Szeged, Hungary

Abstract

Background

Stereotactic targets for thalamotomy are usually derived from population-based coordinates. Individual anatomy is used only to scale the coordinates based on the location of some internal guide points. While on conventional MR imaging the thalamic nuclei are indistinguishable, recently it has become possible to identify individual thalamic nuclei using different connectivity profiles, as defined by MR diffusion tractography.

Methodology and Principal Findings

Here we investigated the inter-individual variation of the location of target nuclei for thalamotomy: the putative ventralis oralis posterior (Vop) and the ventral intermedius (Vim) nucleus as defined by probabilistic tractography. We showed that the mean inter-individual distance of the peak Vop location is 7.33 mm and 7.42 mm for Vim. The mean overlap between individual Vop nuclei was 40.2% and it was 31.8% for Vim nuclei. As a proof of concept, we also present a patient who underwent Vop thalamotomy for untreatable tremor caused by traumatic brain injury and another patient who underwent Vim thalamotomy for essential tremor. The probabilistic tractography indicated that the successful tremor control was achieved with lesions in the Vop and Vim respectively.

Conclusions

Our data call attention to the need for a better appreciation of the individual anatomy when planning stereotactic functional neurosurgery.

Introduction

Thalamotomy was introduced in the treatment of tremor by Hassler in 1954 [1], the selective stereotactic lesioning of the ventralis intermedius nucleus of the thalamus (Vim) was described by Narabayashi [2] and later the electrical stimulation of Vim in 1987 by Benabid [3]. The precise targeting within the brain is of crucial importance for successful surgical intervention. Targeting the desired thalamic nucleus is usually carried out by using stereotactic coordinates specified in relation to a point on the anterior commissure – posterior commissure (AC-PC) line [4], [5], [6]. Other methods try to establish population-based stereotactic coordinates of the target nuclei based on postmortem histological data or intraoperative stimulation techniques [7], [8], [9]. While achieving reasonable results, these methods neglect individual anatomical variations.

Since the identification of thalamic nuclei on conventional imaging modalities is difficult, several novel approaches were proposed to aid the visualisation of the functionally important thalamic nuclei. It is now possible to segregate the major thalamic structures using MR relaxometry, even in a clinically acceptable time [10], [11]. Recently, probabilistic tractography was successfully used to investigate the connectivity profile of two major thalamic target nuclei for functional neurosurgery: ventralis intermedius (Vim) and ventralis oralis posterior (Vop) [12]. In addition, it is possible to segment the thalamic nuclei based on connectivity patterns defined by MR diffusion tractography [13], [14]. This approach has been suggested by others [12], [15], [16], but the benefit of such a method has not been evaluated systematically so far.

In the current investigation we calculated the interindividual variability of the target thalamic nuclei, (Vop and Vim) location in healthy controls. We also present the retrospective DTI identification of the target thalamic nuclei in two patients that underwent stereotactic Vim and Vop thalamotomy.

Methods

Participants

Nine healthy individuals, with no history of neurological or psychiatric diseases were included in the study (mean age: 28.3±7.09, male: 3). Furthermore, two patients were also included, who underwent stereotactic thalamotomy.

Ethics

The study was approved by the Ethics Committee of University of Szeged (authority number: 87/2009), and all subjects provided written consent.

Image acquisition

Imaging was carried out using a 1.5 T GE Signa Excite scanner. High resolution T1 weighted images (3D IR-FSPGR: TR/TE/TI: 10.3/4.2/450ms, flip angle: 15°, ASSET: 2, FOV: 25*25 cm, matrix: 256*256, slice thickness: 1mm) and diffusion-weighted images (DTI: TR/TE: 16000/93.8ms, flip angle: 90°, FOV: 23*23cm, matrix: 96*96, slice thickness: 2.4mm, ASSET: 2, NEX: 2) were acquired. Diffusion weighting was performed along 60 independent directions, with a b-value of 1000 s/mm2. Six reference images, with no diffusion weighting were also obtained. The preoperative scans of the patients were obtained using identical acquisition sequences to those used for control subjects. Postoperative scans were carried out three months after the surgery. High-resolution-T1 weighted images (with parameters identical to the preoperative ones) and sagittal 3D FLAIR images (3D FLAIR: TR/TE/TI: 6000/134.6/1839ms, flip angle: 90°, FOV: 23*23cm, matrix: 256*256, slice thickness: 2mm) were acquired to localise the thalamotomic lesion.

Image processing

Data were processed using the tools from the FMRIB Software Library (FSL, version 5.0; Oxford Centre for Functional MRI of the Brain (FMRIB), UK; www.fmrib.ox.ac.uk/fsl [17]) according to the method reported by Behrens and colleagues [14]. Initially the raw diffusion data was corrected for eddy-currents and motion artefacts with FMRIB's Diffusion Toolbox (FDT [18]). The images were then skull stripped using the Brain Extraction Tool (BET [19]) and the diffusion-weighted images were registered to the high-resolution T1-weighted image with a 6 degree-of-freedom linear registration using the FMRIB Linear Image Registration Tool (FLIRT [20]). Probability distributions of fiber orientation were estimated for each brain voxel in the acquired diffusion space using a multi-fiber extension [21] of the probabilistic tractography available in FDT [18].

Probabilistic diffusion tractography can estimate the pathways that originate at any given seed voxel, as well as the probability that such a pathway will pass through any other voxel in the brain [18]. Binary masks of the thalamus and the cortical targets were drawn manually for each subject. Probabilistic multi-fiber diffusion tractography was initiated from every voxel inside the thalamic mask: 5000 times for each voxel. Counters were increased every time an individual streamline reached the cortical target region. Hence, the values stored in thalamic voxels in one of the resulting images represent the probability of those voxels being connected to the particular target cortical mask assigned to that image [14].

To investigate the inter-subject variability of the target nuclei, a specific distance reserving registration method was used: the high-resolution T1-weighted images were registered to the standard MNI brain with a 6 degree-of-freedom transformation (rotations and translations only) that was based on landmark points in a way that kept the position of the anterior commissure fixed, while aligning the ACPC line and the mid-sagittal plane between images.

The Euclidean distance, given by the equation:

was calculated for each subject between points representing the thalamic voxel with peak connection probability to the premotor cortex (for Vop) and to the motor cortex (for Vim) and the standard space voxel with peak connectivity as indicated by the Oxford Thalamic Connectivity Map [13]. The mean pair-wise distance was also calculated: the distance between the voxel with peak connection probability to the premotor (for Vop) and to the motor (for Vim) cortex was calculated between all possible pair of subjects and this inter-subject distance was averaged for Vop and Vim separately.

The connectivity images for the thalamus were thresholded at 10% of individual connectivity maximum and binarised in order to create masks of the thalamic regions connected to the premotor or motor cortex. To assess the similarity of positions the overlap between each pair of masks (between subjects) were calculated according to the method proposed by Crum and colleagues [22]. Overlap was measured by the Tanimoto Coefficient (TC), which is defined as the ratio of the number of voxels in the intersection of the two regions to the number of voxels in the union:

We also took further steps to exclude possible registration bias between the DTI and T1-weighted images that might originate from the EPI distortions. We manually compared the high-resolution T1-weighted image and the first volume of the diffusion data (no diffusion gradient applied) registered to the T1-weighted image for each subjects. The positions of the following landmarks were compared: AC, PC, the maximal width of the third ventricle at the level AC and the highest point of the corpus callosum in the midsagittal line.

Results

Inter-subject variability of the position of thalamic nuclei

As found with the manual comparison the shape and the position of the thalamus was not affected by EPI distortions. The largest misregistrations were found along the anterior-posterior axis, but even that was minimal (for AC and for PC 0.6±0.5mm). The width of the third ventricle was not different, and the position of the highest point of the corpus callosum differed only in case of a single subject with 1 mm.

The premotor thalamus was consistently localised in all subjects, however, since the individual brain shapes differed, the exact location of the Vop and Vim nucleus varied substantially across subjects (Fig 1). The mean distance of the peak connection probability of Vop from the peak probability indicated by the Oxford Thalamic Connectivity map was 5.08mm. The equivalent distance was 6.26mm for the Vim. The mean pair-wise inter-subject distance of the peak connectivity voxel of Vop was 7.33±3.37mm (range: 0–14.56mm) and 7.42±3.35 mm (range: 2–14.28mm) for Vim. The mean distance of the tractography defined Vop and Vim coordinates from the Stereotactic target point defined by Hyam's method [12] was 7.19±4.36mm (range: 2.45–14.89mm) for Vim and 9.58±4.82mm (range: 3.0–17.12mm) for Vop.

thumbnail
Figure 1. Spatial localisation of peak Vop and Vim voxels as defined by probabilistic tractography.

Individual images were realigned to standard space with a 6 DOF transformation in a way to match the location of the AC and align the AC-PC line and the midsagittal plane. The registered peak connectivity voxels for the primary motor (Vim - in red) and for the premotor cortex (Vop – in blue) are shown on three orthogonal slices.

https://doi.org/10.1371/journal.pone.0029969.g001

The mean pair-wise overlap for Vop as calculated by the Tannimoto Coefficient was 40.2% (range: 15.5%–66.2%) and 31.8% (range: 3.2%–66.2%) for Vim.

Case #1: Vop thalamotomy for posttraumatic tremor

A left-handed, 32 year-old male, after having a motorbike accident that caused brainstem bleeding (at the level of the red nucleus; Fig 2), developed a tremor having resting, postural, as well as intention components, bradykinesia, rigidity and a slight paresis of the left arm. Seven years after the accident, because of ineffective tremor control by various medication regimes, functional neurosurgery was offered. Because the patient has not accepted deep-brain-stimulation implantation stereotactic thalamotomy was carried out on the right side in the region of Vop/Vim border with target coordinates according to Guiot's method [5], [23] that is, given the AC-PC length of 24.9mm, was 8mm anterior of PC and 12.7mm lateral from the mid-sagittal line on the AC-PC plane. During the standard surgical procedure the final position of the lesion along the selected trajectory was adjusted according to the tremor control of the intraoperative stimulation. After surgery the patient had an immediate relief from his left upper extremity tremor (his score on the part A of the Fahn-Tolosa-Marin tremor rating scale (FMT) reduced form 16 to 4 points). Neither the muscle strength, nor the bradykinesia improved. At the three-month checkup these findings had not changed for the left upper extremity, however a small amplitude flexion-extension tremor appeared in the left toes.

thumbnail
Figure 2. Preoperative FLAIR images show the right perirubral mesencephalic traumatic haemorrhage.

Images are presented in radiological convention.

https://doi.org/10.1371/journal.pone.0029969.g002

As shown by the connectivity-based segmentation of the thalamus, regions connected with highest probability to the prefrontal, premotor, motor and sensory, parietal, occipital and temporal cortices were clearly delineated (Fig 3). The operative lesion of the thalamus, identified on the postoperative FLAIR images (registered to the preoperative T1-weighted structural image with 6 DOF linear registration) was situated in the region of the thalamus connected with the highest probability to the premotor cortex (Fig 3).

thumbnail
Figure 3. Connectivity-based segmentation of the right thalamus of the patient operated on for tremor caused by traumatic mesencephalic bleeding.

Results are overlaid on the postoperative FLAIR images. The bright area represents the Stereotactic thalamotomy. The key to the colour coding of the thalamic segmentation is shown in the inset. The colour in each thalamic voxel represents the colour of the cortical area that has the highest connection probability to that voxel.

https://doi.org/10.1371/journal.pone.0029969.g003

Probabilistic multi-fiber tractography initiated from the thalamic lesion showed, on one hand, connection to the medial premotor - prefrontal cortex and, on the other hand, to the brain stem particularly to the rubral/perirubral mesencephalic posthaemorrhagic lesion and some fibres travelling further on to the cerebellum (Fig 4).

thumbnail
Figure 4. Connectivity of the thalamic lesion, evaluated using multi-fiber probabilistic tractography.

The tracts reach the premotor cortex, on one hand, and the nucleus ruber and the perirubral mesencephalon, on the other hand, with some travelling further to the cerebellum. The image is thresholded at 500 particles (10%).

https://doi.org/10.1371/journal.pone.0029969.g004

Case #2: Vim thalamotomy for essential tremor

A right-handed, 50-year-old women had developed tremor ten years earlier. The tremor was most prominent in the hands but also appeared in the head, having postural and action components, and was also present when walking. Her family history was unremarkable. Medication regimes that had been tried over the years were ineffective, and the tremor severely affected the quality of life of the patient. Stereotactic Vim thalamotomy was carried out with target coordinates according to Guiot [5], [23] that correspond well with the coordinates defined by Hyam's method [12]. After surgery the tremor of the right hand improved significantly and on the three-month control examination this effect was found to be enduring (FMT score reduced from 17 to 11 points). No significant change was detected in the head-tremor or in the tremor of the left hand.

The postoperative T1-weighted image was registered to the preoperative T1-weighted image in order to identify the position of the lesion in relation to the diffusion–based connectivity results. The connectivity-based segmentation of the thalamus indicated that the surgical lesion was in the putative Vim nucleus, the region that showed the highest probability of connection to the primary motor cortex (Fig 5). Postoperative probabilistic tractography from the lesion showed that the lesioned thalamus was connected to the primary motor cortex, mainly to the medial structures, and downstream to the ipsilateral cerebellar hemisphere (Fig 6).

thumbnail
Figure 5. Connectivity-based segmentation of the left thalamus of the patient having essential tremor.

Results are overlaid on the postoperative T1-weighted image. White voxels represent the location of the Stereotactic lesion. The colour coding of the thalamic segmentation is the same as that using in Figure 3. The lesion is in the area of the thalamus connected most strongly to the primary motor cortex: the putative Vim.

https://doi.org/10.1371/journal.pone.0029969.g005

thumbnail
Figure 6. Connectivity of the thalamic lesion of the essential tremor patient, evaluated using multi-fiber probabilistic tractography.

The tracts reach the primary motor cortex, on one hand, and travelling to the cerebellum on the other hand. The image is thresholded at 250 particles (5%).

https://doi.org/10.1371/journal.pone.0029969.g006

Discussion

In the current study we investigated the inter-individual variability of the position of the thalamotomic target nuclei Vim and Vop as defined by probabilistic diffusion tractography. We found that, compared to the size of the nucleus, the variability of its position is substantial. Our results call attention to the importance of defining the thalamotomic lesion target individually. We also showed that the lesion in two successful thalamotomic cases was in fact in the target nuclei as defined by tractography based segmentation.

Despite the accurate stereotactic atlases [7], [24] that guide the standard methods for targeting, significant intersubject variability still exists [7]. The differences in gross thalamic morophology is already conspicuous: the medio-lateral aspect as defined by the position of the internal capsule, as well as the height of the thalamus vary between subjects [7]. The position of the pulvinar system varies even between subjects with the same AC-PC distance [7]. Moreover, as emphasized by Morel, these variations are not even homogenous within the thalamus, particularly along the mediolateral axis [7]. Our results, using a different modality, also indicate significant variability of the position of the Vop and Vim nuclei, especially when the size of the nucleus and the surgical lesion is considered. The distance between peak-probability Vop voxels was 7 mm on average and is similar for Vim, but in certain cases the distance was as much as one and a half centimetres. The overlap between Vops on average was around only 40% and varies from 15% to 60%. Similar overlap values were found for Vim. It is also important to emphasize that in the current variability estimation only young, healthy individuals were included. Brain pathology may further increase the variability.

One possible caveat of our analysis that cannot be neglected is the misregistration when aligning the high resolution T1-weighted and diffusion weighted images. The diffusion-weighted images suffer significant distortions because of susceptibility artefacts, however our analysis indicated that this distortions was minimal in the region of thalamus.

The high inter-subject variability should be accounted for during surgery. The first approach is to standardise the position of the lesion according to internal reference distances such as the length of the AC-PC line, the height of the thalamus and the width of the third ventricle [5], [24]. Instead of using the histology-based stereotactic atlas derived from a limited number of brains, recently a population based-probabilistic functional atlas was developed from data gathered from pre-, intra- and postoperative neuroimaging and electrophysiological investigations [8], [9]. The extension of this approach is the physiological control of lesion depth during awaken surgery by monitoring patients complains while trying different stimulation positions.

It follows from the discussion above that, ideally, individual anatomy should be respected when performing stereotactic thalamotomy. Recently, Yamada and colleagues suggested a diffusion-based method for identification of the ventral thalamic nucleus [15]. Their approach incorporates tracking of the cerebello-thalamo-cortical and the spino-thalamic tracts that cross the thalamus, and the identification of anatomical landmarks on short tau inversion recovery (STIR) and fractional anisotropy images [15]. Sedrak and co-workers used the fractional anisotropy weighted color-coded diffusion direction maps to manually segment the thalamus from the surrounding white matter and identify the ventral intermedius nucleus [16] to define target for deep-brain-stimulation (DBS). The connectivity of effective electrode location investigated by tractography showed connections to the sensory-motor cortex and to the cerebellum. Compared to these deterministic approaches our method offers further advantages by utilising probabilistic tractography [14], [18]. This method is able to give statistical evaluation of connection probabilities to cortical targets for every thalamic voxel that - by appropriate selection of the cortical targets – are then able to delineate thalamic regions that correspond well to thalamic nuclei as defined by histology [14], [18], [25]. This approach shows a high intra-subject reproducibility in terms of repeated analysis as well as repeated data acquisitions [25].

Within the ventral thalamic nuclear complex the ventralis oralis anterior/posterior (Voa/Vop) and the ventralis intermedius nucleus (Vim) can be delineated. The true connectivity of these nuclei is debated, and it is hard to mach different nomenclatures [26], but as a rule of thumb it can be said that Vim is the cerebellar receiving area while Vop is the pallidal receiving area of the thalamus [27], [28]. More importantly, from our point of view, cortical connectivity of the putative Vim and Vop was recently investigated using diffusion imaging [12]. It was found that Vop seeded tractograms connected to the supplementary motor area (SMA) and to the dorsolateral prefrontal cortex (DLPFC) with higher probability than Vim-seeded tractograms; while the Vim was more likely to connect to the primary motor cortex (M1) and contralateral cerebellum, although significant cerebellar connections were found for Vim as well as for Vop [12]. Here we provide a confirmation of Hyam's pioneering analysis by showing that the Vop and Vim coordinates selected for thalamotomy were indeed effective. The connectivity profile of our patients' thalamotomic lesions were strikingly similar to those described as Vop and Vim connectivity profiles by Hyam [12].

Lesioning or stimulation of the ventral thalamic nuclei is a frequently considered option for the treatment of tremor of various causes including post-traumatic Holmes tremor [29], [30], [31], [32], [33], [34]. Within the ventral nuclei, Vim is the optimal target in essential tremor [35], [36]. It is less and less frequently used in the treatment of Parkinsonian tremor because of the poor benefits for bradykinesia and rigidity [37]. In multiple sclerosis, beside targeting Vim some authors prefer to target Vop suggesting better control of the ataxic component of the tremor [4], [38]. In tremor cases developed after injury of the superior cerebellar peduncules or the rubrocerebellar connections (e.g. traumatic brain injury, tumour) tremor control usually is difficult. Reports about the stimulation/lesion of the Vim [32], [33], [39] and Voa/Vop [40] indicated positive outcome. Foote and colleagues also reported three posttraumatic tremor cases and one MS tremor case that responded to dual Vim-Voa/Vop stimulation [41]. Based on the connectivity profile, the thalamic lesion in our patient was in the Vop. This lesion resulted in excellent tremor control, however the rigidity and the bradykinesia of the affected extremity did not improve significantly. It also has to be mentioned that probabilistic tractography was recently also successfully used in defining connectivity of the pedunculo-pontine nucleus, a novel target for neurosurgical treatment of Parkinson's disease [42].

Conclusions and limitations

Atlas-based targeting for stereotactic functional neurosurgery is a frequently considered option for the treatment of movement disorders. While the overall outcome is reported to be positive, the failures and side effects might be related to the unaccounted inter-individual variability of the target localisation. To overcome such limitations we suggest incorporating information about individual anatomy from diffusion imaging, using connectivity-based segmentations of the thalamic nuclei. However, it has to be pointed out that in the two cases we reported in current analysis we did not plan the targeting with tractography based segmentation a priori of the surgery. Extension of this study to compare the result of the stereotaxic neurosurgery with and without tragtography based targeting is crucial. Further improvements might be expected if additional imaging modalities (e.g. relaxometry, fMRI) could also be utilised.

Author Contributions

Conceived and designed the experiments: ZTK EV PB LV. Performed the experiments: NS LD ZTK PK IV ZK AK MB. Analyzed the data: ZTK NS AK MJ. Contributed reagents/materials/analysis tools: MJ. Wrote the paper: ZTK NS LV. Carried out surgical intervention: ZK IV.

References

  1. 1. Hassler R, Riechert T (1954) [Indications and localization of stereotactic brain operations.]. Nervenarzt 25: 441–447.
  2. 2. Narabayashi H, Ohye C (1980) Importance of microstereoencephalotomy for tremor alleviation. Appl Neurophysiol 43: 222–227.
  3. 3. Benabid AL, Pollak P, Louveau A, Henry S, de Rougemont J (1987) Combined (thalamotomy and stimulation) stereotactic surgery of the VIM thalamic nucleus for bilateral Parkinson disease. Appl Neurophysiol 50: 344–346.
  4. 4. Bittar RG, Hyam J, Nandi D, Wang S, Liu X, et al. (2005) Thalamotomy versus thalamic stimulation for multiple sclerosis tremor. J Clin Neurosci 12: 638–642.
  5. 5. Dormont D, Cornu P, Pidoux B, Bonnet AM, Biondi A, et al. (1997) Chronic thalamic stimulation with three-dimensional MR stereotactic guidance. AJNR Am J Neuroradiol 18: 1093–1107.
  6. 6. Alusi SH, Aziz TZ, Glickman S, Jahanshahi M, Stein JF, et al. (2001) Stereotactic lesional surgery for the treatment of tremor in multiple sclerosis: a prospective case-controlled study. Brain 124: 1576–1589.
  7. 7. Morel A, Magnin M, Jeanmonod D (1997) Multiarchitectonic and stereotactic atlas of the human thalamus. J Comp Neurol 387: 588–630.
  8. 8. Nowinski WL, Belov D, Thirunavuukarasuu A, Benabid AL (2005) A probabilistic functional atlas of the VIM nucleus constructed from pre-, intra- and postoperative electrophysiological and neuroimaging data acquired during the surgical treatment of Parkinson's disease patients. Stereotact Funct Neurosurg 83: 190–196.
  9. 9. Nowinski WL, Belov D, Benabid AL (2002) A community-centric internet portal for stereotactic and functional neurosurgery with a probabilistic functional atlas. Stereotact Funct Neurosurg 79: 1–12.
  10. 10. Deoni SC, Rutt BK, Parrent AG, Peters TM (2007) Segmentation of thalamic nuclei using a modified k-means clustering algorithm and high-resolution quantitative magnetic resonance imaging at 1.5 T. Neuroimage 34: 117–126.
  11. 11. Deoni SC, Josseau MJ, Rutt BK, Peters TM (2005) Visualization of thalamic nuclei on high resolution, multi-averaged T1 and T2 maps acquired at 1.5 T. Hum Brain Mapp 25: 353–359.
  12. 12. Hyam JA, Owen SL, Kringelbach ML, Phd NJ, Stein JF, et al. (2011) Contrasting Connectivity of the Vim and Vop Nuclei of the Motor Thalamus Demonstrated by Probabilistic Tractography. Neurosurgery.
  13. 13. Johansen-Berg H, Behrens TE, Sillery E, Ciccarelli O, Thompson AJ, et al. (2005) Functional-anatomical validation and individual variation of diffusion tractography-based segmentation of the human thalamus. Cereb Cortex 15: 31–39.
  14. 14. Behrens TE, Johansen-Berg H, Woolrich MW, Smith SM, Wheeler-Kingshott CA, et al. (2003) Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat Neurosci 6: 750–757.
  15. 15. Yamada K, Akazawa K, Yuen S, Goto M, Matsushima S, et al. (2010) MR imaging of ventral thalamic nuclei. AJNR Am J Neuroradiol 31: 732–735.
  16. 16. Sedrak M, Gorgulho A, Frew A, Behnke E, Desalles A, et al. (2011) Diffusion Tensor Imaging (DTI) and Colored Fractional Anisotropy (FA) Mapping of the Ventralis Intermedius (VIM) Nucleus of the Thalamus. Neurosurgery.
  17. 17. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, et al. (2004) Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23: Suppl 1S208–219.
  18. 18. Behrens TE, Woolrich MW, Jenkinson M, Johansen-Berg H, Nunes RG, et al. (2003) Characterization and propagation of uncertainty in diffusion-weighted MR imaging. Magn Reson Med 50: 1077–1088.
  19. 19. Smith SM (2002) Fast robust automated brain extraction. Hum Brain Mapp 17: 143–155.
  20. 20. Jenkinson M, Smith S (2001) A global optimisation method for robust affine registration of brain images. Med Image Anal 5: 143–156.
  21. 21. Behrens TE, Berg HJ, Jbabdi S, Rushworth MF, Woolrich MW (2007) Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? Neuroimage 34: 144–155.
  22. 22. Crum WR, Camara O, Hill DL (2006) Generalized overlap measures for evaluation and validation in medical image analysis. IEEE Trans Med Imaging 25: 1451–1461.
  23. 23. Guiot G, Arfel G, Derome P, Kahn A (1968) [Neurophysiologic control procedures for sterotaxic thalamotomy]. Neurochirurgie 14: 553–566.
  24. 24. Gallay MN, Jeanmonod D, Liu J, Morel A (2008) Human pallidothalamic and cerebellothalamic tracts: anatomical basis for functional stereotactic neurosurgery. Brain Struct Funct 212: 443–463.
  25. 25. Traynor C, Heckemann RA, Hammers A, O'Muircheartaigh J, Crum WR, et al. (2010) Reproducibility of thalamic segmentation based on probabilistic tractography. Neuroimage.
  26. 26. Hamani C, Dostrovsky JO, Lozano AM (2006) The motor thalamus in neurosurgery. Neurosurgery 58: 146–158; discussion 146–158.
  27. 27. Macchi G, Jones EG (1997) Toward an agreement on terminology of nuclear and subnuclear divisions of the motor thalamus. J Neurosurg 86: 670–685.
  28. 28. Molnar GF, Pilliar A, Lozano AM, Dostrovsky JO (2005) Differences in neuronal firing rates in pallidal and cerebellar receiving areas of thalamus in patients with Parkinson's disease, essential tremor, and pain. J Neurophysiol 93: 3094–3101.
  29. 29. Chen H, Zhuang P, Miao SH, Yuan G, Zhang YQ, et al. (2010) Neuronal firing in the ventrolateral thalamus of patients with Parkinson's disease differs from that with essential tremor. Chin Med J (Engl) 123: 695–701.
  30. 30. Franzini A, Cordella R, Messina G, Marras CE, Romito LM, et al. (2011) Deep brain stimulation for movement disorders.Considerations on 276 consecutive patients. J Neural Transm.
  31. 31. Hooper AK, Okun MS, Foote KD, Fernandez HH, Jacobson C, et al. (2008) Clinical cases where lesion therapy was chosen over deep brain stimulation. Stereotact Funct Neurosurg 86: 147–152.
  32. 32. Kim MC, Son BC, Miyagi Y, Kang JK (2002) Vim thalamotomy for Holmes' tremor secondary to midbrain tumour. J Neurol Neurosurg Psychiatry 73: 453–455.
  33. 33. Sanborn MR, Danish SF, Ranalli NJ, Grady MS, Jaggi JL, et al. (2009) Thalamic deep brain stimulation for midbrain tremor secondary to cystic degeneration of the brainstem. Stereotact Funct Neurosurg 87: 128–133.
  34. 34. Blomstedt P, Hariz GM, Hariz MI, Koskinen LO (2007) Thalamic deep brain stimulation in the treatment of essential tremor: a long-term follow-up. Br J Neurosurg 21: 504–509.
  35. 35. Flora ED, Perera CL, Cameron AL, Maddern GJ (2010) Deep brain stimulation for essential tremor: a systematic review. Mov Disord 25: 1550–1559.
  36. 36. Zesiewicz TA, Elble R, Louis ED, Hauser RA, Sullivan KL, et al. (2005) Practice parameter: therapies for essential tremor: report of the Quality Standards Subcommittee of the American Academy of Neurology. Neurology 64: 2008–2020.
  37. 37. Pahwa R, Lyons KE, Wilkinson SB, Simpson RK Jr, Ondo WG, et al. (2006) Long-term evaluation of deep brain stimulation of the thalamus. J Neurosurg 104: 506–512.
  38. 38. Critchley GR, Richardson PL (1998) Vim thalamotomy for the relief of the intention tremor of multiple sclerosis. Br J Neurosurg 12: 559–562.
  39. 39. Broggi G, Brock S, Franzini A, Geminiani G (1993) A case of posttraumatic tremor treated by chronic stimulation of the thalamus. Mov Disord 8: 206–208.
  40. 40. Broggi G, Franzini A, Tringali G, Ferroli P, Marras C, et al. (2006) Deep brain stimulation as a functional scalpel. Acta Neurochir Suppl 99: 13–19.
  41. 41. Foote KD, Seignourel P, Fernandez HH, Romrell J, Whidden E, et al. (2006) Dual electrode thalamic deep brain stimulation for the treatment of posttraumatic and multiple sclerosis tremor. Neurosurgery 58: ONS–280-285; discussion ONS-285-286.
  42. 42. Muthusamy KA, Aravamuthan BR, Kringelbach ML, Jenkinson N, Voets NL, et al. (2007) Connectivity of the human pedunculopontine nucleus region and diffusion tensor imaging in surgical targeting. J Neurosurg 107: 814–820.