Scaling of Motor Output, From Mouse to Humans

Effect of size on expression of bistability in mouse spinal motoneurons

Bistability in spinal motoneurons supports tonic spike activity in the absence of excitatory drive. Earlier work in adult preparations suggested that smaller motoneurons innervating slow antigravity muscle fibers are more likely to generate bistability for postural maintenance. However, whether large motoneurons innervating fast-fatigable muscle fibers display bistability is still controversial. To address this, we examined the relationship between soma size and bistability in lumbar (L4-L5) ventrolateral α-motoneurons of choline acetyltransferase (ChAT)-green fluorescent protein (GFP) and Hb9-GFP mice during the first 4 wk of life. We found that as neuron size increases, the prevalence of bistability rises. Smaller α-motoneurons lack bistability, whereas larger fast α-motoneurons [matrix metalloproteinase-9 (MMP-9)+/Hb9+] with a soma area ≥ 400 µm2 exhibit significantly higher bistability. Ionic currents associated with bistability, including the persistent Nav1.6 current, the thermosensitive Trpm5 Ca2+-activated Na+ current, and the slowly inactivating Kv1.2 current, also scale with cell size. Serotonin evokes full bistability in large motoneurons with partial bistable properties but not in small motoneurons. Our study provides important insights into the neural mechanisms underlying bistability and how motoneuron size correlates with bistability in mice.NEW & NOTEWORTHY Bistability is not a common feature of all mouse spinal motoneurons. It is absent in small, slow motoneurons but present in most large, fast motoneurons. This difference results from differential expression of ionic currents that enable bistability, which are highly expressed in large motoneurons but small or absent in small motoneurons. These results support a possible role for fast motoneurons in maintenance of tonic posture in addition to their known roles in fast movements.

Myomatrix arrays for high-definition muscle recording

Neurons coordinate their activity to produce an astonishing variety of motor behaviors. Our present understanding of motor control has grown rapidly thanks to new methods for recording and analyzing populations of many individual neurons over time. In contrast, current methods for recording the nervous system's actual motor output - the activation of muscle fibers by motor neurons - typically cannot detect the individual electrical events produced by muscle fibers during natural behaviors and scale poorly across species and muscle groups. Here we present a novel class of electrode devices ('Myomatrix arrays') that record muscle activity at unprecedented resolution across muscles and behaviors. High-density, flexible electrode arrays allow for stable recordings from the muscle fibers activated by a single motor neuron, called a 'motor unit,' during natural behaviors in many species, including mice, rats, primates, songbirds, frogs, and insects. This technology therefore allows the nervous system's motor output to be monitored in unprecedented detail during complex behaviors across species and muscle morphologies. We anticipate that this technology will allow rapid advances in understanding the neural control of behavior and identifying pathologies of the motor system.

Motoneuron-driven computational muscle modelling with motor unit resolution and subject-specific musculoskeletal anatomy

The computational simulation of human voluntary muscle contraction is possible with EMG-driven Hill-type models of whole muscles. Despite impactful applications in numerous fields, the neuromechanical information and the physiological accuracy such models provide remain limited because of multiscale simplifications that limit comprehensive description of muscle internal dynamics during contraction. We addressed this limitation by developing a novel motoneuron-driven neuromuscular model, that describes the force-generating dynamics of a population of individual motor units, each of which was described with a Hill-type actuator and controlled by a dedicated experimentally derived motoneuronal control. In forward simulation of human voluntary muscle contraction, the model transforms a vector of motoneuron spike trains decoded from high-density EMG signals into a vector of motor unit forces that sum into the predicted whole muscle force. The motoneuronal control provides comprehensive and separate descriptions of the dynamics of motor unit recruitment and discharge and decodes the subject's intention. The neuromuscular model is subject-specific, muscle-specific, includes an advanced and physiological description of motor unit activation dynamics, and is validated against an experimental muscle force. Accurate force predictions were obtained when the vector of experimental neural controls was representative of the discharge activity of the complete motor unit pool. This was achieved with large and dense grids of EMG electrodes during medium-force contractions or with computational methods that physiologically estimate the discharge activity of the motor units that were not identified experimentally. This neuromuscular model advances the state-of-the-art of neuromuscular modelling, bringing together the fields of motor control and musculoskeletal modelling, and finding applications in neuromuscular control and human-machine interfacing research.

Combining electromyographic and electrical impedance data sets through machine learning: A study in D2-mdx and wild-type mice

INTRODUCTION/AIMS

Needle impedance-electromyography (iEMG) assesses the active and passive electrical properties of muscles concurrently by using a novel needle with six electrodes, two for EMG and four for electrical impedance myography (EIM). Here, we assessed an approach for combining multifrequency EMG and EIM data via machine learning (ML) to discriminate D2-mdx muscular dystrophy and wild-type (WT) mouse skeletal muscle.

METHODS

iEMG data were obtained from quadriceps of D2-mdx mice, a muscular dystrophy model, and WT animals. EIM data were collected with the animals under deep anesthesia and EMG data collected under light anesthesia, allowing for limited spontaneous movement. Fourier transformation was performed on the EMG data to provide power spectra that were sampled across the frequency range using three different approaches. Random forest-based, nested ML was applied to the EIM and EMG data sets separately and then together to assess healthy versus disease category classification using a nested cross-validation procedure.

RESULTS

Data from 20 D2-mdx and 20 WT limbs were analyzed. EIM data fared better than EMG data in differentiating healthy from disease mice with 93.1% versus 75.6% accuracy, respectively. Combining EIM and EMG data sets yielded similar performance as EIM data alone with 92.2% accuracy.

DISCUSSION

We have demonstrated an ML-based approach for combining EIM and EMG data obtained with an iEMG needle. While EIM-EMG in combination fared no better than EIM alone with this data set, the approach used here demonstrates a novel method of combining the two techniques to characterize the full electrical properties of skeletal muscle.

Myomatrix arrays for high-definition muscle recording

Neurons coordinate their activity to produce an astonishing variety of motor behaviors. Our present understanding of motor control has grown rapidly thanks to new methods for recording and analyzing populations of many individual neurons over time. In contrast, current methods for recording the nervous system's actual motor output - the activation of muscle fibers by motor neurons - typically cannot detect the individual electrical events produced by muscle fibers during natural behaviors and scale poorly across species and muscle groups. Here we present a novel class of electrode devices ("Myomatrix arrays") that record muscle activity at unprecedented resolution across muscles and behaviors. High-density, flexible electrode arrays allow for stable recordings from the muscle fibers activated by a single motor neuron, called a "motor unit", during natural behaviors in many species, including mice, rats, primates, songbirds, frogs, and insects. This technology therefore allows the nervous system's motor output to be monitored in unprecedented detail during complex behaviors across species and muscle morphologies. We anticipate that this technology will allow rapid advances in understanding the neural control of behavior and in identifying pathologies of the motor system.

A cellular taxonomy of the adult human spinal cord

The mammalian spinal cord functions as a community of cell types for sensory processing, autonomic control, and movement. While animal models have advanced our understanding of spinal cellular diversity, characterizing human biology directly is important to uncover specialized features of basic function and human pathology. Here, we present a cellular taxonomy of the adult human spinal cord using single-nucleus RNA sequencing with spatial transcriptomics and antibody validation. We identified 29 glial clusters and 35 neuronal clusters, organized principally by anatomical location. To demonstrate the relevance of this resource to human disease, we analyzed spinal motoneurons, which degenerate in amyotrophic lateral sclerosis (ALS) and other diseases. We found that compared with other spinal neurons, human motoneurons are defined by genes related to cell size, cytoskeletal structure, and ALS, suggesting a specialized molecular repertoire underlying their selective vulnerability. We include a web resource to facilitate further investigations into human spinal cord biology.

Local and systemic transcriptomic responses from acute exercise induced muscle damage of the human knee extensors

Skeletal muscle is adaptable to a direct stimulus of exercise-induced muscle damage (EIMD). Local muscle gene networks and systemic circulatory factors respond to EIMD within days, mediating anti-inflammation and cellular proliferation. Here we show in humans that local EIMD of one muscle group is associated with a systemic response of gene networks that regulate muscle structure and cellular development in non-local homologous muscle not directly altered by EIMD. In the non-dominant knee-extensors of seven males, EIMD was induced through voluntary contractions against an electric motor that lengthened muscles. Neuromuscular assessments, vastus lateralis muscle biopsies and blood draws occurred at two days prior, and one and two days post the EIMD intervention. From the muscle and blood plasma samples, RNA-seq measured transcriptome changes of differential expression using bioinformatic analyses.Relative to the time of the EIMD intervention, local muscle that was mechanically damaged had 475 genes differentially expressed, as compared to 33 genes in the non-local homologous muscle. Gene and network analysis showed that activity of the local muscle was related to structural maintenance, repair, and energetic processes, whereas gene and network activity of the non-local muscle (that was not directly modified by the EIMD) were related to muscle cell development, stress response, and structural maintenance. Altered expression of two novel miRNAs related to the EIMD response supported that systemic factors were active. Together, these results indicate that the expression of genes and gene networks that control muscle contractile structure can be modified in response to non-local EIMD in humans.

Distribution of motor unit properties across human muscles

The purpose of our review was to compare the distribution of motor unit properties across human muscles of different sizes and recruitment ranges. Although motor units can be distinguished based on several different attributes, we focused on four key parameters that have a significant influence on the force produced by muscle during voluntary contractions: the number of motor units, average innervation number, the distributions of contractile characteristics, and discharge rates within motor unit pools. Despite relatively few publications on this topic, current data indicate that the most influential factor in the distribution of these motor unit properties between muscles is innervation number. Nonetheless, despite a fivefold difference in innervation number between a hand muscle (first dorsal interosseus) and a lower leg muscle (tibialis anterior), the general organization of their motor unit pools, and the range of discharge rates appear to be relatively similar. These observations provide foundational knowledge for studies on the control of movement and the changes that occur with aging and neurological disorders.

Diversity of Mammalian Motoneurons and Motor Units

Although they share the common function of controlling muscle fiber contraction, spinal motoneurons display a remarkable diversity. Alpha-motoneurons are the "final common pathway", which relay all the information from spinal and supraspinal centers and allow the organism to interact with the outside world by controlling the contraction of muscle fibers in the muscles. On the other hand, gamma-motoneurons are specialized motoneurons that do not generate force and instead specifically innervate muscle fibers inside muscle spindles, which are proprioceptive organs embedded in the muscles. Beta-motoneurons are hybrid motoneurons that innervate both extrafusal and intrafusal muscle fibers. Even among alpha-motoneurons, there exists an exquisite diversity in terms of motoneuron electrical and molecular properties, physiological and structural properties of their neuromuscular junctions, and molecular and contractile properties of the innervated muscle fibers. This diversity, across species, across muscles, and across muscle fibers in a given muscle, underlie the vast repertoire of movements that one individual can perform.

Control of Mammalian Locomotion by Somatosensory Feedback

When animals walk overground, mechanical stimuli activate various receptors located in muscles, joints, and skin. Afferents from these mechanoreceptors project to neuronal networks controlling locomotion in the spinal cord and brain. The dynamic interactions between the control systems at different levels of the neuraxis ensure that locomotion adjusts to its environment and meets task demands. In this article, we describe and discuss the essential contribution of somatosensory feedback to locomotion. We start with a discussion of how biomechanical properties of the body affect somatosensory feedback. We follow with the different types of mechanoreceptors and somatosensory afferents and their activity during locomotion. We then describe central projections to locomotor networks and the modulation of somatosensory feedback during locomotion and its mechanisms. We then discuss experimental approaches and animal models used to investigate the control of locomotion by somatosensory feedback before providing an overview of the different functional roles of somatosensory feedback for locomotion. Lastly, we briefly describe the role of somatosensory feedback in the recovery of locomotion after neurological injury. We highlight the fact that somatosensory feedback is an essential component of a highly integrated system for locomotor control. © 2021 American Physiological Society. Compr Physiol 11:1-71, 2021.

Estimates of persistent inward currents increase with the level of voluntary drive in low-threshold motor units of plantar flexor muscles

This study tested whether estimates of persistent inward currents (PICs) in the human plantar flexors would increase with the level of voluntary drive. High-density surface electromyograms were collected from soleus and gastrocnemius medialis of 21 participants (29.2 ± 2.6 yr) during ramp-shaped isometric contractions to 10%, 20%, and 30% (torque rise and decline of 2%/s and 30-s duration) of each participant's maximal torque. Motor units identified in all the contraction intensities were included in the paired-motor unit analysis to calculate delta frequency (ΔF) and estimate the PICs. ΔF is the difference in discharge rate of the control unit at the time of recruitment and derecruitment of the test unit. Increases in PICs were observed from 10% to 20% [Δ = 0.6 pulse per second (pps); P < 0.001] and from 20% to 30% (Δ = 0.5 pps; P < 0.001) in soleus and from 10% to 20% (Δ = 1.2 pps; P < 0.001) but not from 20% to 30% (Δ = 0.09 pps; P = 0.724) in gastrocnemius medialis. Maximal discharge rate increased for soleus and gastrocnemius medialis from 10% to 20% [Δ = 1.75 pps (P < 0.001) and Δ = 2.43 pps (P < 0.001), respectively] and from 20% to 30% [Δ = 0.80 pps (P < 0.017) and Δ = 0.92 pps (P = 0.002), respectively]. The repeated-measures correlation identified associations between ΔF and increases in maximal discharge rate for soleus (r = 0.64; P < 0.001) and gastrocnemius medialis (r = 0.77; P < 0.001). An increase in voluntary drive tends to increase PIC strength, which has key implications for the control of force but also for comparisons between muscles or studies when relative force levels might be different. Increases in voluntary descending drive amplify PICs in humans and provide an important spinal mechanism for motor unit discharging, and thus force output modulation.NEW & NOTEWORTHY Animal experiments and computational models have shown that motor neurons can amplify the synaptic input they receive via persistent inward currents. Here we show in humans that this amplification varies proportionally to the magnitude of the voluntary drive to the muscle.

Increasing Muscle Speed Drives Changes in the Neuromuscular Transform of Motor Commands during Postnatal Development in Songbirds

Progressive changes in vocal behavior over the course of vocal imitation leaning are often attributed exclusively to developing neural circuits, but the effects of postnatal body changes remain unknown. In songbirds, the syrinx transforms song system motor commands into sound and exhibits changes during song learning. Here we test the hypothesis that the transformation from motor commands to force trajectories by syringeal muscles functionally changes over vocal development in zebra finches. Our data collected in both sexes show that, only in males, muscle speed significantly increases and that supralinear summation occurs and increases with muscle contraction speed. Furthermore, we show that previously reported submillisecond spike timing in the avian cortex can be resolved by superfast syringeal muscles and that the sensitivity to spike timing increases with speed. Because motor neuron and muscle properties are tightly linked, we make predictions on the boundaries of the yet unknown motor code that correspond well with cortical activity. Together, we show that syringeal muscles undergo essential transformations during song learning that drastically change how neural commands are translated into force profiles and thereby acoustic features. We propose that the song system motor code must compensate for these changes to achieve its acoustic targets. Our data thus support the hypothesis that the neuromuscular transformation changes over vocal development and emphasizes the need for an embodied view of song motor learning.SIGNIFICANCE STATEMENT Fine motor skill learning typically occurs in a postnatal period when the brain is learning to control a body that is changing dramatically due to growth and development. How the developing body influences motor code formation and vice versa remains largely unknown. Here we show that vocal muscles in songbirds undergo critical transformations during song learning that drastically change how neural commands are translated into force profiles and thereby acoustic features. We propose that the motor code must compensate for these changes to achieve its acoustic targets. Our data thus support the hypothesis that the neuromuscular transformation changes over vocal development and emphasizes the need for an embodied view of song motor learning.

Spinal cord injury and diaphragm neuromotor control

Introduction: Neuromotor control of diaphragm muscle and the recovery of diaphragm activity following spinal cord injury have been narrowly focused on ventilation. By contrast, the understanding of neuromotor control for non-ventilatory expulsive/straining maneuvers (including coughing, defecation, and parturition) is relatively impoverished. This variety of behaviors are achieved via the recruitment of the diverse array of motor units that comprise the diaphragm muscle.Areas covered: The neuromotor control of ventilatory and non-ventilatory behaviors in health and in the context of spinal cord injury is explored. Particular attention is played to the neuroplasticity of phrenic motor neurons in various models of cervical spinal cord injury.Expert opinion: There is a remarkable paucity in our understanding of neuromotor control of maneuvers in spinal cord injury patients. Dysfunction of these expulsive/straining maneuvers reduces patient quality of life and contributes to severe morbidity and mortality. As spinal cord injury patient life expectancies continue to climb steadily, a nexus of spinal cord injury and age-associated comorbidities are likely to occur. While current research remains concerned only with the minutiae of ventilation, the major functional deficits of this clinical cohort will persist intractably. We posit some future research directions to avoid this scenario.

Linking Motoneuron PIC Location to Motor Function in Closed-Loop Motor Unit System Including Afferent Feedback: A Computational Investigation

The goal of this study is to investigate how the activation location of persistent inward current (PIC) over motoneuron dendrites is linked to motor output in the closed-loop motor unit. Here, a physiologically realistic model of a motor unit including afferent inputs from muscle spindles was comprehensively analyzed under intracellular stimulation at the soma and synaptic inputs over the dendrites during isometric contractions over a full physiological range of muscle lengths. The motor output of the motor unit model was operationally assessed by evaluating the rate of force development, the degree of force potentiation and the capability of self-sustaining force production. Simulations of the model motor unit demonstrated a tendency for a faster rate of force development, a greater degree of force potentiation, and greater capacity for self-sustaining force production under both somatic and dendritic stimulation of the motoneuron as the PIC channels were positioned farther from the soma along the path of motoneuron dendrites. Interestingly, these effects of PIC activation location on force generation significantly differed among different states of muscle length. The rate of force development and the degree of force potentiation were systematically modulated by the variation of PIC channel location for shorter-than-optimal muscles but not for optimal and longer-than-optimal muscles. Similarly, the warm-up behavior of the motor unit depended on the interplay between PIC channel location and muscle length variation. These results suggest that the location of PIC activation over motoneuron dendrites may be distinctively reflected in the motor performance during shortening muscle contractions.

MOTOR NEURON LOSS IN AGING AND AMYOTROPHIC LATERAL SCLEROSIS: DIFFERENT FUSE LENGTHS, SAME EXPLOSION

Advanced age and amyotrophic lateral sclerosis (ALS) are both associated with a loss of motor neurons resulting in muscle fiber atrophy and muscle weakness. Aging associated muscle fiber atrophy and weakening is termed sarcopenia, but the association with motor neuron loss is not as clearly established as in ALS, probably related to the prolonged time course of aging-related changes. Although aging and ALS effects on limb muscle strength and neuromotor performance are serious, such effects on the diaphragm muscle can be life threatening. Converging evidence indicates that larger phrenic motor neurons, innervating more fatigable type IIx and/or IIb diaphragm muscle fibers (fast fatigue intermediate, FInt and fast fatigable, FF motor units) are more susceptible to degeneration with both aging and ALS compared to smaller phrenic motor neurons innervating type I and IIa diaphragm muscle fibers (slow and fast fatigue resistant motor units, respectively). The etiology of ALS and age-related loss of motor neurons appears to involve mitochondrial function and neuroinflammation, both chronic and acute exacerbation. How mitochondrial dysfunction, neuroinflammation and motor neuron size intersect is the focus of continuing investigation.

Size-Dependent Vulnerability of Lumbar Motor Neuron Dendritic Degeneration in SOD1G93A Mice

The motor neuron (MN) soma surface area is correlated with motor unit type. Larger MNs innervate fast fatigue-intermediate (FInt) or fast-fatiguable (FF) muscle fibers in type FInt and FF motor units, respectively. Smaller MNs innervate slow-twitch fatigue-resistant (S) or fast fatigue-resistant (FR) muscle fibers in type S and FR motor units, respectively. In amyotrophic lateral sclerosis (ALS), FInt and FF motor units are more vulnerable, with denervation and MN death occurring for these units before the more resilient S and FR units. Abnormal MN dendritic arbors have been observed in ALS in humans and rodent models. We used a Golgi-Cox impregnation protocol to examine soma size-dependent changes in the dendritic morphology of lumbar MNs in SOD1^G93A mice, a model of ALS, at pre-symptomatic, onset and mid-disease stages. In wildtype control mice, the relationship between MN soma surface area and dendritic length or dendritic spine number was highly linear (i.e., increased MN soma size correlated with increased dendritic length and spines). By contrast, in SOD1^G93A mice, this linear relationship was lost and dendritic length reduction and spine loss were observed in larger MNs, from pre-symptomatic stages onward. These changes correlated with the neuromotor symptoms of ALS in rodent models. At presymptomatic ages, changes were restricted to the larger MNs, likely to comprise vulnerable FInt and FF motor units. Our results suggest morphological changes of MN dendrites and dendritic spines are likely to contribute ALS pathogenesis, not compensate for it. Anat Rec, 303:1455-1471, 2020. © 2019 American Association for Anatomy.