The variability of reflex amplitude estimates in motor unit pools depends on the phenotype distribution and discharge statistics

Motor unit (MU) activity during electrically or mechanically evoked reflexes is used to examine the relationship between neural inputs and MU properties. However, variations in single-MU reflex amplitudes are not fully understood and limit their reliability in determining the input-output relation of motor neurons (MNs). Using experiments and computer simulations, we investigated (i) whether MN discharge statistics and muscle activation explain the variability of reflex amplitude estimates and (ii) whether these variations are reflected differently across distinct reflex amplitude estimation methods. We analyzed MU spike trains extracted from isometric contractions of the tibialis anterior muscle at 10 % and 20 % MVC (maximum voluntary contraction). Estimating reflex amplitudes based on the peristimulus frequencygram (PSF) at 10 % MVC, the linear regression between discharge rate (DR) and reflex amplitude was always positive, with p < 0.05 in 3 out of 6 subjects; however, the linear correlation was inconsistent at 20 % MVC. We thereby observed that inter-subject variability was associated with the coefficient of variation of the interspike intervals. Furthermore, the linear correlation between DR and peristimulus time histogram (PSTH) based reflex amplitudes was inconsistent for both contraction forces. To obtain further insights into the influence of MN properties, we simulated reflexes in a heterogeneous MN population using electrical circuit models and varied MN inputs. The simulations indicate that, besides mean input current and membrane noise, MN properties also contribute to the variability of reflex amplitude estimates. The MN heterogeneity is well captured by PSF-based reflex estimates but not by PSTH-based ones. These results show that variations in amplitude estimates of individual MU reflexes are due to complex interactions between intrinsic and extrinsic factors. As PSF-based reflex amplitude estimates reflect the MN size distribution, tracking PSF-based reflex amplitudes at fixed MVC levels across individual subjects could serve as a marker for investigating spinal adaptations under (patho)physiological conditions. Author summaryMotor neurons are specialized nerve cells that control human movement. Each motor neuron activates a specific set of muscle fibers, and the functional unit consisting of a motor neuron and muscle fibers is called a motor unit. We can observe the activity of motor neurons in humans by decomposing the electrical activity of muscles (the electromyogram) into contributions from individual motor units. Reflex responses of motor units are often used to study the input-output relation of motor neurons in humans. We used a combination of experiments and computer simulations to study the factors that influence the reflex amplitude of motor units during an excitatory reflex. We found that the reflex amplitude is non-linearly influenced by a number of intrinsic and extrinsic factors, e. g., motor neuron size, but also the muscle force. Additionally, we found that these factors have different effects on the results of the two common methods used to calculate the reflex amplitude. These results provide guidance on choosing a suitable evaluation method and on interpreting reflex experiments.