Patterns of motor recruitment can be determined using surface EMG

Children With Fragile X Syndrome Display a Switch Towards Fast Fibres in Their Recruitment Strategy During Gait

BACKGROUND

Fragile X Syndrome (FXS) is a genetic disorder caused by the lack of FMRP, a crucial protein for brain development and function. FMR1 mutations are categorized into premutation and full mutation (FXSFull), with somatic mosaicism (FXSMos) modulating the FXS phenotype. Recent studies identified muscle activity alterations during gait in FXS children. This study aims to explore the relationship between these muscle activity changes and motor fibre recruitment strategies during gait in FXS children.

METHODS

Fifty-four FXS children and fourteen healthy controls participated in the study. Gait trials at self-selected speeds were recorded using four synchronized cameras and a surface electromyography system that captured bilateral activity of Gastrocnemius lateralis, Tibialis anterior, Rectus and Biceps femoris muscles. The continuous wavelet transform, using the 'bump' mother wavelet, provided the percentage distribution of signal energy across nine frequency bands (50-Hz increments within a 450- to 10-Hz spectrum) and the Instantaneous MeaN Frequency (IMNF) time-frequency distribution.

RESULTS

Results indicated that both FXSFull and FXSMos children exhibit a distinct fibre recruitment strategy compared to controls, with a higher percentage of total energy and elevated IMNF (p < 0.05).

CONCLUSIONS

This increased reliance on fast-twitch fibres may contribute to the observed fatigability and exercise intolerance in FXS children.

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 diagnostic model of nerve root compression localization in lower lumbar disc herniation based on random forest algorithm and surface electromyography

Objective

This study aimed to investigate the muscle activation of patients with lumbar disc herniation (LDH) during walking by surface electromyography (SEMG) and establish a diagnostic model based on SEMG parameters using random forest (RF) algorithm for localization diagnosis of compressed nerve root in LDH patients.

Methods

Fifty-eight patients with LDH and thirty healthy subjects were recruited. The SEMG of tibialis anterior (TA) and lateral gastrocnemius (LG) were collected bilaterally during walking. The peak root mean square (RMS-peak), RMS-peak time, mean power frequency (MPF), and median frequency (MF) were analyzed. A diagnostic model based on SEMG parameters using RF algorithm was established to locate compressed nerve root, and repeated reservation experiments were conducted for verification. The study evaluated the diagnostic efficiency of the model using accuracy, precision, recall rate, F1-score, Kappa value, and area under the receiver operating characteristic (ROC) curve.

Results

The results showed that delayed activation of TA and decreased activation of LG were observed in the L5 group, while decreased activation of LG and earlier activation of LG were observed in the S1 group. The RF model based on eight SEMG parameters showed an average accuracy of 84%, with an area under the ROC curve of 0.93. The RMS peak time of TA was identified as the most important SEMG parameter.

Conclusion

These findings suggest that the RF model can assist in the localization diagnosis of compressed nerve roots in LDH patients, and the SEMG parameters can provide further references for optimizing the diagnosis model in the future.

Correlation Between Skin Autofluorescence and Muscle Activities of Lower Limb in Aging Without Disease and Disability

Skin autofluorescence is a useful index to estimate the accumulation of advanced glycation end-products in human tissues. Elderly persons with higher skin autofluorescence have lower muscle mass, muscle strength and muscle power, however, little is known about the relationship between the skin autofluorescence level and each muscle activity. We measured the values of skin autofluorescence from five places on a lower limb, and the signals of surface electromyogram during isometric contractions from five muscles on that, simultaneously. The waveforms of surface electromyogram were analyzed by Daubechies-4 wavelet transformation. The value of skin autofluorescence was increased in the proximal part of the lower limb compared with the value of the distal part. The principal component of surface electromyogram activity in a time-frequency domain was lower in the proximal part compared with that of the distal part. There was a weak negative correlation between the value of skin autofluorescence on the gluteal region and the value of the mean wavelet coefficient of the surface electromyogram signals within the gluteus maximus muscle. The higher accumulation of advanced glycation end-products on the gluteal region might suggest the lower muscle activity in aging without disease and disability.

Inter-muscular networks of synchronous muscle fiber activation

Skeletal muscles continuously coordinate to facilitate a wide range of movements. Muscle fiber composition and timing of activation account for distinct muscle functions and dynamics necessary to fine tune muscle coordination and generate movements. Here we address the fundamental question of how distinct muscle fiber types dynamically synchronize and integrate as a network across muscles with different functions. We uncover that physiological states are characterized by unique inter-muscular network of muscle fiber cross-frequency interactions with hierarchical organization of distinct sub-networks and modules, and a stratification profile of links strength specific for each state. We establish how this network reorganizes with transition from rest to exercise and fatigue-a complex process where network modules follow distinct phase-space trajectories reflecting their functional role in movements and adaptation to fatigue. This opens a new area of research, Network Physiology of Exercise, leading to novel network-based biomarkers of health, fitness and clinical conditions.

Changes in ankle work, foot work, and tibialis anterior activation throughout a long run

BACKGROUND

The ankle and foot together contribute to over half of the positive and negative work performed by the lower limbs during running. Yet, little is known about how foot kinetics change throughout a run. The amount of negative foot work may decrease as tibialis anterior (TA) electromyography (EMG) changes throughout longer-duration runs. Therefore, we examined ankle and foot work as well as TA EMG changes throughout a changing-speed run.

METHODS

Fourteen heel-striking subjects ran on a treadmill for 58 min. We collected ground reaction forces, motion capture, and EMG. Subjects ran at 110%, 100%, and 90% of their 10-km running speed and 2.8 m/s multiple times throughout the run. Foot work was evaluated using the distal rearfoot work, which provides a net estimate of all work contributors within the foot.

RESULTS

Positive foot work increased and positive ankle work decreased throughout the run at all speeds. At the 110% 10-km running speed, negative foot work decreased and TA EMG frequency shifted lower throughout the run. The increase in positive foot work may be attributed to increased foot joint work performed by intrinsic foot muscles. Changes in negative foot work and TA EMG frequency may indicate that the TA plays a role in negative foot work in the early stance of a run.

CONCLUSION

This study is the first to examine how the kinetic contributions of the foot change throughout a run. Future studies should investigate how increases in foot work affect running performance.

Universal spectral profile and dynamic evolution of muscle activation: a hallmark of muscle type and physiological state

The skeletal muscle is an integrated multicomponent system with complex dynamics of continuous myoelectrical activation of various muscle types across time scales to facilitate muscle coordination among units and adaptation to physiological states. To understand the multiscale dynamics of neuromuscular activity, we investigated spectral characteristics of different muscle types across time scales and their evolution with physiological states. We hypothesized that each muscle type is characterized by a specific spectral profile, reflecting muscle composition and function, that remains invariant over time scales and is universal across subjects. Furthermore, we hypothesized that the myoelectrical activation and corresponding spectral profile during certain movements exhibit an evolution path in time that is unique for each muscle type and reflects responses in muscle dynamics to exercise, fatigue, and aging. To probe the multiscale mechanism of neuromuscular regulation, we developed a novel protocol of repeated squat exercise segments, each performed until exhaustion, and we analyzed differentiated spectral power responses over a range of frequency bands for leg and back muscle activation in young and old subjects. We found that leg and back muscle activation is characterized by muscle-specific spectral profiles, with differentiated frequency band contribution, and a muscle-specific evolution path in response to fatigue and aging that is universal across subjects in each age group. The uncovered universality among subjects in the spectral profile of each muscle at a given physiological state, as well as the robustness in the evolution of these profiles over a range of time scales and states, reveals a previously unrecognized multiscale mechanism underlying the differentiated response of distinct muscle types to exercise-induced fatigue and aging.NEW & NOTEWORTHY To understand coordinated function of distinct fibers in a muscle, we investigated spectral dynamics of muscle activation during maximal exercise across a range of frequency bands and time scales of observation. We discovered a spectral profile that is specific for each muscle type, robust at short, intermediate, and large time scales, universal across subjects, and characterized by a muscle-specific evolution path with accumulation of fatigue and aging, indicating a previously unrecognized multiscale mechanism of muscle tone regulation.

Quantifying muscle coactivation in individuals with incomplete spinal cord injury using wavelets

BACKGROUND

Individuals with incomplete spinal cord injury often have decreased gait function and coactivation of antagonistic muscle pairs. Common ways of quantifying coactivation using electromyographic signals do not consider frequency information in the signal. As electromyographic signals from different motor unit types have different frequency components and muscle fiber type can change in individuals with spinal cord injury, it may be beneficial to consider frequency components. The aims were to demonstrate the utility of using a method which considers temporal and frequency components of the electromyographical signal to quantify coactivation in lower extremity muscles in individuals with incomplete spinal cord injury through 1) comparison with able-bodied individuals and 2) comparison before and after body weight supported treadmill training.

METHODS

Frequency decomposition techniques were applied to electromyographical signals to consider the temporal and frequency components of the electromyographical signals to quantify coactivation over a range of frequencies.

RESULTS

Our main findings show that correlation coefficients between total EMG intensities of rectus femoris-biceps femoris and medial gastrocnemius-tibialis anterior were significantly different between able-bodied individuals and those with incomplete spinal cord injury (p = 0006, p = 0.01). The correlation spectra of medial gastrocnemius-tibialis anterior of the spinal cord injury group were substantially different than those the able-bodied group, while the EMG normalcy score was significantly different (p = 0.002). We also found that there was a change in coactivation of ankle muscles after body weight supported treadmill training.

INTERPRETATION

Our findings indicate that there may be frequency specific differences in muscle coactivation between able-bodied individuals and those with incomplete spinal cord injury. Changes in coactivation were also observed before and after body weight supported treadmill training. These differences may reflect the changes in recruitment patterns of different motor unit types.

NIRS-EMG for Clinical Applications: A Systematic Review

In this review, we present an overview of the applications and computed parameters of electromyography (EMG) and near-infrared spectroscopy (NIRS) methods on patients in clinical practice. The eligible studies were those where both techniques were combined in order to assess muscle characteristics from the electrical and hemodynamic points of view. With this aim, a comprehensive screening of the literature based on related keywords in the most-used scientific data bases allowed us to identify 17 papers which met the research criteria. We also present a brief overview of the devices designed specifically for muscular applications with EMG and NIRS sensors (a total of eight papers). A critical analysis of the results of the review suggests that the combined use of EMG and NIRS on muscle has been only partially exploited for assessment and evaluation in clinical practice and, thus, this field shows promises for future developments.

Quantifying the velocity-dependent muscle response during gait of children with Cerebral Palsy

A new method is introduced quantifying the velocity-dependent muscle response during gait in spastic muscles of children with Cerebral Palsy. The velocity-dependent muscle activation Index is calculated during a 3-dimensional gait analysis using segment angular velocity and the Instantaneous Mean Frequency calculated from surface electromyography. Typical developed children (n = 11) and children with hemiplegia (n = 11) aging from 8 to 19 years participated in the study. The rectus femoris and the medial gastrocnemius were assessed by calculating the velocity dependent muscle activation Index and the modified Ashworth Scale. Greater velocity-dependent muscle activation Index values for both medial gastrocnemius and rectus femoris muscles were associated with greater Ashworth Scale. Post hoc analysis revealed significant lower velocity-dependent muscle activation Index means in the Typical developed group compared with Ashworth Scale scores of 1, 2, 3, and 5. In addition, velocity-dependent muscle activation Index for Ashworth Scale 0, 1, and 2 were significantly lower than for Ashworth Scale 3 and 5. The velocity dependent muscle activation Index showed negative low correlation with walking speed and cadence. Findings show that spastic muscles can be quantified during dynamic functional task such as walking. Future studies should investigate the reliability of the velocity-dependent muscle activation Index that may be used for the assessment of spasticity management such as Botulinum toxin A interventions.

How Do the Mechanical Demands of Cycling Affect the Information Content of the EMG?

PURPOSE

The persistence of phase-related information in EMG signals can be quantified by its entropic half-life (EnHL). It has been proposed that the EnHL would increase with the demands of a movement task, and thus increase as the pedaling power increased during cycling. However, simulation work on the properties of EMG signals suggests that the EnHL depends on burst duration and duty cycle in the EMG that may not be related to task demands. This study aimed to distinguish between these alternate hypotheses.

METHODS

The EnHL was characterized for 10 muscles from nine cyclists cycling at a range of powers (35 to 260 W) and cadences (60-140 rpm) for the raw EMG, phase-randomized surrogate EMG, EMG intensity, and the principal components describing the muscle coordination patterns.

RESULTS

There was phase-related information in the raw EMG signals and EMG intensities that was related to the EMG burst duration, duty cycle pedaling cadence, and power. The EnHL for the EMG intensities of the individual muscles (excluding quadriceps) and for the coordination patterns decreased as cycling power and cadence increased.

CONCLUSIONS

The EnHL provide information on the structure of the motor control signals and their constituent motor unit action potentials, both within and between muscles, rather than on the mechanical demands of the cycling task per se.

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.

Energy cost and lower leg muscle activities during erect bipedal locomotion under hyperoxia

Background

Energy cost of transport per unit distance (CoT) against speed shows U-shaped fashion in walking and linear fashion in running, indicating that there exists a specific walking speed minimizing the CoT, being defined as economical speed (ES). Another specific gait speed is the intersection speed between both fashions, being called energetically optimal transition speed (EOTS). We measured the ES, EOTS, and muscle activities during walking and running at the EOTS under hyperoxia (40% fraction of inspired oxygen) on the level and uphill gradients (+ 5%).

Methods

Oxygen consumption \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left(\dot{V}{\mathrm{O}}_2\right) $$\end{document}V˙O2 and carbon dioxide output \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left(\dot{V}{\mathrm{CO}}_2\right) $$\end{document}V˙CO2 were measured to calculate the CoT values at eight walking speeds (2.4–7.3 km h⁻¹) and four running speeds (7.3–9.4 km h^− 1) in 17 young males. Electromyography was recorded from gastrocnemius medialis, gastrocnemius lateralis (GL), and tibialis anterior (TA) to evaluate muscle activities. Mean power frequency (MPF) was obtained to compare motor unit recruitment patterns between walking and running.

Results

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \dot{V}{\mathrm{O}}_2 $$\end{document}V˙O2, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \dot{V}{\mathrm{CO}}_2 $$\end{document}V˙CO2, and CoT values were lower under hyperoxia than normoxia at faster walking speeds and any running speeds. A faster ES on the uphill gradient and slower EOTS on both gradients were observed under hyperoxia than normoxia. GL and TA activities became lower when switching from walking to running at the EOTS under both FiO₂ conditions on both gradients, so did the MPF in the TA.

Conclusions

ES and EOTS were influenced by reduced metabolic demands induced by hyperoxia. GL and TA activities in association with a lower shift of motor unit recruitment patterns in the TA would be related to the gait selection when walking or running at the EOTS.

Trial registration

UMIN000017690 (R000020501). Registered May 26, 2015, before the first trial.

Electronic supplementary material

The online version of this article (10.1186/s40101-018-0177-7) contains supplementary material, which is available to authorized users.

A wavelet based time frequency analysis of electromyograms to group steps of runners into clusters that contain similar muscle activation patterns

PURPOSE

To wavelet transform the electromyograms of the vastii muscles and generate wavelet intensity patterns (WIP) of runners. Test the hypotheses: 1) The WIP of the vastus medialis (VM) and vastus lateralis (VL) of one step are more similar than the WIPs of these two muscles, offset by one step. 2) The WIPs within one muscle differ by having maximal intensities in specific frequency bands and these intensities are not always occurring at the same time after heel strike. 3) The WIPs that were recorded form one muscle for all steps while running can be grouped into clusters with similar WIPs. It is expected that clusters might have distinctly different, cluster specific mean WIPs.

METHODS

The EMG of the vastii muscles from at least 1000 steps from twelve runners were recorded using a bipolar current amplifier and yielded WIPs. Based on the weights obtained after a principal component analysis the dissimilarities (1-correlation) between the WIPs were computed. The dissimilarities were submitted to a hierarchical cluster analysis to search for groups of steps with similar WIPs. The clusters formed by random surrogate WIPs were used to determine whether the groups were likely to be created in a non-random manner.

RESULTS

The steps were grouped in clusters showing similar WIPs. The grouping was based on the frequency bands and their timing showing that they represented defining parts of the WIPs. The correlations between the WIPs of the vastii muscles that were recorded during the same step were higher than the correlations of WPIs that were recorded during consecutive steps, indicating the non-randomness of the WIPs.

CONCLUSIONS

The spectral power of EMGs while running varies during the stance phase in time and frequency, therefore a time averaged power spectrum cannot reflect the timing of events that occur while running. It seems likely that there might be a set of predefined patterns that are used upon demand to stabilize the movement.

Movement Complexity and Neuromechanical Factors Affect the Entropic Half-Life of Myoelectric Signals

Appropriate neuromuscular functioning is essential for survival and features underpinning motor control are present in myoelectric signals recorded from skeletal muscles. One approach to quantify control processes related to function is to assess signal variability using measures such as Sample Entropy. Here we developed a theoretical framework to simulate the effect of variability in burst duration, activation duty cycle, and intensity on the Entropic Half-Life (EnHL) in myoelectric signals. EnHLs were predicted to be <40 ms, and to vary with fluctuations in myoelectric signal amplitude and activation duty cycle. Comparison with myoelectic data from rats walking and running at a range of speeds and inclines confirmed the range of EnHLs, however, the direction of EnHL change in response to altered locomotor demand was not correctly predicted. The discrepancy reflected different associations between the ratio of the standard deviation and mean signal intensity (Ist:It¯) and duty factor in simulated and physiological data, likely reflecting additional information in the signals from the physiological data (e.g., quiescent phase content; variation in action potential shapes). EnHL could have significant value as a novel marker of neuromuscular responses to alterations in perceived locomotor task complexity and intensity.

Alterations in Spectral Attributes of Surface Electromyograms after Utilization of a Foot Drop Stimulator during Post-Stroke Gait

BACKGROUND

A foot drop stimulator (FDS) is a rehabilitation intervention that stimulates the common peroneal nerve to facilitate ankle dorsiflexion at the appropriate time during post-stroke hemiplegic gait. Time-frequency analysis (TFA) of non-stationary surface electromyograms (EMG) and spectral variables such as instantaneous mean frequency (IMNF) can provide valuable information on the long-term effects of FDS intervention in terms of changes in the motor unit (MU) recruitment during gait, secondary to improved dorsiflexion.

OBJECTIVE

The aim of this study was to apply a wavelet-based TFA approach to assess the changes in neuromuscular activation of the tibialis anterior (TA), soleus (SOL), and gastrocnemius (GA) muscles after utilization of an FDS during gait post-stroke.

METHODS

Surface EMG were collected bilaterally from the TA, SOL, and GA muscles from six participants (142.9 ± 103.3 months post-stroke) while walking without the FDS at baseline and 6 months post-FDS utilization. Continuous wavelet transform was performed to get the averaged time-frequency distribution of band pass filtered (20-300 Hz) EMGs during multiple walking trials. IMNFs were computed during normalized gait and were averaged during the stance and swing phases. Percent changes in the energies associated with each frequency band of 25 Hz between 25 and 300 Hz were computed and compared between visits.

RESULTS

Averaged time-frequency representations of the affected TA, SOL, and GA EMG show altered spectral attributes post-FDS utilization during normalized gait. The mean IMNF values for the affected TA were significantly lower than the unaffected TA at baseline (p = 0.026) and follow-up (p = 0.038) during normalized stance. The mean IMNF values significantly increased (p = 0.017) for the affected GA at follow-up during normalized swing. The frequency band of 250-275 Hz significantly increased in the energies post-FDS utilization for all muscles.

CONCLUSION

The application of wavelet-based TFA of EMG and outcome measures (IMNF, energy) extracted from the time-frequency distributions suggest alterations in MU recruitment strategies after the use of FDS in individuals with chronic stroke. This further establishes the efficacy of FDS as a rehabilitation intervention that may promote motor recovery in addition to treating the secondary complications of foot drop due to post-stroke hemiplegia.

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.

Warm-Up Strategies for Sport and Exercise: Mechanisms and Applications

It is widely accepted that warming-up prior to exercise is vital for the attainment of optimum performance. Both passive and active warm-up can evoke temperature, metabolic, neural and psychology-related effects, including increased anaerobic metabolism, elevated oxygen uptake kinetics and post-activation potentiation. Passive warm-up can increase body temperature without depleting energy substrate stores, as occurs during the physical activity associated with active warm-up. While the use of passive warm-up alone is not commonplace, the idea of utilizing passive warming techniques to maintain elevated core and muscle temperature throughout the transition phase (the period between completion of the warm-up and the start of the event) is gaining in popularity. Active warm-up induces greater metabolic changes, leading to increased preparedness for a subsequent exercise task. Until recently, only modest scientific evidence was available supporting the effectiveness of pre-competition warm-ups, with early studies often containing relatively few participants and focusing mostly on physiological rather than performance-related changes. External issues faced by athletes pre-competition, including access to equipment and the length of the transition/marshalling phase, have also frequently been overlooked. Consequently, warm-up strategies have continued to develop largely on a trial-and-error basis, utilizing coach and athlete experiences rather than scientific evidence. However, over the past decade or so, new research has emerged, providing greater insight into how and why warm-up influences subsequent performance. This review identifies potential physiological mechanisms underpinning warm-ups and how they can affect subsequent exercise performance, and provides recommendations for warm-up strategy design for specific individual and team sports.

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.

A survey of sensor devices: use in sports biomechanics

This paper examines the use of sensor devices in sports biomechanics, focusing on current frequency of use of Electromyography (EMG) device preferences. Researchers in the International Society of Biomechanics in Sports were invited to participate in an online survey. Responses on multiple sensor devices highlighting frequency of use, device features and improvements researchers sought in acquisition and analysis methods were obtained via an online questionnaire. Results of the investigation showed that the force platform is the most frequently used device, with inertial measurement units and EMG devices growing in popularity. Wireless functionality and ease of use for both the participant and the practitioner proved to be important features. The main findings of the survey demonstrated need for a simple, low power, multi-channel device which incorporates the various sensors into one single device. Biomechanists showed they were looking for more availability of wireless sensor devices with acquisition and analysis features. The study found there is a need to develop software analysis tools to accompany the multi-channel device, providing all the basic functions while maintaining compatibility with existing systems.

A Finite Element Model Approach to Determine the Influence of Electrode Design and Muscle Architecture on Myoelectric Signal Properties

Introduction

Surface electromyography (sEMG) is the measurement of the electrical activity of the skeletal muscle tissue detected at the skin’s surface. Typically, a bipolar electrode configuration is used. Most muscles have pennate and/or curved fibres, meaning it is not always feasible to align the bipolar electrodes along the fibres direction. Hence, there is a need to explore how different electrode designs can affect sEMG measurements.

Method

A three layer finite element (skin, fat, muscle) muscle model was used to explore different electrode designs. The implemented model used as source signal an experimentally recorded intramuscular EMG taken from the biceps brachii muscle of one healthy male. A wavelet based intensity analysis of the simulated sEMG signal was performed to analyze the power of the signal in the time and frequency domain.

Results

The model showed muscle tissue causing a bandwidth reduction (to 20-92- Hz). The inter-electrode distance (IED) and the electrode orientation relative to the fibres affected the total power but not the frequency filtering response. The effect of significant misalignment between the electrodes and the fibres (60°- 90°) could be reduced by increasing the IED (25–30 mm), which attenuates signal cancellation. When modelling pennated fibres, the muscle tissue started to act as a low pass filter. The effect of different IED seems to be enhanced in the pennated model, while the filtering response is changed considerably only when the electrodes are close to the signal termination within the model. For pennation angle greater than 20°, more than 50% of the source signal was attenuated, which can be compensated by increasing the IED to 25 mm.

Conclusion

Differences in tissue filtering properties, shown in our model, indicates that different electrode designs should be considered for muscle with different geometric properties (i.e. pennated muscles).

One-leg hop kinematics 20 years following anterior cruciate ligament rupture: Data revisited using functional data analysis

BACKGROUND

Despite interventions, anterior cruciate ligament ruptures can cause long-term deficits. To assist in identifying and treating deficiencies, 3D-motion analysis is used for objectivizing data. Conventional statistics are commonly employed to analyze kinematics, reducing continuous data series to discrete variables. Conversely, functional data analysis considers the entire data series.

METHODS

Here, we employ functional data analysis to examine and compare the entire time-domain of knee-kinematic curves from one-leg hops between and within three groups. All subjects (n=95) were part of a long-term follow-up study involving anterior cruciate ligament ruptures treated ~20 years ago conservatively with physiotherapy only or with reconstructive surgery and physiotherapy, and matched knee-healthy controls.

FINDINGS

Between-group differences (injured leg, treated groups; non-dominant leg, controls) were identified during the take-off and landing phases, and in the sagittal (flexion/extension) rather than coronal (abduction/adduction) and transverse (internal/external) planes. Overall, surgical and control groups demonstrated comparable knee-kinematic curves. However, compared to controls, the physiotherapy-only group exhibited less flexion during the take-off (0-55% of the normalized phase) and landing (44-73%) phase. Between-leg differences were absent in controls and the surgically treated group, but observed during the flight (4-22%, injured leg>flexion) and the landing (57-85%, injured leg<internal rotation) phases in the physiotherapy-only group.

INTERPRETATION

Functional data analysis identified specific functional knee-joint deviations from controls persisting 20 years post anterior cruciate ligament rupture, especially when treated conservatively. This approach is suggested as a means for comprehensively analyzing complex movements, adding to previous analyses.

Muscle coordination limits efficiency and power output of human limb movement under a wide range of mechanical demands

This study investigated the influence of cycle frequency and workload on muscle coordination and the ensuing relationship with mechanical efficiency and power output of human limb movement. Eleven trained cyclists completed an array of cycle frequency (cadence)-power output conditions while excitation from 10 leg muscles and power output were recorded. Mechanical efficiency was maximized at increasing cadences for increasing power outputs and corresponded to muscle coordination and muscle fiber type recruitment that minimized both the total muscle excitation across all muscles and the ineffective pedal forces. Also, maximum efficiency was characterized by muscle coordination at the top and bottom of the pedal cycle and progressive excitation through the uniarticulate knee, hip, and ankle muscles. Inefficiencies were characterized by excessive excitation of biarticulate muscles and larger duty cycles. Power output and efficiency were limited by the duration of muscle excitation beyond a critical cadence (120-140 rpm), with larger duty cycles and disproportionate increases in muscle excitation suggesting deteriorating muscle coordination and limitations of the activation-deactivation capabilities. Most muscles displayed systematic phase shifts of the muscle excitation relative to the pedal cycle that were dependent on cadence and, to a lesser extent, power output. Phase shifts were different for each muscle, thereby altering their mechanical contribution to the pedaling action. This study shows that muscle coordination is a key determinant of mechanical efficiency and power output of limb movement across a wide range of mechanical demands and that the excitation and coordination of the muscles is limited at very high cycle frequencies.

Ballistic exercise as a pre-activation stimulus: a review of the literature and practical applications

Post-activation potentiation (PAP) refers to the acute enhancement of muscular function as a direct result of its contractile history. Protocols designed to elicit PAP have commonly employed heavy resistance exercise (HRE) as the pre-activation stimulus; however, a growing body of research suggests that low-load ballistic exercises (BE) may also provide an effective stimulus. The ability to elicit PAP without the need for heavy equipment would make it easier to utilise prior to competition. It is hypothesised that BE can induce PAP given the high recruitment of type II muscle fibres associated with its performance. The literature has reported augmentations in power performance typically ranging from 2 to 5 %. The performance effects of BE are modulated by loading, recovery and physical characteristics. Jumps performed with an additional loading, such as depth jumps or weighted jumps, appear to be the most effective activities for inducing PAP. Whilst the impact of recovery duration on subsequent performance requires further research, durations of 1-6 min have been prescribed successfully in multiple instances. The effect of strength and sex on the PAP response to BE is not yet clear. Direct comparisons of BE and HRE, to date, suggest a tendency for HRE protocols to be more effective; future research should consider that these strategies must be optimised in different ways. The role of acute augmentations in lower limb stiffness is proposed as an additional mechanism that may further explain the PAP response following BE. In summary, BE demonstrates the potential to enhance performance in power tasks such as jumps and sprints. This review provides the reader with some practical recommendations for the application of BE as a pre-activation stimulus.

Peak torque and rate of torque development influence on repeated maximal exercise performance: contractile and neural contributions

Rapid force production is critical to improve performance and prevent injuries. However, changes in rate of force/torque development caused by the repetition of maximal contractions have received little attention. The aim of this study was to determine the relative influence of rate of torque development (RTD) and peak torque (T_peak) on the overall performance (i.e. mean torque, T_mean) decrease during repeated maximal contractions and to investigate the contribution of contractile and neural mechanisms to the alteration of the various mechanical variables. Eleven well-trained men performed 20 sets of 6-s isokinetic maximal knee extensions at 240°·s^-1, beginning every 30 seconds. RTD, T_peak and T_mean as well as the Rate of EMG Rise (RER), peak EMG (EMG_peak) and mean EMG (EMG_mean) of the vastus lateralis were monitored for each contraction. A wavelet transform was also performed on raw EMG signal for instant mean frequency (if_mean) calculation. A neuromuscular testing procedure was carried out before and immediately after the fatiguing protocol including evoked RTD (eRTD) and maximal evoked torque (eT_peak) induced by high frequency doublet (100 Hz). T_mean decrease was correlated to RTD and T_peak decrease (R²=0.62; p<0.001; respectively β=0.62 and β=0.19). RER, eRTD and initial if_mean (0-225 ms) decreased after 20 sets (respectively -21.1±14.1, -25±13%, and ~20%). RTD decrease was correlated to RER decrease (R²=0.36; p<0.05). The eT_peak decreased significantly after 20 sets (24±5%; p<0.05) contrary to EMG_peak (-3.2±19.5 %; p=0.71). Our results show that reductions of RTD explained part of the alterations of the overall performance during repeated moderate velocity maximal exercise. The reductions of RTD were associated to an impairment of the ability of the central nervous system to maximally activate the muscle in the first milliseconds of the contraction.

Regionalizing muscle activity causes changes to the magnitude and direction of the force from whole muscles-a modeling study

Skeletal muscle can contain neuromuscular compartments that are spatially distinct regions that can receive relatively independent levels of activation. This study tested how the magnitude and direction of the force developed by a whole muscle would change when the muscle activity was regionalized within the muscle. A 3D finite element model of a muscle with its bounding aponeurosis was developed for the lateral gastrocnemius, and isometric contractions were simulated for a series of conditions with either a uniform activation pattern, or regionally distinct activation patterns: in all cases the mean activation from all fibers within the muscle reached 10%. The models showed emergent features of the fiber geometry that matched physiological characteristics: with fibers shortening, rotating to greater pennation, adopting curved trajectories in 3D and changes in the thickness and width of the muscle belly. Simulations were repeated for muscle with compliant, normal and stiff aponeurosis and the aponeurosis stiffness affected the changes to the fiber geometry and the resultant muscle force. Changing the regionalization of the activity resulted to changes in the magnitude, direction and center of the force vector from the whole muscle. Regionalizing the muscle activity resulted in greater muscle force than the simulation with uniform activity across the muscle belly. The study shows how the force from a muscle depends on the complex interactions between the muscle fibers and connective tissues and the region of muscle that is active.

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.

Motor strategy patterns study of diabetic neuropathic individuals while walking. A wavelet approach

The aim of this study was to investigate muscle׳s energy patterns and spectral properties of diabetic neuropathic individuals during gait cycle using wavelet approach. Twenty-one diabetic patients diagnosed with peripheral neuropathy, and 21 non-diabetic individuals were assessed during the whole gait cycle. Activation patterns of vastus lateralis, medial gastrocnemius and tibialis anterior were studied by means of bipolar surface EMG. The signal׳s energy and frequency were compared between groups using t-test. The energy was compared in each frequency band (7-542 Hz) using ANOVAs for repeated measures for each group and each muscle. The diabetic individuals displayed lower energies in lower frequency bands for all muscles and higher energies in higher frequency bands for the extensors׳ muscles. They also showed lower total energy of gastrocnemius and a higher total energy of vastus, considering the whole gait cycle. The overall results suggest a change in the neuromuscular strategy of the main extensor muscles of the lower limb of diabetic patients to compensate the ankle extensor deficit to propel the body forward and accomplish the walking task.

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.

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.

Contraction Intensity and Velocity on Vastus Lateralis Semg Power Spectrum and Amplitude

Effect of contraction intensity [100%, 75%, 50%, and 25% maximum voluntary contraction (MVC)] and movement velocity [0° (isometric)], 50°, 100°, 200°, and 400°/sec. [isovelocities]) on root mean square amplitude (SEMG–RMS) and median frequency power spectrum (SEMG–MNF) of vastus lateralis (VL) surface electromyography was investigated with ten healthy female university students. Peak torque (PT), mean torque (MT), SEMG–MNF, and SEMG–RMS, analyzed using separate repeated-measures analyses of variance ( p ≤ .05), indicated: (1) an inverse relation between PT and MT and movement velocity, (2) greater SEMG–MNF values during all isovelocity conditions compared with isometric conditions, with highest values occurring at 50°/sec. and at 100% and 75% MVC, and (3) at all contraction intensities SEMG–RMS values were higher during dynamic movements than isometric movements and highest at 200° / sec. Isovelocity contractions were inferred to facilitate a greater recruitment of fast-twitch fibers (via increased SEMG–MNF), which was intensified at 50°/sec, whereas greater overall muscle activation was found at 200° / sec.

Shoulder muscles recruitment during a power backward giant swing on high bar: a wavelet-EMG-analysis

This study aimed at determining the upper limb muscles coordination during a power backward giant swing (PBGS) and the recruitment pattern of motor units (MU) of co-activated muscles. The wavelet transformation (WT) was applied to the surface electromyographic (EMG) signal of eight shoulder muscles. Total gymnast's body energy and wavelet synergies extracted from the WT-EMG by using a non-negative matrix factorization were analyzed as a function of the body position angle of the gymnast. A cross-correlation analysis of the EMG patterns allowed determining two main groups of co-activated muscles. Two wavelet synergies representing the main spectral features (82% of the variance accounted for) discriminated the recruitment of MU. Although no task-group of MU was found among the muscles, it appeared that a higher proportion of fast MU was recruited within the muscles of the first group during the upper part of the PBGS. The last increase of total body energy before bar release was induced by the recruitment of the muscles of the second group but did not necessitate the recruitment of a higher proportion of fast MU. Such muscle coordination agreed with previous simulations of elements on high bar as well as the findings related to the recruitment of MU.

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 electromyographic and mechanomyographic signals during dynamic concentric and eccentric contractions of the human biceps brachii muscle

The purpose of this study was to describe and examine the variations in recruitment patterns of motor units (MUs) in biceps brachii (BB) through a range of joint motion during dynamic eccentric and concentric contractions. Twelve healthy participants (6 females, 6 males, age=30±8.5 years) performed concentric and eccentric contractions with constant external loading at different levels. Surface electromyography (EMG) and mechanomyography (MMG) were recorded from BB. The EMGs and MMGs were decomposed into their intensities in time-frequency space using a wavelet technique. The EMG and MMG spectra were then compared using principal component analysis. Variations in total intensity, first principal component (PCI), and the angle θ formed by first component (PCI) and second component (PCII) loading scores were explained in terms of MU recruitment patterns and elbow angles. Elbow angle had a significant effect on dynamic concentric and eccentric contractions. The EMG total intensity was greater for concentric than for eccentric contractions in the present study. MMG total intensity, however, was lower during concentric than during eccentric contractions. In addition, there was no significant difference in θ between concentric and eccentric contractions for both EMG and MMG. Selective recruitment of fast MUs from BB muscle during eccentric muscle contractions was not found in the present study.

Movement mechanics as a determinate of muscle structure, recruitment and coordination

During muscle contractions, the muscle fascicles may shorten at a rate different from the muscle-tendon unit, and the ratio of these velocities is its gearing. Appropriate gearing allows fascicles to reduce their shortening velocities and allows them to operate at effective shortening velocities across a range of movements. Gearing of the muscle fascicles within the muscle belly is the result of rotations of the fascicles and bulging of the belly. Variable gearing can also occur as a result of tendon length changes that can be caused by changes in the relative timing of muscle activity for different mechanical tasks. Recruitment patterns of slow and fast fibres are crucial for achieving optimal muscle performance, and coordination between muscles is related to whole limb performance. Poor coordination leads to inefficiencies and loss of power, and optimal coordination is required for high power outputs and high mechanical efficiencies from the limb. This paper summarizes key studies in these areas of neuromuscular mechanics and results from studies where we have tested these phenomena on a cycle ergometer are presented to highlight novel insights. The studies show how muscle structure and neural activation interact to generate smooth and effective motion of the body.

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

Recruitment patterns and activation dynamics of different motor units greatly influence the temporal pattern and magnitude of muscle force development, yet these features are not often considered in muscle models. The purpose of this study was to characterize the recruitment and activation dynamics of slow and fast motor units from electromyographic (EMG) recordings and twitch force profiles recorded directly from animal muscles. EMG and force data from the gastrocnemius muscles of seven goats were recorded during in vivo tendon-tap reflex and in situ nerve stimulation experiments. These experiments elicited EMG signals with significant differences in frequency content (p<0.001). The frequency content was characterized using wavelet and principal components analysis, and optimized wavelets with centre frequencies, 149.94 Hz and 323.13 Hz, were obtained. The optimized wavelets were used to calculate the EMG intensities and, with the reconstructed twitch force profiles, to derive transfer functions for slow and fast motor units that estimate the activation state of the muscle from the EMG signal. The resulting activation-deactivation time constants gave r values of 0.98-0.99 between the activation state and the force profiles. This work establishes a framework for developing improved muscle models that consider the intrinsic properties of slow and fast fibres within a mixed muscle, and that can more accurately predict muscle force output from EMG.

Muscular timing and inter-muscular coordination in healthy females while walking

The dynamic interplay between muscles surrounding the knee joint, the central nervous system and external factors require a control strategy to generate and stabilise the preferred gait pattern. The electromyographic (EMG) signal is a common measure reflecting the neuromuscular control strategies during dynamic tasks. Neuromuscular control mechanisms, found in processed EMG signals, showed a precise pacing with a pacing rhythm and a tight control of muscle activity in running and maximally contracted muscles. The purpose of this study was to provide an insight how muscles get activated during walking. The EMG power, extracted by the wavelet transform (92-395Hz), over a time period encompassing 250ms before and 250ms after heel strike was analysed. The study showed that the wavelet-based analysis of EMG signals was sufficiently sensitive to detect a synchronisation of the activation of thigh muscles while walking. The results within each single subject and within the group consisting of 10 healthy females showed that, although there was a lot of jitter in the locations of the intensity peaks, the muscle activation is controlled, on average, by a neuromuscular activity paced at about 40ms, however with variable amplitudes. Albeit the jitter of the signal, the results resolved the temporal dependency of intensity peaks within muscles surrounding the knee and provided an insight into neural control of locomotion. The methodology to assess the stabilising muscle activation pattern may provide a way to discriminate subjects with normal gait pattern form those with a deteriorated neuromuscular control strategy.

Isometric fatigue patterns in time and time-frequency domains of triceps surae muscle in different knee positions

The occurrence of fatigue in triceps surae (TS) muscles during sustained plantar flexion contraction is investigated by means of the RMS electromyogram (EMG) and the instantaneous median frequency (IMF) of the short time Fourier transform (STFT). Six male subjects realized a 40% maximal plantar flexion isometric voluntary contraction until fatigue in two knee positions. Electrodes were positioned on gastrocnemius medialis, gastrocnemius lateralis and soleus muscles. The torque (TO) and EMG signals were synchronized. The RMS and the median of the IMF values were obtained, respectively, for each 250 ms and 1s windows of signal. Each signal was segmented into 10 epochs, from which the mean values of IMF, RMS and TO were obtained and submitted to linear regressions to determine parameter trends. Friedman test with the Dunn's post hoc were used to test for differences among muscles activation for each knee position and among slopes of regression curves, as well as to observe changes in TS RMS values over time. The results indicate different activation strategies with the knee extended (KE) in contrast to knee flexed (KF). With the KE, the gastrocnemii showed typical fatigue behavior with significant (p<0.05) IMF reductions and RMS increases over time, while soleus showed concomitant RMS and IMF increases (p<0.05) suggesting an increased soleus contribution to the torque production. With KF, the gastrocnemii were under activated, increasing the role of soleus. Thus, time-frequency analysis represented an important tool for TS muscular fatigue evaluation, allowing differentiates the role of soleus muscle.