Search

The EMG filling factor is an index to quantify the degree to which an EMG signal has been filled. Here, we tested the validity of such index to analyse the EMG filling process as contraction force was slowly increased. Surface EMG signals were recorded from the quadriceps muscles of healthy subjects as force was gradually increased from 0 to 40% MVC. The sEMG filling process was analyzed by measuring the EMG filling factor (calculated from the non-central moments of the rectified sEMG). (1) As force was gradually increased, one or two prominent abrupt jumps in sEMG amplitude appeared between 0 and 10% of MVC force in all the vastus lateralis and medialis. (2) The jumps in amplitude were originated when a few large-amplitude MUPs, clearly standing out from previous activity, appeared in the sEMG signal. (3) Every time an abrupt jump in sEMG amplitude occurred, a new stage of sEMG filling was initiated. (4) The sEMG was almost completely filled at 2–12% MVC. (5) The filling factor decreased significantly upon the occurrence of an sEMG amplitude jump, and increased as additional MUPs were added to the sEMG signal. (6) The filling factor curve was highly repeatable across repetitions. It has been validated that the filling factor is a useful, reliable tool to analyse the sEMG filling process. As force was gradually increased in the vastus muscles, the sEMG filling process occurred in one or two stages due to the presence of abrupt jumps in sEMG amplitude. [ABSTRACT FROM AUTHOR]

The motor unit profile, a representation of the trajectories of positive and negative turns of a scanning-EMG signal, is a new way to characterize the motor unit potential. Such characterization allows quantification of the scanning-EMG signal's complexity, which is closely related to the anatomy and physiology of the motor unit. To extract the motor unit profile, an algorithm that detects the turns of the scanning-EMG signal and links them using point-tracking techniques has been developed. The performance of this algorithm is sensitive to three parameters: the turn detection threshold, the maximum tracking interval threshold, and the trajectory purge threshold. Real scanning-EMG signals have been used to analyze the algorithm's behavior and the influence of the algorithm's parameters and to determine which parameter values provide the best performance. [ABSTRACT FROM AUTHOR]

The transient enlargement of the compound muscle action potential (M wave) after a conditioning contraction is referred to as potentiation. It has been recently shown that the potentiation of the first and second phases of a monopolar M wave differed drastically; namely, the first phase remained largely unchanged, whereas the second phase underwent a marked enlargement and shortening. This dissimilar potentiation of the first and second phases has been suggested to be attributed to a transient increase in conduction velocity after the contraction. Here, we present a series of simulations to test if changes in the timing variability between motor unit potentials (MUPs) can be responsible for the unequal potentiation (and shortening) of the first and the second M-wave phases. We found that an increase in the mean motor unit conduction velocity resulted in a marked enlargement and narrowing of both the first and second M-wave phases. The enlargement of the first phase caused by a global increase in motor unit conduction velocities was apparent even for the electrode located over the innervation zone and became more pronounced with increasing distance to the innervation zone, whereas the potentiation of the second phase was largely independent of electrode position. Our simulations indicate that it is unlikely that an increase in motor unit conduction velocities (accompanied or not by changes in their distribution) could account for the experimental observation that only the second phase of a monopolar M wave, but not the first, is enlarged after a brief contraction. However, the combination of an increase in the motor unit conduction velocities and a spreading of the motor unit activation times could potentially explain the asymmetric potentiation of the M-wave phases. [ABSTRACT FROM AUTHOR]

In quantitative electromyography (EMG), the set of potentials that constitute a motor unit action potential (MUAP) train are represented by a single waveform from which various parameters are determined in order to characterize the MUAP for diagnostic analysis. Several methods that extract such a waveform are currently available, and they are, in essence, based on two operations: averaging and selection, which are performed either sample-by-sample or on the whole-potential. We present a new approach that carries out selection and averaging on a local interval basis. We tested our algorithm with a dataset of MUAP records extracted from the tibialis anterioris muscle of healthy subjects and compared it with some of the most relevant state-of-the-art methods considered in a previous work (Malanda et al., J. Electromyogr. Kinesiol., 2015). The comparison covered general purpose signal processing figures of merit and clinically used MUAP waveform parameters. Significantly better results in both sets of figures of merit were obtained with the new approach. In addition, relative to the other algorithms tested, the new approach required fewer potentials from the MUAP set to obtain an accurate representative waveform. [ABSTRACT FROM AUTHOR]

In the context of quantitative electromyography (EMG), it is of major interest to obtain a waveform that faithfully represents the set of potentials that constitute a motor unit action potential (MUAP) train. From this waveform, various parameters can be determined in order to characterize the MUAP for diagnostic analysis. The aim of this work was to conduct a thorough, in-depth review, evaluation and comparison of state-of-the-art methods for composing waveforms representative of MUAP trains. We evaluated nine averaging methods: Ensemble (EA), Median (MA), Weighted (WA), Five-closest (FCA), MultiMUP (MMA), Split-sweep median (SSMA), Sorted (SA), Trimmed (TA) and Robust (RA) in terms of three general-purpose signal processing figures of merit (SPMF) and seven clinically-used MUAP waveform parameters (MWP). The convergence rate of the methods was assessed as the number of potentials per MUAP train (NPM) required to reach a level of performance that was not significantly improved by increasing this number. Test material comprised 78 MUAP trains obtained from the tibialis anterioris of seven healthy subjects. Error measurements related to all SPMF and MWP parameters except MUAP amplitude descended asymptotically with increasing NPM for all methods. MUAP amplitude showed a consistent bias (around 4% for EA and SA and 1–2% for the rest). MA, TA and SSMA had the lowest SPMF and MWP error figures. Therefore, these methods most accurately preserve and represent MUAP physiological information of utility in clinical medical practice. The other methods, particularly WA, performed noticeably worse. Convergence rate was similar for all methods, with NPM values averaged among the nine methods, which ranged from 10 to 40, depending on the waveform parameter evaluated. [ABSTRACT FROM AUTHOR]

Abstract: Muscle fatigue is a common experience in daily life. Many authors have defined it as the incapacity to maintain the required or expected force, and therefore, force, power and torque recordings have been used as direct measurements of muscle fatigue. In addition, the measurement of these variables combined with the measurement of surface electromyography (sEMG) recordings (which can be measured during all types of movements) during exercise may be useful to assess and understand muscle fatigue. Therefore, there is a need to develop muscle fatigue models that relate changes in sEMG variables with muscle fatigue. However, the main issue when using conventional sEMG variables to quantify fatigue is their poor association with direct measures of fatigue. Therefore, using different techniques, several authors have combined sets of sEMG parameters to assess muscle fatigue. The aim of this paper is to serve as a state-of-the-art summary of different sEMG models used to assess muscle fatigue. This paper provides an overview of linear and non-linear sEMG models for estimating muscle fatigue, their ability to assess power loss and their limitations due to neuromuscular changes after a training period. [Copyright &y& Elsevier]

Abstract: In a previous work we studied the ratio between the amplitudes of the second and first phases (which we call PPR, after peak-to-peak ratio) of the single fibre action potential (SFAP) for a collection of fibrillation potentials (FPs) extracted from two pathological muscles. These FPs showed a wider PPR range than the Dimitrov–Dimitrova (D–D) convolutional model could provide. We proposed a modification of the D–D intracellular action potential (IAP) in order to obtain a range of PPRs comparable to that observed in our FPs. This paper extends that study to a large number of SFAPs extracted from the tibialis anterior muscle of normal subjects. The estimation of the average PPR range of non-diseased muscles in non-fatigued conditions is important since it can be used as a reference to establish a comparison with PPR ranges from muscles suffering some disorder or from fibres that are fatigued. Other aspects of the PPR, as its sensitivity with volume conductor parameters or to what extent changes in the SFAP PPR reflects changes in IAP spatial profile are also examined. We found that the PPR of experimental SFAPs ranges from 0.3 to 2.5 in all subjects and that all PPR histograms contain a well-defined single peak around the PPR value 1.0. [ABSTRACT FROM AUTHOR]

Abstract: The present paper describes a global procedure for estimating all the synthesis parameters that generate a single fibre action potential (SFAP) in the Dimitrov–Dimitrova (D–D) convolutional model. We call this inverse problem Identification Procedure, and it is presented in two parts, this paper being the second. The procedure incorporates the candidate pair (CP) method developed in Part I, which provides the values of radial distance r and fibre diameter d of the simulated SFAP that best matches a potential under study. The CP-method required prior knowledge of all the excitation parameters. However, since the Identification Procedure makes no assumption about the excitation, multiple combinations of the synthesis parameters result in very similar SFAPs whose shape is close the signal under study. Analysis of the possible combinations reveals that r and d can be modelled as two jointly Gaussian random variables. The interest of the Identification Procedure is that, for a certain SFAP, it provides estimates of r and d, along with estimates of different parameters that determine the IAP waveform. Moreover, the procedure is able to determine the degree of error that accompanies the estimation of r and d. [Copyright &y& Elsevier]

The Dimitrov–Dimitrova (D–D) model generates a single fibre action potential (SFAP) as the convolution of an excitation function and a filter impulse function. We propose a method to estimate the parameters involved in these functions from a SFAP waveform (inverse problem) and call it Identification Procedure. The Identification Procedure comprises two parts. The present paper, Part I, is centred on the estimation of the radial distance, r, and the fibre diameter, d, of a given SFAP. To this end we develop a technique which we call the candidate pair method (CP-method), and we test it on fibrillation potentials (FPs) as experimental data. We found that the D–D model cannot synthesize all the SFAP waveforms observed in the experimental recordings, but in the cases where it can, the CP-method then provides the values of r and d that are more likely to have synthesized the SFAP. Having a method that provides information about the fibre diameter straight from a SFAP waveform is very desirable as this parameter has clinical and physiological relevance. Moreover, the CP-method plays a major role in the Identification Procedure that is carried out in Part II. [Copyright &y& Elsevier]

Introduction: We compared the recovery of muscle electrical properties after maximal voluntary contractions (MVCs) of 1 and 3 min duration by examining separately the first and second phases of the muscle compound action potential (M wave).Methods: M waves were evoked by supramaximal single shocks to the femoral nerve throughout the 30-min recovery following 1-min and 3-min MVCs. The amplitude, duration, and area of the M-wave first and second phases, along with peak-to-peak amplitude and total area, were measured from the knee extensors.Results: (1) The amplitude of the M-wave first phase increased to the same extent (and had the same time course of recovery) after the 1 and 3-min MVCs, whereas the amplitude of the second phase increased more markedly after the 1-min than after the 3-min MVC (P < 0.05). (2) The first phase remained enlarged for 2 min after exercise, whereas the augmentation of the second phase only lasted for 30 s. (3) After 30 min of recovery, the amplitude, area, and duration of both the first and second phases were decreased compared to control values (P < 0.05).Conclusions: The similar enlargement of the M-wave first phase after the 1 and 3-min MVCs suggests that the extracellular K+ concentration attained after these contractions was similar. The mechanisms responsible for the long-term decreases in M-wave amplitude and duration are unknown at present, but are likely due to a decrease in the amplitude of individual transmembrane potentials and an increase in conduction velocity. [ABSTRACT FROM AUTHOR]

Abstract: 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 40mm, 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. [Copyright &y& Elsevier]

Introduction: In the compound muscle action potential (M wave) recorded using the belly-tendon configuration, the contribution of the tendon electrode is assumed to be negligible compared to the belly electrode. We tested this assumption by placing the reference electrode at a distant (contralateral) site, which allowed separate recording of the belly and tendon contributions.Methods: M waves were recorded at multiple selected sites over the right quadriceps heads and lower leg using two different locations for the reference electrode: the ipsilateral (right) and contralateral (left) patellar tendon. The general parameters of the M wave (amplitude, area, duration, latency, and frequency) were measured.Results: (1) The tendon potential had a small amplitude (<30%) compared to the belly potential; (2) Changing the reference electrode from the ipsilateral to the contralateral patella produced moderate changes in the M wave recorded over the innervation zone, these changes affecting significantly the amplitude of the M-wave second phase (p = 0.006); (3) Using the contralateral reference system allowed recording of short-latency components occurring immediately after the stimulus artefact, which had the same latency and amplitude (p = 0.18 and 0.25, respectively) at all recording sites over the leg.Conclusions: The potential recorded at the "tendon" site after femoral nerve stimulation is small (compared to the belly potential), but not negligible, and makes a significant contribution to the second phase of belly-tendon M wave. Adopting a distant (contralateral) reference allowed recording of far-field components that may aid in the understanding of the electrical formation of the M wave. [ABSTRACT FROM AUTHOR]

Abstract: In situ recording of the intracellular action potential (IAP) of human muscle fibres is not yet feasible, and consequently, knowledge about certain IAP characteristics of these IAPs is still limited. The ratio between the amplitudes of the second and first phases (the so-called peak-to-peak ratio, PPR) of a single fibre action potential (SFAP) is known to be closely related to the IAP profile. The PPR of experimentally recorded SFAPs has been found to be largely independent of changes in the fibre-to-electrode (radial) distance. The main goal of this paper is to analyze the effect of changes in different aspects of the IAP spike on the relationship between PPR and radial distance. Based on this analysis, we hypothesize about the characteristics of IAPs obtained experimentally. It was found that the sensitivity of the SFAP PPR to changes in radial distance is essentially governed by the duration of the IAP spike. Assuming that, for mammals, the duration of the IAP rising phase lies within the range 0.2–0.4ms, we tentatively suggest that the duration of the IAP spike should be over approximately 0.75ms, with the shape of the spike strongly asymmetric. These IAP characteristics broadly coincide with those observed in mammal IAPs. [Copyright &y& Elsevier]

Abstract: In a series of previous works we studied the ratio between the amplitudes of the second and first phases (the peak-to-peak ratio) of single fibre action potential (SFAPs) using the Dimitrov–Dimitrova SFAP convolutional model as a reference. From experimental potentials extracted from both healthy and diseased muscles, we determined typical peak-to-peak ratio (PPR) values and ranges for both normal and pathological conditions. In addition, we investigated the changes observed in the PPR of consecutive potentials recorded at different fibre-to-electrode distances. However, our results were not conclusive due to insufficient data. The objective of the present work was to obtain a more concrete description of the relation between PPR and radial distance. To this end, we recorded 135 sets of consecutive SFAPs from the m. tibialis anterior of four normal subjects. The needle was intentionally moved whilst recording each SFAP set. We found that PPR was largely independent of small changes in electrode position when the electrode was close to the fibre and sufficiently far from the neuromuscular and/or fibre–tendon junctions. In the discussion, we provide evidence that this result is in agreement with the generation of extracellular potentials considering the spatial extension of the intracellular action potential (IAP) along the fibre. [Copyright &y& Elsevier]

Abstract: In situ recording of the intracellular action potential (IAP) of human muscle fibres is not yet possible, and consequently, knowledge concerning certain IAP characteristics is still limited. According to the core-conductor theory, close to a fibre, a single fibre action potential (SFAP) can be assumed to be proportional to the IAP second derivative. Thus, we might expect to be able to derive some characteristics of the IAP, such as the duration of its spike, from the SFAP waveform. However, SFAP properties not only depend on the IAP shape but also on the fibre-to-electrode (radial) distance and other physiological properties of the fibre. In this paper we, first, propose an SFAP parameter (the negative phase duration, NPD) appropriate for estimating the IAP spike duration and, second, show that this parameter is largely independent of changes in radial distance and muscle fibre propagation velocity. Estimation of the IAP spike duration from a direct measurement taken from the SFAP waveform provides a possible way to enhance the accuracy of SFAP models. Because IAP spike duration is known to be sensitive to the effects of fatigue and calcium accumulation, the proposed SFAP parameter, the NPD, has potential value in electrodiagnosis and as an indicator of IAP profile changes due to peripheral fatigue. [Copyright &y& Elsevier]

Abstract: Using the core-conductor theory, a single fibre action potential (SFAP) can be expressed as the convolution of a biolectrical source and a weight function. In the Dimitrov–Dimitrova (D–D) SFAP convolutional model, the first temporal derivative of the intracellular action potential (IAP) is used as the source. The present work evaluates the relationship between the SFAP peak-to-peak amplitude (V pp) and peak-to-peak interval (rise-time, RT) at different fibre-to-electrode distances using simulated signals obtained by the D–D model as well as real recordings. With a single fibre electrode, we recorded 63 sets of consecutive SFAPs from the m. tibialis anterior of four normal subjects. The needle was intentionally moved whilst recording each SFAP set. We used the observed changes in RT and V pp within each SFAP set as a point of reference with which to evaluate how closely the relationship between RT and V pp provided by the D–D model reflects real data. We found that half of the recorded SFAP sets had rise-times higher than those generated by the D–D model. We also showed the influence of the IAP spatial length on the sensitivity of RT and V pp with radial distance. The study reveals some inaccuracies in simulated SFAPs whose origin might be related to the assumptions made in the core-conductor theory. [ABSTRACT FROM AUTHOR]

Abstract: Objective: This study aims to check the theoretical predictions of the small changes in the rise-time of human single muscle fibre action potentials (SFAPs) experimentally when the recording electrode is in the proximity of the fibre. Methods: Using a single-fibre electrode, 93 sets of consecutive SFAPs were recorded under needle movement from the m. tibialis anterior of four normal subjects. In each SFAP set, the excitation source was assumed to remain unchanged and thus only the fibre–electrode distance varied between consecutive discharges. The changes of SFAP rise-time with peak-to-peak amplitude (V pp) within each SFAP set were used to study the relationship between rise-time and radial distance. Results: SFAP rise-time was relatively unchanged when V pp was greater than 1.0–1.5mV. For V pp smaller than 1.0–1.5mV, the sensitivity of rise-time to changes in radial distance increased considerably. Conclusions: Although the SFAP rise-time is practically not influenced by changes in radial distance in the proximity of muscle fibres, future work is necessary to confirm the hypothesis of using the rise-time for estimating changes in the intracellular action potential (IAP) depolarisation phase. Significance: The studies are useful in diagnostics and following the results of treatment in patients suffering from muscular dystrophy or myopathy, they are known to be accompanied by severe changes in the duration of IAP. [Copyright &y& Elsevier]

Objective: We studied movement-related electroencephalographic oscillatory changes in the alpha and beta range during a sequence of two movements in 7 healthy volunteers, in order to investigate the relationship between these changes and each component in the sequence.Methods: The sequence consisted of a wrist active extension-passive flexion followed by a first and second finger pincer. A total of 10.5 s sweeps were recorded using the level of surface electromyographic (EMG) activity in wrist extensors as trigger, including a 7.5 s pre-stimulus. The sweeps were also realigned manually offline using as trigger the end of the first EMG burst, or the beginning of the second movement. An index of the changes in non-phase-locked energy in the 7–37 Hz range was obtained by averaging single-sweep time-frequency transforms.Results: The duration of each of the movements in the sequence and the relationship between them were compatible with the use of two different motor programmes in the sequence. In the beta band, a decrease in energy (event-related desynchronisation, ERD) began 1.5 s before the onset of the first movement, and was sustained until the end of the second movement. No energy increases were observed until the end of the second movement. In the alpha band, the ERD began 0.5 seconds before the first movement and was sustained throughout the recording.Conclusion: These findings suggest that the beta-event-related synchronisation is related to the end of the whole motor process, and not to the end of each motor programme. [Copyright &y& Elsevier]