Population coding of forelimb joint kinematics by peripheral afferents in monkeys

Distribution of Large and Small Dorsal Root Ganglion Neurons in Common Marmosets

The aim of this study was to elucidate the size and distribution of dorsal root ganglion (DRG) neurons in non-human primates and to compare them with those of rodent DRG neurons. By measuring the size of NeuN-, NF200-, and peripherin-positive DRG neurons in the lumbar spinal cord of rats and marmosets, we found that the cell size distribution pattern was comparable in both species, although DRG neurons in marmosets were larger than those of rodents. This is the first demonstration that DRG neurons in marmosets have a bimodal size distribution, which has been well established in rodents and humans.

Information about contact force and surface texture is mixed in the firing rates of cutaneous afferent neurons

Cutaneous mechanoreceptors in our hands gather information about the objects we handle. Tactile fibers encode mixed information about contact events and object properties. Neural coding in tactile afferents is typically studied by varying a single aspect of tactile stimuli, avoiding the confounds of real-world haptic interactions. We instead record responses of small populations of dorsal root ganglia (DRG) neurons to variable tactile stimuli and find that neurons primarily respond to force, though some texture information can be detected. Tactile nerve fibers convey information about many features of haptic interactions, including the force and speed of contact, as well as the texture and shape of the objects being handled. How we perceive these object features is relatively unaffected by the forces and movements we use when interacting with the object. Because signals related to contact events and object properties are mixed in the responses of tactile fibers, our ability to disentangle these different components of our tactile experience implies that they are demultiplexed as they propagate along the neuraxis. To understand how texture and contact mechanics are encoded together by tactile fibers, we studied the activity of multiple neurons recorded simultaneously in the cervical DRG of two anesthetized rhesus monkeys while textured surfaces were applied to the glabrous skin of the fingers and palm using a handheld probe. A transducer at the tip of the textured probe measured contact forces as tactile stimuli were applied at different locations on the finger-pads and palm. We examined how a sample population of DRG neurons encode force and texture and found that firing rates of individual neurons are modulated by both force and texture. In particular, slowly adapting (SA) neurons were more responsive to force than texture, and rapidly adapting (RA) neurons were more responsive to texture than force. Although force could be decoded accurately throughout the entire contact interval, texture signals were most salient during onset and offset phases of the contact interval.NEW & NOTEWORTHY Cutaneous mechanoreceptors in our hands gather information about the objects we handle. Tactile fibers encode mixed information about contact events and object properties. Neural coding in tactile afferents is typically studied by varying a single aspect of tactile stimuli, avoiding the confounds of real-world haptic interactions. We instead record responses of small populations of DRG neurons to variable tactile stimuli and find that neurons primarily respond to force, though some texture information can be detected.

Spatial models of cell distribution in human lumbar dorsal root ganglia

Dorsal root ganglia (DRG), which contain the somata of primary sensory neurons, have increasingly been considered as novel targets for clinical neural interfaces, both for neuroprosthetic and pain applications. Effective use of either neural recording or stimulation technologies requires an appropriate spatial position relative to the target neural element, whether axon or cell body. However, the internal three-dimensional spatial organization of human DRG neural fibers and somata has not been quantitatively described. In this study, we analyzed 202 cross-sectional images across the length of 31 human L4 and L5 DRG from 10 donors. We used a custom semi-automated graphical user interface to identify the locations of neural elements in the images and normalize the output to a consistent spatial reference for direct comparison by spinal level. By applying a recursive partitioning algorithm, we found that the highest density of cell bodies at both spinal levels could be found in the inner 85% of DRG length, the outer-most 25-30% radially, and the dorsal-most 69-76%. While axonal density was fairly homogeneous across the DRG length, there was a distinct low density region in the outer 7-11% radially. These findings are consistent with previous qualitative reports of neural distribution in DRG. The quantitative measurements we provide will enable improved targeting of future neural interface technologies and DRG-focused pharmaceutical therapies, and provide a rigorous anatomical description of the bridge between the central and peripheral nervous systems.

Flexible microelectrode array for interfacing with the surface of neural ganglia

OBJECTIVE

The dorsal root ganglia (DRG) are promising nerve structures for sensory neural interfaces because they provide centralized access to primary afferent cell bodies and spinal reflex circuitry. In order to harness this potential, new electrode technologies are needed which take advantage of the unique properties of DRG, specifically the high density of neural cell bodies at the dorsal surface. Here we report initial in vivo results from the development of a flexible non-penetrating polyimide electrode array interfacing with the surface of ganglia.

APPROACH

Multiple layouts of a 64-channel iridium electrode (420 µm²) array were tested, with pitch as small as 25 µm. The buccal ganglia of invertebrate sea slug Aplysia californica were used to develop handling and recording techniques with ganglionic surface electrode arrays (GSEAs). We also demonstrated the GSEA's capability to record single- and multi-unit activity from feline lumbosacral DRG related to a variety of sensory inputs, including cutaneous brushing, joint flexion, and bladder pressure.

MAIN RESULTS

We recorded action potentials from a variety of Aplysia neurons activated by nerve stimulation, and units were observed firing simultaneously on closely spaced electrode sites. We also recorded single- and multi-unit activity associated with sensory inputs from feline DRG. We utilized spatial oversampling of action potentials on closely-spaced electrode sites to estimate the location of neural sources at between 25 µm and 107 µm below the DRG surface. We also used the high spatial sampling to demonstrate a possible spatial sensory map of one feline's DRG. We obtained activation of sensory fibers with low-amplitude stimulation through individual or groups of GSEA electrode sites.

SIGNIFICANCE

Overall, the GSEA has been shown to provide a variety of information types from ganglia neurons and to have significant potential as a tool for neural mapping and interfacing.

Quantitative models of feline lumbosacral dorsal root ganglia neuronal cell density

BACKGROUND

Dorsal root ganglia (DRG) are spinal root components that contain the cell bodies of converging primary sensory neurons. DRG are becoming a therapeutic target for electrical neural interfaces. Our purpose was to establish methods for quantifying the non-random nature and distribution of neuronal cell bodies within DRG.

NEW METHOD

We identified neuronal cell body locations in 26 feline lumbosacral DRG cross-section histological images and used computational tools to quantify spatial trends. We first analyzed spatial randomness using the nearest-neighbor distance method. Next we overlaid a 6×6 grid, modeling neuronal cellular density in each grid square and comparing regions statistically. Finally we transformed DRG onto a polar map and calculated neuronal cellular density in annular sectors. We used a recursive partition model to determine regions of high and low density, and validated the model statistically.

RESULTS

We found that the arrangement of neuronal cell bodies at the widest point of DRG is distinctly non-random with concentration in particular regions. The grid model suggested a radial trend in density, with increasing density toward the outside of the DRG. The polar transformation model showed that the highest neuronal cellular density is in the outer 23.9% radially and the dorsal ±61.4° angularly.

COMPARISON WITH EXISTING METHODS

To our knowledge, DRG neuronal cell distribution has not been previously quantified.

CONCLUSIONS

These results confirm and expand quantitatively on the existing understanding of DRG anatomy. Our methods can be useful for analyzing the distribution of cellular components of other neural structures or expanding to three-dimensional models.

Muscle spindles in human tibialis anterior encode muscle fascicle length changes

Muscle spindles provide exquisitely sensitive proprioceptive information regarding joint position and movement. Through passively driven length changes in the muscle-tendon unit (MTU), muscle spindles detect joint rotations because of their in-parallel mechanical linkage to muscle fascicles. In human microneurography studies, muscle fascicles are assumed to follow the MTU and, as such, fascicle length is not measured in such studies. However, under certain mechanical conditions, compliant structures can act to decouple the fascicles, and, therefore, the spindles, from the MTU. Such decoupling may reduce the fidelity by which muscle spindles encode joint position and movement. The aim of the present study was to measure, for the first time, both the changes in firing of single muscle spindle afferents and changes in muscle fascicle length in vivo from the tibialis anterior muscle (TA) during passive rotations about the ankle. Unitary recordings were made from 15 muscle spindle afferents supplying TA via a microelectrode inserted into the common peroneal nerve. Ultrasonography was used to measure the length of an individual fascicle of TA. We saw a strong correlation between fascicle length and firing rate during passive ankle rotations of varying rates (0.1-0.5 Hz) and amplitudes (1-9°). In particular, we saw responses observed at relatively small changes in muscle length that highlight the sensitivity of the TA muscle to small length changes. This study is the first to measure spindle firing and fascicle dynamics in vivo and provides an experimental basis for further understanding the link between fascicle length, MTU length, and spindle firing patterns.NEW & NOTEWORTHY Muscle spindles are exquisitely sensitive to changes in muscle length, but recordings from human muscle spindle afferents are usually correlated with joint angle rather than muscle fascicle length. In this study, we monitored both muscle fascicle length and spindle firing from the human tibialis anterior muscle in vivo. Our findings are the first to measure these signals in vivo and provide an experimental basis for exploring this link further.

Modulation of spinal motor output by initial arm postures in anesthetized monkeys

Proper execution of voluntary movement requires a sensorimotor transformation based on the initial limb state. For example, successfully reaching to a stable target requires the recruitment of different muscle groups depending on limb position at movement initiation. To test whether this transformation could occur at the spinal level, we stimulated the cervical spinal cord of anesthetized monkeys while systematically changing initial posture and examined the modulation of the twitch response induced in the upper limb muscles. In three monkeys, a multichannel microelectrode array was implanted into the C6 segment of the spinal cord and electromyographic electrodes were implanted in 12 limb muscles (five hand, four elbow, and three shoulder muscles). The magnitude and onset latency of the evoked response in each electrode-muscle pair were examined by systematically changing the hand position through nine positions in a horizontal plane with the monkey prone. Among 330 electrode-muscle pairs examined, 61% of pairs exhibited significant modulation of either magnitude or latency of twitch responses across different hand/arm configurations (posture dependency). We found that posture dependency occurred preferentially in the distal rather than proximal muscles and was not affected by the location of the electrode within the stimulated spinal segment. Importantly, this posture dependency was not affected by spinalization at the C2 level. These results suggest that excitability in the cervical spinal cord is affected by initial arm posture through spinal reflex pathways. This posture dependency of spinal motor output could affect voluntary arm movement by adjusting descending motor commands relative to the initial arm posture.

Decoding of the spike timing of primary afferents during voluntary arm movements in monkeys

Understanding the mechanisms of encoding forelimb kinematics in the activity of peripheral afferents is essential for developing a somatosensory neuroprosthesis. To investigate whether the spike timing of dorsal root ganglion (DRG) neurons could be estimated from the forelimb kinematics of behaving monkeys, we implanted two multi-electrode arrays chronically in the DRGs at the level of the cervical segments in two monkeys. Neuronal activity during voluntary reach-to-grasp movements were recorded simultaneously with the trajectories of hand/arm movements, which were tracked in three-dimensional space using a motion capture system. Sixteen and 13 neurons, including muscle spindles, skin receptors, and tendon organ afferents, were recorded in the two monkeys, respectively. We were able to reconstruct forelimb joint kinematics from the temporal firing pattern of a subset of DRG neurons using sparse linear regression (SLiR) analysis, suggesting that DRG neuronal ensembles encoded information about joint kinematics. Furthermore, we estimated the spike timing of the DRG neuronal ensembles from joint kinematics using an integrate-and-fire model (IF) incorporating the SLiR algorithm. The temporal change of firing frequency of a subpopulation of neurons was reconstructed precisely from forelimb kinematics using the SLiR. The estimated firing pattern of the DRG neuronal ensembles encoded forelimb joint angles and velocities as precisely as the originally recorded neuronal activity. These results suggest that a simple model can be used to generate an accurate estimate of the spike timing of DRG neuronal ensembles from forelimb joint kinematics, and is useful for designing a proprioceptive decoder in a brain machine interface.

Neuroprosthetic technology for individuals with spinal cord injury

CONTEXT

Spinal cord injury (SCI) results in a loss of function and sensation below the level of the lesion. Neuroprosthetic technology has been developed to help restore motor and autonomic functions as well as to provide sensory feedback.

FINDINGS

This paper provides an overview of neuroprosthetic technology that aims to address the priorities for functional restoration as defined by individuals with SCI. We describe neuroprostheses that are in various stages of preclinical development, clinical testing, and commercialization including functional electrical stimulators, epidural and intraspinal microstimulation, bladder neuroprosthesis, and cortical stimulation for restoring sensation. We also discuss neural recording technologies that may provide command or feedback signals for neuroprosthetic devices.

CONCLUSION/CLINICAL RELEVANCE

Neuroprostheses have begun to address the priorities of individuals with SCI, although there remains room for improvement. In addition to continued technological improvements, closing the loop between the technology and the user may help provide intuitive device control with high levels of performance.

Real-time control of hind limb functional electrical stimulation using feedback from dorsal root ganglia recordings

OBJECTIVE

Functional electrical stimulation (FES) approaches often utilize an open-loop controller to drive state transitions. The addition of sensory feedback may allow for closed-loop control that can respond effectively to perturbations and muscle fatigue.

APPROACH

We evaluated the use of natural sensory nerve signals obtained with penetrating microelectrode arrays in lumbar dorsal root ganglia (DRG) as real-time feedback for closed-loop control of FES-generated hind limb stepping in anesthetized cats.

MAIN RESULTS

Leg position feedback was obtained in near real-time at 50 ms intervals by decoding the firing rates of more than 120 DRG neurons recorded simultaneously. Over 5 m of effective linear distance was traversed during closed-loop stepping trials in each of two cats. The controller compensated effectively for perturbations in the stepping path when DRG sensory feedback was provided. The presence of stimulation artifacts and the quality of DRG unit sorting did not significantly affect the accuracy of leg position feedback obtained from the linear decoding model as long as at least 20 DRG units were included in the model.

SIGNIFICANCE

This work demonstrates the feasibility and utility of closed-loop FES control based on natural neural sensors. Further work is needed to improve the controller and electrode technologies and to evaluate long-term viability.