Open Access
Validation of the filling factor index to study the filling process of the sEMG signal in the quadriceps
(Elsevier, 2023) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Mariscal Aguilar, Cristina; Niazi, Imran Khan; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
Introduction: 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. Methods: 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). Results: (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. Conclusions: 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.
Open Access
M-wave changes caused by brief voluntary and stimulated isometric contractions
(Springer, 2023) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Navallas Irujo, Javier; Place, Nicolas; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Introduction Under isometric conditions, the increase in muscle force is accompanied by a reduction in the fbers’ length. The efects of muscle shortening on the compound muscle action potential (M wave) have so far been investigated only by computer simulation. This study was undertaken to assess experimentally the M-wave changes caused by brief voluntary and stimulated isometric contractions. Methods Two diferent methods of inducing muscle shortening under isometric condition were adopted: (1) applying a brief (1 s) tetanic contraction and (2) performing brief voluntary contractions of diferent intensities. In both methods, supramaximal stimulation was applied to the brachial plexus and femoral nerves to evoke M waves. In the frst method, electrical stimulation (20 Hz) was delivered with the muscle at rest, whereas in the second, stimulation was applied while participants performed 5-s stepwise isometric contractions at 10, 20, 30, 40, 50, 60, 70, and 100% MVC. The amplitude and duration of the frst and second M-wave phases were computed. Results The main fndings were: (1) on application of tetanic stimulation, the amplitude of the M-wave frst phase decreased (~10%, P<0.05), that of the second phase increased (~50%, P<0.05), and the M-wave duration decreased (~20%, P<0.05) across the frst fve M waves of the tetanic train and then plateaued for the subsequent responses; (2) when superimposing a single electrical stimulus on muscle contractions of increasing forces, the amplitude of the M-wave frst phase decreased (~20%, P<0.05), that of the second phase increased (~30%, P<0.05), and M-wave duration decreased (~30%, P<0.05) as force was raised from 0 to 60–70% MVC force. Conclusions The present results will help to identify the adjustments in the M-wave profle caused by muscle shortening and also contribute to diferentiate these adjustments from those caused by muscle fatigue and/or changes in Na+–K+ pump activity.
Open Access
The filling factor of the sEMG signal at low contraction forces in the quadriceps muscles is influenced by the thickness of the subcutaneous layer
(Frontiers Media, 2023) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Mariscal Aguilar, Cristina; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Institute of Smart Cities - ISC; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
Introduction: It has been shown that, for male subjects, the sEMG activity at low contraction forces is normally 'pulsatile', i.e., formed by a few large-amplitude MUPs, coming from the most superficial motor units. The subcutaneous layer thickness, known to be greater in females than males, influences the electrode detection volume. Here, we investigated the influence of the subcutaneous layer thickness on the type of sEMG activity (pulsatile vs. continuous) at low contraction forces. Methods: Voluntary surface EMG signals were recorded from the quadriceps muscles of healthy males and females as force was gradually increased from 0% to 40% MVC. The sEMG filling process was examined by measuring the EMG filling factor, computed from the non-central moments of the rectified sEMG signal. Results: 1) The sEMG activity at low contraction forces was ¿continuous¿ in the VL, VM and RF of females, whereas this sEMG activity was ¿pulsatile¿ in the VL and VM of males. 2) The filling factor at low contraction forces was lower in males than females for the VL (p = 0.003) and VM (p = 0.002), but not for the RF (p = 0.54). 3) The subcutaneous layer was significantly thicker in females than males for the VL (p = 0.001), VM (p = 0.001), and RF (p = 0.003). 4) A significant correlation was found in the vastus muscles between the subcutaneous layer thickness and the filling factor (p < 0.05). Discussion: The present results indicate that the sEMG activity at low contraction forces in the female quadriceps muscles is 'continuous' due to the thick subcutaneous layer of these muscles, which impedes an accurate assessment of the sEMG filling process.
Open Access
Automatic jitter measurement in needle-detected motor unit potential trains
(Elsevier, 2022) Malanda Trigueros, Armando; Stashuk, Daniel W.; Navallas Irujo, Javier; Rodríguez Falces, Javier; Rodríguez Carreño, Ignacio; Valle, César; Garnés Camarena, Óscar; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
In an active motor unit (MU), the time intervals between the firings of its muscle fibers vary across successive MU activations. This variability is called jitter and is increased in pathological processes that affect the neuromuscular junctions or terminal axonal segments of MUs. Traditionally, jitter has been measured using single fiber electrodes (SFEs) and a difficult and subjective manual technique. SFEs are expensive and reused, implying a potential risk of patient infection; so, they are being gradually substituted by safer, disposable, concentric needle electrodes (CNEs). As CNEs are larger, voltage contributions from individual fibers of a MU are more difficult to detect, making jitter measurement more difficult. This paper presents an automatic method to estimate jitter from trains of motor unit potentials (MUPs), for both SFE and CNE records. For a MUP train, segments of MUPs generated by single muscle fibers (SF MUP segments) are found and jitter is measured between pairs of these segments. Segments whose estimated jitter values are not reliable, according to several SF MUP segment characteristics, are excluded. The method has been tested in several simulation studies that use mathematical models of muscle fiber potentials. The results are very satisfactory in terms of jitter estimation error (less than 10% in most of the cases studied) and mean number of valid jitter estimates obtained per simulated train (greater than 1.0 in many of the cases and less than 0.5 only in the most complicated). A preliminary study with real signals was also performed, using 19 MUP trains from 3 neuropathic patients. Jitter measurements obtained by the automatic method were compared with those extracted from a commercial system (Keypoint) and the edition and supervision of an expert electromyographer. From these measurements 63% were taken from equivalent interval pair sites within the time span of the MUP trains and, as such, were considered as compatible measurements. Differences in jitter of these compatible measurements were very low (mean value of 1.3 μs, mean of absolute differences of 2.97 μs, 25% and 75% percentile intervals of − 0.85 and 3.82 μs, respectively). Although new tests with larger number of real recordings are still required, the method seems promising for clinical practice.
Open Access
EMG probability density function: a new way to look at EMG signal filling from single motor unit potential to full interference pattern
(IEEE, 2023) Navallas Irujo, Javier; Eciolaza Ferrando, Adrián; Mariscal Aguilar, Cristina; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
An analytical derivation of the EMG signal's amplitude probability density function (EMG PDF) is presented and used to study how an EMG signal builds-up, or fills, as the degree of muscle contraction increases. The EMG PDF is found to change from a semi-degenerate distribution to a Laplacian-like distribution and finally to a Gaussian-like distribution. We present a measure, the EMG filling factor, to quantify the degree to which an EMG signal has been built-up. This factor is calculated from the ratio of two non-central moments of the rectified EMG signal. The curve of the EMG filling factor as a function of the mean rectified amplitude shows a progressive and mostly linear increase during early recruitment, and saturation is observed when the EMG signal distribution becomes approximately Gaussian. Having presented the analytical tools used to derive the EMG PDF, we demonstrate the usefulness of the EMG filling factor and curve in studies with both simulated signals and real signals obtained from the tibialis anterior muscle of 10 subjects. Both simulated and real EMG filling curves start within the 0.2 to 0.35 range and rapidly rise towards 0.5 (Laplacian) before stabilizing at around 0.637 (Gaussian). Filling curves for the real signals consistently followed this pattern (100% repeatability within trials in 100% of the subjects). The theory of EMG signal filling derived in this work provides (a) an analytically consistent derivation of the EMG PDF as a function of motor unit potentials and motor unit firing patterns; (b) an explanation of the change in the EMG PDF according to degree of muscle contraction; and (c) a way (the EMG filling factor) to quantify the degree to which an EMG signal has been built-up.
Open Access
The first and second phases of the muscle compound action potential in the thumb are differently affected by electrical stimulation trains
(American Physiological Society, 2024) Lanfranchi, Clément; Rodríguez Falces, Javier; Place, Nicolas; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoa eta Telekomunikazio Ingeniaritza; Institute of Smart Cities - ISC
Sarcolemmal membrane excitability is often evaluated by considering the peak-to-peak amplitude of the compound muscle action potential (M wave). However, the first and second M-wave phases represent distinct properties of the muscle action potential, which are differentially affected by sarcolemma properties and other factors such as muscle architecture. Contrasting with previous studies in which voluntary contractions have been used to induce muscle fatigue, we used repeated electrically induced tetanic contractions of the adductor pollicis muscle and assessed the kinetics of M-wave properties during the course of the contractions. Eighteen participants (24 ± 6 yr; means ± SD) underwent 30 electrically evoked tetanic contractions delivered at 30 Hz, each lasting 3 s with 1 s intervals. We recorded the amplitudes of the first and second M-wave phases for each stimulation. During the initial stimulation train, the first and second M-wave phases exhibited distinct kinetics. The first phase amplitude showed a rapid decrease to reach ~59% of its initial value (P < 0.001), whereas the second phase amplitude displayed an initial transient increase of ~19% (P ¼ 0.007). Within subsequent trains, both the first and second phase amplitudes consistently decreased as fatigue developed with a reduction during the last train reaching ~47% of its initial value (P < 0.001). Analyzing the first M wave of each stimulation train unveiled different kinetics for the first and second phases during the initial trains, but these distinctions disappeared as fatigue progressed. These findings underscore the interplay of factors affecting the M wave and emphasize the significance of separately scrutinizing its first and second phases when assessing membrane excitability adjustments during muscle contractions.
Open Access
Understanding EMG PDF changes with motor unit potential amplitudes, firing rates, and noise level through EMG filling curve analysis
(IEEE, 2024-08-30) Navallas Irujo, Javier; Mariscal Aguilar, Cristina; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoa eta Telekomunikazio Ingeniaritza
EMG filling curve characterizes the EMG filling process and EMG probability density function (PDF) shape change for the entire force range of a muscle.We aim to understand the relation between the physiological and recording variables, and the resulting EMG filling curves. We thereby present an analytical and simulation study to explain how the filling curve patterns relate to specific changes in the motor unit potential (MUP) waveforms and motor unit (MU) firing rates, the two main factors affecting the EMG PDF, but also to recording conditions in terms of noise level. We compare the analytical results with simulated cases verifying a perfect agreement with the analytical model. Finally, we present a set of real EMG filling curves with distinct patterns to explain the information about MUP amplitudes, MU firing rates, and noise level that these patterns provide in the light of the analytical study. Our findings reflect that the filling factor increases when firing rate increases or when newly recruited motor unit have potentials of smaller or equal amplitude than the former ones. On the other hand, the filling factor decreases when newly recruited potentials are larger in amplitude than the previous potentials. Filling curves are shown to be consistent under changes of the MUP waveform, and stretched under MUP amplitude scaling. Our findings also show how additive noise affects the filling curve and can even impede to obtain reliable information from the EMG PDF statistics.
Open Access
The contribution of the tendon electrode to M-wave characteristics in the biceps brachii, vastus lateralis and tibialis anterior
(Wiley, 2023) Rodríguez Falces, Javier; Etxaleku, Saioa; Trajano, Gabriel S.; Setuain Chourraut, Igor; Ciencias de la Salud; Ingeniería Eléctrica, Electrónica y de Comunicación; Osasun Zientziak; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
In some compound muscle action potentials (M waves) recorded using the belly–tendon configuration, the tendon electrode makes a noticeable contribution to the M wave. However, this finding has only been demonstrated in some hand and foot muscles. Here, we assessed the contribution of the tendon potential to the amplitude of the vastus lateralis, biceps brachii and tibialis anterior M waves, and we also examined the role of this tendon potential in the shoulder-like feature appearing in most M waves. M waves were recorded separately at the belly and tendon locations of the vastus lateralis, biceps brachii and tibialis anterior from 38 participants by placing the reference electrode at a distant (contralateral) site. The amplitude of the M waves and the latency of their peaks and shoulders were measured. In the vastus lateralis, the tendon potential was markedly smaller in amplitude (∼75%) compared to the belly M wave (P = 0.001), whereas for the biceps brachii and tibialis anterior, the tendon and belly potentials had comparable amplitudes. In the vastus lateralis, the tendon potential showed a small positive peak coinciding in latency with the shoulder of the belly–tendon M wave, whilst in the biceps brachii and tibialis anterior, the tendon potential showed a clear negative peak which coincided in latency with the shoulder. The tendon potential makes a significant contribution to the belly–tendon M waves of the biceps brachii and tibialis anterior muscles, but little contribution to the vastus lateralis M waves. The shoulder observed in the belly–tendon M wave of the vastus lateralis is caused by the belly potential, the shoulder in the biceps brachii M wave is generated by the tendon potential, whereas the shoulder in the tibialis anterior M wave is caused by both the tendon and belly potentials.
Embargo
Modeling the extracellular potential generated by a muscle fiber as the output signal of a convolutional system
(American Physiological Society, 2024-09-01) Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoa eta Telekomunikazio Ingeniaritza; Institute of Smart Cities - ISC
A central topic in Bioelectricity is the generation of the extracellular potential that results from the propagation of a transmembrane action potential along the muscle fiber. However, the way in which the extracellular potential is determined by the propagating action potential is difficult to describe, conceptualize, and visualize. Moreover, traditional quantitative approaches aimed at modeling extracellular potentials involve complex mathematical formulations, which do not allow students to visualize how the extracellular potential is generated around the active fiber. The present study is aimed at presenting a novel pedagogical approach to teaching the generation of extracellular potentials produced by muscle fibers based on the convolution operation. The effectiveness of this convolutional model was tested using a written exam and a satisfaction survey. Most students reported that a great advantage of this model was that it simplifies the problem by dividing it into three distinct components: 1) the input signal (associated with the action potential), 2) the impulse response (linked to the system formed by the fiber and the recording electrode), and 3) the output signal (the extracellular potential). Another key aspect of the present approach was that the input signal was represented by a sequence of electric dipoles, which allowed students to visualize the individual contribution of each dipole to the resulting extracellular potential. The results of the survey indicate that the combination of basic principles of electrical fields and intuitive graphical representations largely improves students' understanding of Bioelectricity concepts and enhances their motivation to complete their studies of biomedical engineering.
Open Access
The probability density function of the surface electromyogram and its dependence on contraction force in the vastus lateralis
(BMC, 2024-10-26) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Mariscal Aguilar, Cristina; Recalde Villamayor, Silvia; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoa eta Telekomunikazio Ingeniaritza; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Introduction: the probability density function (PDF) of the surface electromyogram (sEMG) depends on contraction force. This dependence, however, has so far been investigated by having the subject generate force at a few fixed percentages of MVC. Here, we examined how the shape of the sEMG PDF changes with contraction force when this force was gradually increased from zero. Methods: voluntary surface EMG signals were recorded from the vastus lateralis of healthy subjects as force was increased in a continuous manner vs. in a step-wise fashion. The sEMG filling process was examined by measuring the EMG filling factor, computed from the non-central moments of the rectified sEMG signal. Results: in 84% of the subjects, as contraction force increased from 0 to 10% MVC, the sEMG PDF shape oscillated back and forth between the semi-degenerate and the Gaussian distribution; the PDF–force relation varied greatly among subjects for forces between 0 and ~ 10% MVC, but this variability was largely reduced for forces above 10% MVC; the pooled analysis showed that, as contraction force gradually increased, the sEMG PDF evolved rapidly from the semi-degenerate towards the Laplacian distribution from 0 to 5% MVC, and then more slowly from the Laplacian towards the Gaussian distribution for higher forces. Conclusions: the study demonstrated that the dependence of the sEMG PDF shape on contraction force can only be reliably assessed by gradually increasing force from zero, and not by performing a few constant-force contractions. The study also showed that the PDF–force relation differed greatly among individuals for contraction forces below 10% MVC, but this variability was largely reduced when force increased above 10% MVC.