Reliability of muscle fibre conduction velocity distribution estimation from surface EMG

Most of the neuromuscular diseases induce changes in muscle fibre characteristics. For example, Duchenne dystrophy is characterized by a specific loss of fast fibres, and an increase in small diameter fibres. These morphological changes may lead to large modifications in the distribution of fibre diameters, possibly producing bimodal distributions. It has already been shown that it is possible to reveal these morphological modifications through the distribution of muscle fibre conduction velocity (MFCV) as estimated from needle electromyography (EMG) recordings. In this paper, we investigate whether such changes can be extracted from surface EMG signals.

Simulation allows generation of surface EMG signals in which features are well described especially at a morphological level. Therefore, we generated a database of simulated signals both in voluntary and electrically elicited contraction conditions using a bimodal distribution of muscle fibre diameters. MFCV distributions were computed using two short-term methods based on cross-correlation and peak-to-peak techniques for voluntary contraction signals, and using a deconvolution method in time domain for electrically elicited signals. MFCV distributions were compared with true ones, as generated from modelling.

This work reveals that estimating MFCV distribution through these methods does not appear yet as precise and robust enough to accurately characterize changes in redistribution of various muscle fibre diameters. However, it opens to new experimental protocols that can be explored in order to improve the robustness of MFCV distribution estimation for the follow-up of patients suffering from neuromuscular disorders.     

Distribution of innervation zone and muscle fiber conduction velocity in the biceps brachii muscle

The spatial distributions of muscle innervation zone (IZ) and muscle fiber conduction velocity (CV) were examined in nine healthy young male participants. High-density surface electromyography (EMG) was collected from the biceps brachii muscle when subjects performed isometric elbow flexions at 20% to 80% of the maximal voluntary contraction (MVC). A total of 9498 samples of IZs were identified and CVs were calculated using the Radon transform. The center and width of IZ sample distribution were compared within four different force levels and six medial to lateral electrode column positions using repeated measures ANOVA and multiple comparison tests. Significant shifts of IZ center were observed in the medial columns (Columns 5, 6, and 7) compared with the lateral columns (Columns 3 and 4) (p < 0.05). Similarly, significant differences in the IZ width were found in Column 7 and 8 compared to Column 3 (p < 0.05). In contrast, muscle CV was unaffected by column position. Instead, muscle CV was faster at 40% and 80% MVC compared to 20% MVC (p < 0.05). The findings of this study add further insights into the physiological properties of the biceps brachii muscle.     

Estimation of muscle fiber conduction velocity of doublet discharges

A doublet is defined as two consecutive discharges of a motor unit occurring at short time interval between each other (e.g., <20 ms). In this paper, we propose a method for the estimation of muscle fiber conduction velocity (CV) from two partly overlapping action potentials generated by the same motor units. The method is based on the minimization of the mean square error between time-filtered versions of two surface EMG signals recorded along the direction of muscle fibers. The minimization is performed over the filter parameters that define the two propagation delays. The method was tested on simulated and experimental signals. Simulation results showed that the method is only in some cases superior to the simpler peak approach, due to limitations in the ideal model used for the algorithm development. However, application to experimental signals that mimic doublet motor unit discharges showed a substantial improvement in estimation quality of the new method with respect to the peak method.     

Spatial distribution of active muscle fibre characteristics in the upper trapezius muscle and its dependency on contraction level and duration

The aim of this study was to provide direct in vivo information of the physiological and structural characteristics of active muscle fibres from a large part of the upper trapezius muscle. Two-dimensional (2-D) multi-channel surface electromyography recordings were used, with 13 × 10 electrodes covering 6 × 4.5 cm of the skin’s surface. A previously developed method was applied to detect individual propagating motor unit action potentials and to estimate their corresponding muscle fibre conduction velocity (MFCV) and muscle fibre orientation (MFO). Using these estimates, spatial distributions of MFCV and MFO were examined for five male subjects performing isometric shoulder elevation at different force levels. The main results were: (1) the general relationship between MFCV and force generation was non-systematic, with a positive relationship at the inferior part of the muscle, (2) the spatial distribution of MFCV at different force levels and fatigue was inhomogeneous and (3) the MFO was slightly different (6°) of the muscle fibres with origin superior compared to inferior to the C7 vertebra. These findings provide new information of the MFO of contracting muscle fibres and knowledge of the physiological characteristics of a large part of the upper trapezius muscle that previously was based on observations from human cadavers only.     

Spatio-temporal processing of surface EMG signals from the sternocleidomastoideus muscle to assess effects of radiotherapy on motor unit conduction velocity and firing rate—A pilot study

Radiation therapy causes both muscle and nerve tissue damage. However, the evolution and mechanisms of these damages are not fully understood. Information on the state of active muscle fibres and motoneurons can be obtained by measuring sEMG signals and calculating the conduction velocity (CV) and firing rate of individual motor units, respectively. The aim of this pilot study was to evaluate if the multi-channel surface EMG (sEMG) technique could be applied to the sternocleidomastoideus muscle (SCM) of radiotherapy patients, and to assess if the CV and firing rate are altered as a consequence of the radiation.

Surface EMG signals were recorded from the radiated and healthy SCM muscles of 10 subjects, while subjects performed isometric rotation of the head. CV and firing rate were calculated using two recently proposed methods based on spatio-temporal processing of the sEMG signals. The multi-channel sEMG technique was successfully applied to the SCM muscle and CV and firing rates were obtained. The measurements were fast and simple and comfortable for the patients. Sufficient data quality was obtained from both sides of seven and four subjects for the CV and firing rate analysis, respectively. No differences in CV or firing rate were found between the radiated and non-radiated sides (p = 0.13 and p = 0.20, respectively). Firing rate and CV were also obtained from a myokymic discharge pattern. It was found that the CV decreased significantly (p = 0.01) during the bursts.     

Dependence of average muscle fibre conduction velocity on voluntary contraction force

Average muscle fibre conduction velocity (CV) measured with multichannel surface electrodes decreases with time during sustained isometric contraction. Based on this property, CV is considered a candidate for an objective index to localized muscular fatigue. CV, however, also depends on many other factors that include muscle temperature and voluntary contraction force. In this paper, the effect of contraction force on CV was studied by defining not only the target force level but also the whole force trajectory. The contraction was isometric and lasted 14 s. The target force was set at four levels from 30% to 90% of the maximal voluntary contraction (MVC). Three typical muscles were studied in seven healthy male subjects. In the vastus lateralis, CV increased with contraction force in many cases. In the biceps brachii, CV decreased rapidly with time before the contraction force reached the target levels of 70% or 90% MVC. At these force levels, CV was smaller than that at 50% MVC. CV in the biceps consequently showed no apparent dependence on the contraction force. The tibialis anterior showed intermediate change in CV between the vastus lateralis and the biceps brachii. These results indicate that CV basically increases with contraction force, but this relationship becomes unclear when CV decreases rapidly with time.     

Low back muscle action potential conduction velocity estimated using high-density electromyography

While a decreasing spectral content of surface electromyography reflects low back muscle fatigue development, reliability of these decreases may be insufficient. Decreasing frequency content is largely determined by decreasing average motor unit action potential conduction velocities (CV), which is considered a more direct measure of muscle fatigue development. However, for the low back muscles it has been proven difficult to identify propagating potentials and consequently estimate the CV. The aim of this study was to estimate the low back muscle CV from high-density multi-channel electromyography by using peak-delay and cross-correlation methods. Fourteen healthy male participants without a history of low-back pain performed a 30 degrees lumbar flexion trial until exhaustion while standing. For 10 out of the 14 participants (118 out of 560 sites) realistic CV estimates were obtained using both methods, the majority likely over the iliocostalis lumborum muscle. Between-method CV differences appeared to be small. Close to the spine a considerable number of sites (79) yielded systematically overestimated low back muscle CV values. Estimating low back muscle CV may allow additional insight into low back muscle fatigue development and potentially improve its monitoring using (high-density) surface electromyography.     

The effect of exercise induced hyperthermia on muscle fibre conduction velocity during sustained isometric contraction

This study investigated the effect of dynamic exercise in a hot environment on muscle fibre conduction velocity (MFCV) of the knee extensors during a sustained isometric contraction. Seven trained male cyclists (mean [±SD], age, and were 35 ± 9.9 and 57.4 ± 6.6 ml kg⁻¹ min⁻¹) cycled for 50 min at 60% of peak power output in either: (1) 40 °C (HOT); or (2) 19 °C (NEUTRO); and (3) remained passive in 40 °C (PASS). Post-intervention a 100 s maximal sustained isometric contraction (SMC) of the knee extensors was performed. Rectal temperature increased (p < 0.01) for both HOT and NEUTRO with PASS unchanged and with HOT rising higher (p < 0.01) than NEUTRO (38.6 ± 0.4 vs. 37.6 ± 0.4 °C). Muscle temperature increased (p < 0.01) for all three conditions with HOT rising the highest (p < 0.01) (40.3 ± 0.5 vs. 38.3 ± 0.3 and 37.6 ± 1.3 °C for NEUTRO and PASS, respectively). Lactate showed higher accumulation (p < 0.01) for HOT than NEUTRO (6.9 ± 2.3 vs. 4.2 ± 2.1 mmol l⁻¹). During SMC the torque, electromyography root mean squared (RMS) and MFCV all significantly (p < 0.01) declined. Only in HOT did MFCV decline significantly (p < 0.01) less than torque and RMS (9.9 ± 6.2% vs. 37.5 ± 17.8% and 37.6 ± 21.4%, respectively). In conclusion, during exercise induced hyperthermia, reduced motor unit recruitment as opposed to slower conducting properties of the muscle fibre appears to be responsible for the greater reduction in torque output.     

Muscle fiber conduction velocity estimation by using normalized peak-averaging technique

In order to compute the muscle fiber conduction velocity (MFCV) and to clarify how action potentials are conducted, the normalized peak-averaging technique (NPAT) was newly employed together with computer softwares. Twelve pairs of surface electromyograms were selected from biceps brachii muscles during contraction at a level of 50% of the maximum voluntary isometric contraction in seven healthy volunteers. The techniques to compute MFCV from the time delay of the peaks (P-NPAT) and from the cross correlation (CC-NPAT) of averaged pulses were compared to the cross-correlation technique (CCT). The spread rate of averaged pulses was computed to estimate the spread of MFCVs in different motor units. Tri-phasic averaged pulses were obtained clearly by averaging more than 500 detected pulses. The P-NPAT and CC-NPAT highly correlated with the CCT in the computed MFCVs. The MFCVs obtained by P-NPAT were generally larger than those obtained by CCT, and the spread rates had in the definite values. These results suggest that the MFCV could be computed and the spread of MFCVs would be estimated from averaged pulses. The MFCV of a patient with myotonic dystrophy was also studied, and it was suggested that the NPAT would be clinically useful.     

Clinical applications of high-density surface EMG: A systematic review

High density-surface EMG (HD-sEMG) is a non-invasive technique to measure electrical muscle activity with multiple (more than two) closely spaced electrodes overlying a restricted area of the skin. Besides temporal activity HD-sEMG also allows spatial EMG activity to be recorded, thus expanding the possibilities to detect new muscle characteristics. Especially muscle fiber conduction velocity (MFCV) measurements and the evaluation of single motor unit (MU) characteristics come into view. This systematic review of the literature evaluates the clinical applications of HD-sEMG. Although beyond the scope of the present review, the search yielded a large number of “non-clinical” papers demonstrating that a considarable amount of work has been done and that significant technical progress has been made concerning the feasibility and optimization of HD-sEMG techniques. Twenty-nine clinical studies and four reviews of clinical applications of HD-sEMG were considered. The clinical studies concerned muscle fatigue, motor neuron diseases (MND), neuropathies, myopathies (mainly in patients with channelopathies), spontaneous muscle activity and MU firing rates. In principle, HD-sEMG allows pathological changes at the MU level to be detected, especially changes in neurogenic disorders and channelopathies. We additionally discuss several bioengineering aspects and future clinical applications of the technique and provide recommendations for further development and implementation of HD-sEMG as a clinical diagnostic tool.     

Muscle fiber conduction velocity in different gait phases of early and late-stage diabetic neuropathy

We investigated the muscle fiber conduction velocity (MFCV) during gait phases of the lower limb muscles in individuals with various degrees of diabetic peripheral neuropathy (DPN). Forty-five patients were classified into severity degrees of DPN by a fuzzy model. The stages were absent (n = 11), mild (n = 14), moderate (n = 11) and severe (n = 9), with 10 matched healthy controls. While walking, all subjects had their sEMG (4 linear electrode arrays) recorded for tibialis anterior (TA), gastrocnemius medialis (GM), vastus lateralis (VL) and biceps femoris (BF). MFCV was calculated using a maximum likelihood algorithm with 30 ms standard deviation Gaussian windows. In general, individuals in the earlier stages of DPN showed lower MFCV of TA, GM and BF, whilst individuals with severe DPN presented higher MFCV of the same muscles. We observed that mild patients already showed lower MFCV of TA at early stance and swing, and lower MFCV of BF at swing. All diabetic groups showed a markedly reduction in MFCV of VL, irrespective of DPN. Severe patients presented higher MFCV mainly in distal muscles, TA at early and swing phases and GM at propulsion and midstance. The absent group already showed MFCV of VL and GM reductions at the propulsion phase and of VL at early stance. Although MFCV changes were not as progressive as the DPN was, we clearly distinguished diabetic patients from controls, and severe patients from all others.     

The effect of temperature on conduction velocity in human muscle fibers

The effects of variation of intramuscular temperature (T) on conduction velocity (CV) of the action potential along single human muscle fibers of the biceps brachii was studied in situ in 15 normal volunteers (mean age 39 years, range 21–62 years). Cooling was obtained by direct application of ice over a rectangular skin region including the stimulating and recording area. The intramuscular T was monitored by a needle thermocouple (copperconstantane). In all the 24 muscle fibers studied, a linear relationship was observed between CV and T. The slopes of the regression lines, ranging between 0.190 and 0.079 m/s, were positively correlated with the starting CV at 36°C ranging between 2.2 and 5.2 m/s. If conduction changes are expressed as a percentage of the basal CV at 36°C, the CV/T coefficient is the same for all the fibers and independent of the individual CV: 3.4% of CV/°C.     

Task-related differences in peroneus longus muscle fiber conduction velocity

It has been identified that the peroneus longus presents a regional activity. Specifically, a greater activation of the anterior and posterior compartments has been observed during eversion, whereas a lower activation of the posterior compartment has been reported during plantarflexion. In addition to myoelectrical amplitude, motor unit recruitment can be inferred indirectly from muscle fiber conduction velocity (MFCV). However, there are few reports of MFCV of the regions that make up a muscle, and even less, MFCV of the peroneus longus compartments. This study aimed to analyze the MFCV of peroneus longus compartments during eversion and plantarflexion. Twenty-one healthy individuals were assessed. High-density surface electromyography was recorded from the peroneus longus during eversion and plantarflexion at 10%, 30%, 50%, and 70% of maximal voluntary isometric contraction. The posterior compartment presented a lower MFCV than the anterior compartment during plantarflexion, and both compartments did not show differences in MFCV during eversion; however, the posterior compartment showed an increase in MFCV during eversion compared to plantarflexion. Differences observed in the MFCV of the peroneus longus compartments could support a regional activation strategy and, to some extent, explain different motor unit recruitment strategies of the peroneus longus during ankle movements.     

Multichannel EMG-based estimation of fiber conduction velocity during isometric contraction of patients with different stages of diabetic neuropathy

This study compares muscle fiber conduction velocities estimated using surface electromyography during isometric maximal voluntary contraction in different stages of diabetic neuropathy. Eighty-five adults were studied: 16 non-diabetic individuals and 69 diabetic patients classified into four neuropathy stages, defined by a fuzzy expert system: absent (n = 26), mild (n = 21), moderate (n = 11) and severe (n = 11). Average muscle fiber conduction velocities of gastrocnemius medialis, tibialis anterior, vastus lateralis and biceps femoris were assessed using linear array electrodes, and were compared by ANOVA. Conduction velocities were significantly decreased in the moderate neuropathy group for the vastus lateralis compared to other groups (from 18% to 21% decrease), and were also decreased in all diabetic groups for the tibialis anterior (from 15% to 20% from control group). Not only the distal anatomical localization of the muscle affects the conduction velocity, but also the proportion of muscle fiber type, where the tibialis anterior with greater type I fiber proportion is affected earlier while the vastus lateralis with greater type II fiber proportion is affected in later stages of the disease. Generally, the muscles of the lower limb have different responsiveness to the effects of diabetes mellitus and show a reduction in the conduction velocity as neuropathy progresses.     

Reproducibility of muscle fibre conduction velocity during linearly increasing force contractions

Muscle fibre conduction velocity (MFCV) is a basic physiological parameter biophysically related to the diameter of muscle fibres and properties of the sarcolemma. The aim of this study was to assess the intersession reproducibility of the relation between voluntary force and estimates of average muscle fibre conduction velocity (MFCV) from multichannel high-density surface electromyographic recordings (HDsEMG). Ten healthy men performed six linearly increasing isometric ankle dorsiflexions on two separate experimental sessions, 4 weeks apart. Each session involved the recordings of voluntary force during maximal isometric (MViF) and submaximal ramp contractions at 35–50–70% of MViF. Concurrently, the HDsEMG activity was detected from the tibialis anterior muscle and MFCV estimates were derived in 250-ms epochs. Absolute and relative reproducibility of MFCV initial value (intercept) and rate of change (regression slope) as a function of force were assessed by within-subject coefficient of correlation (CV_w) and with intraclass correlation coefficient (ICC). MFCV was positively correlated with voluntary force (R² = 0.75 ± 0.12) in all individuals and test conditions (P < 0.001). Average CV_w for MFCV intercept and slope were of 2.6 ± 2.0% and 11.9 ± 3.2% and ICC values of 0.96 and 0.94, respectively.Overall, MFCV regression coefficients showed a high degree of intersession reproducibility in both absolute and relative terms. These results may have important practical implications in the tracking of training-induced neuromuscular changes and/or in the monitoring of the progress of neuromuscular disorders when a full sEMG signal decomposition is problematic or not possible.     

Character of the effect of microwave on conduction velocity of frog ventricular muscle

Thirty-two frog hearts were divided into four groups and placed individually in temperature-controlled waveguides filled with Ringer's solution. The pacemaker was removed, and stimulation was provided at 0.3 Hz by three carbon-loaded Teflon electrodes located on the aorta and the ventricular muscle. Conduction velocity was measured from the difference between two action potentials. One group served as control; the three other groups were exposed for 2 h to pulsed 2,450 MHz microwave fields (10 μs, 0.001 duty cycle, 16 Hz modulation) at specific absorption rates (SARs) of 0.003, 2, and 6 W/kg, respectively. No significant difference in conduction velocity was found between the control and the exposed groups. © 1994 Wiley-Liss, Inc.     

Comparison of the duration and power spectral changes of monopolar and bipolar M waves caused by alterations in muscle fibre conduction velocity

The muscle compound action potential (M wave) recorded under monopolar configuration reflects both the propagation of the action potentials along the muscle fibres and their extinction at the tendon. M waves recorded under a bipolar configuration contain less cross talk and noise than monopolar M waves, but they do not contain the entire informative content of the propagating potential. The objective of this study was to compare the effect of changes in muscle fibre conduction velocity (MFCV) on monopolar and bipolar M waves and how this effect depends on the distance between the recording electrodes and tendon. The study was based on a simulation approach and on an experimental investigation of the characteristics of surface M waves evoked in the vastus lateralis during 4-s step-wise isometric contractions in knee extension at 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90% MVC. The peak-to-peak duration (Durpp) and median frequency (Fmedian) of the M waves were calculated. For monopolar M waves, changes in Durpp and Fmedian produced by MFCV depended on the distance from the electrode to the tendon, whereas, for bipolar M waves, changes in Durpp and Fmedian were largely independent of the electrode-to-tendon distance. When the distance between the detection point and tendon lay between approximately 15 and 40 mm, changes in Durpp of bipolar M waves were more pronounced than those of distal monopolar M waves but less marked than those of proximal monopolar M waves, and the opposite occurred for Fmedian. Since, for bipolar M waves, changes in duration and power spectral features produced by alterations in MFCV are not influenced by the electrode-to-tendon distance, the bipolar electrode configuration is a preferable choice over monopolar arrangements to estimate changes in conduction velocity.     

Separation of propagating and non propagating components in surface EMG

Surface electromyogram (EMG) detected by electrode arrays along the muscle fibre direction can be approximated by the sum of propagating and non propagating components. A technique to separate propagating and non propagating components in surface EMG signals is developed. The first step is an adaptive filter, which allows obtaining an estimation of the delay between signals detected at different channels and a first estimate of propagating and non propagating components; the second step is used to optimise the estimation of the two components. The method is applicable to signals with one propagating and one non propagating component. It was optimised on simulated signals, and then applied to single motor unit action potentials (MUAP) and to electrically elicited EMG (M-waves).

The new method was first tested on phenomenological signals constituted by the sum of a propagating and a non propagating signal and then applied to simulated and experimental EMG signals. Simulated signals were generated by a cylindrical, layered volume conductor model. Experimental signals were monopolar surface EMG signals collected from the abductor pollicis brevis muscle and M-waves recorded during transcutaneous electrical stimulation of the biceps muscle. The technique may find different applications: in single motor unit (MU) studies (a) for decreasing the variability and bias of CV estimates due to the presence of the non propagating components, (b) for estimating automatically the length of the muscle fibres from only three detected channels and (c) for removal of the stimulation artifact M-waves.     

Network modeling and analysis of lumbar muscle surface EMG signals during flexion–extension in individuals with and without low back pain

In this paper, we propose modeling the activity coordination network between lumbar muscles using surface electromyography (sEMG) signals and performing the network analysis to compare the lumbar muscle coordination patterns between patients with low back pain (LBP) and healthy control subjects. Ten healthy subjects and eleven LBP patients were asked to perform flexion–extension task, and the sEMG signals were recorded. Both the subject-level and the group-level PC_fdr algorithms are applied to learn the sEMG coordination networks with the error-rate being controlled. The network features are further characterized in terms of network symmetry, global efficiency, clustering coefficient and graph modules. The results indicate that the networks representing the normal group are much closer to the order networks and clearly exhibit globally symmetric patterns between the left and right sEMG channels. While the coordination activities between sEMG channels for the patient group are more likely to cluster locally and the group network shows the loss of global symmetric patterns. As a complementary tool to the physical and anatomical analysis, the proposed network analysis approach allows the visualization of the muscle coordination activities and the extraction of more informative features from the sEMG data for low back pain studies.     

Spatial EMG potential distribution pattern of vastus lateralis muscle during isometric knee extension in young and elderly men

The aim of the present study was to compare spatial electromyographic (EMG) potential distribution during force production between elderly and young individuals using multi-channel surface EMG (SEMG). Thirteen elderly (72-79 years) and 13 young (21-27 years) healthy male volunteers performed ramp submaximal contraction during isometric knee extension from 0% to 65% of maximal voluntary contraction. During contraction, multi-channel EMG was recorded from the vastus lateralis muscle. To evaluate alteration in heterogeneity and pattern in spatial EMG potential distribution, coefficient of variation (CoV), modified entropy and correlation coefficients with initial torque level were calculated from multi-channel SEMG at 5% force increment. Increase in CoV and decrease in modified entropy of RMS with increase of exerted torque were significantly smaller in elderly group (p < 0.05) and correlation coefficients with initial torque level were significantly higher in elderly group than in young group at moderate torque levels (p < 0.05). These data suggest that the increase of heterogeneity and the change in the activation pattern are smaller in elderly individuals than in young individuals. We speculated that multi-channel SEMG pattern in elderly individual reflects neuromuscular activation strategy regulated predominantly by clustering of similar type of muscle fibers in aged muscle.