QUICK SEARCH:   [advanced] Author: Keyword(s):
Home     Help     Feedback     Subscriptions     Archive     Search     Table of Contents

First published online August 17, 2006
Journal of Experimental Biology 209, 3345-3357 (2006)
Published by The Company of Biologists 2006
doi: 10.1242/jeb.02340

This Article Summary Figures Only Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Arampatzis, A. Articles by Brüggemann, G.-P. Search for Related Content PubMed PubMed Citation Articles by Arampatzis, A. Articles by Brüggemann, G.-P.

Influence of the muscle-tendon unit's mechanical and morphological properties on running economy

Adamantios Arampatzis^*, Gianpiero De Monte, Kiros Karamanidis, Gaspar Morey-Klapsing, Savvas Stafilidis and Gert-Peter Brüggemann

Adamantios Arampatzis, German Sport University of Cologne, Institute of Biomechanics and Orthopaedics, Carl-Diem-Weg 6, 50933 Cologne, Germany

^* Author for correspondence (e-mail: Arampatzis{at}dshs-koeln.de)

Accepted 18 May 2006

	Summary

TOP Summary Introduction Materials and methods Results Discussion References

 
The purpose of this study was to test the hypothesis that runners having different running economies show differences in the mechanical and morphological properties of their muscle-tendon units (MTU) in the lower extremities. Twenty eight long-distance runners (body mass: 76.8±6.7 kg, height: 182±6 cm, age: 28.1±4.5 years) participated in the study. The subjects ran on a treadmill at three velocities (3.0, 3.5 and 4.0 m s^-1) for 15 min each. The consumption was measured by spirometry. At all three examined velocities the kinematics of the left leg were captured whilst running on the treadmill using a high-speed digital video camera operating at 250 Hz. Furthermore the runners performed isometric maximal voluntary plantarflexion and knee extension contractions at eleven different MTU lengths with their left leg on a dynamometer. The distal aponeuroses of the gastrocnemius medialis (GM) and vastus lateralis (VL) were visualised by ultrasound during plantarflexion and knee extension, respectively. The morphological properties of the GM and VL (fascicle length, angle of pennation, and thickness) were determined at three different lengths for each MTU. A cluster analysis was used to classify the subjects into three groups according to their consumption at all three velocities (high running economy, N=10; moderate running economy, N=12; low running economy, N=6). Neither the kinematic parameters nor the morphological properties of the GM and VL showed significant differences between groups. The most economical runners showed a higher contractile strength and a higher normalised tendon stiffness (relationship between tendon force and tendon strain) in the triceps surae MTU and a higher compliance of the quadriceps tendon and aponeurosis at low level tendon forces. It is suggested that at low level forces the more compliant quadriceps tendon and aponeurosis will increase the force potential of the muscle while running and therefore the volume of active muscle at a given force generation will decrease.

Key words: tendon elasticity, tendon stiffness, running economy, ultrasonography, running kinematics, energy exchange, skeletal muscle

	Introduction

TOP Summary Introduction Materials and methods Results Discussion References

 
The literature defines running economy as the rate of oxygen consumption per unit body mass when running at a constant pace (Daniels et al., 1978; Cavanagh and Kram, 1985; Williams and Cavanagh, 1987). It has been shown that distance runners demonstrate significant differences in the rate of oxygen consumption while running at the same velocity (Williams and Cavanagh, 1987; Daniels and Daniels, 1992). Furthermore, at a wide range of velocities there is a close relationship between metabolic energy cost and mechanical power (Bijker et al., 2001). As a consequence, many studies on running economy have been motivated by the suggestion that biomechanical factors might explain the differences of running economy between individuals (Cavanagh and Williams, 1982; Williams and Cavanagh, 1987; Kyröläinen et al., 2001). However, in general the relationship observed between biomechanical factors and running economy is weak and it has been concluded that descriptive kinematic and kinetic parameters alone cannot explain the complexity of running economy (Williams and Cavanagh, 1987; Martin and Morgan, 1992; Kyröläinen et al., 2001).

It has been suggested that variables that describe muscle force production (i.e. force-length-velocity relationship and activation) are probably more suitable for explaining running economy (Martin and Morgan, 1992). From a mechanical point of view there are two main issues that can affect the force-length-velocity relationship and the activation of the muscles while running. The mechanical advantages of the muscles (ratio of an agonist muscle group moment arm to that of the ground reaction force acting about a joint) may affect the force production in relation to the active muscle volume. For example, it is well accepted in the literature that small mammals show lower effective mechanical advantages during running than larger mammals (Biewener, 1989; Biewener, 1990), which lead to a decrease in the force production per active muscle volume. Recently Biewener et al. reported that differences in the effective mechanical advantages between walking and running explain the higher energy transport during running compared to walking (Biewener et al., 2004). A second issue that can influence the force-length-velocity relationship and activation of the muscles is the non-rigidity of the tendon and aponeurosis (Bobbert, 2001; Hof et al., 2002; Roberts, 2002). A higher compliance will allow the muscle fibres to contract at lower shortening velocities than the whole muscle-tendon unit (MTU) (Ettema et al., 1990a; Ettema et al., 1990b), and as a consequence of the force-velocity relationship their force-generating potential will be higher (Hof et al., 1983; Hof et al., 2002; Bobbert, 2001). Furthermore, due to the non-rigidity of the tendon and aponeurosis, when the MTU is elongated, strain energy can be stored that is independent of metabolic processes (Roberts, 2002). This way the whole mechanical energy produced during the shortening of the MTU can be enhanced (Alexander and Bennet-Clark, 1977; Ker et al., 1987; de Haan et al., 1989; Ettema, 1996; Roberts et al., 1997).

Despite these phenomena being well known, there is currently no study showing the influence of the mechanical properties of the MTU on running economy in humans, nor has the role of muscle architecture in enhancing running economy been studied. Recent in vivo studies investigating the mechanical and morphological properties of the MTU have demonstrated differences between athletes pertaining to different disciplines (Kawakami et al., 1993; Kawakami et al., 1995; Abe et al., 2000). Abe et al., for example, compared sprinters with long-distance runners, and found that the sprinters had longer fascicles and lower pennation angles in the mm. vastus lateralis and gastrocnemius (Abe et al., 2000). Longer muscle fascicles can exhibit higher shortening velocities and mechanical powers than shorter fascicles. In general, literature reports recognised a significant correlation between the fascicle lengths of the lower extremity muscles (vastus lateralis, gastrocnemii) and sprint performance (Kumagai et al., 2000; Abe et al., 2001).

While running the muscles acting around the ankle and knee joints (i.e. triceps surae and quadriceps femoris) contribute more than 70% of the total mechanical work (Winter, 1983; Sasaki and Neptune, 2005). Therefore it can be argued that they belong to the main muscles expending energy during sub-maximal running. Furthermore Sasaki and Neptune (Sasaki and Neptune, 2005) reported that during sub-maximal running the energy stored in the tendon and aponeurosis of the triceps surae and quadriceps femoris MTU is about 75% of the energy stored in all tendons of the muscoloskeletal system. Reports on the influence of the non-rigidity of the tendon and aponeurosis on the effectivity of muscle force production (Ettema et al., 1990a; Ettema et al., 1990b; Roberts et al., 1997; Hof et al., 2002) and the effect of fascicle lengths on performance of sport activities (Kumagai et al., 2000; Abe et al., 2001) reveal the expectation that running economy may be affected by the mechanical and morphological properties of the triceps surae and quadriceps femoris MTUs. Basing on the above expectation it can be hypothesised that runners having different running economy would show differences in the mechanical and morphological properties of their MTUs in the lower extremities. Therefore we examined the mechanical properties and the architecture of the MTUs of the lower extremities from runners displaying different running economies, together with their running kinematics.

	Materials and methods

TOP Summary Introduction Materials and methods Results Discussion References

 
Subjects
Twenty eight male long-distance runners (body mass 76.8±6.7 kg, height 182±6 cm, age 28.1±4.5 years), all regularly participating in running competitions locally, took part in the study. The runners gave their informed written consent to the experimental procedure complying with the rules of the local scientific board. All subjects performed endurance running training between 4 and 9 times per week. The training volume ranged from 40 to 120 km week^-1. None of the subjects had a history of neuromuscular or musculoskeletal impairments at the time of the study that could affect their running technique.

Oxygen consumption
After a warm-up period of 5 min at a running velocity of 3.0 m s^-1 the subjects ran three different velocities for 15 min on a treadmill in the same order (3.0, 3.5 and 4.0 m s^-1) wearing their own running shoes. Oxygen consumption (; ml kg^-1 min^-1) was measured during this 15 min period using a breath-by-breath spirometer (Jaeger Oxycon , Hoechberg, Germany). The spirometer was calibrated before each session by means of a two-point calibration using environment air and a gas mixture (5.5% CO₂, 0% O₂, balance N₂). The volume sensor was calibrated by means of a manual 2 litre syringe. The accuracy values provided by the manufacturer were 0.01% for O₂ and CO₂ with a drift of 0.02% per hour, and 0.02% for volume. For each velocity the average value of the was calculated from 4 min of running at steady state (Fig. 1, min 10-14). There were 10 min rests between running at each velocity test. Blood samples were taken from the earlobe directly after finishing each velocity test within the first 30 s of the rest to determine blood lactate concentration, which helps to identify differences in the anaerobic energy cost between the examined subjects that might occur.

View larger version (22K): [in this window] [in a new window]   Fig. 1. Oxygen consumption () of a runner at all three velocities. The two vertical lines delimit the period during which oxygen consumption was determined.

Measurement of maximal isometric ankle and knee joint moment
The subjects performed isometric maximal voluntary ankle plantarflexion and knee extension contractions (MVC) of their left leg on two separate test days. The warm-up consisted of 2-3 min performing submaximal isometric contractions and three MVCs. Afterwards the subjects performed isometric maximal voluntary ankle plantarflexion or knee extension contractions at eleven different ankle-knee and knee-hip joint angle configurations, respectively (Table 1), on a dynamometer (Biodex Medical Systems. Inc., Shirley, NY, USA). Different joint angle configurations were chosen in order to examine triceps surae and quadriceps femoris muscle strength potential over the whole range of achievable MTU lengths. The different joint angle configurations were applied in random order. 3 min rest between contractions were allowed. The subjects were instructed and encouraged to produce a maximal isometric moment and to hold it for about 2-3 s.

Before each MVC the axis of rotation of the dynamometer was carefully aligned with the axis of rotation of the ankle and knee joints. The axis of rotation of the ankle joint was defined to be parallel to the axis of the dynamometer and passing through the midpoint of the line connecting both malleoli. In the same way the axis of rotation of the knee joint was defined to be parallel to the axis of the dynamometer and passing through the midpoint of the line connecting the lateral and medial femoral condyles. During the contraction the axes clearly shifted away from each other. Therefore, kinematic data were recorded using a Vicon 624 system (Vicon Motion Systems, Oxford, UK) with eight cameras operating at 120 Hz to calculate the resultant joint moments. To calculate the lever arm of the ankle joint during ankle plantarflexion the centre of pressure under the foot was determined by means of a flexible pressure distribution insole (Pedar, Novel GmbH, Munich, Germany) operating at 99 Hz. The compensation of moments due to gravitational forces was done for all subjects before each ankle plantarflexion or knee extension contraction. The exact method for calculating the resultant joint moments has been previously described (Arampatzis et al., 2004; Arampatzis et al., 2005b).

The moments arising from antagonistic coactivation during the ankle plantarflexion and knee extension efforts were quantified by assuming a linear relationship between surface electromyography (EMG) amplitude of the ankle dorsiflexor or knee flexor muscles and moment (Baratta et al., 1988). This was established by measuring EMG and moment during one relaxed condition and two submaximal ankle dorsiflexion or knee flexion contractions at each joint angle configuration (Mademli et al., 2004). Therefore, in the text below, maximal knee and ankle joint moments refer to the maximal joint moment values considering the effect of gravitational forces, the effect of the joint axis alignment relative to the dynamometer axis and the effect of the antagonistic moment on the moment measured at the dynamometer.

Measurement of EMG-activity during isometric contractions
Bipolar EMG lead-offs with pre-amplification (analogue RC-filter 10-500 Hz bandwidth, Biovision, Wehrheim, Germany) and adhesive surface electrodes (blue sensor; Medicotest, Ballerup, Denmark) were used to analyse muscle activity. Before placing the electrodes the skin was carefully prepared (shaved and cleaned with alcohol) to reduce skin impedance. The electrodes were positioned above the midpoint of the muscle belly as assessed by palpation, parallel to the presumed direction of the muscle fibres. The inter-electrode distance was 2 cm. The activation of the triceps surae muscle was assessed from the EMGs of the gastrocnemius medialis (GM), gastrocnemius lateralis (GL) and soleus (SOL). During knee extension the EMG-activities of the vastus lateralis (VL), vastus medialis (VM) and rectus femoris (RF) were analysed. The EMG signals were recorded at 1080 Hz by the Vicon system. Before starting the experiment, tests including submaximal and maximal isometric contractions for each muscle group were undertaken to determine whether an adequate signal was obtained from each muscle and to adjust the amplifier gains. The EMG signal from each muscle was checked online for artefacts due to mechanical causes by passively shaking the leg. Additionally, several functional tests (i.e. hopping in place) were undertaken to determine whether a good signal was obtained from each muscle. The preparation was renewed when such artefacts were observed. All isometric contractions at the knee or the ankle joint were performed within one testing session. No electrode replacement or re-adjusting of the EMG pre-amplification gain was done during the measurements.

The EMG-activity is described by the root mean square (RMS) of the raw signals for a time interval of 1000 ms at peak joint moment. The RMS from each muscle was normalised to the individual maximal RMS value of each muscle for each subject during the eleven isometric contractions. In order to determine the EMG-activity of the ankle plantarflexor and knee extensor muscles, the normalised RMS of the examined muscles were averaged and weighted by their physiological cross sectional areas (PCSA). For the TS, a PCSA ratio of 6:2:1 for the SOL, GM, GL (Out et al., 1996) and for the QF, a PCSA ratio of 0.92, 1.00 and 0.72 for the RF, VL and VM (Herzog et al., 1990) were assumed.

View larger version (66K): [in this window] [in a new window]   Fig. 2. Muscle ultrasound images at rest. Left, gastrocnemius medialis (GM); right, vastus lateralis (VL). L, fascicle length; d, muscle thickness; , pennation angle; digitalised cross-point, insertion of the fascicle into the deeper aponeurosis.

The tendon force was calculated by dividing the ankle or knee joint moment by the corresponding tendon moment arm. The tendon moment arms of the Achilles tendon and the patellar tendon were calculated using the data provided by Maganaris et al. (Maganaris et al., 1998) and Herzog and Read (Herzog and Read, 1993), respectively. The stiffness of the tendon and aponeurosis was calculated by means of linear regression equations. The stiffness (normalised) represented the relationship between the tendon force and the strain of the tendon and aponeurosis between 45% and 100% of the maximal tendon force. We used the normalised stiffness because the amount of elongation of a tendon at a given exerted force depends on the rest length of the tendon. In vitro studies examining the elongation of the tendon used the same rest lengths. In vivo it is very difficult to do this, because of the differences in the anthropometrical characteristics of the subjects and also the differences in the localisation of the ultrasound probe (it is practically impossible to place the ultrasound probe exactly at the same position at all experiments). Differences in the rest length would influence the calculated stiffness (relationship between tendon force and elongation) of the tendon and aponeurosis (Rack and Westbury, 1984; Muramatsu et al., 2001; Arampatzis et al., 2005a; Stafilidis et al., 2005). This fact makes it difficult to compare the stiffness between different subjects or groups. Recently it has been reported that the strain measured at the myotendinous junction and at the muscle belly is similar (Muramatsu et al., 2001; Arampatzis et al., 2005a; Stafilidis et al., 2005), so the choice of the cross-point does not effect the calculated strain. Therefore in the present study we used the normalised stiffness. The linearity between tendon force and strain was checked using the coefficient of determination (r²). The coefficients of determination were reasonably high (r²=0.98 to 0.99). The energy storage capacity of the tendon and aponeurosis during the maximal voluntary contraction was calculated as the integral of the tendon force over the tendon strain.

Measurement of muscle architecture
The muscle architecture of the GM and VL (fascicle length, angle of pennation, and thickness) was determined at three different lengths for each muscle-tendon unit. Following ankle-knee joint angle combinations were chosen for the GM: Position 1: ankle angle 90°, knee angle 180°; Position 2: ankle angle 110°, knee angle 160°; Position 3: ankle angle 120°, knee angle 110°. The corresponding knee-hip joint angle combinations for the VL were: Position 1: knee angle 80°, hip angle 140°; Position 2: knee angle 115°, hip angle 140°; Position 3: knee angle 170°, hip angle 140°. All measurements were done on the relaxed muscle at the cited positions. The pennation angles of the GM and VL were measured as the angle of insertion of the muscle fascicles into the deep aponeurosis. The fascicle length was defined as the length of the fascicular path between the insertions of the fascicle into the upper and deeper aponeurosis (Fig. 2). The muscle thickness was defined as the distance between the deeper and upper aponeurosis.

Statistics
The examined runners were divided into groups by means of a cluster analysis, based on their oxygen consumption (ml kg^-1 min^-1) at all three velocities (3.0, 3.5 and 4.0 m s^-1). The cluster analysis revealed three relatively homogeneous groups (group 1: high running economy, N=10; group 2: moderate running economy, N=12; group 3: low running economy, N=6). All parameters (running kinematics, mechanical and morphological properties of the lower extremity MTUs) were checked for differences between groups using one-way analysis of variance (ANOVA) and Tukey post-hoc comparisons. The level of significance was set to P=0.05.

	Results

TOP Summary Introduction Materials and methods Results Discussion References

 
Metabolic parameters
The runners of the three groups showed statistically significant differences (P<0.05) in their at all examined velocities (Table 2). Group 1 was the most economical, followed by groups 2 and 3. No differences between groups in body mass, body height or lactate concentration at each running velocity were found (Table 2).

Running kinematics
Fig. 3 shows the average ankle, knee, and hip angles at 3.0, 3.5 and 4.0 m s^-1 for all three groups. The shape of the curves, as well as the angular and temporal parameters at all velocities, showed no statistically significant (P>0.05) differences between groups (Table 3). For all three velocities the average duty factor ranged from 36% to 39% and did not show any statistically significant (P>0.05) difference between groups (Table 3).

View larger version (31K): [in this window] [in a new window]   Fig. 3. Average values of the ankle (AA) knee (KA) and hip (HA) angles during running for all velocities and groups. The x-axis is normalised to the stride cycle (from heel strike to the next heel strike of the same foot). The vertical lines separate the contact and swing phases (mean ± 1 s.d.). Group 1, high running economy (N=10). Group 2, moderate running economy (N=12). Group 3, low running economy (N=6). (A) Velocity 1 (3.0 m s-1), (B) velocity 2 (3.5 m s-1), (C) velocity 3 (4.0 m s-1).

Triceps surae and quadriceps femoris muscle-tendon units
The maximal calculated tendon force, maximal plantarflexion moment, normalised stiffness and energy storage capacity during the MVC of the triceps surae tendon and aponeurosis were highest (P<0.05) for group 1 (Table 4). In contrast, the maximal strain of the tendon and aponeurosis showed no significant (P>0.05) differences between the groups (Table 4). The three examined groups showed significant differences (P<0.05) in their mechanical properties of the quadriceps femoris tendon and aponeurosis (Table 5). Although the maximal knee extension moment, the maximal calculated tendon force, and the normalised stiffness showed no significant differences (P>0.05), the maximal strain and the energy storage capacity of the tendon and aponeurosis during the MVC were highest (P<0.05) for group 1 (Table 5). Furthermore, the normalised stiffness of the quadriceps tendon and aponeurosis of the group 1 runners was significantly lower (P<0.05) at low level forces (up to 45% MVC) compared to the other two groups (Fig. 4). The calculated stiffness of the quadriceps tendon and aponeurosis from 0 to 45% MVC for all the three groups was: 29.5±4.8 kN strain^-1, 42.5± 6.2 kN strain^-1 and 38.1±5.9 kN strain^-1, respectively.

View larger version (22K): [in this window] [in a new window]   Fig. 4. Force-strain relationship of the quadriceps femoris tendon and aponeurosis during maximal knee extension effort for the three groups (mean ± 1 s.d.). Group 1, high running economy; Group 2, moderate running economy; Group 3, low running economy. *Statistically significant differences in maximal tendon strain between group 1 and the other two groups (P<0.05).

As expected, on account of the length of the GM MTU at the examined ankle-knee angle combinations, the fascicle length of the GM decreased and the angle of pennation increased from position 1 to position 3 (Table 6). The thickness showed nearly constant values for all three examined positions. The morphological parameters: fascicle length, ratio (fascicle length/tibia length), angle of pennation and thickness (Table 6) showed no statistically significant differences between groups at any position (P>0.05). Similar to the results obtained for the GM, while the fascicle length of the VL decreased and its angle of pennation increased from position 1 to position 3, its thickness showed nearly constant values (Table 7). The comparison of all three groups, as for the GM, revealed no significant (P>0.05) differences in fascicle length, ratio (fascicle length/femur length), angle of pennation or thickness at any examined position (Table 7).

The maximal plantarflexion moment was significantly higher (P<0.05) for group 1 at seven of the eleven examined positions (Fig. 5). However at four positions (i.e. pronounced ankle plantarflexion and knee flexion combinations) there were no differences between the three groups. Comparison of the EMG-activity during the maximal plantarflexion effort revealed no significant differences between groups at any ankle/knee joint combination (P>0.05) (Fig. 5). All three runner groups exhibited similar knee extension moments and EMG-activity in the quadriceps femoris muscles at all knee/hip angle combinations (P>0.05) (Fig. 6).

View larger version (55K): [in this window] [in a new window]   Fig. 5. Normalised RMS values of the EMG signal (A) for the triceps surae muscle (gastrocnemius medialis, gastrocnemius lateralis and soleus) and ankle joint moment (B) during isometric maximal voluntary ankle plantarflexion contraction at 11 different joint angle configurations for the examined groups (means ±1 s.d.). RMS values from each subject were normalised to the highest RMS value measured over all joint angle configurations. Group 1, high running economy; Group 2, moderate running economy; Group 3, low running economy. *Statistically significant differences between group 1 and the other two groups (P<0.05).

View larger version (57K): [in this window] [in a new window]   Fig. 6. Normalised RMS values of the EMG signal (A) for the quadriceps femoris muscles (vastus lateralis, vastus medialis and rectus femoris) and knee joint moment (B) during isometric maximal voluntary knee extension contraction at 11 different joint angle configurations for the examined groups (means ± 1 s.d.). RMS values from each subject were normalised to the highest RMS value measured over all joint angle configurations. Group 1, high running economy; Group 2, moderate running economy; Group 3, low running economy. *Statistically significant differences between group 1 and the other two groups (P<0.05).

	Discussion

TOP Summary Introduction Materials and methods Results Discussion References

 
The main findings of this study were that the most economical runners showed (i) a higher contractile strength of the triceps surae muscle tendon unit, (ii) a higher compliance of the quadriceps tendon and aponeurosis at low level forces, (iii) a higher energy storage capacity during MVC in both the triceps surae and the quadriceps muscle-tendon units, (iv) no differences in the muscle architecture of the GM and VL and (v) no differences in the kinematic characteristics of running.

It is noteworthy, however, that although groups 2 and 3 (moderate and low economy runners) showed clearly different oxygen consumptions at all three examined velocities they did not show differences in the mechanical (maximal joint moment, maximal calculated tendon force, maximal strain, normalised stiffness, energy storage capacity) properties of either the triceps surae or the quadriceps femoris MTUs. Thus the results demonstrate that the differences found in the mechanical properties alone between the three groups, cannot entirely explain the differences in running economy for all examined subjects and support the opinion (Williams and Cavanagh, 1987; Lake and Cavanagh, 1996; Kyrölainen et al., 2001) that a global explanation of running economy is very complex; many variables may affect the existing differences in running economy between individuals. In the present study, however, the most economical runners showed differing values in important mechanical properties at the triceps surae as well as in the quadriceps femoris MTUs. The different mechanical properties found in the lower extremities of the most economical runners may affect intrinsic muscle properties for force and energy production (i.e. force-length-velocity relationship, activation and energy storage and return during the running task) and may contribute to running economy.

Another point that should be discussed is that the examined runners showed a wide range of training volumes (4-9 times per week, 40-120 km week^-1). Although some studies reported no differences in running economy between trained and untrained subjects (Dolgener, 1982), or no relationship between running economy and training volume (Pate et al., 1992; Weston et al., 2000), others revealed that trained subjects are more economical than untrained subjects and elite runners more economical than sub-elite runners (Morgan et al., 1995). Therefore to examine whether the training volume of the examined runners could have an effect on the differences found in running economy we compared the training volumes (km week^-1) among the runner groups. No significant differences in training volume between the groups could be found (P>0.05) (group 1: 76.0±28.3 km week^-1, group 2: 74.1±13.4 km week^-1, group 3: 72.5±20.6 km week^-1) indicating that training volume was not responsible for the differences in running economy observed at our subject groups. Furthermore there was no relationship between training volume and _O 2 at all examined velocities (r²=0.08-0.09, P>0.05).

Kinematic characteristics
There were no detectable differences between groups in contact time, swing time, duty factor or stride frequency at any running velocity. Furthermore, the range of motion of the lower extremity's joints during running at all three examined velocities showed no statistically significant differences between groups. These data suggest that the differences in running economy between the three groups are not related to different stride frequencies or to different kinematic characteristics. Similar conclusions were reported in earlier studies examining kinematic parameters in runners with different running economies (Williams and Cavanagh, 1987; Kyrölainen et al., 2001). In summary, the findings of the kinematic analysis show that in a homogeneous group of runners, differences in running economy are not related to their kinematic parameters. As the effective mechanical advantage was not examined in the present study it is not possible to glean any information regarding this parameter for the three groups. However, it is unlikely that groups of runners showing no differences in running kinematics and body mass would show differences in the effective mechanical advantage during running.

Higher compliance of tendon at low level forces in the quadriceps femoris
The maximal strain of the quadriceps tendon and aponeurosis was higher (about 38%) in the high running economy group (group 1) compared to the other two groups. However at higher exerted forces the normalised stiffness showed no significant differences between the groups (45-100% of the MVC). This means that the quadriceps tendon and aponeurosis of the group 1 runners are more compliant at low level forces (up to 45% MVC) in comparison to the other two groups (Fig. 4).

No significant differences in joint kinematics between the three groups could be identified at any running velocity. Thus it can be argued that the velocities of the triceps surae and quadriceps femoris MTUs would not show significant differences between groups either. While running, during the first part of the contact phase the triceps surae and quadriceps femoris MTU are lengthening and so does the tendon and aponeurosis due to the developed forces (Hof et al., 2002; Sasaki and Neptune, 2005). During the second part of the contact phase the shortening velocity of the muscle fibres (contractile element, CE) is lower than the shortening velocity of the MTU due to the additional shortening of the tendon and aponeurosis (Roberts, 2002; Hof et al., 2002). Therefore at submaximal running intensities (velocities from 3.0 to 4.0 m s^-1) a more compliant quadriceps tendon and aponeurosis at lower force levels will increase the elongation of the series elastic element (SEE) during the first part of the contact phase and consequently decrease the shortening velocity of the CE during the second part of the contact phase (Biewener and Roberts, 2000; Bobbert, 2001). This way the CE increases its force potential due to the force-velocity relationship. A higher force potential of the CE would decrease the volume of active muscle at a given force or a given rate of force generation and consequently would decrease the cost of force production as well (Crow and Kushmerick, 1982; Roberts et al., 1997; Roberts et al., 1998a; Roberts et al., 1998b).

The VL muscle morphology (fascicle length, angle of pennation, muscle thickness) at rest did not differ between the runner groups. This suggests that the working range (width of the force-length relationship) of the VL is similar for all groups. Furthermore the maximal knee extension joint moments were similar at all analysed knee/hip joint angle configurations for the three groups of runners. The activation level during the maximal voluntary knee extension seems not to influence the above findings regarding the knee extension moments because the EMG-activities of the quadriceps femoris muscles did not differ between groups. The lack of differences in VL muscle architecture as well as in the maximal joint moments at different lengths of the quadriceps femoris MTU therefore supports the tenet that the other heads of the quadriceps femoris will also display similar fascicle lengths at all three groups. Thus we can conclude that the morphological properties of the quadriceps femoris muscles were not responsible for the differences in running economy in the examined runners.

It is well known from animal experiments that tendons reduce the mechanical work done by their muscle fibres in each step as a result of their elasticity (Biewener and Baudinette, 1995; Biewener et al., 1998; Biewener and Gillis, 1999) and reduce the metabolic cost of locomotion. During the loading phase, part of the mechanical energy coming from the mammal's body is stored as strain energy in the tendon and converted again into mechanical energy of the mammal's body in the following shortening phase (Biewener and Baudinette, 1995; Biewener and Roberts, 2000). Sasaki and Neptune reported that during submaximal running the energy stored in the SEE of the quadriceps femoris MTU is about 15% of the energy stored in all SEEs of the musculoskeletal system together (Sasaki and Neptune, 2005). Further the total mechanical work done (negative plus positive) by the quadriceps femoris CEs is about 42% of the total mechanical work done by all other CEs of the human system (Sasaki and Neptune, 2005). From this, 42% of the total mechanical work of the quadriceps femoris CE, approximately 70% is negative mechanical work (Sasaki and Neptune, 2005).

The above reports demonstrate that the mechanical work done by the CE and the energy conversion in the SEE of the quadriceps femoris MTU can contribute considerably to the metabolic energy consumption during running. The higher strain at similar exerted tendon force results into a higher energy storage in the tendon of the quadriceps femoris MTU in group 1. Fig. 7 shows the energy storage capacity of the tendon and aponeurosis and the percentage differences between the low/moderate economy and the most economical runners as a function of the tendon force. The values are calculated by fitting a second order polynomial equation to the experimental data in Fig. 4. Concerning the force-strain relationship of the tendon and aponeurosis at the quadriceps femoris we found that the groups 2 and 3 showed similar values. Therefore we considered these two groups together. For comparison of the energy storage capacity of the tendon and aponeurosis (Fig. 7), we presented the values of the high running economy versus the moderate and low running economy runners. The percentage differences between the high economy group of runners and the other two groups in the energy storage capacity at medium and lower tendon forces (i.e. submaximal contractions) are higher than those at maximal forces. The higher ability of the most economical runners to store energy in their tendon and aponeurosis may increase the energy conversion in the SEE and at the same time decrease the negative and positive mechanical work done by the CE of the quadriceps femoris MTU. Therefore it is reasonable to assume that the ability of the most economical runners to achieve a higher energy conversion in the SEE could be an important issue causing the differences in running economy.

View larger version (21K): [in this window] [in a new window]   Fig. 7. Energy storage capacity of the tendon and aponeurosis as function of the tendon force between the runners of group 1 (high running economy) and the runners of groups 2 and 3 (moderate and low running economies), and percentage differences relative to group 1 in the quadriceps femoris muscle-tendon unit (as statistically significant differences in the force-strain relationship were found only between group 1 and both other groups, groups 2 and 3 have been considered together).

Higher contractile strength and tendon stiffness at the triceps surae
At the triceps surae MTU, group 1 shows a higher maximal contractile strength (about 36% higher ankle joint moments and calculated tendon forces) as well as a higher normalised stiffness relative to the other two groups. The higher plantarflexion moments of group 1 were observed at most of the analysed joint angle configurations (Fig. 5) indicating that the contributions of the different heads of the triceps surae to the total moment were similar for the three studied groups. The lack of differences in maximal ankle joint moments at short lengths of the triceps surae MTU (i.e. ankle plantarflexed and at the same time knee flexed positions) can be explained as follows. (1) Because of the parabolic curve of the force-length relationship of the triceps surae muscles (SOL, GM, GL) the group-related differences in plantarflexion moment would be reduced at short fascicle lengths. (2) The active insufficiency of the gastrocnemii at pronounced knee flexed positions, i.e. at short muscle lengths (Cresswell et al., 1995; Miaki et al., 1999) would also reduce the group-related differences in the joint moments. At flexed knee-joint positions the gastrocnemii reach a critical shortened length at which, due to the force-length relationship, the torque output cannot be increased and therefore the gastrocnemii decrease their activation level (Cresswell et al., 1995; Kenedy and Cresswell, 2001). In accordance with these findings the examined runners decrease the EMG-activity of the triceps surae muscles at the more pronounced plantarflexion and knee flexion angles (Fig. 5).

The higher stiffness of group 1 (most economical runners) does not favour the muscle force generation based on a reduction of the shortening velocity of the muscle fibres nor the energy storage and return in the tendon at a given tendon force. However, a higher force would contribute to a higher energy storage and return at identical stiffness values of the tendons. This means that the energy storage capacity is oppositely affected by these two factors. Therefore from the present study, it is difficult to isolate the influence of the higher stiffness and higher maximal contractile strength of the triceps surae MTU on differences in running economy observed between the examined groups. Hof et al. used an inverse dynamics approach to study the triceps surae MTU at submaximal velocities (3.63 and 3.93 m s^-1) comparable to those in the present study (Hof et al., 2002). They showed that during running the subjects adapted to their own tendon stiffness values. Although the range of ankle joint motion was similar, the maximal muscle forces at the triceps surae during running were higher for the subject with the stiffer tendon (Hof et al., 2002). Unfortunately we don't have any data describing the ankle joint moments during running from our subjects. Therefore we can only speculate about possible adjustments in joint kinetics. From the results of the present study (higher contractile strength and higher tendon stiffness at the triceps surae for the most economical runners) and the reports from the literature (Hof et al., 2002) it seems that the functionality of the MTU at submaximal running is not only dependent on the stiffness of the SEE but also on the maximal strength of the CE.

Limitations
We measured the oxygen consumption of the whole body during running at three different velocities. We examined the mechanical and morphological properties of the triceps surae and quadriceps femoris MTUs, however, and argued that the differences found in the mechanical properties of these two MTUs may be responsible for the differences in running economy between the examined runners. This means that we assumed that the triceps surae and quadriceps femoris muscles are the main contributors to the energy expenditure while running. Earlier studies (Winter, 1983; Arampatzis et al., 2000) analysing submaximal running on the basis of inverse dynamics reported that the main contributors to the total mechanical work during running are the muscles acting around the ankle and knee joints (>70%). Similar studies based on forward dynamic simulations of running using a musculoskeletal model (Sasaki and Neptune, 2005) found that the mechanical work done by the CEs of the triceps surae and quadriceps femoris muscles is about 68% of the mechanical work of all muscles CEs. In addition they reported (Sasaki and Neptune, 2005) that the energy storage in the SEEs of the triceps surae and quadriceps femoris MTUs is about 75% of the energy stored in the SEEs of the whole body. Based on all these studies, which rely on inverse dynamic analyses as well as on forward simulations, it is reasonable to assume that the triceps surae and quadriceps femoris MTUs may be representative of the energy expenditure of submaximal running.

In the present study we calculated tendon forces using tendon-moment arm data taken from the literature (Herzog and Read, 1993; Maganaris et al., 1998). It cannot be excluded that individual differences in the anatomical moment arms between groups could exist and could influence the calculated tendon forces. To estimate this potential source of error we analysed the ratio between tendon and aponeurosis displacement and ankle/knee joint angular rotation during a passive condition.

In order to do this we utilised the values from the correction of the tendon and aponeurosis elongation due to joint rotation (see Materials and methods). The average ratios were 0.55-0.69 mm deg.^-1 for the ankle joint and 0.42-0.54 mm deg.^-1 for the knee joint. There were no differences between groups in this ratios. Consequently it is likely that there were no differences between groups in the moment arms either. So a major influence of possible differences in lever arms between groups can be excluded.

Conclusions
All groups of runners (high, moderate and low running economy) showed similar kinematic characteristics during submaximal running, supporting earlier reports, which found that kinematic parameters cannot explain the complexity of running economy (Williams and Cavanagh, 1987; Martin and Morgan, 1992; Kyrölainen et al., 2001). The mechanical properties of the tendon and aponeurosis at the triceps surae and quadriceps femoris MTUs of the most economical runners showed clear differences from those of the moderate and low running economy groups. The quadriceps femoris of the most economical runners had a more compliant tendon at low force levels (< 45% of the MVC), whereas the triceps surae had a higher contractile strength and a higher tendon stiffness from 45% to 100% MVC. It is suggested that the more compliant quadriceps tendon and aponeurosis at low level forces will increase the force potential of the muscle while submaximal running and therefore would decrease the volume of active muscle at a given force generation. Further, we suggest that the efficiency of the triceps surae muscle contraction at submaximal running not only depends on the stiffness of the tendon and aponeurosis, but also on the maximal muscle strength.

List of symbols and abbreviations

CE

contractile element

EMG

electromyography

GL

gastrocnemius lateralis

GM

gastrocnemius medialis

MTU

muscle-tendon unit

MVC

maximal voluntary contraction

PCSA

physiological cross-sectional area

QF

quadriceps femoris

RF

rectus femoris

RMS

root mean square

SEE

series elastic element

SOL

soleus

TS

triceps surae

VL

vastus lateralis

VM

vastus medialis

rate of oxygen consumption

	References

TOP Summary Introduction Materials and methods Results Discussion References

 

Abe, T., Kumagai, K. and Brechue, W. F. (2000). Fascicle length of leg muscles is greater in sprinters than distance runners. Med. Sci. Sports Exerc. 32,1125 -1129.

Abe, T., Fukashiro, S., Harada, Y. and Kawamoto, K. (2001). Relationship between sprint performance and muscle fascicle length in female sprinters. J. Physiol. Anthropol. Appl. Human Sci. 20,141 -147.[CrossRef][Medline]

Alexander, R. McN. and Bennet-Clark, H. C. (1977). Storage of elastic strain energy in muscle and other tissues. Nature 265,114 -117.[CrossRef][Medline]

Arampatzis, A., Knicker, A., Metzler, V. and Brüggemann, G.-P. (2000). Mechanical power in running: a comparison of different approaches. J. Biomech. 33,457 -463.[CrossRef][Medline]

Arampatzis, A., Karamanidis, K., De Monte, G., Stafilidis, S., Morey-Klapsing, G. and Brüggemann, G.-P. (2004). Differences between measured and resultant joint moments during voluntary and artificially elicited isometric knee extension contractions. Clin. Biomech. 19,277 -283.[CrossRef][Medline]

Arampatzis, A., Stafilidis, S., De Monte, G., Karamanidis, K., Morey-Klapsing, G. and Brüggemann, G.-P. (2005a). Strain and elongation of the gastrocnemius tendon and aponeurosis during maximal plantarflexion effort. J. Biomech. 38,833 -841.[CrossRef][Medline]

Arampatzis, A., Morey-Klapsing, G., Karamanidis, K., De Monte, G., Stafilidis, S. and Brüggemann, G.-P. (2005b). Differences between measured and resultant joint moments during isometric contractions at the ankle joint. J. Biomech. 38,885 -892.[CrossRef][Medline]

Baratta, R., Solomonow, M., Zhou, B., Letson, D. and D'Ambrosia, R. (1988). The role of antagonistic musculature in maintaining knee stability. Am. J. Sports Medic. 16,113 -122.

Biewener, A. A. (1989). Scaling body support in mammals: limb posture and muscle mechanics. Science 245, 45-48.[Abstract/Free Full Text]

Biewener, A. A. (1990). Biomechanics of mamalian terrestrial locomotion. Science 250,1097 -1103.[Abstract/Free Full Text]

Biewener, A. A. and Baudinette, R. V. (1995). In vivo muscle force and elastic energy storage during steady-speed hopping of tammar wallabies (Macropus eugenii). J. Exp. Biol. 198,1829 -1841.

Biewener, A. A. and Gillis, G. B. (1999). Dynamics of muscle function during locomotion: accommodating variable conditions. J. Exp. Biol. 202,3387 -3396.[Abstract]

Biewener, A. A. and Roberts, T. J. (2000). Muscle and tendon contributions to force, work, and elastic energy savings: a comparative perspective. Exerc. Sports Sci. Rev. 28, 99-107.

Biewener, A. A., Konieczynski, D. D. and Baudinette, R. V. (1998). In vivo muscle force-length behavior during steady-speed hopping in tammar wallabies. J. Exp. Biol. 201,1681 -1694.[Abstract]

Biewener, A. A., Farley, C. T., Roberts, T. J. and Temaner, M. (2004). Muscle mechanical advantage of human walking and running: implications for energy cost. J. Appl. Physiol. 97,2266 -2274.[Abstract/Free Full Text]

Bijker, K. E., De Groot, G. and Hollander, A. P. (2001). Delta efficiencies of running and cycling. Med. Sci. Sports Exerc. 33,1546 -1551.[Medline]

Bobbert, M. F. (2001). Dependence of human squat jump performance on the series elastic compliance of the triceps surae: a simulation study. J. Exp. Biol. 204,533 -542.[Abstract]

Bojsen-Moller, J., Hansen, P., Aagaard, P. and Magnusson, S. P. (2003). Measuring mechanical properties of the vastus lateralis tendon-aponeurosis complex in vivo by ultrasound imaging. Scand. J. Med. Sci. Sports 13,259 -265.[CrossRef][Medline]

Cavanagh, P. R. and Kram, R. (1985). The efficiency of human movement - a statement of the problem. Med. Sci. Sports Exerc. 17,304 -308.[Medline]

Cavanagh, P. R. and Williams, K. R. (1982). The effect of stride length variation on oxygen uptake during distance running. Med. Sci. Sports Exerc. 14, 30-35.[Medline]

Cresswell, A. G., Löscher, W. N. and Thorstensson, A. (1995). Influence of gastrocnemius muscle length on triceps surae torque development and electromyographic activity in man. Exp. Brain Res. 105,283 -290.[Medline]

Crow, M. T. and Kushmerick, M. J. (1982). Chemical energetics of slow- and fast-twitch muscles of the mouse. J. Gen. Physiol. 79,147 -166.[Abstract/Free Full Text]

Daniels, J. and Daniels, N. (1992). Running economy of elite male and elite female runners. Med. Sci. Sports Exerc. 24,483 -489.[Medline]

Daniels, J. T., Yarbrough, R. A. and Foster, C. (1978). Changes in _Omax and running performance with training. Eur. J. Physiol. 39,249 -254.

de Haan, A., van Ingen Schenau, G. J., Ettema, G. J., Huijing, P. A. and Lodder, A. N. (1989). Efficiency of rat medial gastrocnemius muscle in contractions with and without an active prestretch. J. Exp. Biol. 141,327 -341.[Abstract/Free Full Text]

Dolgener, F. (1982). Oxygen cost of walking and running in untrained, sprint trained, and endurance trained females. J. Sports Med. Phys. Fitness 22, 60-65.[Medline]

Ettema, G. J. C. (1996). Mechanical efficiency and efficiency of storage and release of series elastic energy in skeletal muscle during stretch-shorten cycles. J. Exp. Biol. 199,1983 -1997.[Abstract]

Ettema, G. J. C., Huijing, P. A., van Ingen Schenau, G. J. and de Haan, A. (1990a). Effects of prestretch at the onset of stimulation on mechanical work output of rat medial gastrocnemius muscle-tendon complex. J. Exp. Biol. 152,333 -351.[Abstract/Free Full Text]

Ettema, G. J. C., van Soest, A. J. and Huijing, P. A. (1990b). The role of series elastic structures in prestretch-induced work enhancement during isotonic and isokinetic contractions. J. Exp. Biol. 154,121 -136.[Abstract/Free Full Text]

Herzog, W. and Read, L. J. (1993). Lines of action and moment arms of the major orce-carrying structures crossing the human knee joint. J. Anat. 182,213 -230.

Herzog, W., Abrahamse, S. K. and ter Keurs, H. E. D. J. (1990). Theoretical determination of force-length relations of intact human skeletal muscles using the cross-bridge model. Pflügers Arch. Eur. J. Physiol. 416,113 -119.[CrossRef][Medline]

Hof, A. L., Geelen, B. A. and van den Berg, J. W. (1983). Calf muscle moment, work and efficiency in level walking; role of series elasticity. J. Biomech. 16,523 -537.[CrossRef][Medline]

Hof, A. L., van Zandwijk, J. P. and Bobbert, M. F. (2002). Mechanics of human triceps surae muscle in walking, running and jumping. Acta Physiol. Scand. 174, 17-30.[CrossRef][Medline]

Karamanidis, K., Arampatzis, A. and Brüggemann, G.-P. (2003). Symmetry and reproducilility of kinematic parameters during various running techniques. Med. Sci. Sports Exerc. 35,1009 -1016.[CrossRef][Medline]

Kawakami, Y., Abe, T. and Fukunaga, T. (1993). Muscle-fiber pennation angles are greater in hypertrophied than in normal muscles. J. Appl. Physiol. 74,2740 -2744.[Abstract/Free Full Text]

Kawakami, Y., Abe, T., Kuno, S. Y. and Fukunaga, T. (1995). Training-induced changes in muscle architecture and specific tension. Eur. J. Appl. Physiol. 72, 37-43.[CrossRef]

Kennedy, P. M. and Cresswell, A. G. (2001). The effect of muscle length on motor-unit recruitment during isometric plantar flexion in humans. Exp. Brain Res. 137, 58-64.[CrossRef][Medline]

Ker, R. F., Bennet, M. B., Bibby, S. R., Kester, R. C. and Alexander, R. McN. (1987). The spring in the arch of the human foot. Nature 325,147 -149.[CrossRef][Medline]

Kumagai, K., Abe, T., Brechue, W. F., Ryushi, T., Takano, S. and Mizuno, M. (2000). Sprint performance is related to muscle fascicle length in male 100-m sprinters. J. Appl. Physiol. 88,811 -816.[Abstract/Free Full Text]

Kyröläinen, H., Belli, A. and Komi, P. V. (2001). Biomechanical factors affecting running economy. Med. Sci. Sports Exerc. 8,1330 -1337.

Lake, M. and Cavanagh, P. R. (1996). Six weeks of training does not change running mechanics or improve running economy. Med. Sci. Sports Exerc. 28,860 -869.

Mademli, L., Arampatzis, A., Morey-Klapsing, G. and Brügemann, G.-P. (2004). Effect of ankle joint position and electrode placement on the estimation of the antagonistic moment during maximal plantarflexion. J. Electromyogr. Kinesiol. 14,591 -597.[CrossRef][Medline]

Maganaris, C. N., Baltzopoulos, V. and Sargeant, A. (1998). Changes in Achilles tendon moment arm from rest to maximum isometric plantarflexion: in vivo observations in man. J. Physiol. 510,977 -985.[Abstract/Free Full Text]

Magnusson, S. P., Aagaard, P., Rosager, S., Dyhre-Poulsen, P. and Kjaer, M. (2001). Load-displacement properties of the human triceps surae aponeurosis in vivo. J. Physiol. 531,277 -288.[Abstract/Free Full Text]

Martin, P. E. and Morgan, D. W. (1992). Biomechanical considerations for economical walking and running. Med. Sci. Sports Exerc. 24,467 -474.

McMahon, T. A. (1985). The role of compliance in mammalian running gaits. J. Exp. Biol. 115,263 -282.[Abstract/Free&