- Open Access
- Open Peer Review
Sex-related differences in coordination and variability among foot joints during running
Journal of Foot and Ankle Research volume 11, Article number: 53 (2018)
Women, as compared with men, have a higher proportion of injuries in the ankle/foot region. However, the reason for this sex-related difference in foot injuries remains unclear. Recently, joint coordination and variability of coordination have been suggested to be a critical index for defining both the state of injury and the potential risk of injury. The purpose of this study was to investigate sex-related differences in coordination and variability among the foot joints during running.
Twelve healthy men and 12 healthy women ran on a treadmill. A modified vector coding technique was used to identify coordination and variability among foot joints involving the shank, rearfoot, midfoot, and forefoot segments, and categorized into the following four coordination patterns: in-phase with proximal dominancy, in-phase with distal dominancy, anti-phase with proximal dominancy, and anti-phase with distal dominancy.
There were no differences in all spatiotemporal parameters and in the foot strike angle between men and women. Coordination of variability of the foot joints during running was similar between men and women, but the anti-phase with proximal dominancy in proportion of frontal rearfoot-shank vs. midfoot-rearfoot couple (men; 7.2%, women; 13.9%) and midfoot-rearfoot vs. forefoot-midfoot couple (men; 18.6%, women; 39.8%) in women was significantly increased compared to that in men. Other all coordination of the foot joints during running differed between men and women, and effect sizes of these parameters were all large.
The results may be useful for understanding the underlying mechanism contributing to differences in injury risk in men and women, and may provide novel data on foot joint coordination and variability that could be used as reference data for both biomechanical and clinical running studies.
The increased awareness of exercise as a necessity for maintaining a healthy lifestyle has made running more popular than ever. While there are many beneficial effects of running, the incidence of running injuries is alarmingly high. One study  has suggested that 19.4–79.3% of runners are injured each year. Previous studies [2, 3] have reported a difference in the incidence of ankle/foot injuries between men and women. Frisch et al.  investigated the incidence of sports-related injuries in young athletes and found that ankle/foot injuries were the most common in men and women. Furthermore, women had a significantly higher proportion of injuries in the ankle/foot area than men (34.8% and 16.8%, respectively). In particular, overuse injuries of the foot are common among both elite and recreational runners . Among overuse injuries of the foot, plantar fasciitis, which is a common foot injury in runners, was reported in previous study  to have an increased prevalence in women. Other researchers  have also reported that the incidence of metatarsal stress fractures in women is more than 3.5 times that of men.
Many previous studies [7,8,9,10] have investigated traditional kinematic parameters such as range of motion and peak value during dynamic tasks using single joint analysis to identify the reason for the difference between the sexes with respect to the incidence of running injuries. For example, Takabayashi et al.  found that women had a greater range of motion and peak values in rearfoot and midfoot motions during running compared to men and reported that this result may provide the reason for the high incidence of foot injuries in women. While potential changes of information on the running injury mechanism have been demonstrated, the motor control mechanism behind these changes is not clear . Additionally, abnormal motion in one joint may influence the kinematics in other joints. However, traditional kinematic parameters cannot simultaneously detect two-joint coordinated relationships.
Other researchers have investigated joint or segment coordination based on vector coding technique, which is a dynamic system approaches and reported that coordination information allows a more sensitive measure of joint mechanism . Altered joint or segment coordination results from a change in either the relative timing or amplitude of motion, and this has been suggested to be a cause of running injuries . For example, compared with healthy runners, injured runners demonstrate altered coordination between the thigh and shank . Moreover, coordination between the shank and rearfoot differs between men and women runners . Taking this into consideration, coordination among the foot joints may be altered between men and women during running, since the incidence of ankle/foot injuries between men and women is different. However, the influence of sex-related differences on coordination among the foot joints remains unclear.
Moreover, variability of coordination has been suggested to be a critical index for defining both the state of injury and the potential risk of injury . In healthy able-bodied individuals, there are many available individual degrees of freedom that can be combined or coordinated to achieve a movement task. Hence, low variability means less flexibility of the body to adapt to changing situations, which may lead to greater repetitive stresses on the lower extremity joints . As the frequently used dynamic system approach for quantifying coordination variability, there is a continuous relative phase in addition to vector coding technique. While the continuous relative phase provides spatio-temporal information from the position and velocity signals, this technique has limitations in quantifying non-sinusoidal couplings and may be not appropriate for most lower extremity couplings during the gait . Thus, quantifying coordination variability is often performed using the vector coding technique in many studies, for example runners with patellofemoral syndrome exhibit less coordination variability than healthy runners by vector coding technique . Moreover, women have decreased variability of coordination at the intra-lower extremity (e.g. between thigh and shank) compared to men . While most of the previous studies have focused on variability of coordination in the intra-lower extremity, it is unknown if or how foot joint coordination variability changes in men and women. By conducting this study, we hope to gain a greater understanding of the underlying mechanism contributing to differences of injury risk in men and women.
The purpose of this study was to investigate sex-related differences in coordination and variability among the foot joints during running. Owing to increased foot injury rates in women [2, 3], we hypothesized that there would be a difference in coordination between men and women, and that its variability would be decreased in women compared to men.
Twelve healthy men (age = 20.7 (1.2) years; height = 1.71 (0.05) m; weight = 64.8 (9.0) kg) and twelve women (age = 20.4 (1.0) years; height = 1.60 (0.06) m; weight = 52.8 (8.5) kg) participated in this study. Inclusion criteria were as follows: no history of or currently present with lower limb injuries and normal foot based on the arch height index during 90% weight bearing. The reference value for the normal foot is considered to be 1.5 standard deviation (SD) above or below the mean arch ratio measurement of 0.292 (SD 0.027), based on a previously examined sample population of 102 ft . In this study, all participants provided written informed consent before participation. This study was approved by the ethics committee of our institution (No. 17724–160902).
Reflective markers of 9 mm in diameter were fixed to the right shank and foot at the tibial tuberosity, fibula head, medial malleolus, lateral malleolus, Achilles tendon attachment, sustentaculum tali, peroneal tubercle, navicular bone, cuboid, first metatarsal base, first metatarsal head, second metatarsal base, second metatarsal head, fifth metatarsal head, and head of the proximal phalanx of the hallux (Fig. 1), according to the Rizzoli foot model . The reproducibility of this model has been confirmed in a previous study . The marker was also attached at the right posterior superior iliac spine to determine the stance phase.
Prior to the acquisition of running data, static standing data were measured for each participant. These data were used to calculate the offset values for all joint rotations in each participant . After the measurement, the participants performed barefoot running (individual’s jogging speed and preferred cadence) on a treadmill (Auto Runner AR-100: Minato Medical Science, Japan). Since women are generally smaller in body size as compared to men, women are likely to increase their step length and cadence parameters. Hence, all spatiotemporal parameters were normalized by the following: speed normalized by the square root of the acceleration of gravity × leg length; cadence normalized by the square root of the leg length/acceleration of gravity; and step length normalized by leg length . Such analysis is often used to investigate sex differences in lower extremity kinematics in previous studies [7, 24].
The participants were allowed to practice the tasks repeatedly for more than one minute on the treadmill. The examiner confirmed with the participants through questioning whether they were accustomed to the motion. After all the subjects verbally reported feeling comfortable running on the treadmill, the task was measured for 30 s. From the strides measured in 30 s, five strides were randomly extracted for each participant. The examiner confirmed the question that all participants had no fatigue prior to measurement of the task or ending the task. The static data and task were measured using a three-dimensional motion analysis system (Vicon, Oxford, United Kingdom) that included 13 infrared cameras at a sampling rate of 250 Hz.
Raw marker trajectory data during running were filtered using a second-order, zero-lag Butterworth low-pass filter with a 12-Hz cut-off frequency . The following four segments were defined in the kinematics model: the shank comprising the tibia and fibula; the rearfoot (i.e., calcaneus); the midfoot comprising the navicular, cuneiform, and cuboid bones; and the forefoot comprising the first to fifth metatarsal bones. In this study, the three-dimensional joint angles were calculated at the distal segment and expressed relative to the adjacent proximal segment using a right-handed orthogonal Cardan Xyz sequence of rotations (a sequence of plantarflexion/dorsiflexion, eversion/inversion, and abduction/adduction) , which was selected to be equivalent of the joint coordinate system . Hence, the joint angles were calculated as plantarflexion/dorsiflexion, eversion/inversion, and abduction/adduction of the rearfoot with respect to the shank (RF-SH), midfoot with respect to the rearfoot (MF-RF), and forefoot with respect to the midfoot (FF-MF).
After calculation of the joint angles during the task, the data were time-normalized to the stance phase (100 data points). The stance phase was determined from the marker trajectory data of the calcaneus, second metatarsal head, and posterior superior iliac spine using Smith’s custom designed algorithm . Briefly, foot strike was defined to occur at the point of maximum in vertical displacement between the calcaneus and posterior superior iliac spine markers. Similarly, toe-off was defined to occur at the point of maximum in vertical displacement between the second metatarsal head and posterior superior iliac spine markers. This method can provide an accurate estimation of foot strike and toe-off. The foot strike angle  was also calculated because differences in foot strike angles affect intra-foot kinematics and coupling of the foot joints .
Calculation of the coupling angle and variability
where 0° ≤ γ ≤ 360°, xi, and yi represent the proximal and distal joint angles, respectively. In addition, i represents the percent stance of the jth stride. To determine the coupling angle for a single participant (i.e., five strides) and among the participants, the mean coupling angle was calculated from the mean xi and the mean yi at each percentage of stance. The calculations were performed according to circular statistics . Coordination variability was calculated as the circular SD of the coupling angle for each individual across five strides of data.
The foot joint couplings of interest included the following: (1) sagittal RF-SH vs. sagittal MF-RF, (2) frontal RF-SH vs. frontal MF-RF, (3) sagittal MF-RF vs. sagittal FF-MF, and (4) frontal MF-RF vs. frontal FF-MF. These pairs of joint rotations have been investigated in previous studies . They were chosen to show stronger kinematic coupling during running . The coupling angle represents an instantaneous spatial relationship from which four coordination patterns (Fig. 2b) can be identified: (i) in-phase with proximal dominancy (the same direction and greater angular amplitude of the proximal joint), (ii) in-phase with distal dominancy (the same direction and greater angular amplitude of the distal joint), (iii) anti-phase with proximal dominancy (the opposite direction and greater angular amplitude of the proximal joint), and (iv) anti-phase with distal dominancy (the opposite direction and greater angular amplitude of the distal joint) . In the present study, the positive direction of joint rotation was defined as dorsiflexion and inversion.
There are two functional phases, absorption and propulsion, during the stance phase of running, and the two phases switch at the moment when rearfoot eversion changes to rearfoot inversion . The two phases were divided based on the peak value of rearfoot eversion in each participant. The mean coupling angles were categorized into one of the four coordination patterns at each phase, as shown in Fig. 3. All computations were performed in Scilab version 6.0.0 (Enterprises, Versailles, France).
An independent t-test was used to detect sex-related differences among normalized spatiotemporal parameters, foot strike angles, and foot joint coordination patterns in each functional phase. The alpha level for all tests was set at 0.05. Statistics testing was performed using R version 3.3.2 (The R Foundation for Statistical Computing, Austria). Effect sizes (ES) were also calculated by using Cohen’s d statistics. The ES were evaluated on the following criteria: trivial (0–0.19), small (0.20–0.49), medium (0.50–0.79), and large (> 0.80).
While speed (men: 2.17 (0.12) ms− 1, women: 2.04 (0.13) ms− 1) and step length (men: 0.78 (0.04) m, women: 0.71 (0.06) m) significantly differed between the men and women (P < 0.05), there were no differences in the normalized speed (men: 0.75 (0.04), women: 0.73 (0.04), P > 0.05), normalized cadence (men: 49.1 (2.4), women: 49.3 (2.9), P > 0.05), normalized step length (men: 0.92 (0.05), women: 0.88 (0.07), P > 0.05), and foot strike angle (men: 16.8° (5.1°), women: 15.6° (3.2°), P > 0.05) between men and women.
In the sagittal plane, RF-SH vs. MF-RF couple, during the absorption phase, women showed a significant increase in the proportion of in-phase with distal dominancy (the same direction and greater angular amplitude of MF-RF) compared with men (men: 11% (6.1%), women; 25.4% (12.3%), Fig. 4a). During the absorption phase, in the sagittal MF-RF vs. FF-MF couple, men showed a significant increase in the proportion of anti-phase with distal dominancy (the opposite direction and greater angular amplitude of MF-RF) compared with men (men; 17.0% (7.4%), women; 9.1% (5.7%), Fig. 4b). In the sagittal MF-RF vs. FF-MF couple during the propulsion phase, women showed a significant increase in the proportion of in-phase with proximal dominancy (the same direction and greater angular amplitude of MF-RF) compared to men (men; 43.6% (21.6%), women; 60.6% (12.0%), Fig. 4b). In the frontal RF-SH vs. frontal MF-RF couple during the propulsion phase, women showed a significant increase in the proportion of anti-phase with proximal dominancy (the opposite direction and greater angular amplitude of RF-SH) compared to men (men; 7.2% (7.8%), women; 13.9% (7.5%), Fig. 4c). In the frontal MF-RF vs. frontal FF-MF couple during the absorption phase, women showed a significant increase in the proportion of anti-phase with proximal dominancy (the opposite direction and greater angular amplitude of MF-RF) compared to men (men; 18.6% (14.3%), women; 39.8% (24.1%), Fig. 4d). The ES of these significant parameters were all large.
Variabilities in coordination between the foot joints are shown in Fig. 5. Coordination variabilities for sagittal RF-SH vs. MF-RF, frontal RF-SH vs. MF-RF, sagittal MF-RF vs. FF-MF and frontal MF-RF vs. FF-MF were similar between men and women, and all coordination variabilities were not significantly different (Fig. 5a-d). Moreover, the ES of these parameters were trivial to small. Complete results of coordination and variability are displayed in the Additional file 1: Table S1.
The purpose of this study was to investigate sex-related differences in coordination and variability between the foot joints during running. The results of this study support our first hypothesis that foot joint coordination during running differed between men and women. Contrary to our second hypothesis that coordination variability between foot joints would be lower in women than in men, variabilities in joint coordination were similar between groups.
Some of the coordination patterns were significantly different between men and women (e.g. in-phase with distal dominancy in sagittal RF-SH vs. MF-RF couple). As a rationale for the difference in coordination patterns by sex, DeLeo et al.  reported that altered joint or segment coordination results from a change in either the relative timing or amplitude of motion. The coordination changes in this study could also be a result of alternating joint motion amplitudes, alternating joint motion timing, or a combination of both .
In many previous studies [11,12,13,14, 35], altered joint coordination patterns have been suggested to be a cause or marker of running injury. Chang et al.  reported that anti-phase motion in the frontal plane leads to torsion of the soft tissues. It is thought that this motion may overload the soft tissues or joints. Actually, the coordination pattern between trunk and pelvis rotation in the transverse plane in patients with low back pain decreases in the anti-phase and increases in the in-phase compared to controls. This is because patients adjust to avoid or minimize pain [35, 36]. Additionally, sex-related differences in excessive motions of the lower extremities during running have been suggested to be a contributing factor for running injuries [7, 10, 37], for example women showed excessive midfoot movements in the sagittal plane compared to men . The proportion of frontal RF-SH vs. MF-RF and MF-RF vs. FF-MF couples in women was significantly increased in the anti-phase with proximal dominancy. Therefore, the results of this analysis contribute to the factor that women are more likely to have ankle/foot injuries than men. However, this suggestion will require investigation in a longitudinal study.
In this study, all coordination variabilities were not significantly different. Boyer et al.  showed that coordination variability in transverse thigh and shank rotation was reduced in women compared to men during running, but those results were inconsistent with the results of this study. However, a previous study  reported that the coordination variabilities of the intra-lower extremity were similar between young people and healthy elderly people who have high injury rates. This is due to the fact that those elderly people were healthy, which potentially prevented the injury. Similarly, changes in coordination variability in men and women in this study may have not shown because the target participants were healthy and had no history of disability. Additionally, it was thought to be due to differences in spatiotemporal parameters. Previous studies [14, 39] set the same running speed between men and women. However, because women are generally smaller in body size than men, women are more likely to increase the step length and cadence parameters as compared to men, which may have affected the preferred movement patterns of women. Recent walking and running analysis study [7, 24] has addressed an analysis of sex-related difference by removing the influence of those spatiotemporal parameters. Since the influence of the spatiotemporal parameters was eliminated in this study, coordination variabilities were not significantly different between men and women. Further research into the influence of such factors on variability measures is warranted.
The present study has certain limitations. First, this study was performed with all participants running at a jogging speed because there is a possibility that pain will appear in the foot when barefoot fast running. However, as running kinematics may change when increasing the running speed, the results of this study could have been different under faster running conditions. Second, the present study includes the use of skin-mounted markers to track the underlying skeletal structures. While the Rizzoli foot model used in this study has been validated in an in vitro study , this model has not been validated in an in vivo study. Thus, skin markers mounted on externally identifiable bony landmarks in the foot may not follow the underlying individual skeletal segments to be properly evaluated during running. Finally, since this study used a foot model composed of the rearfoot, midfoot, and forefoot, the motion of each foot bone, such as the metatarsals, cannot be investigated. For example, while the forefoot segment in this study comprises the first to fifth metatarsal bones, previous researches on bone pins have reported that there is considerable movement in each metatarsal bone during running and walking . These limitations will need to be addressed in future studies.
In this study, coordination of variability of the foot joints during running was similar between men and women, but coordination of the foot joints during running differed between men and women. The results of this analysis may be useful for understanding the underlying mechanism contributing to differences in injury risk in men and women. Furthermore, these results provide novel data on foot joint coordination and variability that could be used as reference data for both biomechanical and clinical running studies.
Forefoot with respect to the midfoot
Midfoot with respect to the rearfoot
Rearfoot with respect to the shank
van Gent RN, Siem D, van Middelkoop M, van Os AG, Bierma-Zeinstra SM, Koes BW. Incidence and determinants of lower extremity running injuries in long distance runners: a systematic review. Br J Sports Med 2007; 41:469–480;
Frisch A, Seil R, Urhausen A, Croisier JL, Lair ML, Theisen D. Analysis of sex-specific injury patterns and risk factors in young high-level athletes. Scand J Med Sci Sportss. 2009;19:834–41.
Orr RM, Pope R. Gender differences in load carriage injuries of Australian army soldiers. BMC Musculoskelet Disord. 2016;17:488.
Kindred J, Trubey C, Simons SM. Foot injuries in runners. Curr Sports Med Rep. 2011;10:249–54.
Rano JA, Fallat LM, Savoy-Moore RT. Correlation of heel pain with body mass index and other characteristics of heel pain. J Foot Ankle Surg. 2001;40:351–6 Davis PF, Severud E, Baxter DE. Painful heel syndrome: results of nonoperative treatment. Foot Ankle Int. 1994; 15:531–5.
Ohta-Fukushima M, Mutoh Y, Takasugi S, Iwata H, Ishii S. Characteristics of stress fractures in young athletes under 20 years. J Sports Med Phys Fitness. 2002;42:198–206.
Takabayashi T, Edama M, Nakamura M, Nakamura E, Inai T, Kubo M. Gender differences associated with rearfoot, midfoot, and forefoot kinematics during running. Eur J Sport Sci. 2017;17:1289–96.
Sinclair J, Taylor PJ. Sex differences in tibiocalcaneal kinematics. Human Movement. 2014;15:105–9.
Almonroeder TG, Benson LC. Sex differences in lower extremity kinematics and patellofemoral kinetics during running. J Sports Sci. 2016:1–7.
Ferber R, Davis IM, Williams DS 3rd. Gender differences in lower extremity mechanics during running. Clin Biomech (Bristol, Avon). 2003;18:350–7.
Hafer JF, Freedman Silvernail J, Hillstrom HJ, Boyer KA. Changes in coordination and its variability with an increase in running cadence. J Sports Sci. 2016;34:1388–95.
Hamill J, van Emmerik RE, Heiderscheit BC, Li L. A dynamical systems approach to lower extremity running injuries. Clin Biomech (Bristol, Avon). 1999;14:297–308.
DeLeo AT, Dierks TA, Ferber R, Davis IS. Lower extremity joint coupling during running: a current update. Clin Biomech (Bristol, Avon). 2004;19:983–91.
Rodrigues P, Chang R, TenBroek T, van Emmerik R, Hamill J. Evaluating the coupling between foot pronation and tibial internal rotation continuously using vector coding. J Appl Biomech. 2015;31:88–94.
Hamill J, Palmer C, Van Emmerik RE. Coordinative variability and overuse injury. Sports Med Arthrosc Rehabil Ther Technol. 2012;4:45.
Silvernail JF, Boyer K, Rohr E, Bruggemann GP, Hamill J. Running mechanics and variability with aging. Med Sci Sports Exerc. 2015;47:2175–80.
Peters BT, Haddad JM, Heiderscheit BC, Van Emmerik RE, Hamill J. Limitations in the use and interpretation of continuous relative phase. J Biomech. 2003;36:271–4.
Heiderscheit BC. Variability of stride characteristics and joint coordination among individuals. J Appl Biomech. 2002;18:110–21.
Pollard CD, Heiderscheit BC, van Emmerik RE, Hamill J. Gender differences in lower extremity coupling variability during an unanticipated cutting maneuver. J Appl Biomech. 2005;21:143–52.
Williams DS, McClay IS. Measurements used to characterize the foot and the medial longitudinal arch: reliability and validity. Phys Ther. 2000;80:864–71.
Leardini A, Benedetti MG, Berti L, Bettinelli D, Nativo R, Giannini S. Rear-foot, mid-foot and fore-foot motion during the stance phase of gait. Gait Posture. 2007;25:453–62.
Seo SG, Lee DY, Moon HJ, Kim SJ, Kim J, Lee KM, Chung CY, Choi IH. Repeatability of a multi-segment foot model with a 15-marker set in healthy adults. J Foot Ankle Res. 2014;7:24.
Hof AL. Scaling gait data to body size. Gait & Posture. 1996;(3):222–3.
Pinzone O, Schwartz MH, Baker R. Comprehensive non-dimensional normalization of gait data. Gait & Posture. 2016;44:68–73.
Chang R, Rodrigues PA, Van Emmerik RE, Hamill J. Multi-segment foot kinematics and ground reaction forces during gait of individuals with plantar fasciitis. J Biomech. 2014;47:2571–7.
Grood ES, Suntay WJ. A joint coordinate system for the clinical description of three-dimensional motions: application to the knee. J Biomech Eng. 1983;105:136–44.
Smith L, Preece S, Mason D, Bramah C. A comparison of kinematic algorithms to estimate gait events during overground running. Gait Posture. 2015;41:39–43.
Altman AR, Davis IS. A kinematic method for footstrike pattern detection in barefoot and shod runners. Gait Posture. 2012;35:298–300.
Pohl MB, Buckley JG. Changes in foot and shank coupling due to alterations in foot strike pattern during running. Clin Biomech (Bristol, Avon). 2008;23:334–41.
Chang R, Van Emmerik R, Hamill J. Quantifying rearfoot-forefoot coordination in human walking. J Biomech. 2008;41:3101–5.
Batschelet E. Circular statistics in biology: Academic press London; 1981.
Takabayashi T, Edama M, Yokoyama E, Kanaya C, Kubo M. Quantifying coordination among the rearfoot, midfoot, and forefoot segments during running. Sports Biomech. 2018;17:18–32.
Takabayashi T, Edama M, Yokoyama E, Chiaki K, Inai T, Tokunaga Y, Kubo M. Changes in kinematic coupling among the rearfoot, midfoot, and forefoot segments between running and walking. J Am Podiatr Med Assoc. 2016;108(1):45–51.
Needham RA, Naemi R, Chockalingam N. A new coordination pattern classification to assess gait kinematics when utilising a modified vector coding technique. J Biomech. 2015.
Seay JF, Van Emmerik RE, Hamill J. Low back pain status affects pelvis-trunk coordination and variability during walking and running. Clin Biomech (Bristol, Avon). 2011;26:572–8.
Seay JF, Van Emmerik RE, Hamill J. Influence of low back pain status on pelvis-trunk coordination during walking and running. Spine 2011; 36:E1070-E10E9.
Wearing SC, Smeathers JE, Urry SR, Hennig EM, Hills AP. The pathomechanics of plantar fasciitis. Sports Med. 2006;36:585–611.
Boyer KA, Freedman Silvernail J, Hamill J. Age and sex influences on running mechanics and coordination variability. J Sports Sci. 2017;35:2225–31.
Hafer JF, Boyer KA. Variability of segment coordination using a vector coding technique: reliability analysis for treadmill walking and running. Gait Posture. 2017;51:222–7.
Liu A, Nester C, Ward E, Howard D, Derrick T, Cocheba J. In vitro study of foot kinematics using a walking simulator. Proceeding of biomechanics of the lower limb in health, disease and rehabilitation. 2005:222–3.
Arndt A, Wolf P, Liu A, Nester C, Stacoff A, Jones R, Lundgren P, Lundberg A. Intrinsic foot kinematics measured in vivo during the stance phase of slow running. J Biomech. 2007;40:2672–8.
We thank the members of the Institute for Human Movement and Medical Science, Niigata University of Health and Welfare for their assistance, and Editage (www.editage.jp) for English language editing.
This study was supported by a Research Activity Start-up grant (16H07304) from the Japan Society for the Promotion of Science and Grant-in-Aid for Exploratory Research (H29B17) in our affiliation.
Availability of data and materials
Data may be made available by the authors upon reasonable request.
Ethics approval and consent to participate
All participants provided informed consent prior to participation. The present study was reviewed and approved by the ethical committee in our institution (No. 17724–160902).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.