Spherical statistics for characterizing the spatial distribution of deep brain stimulation effects on neuronal activity

Abstract

Background Computational models of deep brain stimulation (DBS) have played a key role in understanding its physiological mechanisms. By estimating a volume of tissue directly modulated by DBS, one can relate the neuronal pathways within those volumes to the therapeutic efficacy of a particular DBS setting. New Method A spherical statistical framework is described to quantify and determine salient features of such morphologies using visualization techniques, empirical shape analysis, and formal hypothesis testing. This framework is shown using a 3D model of thalamocortical neurons surrounding a radially-segmented DBS array. Results We show that neuronal population volumes modulated by various DBS electrode configurations can be characterized by parametric distribution models, such as Kent and Watson girdle models. Distribution parameters were found to change with stimulus settings, including amplitude and radial distance from the DBS array. Increasing stimulation amplitude through a single electrode resulted in more diffuse neuronal activation and increased rotational symmetry about the mean direction of the activated population. When stimulation amplitude was held constant, the activated neuronal population distribution was more concentrated with distance from the DBS array and was also more rotationally asymmetric. We also show how data representation (e.g. stimulus-entrained cell body vs. axon node) can significantly alter model distribution shape. Comparison to Existing Methods This statistical framework provides a quantitative method to analyze the spatial morphologies of DBS-induced effects on neuronal activity. Conclusions The application of spherical statistics to assess spatial distributions of neuronal activity has potential usefulness for numerous other recording, labeling, and stimulation modalities.

On Optimal Tests for Rotational Symmetry Against New Classes of Hyperspherical Distributions

“…They used mixtures of vMF to model their concept of selectivity profiles, each component of the mixture corresponding to a functional system concentrated about a favored direction on the sphere. Both rotationally and non-rotationally symmetric models are used in Xiao and Johnson (2015) in the study of neuronal population volumes. The objective was to fit the distribution of neuronal spike activity around one or more active electrodes.…”

Multimodal 7T Imaging of Thalamic Nuclei for Preclinical Deep Brain Stimulation Applications

“…In this way, the combined use of imaging techniques described in this study can assist in neurosurgical navigation of DBS targets in a given subject (Kamiryo and Laws, 1996 ; Dormont et al, 1997 ; Dipierro et al, 1999 ). Additionally, the segmented nuclei reconstructions can also aid in the development of more accurate computational models of DBS (Mcintyre et al, 2004 ; Kuncel et al, 2008 ; Keane et al, 2012 ; Zitella et al, 2013 ) to retrospectively quantify the neural pathways modulated by thalamic DBS therapy (Xiao and Johnson, 2015 ) or prospectively predict the stimulation settings necessary to target those pathways on a subject-specific basis (Xiao et al, 2016 ).…”