Experimental and modeling investigation of spectral compression of biceps brachii SEMG activity with increasing force levels

Teaching Essential EMG Theory to Kinesiologists and Physical Therapists Using Analogies Visual Descriptions, and Qualitative Analysis of Biophysical Concepts

Electromyography (EMG) is a multidisciplinary field that brings together allied health (kinesiology and physical therapy) and the engineering sciences (biomedical and electrical). Since the physical sciences are used in the measurement of a biological process, the presentation of the theoretical foundations of EMG is most conveniently conducted using math and physics. However, given the multidisciplinary nature of EMG, a course will most likely include students from diverse backgrounds, with varying levels of math and physics. This is a pedagogical paper that outlines an approach for teaching foundational concepts in EMG to kinesiologists and physical therapists that uses a combination of analogies, visual descriptions, and qualitative analysis of biophysical concepts to develop an intuitive understanding for those who are new to surface EMG. The approach focuses on muscle fiber action potentials (MFAPs), motor unit action potentials (MUAPs), and compound muscle action potentials (CMAPs) because changes in these waveforms are much easier to identify and describe in comparison to the surface EMG interference pattern (IP).

Muscle Fiber Conduction Velocity Correlates With the Age at Onset in Mild FSHD Cases

A majority of patients with facioscapulohumeral muscular dystrophy (FSHD) report severe fatigue. The aim of this study was to explore whether fatigability during a performance task is related to the main clinical features of the disease in mildly affected patients. A total of 19 individuals with a molecular genetic-based diagnosis of FSHD (median D4Z4 deletion length of 27 kb) performed two isometric flexions of the dominant biceps brachii at 20% of their maximal voluntary contraction (MVC) for 2 min, and then at 60% MVC until exhaustion. Fatigability indices (average rectified value, mean frequency, conduction velocity, and fractal dimension) were extracted from the surface electromyogram (sEMG) signal, and their correlations with age, age at onset, disease duration, D4Z4 contraction length, perceived fatigability, and clinical disability score were analyzed. The conduction velocity during the low level contraction showed a significant negative correlation with the age at onset (p < 0.05). This finding suggest the assessment of conduction velocity at low isometric contraction intensities, as a potential useful tool to highlight differences in muscle involvement in FSHD patients.

Increased resistance towards fatigability in patients with facioscapulohumeral muscular dystrophy

PURPOSE

In facioscapulohumeral muscular dystrophy (FSHD) fatigue is a major complaint. We aimed to investigate whether during isometric sustained elbow flexions, performance fatigability indexes differ in patients with FSHD with respect to healthy controls.

METHODS

Seventeen patients with FSHD and seventeen healthy controls performed two isometric flexions of the dominant biceps brachii at 20% of their maximal voluntary contraction (MVC) for 2 min and then at 60% MVC until exhaustion. Muscle weakness was characterized as a percentage of predicted values. Maximal voluntary strength, endurance time and performance fatigability indices (mean frequency of the power spectrum (MNF), muscle fiber conduction velocity (CV) and fractal dimension (FD)), extracted from the surface electromyogram signal (sEMG) were compared between the two groups.

RESULTS

In patients with FSHD, maximal voluntary strength was 68.7% of predicted value (p < 0.01). Compared to healthy controls, FSHD patients showed reduced MVC (p < 0.001; r = 0.62) and lower levels of performance fatigability, characterized by reduced rate of changes in MNF (p < 0.01; r = 0.56), CV (p < 0.05; 0.37) and FD (p < 0.001; r = 0.51) and increased endurance time (p < 0.001; r = 0.63), during the isometric contraction at 60% MVC.

CONCLUSION

A decreased reduction in the slopes of all the considered sEMG parameters during sustained isometric elbow flexions suggests that patients with FSHD experience lower levels of performance fatigability compared to healthy controls.

Spatial filtering for enhanced high-density surface electromyographic examination of neuromuscular changes and its application to spinal cord injury

Background

Spatial filtering of multi-channel signals is considered to be an effective pre-processing approach for improving signal-to-noise ratio. The use of spatial filtering for preprocessing high-density (HD) surface electromyogram (sEMG) helps to extract critical spatial information, but its application to non-invasive examination of neuromuscular changes have not been well investigated.

Methods

Aimed at evaluating how spatial filtering can facilitate examination of muscle paralysis, three different spatial filtering methods are presented using principle component analysis (PCA) algorithm, non-negative matrix factorization (NMF) algorithm, and both combination, respectively. Their performance was evaluated in terms of diagnostic power, through HD-sEMG clustering index (CI) analysis of neuromuscular changes in paralyzed muscles following spinal cord injury (SCI).

Results

The experimental results showed that: (1) The CI analysis of conventional single-channel sEMG can reveal complex neuromuscular changes in paralyzed muscles following SCI, and its diagnostic power has been confirmed to be characterized by the variance of Z scores; (2) the diagnostic power was highly dependent on the location of sEMG recording channel. Directly averaging the CI diagnostic indicators over channels just reached a medium level of the diagnostic power; (3) the use of either PCA-based or NMF-based filtering method yielded a greater diagnostic power, and their combination could even enhance the diagnostic power significantly.

Conclusions

This study not only presents an essential preprocessing approach for improving diagnostic power of HD-sEMG, but also helps to develop a standard sEMG preprocessing pipeline, thus promoting its widespread application.

Alterations of Muscle Activation Pattern in Stroke Survivors during Obstacle Crossing

OBJECTIVE

This study investigates changes in the neuromuscular activation pattern of the lower limb muscles in stroke survivors when crossing obstacles of three different heights.

METHODS

Eight stroke survivors and eight age-, height-, and gender-matched healthy controls were recruited and instructed to cross over obstacles with heights of 10, 20, and 30% leg length. Surface electromyography (EMG) signals were recorded from the rectus femoris (RF), biceps femoris (BF), tibialis anterior (TA), and medial gastrocnemius (MG) of both limbs. Muscle activation signals were normalized to maximum voluntary contraction. Differences between groups and heights were compared using the root mean square of EMG, the cocontraction index of agonist and antagonist muscles, and power spectral analysis based on the mean power frequency (MPF). The correlations between the calculated variables and clinical scales such as Berg Balance Scale and Fugl-Meyer assessment (FMA) were also examined.

RESULTS

During the leading limb swing phase, the activation levels of all four muscles were greater in the stroke group than the healthy controls (p < 0.05), and the TA showed increased activation level with increasing obstacle height in both groups (p < 0.05). Cocontraction between the TA and MG was higher in the stroke group during the swing phase of the leading limb and between the RF and BF during the stance phase (p < 0.05). Similarly, for the trailing limb, increased cocontractions between the two pairs of agonist and antagonist muscles were found during the stance phase in the stroke group (p < 0.05). During the crossing stride, the frequency analysis showed significantly smaller MPF values in all four lower limb muscles in the leading limb of stroke survivors compared with healthy controls (p < 0.05). Moreover, significant correlations were found between the FMA scores and the BF and TA activations in the leading limb during the swing phase (p < 0.05).

CONCLUSION

Greater activation levels of the lower limb muscles resulted in higher muscular demands for stroke survivors, which might lead to greater difficulty in maintaining balance. The increased cocontraction during obstacle crossing might be compensation for the affected stability and enable safe crossing for stroke survivors. The reduced MPF in the affected limb of the stroke group might be due to impairments in motor units or other complex neuromuscular alterations.

The Isokinetic and Electromyographic Assessment of Knee Muscles Strength in the Short- and Long-Term Type 2 Diabetes

BACKGROUND

Type 2 diabetes (T2DM) patients are subject to muscle weakness.

OBJECTIVES

The aim of the study was an assessment of electromyographic (EMG) activity of knee muscles during isometric maximal voluntary contraction in the different disease durations of T2DM.

METHODS

Eighteen patients with less than 10 years and twelve patients with more than 10 years of T2DM were compared with nineteen matched healthy control subjects. EMG of flexor and extensor muscles of knee concurrently with isometric maximal peak torque of knee flexion and extension at 75 degrees of knee flexion were recorded in three groups.

RESULTS

Isometric maximal peak torque of extension and root mean squared (RMS) of vastus lateralis and medial hamstring in the healthy control was significantly higher than both patient groups. Whenever the maximal isometric peak flexion torque was not significantly different between groups, the mean power frequency (MPF) of flexor muscles especially medial hamstrings were higher in the short-term T2DM than healthy control groups. The two factors, gender and age, had significant effect on maximal peak torque and RMS of knee muscles.

CONCLUSIONS

EMG could show the effect of T2DM, gender and age on knee muscles activity. It seems that the medial hamstring was the most sensitive muscle of knee compartment to show the effect of T2DM and difference of short and long-term T2DM in EMG study.

The influence of confounding factors on the relationship between muscle contraction level and MF and MPF values of EMG signal: a review

The purpose of this article is to gather results of studies on the relationship between median frequency (MF) and mean power frequency (MPF) and the level of muscle contraction, and to use those results to discuss the differences in the trends according to factors related to measurement technique and subject. Twenty-one studies with 63 cases for upper limb muscles and nine studies with 31 cases for lower limb muscles were analysed. Most results showed an increase in parameters with an increased level of muscle contraction, only some studies showed a decrease. The influence on parameters of the level of muscle contraction and factors such as subjects, type of contraction, muscle length and electrodes was analysed for each muscle. It was concluded that when analysing the influence of different factors on MF and MPF, because those factors interact they should be considered together, not separately.

Changes in motor unit behavior following isometric fatigue of the first dorsal interosseous muscle

The neuromuscular strategies employed to compensate for fatigue-induced muscle force deficits are not clearly understood. This study utilizes surface electromyography (sEMG) together with recordings of a population of individual motor unit action potentials (MUAPs) to investigate potential compensatory alterations in motor unit (MU) behavior immediately following a sustained fatiguing contraction and after a recovery period. EMG activity was recorded during abduction of the first dorsal interosseous in 12 subjects at 20% maximum voluntary contraction (MVC), before and directly after a 30% MVC fatiguing contraction to task failure, with additional 20% MVC contractions following a 10-min rest. The amplitude, duration and mean firing rate (MFR) of MUAPs extracted with a sEMG decomposition system were analyzed, together with sEMG root-mean-square (RMS) amplitude and median frequency (MPF). MUAP duration and amplitude increased immediately postfatigue and were correlated with changes to sEMG MPF and RMS, respectively. After 10 min, MUAP duration and sEMG MPF recovered to prefatigue values but MUAP amplitude and sEMG RMS remained elevated. MU MFR and recruitment thresholds decreased postfatigue and recovered following rest. The increase in MUAP and sEMG amplitude likely reflects recruitment of larger MUs, while recruitment compression is an additional compensatory strategy directly postfatigue. Recovery of MU MFR in parallel with MUAP duration suggests a possible role for metabolically sensitive afferents in MFR depression postfatigue. This study provides insight into fatigue-induced neuromuscular changes by examining the properties of a large population of concurrently recorded single MUs and outlines possible compensatory strategies involving alterations in MU recruitment and MFR.

Segment lengths influence hill walking strategies

Segment lengths are known to influence walking kinematics and muscle activity patterns. During level walking at the same speed, taller individuals take longer, slower strides than shorter individuals. Based on this, we sought to determine if segment lengths also influenced hill walking strategies. We hypothesized that individuals with longer segments would display more joint flexion going uphill and more extension going downhill as well as greater lateral gastrocnemius and vastus lateralis activity in both directions. Twenty young adults of varying heights (below 155 cm to above 188 cm) walked at 1.25 m/s on a level treadmill as well as 6° and 12° up and downhill slopes while we collected kinematic and muscle activity data. Subsequently, we ran linear regressions for each of the variables with height, leg, thigh, and shank length. Despite our population having twice the anthropometric variability, the level and hill walking patterns matched closely with previous studies. While there were significant differences between level and hill walking, there were few hill walking variables that were correlated with segment length. In support of our hypothesis, taller individuals had greater knee and ankle flexion during uphill walking. However, the majority of the correlations were between tibialis anterior and lateral gastrocnemius activities and shank length. Contrary to our hypothesis, relative step length and muscle activity decreased with segment length, specifically shank length. In summary, it appears that individuals with shorter segments require greater propulsion and toe clearance during uphill walking as well as greater braking and stability during downhill walking.

Quantifying the familiarization period for maximal resistive exercise

Resistive exercise is used in the assessment of musculoskeletal health, performance, training interventions, and population differences (i.e., gender, age, training status). There is a need to determine the amount of familiarization required to stabilize performance prior to testing. Fifteen males completed a familiarization session consisting of 3 blocks of 5 maximal isometric dorsiflexion contractions, followed by a retention test (an additional block of 5 contractions) performed 3 days later. Mean force and surface electromyography (sEMG) from both the agonist and antagonist muscles were collected. A variance ratio, representing the stability between trials, was calculated for each of the 4 blocks of 5 contractions for both force and sEMG. The variance ratio for both force and agonist sEMG decreased significantly within the first 10 trials and remained stable during the retention test. The variance ratio for antagonist sEMG was stable across the 3 blocks of familiarization and significantly decreased during the retention test. The magnitude variables all remained stable across the 3 familiarization blocks. However, an 11% increase in mean force was seen during the retention test while both agonist and antagonist sEMG remained stable. Although slight changes occurred in the magnitude variables during the retention test, the stabilization of the force and agonist sEMG variance ratios suggest that familiarization to the task was achieved within the first 10 contractions and was sustained over a 3-day period.

A simulation study on the relation between muscle motor unit numbers and the non-Gaussianity/non-linearity levels of surface electromyography

Recent research has demonstrated that surface electromyography (sEMG) signals have non-Gaussianity and non-linearity properties. It is known that more muscle motor units are recruited and firing rates (FRs) increase as exertion increases. A hypothesis was proposed that the Gaussianity test (S (g)) and linearity test (S (ℓ)) levels of sEMG signals are associated with the number of active motor units (nMUs) and the FR. The hypothesis has only been preliminarily discussed in experimental studies. We used a simulation sEMG model involving spatial (active MUs) and temporal (three FRs) information to test the hypothesis. Higher-order statistics (HOS) from the bi-frequency domain were used to perform S (g) and S (ℓ). Multivariate covariance analysis and a correlation test were employed to determine the nMUs-S (g) relationship and the nMUs-S (ℓ) relationship. Results showed that nMUs, the FR, and the interaction of nMUs and the FR all influenced the S (g) and S (ℓ) values. The nMUs negatively correlated to both the S (g) and S (ℓ) values. That is, at the three FRs, sEMG signals tended to a more Gaussian and linear distribution as exertion and nMUs increased. The study limited experiment factors to the sEMG non-Gaussianity and non-linearity levels. The study quantitatively described nMUs and the FR of muscle that are not directly available from experiments. Our finding has guiding significance for muscle capability assessment and prosthetic control.

IMPACT OF FINGERTIP ACTIONS ON TOTAL POWER OF SURFACE ELECTROMYOGRAPHY FROM EXTRINSIC HAND MUSCLES

Quantizing the relationship between finger force and multitendoned extrinsic hand muscles could be useful for understanding the control strategies that underlie the coordination of finger movements and forces. The objective of this study is to explore the relationship of fingertip force production and total power of surface electromyography (sEMG) recorded on extrinsic hand muscles under isometric voluntary contraction. Thirteen healthy volunteers were recruited to participate in this study. In the designed force-tracking tasks, all volunteers were required to produce a certain force with either index finger or middle finger to match the target force for 5 s. Meanwhile, the sEMG signals were acquired from two extrinsic hand muscles: extensor digitorum (ED) and flexor digitorum superficialis (FDS). For each trial, sEMG of the effective force segment was extracted; then, the power spectrum was estimated based on autoregressive (AR) model and from which the corresponding total power of sEMG was computed. The experimental results reveal that the total power of sEMG linearly increases with force level regardless of the task finger and extrinsic hand muscle. It is also found that the total power obtained from index finger is significantly less than that of middle finger for FDS at the same force level (p < 0.05), while this kind of statistical significance cannot be found for ED. However, with respect to the measurement of total power, the type of extrinsic hand muscle has not exhibited significantly different contribution to the task finger under a certain fingertip force level. The findings of this study indicate that the total power of the extrinsic hand muscle's sEMG can be used to characterize finger's activities.

Variability of neural activation during walking in humans: short heels and big calves

People come in different shapes and sizes. In particular, calf muscle size in humans varies considerably. One possible cause for the different shapes of calf muscles is the inherent difference in neural signals sent to these muscles during walking. In sedentary adults, the variability in neural control of the calf muscles was examined with muscle size, walking kinematics and limb morphometrics. Half the subjects walked while activating their medial gastrocnemius (MG) muscles more strongly than their lateral gastrocnemius (LG) muscles during most walking speeds ('MG-biased'). The other subjects walked while activating their MG and LG muscles nearly equally ('unbiased'). Those who walked with an MG-biased recruitment pattern also had thicker MG muscles and shorter heel lengths, or MG muscle moment arms, than unbiased walkers, but were similar in height, weight, lower limb length, foot length, and exhibited similar walking kinematics. The relatively less plastic skeletal system may drive calf muscle size and motor recruitment patterns of walking in humans.

Experimental and modelling investigation of surface EMG spike analysis

A pattern classification method based on five measures extracted from the surface electromyographic (sEMG) signal is used to provide a unique characterization of the interference pattern for different motor unit behaviours. This study investigated the sensitivity of the five sEMG measures during the force gradation process. Tissue and electrode filtering effects were further evaluated using a sEMG model. Subjects (N=8) performed isometric elbow flexion contractions from 0 to 100% MVC. The sEMG signals from the biceps brachii were recorded simultaneously with force. The basic building block of the sEMG model was the detection of single fibre action potentials (SFAPs) through a homogeneous, equivalent isotropic, infinite volume conduction medium. The SFAPs were summed to generate single motor unit action potentials. The physiologic properties from a well-known muscle model and motor unit recruitment and firing rate schemes were combined to generate synthetic sEMG signals. The following pattern classification measures were calculated: mean spike amplitude, mean spike frequency, mean spike slope, mean spike duration, and the mean number of peaks per spike. Root-mean-square amplitude and mean power frequency were also calculated. Taken together, the experimental data and modelling analysis showed that below 50% MVC, the pattern classification measures were more sensitive to changes in force than traditional time and frequency measures. However, there are additional limitations associated with electrode distance from the source that must be explored further. Future experimental work should ensure that the inter-electrode distance is no greater than 1cm to mitigate the effects of tissue filtering.

Relationships between surface EMG variables and motor unit firing rates

Although surface electromyography (sEMG) is a widely used electrophysiological technique, its physiological interpretation remains somewhat controversial. This study examined the relationship between motor unit firing rates (MUFR) and the root mean square (RMS) amplitude and mean power frequency (MPF) of the sEMG signal in the biceps brachii. Eleven subjects performed maximal isometric elbow flexion while indwelling and sEMG recordings were obtained from the biceps. The RMS amplitude and MPF of the surface signal, and the mean MUFR from the indwelling signal, were calculated over 500 ms epochs. Group means showed a strong MUFR-RMS amplitude relationship (r (2) = 0.91), but a weak MUFR-MPF relationship (r (2) = 0.20). Using all trials, the MUFR-RMS amplitude (r (2) = 0.19) and MUFR-MPF (r (2) = 0.0037) relationships were much weaker. Within individual subjects, the MUFR-RMS amplitude (mean r (2) = 0.13 +/- 0.17) and the MUFR-MPF (mean r (2) = 0.040 +/- 0.041) relationships were also weak. These results suggest that MUFR cannot be predicted from the characteristics of the sEMG signal.

Motor unit synchronization in FDI and biceps brachii muscles of strength-trained males

Motor unit (MU) synchronization is the simultaneous or near-simultaneous firing of two MUs which occurs more often than would be expected by chance. The present study sought to investigate the effects of exercise training, muscle group, and force level, by comparing the magnitude of synchronization in the biceps brachii (BB) and first dorsal interosseous (FDI) muscles of untrained and strength-trained college-aged males at two force levels, 30% of maximal voluntary contraction (MVC) and 80% MVC. MU action potentials were recorded directly via an intramuscular needle electrode. The magnitude of synchronization was assessed using previously-reported synchronization indices: k', E, and CIS. Synchronization was significantly higher in the FDI than in the BB. Greater synchronization was observed in the strength-trained group with CIS, but not with E or k'. Also, synchronization was significantly greater at 80% MVC than at 30% MVC with E, but only moderately greater with CIS and there was no force difference with k'. Synchronization prevalence was found to be greater in the BB (80.1%) than in the FDI (71.5%). Thus, although the evidence is a bit equivocal, it appears that MU synchronization is greater at higher forces, and greater in strength-trained individuals than in untrained subjects.