Brain Advance Access originally published online on August 2, 2004
Brain 2004 127(10):2232-2246; doi:10.1093/brain/awh252
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Brain Vol. 127 No. 10 © Guarantors of Brain 2004; all rights reserved
The human spinal cord interprets velocity-dependent afferent input during stepping
1 Department of Neurology and 2 Brain Research Institute, University of California, Los Angeles, CA, USA
Correspondence to: Susan J. Harkema, Department of Neurology, Brain Research Institute, David Geffen School of Medicine at UCLA, Human Locomotion Research Center, Rehabilitation Building, 1000 Veteran Avenue Suite A386, Los Angeles, CA 90095-7147 USA E-mail: sharkema{at}mednet.ucla.edu
Received July 28, 2003. Revised March 15, 2004. Second revision on May 7, 2004. Accepted May 25, 2004.
| Summary |
|---|
|
|
|---|
We studied the motor response to modifying the rate of application of sensory input to the human spinal cord during stepping. We measured the electromyographic (EMG), kinematic and kinetic patterns of the legs during manually assisted or unassisted stepping using body weight support on a treadmill (BWST) in eight individuals with spinal cord injury (SCI). At various treadmill speeds (0.271.52 m/s), we measured the EMG activity of the soleus (SOL), medial gastrocnemius (MG), tibialis anterior (TA), medial hamstrings (MH), vastus lateralis (VL), rectus femoris (RF) and iliopsoas (ILIO); the hip, knee and ankle joint angles; the amount of body weight support (BWS); and lower limb loading. The EMG amplitude and burst duration of the SOL, MG, TA, MH, VL, RF and ILIO were related to the step cycle duration during stepping using BWST. EMG mean amplitudes increased at faster treadmill speeds, and EMG burst durations shortened with decreased step cycle durations. Muscle stretch of an individual muscle could not account for the EMG amplitude modulation in response to stepping speed. The effects on the EMG amplitude and burst duration were similar in subjects with partial and no detectable supraspinal input. We propose that the human spinal cord can interpret complex step-related, velocity-dependent afferent information to contribute to the neural control of stepping.
Key Words: spinal cord injury; locomotion; sensory processing; speed; locomotor training
Abbreviations: ASIA = American Spinal Injury Association; BWL = body weight load; BWST = body weight support on a treadmill; EMG = electromyography; ILIO = iliopsoas; LL = limb load; MG = medial gastrocnemius; MH = medial hamstrings; MTL = muscletendon length; RF = rectus femoris; SCI = spinal cord injury; SL = shank length; SOL = soleus; TA = tibialis anterior; VL = vastus lateralis; VMTL = velocity of muscletendon length
| Introduction |
|---|
|
|
|---|
Peripheral feedback interacting with spinal neural centres allows spinalized animals to adapt to changes in treadmill speed during stepping and execute successful hindlimb stepping (Forssberg and Grillner, 1973
The human nervous system can respond to specific sensory cues related to lower extremity weight-bearing and the kinematics of stepping when supraspinal input is severely compromised (Maegele et al., 2002
). These results suggest that the recovery of walking after spinal cord injury (SCI) may be more dependent on sensory processing than previously considered in humans (Crozier et al., 1992
; Waters et al., 1993
). During stepping, the human spinal cord may be able to generate effective activity in flexors and extensors by integrating complex afferent information related to stepping with remaining supraspinal input after injury. We have previously shown that step-associated proprioception related to alternating leg loading and kinematics can activate and coordinate spinal motor pools in an effective temporal sequence to generate stepping even when an individual is not able to voluntarily execute similar motor patterns in resting positions (Maegele et al., 2002
). Presumably, manually assisted stepping provided an ensemble of sensory inputs derived from the legs that are associated with walking, which facilitated the ability of clinically complete and incomplete SCI subjects to activate leg muscles. For example, the level of limb load (Harkema et al., 1997
) and information from the other leg (Ferris et al., 2004
) were shown to modulate EMG amplitudes during stepping. These modulations were attributed to sensory processing by the human spinal cord.
The aim of this study was to assess whether velocity-dependent afferent input can modulate efferent patterns when supraspinal control is compromised or absent during manually assisted stepping in humans. We varied the stepping velocity by changing treadmill speed, and assessed the modulation of step cycle characteristics as well as leg muscle EMG burst durations and mean amplitudes from individuals with clinically complete and incomplete SCI. We studied these parameters to assess whether velocity-dependent afferent information from the legs during stepping can provide important information to facilitate the generation and modulation of locomotor efferent patterns by the human spinal cord.
| Material and methods |
|---|
|
|
|---|
Subject population
Eight healthy subjects with a thoracic-level SCI resulting from trauma or ischaemia volunteered for this study (Table 1). None of the subjects were taking any antispasticity medication during the study. A clinician assessed the level and extent of SCI according to the American Spinal Injury Association (ASIA) impairment scale (Maynard et al., 1997
|
Experimental design
In all subjects we measured EMG activity, hip, knee and ankle joint angles, level of body weight support (BWS), and amount of lower leg loading during stepping at various treadmill speeds. All subjects wore a harness and were suspended by an overhead, motorized lift. Trainers provided manual assistance as necessary during stepping. Hand placement distal to the patella facilitated extension during the stance phase, and at the popliteal crease for flexion during swing. Hand placement proximal to the ankle assisted proper foot placement and foot clearance at lift-off. A trainer positioned behind the subject aided in hip stabilization and weight shifting between legs during stepping as needed. Data were recorded during stepping of SCI subjects at treadmill speeds between 0.27 and 1.52 m/s while at the same BWS. Speed trials were randomized and the treadmill belt was stopped between each trial.
Data acquisition
We collected EMG, joint angle, footswitch and BWS data at 1 kHz using a 24-channel hard-wired analogue-to-digital board and a custom-written Labview software acquisition program (National Instruments, Austin, TX, USA). EMG data were sampled and AC-coupled into a differential amplifier (Konigsberg Instruments, Pasedena, CA, USA). We measured EMG activity bilaterally from the soleus (SOL), medial gastrocnemius (MG), tibialis anterior (TA), medial hamstrings (MH), vastus lateralis (VL) and rectus femoris (RF) using bipolar surface electrodes. The electrode placements have been described previously (Harkema et al., 1997
; Beres-Jones et al., 2003
). We verified the efficacy of electrode placement in order to avoid contamination by crosstalk among muscles by eliciting and recording EMG activity in an individual muscle while recording inactivity in the agonists and antagonists. The iliopsoas (ILIO) was recorded in five subjects (A1, A7, A8, C4 and C7) using a fine wire electrode placed medial to the anterior superior iliac spine. Eliciting a flexion response and recording EMG through the electrode verified accurate placement. Joint angles were measured using electrogoniometers placed bilaterally at the ankle, knee and hip. BWS was recorded from a loadcell in series with the cable attached to the harness.
Individual leg load was measured during each speed condition by recording pressure distribution of the foot's plantar surface using pressure sensor shoe inserts (Tekscan, Boston, MA, USA). We used Tekscan software to digitally sample the pressure data (50 Hz), which were then converted to vertical ground reaction force data. We dynamically calibrated the Tekscan shoe inserts with a force platform (Kistler, Amherst, NY, USA) using non-disabled subjects of comparable weight and shoe size to the SCI subjects. A previous study describes these procedures in more detail (Harkema et al., 1997
). We used a customized Labview program to interpolate the ground reaction force data acquired from Tekscan software from 50 to 1000 Hz in order to continuously synchronize the load signal with the EMG signals. We also detected foot contact using the footswitch data acquired directly with the Labview acquisition system so we could confirm accurate synchronization of the ground reaction force data. The EMG data were rectified and high-pass filtered at 32 Hz using Labview software customized by our laboratory.
Data analyses
We manipulated the treadmill speed (0.271.52 m/s) to alter the step cycle characteristics, and thus the rate of the kinematics and kinetics during stepping. Figure 1A and B illustrates a typical recording of BWS (top panel) and left and right limb loading calculated from the vertical ground reaction force data (middle and bottom panels) during stepping using BWST and manual assistance. BWS was calculated as the mean force from the loadcell during the series of steps. Step cycle characteristics were determined from the vertical ground reaction force data: step cycle duration (s) as the time from foot contact to next foot contact; stance duration (s) as the time from foot contact to lift-off; and swing duration (s) as the time from lift-off to the subsequent foot contact. The step cycle characteristics were consistently modulated by changing treadmill speeds, as evidenced by exemplary data from these two SCI subjects (Fig. 1C and D).
|
Data from approximately 713 steps from the eight subjects were analysed. Onset of an EMG burst was defined as the time when the signal amplitude remained above the threshold (mean of the baseline+3 standard deviations) for 30 ms. The end of the EMG burst was defined as the time when the signal amplitude remained below the threshold level for 50 ms. Burst duration was calculated as the time between the onset and the end of the EMG burst. Mean EMG amplitude was calculated by dividing the sum of the amplitudes of each burst by the burst duration.
To understand whether the EMG activity modulation was the result of changes in the level of leg loading or immediate afferent feedback from biomechanical events, we assessed the relationships among the EMG activity, limb load (LL), and changes in muscletendon length (MTL) for each muscle. LL was calculated for each leg as the mean ground reaction force during stance. Percentage body weight load (% BWL) was LL divided by the total body weight. MTL was calculated with the acquired joint angle measurements using regression equations that were derived from correlations of joint angle and direct muscletendon complex length measurements from anthropometric specimens (Hawkins and Hull, 1990
). The velocity of muscletendon length change (VMTL) was the derivation of the MTL curve. These values were normalized to the shank length (SL). The relationship among mean EMG activity, MTL and VTML was compared in two ways. First, the mean EMG amplitude was correlated with the muscletendon stretch, defined as the maximum increase in MTL during the entire step cycle. Secondly, to assess the immediate response of the EMG activity to muscletendon stretch, the mean EMG amplitude was correlated to the muscletendon stretch coinciding with the EMG activity, incorporating a 70 ms delay for signal conduction velocity. Onset latencies of EMG activity from biomechanical events (LL, muscletendon stretch) were calculated as the time from the onset of the event to the onset of EMG activity. End latencies of EMG activity from biomechanical events (LL, muscletendon stretch) were calculated as the time from the end of the event to the end of the EMG activity.
Statistical analyses
We used the statistical method for repeated measures for the following relationships for each leg of eight subjects: (i) step cycle duration, stance duration, swing duration and LL versus treadmill speed; (ii) EMG amplitude and burst duration for each muscle versus step cycle duration (load and leg were included in the model below as covariates); and (iii) MTL and VMTL for each muscle versus EMG amplitude. Group analyses (all SCI subjects) and subgroup analyses (clinically complete and incomplete SCI subjects analysed separately) were also performed for each muscle EMG amplitude and burst duration. Stepping tests were carried out 713 times at four to eight different speeds. These variables were repeatedly measured for 610 steps for the seven muscles on both legs for each test and subject. We used a mixed linear model method as a statistical method that incorporates analyses for random effects, repeated measurements, random coefficients, and other statistical problems. For the same muscle, several statistical tests were performed, each with a different treadmill speed. There are several factors that may have influenced the outcome: load, stepping rate, muscle and/or leg. These factors are included in the model as fixed effects. The sequences of the tests were considered as random effects to allow the variation among the tests on the same muscle of a subject. In addition, during each test, the outcome was repeatedly measured for 612 steps; thus, the measurements for the steps in a test might be correlated (e.g. not independent). Therefore, we assume that the error terms are correlated and the covariance matrix of the error term is block diagonal, with each block corresponding to a single test. The structure for the block is assumed to be first-order autoregressive:
![]() |
For all subjects, the model is:
![]() |
is the intercept, Xijk is the vector of fixed effects, such as speed, load, side and step cycle duration, ß is the vector of the regression coefficients for the fixed effects, tij is a random effect of the jth test for subject i; tij
iid N(0,
), and
ijk is the error term. The covariance matrix of the error term is block diagonal, each block corresponding to a single test for subject i. The structure for the block is assumed to be autoregressive order 1. For an individual leg, the same model is used except that there is no i in the subscript. The statistical software SAS was used to carry out all analyses. All tests were two-sided with significance level 0.05.
| Results |
|---|
|
|
|---|
In this study we will present results relating the speed of stepping to step cycle and EMG burst characteristics, including the effect of stepping velocity on step cycle, stance and swing durations and level of loading; the effect of stepping velocity on leg EMG mean amplitudes and burst durations; and the characteristics of the EMG burst pattern in relation to the dynamics of the step cycle.
Stepping speed modulated step cycle, stance and swing durations
Step cycle durations were lower at higher stepping speeds, primarily attributed to a decrease in the duration of the stance phase (Fig. 1AD). Step cycle duration and stance duration decreased in the clinically complete (Fig. 1A and C) and incomplete (Fig. 1B and D) SCI subjects with increasing treadmill speeds. In all SCI subjects, step cycle and stance duration decreased significantly (P < 0.05) in both legs, and swing duration (P < 0.05) decreased significantly in 14/16 legs with higher treadmill speeds (Fig. 2A). However, in all subjects, the reduction in swing duration (0.70 ± 0.40 s) was considerably less than the reduction in stance duration (1.33 ± 0.40 s). Stance duration was longer than swing duration in all SCI subjects except in one clinically complete SCI subject (A14). Step cycle, stance and swing durations were similar between the left and right legs.
|
There was no predictable or consistent trend between load and speed in the SCI subjects. Exemplary data from a clinically complete (Fig. 1E) and incomplete (Fig. 1F) SCI subject show the variation in leg load across the treadmill speeds and also illustrate that while changes may occur from speed to speed, there was no predictable or consistent trend between load and speed. The correlation between leg load and treadmill speed was significant in only 6/16 legs. Furthermore, the leg load changed by less than 13% in 14/16 legs from the slowest to fastest treadmill speeds and actually decreased in 11/16 legs (Fig. 2B). In one subject, there was a significant decrease of leg load (21%) with higher treadmill speeds. We have previously shown that increasing leg loading during stepping can increase EMG mean amplitudes when stepping at a constant speed (Harkema et al., 1997
Stepping velocity modulated leg EMG mean amplitudes and burst durations
Faster stepping speeds resulted in higher motor pool output and shorter burst durations (Table 2). The percentage change in EMG burst duration ranged from 41 to 60% and EMG mean amplitude ranged from 28 to 70% in muscles from the slowest to fastest stepping speed in all SCI subjects. EMG burst durations were significantly (P < 0.05) shorter and EMG mean amplitudes were significantly (P < 0.05) higher in all muscles by group analyses in all SCI subjects. In subgroup analyses (clinically complete and incomplete SCI subjects analysed separately) similar results were found (Table 2). In both SCI groups, EMG burst durations were significantly (P < 0.05) shorter at faster speeds in all muscles. EMG amplitudes were higher in the clinically incomplete than in the complete SCI subjects in all muscles except the ILIO. In clinically complete SCI subjects (n = 5), SOL, MG, TA, MH and ILIO EMG amplitudes were significantly higher at the faster speeds (P < 0.05); however, VL and RF were not. In clinically incomplete SCI subjects (n = 3), SOL, MG, MH, VL and RF were significantly higher at the faster speeds (P < 0.05), but the TA was not. In individual legs analyses, significance in the increase of EMG amplitudes at faster speeds was observed from clinically complete SCI subjects of VL (4/9 legs) and RF (5/9 legs) and from clinically incomplete SCI subjects from the TA (4/6 legs). Thus, non-significant results of the subgroup analyses are probably attributable to individual subject amplitude differences and the lower number of subjects within each group, rather than to subgroup differences.
|
This increase in EMG activity was reflected in a variety of ways. For example, in muscles that displayed little or no EMG burst activity at the slowest speeds, activity often emerged at faster stepping speeds (Fig. 3A, TA and VL muscles; 0.27 and 0.54 m/s). Furthermore, in muscles in which EMG burst activity was already present at a lower speed, amplitudes were even greater at faster velocities (Fig. 3A). This finding of higher EMG amplitude at faster speeds is consistent across steps, as shown by comparing the average EMG mean amplitude from approximately ten steps taken at each speed (Fig. 3B). Burst durations were shorter at faster stepping speeds, as shown by the difference in duration of the EMG burst activity at the three speeds (Fig. 3AC). For example, at 0.27 m/s only one burst of activity occurred in 4 s, whereas two bursts were present at 0.54 m/s and three occurred at 0.89 m/s (Fig. 3A). The increase in EMG amplitude was not attributed to an increase in leg load since there was little change in loading, with even a small decrease at the fastest speed (Fig. 3D). Representative data from six clinically complete and incomplete SCI subjects show the relationships among EMG burst duration (Fig. 4), EMG mean amplitude (Fig. 5) and step cycle duration across the treadmill speeds.
|
|
|
Characteristics of the EMG burst pattern in relation to afferent feedback
The EMG amplitude changes could not be solely attributed to the immediate afferent feedback from stretching and the rate of stretching of the muscletendon complex. As stepping velocity increases, presumably so too does the changes in length of the muscletendon complex, and the velocity at which it lengthens. However, the relationship between EMG amplitude and step cycle duration remained significant (P < 0.05) in all muscles when including the effect of muscletendon length (MTL) and velocity of muscletendon length (VMTL) as covariates with group analyses and subgroup analyses, both when considering muscle stretch during the entire step cycle and only during muscle activation (for details of analyses see Material and methods, and Harkema et al., 1997
The muscle stretch during the entire step cycle could not account for EMG amplitude changes, as it did not increase at the faster stepping speeds. In fact, in all muscles the MTL significantly (P < 0.05) decreased both with group and subgroup analyses. This can be illustrated by examining the changes in EMG amplitude and SOL MTL from a clinically complete SCI subject at three stepping speeds (Fig. 6). The EMG amplitude is higher at the faster stepping speeds (Fig. 6A and D). However, the absolute length change is greater at the slowest speed (0.36 m/s) than at the two faster speeds (0.63 and 0.89 m/s; Fig. 6B). This is also reflected when examining the relationship between EMG amplitude and the extent of SOL muscle stretch at all the stepping speeds; the greater the EMG amplitude the lower the SOL MTL (Fig. 6E). The rate of muscle stretch during the entire step cycle was higher at the faster stepping speed, as expected (Fig. 6C and F); however the correlation between EMG amplitude and VTML was much lower (Fig. 6F; r2 = 0.04) than the correlation between EMG amplitude and step cycle duration (Fig. 6D; r2 = 0.30).
|
The muscle stretch that occurred during the EMG activity could not account for the EMG amplitude modulation with changes in step cycle duration. Several muscles (TA, VL, RF and ILIO) would routinely shorten during EMG activation; thus, it is unlikely that immediate afferent input from muscle stretch solely contributed to the EMG amplitude modulation with speed of stepping. Even in those muscles that primarily lengthened during EMG activation, muscle stretch did not solely account for the EMG modulation with stepping speed. This can be shown by examining these relationships in the SOL muscle, where it has been demonstrated that homonymous muscle stretch can influence the EMG amplitude during stepping (Yang et al., 1991
In addition the SOL muscle burst activity was not phase-locked to either the stretch of the muscletendon complex or the leg loading, indicating that the activity was not dependent only on immediate afferent input related to loading during stance or muscle stretch (Fig. 7). The SOL burst should occur immediately after loading if this stimulus were the primary afferent input generating the activity. However the SOL onset always preceded the onset of the loading in this subject (Fig. 7A, C, D). Furthermore, the end of the SOL burst preceded the end of the loading at all speeds studied (Fig. 7A, C, D). We also measured these latencies of the burst onset and end of loading across all eight subjects in all muscles (15 legs). The timing of the burst pattern in relation to loading did not show that the immediate feedback from loading consistently resulted in the muscle activity. Likewise, if SOL burst activity was dependent only on immediate afferent input related to muscletendon stretch, then the onset of the SOL burst should always have occurred after the initiation of stretch. However, the SOL burst onset preceded the initiation of stretch in the SOL muscle across all speeds (Fig. 7A, B, E). When we measured the latency of the SOL burst onset to the onset of muscletendon stretch in all eight subjects (15 legs) during EMG activity, the timing of the burst pattern in relation to the onset of muscletendon stretch did not support the idea that immediate feedback from stretch consistently resulted in the muscle activity.
|
We also compared the relationships among step cycle duration, EMG amplitude, hip extension position and hip extension velocity (Fig. 8). As exemplified by data from a clinically complete SCI subject, there were significant correlations between hip extension velocity and step cycle duration (r2 = 0.77), EMG amplitude and hip extension velocity (r2 = 0.61), and EMG amplitude and hip extension position (r2 = 0.15) across stepping speeds.
|
| Discussion |
|---|
|
|
|---|
The present data demonstrate that the rate of application of afferent input provides important sensory information that modulates leg motor pool activity in humans with compromised or no detectable supraspinal input during stepping. This efferent modulation cannot be explained by the response to immediate afferent stimuli from a specific stimulus, but rather is the result of an ensemble of stimuli dependent on the rate of application of the sensory input.
Sensory processing by the human spinal cord
We have demonstrated that velocity-dependent modulation of locomotor activity during stepping in human subjects occurs without detectable supraspinal input. In subjects with clinically complete SCI, step cycle and stance phase durations shortened at faster stepping velocities, while swing phase duration showed less modulation (Figs 1 and 2). EMG mean amplitudes were higher and burst durations were shorter at faster stepping velocities in the muscles studied in all SCI subjects (Figs 3
5; Table 2). Since these modulations occurred with no or compromised supraspinal input available, a significant level of neural control of stepping can be attributed to the interpretation of particular patterns of afferent information by the human spinal cord. These data support the idea that the human spinal cord can interpret complex velocity-dependent afferent information during stepping to facilitate the generation and modulation of locomotor efferent patterns.
Spinal neural centres responding to peripheral feedback are essential for spinalized animals to adapt to changes in treadmill speed during stepping (Forssberg and Grillner, 1973
; Forssberg et al., 1980a
; Andersson et al., 1981
; Pierotti et al., 1989
; de Guzman et al., 1991
; Roy et al., 1991
; de Leon et al., 1998
). These locomotor patterns observed at various stepping speeds are similar to those observed in intact cats (Goslow et al., 1973
; Halbertsma et al., 1976
). Furthermore, studies of fictive locomotion report that cycle period modulates burst duration to a greater degree in extensors than flexors (Andersson and Grillner, 1983
; Pearson and Rossignol, 1991
) and more consistently after repetitive step training (Baker et al., 1984
). In our study, the extensor muscles SOL, MG and VL showed the most consistent burst duration response to stepping velocity, as burst durations were shorter at faster speeds in nearly every case (Fig. 4). The changes in step cycle characteristics in SCI subjects observed in this study are also similar to reports from intact animals (Gardiner et al., 1982
; Pierotti et al., 1989
; Roy et al., 1991
), non-disabled humans during walking (Murray et al., 1966
; Grillner et al., 1979
; Shiavi et al., 1981
, 1987
; Nilsson et al., 1985
; Nilsson and Thorstensson, 1987
; Dietz et al., 1994
; Andersson et al., 1997
) and ambulatory clinically incomplete SCI subjects (Pepin et al., 2003a
, b
). However, in these SCI subjects stepping independently on a treadmill, hip, knee and ankle angles and stride length were significantly different from non-disabled subjects during similar slower walking speeds (Visintin and Barbeau, 1994
; Pepin et al., 2003a
, b
).
Afferent modulation of locomotor patterns
Load-related sensory pathways have been shown to modulate the locomotor networks in spinalized animals. The locomotor rhythm can be entrained by stimulation of group I afferents (Conway et al., 1987
; Pearson et al., 1992
) and limb loading during stance activates positive feedback pathways that contribute to ipsilateral extensor muscle activation in animals (Duysens and Pearson, 1980
; Gossard et al., 1994
; Guertin et al., 1995
; McCrea et al., 1995
; Whelan et al., 1995
; Whelan and Pearson, 1997
). Leg loading during the stance phase has also been shown to modulate extensor efferent output during walking in humans (Bastiaanse et al., 2000
; Sinkjaer et al., 2000
; Stephens and Yang, 1999
) and to be mediated at the level of the human spinal cord (Harkema et al., 1997
). Furthermore, it has been reported in humans (Sinkjaer et al., 2000
; Stephens and Yang, 1999
) and animals (Hiebert and Pearson, 1999
) that extensor muscle group I afferent feedback comprises more than 50% of the extensor output during walking. In this study there was either no significant change or a decrease in the amount of weight-bearing at faster speeds; thus, the EMG amplitude modulation with step cycle duration could not be attributed to changes in the level of loading (Figs 1 and 2). In some subjects EMG activity emerged only at the highest stepping speeds even with a lower level of load (Fig. 3). Although the amount of load did not account for the EMG amplitude modulations with stepping speed, the rate of the Ib afferent input could have contributed to the modulation. In addition, the length of time during which load-related information was presented to the spinal networks may have been an important factor in determining the EMG burst durations (Andersson et al., 1981
; Kriellaars et al., 1994
; Whelan and Pearson, 1997
).
Stretch-related afferent information has been demonstrated to contribute 3060% of SOL EMG amplitude during human walking (Yang et al., 1991
). In our study, when MTL and VMTL were considered as covariates, the relationship between SOL EMG amplitude and step cycle remained significant. Furthermore, the onset and end of the SOL burst could precede the onset and end of either the stance phase or the stretch of the SOL muscletendon complex (Fig. 7). In all muscles, the increase in EMG mean amplitude was not highly correlated to the amplitude of muscletendon stretch or rate of stretch in SCI subjects (Fig. 6). Although stretch-related afferent information is likely to be important (Duysens et al., 1991
), especially in the triceps surae, in generating EMG activity, this afferent feedback alone cannot explain the EMG amplitude modulations observed at different stepping velocities in these studies. These results indicate that the additional sensory input derived from stepping at faster velocities also modulates EMG amplitude and cannot be attributed only to immediate input from just one afferent to the same muscle.
Evidence that hip-associated velocity-dependent sensory information can influence central spinal locomotor networks includes the entraining of the bursting cycle by the frequency of hip movements in the fictive cat preparation, which could not be explained by simple immediate afferent reflex (Andersson and Grillner, 1983
). In this preparation, fictive locomotion was pharmacologically induced in a paralysed and spinalized cat that had been completely deafferented except for the hip joint and the surrounding small muscles. They found that varying the frequency of passive sinusoidal hip movement could entrain the fictive locomotor rhythm in both the ipsilateral and contralateral limb in hip, knee and ankle muscles. EMG burst durations shortened and EMG mean amplitudes increased, with extensor burst durations more affected than flexor burst durations. The increase in amplitude was attributed to velocity-dependent response of hip joint and muscle receptors since total angular excursion was the same across frequencies and there was no loading or stance phase. In this study, there was a significant correlation between hip extension velocity and EMG amplitude; however, based on these experiments alone we cannot tell whether the rate of hip extension is the primary afferent stimulus that resulted in the EMG modulation. For example, flexor afferents from other muscles may have contributed to shortening of the extensor phase and onset of flexion phase, as observed in animal studies (Hiebert et al., 1996
), or the rate of load and stretch-related information may also contribute (Yang et al., 1991
).
Other studies have also shown that sensory pathways related to the hip joint and flexor muscles can effectively modulate locomotor patterns (Kriellaars et al., 1994
; Hiebert et al., 1996
; Lam and Pearson, 2001
). Furthermore, the combination of hip extension (Grillner and Rossignol, 1978
; Andersson and Grillner, 1981
) and unloading of one limb at terminal stance was shown to be closely linked to the initiation of swing in the spinal cat (Duysens and Pearson, 1980
; Conway et al., 1987
). Human studies have also identified hip extension as an important sensory cue in infant stepping (Pang and Yang, 2000
) and after human spinal cord injury (Dietz et al., 2002
). In this study, hip extension position was not significantly different across stepping speeds. However, it is interesting to note that at the same hip position, at the slower rates flexor EMG activity was not elicited as it was at the higher rates (Fig. 8). Future studies should be designed to individually manipulate the rate of joint kinematics independent of stepping rate and levels of loading to distinguish the relative contributions of velocity components of hip and other flexor afferents in generating human locomotion.
Response to the rate of application of sensory input in SCI subjects
We propose that the functionally isolated human spinal cord does not just respond to one simple afferent signal, but that it can process and interpret an ensemble of velocity- and load-dependent afferent information during stepping to effectively modulate locomotor EMG patterns. Immediate group I afferent responses (e.g. from rate of stretch or from increased limb load amplitude at faster speeds) may not be solely responsible for the efferent modulations observed in this study, but that is not to say that these afferents do not play a role in modulation of EMG activity during locomotion. However, the rate at which afferent input is delivered to the CNS may be important in more complex regulation of efferent patterns. These results support the idea that the human spinal cord plays a key role in the complex neural control of locomotion, as described in other mammals, incorporating a myriad of sensory input into coherent locomotor output.
Clinical implications
During conventional gait rehabilitation (Ford and Duckworth, 1974
; Somers, 1992
), locomotor training for recovery of walking (Stewart et al., 1991
; Wernig and Müller, 1991
; Barbeau and Fung, 1992
; Dietz et al., 1995
; Wernig et al., 1995
; Colombo et al., 1998
; Field-Fote, 2001
) and studies of locomotion after SCI (Fung and Barbeau, 1994
; Visintin and Barbeau, 1994
; Harkema et al., 1997
; Maegele et al., 2002
), slower than normal walking speeds are primarily used. Furthermore, normative speeds of stepping have been unattainable in ambulatory SCI individuals during independent stepping on a treadmill (Pepin et al., 2003a
, b
).
The results from this study demonstrate that stepping at more normative treadmill speeds results in higher leg EMG output during stepping when compromised or when no supraspinal input is available to the spinal cord. The level of leg loading has also been shown to modulate EMG amplitude in SCI subjects stepping at a constant treadmill speed. From a clinical perspective it will be important in future studies to assess the interactive effect of the level of loading and stepping speed in generating the most effective locomotor patterns in individuals following SCI. Furthermore, to maximize recovery of walking after neurological injury, rehabilitative strategies should include approaches that repetitively provide the appropriate sensory cues associated with locomotion by including normative stepping speeds while maximizing loading on the legs.
| Acknowledgements |
|---|
This research was supported by NIH grants NS16333, NS38654, HD07416 and MO1-RR00865. The authors would like to thank the research volunteers for their valuable contribution to this study and UCLA Human Locomotion Research Center members for their assistance in data collection and management. The authors would like to acknowledge the collaboration with Jennifer Foster and Uday Patel for their earlier contributions to this study and HeJing Wang for statistical support. The authors also acknowledge the personnel of Rehabilitation Clinics, Langensteinbach, Karlsbad, Germany, including Dr Anton Wernig and Sabina Müller, for locomotor training (Laufband therapy) of two participants in this study.
| References |
|---|
|
|
|---|
Andersson O, Grillner S. Peripheral control of the cat's step cycle. I. Phase dependent effects of ramp-movements of the hip during fictive locomotion. Acta Physiol Scand 1981; 113: 89101.[ISI][Medline]
Andersson O, Grillner S. Peripheral control of the cat's step cycle. II. Entrainment of the central pattern generators for locomotion by sinusoidal hip movements during fictive locomotion. Acta Physiol Scand 1983; 118: 22939.[ISI][Medline]
Andersson O, Forssberg H, Grillner S, Wallen P. Peripheral feedback mechanisms acting on the central pattern generators for locomotion in fish and cat. Can J Physiol Pharmacol 1981; 59: 71326.[ISI][Medline]
Andersson EA, Nilsson J, Thorstensson A. Intramuscular EMG from the hip flexor muscles during human locomotion. Acta Physiol Scand 1997; 161: 36170.[CrossRef][ISI][Medline]
Baker LL, Chandler SH, Goldberg LJ. L-Dopa-induced locomotor-like activity in ankle flexor and extensor nerves of chronic and acute spinal cats. Exp Neurol 1984; 86: 51526.[CrossRef][ISI][Medline]
Barbeau H, Fung J. New experimental approaches in the treatment of spastic gait disorders. Med Sci Sports 1992; 36: 23446.
Bastiaanse CM, Duysens J, Dietz V. Modulation of cutaneous reflexes by load receptor input during human walking. Exp Brain Res 2000; 135: 18998.[CrossRef][ISI][Medline]
Beres-Jones JA, Johnson TD, Harkema SJ. Clonus after human spinal cord injury cannot be attributed solely to recurrent muscletendon stretch. Exp Brain Res 2003; 149: 22236.[ISI][Medline]
Colombo G, Wirz M, Dietz V. Effect of locomotor training related to clinical and electrophysiological examinations in spinal cord injured humans. Ann NY Acad Sci 1998; 860: 5368.
Conway BA, Hultborn H, Kiehn O. Proprioceptive input resets central locomotor rhythm in the spinal cat. Exp Brain Res 1987; 68: 64356.[ISI][Medline]
Crozier KS, Cheng LL, Graziani V, Zorn G, Herbison GJ, Ditunno JF. Spinal cord injury: prognosis for ambulation based on quadriceps recovery. Paraplegia 1992; 30: 7627.[ISI][Medline]
de Guzman CP, Roy RR, Hodgson JA, Edgerton VR. Coordination of motor pools controlling the ankle musculature in adult spinal cats during treadmill walking. Brain Res 1991; 555: 20214.[CrossRef][ISI][Medline]
de Leon RD, Hodgson JA, Roy RR, Edgerton VR. Locomotor capacity attributable to step training versus spontaneous recovery after spinalization in adult cats. J Neurophysiol 1998; 79: 132940.
Dietz V, Zijlstra W, Duysens J. Human neuronal interlimb coordination during split-belt locomotion. Exp Brain Res 1994; 101: 51320.[ISI][Medline]
Dietz V, Colombo G, Jensen L, Baumgartner L. Locomotor capacity of spinal cord in paraplegic patients. Ann Neurol 1995; 37: 57482.[CrossRef][ISI][Medline]
Dietz V, Muller R, Colombo G. Locomotor activity in spinal man: significance of afferent input from joint and load receptors. Brain 2002; 125: 262634.
Duysens J, Pearson KG. Inhibition of flexor burst generation by loading ankle extensor muscle in walking cats. Brain Res 1980; 187: 32132.[CrossRef][ISI][Medline]
Duysens J, Tax AA, van der Doelen B, Trippel M, Dietz V. Selective activation of human soleus or gastrocnemius in reflex responses during walking and running. Exp Brain Res 1991; 87: 193204.[ISI][Medline]
Ferris DP, Gordon KE, Beres-Jones JA, Harkema SJ. Muscle activation during unilateral stepping occurs in the nonstepping limb of humans with clinically complete spinal cord injury. Spinal Cord 2004; 42: 1423.[CrossRef][ISI][Medline]
Field-Fote EC. Combined use of body weight support, functional electric stimulation, and treadmill training to improve walking ability in individuals with chronic incomplete spinal cord injury. Arch Phys Med Rehabil 2001; 82: 81824.[CrossRef][ISI][Medline]
Ford JR, Duckworth B. Physical management for the quadriplegic patient. Philadelphia: F.A. Davis, 1974.
Forssberg H, Grillner S. The locomotion of the acute spinal cat injected with clonidine i.v. Brain Res 1973; 50: 1846.[CrossRef][ISI][Medline]
Forssberg H, Grillner S, Halbertsma J. The locomotion of the low spinal cat. I. Coordination within a hindlimb. Acta Physiol Scand 1980a; 108: 26981.[ISI][Medline]
Forssberg H, Grillner S, Halbertsma J, Rossignol S. The locomotion of the low spinal cat. II. Interlimb coordination. Acta Physiol Scand 1980b; 108: 28395.[ISI][Medline]
Fung J, Barbeau H. Effects of conditioning cutaneomuscular stimulation on the soleus H-reflex in normal and spastic paretic subjects during walking and standing. J Neurophysiol 1994; 72: 2090104.
Gardiner KR, Gardiner PF, Edgerton VR. Guinea pig soleus and gastrocnemius electromyograms at varying speeds, grades, and loads. J Appl Physiol 1982; 52: 4517.
Goslow GE, Stauffer EK, Nemeth WC, Stuart DG. The cat step cycle: responses of muscle spindles and tendon organs to passive stretch within the locomotor range. Brain Res 1973; 60: 3554.[CrossRef][ISI][Medline]
Gossard J-P, Brownstone R, Barajon I, Hultborn H. Transmission in a locomotor-related group Ib pathway from hindlimb extensor muscles in the cat. Exp Brain Res 1994; 98: 21328.[ISI][Medline]
Grillner S, Rossignol S. On the initiation of the swing phase of locomotion in chronic spinal cats. Brain Res 1978; 146: 26977.[CrossRef][ISI][Medline]
Grillner S, Halbertsma J, Nilsson J, Thorstensson A. The adaptation to speed in human locomotion. Brain Res 1979; 165: 17782.[CrossRef][ISI][Medline]
Guertin P, Angel MJ, Perreault MC, McCrea DA. Ankle extensor group I afferents excite extensors throughout the hindlimb during fictive locomotion in the cat. J Physiol 1995; 487.1: 197209.
Halbertsma J, Miller S, van der Meche FGA. Basic phasing for the flexion and extension movements of the limbs during locomotion. In: Herman RM, editor. Neural control of locomotion. New York: Plenum Press, 1976. p. 489518.
Harkema SJ, Hurley SL, Patel UK, Requejo PS, Dobkin BH, Edgerton VR. Human lumbosacral spinal cord interprets loading during stepping. J Neurophysiol 1997; 77: 797811.
Hawkins D, Hull ML. A method for determining lower extremity muscletendon lengths during flexion/extension movements. J Biomech 1990; 23: 48794.[CrossRef][ISI][Medline]
Hiebert GW, Pearson KG. Contribution of sensory feedback to the generation of extensor activity during walking in the decerebrate cat. J Neurophysiol 1999; 81: 75870.
Hiebert GW, Whelan PJ, Prochazka A, Pearson KG. Contribution of hind limb flexor muscle afferents to the timing of phase transitions in the cat step cycle. J Neurophysiol 1996; 75: 112637.
Kriellaars DJ, Brownstone R, Noga BR, Jordan LM. Mechanical entrainment of fictive locomotion in the decerebrate cat. J Neurophysiol 1994; 71: 207486.
Lam T, Pearson KG. Proprioceptive modulation of hip flexor activity during the swing phase of locomotion in decerebrate cats. J Neurophysiol 2001; 86: 132132.
Maegele M, Muller S, Wernig A, Edgerton VR, Harkema SJ. Recruitment of spinal motor pools during voluntary movements versus stepping after human spinal cord injury. J Neurotrauma 2002; 19: 121729.[CrossRef][ISI][Medline]
Maynard FM, Bracken MB, Creasey G, Ditunno JF, Donovan WH, Ducker TB, et al. International Standards for Neurological and Functional Classification of Spinal Cord Injury. Spinal Cord 1997; 35: 26674.[CrossRef][ISI][Medline]
McCrea DA, Shefchyk SJ, Stephens MJ, Pearson KG. Disynaptic group I excitation of synergist ankle extensor motoneurones during fictive locomotion in the cat. J Physiol 1995; 487: 52739.[ISI][Medline]
Murray MP, Kory RC, Clarkson BH, Sepic SB. Comparison of free and fast speed walking patterns of normal men. Am J Phys Med 1966; 45: 823.[Medline]
Nilsson J, Thorstensson A. Adaptability in frequency and amplitude of leg movements during human locomotion at different speeds. Acta Physiol Scand 1987; 129: 10714.[ISI][Medline]
Nilsson J, Thorstensson A, Halbertsma J. Changes in leg movements and muscle activity with speed of locomotion and mode of progression in humans. Acta Physiol Scand 1985; 123: 45775.[ISI][Medline]
Pang MY, Yang JF. The initiation of the swing phase in human infant stepping: importance of hip position and leg loading. J Physiol 2000; 528: 389404.
Pearson KG, Rossignol S. Fictive motor patterns in chronic spinal cats. J Neurophysiol 1991; 66: 187487.
Pearson KG, Ramirez JM, Jiang W. Entrainment of the locomotor rhythm by group Ib afferents from ankle extensor muscles in spinal cats. Exp Brain Res 1992; 90: 55766.[ISI][Medline]
Pepin A, Ladouceur M, Barbeau H. Treadmill walking in incomplete spinal-cord-injured subjects: 2. Factors limiting the maximal speed. Spinal Cord 2003a; 41: 2719.[CrossRef][ISI][Medline]
Pepin A, Norman KE, Barbeau H. Treadmill walking in incomplete spinal-cord-injured subjects: 1. Adaptation to changes in speed. Spinal Cord 2003b; 41: 25770.[CrossRef][ISI][Medline]
Pierotti DJ, Roy RR, Gregor R, Edgerton VR. Electromyographic activity of cat hindlimb flexors and extensors during locomotion at varying speeds and inclines. Brain Res 1989; 481: 5766.[CrossRef]









