 |
||||
|
||||
|
||||
|
||||
|
||||
Previous Article | Next Article
The Journal of Neuroscience, January 1, 2000, 20(1):409-426
Rapid Correction of Aimed Movements by Summation of Force-Field Primitives
William J. Kargo and Simon F. Giszter Neurobiology Department, Medical College of Pennsylvania/Hahnemann University, Philadelphia, Pennsylvania 19129
ABSTRACT
TOP
ABSTRACT
INTRODUCTION
Spinal circuits form building blocks for movement construction. In the frog, such building blocks have been described as isometric force fields. Microstimulation studies showed that individual force fields can be combined by vector summation. Summation and scaling of a few force-field types can, in theory, produce a large range of dynamic force-field structures associated with limb behaviors. We tested for the first time whether force-field summation underlies the construction of real limb behavior in the frog. We examined the organization of correction responses that circumvent path obstacles during hindlimb wiping trajectories. Correction responses were triggered on-line during wiping by cutaneous feedback signaling obstacle collision. The correction response activated a force field that summed with an ongoing sequence of force fields activated during wiping. Both impact force and time of impact within the wiping motor pattern scaled the evoked correction response amplitude. However, the duration of the correction response was constant and similar to the duration of other muscles activated in different phases of wiping. Thus, our results confirm that both force-field summation and scaling occur during real limb behavior, that force fields represent fixed-timing motor elements, and that these motor elements are combined in chains and in combination contingent on the interaction of feedback and central motor programs.
Key words: force field; spinal cord; reflex; superposition; primitive; motor control
INTRODUCTION
The spinal cord of lower vertebrates and mammals may be organized into movement control modules (Grillner and Wallen, 1985
; Loeb, 1985
Spinal microstimulation provided the initial evidence that the frog spinal cord may be organized into modules that produce force-field primitives (Bizzi et al., 1991
Descending and segmental systems might recruit, scale, and sum force-field primitives in a manner similar to microstimulation (Bizzi et al., 1991
We therefore examined in this paper how force-field mechanisms might be used to generate and control limb behavior. Specifically, we examined how the targeted wiping reflex of the spinal frog is organized and how on-line trajectory corrections are achieved. We found that force-field primitives were recruited in chains and in combination to produce wiping trajectories and their on-line adjustments. Furthermore, each primitive exhibited a common activation waveform of constant duration. Taken together, our data suggest that constant-duration force-field primitives form building blocks for the construction of certain limb behaviors.
MATERIALS AND METHODS
Surgery
Detailed surgical procedures have been outlined in previously published material (Giszter et al., 1993Experimental protocol and design
Fourteen frogs were used in these experiments. We tested the hindlimb-to-hindimb form of wiping. In this form, the wiping limb and target limb both move such that the stimulus and the ankle of the wiping limb (i.e., the wiping tool in this form) meet near the midline of the frog, just caudal to the cloaca [for detailed kinematic descriptions, see Giszter et al. (1989)Data collection
Limb kinematics, EMGs, and contact forces with the rod were recorded during free limb trials. Movements of the wiping and target limbs were videotaped at 60 fields per second perpendicular to the horizontal plane of the wiping movement. A light-emitting diode was triggered at the initiation of the wiping stimulus and allowed movements and EMGs to be synchronized within 1-16.67 msec. EMGs were amplified 1000× and bandpass-filtered (cutoffs 100 Hz and 10 kHz) using A-M Systems amplifiers (A-M Systems, Everett, WA). EMGs were analog-to-digital converted using a DAS16 A/D converter and stored on computer using custom-made software. Forces were sampled at 250 Hz, and EMGs were sampled at 1000 Hz. During isometric trials, forces generated at the ankle were recorded by using the same force/torque transducer to immobilize the ankle. Figure 1 shows the experimental conditions (with and without local cutaneous deafferentation) and the data types collected (free limb kinematics and isometric forces; EMGs were collected in both restrained and unrestrained trials).
View larger version (29K):[in this window]
[in a new window]
Figure 1. The hindlimb to hindlimb wiping trials described here were evoked in spinalized frogs in two conditions: with sensory feedback intact and with a local cutaneous deafferentation of the calf/foot of the wiping limb (shaded area). Wiping was evoked in all trials by electrical stimulation of the target limb (see Materials and Methods). Two data types were collected: isometric forces and free limb kinematics. Joint angles were measured as shown and converted to engineering format as needed (e.g., for torque calculations). In trials in which the wiping limb was free to move, kinematic data were recorded after the wiping was evoked. In trials in which the wiping limb was restrained, isometric forces were recorded at the ankle. Hindlimb muscle EMGs were recorded during both trial types (free limb and isometric) and during both conditions (feedback intact vs cutaneous deafferentation).
Data analysis
Kinematic analysis. We examined changes in the free limb kinematics that occurred during obstacle collision trials. The hip and knee joint locations of the wiping and target limbs (in the horizontal plane of the limb) were digitized from video recordings using a custom-designed program. The joint coordinates from the digitization procedure were imported into the S-Plus Statistics and Display package (Statistical Science, Seattle, WA). These were used to reconstruct movement of the thigh and calf and to derive the hip and knee joint angles over time. Hip angle was measured as the angle between the longitudinal axis of the frog and the thigh (Fig. 1). Knee angle was measured as the internal angle between the thigh and calf (Fig. 1). Because there was little variation in trial-to-trial duration of the wiping movement for an individual frog, time normalization was not necessary. Thus, averages and SDs of the hip and knee joint angles over time were constructed after simply aligning individual trials at the onset of movement for each condition. Force analysis (single limb position). We examined differences in the endpoint isometric forces generated during wiping trials and between conditions, i.e., when sensory feedback was intact and when skin regions contacting the ankle restraint were deafferented. We first examined differences in the direction and magnitude of force vectors generated at a single limb configuration. The direction and magnitude of the wiping force vectors were calculated for a 3 sec period after the reflex was evoked (wiping motor pattern typically lasted 600-1000 msec). Endpoint polar force direction was calculated in the plane of the limb (positive x-axis at 0°, longitudinal axis of frog at 90°) from the measured endpoint force at time (t) by:
where Fy and Fx are the recorded horizontal force components over time. The magnitude of the force vectors was calculated by:
(1)
Averages and SDs of the force direction and magnitude at each time point after initiation of wiping were constructed by aligning individual trials at the onset of force. The onset of force for each trial was determined after the averaged baseline force (200 msec before wiping was initiated) was subtracted from the force record. Both force direction and magnitude were determined for active forces, i.e., when resting, baseline forces were subtracted. On the basis of analysis of EMG patterns (see EMG analysis), we believed that feedback caused by skin contact with the ankle restraint initiated a corrective response during isometric wiping trials. We tested whether correction responses recruited an isometric force-field primitive and whether this primitive was simply superimposed on the underlying wiping forces. According to the definition used in this paper, a force-field primitive is observed as a structurally invariant force field over time. In its most general form, a force-field primitive would be a function of both position and velocity. In this paper, we confine ourselves to isometric force measurements at which velocity is always zero. Other experiments [for preliminary description, see Giszter et al. (2000)
(2) (r) (Eq. 3). That is to say, a constant structure force field is scaled in amplitude through time in a similar way at all locations. To derive the endpoint forces generated by the evoked correction response under the hypothesis that force-field primitives comprised the wiping and corrective patterns, we therefore proceeded as follows. We subtracted the averaged force vectors generated during wiping trials alone, i.e., after cutaneous deafferentation, from the averaged force vectors generated during wiping trials with an evoked correction response, i.e., with cutaneous feedback intact:
(3)
(4)
where FC(r,t) represents the total forces generated in a corrected wipe over time (t) at configuration (r) with cutaneous feedback intact. This field is minimally composed of a sum of primitives
(5) (r) and
where lt and lc are lengths of the thigh and calf and
(6) h and
where "
(7) " represents the inner product of two vectors. Ft and Ft+m denote the two resultant force vectors that are compared at times t and t + m at N locations x1, x2, ... xN. The cosine of the angle between the two sampled fields is calculated:
where:
(8)
represents the norm of a sampled force field. If the correction is generated by a field, which can be expressed in the form of Equation 5 above, then at two chosen time points (t) and (t + m):
(9)
(10)
because the scalars b(t) and b(t + m) determining force magnitude do not affect the result of the cosine operator in Equation 8. The value of 1 indicates that the two fields tested are simply scaled versions of one another. If the value of the cosine measure from Equation 8 was 0.9 for two test fields, these were considered to have a similar force-field structure, as described in Mussa-Ivaldi et al. (1994)
RESULTS
Kinematic observations of the trajectory correction response
We first analyzed the kinematic features of corrective responses during the hindlimb-to-hindlimb form of wiping. These corrections have never been described in any detail. This kinematic analysis will form the basis of subsequent analyses. In this form of wiping, when an irritant is applied to the foot of the target limb, the wiping limb attempts to reach and remove this stimulus. Figure 2A shows the normal, free limb kinematics of the wiping and target limbs (i.e., the thigh and calf of each limb) after we evoked this reflex behavior using our standard electrical stimulation (500 msec duration, 2 msec biphasic pulses, 33 Hz, 5 V). The ankle of the wiping limb followed a fairly direct trajectory to a region caudal to the cloaca where target limb contact occurred. In the trial shown, target limb contact occurred 233 msec after the onset of wiping limb movement. After contact with the stimulus on the target limb, the ankle of the wiping limb followed a trajectory that was directed away from the body of the frog. This final portion of the trajectory ["whisk"; see Berkinblit et al. (1986)
View larger version (32K):[in this window]
[in a new window]
Figure 2. Path obstacles evoked trajectory correction responses during hindlimb wiping movements. A-D, Digitized stick figures of the wiping and target limbs (thigh and calf) and the ankle path of both limbs (traced by a bold line; arrows mark direction of movement along path) are shown after initiation of the wiping reflex in four experimental manipulations. Each stick figure is separated in time from the next by 16.67 msec ( sec). Bold stick figures mark the initial limb configurations and their configuration at target limb contact. The frog's body would be above each stick figure. A, The unobstructed kinematics of the wiping limb with feedback intact consisted of a knee flexion and hip extension that resulted in the ankle moving along a trajectory directed to a midline region of stimulus contact. A final knee extension ("whisk") occurred on stimulus contact. B, When the wiping limb collided with a path obstacle en route to stimulus contact, a hip flexion was evoked that enabled the limb to clear the obstacle. Once clear there was an enhanced knee flexion. Target limb contact occurred at a similar midline site and at a time similar to unobstructed trials (see Results). C, After local removal of cutaneous feedback from the wiping limb, unobstructed free limb kinematics showed no apparent differences from the kinematics seen with feedback intact (compare A, C). D, When obstacles were placed in the path of the locally deafferented limb, corrective responses were no longer evoked, and the limb was driven further into the obstacle (compare B, D).
Spinal frogs exhibited trajectory correction responses when obstacles were placed in the movement path of the wiping limb. In Figure 2B such a trial is shown. The trajectory of the ankle, after contact with the obstacle, was redirected toward the head of the frog. This rostral redirection of the ankle trajectory was caused by a hip flexion, which was not seen in the unimpeded movement. Once clear of the obstacle, there was an enhanced knee flexion, and the ankle trajectory exhibited an increased velocity toward the target region (compare distance between digitized ankle positions). Interestingly, the time of target limb contact as determined from the video record was similar to the unimpeded trial. In the frog shown in Figure 2, the averaged difference in the time of target limb contact between unimpeded and impeded wiping trials (in successful obstacle avoidance) was only 30 msec.
In this portion of the study we examined nine frogs that all exhibited successful correction responses when path obstacles were placed close to the ankle of the wiping limb, i.e., with obstacles placed ~2-3 cm or less from the initial position of the ankle. The straight-line distance from the initial position of the ankle to the site of target limb contact ranged from 10 to 13 cm for the nine frogs. In each frog, two common kinematic features characterized successful correction responses (Figs. 2, 3A,B, bottom panels): (1) an evoked hip flexion that allowed the ankle to clear the obstacle and (2) an enhanced knee flexion once the obstacle was cleared. The enhanced knee flexion allowed the limb to extend at the hip without further collision and also resulted in an increased ankle velocity toward the target region. Taken together, these effects acted to maintain target limb contact times that were close to contact times in unimpeded wiping trials. In five of nine frogs there was only a small increase of ~30 msec (significant at p < 0.05) in the mean latency to target limb contact (mean of 283 ± 20 msec SD with path obstacles compared with 250 ± 13 msec SD without obstacles) (Table 1). Thus, although there was a small and significant difference, these frogs were able to maintain contact times within the extremes of the range of variations seen during unobstructed trials. However, in the other four frogs there were larger increases (p < 0.01) in the mean latency to target limb contact after movement initiation (mean 316 ± 25 msec SD with path obstacles compared with 245 ± 12 msec SD without). The increased latency to target limb contact in these four frogs was mainly attributable to an increase in the time to clear the obstacle (mean duration of 125 msec in contact with the obstacle in this group compared with a mean duration of 66 msec in the five frogs with almost normal movement durations).
View larger version (43K):[in this window]
[in a new window]
Figure 3. Wiping motor patterns and evoked correction responses were similar for free limb and isometric trials. Rectified EMGs are shown for six wiping limb muscles after the wiping reflex was evoked in four separate conditions: (A) unobstructed free limb trial, (B) free limb trial with obstacle collision, (C) isometric trial, and (D) isometric trial after local cutaneous deafferentation of the wiping limb. Below the EMG panels in A and B are shown obstacle contact forces and hip and knee joint kinematics (in A, there was no obstacle collision). First and second vertical lines mark movement onset and target limb contact, respectively. Below the EMG panels in C and D are shown the recorded ankle force components (Fx, Fy, and Fz) and the derived hip and knee torques. A, The wiping motor pattern consisted of three phases: initial knee flexor (ST), subsequent hip extensor (RI), and final knee extensor (VE) activation. Hip flexors (IL) exhibited phasic activity near the time of stimulus contact. B, The timing of wiping motor phases was unaltered by obstacle collision. Hip flexor muscles (IL, GL, and SA) showed a strong synchronous activation after obstacle contact (see contact force) that ended near the onset of RI. Their activation was associated with a reversal of hip joint motion to flexion. After obstacle clearance there was an enhanced knee flexion and resumed hip extension. C, The time of onset of knee flexor (ST) and hip extensor phases (RI) was unaltered during isometric trials. The knee extensor phase (VE) was absent. Hip flexor muscles showed a strong activation after the onset of force, which represented the correction response that ended near the onset of RI. The time of hip flexor activation was associated with production of a hip flexor torque (bolded line below), whereas RI activation was associated with the onset of hip extensor torques. D, The time of onset of knee flexor and hip extensor phases was unaltered after local cutaneous deafferentation. Hip flexor activity and hip flexor torque associated with the correction response were eliminated by the deafferentation.
View this table: [in this window]
[in a new window]
Table 1. Target-to-target limb contact The movement path of the contralateral (i.e., target) limb did not exhibit any major differences during wiping trials with and without the path obstacle (Fig. 2), and in this study, we focus only on corrective changes within the wiping limb.
Trajectory correction responses are evoked by cutaneous feedback
The sensory basis of the evoked correction response in frogs has not been described. Muscle, joint, or skin feedback could all contribute to initiating and controlling the response. We determined that the correction response in the wiping limb was initiated primarily by feedback from cutaneous receptors. We could demonstrate this by local removal of cutaneous feedback from the calf and foot of the wiping limb alone (see Materials and Methods). The successful correction responses depended on cutaneous feedback stemming from the obstacle collision. In the frog shown in Figure 2, after cutaneous deafferentation, the wiping limb did not exhibit a successful corrective response and was in fact driven further into the path obstacle by the extending hip.
Cutaneous feedback, although necessary for evoking successful correction responses, was not necessary for controlling the unimpeded free limb trajectory. In Figure 2C, the unimpeded kinematics of the wiping limb are shown after the local cutaneous deafferentation. The kinematics of the wiping limb did not exhibit any major differences, e.g., in the latency to target limb contact (233 msec), in its position at target limb contact (knee angle 50°, hip angle 42°), and in the shape of the endpoint path, from the trials in which cutaneous feedback was intact [see also Kargo and Giszter (2000)
To show that loss of corrections was not attributable to nonspecific effects of the loss of cutaneous feedback, we performed a second series of experiments. Rather than using cutaneous deafferentation, we used bone pins placed in the tibia to immobilize the limb during the wiping reflex just as an obstacle would, but without skin contact. At the same time, skin sensory fields were largely undisturbed. We were unable to initiate corrective EMGs during isometric holds (see below) in this arrangement in the absence of cutaneous contact. However, true collisions evoked corrections. These approaches allowed us to manipulate wiping and corrective responses in isometric experiments.
Basis for isometric examination of wiping and correction responses
In the following sections, we tested whether the observed correction response forces FC(r,t) were generated by the on-line superposition of a structurally invariant, corrective force field with an underlying wiping force field. Specifically:
where FN represents the normal pattern of forces over time (t) at position (r) activated by the skin stimulus on the contralateral limb and FR represents the force-field primitive activated by ipsilateral obstacle contact. To test this hypothesis we recorded the endpoint isometric forces generated during wiping and trajectory correction responses. Previous studies have shown that forces and EMGs appropriate for wiping can be evoked in the isometric limb (Giszter et al., 1993
(11) Free limb EMG patterns generated during wiping and trajectory corrections
We first examined the basic motor pattern that occurs during free limb wiping trials. The normal wiping motor pattern consisted of three periods of muscle activity (for detailed analysis, also see Kargo and Giszter (2000)
A specific group of hindlimb muscles was commonly activated in response to obstacle collision during free limb wiping trials. With obstacles placed 2-3 cm from the initial position of the ankle, collisions occurred during the initial knee flexor phase (see contact force below EMGs in Fig. 3B). In the frog shown in Figure 3A,B, a group of hip flexor-related muscles (IL and GL, monoarticular hip flexors, and SA, a hip/knee flexor) were reflexively activated ~15 msec after obstacle contact. The duration of their activation was ~150 msec. The reflex activation of these muscles did not appear to interfere with or cause variation of the onset times of subsequent wiping motor phases (see section entitled Correction responses are generated by addition of a muscle synergy to the basic pattern and Table 2). None of the comparisons of motor pattern latency differed significantly between correction and free limb, unobstructed trials.
View this table: [in this window]
[in a new window]
Table 2. Latency in milliseconds of the first two motor phases of wiping in four different conditions Motor patterns are similar under isometric conditions
After evoking a series of free limb wiping trials in each frog, we immobilized the wiping limb using our ankle restraint. This apparatus allowed us to examine isometric forces and force-field patterns. We again recorded EMGs and also isometric forces generated at the ankle after evoking the wiping reflex. We found that the initial knee flexor phase and the subsequent hip extensor phases of the wiping motor pattern were initiated in all frogs and remained similar to free limb patterns (see Fig. 3C, restrained trial, and compare with Fig. 3B, free limb trial with correction). The hip extensor phase maintained a similar time of onset relative to knee flexor onset in eight of nine frogs (Table 2). With the exception of bf74, the ST to RI latency did not differ significantly between free limb, correction, and isometric tests. However, in the majority of frogs (eight of nine), the usual, final phase of knee extensor activation was absent (Fig. 3C,D). In the one exception, the knee extensor phase was much reduced, although present. Loss of knee extensor activation was often associated with a more abrupt termination of knee flexor and hip extensor EMG activity (compare ST and RI EMGs in Fig. 3C,D, isometric, with Fig. 3A,B, free limb). The lack of the knee extensor burst representing initiation of the final motor phase suggests that either triggering or aborting of this phase of the motor program may occur based on successful or unsuccessful achievement of target limb contact and/or a specific wiping limb configuration. These observations are compatible with the kinematic observations reported in Giszter et al. (1989)
We also examined the EMGs recorded during isometric wiping trials to assess whether the ankle restraint evoked a correction response during wiping. The ankle restraint was placed in contact with the same skin areas (~2-3 mm above the ankle) upon which obstacle collision occurred during the free limb trials. We found that the same hip flexor-related muscles that were activated after obstacle collision during free limb wiping were activated during isometric wiping in our tested frogs. In the frog shown in Figure 3C,D, these hip flexor muscles were activated during the initial knee flexor phase, ~90 msec after ST onset and 65 msec after force onset (average onset time of 55 ± 11 msec relative to first measured force onset for the nine frogs taken together). The evoked EMG activity always ended close to the onset of the hip extensor phase. In fact, in eight of nine frogs of this study, the duration of activation of these muscles during isometric trials was similar to their duration of activation during collisions in free limb trials, e.g., the averaged duration of the evoked IL EMG during free limb trials was 154 ± 16 msec, whereas in isometric trials it was 165 ± 20 msec (Table 2) (p > 0.1).
In addition to the EMG similarities, correction responses initiated during isometric wiping could be abolished by cutaneous deafferentation (Fig. 3D), similar to free limb trials. The deafferentation eliminated pertinent skin feedback caused by contact with the ankle restraint. In the frog shown in Figure 3D, the IL, SA, and GL EMGs no longer exhibited a strong, phasic activation after cutaneous deafferentation. However, the knee flexor and hip extensor phases of the wiping motor pattern were clearly unaffected and had similar onset and offset times. An examination of the derived hip and knee torques generated during isometric trials showed the torque patterns produced when cutaneous feedback was intact and when this feedback was eliminated by the deafferentation procedure (Fig. 3C,D, bottom panels). With feedback intact, a hip flexor torque was initiated soon after the initial knee flexor torque and was associated with the onset of IL, GL, and SA. After removal of cutaneous feedback, this strong hip flexor torque was absent, as were the enhanced activities of IL, GL, and SA. The wiping motor pattern also involved a substantial z-direction force (Fig. 3C,D). Variations in the z-force associated with corrections were small (Fig. 3, compare C and D, Fz).
In summary, our analysis showed that wiping motor patterns can be initiated in the isometric limb, and these differ primarily only in amplitude from free limb EMGs. The first two phases of the motor pattern, the initial knee flexor (ST and BI) and subsequent hip extensor (RI, SM and AD) phases, had similar onset times as in the free limb. A different report (Giszter et al., 2000
In the following sections we exploit our observation that wiping and correction motor patterns appear conserved in the isometric limb. We used this to test the hypothesis that on-line superposition of force fields underlies the generation of trajectory corrections during wiping. We test this first at a single limb position and then at multiple limb positions.
Superposition of isometric endpoint forces underlies trajectory corrections
We examined the endpoint forces generated during wiping trials with cutaneous feedback intact, i.e., when correction responses were evoked, and after cutaneous deafferentation at a single, standard limb position (hip 90°, knee 105°) in nine frogs. In each frog, when cutaneous feedback was intact we observed a temporal pattern of endpoint forces similar to those shown in Figure 4B. We plotted both the angular directions and amplitudes of force vectors over time. Figure 4C shows the directions, expressed as an angle, over time in each condition. Figure 4D shows the magnitudes. We also expressed the angular variations of force direction as distributions in circular histograms (Fig. 5). In the frog shown in Figure 4B, the force vectors generated by the wiping limb after evoking the wiping reflex were initially directed toward the midline (~160°), where target limb contact normally occurs, and remained in this direction for ~75 msec (Fig. 4B, WCF). The forces then began to rotate toward the head of the frog over the next 150 msec. This continuous rostral rotation immediately followed the activation of IL, a hip flexor, and SA, a hip/knee flexor (averaged period of activations shown in Fig. 4E as bars). At the onset of the hip extensor burst (RI), the endpoint forces rotated in an opposite direction, away from the frog, and back toward and through the target region. They stabilized at a final direction of ~260° (measured counterclockwise and relative to horizontal, 0°) (Fig. 4C).
View larger version (25K):[in this window]
[in a new window]
Figure 4. Cutaneous reflex pathways recruited a fixed-direction force component that summed with an underlying sequence of forces generated during wiping. In this figure, time is referenced to stimulus onset. A, Force vectors in B are plotted with respect to the frog shown. Thus, vectors pointing to the left of the page are pointing to the midline of the frog. The force vector shown here represents the scale for B. Vector directions in C and Figures 5 and 6 are referenced relative to the axes in this panel. B, Isometric force vectors are shown at 20 msec intervals after initiation of wiping for two conditions (top two rows; WCF, with feedback intact; NCF, no cutaneous feedback). Forces in each condition represent an average of seven trials. In a correction response (WCF), force vectors were initially directed to the midline, rotated rostrally toward the hip, and then rotated back to and through the midline to caudal directions. In the absence of a correction (NCF), force vectors were initially directed to the midline for 200 msec and then rotated caudally. In the third row, vectors, C, produced by a subtraction of NCF vectors from WCF vectors are shown. These vectors represent the correction response. Note that these appear uniformly directed. C, Angular directions (see A) of the force vectors for the two conditions and for the correction vectors are shown at 4 msec intervals. In all trials (WCF and NCF), force vectors were initially similarly directed to 180°. In corrections (WCF), force vectors rotated to ~140°, whereas NCF vectors remained fixed in direction. In both cases (WCF and NCF), force vectors then rotated caudally at around the same time to a direction of 260°. Corrective, C, vectors as noted in B remained stable in time. D, Vector magnitudes are shown over time. With no correction (NCF), force vectors exhibited a biphasic pattern of growth. Each phase is associated with a period of fixed-direction force production as shown in B. In corrections (WCF), initial phase of vector growth is increased because of the summation of C vectors along a similar direction. The second phase is decreased because correction and wiping vectors are oppositely directed. C vectors show a monotonic increase and decline. E, Bar plots representing the averaged period of activation of ST, RI, IL, and SA EMGs are shown. The horizontal bars at each end represent ±1 SD. With feedback intact, the onset of IL and SA activation was associated with the onset of the rostral force rotation. After cutaneous deafferentation, IL and SA activation were absent, as was the rostral rotation of forces. RI onset was similar in both conditions and associated with the onset of caudally rotated forces.
View larger version (37K):[in this window]
[in a new window]
Figure 5. Corrective force components are directionally stable. A-C, Circular histograms of the directions of force vectors generated over wiping in five trials (600 msec duration trials, 250 Hz sampling rate; thus ~750 samples total) are shown for each condition: with feedback and correction response (A) and without cutaneous feedback or a correction response (B). The derived corrective vectors are shown in C. Each bin represents 12°. See D for orientation of the circular histograms to the frog's body and wiping limb. A, With feedback intact, vector directions exhibited two clusters: an initial cluster directed toward the midline and a second cluster directed away from the frog. The first cluster was broadly spread and spanned four bins. B, Without cutaneous feedback, the initial cluster of vector directions was more focally distributed to the midline, i.e., mainly in one bin (78-90°), whereas the second cluster of forces was similar to that with feedback intact (252-264°). C, Corrective vectors were obtained by subtracting time-matched forces in B from A. The directions of the corrective vectors clustered tightly in a single bin from 108 to 120° and were directed to the hip joint. Thus, the corrective vectors summed with the initial set of wiping vectors directed to the midline and resulted in the rostral spread of these midline-directed force clusters when feedback was intact.
After local removal of cutaneous feedback, the pattern of forces generated by the wiping limb consisted of only two sequentially activated and uniformly directed force components (Fig. 4B, NCF). These could correspond to stable force fields in each phase. In the absence of a corrective response, the initial forces after the wiping reflex was evoked grew in magnitude along a uniform direction (~160°) for ~200 msec (Fig. 4B,C). As the hip extensor burst was activated, the forces rotated away from the frog and reached a second set of uniformly directed forces (~260°) for the remainder of the period of RI activation. There was no period of IL/SA activation and no corresponding rostral rotation of endpoint forces without cutaneous feedback.
We analyzed the distribution of the force directions expressed as angles (as in Fig. 4C) that were generated during wiping trials both with and without cutaneous feedback intact. After cutaneous deafferentation, i.e., no correction response, the force directions clustered into two discrete subgroups (Fig. 5B). The forces in the initial (midline-directed) subgroup were contained mainly within a single 12° bin directed to 168-180° (mean direction of 175 ± 5 SD). Forces in the second (caudally directed) subgroup were mainly contained within two bins directed to 252-276° (mean direction of 261 ± 7 SD). When feedback was intact (Fig. 5A), the initial subgroup of forces directed toward the midline was more broadly distributed than the corresponding initial subgroup of forces produced when cutaneous feedback was removed (average direction of 156 ± 13 SD) (Fig. 5A). The subsequent caudal subgroup of forces with feedback intact was similarly directed to the corresponding caudal subgroup produced with cutaneous feedback removed, although the variance was greater (mean direction of 258 ± 13 SD).
We hypothesized (Eq. 10; also see Eq. 4 and 5) that differences in the variance and direction of the initial subgroup of forces produced with and without cutaneous feedback intact were caused by differences in the contribution of two fixed-direction force components (i.e., of low directional variance). When cutaneous feedback was intact, a corrective force component was activated and summed through time with an initial wiping component. This summation resulted in the rostral rotation and increased variance of the midline-directed force vectors. To test this, we subtracted the time series of force vectors that were generated when cutaneous feedback was removed from the force vectors that were generated at the corresponding times with feedback intact. We examined the time-varying pattern of resultant (corrective) forces, i.e., direction and magnitude. These are shown in Figure 4B-D (series marked C). In the frog shown in Figures 4 and 5, the corrective force vectors exhibited a monophasic increase and decrease in amplitude and were uniformly directed to ~115° during the period of IL/SA activation (Fig. 4C,D). The distribution of the corrective force directions over time showed low variance (Fig. 5C). Eighty-one percent of the corrective vectors (n = 675 vectors; 5 trials × 135 force samples; 540 msec long trials) were contained in a single bin (118-130°). The mean direction was 120°, and the SD was only 3°.
To summarize, cutaneous feedback triggered the activation of a fixed-direction force component associated with a correction response. This corrective component summed with an ongoing force component produced during the initial phase of wiping. The subsequent onset of the hip extensor phase of wiping was associated with the production of a final, fixed-direction force component. Thus these data indicate that three fixed-direction force components are activated both synchronously and sequentially to generate correction responses during wiping. This applies to isometric conditions and when testing at a single configuration.
We examined how the force component recruited during correction responses varied among frogs. Polar plots of the distributions of corrective force directions are shown for four frogs in Figure 6A. The limb posture at which the wiping trials were evoked is shown in the inset of Figure 6A (90° hip and 105° knee angles) and was the same for each frog. For six of nine frogs tested, a large percentage (>80%) of the corrective forces over time (and for multiple trials) were contained in a single bin. For each of these frogs, the variance of the corrective force directions was <12°. In addition to the uniformity of corrective force directions over time for each frog, the corrective vectors were similarly directed among frogs (the modal bin for all frogs lay within a 30° segment).
View larger version (25K):[in this window]
[in a new window]
Figure 6. Cutaneous reflex pathways that signal obstacle contact recruited similar corrective force components among frogs. A, Circular histograms of the directions of corrective forces over time are shown for four frogs (obtained as described in Figs. 4 and 5 and Materials and Methods). For each frog, the collection of corrective force directions clustered mainly within one bin and the modal bins for each frog in the group were localized within a 36° segment (from 84 to 120°). B, The averaged corrective hip and knee torques are shown for each of the four frogs shown in A. Knee and hip torques are normalized to the maximal hip torque produced for each frog. For each frog, the corrective forces were produced by a multi-joint torque response. In addition, in each frog the knee to hip torque ratios were constant over time (linear up to the peak corrective force). The small variability among frogs in the direction of corrective forces was related to the magnitude of the ratio of knee flexor torque to hip flexor torque, and in some frogs (bf65) also to the direction of the knee torque component, e.g., extensor versus flexor.
The other three frogs showed deletions of the extensor phase. Extensor phase deletions will not be further pursued here but are reported in more detail in Giszter and Kargo (2000)
A hip flexion torque was primarily responsible for generating the corrective endpoint forces. Figure 6B shows the balance of corrective knee torque versus hip torque, displayed from the onset of corrective force up to the peak corrective force. In each frog, there was an almost constant ratio of hip flexion torque to knee torque over time. This constant joint torque ratio resulted in the uniformly directed corrective forces. However, the direction (flexor vs extensor) and relative magnitude of the corrective knee torque exhibited some variability among frogs (knee torques are normalized to peak hip torque in Fig. 6B). In four of six frogs, a corrective knee flexion torque accompanied the hip flexion torque. In the other two frogs, there was a small knee extensor torque component. This variability may be related to differences among frogs in the coactivation of other muscles in addition to monoarticular hip flexors (IL and GL) during correction responses. For instance, SA, a hip flexor/knee flexor, and RA, a hip flexor/knee extensor, were both commonly activated during corrective responses. The balance of these muscles' activations would then determine the resulting knee torque. We examine in more detail the EMGs underlying corrective responses (in the sections entitled Correction responses are generated by addition of a muscle synergy to the basic pattern and Muscles recruited during correction responses exhibit covariation).
Superposition of isometric force fields underlies trajectory corrections
Although the correction responses at a single location showed clear superposition of forces, this guarantees nothing about the forces generated at other configurations or the stability over time of a force-field primitive. We therefore examined the force fields generated during wiping trials with and without cutaneous feedback intact, i.e., with and without evoked correction responses. Force fields were constructed by measuring the isometric forces generated at a range of initial limb positions. Thus, force-field measurements describe the configuration-dependent effects of both muscle actions and changes in muscle activation (caused by changes in proprioceptive feedback). Conceivably, responses might differ radically across the workspace, and field structure might be complex or divergent.
In Figure 7, a time series of force fields generated by the wiping limb with cutaneous feedback intact (Fig. 7A) and after cutaneous deafferentation (Fig. 7B) are shown for one frog. The force fields were initially similar between conditions. Approximately 90 msec after the onset of force, the force fields began to exhibit apparent differences. With feedback, intact force fields rotated rostrally toward the hip joint, whereas force fields produced with no cutaneous feedback continued to grow in magnitude and remained directed toward the midline region of stimulus contact. The time points at which the force fields are shown, i.e., 130-330 msec at 40 msec intervals, are marked as successive vertical lines on an EMG record collected at a single, central location in the workspace (Fig. 7; EMG record to the right of each force-field series). The EMG record includes four muscles that represent the knee flexor phase (BI), the hip extensor phase (RI), and the corrective response (IL and RA). The alteration in force-field structure between conditions was associated with the reflex activation of IL and RA when feedback was intact, i.e., a correction response.
View larger version (56K):[in this window]
[in a new window]
Figure 7. Cutaneous reflex pathways signaling obstacle contact recruited a corrective force field that summed with an underlying force-field sequence produced during wiping. Time series of force fields are shown after initiation of wiping both with feedback intact and after cutaneous deafferentation. A typical EMG pattern of representative muscles associated with field generation is shown to the right of each time series. For the fields, a single vector represents the force generated at that position at the specified latency after wiping initiation. The black dot marks the location of the hip joint. The lightly shaded ellipse indicates the region of contact with the target limb in free limb trials. A, With feedback intact, force fields were initially directed to the midline where target limb contact normally occurred (bold line overlaid on the 130 msec field shows the unrestrained ankle path up to limb contact. At 210 msec up to 330 msec, force fields rotated rostrally toward the hip and away from the target region. At 330 msec, force fields began to rotate back to the target region. The hip extensor burst began around this time (see EMG record to the right; rectified EMGs from BI knee flexor, RI hip extensor, and IL/RA-hip flexors are shown for a single trial taken from the middle of the workspace). B, After removal of cutaneous feedback, initial force fields were directed to the midline. Subsequent fields (210-330 msec) converged to the target region and did not rotate to the hip joint. The strong IL/RA activation was absent without cutaneous feedback (see EMGs to the right). At 330 msec, the field began to rotate caudally because of hip extensor activation. C, Procedure used to produce corrective force fields. Time-matched fields in B were subtracted from corresponding fields in A. The procedure was performed at each time point after wiping initiation. D, Corrective fields are shown at three times, at 40 msec intervals. Corrective fields were structurally similar (similarity measures > 0.95) over a 100 msec period starting shortly after IL EMG onset. These rostrally directed fields summed with the underlying wiping force fields.
We hypothesized that the alteration in force-field structure was caused by the reflex activation and superposition of a structurally invariant corrective field on a wiping field. To examine the time-varying structure of the corrective force fields, we subtracted vector fields between conditions (feedback intact fields
feedback removed fields; procedure graphically represented in Fig. 7C). In Figure 7D, corrective fields are shown at three successive time points (130, 170, and 210 msec after force onset). We used a correlation method to examine the similarity of corrective force fields over a 200 msec interval (IL and RA were activated for ~200 msec) starting at 90 msec after force onset [see Materials and Methods, Mussa-Ivaldi et al. (1994)
In summary, we addressed configuration-based variations in forces and corrective responses by measuring entire force fields. Examination of these force fields showed that trajectory corrections were generated by a dynamically evolving force-field structure. Our subtraction procedure showed that this dynamic force-field evolution was generated by combining the force-field primitive evoked in the correction response with the force-field primitives associated with the wiping reflex.
Correction responses are generated by addition of a muscle synergy to the basic pattern
The data above suggest that cutaneous feedback triggered the activation of a force-field primitive. This corrective force field sums with and does not disrupt the underlying wiping force-field sequence. In this section, we examine directly whether muscle activation triggered by obstacle collision sums with and does not disrupt the ongoing wiping motor pattern.
We found that two to three hip flexor-related muscles were activated in response to obstacle collision in each of the nine frogs tested in free limb conditions. The activation of these muscles was added to and did not interfere with the underlying motor pattern (Table 2; see list of reflexly activated muscles). In Figure 8A,B, rectified, filtered, and averaged EMGs (n = 11 trials per condition) are shown for one frog during unobstructed, free limb wiping trials (bold lines, thin lines, 1 SD above and below average, respectively) and for obstacle collision trials, respectively (similar convention). The obstacle was placed next to but not touching the wiping limb. Collision occurred immediately after movement onset. Both movement onset and obstacle collision are indicated by the dotted vertical line. SA, GL, and RA muscles that all generate a hip flexor torque (Loeb et al., 2000) were activated ~20 msec after obstacle contact (marked by the dotted vertical line). Reflex activation of these muscles during the correction appeared not to interfere with the timing or magnitude of the knee flexor phase (characterized by ST), hip extensor phase (characterized by RI), or knee extensor phase (characterized by VE) that comprised the wiping motor pattern (Fig. 8A-C). We directly examined this by subtracting the individual muscles' averaged EMG waveforms in trials without obstacle collision from averaged EMGs obtained in trials with collision. The EMG differences are shown in Figure 8D. Only SA, RA, and GL showed significant EMG changes during obstacle collision trials. Horizontal lines mark the maximum calculated SD for the difference for each muscle between the two conditions. EMG differences at peak for SA, GL, and RA were significant (p < 0.05).
View larger version (47K):[in this window]
[in a new window]
Figure 8. Evoked muscle activities during correction responses sum with the underlying wiping motor pattern. On obstacle trials, the obstacle was placed next to but not touching the wiping limb. Collision occurred immediately after movement onset. A, Rectified, filtered, and averaged EMGs for unobstructed wiping trials. Averages were constructed from 10 trials in each condition. Lines indicating 1 SD above and below the mean value are shown (fine lines). The first vertical line (dotted) marks the averaged time of movement onset and obstacle contact, and the second vertical line (solid) marks the averaged time of target limb contact for both conditions. B, EMGs for obstacle collision trials are shown. EMGs were calculated as in A. C, The EMGs for both conditions are overlaid (bold lines/no shading: unobstructed trials; hairlines/shaded envelopes: obstacle collision trials). D, The EMG differences obtained by subtracting the unobstructed EMGs from the collision trial EMGs are shown. Horizontal lines mark the maximum SDs for the EMG differences between conditions (see Materials and Methods). Obstacle contact evoked the synchronous activation of hip flexor-related muscles SA, GL, and RA without significantly disrupting other EMG bursts of the wiping motor pattern. E, The peak EMGs of GL, SA, and RA muscles recruited during the correction response are plotted versus one another for multiple correction trials (n = 10). Peak EMGs were significantly correlated (all p < 0.01) over multiple trials for each of the three muscles. F, The normalized SA EMG is plotted versus the normalized GL EMG over time from the onset of the correction response up to peak activation for multiple trials (n = 10). The EMG activation waveforms up to peak EMG were highly correlated (p < 0.001).
Muscles recruited during correction responses exhibit covariation
The extracted corrective fields described earlier (see sections entitled Superposition of isometric endpoint forces underlies trajectory corrections and Superposition of isometric force fields underlies trajectory corrections) were structurally invariant over time. Structurally invariant fields are produced either by activation of single muscles or by the synchronous and linearly covarying activation of a group of muscles. We examined EMG patterns related to the correction. These were extracted as described in the preceding section (Fig. 8D). We analyzed these to confirm that muscle activation covaried as expected. Such a pattern of activation could result from a common motor drive that is distributed to a group of motor pools. In the following section, we examined whether muscles recruited during correction responses exhibited correlated activation, which could be associated with a common motor drive.
We correlated the SA and GL EMG waveforms from the onset of the corrective response up to their peak amplitude for the 11 correction response trials (Fig. 8F). The EMGs were normalized so the peak amplitude for each trial and for each muscle was equal to one. The GL to SA EMG correlation was highly significant (p < 0.001), and the r2 value was 0.85. This high correlation over time shows that the recruitment of muscles comprising the correction response was consistent with a common motor drive. SA and GL exhibited quite different EMG waveforms during unobstructed trials, and the muscles were recorded at distant sites (ventral/distal thigh, SA; dorsal/proximal thigh, GL). The high correlation of these muscles in correction trials could not therefore be attributed to EMG cross-talk or to their very similar mechanical actions.
Our normalization procedure in the preceding analysis showed that muscles recruited in a correction response exhibited within-trial coupling. This does not guarantee an invariant gain or coupling ratio of muscle activations between trials. However, we found that the peak amplitude of the non-normalized EMG of recruited muscles covaried between correction trials. This indicated that there was little trial-to-trial variation in the coupling of muscles recruited in a correction response. The peak EMG amplitudes of SA, GL, and RA, the muscles activated during the correction response, are plotted versus one another in three separate panels in Figure 8E. The peak EMG (within a 100 msec time window after obstacle contact; thus the knee extensor burst in RA was not considered in calculation of the peak) of an individual muscle varied from trial to trial. However, the peak EMGs of SA, GL, and RA covaried together and were significantly correlated (p < 0.01; r2 values ranged from 0.80 to 0.90).
In summary, EMG analyses showed that the several muscles that were recruited during the correction response linearly covaried through time and behaved as a regulated group. Thus, in the frogs tested here, muscle groups and not single hindlimb muscles produced the fixed-structure corrective fields.
Correction response amplitude is regulated based on phase and impact force
We examined what factors might regulate the peak amplitude and duration of the correction response from trial to trial. We found that the peak amplitude of the corrective motor response depended on the magnitude of obstacle collision, i.e., the contact force. In Figure 9A, rectified and filtered SA EMGs evoked during successive obstacle collision trials are shown plotted on top of one another. Each wiping trial started from the same initial limb position with the path obstacle at the same distance (~2 mm) from the ankle. Obstacle collision occurred at a similar time for each trial, within 16.67 msec of movement onset, and collisions evoked correction responses at a similar time. There was a 15 msec range of variation in the onset of the initial EMG response. We correlated the peak of the SA EMG evoked by obstacle contact with the highest initial contact force exerted against the obstacle before the onset of the SA EMG. For most frogs (five of six) as shown for the frog in Figure 9B, there was a significant (p < 0.05) correlation between initial rise of contact force before the corrective response was evoked and the peak of the ensuing corrective response (here represented by the SA EMG). However, the duration of the corrective motor response showed no correlation with contact force. The total range of variation in the duration of the corrective response was ±5 msec. Variations in corrective response amplitude could be caused by variations in the magnitude of the feedback signal related to strength of impact. However, we expect that contact force in this experiment will vary with the magnitude of the initial knee flexor burst. Fluctuations in the overall strength of the motor program could cause variations in the strength of this burst. If corrective response strength covaried with overall motor program strength, this could therefore also account for the covariation of corrective response amplitude and contact force.
View larger version (39K):[in this window]
[in a new window]
Figure 9. Correction response amplitude was regulated both by motor phase and by feedback arising from obstacle contact. A, Rectified and filtered SA EMGs are overlaid for multiple trials (n = 11) where obstacle collision occurred at a fixed latency. For each trial, obstacle collision occurred within 16.67 msec after movement onset (the mean time of movement onset and the SDs from the mean are shown as bold and dotted vertical lines, respectively). Collision occurred slightly after the second dotted line. Evoked responses represent scaled versions of one another. B, The peak magnitude of the evoked SA EMG is plotted versus the highest contact force before SA onset as the limb collided with the obstacle. The peak SA EMG was significantly correlated with the highest preceding contact force (p < 0.05). This relationship could indicate regulation of corrective strength by contact force or covariation of corrective strength with wiping motor pattern strength. C, SA EMG duration at half peak showed no correlation with the magnitude of obstacle collision. D, In a second frog held under isometric conditions, we generated impacts of varying strength using a stepping motor to accelerate a mass into the normal site of obstacle contact. Rectified and filtered peak IL EMG (n = 34 trials) significantly covaried with the applied force (p < 0.01) independently of motor program strength, strength of correction, and time of collision. E, The rectified and filtered EMGs of evoked IL response are shown for multiple trials (n = 12) in which obstacle collision occurred at different distances (and therefore times) along the movement path of the wiping limb. The IL EMG was evoked at progressively longer latencies from movement onset because contact occurred at later times. F, The peak magnitude of the evoked IL EMG is plotted versus the time of obstacle contact. Successful responses ( ; see Results). Evoked IL responses had reduced amplitudes in these trials, departing from the linear relationship with impact force. G, The duration of evoked IL responses showed no correlation with the time of obstacle contact and varied by <10 msec. H, The contact force showed a significant linear correlation with the time of obstacle contact (p < 0.001).
In three frogs, we directly tested whether the magnitude of contact force affected the amplitude of the ensuing correction response by varying contact force independently of the motor program. In these frogs, we held the wiping limb immobile by placing bone pins in the distal portion of the tibia and attaching this arrangement to the force transducer (see Materials and Methods). Thus, skin contact did not occur with any restraint device during the wiping reflex, yet the limb was held immobile. To evoke correction responses in this condition, we used a stepper motor to drive a mass into the heel of the wiping limb with a known velocity. We varied the speed of the collision with the limb and measured the force transient produced by the impact. The times at which the mass collided with the limb were randomly varied within a 100 msec time window starting at the onset of the wiping motor pattern. We found a significant correlation between the peak impact force produced by the accelerated mass and the ensuing peak correction response. In Figure 9D, data are shown for one frog (n = 34 trials). Peak impact force varied between 0.5 and 1.0 N, and the peak-evoked IL EMG was significantly (p < 0.01) correlated with this force. The correlation between the peak IL EMG and the time of mass impact was not significant in this experiment. However, we suspect that the correlation between amplitude of the correction response and peak impact force would be higher if the time of mass impact was also taken into account. In the next section we show that motor phase, in addition to impact force, does in fact regulate the amplitude of the correction response.
In five frogs, we systematically varied the position of the path obstacle so that collisions occurred at different times in the motor pattern. Collision times varied from just after movement onset, i.e., near the onset of the knee flexor burst, to just before target limb contact, i.e., near the middle or end of the hip extensor burst and onset of the knee extensor burst. The onset of limb motion was typically 45-65 msec after ST onset, and limb contact occurred ~200-250 msec after movement onset. We found in each of these frogs that there was a time window (100-150 msec in duration) in which correction responses could be evoked. This window spanned from the time of movement onset up to the time of RI onset, i.e., the hip extensor burst.
Figure 9E shows a set of rectified and filtered IL EMGs evoked in this type of experiment. Correction responses were evoked by placing obstacles at successively further distances from the initial position of the ankle. The evoked corrective responses shown in Figure 9E had a duration that varied as little as 10 msec [duration measured at half peak (Fig. 9G)]. Thus, the IL EMGs shown in Figure 9E may represent time-shifted and similar duration (~120 msec) EMG waveforms. Figure 9F shows that these IL EMGs were evoked when obstacle contact occurred between the time of movement onset (the first set of vertical lines denote average and SDs of movement onset) and RI onset (second solid vertical line), a time window of ~100 msec. However, not all of the evoked correction responses successfully cleared the obstacle. In particular, when obstacle contact occurred close to RI onset, the corrective response often failed (Fig. 9F-H, circles). The IL EMG magnitude began to decline in these unsuccessful trials, independently of contact force, which continued to increase (Fig. 9H). Presumably in addition to this decline, the corrective hip flexion was made ineffective by the onset of the relatively powerful hip extensors and increased opposing friction at the obstacle [caused by increased obstacle contact forces (Fig. 9H)]. Once again, the amplitude of successful correction responses was significantly correlated with the peak contact force (Fig. 9F,
In summary, the corrective response is produced by a constant duration (and possibly common) motor drive to multiple motor pools. Both cutaneous feedback and motor phase can regulate the amplitude of the corrective response, but neither affects the duration. The invariance of response duration is explored in more detail in the following section.
Constant timing of corrective response and other EMG elements of the wiping motor pattern
In Figure 9, we showed that corrective responses could be evoked at different phases in the wiping motor program with little variation in the duration of the response. Despite constant duration of the EMG measured at half peak, it remained possible that the shape of EMG waveforms underlying correction responses varied depending on phase and amplitude of response. To examine this, we aligned the peaks of individual corrective EMG responses [shown in Fig. 10A (IL EMG)] and normalized their magnitudes at this peak. We then obtained the averaged EMG waveforms and their variances [for example, Fig. 10B (IL waveform)]. We found that when the waveform of activation was scaled in this way in amplitude alone, EMG waveforms of different corrective responses were very similar in shape and displayed low variance. We examined our filter characteristics (phase and frequency response) and averaging procedures, and these cannot account for the low variance of the waveform shape and duration reported here.
We tested whether the conserved shape of the activation waveform was a unique feature of the corrective response. To do this, we applied the same method of time alignment of EMG peaks and amplitude normalization to other muscles' EMGs. We chose muscles that were activated in different phases of the wiping motor pattern and not associated with corrections (for example, RI and ST). We found that these other activation waveforms also exhibited shapes (Fig. 10B, RI and ST) very similar to those observed in correction responses. When the averaged waveforms of different muscles were overlaid (Fig. 10B, last row), it could be seen that they all lay within one SD of the IL corrective waveform. We repeated this analysis procedure in a series of frogs. A second example is shown in Figure 10C. Activation waveforms of different muscles in different frogs showed similar shapes and remarkable invariance [Fig. 10D (EMG overlay for frogs bf89 and bf64)]. For certain muscles, [for example in Fig. 10, B (RI) and C (ST)], we were forced to infer inactivation shape because of overlapping phases of activity (for example, ST was commonly activated in both the knee flexor and hip extensor phases of wiping). However, in different frogs, the waveform of these same muscles could be observed directly (Fig. 10C, RI). In general, these data support the idea that different coactivated muscle groups showed very similar patterns of activation through time. We speculate that the conserved activation dynamic described here may be an important feature of the operation of motor primitives recruited in reflex behaviors and pattern generation in the frog.
View larger version (51K):[in this window]
[in a new window]
Figure 10. The activation waveform of EMGs comprising the correction response had a shape and duration that were similar to the activation waveforms of other muscles' EMGs in different phases of the wiping motor pattern. A, In the top two rows, the averaged EMG waveforms of ST and RI from movement onset up to target contact are shown. In the bottom row, the individually evoked IL responses are shown (from Fig. 9). IL responses were evoked by obstacle contact that occurred at different times in the interval between movement onset and RI onset. This was a 100 msec time window. The averaged time of movement onset is marked by the first set of vertical lines; dotted lines represent 1 SD about this. B, The averaged EMG waveforms of ST, RI, and the corrective IL burst during trajectory correction trials (n = 11) for the data in A are shown in the top three rows. Averages were constructed by aligning individual EMGs from each trial (n = 11) at its initial peak value (ST and RI exhibit multiple bursts in the wiping motor pattern) and normalizing this peak to a value of 1 (see Results for more detail) before averaging the resultant waveforms. SDs are shown around the averaged EMGs. The dotted lines indicate the inferred relaxation of the activation waveform in the absence of the later, overlapping phases (note however that in frog C, the relaxation of RI can be directly observed). These were obtained by assuming waveform symmetry and reflecting the rising phase of the EMG. In the bottom row, the three averaged EMG waveforms for RI, ST, and IL are overlaid. These waveforms have a similar shape and duration despite occurring at different phases of the motor program, and all lie within 1 SD of the IL waveform. C, The same averaging procedure was performed in a second frog. Individual EMG waveform averages were similar and showed low variance. Overlaid EMG waveforms have a similar shape and duration. D, The ST, RI, SA, and IL EMG waveforms for the two frogs shown in B and C are overlaid. The EMG waveforms in these spinal frogs all exhibit similar shape and duration. The data suggest that different motor primitives may all have common activation dynamics even when recruited in different phases and combinations during reflex behaviors.
DISCUSSION
Overview
Construction of limb movements by summation of force-field primitives was suggested by previous work. This paper examines this mechanism of movement generation and provides the first analysis of trajectory correction responses during hindlimb wiping in spinal frogs.
We showed that cutaneous feedback critical to obstacle detection acted on-line to trigger activation of a correction response during free limb wiping. We exploited this property to test in the immobilized limb whether correction responses were produced by on-line summation of force-field primitives. A force-field primitive was observed here as an isometric force field that exhibited a conserved structure over time (Giszter et al., 1993
Thus, our data are consistent with the view that lower vertebrates construct and adjust certain hindlimb movements by sequencing and summing force-field primitives that each exhibits a common activation waveform. In addition, our data indicate that vector summation, which was demonstrated in dual microstimulation experiments, occurs naturally and describes the effect of combining force-field primitives during real behavior. We discuss these findings in greater detail below.
Spinally organized correction responses
The spinal cord of lower vertebrates can generate targeted scratching or wiping movements to different skin areas (Stein, 1983
Value and limitations of force-field descriptions
Primitives were observed as isometric force fields that exhibited a conserved structure over time. Force fields map the limb's configuration into isometric forces at the ankle. Force fields were measured by holding the limb immobile in some location and recording the isometric forces produced after evoking wiping. This measurement was repeated at each of an array of positions, for the same stimulation. In this way a sampled map of force-position pairs was constructed. This map is a summary of the endpoint effects of the mechanical properties of the activated muscles in isometric conditions, limb configuration, and the effects of neural feedback pathways operating under isometric conditions. The isometric description summarizes the output of the primitive under defined conditions. In part this may help to predict the motion of the leg and its interactions with the environment. However, in free-limb conditions, limb velocity and motion-related feedback must also be considered, position feedback is altered, muscles operate in different ranges of their force-velocity curves, and inertial interactions among limb segments must be accounted for. In the frog, data support an important role of proprioceptive feedback in the organization and control of wiping trajectories (Kargo and Giszter, 2000
We focused on an isometric description of primitives to test the hypothesis that force-field summation underlies trajectory adjustments. As outlined in the preceding paragraph, the isometric description provides only part of what is needed to fully understand the structure and regulation of field primitives and their superposition in trajectory generation. However, in the behavior observed here, the isometric description was directly relevant to the function of the correction response. During free limb trials, collision with the path obstacle brings the limb to rest and initially prevents further limb movement. Thus, as a result of the collision, the limb enters isometric conditions. The isometric measurement therefore matches the initial mechanical effects of the obstacle and directly represents the forces generated by the frog that serve to accelerate the arrested limb in a new direction to clear the obstacle.
Force-field summation during real behavior
In past work, microstimulation of premotor regions of the spinal cord elicited only a few types of force fields (Bizzi et al., 1991
Combination of primitives does not alter the normal motor program
It was unexpected that despite the addition of the corrective response, the timing of phases comprising the wiping motor pattern was not affected in obstacle collision trials when compared with unperturbed trials. This was true when the motor program was examined at the level of free limb EMGs, isometric EMGs, or force production. Despite the lack of adjustment in the timing of the wiping motor pattern, the time to target limb contact was preserved and similar to the times to contact seen in unimpeded trials. Kargo and Giszter (2000)
A unit of fixed temporal dynamics
How force-field primitives are recruited, modulated, and combined to construct movements could occur in several ways. First, a primitive could represent a simple feedforward synergy that passively follows arbitrary activation patterns. The tasks of planning and execution systems then involve determining the activation waveforms and feedback controls (Mussa-Ivaldi and Giszter 1992
Voluntary human movement exhibits a wide temporal range, which may vary depending on task. However, several investigators have provided support for schemes of voluntary trajectory generation involving superposition of pulsed or fixed timing elements (Vicario and Ghez, 1984
FOOTNOTES Received Aug. 19, 1999; revised Oct. 25, 1999; accepted Oct. 28, 1999.
This work was supported by National Institutes of Health (NIH) Grant NS43460 to S.F.G., and NIH Grant 5T32 HD07467. We thank Drs. Donald Faber, Marion Murray, Alberto Pereda, Richard Poppele, Terence Sanger, and John Soechting for comments on earlier versions of this work. Nicholas Auyong designed the driver circuit for the stepper motor system used in part of this work.
Correspondence should be addressed to Simon Giszter, Department of Neurobiology and Anatomy, Medical College of Pennsylvania/Hahnemann Medical School, 3200 Henry Avenue, Philadelphia, PA 19129. E-mail: simon{at}swampthing.auhs.edu.
REFERENCES
- Asada H, Slotine JJE (1986) In: Robot analysis and control. New York: Wiley.
- Basmajian JV, DeLuca CJ (1985) In: Muscles alive: their functions revealed through electromyography, Ed 5. Baltimore: Williams & Wilkins.
- Berkinblit MB, Feldman AG, Fukson OI (1986) Adaptability of innate motor patterns and motor control mechanisms. Behav Brain Sci 9:585-638[Web of Science].
- Bizzi E, Giszter SF, Mussa-Ivaldi FA (1991) Computations underlying the execution of movement: a novel biological perspective. Science 253:287-291
[Abstract/Free Full Text] .- Bizzi E, Giszter SF, Loeb E, Mussa-Ivaldi FA, Saltiel P (1995) Modular organization of motor behavior in the frog's spinal cord. Trends Neurosci 18:442-446[Web of Science][Medline].
- Buford JA, Smith JL (1993) Adaptive control for backward quadrupedal walking. III. Stumbling corrective reactions and cutaneous reflex sensitivity. J Neurophysiol 70:1102-1114
[Abstract/Free Full Text] .- Burdet E, Milner TE (1998) Quantization of human motions and learning of accurate movements. Biol Cybern 78:307-318[Web of Science][Medline].
- Cannon M, Slotine JJE (1995) Space frequency localized basis networks for nonlinear estimation and control. Neurocomputing 9:000-000.
- d'Avella A, Bizzi E (1998) Low dimensionality of supraspinally induced force fields. Proc Natl Acad Sci USA 95:7711-7714
[Abstract/Free Full Text] .- Doeringer JA, Hogan N (1998) Intermittancy in preplanned elbow movements persist in the absence of visual feedback. J Neurophysiol 80:1787-1799
[Abstract/Free Full Text] .- Drew T, Rossignol S (1987) A kinematic and electromyographic study of cutaneous reflexes evoked from the forelimb of unrestrained walking cats. J Neurophysiol 57:1160-1184
[Abstract/Free Full Text] .- Forssberg H (1979) Stumbling corrective reaction: a phase-dependent compensatory reaction during locomotion. J Neurophysiol 42:936-953
[Abstract/Free Full Text] .- Giszter SF, Kargo WJ (2000) Conserved temporal dynamics and vector superposition of primitives in frog wiping reflexes during spontaneous extensor deletions. Neurocomputing, in press.
- Giszter SF, McIntyre J, Bizzi E (1989) Kinematic strategies and sensorimotor transformations in the wiping movements of frogs. J Neurophysiol 62:750-767
[Abstract/Free Full Text] .- Giszter SF, Mussa-Ivaldi FA, Bizzi E (1993) Convergent force field organized in the frog's spinal cord. J Neurosci 13:467-491[Abstract].
- Giszter SF, Kargo WJ, Davies MR (2000) Augmenting postural primitives in spinal cord: dynamic force-field structures used in trajectory generation. In: Biomechanics and neural control of movement (Winter J, Crago PJ, eds). New York: Springer, in press.
- Gottlieb GL (1998) Muscle activation patterns during two types of voluntary single-joint movement. J Neurophysiol 80:1860-1867
[Abstract/Free Full Text] .- Grillner S, Wallen P (1985) Central pattern generators for locomotion, with special reference to vertebrates. Annu Rev Neurosci 8:233-261[Web of Science][Medline].
- Kargo WJ, Giszter SF (2000) Afferent roles in hindlimb wipe reflex trajectories: free limb kinematics and motor patterns. J Neurophysiol, in press.
- Kargo WJ, Davies MR, Giszter SF (1998) Segmental afferent regulation of hindlimb wiping in the spinal frog. Ann NY Acad Sci 860:456-457[Web of Science][Medline].
- Kiehn O, Hounsgaard J, Sillar K (1997) Basic building blocks of vertebrate spinal central pattern generators. In: Neurons, networks, and motor behavior (Stein P, Grillner S, Selverston A, Stuart D, eds), pp 47-60. Boston: MIT.
- Krebs HI, Aisen ML, Volpe BT, Hogan N (1999) Quantization of continuous arm movements in humans with brain injury. Proc Natl Acad Sci USA 96:4645-4649
[Abstract/Free Full Text] .- Loeb EP, Giszter SF, Saltiel P, Mussa-Ivaldi FA, Bizzi E (1999) Output units of motor behavior: an experimental and modeling study. J Cognit Neurosci 12:1-20[Web of Science][Medline].
- Loeb GE (1985) Motoneurone task groups: coping with kinematic heterogeneity. J Exp Biol 115:137-146
[Abstract/Free Full Text] .- Martin JH, Cooper SE, Ghez C (1995) Kinematic analysis of reaching in the cat. Exp Brain Res 102:379-392[Web of Science][Medline].
- Massaquoi SG, Slotine JE (1996) The intermediate cerebellum may function as a wave-variable processor. Neurosci Lett 215:60-64[Web of Science][Medline].
- Mataric MJ, Williamson M, Demiris J, Mohan A (1998) Behavior based primitives for articulated control. In: Proceedings from animal to animats 5, Fifth International Conference on Simulation of Adaptive Behavior (Pfeifer R, Blumberg B, Meyer J, Wilson S, eds), pp 165-170. Boston: MIT.
- Mortin LI, Stein PSG (1989) Spinal cord segments containing key elements of the central pattern generators for three forms of the scratch reflex in turtles. J Neurosci 9:2285-2296[Abstract].
- Mussa-Ivaldi FA (1992) From basis functions to basis fields: vector field approximation from sparse data. Biol Cybern 67:479-489[Web of Science][Medline].
- Mussa-Ivaldi FA (1997) Nonlinear force fields: a distributed system of control primitives for representing and learning movements. Proc IEEE Int Symp Comput Intell Robot Autom 84-90.
- Mussa-Ivaldi FA, Giszter SF (1992) Vector field approximation: a computational paradigm for motor control and learning. Biol Cybern 67:491-500[Web of Science][Medline].
- Mussa-Ivaldi FA, Giszter SF, Bizzi E (1994) Linear combinations of primitives in vertebrate motor control. Proc Natl Acad Sci USA 91:7534-7538
[Abstract/Free Full Text] .- Nichols TR (1994) A biomechanical perspective on spinal mechanisms of coordinated muscular action: an architecture principle. Acta Anat 151:1-13[Web of Science][Medline].
- Pellegrini JJ, Flanders M (1996) Force path curvature and conserved features of muscle activation. Exp Brain Res 110:80-90[Web of Science][Medline].
- Prochazka A, Gillard D (1997) Sensory control of locomotion. Proc IEEE Am Control Conf 5:2846-2850.
- Rossignol S (1996) Neural control of stereotypic limb movements. In: Handbook of physiology, Sec 12 (Rowell L, Shepherd J, eds), pp 173-216. New York: Oxford UP.
- Saltiel P, Tresch MC, Bizzi E (1998) Spinal cord modular organization and rhythm generation: an NMDA iontophoretic study in the frog. J Neurophysiol 80:2323-2339
[Abstract/Free Full Text] .- Sanner R, Slotine JJE (1995) Stable robotic learning using "neural networks." Neural Comput 7:753-790[Web of Science].
- Schotland JL, Rymer WZ (1993) Wipe and flexion reflexes of the frog. I. Kinematics and EMG patterns. J Neurophysiol 69:1725-1735
[Abstract/Free Full Text] .- Sergio LE, Ostry DJ (1993) Three-dimensional kinematic analysis of frog hindlimb movement in reflex wiping. Exp Brain Res 94:53-64[Web of Science][Medline].
- Simpson J (1969) On why a frog is not a cat. In: PhD thesis MIT.
- Smith JL, Stein PSG (1997) Neural and Biomechanical control strategies for different forms of vertebrate hindlimb motor tasks. In: Neurons, networks, and motor behavior (Stein P, Grillner S, Selverston A, Stuart D, eds), pp 47-60. Boston: MIT.
- Stein PSG (1983) The vertebrate scratch reflex. Symp Soc Exp Biol 37:383-403[Medline].
- Tresch MC (1997) Discreteness in spinal motor systems in the rat and frog. In: PhD thesis MIT.
- Tresch MC, Saltiel P, Bizzi E (1999) The construction of movement by the spinal cord. Nat Neurosci 2:162-167[Web of Science][Medline].
- Vicario DS, Ghez C (1984) The control of rapid limb movement in the cat. IV. Updating of ongoing isometric responses. Exp Brain Res 55:134-144[Web of Science][Medline].
- Wand P, Prochazka A, Sontag KH (1980) Neuromuscular responses to gait perturbations in freely moving cats. Exp Brain Res 38:109-114[Web of Science][Medline].
- Williamson MM (1996) In: Postural primitives: interactive behavior for a humanoid robot arm. Proceedings of the Fourth International Conference on Simulation of Adaptive Behavior (Maes P, Mataric M, Meyer J, Pollack J, Wilson S, eds), pp 124-135. Boston: MIT.
Copyright © 2000 Society for Neuroscience 0270-6474/0/201409-18$05.00/0
This article has been cited by other articles:
M. Berniker, A. Jarc, E. Bizzi, and M. C. Tresch
Simplified and effective motor control based on muscle synergies to exploit musculoskeletal dynamics
PNAS, May 5, 2009; 106(18): 7601 - 7606.
[Abstract] [Full Text] [PDF]
V. C. K. Cheung, A. d'Avella, and E. Bizzi
Adjustments of Motor Pattern for Load Compensation Via Modulated Activations of Muscle Synergies During Natural Behaviors
J Neurophysiol, March 1, 2009; 101(3): 1235 - 1257.
[Abstract] [Full Text] [PDF]
V. S. Boyce and M. A. Lemay
Modularity of Endpoint Force Patterns Evoked Using Intraspinal Microstimulation in Treadmill Trained and/or Neurotrophin-Treated Chronic Spinal Cats
J Neurophysiol, March 1, 2009; 101(3): 1309 - 1320.
[Abstract] [Full Text] [PDF]
J. J. Kutch, A. D. Kuo, A. M. Bloch, and W. Z. Rymer
Endpoint Force Fluctuations Reveal Flexible Rather Than Synergistic Patterns of Muscle Cooperation
J Neurophysiol, November 1, 2008; 100(5): 2455 - 2471.
[Abstract] [Full Text] [PDF]
E. J. Perreault, K. Chen, R. D. Trumbower, and G. Lewis
Interactions With Compliant Loads Alter Stretch Reflex Gains But Not Intermuscular Coordination
J Neurophysiol, May 1, 2008; 99(5): 2101 - 2113.
[Abstract] [Full Text] [PDF]
W. J. Kargo and S. F. Giszter
Individual Premotor Drive Pulses, Not Time-Varying Synergies, Are the Units of Adjustment for Limb Trajectories Constructed in Spinal Cord
J. Neurosci., March 5, 2008; 28(10): 2409 - 2425.
[Abstract] [Full Text] [PDF]
Y. P. Ivanenko, G. Cappellini, N. Dominici, R. E. Poppele, and F. Lacquaniti
Modular Control of Limb Movements during Human Locomotion
J. Neurosci., October 10, 2007; 27(41): 11149 - 11161.
[Abstract] [Full Text] [PDF]
A. d'Avella, A. Portone, L. Fernandez, and F. Lacquaniti
Control of fast-reaching movements by muscle synergy combinations.
J. Neurosci., July 26, 2006; 26(30): 7791 - 7810.
[Abstract] [Full Text] [PDF]
G. Cappellini, Y. P. Ivanenko, R. E. Poppele, and F. Lacquaniti
Motor Patterns in Human Walking and Running
J Neurophysiol, June 1, 2006; 95(6): 3426 - 3437.
[Abstract] [Full Text] [PDF]
M. G. Sirota, G. A. Pavlova, and I. N. Beloozerova
Activity of the Motor Cortex During Scratching
J Neurophysiol, February 1, 2006; 95(2): 753 - 765.
[Abstract] [Full Text] [PDF]
V. S. Chib, J. L. Patton, K. M. Lynch, and F. A. Mussa-Ivaldi
Haptic Identification of Surfaces as Fields of Force
J Neurophysiol, February 1, 2006; 95(2): 1068 - 1077.
[Abstract] [Full Text] [PDF]
Y. P. Ivanenko, G. Cappellini, N. Dominici, R. E. Poppele, and F. Lacquaniti
Coordination of Locomotion with Voluntary Movements in Humans
J. Neurosci., August 3, 2005; 25(31): 7238 - 7253.
[Abstract] [Full Text] [PDF]
V. C. K. Cheung, A. d'Avella, M. C. Tresch, and E. Bizzi
Central and Sensory Contributions to the Activation and Organization of Muscle Synergies during Natural Motor Behaviors
J. Neurosci., July 6, 2005; 25(27): 6419 - 6434.
[Abstract] [Full Text] [PDF]
Y. Zhurov, K. R. Weiss, and V. Brezina
Tight or Loose Coupling Between Components of the Feeding Neuromusculature of Aplysia?
J Neurophysiol, July 1, 2005; 94(1): 531 - 549.
[Abstract] [Full Text] [PDF]
A. G. Richardson, J.-J. E. Slotine, E. Bizzi, and M. C. Tresch
Intrinsic Musculoskeletal Properties Stabilize Wiping Movements in the Spinalized Frog
J. Neurosci., March 23, 2005; 25(12): 3181 - 3191.
[Abstract] [Full Text] [PDF]
A. d'Avella and E. Bizzi
Shared and specific muscle synergies in natural motor behaviors
PNAS, February 22, 2005; 102(8): 3076 - 3081.
[Abstract] [Full Text] [PDF]
J.-F. Perrier and M. C. Tresch
Recruitment of motor neuronal persistent inward currents shapes withdrawal reflexes in the frog
J. Physiol., January 15, 2005; 562(2): 507 - 520.
[Abstract] [Full Text] [PDF]
J. Jing, E. C. Cropper, I. Hurwitz, and K. R. Weiss
The Construction of Movement with Behavior-Specific and Behavior-Independent Modules
J. Neurosci., July 14, 2004; 24(28): 6315 - 6325.
[Abstract] [Full Text] [PDF]
C. B. Hart and S. F. Giszter
Modular Premotor Drives and Unit Bursts as Primitives for Frog Motor Behaviors
J. Neurosci., June 2, 2004; 24(22): 5269 - 5282.
[Abstract] [Full Text] [PDF]
Y. P. Ivanenko, R. E. Poppele, and F. Lacquaniti
Five basic muscle activation patterns account for muscle activity during human locomotion
J. Physiol., April 1, 2004; 556(1): 267 - 282.
[Abstract] [Full Text] [PDF]
W. J. Kargo and D. A. Nitz
Early Skill Learning Is Expressed through Selection and Tuning of Cortically Represented Muscle Synergies
J. Neurosci., December 3, 2003; 23(35): 11255 - 11269.
[Abstract] [Full Text] [PDF]
E. Bizzi, A. D'Avella, P. Saltiel, and M. Tresch
Book Review: Modular Organization of Spinal Motor Systems
Neuroscientist, October 1, 2002; 8(5): 437 - 442.
[Abstract] [PDF]
B. Rohrer, S. Fasoli, H. I. Krebs, R. Hughes, B. Volpe, W. R. Frontera, J. Stein, and N. Hogan
Movement Smoothness Changes during Stroke Recovery
J. Neurosci., September 15, 2002; 22(18): 8297 - 8304.
[Abstract] [Full Text] [PDF]
W. J. Kargo and L. C. Rome
Functional morphology of proximal hindlimb muscles in the frog Rana pipiens
J. Exp. Biol., July 15, 2002; 205(14): 1987 - 2004.
[Abstract] [Full Text] [PDF]
R. E. Poppele, G. Bosco, and A. M. Rankin
Independent Representations of Limb Axis Length and Orientation in Spinocerebellar Response Components
J Neurophysiol, January 1, 2002; 87(1): 409 - 422.
[Abstract] [Full Text] [PDF]
G. Bosco and R. E. Poppele
Proprioception From a Spinocerebellar Perspective
Physiol Rev, April 1, 2001; 81(2): 539 - 568.
[Abstract] [Full Text] [PDF]
P. Saltiel, K. Wyler-Duda, A. D'Avella, M. C. Tresch, and E. Bizzi
Muscle Synergies Encoded Within the Spinal Cord: Evidence From Focal Intraspinal NMDA Iontophoresis in the Frog
J Neurophysiol, February 1, 2001; 85(2): 605 - 619.
[Abstract] [Full Text] [PDF]
W. J. Kargo and S. F. Giszter
Afferent Roles in Hindlimb Wipe-Reflex Trajectories: Free-Limb Kinematics and Motor Patterns
J Neurophysiol, March 1, 2000; 83(3): 1480 - 1501.
[Abstract] [Full Text] [PDF]
This Article Abstract 
Full Text (PDF) Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager 
Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (57) Citing Articles via Google Scholar Google Scholar Articles by Kargo, W. J. Articles by Giszter, S. F. Search for Related Content PubMed PubMed Citation Articles by Kargo, W. J. Articles by Giszter, S. F.
|
|
|
|
 |
|