|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Laboratorio di Ingegneria del Sistema Neuromusculore, Dipartimento di Elettronica, Politecnico di Torino, Torino, Italy; and 2Center for Sensory-Motor Interaction (SMI), Department of Health Science and Technology, Aalborg University, Aalborg, Denmark
Submitted 16 November 2004; accepted in final form 6 February 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Multi-channel surface EMG recordings can be applied to analyze the membrane properties of muscle fibers both at global and single motor unit level. However, surface EMG has poor spatial selectivity that hinders the identification of single motor units, especially at medium/high contraction levels. To overcome this problem, we recently proposed the use of a visual feedback on multi-channel surface EMG signals to identify and follow the activity of single motor units (Farina et al. 2004b
). With this method, it was possible to analyze each action potential generated by a target motor unit in normal and ischemic conditions (Farina et al. 2004a
,b
). The possibility of identifying and recalling the same motor unit with surface EMG feedback allows the analysis of changes in motor unit conduction velocity modulated by sustained contraction, ischemia, or other conditions.
The aim of this study was to test the hypothesis that a sustained contraction determines changes in the conduction velocity of motor units, independently of their activation. Thus single motor units were analyzed with multi-channel surface EMG feedback before and after a low-force sustained contraction during which they were not activated. To induce appreciable changes in metabolite concentration during relatively short and low-force contractions, the muscle was studied in ischemic conditions. This resembled the situation at higher force contraction levels when the intramuscular pressure overcomes the systolic pressure and the muscle becomes ischemic.
| METHODS |
|---|
|
|
|---|
Ten healthy male subjects [age, 26.3 ± 2.1 (SD) yr; height, 1.82 ± 0.06 m; weight, 72.4 ± 6.9 kg) participated in the study. The local ethics committee approved the study, and all subjects signed an informed consent form before participation.
Surface EMG recordings
Surface EMG signals were detected with a two-dimensional array of 61 silver electrodes (1 mm diam, 3-mm interelectrode distance, 5 columns and 13 rows without the 4 corner electrodes; Fig. 1A) from the abductor pollicis brevis of the dominant hand. The small electrode size and interelectrode distance allowed high spatial selectivity (Reucher et al. 1987
). The EMG signals were amplified with four 16-channel EMG amplifiers (LISiN-Prima Biomedical & Sport, Treviso, Italy), band-pass filtered (3-dB bandwidth, 10500 Hz), sampled at 1,650 Hz, converted to digital form by 12-bit A/D converters, and displayed in real time as bipolar derivations to the subject.
|
General procedures
A custom-designed brace was used to measure abduction force (Fig. 1B). The subject's wrist was fixed in a padded wood support with the head of the thumb phalanx in touch with a load cell (model 8523-50N, Burster, Gernsbach, Germany). The force signal was amplified (Force Amplifier, MISO-II, LISiN) and recorded in parallel with the EMG signals.
The subjects performed three maximal voluntary contractions (MVCs) separated by 2-min rest, after which the electrode grid was located over the abductor pollicis brevis. The subject was asked to identify the activity of a single motor unit (target motor unit) in the surface EMG recordings. To do so, he varied the force level until the activity of a single motor unit could be visually identified from the surface recordings. The feedback consisted of the display of the bipolar surface EMG signals in segments of 500 ms, as previously described (Farina et al. 2004b
). The subject selected the column of the matrix that provided the best feedback. The amplification factor was adjusted to optimize the feedback.
The training phase lasted
20 min, after which ischemia was induced in the hand with a cuff inflated around the forearm at 180 mmHg pressure. The blocking of blood flow was confirmed by the absence of a palpable peripheral pulse. The subject rested for 5 min in ischemic conditions and then activated the target motor unit for 10 s at
12 pulses per second (pps; contraction C1; Fig. 1C) using the surface EMG feedback. Immediately after contraction C1, he performed a 3-min contraction at a force level just below the recruitment threshold of the target motor unit (contraction C2). To avoid recruitment of the target motor unit during the entire duration of contraction C2 (subthreshold contraction), the subject used the surface EMG feedback and continuously checked for the absence of the target motor unit action potentials. No feedback on the exerted force was provided. After contraction C2, the subject repeated a contraction identical to C1, with the target motor unit activated at
12 pps for 10 s (contraction C3). Comparison of the target motor unit conduction velocity between contractions C1 and C3 allowed us to determine if the sustained contraction C2 affected the properties of nonactive muscle fibers.
Eight of the 10 subjects participated in a second experimental session, at least 2 days after the first. The second session was identical to the first, with the only difference that, between contractions C1 and C3, instead of contraction C2, the subject rested for 3 min in ischemic conditions (Fig. 1C). This session served as a control to assure that ischemia by itself did not induce modifications in muscle fiber conduction velocity.
Signal analysis
The target motor unit action potentials in C1 and C3 were detected off-line from the interference EMG signals with a segmentation-classification algorithm previously described (Gazzoni et al. 2004
). The column of the electrode grid that detected action potentials of the target motor unit with the highest amplitude was used for conduction velocity estimation. From the double differential signals computed from the selected column, those distal to the innervation zone showing clear propagation and small shape changes of the potential waveform were selected for conduction velocity estimation. Conduction velocity was computed from each identified action potential of the target motor unit with a multi-channel algorithm previously described (Farina et al. 2001
). Peak-to-peak amplitude, duration (area divided by the peak-to-peak amplitude; Nandedkar et al. 1986
), and mean power spectral frequency of the target motor unit action potentials were estimated from the central channel of those used for conduction velocity estimation. Discharge rate was computed as the inverse of the time interval between subsequent detected discharges of the target motor unit.
For the subject group that performed contraction C2 (C2 group), global conduction velocity, average rectified value, and mean power spectral frequency were estimated during C2 from consecutive and nonoverlapping EMG signal portions of 1 s, using the same channels as for the analysis of the target motor unit action potentials. Global conduction velocity was computed with the same algorithm used for single motor unit conduction velocity estimation (Farina et al. 2001
). Global variables reflected the properties of all active motor units in the C2 contraction.
Data were fit with a regression line, which defined the initial value (intercept of the regression line at time t = 0) and the rate of change over time (slope of the regression line) of the computed variables. The percent decrease was defined as the difference between the value in the end and in the beginning of the contraction, divided by the initial value and expressed in percentage.
Statistical analysis
Data are presented as means ± SD. The target motor unit properties were analyzed using two-way mixed model ANOVA. The repeated measure was the contraction (C1 and C3), and the between-group factor was the subject group (C2 group and control group). Exerted force in the three contractions of the C2 group was analyzed with one-way repeated measures ANOVA, with the contraction as the factor (C1, C2, and C3). Significances revealed by ANOVA were followed by posthoc Student-Newman-Keuls (SNK) pair-wise comparisons. Global surface EMG variables were compared in the beginning and end of C2 with the Wilcoxon matched pairs test. Significance was accepted for P < 0.05.
| RESULTS |
|---|
|
|
|---|
|
The force exerted in C2 (6.5 ± 4.1% MVC) was lower than that in C1 (ANOVA: F = 5.14, P < 0.05; SNK: P < 0.05). Global conduction velocity (initial value, 3.45 ± 0.66 m/s), average rectified value (147 ± 79 µV), and mean power spectral frequency (120 ± 14 Hz) significantly decreased in the 3-min C2 contraction (9.6 ± 5.4, 30.2 ± 30.3, and 6.3 ± 5.6%, respectively; P < 0.01).
Contractions C1 and C3
None of the computed variables significantly changed over time during C1 and C3 in both subject groups. Thus the initial variable values (intercept of the regression lines) were assumed as representative of the contraction and used for further statistical analysis.
FORCE AND DISCHARGE RATE. Exerted force (C2 group: 19.3 ± 17.4 and 12.3 ± 5.7% MVC, for C1 and C3; control group: 16.3 ± 21.2 and 14.9 ± 18.7% MVC, respectively) and discharge rate (C2 group: 14.0 ± 2.8 and 12.3 ± 2.2 pps for C1 and C3, respectively; control group: 15.3 ± 2.0 and 13.7 ± 2.3 pps) were not significantly different between the two subject groups and between C1 and C3.
TARGET MOTOR UNIT ACTION POTENTIAL. Figure 3 shows examples of single motor unit action potentials detected with the two-dimensional array of electrodes. The detection of motor unit action potentials with the array allowed the analysis of the EMG potential distribution over the skin due to single motor unit activity. The column of the grid with the maximum amplitude potentials corresponded to the location of the motor unit over the skin plane (Fig. 3). Increasing distance from the source in the transverse direction decreased the amplitude of the potentials. Potentials detected along fiber direction (i.e., by electrodes in a column) were similar in shape and delayed between each other. The estimation of the delay between action potential waveforms is inversely related to conduction velocity. The innervation zone was identified as the point of inversion of propagation of the action potentials along the columns.
|
|
| DISCUSSION |
|---|
|
|
|---|
Identification of single motor units with surface EMG feedback
We have previously shown that multi-channel surface EMG is an effective feedback for controlling the activity of single motor units at low contraction levels (25% MVC) in normal and ischemic conditions (Farina et al. 2004a
,b
). This study indicates that the same technique can be applied for the identification of motor units at contraction levels of
20% MVC on average, in ischemic conditions. To enhance the feedback in a condition of relatively high background activity of other motor units (Fig. 2), we used a two-dimensional surface EMG recording system instead of the one-dimensional array applied in previous studies. This increased the likelihood that the subject could identify specific EMG channels in which the discrimination of the target motor unit was possible (Fig. 2). All subjects were able to recruit and derecruit the target motor unit over the entire experimental session.
Membrane properties of quiescent muscle fibers
The properties of the target motor unit were analyzed before and after a sustained low force contraction during which it was not active. To enhance the modifications in electrolyte concentrations during low force tasks, the contractions were performed with blood flow occlusion, which avoided metabolite washout. This condition resembles contractions at high force levels, when the intramuscular pressure exceeds the systolic pressure.
The decrease in motor unit conduction velocity between C1 and C3 could not be due to the presence of ischemia by itself because it was not observed in the control group. Moreover, in similar experimental conditions and in the same muscle, conduction velocity of nonactive motor units was not affected by ischemia for
13 min (Farina et al. 2004a
). Conduction velocity depends on the discharge rate (Farina et al. 2004b
; Nishizono et al. 1989
); thus a different discharge rate in C1 and C3 could have influenced the results. Discharge rate in C3 was <2 pps smaller than in C1 for both subject groups, and this difference was not significant. In the same muscle, we have previously observed a sensitivity of conduction velocity to discharge rate of
1%/pps (Farina et al. 2004b
). Thus the eventual decrease in conduction velocity attributable to discharge rate was <2%, much lower than the observed change. In addition, although the decrease in mean discharge rate was similar in the C2 and control groups, no decrease in conduction velocity was detected in the control group.
In the abductor pollicis, we previously observed an
2% decrease in conduction velocity of motor units active for 5 min, with normal blood circulation (Farina et al. 2004b
). In this study, conduction velocity of nonactive motor units decreased three times more in a shorter period of time (3 min) with occluded blood flow. The observed changes in conduction velocity were similar to those reported with motor unit activation after prolonged ischemia (>13 min) (Farina et al. 2004a
). These findings can be explained by the changes that occurred in the extracellular environment during the 3-min subthreshold contraction and indicate that muscle fiber electrophysiological properties are affected by the overall muscle activity, partly independently of the specific fiber activity. This effect was previously advocated as one of the potential mechanisms underlying the modifications in M-wave observed in parts of the gastrocnemius muscle of the cat not elicited by a conditioning fatigue stimulation (Kostyukov et al. 2002
). This is the first time that an effect of active muscle fibers on membrane properties of nonactive ones is shown in vivo at the level of single motor units.
The muscle fiber membrane properties depend on the electrolyte concentration within the muscle. In particular, Na+ and K+ concentration gradients across the fiber membrane influence muscle fiber properties. A flux of Na+ and K+ across the fiber membrane is associated to the action potential generation (Fenn and Cobb 1936
) and is only in part counteracted by the Na+/K+ active pump. As a consequence, K+ concentration in the extracellular environment is modified by the activity of muscle fibers. K+ concentration affects fiber contractile properties and membrane excitability and is an important determinant of muscle fatigue (Jorgensen et al. 1988
; Nielsen et al. 2004
; Nordsborg et al. 2003
; Sejersted and Sjogaard 2000
), in both fast- and slow-twitch muscles (Cairns et al. 1997
).
Increased extracellular K+ concentration decreases action potential conduction velocity (Kossler et al. 1991
). In absence of blood flow, K+ accumulates at a faster rate than with normal circulation (Barcroft and Millen 1939
). Blood flow has indeed a main role in the maintenance of the force level during muscle contraction, as shown by the comparison of sustained and intermittent maximal exercises (Pitcher and Miles 1997
). Accordingly, larger K+ extracellular concentrations were observed after high-force isometric contractions with respect to intermittent contractions (Vyskocil et al. 1983
). In the ischemic conditions studied, the accumulation of K+ thus progressed at a faster rate than with normal blood circulation.
Because the fibers of different motor units are intermingled (Brandstater and Lambert 1973
; Kugelberg and Edstrom 1968
), production of K+ due to the activity of recruited muscle fibers may influence the membrane properties of quiescent ones. The results of this study provide evidence for this hypothesis in single motor units analyzed in vivo and quantify the effect of muscle activity on nonactive muscle fibers. The percent decrease in conduction velocity of the target motor unit between C1 and C3 was not significantly different from the percent decrease in global conduction velocity during C2. Thus in the conditions analyzed, the main determinant of motor unit conduction velocity decrease was not the rate of activation but the modification of the extracellular environment induced by muscle contraction. It is expected that similar changes in conduction velocity could be detected during contractions with normal blood circulation by increasing the contraction level or the duration of the sustained contraction. Ischemia increased the magnitude of the change and allowed us to appreciate it with relative low contraction force and duration, necessary for a successful surface EMG feedback.
The percent decrease in conduction velocity had a large variability among subjects (Table 1). This was probably due to many factors that could not be controlled. The force level was not fixed but corresponded to the subject self-selected force for optimal feedback. The force exerted was indeed very different among different subjects. Accordingly, there was a large variability in the amplitude of the motor unit action potentials among subjects. During C2, the subjects varied the exerted force to avoid the recruitment of the target motor unit. This was probably the reason why EMG average rectified value decreased during C2. Moreover, it is expected that the effect of the activity of muscle fibers on quiescent ones depends on the relative location of the active motor units with respect to the target one. Due to the low contraction level, there was probably a large variability among subjects in the geometrical arrangement of the active motor units, which contributed to the spread of the values reported in Table 1. For higher force levels, this variability would probably be reduced. Although the percent decrease in conduction velocity in the C2 group had large intersubject variability, the muscle activation had a clear effect on the quiescent fibers, which was not observed in the control group.
The change in conduction velocity induced similar relative changes in action potential mean frequency as those observed in previous work (Farina et al. 2004a
). However, although conduction velocity decreased, the duration of the action potential did not significantly change, whereas it was expected to increase (Lindstrom and Magnusson 1977
). This may be due to the variability in the estimation of duration (related to the estimation method) that may have masked the small changes expected (of the order of 6%).
Interpretation of surface EMG during sustained contractions
These findings also have relevance for the interpretation of global surface EMG variables during sustained contraction. Previous studies have shown that average muscle fiber conduction velocity (representative of the mean conduction velocity of all active motor units) initially decreases and then increases during endurance contractions (Gazzoni et al. 2001
; Houtman et al. 2003
). This was interpreted as due to a change in the recruited motor unit pool. Houtman et al. (2003)
hypothesized that the initial conduction velocity decrease was due to the slowing of conduction velocity of the motor units active since the beginning of the contraction, while the subsequent increase revealed progressive recruitment of additional motor units due to fatigue. The results of this study suggest that the increase in average conduction velocity due to additional motor unit recruitment is probably smaller than that expected under the assumption of recruitment of fresh motor units. Thus the quantification of changes in the active motor unit population during sustained contraction, on the basis of the analysis of average conduction velocity time-course, is complicated by the membrane property modifications of the whole motor unit pool.
The change in conduction velocity of quiescent fibers in a submaximal contraction also complicates the interpretation of surface EMG amplitude with fatigue. It has been previously reported that, although EMG amplitude increases during submaximal fatiguing contractions, the amplitude of the surface EMG is significantly less than maximum at the endurance limit (Fuglevand et al. 1993
). It was proposed that the reduced EMG amplitude was due to a deficit in the ability to activate the muscle maximally. However, these results underline at least two additional factors in the interpretation of EMG amplitude at the endurance point. Changed membrane properties of all fibers in the muscle may result in modifications of the twitch torque (increased contraction time) with concomitant decrease in discharge rates of additionally recruited motor units for maintaining tetanic fusion with respect to the condition of fresh fibers (Bigland-Ritchie 1981
). Moreover, decreased conduction velocity increases surface EMG amplitude cancellation (Farina et al. 2004c
; Keenan et al. 2005
).
In summary, this study showed that conduction velocity of quiescent muscle fibers decreases during muscle contraction as a consequence of the alteration of the extracellular environment. Thus the electrophysiological properties of muscle fibers change during muscle contraction, partly independently of their activation.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: D. Farina, Center for Sensory-Motor Interaction, Dept. of Health Science and Technology, Aalborg Univ., Fredrik Bajers Vej 7 D-3, DK-9220 Aalborg, Denmark (E-mail: df{at}hst.aau.dk)
| REFERENCES |
|---|
|
|
|---|
Bigland-Ritchie B. EMG and fatigue of human voluntary and stimulated contractions. Ciba Found Symp 82: 130156, 1981.[Medline]
Brandstater M and Lambert E. Motor unit anatomy. In: New Developments in Electromyography and Clinical Neurophysiology, edited by Desmedt JE. Basel: Karger, 1973, p. 1422.
Cairns SP, Hing WA, Slack JR, Mills RG, and Loiselle DS. Different effects of raised [K+]o on membrane potential and contraction in mouse fast- and slow-twitch muscle. Am J Physiol 273: C598C611, 1997.
Farina D, Gazzoni M, and Camelia F. Conduction velocity of low-threshold motor units during ischemic contractions performed with surface EMG feedback. J Appl Physiol 98: 14871494, 2004a.
Farina D, Gazzoni M, and Camelia F. Low-threshold motor unit membrane properties vary with contraction intensity during sustained activation with surface EMG visual feedback. J Appl Physiol 96: 15051515, 2004b.
Farina D and Merletti R. Methods for estimating muscle fibre conduction velocity from surface electromyographic signals. Med Biol Eng Comput 42: 432445, 2004.[CrossRef][ISI][Medline]
Farina D, Merletti R, and Enoka RM. The extraction of neural strategies from the surface EMG. J Appl Physiol 96: 14861495, 2004c.
Farina D, Muhammad W, Fortunato E, Meste O, Merletti R, and Rix H. Estimation of single motor unit conduction velocity from surface electromyogram signals detected with linear electrode arrays. Med Biol Eng Comput 39: 225236, 2001.[CrossRef][ISI][Medline]
Fenn W and Cobb D. Electrolyte changes in muscle during activity. Am J Physiol 115: 345356, 1936.
Fuglevand AJ, Zackowski KM, Huey KA, and Enoka RM. Impairment of neuromuscular propagation during human fatiguing contractions at submaximal forces. J Physiol 460: 549572, 1993.
Gazzoni M, Farina D, and Merletti R. Motor unit recruitment during constant low force and long duration muscle contractions investigated with surface electromyography. Acta Physiol Pharmacol Bulg 26: 6771, 2001.[Medline]
Gazzoni M, Farina D, and Merletti R. A new method for the extraction and classification of single motor unit action potentials from surface EMG signals. J Neurosci Methods 136: 165177, 2004.[CrossRef][ISI][Medline]
Houtman CJ, Stegeman DF, Van Dijk JP, and Zwarts MJ. Changes in muscle fiber conduction velocity indicate recruitment of distinct motor unit populations. J Appl Physiol 95: 10451054, 2003.
Jorgensen K, Fallentin N, Krogh-Lund C, and Jensen B. Electromyography and fatigue during prolonged, low-level static contractions. Eur J Appl Physiol Occup Physiol 57: 316321, 1988.[Medline]
Juel C. Muscle action potential propagation velocity changes during activity. Muscle Nerve 11: 714719, 1988.[CrossRef][ISI][Medline]
Keenan KG, Farina D, Maluf KS, Merletti R, and Enoka RM. The influence of amplitude cancellation on the simulated surface electromyogram. J Appl Physiol 98: 120131, 2005.
Kossler F, Lange F, Caffier G, and Kuchler G. External potassium and action potential propagation in rat fast and slow twitch muscles. Gen Physiol Biophys 10: 485498, 1991.[ISI][Medline]
Kostyukov AI, Kalezic I, Serenko SG, Ljubisavljevic M, Windhorst U, and Johansson H. Spreading of fatigue-related effects from active to inactive parts in the medial gastrocnemius muscle of the cat. Eur J Appl Physiol 86: 295307, 2002.[CrossRef][Medline]
Kugelberg E and Edstrom L. Differential istochemical effects of muscle contraction on phosphorylase and glycogen in various types of fibers: relation to fatigue. J Neurol Neurosurg Psychiatr 31: 415423, 1968.
Lindstrom L and Magnusson R. Interpretation of myoelectric power spectra: a model and its applications. Proc IEEE 65: 65362, 1977.
McKenna MJ. The roles of ionic processes in muscular fatigue during intense exercise. Sports Med 13: 134145, 1992.[ISI][Medline]
Nandedkar SD, Sanders DB, and Stalberg EV. Automatic analysis of the electromyographic interference pattern. Part I: Development of quantitative features. Muscle Nerve 9: 431439, 1986.[CrossRef][ISI][Medline]
Nielsen JJ, Mohr M, Klarskov C, Kristensen M, Krustrup P, Juel C, and Bangsbo J. Effects of high-intensity intermittent training on potassium kinetics and performance in human skeletal muscle. J Physiol 554: 857870, 2004.
Nishizono H, Kurata H, and Miyashita M. Muscle fiber conduction velocity related to stimulation rate. Electroencephalogr Clin Neurophysiol 72: 52934, 1989.[CrossRef][ISI][Medline]
Nordsborg N, Mohr M, Pedersen LD, Nielsen JJ, Langberg H, and Bangsbo J. Muscle interstitial potassium kinetics during intense exhaustive exercise: effect of previous arm exercise. Am J Physiol Regul Integr Comp Physiol 285: R143R148, 2003.
Pitcher JB and Miles TS. Influence of muscle blood flow on fatigue during intermittent human hand-grip exercise and recovery. Clin Exp Pharmacol Physiol 24: 471476, 1997.[ISI][Medline]
Reucher H, Rau G, and Silny J. Spatial filtering of noninvasive multielectrode EMG: Part Iintroduction to measuring technique and applications. IEEE Trans Biomed Eng 34: 98105, 1987.[ISI][Medline]
Sejersted OM and Sjogaard G. Dynamics and consequences of potassium shifts in skeletal muscle and heart during exercise. Physiol Rev 80: 14111481, 2000.
Sjogaard G and McComas AJ. Role of interstitial potassium. Adv Exp Med Biol 384: 6980, 1995.[Medline]
Vyskocil F, Hnik P, Rehfeldt H, Vejsada R, and Ujec E. The measurement of K+e concentration changes in human muscles during volitional contractions. Pfluegers 399: 235237, 1983.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |