EMG analysis tuned for determining the timing and level of activation in different motor units

Surface electromyographic frequency characteristics of the quadriceps differ between continuous high- and low-torque isometric knee extension to momentary failure

Surface EMG (sEMG) has been used to compare loading conditions during exercise. Studies often explore mean/median frequencies. This potentially misses more nuanced electrophysiological differences between exercise tasks. Therefore, wavelet-based analysis was used to evaluate electrophysiological characteristics in the sEMG signal of the quadriceps under both higher- and lower-torque (70 % and 30 % of MVC, respectively) isometric knee extension performed to momentary failure. Ten recreationally active adult males with previous resistance training experience were recruited. Using a within-session, repeated-measures, randomised crossover design, participants performed isometric knee extension whilst sEMG was collected from the vastus medialis (VM), rectus femoris (RF) and vastus lateralis (VL). Mean signal frequency showed similar characteristics in each condition at momentary failure. However, individual wavelets revealed different frequency component changes between the conditions. All frequency components increased during the low-torque condition. But low-frequency components increased, and high-frequency components decreased, in intensity throughout the high-torque condition. This resulted in convergence of the low-torque and high-torque trial wavelet characteristics towards the end of the low-torque trial. Our results demonstrate a convergence of myoelectric signal properties between low- and high-torque efforts with fatigue via divergent signal adaptations. Further work should disentangle factors influencing frequency characteristics during exercise tasks.

A review of the efforts to develop muscle and musculoskeletal models for biomechanics in the last 50 years

Both the Hill and the Huxley muscle models had already been described by the time the International Society of Biomechanics was founded 50 years ago, but had seen little use before the 1970s due to the lack of computing. As computers and computational methods became available in the 1970s, the field of musculoskeletal modeling developed and Hill type muscle models were adopted by biomechanists due to their relative computational simplicity as compared to Huxley type muscle models. Muscle forces computed by Hill type muscle models provide good agreement in conditions similar to the initial studies, i.e. for small muscles contracting under steady and controlled conditions. However, more recent validation studies have identified that Hill type muscle models are least accurate for natural in vivo locomotor behaviours at submaximal activations, fast speeds and for larger muscles, and thus need to be improved for their use in understanding human movements. Developments in muscle modelling have tackled these shortcomings. However, over the last 50 years musculoskeletal simulations have been largely based on traditional Hill type muscle models or even simplifications of this model that neglected the interaction of the muscle with a compliant tendon. The introduction of direct collocation in musculoskeletal simulations about 15 years ago along with further improvements in computational power and numerical methods enabled the use of more complex muscle models in simulations of whole-body movement. Whereas Hill type models are still the norm, we may finally be ready to adopt more complex muscle models into musculoskeletal simulations of human movement.

The impact of age-related increase in passive muscle stiffness on simulated upper limb reaching

Ageing changes the musculoskeletal and neural systems, potentially affecting a person's ability to perform daily living activities. One of these changes is increased passive stiffness of muscles, but its contribution to performance is difficult to separate experimentally from other ageing effects such as loss of muscle strength or cognitive function. A computational upper limb model was used to study the effects of increasing passive muscle stiffness on reaching performance across the model's workspace (all points reachable with a given model geometry). The simulations indicated that increased muscle stiffness alone caused deterioration of reaching accuracy, starting from the edges of the workspace. Re-tuning the model's control parameters to match the ageing muscle properties does not fully reverse ageing effects but can improve accuracy in selected regions of the workspace. The results suggest that age-related muscle stiffening, isolated from other ageing effects, impairs reaching performance. The model also exhibited oscillatory instability in a few simulations when the controller was tuned to the presence of passive muscle stiffness. This instability is not observed in humans, implying the presence of natural stabilizing strategies, thus pointing to the adaptive capacity of neural control systems as a potential area of future investigation in age-related muscle stiffening.

Interfacing Motor Units in Nonhuman Primates Identifies a Principal Neural Component for Force Control Constrained by the Size Principle

Motor units convert the last neural code of movement into muscle forces. The classic view of motor unit control is that the CNS sends common synaptic inputs to motoneuron pools and that motoneurons respond in an orderly fashion dictated by the size principle. This view, however, is in contrast with the large number of dimensions observed in motor cortex, which may allow individual and flexible control of motor units. Evidence for flexible control of motor units may be obtained by tracking motor units longitudinally during tasks with some level of behavioral variability. Here we identified and tracked populations of motor units in the brachioradialis muscle of two macaque monkeys during 10 sessions spanning >1 month with a broad range of rate of force development (1.8-38.6 N · m · s-1). We found a very stable recruitment order and discharge characteristics of the motor units over sessions and contraction trials. The small deviations from orderly recruitment were fully predicted by the motor unit recruitment intervals, so that small shifts in recruitment thresholds happened only during contractions at a high rate of force development. Moreover, we also found that one component explained more than ∼50% of the motor unit discharge rate variance, and that the remaining components represented a time-shifted version of the first. In conclusion, our results show that the recruitment of motoneurons is determined by the interplay of the size principle and common input and that this recruitment scheme is not violated over time or by the speed of the contractions.SIGNIFICANCE STATEMENT With a new noninvasive high-density electromyographic framework, we show the activity of motor unit ensembles in macaques during voluntary contractions. The discharge characteristics of brachioradialis motor units revealed a relatively fixed recruitment order and discharge characteristics across days and rate of force developments. These results were further confirmed through invasive axonal stimulation and recordings of intramuscular electromyographic activity from 16 arm muscles. The study shows for the first time the feasibility of longitudinal noninvasive motor unit interfacing and tracking of the same motor units in nonhuman primates.

The impact of intrinsic muscle properties on simulated reaching performance

Musculoskeletal modelling is used widely for studying limb motion and its control, but simulation outcomes may depend heavily on the underlying muscle model used. The aim of this study was to investigate how intrinsic muscle properties affect reaching movements in a simple upper limb model. The simulations suggest that more realistic, higher-order activation dynamics requires longer prediction from a forward model and gives rise to a higher level of unplanned co-contraction than simple activation models. Consistent with prior work, muscle force-length-velocity properties stabilised and smoothed limb movements and furthermore helped promote accurate reaching performance with the high-order activation model.

Modeling muscle function using experimentally determined subject-specific muscle properties

Muscle models are commonly based on intrinsic properties pooled across a number of individuals, often from a different species, and rarely validated against directly measured muscle forces. Here we use a rich data set of rat medial gastrocnemius muscle forces recorded during in-situ and in-vivo isometric, isotonic, and cyclic contractions to test the accuracy of forces predicted using Hill-type muscle models. We identified force-length and force-velocity parameters for each individual, and used either these subject-specific intrinsic properties, or population-averaged properties within the models. The modeled forces for cyclic in-vivo and in-situ contractions matched with measured muscle-tendon forces with r² between 0.70 and 0.86, and root-mean square errors (RMSE) of 0.10 to 0.13 (values normalized to the maximum isometric force). The modeled forces were least accurate at the highest movement and cycle frequencies and did not show an improvement in r² when subject-specific intrinsic properties were used; however, there was a reduction in the RMSE with fewer predictions having higher errors. We additionally recorded and tested muscle models specific to proximal and distal regions of the muscle and compared them to measures and models from the whole muscle belly: there was no improvement in model performance when using data from specific anatomical regions. These results show that Hill-type muscle models can yield very good performance for cyclic contractions typical of locomotion, with small reductions in errors when subject-specific intrinsic properties are used.

Effect of Sudden Deprivation of Sensory Inputs From Periodontium on Mastication

Objective

To investigate the effect of sudden deprivation of sensory inputs from the periodontium on jaw kinematics and time-varying activation profile of the masseter muscle.

Methods

Fourteen (age range: 22–26 years; four men) healthy and natural dentate volunteers participated in a single experimental session. During the experiment, the participants were asked to eat six hard visco-elastic test food models, three each before and after an anesthetic intervention. The movements of the jaw in three dimensions and electromyographic (EMG) activity of the masseter muscle on the chewing side were recorded.

Results

The results of the study showed no significant differences in the number of chewing cycles (P = 0.233) and the duration of chewing sequence (P = 0.198) due to sudden deprivation of sensory inputs from the periodontium. However, there was a significant increase in the jaw opening velocity (P = 0.030) and a significant increase in the duration of occlusal phase (P = 0.004) during the anesthetized condition. The EMG activity of the jaw closing phase was significantly higher during the control condition [116.5 arbitrary units (AU)] than anesthetized condition (93.9 AU). The temporal profile of the masseter muscle showed a biphasic increase in the excitatory muscle drive in the control condition but this increase was virtually absent during the anesthetized condition.

Conclusion

Sudden deprivation of sensory inputs from the periodontium affects the jaw kinematics and jaw muscle activity, with a clear difference in the time-varying activation profile of the masseter muscle. The activation profile of the masseter muscle shows that periodontal mechanoreceptors contribute to approximately 20% of the EMG activity during the jaw closing phase.

Investigation into pathophysiology of naturally occurring palatal instability and intermittent dorsal displacement of the soft palate (DDSP) in racehorses: Thyro-hyoid muscles fatigue during exercise

Exercise induced intermittent dorsal displacement of the soft palate (DDSP) is a common cause of airway obstruction and poor performance in racehorses. The definite etiology is still unclear, but through an experimental model, a role in the development of this condition was identified in the dysfunction of the thyro-hyoid muscles. The present study aimed to elucidate the nature of this dysfunction by investigating the spontaneous response to exercise of the thyro-hyoid muscles in racehorses with naturally occurring DDSP. Intramuscular electrodes were implanted in the thyro-hyoid muscles of nine racehorses, and connected to a telemetric unit for electromyographic monitoring implanted subcutaneously. The horses were recruited based on upper airway function evaluated through wireless endoscopy during exercise. Five horses, with normal function, were used as control; four horses were diagnosed as DDSP-affected horses based on repeated episodes of intermittent dorsal displacement of the soft palate. The electromyographic activity of the thyro-hyoid muscles recorded during incremental exercise tests on a high-speed treadmill was analyzed to measure the mean electrical activity and the median frequency of the power spectrum, thereafter subjected to wavelet decomposition. The affected horses had palatal instability with displacement on repeated exams prior to surgical implantation. Although palatal instability persisted after surgery, only two of these horses displaced the palate after instrumentation. The electromyographic traces from this group of four horses showed, at highest exercise intensity, a decrease in mean electrical activity and median power frequency, with progressive decrease in the contribution of the high frequency wavelets, consistent with development of thyro-hyoid muscle fatigue. The results of this study identified fatigue as the main factor leading to exercise induced palatal instability and DDSP in a group of racehorses. Further studies are required to evaluate the fiber type composition and metabolic characteristics of the thyro-hyoid muscles that could predispose to fatigue.

Adaptation of Fatigue Affected Changes in Muscle EMG Frequency Characteristics for the Determination of Training Load in Physical Therapy for Cancer Patients

Cancer patients often experience loss in body weight and also a decrease in muscle mass, which results in the reduction of physical activity and mobilization of the patient. To decelerate the loss of muscle mass, as part of the cancer treatment patients frequently undergo physical therapy and considering the physical capabilities of the patients, with moderate loads. Moreover, frequent studies also observed for cancer patients, together with the decrease in muscle mass a shift into fast-twitch muscle fibers from slow-twitch fibers. The aim of our study therefore was to determine how motor fibers behave under moderate isometric load executed until total exhaustion. 11 university students (G1), and 14 elite athletes (G2) participated in the study. 65% of the maximal voluntary contraction (MVC) was determined for the biceps brachii muscle, and with this load holding a weight, participants had to sustain a 90 deg. isometric elbow flexion in a standing posture until complete fatigue occurred. EMG activity for the biceps brachii muscle was measured and frequency analysis was performed. 3 windows were determined in the fatiguing protocol: the first (W1), middle (W2), and last (W3) 5 s, and also frequency analysis for MVC was performed (MAX) between 0 and 260 Hz with 20 Hz wide frequency bands. The results indicate, that as the protocol progressed in time and the effect of fatigue increased (from W1 to W3) the activity of low frequency muscle fibers significantly increased (0-40 Hz) while activity of high frequency muscle fibers (60-260 Hz) significantly decreased for G1 and G2 groups identically. We can conclude, that training applied with constant moderate tension as fatigue increases will result in the increased activation of the lower frequency slow-twitch muscle fibers, but the increase of fatigue in the lower frequency fibers will not result in the increase in the activation level of the higher frequency fast-twitch fibers. Consequently, because as slow-twitch fibers are being used at moderate loads and even when fatigue occurs in these fibers the fast-twitch fibers will not work, higher muscle loads are needed if the aim is to activate fast-twitch fibers. Considering the shift into fast-twitch muscle fibers from slow-twitch fibers for cancer patients, in some cases if the patient's age and physical status allows during the physical treatment, higher loads and consequently higher levels of activation might be beneficial for the retardment of loss concerning the fast-twitch fiber mass.

Assessing aberrant muscle activity patterns via the analysis of surface EMG data collected during a functional evaluation

BACKGROUND

Surface electromyographic (EMG) recordings collected during the performance of functional evaluations allow clinicians to assess aberrant patterns of muscle activity associated with musculoskeletal disorders. This assessment is typically achieved via visual inspection of the surface EMG data. This approach is time-consuming and leads to accurate results only when the assessment is carried out by an EMG expert.

METHODS

A set of algorithms was developed to automatically evaluate aberrant patterns of muscle activity. EMG recordings collected during the performance of functional evaluations in 62 subjects (22 to 61 years old) were used to develop and characterize the algorithms. Clinical scores were generated via visual inspection by an EMG expert using an ordinal scale capturing the severity of aberrant patterns of muscle activity. The algorithms were used in a case study (i.e. the evaluation of a subject with persistent back pain following instrumented lumbar fusion who underwent lumbar hardware removal) to assess the clinical suitability of the proposed technique.

RESULTS

The EMG-based algorithms produced accurate estimates of the clinical scores. Results were primarily obtained using a linear regression approach. However, when the results were not satisfactory, a regression implementation of a Random Forest was utilized, and the results compared with those obtained using a linear regression approach. The root-mean-square error of the clinical score estimates produced by the algorithms was a small fraction of the ordinal scale used to rate the severity of the aberrant patterns of muscle activity. Regression coefficients and associated 95% confidence intervals showed that the EMG-based estimates fit well the clinical scores generated by the EMG expert. When applied to the clinical case study, the algorithms appeared to capture the characteristics of the muscle activity patterns associated with persistent back pain following instrumented lumbar fusion.

CONCLUSIONS

The proposed approach relies on EMG-based measures to generate accurate estimates of the severity of aberrant patterns of muscle activity. The results obtained in the case study suggest that the proposed technique is suitable to derive clinically-relevant information from EMG data collected during functional evaluations.

Shifting gears: dynamic muscle shape changes and force-velocity behavior in the medial gastrocnemius

When muscles contract, they bulge in thickness or in width to maintain a (nearly) constant volume. These dynamic shape changes are tightly linked to the internal constraints placed on individual muscle fibers and play a key functional role in modulating the mechanical performance of skeletal muscle by increasing its range of operating velocities. Yet to date we have a limited understanding of the nature and functional implications of in vivo dynamic muscle shape change under submaximal conditions. This study determined how the in vivo changes in medial gastrocnemius (MG) fascicle velocity, pennation angle, muscle thickness, and subsequent muscle gearing varied as a function of force and velocity. To do this, we obtained recordings of MG tendon length, fascicle length, pennation angle, and thickness using B-mode ultrasound and muscle activation using surface electromyography during cycling at a range of cadences and loads. We found that that increases in contractile force were accompanied by reduced bulging in muscle thickness, reduced increases in pennation angle, and faster fascicle shortening. Although the force and velocity of a muscle contraction are inversely related due to the force-velocity effect, this study has shown how dynamic muscle shape changes are influenced by force and not influenced by velocity.NEW & NOTEWORTHY During movement, skeletal muscles contract and bulge in thickness or width. These shape changes play a key role in modulating the performance of skeletal muscle by increasing its range of operating velocities. Yet to date the underlying mechanisms associated with muscle shape change remain largely unexplored. This study identified muscle force, and not velocity, as the mechanistic driving factor to allow for muscle gearing to vary depending on the contractile conditions during human cycling.

Associations between motor unit action potential parameters and surface EMG features

The surface interference EMG signal provides some information on the neural drive to muscles. However, the association between neural drive to muscle and muscle activation has long been debated with controversial indications due to the unavailability of motor unit population data. In this study, we clarify the potential and limitations of interference EMG analysis to infer motor unit recruitment strategies with an experimental investigation of several concurrently active motor units and of the associated features of the surface EMG. For this purpose, we recorded high-density surface EMG signals during linearly increasing force contractions of the tibialis anterior muscle, up to 70% of maximal force. The recruitment threshold (RT), conduction velocity (MUCV), median frequency (MDF_MU), and amplitude (RMS_MU) of action potentials of 587 motor units from 13 individuals were assessed and associated with features of the interference EMG. MUCV was positively associated with RT (R² = 0.64 ± 0.14), whereas MDF_MU and RMS_MU showed a weaker relation with RT (R² = 0.11 ± 0.11 and 0.39 ± 0.24, respectively). Moreover, the changes in average conduction velocity estimated from the interference EMG predicted well the changes in MUCV (R² = 0.71), with a strong association to ankle dorsiflexion force (R² = 0.81 ± 0.12). Conversely, both the average EMG MDF and RMS were poorly associated with motor unit recruitment. These results clarify the limitations of EMG spectral and amplitude analysis in inferring the neural strategies of muscle control and indicate that, conversely, the average conduction velocity could provide relevant information on these strategies.NEW & NOTEWORTHY The surface EMG provides information on the neural drive to muscles. However, the associations between EMG features and neural drive have been long debated due to unavailability of motor unit population data. Here, by using novel highly accurate decomposition of the EMG, we related motor unit population behavior to a wide range of voluntary forces. The results fully clarify the potential and limitation of the surface EMG to provide estimates of the neural drive to muscles.

Comparison of human gastrocnemius forces predicted by Hill-type muscle models and estimated from ultrasound images

Hill-type models are ubiquitous in the field of biomechanics, providing estimates of a muscle's force as a function of its activation state and its assumed force-length and force-velocity properties. However, despite their routine use, the accuracy with which Hill-type models predict the forces generated by muscles during submaximal, dynamic tasks remains largely unknown. This study compared human gastrocnemius forces predicted by Hill-type models with the forces estimated from ultrasound-based measures of tendon length changes and stiffness during cycling, over a range of loads and cadences. We tested both a traditional model, with one contractile element, and a differential model, with two contractile elements that accounted for independent contributions of slow and fast muscle fibres. Both models were driven by subject-specific, ultrasound-based measures of fascicle lengths, velocities and pennation angles and by activation patterns of slow and fast muscle fibres derived from surface electromyographic recordings. The models predicted, on average, 54% of the time-varying gastrocnemius forces estimated from the ultrasound-based methods. However, differences between predicted and estimated forces were smaller under low speed-high activation conditions, with models able to predict nearly 80% of the gastrocnemius force over a complete pedal cycle. Additionally, the predictions from the Hill-type muscle models tested here showed that a similar pattern of force production could be achieved for most conditions with and without accounting for the independent contributions of different muscle fibre types.

Increased intensity and reduced frequency of EMG signals from feline self-reinnervated ankle extensors during walking do not normalize excessive lengthening

Kinematics of cat level walking recover after elimination of length-dependent sensory feedback from the major ankle extensor muscles induced by self-reinnervation. Little is known, however, about changes in locomotor myoelectric activity of self-reinnervated muscles. We examined the myoelectric activity of self-reinnervated muscles and intact synergists to determine the extent to which patterns of muscle activity change as almost normal walking is restored following muscle self-reinnervation. Nerves to soleus (SO) and lateral gastrocnemius (LG) of six adult cats were surgically transected and repaired. Intramuscular myoelectric signals of SO, LG, medial gastrocnemius (MG), and plantaris (PL), muscle fascicle length of SO and MG, and hindlimb mechanics were recorded during level and slope (±27°) walking before and after (10-12 wk postsurgery) self-reinnervation of LG and SO. Mean myoelectric signal intensity and frequency were determined using wavelet analysis. Following SO and LG self-reinnervation, mean myoelectric signal intensity increased and frequency decreased in most conditions for SO and LG as well as for intact synergist MG (P < 0.05). Greater elongation of SO muscle-tendon unit during downslope and unchanged magnitudes of ankle extensor moment during the stance phase in all walking conditions suggested a functional deficiency of ankle extensors after self-reinnervation. Possible effects of morphological reorganization of motor units of ankle extensors and altered sensory and central inputs on the changes in myoelectric activity of self-reinnervated SO and LG are discussed.

Validation of Hill-type muscle models in relation to neuromuscular recruitment and force-velocity properties: predicting patterns of in vivo muscle force

We review here the use and reliability of Hill-type muscle models to predict muscle performance under varying conditions, ranging from in situ production of isometric force to in vivo dynamics of muscle length change and force in response to activation. Muscle models are frequently used in musculoskeletal simulations of movement, particularly when applied to studies of human motor performance in which surgically implanted transducers have limited use. Musculoskeletal simulations of different animal species also are being developed to evaluate comparative and evolutionary aspects of locomotor performance. However, such models are rarely validated against direct measures of fascicle strain or recordings of muscle-tendon force. Historically, Hill-type models simplify properties of whole muscle by scaling salient properties of single fibers to whole muscles, typically accounting for a muscle's architecture and series elasticity. Activation of the model's single contractile element (assigned the properties of homogenous fibers) is also simplified and is often based on temporal features of myoelectric (EMG) activation recorded from the muscle. Comparison of standard one-element models with a novel two-element model and with in situ and in vivo measures of EMG, fascicle strain, and force recorded from the gastrocnemius muscles of goats shows that a two-element Hill-type model, which allows independent recruitment of slow and fast units, better predicts temporal patterns of in situ and in vivo force. Recruitment patterns of slow/fast units based on wavelet decomposition of EMG activity in frequency-time space are generally correlated with the intensity spectra of the EMG signals, the strain rates of the fascicles, and the muscle-tendon forces measured in vivo, with faster units linked to greater strain rates and to more rapid forces. Using direct measures of muscle performance to further test Hill-type models, whether traditional or more complex, remains critical for establishing their accuracy and essential for verifying their applicability to scientific and clinical studies of musculoskeletal function.

The effect of fast and slow motor unit activation on whole-muscle mechanical performance: the size principle may not pose a mechanical paradox

The output of skeletal muscle can be varied by selectively recruiting different motor units. However, our knowledge of muscle function is largely derived from muscle in which all motor units are activated. This discrepancy may limit our understanding of in vivo muscle function. Hence, this study aimed to characterize the mechanical properties of muscle with different motor unit activation. We determined the isometric properties and isotonic force-velocity relationship of rat plantaris muscles in situ with all of the muscle active, 30% of the muscle containing predominately slower motor units active or 20% of the muscle containing predominately faster motor units active. There was a significant effect of active motor unit type on isometric force rise time (p < 0.001) and the force-velocity relationship (p < 0.001). Surprisingly, force rise time was longer and maximum shortening velocity higher when all motor units were active than when either fast or slow motor units were selectively activated. We propose this is due to the greater relative effects of factors such as series compliance and muscle resistance to shortening during sub-maximal contractions. The findings presented here suggest that recruitment according to the size principle, where slow motor units are activated first and faster ones recruited as demand increases, may not pose a mechanical paradox, as has been previously suggested.

Early deactivation of slower muscle fibres at high movement frequencies

Animals produce rapid movements using fast cyclical muscle contractions. These types of movements are better suited to faster muscle fibres within muscles of mixed fibre types as they can shorten at faster velocities and achieve higher activation-deactivation rates than their slower counterparts. Preferential recruitment of faster muscle fibres has previously been shown during high velocity contractions. Additionally, muscle deactivation takes longer than activation and therefore may pose a limitation to fast cyclical contractions. It has been speculated that slower fibres maybe deactivated before faster fibres to accommodate their longer deactivation time. This study aimed to test whether shifts in muscle fibre recruitment occur with derecruitment of slow fibres before the faster fibres at high cycle frequencies. Electromyographic (EMG) signals were collected from the medial gastrocnemius at an extreme range of cycle frequencies and workloads. Wavelets were used to resolve the EMG signals into time and frequency space and the primary sources of variability within the EMG frequency spectra were identified through principal component analysis. A general early derecruitment of slower fibres was evident at the end of muscle excitation for the higher cycle frequencies, and additional slower fibre recruitment was present at the highest cycle frequency. The duration of muscle excitation reached a minimum of about 150 ms and did not change for the three highest cycle frequencies suggesting a duration limit for the medial gastrocnemius. This study provides further evidence of modifications of muscle fibre recruitment strategies to meet the mechanical demands of movement.

Accuracy of gastrocnemius muscles forces in walking and running goats predicted by one-element and two-element Hill-type models

Hill-type models are commonly used to estimate muscle forces during human and animal movement-yet the accuracy of the forces estimated during walking, running, and other tasks remains largely unknown. Further, most Hill-type models assume a single contractile element, despite evidence that faster and slower motor units, which have different activation-deactivation dynamics, may be independently or collectively excited. This study evaluated a novel, two-element Hill-type model with "differential" activation of fast and slow contractile elements. Model performance was assessed using a comprehensive data set (including measures of EMG intensity, fascicle length, and tendon force) collected from the gastrocnemius muscles of goats during locomotor experiments. Muscle forces predicted by the new two-element model were compared to the forces estimated using traditional one-element models and to the forces measured in vivo using tendon buckle transducers. Overall, the two-element model resulted in the best predictions of in vivo gastrocnemius force. The coefficient of determination, r(2), was up to 26.9% higher and the root mean square error, RMSE, was up to 37.4% lower for the two-element model than for the one-element models tested. All models captured salient features of the measured muscle force during walking, trotting, and galloping (r(2)=0.26-0.51), and all exhibited some errors (RMSE=9.63-32.2% of the maximum in vivo force). These comparisons provide important insight into the accuracy of Hill-type models. The results also show that incorporation of fast and slow contractile elements within muscle models can improve estimates of time-varying, whole muscle force during locomotor tasks.

Built for rowing: frog muscle is tuned to limb morphology to power swimming

Rowing is demanding, in part, because drag on the oars increases as the square of their speed. Hence, as muscles shorten faster, their force capacity falls, whereas drag rises. How do frogs resolve this dilemma to swim rapidly? We predicted that shortening velocity cannot exceed a terminal velocity where muscle and fluid torques balance. This terminal velocity, which is below Vmax, depends on gear ratio (GR = outlever/inlever) and webbed foot area. Perhaps such properties of swimmers are 'tuned', enabling shortening speeds of approximately 0.3Vmax for maximal power. Predictions were tested using a 'musculo-robotic' Xenopus laevis foot driven either by a living in vitro or computational in silico plantaris longus muscle. Experiments verified predictions. Our principle finding is that GR ranges from 11.5 to 20 near the predicted optimum for rowing (GR ≈ 11). However, gearing influences muscle power more strongly than foot area. No single morphology is optimal for producing muscle power. Rather, the 'optimal' GR decreases with foot size, implying that rowing ability need not compromise jumping (and vice versa). Thus, despite our neglect of additional forces (e.g. added mass), our model predicts pairings of physiological and morphological properties to confer effective rowing. Beyond frogs, the model may apply across a range of size and complexity from aquatic insects to human-powered rowing.

Estimating changes in metabolic power from EMG

Metabolic rates can increase 21 times above resting levels during cycling with the majority attributed to muscular contractions. Metabolic estimates attained through gas exchange parameters are limited by the respiration rate and time delay with respect to these contractions. In contrast surface electromyography (EMG) contains instantaneous muscle contraction information at higher temporal resolutions. An adequate metabolic power-EMG relationship has not been established to use EMG as a metabolic estimate during dynamic activities. The purpose of this study was to establish a metabolic power-EMG relationship during non steady-state conditions. Participants cycled at workloads between 25 and 90% O2max while EMG and gas exchange were monitored. The EMG was resolved into intensities and total EMG intensity was calculated as the sum of intensities across all muscles for each pedal cycle. Metabolic power was estimated from gas exchange parameters and the mean total EMG intensity between breaths was calculated and used as breath-by-breath values. Comparisons were made between breath-by-breath resolutions of metabolic power and total EMG intensity. Different weighting coefficients were also applied to the EMG for each muscle to analyze the effects of different muscle weightings on metabolic power estimations. There was a significant correlation (r = 0.91) between estimates of metabolic power from EMG and gas exchange. Muscle weighting had a significant effect on metabolic power determination with the highest and lowest correlated estimates having the largest weightings on muscles proximal and distal to the knee respectively. This study demonstrates that EMG contains important information about the metabolic costs of muscle contractions and provides good predictions of metabolic changes during non steady-state conditions. Also, the importance of each muscle is workload dependent with inappropriate weightings reducing metabolic estimations. These findings have implications for future EMG applications as they provide more immediate, higher temporal resolution predictions of changes in metabolic power.

Recruitment of faster motor units is associated with greater rates of fascicle strain and rapid changes in muscle force during locomotion

Animals modulate the power output needed for different locomotor tasks by changing muscle forces and fascicle strain rates. To generate the necessary forces, appropriate motor units must be recruited. Faster motor units have faster activation-deactivation rates than slower motor units, and they contract at higher strain rates; therefore, recruitment of faster motor units may be advantageous for tasks that involve rapid movements or high rates of work. This study identified motor unit recruitment patterns in the gastrocnemii muscles of goats and examined whether faster motor units are recruited when locomotor speed is increased. The study also examined whether locomotor tasks that elicit faster (or slower) motor units are associated with increased (or decreased) in vivo tendon forces, force rise and relaxation rates, fascicle strains and/or strain rates. Electromyography (EMG), sonomicrometry and muscle-tendon force data were collected from the lateral and medial gastrocnemius muscles of goats during level walking, trotting and galloping and during inclined walking and trotting. EMG signals were analyzed using wavelet and principal component analyses to quantify changes in the EMG frequency spectra across the different locomotor conditions. Fascicle strain and strain rate were calculated from the sonomicrometric data, and force rise and relaxation rates were determined from the tendon force data. The results of this study showed that faster motor units were recruited as goats increased their locomotor speeds from level walking to galloping. Slow inclined walking elicited EMG intensities similar to those of fast level galloping but different EMG frequency spectra, indicating that recruitment of the different motor unit types depended, in part, on characteristics of the task. For the locomotor tasks and muscles analyzed here, recruitment patterns were generally associated with in vivo fascicle strain rates, EMG intensity and tendon force. Together, these data provide new evidence that changes in motor unit recruitment have an underlying mechanical basis, at least for certain locomotor tasks.

Task-dependent activity of motor unit populations in feline ankle extensor muscles

Understanding the functional significance of the morphological diversity of mammalian skeletal muscles is limited by technical difficulties of estimating the contribution of motor units with different properties to unconstrained motor behaviours. Recently developed wavelet and principal components analysis of intramuscular myoelectric signals has linked signals with lower and higher frequency contents to the use of slower and faster motor unit populations. In this study we estimated the relative contributions of lower and higher frequency signals of cat ankle extensors (soleus, medial and lateral gastrocnemii, plantaris) during level, downslope and upslope walking and the paw-shake response. This was done using the first two myoelectric signal principal components (PCI, PCII), explaining over 90% of the signal, and an angle θ, a function of PCI/PCII, indicating the relative contribution of slower and faster motor unit populations. Mean myoelectric frequencies in all walking conditions were lowest for slow soleus (234 Hz) and highest for fast gastrocnemii (307 and 330 Hz) muscles. Motor unit populations within and across the studied muscles that demonstrated lower myoelectric frequency (suggesting slower populations) were recruited during tasks and movement phases with lower mechanical demands on the ankle extensors--during downslope and level walking and in early walking stance and paw-shake phases. With increasing mechanical demands (upslope walking, mid-phase of paw-shake cycles), motor unit populations generating higher frequency signals (suggesting faster populations) contributed progressively more. We conclude that the myoelectric frequency contents within and between feline ankle extensors vary across studied motor behaviours, with patterns that are generally consistent with muscle fibre-type composition.

Effects of homosynaptic depression on spectral properties of H-reflex recordings

The purpose of this study was to determine the effects of homosynaptic depression (HD) on spectral properties of the soleus (SOL) H-reflex. Paired stimulations, separated by 100 ms, were used to elicit an unconditioned and conditioned H-reflex in the SOL muscle of 20 participants during quiet standing. Wavelet and principal component analyses were used to analyze features of the time-varying spectral properties of the unconditioned and conditioned H-reflex. The effects of HD on spectral properties of the H-reflex signal were quantified by comparing extracted principal component scores. The analysis extracted two principal components: one associated with the intensity of the spectra and one associated with its frequency. The scores for both principal components were smaller for the conditioned H-reflex. HD decreases the spectral intensity and changes the spectral frequency of H-reflexes. These results suggest that HD changes the recruitment pattern of the motor units evoked during H-reflex stimulations, in that it not only decreases the intensity, but also changes the types of motor units that contribute to the H-reflex signal.

A muscle's force depends on the recruitment patterns of its fibers

Biomechanical models of whole muscles commonly used in simulations of musculoskeletal function and movement typically assume that the muscle generates force as a scaled-up muscle fiber. However, muscles are comprised of motor units that have different intrinsic properties and that can be activated at different times. This study tested whether a muscle model comprised of motor units that could be independently activated resulted in more accurate predictions of force than traditional Hill-type models. Forces predicted by the models were evaluated by direct comparison with the muscle forces measured in situ from the gastrocnemii in goats. The muscle was stimulated tetanically at a range of frequencies, muscle fiber strains were measured using sonomicrometry, and the activation patterns of the different types of motor unit were calculated from electromyographic recordings. Activation patterns were input into five different muscle models. Four models were traditional Hill-type models with different intrinsic speeds and fiber-type properties. The fifth model incorporated differential groups of fast and slow motor units. For all goats, muscles and stimulation frequencies the differential model resulted in the best predictions of muscle force. The in situ muscle output was shown to depend on the recruitment of different motor units within the muscle.

Spectral properties of H-reflex recordings after an acute bout of whole-body vibration

Although research supports the use of whole-body vibration (WBV) to improve neuromuscular performance, the mechanisms for these improvements remain unclear. The purpose of this study was to identify the effect of WBV on the spectral properties of electrically evoked H-reflex recordings in the soleus (SOL) muscle. The H-reflex recordings were measured in the SOL muscle of 20 participants before and after a bout of WBV. The H-reflexes were evoked every 15 seconds for 150 seconds after WBV. A wavelet procedure was used to extract spectral data, which were then quantified with a principle components analysis. Resultant principle component scores were used for statistical analysis. The analysis extracted 1 principle component associated with the intensity of the myoelectric spectra and 1 principle component associated with the frequency. The scores of the principle component that were related to the myoelectric intensity were smaller at 30 and 60 milliseconds after WBV than before WBV. The WBV transiently decreased the intensity of myoelectric spectra during electrically evoked contractions, but it did not influence the frequency of the spectra. The decrease in intensity likely indicates a smaller electrically evoked muscle twitch response, whereas the lack of change in frequency would indicate a similar recruitment pattern of motor units before and after WBV.