Javier Navallas

Motor unit action potential duration, I: variability of manual and automatic measurements.

The aim of this work is to analyze the variability in manual measurements of motor unit action potential (MUAP) duration and to evaluate the effectiveness of well-known algorithms for automatic measurement. Two electromyographists carried out three independent duration measurements of a set of 240 MUAPs. The intraexaminer and interexaminer variabilities were analyzed by means of the Gage Reproducibility and Repeatability method. The mean of the three closest manually marked positions was considered the gold standard of the duration markers positions (GSP). The results of four well-known automatic methods for estimating MUAP duration were compared to the GSP. Manual measurements of duration showed a lot of variability, with the combined intraoperator and interoperator variability greater than 30%. The greatest difference between manual positions was 11.2 ms. The mean differences between the GSP and those obtained with the four automatic methods ranged between 0.6 and 8.5 ms. Both manual and automatic measurements of MUAP duration show a high degree of variability. More precise methods are needed to improve the accuracy and reliability of the estimates of this parameter.

Estimating the duration of intracellular action potentials in muscle fibres from single-fibre extracellular potentials.

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.

Studying motor end-plate topography by means of scanning-electromyography

OBJECTIVE This study investigates the anatomical and physiological causes of the temporal jumps between motor unit fractions observed in scanning-EMG recordings; i.e., temporal gaps in the distribution of the electrical activity of the motor unit within a transverse section of its territory. METHODS Scanning-EMG signals were extracted from the biceps brachii. Spatial and temporal aspects of motor unit fractions were characterized. Simulated motor unit potentials (MUPs) were generated under different distributions of the anatomical and physiological parameters. Complexity and latency of the MUPs were calculated. RESULTS To obtain realistic simulated MUPs it was necessary to reduce the width of the motor end-plate zone and variability of the muscle fiber conduction velocity compared to the overall data available for the entire muscle. To simulate temporal jumps between fractions, motor end-plates must be distributed in sub-bands of narrow width with different mean positions, innervated by separate axonal branches. CONCLUSIONS Differences in the motor end-plates positions, axonal branches lengths and conduction velocities explain the temporal jumps between fractions. SIGNIFICANCE Characterization of the motor unit fractions, observable by means of scanning-EMG, can be used to obtain valuable information about the motor end-plate topography and branching pattern of the axon innervating the motor unit.

Motor unit action potential duration, II: a new automatic measurement method based on the wavelet transform.

The aim of this work is 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 to the gold standard of duration marker positions obtained by manual measurement. The new method was also compared to 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.

Optimizing interoperability between video-oculographic and electromyographic systems.

A new system is presented that enhances the interoperability between a video-oculographic (VOG) system for mouse movement control and an electromyographic (EMG) system for mouse click detection. The proposed VOG-EMG system combines gaze and muscle information to minimize the number of undesired clicks due to involuntary activations and environmental noise. We tested the system with 24 subjects, comparing three different configurations: one in which the VOG and EMG systems worked independently and two in which we used VOG gaze information to improve the EMG click detection. Results show that the number of false-positive click detections can be reduced when VOG and EMG information is combined. In addition, the third configuration, including extra processing, can reduce the activation delay produced because of the combined use of the VOG and EMG systems. The new VOG-EMG system is meant to be used in noisy environments in which the number of false clicks may impeach a reliable human-computer interaction.

Averaging methods for extracting representative waveforms from motor unit action potential trains

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.

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 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.

Effects of changes in the shape of the intracellular action potential on the peak-to-peak ratio of single muscle fibre potentials

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.

A muscle architecture model offering control over motor unit fiber density distributions

The aim of this study was to develop a muscle architecture model able to account for the observed distributions of innervation ratios and fiber densities of different types of motor units in a muscle. A model algorithm is proposed and mathematically analyzed in order to obtain an inverse procedure that allows, by modification of input parameters, control over the output distributions of motor unit fiber densities. The model’s performance was tested with independent data from a glycogen depletion study of the medial gastrocnemius of the rat. Results show that the model accurately reproduces the observed physiological distributions of innervation ratios and fiber densities and their relationships. The reliability and accuracy of the new muscle architecture model developed here can provide more accurate models for the simulation of different electromyographic signals.

Comparative evaluation of motor unit architecture models

We present a statistical evaluation and comparison of the simulation outcomes of nine different motor unit architecture modeling approaches, which derive from combinations of four motor unit territory placement algorithms with two innervation pattern algorithms (one of the combinations allows for a double approach). We test how well the outcomes of these models agree with well-established physiological principles. Our results show that algorithms based on independent and uniformly distributed territory placement always lead to an unwanted edge effect consisting in a decay of the number of overlapping motor unit territories toward the edge of the muscle cross section. This edge effect interferes with the properties of the simulated motor units. On the other hand, controlled placement of motor units, so as to minimize the spatial variance of muscle fiber density (MFD), can result in simulated motor units which better reflect the empirical evidence currently available about motor unit architecture.