Selective direct motor cortical influence during naturalistic climbing.

It remains poorly resolved when and how motor cortical output directly influences limb muscle activity through descending projections, which impedes mechanistic understanding of cortical movement control. Here we addressed this in mice performing an ethologically inspired all-limb climbing behavior. We quantified the direct influence of forelimb primary motor cortex (caudal forelimb area, CFA) on muscle activity comprehensively across the muscle activity states that occur during climbing. We found that CFA informs muscle activity pattern, mainly by selectively activating certain muscles while exerting much smaller, bidirectional effects on their antagonists. From Neuropixel recordings, we identified linear combinations (components) of motor cortical activity that covary with these effects, finding that these components differ from those that covary with muscle activity or kinematics. Collectively, our results reveal an instructive direct motor cortical influence on limb muscles that is selective within a motor behavior and reliant on a new type of neural activity subspace.

Hand-jaw coordination as mice handle food is organized around intrinsic structure-function relationships.

Rodent jaws evolved structurally to support dual functionality, for either biting or chewing food. Rodent hands also function dually during food handling, for actively manipulating or statically holding food. How are these oral and manual functions coordinated? We combined electrophysiological recording of muscle activity and kilohertz kinematic tracking to analyze masseter and hand actions as mice of both sexes handled food. Masseter activity was organized into two modes synchronized to hand movement modes. In holding/chewing mode, mastication occurred as rhythmic (∼5 Hz) masseter activity while the hands held food below the mouth. In oromanual/ingestion mode, bites occurred as lower-amplitude aperiodic masseter events that were precisely timed to follow regrips (by ∼200 ms). Thus, jaw and hand movements are flexibly coordinated during food handling: uncoupled in holding/chewing mode, and tightly coordinated in oromanual/ingestion mode as regrip-bite sequences. Key features of this coordination were captured in a simple model of hierarchically orchestrated mode-switching and intra-mode action sequencing. We serendipitously detected an additional masseter-related action, tooth-sharpening, identified as bouts of higher-frequency (∼13 Hz) rhythmic masseter activity, which was accompanied by eye displacement, including rhythmic proptosis, attributable to masseter contractions. Collectively, the findings demonstrate how a natural, complex, and goal-oriented activity is organized as an assemblage of distinct modes and complex actions, adapted for the divisions of function arising from anatomical structure. These results reveal intricate, high-speed coordination of disparate effectors, and show how natural forms of dexterity can serve as a model for understanding the behavioral neurobiology of multi-body-part coordination in general.Significance statement Survival hinges on efficiently handling and ingesting food. During food handling, mice switch between statically holding and actively manipulating food with the hands, and also between biting with the incisors or chewing with the molars. Using masseter electromyography and kilohertz frame rate video, we show that the two modes of hand movement and of jaw movement are tightly linked, with holding co-occurring with chewing and handling with biting. On faster timescales, biting is precisely timed to occur after rapid sequences of manipulative hand movements. Our findings reveal fast, intricate, hierarchical coordination of hand and jaw movements, show how morphology influences the organization of complex behavior, and establish food handling as a model for neurobiological studies of multi-body-part coordination.

Hierarchy between forelimb premotor and primary motor cortices and its manifestation in their firing patterns.

Though hierarchy is commonly invoked in descriptions of motor cortical function, its presence and manifestation in firing patterns remain poorly resolved. Here we use optogenetic inactivation to demonstrate that short-latency influence between forelimb premotor and primary motor cortices is asymmetric during reaching in mice, demonstrating a partial hierarchy between the endogenous activity in each region. Multi-region recordings revealed that some activity is captured by similar but delayed patterns where either region's activity leads, with premotor activity leading more. Yet firing in each region is dominated by patterns shared between regions and is equally predictive of firing in the other region at the single-neuron level. In dual-region network models fit to data, regions differed in their dependence on across-region input, rather than the amount of such input they received. Our results indicate that motor cortical hierarchy, while present, may not be exposed when inferring interactions between populations from firing patterns alone.

Recording Forelimb Muscle Activity in Head-Fixed Mice with Chronically Implanted EMG Electrodes.

Powerful genetic and molecular tools available in mouse systems neuroscience research have enabled researchers to interrogate motor system function with unprecedented precision in head-fixed mice performing a variety of tasks. The small size of the mouse makes the measurement of motor output difficult, as the traditional method of electromyographic (EMG) recording of muscle activity was designed for larger animals like cats and primates. Pending commercially available EMG electrodes for mice, the current gold-standard method for recording muscle activity in mice is to make electrode sets in-house. This article describes a refinement of established procedures for hand fabrication of an electrode set, implantation of electrodes in the same surgery as headplate implantation, fixation of a connector on the headplate, and post-operative recovery care. Following recovery, millisecond-resolution EMG recordings can be obtained during head-fixed behavior for several weeks without noticeable changes in signal quality. These recordings enable precise measurement of forelimb muscle activity alongside in vivo neural recording and/or perturbation to probe mechanisms of motor control in mice.

Motor cortical influence relies on task-specific activity covariation.

During limb movement, spinal circuits facilitate the alternating activation of antagonistic flexor and extensor muscles. Yet antagonist cocontraction is often required to stabilize joints, like when loads are handled. Previous results suggest that these different muscle activation patterns are mediated by separate flexion- and extension-related motor cortical output populations, while others suggest recruitment of task-specific populations. To distinguish between hypotheses, we developed a paradigm in which mice toggle between forelimb tasks requiring antagonist alternation or cocontraction and measured activity in motor cortical layer 5b. Our results conform to neither hypothesis: consistent flexion- and extension-related activity is not observed across tasks, and no task-specific populations are observed. Instead, activity covariation among motor cortical neurons dramatically changes between tasks, thereby altering the relation between neural and muscle activity. This is also observed specifically for corticospinal neurons. Collectively, our findings indicate that motor cortex drives different muscle activation patterns via task-specific activity covariation.

Oculomotor plant and neural dynamics suggest gaze control requires integration on distributed timescales.

A fundamental principle of biological motor control is that the neural commands driving movement must conform to the response properties of the motor plants they control. In the oculomotor system, characterizations of oculomotor plant dynamics traditionally supported models in which the plant responds to neural drive to extraocular muscles on exclusively short, subsecond timescales. These models predict that the stabilization of gaze during fixations between saccades requires neural drive that approximates eye position on longer timescales and is generated through the temporal integration of brief eye velocity-encoding signals that cause saccades. However, recent measurements of oculomotor plant behaviour have revealed responses on longer timescales. Furthermore, measurements of firing patterns in the oculomotor integrator have revealed a more complex encoding of eye movement dynamics. Yet, the link between these observations has remained unclear. Here we use measurements from the larval zebrafish to link dynamics in the oculomotor plant to dynamics in the neural integrator. The oculomotor plant in both anaesthetized and awake larval zebrafish was characterized by a broad distribution of response timescales, including those much longer than 1 s. Analysis of the firing patterns of oculomotor integrator neurons, which exhibited a broadly distributed range of decay time constants, demonstrates the sufficiency of this activity for stabilizing gaze given an oculomotor plant with distributed response timescales. This work suggests that leaky integration on multiple, distributed timescales by the oculomotor integrator reflects an inverse model for generating oculomotor commands, and that multi-timescale dynamics may be a general feature of motor circuitry. KEY POINTS: Recent observations of oculomotor plant response properties and neural activity across the oculomotor system have called into question classical formulations of both the oculomotor plant and the oculomotor integrator. Here we use measurements from new and published experiments in the larval zebrafish together with modelling to reconcile recent oculomotor plant observations with oculomotor integrator function. We developed computational techniques to characterize oculomotor plant responses over several seconds in awake animals, demonstrating that long timescale responses seen in anaesthetized animals extend to the awake state. Analysis of firing patterns of oculomotor integrator neurons demonstrates the sufficiency of this activity for stabilizing gaze given an oculomotor plant with multiple, distributed response timescales. Our results support a formulation of gaze stabilization by the oculomotor system in which commands for stabilizing gaze are generated through integration on multiple, distributed timescales.

Towards Cell and Subtype Resolved Functional Organization: Mouse as a Model for the Cortical Control of Movement.

Despite a long history of interrogation, the functional organization of motor cortex remains obscure. A major barrier has been the inability to measure and perturb activity with sufficient resolution to reveal clear functional elements within motor cortex and its associated circuits. Increasingly, the mouse has been employed as a model to facilitate application of contemporary approaches with the potential to surmount this barrier. In this brief essay, we consider these approaches and their use in the context of studies aimed at resolving the logic of motor cortical operation.

Motor Cortex Embeds Muscle-like Commands in an Untangled Population Response.

Primate motor cortex projects to spinal interneurons and motoneurons, suggesting that motor cortex activity may be dominated by muscle-like commands. Observations during reaching lend support to this view, but evidence remains ambiguous and much debated. To provide a different perspective, we employed a novel behavioral paradigm that facilitates comparison between time-evolving neural and muscle activity. We found that single motor cortex neurons displayed many muscle-like properties, but the structure of population activity was not muscle-like. Unlike muscle activity, neural activity was structured to avoid "tangling": moments where similar activity patterns led to dissimilar future patterns. Avoidance of tangling was present across tasks and species. Network models revealed a potential reason for this consistent feature: low tangling confers noise robustness. Finally, we were able to predict motor cortex activity from muscle activity by leveraging the hypothesis that muscle-like commands are embedded in additional structure that yields low tangling.

Behaviorally Selective Engagement of Short-Latency Effector Pathways by Motor Cortex.

Blocking motor cortical output with lesions or pharmacological inactivation has identified movements that require motor cortex. Yet, when and how motor cortex influences muscle activity during movement execution remains unresolved. We addressed this ambiguity using measurement and perturbation of motor cortical activity together with electromyography in mice during two forelimb movements that differ in their requirement for cortical involvement. Rapid optogenetic silencing and electrical stimulation indicated that short-latency pathways linking motor cortex with spinal motor neurons are selectively activated during one behavior. Analysis of motor cortical activity revealed a dramatic change between behaviors in the coordination of firing patterns across neurons that could account for this differential influence. Thus, our results suggest that changes in motor cortical output patterns enable a behaviorally selective engagement of short-latency effector pathways. The model of motor cortical influence implied by our findings helps reconcile previous observations on the function of motor cortex.

Spinal Inhibitory Interneuron Diversity Delineates Variant Motor Microcircuits.

Animals generate movement by engaging spinal circuits that direct precise sequences of muscle contraction, but the identity and organizational logic of local interneurons that lie at the core of these circuits remain unresolved. Here, we show that V1 interneurons, a major inhibitory population that controls motor output, fractionate into highly diverse subsets on the basis of the expression of 19 transcription factors. Transcriptionally defined V1 subsets exhibit distinct physiological signatures and highly structured spatial distributions with mediolateral and dorsoventral positional biases. These positional distinctions constrain patterns of input from sensory and motor neurons and, as such, suggest that interneuron position is a determinant of microcircuit organization. Moreover, V1 diversity indicates that different inhibitory microcircuits exist for motor pools controlling hip, ankle, and foot muscles, revealing a variable circuit architecture for interneurons that control limb movement.

Primacy of Flexor Locomotor Pattern Revealed by Ancestral Reversion of Motor Neuron Identity.

Spinal circuits can generate locomotor output in the absence of sensory or descending input, but the principles of locomotor circuit organization remain unclear. We sought insight into these principles by considering the elaboration of locomotor circuits across evolution. The identity of limb-innervating motor neurons was reverted to a state resembling that of motor neurons that direct undulatory swimming in primitive aquatic vertebrates, permitting assessment of the role of motor neuron identity in determining locomotor pattern. Two-photon imaging was coupled with spike inference to measure locomotor firing in hundreds of motor neurons in isolated mouse spinal cords. In wild-type preparations, we observed sequential recruitment of motor neurons innervating flexor muscles controlling progressively more distal joints. Strikingly, after reversion of motor neuron identity, virtually all firing patterns became distinctly flexor like. Our findings show that motor neuron identity directs locomotor circuit wiring and indicate the evolutionary primacy of flexor pattern generation.

Prolonged, brain-wide expression of nuclear-localized GCaMP3 for functional circuit mapping.

Larval zebrafish offer the potential for large-scale optical imaging of neural activity throughout the central nervous system; however, several barriers challenge their utility. First, ~panneuronal probe expression has to date only been demonstrated at early larval stages up to 7 days post-fertilization (dpf), precluding imaging at later time points when circuits are more mature. Second, nuclear exclusion of genetically-encoded calcium indicators (GECIs) limits the resolution of functional fluorescence signals collected during imaging. Here, we report the creation of transgenic zebrafish strains exhibiting robust, nuclearly targeted expression of GCaMP3 across the brain up to at least 14 dpf utilizing a previously described optimized Gal4-UAS system. We confirmed both nuclear targeting and functionality of the modified probe in vitro and measured its kinetics in response to action potentials (APs). We then demonstrated in vivo functionality of nuclear-localized GCaMP3 in transgenic zebrafish strains by identifying eye position-sensitive fluorescence fluctuations in caudal hindbrain neurons during spontaneous eye movements. Our methodological approach will facilitate studies of larval zebrafish circuitry by both improving resolution of functional Ca(2+) signals and by allowing brain-wide expression of improved GECIs, or potentially any probe, further into development.

Edging toward entelechy in motor control.

The organization and functional logic of corticospinal motor neurons and their target connections remains unclear, despite their evident influence on movement. Spinal interneurons mediate much of this influence, yet we know little about the way in which corticospinal neurons engage spinal interneurons. This is perhaps not surprising given that the principles of organization of local spinal microcircuits remain elusive--we have glimpses of an underlying order but lack a comprehensive view of their functional architecture. In this brief essay we make a case that a new focus on the intersection of cortical and spinal circuits may provide clarity to the interpretation of corticospinal motor neuron firing patterns and help specify the logic of corticospinal motor neuronal function.

A dynamical feedback model for adaptation in the olfactory transduction pathway.

Olfactory transduction exhibits two distinct types of adaptation, which we denote multipulse and step adaptation. In terms of measured transduction current, multipulse adaptation appears as a decrease in the amplitude of the second of two consecutive responses when the olfactory neuron is stimulated with two brief pulses. Step adaptation occurs in response to a sustained steplike stimulation and is characterized by a return to a steady-state current amplitude close to the prestimulus value, after a transient peak. In this article, we formulate a dynamical model of the olfactory transduction pathway, which includes the kinetics of the CNG channels, the concentration of Ca ions flowing through them, and the Ca-complexes responsible for the regulation. Based on this model, a common dynamical explanation for the two types of adaptation is suggested. We show that both forms of adaptation can be well described using different time constants for the kinetics of Ca ions (faster) and the kinetics of the feedback mechanisms (slower). The model is validated on experimental data collected in voltage-clamp conditions using different techniques and animal species.

Spatial gradients and multidimensional dynamics in a neural integrator circuit.

In a neural integrator, the variability and topographical organization of neuronal firing-rate persistence can provide information about the circuit's functional architecture. We used optical recording to measure the time constant of decay of persistent firing (persistence time) across a population of neurons comprising the larval zebrafish oculomotor velocity-to-position neural integrator. We found extensive persistence time variation (tenfold; coefficients of variation = 0.58-1.20) across cells in individual larvae. We also found that the similarity in firing between two neurons decreased as the distance between them increased and that a gradient in persistence time was mapped along the rostrocaudal and dorsoventral axes. This topography is consistent with the emergence of persistence time heterogeneity from a circuit architecture in which nearby neurons are more strongly interconnected than distant ones. Integrator circuit models characterized by multiple dimensions of slow firing-rate dynamics can account for our results.

Regression-based identification of behavior-encoding neurons during large-scale optical imaging of neural activity at cellular resolution.

The advent of methods for optical imaging of large-scale neural activity at cellular resolution in behaving animals presents the problem of identifying behavior-encoding cells within the resulting image time series. Rapid and precise identification of cells with particular neural encoding would facilitate targeted activity measurements and perturbations useful in characterizing the operating principles of neural circuits. Here we report a regression-based approach to semiautomatically identify neurons that is based on the correlation of fluorescence time series with quantitative measurements of behavior. The approach is illustrated with a novel preparation allowing synchronous eye tracking and two-photon laser scanning fluorescence imaging of calcium changes in populations of hindbrain neurons during spontaneous eye movement in the larval zebrafish. Putative velocity-to-position oculomotor integrator neurons were identified that showed a broad spatial distribution and diversity of encoding. Optical identification of integrator neurons was confirmed with targeted loose-patch electrical recording and laser ablation. The general regression-based approach we demonstrate should be widely applicable to calcium imaging time series in behaving animals.

All-trans-retinal is a closed-state inhibitor of rod cyclic nucleotide-gated ion channels.

Rod vision begins when 11-cis-retinal absorbs a photon and isomerizes to all-trans-retinal (ATR) within the photopigment, rhodopsin. Photoactivated rhodopsin triggers an enzyme cascade that lowers the concentration of cGMP, thereby closing cyclic nucleotide-gated (CNG) ion channels. After isomerization, ATR dissociates from rhodopsin, and after a bright light, this release is expected to produce a large surge of ATR near the CNG channels. Using excised patches from Xenopus oocytes, we recently showed that ATR shuts down cloned rod CNG channels, and that this inhibition occurs in the nanomolar range (aqueous concentration) at near-physiological concentrations of cGMP. Here we further characterize the ATR effect and present mechanistic information. ATR was found to decrease the apparent cGMP affinity, as well as the maximum current at saturating cGMP. When ATR was applied to outside-out patches, inhibition was much slower and less effective than when it was applied to inside-out patches, suggesting that ATR requires access to the intracellular surface of the channel or membrane. The apparent ATR affinity and maximal inhibition of heteromeric (CNGA1/CNGB1) channels was similar to that of homomeric (CNGA1) channels. Single-channel and multichannel data suggest that channel inhibition by ATR is reversible. Inhibition by ATR was not voltage dependent, and the form of its dose-response relation suggested multiple ATR molecules interacting per channel. Modeling of the data obtained with cAMP and cGMP suggests that ATR acts by interfering with the allosteric opening transition of the channel and that it prefers closed, unliganded channels. It remains to be determined whether ATR acts directly on the channel protein or instead alters channel-bilayer interactions.

All-trans-retinal shuts down rod cyclic nucleotide-gated ion channels: a novel role for photoreceptor retinoids in the response to bright light?

In retinal rods, light-induced isomerization of 11-cis-retinal to all-trans-retinal within rhodopsin triggers an enzyme cascade that lowers the concentration of cGMP. Consequently, cyclic nucleotide-gated (CNG) ion channels close, generating the first electrical response to light. After isomerization, all-trans-retinal dissociates from rhodopsin. We now show that all-trans-retinal directly and markedly inhibits cloned rod CNG channels in excised patches. 11-cis-retinal and all-trans-retinol also inhibited the channels, but at somewhat higher concentrations. Single-channel analysis suggests that all-trans-retinal reduces average open probability of rod CNG channels by inactivating channels for seconds at a time. At physiological cGMP levels, all-trans-retinal inhibited in the nanomolar range. Our results suggest that all-trans-retinal may be a potent regulator of the channel in rods during the response to bright light, when there is a large surge in the concentration of all-trans-retinal.