|
|
||||||||
REPORT
Neurological Sciences Institute, Oregon Health and Sciences University, Beaverton, Oregon; and Department of Biomedical Engineering, Emory University and Georgia Institute of Technology, Atlanta, Georgia
Submitted 6 July 2004; accepted in final form 25 August 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
If synergies exist, what is the nature of the command signal, and furthermore, what is the biomechanical significance of each synergy? We hypothesized that the modularity in control afforded by a synergy organization must be matched by a similar reduction in degrees of freedom at the biomechanical level. To this end, we studied how synergy activation levels and biomechanical signals varied during the automatic postural response in cats to linear translation of the support surface in 16 directions in the horizontal plane.
Feline postural responses have been well characterized and are sufficiently complex for investigating how synergies can control natural behaviors. For each direction of translation, a different pattern of EMG activity and active force is elicited to stabilize balance (Macpherson 1988a, b
). The purpose of this report is to examine the synergy organization of this natural behavior and to test whether the synergies were related to the functional, biomechanical variable of endpoint force between the limb and the support surface.
| METHODS |
|---|
|
|
|---|
|
EMG responses were analyzed using a nonnegative optimization procedure similar to that described by Tresch et al. (1999)
. It was assumed that 1) any given muscle can belong to more than one synergy and 2) the muscles within a given synergy have a fixed proportional activation. For example, given three muscles, each synergy, Wi, is represented as a 3 x 1 vector (with constant elements wi1, wi2, wi3) that specifies the relative activation level of each muscle between 0 and 1 (Fig. 1A, red and blue boxes). When a synergy is recruited during a motor task, each element in the synergy vector Wi is multiplied by a scaling coefficient, ci, yielding a specific pattern of muscle activity ci·Wi (Fig. 1A, red or blue lines). Because several synergies may act on a given muscle, the net activation of that muscle is the sum of activations due to each synergy. In Fig. 1A, the total activity of muscle 1 is c1·w11 + c2·w21. In general, the ensemble EMG pattern can be expressed as the weighted sum of synergies, E
G =
Wi · ci. The pattern of each synergy's activation level over the 16 perturbation directions is expressed as a 16 x 1 vector Ci (with elements ci1, ci2. ...ci16). Ci represents the command to synergy Wi and expresses how the synergy activation varies as a function of direction (synergy tuning curves).
In the optimization procedure, the number of synergies is first specified by the user, and the elements of each synergy vector are randomly assigned initial values. A gradient search algorithm (MATLAB function fmincon.m) varies the elements of each synergy, wi1, wi2..., to minimize the difference between the actual data and the reconstructed muscle activation pattern, using additive linear combinations of the synergies. At each iteration, the reconstruction coefficients cij are found using a nonnegative least squares fit (MATLAB lsnonneg.m) of the synergies to the data. This process is continued until the least-squared errors between EMG data and the reconstructed EMG values are minimized. The same synergy vectors were found over repeated optimizations with random initial conditions.
Optimizations for each cat were run using one to six synergies. The number of synergies that best characterized the data were determined using two criteria: 1) that a large percentage (>95%) of the total variability was accounted for (VAF) over the entire run and 2) that the distribution of errors as a function of perturbation direction was roughly uniform. These two criteria ensured that each synergy chosen added a component that was necessary so that muscle activation patterns in each perturbation direction were adequately reconstructed. Whereas the addition of a synergy might increase the global VAF by only a few percentage points, synergy was included if it improved the reconstruction of specific muscle activation patterns in response to particular perturbation directions. If adding a synergy improved the reconstruction only incrementally and roughly equally across all directions, the synergy probably accounted for random variations and was not added. At least five optimizations were performed for each set of synergies, and the VAF varied by <0.01% across all optimizations with the same number of synergies.
Decomposition of ground reaction forces
To test the hypothesis that activation of each muscle synergy generates a unique direction of endpoint force, we decomposed the change in applied force at the left hindpaw into basis force vectors using the same nonnegative factorization technique describe above. The basis force vectors are analogous to muscle synergies: a set of nonorthogonal basis force vectors that are scaled by activation coefficients are extracted. Thus the nonnegative routine removes any dependence on the defined Cartesian coordinate system and essentially determines a new, nonorthogonal, force coordinate system that can reconstruct the entire force response data for all directions of translation. This analysis also allowed us to determine the number of degrees of freedom in the biomechanical output variable of force applied by the hindlimb to the ground. First, the change in active force response was computed for each x, y, and z component, by subtracting the mean background levels prior to the perturbation. Force changes of opposite sign elicit responses in different muscles; for example, loading of the limb (+Fz) elicits an active extensor response, whereas unloading (Fz) elicits an active flexor response (Macpherson 1988a, b
). Thus the positive and negative changes in force were separated, resulting in six different vectors representing the positive and negative changes in Fx, Fy, and Fz applied forces. If the change in force was positive for a given direction, the corresponding negative component was set to zero. Each force direction was normalized to its own maximum value to eliminate the magnitude differences between vertical and horizontal force changes. The absolute value of these six values for each of the 16 directions of translation were input to the nonnegative factorization. As in the EMG analysis, the number of basis force vectors was determined by optimizing how well a set of nonorthogonal basis force vectors could reconstruct the changes in applied force over all perturbation directions. These basis force vectors were converted back to force units, and the activation levels of muscle synergies were correlated to the activation levels of the basis force vectors across perturbation direction.
| RESULTS |
|---|
|
|
|---|
Number of synergies
Using both global and local criteria, four synergies were found to reliably reproduce >95% variability in EMG patterns across all perturbation directions (Fig. 1B) and to reproduce >90% variability in any single perturbation direction (Fig. 1C). Although the addition of a fifth or higher synergy increased the overall VAF, the improvement was incremental and evenly spread across all perturbation directions. In contrast, with three synergies, the VAF for some directions dropped <90% even though overall VAF exceeded that value in all cats.
Muscle synergy composition
Each synergy was characterized by a dominant subset of muscles with relatively high weights and the remainder with low or zero weights (Fig. 2A). Most muscles were activated by more than one synergy, but usually predominated in only one synergy and exhibited lower weights in the others. Moreover, each synergy was comprised of muscles crossing all of the hindlimb joints including both uni- and biarticular muscles. For example in cat Kn, gluteus medius (GLUT) was activated by synergies W1 and W3 (Fig. 2A), middle biceps femoris (BFMM) was primarily activated by W3, and posterior semimembranosus (SEMP) was activated by synergies W1, W2, and W3. Furthermore, W1 was dominated by extensors of the hip, knee, and ankle and W2 by flexors, whereas W3 and W4 tended toward higher activations in biarticular thigh muscles. Six of the recorded muscles were common across all three cats, allowing us to compare the synergy composition across subjects. The weightings (wi) for those six muscles were correlated across cats for each of the four synergies. Coefficients of determination (r2) between cat An and cats Kn and Wo, respectively, for the first synergy (W1) were 0.96 and 0.89. For W2, they were 0.71 and 0.72, for W3, they were 0.90 and 0.80, and for W4, they were 0.52 and 0.66. Thus synergy composition is robust across subjects and reflects the consistency in the EMG tuning curves of the postural response (Macpherson 1988b
).
|
The activation level of each muscle synergy (coefficients Ci) varied with perturbation direction with a single, monophasic peak at a characteristic direction. (Fig. 2B). This "synergy tuning curve" may be considered as the neural command signal to the synergy, which activates a pattern of muscles appropriate to the postural response. For each perturbation direction, a unique pattern of synergies is activated. For example, for a rightward perturbation at 270°, W1 and W4 are activated, but W2 and W3 are not. The directions and shapes of the synergy tuning curves were quite similar across cats. Peak activation for W1 occurred at 202, 180, and 225° for cats An, Kn, and Wo, respectively. Peak activation for W2 occurred at 247, 270, and 337°; peak activation for W3 occurred at 90, 90, and 135°.
Reconstruction of EMG tuning curves
Because the EMG data were reconstructed as a weighted sum of synergies at each perturbation direction (Fig. 2C), the tuning curves of each muscle could only take the shape of additive linear combinations of the synergy tuning curves, Ci. For example, the tuning curve of GLUT results from the synergy tuning curve C1 multiplied by the value of the GLUT weighting in W1 (Fig. 2, A and B; red), plus the synergy tuning curve C3 multiplied by the value of the GLUT weighting in W3 (Fig. 2, A and B; green). BFMM is activated by W3 only, and therefore, has the same single peak as C3, whereas SEMP has two peaks due to activation by both W2 and W3. Thus EMG tuning curves were partitioned according to each synergy's contribution to the activity. The reconstructed EMG tuning curves closely matched the original data (compare solid black lines and dashed black lines in Fig. 2C).
Relationship between muscle synergies and hindlimb force
The synergy activation tuning curves, Ci, were significantly correlated with the activation of individual basis force vectors (Fig. 3), suggesting a functional association between muscle synergies and limb endpoint force. The change in force applied by the limb to the ground following perturbations (Fig. 3A) was decomposed into four nonorthogonal basis force vectors that each pointed in a different direction (Fig. 3B). The magnitude of each basis force vector varied as a function of perturbation direction similar to the synergy tuning curves, and each matched-pair of synergy and force tuning curves were significantly correlated (Fig. 3C). The extensor synergy (W1) was associated with a downward and outward force, the flexor synergy (W2) with an upward and inward force, and synergy W3 primarily with a lateral force. Synergy W4 showed the lowest correlation with endpoint force vector.
|
| DISCUSSION |
|---|
|
|
|---|
The automatic postural response is highly stereotyped and repeatable (Macpherson 1988a, b
), and these features are reflected in the robustness of the synergy composition and tuning across days and across cats. We chose to examine the relation of synergies to endpoint force because maintaining balance is inherently a force control taskforces applied against the support surface act to accelerate the center of mass. Moreover, our recent study suggested that paw afferents encoding change in endpoint force at the support surface provide the sensory cues for selection of the appropriate muscle pattern, or synergy activation, in the postural response (Ting and Macpherson 2004
). The first two synergies, W1 and W2 were dominated by the antigravity-related functions of loading and unloading of the limb. It is likely that the third synergy fine-tuned the horizontal plane forces, contributing to the previously reported force constraint strategy in the response to translation, in which the applied forces tend to group around two main directions in the horizontal plane, regardless of the direction of translation (Macpherson 1988a
). The fourth synergy, W4, had the lowest correlation to active force production at the ground, with a peak activation at 337°. In this direction, the left hindlimb generated very little change in force, while the contralateral, right hindlimb played a large role in stabilizing the body by applying an unloading force directed forward and medially. Thus, rather than generating a specific endpoint force, W4 may have acted to stiffen the leg, perhaps to minimize unwanted effects from the strong flexor activation of the contralateral leg and to stabilize the pelvis. Furthermore, the mechanical interactions between the limbs could also account for the portions of the ground reaction forces that are not matched to synergy activations.
The correspondence between synergies and whole limb force is consistent not only with balance control but also with pedaling for which muscle groupings correlate with accelerations of the toe relative to the hip (Ting et al. 1999
) and vary in amplitude as a function of the sensorimotor state of the task (Ting et al. 2000
). We expect that other variableskinematic, kinetic, and impedancemay also be controlled by a muscle synergy organization, depending on the task (Ivanenko et al. 2003
; Krishnamoorthy et al. 2003a, b
). By including these variables, the number of possible degrees of freedom in biomechanical outputs is vastly increased. For example, Saltiel et al. (2001)
showed that similar kinematics, such as leg extension, could be produced by synergies with different muscle activations, although the movements might be different under different loading conditions. We anticipate that there are many synergies involved in the performance of different types of tasks.
Advantages of our technique include the extraction of physiologically interpretable components and the fact that synergy organization can be identified without a priori knowledge of function (cf. Jacobs and Macpherson 1996
; Raasch and Zajac 1999
). Nonnegative methods are more successful than strictly orthogonal methods like principal components analysis (PCA) at breaking complex objects into meaningful parts, particularly in such inherently positive-valued data as neural spike trains or muscle activations (Lee and Seung 1999
). Basis vectors only add and never subtract features, and they are nonorthogonal. Therefore each synergy's activation can be analyzed independently, whereas in PCA, many features of the components cancel each other out on data reconstruction, which means that their contribution to the final output pattern depends on the activity of other synergies. In contrast, the extracted components in nonnegative techniques are based on the compositionality of the data set, rather than holistic features. For example, when applied to images of faces, a nonnegative extraction routine generates basis vectors representing noses, ears, and eyes, whereas PCA generates components that all tend to look roughly like an entire face (Lee and Seung 1999
). A limitation of our technique is that proper identification of muscle synergies requires a high degree of complexity in the sampling of muscle coordination patterns. If, for example, only two perturbation directions had been studied, only two muscle synergies would have been identified. For our paradigm, similar synergies could be extracted from as few as 8 of the 16 perturbations directions.
Our approach may be useful for understanding neural feedback control and biomechanical functions of muscle synergies in natural behaviors. The idea that muscle synergies produce whole limb functions is consistent with the findings that sensory signals arising from the spinal cord can represent integrative variables such as limb length and limb orientation (Bosco et al. 2000
). It is also consistent with recent studies suggesting that focal stimulation of the spinal cord may produce consistent patterns of muscle activation (Ivanenko et al. 2003
; Lemay and Grill 2004
; Saltiel et al. 2001
; Tresch et al. 2002
). Thus muscle synergy activations could be modulated directly by such sensory signals to control overall limb function, allowing higher CNS levels to encode task-level variables without regard for individual muscle activations. Our synergy analysis may provide a pattern that varies more closely with neural signals than biomechanics or individual muscle activation patterns. It remains to be seen whether the synergy organization can be used to control a wide range of activities, and the extent to which synergies are shared across different types of tasks.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: L. H. Ting, Wallace H. Coulter Dept. of Biomedical Engineering at Georgia Tech and Emory Univ., 313 Ferst Dr., Atlanta, GA 30322-0535 (E-mail: lting{at}emory.edu)
| REFERENCES |
|---|
|
|
|---|
Bosco G, Poppele RE, and Eian J. Reference frames for spinal proprioception: limb endpoint based or joint-level based? J Neurophysiol 83: 29312945, 2000.
Buchanan TS, Almdale DPJ, Lewis JL, and Rymer WZ. Characteristics of synergic relations during isometric contractions of human elbow muscles. J Neurophysiol 56: 12251241, 1986.
Ivanenko YP, Grasso R, Zago M, Molinari M, Scivoletto G, Castellano V, Macellari V, and Lacquaniti F. Temporal components of the motor patterns expressed by the human spinal cord reflect foot kinematics. J Neurophysiol 90: 35553565, 2003.
Jacobs R and Macpherson JM. Two functional muscle groupings during postural equilibrium tasks in standing cats. J Neurophysiol 76: 24022411, 1996.
Krishnamoorthy V, Goodman S, Zatsiorsky V, and Latash ML. Muscle synergies during shifts of the center of pressure by standing persons: identification of muscle modes. Biol Cybern 89: 152161, 2003a.[CrossRef][ISI][Medline]
Krishnamoorthy V, Latash ML, Scholz JP, and Zatsiorsky VM. Muscle synergies during shifts of the center of pressure by standing persons. Exp Brain Res 152: 281292, 2003b.[CrossRef][ISI][Medline]
Lee DD and Seung HS. Learning the parts of objects by non-negative matrix factorization. Nature 401: 788791, 1999.[CrossRef][Medline]
Lee WA. Neuromotor synergies as a basis for coordinated intentional action. J Motor Behav 16: 135170, 1984.[ISI][Medline]
Lemay MA and Grill WM. Modularity of motor output evoked by intraspinal microstimulation in cats. J Neurophysiol 91: 502514, 2004.
Macpherson JM. Strategies that simplify the control of quadrupedal stance. I. Forces at the ground. J Neurophysiol 60: 204217, 1988a.
Macpherson JM. Strategies that simplify the control of quadrupedal stance. II. Electromyographic activity. J Neurophysiol 60: 218231, 1988b.
Macpherson JM. How flexible are muscle synergies? In: Motor Control: Concepts and Issues, edited by Humphrey DR and Freund H-J. New York: Wiley Press, 1991, p. 3347.
Raasch CC and Zajac FE. Locomotor strategy for pedaling: muscle groups and biomechanical functions. J Neurophysiol 82: 515525, 1999.
Saltiel P, Wyler-Duda K, D'Avella A, Tresch MC, and Bizzi E. Muscle synergies encoded within the spinal cord: evidence from focal intraspinal NMDA iontophoresis in the frog. J Neurophysiol 85: 605619, 2001.
Soechting JF and Lacquaniti F. An assessment of the existence of muscle synergies during load perturbations and intentional movements of the human arm. Exp Brain Res 74: 535548, 1989.[ISI][Medline]
Ting LH, Kautz SA, Brown DA, and Zajac FE. Phase reversal of biomechanical function pairs and muscle activity in backward pedaling. J Neurophysiol 81: 544551, 1999.
Ting LH, Kautz SA, Brown DA, and Zajac FE. Contralateral movement and extensor force generation alter flexion phase muscle coordination in pedaling. J Neurophysiol 83: 33513365, 2000.
Ting LH and Macpherson JM. Ratio of shear to load ground-reaction force may underlie the directional tuning of the automatic postural response to rotation and translation. J Neurophysiol 92: 808823, 2004.
Tresch MC, Saltiel P, and Bizzi E. The construction of movement by the spinal cord. Nat Neurosci 2: 162167, 1999.[CrossRef][ISI][Medline]
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
S. A. Winges, S. E. Eonta, J. F. Soechting, and M. Flanders Multi-Digit Control of Contact Forces During Rotation of a Hand-Held Object J Neurophysiol, April 1, 2008; 99(4): 1846 - 1856. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Drew, J. Kalaska, and N. Krouchev Muscle synergies during locomotion in the cat: a model for motor cortex control J. Physiol., March 1, 2008; 586(5): 1239 - 1245. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Overduin, A. d'Avella, J. Roh, and E. Bizzi Modulation of Muscle Synergy Recruitment in Primate Grasping J. Neurosci., January 23, 2008; 28(4): 880 - 892. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Torres-Oviedo and L. H. Ting Muscle Synergies Characterizing Human Postural Responses J Neurophysiol, October 1, 2007; 98(4): 2144 - 2156. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nishikawa, A. A. Biewener, P. Aerts, A. N. Ahn, H. J. Chiel, M. A. Daley, T. L. Daniel, R. J. Full, M. E. Hale, T. L. Hedrick, et al. Neuromechanics: an integrative approach for understanding motor control Integr. Comp. Biol., July 1, 2007; 47(1): 16 - 54. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. M. Herter, I. Kurtzer, D. W. Cabel, K. A. Haunts, and S. H. Scott Characterization of Torque-Related Activity in Primary Motor Cortex During a Multijoint Postural Task J Neurophysiol, April 1, 2007; 97(4): 2887 - 2899. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. F. Giszter, M. R. Davies, and V. Graziani Motor Strategies Used by Rats Spinalized at Birth to Maintain Stance in Response to Imposed Perturbations J Neurophysiol, April 1, 2007; 97(4): 2663 - 2675. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. G. Larson and J. T. Stern Jr Humeral retractor EMG during quadrupedal walking in primates J. Exp. Biol., April 1, 2007; 210(7): 1204 - 1215. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Klein Breteler, K. J. Simura, and M. Flanders Timing of Muscle Activation in a Hand Movement Sequence Cereb Cortex, April 1, 2007; 17(4): 803 - 815. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Kurtzer, T. M. Herter, and S. H. Scott Nonuniform Distribution of Reach-Related and Torque-Related Activity in Upper Arm Muscles and Neurons of Primary Motor Cortex J Neurophysiol, December 1, 2006; 96(6): 3220 - 3230. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Schepens and T. Drew Descending Signals From the Pontomedullary Reticular Formation Are Bilateral, Asymmetric, and Gated During Reaching Movements in the Cat J Neurophysiol, November 1, 2006; 96(5): 2229 - 2252. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Daley and A. A. Biewener Running over rough terrain reveals limb control for intrinsic stability PNAS, October 17, 2006; 103(42): 15681 - 15686. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Krouchev, J. F. Kalaska, and T. Drew Sequential Activation of Muscle Synergies During Locomotion in the Intact Cat as Revealed by Cluster Analysis and Direct Decomposition J Neurophysiol, October 1, 2006; 96(4): 1991 - 2010. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Torres-Oviedo, J. M. Macpherson, and L. H. Ting Muscle Synergy Organization Is Robust Across a Variety of Postural Perturbations J Neurophysiol, September 1, 2006; 96(3): 1530 - 1546. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. P. Ivanenko, R. E. Poppele, and F. Lacquaniti Motor Control Programs and Walking Neuroscientist, August 1, 2006; 12(4): 339 - 348. [Abstract] [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. C. Tresch, V. C. K. Cheung, and A. d'Avella Matrix Factorization Algorithms for the Identification of Muscle Synergies: Evaluation on Simulated and Experimental Data Sets J Neurophysiol, April 1, 2006; 95(4): 2199 - 2212. [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] |
||||
![]() |
I. Kurtzer, J. A. Pruszynski, T. M. Herter, and S. H. Scott Primate Upper Limb Muscles Exhibit Activity Patterns That Differ From Their Anatomical Action During a Postural Task J Neurophysiol, January 1, 2006; 95(1): 493 - 504. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Gruneberg, J. Duysens, F. Honegger, and J.H.J. Allum Spatio-Temporal Separation of Roll and Pitch Balance-Correcting Commands in Humans J Neurophysiol, November 1, 2005; 94(5): 3143 - 3158. [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] |
||||
![]() |
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] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||