Motor unit conduction velocity distribution estimation: assessment of two short-term processing methods

Central and peripheral components of exercise-related fatigability in myotonic dystrophy type 1

BACKGROUND

Fatigue frequently occurs in myotonic dystrophy type 1 (DM1), but its pathophysiology remains unclear. This study assessed central and peripheral components of exercise-related fatigability in patients with DM1, compared to controls.

METHODS

Examinations were performed before and after a contraction of the abductor digiti minimi (ADM) muscle sustained for 45 s at 60% of maximal voluntary contraction (MVC). Myoelectric activity was recorded using high spatial resolution surface EMG during twitch stimulations and MVC and was characterized by root mean square, mean power frequency (MPF), and muscle fiber conduction velocity (MFCV). Peripheral nerve excitability was assessed by stimulating the ulnar nerve at the wrist with ADM recordings. Motor cortex excitability testing to transcranial magnetic stimulation included measures of intracortical facilitation and inhibition of motor evoked potentials (MEPs) in ADM muscle.

RESULTS

At baseline, patients with DM1 showed altered peripheral nerve and cortical excitability (reduced intracortical facilitation) associated with impaired myoelectric properties. During the fatiguing exercise, the force remained stable, while MPF and MFCV decreased in both DM1 and control groups. After exercise, only refractoriness was reduced in patients with DM1, whereas controls showed marked neuromuscular and cortical changes.

CONCLUSION

Patients with DM1 showed altered excitability of various cortical and neuromuscular components at baseline. However, most of excitability parameters did not vary after exercise in patients with DM1, in contrast to controls. This suggests that excitability properties, frankly altered at baseline, were not prone to be affected further after exercise in patients with DM1.

Simulation of Surface EMG Signals for a Multilayer Volume Conductor With Triangular Model of the Muscle Tissue

This study analytically describes surface electromyogram (sEMG) signals generated by a model of a triangular muscle, i.e., a muscle with fibers arranged in a fan shape. Examples of triangular muscles in the human body are the deltoid, the pectoralis major, the trapezius, the adductor pollicis. A model of triangular muscle is proposed. It is a sector of a cylindrical volume conductor (with the fibers directed along the radial coordinate) bounded at the muscle/fat interface. The muscle conductivity tensor reflects the fan anisotropy. Edge effects have been neglected. A solution of the nonspace invariant problem for a triangular muscle is provided in the Fourier domain. An approximate analytical solution for a two plane layer volume conductor model is obtained by introducing a homogeneous layer (modeling the fat) over the triangular muscle. The results are implemented in a complete sEMG generation model (including the finite length of the fibers), simulating single fiber action potentials. The model is not space invariant due to the changes of the volume conductor along the direction of action potential propagation. Thus the detected potentials at the skin surface change shape as they propagate. This determines problems in the extraction and interpretation of parameters. As a representative example of application of the simulation model, the influence of the inhomogeneity of the volume conductor in conduction velocity (CV) estimation is addressed (for two channels; maximum likelihood and reference point methods). Different fiber depths, electrode placements and small misalignments of the detection system with respect to the fiber have been simulated. The error in CV estimation is large when the depth of the fiber increases, when the detection system is not aligned with the fiber and close to the innervation point and to the tendons.

Motor unit conduction velocity distribution estimation from evoked motor responses

Action potentials travel along the muscle fibers with a specific conduction velocity that depends on their structural and functional properties. Only the estimation of muscle conduction velocity distribution (MCVD) may be able to depict this propagation heterogeneity. Based on the method proposed by Cummins et al. (Electroenceph Clin Neurophysiol, 46:647-658, 1979) to estimate nerve conduction velocity distribution (NCVD), the present paper proposes a method that modifies the Cummins' approach to make it suitable for MCVD estimation from electrically evoked motor responses. The MCVD estimation algorithm was first assessed by means of simulated signals in order to control all signal features during the optimization process. Simulations showed that estimated distributions were very close to the true ones when taking into account the specificities of the muscle action potential, due to its generation and extinction (MSE divided by 5 on distribution standard deviation). This method was then applied to real signals. Elicited motor responses were recorded on the biceps brachii of healthy subjects either during repeated maximal stimulations at 20 Hz or during increasing intensity stimulations at 1 Hz. MCVD estimates were used to analyze fatigue and motor unit recruitment processes, respectively.

Tracking motor unit action potentials in the tibialis anterior during fatigue

New surface electromyogram (SEMG) techniques offer the potential to advance knowledge of healthy and diseased motor units. Conduction velocity (CV) estimates, obtained from indwelling electrodes, may provide diagnostic information, but the standard method of CV estimation from SEMG may be of only limited value. We developed a motor unit (MU) tracking algorithm to extract motor unit conduction velocity (MUCV) and motor unit action potential (MUAP) amplitude estimates from SEMG. The technique is designed to provide a noninvasive means of accessing fatigue and recruitment behavior of individual MUs. We have applied this MU tracking algorithm to SEMG data recorded during isometric fatiguing contractions of the tibialis anterior (TA) muscle in nine healthy subjects, at 30%-40% maximum voluntary contraction (MVC). The results reveal that MUCVs and MUAP amplitudes of individual MUs can be estimated and tracked across time. Time-related changes in the MU population may also be monitored. Thus, the SEMG technique employed provides insight into the behavior of the underlying muscle at the MU level by noninvasive means.

Clinical applications of surface electromyography in neuromuscular disorders

Surface electromyography (SEMG) is still rarely used in clinical settings for the detection and analysis of myoelectric signals. The electromyographic signal detected on the skin surface includes information from a greater proportion of the muscle of interest than conventional clinical EMG, acquired using needle electrodes. SEMG is therefore more representative than the localised, and thus very selective needle EMG signals currently used. However, both reliability and interpretation of surface EMG need to be questioned. This review looks at the studies concerned with the characterisation of neuromuscular pathologies using EMG parameters. After introducing principles and limitations of surface EMG, an overview of the main results obtained in clinical settings is presented and discussed. There is a particular focus on high spatial resolution surface EMG as it is currently the best compromise between the selectivity of needle EMG and the representative nature of classical SEMG. Several perspectives are proposed that underline the fact that surface EMG is an evolving discipline and should be worthy of a place in routine clinical examinations.

A Technique to Track Individual Motor Unit Action Potentials in Surface EMG by Monitoring Their Conduction Velocities and Amplitudes

The speed of propagation of an action potential along a muscle fiber, its conduction velocity (CV), can be used as an indication of the physiological or pathological state of the muscle fiber membrane. The motor unit action potential (MUAP), the waveform resulting from the spatial and temporal summation of the individual muscle fiber action potentials of that motor unit (MU), propagates with a speed referred to as the motor unit conduction velocity (MUCV). This paper introduces a new algorithm, the MU tracking algorithm, which estimates MUCVs and MUAP amplitudes for individual MUs in a localized MU population using SEMG signals. By tracking these values across time, the electrical activity of the localized MU pool can be monitored. An assessment of the performance of the algorithm has been achieved using simulated SEMG signals. It is concluded that this analysis technique enhances the suitability of SEMG for clinical applications and points toward a future of noninvasive diagnosis and assessment of neuromuscular disorders.

Simultaneous estimation of muscle fibre conduction velocity and muscle fibre orientation using 2D multichannel surface electromyogram

The paper presents a new approach for simultaneous estimation of muscle fibre conduction velocity (MFCV) and muscle fibre orientation (MFO) for motor units (MUs) in two-dimensional (2D) multichannel surface electromyography recordings. This is an important tool for detecting changes and abnormalities in muscle function and structure. In addition, simultaneous estimation of MFO and MFCV avoids the necessity of manual electrode alignment. The proposed method detected propagating MU action potentials (MUAPs) in a running time window as moving components in amplitude maps. Thereafter, estimations were obtained by fitting a three-dimensional function to these maps. The performance was evaluated using synthetic MU signals at 10 dB SNR and authentic biceps brachii measurements. Results demonstrated MFCV and MFO estimates with standard deviations of less than 0.05 m s(-1) and 1 degrees for simulated signals, and less than 0.2 m s(-1) and 4 degrees for experimental data. However, standard deviations as low as 0.12 m s(-1) and 1.6 degrees from real signals were demonstrated. It was concluded that the method performs as well as, or better than, linear array multichannel methods when individual propagating MUAPs can be identified, even if electrodes are not aligned with fibre direction.

The effects of bipolar electrode montage on conduction velocity estimation from the surface electromyogram

This study examines the influence of the bipolar electrode montage on conduction velocity (CV) estimation. Electrode montage refers to the combination of two parameters, the inter-electrode distance (IED), the distance between the two electrodes of a bipolar pair, and the inter-signal distance (ISD), the distance between two bipolar signals used to calculate CV. Data from the biceps brachii (BB) and tibialis anterior (TA) of healthy subjects are analysed. Two approaches are used for CV estimation. The first returns a single value per epoch. The second is based on finding velocity values from individual peaks in the signal and results in a peak velocity (PV) distribution being generated per epoch. It is concluded that CV estimation is significantly dependent on the choice of the (IED, ISD) electrode montage. The main results are that the electrode montage affects (1) the mean PV and CV estimates, (typically P < 0.001), (2) the degree of spatial variability, and (3) the width of the PV distributions. The combination of a small IED and large an ISD is recommended.

Methods for estimating muscle fibre conduction velocity from surface electromyographic signals

The review focuses on the methods currently available for estimating muscle fibre conduction velocity (CV) from surface electromyographic (EMG) signals. The basic concepts behind the issue of estimating CV from EMG signals are discussed. As the action potentials detected at the skin surface along the muscle fibres are, in practice, not equal in shape, the estimation of the delay of propagation (and thus of CV) is not a trivial task. Indeed, a strictly unique definition of delay does not apply in these cases. Methods for estimating CV can thus be seen as corresponding to specific definitions of the delay of propagation between signals of unequal shape. The most commonly used methods for CV estimation are then reviewed. Together with classic methods, recent approaches are presented. The techniques are described with common notations to underline their relationships and to highlight when an approach is a generalisation of a previous one or when it is based on new concepts. The review identifies the difficulties of CV estimation and underlines the issues that should be considered by the investigator when selecting a particular method and detection system for assessing muscle fibre CV. The many open issues in CV estimation are also presented.

Simulation analysis of the ability to estimate motor unit propagation velocity non-invasively by different two-channel methods and types of multi-electrodes

Ability to estimate motor unit propagation velocity correctly using different two-channel methods for delay estimation and different non-invasive spatial filters was analysed by simulation. It was established that longitudinal double difference electrodes could be not a better choice than simple bipolar parallel electrodes. Spatial filtration with a new multi-electrode (performing difference between signals detected by two transversal double difference electrodes positioned along the muscle fibres) promises to give the best estimate. Delay estimation between reference points is more preferable than that based on the cross-correlation technique, which is considerably sensitive to the fundamental properties of the muscle fibre extracellular fields. Preliminary averaging and approximation of the appropriate parts of the signals around chosen reference points could reduce the larger noise sensitivity and the effects of local tissue inhomogeneities as well as eliminate the sampling problem. A correct estimate of the propagation velocity could be impossible, even in the case of not very deep motor units (15 or 10 mm, depending on the spatial filter used) with relatively long (about 120 mm) muscle fibres. In the case of fibres with asymmetrical location of the end-plates in respect to the fibre ends, the propagation velocity estimates could be additionally biased above the longer semilength of the motor unit fibres.