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 <i>p</i> < 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.