skip to main content


This content will become publicly available on June 12, 2024

Title: An information theoretic method to resolve millisecond-scale spike timing precision in a comprehensive motor program
Sensory inputs in nervous systems are often encoded at the millisecond scale in a precise spike timing code. There is now growing evidence in behaviors ranging from slow breathing to rapid flight for the prevalence of precise timing encoding in motor systems. Despite this, we largely do not know at what scale timing matters in these circuits due to the difficulty of recording a complete set of spike-resolved motor signals and assessing spike timing precision for encoding continuous motor signals. We also do not know if the precision scale varies depending on the functional role of different motor units. We introduce a method to estimate spike timing precision in motor circuits using continuous MI estimation at increasing levels of added uniform noise. This method can assess spike timing precision at fine scales for encoding rich motor output variation. We demonstrate the advantages of this approach compared to a previously established discrete information theoretic method of assessing spike timing precision. We use this method to analyze the precision in a nearly complete, spike resolved recording of the 10 primary wing muscles control flight in an agile hawk moth, Manduca sexta . Tethered moths visually tracked a robotic flower producing a range of turning (yaw) torques. We know that all 10 muscles in this motor program encode the majority of information about yaw torque in spike timings, but we do not know whether individual muscles encode motor information at different levels of precision. We demonstrate that the scale of temporal precision in all motor units in this insect flight circuit is at the sub-millisecond or millisecond-scale, with variation in precision scale present between muscle types. This method can be applied broadly to estimate spike timing precision in sensory and motor circuits in both invertebrates and vertebrates.  more » « less
Award ID(s):
1806833
NSF-PAR ID:
10459089
Author(s) / Creator(s):
; ; ; ;
Editor(s):
Latham, Peter E.
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
19
Issue:
6
ISSN:
1553-7358
Page Range / eLocation ID:
e1011170
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Traditional models of motor control typically operate in the domain of continuous signals such as spike rates, forces, and kinematics. However, there is growing evidence that precise spike timings encode significant information that coordinates and causally influences motor control. Some existing neural network models incorporate spike timing precision but they neither predict motor spikes coordinated across multiple motor units nor capture sensory-driven modulation of agile locomotor control. In this paper, we propose a visual encoder and model of a sensorimotor system based on a recurrent neural network (RNN) that utilizes spike timing encoding during smooth pursuit target tracking. We use this to predict a nearly complete, spike-resolved motor program of a hawkmoth that requires coordinated millisecond precision across 10 major flight motor units. Each motor unit enervates one muscle and utilizes both rate and timing encoding. Our model includes a motion detection mechanism inspired by the hawkmoth's compound eye, a convolutional encoder that compresses the sensory input, and a simple RNN that is sufficient to sequentially predict wingstroke-to-wingstroke modulation in millisecond-precise spike timings. The two-layer output architecture of the RNN separately predicts the occurrence and timing of each spike in the motor program. The dataset includes spikes recorded from all motor units during a tethered flight where the hawkmoth attends to a moving robotic flower, with a total of roughly 7000 wingstrokes from 16 trials on 5 hawkmoth subjects. Intra-trial and same-subject inter-trial predictions on the test data show that nearly every spike can be predicted within 2 ms of its known spike timing precision values. Whereas, spike occurrence prediction accuracy is about 90%. Overall, our model can predict the precise spike timing of a nearly complete motor program for hawkmoth flight with a precision comparable to that seen in agile flying insects. Such an encoding framework that captures visually-modulated precise spike timing codes and coordination can reveal how organisms process visual cues for agile movements. It can also drive the next generation of neuromorphic controllers for navigation in complex environments. 
    more » « less
  2. Sequences of action potentials, or spikes, carry information in the number of spikes and their timing. Spike timing codes are critical in many sensory systems, but there is now growing evidence that millisecond-scale changes in timing also carry information in motor brain regions, descending decision-making circuits, and individual motor units. Across all of the many signals that control a behavior, how ubiquitous, consistent, and coordinated are spike timing codes? Assessing these open questions ideally involves recording across the whole motor program with spike-level resolution. To do this, we took advantage of the relatively few motor units controlling the wings of a hawk moth, Manduca sexta . We simultaneously recorded nearly every action potential from all major wing muscles and the resulting forces in tethered flight. We found that timing encodes more information about turning behavior than spike count in every motor unit, even though there is sufficient variation in count alone. Flight muscles vary broadly in function as well as in the number and timing of spikes. Nonetheless, each muscle with multiple spikes consistently blends spike timing and count information in a 3:1 ratio. Coding strategies are consistent. Finally, we assess the coordination of muscles using pairwise redundancy measured through interaction information. Surprisingly, not only are all muscle pairs coordinated, but all coordination is accomplished almost exclusively through spike timing, not spike count. Spike timing codes are ubiquitous, consistent, and essential for coordination. 
    more » « less
  3. null (Ed.)
    Animals rapidly collect and act on incoming information to navigate complex environments, making the precise timing of sensory feedback critical in the context of neural circuit function. Moreover, the timing of sensory input determines the biomechanical properties of muscles that undergo cyclic length changes, as during locomotion. Both of these issues come to a head in the case of flying insects, as these animals execute steering manoeuvres at timescales approaching the upper limits of performance for neuromechanical systems. Among insects, flies stand out as especially adept given their ability to execute manoeuvres that require sub-millisecond control of steering muscles. Although vision is critical, here I review the role of rapid, wingbeat-synchronous mechanosensory feedback from the wings and structures unique to flies, the halteres. The visual system and descending interneurons of the brain employ a spike rate coding scheme to relay commands to the wing steering system. By contrast, mechanosensory feedback operates at faster timescales and in the language of motor neurons, i.e. spike timing, allowing wing and haltere input to dynamically structure the output of the wing steering system. Although the halteres have been long known to provide essential input to the wing steering system as gyroscopic sensors, recent evidence suggests that the feedback from these vestigial hindwings is under active control. Thus, flies may accomplish manoeuvres through a conserved hindwing circuit, regulating the firing phase—and thus, the mechanical power output—of the wing steering muscles. 
    more » « less
  4. Morrison, Abigail (Ed.)
    Assessing directional influences between neurons is instrumental to understand how brain circuits process information. To this end, Granger causality, a technique originally developed for time-continuous signals, has been extended to discrete spike trains. A fundamental assumption of this technique is that the temporal evolution of neuronal responses must be due only to endogenous interactions between recorded units, including self-interactions. This assumption is however rarely met in neurophysiological studies, where the response of each neuron is modulated by other exogenous causes such as, for example, other unobserved units or slow adaptation processes. Here, we propose a novel point-process Granger causality technique that is robust with respect to the two most common exogenous modulations observed in real neuronal responses: within-trial temporal variations in spiking rate and between-trial variability in their magnitudes. This novel method works by explicitly including both types of modulations into the generalized linear model of the neuronal conditional intensity function (CIF). We then assess the causal influence of neuron i onto neuron j by measuring the relative reduction of neuron j ’s point process likelihood obtained considering or removing neuron i . CIF’s hyper-parameters are set on a per-neuron basis by minimizing Akaike’s information criterion. In synthetic data sets, generated by means of random processes or networks of integrate-and-fire units, the proposed method recovered with high accuracy, sensitivity and robustness the underlying ground-truth connectivity pattern. Application of presently available point-process Granger causality techniques produced instead a significant number of false positive connections. In real spiking responses recorded from neurons in the monkey pre-motor cortex (area F5), our method revealed many causal relationships between neurons as well as the temporal structure of their interactions. Given its robustness our method can be effectively applied to real neuronal data. Furthermore, its explicit estimate of the effects of unobserved causes on the recorded neuronal firing patterns can help decomposing their temporal variations into endogenous and exogenous components. 
    more » « less
  5. Abstract

    Objective.Characterizing the relationship between neuron spiking and the signals that electrodes record is vital to defining the neural circuits driving brain function and informing clinical brain-machine interface design. However, high electrode biocompatibility and precisely localizing neurons around the electrodes are critical to defining this relationship.Approach.Here, we demonstrate consistent localization of the recording site tips of subcellular-scale (6.8µm diameter) carbon fiber electrodes and the positions of surrounding neurons. We implanted male rats with carbon fiber electrode arrays for 6 or 12+ weeks targeting layer V motor cortex. After explanting the arrays, we immunostained the implant site and localized putative recording site tips with subcellular-cellular resolution. We then 3D segmented neuron somata within a 50µm radius from implanted tips to measure neuron positions and health and compare to healthy cortex with symmetric stereotaxic coordinates.Main results.Immunostaining of astrocyte, microglia, and neuron markers confirmed that overall tissue health was indicative of high biocompatibility near the tips. While neurons near implanted carbon fibers were stretched, their number and distribution were similar to hypothetical fibers placed in healthy contralateral brain. Such similar neuron distributions suggest that these minimally invasive electrodes demonstrate the potential to sample naturalistic neural populations. This motivated the prediction of spikes produced by nearby neurons using a simple point source model fit using recorded electrophysiology and the mean positions of the nearest neurons observed in histology. Comparing spike amplitudes suggests that the radius at which single units can be distinguished from others is near the fourth closest neuron (30.7 ± 4.6µm,Xˉ± S) in layer V motor cortex.Significance.Collectively, these data and simulations provide the first direct evidence that neuron placement in the immediate vicinity of the recording site influences how many spike clusters can be reliably identified by spike sorting.

     
    more » « less