Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Anatomically Asymmetrical Runners Move More Asymmetrically at the Same Metabolic Cost

  • Elena Seminati ,

    elena.seminati@unimi.it

    Affiliation Department of Pathophysiology and Transplantation, Faculty of Medicine, University of Milan, Milan, Italy

  • Francesca Nardello,

    Affiliation Department of Neurological and Movement Sciences, School of Exercise and Sport Sciences, University of Verona, Verona, Italy

  • Paola Zamparo,

    Affiliation Department of Neurological and Movement Sciences, School of Exercise and Sport Sciences, University of Verona, Verona, Italy

  • Luca P. Ardigò,

    Affiliation Department of Neurological and Movement Sciences, School of Exercise and Sport Sciences, University of Verona, Verona, Italy

  • Niccolò Faccioli,

    Affiliation Department of Pathology and Diagnostics, Section of Radiology, University of Verona, Verona, Italy

  • Alberto E. Minetti

    Affiliation Department of Pathophysiology and Transplantation, Faculty of Medicine, University of Milan, Milan, Italy

Abstract

We hypothesized that, as occurring in cars, body structural asymmetries could generate asymmetry in the kinematics/dynamics of locomotion, ending up in a higher metabolic cost of transport, i.e. more ‘fuel’ needed to travel a given distance. Previous studies found the asymmetries in horses’ body negatively correlated with galloping performance. In this investigation, we analyzed anatomical differences between the left and right lower limbs as a whole by performing 3D cross-correlation of Magnetic Resonance Images of 19 male runners, clustered as Untrained Runners, Occasional Runners and Skilled Runners. Running kinematics of their body centre of mass were obtained from the body segments coordinates measured by a 3D motion capture system at incremental running velocities on a treadmill. A recent mathematical procedure quantified the asymmetry of the body centre of mass trajectory between the left and right steps. During the same sessions, runners’ metabolic consumption was measured and the cost of transport was calculated. No correlations were found between anatomical/kinematic variables and the metabolic cost of transport, regardless of the training experience. However, anatomical symmetry significant correlated to the kinematic symmetry, and the most trained subjects showed the highest level of kinematic symmetry during running. Results suggest that despite the significant effects of anatomical asymmetry on kinematics, either those changes are too small to affect economy or some plastic compensation in the locomotor system mitigates the hypothesized change in energy expenditure of running.

Introduction

The symmetry between the left and right sides of the body plays an important role in legged locomotion. The symmetrical behaviour of lower limbs during gait has often been taken for granted, mainly for simplicity in data collection and analysis, while the lack of it was frequently considered as an indicator of gait pathology [1]. Differently from what expected, healthy human gait is rather asymmetrical [2], [3]. This seems to reflect a functional difference inherently associated to the laterality of the dominant side characterising each individual [4], [5]. This topic was introduced more than 80 years ago by Lund [6] who showed the effects of structural/anatomical asymmetry on lateral drift in human locomotion. The same experiments were recently repeated and supported the hypothesis of a relationship between leg length inequality and asymmetry in locomotion [7][9].

Body symmetry can be further modulated in sports: depending on the discipline, relevant muscles become asymmetrically different (tennis, fencing, throwing, etc.), or they are required to reach similar hypertrophy (ice-skating, downhill skiing, front crawl, etc.) on the two sides of the sagittal plane. Thus, body changes towards or from symmetry are not just the consequence of genetics and laterality, being also caused by specific training protocols.

As the concept of symmetry has an important influence in human locomotion, it plays a key role in the design and maintenance of vehicles, which are periodically inspected and serviced to guarantee wheel balance and homogeneous tyre wearing, in order to reduce fuel consumption and ensure a safe drive. Would it be the same for human running? Can an anatomical/structural asymmetry of the human body cause kinematic/dynamic asymmetry of locomotion? Also, can structural or functional asymmetries be related to some increase of the metabolic cost of transport?

Several authors studied symmetry in locomotion in humans [1][4], [10], [11] and also in animals [12], but only few of them investigated the possible interaction between symmetry and energy saving. Manning and collaborators found negative correlations between anatomical symmetry and race time during competitions, both in human running and in galloping horses [13], [14]. These preliminary findings encouraged us to study the possible interactions between different kinds of symmetry (anatomical and dynamical) and the human running performance, not only in term of race time, but also of energy saving. In the present study, we investigate the relationship between the cost of transport (C) while running at different increasing velocities and individual anatomical and dynamical symmetries in three differently trained groups of subjects, with the idea that ‘race cars’ should more strongly rely on symmetry than ordinary ‘automobiles’.

Materials and Methods

Subjects

Nineteen healthy male subjects volunteered to participate this investigation. Exclusion criteria included neurological or musculoskeletal pathologies affecting running ability. The institutional ethics committee of the University of Milano had approved all methods and procedures, and subjects gave their written informed consent (approved by the same committee) prior to the start of testing. We clustered participants into three different groups, based on their specific running ability:

  • group 1, (n = 7): Untrained Runners (UR), who practiced sport (not specifically running) 3 times per week (less than 2 hours per week)
  • group 2, (n = 7): Occasional Runners (OR), fit athletes, who trained more than 3 times per week, (between 2 and 6 hours per week). Each of them had previously participated in a national competition (half marathon or 10 km competition)
  • group 3, (n = 5): Skilled Runners (SR), master athletes who trained more than 3 times per week (at least 6 hours per week); they were marathon runners, with a mean performance time of 2 h 44 min 24 s ±10 min 12 s standard deviation (SD).

Anthropometric characteristic of the different subject groups are shown in Table 1.

MR Dataset and 3D Images Processing

In order to evaluate the anatomical symmetries, each participant underwent Magnetic Resonance (MR) imaging. Subjects were adjusted in a supine position as to preserve the maximal body symmetry in the sagittal plane.

MR scans were performed with a 1.5-T superconductive magnet (Siemens, Erlangen, Germany). In all subjects multiplanar T1-weighted Spin-echo sequences were obtained (TE 11, TR 565, flip angle 90°), on a coronal plane for three different anatomical districts: Pelvis district (PD), Upper-Leg district (UD), including thigh and knee, Lower-Leg district (LD), including calf and ankle, with slice thickness of 4 mm. The matrix was 320×320 and the field of view (FOV) was 460×460. Total examination time was less than 7 minutes (36 coronal slices for each district).

All the recorded images (saved in DICOM format) were subsequently analyzed with a custom, ad hoc program written in LabVIEW 8.6 (National Instrument, Austin, Texas, USA). The procedure we implemented exports, for each districts, 36 MR images (slices) as two-dimensional matrix of 320×320 pixels, each of which 1.44×1.44 mm, and includes several post-processing steps, as shown in Figure 1.

thumbnail
Figure 1. Principal steps involved in the 3D cross-correlation algorithm.

a) The 36 slices of the MR sequence, create a 3D volume whose sizes are laterally indicated, b) right volume (Rv) and left volumes (Lv) separated, c) left reflected volume (Lrv) on the sagittal plane (in the mirror), d) zero-padding operation around right volume, e) Lrv superimposed to Rv in order to find the position that maximize the cross-correlation value.

https://doi.org/10.1371/journal.pone.0074134.g001

The 36 coronal slices, for every district, assembled together, re-create a three-dimensional (3D) volume, whose elements (voxel) are values corresponding to a grey level intensity (8 bit scale), reflecting proton density, (Figure 1a). In order to compare the subject’s left lower limb with the right one, firstly, the initial 3D volume has to be split in two separated volumes, right volume (Rv) and left volume (Lv), (Figure 1b). Successively the Lv is specularly reflected, with respect to the sagittal plane (Figure 1c), whilst the Rv is bordered by zero intensity voxel (Figure 1d), through a zero-padding operation, so that the left reflected volume (Lrv) can be virtually superimposed on the Rv (Figure 1e), and moved along the three axes in order to find the best matching overlap and to evaluate the ‘overall’ similarity (i.e. symmetry) between the two limbs. To achieve this aim the algorithm performs a 3D correlation between the contents of the two respective anatomical volumes.

The correlation between two signals (cross-correlation) is a standard approach for signal processing and it has been recently designed in 3D in order to consider simultaneously the full anatomical volume information, to assist radiologists in providing correct diagnosis of metastases within the lungs [15], [16] or brain [17], for instance.

Following Lewis’ approach [18], a normalised cross-correlation coefficient (ri,j,k), was adopted to identify the symmetry degree between the 3D split volumes:where and are the voxel mean value of the left reflected volume and the right volume, respectively. The two volumes are virtually superimposed at coordinates i, j and k, and calculations are performed for all pairs of corresponding voxels along x, y and z axes.

For every subject and each anatomical district we evaluated the maximal cross correlation value (rmax) (i.e. the value corresponding to the best overlap between right and left reflected volumes). This coefficient can assume a range of values between −1 and 1, depending upon the similarity of the 3D analyzed volumes, where a value of 1 indicates an exact matching of the Lrv with the Rv, a value of −1 indicates opposite grey values for voxels in Lrv with respect to Rv, and a value of 0 indicates no correlation between the two volumes.

Software reliability and accuracy were validated by comparing two identical bottles filled up with water (see Figure 2a), resulting in a maximal cross correlation value of rmax = 0.99. The algorithm provided a value of rmax = 1 only when the right volume of a specific subject was compared with itself, while the lowest value of rmax was obtained when the right volume of a specific subject was compared with the left reflected volume of an other different subject (rmax = 0.51, as shown in Figure 2b).

thumbnail
Figure 2. Examples of obtained cross-correlation values plotted versus iterations.

Cross correlation values (r) of all iterations (134,400 overlap positions = 28 (i)×60 (j)×80 (k), between right volume and left reflected volume); a) comparison between two bottles, filled with the same volume of water (rmax = 0.99); b) comparison between the right upper leg of a subject and the left upper leg of a different subject (rmax = 0.51); c) comparison between right and left upper legs of the same subject. Inset: enlargement of cross correlation pattern showing the inner processing loop (i coordinates).

https://doi.org/10.1371/journal.pone.0074134.g002

We evaluated for every subject a single maximal cross correlation value (rmax) for each district, (rmax(PD) for Pelvis district, rmax(UD) for Upper-Leg district and rmax (LD) for Lower-Leg district) and secondly a ‘global’ anatomical cross correlation value () as the mean of the three districts:

Kinematics

In order to capture kinematic functional symmetries on many steps, all the subjects performed level shod running on a treadmill (h/p/Cosmos Saturn 4.0, Germany).

Human body has been modelled as a series of linked, rigid segments: 18 reflective markers were placed bilaterally on anatomical landmark points (immediately anterior to ear tragus, shoulder, elbow, wrist, greater trochanter, lateral epicondyle of femur, lateral malleolus, calcaneus, and 5th metatarsal head) and their 3D position was captured at 100 Hz, using an eight-camera Vicon MX optoelectronic system (Vicon, Oxford, UK). In this way, 12 body segments were defined [19].

After a brief period of familiarization on the treadmill, each subject ran at six different incremental speeds: from 2.22 m/s to 5 m/s, step 0.56 m/s. Each speed was maintained for at least 5 min, with a rest period of at least 5 min between successive trials.

The 3D recorded coordinates of the 12 segments, together with the anthropometric tables [20], [21], were used to compute the experimental trajectory of the Body Centre of Mass (BCOM). Successively, we adopted a recent mathematical method [22], [23] simultaneously capturing the spatial and dynamical features of that 3D BCoM trajectory, which allows to quantify dynamical symmetry indices of locomotion in the 3 spatial axes; by having sampled the body motion on a treadmill, the trajectory of the BCOM can be represented by a closed 3D loops (Lissajous contours), representing its displacement with respect to the average position. The 3D trajectory is mathematically defined by a 6-harmonic Fourier series, whose coefficients are used to calculate the Dynamical Symmetry Indices (for progression axis), (for vertical axis) and (for lateral axis). The motion of the BCOM is expected to exhibit perfect right–left symmetry if it contained just even harmonics in the progression and y vertical directions, and just odd harmonics in the lateral direction, as within a stride it oscillates twice in the sagittal (y–x) plane and only once in the horizontal (x–z) plane. Dynamical Symmetry indices are then averaged among the strides number (n) as to obtain for each velocity and each subject:

(SI, 0: no symmetry between right and left steps, 1: complete symmetry).

Successively, the three mean dynamic indices (, , ) are weighted according to the ‘real’ maximum displacement range of the BCOM, i.e. dx ( = running speedXstride frequency), dy and dz, respectively, and a Global symmetry Index (GI) is calculated as

(GI, 0: no symmetry between right and left steps, 1: complete symmetry).

Energy Cost Measurement

Oxygen consumption () of running was measured with a breath-by-breath gas analyzer (Cosmed K4b2, Rome, Italy). Data, including heart rate (HR), were recorded at each progression speed, after the metabolic steady state had been achieved (3 min), for further 2 minutes. 5 minutes of testing was performed at each speed. Resting was measured while standing. Respiratory Exchange Ratio (RER) was monitored in order to check for aerobic conditions (RER<1). We expressed the metabolic Cost of Transport (C), i.e. the oxygen consumed to move 1 kg of body mass 1 m distance, in J (kg m)−1 by dividing the net [measured - resting, [ml O2 (Kg min)−1] by the progression speed (m min−1), and by assuming an energy equivalent of 20.9 J ml O2−1.

Statistical Analysis

Relationships between variable pairs were investigated using Pearson’s correlation coefficient. To compare speed dependent variables (C, HR, , , and GI), differences were analyzed using a two-ways ANOVA (groupxspeed) (with a post-hoc Bonferroni correction). For speed independent variables (rmax(PD), rmax(UD) and rmax(LD)), we performed a one-way ANOVA for repeated measures in order to detect difference among districts. Furthermore, Principal Component Analysis (PCA) was performed on the three anatomical indices, in order to estimate their relative contribution to the total variance.

Statistical significance was accepted when p<0.05.

Results

Since only five OR and five SR subjects were able to complete all the running protocols up to 5.0 m/s, and UR subjects stopped at the speed of 4.44 m/s, we did not consider in the statistical analysis the highest speed level.

Anatomical Symmetries

Anatomical symmetries are described by the maximal cross-correlation value for each district (rmax(PD), rmax(UD) and rmax(LD)), and by the global anatomical cross correlation value (). These values are limited to only 17 subjects, because two MR tests (one for the UR and one for the SR) had to be discarded due to technical problems.

One-way ANOVA between the three groups of subject didn’t show any difference between UR, OR and SR for the cross-correlation values, while we found significantly lower values of anatomical symmetry for pelvis district, compared to the upper (p<0.05) and lower leg district (p<0.01) (rmax(PD) = 0.77±0.09, rmax(UD) = 0.82±0.05 and rmax(LD) = 0.83±0.05). PCA showed that 65.8% of the total variance was explained by the first principal component, where the three considered parameters (rmax(LD), rmax(UD) and rmax(PD)) had almost the same weight. However UD seems to give the greatest contribution to the first principal component, with respect to the other two districts.

Results regarding pairwise correlations between variables are summarized in Table 2: rmax(UD) is significantly correlated with rmax(PD) and rmax(LD) (p<0.05), also rmax(PD) and rmax(LD) seem to be positively correlated even if not significantly (p = 0.087). Significant results were found also between anatomical symmetries and kinematics (mean values for the Global Symmetry Index () were evaluated starting from the single values of , and and averaged within each group of speeds for each subject); in particular rmax(PD) is positively correlated with (p<0.01) and (p<0.05) and also is positively and significantly correlated with (p<0.05), while we observed a positive trend between and GI, even if not significantly (p = 0.055) (see Figure 3).

thumbnail
Figure 3. Regression of the mean dynamic Global Symmetry Index (GI) versus the global anatomical cross correlation value ().

Each point represents the mean Global Symmetry Index averaged among the different running speeds for each subject; Untrained Runners (UR), Occasional Runners (OR) and Skilled Runners (SR) (r = 0.473; p = 0.055).

https://doi.org/10.1371/journal.pone.0074134.g003

thumbnail
Table 2. Statistical correlation matrix results between variable pairs.

https://doi.org/10.1371/journal.pone.0074134.t002

Kinematics

Mean values for the Global Symmetry Index (), evaluated starting from the single values of , and and averaged within each group of subjects, ± SD, are shown in Figure 4. We performed a two-ways ANOVA, where independent variables were running speed and subject group and the dependent variable was GI. Results show that UR have a GI significantly lower than both OR and SR at each velocity (p<0.01). Also, GI for UR seems to decrease with increasing running velocity, even if not significantly. Statistical analysis did not show any difference between UR, OR and SR for the single kinematic symmetry indices, while one-way ANOVA for repeated measure shown significantly lower values for , (0.72±0.06) compared to (0.88±0.04) and (0.85±0.05), (p<0.01).

thumbnail
Figure 4. Mean values for the dynamic Global Symmetry Index () plotted against running speed.

Mean values for the dynamic Global Symmetry Index () are evaluated starting from the single values of , and and averaged within each group of subjects, ± SD, in untrained runners (UR), occasional runners (OR) and skilled runners (SR). Two-way ANOVA (group×running speed) show that the group of UR had a mean GI always lower compared to the OR and SR, (* = p<0.01), independently from the running speed.

https://doi.org/10.1371/journal.pone.0074134.g004

Cost of Transport

Results for the metabolic cost C and HR are presented in Figure 5. C is confirmed to be independent of speed, with no differences among running groups. At the same speed, HR decreased as runners’ ability increased, with values for SR significantly lower than for OR and UR. No significant correlation was found between the C and the previously analysed parameters, both for kinematics and for anatomical values (Table 2).

thumbnail
Figure 5. Mean values ± SD for the cost of transport (C) (lower curves) and for the heart rate (HR) (upper curves).

C and HR are plotted against running speed for untrained runners (UR), occasional runners (OR) and skilled runners (SR). Results obtained with the two-way ANOVA (group×running speed) show no significant difference among groups of subjects across velocity for C, which results to be independent of the running speed. HR increased significantly with the running speed for all the three group of subjects. Furthermore we obtained significantly higher HR values for UR compared to OR and SR (* = p<0.01).

https://doi.org/10.1371/journal.pone.0074134.g005

Discussion

The main aim of this project was to investigate the relationship among the anatomical/structural symmetry of the lower limbs, the dynamical symmetry of the 3D BCOM displacement and the metabolic cost of human running. C has been considered as an indirect index of running performance: at the same sustainable fraction of maximal , the lower the cost the higher the average speed [24]. While being aware of the speed and training level independency of C, as debated and reported in the literature [25][29], our hypothesis was that more asymmetrical limbs, in subjects committed to run with symmetrical steps, would have involved a higher C. In other words, part of the inter-subject C variance could have been explained by different level of anatomical asymmetry.

Differently from previous studies dealing with gross morphological features (bones length [5], [8], human face [13] and horse muzzle [14] landmarks) and isolated gait parameters (stride length and frequency [30], [31], joint angles [3] and ground reaction forces [10]), we analysed the symmetry of the ‘whole’ (left and right) lower limb anatomy and of the global running kinematics (3D trajectory of BCOM), in three groups of differently trained athletes. Our hypothesis, inspired by the engineering of motor vehicles, was not completely verified. C was not significantly correlated either with anatomical symmetries or with dynamical symmetries in running, while we found significant correlations between the anatomical and dynamical symmetries indices (Table 2). This indicates that the more anatomically symmetrical are the subjects, the more symmetrical is their running gait (especially in the forward (x) direction).

This finding is in accordance with the recent literature, where high level of leg length discrepancy (LLD) is correlated with low symmetrical gait coefficients [7] in walking. In our work, individual LLD was always lower than 2 cm (Table 1), and had no effect on C, according to the studies of Gurney [32].

It is possible that some physiological adaptations of the human machinery compensate for small asymmetries typical of the mechanics of our legged system [1], [2], with no influence on C. Rather, larger anatomical discrepancies, like a LLD higher than 2 cm [32] or a body mass not uniformly distributed [33], [11], could influence economy. Similar adaptations behaviours might have occurred in runners wearing new and worn shoes [34], or on surfaces of different stiffness [35]. Despite of the changed properties of materials, runners modified their motion pattern as to retain their original dynamics of running.

This could occur also in the subjects of this study, who seem to compensate their anatomical body asymmetries and minimize C, a strategy frequently adopted by animals [36]. With the main propulsive muscles operating close to isometric in running [37], tendons can store (stretching) and release (shortening) variable amounts of elastic energy during each step, in the attempt to adapt to different anatomical asymmetries. In this way the metabolic cost can be potentially kept unchanged.

In addition, although HR results (Figure 5) witness the appropriateness of clustering subjects according to the different training status (most skilled runners reported the lowest HR, at the same speed, p<0.01), the almost speed-independent C values seem not to be influenced by the different fitness level, as also found by other investigators [27], [29].

As also indicated in previous studies, training and experience seem to be important elements in the lower limb joint angle symmetry and in the stride variability of running, even at no apparent metabolic benefit [29][31]. The most experienced and high performance athletes can maintain, even at high velocities, higher dynamical symmetry than untrained runners (Figure 4). As step frequency and muscles effort increase, the higher physical demand and peripheral fatigue could impair the maintenance of a symmetrical gait and a consistent locomotion pattern, as seen for the UR group.

Furthermore, MRI measurements showed that the anatomical symmetry does not depend on the investigated district. PCA and correlation among lower limb districts could have been caused by misalignments of the two limbs during MRI test. However, due to the use of alignment tools during the tests, we feel confident that the intra-subject symmetry correlation among districts is not a measurement artefact. Similar eigenvalues from PCA suggest that the total variance of symmetry is equally explained by the three districts.

This work brings developments in the study of locomotion symmetry, also by means of newly introduced methodologies (BCOM 3D trajectory analysis and 3D cross-correlation between ‘whole’ limb MRI voxels). Differently from the original hypothesis, asymmetrical limbs generate asymmetrical body running at no apparent additional metabolic cost. This suggests some plasticity of the human body in coping with structural changes, with the final result of preserving locomotion economy. Deeper insights have been obtained regarding the relationship between the symmetries correlation residuals and the cost of transport, with the idea that subjects would be less economic when their anatomical and dynamical symmetry values do not match. Supplemental analysis and discussion regarding this hypothesis have been reported in the Appendix S1. Even if statistical results in this perspective are weak, possibly due to the relatively small sample size and low asymmetry level, there are some hints suggesting that only the runners who fail to match their anatomy and dynamics features have an increased cost of locomotion. Therefore, the initial hypothesis embedded in the title “anatomically asymmetrical runners move more asymmetrically at the same metabolic cost” is still valid (i.e. the cost would increase when an anatomically asymmetrical runner attempts to move in a symmetrical way). Further studies focusing on adaptations of the muscle-tendon interplay could reveal how human machine compensate the small structural asymmetries that characterize our legged system. The anatomical asymmetry threshold, above which the now expected asymmetrical gait will also involve an increase in running cost, is the challenge for future investigations.

Acknowledgments

The authors would like to thank all the subjects for their participation in the study, and the Technician Lauro Dalla Chiara for his help during the MR scans performed at the University Hospital Polyclinic “Borgo Roma” in Verona (Italy). Statistical support from Carlo M. Biancardi is also greatly appreciated.

Supporting Information

Figure S1.

Examples of univariate and bivariate regressions. Four different types of linear regressions are presented as examples of correlation between Dynamical Symmetry index in forward direction () and maximal cross-correlation value for Pelvis District (): a) Univariate regression, b) Univariate regression with intercept forced to be equal to 0, c) Bivariate regression, d) Bivariate regression with intercept forced to be equal to 0. Untrained Runners (UR), Occasional Runners (OR) and Skilled Runners (SR) symbols as in Figure 3. N.B. The determination coefficient in regressions lines forced through the origin, differently from the general model, does not reflect the fraction of the variability in the dependent variable explained by the independent variable. This makes R2 values unrealistically high and not comparable with the ones obtained in the general models.

https://doi.org/10.1371/journal.pone.0074134.s001

(TIF)

Author Contributions

Conceived and designed the experiments: AEM LPA PZ. Performed the experiments: ES FN NF LPA PZ. Analyzed the data: ES FN AEM. Contributed reagents/materials/analysis tools: NF. Wrote the paper: ES AEM. Designed the software used in analysis: AEM ES. Final approval of the paper version to be published: ES FN PZ LPA NF AEM.

References

  1. 1. Sadeghi H, Allard P, Prince F, Labelle H (2000) Symmetry and limb dominance in able-bodied gait: a review. Gait & Posture 12: 34–45.
  2. 2. Nardello F, Ardigò LP, Minetti AE (2009) Human locomotion: Right/left symmetry in 3D trajectory of body centre of mass. Gait & Posture 30(S): S80–S81.
  3. 3. Forczek W, Staszkiewicz R (2012) An evaluation of symmetry in the lower limb joint during the able-bodied gait of women and men. Journal of Human Kinetics 35: 47–57.
  4. 4. Maupas E, Paysant J, Martinet N, Andre J (1999) Asymmetric leg activity in healthy subjects during walking, detected by electrogoniometry. Clin Biomech (Bristol, Avon) 14: 403–411.
  5. 5. Cuk T, Leben-Seljak P, Stefancic M (2001) Lateral asymmetry of human long bones. Variability and Evolution 9: 19–32.
  6. 6. Lund FH (1930) Physical asymmetries and disorientation. The American journal of Psychology 42(1): 51–62.
  7. 7. Gurney B, Mermier C, Robergs R, Gibson A, Rivero D (2001) Effects of limb-length discrepancy on gait economy and lower-extremity muscle activity in older adults. J Bone Joint Surg Am 83-A: 907–915.
  8. 8. Seeley MK, Umberger BR, Clasey JL, Shapiro R (2010) The relation between mild leg-lenght inequality and able-bodied gait asymmetry. J Sports Sci Med 9(4): 572–579.
  9. 9. Souman JL, Frissen I, Sreenivasa MN, Ernst MO (2009) Walking straight into circles Curr Biol. 19: 1538–1542.
  10. 10. Herzog W, Nigg BM, Read LJ, Olsson E (1989) Asymmetries in ground reaction force patterns in normal human gait. Med Sci Sports Exerc 21: 110–114.
  11. 11. Mattes SJ, Martin PE, Royer TD (2000) Walking symmetry and energy cost in persons with unilateral transtibial amputations: matching prosthetic and intact limb inertial properties. Arch Phys Med Rehabil 81: 561–568.
  12. 12. Halling Thomsen M, Tolver Jensen A, Sorensen H, Lindegaard C, Haubro Andersen P (2010) Symmetry indices based on accelerometric data in trotting horses. J Biomech 43: 2608–2612.
  13. 13. Manning JT, Pickup LJ (1998) Symmetry and performance in middle distance runners. Int J Sports Med 19: 205–209.
  14. 14. Manning J, Ockenden L (1994) Fluctuating asymmetry in racehorses. Nature 370: 185–186.
  15. 15. Lee Y, Hara T, Fujita H, Itoh S, Ishigaki T (2001) Automated detection of pulmonary nodules in helical CT images based on an improved template-matching technique. IEEE Trans Med Imaging 20: 595–604.
  16. 16. Wang P, DeNunzio A, Okunieff P, O′Dell WG (2007) Lung metastases detection in CT images using 3D template matching. Med Phys 34: 915–922.
  17. 17. Ambrosini R, Wang P, O′Dell W (2010) Computer-aided detection of metastatic brain tumors using automated three-dimensional template matching. J Magn Reson Imaging 31: 85–93.
  18. 18. Lewis J (1996) Fast normalized cross-correlation. Industrial Light and Magic.
  19. 19. Minetti AE, Ardigò LP, Saibene F (1994) Mechanical determinants of the minimum energy cost of gradient running in humans. J exp Biol 195: 211–225.
  20. 20. Dempster WT, Gabel WC, Felts WJ (1959) The anthropometry of the manual work space for the seated subject. Am J Phys Anthropol 17: 289–317.
  21. 21. Winter D (2005) Biomechanics and motor control of human movement. New York: John Wiley and Sons, Inc, Third Edition.
  22. 22. Minetti AE (2009) The mathematical description (Lissajous contour) of the 3D trajectory of the body centre of mass: A locomotor ‘signature’ for the physiology, biomechanics and pathology of human and animal gaits. Gait & Posture 30(S): S153–S153.
  23. 23. Minetti AE, Cisotti C, Mian OS (2011) The mathematical description of the body centre of mass 3D path in human and animal locomotion. J Biomech 44: 1471–1477.
  24. 24. Di Prampero PE (1986) The energy cost of human locomotion on land and in water. Int J Sports Med 7(2): 55–72.
  25. 25. Margaria R, Cerretelli P, Aghemo P, Sassi G (1963) Energy cost of running. J Appl Physiol 18: 367–370.
  26. 26. Daniels JT (1985) A physiologist’s view of running economy. Med Sci Sports Exerc 17: 332–338.
  27. 27. Slawinski JS, Billat VL (2004) Difference in mechanical and energy cost between highly, well, and nontrained runners. Med Sci Sports Exerc 36: 1440–1446.
  28. 28. Beneke R, Hutler M (2005) The effect of training on running economy and performance in recreational athletes. Med Sci Sports Exerc 37: 1794–1799.
  29. 29. McGregor SJ, Busa MA, Yaggie JA, Bollt EM (2009) High resolution MEMS accelerometers to estimate VO2 and compare running mechanics between highly trained inter-collegiate and untrained runners. PLoS One 4: e7355.
  30. 30. Cavanagh PR, Pollock ML, Landa J (1977) A biomechanical comparison of elite and good distance runners. Ann N Y Acad Sci 301: 328–45.
  31. 31. Nakayama Y, Kudo K, Ohtsuki T (2010) Variability and fluctuation in running gait cycle of trained runners and non-runners. Gait & Posture 31: 331–335.
  32. 32. Gurney B (2002) Leg length discrepancy. Gait & Posture 15(2): 195–206.
  33. 33. Saibene F, Minetti AE (2003) Biomechanical and physiological aspects of legged locomotion in humans. Eur J Appl Physiol 88: 297–316.
  34. 34. Kong PW, Candelaria NG, Smith DR (2009) Running in new and worn shoes: a comparison of three types of cushioning footwear. Br J Sports Med 43(10): 745–9.
  35. 35. Hardin EC, van den Bogert AJ, Hamill J (2004) Kinematic adaptations during running: effects of footwear, surface, and duration. Med Sci Sports Exerc 36(5): 838–44.
  36. 36. Alexander RM (1989) Optimization and gaits in the locomotion of vertebrates. Physiol Rev 69: 1199–1227.
  37. 37. Srinivasan M (2011) Fifteen observations on the structure of energy-minimizing gaits in many simple biped models. J R Soc Interface 8: 74–98.