Objective To evaluate the effects of bilateral caudal zona incerta (cZi) deep brain stimulation (DBS) for Parkinson's disease (PD) one year after surgery and to create anatomical improvement maps based on patient-specific simulation of the electric field. Materials and Methods We report the one-year results of bilateral cZi-DBS in 15 patients with PD. Patients were evaluated on/off medication and stimulation using the Unified Parkinson's Disease Rating Scale (UPDRS). Main outcomes were changes in motor symptoms (UPDRS-III) and quality of life according to Parkinson's Disease Questionnaire-39 (PDQ-39). Secondary outcomes included efficacy profile according to sub-items of UPDRS-III and simulation of the electric field distribution around the DBS lead using the finite element method. Simulations from all patients were transformed to one common magnetic resonance imaging template space for the creation of improvement maps and anatomical evaluation. Results Median UPDRS-III score off medication improved from 40 at baseline to 21 on stimulation at one-year follow-up (48%, p < 0.0005). PDQ-39 summary index did not change, but the subdomain activities of daily living (ADL) and stigma improved (25%, p < 0.03 and 75%, p < 0.01), whereas communication worsened (p < 0.03). For UPDRS-III sub-items, stimulation alone reduced median tremor score by 9 points, akinesia by 3, and rigidity by 2 points at one-year follow-up in comparison to baseline (90%, 25%, and 29%, respectively, p < 0.01). Visual analysis of the anatomical improvement maps based on simulated electrical fields showed no evident relation with the degree of symptom improvement and neither did statistical analysis show any significant correlation. Conclusions Bilateral cZi-DBS alleviates motor symptoms, especially tremor, and improves ADL and stigma in PD patients one year after surgery. Improvement maps may be a useful tool for visualizing the spread of the electric field. However, there was no clear-cut relation between anatomical location of the electric field and the degree of symptom relief.

Deep brain stimulation (DBS) in zona incerta (Zi) is used for symptom alleviation in essential tremor (ET). Zi is positioned along the dentato-rubro-thalamic tract (DRT). Electric field simulations with the finite element method (FEM) can be used for estimation of a volume where the stimulation affects the tissue by applying a fixed isolevel (VDBS). This work aims to develop a workflow for combined patient-specific electric field simulation and white matter tracing of the DRT, and to investigate the influence on the VDBS from different brain tissue models, lead design and stimulation modes. The novelty of this work lies in the combination of all these components. Patients with ET were implanted in Zi (lead 3389, n = 3, voltage mode; directional lead 6172, n = 1, current mode). Probabilistic reconstruction from diffusion MRI (dMRI) of the DRT (n = 8) was computed with FSL Toolbox. Brain tissue models were created for each patient (two homogenous, one heterogenous isotropic, one heterogenous anisotropic) and the respective VDBS (n = 48) calculated from the Comsol Multiphysics FEM simulations. The DRT and VDBS were visualized with 3DSlicer and superimposed on the preoperative T2 MRI, and the common volumes calculated. Dice Coefficient (DC) and level of anisotropy were used to evaluate and compare the brain models. Combined patient-specific tractography and electric field simulation was designed and evaluated, and all patients showed benefit from DBS. All VDBS overlapped the reconstructed DRT. Current stimulation showed prominent difference between the tissue models, where the homogenous grey matter deviated most (67 < DC < 69). Result from heterogenous isotropic and anisotropic models were similar (DC > 0.95), however the anisotropic model consistently generated larger volumes related to a greater extension of the electric field along the DBS lead. Independent of tissue model, the steering effect of the directional lead was evident and consistent. A workflow for patient-specific electric field simulations in combination with reconstruction of DRT was successfully implemented. Accurate tissue classification is essential for electric field simulations, especially when using the current control stimulation. With an accurate targeting and tractography reconstruction, directional leads have the potential to tailor the electric field into the desired region.

Group analysis of patients with deep brain stimulation (DBS) has the potential to help understand and optimize the treatment of patients with movement disorders. Probabilistic stimulation maps (PSM) are commonly used to analyze the correlation between tissue stimulation and symptomatic effect but are applied with different methodological variations. To compute a group-specific MRI template and PSMs for investigating the impact of PSM model parameters. Improvement and occurrence of dizziness in 68 essential tremor patients implanted in caudal zona incerta were analyzed. The input data includes the best parameters for each electrode contact (screening), and the clinically used settings. Patient-specific electric field simulations (n = 488) were computed for all DBS settings. The electric fields were transformed to a group-specific MRI template for analysis and visualization. The different comparisons were based on PSMs representing occurrence (N-map), mean improvement (M-map), weighted mean improvement (wM-map), and voxel-wise t-statistics (p-map). These maps were used to investigate the impact from input data (clinical/screening settings), clustering methods, sampling resolution, and weighting function. Screening or clinical settings showed the largest impacts on the PSMs. The average differences of wM-maps were 12.4 and 18.2% points for the left and right sides respectively. Extracting clusters based on wM-map or p-map showed notable variation in volumes, while positioning was similar. The impact on the PSMs was small from weighting functions, except for a clear shift in the positioning of the wM-map clusters. The distribution of the input data and the clustering method are most important to consider when creating PSMs for studying the relationship between anatomy and DBS outcome.

Deep brain stimulation (DBS) is a therapy for Parkinson’s disease (PD) and essential tremor (ET). The mechanism of action of DBS is still incompletely understood. Retrospective group analysis of intra-operative data recorded from ET patients implanted in the ventral intermediate nucleus of the thalamus (Vim) is rare. Intra-operative stimulation tests generate rich data and their use in group analysis has not yet been explored. Objective. To implement, evaluate, and apply a group analysis workflow to generate probabilistic stimulation maps (PSMs) using intra-operative stimulation data from ET patients implanted in Vim. Approach. A group-specific anatomical template was constructed based on the magnetic resonance imaging scans of 6 ET patients and 13 PD patients. Intra-operative test data (total: n = 1821) from the 6 ET patients was analyzed: patient-specific electric field simulations together with tremor assessments obtained by a wrist-based acceleration sensor were transferred to this template. Occurrence and weighted mean maps were generated. Voxels associated with symptomatic response were identified through a linear mixed model approach to form a PSM. Improvements predicted by the PSM were compared to those clinically assessed. Finally, the PSM clusters were compared to those obtained in a multicenter study using data from chronic stimulation effects in ET. Main results. Regions responsible for improvement identified on the PSM were in the posterior sub-thalamic area (PSA) and at the border between the Vim and ventro-oral nucleus of the thalamus (VO). The comparison with literature revealed a center-to-center distance of less than 5 mm and an overlap score (Dice) of 0.4 between the significant clusters. Our workflow and intra-operative test data from 6 ET-Vim patients identified effective stimulation areas in PSA and around Vim and VO, affirming existing medical literature. Significance. This study supports the potential of probabilistic analysis of intra-operative stimulation test data to reveal DBS’s action mechanisms and to assist surgical planning.

Deep Brain Stimulation (DBS) is an established therapy for movement and neuropsychiatric disorders. Identifying brain regions (Probabilistic Sweet Spots, PSS) linked with the greatest symptom improvement is crucial for refining pre-operative targeting and post-operative programming. Probabilistic stimulation mapping is a powerful data-driven tool to delineate these regions. However, the chosen statistical methods can influence the identified PSS. A comprehensive evaluation of their impact is lacking in DBS research. The present study compares the PSS generated with four voxel-wise statistical approaches - t-test, Wilcoxon test, Linear Mixed Model, and Bayesian t-test - with the aim of assessing their influence on computed results on the same dataset. Intra-operative stimulation test data of 23 Essential Tremor (ET) patients was used to run patient-specific electric field simulations and to generate PSS in a group-specific anatomical template space. The PSS for the different statistical tests were first compared in terms of size and topography. Then, their correlation with clinical improvement was calculated in a leave-one-out cross-validation scheme and PSS consistency across datasets with different compositions was assessed. Our findings emphasize the impact of statistical test selection on both the anatomical location and volume of the extracted PSS, highlighting the importance of careful methodological choices in future DBS mapping studies. The Bayesian t-test and a voxel-wise application of nonparametric permutation testing, introduced for the first time in DBS research, showed promising results in identifying PSS representative of improvement and exhibited robustness to variations in the dataset.

Probabilistic stimulation maps of deep brain stimulation (DBS) effect based on voxel-wise statistics (p-maps) have increased in literature over the last decade. These p-maps require correction for Type-1 errors due to multiple testing based on the same data. Some analyses do not reach overall significance, and this study aims to evaluate the impact of sample size on p-map computation. A dataset of 61 essential tremor patients treated with DBS was used for the investigation. Each patient contributed with four stimulation settings, one for each contact. From the dataset, 5 to 61 patients were randomly sampled with replacement for computation of p-maps and extraction of high- and low-improvement volumes. For each sample size, the process was iterated 20 times with new samples generating in total 1140 maps. The overall p-value corrected for multiple comparisons, significance volumes, and dice coefficients (DC) of the volumes within each sample size were evaluated. With less than 30 patients (120 simulations) in the sample, the variation in overall significance was larger and the median significance volumes increased with sample size. Above 120 simulations, the trends stabilize but present some variations in cluster location, with a highest median DC of 0.73 for n = 57. The variation in location was mainly related to the region between the high- and low-improvement clusters. In conclusion, p-maps created with small sample sizes should be evaluated with caution, and above 120 simulations in single-center studies are probably required for stable results.