- 1School of Biomedical Engineering, Med-X Research Institute, Shanghai Jiao Tong University, Shanghai, China
- 2School of Technology and Health, Royal Institute of Technology, Stockholm, Sweden
- 3Department of Rehabilitation, Ruijin Hospital of School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 4Division of Biokinesiology and Physical Therapy, University of Southern California, Los Angeles, CA, USA
Mounting evidence suggests that both α and γ motoneurons are active during movement and posture, but how does the central motor system coordinate the α-γ controls in these tasks remains sketchy due to lack of in vivo data. Here a computational model of α-γ control of muscles and spindles was used to investigate α-γ integration and coordination for movement and posture. The model comprised physiologically realistic spinal circuitry, muscles, proprioceptors, and skeletal biomechanics. In the model, we divided the cortical descending commands into static and dynamic sets, where static commands (αs and γs) were for posture maintenance and dynamic commands (αd and γd) were responsible for movement. We matched our model to human reaching movement data by straightforward adjustments of descending commands derived from either minimal-jerk trajectories or human EMGs. The matched movement showed smooth reach-to-hold trajectories qualitatively close to human behaviors, and the reproduced EMGs showed the classic tri-phasic patterns. In particular, the function of γd was to gate the αd command at the propriospinal neurons (PN) such that antagonistic muscles can accelerate or decelerate the limb with proper timing. Independent control of joint position and stiffness could be achieved by adjusting static commands. Deefferentation in the model indicated that accurate static commands of αs and γs are essential to achieve stable terminal posture precisely, and that the γd command is as important as the αd command in controlling antagonistic muscles for desired movements. Deafferentation in the model showed that losing proprioceptive afferents mainly affected the terminal position of movement, similar to the abnormal behaviors observed in human and animals. Our results illustrated that tuning the simple forms of α-γ commands can reproduce a range of human reach-to-hold movements, and it is necessary to coordinate the set of α-γ descending commands for accurate and stable control of movement and posture.
Introduction
The physiological system of human motor control is not only highly redundant (Bernstein, 1967; Martin et al., 2009), but also endowed with intricate dual α and γ sensorimotor control (Granit, 1975; Pierrot-Deseilligny and Burke, 2005; Lan and He, 2012; Prochazka and Ellaway, 2012). Our understanding about how movements are organized and how muscles are coordinated in performing different tasks remains incomplete due to lack of in vivo data during behaviors. One of the remaining issues in sensorimotor control is to account for the role of γ motor system in movement and posture. In spite of thorough elucidation of peripheral efferent and afferent innervations of the spindle organ (Matthews, 1964; Boyd, 1980; Hulliger, 1984; Prochazka, 1999), the importance of γ motor system in motor control is still not well understood.
A consistently observed phenomenon during both movement and posture control is α-γ co-activation (Vallbo, 1971; Taylor et al., 2004, 2006; Prochazka and Ellaway, 2012). Direct recording by Taylor et al. (2004, 2006) revealed a co-varying pattern of gamma-static and dynamic firings with joint angle during locomotion. There was plenty evidence of independent control of γ-motoneurons during movement (Prochazka et al., 1985; Dimitriou and Edin, 2008). The co-activation of γ motoneurons with α motor activity was generally viewed to compensate for the unloading effects of muscle contraction to spindle sensitivity. Using a realistic virtual arm (VA) model (He et al., 2013), Lan and He (2012) suggested a plausible function of γs fusimotor control to convey centrally encoded joint angle information to the periphery through regulating spindle sensitivity, so that the Ia signaling from spindle afferents is kept faithfully proportional to the joint angle during movement and muscle contraction. This not only explains γ co-activation with α activity, but also supports the hypothesis that the γs command reinforces the centrally planned kinematics of movement and posture by way of spinal circuits.
Neurophysiological studies have identified separate spinal pathways and circuits of sensorimotor system in details, where the α and γ commands interact with each other to produce sensory and motor outputs (Baldissera et al., 1981; Lemon et al., 2004; Pierrot-Deseilligny and Burke, 2005). In addition to the mono-synaptic cortico-motoneuronal pathway (Lawrence and Kuypers, 1968; Lemon, 1999; Quallo et al., 2012), disynaptic excitatory and inhibitory cortico-motoneuronal pathways via PNs in C3-C4 were found to exist in cats and in macaque monkeys, as well as in humans (Malmgren and Pierrot-Deseilligny, 1988; Gracies et al., 1994; Alstermark et al., 2007; Isa et al., 2007). The PN was shown to play an important role in reaching movement of upper limb (Alstermark et al., 1981; Alstermark and Isa, 2012). In a computational analysis to emulate involuntary oscillatory movements in human upper extremity (Hao et al., 2013), it is postulated that movement signals are transmitted via the disynaptic pathway of the PN network, where the γ-dynamic (γd) command is integrated with the α-dynamic (αd) command to produce pre-motoneuronal outputs. The γd command encodes kinetic information of joint acceleration (or deceleration), and gates the αd command of double frequency at the PN network to determine the timing of activation for a pair of flexor and extensor muscles during oscillatory movements.
These studies suggest that while the α motor system provides the main drive for muscles, the γ motor system executes a more subtle control for movement dynamics and the maintenance of posture. In this paper, we extend the α-γ model (Hao et al., 2013) to investigate the modular control of voluntary movement and posture, and to demonstrate coordination of a set of α-γ descending commands in the control of movement and posture. Human reach-and-hold movements and muscle EMG activities were recorded and analyzed to guide the specification of the central descending commands. The model behavior was matched to the data of human movement and posture in all subjects. A variety of modifications were introduced in model structures to assess the functional significance of α-γ coordination by “deafferenting” or “deefferenting” the model. Results of this study provided a quantitative evaluation of the relative importance that each proprioceptive afferent and descending α-γ command may have on movement and posture control. Agreement between model predictions and human movement data further corroborated the α-γ coordination as a rudimentary means of sensorimotor control. The results argue that coordinated γ control with α activation is essential for accurate and stable control of movement and posture. Preliminary analysis of this study appeared in a conference proceeding (Li et al., 2014).
Materials and Methods
Corticospinal Virtual Arm (CS-VA) Model
The module of posture and movement control based on α-γ dual control system was implemented in the computational CS-VA model (Figure 1). This model is based on realistic physiological studies (Cheng et al., 2000; Mileusnic and Loeb, 2006; Mileusnic et al., 2006; Taylor et al., 2006; Alstermark et al., 2007; Isa et al., 2007) and has been validated to capture the realistic neuromechanical properties of human upper limb in previous work (Song et al., 2008a,b; Lan and He, 2012; He et al., 2013). It consists of four parts, the primary motor and sensory cortex, PN network, spinal cord circuitry and virtual arm (VA) model. The spinal cord network with PNs innervating a pair of antagonistic muscles (Figure 2) has been validated in the mechanism study of Parkinsonian tremor (Hao et al., 2013). The VA has two joints including shoulder and elbow in the horizontal plane, two degrees of freedom, namely shoulder flexion/extension and elbow flexion/extension. The six dominating muscles includes two pairs of monoarticular muscles, shoulder flexor Pectoralis Major Clavicle (PC), extensor Deltoid Posterior (DP), elbow flexor Brachialis (BS), extensor Triceps Lateral (Tlt) and one pair of biarticular muscles, flexor Biceps Short Head (Bsh), extensor Triceps Long Head (Tlh). Muscle spindle and Golgi Tendon Organ (GTO) are implemented within the VA to provide proprioceptive feedback. Three parts of the CS-VA model, including PN network, spinal cord circuitry and VA, have been integrated in SIMULINK/MATLAB platform. In this study, bi-articular muscles of Bsh and Tlh were mainly used to realized movement and posture.
Figure 1. Corticospinal virtual arm (CS-VA) model. Descending α and γ commands from motor cortex were processed in PN network and spinal cord circuitry, the coordinated output of α and γ motor neurons are transmitted to virtual muscles and spindles, the virtual arm is then activated. PN represents propriospinal interneuron. The subscripts “s” and “d” in α and γ motor signals stand for static and dynamic commands, respectively. Um is muscle input, which is equivalent to human EMG; Us is spindle input; Lmt and Fm represent muscle tendon length and muscle force computed form virtual arm model; Lce represents fascicle length; Ia and II are sensory feedback of primary and secondary afferents from muscle spindle; and Ib is the feedback of Golgi Tendon Organ (GTO). The virtual arm has two joints including shoulder and elbow in the horizontal plane, with two degrees of freedom, i.e., shoulder flexion/extension and elbow flexion/extension. Three pairs of antagonistic muscles are included in the model: shoulder flexor Pectoralis Major Clavicle (PC), extensor Deltoid Posterior (DP), elbow flexor Brachialis (BS), extensor Triceps Lateral (Tlt), and biarticular flexor Biceps Short Head (Bsh), extensor Triceps Long Head (Tlh). Simulated joint kinematics, especially elbow angle is compared to elbow angle (θel) recorded in human experiment. Modified from Hao et al. (2013) with permission.
Figure 2. α-γ control model in spinal cord circuitry for a pair of antagonistic muscles. α and γ commands are coordinated by PN network and reflex circuitry within spinal cord. The subscripts “e” and “f” stand for extensor and flexor; MN represents motor neuron pool; U represents muscle input; de and df are gains of γd inhibition on PN(γ-PN); pe and pf are gains of PN to reciprocal inhibition; ge and gf are gains of recurrent inhibition from Renshaw Cell (RC); be and bf represent inhibition gains on α motor neurons from Ib afferents (Ib(−)); se and sf are stretch reflex gains from Ia afferents ((Ia(+)); re and rf are reciprocal inhibition gains on α motor neurons (Ia(−)); ae and af are Ia afferent gains on PN(Ia-PN). Excitatory synapse at neurons is indicated by a “y” termination, and inhibitory synapse is indicated by a filled dot termination. The function of the feedback from Ia afferents in dashed lines are discussed later.
Within the CS-VA model, central motor commands from motor cortex are transmitted to spinal α and γ motoneurons (αMN, γMN) in two pathways, mono-synaptic pathway carrying static α and γ motor commands to corresponding motoneurons for posture control, and multi-synaptic pathway transmitting dynamic α and γ motor commands to MNs via PN network for movement control. Muscles and spindles are activated by α and γ motor signals respectively after the regulation of PN and reflex network within spinal cord. The muscular dynamics computed by SIMM model of arm give a real-time muscle tendon length to virtual muscle, at the same time, the dynamics are reflected by Ia and II afferents from muscle spindle, and Ib afferent from GTO. Feedback from proprioceptive afferents participates in signal processing within spinal cord network and cortex simultaneously. In this framework, central inputs of static and dynamic α and γ motor commands are specialized in Section Specification of Central Descending Commands, the input of virtual muscle (Um) is equivalent to recorded EMG, and the output joint kinematics gives a behavioral indicator to compare with elbow angle (θel) in human experiment.
PNs in C3-C4 receive wide descending control from cortico-, rubro-, reticulo and tectospinal tracts, and project inhibition or excitation to downstream αMN (Alstermark et al., 2007), at the same time receive feedback from Ia and cutaneous afferents (Alstermark et al., 1984a,b,c). Based on neurophysiology structures of PN network, αMN, γMN, and the spinal reflexes, a model of spinal circuitry was implemented within the CS-VA model as is shown in Figure 2 (Hao et al., 2013; Li et al., 2014). The central descending commands in the model were passed down to αMN and γMN through mono-synaptic and multi-synaptic PN pathways, respectively (Isa et al., 2007). There was experimental evidence indicating that PNs mediate motor commands for reaching movement (Alstermark et al., 2007; Alstermark and Isa, 2012). Computational analysis (Lan and He, 2012) suggested that γs conveys kinematic information of joint angle. Therefore, in this model, we assume that the descending commands can be divided into static set (αs and γs) for posture and dynamic set (αd and γd) for movement. This division is consistent to the evidence of individual posture and movement modules in central motor control system (Kurtzer et al., 2005).
As illustrated in Figure 2, α and γ commands control muscles and spindles through a set of spinal circuits. αd is integrated at the PN with γd in a pattern of mirror inhibition (γd or 1-γd). The recording from dynamic γMN in Taylor et al. (2006) revealed a reciprocal change of γd with joint angle. The PN also receives excitation from autogenic Ia afferent (Ia-PN) (Alstermark et al., 1984b; Malmgren and Pierrot-Deseilligny, 1988). After the PN, αs and αd commands converge at αMN, and its output is further regulated by reciprocal inhibition (Ia(−)) from antagonistic muscle, recurrent inhibition of Renshaw Cell (RC), autogenic stretch reflexes from Ia (Ia(+)), and Ib afferents (Ib(−)) (Eccles et al., 1957, 1960; Windhorst, 2007). The output of PN (YPN) and αMN (YαMN) are given in Equation (1), in which subscript “e” and “f” represent extensor and flexor,
N(t) is the sum of descending commands to αMN, C(t) represents background activation of αMN pools, τ is the time constant of excitation with the value 0.029 (sec), σ(t) represents the signal dependent noise, which is a Gaussian distributed random signal (Fuglevand et al., 1993; Jones et al., 2002). C(t) was multiplied by σ(t) to yield the noise corrupted background activations (C′(t)) in Equations (1g) and (1h). The output of αMN (Equations 1i and 1j) is the sum of all excitatory and inhibition inputs; here, the values of υ′ and φ′ are proportional to Ia and Ib afferent discharge frequencies, respectively; and g, s, r, b are reflex gains of RC, Ia(+), Ia(−), Ib(−). The values of spinal reflex gains are listed in Table 1.
Human Reach and Hold Experiment
Subjects and Experiments
Seven healthy adults participated in this elbow joint's reach and hold experiment. The human subject study was approved by the Internal Review Board (IRB) of University of Southern California (USC). All of them obtained a brief explanation of this study before the experiment, and signed the informed consent.
During the experiment, the subject sat comfortably with the upper arm maintained perpendicular to the trunk in the horizontal plane (Figure 3A), the hand and forearm were secured on the manipulation to make single joint movements around the elbow. All subjects initially held their elbow at 90°, and performed successive reach and hold movements triggered by the visual cue of LED light, following three blocks. In Block 1, the elbow was extended with the range of 30 degrees, and after three successive extensions, the elbow stopped at 0°, the reversal flexion with the same movement range and holding posture was performed after a 5 s holding period at 0°. In Block 2, the range of extension and flexion was 45°, thus one holding period at the mid-posture of 45° was performed, and in Block 3 the elbow had direct extension and flexion movement between 90° and 0°. During the experiments, the subjects moved as fast as them could between postures, and a 5 s holding was required at each posture. Each block was repeated for five trials, and between blocks, the subjects had more than 10 min to rest. During the movement, the elbow angle was recorded and EMGs of biceps short head (Bsh) and triceps long head (Tlh) were collected using bi-polar surface electrodes, the EMG signals were pre-amplified at the gain of 1000, band-pass filtered with cut-off frequencies at 5 Hz and 1000 Hz, and then sampled at 2000 Hz for post-processing.
Figure 3. Human reach and hold experiments, and analysis of kinematic and EMG data during posture and movement. (A) Experimental set up: the subject is seated at the table to make fast reaching movement of elbow joint in the horizontal plane. Subjects' arm and hand were secured on a manipulandum to keep the upper arm stable and ensure smooth extension and flexion of the forearm. Elbow angle (θel) is defined as the included angle between forearm and the extending line of upper arm. Movements with a range of 30, 45, and 90 degrees are performed in Block 1, 2, and 3, respectively. During extension, the elbow angle (θel) changes from initial posture of 90° to terminal posture of 0°, the filled dot indicates the holding posture between reaching movements. The reversal flexion from 0° to 90° has the same movement range and holding postures. The shoulder angle (θsh) kept at 90° during all experiments. (B) Analysis of kinematic and EMG data during posture and movement. During movement, the elbow moves from initial posture to terminal posture, and the movement onset (t0) and offset (t1) are defined as the time at which the velocity change (increases or decreases) is 10% of the peak velocity (Vpeak). Initial posture period starts at time (t0 – 1.5) and terminal posture is defined from time (t0 + 4) after the angle stabilizes. The time window is 1 s. The tri-phasic EMG pattern of extensor Triceps long head (Tlh) and flexor Biceps short head (Bsh with inverse value) is shown in the bottom (low passed at a cut off frequency of 50 Hz with digital Butterworth filter).
Analysis of Kinematics and EMG Data
Recorded joint trajectory and EMG were analyzed to provide guidance for simulation. Since Subject 3 couldn't follow the experimental protocol, the data was excluded from result analysis. The joint angles were low-pass filtered with a cut-off frequency of 10 Hz to remove high-frequency noise, and differentiated to obtain velocity. The collected EMG of biceps and triceps were band-pass filtered with a cut-off frequency between 20 and 500 Hz to remove motion artifacts and high-frequency noise, rectified, and then low passed at the cut-off frequency 50 Hz for further analysis. As shown in Figure 3, the time for movement onset (t0) and offset (t1) was defined as the time where velocity increased and decreased to 10% of peak velocity (Vpeak) respectively (Atkeson and Hollerbach, 1985). The static angles and EMG during initial posture were then defined as the averaged value from time (t0−1.5) (sec) to time (t0−0.5) (sec) before movement onset. Terminal posture phase started from time (t1 + 4) (sec), and steady state posture angle was obtained in a window of 1 (sec). The onset of muscle firing was calculated as the time at which EMG amplitude exceeded static EMG for three times of standard deviation (3SD) of it, and maintained higher than it for at least 25 ms (Hodges and Bui, 1996; Takatoku and Fujiwara, 2010), the offset of muscle firing was calculated as the time EMG decreased lower than 3 SD below static EMG. According to this threshold criterion, the duration of muscle firing was obtained, and a low passed EMG at the cut off frequency of 6 Hz (Steele, 1994) was used to detect the amplitude of EMG. Off-line Digital Butterworth filters were used in this study with forward and reverse passes to avoid phase shift. The static and dynamic features of movements and EMGs were used to tune parameters of descending commands for corresponding posture and movement.
Specification of Central Descending Commands
To fit the behaviors of the model in Section Corticospinal Virtual Arm (CS-VA) Model to the experimental data in Section Human Reach and Hold Experiment, the parameters of the model are fixed throughout the simulation. Only parameters of central descending commands are adjusted to capture the realistic feature of human movements and postures, as described in the following sub-sections.
Determining Static Command set (αs, γs) for Posture
Posture was realized by adjusting static α and γ commands in this model. Earlier results (Lestienne et al., 1981) indicated that terminal posture could be coded as the ratio of activation levels between antagonistic muscles of a joint, thus αs of Bsh and Tlh in the simulation was set according to the relationship between ratio of static EMG (Bsh/Tlh) (RS) and static elbow angle (θel) of human experiment (Equation 2). The initial and terminal αs of Bsh and Tlh were then decided by corresponding RS at angle θel from the experiment. The value of αs was adjusted in the model to keep a low activation on muscles during postures. The transition of αs between initial and terminal postures was assumed as a ramped changing pattern. To ensure a single joint movement in the simulation, the shoulder was fixed through adding a constant high co-activation level on αs of PC and DP.
γs is related to centrally planned kinematics. In this study, we adopt an experimentally approved criterion for the central planned angle trajectory, the minimal-jerk criterion (Hogan, 1984; Flash and Hogan, 1985). The objective function is given in Equation (3a), and the minimal-jerk trajectory (MJT) is obtained by integration from t0 and t1, which are times of movement onset and end. The outcome is a smooth angle trajectory in Equation (3b), as follows.
The γs input for muscle spindles was calculated from the MJT during posture and movement, according to the quadratic relationship between γs and joint angle obtained in Lan and He (2012). The γs inputs to Bsh and Tlh were determined in Equations (4a) and (4b) in the following:
where the shoulder angle (θsh) was fixed, and the elbow angle trajectory (θel) was the MJT from Equation (3b). The unit of joint angles in Equation (4) is degree.
Determining Dynamic Command Set (αd, γd) for Movement
According Taylor's recording (Taylor et al., 2000, 2004, 2006), γd efferent had a constant baseline of firing rate during both posture and movement, and its firing rate sensitively changed with onset of muscle lengthening, thus presented a leading phase before movement. In the previous analysis of involuntary oscillatory movements (Hao et al., 2013), it was hypothesis that γd represented joint acceleration, and was used to steer activations of antagonistic muscles for joint acceleration and deceleration. This reproduced all features of involuntary oscillatory movements such as Parkinsionian tremor. Thus in the present study, we also adopted this hypothesis, and determined the γd command as the acceleration of the MJT as follows:
In which θel(t) was adopted from Equation (3b). The constant bias was set at 0.5 (normalized) to give a balanced (or mirror) inhibition on antagonistic PNs (see Figure 2). Parameter k was chosen so that the value of γd was within 0 and 1.
Movement was performed through integrating dynamic α and γ commands at the PN network. It has been proposed that, the motor system modulated the pulse amplitude and duration of muscle activations to produce movement with different velocities and ranges (Corcos et al., 1989; Gottlieb et al., 1989), and the pulse strategy has successfully generated scaled movements (Lan, 1997; Lan et al., 2005). Thus, a pair of pulses was utilized on αd to act as agonistic acceleration (AG1) and antagonistic deceleration (ANT) phase respectively. The pulse waveform was given in Equations (6a) and (6b):
where ampAG1 and ampANT were pulse amplitudes of AG1 and ANT (for extension movement, AG1 is on Tlh and ANT is on Bsh respectively), and their values were based on the ratio of peak EMG (ANT/AG1) (RP) obtained in human movement. pwAG1 and pwANT were the pulse widths for AG1 and ANT, which were set as the durations of experimental EMGs according to the 3SD threshold criterion in Section Analysis of Kinematics and EMG Data. tAG1 andtANT were the starting times for AG1 phase and ANT phase, tsim was the terminal time of simulation. Therefore, by adjusting the amplitude and width of the bi-phasic pulses, movements with different ranges and velocities could be achieved by the model.
Determining the Stabilization Pulse
A third pulse with a ramped decreasing form was necessary to stabilize the joint after movement. It was implemented in αs of agonist muscle (second agonistic muscle burst, AG2), since tri-phasic patterned EMG is found a common feature in fast reaching movement (Ghez and Martin, 1982; Flanders et al., 1994; Berardelli et al., 1996), and AG2 has been widely accepted as a phase to help stabilize the elbow at terminal posture after movement (Hannaford and Stark, 1985; Takatoku and Fujiwara, 2010). The amplitude and duration of AG2 were adjusted in the simulation to stabilize the joint after movement.
Simulation with Intact Feedforward and Feedback Controls
Normal condition with intact feedforward and feedback control was simulated to match typical trials in human experiment. For each subject, different postures were simulated and compared to the experimental relationship between RS and steady state angle. Stiffness control by varying co-activation level of antagonist muscles at different postures was demonstrated. To determine joint stiffness, a force perturbation was added directly on single joint muscle BS, and the stiffness was calculated as the ratio of changes in elbow torque with respect to joint angle (He et al., 2013). For each subject, the extension movement from 60° to 30° (Movement 1) was specifically simulated and matched to experimental data. More simulations were performed to fit experimental movements of larger range and in opposite direction for Subject 1, Simulation of 11 s was usually performed. After the system had converged at a steady initial state, a random, signal dependent noise (Jones et al., 2002) was added on Um, the elbow was driven at time 6 (sec). Reflex gains used in this study were set within the range that kept the system stable (Table 1) (He et al., 2013). All α and γ commands used in simulation were normalized to values between 0 and 1. Nominal parameter values of spinal circuits were adopted from Hao et al. (2013), and were listed in Table 1.
Simulation with Abnormal Feedforward and Feedback Controls
In this study, we examined the effects of altering model structure on movement and posture by deafferenting and deefferenting the model. Deafferented conditions were modeled by assigning the related gains of Ia, and Ib to zero. Deefferented conditions of the model were obtained by setting one of the descending commands to zero at a time, while keeping others intact. The effect of abnormal ratio (Bsh/Tlh) of αs on terminal postures was also examined by assigning static ratio out of normal ranges of experiments. Note that γd was always maintained a bias of 0.5 in the analysis of deefferentation study. The abnormal conditions were applied to the model 0.1 (sec) before movement initiation in simulations.
The difference between experimental and simulated movements was used to quantify the effects of abnormal control of movement and posture. Angular errors for movements and terminal steady state posture were evaluated as in Equations (7a) and (7b):
where ErrorP and ErrorM gave the averaged differences in terminal steady state posture and movement between simulated and experiment respectively (Figure 3), W and V were the number of resampled data points during terminal steady state posture and movement phases.
Results
Features in Kinematics and EMGs of Human Posture and Movement
All subjects presented a common pattern in movement and posture at the elbow joint. The joint angle was stabilized after a slight overshoot, and the velocity displayed a bell-shaped profile. Tri-phasic EMG was observed in fast movements, with AG1 firing at accelerating phase, ANT at deceleration phase. The AG2 of agonist muscle appeared during movement offset and lasted until the joint was stabilized at terminal posture. When the elbow was stabilized after movement, the EMG of Tlh and Bsh were maintained at different steady state levels, which varied with elbow angle.
An example of the relations between EMG ratio and joint angle during posture and movement was presented in Figure 4. The EMG of Tlh and Bsh were found to vary reciprocally with the elbow angle from 0° to 110°. The ratio of static EMG (Bsh/Tlh) (RS) was well fitted against the elbow angle in a linear manner (Figure 4A, P < 0.01). During movement although the peak velocity increased with movement range and speed, The relationship between peak velocity and ratio of peak EMG (ANT/AG1) (RP) was almost flat (P > 0.05), as shown in Figure 4B, and the ratio of peak EMG had an average value of 0.67 (±0.29).
Figure 4. Experimental kinematic and EMG during posture and movement of Subject 1. (A) Relationship between postural EMG ratio of Biceps and Triceps and elbow angle, and a linear regression is fitted with P < 0.01, R2 = 0.8715. (B) Relationship between the ratio of peak EMG (ANT/AG1, Bsh/Tlh for extension, Tlh/Bsh for flexion) and peak velocity for fast reaching. The linear regression is rejected (intercept P > 0.05).
These features were present in all subjects. A summary of experimental relation between the ratio of static EMGs and elbow angle for all 6 subjects was given in Table 2. This indicates that the central module for posture control purposefully varies the static levels of antagonistic muscle activation in such particular linear manner to compensate for the change in moment arms of muscles at different joint angles. This experimental linear relationship is sufficient to guide the specification of αs descending commands in the model for posture maintenance.
Table 2. Linear relationship between ratio of static EMG (RS) and elbow angle (θel) (Experiment (Exp.) vs. Simulation (Sim.)).
Simulated Posture Control
Posture control was then simulated by tuning αs inputs to Bsh and Tlh for each subject, and by setting γs according to the static angle θel. Using the linear relation to guide the specification of αs commands for Bsh and Tlh muscles in the model, a comparable relationship was obtained in simulated postures for all subjects. The results in Table 2 showed that the simulated linear equations for each subject matched to the experimental relationship closely (P < 0.05). Overall an average linear relationship existed in experiment and simulated groups for all six subjects (P < 0.05), and they had a similar slope and intercept. This indicates that posture control can be achieved by tuning αs, as well as setting γs inputs to relevant antagonistic muscles.
Thus, the posture module of motor system can control joint stiffness while maintaining the same joint posture. Figure 5 demonstrated such a strategy by increasing the co-activation levels of antagonistic muscles for Bsh and Tlh, while keeping the γs commands at the values corresponding to elbow angles of 30°, 50°, 65° respectively. Figure 5A showed that the simulated posture angles against the experimental linear relation. Since multiple activation levels of muscles can achieve the same elbow posture, joint stiffness can be controlled to desired levels for different motor tasks. We calculated joint stiffness by applying force perturbations on BS and measured the ratio between change of joint torque and angle. Figure 5B showed that joint stiffness increased with the muscle activation, while the same posture was maintained. The results indicate that the central motor system can control joint stiffness and posture independently by tuning the levels of αs commands with programmed γs control for postures.
Figure 5. Simulated posture control. (A) Muscle activation ratio of simulated postures in comparison with experiment's linear regression line. Different postures were maintained while keeping the ratio of muscle activation of Bsh and Tlh within the range of experimental value. (B) Postures maintained with increasing joint stiffness. The elbow angles were maintained at 30°, 50°, and 65° by increasing muscles activation on Bsh, but keeping the ratio of Bsh and Tlh within experimental range, independently controlling joint stiffness was realized while maintaining the same posture.
Simulated Movement Control
The extension movement from 60° to 30° (Movement 1) was simulated for six subjects individually. An extension movement from 90° to 45° (Movement 2), a flexion movement from 30° to 60° (Movement 3) were simulated for Subject 1. Figures 6A–C depict the descending commands tuned for Movement 1, 2, and 3 respectively. During commands assignment, γs and γd commands were calculated from the MJT of target trajectory (Equations 4 and 5). And the initial and terminal postures were maintained by setting the ratio of αs of Bsh and Tlh as experimental values. Thus, movements between postures were simulated with tuning of the pulse height and width of the αd command. The amplitudes and pulse widths were adopted from experimental EMG (Equation 6). The third pulse in the αs of Tlh acted to help stabilize the elbow after fast movement. The posture commands were also specified for the other four muscles in the model to keep the simulation running properly.
Figure 6. Descending commands of Simulated Movement 1, 2, and 3 for Subject 1. (A) Descending commands of Movement 1. The joint was held at initial posture of 63.36° and terminal posture of 29.67° by posture commands, the amplitude and pulse width for AG1 (Tlh) and ANT (Bsh) were 0.60, 0.27 (sec) and 0.25, 0.18 (sec) respectively. (B,C) Present the descending commands for movements 2 and 3 respectively. The amplitudes and pulse widths of αd for Movement 2 were ampAG1(0.69), pwAG1(0.28(sec)), ampANT(0.30), pwANT(0.20(sec)), in Movement 3, they were ampAG1(0.34), pwAG1(0.26(sec)), ampANT(0.26), pwANT(0.19 (sec)).
The simulated angle, velocity and muscle activation of Movement 1, 2, and 3 were compared with the experimental movements (Figures 7A–C). The initial and terminal angles could finely match the experimental postures, and the three movements displayed the signature bell-shaped velocity profile, and tri-phasic firing pattern of Bsh and Tlh. Results showed that the elbow could be maintained closely at terminal postures, and the errors during movements were relatively small compared to the range of movements.
Figure 7. Simulated Movement 1, 2, and 3, in comparison with experimental data, and they are presented in (A–C) respectively. The angle and velocity profile were low-passed filtered with the cut off frequency 10 Hz, and the simulated and experimental muscle activation were low-passed filtered with the cut off frequency 50 Hz (Bsh with inverse value). The ErrorP of Movement 1, 2, and 3 were 0.19 (deg), −0.04 (deg) and −0.25 (deg), and their ErrorM were −3.99 (deg), −5.66 (deg), and −1.44 (deg).
Furthermore, all simulated movements matched well to those of experimental trajectories and EMGs in all subjects by tuning the αd pulse command. The errors of trajectory fitting for Movement 1 were calculated according to Equation (7) and listed in Table 3. It is demonstrated that the movement module of the central motor system could control movements of a range of angles and velocities by tuning the αd pulse command while coordinating posture commands.
Simulation with Abnormal Afferent and Efferent Controls
Movements under abnormal afferent and efferent controls for Movement 1 were simulated to assess the individual contribution of each descending command to movement and posture. Figure 8 illustrated the various effects for Subject 1. Figure 8A showed that deafferentation after Ia(−) was removed overshot the target of movement with a smaller steady state angle; while removing Ia(+) and Ia-PN undershot the target of movement, note that removing Ia(+) resulted in a larger steady state angle, while Ia-PN had no effect on final posture. When all Ia afferents to αMN were removed in deafferented state, the joint undershot the target of movement, and stabilized gradually to a slightly larger steady state posture.
Figure 8. Simulated movement with abnormal afferent and efferent controls. The simulations were based on Movement 1 of Subject 1. (A) Simulated elbow angle trajectory and velocity under conditions of abnormal afferent feedback. Each feedback gain of Ia afferent was set to zero prior to movement initiation in separated trials. The steady state angle of Ia (+), Ia (−), Ia-PN removed and deafferented conditions were 33.20°, 27.53°, 29.71°, and 31.03°. (B) Simulated elbow trajectory and velocity with abnormal static and dynamic efferent commands. The joint terminated at 30.58° under doubled αs. After γs was removed, the terminal angle was 32.45 (±0.42)°, and the standard variation was 156% of that under normal condition.
The results of abnormal efferent commands were shown in Figure 8B. When the ratio of αs was doubled, the movement and steady state posture were slightly affected compared to those in normal condition. When γs was removed, the joint undershot the target of movement similar to that of deafferentation. After αd was removed or γd kept a constant inhibition on PNs, the elbow couldn't make the fast movement as in normal condition. The slow approach to the target posture was brought about by the action of spinal reflexes. Thus, feedforward control of descending commands (αd and γd) is absolutely essential for fast movements.
This general pattern of behavior in all subjects is summarized in the phase diagram of errors in Figure 9. As illustrated in Figure 9A, removing Ia(+), Ia(−) and all Ia afferents generally impacted the terminal postures, these results were consistent with the behaviors observed in humans and animals (Polit and Bizzi, 1979; Gentilucci et al., 1994; Gordon et al., 1995). But removing Ia-PN seemed to reduce errors in both movement and posture. Static and dynamic descending commands showed distinct impacts on posture and movement (Figure 9B). The deviated ratios of αs for Bsh and Tlh from normal values yielded a wider departure from the target posture, which was different for individual subjects. The absence of γs resulted in an error primarily in posture. The importance of dynamic commands for movement was clearly seen from the large ErrorM in Figure 9B. This is consistent with the observation in cats after the PNs innervating the upper extremity muscles were removed (Alstermark et al., 1981; Alstermark and Isa, 2012).
Figure 9. Phase diagram of the effect of abnormal afferent feedback and efferent commands on posture and movement for six subjects. All simulations were based on the Movement 1, shown in Table 3. (A) The errors of terminal posture and movement comparing with experimental trials, without Ia(+), Ia(−), Ia-PN, and all of them, respectively. ErrorP and ErrorM are the mean difference of θel between simulated and experimental trial, during terminal posture and movement, respectively. The dashed lines represented no posture or movement error. (B) The errors of terminal posture and movement comparing with experimental trials, with abnormal αs, blocked γs, blocked αd and constant γd, respectively. The average of ErrorP and ErrorM was 0.21(deg) and 14.75 (deg) for condition of αd removed, and 0.20 (deg) and 14.92 (deg) for condition of γd kept constant.
Discussion
An important issue in motor control has been to differentiate the role of the γ motor system from that of the α motor system (Stein, 1974; Houk and Rymer, 1981; Bizzi et al., 1991). Although increasing evidence revealed α-γ co-activation during movement and posture control (Vallbo, 1971; Taylor et al., 2006; Prochazka and Ellaway, 2012), the dominant effects of the α motor system tended to diminish the function significance of the γ motor system. Using a computational virtual arm (VA) model (He et al., 2013), Lan and He (2012) re-interpreted a set of experiment data from Stein et al. (2004), Cordo et al. (2002), Taylor et al. (2006), and suggested that γs may encode centrally planned information of joint angle and reinforce the planned joint angle though regulating spindle sensitivity. By coupling the VA model with a corticospinal (CS) network of PN in analyzing involuntary oscillatory movements (Hao et al., 2013), we further hypothesized that γd represents centrally planned joint acceleration. In this paper, the combined CS-VA model was used to implicate the necessity of α-γ coordination in movement and posture control. Only the set of α-γ descending commands were adjusted to fit human movement data. Close match of model behaviors to those observed in human experiments as demonstrated in the results (Figure 7, Table 3) established the validity of the model, therefore, providing a neurophysiologically realistic, multi-scale computational model to evaluate the contribution of various components of descending control in sensorimotor functions. In particular, this model will also be valuable to understand sensorimotor dysfunctions (Hao et al., 2013) and to design novel rehabilitation strategies for motor relearning (Zhuang et al., 2015).
A distinct feature of this computational model is the division of sensorimotor control into movement module and posture module by the dynamic and static descending command sets (Bizzi et al., 2008; Diedrichsen and Classen, 2012). Results of this study indicate that it is possible to coordinate the two sets of descending α-γ commands to achieve accurate control of movement dynamics and stable maintenance of final posture. Methods adopted in this study to determine descending commands are novel, in that these descending commands are specified according to proven rules of central planning for kinematics (Hogan, 1984; Flash and Hogan, 1985) and EMG data collected in human subjects performing reach-and-hold tasks. Assumptions are also adopted from previous analytical studies regarding γs encoding of joint angle (Lan and He, 2012), and γd representation of joint acceleration of planned movement trajectory (Taylor et al., 2006; Hao et al., 2013). Results in Figure 5 and Table 2 indicate that αs and αd commands can be tuned based on EMG signals of human data at steady state and during movement. Figure 5 shows that tuning αs based on the linear ratio of antagonistic muscles in Figure 4 is necessary to achieve independent control of joint posture and stiffness, which is an important aspect of regulation of motor functions by the central motor system (Mussa-Ivaldi et al., 1985). For movement control, the αd command is integrated with γd command at the PN to distribute properly the activation to flexor or extensor acting at the joint, and its pulse amplitude and width can be tuned according to the speed of movement and duration of EMG bursts. We illustrated that adjusting these descending commands can fit reach-and-hold movements for a range of amplitude and direction in different subjects. The ability of the model to fit experiment movements suggests that the computational model captures the neural mechanism of corticospinal computation, as well as the modular nature of organization and coordination of descending α-γ commands by the central motor control system (Ghez et al., 2007; Scheidt and Ghez, 2007; Scheidt et al., 2011; Poston et al., 2013).
The model predicted the contribution of descending α-γ commands to movement and posture by deefferenting the CS-VA model in simulation. With abnormal αs command out of the range of experimental ratio (Figures 8B, 9B), the terminal angle deviated from its targeted position. The large variation shown at the terminal position suggests that accurate γs command is essential for posture maintenance. The tardy movements under abnormal αd and γd (Figure 9B) indicated that fast reaching movement must be performed with proper coordination of αd and γd commands. This is consistent with the observation that cats were unable to carry out skilled reaching movement without PN (Alstermark et al., 1981; Alstermark and Isa, 2012).
Proprioceptive afferents from muscle spindle are important for motor learning (Jeannerod, 1988; Schmidt and Lee, 2011), but are not found indispensible for control of movements, since deafferentation in human patients and animals did not entirely disable their movement execution. Early study in deafferented patients indicated obvious motor dysfunction only with larger errors in terminal position acquisition (Polit and Bizzi, 1979; Gentilucci et al., 1994; Gordon et al., 1995) and lower joint stiffness during posture maintenance (Bizzi et al., 1984), but changes in movement kinematic and muscle firing pattern were not obvious (Taub et al., 1966, 1975; Vaughan et al., 1970; Rothwell et al., 1982; Bizzi et al., 1984; Gordon et al., 1995). This implied that proprioceptive afferent played an important role in posture maintenance and contributed to fine control of movement.
This functional role of proprioceptive afferents is reiterated by deafferenting the CS-VA model (Figures 8A, 9A). It was shown that, when Ia(+) and Ia-PN was removed respectively, the movement slowed down, and the terminal posture shifted, and after Ia(−) was removed, the movement speeded up, and the terminal posture targeted at a lower angle. This confirmed the positive feedback of Ia(+) and Ia-PN, and inhibition role of Ia(−). However, despite the difference of peak velocity, these movements presented similar bell-shaped velocity profile, only settled at different terminal angles. Deafferented simulation showed slowed movement and increased errors in terminal angle. This was similar to the behaviors observed in deafferented primates (Polit and Bizzi, 1979), which showed that the deafferented primate after intensive training was able to control pointing movement with normal like kinematics and muscle activation, but inaccurate positions.
Conclusion
A corticospinal computational model based on the modular organization of movement and posture was validated in this study by fitting the model to experimental human data. Analysis of simulated movement and posture with intact and altered model structures demonstrated that it is necessary to coordinate the set of α-γ descending commands in order to achieve effective control of accurate movement dynamics and stable postures. Results suggest that the central commands of posture module are mediated via a mono-synaptic corticospinal pathway, while those of movement module are transmitted to spinal motoneruons through a multi-synaptic corticospinal pathway involving the propriospinal neurons (PN). The PN network plays the pivotal role to integrate the αd and γd commands for movement generation. The model is able to capture many essential aspects of motor behaviors, such as independent regulation of joint angle and stiffness and the signature temporal pattern of EMGs, by simply tuning the αs and αd commands. This study suggests a plausible neural computational mechanism for the central motor system to control movement and posture. The model will be useful as a complementary tool to understand neural control of movements, as well as a valuable platform to aid design of novel rehabilitation strategy for motor disabilities.
Author Contributions
SL performed model simulation, analyzed experimental data, prepared figures and tables and drafted the manuscript; XH and MH contributed to set up the computational model; JM and CZ contributed in the analysis of human experiment data. CZ also carried out part of simulation work. CN offered intellectual suggestions and edited the manuscript. NL conceived the computational approach, designed the human subject experiment, proposed analytical method and edited the final version of the manuscript.
Conflict of Interest Statement
The Chief Editor Prof Wu declares that, despite having sharing affiliation with the Reviewer Dr. Wang, the review process was handled objectively and no conflict of interest exists. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work is supported in part by grants from the Natural Science Foundation of China (No. BC0820041 and No. BC0820010), and the National Basic Research Program of Project 973 of the Ministry of Science and Technology of China (2011CB013304).
Abbreviations
α, γ, Alpha, gamma motor nerves and neurons; PN, Propriospinal neuron; EMG, Electromyography; GTO, Golgi Ten Organ; CS-VA, Corticospinal virtual arm; VA, Virtual arm; αs, γs, Static alpha, gamma command; αd, γd, Dynamic alpha, gamma command; αMN, γMN, Alpha, gamma motoneuron; Um, Muscle input in the CS-VA model; Us, Spindle input in the CS-VA model; Fm, Muscle force; Lce, Muscle fascicle length; Lmt, Musculo-tendon length; Ia, Primaryafferent from spindle; II, Secondary afferent from spindle; Ib, Afferent from GTO; PC, Pectoralis Major Clavicle; DP, Deltoid Posterior; BS, Brachialis; Tlt, Triceps Lateral; Bsh, Biceps Short Head; Tlh, Triceps Long Head; θel, Elbow angle; θsh, Shoulder angle; γ-PN, Dynamic gamma inhibition on PN; Ia(−), Ia reciprocal inhibition on αMN; Ia-PN, Ia afferent gain on PN; Ia(+), Stretch reflex; Ib(−), Ib inhibition from GTO; RC, Renshaw Cell; ErrorP, ErrorM, mean difference between simulation and experiment of terminal posture, movement; Vpeak, Peak velocity; 3SD, Three times of standard deviation of static EMG; AG1, First agonistic muscle burst phase; AG2, Second agonistic muscle burst phase; ANT, Antagonistic muscle burst phase; RS, Ratio of static EMG (Bsh/Tlh); RP, Ratio of peak EMG (ANT/AG1); ampAG1 ampANT, Amplitude of AG1 and ANT phase on αd; pwAG1,pwANT, Pulse width of AG1 and ANT phase on αd; MJT, Minimum jerk trajectory.
References
Alstermark, B., and Isa, T. (2012). Circuits for skilled reaching and grasping. Annu. Rev. Neurosci. 35, 559–578. doi: 10.1146/annurev-neuro-062111-150527
Alstermark, B., Isa, T., Pettersson, L. G., and Sasaki, S. (2007). The C3–C4 propriospinal system in the cat and monkey: a spinal pre-motoneuronal centre for voluntary motor control. Acta Physiol. 189, 123–140. doi: 10.1111/j.1748-1716.2006.01655.x
Alstermark, B., Lundberg, A., Norrsell, U., and Sybirska, E. (1981). Integration in descending motor pathways controlling the forelimb in the cat - 9. Differential behavioural defects after spinal cord lesions interrupting defined pathways from higher centres to motoneurones. Exp. Brain Res. 42, 299–318. doi: 10.1007/BF00237496
Alstermark, B., Lundberg, A., and Sasaki, S. (1984a). Integration in descending motor pathways controlling the forelimb in the cat - 10. Inhibitory pathways to forelimb motoneurons via C3–C4 propriospinal neurones. Exp. Brain Res. 56, 279–292. doi: 10.1007/BF00236284
Alstermark, B., Lundberg, A., and Sasaki, S. (1984b). Integration in descending motor pathways controlling the forelimb in the cat - 11. Inhibitory pathways from higher motor centres and forelimb afferents to C3–C4 propriospinal neurones. Exp. Brain Res. 56, 293–307. doi: 10.1007/BF00236285
Alstermark, B., Lundberg, A., and Sasaki, S. (1984c). Integration in descending motor pathways controlling the forelimb in the cat. 12. Interneurones which may mediate descending feed-forward inhibition and feed-back inhibition from the forelimb to C3–C4 propriospinal neurones. Exp. Brain Res. 56, 308–322. doi: 10.1007/BF00236286
Atkeson, C. G., and Hollerbach, J. M. (1985). Kinematic features of unrestrained vertical arm movements. J. Neurosci. 5, 2318–2330.
Baldissera, F., Hultborn, H., and Illert, M. (1981). “Integration in spinal neuronal systems,” in Handbook of Physiology, Section I, The Nervous System, Vol II, Motor Control, Part I, ed V. B. Brooks (Bethesda, MD: American Physiological Society), 509–595.
Berardelli, A., Hallett, M., Rothwell, J. C., Agostino, R., Manfredi, M., Thompson, P. D., et al. (1996). Single-joint rapid arm movements in normal subjects and in patients with motor disorders. Brain 119, 661–674. doi: 10.1093/brain/119.2.661
Bizzi, E., Accornero, N., Chapple, W., and Hogan, N. (1984). Posture control and trajectory formation during arm movement. J. Neurosci. 4, 2738–2744.
Bizzi, E., Cheung, V. C. K., d'Avella, A., Saltiel, P., and Tresch, M. (2008). Combining modules for movement. Brain Res. Rev. 57, 125–133. doi: 10.1016/j.brainresrev.2007.08.004
Bizzi, E., Mussa-Ivaldi, F. A., and Giszter, S. (1991). Computations underlying the execution of movement: a biological perspective. Science 253, 287–291. doi: 10.1126/science.1857964
Boyd, I. A. (1980). The isolated mammalian muscle spindle. Trends Neurosci. 3, 258–265. doi: 10.1016/0166-2236(80)90096-X
Cheng, E. J., Brown, I. E., and Loeb, G. E. (2000). Virtual muscle: a computational approach to understanding the effects of muscle properties on motor control. J. Neurosci. Methods 101, 117–130. doi: 10.1016/S0165-0270(00)00258-2
Corcos, D. M., Gottlieb, G. L., and Agarwal, G. C. (1989). Organizing principles for single-joint movements. II. A speed-sensitive strategy. J. Neurophysiol. 62, 358–368.
Cordo, P. J., Flores-Vieira, C., Verschueren, S. M., Inglis, J. T., and Gurfinkel, V. (2002). Position sensitivity of human muscle spindles: single afferent and population representations. J. Neurophysiol. 87, 1186–1195. doi: 10.1152/jn.00393.2001
Diedrichsen, J., and Classen, J. (2012). Stimulating news about modular motor control. Neuron 76, 1043–1045. doi: 10.1016/j.neuron.2012.12.001
Dimitriou, M., and Edin, B. B. (2008). Discharges in human muscle spindle afferents during a key-pressing task. J. Physiol. 586, 5455–5470. doi: 10.1113/jphysiol.2008.160036
Eccles, J. C., Eccles, R. M., and Lundberg, A. (1957). The convergence of monosynaptic excitatory afferents on to many different species of alpha motoneurones. J. Physiol. 137, 22–50. doi: 10.1113/jphysiol.1957.sp005794
Eccles, J. C., Eccles, R. M., and Magni, F. (1960). Monosynaptic excitatory action on motoneurones regenerated to antagonistic muscles. J. Physiol. 154, 68–88. doi: 10.1113/jphysiol.1960.sp006565
Flanders, M., Pellegrini, J. J., and Soechting, J. F. (1994). Spatial/temporal characteristics of a motor pattern for reaching. J. Neurophysiol. 71, 811–813.
Flash, T., and Hogan, N. (1985). The coordination of arm movements: an experimentally confirmed mathematical model. J. Neurosci. 5, 1688–1703.
Fuglevand, A. J., Winter, D. A., and Patla, A. E. (1993). Models of recruitment and rate coding organization in motor-unit pools. J. Neurophysiol. 70, 2470–2488.
Gentilucci, M., Toni, I., Chieffi, S., and Pavesi, G. (1994). The role of proprioception in the control of prehension movements: a kinematic study in a peripherally deafferented patient and in normal subjects. Exp. Brain Res. 99, 483–500. doi: 10.1007/BF00228985
Ghez, C., and Martin, J. H. (1982). The control of rapid limb movement in the cat. III. Agonist-antagonist coupling. Exp. Brain Res. 45, 115–125. doi: 10.1007/bf00235770
Ghez, C., Scheidt, R., and Heijink, H. (2007). Different learned coordinate frames for planning trajectories and final positions in reaching. J. Neurophysiol. 98, 3614–3626. doi: 10.1152/jn.00652.2007
Gordon, J., Ghilardi, M. F., and Ghez, C. (1995). Impairments of reaching movements in patients without proprioception. I. Spatial errors. J. Neurophysiol. 73, 347–360.
Gottlieb, G. L., Corcos, D. M., and Agarwal, G. C. (1989). Organizing principles for single-joint movements. I. A speed-insensitive strategy. J. Neurophysiol. 62, 342–357.
Gracies, J. M., Meunier, S., and Pierrot-Deseilligny, E. (1994). Evidence for corticospinal excitation of presumed propriospinal neurones in man. J. Physiol. 475, 509–518. doi: 10.1113/jphysiol.1994.sp020089
Granit, R. (1975). The functional role of the muscle spindles: facts and hypotheses. Brain 98, 531–556. doi: 10.1093/brain/98.4.531
Hannaford, B., and Stark, L. (1985). Roles of the elements of the triphasic control signal. Exp. Neurol. 90, 619–634. doi: 10.1016/0014-4886(85)90160-8
Hao, M., He, X., Xiao, Q., Alstermark, B., and Lan, N. (2013). Corticomuscular transmission of tremor signals by propriospinal neurons in Parkinson's disease. PLoS ONE 8:e79829. doi: 10.1371/journal.pone.0079829
He, X., Du, Y. F., and Lan, N. (2013). Evaluation of feedforward and feedback contributions to hand stiffness and variability in multijoint Arm control. Neural Syst. Rehabil. Eng. IEEE Trans. 21, 634–647. doi: 10.1109/TNSRE.2012.2234479
Hodges, P. W., and Bui, B. H. (1996). A comparison of computer-based methods for the determination of onset of muscle contraction using electromyography. Electroencephalogr. Clin. Neurophysiol. 101, 511–519. doi: 10.1016/S0921-884X(96)95190-5
Hogan, N. (1984). An organizing principle for a class of voluntary movements. J. Neurosci. 4, 2745–2754.
Houk, J. C., and Rymer, W. Z. (1981). “Neural control of muscle length and tension,” in Handbook of Physiology. Sect. I. The Nervous System. Vol. II, Motor Control, Part 1, ed V. B. Brooks (Bethesda, MD: American Physiological Society), 257–323.
Hulliger, M. (1984). The mammalian muscle spindle and its central control. Rev. Physiol. Biochem. Pharmacol. 101, 1–110. doi: 10.1007/BFb0027694
Isa, T., Ohki, Y., Alstermark, B., Pettersson, L. G., and Sasaki, S. (2007). Direct and indirect cortico-motoneuronal pathways and control of hand/arm movements. Physiology 22, 145–152. doi: 10.1152/physiol.00045.2006
Jeannerod, M. (1988). The Neural and Behavioural Organization of Goal-directed Movements. New York, NY: Clarendon Press/Oxford University Press.
Jones, K. E., Hamilton, A. F. D. C., and Wolpert, D. M. (2002). Sources of signal-dependent noise during isometric force production. J. Neurophysiol. 88, 1533–1544. doi: 10.1152/jn.00985.2001
Kurtzer, I., Herter, T. M., and Scott, S. H. (2005). Random change in cortical load representation suggests distinct control of posture and movement. Nat. Neurosci. 8, 498–504. doi: 10.1038/nn1420
Lan, N. (1997). Analysis of an optimal control model of multi-joint arm movements. Biol. Cybern. 76, 107–117. doi: 10.1007/s004220050325
Lan, N., and He, X. (2012). Fusimotor control of spindle sensitivity regulates central and peripheral coding of joint angles. Front. Comput. Neurosci. 6:66. doi: 10.3389/fncom.2012.00066
Lan, N., Li, Y., Sun, Y., and Yang, F. S. (2005). Reflex regulation of antagonist muscles for control of joint equilibrium position. Neural Syst. Rehabil. Eng. IEEE Trans. 13, 60–71. doi: 10.1109/TNSRE.2004.841882
Lawrence, D. G., and Kuypers, H. G. J. M. (1968). The functional organization of the motor system in the monkey: I. The effects of bilateral pyramidal lesions. Brain 91, 1–14. doi: 10.1093/brain/91.1.1
Lemon, R. N. (1999). Neural control of dexterity: what has been achieved? Exp. Brain Res. 128, 6–12. doi: 10.1007/s002210050811
Lemon, R. N., Kirkwood, P. A., Maier, M. A., Nakajima, K., and Nathan, P. (2004). Direct and indirect pathways for corticospinal control of upper limb motoneurons in the primate. Prog. Brain Res. 143, 263–279. doi: 10.1016/S0079-6123(03)43026-4
Lestienne, F., Polit, A., and Bizzi, E. (1981). Functional organization of the motor process underlying the transition from movement to posture. Brain Res. 230, 121–131. doi: 10.1016/0006-8993(81)90396-6
Li, S., He, X., and Lan, N. (2014). Modular control of movement and posture by the corticospinal alpha-gamma motor systems. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2014, 4079–4082. doi: 10.1109/EMBC.2014.6944520
Malmgren, K., and Pierrot-Deseilligny, E. (1988). Evidence for non-monosynaptic Ia excitation of human wrist flexor motoneurones, possibly via propriospinal neurones. J. Physiol. 405, 747–764. doi: 10.1113/jphysiol.1988.sp017359
Martin, V., Scholz, J. P., and Schöner, G. (2009). Redundancy, self-motion, and motor control. Neural Comput. 21, 1371–1414. doi: 10.1162/neco.2008.01-08-698
Mileusnic, M. P., Brown, I. E., Lan, N., and Loeb, G. E. (2006). Mathematical models of proprioceptors. I. Control and transduction in the muscle spindle. J. Neurophysiol. 96, 1772–1788. doi: 10.1152/jn.00868.2005
Mileusnic, M. P., and Loeb, G. E. (2006). Mathematical models of proprioceptors. II. Structure and function of the Golgi tendon organ. J. Neurophysiol. 96, 1789–1802. doi: 10.1152/jn.00869.2005
Mussa-Ivaldi, F. A., Hogan, N., and Bizzi, E. (1985). Neural, mechanical, and geometric factors subserving arm posture in humans. J. Neurosci. 5, 2732–2743.
Pierrot-Deseilligny, E., and Burke, D. (2005). The Circuitry of the Human Spinal Cord. Its Role in Motor Control and Movement Disorders. New York, NY: Cambridge University Press, 642.
Polit, A., and Bizzi, E. (1979). Characteristics of motor programs underlying arm movements in monkeys. J. Neurophysiol. 42, 183–194.
Poston, B., Van Gemmert, A. W. A., Sharma, S., Chakrabarti, S., Zavaremi, S. H., and Stelmach, G. (2013). Movement trajectory smoothness is not associated with the endpoint accuracy of rapid multi-joint arm movements in young and older adults. Acta Psychol. 143, 157–167. doi: 10.1016/j.actpsy.2013.02.011
Prochazka, A. (1999). Quantifying proprioception. Prog. Brain Res. 123, 133–142. doi: 10.1016/S0079-6123(08)62850-2
Prochazka, A., and Ellaway, P. (2012). Sensory systems in the control of movement. Compr. Physiol. 2, 2615–2627. doi: 10.1002/cphy.c100086
Prochazka, A., Hulliger, M., Zangger, P., and Appenteng, K. (1985). ‘Fusimotor set’: new evidence for α-independent control of γ-motoneurones during movement in the awake cat. Brain Res. 339, 136–140. doi: 10.1016/0006-8993(85)90632-8
Quallo, M. M., Kraskov, A., and Lemon, R. N. (2012). The activity of primary motor cortex corticospinal neurons during tool use by macaque monkeys. J. Neurosci. 32, 17351–17364. doi: 10.1523/jneurosci.1009-12.2012
Rothwell, J. C., Traub, M. M., Day, B. L., Obeso, J. A., Thomas, P. K., and Marsden, C. D. (1982). Manual motor performance in a deafferented man. Brain 105, 515–542. doi: 10.1093/brain/105.3.515
Scheidt, R. A., and Ghez, C. (2007). Separate adaptive mechanisms for controlling trajectory and final position in reaching. J. Neurophysiol. 98, 3600–3613. doi: 10.1152/jn.00121.2007
Scheidt, R. A., Ghez, C., and Asnani, S. (2011). Patterns of hypermetria and terminal cocontraction during point-to-point movements demonstrate independent action of trajectory and postural controllers. J. Neurophysiol. 106, 2368–2382. doi: 10.1152/jn.00763.2010
Schmidt, R. A., and Lee, T. D. (2011). Motor Control and Learning: A Behavioral Emphasis. Champaign, IL: Human Kinetics.
Song, D., Lan, N., Loeb, G. E., and Gordon, J. (2008a). Model-based sensorimotor integration for multi-joint control: development of a virtual arm model. Ann. Biomed. Eng. 36, 1033–1048. doi: 10.1007/s10439-008-9461-8
Song, D., Raphael, G., Lan, N., and Loeb, G. E. (2008b). Computationally efficient models of neuromuscular recruitment and mechanics. J. Neural Eng. 5, 175–184. doi: 10.1088/1741-2560/5/2/008
Steele, C. (1994). Anticipatory Postural Reactions in Clumsy Children: Changes after Physiotherapy Treatment. Master of Physiotherapy Thesis, University of Queensland.
Stein, R. B., Weber, D. J., Aoyagi, Y., Prochazka, A., Wagenaar, J. B. M., Shoham, S., et al. (2004). Coding of position by simultaneously recorded sensory neurones in the cat dorsal root ganglion. J. Physiol. 560, 883–896. doi: 10.1113/jphysiol.2004.068668
Takatoku, N., and Fujiwara, M. (2010). Muscle activity patterns during quick increase of movement amplitude in rapid elbow extensions. J. Electromyogr. Kinesiol. 20, 290–297. doi: 10.1016/j.jelekin.2009.03.007
Taub, E., Ellman, S. J., and Berman, A. J. (1966). Deafferentation in monkeys: effect on conditioned grasp response. Science 151, 593–594. doi: 10.1126/science.151.3710.593
Taub, E., Goldberg, I. A., and Taub, P. (1975). Deafferentation in monkeys: pointing at a target without visual feedback. Exp. Neurol. 46, 178–186. doi: 10.1016/0014-4886(75)90040-0
Taylor, A., Durbaba, R., and Ellaway, P. H. (2004). Direct and indirect assessment of γ-motor firing patterns. Can. J. Physiol. Pharmacol. 82, 793–802. doi: 10.1139/y04-053
Taylor, A., Durbaba, R., Ellaway, P. H., and Rawlinson, S. (2006). Static and dynamic γ-motor output to ankle flexor muscles during locomotion in the decerebrate cat. J. Physiol. 571, 711–723. doi: 10.1113/jphysiol.2005.101634
Taylor, A., Ellaway, P. H., Durbaba, R., and Rawlinson, S. (2000). Distinctive patterns of static and dynamic gamma motor activity during locomotion in the decerebrate cat. J. Physiol. 529, 825–836. doi: 10.1111/j.1469-7793.2000.00825.x
Vallbo, A. B. (1971). Muscle spindle response at the onset of isometric voluntary contractions in man. Time difference between fusimotor and skeletomotor effects. J. Physiol. 218, 405–431. doi: 10.1113/jphysiol.1971.sp009625
Vaughan, H. G. Jr., Gross, E. G., and Bossom, J. (1970). Cortical motor potential in monkeys before and after upper limb deafferentation. Exp. Neurol. 26, 253–262. doi: 10.1016/0014-4886(70)90123-8
Windhorst, U. (2007). Muscle proprioceptive feedback and spinal networks. Brain Res. Bull. 73, 155–202. doi: 10.1016/j.brainresbull.2007.03.010
Keywords: α-γ motor system, propriospinal neurons, spinal circuits, muscle and spindle, computational modeling, simulation, movement and posture
Citation: Li S, Zhuang C, Hao M, He X, Marquez JC, Niu CM and Lan N (2015) Coordinated alpha and gamma control of muscles and spindles in movement and posture. Front. Comput. Neurosci. 9:122. doi: 10.3389/fncom.2015.00122
Received: 27 April 2015; Accepted: 14 September 2015;
Published: 09 October 2015.
Edited by:
Si Wu, Beijing Normal University, ChinaReviewed by:
Da-Hui Wang, Beijing Normal University, ChinaDanke Zhang, Hangzhou Dianzi University, China
Copyright © 2015 Li, Zhuang, Hao, He, Marquez, Niu and Lan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ning Lan, School of Biomedical Engineering, Med-X Research Institute, Shanghai Jiao Tong University, 1954 Hua Shan Road, Shanghai 200030, China, ninglan@sjtu.edu.cn