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
Evaluation of the electromyography test for the analysis of the aerobic-anaerobic transition in elite cyclists during incremental exercise
(MDPI, 2019) Latasa Zudaire, Iban Alexander; Córdova Martínez, Alfredo; Quintana Ortí, Gregorio; Lavilla Oiz, Ana; Navallas Irujo, Javier; Rodríguez Falces, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
(1) Background: The aim of this study was to investigate the validity and reliability of surface electromyography (EMG) for automatic detection of the aerobic and anaerobic thresholds during an incremental continuous cycling test using 1 min exercise periods in elite cyclists. (2) Methods: Sixteen well-trained cyclists completed an incremental exercise test (25 W/1 min) to exhaustion. Surface bipolar EMG signals were recorded from the vastus lateralis, vastus medialis, biceps femoris, and gluteus maximus, and the root mean square (RMS) were assessed. The multi-segment linear regression method was used to calculate the first and second EMG thresholds (EMG(T1) and EMG(T2)). During the test, gas exchange data were collected to determine the first and second ventilatory thresholds (VT1 and VT2). (3) Results: Two breakpoints (thresholds) were identified in the RMS EMG vs. time curve for all muscles in 75% of participants. The two breakpoints, EMG(T1) and EMG(T2) , were detected at around 70%-80% and 90%-95% of VO2MAX, respectively. No significant differences were found between the means of VT(1 )and EMGT(1) for the vastii and biceps femoris muscles (p > 0.05). There were no significant differences between means of EMG(T2) and VT2 (p > 0.05). (4) Conclusions: It is concluded that the multi-segment linear regression algorithm is a valid non-invasive method for analyzing the aerobic-anaerobic transition during incremental tests with 1 min stage durations.
Open Access
End-of-fiber signals strongly influence the first and second phases of the M wave in the vastus lateralis: implications for the study of muscle excitability
(Frontiers Media, 2018) Rodríguez Falces, Javier; Place, Nicolas; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
It has been recurrently observed that, for compound muscle action potentials (M wave) recorded over the innervation zone of the vastus lateralis, the descending portion of the first phase generally shows an 'inflection' or 'shoulder'. We sought to clarify the electrical origin of this shoulder-like feature and examine its implications. M waves evoked by maximal single shocks to the femoral nerve were recorded in monopolar and bipolar configurations from 126 individuals using classical (10-mm recording diameter, 20-mm inter-electrode distance) electrodes and from eight individuals using small electrodes arranged in a linear array. The changes of the M-wave waveform at different positions along the muscle fibers' direction were examined. The shoulder was identified more frequently in monopolar (97%) than in bipolar (46%) M waves. The shoulder of M waves recorded at different distances from the innervation zone had the same latency. Furthermore, the shoulder of the M wave recorded over the innervation zone coincided in latency with the positive peak of that recorded beyond the muscle. The positive phase of the M wave detected 20 mm away from the innervation zone was essentially composed of non-propagating components. The shoulder-like feature in monopolar and bipolar M waves results from the termination of action potentials at the superficial aponeurosis of the vastus lateralis. We conclude that, only the amplitude of the first phase, and not the second, of M waves recorded monopolarly and/or bipolarly in close proximity to the innervation zone can be used reliably to monitor possible changes in muscle membrane excitability.
Open Access
Modelling fibrillation potentials: analysis of time parameters in the muscle intracellular action potential
(IEEE, 2007) Rodríguez Falces, Javier; Malanda Trigueros, Armando; Gila Useros, Luis; Rodríguez Carreño, Ignacio; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
A single fibre action potential (SFAP) can be modelled as the convolution of a biolectrical source and a filter impulse response. In the Dimitrov-Dimitrova (D-D) convolutional model, the first temporal derivative of the intracellular action potential (IAP) is used as the source, and T spl is a time parameter related to the duration of the IAP waveform. Our work is centred on the relation between Tspl and the main spike duration (MSD), defined as the time interval between the first and third phases of the SFAP. We show that Tspl essentially determines the MSD parameter. As experimental data, we used fibrillation potentials (FPs) of two different muscles to study the D-D model. We found that T spl should have a certain statistical variability in order to explain the variability in the MSD of our FPs. In addition, we present a method to estimate the T spl values corresponding to a given SFAP from its measured MSD.
Open Access
A masked least-squares smoothing procedure for artifact reduction in scanning-EMG recordings
(Springer, 2018) Corera Orzanco, Íñigo; Eciolaza Ferrando, Adrián; Rubio Zamora, Oliver; Malanda Trigueros, Armando; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica, Electrónica y de Comunicación; Ingeniaritza Elektrikoa, Elektronikoaren eta Telekomunikazio Ingeniaritzaren
Scanning-EMG is an electrophysiological technique in which the electrical activity of the motor unit is recorded at multiple points along a corridor crossing the motor unit territory. Correct analysis of the scanning-EMG signal requires prior elimination of interference from nearby motor units. Although the traditional processing based on the median filtering is effective in removing such interference, it distorts the physiological waveform of the scanning-EMG signal. In this study, we describe a new scanning-EMG signal processing algorithm that preserves the physiological signal waveform while effectively removing interference from other motor units. To obtain a cleaned-up version of the scanning signal, the masked least-squares smoothing (MLSS) algorithm recalculates and replaces each sample value of the signal using a least-squares smoothing in the spatial dimension, taking into account the information of only those samples that are not contaminated with activity of other motor units. The performance of the new algorithm with simulated scanning-EMG signals is studied and compared with the performance of the median algorithm and tested with real scanning signals. Results show that the MLSS algorithm distorts the waveform of the scanning-EMG signal much less than the median algorithm (approximately 3.5 dB gain), being at the same time very effective at removing interference components.
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 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.
Open Access
Motor unit action potential duration, II: a new automatic measurement method based on the wavelet transform
(Lippincott, Williams & Wilkins, 2007) Rodríguez Carreño, Ignacio; Gila Useros, Luis; Malanda Trigueros, Armando; García Gurtubay, Ignacio; Mallor Giménez, Fermín; Gómez Elvira, Sagrario; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Estadística e Investigación Operativa; Ingeniaritza Elektrikoa eta Elektronikoa; Estatistika eta Ikerketa Operatiboa
To present and evaluate a new algorithm, based on the wavelet transform, for the automatic measurement of motor unit action potential (MUAP) duration. A total of 240 MUAPs were studied. The waveform of each MUAP was wavelet-transformed, and the start and end points were estimated by regarding the maxima and minima points in a particular scale of the wavelet transform. The results of the new method were compared with the gold standard of duration marker positions obtained by manual measurement. The new method was also compared with a conventional algorithm, which we had found to be best in a previous comparative study. To evaluate the new method against manual measurements, the dispersion of automatic and manual duration markers were analyzed in a set of 19 repeatedly recorded MUAPs. The differences between the new algorithm’s marker positions and the gold standard of duration marker positions were smaller than those observed with the conventional method. The dispersion of the new algorithm’s marker positions was slightly less than that of the manual one. Our new method for automatic measurement of MUAP duration is more accurate than other available algorithms and more consistent than manual measurements.
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
Correlation between discharge timings of pairs of motor units reveals the presence but not the proportion of common synaptic input to motor neurons
(American Physiological Society, 2017) Rodríguez Falces, Javier; Negro, Francesco; Farina, Dario; Ingeniería Eléctrica y Electrónica; Ingeniaritza Elektrikoa eta Elektronikoa
We investigated whether correlation measures derived from pairs of motor unit (MU) spike trains are reliable indicators of the degree of common synaptic input to motor neurons. Several 50-s isometric contractions of the biceps brachii muscle were performed at different target forces ranging from 10 to 30% of the maximal voluntary contraction relying on force feedback. Forty-eight pairs of MUs were examined at various force levels. Motor unit synchrony was assessed by cross-correlation analysis using three indexes: the output correlation as the peak of the cross-histogram (ρ) and the number of synchronous spikes per second (CIS) and per trigger (E). Individual analysis of MU pairs revealed that ρ, CIS, and E were most often positively associated with discharge rate (87, 85, and 76% of the MU pairs, respectively) and negatively with interspike interval variability (69, 65, and 62% of the MU pairs, respectively). Moreover, the behavior of synchronization indexes with discharge rate (and interspike interval variability) varied greatly among the MU pairs. These results were consistent with theoretical predictions, which showed that the output correlation between pairs of spike trains depends on the statistics of the input current and motor neuron intrinsic properties that differ for different motor neuron pairs. In conclusion, the synchronization between MU firing trains is necessarily caused by the (functional) common input to motor neurons, but it is not possible to infer the degree of shared common input to a pair of motor neurons on the basis of correlation measures of their output spike trains.