Comparison of MRI-based automated segmentation methods and functional neurosurgery targeting with direct visualization of the Ventro-intermediate thalamic nucleus at 7T

The ventro-intermediate nucleus (Vim), as part of the motor thalamic nuclei, is a commonly used target in functional stereotactic neurosurgery for treatment of drug-resistant tremor. As it cannot be directly visualized on routinely used magnetic resonance imaging (MRI), its clinical targeting is performed using indirect methods. Recent literature suggests that the Vim can be directly visualized on susceptibility-weighted imaging (SWI) acquired at 7 T. Our work aims to assess the distinguishable Vim on 7 T SWI in both healthy-population and patients and, using it as a reference, to compare it with: (1) The clinical targeting, (2) The automated parcellation of thalamic subparts based on 3 T diffusion MRI (dMRI), and (3) The multi-atlas segmentation techniques. In 95.2% of the data, the manual outline was adjacent to the inferior lateral border of the dMRI-based motor-nuclei group, while in 77.8% of the involved cases, its ventral part enclosed the Guiot points. Moreover, the late MRI signature in the patients was always observed in the anterior part of the manual delineation and it overlapped with the multi-atlas outline. Overall, our study provides new insight on Vim discrimination through MRI and imply novel strategies for its automated segmentation, thereby opening new perspectives for standardizing the clinical targeting.

The Vim cannot be directly visualized on current magnetic resonance imaging (MRI) sequences that are routinely used in clinical practice. The targeting methods are therefore indirect, employing either stereotactic coordinates or the quadrilatere of Guiot 8,9,[13][14][15] . The quadrilatere of Guiot has been established based upon the Vim definition by Guiot and Albe-Fesard while performing electrophysiological recordings 16 . It is a robust and reliable method and is used both in the field of the reference technique (e.g. DBS) or in RS.
In the RS treatment settings, there is no intraoperative confirmation of the target. However, clinical results have shown a reduction in tremor amplitude, comparable to the reference technique 14 . In fact, good clinical outcome was achieved in more than 70% of the cases 14,17 , using the indirect targeting of the quadrilatere of Guiot. Better clinical outcome is usually correlated with higher 1-year MR signature volumes 18,19 . Radiological hyporesponders are rather exceptional, usually accounting for approximately 20% of patients and are frequently, but not always, related to lower clinical improvement.
Nevertheless, the used targeting techniques, including the quadrilatere of Giout, are built upon stereotactic standard landmarks of the brain, not necessarily in a well-established relationship with the underlying thalamic morphology 2 and therefore, are not sensitive to inter-subject variability 20 . Additionally, the potential of variability increases with the lack of consensus for a "gold standard" targeting method among the centres. Such a drawback is particularly crucial for Vim RS, which cannot rely on intraoperative target confirmation. Hence, there is a need for an improvement and standardization the targeting procedure.
Advanced MRI techniques applied for acquisitions at ultra-high field could help to address the limitation in Vim-targeting. Recently, Abosch et al. 21 suggested the possibility of a direct Vim visualization by using susceptibility weighted imaging (SWI) acquired at 7 T that showing an enhanced image contrast inside the thalamic area 21,22 (see Fig. 1).
Susceptibility-weighted imaging (SWI) is a MRI technique particularly sensitive to magnetic susceptibility effects 23,24 . For instance, de-oxygenated blood vessels, as well as the calcium and iron-loaded tissue demonstrate, with respect to the surroundings, strong magnetic susceptibility contrast. Therefore, in such regions, SWI provides an enhanced contrast when compared to standard T1, T2, or T2*-weighted imaging. Aside the increased signal-to-noise ratio that can be traded for improved spatial resolution, ultra-high magnetic field strength, such as 7T, has a superior sensitivity to magnetic susceptibility-related contrasts.
SWI acquired at 7T has already been used for the study of other deep brain structures such as the subthalamic nucleus (STN) and substantia nigra (SN) 21,22,25 . A related technique often used in the context of targeting is the quantitative susceptibility mapping (QSM), where the phase of the MR signal is used to compute a map of the susceptibility distribution. QSM has been mainly explored for a direct visualization of the STN 26-28 , but its contrast variation inside the thalamic area is however limited for discriminating the Vim 27,28 . The existing approach for improving the QSM potential for the visualization of the Vim 29 requires data to be acquired with subjects' head in various positions in respect to the magnetic field, which is very demanding in terms of clinical feasibility. Moreover, until now, other than a descriptive study performed by Abosch et al., no investigation based on 7T MRI, involved a profound analysis of the thalamic subparts or the Vim in particular.
Other advanced high field MRI techniques have been explored for thalamic nuclei identification. Particularly, diffusion MRI (dMRI) has gained interest in this domain thanks to its ability to depict exquisite microstructural details related to the orientation of the white matter fibres, their coherence and the diffusivity affected by the grey Figure 1. Illustration of the visible structures corresponding to the thalamic area in axial view on the SWI acquired at 7T, based upon the illustration from the pioneering work of Abosch et al. 21 . In panel A the SWI features are compared with the Schalterbrand atlas (plate 53 Hd + 3.5) superimposed on the right thalamus. The arrows and the respective color contours indicate the nuclei: Vim, Pulvinar (Pu), ventro-caudalus (Vc), ventroodalis (Vo), the medio-dorsal group (MD) and the internal capsule (ic). The shown SWI image is part of the dataset used in this study. Panel B gives a corresponding axial plate of the Morel's atlas where the (same) color (shade) matches appropriate regions of the Shalternbrand's atlas, while keeping the same nomenclature used in each one of them. In fact, considering the Morel's atlas 59 nomenclature, Vim is part of the Ventro-Lateralposterior nuclei, which furthermore, together with ventro-lateral anterior and ventro-posterior nuclei form the ventral latero-ventral (VLV) group of nuclei. matter (e.g. size of neurons or their arrangement), which are different for each thalamic nucleus. However, with the disadvantage of relatively low spatial resolution, the present standard dMRI techniques are limited in distinguishing small thalamic regions such as the Vim.
In this study we aim to take a step forward in extended and quantitative analysis of the advantages brought by 7T SWI. Consequently, our goal was a better direct characterization of the Vim. To achieve this, we studied two different types of subjects: a healthy population, in the frame of a preliminary analysis, and a drug-resistant tremor patients' population, to be able to evaluate the applicability of the previous on real clinical cases. Additionally, in order to assess the performance of the known approaches, the visually distinguishable Vim on 7T SWI will be also compared with three different state-of-the-art methods, namely: (1) the quadrilatere of Guiot, used routinely as clinical targeting and built on 3T data; (2) The automated 3T dMRI-based thalamic subdivision; and (3) The atlas-based Vim outline employing both 3T and 7T data. Each one of these methods is defining the Vim in a different mode. The first method represents the Vim in a functional manner based on statistical findings from previous electrophysiological recordings done on a group of tremor patients. Its outcome is a single point only, that is expected to fall inside the Vim 16 . The remaining two methods, as complementary to the first one, are part of computer-assisted image-analysis techniques and reveal the anatomical structure of the Vim, but differ between them in terms of size of the provided outline. Figure 2 shows an illustrative overview of the performed comparisons.

Results
Nine healthy subjects without particular brain disease and/or deformation caused by intracranial lesions were scanned at both 3T and 7T. Five of them were relatively young (25+/−2 years old, 2 females, will be referred further as YS1-YS5 respectively) and the remaining four represented the older population (67.2+/−9.5 years old, 1 female, accordingly referred further as ES1-ES4). Please refer to Methods for the acquisition details.
The four approaches investigated for defining the Vim provide outcomes at different sizes: from a single point (e.g. the quadrilatere of Guiot) to a region that encloses a group of nuclei (segmentation from diffusion data). To better assess the spatial distribution, we further divided each region of interest (ROI), either manually or automatically obtained, into eight geometrical sub-regions (refer to Fig. S1 in the supplementary material). More specifically, for the ROI, we calculated the smallest rectangular cuboid containing all non-zero voxels whose mid-plane isolated the superior from the inferior part, while the in-plane diagonals isolated the anterior, lateral, posterior, and ventral parts of the ROI.
The quantitative evaluation of the results included: • Estimating the intra-subject reproducibility of the Guiot targeting procedures per hemisphere for each subject (each repeated 6 times), by employing the Euclidean distance between the respective points, taking as reference the first targeted point; Manual Vim delineation on 7T SWI. The Vims were identified and manually delineated in both hemispheres for eight out of nine subjects. In the ninth case, only the right Vim area was distinguishable, while on the left side, a blood vessel made it difficult to discriminate any contrast difference that could have corresponded to the Vim (Figs 3 and 4). Hence, in total, we had 17 manually delineated Vims. The manual delineation was nonetheless challenging and time-consuming. In fact, the area of the Vim, appearing as a hyper-intensity structure, was directly detectable, but its peripheral borders were not always unquestionably discernible, due either to the lack of a consistent image intensity contrast or the presence of small blood vessels that sometimes could have been confounded with the surrounding nuclei. Moreover, the image intensity and its contrast variations differed between subjects and different age ranges (Fig. 4).
With the aim of providing reliable outlines, the rater chose not to include in the delineation the surrounding, more ambiguous darker regions. These outlines were then considered as a reference for the subsequently performed analyses.
The manually delineated volumes were in the range [35.5, 94.5] mm 3 and occupied 0.4-1.5% of the total thalamic volume (Table 1).
Comparison between the functional and the structural Vim definition. The obtained points from the individually repeated targeting with the quadrilatere of Guiot, a task performed only on the young cohort, were either overlapping or differing by one voxel for all subjects involved in this study (YS1-YS5). The maximum  . Illustration of the difficulties encountered for the manual delineation of the Vim regarding the image contrast on SWI acquired at 7T and the presence of blood vessels in the targeted area. We can observe that the contrast varies between subjects, but also between the two groups, young (here YS1 and YS2) versus the elderly (ES2 and ES3). The arrow for Subject YS1, illustrates the relatively big blood vessel passing through the left targeted thalamic region that prevented manual discrimination of the Vim. The presence of small vessels surrounding the Vim can be observed in each panel.
Scientific RepoRts | (2019) 9:1119 | https://doi.org/10.1038/s41598-018-37825-8 intra-subject difference was less than 1.3 mm, confirming the expected high reproducibility of this targeting (see the boxplots in Fig. 5). The Guiot points were always inside or on the border of the manual delineation. Considering the geometrical subdivision (Fig. S1), for seven out of nine cases the Guiot points were in the ventral part of the manual delineation.
Comparison with the automated methods. Automated 3T dMRI-based segmentation. For all nine subjects, the thalamic parcellation showed a similar segmentation pattern as previously reported for a larger cohort of 35 healthy subjects 31 (Fig. 3). The segmented ventral latero-ventral (VLV) cluster, which is the outline of interest since it encloses the motor-related nuclei including the Vim 31 ( Fig. 1), represented the expected spatial extent and its volume was in the interval [0.74, 1.4] cm 3 or 13-19% of the corresponding thalamic volume. Detailed values are reported in Table 1.
The VLV cluster consistently included the manual Vim outline. The size of the VLV was on average 15 times bigger than the manual outline (see Fig. 3). With respect to the spatial localization assessed via the geometrical subdivision ( Fig. S1), for 16 out of 17 cases the manual Vim outline was found in the inferior VLV, in the last case being around the mid-plane. Furthermore, the manual outline was always next to the lateral VLV border and, in 15 cases, it was situated in the anterior/anterior-lateral part of this cluster. For the remaining two cases, the Vim outline was around the respective in-plane diagonal.
Multi-atlas segmentation. To assure an equal contribution from both hemispheres, the subject with no manual outline of the Vim (YS1) was not considered for building the multi-atlas framework, leading to eight subjects in total for this analysis. Hence, as we applied the leave-one-out (LOO) cross-validation, seven atlases were used per hemisphere.
The multi-atlas segmentation outlines were fairly comparable to the manual delineation ( Fig. 6), with the exception of one subject (ES3) right hemisphere where this methodology failed to provide an outline overlapping with the manual delineation. The calculated Dice overlap was 45.1+/−12.9% for the left and 45.0+/−26.4% for the right hemisphere. The individual Dice coefficients are given in Fig. 6 together with the visual representation of the outcome.
Patients treated with Vim radiosurgery for essential tremor. We applied the same four approaches for defining the Vim, as performed in the healthy population, on two patients, with drug-resistant tremor, treated  with left Vim RS for right-sided drug-resistant essential tremor. One was 86-years old and one 54-years old, referred to as P1 and P2, respectively. Indirect targeting was performed for both using the quadrilatere of Guiot.
A dose of 130 Gy at the 100% isodose line was prescribed in each patient.
In the particular case of these patients, we observed that the Vim is also distinguishable on 7 Tesla imaging, fairly comparable as for the explored healthy population, previously described. The performed manual delineation presented volumes between 36.2 and 43.6 mm 3 . Furthermore, the thalamus parcellation provided again the expected segmentation pattern and the volumes of the VLV clusters were in the interval [0.77,1] cm 3 . Moreover, all four VLVs (left and right side of each patient) enclosed the manual delineation on the 7T SWI in their inferior anterio-lateral segment. Additionally, for three out of four Vims, the multi-atlas based segmentation gave outlines of volume between 19.5 and 40.8 mm 3 and overlap with the corresponding manual delineation in the range between 31.4 and 51.6%. For the right side of P1, the outline enclosed only a few voxels and did not overlap the manual delineation.
Moreover, the MR signature usually seen on the follow-up images after Vim RS (here at 3 and 9 months for P1 and P2 respectively), due to its delayed effect 32,33 , showed the classical "cockade aspect", as in previously published reports 32,34 . The "cockade aspect", as referred in the literature, contains both the peripheral ring of the contrast enhancement as well as the hypo-intensity, central area on T1 weighted images, all corresponding to the 90 Gy isodose line 34 (refer to Fig. 7).
The MR signature overlapped with the anterior half of the manually delineated Vim. Furthermore, the centroid of this signature, which should represent the targeted point during the treatment, for P1, was in the ventral, while for P2 in the anterior Vim's segment delineated geometrically (conforming to the Fig. S1 from the supplementary material). Additionally, these centroids were, in both cases, inside the multi-atlas segmentation outline. Overall results of the patients' data are illustrated in Fig. 7.

Discussion
This study presents extensive Vim-related analyses of 7T SWI whose enhanced image contrast in the central brain allows a direct visualization of the thalamic subparts. This advantage of directly distinguishing the Vim was demonstrated here for data coming from both healthy volunteers representing broad age intervals and drug-resistant essential tremor patients. Additionally, having the manual delineation of the Vim done on 7T SWI, the study also compares the current indirect clinical targeting and the existing automated image-segmentation techniques against this outline. Hence, the presented work provides a reinforced knowledge characterizing the Vim.
In the young healthy-population data, we observed that the Guiot's targeting points, are for 77.8% of the cases in the geometrically delineated ventral part of the manual Vim delineation based on 7T SWI. These observations were in line with the findings from the patients' data, where the MRI signature appearing on the follow-up images was enclosed in the anterior half of the manual outline. Although the Guiot targeting was not performed on the healthy elderly population, observing the high robustness of the outcome in young subjects and patients, we assume that the same findings would be found in any additional data.
The Guiot targeting allocates only a single point and, as a consequence, is limited in providing a full structural discrimination of the Vim. In this context, our findings, obtained from 11 Vims in total, indicate a confined localization of the clinically aimed points inside a structural outline. Hence, the ventral part of the Vim should be further investigated in terms of connectivity and function for potentially obtaining advanced target-related insights. Moreover, an extension of this study on a larger dataset could possibly lead to a robust prediction of the neurosurgical target emplacement inside the Vim's structural borders. Our results also suggested that the manual delineation, 15 times smaller than the size of the VLV, is mainly enclosed (95.2% of the studied cases) in the inferior anterior-lateral division of the VLV cluster, which is anatomically relevant. Moreover, such findings reveal a constrained area that most likely corresponds to the Vim within a region automatically delineated by computer-assisted image-analysis techniques, which has not been demonstrated previously. Future developments could exploit this new spatial prior knowledge as a constraint for further subdivision of the VLV cluster, for instance based again on dMRI or in combination, if available, with 7T SWI features.
The proposed multi-atlas segmentation framework shows promising potential of delineating the Vim area in automated manner. Although the outcome tends to slightly underestimate the volume of the Vim, in both patients' data it overlapped with the late MRI signature enclosing the centroid of this region. In the healthy population data, as well, the multi-atlas outline is relatively small but, however, fairly acceptable, mainly in terms of identifying the Vim area. With regard to the quantitative estimation of the overlap, a Dice Coefficient greater than 70%, is commonly considered as a good match. Nevertheless, this value has been established for evaluating the segmentation accuracy for relatively big structures. In fact, Dice is more sensitive to the variation of edge voxels than the inner ones and therefore, for small and elongated structures as the Vim, one could argue that such measure could bias the estimation of the overlap 35 .
We assume that the presented promising potential of the proposed multi-atlas segmentation is closely related to the use of SWI-intensity instead of T1w or T2w, lacking in intrathalamic contrast variation. Consequently, when 7T SWI is available, this approach could be used as an alternative or as an initialization of a subject-related segmentation of the Vim, either manual or automated one. Nevertheless, we expect improvement of the presented findings with an extension/adaptation of the number of employed atlases to its optimum 36 .
Although our main aim focuses on the Vim, this study also provides validation of the automated dMRI-based clustering against the directly visible SWI regions inside the thalamus. For example, in SWI the Pulvinar appears as the most-posterior distinguishable thalamic feature very dark laterally and brighter next to the ventricles (Fig. 3). Accordingly to previously reported findings 31 , in Fig. 3 we see that the green Pu-cluster contour delineates mainly the darker nuclei part and a portion of medial Pulvinar (the brighter part) is enclosed in the cluster representing the centro-lateral and the lateral posterior nuclei along with the medial part of the Pulvinar (CL-LP-PuM). Furthermore, the spatial extent of the medio-dorsal (MD) cluster coincides in general to the equivalent group of nuclei. All these findings support, above all, the spatial distribution of the VLV and thereupon, our findings related to this cluster. The manual delineation on the 7T SWI data is however a time-consuming task. In fact, the region corresponding to the Vim is, as shown by our data, visually directly distinguishable both in healthy controls and in patients. However, the difficulties arise in the presence of blood vessels or a low local image contrast, both of them making more challenging the precise discrimination of the peripheral borders of this nucleus. Hence, the provided outline encloses only the region uncontestably corresponding to the Vim. Or, in other words, it represents well the Vim, but tends towards a slight underestimation of its extent. This could also be deduced from the calculated volumes (Table 1), which are near the lower limits of the known range 3 .
Moreover, despite the fact that no differences were noticed between the SWI acquired from the patients and the elderly healthy volunteers, we however observed inconsistency in terms of the contrast variation between the cohorts representing different age ranges. Namely, the elderly subjects presented poorer contrast, which increased the difficulty in immediate identification of the Vim in those cases (Fig. 4). Nevertheless, the intrathalamic contrast remains globally similar and adequate for outlining the supposed Vim area.
Supplementary SWI analysis indicated that the image intensity contrast could also depend on the head orientation inside the scanner. This could be caused by anisotropic contributions to the local tissue susceptibility. For instance, white matter fibres have significantly anisotropic susceptibility 37 . Depending on the head orientation relative to the static MR field, the local susceptibility-based perturbations of the field may differ, and the resulting effects may thereby show significant changes in contrast. Anisotropic voxels may also add different partial volume effects depending on the head orientation.
The multi-atlas segmentation could also be troubled by the mentioned contrast-related drawbacks. We think that the relatively small volumes obtained from the multi-atlas segmentation are related on one side to eventual minor registration errors and on the other side to the Joint Label Fusion method we apply a posteriori, as both registration and fusion can be disturbed by the random presence of vessels. Moreover, the stated drawbacks could also be associated with the non-symmetrical tendency of the multi-atlas outcome. We assume as well that using a bigger cohort will not only allow finding the optimum number of atlases, but it could also decrease the left-right asymmetry between the multi-atlas outlines.
Regarding the stated limitations, further improvements including contrast standardization and vessel removal 38 should be addressed in order to make 7T SWI an even more powerful and robust tool for Vim discrimination. Additionally, such standardization could be beneficial for further investigating if the 7T SWI is capable to capture eventual anatomical differences between the controls and the patients' data.
The SWI sequences are also prone to spatial distortions, but as Duchin et al. have already demonstrated 39 , the central brain area and therefore the thalamus as well, are negligibly affected by those distortions. However, to compensate for the peripheral spatial distortions observed on the 7T data, an affine transform was employed to match them to the corresponding 3T images. Additionally, considering the motion, even though the tremor patients are more prone to such artefacts, there were no particular differences between the healthy subjects and the patients in terms of the acquired data and the obtained findings.
Although the number of included clinical cases is low, the findings observed in the patients are in line with those from the explored healthy population. Moreover, as both of the patients are presenting a good clinical outcome, their contribution to this work becomes even more valuable.
Nevertheless, in line with the pioneer work of Abosch et al. 21 , which is mainly in the scope of defining the subthalamic nucleus and other basal ganglia structures on 7T SWI, by exploring a larger dataset from healthy subjects and, moreover, extending it to clinical cases, the present study demonstrates even further the high potential of SWI acquired at 7T to provide more precise characterization of the Vim as a structure than the other existing methods currently used to define this nucleus. Furthermore, our comparison analyses, where the manual Vim delineation done on 7T SWI was taken as a reference, indicate two new promising approaches for an automated localization and delineation of the Vim: 1) the localization likelihood inside the cluster obtained from the robust and reproducible 3T dMRI-based thalamic parcellation and 2) the multi-atlas segmentation exploring the 7T SWI information.
Further studies on a larger dataset and including more raters for the manual delineation should confirm the reported findings. Improvements could also include multi-atlas segmentation based on 7T SWI only. Moreover, given the potential of the automated dMRI-based segmentation, we assume that combining it with the complementary 7T SWI information would further lead to a more accurate Vim outline based only on the individual anatomy.

Methods
Local institutional review board named Commission cantonale d'éthique de la recherche sur l'être humain (CER-VD) approved the study and all participants gave written informed consent. All the analyses performed in the presented work were done in accordance with the relevant guidelines and regulation stated in the obtained approval from CER-VD with protocol number 165/15.  Table 2.
As an added value, to evaluate the appearance of the MR signature after Vim RS, Gadolinium-enhanced MPRAGE was acquired every three months for the two patients, with drug-resistant essential tremor. The acquisition was done at 3T with similar parameters, as for the preoperational 3T MPRAGE. The first patient (P1) showed decreased tremor amplitude on the right, contralateral hand at four months, with further complete tremor arrest at six months, while the second patient (P2), had major tremor decrease at 3 months after the procedure. This type of clinical answer is common after Vim RS for tremor 18 . The subjects were carefully immobilized with head positioning pads to minimize motion and were instructed to avoid moving during the individual image acquisitions.
Common image space. The analyses, for each subject respectively, were performed in the individual anterior commissure -posterior commissure (AC-PC) image space.
For the young population, we first transformed the T2w into the AC-PC space by employing 3D Slicer 41 and choosing manually 30 brain mid-plane points. The resulting image was then used as a reference for the AC-PC alignment of the remaining sequences. More precisely, for each 3T contrast, with the exception of the DWI, we performed a rigid-body transformation (six degrees of freedom, DOF).
For the dMRI data, we registered the fractional anisotropy (FA) map to the MPRAGE with a non-linear transform using FSL FNIRT 42,43 , to simultaneously account for the echo planar imaging (EPI) distortions 44 .
Moreover, the observed distortions in the frontal and the parietal cortex of the 7T data were compensated with a 12-DOF linear transform between both skull-stripped 7T MP2RAGE and 3T MPRAGE. This transform was applied to the SWI rigidly (6-DOF) aligned beforehand to the corresponding MP2RAGE.
All registrations 45 were performed with 100,000 voxel samples and Mattes Mutual Information as the cost function. The outcome quality was assessed by visual inspection of the matching between the ventricles. Figure S2 in the supplementary material summarizes the registration process. After being transformed to the AC-PC space, all images were resampled to the T2 CISS spatial resolution (0.4 × 0.4 × 0.4 mm 3 ) using linear interpolation.
Since the 3T images of the elderly cohort were acquired directly in AC-PC space, the AC-PC transformation was not required in these cases, but correspondence between the respective images was however performed according to the previously described flow for the young cohort. After being transformed, all the images were interpolated to the voxel size of 0.5 × 0.5 × 0.5 mm 3 defined by the T1 template in MNI space 46 . Manual Vim delineation on 7T SWI. SWI combines contrast from differences in effective transverse signal relaxation (T2*) and in MR signal phase, both strongly affected by the magnetic susceptibility and the underlying tissue geometry 23,24 . Accordingly, SWI has shown increased sensitivity with respect to standard T1-or T2*-weighted imaging for visualizing deoxygenated blood vessels, vascular and axonal lesions, as well as calcium and iron depositions 47 and it allows better visualization of certain deep brain structures, including the thalamic subparts. Particularly, the Vim appears as a well-distinguishable hyper-intensity structure surrounded by darker regions 21 (Fig. 1a). This is especially evident at higher magnetic field strengths, such as 7T, where the susceptibility-dependent contrast is more accentuated, and moreover, sub-millimeter spatial resolution can be easily achieved.
The manual delineation of the Vim on 7T SWI was done by an experienced neurosurgeon (CT) and was primarily based upon the previous observations of Abosch et al. 21 and the Schaltenbrand and Wahren atlas 48 also used as reference in the mentioned study.
The thalamic lateral border and the internal capsule are well distinguishable on the 7T SWI. Referring to the Shaltenbrand's atlas, the Pulvinar as well is directly recognizable in the axial SWI plane (in green, Fig. 1a). Superior to it, again in axial plane, is the ventro-caudalus nucleus (Vc) appearing as a dark region (Fig. 1a, bright-violet) immediately posterior to a narrower zone with brighter image intensity that is considered as Vim (Fig. 1a, yellow). Above it, ventro-odalis (Vo) appears (Fig. 1a, bright-  As the targeting procedure was done outside of regular treatment-planning software, the quadrilatere of Guiot was built independently six times for each subject, bilaterally and blindly with the aim of assessing its reproducibility. A sphere with 2 mm radius was drawn around each Guiot point to simulate both the 90 Gy isodose line and the contrast-enhanced area 49 visualized on follow-up MRI if a radiosurgical thalamotomy 50 would have been applied in the presented cases. Automated 3T dMRI-based segmentation. Several research groups explored the dMRI advantages for automated thalamic subdivision 31,[51][52][53][54][55][56][57][58] . The initial approaches are mainly based on diffusion tensor images i.e. a coarse diffusion representation. Recently, our group developed an approach 31 that relies on a more detailed diffusion feature -the Spherical Harmonic (SH) representation of the Orientation Distribution Functions (ODFs) and outperforms the state-of-the-art methods by providing robust and reproducible segmentation pattern for a large dataset and different diffusion sequences. Our approach subdivides the thalamus into seven groups of nuclei corresponding to the anatomy in Morel's atlas 59 . One of the segmented groups of nuclei is the ventral latero-ventral cluster (VLV) enclosing the motor-related nuclei including the Vim (Fig. 1).
The pre-processing of the dMRI data included data denoising [60][61][62] , bias field 43,63 , motion 64 and eddy current correction 65 . The SH coefficients of the constant solid angle ODFs of maximum order 6 were calculated using the FSL's qboot function 31 .
The thalamic masks were initially obtained from the Freesurfer parcellation 66,67 on each MPRAGE and redefined, as described by Battistella et al. 31 , by eliminating the voxels having FA value greater than 0.55 or exceeding 5% probability to belong to the cerebrospinal fluid. The refined masks were further compared to the thalamic borders appearing more evident in SWI and consequently, for providing even more precise thalamic outline, several voxels were manually refined.
The thalamic subdivision was completed in the diffusion space of each subject and the resulting clusters were brought in the AC-PC space.
Multi-atlas segmentation. One widely used state-of-the-art segmentation technique is the atlas-based registration, if an atlas is available. Having achieved the manual Vim outline, we considered each one of them as an atlas within a multi-atlas segmentation framework.
The framework was built in leave-one-out scenario, meaning that for each subject-target its own Vim outlines were not considered as atlases, but only those from the remaining subjects. The atlas registration was done in 2 steps.
First, with an affine transform 45 we matched the respective MPRAGEs and we applied it to the corresponding SWI. Second, we performed a non-linear registration for a local correspondence using only the SWI as it gives more distinguishable thalamic features than the T1w, thus potentially allowing more accurate matching between the thalamic nuclei.
To this end, we defined a volume of interest (VOI) surrounding the thalamus in which, in order to standardize the SWI intensities among all the subjects, we performed a histogram matching proposed by Nyiul et al. 68 . The standard scale of this intensity matching method was built upon the histogram deciles from the SWI-VOI of the elderly population and was then used as the reference scale for the one-to-one intensity mapping of all the SWI-VOIs.
Non-rigid BSpline transform using maximum displacement of 1 mm 45 was computed in between the respective SWI-VOIs. The combination of both linear and non-linear transform was applied to the Vim outlines respectively and the final multi-atlas segmentation outcome was obtained by employing the Joint Label Fusion method with corrective learning 69 with the default parameters.

Data Availability
All the original images used in the presented study and the corresponding manual delineation of both left and right Vim are uploaded in Zenodo (https://doi.org/10.5281/zenodo.1438358).