Characteristics of surface mechanomyogram are dependent on development of fusion of motor units in humans

Accurate Identification of Motoneuron Discharges From Ultrasound Images Across the Full Muscle Cross-Section

OBJECTIVE

Non-invasive identification of motoneuron (MN) activity commonly uses electromyography (EMG). However, surface EMG (sEMG) detects only superficial sources, at less than approximately 10-mm depth. Intramuscular EMG can detect deep sources, but it is limited to sources within a few mm of the detection site. Conversely, ultrasound (US) images have high spatial resolution across the whole muscle cross-section. The activity of MNs can be extracted from US images due to the movements that MN activation generates in the innervated muscle fibers. Current US-based decomposition methods can accurately identify the location and average twitch induced by MN activity. However, they cannot accurately detect MN discharge times.

METHODS

Here, we present a method based on the convolutive blind source separation of US images to estimate MN discharge times with high accuracy. The method was validated across Ten participants using concomitant sEMG decomposition as the ground truth.

RESULTS

140 unique MN spike trains were identified from US images, with a rate of agreement (RoA) with sEMG decomposition of 87.4 ± 10.3%. Over 50% of these MN spike trains had a RoA greater than 90%. Furthermore, with US, we identified additional MUs well beyond the sEMG detection volume, at up to >30 mm below the skin.

CONCLUSION

The proposed method can identify discharges of MNs innervating muscle fibers in a large range of depths within the muscle from US images.

SIGNIFICANCE

The proposed methodology can non-invasively interface with the outer layers of the central nervous system innervating muscles across the full cross-section.

Non-Linearity in Motor Unit Velocity Twitch Dynamics: Implications for Ultrafast Ultrasound Source Separation

Ultrasound (US) muscle image series can be used for peripheral human-machine interfacing based on global features, or even on the decomposition of US images into the contributions of individual motor units (MUs). With respect to state-of-the-art surface electromyography (sEMG), US provides higher spatial resolution and deeper penetration depth. However, the accuracy of current methods for direct US decomposition, even at low forces, is relatively poor. These methods are based on linear mathematical models of the contributions of MUs to US images. Here, we test the hypothesis of linearity by comparing the average velocity twitch profiles of MUs when varying the number of other concomitantly active units. We observe that the velocity twitch profile has a decreasing peak-to-peak amplitude when tracking the same target motor unit at progressively increasing contraction force levels, thus with an increasing number of concomitantly active units. This observation indicates non-linear factors in the generation model. Furthermore, we directly studied the impact of one MU on a neighboring MU, finding that the effect of one source on the other is not symmetrical and may be related to unit size. We conclude that a linear approximation is partly limiting the decomposition methods to decompose full velocity twitch trains from velocity images, highlighting the need for more advanced models and methods for US decomposition than those currently employed.

A non-volitional skeletal muscle endurance test measures functional changes associated with impaired blood flow

Introduction: An electrically stimulated intermittent fatigue test using mechanomyography was recently proposed as a possible tool for detecting clinically relevant changes in muscle function. This study was designed to determine whether the proposed test can detect additional fatigue when it should be present. Methods: Subjects (n = 10) underwent two trials each (occluded and normal blood flow) with a standardized fatigue protocol on the Ankle Dorsiflexors (AD) and Wrist Extensors (WE) using a clinical electrical stimulator. Results: Mean normalized twitch acceleration was strongly predictive of mean normalized torque (R ² = 0.828). The WE experienced lower twitch magnitudes throughout the tourniquet trial (10.81 ± 1.25 m/s²) compared to normal blood flow (18.05 ± 1.06 m/s²). The AD twitches were overall reduced in the tourniquet trial (3.87 ± 0.48 m/s²) compared with the control trial (8.57 ± 0.91 m/s²). Conclusion: Occluding blood flow to a muscle should cause greater muscle fatigue. The ability to detect reduced contraction magnitudes during an electrically stimulated fatigue protocol resulting from low blood flow suggests the proposed test may be capable of detecting clinically relevant muscle deficits.

Estimation of contractile parameters of successive twitches in unfused tetanic contractions of single motor units - A proof-of-concept study using ultrafast ultrasound imaging in vivo

During a voluntary contraction, motor units (MUs) fire a train of action potentials, causing summation of the twitch forces, resulting in fused or unfused tetanus. Twitches have been important in studying whole-muscle contractile properties and differentiation between MU types. However, there are still knowledge gaps concerning the voluntary force generation mechanisms. Current methods rely on the spike-triggered averaging technique, which cannot track changes in successive twitches' properties in response to individual neural firings. This study proposes a method that estimates successive twitches contractile parameters of single MUs during low force voluntary isometric contractions in human biceps brachii. We used a previously developed ultrafast ultrasound imaging method to estimate unfused tetanic activity signals of single MUs. A twitch decomposition model was used to decompose unfused tetanic activity signals into individual twitches. This study found that the contractile parameters varied within and across MUs. There was an association between the inter-spike interval and the contraction time (r = 0.49,p < 0.001) and the half-relaxation time (r = 0.58,p < 0.001), respectively. The method shows the proof-of-concept to study MU contractile properties of individual twitches in vivo, which can provide further insights into the force generation mechanisms of voluntary contractions and response to individual neural discharges.

Physical and electrophysiological motor unit characteristics are revealed with simultaneous high-density electromyography and ultrafast ultrasound imaging

Electromyography and ultrasonography provide complementary information about electrophysiological and physical (i.e. anatomical and mechanical) muscle properties. In this study, we propose a method to assess the electrical and physical properties of single motor units (MUs) by combining High-Density surface Electromyography (HDsEMG) and ultrafast ultrasonography (US). Individual MU firings extracted from HDsEMG were used to identify the corresponding region of muscle tissue displacement in US videos. The time evolution of the tissue velocity in the identified region was regarded as the MU tissue displacement velocity. The method was tested in simulated conditions and applied to experimental signals to study the local association between the amplitude distribution of single MU action potentials and the identified displacement area. We were able to identify the location of simulated MUs in the muscle cross-section within a 2 mm error and to reconstruct the simulated MU displacement velocity (cc > 0.85). Multiple regression analysis of 180 experimental MUs detected during isometric contractions of the biceps brachii revealed a significant association between the identified location of MU displacement areas and the centroid of the EMG amplitude distribution. The proposed approach has the potential to enable non-invasive assessment of the electrical, anatomical, and mechanical properties of single MUs in voluntary contractions.

A Wearable All-Gel Multimodal Cutaneous Sensor Enabling Simultaneous Single-Site Monitoring of Cardiac-Related Biophysical Signals

The human cutaneous sensory organ is a highly evolved biosensor that is efficient, sensitive, selective, and adaptable. Recently, with the development of various materials and structures inspired by sensory organs, artificial cutaneous sensors have been widely studied. In this study, the acquisition of biophysical signals is demonstrated at one point on the body using a wearable all-gel-integrated multimodal sensor composed of four element sensors, inspired by the slow/rapid adapting functions of the skin sensory receptors. The gel-type sensors ensure flexibility, compactness, portability, adherence, and integrity. The wearable all-gel multimodal sensor is easily attached to the wrist and simultaneously gathers blood pressure (BP), electrocardiogram (ECG), electromyogram (EMG), and mechanomyogram (MMG) signals related to cardiac and muscle health. Human activity causes muscle contraction, which affects blood flow; therefore, the relationship between the muscle and heart is crucial for screening and predicting heart health. Cardiac health is monitored by obtaining the two types of phase time differences (i.e., Δt_be : BP and ECG, Δt_em : ECG and MMG) generated during muscle movement. The suggested multimodal sensor has potential applicability in monitoring biophysical conditions and diagnosing cardiac-related health problems.

Identification of single motor units in skeletal muscle under low force isometric voluntary contractions using ultrafast ultrasound

The central nervous system (CNS) controls skeletal muscles by the recruitment of motor units (MUs). Understanding MU function is critical in the diagnosis of neuromuscular diseases, exercise physiology and sports, and rehabilitation medicine. Recording and analyzing the MUs’ electrical depolarization is the basis for state-of-the-art methods. Ultrafast ultrasound is a method that has the potential to study MUs because of the electrical depolarizations and consequent mechanical twitches. In this study, we evaluate if single MUs and their mechanical twitches can be identified using ultrafast ultrasound imaging of voluntary contractions. We compared decomposed spatio-temporal components of ultrasound image sequences against the gold standard needle electromyography. We found that 31% of the MUs could be successfully located and their firing pattern extracted. This method allows new non-invasive opportunities to study mechanical properties of MUs and the CNS control in neuromuscular physiology.

Mechanotendography in Achillodynia shows reduced oscillation variability of pre-loaded Achilles tendon: a pilot study

The present study focuses on an innovative approach in measuring the mechanical oscillations of pre-loaded Achilles tendon by using Mechanotendography (MTG) during application of a short yet powerful mechanical pressure impact. This was applied on the forefoot from the plantar side in direction of dorsiflexion, while the subject stood on the ball of the forefoot on one leg. Participants with Achilles tendinopathy (AT; n = 10) were compared to healthy controls (Con; n = 10). Five trials were performed on each side of the body. For evaluation, two intervals after the impulse began (0-100ms; 30-100ms) were cut from the MTG and pressure raw signals. The intrapersonal variability between the five trials in both intervals were evaluated using the arithmetic mean and coefficient of variation of the mean correlation (Spearman rank correlation) and the normalized averaged mean distances, respectively. The AT-group showed a significantly reduced variability in MTG compared to the Con-group (from p = 0.006 to p = 0.028 for different parameters). The 95% confidence intervals (CI) of MTG results were disjoint, whereas the 95% CIs of the pressure signals were similar (p = 0.192 to p = 0.601). We suggest from this work that the variability of mechanical tendon oscillations could be an indicative parameter of an altered Achilles tendon functionality.

Mechanotendography in Achillodynia shows reduced oscillation variability of pre-loaded Achilles tendon: a pilot study

The present study focuses on an innovative approach in measuring the mechanical oscillations of pre-loaded Achilles tendon by using Mechanotendography (MTG) during application of a short yet powerful mechanical pressure impact. This was applied on the forefoot from the plantar side in direction of dorsiflexion, while the subject stood on the ball of the forefoot on one leg. Participants with Achilles tendinopathy (AT; n = 10) were compared to healthy controls (Con; n = 10). Five trials were performed on each side of the body. For evaluation, two intervals after the impulse began (0-100ms; 30-100ms) were cut from the MTG and pressure raw signals. The intrapersonal variability between the five trials in both intervals were evaluated using the arithmetic mean and coefficient of variation of the mean correlation (Spearman rank correlation) and the normalized averaged mean distances, respectively. The AT-group showed a significantly reduced variability in MTG compared to the Con-group (from p = 0.006 to p = 0.028 for different parameters). The 95% confidence intervals (CI) of MTG results were disjoint, whereas the 95% CIs of the pressure signals were similar (p = 0.192 to p = 0.601). We suggest from this work that the variability of mechanical tendon oscillations could be an indicative parameter of an altered Achilles tendon functionality.

Normalized maximal rate of torque development during voluntary and stimulated static contraction in human tibialis anterior: Influence of age

The risk of falling in older adults has been related, among other factors, to the reduction of the rate of torque development (RTD) with age. It is well known that both structural/peripheral and neural factors can influence the RTD. The purpose of this study was to compare the normalized RTD in young and older participants obtained during a) rapid voluntary tension production and b) neuromuscular electrical stimulation. The tibialis anterior of 19 young subjects (10 males and 9 females; age 21-33 years old) and 19 older participants (10 males and 9 females; age 65-80 years old) was studied. The subjects performed a series of maximal isometric explosive dorsiflexions and underwent trains of supra-maximal electrical stimulations (35 Hz) on the tibialis anterior motor point. Muscle shortening was indirectly measured using a laser (surface mechanomyogram, MMG). Both torque and MMG were normalized to their maximum value. Using a 20 ms sliding window on the normalized torque signal, the normalized maximum RTD was calculated for both voluntary and stimulated contractions. Active stiffness of the muscle- tendon unit was calculated as the area of the normalized torque with respect to the normalized MMG. Normalized maximum RTD was found significantly lower in older adults during voluntary activity (young: 751.9 ± 216.3%/s and old: 513.9 ± 173.9%/s; P < .001), and higher during stimulated contractions (young: 753.1 ± 225.9%/s and old: 890.1 ± 221.3%/s; P = .009). Interestingly, active stiffness was also higher in older adults (young: 3524.6 ± 984.6‰ and old 4144.6 ± 816.6‰; P = .041) and significantly correlated to the normalized maximum RTD during stimulated contractions. This dichotomy suggests that modifications in the structural/peripheral muscle properties are not sufficient to counteract the age-related decrease in neural drive to the muscle during voluntary isometric contractions in aged participants.

Assessment of Skeletal Muscle Contractile Properties by Radial Displacement: The Case for Tensiomyography

Skeletal muscle operates as a near-constant volume system; as such muscle shortening during contraction is transversely linked to radial deformation. Therefore, to assess contractile properties of skeletal muscle, radial displacement can be evoked and measured. Mechanomyography measures muscle radial displacement and during the last 20 years, tensiomyography has become the most commonly used and widely reported technique among the various methodologies of mechanomyography. Tensiomyography has been demonstrated to reliably measure peak radial displacement during evoked muscle twitch, as well as muscle twitch speed. A number of parameters can be extracted from the tensiomyography displacement/time curve and the most commonly used and reliable appear to be peak radial displacement and contraction time. The latter has been described as a valid non-invasive means of characterising skeletal muscle, based on fibre-type composition. Over recent years, applications of tensiomyography measurement within sport and exercise have appeared, with applications relating to injury, recovery and performance. Within the present review, we evaluate the perceived strengths and weaknesses of tensiomyography with regard to its efficacy within applied sports medicine settings. We also highlight future tensiomyography areas that require further investigation. Therefore, the purpose of this review is to critically examine the existing evidence surrounding tensiomyography as a tool within the field of sports medicine.

Mechanomyography responses characterize altered muscle function during electrical stimulation-evoked cycling in individuals with spinal cord injury

BACKGROUND

Investigation of muscle fatigue during functional electrical stimulation (FES)-evoked exercise in individuals with spinal cord injury using dynamometry has limited capability to characterize the fatigue state of individual muscles. Mechanomyography has the potential to represent the state of muscle function at the muscle level. This study sought to investigate surface mechanomyographic responses evoked from quadriceps muscles during FES-cycling, and to quantify its changes between pre- and post-fatiguing conditions in individuals with spinal cord injury.

METHODS

Six individuals with chronic motor-complete spinal cord injury performed 30-min of sustained FES-leg cycling exercise on two days to induce muscle fatigue. Each participant performed maximum FES-evoked isometric knee extensions before and after the 30-min cycling to determine pre- and post- extension peak torque concomitant with mechanomyography changes.

FINDINGS

Similar to extension peak torque, normalized root mean squared (RMS) and mean power frequency (MPF) of the mechanomyography signal significantly differed in muscle activities between pre- and post-FES-cycling for each quadriceps muscle (extension peak torque up to 69%; RMS up to 80%, and MPF up to 19%). Mechanomyographic-RMS showed significant reduction during cycling with acceptable between-days consistency (intra-class correlation coefficients, ICC = 0.51-0.91). The normalized MPF showed a weak association with FES-cycling duration (ICC = 0.08-0.23). During FES-cycling, the mechanomyographic-RMS revealed greater fatigue rate for rectus femoris and greater fatigue resistance for vastus medialis in spinal cord injured individuals.

INTERPRETATION

Mechanomyographic-RMS may be a useful tool for examining real time muscle function of specific muscles during FES-evoked cycling in individuals with spinal cord injury.

Investigation of the Relationship Between Electrical Stimulation Frequency and Muscle Frequency Response Under Submaximal Contractions

Neuromuscular electrical stimulation (NMES) is a common tool that is used in clinical and laboratory experiments and can be combined with mechanomyography (MMG) for biofeedback in neuroprostheses. However, it is not clear if the electrical current applied to neuromuscular tissues influences the MMG signal in submaximal contractions. The objective of this study is to investigate whether the electrical stimulation frequency influences the mechanomyographic frequency response of the rectus femoris muscle during submaximal contractions. Thirteen male participants performed three maximal voluntary isometric contractions (MVIC) recorded in isometric conditions to determine the maximal force of knee extensors. This was followed by the application of nine modulated NMES frequencies (20, 25, 30, 35, 40, 45, 50, 75, and 100 Hz) to evoke 5% MVIC. Muscle behavior was monitored by the analysis of MMG signals, which were decomposed into frequency bands by using a Cauchy wavelet transform. For each applied electrical stimulus frequency, the mean MMG spectral/frequency response was estimated for each axis (X, Y, and Z axes) of the MMG sensor with the values of the frequency bands used as weights (weighted mean). Only with respect to the Z (perpendicular) axis of the MMG signal, the stimulus frequency of 20 Hz did not exhibit any difference with the weighted mean (P = 0.666). For the frequencies of 20 and 25 Hz, the MMG signal displayed the bands between 12 and 16 Hz in the three axes (P < 0.050). In the frequencies from 30 to 100 Hz, the muscle presented a higher concentration of the MMG signal between the 22 and 29 Hz bands for the X and Z axes, and between 16 and 34 Hz bands for the Y axis (P < 0.050 for all cases). We observed that MMG signals are not dependent on the applied NMES frequency, because their frequency contents tend to mainly remain between the 20- and 25-Hz bands. Hence, NMES does not interfere with the use of MMG in neuroprosthesis.

Mechanomyographic responses for the biceps brachii are associated with failure times during isometric force tasks

In order to characterize the physiological adjustments within the neuromuscular system that contribute to task failure, this study examined the surface mechanomyographic (MMG) response during maximal and submaximal isometric force tasks of the elbow flexors sustained to failure. The time and frequency components of the MMG signal have shown to be influenced by motor unit activation patterns as well as tetanus. Therefore, it was hypothesized that the rate of change for the MMG response would associate with failure times and would be reduced to a similar degree between the two tasks. The isometric force tasks were performed by the dominant elbow flexors of twenty healthy males (age: 25 ± 4 years) and MMG was collected from the biceps brachii. Regression analyses were used to model the relationships between the rates of change for MMG versus failure times. There were high levels of interindividual variability in the response patterns, yet the models demonstrated significant negative associations between the rate of change for the MMG responses and failure times during both tasks (R2  = 0.41-0.72, P < 0.05). Similarly, the mean MMG amplitude and frequency values were reduced to comparable levels at the failure point of the two tasks. The results of this study demonstrated that force failure is associated with the rate of diminution in the properties of the muscle force twitch.

Comparing electro- and mechano-myographic muscle activation patterns in self-paced pediatric gait

Electromyography (EMG) is the standard modality for measuring muscle activity. However, the convenience and availability of low-cost accelerometer-based wearables makes mechanomyography (MMG) an increasingly attractive alternative modality for clinical applications. Literature to date has demonstrated a strong association between EMG and MMG temporal alignment in isometric and isokinetic contractions. However, the EMG-MMG relationship has not been studied in gait. In this study, the concurrence of EMG- and MMG-detected contractions in the tibialis anterior, lateral gastrocnemius, vastus lateralis, and biceps femoris muscles were investigated in children during self-paced gait. Furthermore, the distribution of signal power over the gait cycle was statistically compared between EMG-MMG modalities. With EMG as the reference, muscular contractions were detected based on MMG with balanced accuracies between 88 and 94% for all muscles except the gastrocnemius. MMG signal power differed from that of EMG during certain phases of the gait cycle in all muscles except the biceps femoris. These timing and power distribution differences between the two modalities may in part be related to muscle fascicle length changes that are unique to muscle motion during gait. Our findings suggest that the relationship between EMG and MMG appears to be more complex during gait than in isometric and isokinetic contractions.

Age and Sex Effects on the Active Stiffness of Vastus Intermedius under Isometric Contraction

Previously, a novel technique was proposed to quantify the relationship between the muscle stiffness and its nonfatigue contraction intensity. The method extended the measured range of isometric contraction to 100% maximum voluntary contraction (MVC) using an ultrasonic shear wave measurement setup. Yet, it has not been revealed how this relationship could be affected by factors like age or sex. To clarify these questions, vastus intermedius (VI) stiffness of 40 healthy subjects was assessed under 11 step levels of isometric contraction. The subjects were divided into four groups: young males, young females, elderly males, and elderly females (n = 10 for each). In a relaxed state, no significant difference was observed between the male and female subjects (p = 0.156) nor between the young and elderly subjects (p = 0.221). However, when performing isometric contraction, the VI stiffness of males was found to be significantly higher than that of females at the same level (p < 0.001), and that of the young was higher than the elderly (p < 0.001). Meanwhile, for two knee joint angles used, the stiffness measured at a 90° knee joint angle was always significantly larger than that measured at 60° (p < 0.001). Recognizing the active muscle stiffness of VI contributes to body stability, and these results may provide insight into the age and sex bias in musculoskeletal studies, such as those on fall risks.

Estimating the minimum stimulation frequency necessary to evoke tetanic progression based on muscle twitch parameters

The summation of the muscle force caused by an increase in the firing rate is named a tetanic contraction (tetanus), and the minimum stimulation frequency necessary to evoke an unfused/fused tetanus is related to the contraction time (CT) and relaxation time (RT) of the twitch. In particular, the fusion index (FI) is a very useful indicator, and it is used to evaluate the change in the muscle fiber component ratio. However, the measurement of the FI is invasive, because most patients experience pain during the electrical stimulation for tetanus. We expect that the twitch parameters CT and RT can substitute for the FI in the future. We found that the minimum stimulation frequency necessary to evoke the unfused/fused tetanus can be estimated from the twitch parameters as a first step. The results showed that (1) the minimum stimulation frequencies calculated from twitch parameters during unfused/fused tetanus were very similar to those calculated from FI parameters, and (2) they were also strongly correlated with FI parameters regardless of fiber components. The basic characteristics of tetanic progression in different fiber types could be estimated from twitch parameters.

Torque and mechanomyogram relationships during electrically-evoked isometric quadriceps contractions in persons with spinal cord injury

The interaction between muscle contractions and joint loading produces torques necessary for movements during activities of daily living. However, during neuromuscular electrical stimulation (NMES)-evoked contractions in persons with spinal cord injury (SCI), a simple and reliable proxy of torque at the muscle level has been minimally investigated. Thus, the purpose of this study was to investigate the relationships between muscle mechanomyographic (MMG) characteristics and NMES-evoked isometric quadriceps torques in persons with motor complete SCI. Six SCI participants with lesion levels below C4 [(mean (SD) age, 39.2 (7.9) year; stature, 1.71 (0.05) m; and body mass, 69.3 (12.9) kg)] performed randomly ordered NMES-evoked isometric leg muscle contractions at 30°, 60° and 90° knee flexion angles on an isokinetic dynamometer. MMG signals were detected by an accelerometer-based vibromyographic sensor placed over the belly of rectus femoris muscle. The relationship between MMG root mean square (MMG-RMS) and NMES-evoked torque revealed a very high association (R(2)=0.91 at 30°; R(2)=0.98 at 60°; and R(2)=0.97 at 90° knee angles; P<0.001). MMG peak-to-peak (MMG-PTP) and stimulation intensity were less well related (R(2)=0.63 at 30°; R(2)=0.67 at 60°; and R(2)=0.45 at 90° knee angles), although were still significantly associated (P≤0.006). Test-retest interclass correlation coefficients (ICC) for the dependent variables ranged from 0.82 to 0.97 for NMES-evoked torque, between 0.65 and 0.79 for MMG-RMS, and from 0.67 to 0.73 for MMG-PTP. Their standard error of measurements (SEM) ranged between 10.1% and 31.6% (of mean values) for torque, MMG-RMS and MMG-PTP. The MMG peak frequency (MMG-PF) of 30Hz approximated the stimulation frequency, indicating NMES-evoked motor unit firing rate. The results demonstrated knee angle differences in the MMG-RMS versus NMES-isometric torque relationship, but a similar torque related pattern for MMG-PF. These findings suggested that MMG was well associated with torque production, reliably tracking the motor unit recruitment pattern during NMES-evoked muscle contractions. The strong positive relationship between MMG signal and NMES-evoked torque production suggested that the MMG might be deployed as a direct proxy for muscle torque or fatigue measurement during leg exercise and functional movements in the SCI population.

The Skin Acts to Maintain Muscle Shear Modulus

It is not clear how the tissues covering the skeletal muscles affect the muscles' mechanical properties. The main purpose of this study was to examine changes in muscle shear modulus as a representative mechanical property of muscle with and without the covering tissues of skin and epimysium (fascia). Shear modulus of the medial gastrocnemius (MG) muscle was determined using ultrasound shear-wave elastography in the Thiel's embalmed cadavers under three different conditions: original (intact cadavers), removal of the skin on the MG and subsequent removal of the epimysium. Muscle shear modulus significantly decreased by 50% after removal of the skin, whereas no additional changes in shear modulus were observed after subsequent removal of the epimysium. This study suggests that the skin is a main contributor for maintaining the muscle mechanical properties among tissues covering the skeletal muscle.

Mechanomyographic parameter extraction methods: an appraisal for clinical applications

The research conducted in the last three decades has collectively demonstrated that the skeletal muscle performance can be alternatively assessed by mechanomyographic signal (MMG) parameters. Indices of muscle performance, not limited to force, power, work, endurance and the related physiological processes underlying muscle activities during contraction have been evaluated in the light of the signal features. As a non-stationary signal that reflects several distinctive patterns of muscle actions, the illustrations obtained from the literature support the reliability of MMG in the analysis of muscles under voluntary and stimulus evoked contractions. An appraisal of the standard practice including the measurement theories of the methods used to extract parameters of the signal is vital to the application of the signal during experimental and clinical practices, especially in areas where electromyograms are contraindicated or have limited application. As we highlight the underpinning technical guidelines and domains where each method is well-suited, the limitations of the methods are also presented to position the state of the art in MMG parameters extraction, thus providing the theoretical framework for improvement on the current practices to widen the opportunity for new insights and discoveries. Since the signal modality has not been widely deployed due partly to the limited information extractable from the signals when compared with other classical techniques used to assess muscle performance, this survey is particularly relevant to the projected future of MMG applications in the realm of musculoskeletal assessments and in the real time detection of muscle activity.

Mechanomyography and muscle function assessment: a review of current state and prospects

Previous studies have explored to saturation the efficacy of the conventional signal (such as electromyogram) for muscle function assessment and found its clinical impact limited. Increasing demand for reliable muscle function assessment modalities continues to prompt further investigation into other complementary alternatives. Application of mechanomyographic signal to quantify muscle performance has been proposed due to its inherent mechanical nature and ability to assess muscle function non-invasively while preserving muscular neurophysiologic information. Mechanomyogram is gaining accelerated applications in evaluating the properties of muscle under voluntary and evoked muscle contraction with prospects in clinical practices. As a complementary modality and the mechanical counterpart to electromyogram; mechanomyogram has gained significant acceptance in analysis of isometric and dynamic muscle actions. Substantial studies have also documented the effectiveness of mechanomyographic signal to assess muscle performance but none involved comprehensive appraisal of the state of the art applications with highlights on the future prospect and potential integration into the clinical practices. Motivated by the dearth of such critical review, we assessed the literature to investigate its principle of acquisition, current applications, challenges and future directions. Based on our findings, the importance of rigorous scientific and clinical validation of the signal is highlighted. It is also evident that as a robust complement to electromyogram, mechanomyographic signal may possess unprecedented potentials and further investigation will be enlightening.

Laser Doppler vibrometry measurement of the mechanical myogram

Contracting muscles show complex dimensional changes that include lateral expansion. Because this expansion process is intrinsically vibrational, driven by repetitive actions of multiple motor units, it can be sensed and quantified using the method of Laser Doppler Vibrometry (LDV). LDV has a number of advantages over more traditional mechanical methods based on microphones and accelerometers. The LDV mechanical myogram from a small hand muscle (the first dorsal interosseous) was studied under conditions of elastic loading applied to the tip of the abducted index finger. The LDV signal was shown to be related systematically to the level of force production, and to compare favorably with conventional methods for sensing the mechanical and electrical aspects of muscle contraction.

MUSCLES WITHIN MUSCLES: A MECHANOMYOGRAPHIC ANALYSIS OF MUSCLE SEGMENT CONTRACTILE PROPERTIES WITHIN HUMAN GLUTEUS MAXIMUS

The aim of this investigation was to determine the contractile properties of motor units within 3 segments of the gluteus maximus utilizing a laser-based mechanomyographic (MMG) technique. The intention was to determine whether there were segmental differences in motor unit contractile properties and whether these differences may be related to the muscle segment's function and its fibre type composition.

Ten subjects were recruited from the student population at the University of Wollongong. Maximal percutaneous neuromuscular stimulation (PNS) was delivered to the medial and lateral portions of three (cranial, middle, caudal) muscle segments of the gluteus maximus by an MMG stimulator. An MMG laser sensor measured the lateral displacement of the muscle segment belly resulting from the development of maximal isometric tension. Parameters characterizing the MMG waveforms were statistically compared to determine variations in contractile properties both within (medial to lateral) and between segments.

Our results indicated that the contractile properties of motor units varied significantly (p < 0.05) between, but not within (medial to lateral), the three segments of the gluteus maximus. Most the gluteus maximus. Most notably, segment contraction time (tc) decreased significantly (p < 0.05) in a cranio to caudal direction suggesting a variation in muscle fibre type composition within the three segments of the muscle. Even when corrected for differences in muscle belly displacement between subjects, the cranial segment was found to have a significantly (p < 0.05) longer contraction time than the two more caudal segments. The results suggest that the gluteus maximus was composed of muscle segments that were physiologically, as well as anatomically, designed to fulfil particular roles during everyday motor tasks. Based upon these results, the MMG technique appears to have considerable utility for the non-invasive assessment of muscle segment contractile properties for both laboratory and clinical applications.

Mechanomyographic response during FES in healthy and paraplegic subjects

Mechanomyography (MMG) registers lateral oscillations of contracting muscles. Functional electrical stimulation (FES) improves the rehabilitation of paraplegic subjects and can be used in neuroprosthesis control. During FES application, muscular contraction responses may vary, possibly due to fatigue or adaptation of nerve cells face to electrical stimuli. This study measured the differences in MMG RMS and median frequency (MF) features between healthy (HV) and spinal cord injury (SCI) volunteers. Ten HV and three SCI participated in the research. FES waveform consisted of a monophasic square wave, 1kHz pulse frequency, 100us active pulse period and 3ms active burst period with burst frequency of 70Hz. For each stimulation series, three analysis windows were inspected. RMS and MF variations were inversely related. The obtained results may help to create new strategies of muscular closed-loop control.

Spectral properties of electromyographic and mechanomyographic signals during isometric ramp and step contractions in biceps brachii

The purposes of this study were: (1) to apply wavelet and principal component analysis to quantify the spectral properties of the surface EMG and MMG signals from biceps brachii during isometric ramp and step muscle contractions when the motor units are recruited in an orderly manner, and (2) to compare the recruitment patterns of motor unit during isometric ramp and step muscle contractions. Twenty healthy participants (age = 34 ± 10.7 years) performed step and ramped isometric contractions. Surface EMG and MMG were recorded from biceps brachii. 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 (PCA) and ANCOVA. Wavelet combined PCA offers a quantitative measure of the contribution of high and low frequency content within the EMG and MMG. The ANCOVA indicated that there was no significant difference in EMG total intensity, EMG(MPF), first and second principal component loading scores (PCI and PCII) between ramp and step contractions, whereas the MMG(MPF) and MMG PCI loading scores were significantly higher during ramp contractions than during step contractions. These findings suggested that EMG and MMG may offer complimentary information regarding the interactions between motor unit recruitment and firing rate that control muscle force production. In addition, our results support the hypothesis that different motor unit recruitment strategy was used by the muscle when contracting under different conditions.

The design and testing of a novel mechanomyogram-driven switch controlled by small eyebrow movements

Background

Individuals with severe physical disabilities and minimal motor behaviour may be unable to use conventional mechanical switches for access. These persons may benefit from access technologies that harness the volitional activity of muscles. In this study, we describe the design and demonstrate the performance of a binary switch controlled by mechanomyogram (MMG) signals recorded from the frontalis muscle during eyebrow movements.

Methods

Muscle contractions, detected in real-time with a continuous wavelet transform algorithm, were used to control a binary switch for computer access. The automatic selection of scale-specific thresholds reduced the effect of artefact, such as eye blinks and head movement, on the performance of the switch. Switch performance was estimated by cued response-tests performed by eleven participants (one with severe physical disabilities).

Results

The average sensitivity and specificity of the switch was 99.7 ± 0.4% and 99.9 ± 0.1%, respectively. The algorithm performance was robust against typical participant movement.

Conclusions

The results suggest that the frontalis muscle is a suitable site for controlling the MMG-driven switch. The high accuracies combined with the minimal requisite effort and training show that MMG is a promising binary control signal. Further investigation of the potential benefits of MMG-control for the target population is warranted.

Two dimensional spatial coherence of the natural vibrations of the biceps brachii muscle generated during voluntary contractions

Surface mechanomyograms (S-MMGs) are recorded from low frequency (〈100 Hz) mechanical oscillations that are naturally generated by skeletal muscle during voluntary contractions. This study investigates a method to determine the propagation directionality of the S-MMG waves. A 3×5 grid of skin mounted accelerometers was mounted on the biceps brachii muscle during submaximal voluntary contractions. This method resulted in findings that the propagation directionality of the S-MMGs are frequency dependent. At high frequencies (>25 Hz), high spatial coherence values were only measured for sensor pairs aligned along the proximal to distal (i.e. longitudinal) orientation, thus indicating that coherent S-MMG were mainly propagating along the muscle fibers direction of the biceps brachii at those frequencies. On the other hand, at lower frequencies (〈25 Hz), the S-MMG spatial coherence values did not exhibit a specific directionality. This method provides results that have an important implication of finding the average phase velocity of the propagating S-MMG wave, which can be used to determine viscoelastic properties of skeletal muscles.

Mechanomyogram for identifying muscle activity and fatigue

Mechanomyogram is the recording of the acoustic activity associated with the muscle contraction. While discovered nearly a decade ago with the intention of providing an alternate to the surface electromyogram, it has not yet been investigated thoroughly and there are no current applications associated with MMG. This paper reports an experimental study of MMG against force of contraction and muscle fatigue during cyclic contraction. The results indicate that there is a relationship between the intensity of the MMG recording and force of contraction. A change in the intensity of MMG is also observed with the onset of muscle fatigue. However, the inter-subject variation is very large. The results also indicate that the spectrum of the MMG is very inconsistent and not a useful feature of the signal.

Use of MMG signals for the control of powered orthotic devices: development of a rectus femoris measurement protocol

A test protocol is defined for the purpose of measuring rectus femoris mechanomyographic (MMG) signals. The protocol is specified in terms of the following: measurement equipment, signal processing requirements, human postural requirements, test rig, sensor placement, sensor dermal fixation, and test procedure. Preliminary tests of the statistical nature of rectus femoris MMG signals were performed, and Gaussianity was evaluated by means of a two-sided Kolmogorov-Smirnov test. For all 100 MMG data sets obtained from the testing of two volunteers, the null hypothesis of Gaussianity was rejected at the 1%, 5%, and 10% significance levels. Most skewness values were found to be greater than 0.0, while all kurtosis values were found to be greater than 3.0. A statistical convergence analysis also performed on the same 100 MMG data sets suggested that 25 MMG acquisitions should prove sufficient to statistically characterize rectus femoris MMG. This conclusion is supported by the qualitative characteristics of the mean rectus femoris MMG power spectral densities obtained using 25 averages.

Assessment of the mechanical properties of the musculoskeletal system using 2-D and 3-D very high frame rate ultrasound

One of the great challenges for understanding muscular diseases is to assess noninvasively the active and passive mechanical properties of the musculoskeletal system. In this paper we report the use of ultrafast ultrasound imaging to explore with a submillimeter resolution the behavior of the contracting tissues in vivo (biceps brachii). To image the contraction, which is a very brief phenomenon (100 ms), a recently designed ultrasound scanner prototype able to take up to 6000 frames/s was used. A very high frame rate from 1000 to 2500 frames/s was used to image the cross section plane of the muscle (transverse to fibers) enabling us to catch in real time the muscle contraction during a transient electrostimulation. Tissue velocities were obtained from radiofrequency based speckle tracking techniques and their profiles are discussed with respect to electrostimulation intensities and pulse repetition frequencies for different volunteers. Three-dimensional (3-D) very high frame rate movies were also acquired by repeating the experiment for different acquisition planes while triggering the imaging system with the electrostimulation device. The reconstructed 3-D velocity field allows the full localization of the contracting fibers bundle. This ultrasound technique, referred to as echo mechanomyography, offers new perspectives for in vivo and in situ noninvasive muscle diagnosis of an active contractile tissue.

Laser-detected lateral muscle displacement is correlated with force fluctuations during voluntary contractions in humans

Fluctuations in muscle force during steady voluntary contractions result from the summation of twitch forces produced by asynchronous activation of multiple motor units. We hypothesized that oscillatory lateral muscle displacement, measured with a non-contact high-resolution laser displacement sensor, is correlated with force fluctuations during steady, voluntary contractions with a human muscle. Eight healthy young adults (20-33 yrs) performed steady isometric contractions with the first dorsal interosseus muscle. Contraction intensity ranged from 2.5% to 60% of the maximal voluntary contraction force. Oscillatory lateral displacement of the muscle surface was measured with a high-resolution laser displacement sensor (0.5 microm resolution), concurrently with abduction force of the index finger. In the time-domain analysis, there was a significant positive peak in the cross-correlation function between lateral muscle displacement and force fluctuations. In addition, the amplitude increased linearly with contraction intensity in both signals. In the frequency-domain analysis, frequency content was similar in both signals, and there was significant coherence between signals for the major frequency range of the signals (<5 Hz). In conclusion, laser-detected lateral displacement of a hand muscle is correlated with force fluctuations across a wide range of contraction intensity during steady voluntary contractions in humans.

Efeitos de dois parâmetros antropométricos no comportamento do sinal mecanomiográfico em testes de força muscular

Apesar da mecanomiografia (MMG) ser uma técnica com crescente destaque em investigações relativas ao fenômeno da contração muscular, poucos trabalhos se dedicaram a entender os possíveis efeitos de variáveis antropométricas no sinal de MMG. Portanto, o objetivo do estudo foi avaliar os efeitos da dobra cutânea e do perímetro do braço no comportamento do sinal de MMG nos domínios do tempo e da freqüência em diferentes níveis de contração isométrica. Participaram do estudo 21 indivíduos do sexo masculino (24,9 ± 7,8 anos) e 21 do feminino (20,7 ± 2,5 anos), todos destros. O protocolo experimental constou de um teste de contração voluntária máxima (CVM) que permitiu determinar cinco cargas percentuais administradas durante os testes de força (20%, 40%, 60%, 80% e 100% da CVM). Um acelerômetro biaxial foi colocado sobre o ventre muscular do bíceps braquial direito. O sinal de MMG foi analisado nos domínios do tempo, mediante o comportamento da amplitude do sinal, por meio do cálculo da raiz média quadrática (valor RMS), e da freqüência, através da freqüência mediana (FM) calculada a partir do seu espectro de potência. Estes parâmetros foram extraídos a partir dos sinais gerados na direção X, perpendicular às fibras musculares. Foram também medidos os parâmetros antropométricos dobra cutânea biciptal (BD_DC_B) e perímetros de braço relaxado (BD_PR) e contraído (BD_PC) no terço médio do braço direito. Ambos os grupos apresentaram um comportamento decrescente da FM com o nível de contração. O valor RMS apresentou comportamento crescente com a carga para ambos os grupos. Houve diferenças estatisticamente significativas (p < 0,05) entre as variáveis antropométricas quando comparados os grupos. Sugere-se que o valor RMS, ao contrário da FM, que poderá sofrer atenuações por parte dos tecidos presentes na interface entre o músculo e o transdutor, seja mais consistente na discriminação da força muscular.

Fatigue and the electromechanical efficiency of the vastus medialis and vastus lateralis muscles

CONTEXT

The relationship between the amplitudes of the mechanomyographic (MMG) and electromyographic (EMG) signals has been used to examine the "electromechanical efficiency" (EME) of normal and diseased muscle. The EME may help us to better understand the neuromuscular relationship between the vastus medialis and vastus lateralis muscles.

OBJECTIVE

To examine the EME of the vastus medialis and vastus lateralis muscles during a fatiguing task.

DESIGN

Repeated-measures design.

SETTING

Research laboratory.

PATIENTS OR OTHER PARTICIPANTS

Ten healthy males (age = 23.2 +/- 1.2 years) with no history of knee injury.

INTERVENTION(S)

Seventy-five consecutive, maximal concentric isokinetic leg extensions at a velocity of 180 degrees /s.

MAIN OUTCOMES MEASURE(S)

Bipolar surface EMG electrodes were placed over the vastus medialis and vastus lateralis muscles, with an MMG contact sensor placed adjacent to the superior EMG electrode on each muscle. The MMG and EMG amplitude values (root mean squares) were calculated for each of the 75 repetitions and normalized to the highest value from the 75 repetitions. The EME was expressed as the ratio of the log-transformed normalized MMG amplitude to the normalized EMG amplitude. For each muscle, the linear relationship for the normalized-group mean EME was determined across the 75 repetitions.

RESULTS

Linear regression indicated decreases in torque (R(2) = .96), vastus medialis EME (R(2) = .73), and vastus lateralis EME (R(2) = .73). The slopes for the vastus medialis and vastus lateralis EME were not different (P > .10).

CONCLUSIONS

The similarities in the fatigue-induced decreases in EME for the vastus medialis and vastus lateralis muscles suggested that symmetry was present between the muscles in the electric and mechanical responses to repeated, maximal muscle actions. The EME measurements may provide a unique insight into the influence of fatigue on the contractile properties of skeletal muscle, including alterations that occur to the intrinsic electric and mechanical components. The EME may be useful in assessing and quantifying clinically relevant asymmetries in vastus medialis and vastus lateralis muscle function in those with knee injuries.

In vivo behaviour of human muscle architecture and mechanomyographic response using the interpolated twitch technique

This study investigated the origin of curvilinear change in the superimposed mechanomyogram (MMG) amplitude of the human medial gastrocnemius muscle (MG) with increasing contraction intensity. The superimposed twitch amplitude, the superimposed MMG amplitude and the extent of fascicle shortening were measured using ultrasonic images of electrical stimulation during isometric plantar flexions at levels 20%, 40%, 60%, 80%, and 100% of the maximal voluntary contraction (MVC). The superimposed twitch amplitude, the superimposed MMG amplitude and the extent of fascicle shortening decreased with increasing contraction intensity. The superimposed MMG amplitude and the extent of fascicle shortening showed a curvilinear decrease, while the superimposed twitch amplitude showed a linear decrease at levels up to 80% of the MVC. There was a linear relationship between the superimposed MMG amplitude and the extent of fascicle shortening at different contraction intensities. These results indicate that the superimposed MMG amplitude reflects changes in the extent of fascicle shortening at different contraction intensities better than the superimposed twitch amplitude. Our study suggests that the origin of the curvilinear decrease of superimposed MMG amplitude is associated with a curvilinear decrease of the extent of fascicle shortening with increasing contraction intensity in the human MG.

In vivo myograph measurement of muscle contraction at optimal length

BACKGROUND

Current devices for measuring muscle contraction in vivo have limited accuracy in establishing and re-establishing the optimum muscle length. They are variable in the reproducibility to determine the muscle contraction at this length, and often do not maintain precise conditions during the examination. Consequently, for clinical testing only semi-quantitative methods have been used.

METHODS

We present a newly developed myograph, an accurate measuring device for muscle contraction, consisting of three elements. Firstly, an element for adjusting the axle of the device and the physiological axis of muscle contraction; secondly, an element to accurately position and reposition the extremity of the muscle; and thirdly, an element for the progressive pre-stretching and isometric locking of the target muscle. Thus it is possible to examine individual in vivo muscles in every pre-stretched, specified position, to maintain constant muscle-length conditions, and to accurately re-establish the conditions of the measurement process at later sessions.

RESULTS

In a sequence of experiments the force of contraction of the muscle at differing stretching lengths were recorded and the forces determined. The optimum muscle length for maximal force of contraction was established. In a following sequence of experiments with smaller graduations around this optimal stretching length an increasingly accurate optimum muscle length for maximal force of contraction was determined. This optimum length was also accurately re-established at later sessions.

CONCLUSION

We have introduced a new technical solution for valid, reproducible in vivo force measurements on every possible point of the stretching curve. Thus it should be possible to study the muscle contraction in vivo to the same level of accuracy as is achieved in tests with in vitro organ preparations.

Two-dimensional spatial distribution of surface mechanomyographical response to single motor unit activity

In order to better understand the mechanisms of generation of mechanomyography (MMG) signals, the two-dimensional distribution of surface MMG produced by the activity of single motor units was analyzed by a novel two-dimensional recording method. Motor unit action potentials were identified from intramuscular electromyographic (EMG) signals and used to trigger the averaging of MMG signals detected over the tibialis anterior muscle of 11 volunteers with a grid of 5x3 accelerometers (20-mm inter-accelerometer distance). The intramuscular wires were inserted between the first and second accelerometer in the middle column of the grid, proximal to the innervation zone. The subjects performed three contractions with visual feedback of the intramuscular EMG signals. In each contraction, a new motor unit was recruited at the minimum stable discharge rate (mean+/-S.D., N = 11 subjects, 7.3+/-2.3 pulse/s), resulting in torque of 2.4+/-2.8% of the maximal voluntary contraction (MVC), 4.6+/-2.7% MVC, and 6.3+/-3.1% MVC (all different, P < 0.01). For 23 out of 33 detected motor units, it was possible to extract the motor unit surface acceleration map (MUAM). A negative MUAM peak (-2.7+/-2.2 mm/s2) was detected laterally and a positive MUAM peak (4.1+/-2.4 mm/s2) medially (P < 0.001). The time-to-peak was shorter in the medial part of the muscle (2.9+/-0.4 ms) than in the other locations (3.4+/-0.5 ms, P < 0.001). The double integrated signals (muscle displacement) indicated negative deflection in the lateral part and inflation close to the tibia bone. The maps of acceleration showed spatial dependency in single motor unit MMG activities. The technique provides a new insight into motor unit contractile properties.

Modulation of muscle activity and force fluctuations in the plantarflexors after bedrest depends on knee position

Force fluctuations in leg muscles increase after bedrest, perhaps due to modulation of the neural strategy that is specific to a muscle or common to agonist muscles. The purpose of this study was to examine the modulation of muscle activity and force fluctuations during steady contractions with variable involvement of plantarflexor muscles after bedrest at knee-flexed (FLX) and extended (EXT) positions. Before and after 20-day bedrest, plantarflexion force and surface electromyogram (EMG) in the medial gastrocnemius (MG), lateral gastrocnemius, and soleus muscles were measured during steady isometric contractions in five young men. In EXT, power <or=10 HZ in the rectified EMG of MG increased significantly after bedrest. This low-frequency modulation of muscle activity in MG accompanied a 29% increase in the standard deviation of force. There was no change in EMG in other muscles. In FLX, there was no adjustment in EMG or force fluctuations. These results suggest that low-frequency modulation of MG plays a role in increasing force fluctuations during steady plantarflexion in EXT after bedrest. The findings indicate task/muscle specificity in the modulation of neural strategy during steady contractions after bedrest and underscore the importance of designing a specific training regimen targeted to particular tasks/muscles with regard to force fluctuations in multiple-agonist systems.

Spatial and force dependency of mechanomyographic signal features

The aim was to investigate with a novel technique the spatial inhomogeneity in surface mechanomyographic (MMG) response to muscle contraction at varying force levels. MMG signals were detected over the dominant tibialis anterior muscle of 10 volunteers using a 5 x 3 grid of accelerometers. The subjects performed 3 s long isometric contractions at forces ranging from 0% to 100% of the maximal force (10% increment) in a randomised order. From the two-dimensional MMG recordings, maps of absolute and normalized temporal and spectral MMG descriptors were obtained. The centroid and entropy of these maps were computed to describe the spatial centre of activity and degree of homogeneity, respectively. MMG absolute amplitude did not depend on location over the muscle while normalized amplitude did and the centroid shifted with increasing force. Amplitude increased with force and its entropy decreased. Absolute and normalized spectral variables depended on location over the muscle and their centroid shifted with increasing force. In addition, the dependency of absolute and normalized spectral variables on force was affected by location. These results highlight limitations when using single-channel MMG features for the assessment of motor unit control strategies, due to a substantial effect of position on the relation between force and MMG characteristics.

Mechanomyographic and electromyographic responses to repeated concentric muscle actions of the quadriceps femoris

In comparison to isometric muscle action models, little is known about the electromyographic (EMG) and mechanomyographic (MMG) amplitude and mean power frequency (MPF) responses to fatiguing dynamic muscle actions. Simultaneous examination of the EMG and MMG amplitude and MPF may provide additional insight with regard to the motor control strategies utilized by the superficial muscles of the quadriceps femoris during a concentric fatiguing task. Thus, the purpose of this study was to examine the EMG and MMG amplitude and MPF responses of the vastus lateralis (VL), rectus femoris (RF), and vastus medialis (VM) during repeated, concentric muscle actions of the dominant leg. Seventeen adults (21.8+/-1.7 yr) performed 50 consecutive, maximal concentric muscle actions of the dominant leg extensors on a Biodex System 3 Dynamometer at velocities of 60 degrees s(-1) and 300 degrees s(-1). Bipolar surface electrode arrangements were placed over the mid portion of the VL, RF, and VM muscles with a MMG contact sensor placed adjacent to the superior EMG electrode on each muscle. Torque, MMG and EMG amplitude and MPF values were calculated for each of the 50 repetitions. All values were normalized to the value recorded during the first repetition and then averaged across all subjects. The cubic decreases in torque at 60 degrees s(-1) (R2 = 0.972) and 300 degrees s(-1) (R2 = 0.931) was associated with a decline in torque of 59+/-24% and 53+/-11%, respectively. The muscle and velocity specific responses for the MMG amplitude and MPF demonstrated that each of the superficial muscles of the quadriceps femoris uniquely contributed to the control of force output across the 50 repetitions. These results suggested that the MMG responses for the VL, RF, VM during a fatiguing task may be influenced by a number of factors such as fiber type differences, alterations in activation strategy including motor unit recruitment and firing rate and possibly muscle wisdom.