A NEW POSTURAL MOTOR RESPONSE TO SPINAL CORD STIMULATION: POST-STIMULATION REBOUND EXTENSION.

Spinal cord stimulation (SCS) has emerged as a therapeutic tool for improving motor function following spinal cord injury. While many studies focus on restoring locomotion, little attention is paid to enabling standing which is a prerequisite of walking. In this study, we fully characterize a new type of response to SCS, a long extension activated post-stimulation (LEAP). LEAP is primarily directed to ankle extensors and hence has great clinical potential to assist postural movements. To characterize this new response, we used the decerebrate cat model to avoid the suppressive effects of anesthesia, and combined EMG and force measurement in the hindlimb with intracellular recordings in the lumbar spinal cord. Stimulation was delivered as five-second trains via bipolar electrodes placed on the cord surface, and multiple combinations of stimulation locations (L4 to S2), amplitudes (50-600 uA), and frequencies (10-40 Hz) were tested. While the optimum stimulation location and frequency differed slightly among animals, the stimulation amplitude was key for controlling LEAP duration and amplitude. To study the mechanism of LEAP, we performed in vivo intracellular recordings of motoneurons. In 70% of motoneurons, LEAP increased at hyperpolarized membrane potentials indicating a synaptic origin. Furthermore, spinal interneurons exhibited changes in firing during LEAP, confirming the circuit origin of this behavior. Finally, to identify the type of afferents involved in generating LEAP, we used shorter stimulation pulses (more selective for proprioceptive afferents), as well as peripheral stimulation of the sural nerve (cutaneous afferents). The data indicates that LEAP primarily relies on proprioceptive afferents and has major differences from pain or withdrawal reflexes mediated by cutaneous afferents. Our study has thus identified and characterized a novel postural motor response to SCS which has the potential to expand the applications of SCS for patients with motor disorders.

Voluntary co-contraction of ankle muscles alters motor unit discharge characteristics and reduces estimates of persistent inward currents.

Motoneuronal persistent inward currents (PICs) are both facilitated by neuromodulatory inputs and highly sensitive to local inhibitory circuits (e.g., Ia reciprocal inhibition). Methods aimed to increase group Ia reciprocal inhibition from the antagonistic muscle have been successful in decreasing PICs, and the diffuse actions of neuromodulators released during activation of remote muscles have increased PICs. However, it remains unknown how motoneurons function in the presence of simultaneous excitatory and inhibitory commands. To probe this topic, we investigated motor unit (MU) discharge patterns and estimated PICs during voluntary co-contraction of ankle muscles, which simultaneously demands the contraction of agonist-antagonist pairs. Twenty young adults randomly performed triangular ramps (10s up and down) of both co-contraction (simultaneous dorsiflexion and plantarflexion) and isometric dorsiflexion to a peak of 30% of their maximum muscle activity from a maximal voluntary contraction. Motor unit spike trains were decomposed from high-density surface electromyography recorded over the tibialis anterior (TA) using blind source separation algorithms. Voluntary co-contraction altered motor unit discharge rate characteristics, decreasing estimates of PICs by 20% (4.47 pulses per second (pps) vs 5.57 pps during isometric dorsiflexion). These findings suggest that, during voluntary co-contraction, the inhibitory input from the antagonist muscle overcomes the additional excitatory and neuromodulatory drive that may occur due to the co-contraction of the antagonist muscle, which constrains PIC behavior.

Muscle Spasms after Spinal Cord Injury Stem from Changes in Motoneuron Excitability and Synaptic Inhibition, Not Synaptic Excitation.

Muscle spasms are common in chronic spinal cord injury (SCI), posing challenges to rehabilitation and daily activities. Pharmacological management of spasms mostly targets suppression of excitatory inputs, an approach known to hinder motor recovery. To identify better targets, we investigated changes in inhibitory and excitatory synaptic inputs to motoneurons as well as motoneuron excitability in chronic SCI. We induced either a complete or incomplete SCI in adult mice of either sex and divided those with incomplete injury into low or high functional recovery groups. Their sacrocaudal spinal cords were then extracted and used to study plasticity below injury, with tissue from naive animals as a control. Electrical stimulation of the dorsal roots elicited spasm-like activity in preparations of chronic severe SCI but not in the control. To evaluate overall synaptic inhibition activated by sensory stimulation, we measured the rate-dependent depression of spinal root reflexes. We found inhibitory inputs to be impaired in chronic injury models. When synaptic inhibition was blocked pharmacologically, all preparations became clearly spastic, even the control. However, preparations with chronic injuries generated longer spasms than control. We then measured excitatory postsynaptic currents (EPSCs) in motoneurons during sensory-evoked spasms. The data showed no difference in the amplitude of EPSCs or their conductance among animal groups. Nonetheless, we found that motoneuron persistent inward currents activated by the EPSCs were increased in chronic SCI. These findings suggest that changes in motoneuron excitability and synaptic inhibition, rather than excitation, contribute to spasms and are better suited for more effective therapeutic interventions.Significance Statement Neural plasticity following spinal cord injury is crucial for recovery of motor function. Unfortunately, this process is blemished by maladaptive changes that can cause muscle spasms. Pharmacological alleviation of spasms without compromising the recovery of motor function has proven to be challenging. Here, we investigated changes in fundamental spinal mechanisms that can cause spasms post-injury. Our data suggest that the current management strategy for spasms is misdirected toward suppressing excitatory inputs, a mechanism that we found unaltered after injury, which can lead to further motor weakness. Instead, this study shows that more promising approaches might involve restoring synaptic inhibition or modulating motoneuron excitability.

Motor unit firing rate modulation is more impaired during flexion synergy-driven contractions of the biceps brachii in chronic stroke.

Following a hemiparetic stroke, individuals exhibit altered motor unit firing patterns during voluntary muscle contractions, including impairments in firing rate modulation and recruitment. These individuals also exhibit abnormal muscle coactivation through multi-joint synergies (e.g., flexion synergy). Here, we investigate whether motor unit firing activity during flexion synergy-driven contractions of the paretic biceps brachii differs from that of voluntary contractions and use these differences to predict changes in descending motor commands. To accomplish this, we characterized motor unit firing patterns of the biceps brachii in individuals with chronic hemiparetic stroke during voluntary isometric elbow flexion contractions in the paretic and non-paretic limbs, as well as during contractions driven by voluntary effort and by flexion synergy expression in the paretic limb. We observed significant reductions in motor unit firing rate modulation from the non-paretic to paretic limb (non-paretic - paretic: 0.14 pps/%MVT, 95% CI: [0.09 0.19]) that were further reduced during synergy-driven contractions (voluntary paretic - synergy driven: 0.19 pps/%MVT, 95% CI: [0.14 0.25]). Moreover, using recently developed metrics, we evaluated how a stroke-induced reliance on indirect motor pathways alters the inputs that motor units receive and revealed progressive increases in neuromodulatory and inhibitory drive to the motor pool in the paretic limb, with the changes greatest during synergy-driven contractions. These findings suggest that an interplay between heightened neuromodulatory drive and alterations in inhibitory command structure may account for the observed motor unit impairments, further illuminating underlying neural mechanisms involved in the flexion synergy and its impact on motor unit firing patterns post-stroke.

Estimates of persistent inward currents in lower limb motoneurons are larger in females than in males.

Noninvasive recordings of motor unit (MU) spike trains help us understand how the nervous system controls movement and how it adapts to various physiological conditions. The majority of participants in human and nonhuman animal physiology studies are male, and it is assumed that mechanisms uncovered in these studies are shared between males and females. However, sex differences in neurological impairment and physical performance warrant the study of sex as a biological variable in human physiology and performance. To begin addressing this gap in the study of biophysical properties of human motoneurons, we quantified MU discharge rates and estimates of persistent inward current (PIC) magnitude in both sexes. We decomposed MU spike trains from the tibialis anterior (TA), medial gastrocnemius (MG), and soleus (SOL) using high-density surface electromyography and blind source separation algorithms. Ten participants of each sex performed slow triangular (10 s up and down) isometric contractions to a peak of 30% of their maximum voluntary contraction. We then used linear mixed-effects models to determine if peak discharge rate and estimates of PICs were predicted by the fixed effects of sex, muscle, and their interaction. Despite a lack of sex-differences in peak discharge rates across all muscles, estimates of PICs were larger [χ2(1) = 6.26, P = 0.012] in females [4.73 ± 0.242 pulses per second (pps)] than in males (3.81 ± 0.240 pps). These findings suggest that neuromodulatory drive, inhibitory input, and/or biophysical properties of motoneurons differ between the sexes and may contribute to differences in MU discharge patterns.NEW & NOTEWORTHY Sex-related differences in motoneuron analyses have emerged with greater inclusion of female participants, however, mechanisms for these differences remain unclear. Estimates of persistent inward currents (i.e., ΔF) in motoneurons of the lower limb muscles were larger in females than in males. This suggests neuromodulatory drive, monoaminergic signaling, intrinsic motoneuron properties, and/or descending motor commands may differ between the sexes, which provides a potential mechanism underlying previously reported sex-related differences in motoneuron discharge patterns.

A geometric approach to quantifying the neuromodulatory effects of persistent inward currents on individual motor unit discharge patterns.

Objective.All motor commands flow through motoneurons, which entrain control of their innervated muscle fibers, forming a motor unit (MU). Owing to the high fidelity of action potentials within MUs, their discharge profiles detail the organization of ionotropic excitatory/inhibitory as well as metabotropic neuromodulatory commands to motoneurons. Neuromodulatory inputs (e.g. norepinephrine, serotonin) enhance motoneuron excitability and facilitate persistent inward currents (PICs). PICs introduce quantifiable properties in MU discharge profiles by augmenting depolarizing currents upon activation (i.e. PIC amplification) and facilitating discharge at lower levels of excitatory input than required for recruitment (i.e. PIC prolongation).Approach. Here, we introduce a novel geometric approach to estimate neuromodulatory and inhibitory contributions to MU discharge by exploiting discharge non-linearities introduced by PIC amplification during time-varying linear tasks. In specific, we quantify the deviation from linear discharge ('brace height') and the rate of change in discharge (i.e. acceleration slope, attenuation slope, angle). We further characterize these metrics on a simulated motoneuron pool with known excitatory, inhibitory, and neuromodulatory inputs and on human MUs (number of MUs; Tibialis Anterior: 1448, Medial Gastrocnemius: 2100, Soleus: 1062, First Dorsal Interosseus: 2296).Main results. In the simulated motor pool, we found brace height and attenuation slope to consistently indicate changes in neuromodulation and the pattern of inhibition (excitation-inhibition coupling), respectively, whereas the paired MU analysis (ΔF) was dependent on both neuromodulation and inhibition pattern. Furthermore, we provide estimates of these metrics in human MUs and show comparable variability in ΔFand brace height measures for MUs matched across multiple trials.Significance. Spanning both datasets, we found brace height quantification to provide an intuitive method for achieving graded estimates of neuromodulatory and inhibitory drive to individual MUs. This complements common techniques and provides an avenue for decoupling changes in the level of neuromodulatory and pattern of inhibitory motor commands.

GABA Increases Sensory Transmission In Monkeys.

Sensory input flow is central to voluntary movements. For almost a century, GABA was believed to modulate this flow by inhibiting sensory axons in the spinal cord to sculpt neural inputs into skilled motor output. Instead, here we show that GABA can also facilitate sensory transmission in monkeys and consequently increase spinal and cortical neural responses to sensory inputs challenging our understanding of generation and perception of movement.

Biomechanical and Sensory Feedback Regularize the Behavior of Different Locomotor Central Pattern Generators.

This work presents an in-depth numerical investigation into a hypothesized two-layer central pattern generator (CPG) that controls mammalian walking and how different parameter choices might affect the stepping of a simulated neuromechanical model. Particular attention is paid to the functional role of features that have not received a great deal of attention in previous work: the weak cross-excitatory connectivity within the rhythm generator and the synapse strength between the two layers. Sensitivity evaluations of deafferented CPG models and the combined neuromechanical model are performed. Locomotion frequency is increased in two different ways for both models to investigate whether the model's stability can be predicted by trends in the CPG's phase response curves (PRCs). Our results show that the weak cross-excitatory connection can make the CPG more sensitive to perturbations and that increasing the synaptic strength between the two layers results in a trade-off between forced phase locking and the amount of phase delay that can exist between the two layers. Additionally, although the models exhibit these differences in behavior when disconnected from the biomechanical model, these differences seem to disappear with the full neuromechanical model and result in similar behavior despite a variety of parameter combinations. This indicates that the neural variables do not have to be fixed precisely for stable walking; the biomechanical entrainment and sensory feedback may cancel out the strengths of excitatory connectivity in the neural circuit and play a critical role in shaping locomotor behavior. Our results support the importance of including biomechanical models in the development of computational neuroscience models that control mammalian locomotion.

The Cellular Basis for the Generation of Firing Patterns in Human Motor Units.

Motor units, which comprise a motoneuron and the set of muscle fibers it innervates, are the fundamental neuromuscular transducers for all motor commands. The one to one relationship between a motoneuron and its innervated muscle fibers allow motoneuron firing patterns to be readily measured in humans. In this chapter, we summarize the current understanding of the cellular basis for the generation of firing patterns in human motor units. We provide a brief review of landmark insights from classic studies and then proceed to consider the features of motor unit firing patterns that are most likely to be sensitive estimators of motoneuron inputs and properties. In addition, we discuss recent advances in technology for recording human motor unit firing patterns and highly realistic computer simulations of motoneurons. The final section presents our recent efforts to use the power of supercomputers for implementation of the motoneuron models, with a goal of achieving a true "reverse engineering" approach that maximizes the insights from motor unit firing patterns into the synaptic structure of motor commands.

Motor Unit Discharge Patterns in Response to Focal Tendon Vibration of the Lower Limb in Cats and Humans.

High-frequency vibration of the tendon provides potent activation of Ia afferents time-locked to the stimulation frequency and provides excitatory ionotropic activation of homonymous motor pools. In cats, the evoked motor unit discharge is constrained to discharge at integer multiples of the vibration frequency, resulting in a probability of discharge that is highly punctuated. Here we quantify the robustness of this punctuated response in the cat and evaluate whether it is present in the human. Soleus electromyography (EMG) was collected from eight cats using 64 channel electrodes during three modes of motoneuron activation. First, tendon vibration parameters were modified. Second, secondary reflex inputs are applied concurrently with tendon vibration. Third, the state of the spinal cord was altered through pharmacological or surgical manipulations. Analogous surface high-density EMG was collected from the lower leg of six humans during both vibration evoked and matched volitional contractions. Array EMG signals from both the cat and human were decomposed into corresponding motor unit action potential spike trains, and the punctuation in discharge was quantified. In the cat, regardless of vibration parameters, secondary synaptic drive, and state of spinal circuitry, focal tendon vibration evoked punctuated motor unit discharge. However, in the human lower limb, the vibration-evoked contractions do not produce punctuated motor unit discharge.

A single-arm pilot of a web-based intervention to improve itch-related quality of life.

BACKGROUND: Many patients with chronic itch and atopic dermatitis (AD) or psoriasis do not receive/use available medical and psychosocial treatments properly due to system, provider and/or patient factors. OBJECTIVE: An educational website (ITCH-RELIEF) to improve itch-related quality of life (QoL) for adults with AD or psoriasis and chronic itch was developed and assessed. ITCH RELIEF stands for Interactive Toolbox of Comprehensive Health Resources to Enhance Living with Itch - Educational Facilitation (for Adults). METHODS: Single-arm pre- and post-test design with 1-month follow-up (N = 137 at baseline). RESULTS: There was statistically and clinically significant improvement in the primary outcome of itch-related QoL impairment as assessed by the ItchyQoL from baseline [M = 78.9, 95% confidence interval (CI) = 75.9, 81.9] to follow up (M = 75.4, CI = 72.4, 78.5), P = 0.007, as well as statistically significant improvement in several itch-related secondary outcomes (all Ps < 0.05). CONCLUSIONS: This study demonstrated initial effectiveness of an online intervention to improve itch-related QoL among individuals with AD or psoriasis and chronic itch. Future studies should address limitations by randomizing more heterogeneous participants, utilizing a longer follow-up and assessing medication use.

A computational approach for generating continuous estimates of motor unit discharge rates and visualizing population discharge characteristics.

Objective. Successive improvements in high density surface electromyography and decomposition techniques have facilitated an increasing yield in decomposed motor unit (MU) spike times. Though these advancements enhance the generalizability of findings and promote the application of MU discharge characteristics to inform the neural control of motor output, limitations remain. Specifically, (1) common approaches for generating smooth estimates of MU discharge rates introduce artifacts in quantification, which may bias findings, and (2) discharge characteristics of large MU populations are often difficult to visualize.Approach. In the present study, we propose support vector regression (SVR) as an improved approach for generating smooth continuous estimates of discharge rate and compare the fit characteristics of SVR to traditionally used methods, including Hanning window filtering and polynomial regression. Furthermore, we introduce ensembles as a method to visualize the discharge characteristics of large MU populations. We define ensembles as the average discharge profile of a subpopulation of MUs, composed of a time normalized ensemble average of all units within this subpopulation. Analysis was conducted with MUs decomposed from the tibialis anterior (N= 2128), medial gastrocnemius (N= 2673), and soleus (N= 1190) during isometric plantarflexion and dorsiflexion contractions.Main result. Compared to traditional approaches, we found SVR to alleviate commonly observed inaccuracies and produce significantly less absolute fit error in the initial phase of MU discharge and throughout the entire duration of discharge. Additionally, we found the visualization of MU populations as ensembles to intuitively represent population discharge characteristics with appropriate accuracy for visualization.Significance. The results and methods outlined here provide an improved method for generating estimates of MU discharge rate with SVR and present a unique approach to visualizing MU populations with ensembles. In combination, the use of SVR and generation of ensembles represent an efficient method for rendering population discharge characteristics.

Estimates of persistent inward currents are reduced in upper limb motor units of older adults.

Ageing is a natural process causing alterations in the neuromuscular system, which contributes to reduced quality of life. Motor unit (MU) contributes to weakness, but the mechanisms underlying reduced firing rates are unclear. Persistent inward currents (PICs) are crucial for initiation, gain control and maintenance of motoneuron firing, and are directly proportional to the level of monoaminergic input. Since concentrations of monoamines (i.e. serotonin and noradrenaline) are reduced with age, we sought to determine if estimates of PICs are reduced in older (>60 years old) compared to younger adults (<35 years old). We decomposed MU spike trains from high-density surface electromyography over the biceps and triceps brachii during isometric ramp contractions to 20% of maximum. Estimates of PICs (ΔFrequency; or simply ΔF) were computed using the paired MU analysis technique. Regardless of the muscle, peak firing rates of older adults were reduced by ∼1.6 pulses per second (pps) (P = 0.0292), and ΔF was reduced by ∼1.9 pps (P < 0.0001), compared to younger adults. We further found that age predicted ΔF in older adults (P = 0.0261), resulting in a reduction of ∼1 pps per decade, but there was no relationship in younger adults (P = 0.9637). These findings suggest that PICs are reduced in the upper limbs of older adults during submaximal isometric contractions. Reduced PIC magnitude represents one plausible mechanism for reduced firing rates and function in older individuals, but further work is required to understand the implications in other muscles and during a variety of motor tasks. KEY POINTS: Persistent inward currents play an important role in the neural control of human movement and are influenced by neuromodulation via monoamines originating in the brainstem. During ageing, motor unit firing rates are reduced, and there is deterioration of brainstem nuclei, which may reduce persistent inward currents in alpha motoneurons. Here we show that estimates of persistent inward currents (ΔF) of both elbow flexor and extensor motor units are reduced in older adults. Estimates of persistent inward currents have a negative relationship with age in the older adults, but not in the young. This novel mechanism may play a role in the alteration of motor firing rates that occurs with ageing, which may have consequences for motor control.

Estimates of persistent inward currents in tibialis anterior motor units during standing ramped contraction tasks in humans.

Persistent inward currents (PICs) play an essential role in setting motor neuron gain and shaping motor unit firing patterns. Estimates of PICs in humans can be made using the paired motor unit analysis technique, which quantifies the difference in discharge rate of a lower threshold motor unit at the recruitment onset and offset of a higher threshold motor unit (ΔF). Because PICs are highly dependent on the level of neuromodulatory drive, ΔF represents an estimate of level of neuromodulation at the level of the spinal cord. Most of the estimates of ΔF are performed under constrained, isometric, seated conditions. In the present study, we used high-density surface EMG arrays to discriminate motor unit firing patterns during isometric seated conditions with torque or EMG visual feedback and during unconstrained standing anterior-to-posterior movements with root mean square EMG visual feedback. We were able to apply the paired motor unit analysis technique to the decomposed motor units in each of the three conditions. We hypothesized that ΔF would be higher during unconstrained standing anterior-to-posterior movements compared with the seated conditions, reflecting an increase in the synaptic input to motoneurons drive while standing. In agreement with previous work, we found that there was no evidence of a difference in ΔF between the seated and standing postures, although slight differences in the initial and peak discharge rates were observed. Taken together, our results suggest that both the standing and seated postures are likely not sufficiently different, both being "upright" postures, to result in large changes in neuromodulatory drive.NEW & NOTEWORTHY In the present study, we show that the discharge rate of a lower threshold motor unit at the recruitment onset and offset of a higher threshold motor unit (ΔF) is similar between standing and seated conditions in human tibialis anterior motor units, suggesting that at least for these two upright postures neuromodulatory drive is similar. We also highlight a proposed technological development in using high-density EMG arrays for real-time muscle activity feedback to accomplish standing ramped contraction tasks and demonstrate the validity of the paired motor unit analysis technique during these conditions.

Slowly activating outward membrane currents generate input-output sub-harmonic cross frequency coupling in neurons.

A major challenge in understanding spike-time dependent information encoding in the neural system is the non-linear firing response to inputs of the individual neurons. Hence, quantitative exploration of the putative mechanisms of this non-linear behavior is fundamental to formulating the theory of information transfer in the neural system. The objective of this simulation study was to evaluate and quantify the effect of slowly activating outward membrane current, on the non-linearity in the output of a one-compartment Hodgkin-Huxley styled neuron. To evaluate this effect, the peak conductance of the slow potassium channel (gK-slow) was varied from 0% to 200% of its normal value in steps of 33%. Both cross- and iso-frequency coupling between the input and the output of the simulated neuron was computed using a generalized coherence measure, i.e., n:m coherence. With increasing gK-slow, the amount of sub-harmonic cross-frequency coupling, where the output frequencies (1-8 Hz) are lower than the input frequencies (15-35 Hz), increased progressively whereas no change in iso-frequency coupling was observed. Power spectral and phase-space analysis of the neuronal membrane voltage vs. slow potassium channel activation variable showed that the interaction of the slow channel dynamics with the fast membrane voltage dynamics generates the observed sub-harmonic coupling. This study provides quantitative insights into the role of an important membrane mechanism i.e. the slowly activating outward current in generating non-linearities in the output of a neuron.

Inability to increase the neural drive to muscle is associated with task failure during submaximal contractions.

We investigated changes in motor unit (MU) behavior and vasti-muscle contractile properties during sustained submaximal fatiguing contractions with a new time-domain tracking technique to understand the mechanisms responsible for task failure. Sixteen participants performed a nonfatiguing 15-s isometric knee extension at 50% of the maximum voluntary (MVC) torque, followed by a 30% MVC sustained contraction until exhaustion. Two grids of 64 surface electromyography electrodes were placed over vastus medialis and lateralis. Signals were decomposed into MU discharge times and the MUs from the 30% MVC sustained contraction were followed until task failure by overlapping decomposition intervals. These MUs were then tracked between 50% and 30% MVC. During the sustained fatiguing contraction, MUs of the two muscles decreased their discharge rate until ∼40% of the endurance time, referred to as the reversal time, and then increased their discharge rate until task failure. This reversal in firing behavior predicted total endurance time and was matched by opposite changes in twitch force (increase followed by a decrease). Despite the later increase in MU firing rates, peak discharge rates at task failure did not reach the frequency attained during a nonfatiguing 50% MVC contraction. These results show that changes in MU firing properties are influenced by adjustments in contractile properties during the course of the contraction, allowing the identification of two phases. Nevertheless, the contraction cannot be sustained, possibly because of progressive motoneuron inhibition/decreased excitability, as the later increase in firing rate saturates at a much lower frequency compared with a higher-force nonfatiguing contraction.NEW & NOTEWORTHY Motor unit firing and contractile properties during a submaximal contraction until failure were assessed with a new tracking technique. Two distinct phases in firing behavior were observed, which compensated for changes in twitch area and predicted time to failure. However, the late increase in firing rate was below the rates attained in absence of fatigue, which points to an inability of the central nervous system to sufficiently increase the neural drive to muscle with fatigue.

Motor Unit Discharge Variability Is Increased in Mild-To-Moderate Parkinson's Disease.

Individuals with Parkinson's disease (PD) demonstrate deficits in muscle activation such as decreased amplitude and inappropriate bursting. There is evidence that some of these disturbances are more pronounced in extensor vs. flexor muscles. Surface EMG has been used widely to quantify muscle activation deficits in PD, but analysis of discharge of the underlying motor units may provide greater insight and be more sensitive to changes early in the disease. Of the few studies that have examined motor unit discharge in PD, the majority were conducted in the first dorsal interosseous, and no studies have measured motor units from extensor and flexor muscles within the same cohort. The objective of this study was to characterize the firing behavior of single motor units in the elbow flexor and extensor muscles during isometric contractions in people with mild-to-moderate PD. Ten individuals with PD (off-medication) and nine healthy controls were tested. Motor unit spike times were recorded via intramuscular EMG from the biceps and triceps brachii muscles during 30-s isometric contractions at 10% maximum voluntary elbow flexion and elbow extension torque, respectively. We selected variables of mean motor unit discharge rate, discharge variability, and torque variability to evaluate motor abnormalities in the PD group. The effects of group, muscle, and group-by-muscle on each variable were determined using separate linear mixed models. Discharge rate and torque variability were not different between groups, but discharge variability was significantly higher in the PD group for both muscles combined (p < 0.0001). We also evaluated the asymmetry in these motor variables between the triceps and biceps for each individual participant with PD to evaluate whether there was an association with disease severity. The difference in torque variability between elbow flexion and extension was significantly correlated with both the Hoehn and Yahr scale (rho = 0.71) and UPDRS (rho = 0.62). Our findings demonstrate that variability in motor output, rather than decreased discharge rates, may contribute to motor dysfunction in people with mild-to-moderate PD. Our findings provide insight into altered neural control of movement in PD and demonstrate the importance of measuring from multiple muscles within the same cohort.