Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 28 June 2022
Sec. Exercise Physiology

Not Lower-Limb Joint Strength and Stiffness but Vertical Stiffness and Isometric Force-Time Characteristics Correlate With Running Economy in Recreational Male Runners

  • 1School of Physical Education and Sport Training, Shanghai University of Sport, Shanghai, China
  • 2Physical Education Department, College of Education, United Arab Emirates University, Al Ain, United Arab Emirates
  • 3Department of Sports Science and Clinical Biomechanics, SDU Sport and Health Sciences Cluster (SHSC), University of Southern Denmark, Odense, Denmark
  • 4School of Kinesiology, Shanghai University of Sport, Shanghai, China

Neuromuscular characteristics, such as lower-limb joint strength, the ability to reuse elastic energy, and to generate force are essential factors influencing running performance. However, their relationship with running economy (RE) remains unclear. The aim of this study was to evaluate the correlations between isokinetic lower-limb joint peak torque (PT), lower-limb stiffness, isometric force-time characteristics and RE among recreational-trained male runners. Thirty male collegiate runners (aged 20–22 years, VO2max: 54.02 ± 4.67 ml·kg−1·min−1) participated in test sessions on four separate days. In the first session, the body composition and RE at 10 km·h−1 were determined. In the second session, leg and vertical stiffness (Kleg and Kvert), knee and ankle stiffness (Kknee and Kankle) were evaluated. In the third session, isokinetic knee and ankle joint PT at velocity of 60°s−1 were tested. The force-time characteristics of isometric mid-thigh pull (IMTP) were evaluated in the final session. The Pearson’s product-moment correlations analysis shows that there were no significant relationships between knee and ankle joint concentric and eccentric PT, Kknee and Kankle, Kleg, and RE at 10 km·h−1. However, Kvert (r = −0.449, p < 0.05) and time-specific rate of force development (RFD) for IMTP from 0 to 50 to 0–300 ms (r = −0.434 to −0.534, p < 0.05) were significantly associated with RE. Therefore, superior RE in recreational runners may not be related to knee and ankle joint strength and stiffness. It seems to be associated with vertical stiffness and the capacity to rapidly produce force within 50–300 ms throughout the lower limb.

Introduction

The running economy (RE) is defined as the steady-state oxygen or energy expenditure at a given running speed and is an essential physiological contributor to distance running success (Blagrove et al., 2018). RE reflects the energy demand during constant submaximal running and is cited as a stronger indicator of endurance performance than maximum oxygen uptake (VO2max) among a homogenous group of runners with similar aerobic capacity (Conley and Krahenbuhl, 1980; Morgan et al., 1989). Thus, a superior RE means that runners use less energy to sustain a given steady-state speed. Given the critical contribution of RE to running performance, exploring the principal factors of RE appears essential. Neuromuscular strength is an important factor in the determination of RE, and many studies have indicated that improvement in RE following strength training interventions is attributed to neuromuscular adaptation, such as a consequence of changes in lower-limb muscle strength, stiffness and power (Barnes and Kilding, 2015; Blagrove et al., 2018; Trowell et al., 2020).

Lower limb muscles, especially those at the knee and ankle play an important role in running (Hamner et al., 2010; Hamner and Delp, 2013). For example, during the stance phase, the knee extensor muscles are the main contributors to decelerate and support the center of mass (COM) at the braking phase (Hamner et al., 2010; Hamner and Delp, 2013) and allow the force production in the knee flexor muscles during the propulsive phase (Liu et al., 2008; Monte et al., 2020). The ankle plantar flexor muscles are also activated during running, generating force more than 12 times the body weight (Komi, 1990), and are the most dominant contributors to force during the propulsion phase (Hamner et al., 2010; Monte et al., 2020). Meanwhile, ankle plantar flexor muscles perform quasi-isometric contractions while generating muscular contraction forces during the braking phase and approach their optimal contraction length, thus facilitating greater displacement of the Achilles tendon for stretch and recoil, which facilitates the storage and recovery of elastic potential energy (Farris and Sawicki, 2012; Lai et al., 2014; Monte et al., 2020). Similarly, the knee extensor muscles also perform quasi-isometric behavior and near-optimal length during the propulsion phase (Bohm et al., 2018; Monte et al., 2020), which allows for more economical production of muscle force (Fletcher and MacIntosh, 2017; Monte et al., 2020). Therefore, neuromuscular characteristics, such as knee and ankle joint muscle strength, may play an important role in improving RE.

However, the relationship between lower-limb joint strength and RE remains unclear. For instance, Sundby and Gorelick. (2014) found no significant relationship between knee flexor, extensor muscles peak torque (PT) and RE in female runners, whereas Westblad et al. (1996) found a modest correlation between eccentric knee extensor strength and RE in male distance runners. In the aspect of ankle joint strength, Bohm et al. (2021) reported that an increase in plantar flexor muscles (i.e., soleus) strength reduces the energy cost of running. To our knowledge, none of the studies have directly investigated the association between ankle joint strength and RE. Therefore, the scientific literature requires a comprehensive investigation of the correlation between lower-limb joint strength and RE.

In athletic performance tasks, lower-limb stiffness is commonly reflected by vertical stiffness (Kvert), leg stiffness (Kleg) and joint stiffness (Maloney and Fletcher, 2021). These measures are identified by different methods, including the spring-mass (Kvert and Kleg) and torsional spring model (joint stiffness). Kvert and Kleg were considered to be contributors to the improvement of RE by reducing energy cost in vertical movements (Heise and Martin, 2001) and increasing the elastic energy storage capacity of the leg muscle-tendon units (Dalleau et al., 1998). However, the results of the existing studies seem contradictory. For example, Man et al. (2016) indicated that both Kvert and Kleg were significantly correlated with RE, whereas Heise and Martin. (1998) found that Kvert was moderately correlated with RE, while Kleg was not. Furthermore, the lower-limb joints constitute a multi-spring system with different elasticities and viscosities (Kuitunen et al., 2002). The torsional spring model provides a different perspective on spring-mass model, while knee or ankle stiffness (Kknee and Kankle) has the greatest effect on leg-spring stiffness during running (Struzik et al., 2021). However, few studies have examined the associations of Kknee and Kankle with RE. To date, only Tam et al. (2018) reported that greater Kknee and lesser Kankle were significantly associated with better RE in well-trained runners. Therefore, the comprehensive investigation of the relationship between Kvert, Kleg and joint stiffness, and RE were essential to facilitate the running performance.

The runner’s capabilities to develop force and power within the transition period during landing could improve RE (Lum et al., 2020). The isometric mid-thigh pull (IMTP) is an accurate and reliable test that measures force-time characteristics, such as peak force (PF) generation and rate of force development (RFD) across various specific times (Wang et al., 2016; Brady et al., 2019). Running involves multi-joint movements and muscle activation (Novacheck, 1998), and IMTP enables the hip, knee and ankle joints to be held at a relatively specific angle during running. Therefore, it is essential to determine the association between the whole lower limb force-generating ability and RE. Currently, only Lum et al. (2020) have reported a relationship between IMTP characteristics and RE. They found that RFD of 0–100 to 0–200 milliseconds (ms) were significantly correlated with RE at 12 km·h−1 in recreational runners. Noticeably, the time that muscles are activated during the contact phase is 200–350 ms at submaximal speeds (10–16 km·h−1) (Gómez-Molina et al., 2017; García-Pinillos et al., 2019a). Therefore, it is important to investigate the RFD over 200 ms in order to have a better understanding of the relationship between force-time characteristics and RE.

To summary, regarding recreational runners, the benefits of neuromuscular factors on the RE have not been sufficiently studied (Sundby and Gorelick, 2014; Silva et al., 2018; Lum et al., 2020). As the increasing number of recreational runners worldwide, identifying these key neuromuscular profiles may benefit runners as this knowledge could be applied to improve RE and running performance through a more optimal strength training program design. The aim of this study was to evaluate the relationship between neuromuscular characteristics and RE in recreational male runners. We hypothesised that knee and ankle joint strength, lower-limb stiffness (Kvert, Kleg, Kknee, and Kankle), and IMTP force-time characteristics (PF and RFD at 0–50 to 0–350 ms) would be significantly associated with RE.

Materials and Methods

Subjects

Thirty recreational-trained male runners (aged 20–22 years, VO2max: 54.02 ± 4.67 ml·kg−1·min−1) from the collegiate running club, who had previously competed in 5- to 21-km races at the collegiate level and had a minimum of 2 years of distance running experience, volunteered for this study. The sample size also satisfied the power requirement of correlation bivariate normal model with G*power 3.1 software (Universitat Dusseldorf, Dusseldorf, Germany), using the setting two-sided test, correlation ρ H1 = 0.5, α err prob = 0.05, Power (1-β err prob) = 0.80, Correlation ρ H0 = 0 (Cohen, 1988; Faul et al., 2009). The subjects’ basic information and physical characteristics are presented in Tables 1. Before the study, all subjects completed the Physical Activity Readiness Questionnaire (PAR-Q), and they respond “no” to all questions. Moreover, they received at least 1 year of training in long-distance running (5–21 km) and ran 20–30 km per week for 3 months prior to the study. Each subject was fully informed about the potential risks and procedures of the experiment and signed an informed consent document. This study was approved by Ethics Committee of Shanghai University of Sport, China (ID number: 2,017,047).

TABLE 1
www.frontiersin.org

TABLE 1. Physical, physiological, and neuromuscular characteristics of the participants (n = 30).

The Experimental Approach to the Problem

Each participant performed four separate test sessions in the laboratory with a rest interval of at least 48 h. The tests flowchart is shown in Figure 1. Before each session, they were instructed to achieve a non-fatigue state (more than 8 h of adequate sleep, no vigorous exercise for 24 h, and no muscle soreness and fatigue). The researcher examined the participants’ commitment to the instructions with a questionnaire when they reported to the laboratory for testing. In the first session, body composition and RE at the speeds of 10 km·h−1 were determined. In the second session, the lower-limb stiffness (Kvert, Kleg, Kknee, and Kankle) at a speed of 10 km·h−1 was calculated. In the third session, the isokinetic knee and ankle joint concentric and eccentric PT at 60°s−1 velocity were tested, while the force-time characteristics of the IMTP (PF and RFD at 0–50 to 0–350 ms) were evaluated in the final session. Pearson’s correlation coefficients were used to investigate the possible relationships.

FIGURE 1
www.frontiersin.org

FIGURE 1. Tests flow diagram. RE, running economy; RFD, rate of force development; IMTP, isometric mid-thigh pull.

Procedures

Body Composition Test

On the measurement day, participants were required to avoid drinking and food intake during the 3 h before testing. The weight, fat mass and fat-free mass were determined by using bioimpedance analyser (X-scan Plus II; Jawon, Korea). Standing height was recorded with a wall-mounted measuring device (Butterfly, Shanghai, China) and the body mass index (BMI) was calculated.

Running Economy Test

Running economy was determined using a treadmill protocol, which showed high intraclass correlation coefficient values (ICC), ranging from 0.92 to 0.94 in our laboratory (Li et al., 2021). The participants warmed up on the treadmill at a running speed of 8 km·h−1 for 4 min. After a 5-min rest, they ran at 10 km·h−1 for 4 min to determine RE, which was defined as the mean oxygen uptake (ml·kg−1·min−1) during the last minute, and the steady-state condition is verified by the respiratory exchange ratio (RER) is <1 of present subjects (Barnes and Kilding, 2015). The submaximal speed is set at 10 km·h−1 because this pace is similar to that used in previous study and reflects the runners’ ability to run at submaximal speeds (Piacentini et al., 2013). Oxygen uptake and heart rate (HR) were continuously monitored using a portable metabolic analyser (K5, Cosmed, Italy) and a HR monitor belt (Garmin, Olathe, United States).

Lower-Limb Stiffness Test

Each participant wore identical running shoes and tight pants provided by the research team before the test, avoiding the effect of this variable on the lower-limb stiffness (Kulmala et al., 2018). Kinematic data were captured using an 8-camera Vicon T40 motion analysis system (Oxford Metrics, Oxford, United Kingdom). Ground reaction force (GRF) data was collected using two 90 cm × 60 cm × 10 cm force platforms (9287 B, Kistler Corporation) with a sampling frequency of 1,000 Hz, synchronised with the motion analysis system. The force platforms are located underneath the treadmill belt. As shown in Figure 2, 36 retroreflective markers were placed in the pelvis and lower limbs to define the foot, shank, thigh and pelvic segments. A standing calibration was recorded to identify the length of each segment and leg, the local coordinate system, and the position of the joint centre for each participant. Participants were then instructed to run on a treadmill at 8 km·h−1 for 4 min to warm up. After warming up, the subjects performed running for 4 min at a speed consistent with the RE test (10 km·h−1), during which biomechanical data collected were more valid due to being in a steady state (RER <1) (Tam et al., 2018). According to previous studies, for healthy adults, recording more than six strides is sufficient to obtain representative data, which is defined as a 95% confidence interval and an error within 5% (Besser et al., 1999). Our study took data from 10 consecutive steps starting at the third minute of running at 10 km·h−1 and averages were calculated for further analysis. The start and end of the support period is determined using a vertical force signal of 50 N (Kyröläinen et al., 2001).

FIGURE 2
www.frontiersin.org

FIGURE 2. Placement of reflective markers and experimental site. The reflective markers are placed at the anterior superior iliac spine, superior margin of the iliac crest, posterior superior iliac spine, medial and lateral femoral condyles, medial and lateral malleolus, first and fifth metatarsal heads, toe and heel to identify the hip, knee and ankle joints. The markers on the T stand are used to track the trajectory of the thigh and shank.

Biomechanical data were processed by Visual 3-dimensional (3D) gait analysis software (v5, C-Motion, Inc., Germantown, MD, United States). Kvert and Kleg measurements using the sine-wave method on a treadmill showed high ICC (0.99 and 0.86) (Pappas et al., 2014). They were calculated according to the recommendations of Morin et al. (2005). Kvert is defined as the ratio of the maximum vertical GRF to the COM vertical displacement (Farley and González, 1996), as in Eq. 1. The Kleg is defined as the ratio of the maximum vertical GRF to the change in leg length (Farley and González, 1996), as in Eq. 2. In addition, Lorimer et al. (2018) confirmed that Kknee and Kankle obtained on a treadmill exhibiting an ICC between 0.75 and 0.90. Based on previous researches, joint stiffness (knee and ankle) is interpreted as the ratio of the change in moment to the angular displacement from touchdown to when the joint is flexed to its maximum angle (Hamill et al., 2014; Tam et al., 2018), as in Eq. 3.

Kvert=FmaxΔy(1)

where Fmax is the maximum vertical GRF, and ∆y denotes the COM vertical displacement during stance phase.

Kleg=FmaxΔL,ΔL=LL2(vtc2)2+Δy(2)

where Fmax is the maximum vertical GRF, and ΔL denotes the change in leg length during ground contact. L is leg length in static stance (the distance from the greater trochanter of the femur to the ankle), ∆y is the COM vertical displacement during stance phase, v is the running speed (m·s−1), and tc denotes the ground contact time (s).

jointstiffness=ΔMΔθ(3)

where ΔM is the amount of change in joint moment, and Δθ denotes the articular angular displacement from touchdown to maximum flexion. Kknee is calculated as the ratio of the change in knee moment to the angular displacement from touchdown to maximum knee flexion, and Kankle is calculated as the ratio of the change in ankle moment to the angular displacement from touchdown to maximum ankle dorsiflexion in midstance.

Isokinetic Strength Test

Knee joint flexor/extensor muscles and ankle joint plantar flexor/dorsiflexor muscles PT were measured using a motor-driven dynamometer (D&R Ferstl GmbH, Hemau, Germany) (Figure 3A) as this test showed high ICC (knee: 0.90 to 0.96; ankle: 0.77–0.98) in previous study (Gonosova et al., 2018; Andrade et al., 2021). The subjects performed two trials to access the right lower-limb knee joint flexion/extension and ankle joint plantarflexion/dorsiflexion PT. Before each trial, the instrument was calibrated according to the manufacturer’s manual. The operator needs to place the shaft of the dynamometer in a substantially horizontal position to allow automatic gravity compensation and not to touch the dynamometer with the hands during weighing. The subject remains relaxed and free from muscle activity to eliminate the effects of gravity. All participants performed three sub-maximal trials to familiarise themselves with the test. As the most commonly used angular speed during isokinetic testing procedures is 60°s−1, which has been suggested to reliably and accurately evaluate maximal capacities of the muscles to produce force (Zapparoli and Riberto., 2017; Silva et al., 2018), five repetitions of concentric and eccentric isokinetic PT were evaluated on the knee and ankle joints at 60° s−1 in the present study (Luna et al., 2012; Andrade et al., 2021). There was a rest period of 1 min between contraction type and a 10-min break between trials. In addition, the subjects were offered encouragement to exert their maximum strength during the trial. The PT was calculated as the relative greatest torque value (maximal PT/body weight) during isokinetic concentric and eccentric phases and was collected for further analysis.

FIGURE 3
www.frontiersin.org

FIGURE 3. Strength test. Interpretation of (A) knee and ankle isokinetic strength test (B) Isometric mid-thigh pull test. IMTP, isometric mid-thigh pull.

As shown in Figure 3A, In the first trial, knee strength was tested in a seated position, with hips at approximately 85° flexion. At the same time, the backrest was adjusted so that the participant could easily flex and extend the knee. The shoulders, hips, and distal femur were then immobilised with instruments and a safety belt. The right knee range of motion was programmed to start from full extension to 90° of flexion (Sundby and Gorelick, 2014). During the test, the subjects were asked to hold onto the handles to keep their posture stable. The knee flexor muscles concentric (Kflex-con), eccentric (Kflex-ecc), knee extensor muscles concentric (Kex-con), and eccentric (Kex-ecc) PT at 60°s−1were measured.

Ankle strength was measured in the second trial. The participants were tested in a supine position on the dynamometer seat with the knees and hips fully extended. The right foot was strapped to the ankle adapter of the dynamometer and fixed with a safety belt. The positions of thighs, hips, and shoulders were fixed. The shaft of the dynamometer was aligned with the axis of the lateral malleolus. The ankle neutral position was programmed to 0°, and the movement range began from 15° dorsiflexion to 40° plantar flexion (Tsiokanos et al., 2002). The ankle dorsiflexor muscles concentric (Adors-con), eccentric (Adors-ecc), plantar flexor muscles concentric (Aplan-con), and eccentric (Aplan-ecc) PT at 60° s−1were measured.

Isometric Mid-Thigh Pull Test

The IMTP test was performed on two force platforms sampling at 1000 HZ (9290AA; Kistler, Winterthur, Switzerland) and a portable test rack (Figure 3B). This testing has consistently been shown to be highly reliable with ICCs ranging from 0.80 to 0.97 (Mcguigan, 2019). Before the test, the subject was instructed to place and mark the crossbar in the clean second pull position, which was defined as 140–150° of hip flexion and 125–145° of knee flexion (Chavda et al., 2020). All participants were asked to use an overhand grip while secured to the barbell with lifting straps to eliminate the influence of grip strength. Four warm-up trials were performed using 50, 70, 80, and 90% of their maximum effort with a rest period of 60 s (Beattie et al., 2017). Before the trial, subjects were instructed to “Push the ground fast and hard with maximum effort” to perform to their best, followed by completing three maximal effort IMTP tests to obtain 3 force-time curves of 2 s standing still and 5 s pulling with full force. The force-time characteristics with the highest PF was included in the statistical analysis. There was a 3-min interval between trials, and the subjects were greatly encouraged in each trial.

The highest force generated during IMTP is reported as PF expressed relative to body weight. In addition, the time-specific force (∆Force) at 50–350 ms from the initiation of the pull was collected. RFD was calculated using the following equation:

RFD=ΔForce/ΔTime(4)

RFD was applied to predetermined specific time intervals: 0–50 to 0–350 ms (RFD0–50, RFD0–100, RFD0–150, RFD0–200, RFD0–250, RFD0–300, and RFD0–350). These time intervals were selected based on the ground-contact times during various submaximal running speeds reported earlier (Gómez-Molina et al., 2017; García-Pinillos et al., 2019a).

Statistical Analyses

The Shapiro-Wilk test was used to test data normality. Values are expressed as mean ± SD. The correlations between the isokinetic lower-limb joint PT (knee and ankle concentric and eccentric PT), lower-limb stiffness (Kvert, Kleg, Kknee, and Kankle), and IMTP force-time characteristics (PF and specific time RFD at 0–50 to 0–350 ms) and RE (10 km·h−1) were analysed using Pearson’s correlation coefficients. Correlations were classified as “small” when r = 0.1–0.3, “moderate”, when r = 0.3–0.5, “large”, when r = 0.5–0.7, “very large”, when r = 0.7–0.9, and “extremely large”, when r = 0.9–1.0 (Hopkins et al., 2009). All reported p values were corrected using the Benjamini–Hochberg procedure, a method used to correct for false-discovery rates arising from statistical multiple comparisons (Benjamini and Hochberg, 1995). Statistical analyses were performed using SPSS (version 22.0; IBM Corp. Armonk, NY, United States) and the statistical programming language R (www.r-project.org). Statistical significance was set at an alpha level of p < 0.05, for Pearson’s correlation coefficients.

Results

All physical and physiological characteristics, lower-limb stiffness, isokinetic strength, and IMTP force-time characteristics values are shown in Table 1. The correlation coefficients between the neuromuscular indicators and RE are demonstrated in Table 2.

TABLE 2
www.frontiersin.org

TABLE 2. Correlations between neuromuscular characteristics and RE at 10 km·h−1 (n = 30).

For isokinetic strength testing, there were non-significant, small correlations between the knee (r = −0.070 to 0.132, p = 0.881–0.937) and ankle joint strength (r = −0.061 to 0.096, p = 0.881–0.955) with RE at 10 km·h−1.

With regard to lower-limb stiffness, we found that the Kvert (r = −0.449, p = 0.049) moderately correlated with RE at 10 km·h−1 (Figure 4). Meanwhile, non-significant, small correlations were found between Kleg (r = −0.100, p = 0.881), Kknee (r = −0.197, p = 0.594), Kankle (r = −0.394, p = 0.073), and RE after Benjamini–Hochberg adjustment.

FIGURE 4
www.frontiersin.org

FIGURE 4. Correlations of Kvert, IMTP force-time characteristics with 10 km·h−1 (A–G) running economy. Kvert, vertical stiffness; RFD, rate of force development; RFD0–50, RFD 0–50 ms; RFD0–100, RFD 0–100 ms; RFD0–150, RFD 0–150 ms; RFD0–200, RFD 0–200 ms; RFD0–250, RFD 0–250 ms; RFD0–300, RFD 0–300 ms; RFD0–350, RFD 0–350 ms; IMTP, isometric mid-thigh pull.

For the IMTP test, we found moderate to large relationships between the specific time RFD at 0–50 to 0–300 ms (r = −0.434 to −0.534, p < 0.05) and RE at 10 km·h−1 after false-discovery rate corrections (Figure 4). There were no significant relationships between the PF (r = 0.070, p = 0.881) and RE.

Discussion

This study aimed to investigate the relationship between isokinetic knee and ankle joint strength, lower-limb stiffness, isometric force-time characteristics, and RE in recreationally trained male runners. Our primary findings were that knee and ankle joint strength, and stiffness were not significantly correlated with RE. The Kvert, but not Kleg, was significantly associated with RE. Furthermore, the time-specific RFD at 0–50 to 0–300 ms were found to be significantly correlated with RE. Non-significant, small-sized relationships were found between IMTP PF and RE at 10 km·h−1.

Isokinetic Lower-Limb Joint Strength and Running Economy

The muscles around the knee and ankle joints play an important role in running. However, the results of this study reject our initial hypothesis and do not show any significant correlation between isokinetic knee and ankle muscle strength and RE. This finding is in line with that the studies of Sundby and Gorelick. (2014) and Andrade et al. (2021), who reported no correlation between concentric or eccentric knee flexion and extension PT and RE in distance runners. It is worth noting that the behavior of the knee/ankle flexor and extensor muscles during running was different from that in the isokinetic strength test. During the running support phase, muscle fascicles length changes are decoupled from tendon length changes (Farris and Sawicki, 2012). The knee extensor muscle fascicles contract quasi-isometrically while the tendon is stretched, optimising the storage and return of elastic strain energy in the tendon (Monte et al., 2020). Similarly, the ankle plantar flexor (i.e., gastrocnemius medialis) and dorsiflexor muscles fascicles (i.e., tibialis anterior) are co-activated and contracted isometrically to develop force, and absorb energy through stretching of the tendinous tissue at foot contact and returning it to the body in subsequent strides (Tam et al., 2018; Maharaj et al., 2019; Monte et al., 2020). For the isokinetic strength test, the knee flexor/extensor and ankle plantar flexor/dorsiflexor muscles are purely concentric or eccentric, and the storage of elastic strain energy is not possible because the tendon does not undergo a preceding stretch (Roberts, 2002).

Apart from the intrinsic muscle-tendon behavior, the isokinetic test performs an open chain movement, while running is a sport that combines closed and open chain movements (Andrade et al., 2021). Stensdotter et al. (2003) found differences in the EMG onset and amplitude of the different knee extensor muscles during the open and closed kinetic chain test. For example, during closed chain knee extension, the quadriceps muscles onset was almost simultaneous, whereas in the open chain the rectus femoris was activated first and the vastus medialis obliquus was activated last. Li et al. (2019) reported that isokinetic leg press tests (closed chain movement) were largely correlated with RE, indicating that the effect of movement pattern choosing in test on predicting RE. Worth to note that the angular velocity generated during running is much greater than the present isokinetic strength test (60°s−1). The maximum angular velocity generated by the kinematics of the joints (knee and ankle) during the human running gait cycle is between 400 and 600°s−1 (Grimmer and Seyfarth, 2014). Meanwhile, the absence of angular acceleration in the isokinetic test is also very distinct from the functional joint movements during running (Andrade et al., 2021). Therefore, the muscle contraction behavior, movement pattern and angular velocity differences may limit the association between isokinetic strength testing and practical running.

Lower-Limb Stiffness and Running Economy

Kvert and Kleg as global stiffnesses describe the ability of the entire lower limb to reduce energy consumption and to utilise elastic energy during vertical and horizontal movements (Dalleau et al., 1998; Heise and Martin, 2001; Struzik et al., 2021). In our study, the Kleg values (13.04 kN·m−1) and Kvert values (25.07 kN·m−1) reported are similar to those of previous studies at similar speeds among recreational runners (Morin et al., 2005; García-Pinillos et al., 2019b). We found that Kvert, but not Kleg, was significantly negatively associated with RE at 10 km·h−1. This finding is consistent with Heise and Martrin. (1998) and partially supports our hypothesis. Kvert is mainly determined by the COM vertical displacement, a higher Kvert allowing for more economical running tasks (with less COM vertical displacement) and with superior performance through greater potential elastic energy return from the tendon structure (Struzik et al., 2021). Folland et al. (2017) reported that vertical displacement of COM (the standard measure of vertical oscillation) explained 28% of the inter-individual variation in RE at 10 and 12 km·h−1 across diverse competitive runners. Therefore, a higher Kvert may be an important contributor to regulate RE by reducing vertical oscillation during the stance phase.

Li et al. (2019) recently reported that Kleg has a large correlation with RE at 12–16 km·h−1, which indicates that Kleg is an essential determinant of RE in well-trained runners. However, this is not consistent with our findings. The difference in Kleg may be due to the experience of endurance running training. For runners, their energy consumption has a curvilinear U-shaped relationship with Kleg rather than a linear one, and there is an optimal Kleg to minimise energy costs (Hunter and Smith, 2007; Moore et al., 2019). Highly trained runners are able to quickly adapt to the optimal regulation of Kleg through prolonged running experience, such as changing gait characteristics (e.g., stride length, stride frequency and contact time) to improving RE (Hunter and Smith, 2007; Moore et al., 2019). In contrast, no optimisation of the above Kleg or gait characteristics was found in recreational runners, who were further away from their optimal stiffness than trained runners (De Ruiter et al., 2014; Bitchell et al., 2019). Moreover, the percentage contribution of elastic strain energy to positive work was shown to be reduced at lower running speeds (Lai et al., 2014). Therefore, it may not be the most advantageous condition for our recreational subjects to improve RE through the reuse of elastic energy in the lower limbs at slow running speeds (i.e., 10 km·h−1).

Joint stiffness reflects the elastic properties of the musculotendinous tissue surrounding individual joint (Kuitunen et al., 2002). Tam et al. (2018) reported that knee and ankle stiffness were significantly associated with RE in trained runners. However, contrary to our initial hypothesis, our results do not exhibit this significant relationship in recreational runners. Many studies indicate that changes in joint stiffness partly be related to foot strike pattern during landing phase (Butler et al., 2003). Compared to a rearfoot strike pattern, knee stiffness was greater and ankle stiffness was lower in the forefoot strike pattern (Hamill et al., 2014). In our study, all subjects were rearfoot strike pattern, whereas Tam et al. (2018) did not report the foot strike pattern of their subjects landed. Considering that the highly trained runners may possibly use the midfoot or forefoot landing pattern (Hasegawa et al., 2007; Mo et al., 2021), thus resulting in the different values of joint stiffness. In addition, the running speed affects joint stiffness. Kuitunen et al. (2002) reported increasing running velocity was associated with increases in knee stiffness. Specifically, Tam et al. (2018) determined KKnee and Kankle using the running speed of 3.28 m·s−1 while our study using of 2.78 m·s−1. Future studies should ascertain if ankle and knee stiffness would be influence at different running speeds.

Although we did not observe Kleg, Kknee and Kankle association with RE at 10 km·h−1, previous studies indicate the significant contributions to running performance especially in well-trained runners, as elastic energy accumulation is an important factor in reducing energy costs (Struzik et al., 2021). Therefore, strategies to improve RE through the reuse of elastic energy may be more advantageous for recreational runners at higher running speeds. At same time, lower-limb muscles also need to shorten at faster velocities and recruit more motor units to generate the necessary high forces in a very short contact time (Fletcher and MacIntosh, 2017).

Isometric Force-Time Characteristics and Running Economy

IMTP is performed to measure the force-time characteristics (i.e., PF and RFD) of various sports (Mcguigan, 2019), and allows the hip, knee, and ankle joints to be at a relatively biomechanical angle during running (Lum et al., 2020). Lum et al. (2020) recently reported that the IMTP time-specific RFD (0–100 to 0–200 ms) significantly correlate with RE at submaximal speeds, indicating that runners with higher force production exhibit better RE. However, the authors did not evaluate forces above 200 ms in their study. Our data (252 ms), as well as previous studies (Gómez-Molina et al., 2017; García-Pinillos et al., 2019a), suggest that ground contact time is greater than 200 ms at running speeds of 10 km·h−1. Therefore, it is essential to evaluate forces above 200 ms for recreational runners. As hypothesised, RFD within the running contact time was significantly correlated with RE. In the present study, we found that time-specific RFD (0–50 to 0–300 ms) with 10 km·h−1 RE were moderate to largely inversely correlated. The reasons for this phenomenon can be explained as follows. Firstly, from the behavior of the lower limb muscle during the support period, a larger RFD may allow the lower-limb muscle to activate rapidly or generate higher forces during a shorter contact time (Lum et al., 2020), promoting favourable muscles conditions (quasi-isometric contraction and near-optimal length), thus optimising the storage and return of elastic energy and reducing the extra work done by the muscles (Lai et al., 2014; Fletcher and MacIntosh, 2017; Bohm et al., 2019; Monte et al., 2020). Secondly, the greater RFD theoretically enables the runners to rapidly push off the ground, thereby decreasing ground contact and muscle contraction time during the stance phase, which results in a more rapid transition of the running gait from the braking to the propulsive phase (Lum et al., 2020). This may potentially reduce vertical oscillation and metabolic demands against gravity (Saunders et al., 2004; Lima and Blagrove, 2020).

Muscles contracting to produce high force theoretically allow the runners to maintain a given speed or perform every running action at a relatively lower intensity (Fletcher and MacIntosh, 2017). However, contrary to our initial hypothesis, we did not find a significant relationship between IMTP PF and RE. This finding is consistent with those of Lum et al. (2020), who reported that PF is not significantly related to RE or running performance in male runners. This can be attributed to the fact that the muscular behavioral processes that drive energy savings need to be timed precisely during the support phase of running (Barnes and Kilding, 2015), whereas the IMTP test require a relatively long time for runners to reach peak force (Chavda et al., 2020). Therefore, compared to the time-specific RFD, the PF is not specific in predicting the energetic cost of running.

This study has certain limitations that must be acknowledged. Firstly, the current findings indicate that the neuromuscular characteristics mentioned above is an acceptable predictor of RE in recreational male runners. Further studies are needed to investigate the potential relationships among elite runners and female runners’ cohort. Secondly, we only test isokinetic joint strength at 60°s−1 velocity, as kinematics of the joints (knee and ankle) produces maximum angular velocities between 400 and 600°s−1 during the human running gait cycle, it is valuable to investigate the relationship between joint strength and RE at higher angular velocities. Finally, a large correlation does not imply a cause-effect relationship. Hence, in the future, studies with longitudinal study design should be performed to confirm the effects of isometric force-time characteristics, lower-limb stiffness and their manipulation through intervention on running economy.

Conclusion

In summary, we did not find significant correlations between knee and ankle joint strength, stiffness, and RE. However, for the whole lower-limb neuromuscular characteristics, IMTP time-specific RFD at 0–50 to 0–300 ms and vertical stiffness significantly correlated with RE at 10 km·h−1 in recreational runners.

Although isokinetic testing is the gold standard for measuring muscle strength (Kambič et al., 2020), given the specific movement of running, the joint strength at higher angular velocity (i.e., 240°s−1) and multi-joint closed-chain movement (i.e., leg press) need be considered when measuring the strength capacity of runners. In addition, since the lower limb muscles preform isometric contractions during the running contact phase, the isometric contraction mode should be utilized in the isokinetic strength test. When evaluating lower limb stiffness, the global stiffness (i.e., vertical or leg stiffness) takes precedence over local stiffness (i.e., ankle or knee joint stiffness) to evaluate runners’ elastic energy utilization. Finally, the results also suggest that runner should focus on ability to rapidly generate force which correspond to the foot contact time (<300 ms).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of Shanghai University of Sport, China (ID number: 102772020RT113). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

QZ and SC conceived and designed the experiments, performed the experiments, analyzed the data, contributed reagents, materials, and analysis tools, prepared figures and/or tables, authored or reviewed drafts of the manuscript, and approved the final draft. GN conceived and designed the experiments, and authored or reviewed drafts of the manuscript. YS performed the experiments, and prepared figures and/or tables. FL conceived and designed the experiments, and approved the final draft. All authors contributed to the article and approved the submitted version.

Funding

This study was sponsored by Shanghai Sailing Program (Fund number: 20YF1446500), Shanghai Key Lab of Human Performance, Shanghai University of Sport (No. 11DZ2261100).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors would like to thank all the runners who were involved in this research and all the members and researchers at the Chinese Marathon College.

References

Andrade M. S., Silva W. A., de Lira C. A. B., Mascarin N. C., Vancini R. L., Nikolaidis P. T., et al. (2021). Isokinetic Muscular Strength and Aerobic Physical Fitness in Recreational Long-Distance Runners: a Cross-Sectional Study. J. Strength Cond. Res. Adv. online Publ. 36, e73–e80. doi:10.1519/JSC.0000000000003837

CrossRef Full Text | Google Scholar

Barnes K. R., Kilding A. E. (2015). Running Economy: Measurement, Norms, and Determining Factors. Sports Med. - Open 1, 8. doi:10.1186/s40798-015-0007-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Beattie K., Carson B. P., Lyons M., Kenny I. C. (2017). The Effect of Maximal- and Explosive-Strength Training on Performance Indicators in Cyclists. Int. J. Sports. Physiol. Perform. 12, 470–480. doi:10.1123/ijspp.2016-0015

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini Y., Hochberg Y. (1995). Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Besser M. P., Kmieczak K., Schwartz L., Snyderman M., Wasko J., Selby- Silverstein L. (1999). Representation of Temporal Spatial Gait Parameters Using Means in Adults without Impairment. Gait Posture 9, 113.

Google Scholar

Bitchell C. L., McCarthy-Ryan M., Goom T., Moore I. S. (2019). Spring-mass Characteristics during Human Locomotion: Running Experience and Physiological Considerations of Blood Lactate Accumulation. Eur. J. Sport Sci. 19, 1328–1335. doi:10.1080/17461391.2019.1609095

PubMed Abstract | CrossRef Full Text | Google Scholar

Blagrove R. C., Howatson G., Hayes P. R. (2018). Effects of Strength Training on the Physiological Determinants of Middle- and Long-Distance Running Performance: a Systematic Review. Sports Med. 48, 1117–1149. doi:10.1007/s40279-017-0835-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohm S., Marzilger R., Mersmann F., Santuz A., Arampatzis A. (2018). Operating Length and Velocity of Human Vastus Lateralis Muscle during Walking and Running. Sci. Rep. 8, 5066. doi:10.1038/s41598-018-23376-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohm S., Mersmann F., Santuz A., Arampatzis A. (2021). Enthalpy Efficiency of the Soleus Muscle Contributes to Improvements in Running Economy. Proc. R. Soc. B 288, 20202784. doi:10.1098/rspb.2020.2784

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohm S., Mersmann F., Santuz A., Arampatzis A. (2019). The Force-Length-Velocity Potential of the Human Soleus Muscle Is Related to the Energetic Cost of Running. Proc. R. Soc. B 286, 20192560. doi:10.1098/rspb.2019.2560

PubMed Abstract | CrossRef Full Text | Google Scholar

Brady C. J., Harrison A. J., Flanagan E. P., Haff G. G., Comyns T. M. (2020). The Relationship between Isometric Strength and Sprint Acceleration in Sprinters. Int. J. Sports. Physiol. Perform. 15, 38–45. doi:10.1123/ijspp.2019-0151

CrossRef Full Text | Google Scholar

Butler R. J., Crowell H. P., Davis I. M. (2003). Lower Extremity Stiffness: Implications for Performance and Injury. Clin. Biomech. 18, 511–517. doi:10.1016/s0268-0033(03)00071-8

CrossRef Full Text | Google Scholar

Chavda S., Turner A. N., Comfort P., Haff G. G., Williams S., Bishop C., et al. (2020). A Practical Guide to Analyzing the Force-Time Curve of Isometric Tasks in Excel. Strength Cond. J. 42, 26–37. doi:10.1519/SSC.0000000000000507

CrossRef Full Text | Google Scholar

Cohen J. (1988). Statistical Power Analysis for the Behavioral Sciences. 2nd ed. New York: L. Erlbaum Associates.

Google Scholar

Conley D. L., Krahenbuhl G. S. (1980). Running Economy and Distance Running Performance of Highly Trained Athletes. Med. Sci. Sports Exerc. 12, 357–360. doi:10.1249/00005768-198012050-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalleau G., Belli A., Bourdin M., Lacour J.-R. (1998). The Spring-Mass Model and the Energy Cost of Treadmill Running. Eur. J. Appl. Physiology 77, 257–263. doi:10.1007/s004210050330

CrossRef Full Text | Google Scholar

De Ruiter C. J., Verdijk P. W. L., Werker W., Zuidema M. J., de Haan A. (2014). Stride Frequency in Relation to Oxygen Consumption in Experienced and Novice Runners. Eur. J. Sport Sci. 14, 251–258. doi:10.1080/17461391.2013.783627

PubMed Abstract | CrossRef Full Text | Google Scholar

Farley C. T., González O. (1996). Leg Stiffness and Stride Frequency in Human Running. J. Biomechanics 29, 181–186. doi:10.1016/0021-9290(95)00029-1

CrossRef Full Text | Google Scholar

Farris D. J., Sawicki G. S. (2012). Human Medial Gastrocnemius Force-Velocity Behavior Shifts with Locomotion Speed and Gait. Proc. Natl. Acad. Sci. U.S.A. 109, 977–982. doi:10.1073/pnas.1107972109

PubMed Abstract | CrossRef Full Text | Google Scholar

Faul F., Erdfelder E., Buchner A., Lang A.-G. (2009). Statistical Power Analyses Using G*Power 3.1: Tests for Correlation and Regression Analyses. Behav. Res. Methods 41, 1149–1160. doi:10.3758/BRM.41.4.1149

PubMed Abstract | CrossRef Full Text | Google Scholar

Fletcher J. R., MacIntosh B. R. (2017). Running Economy from a Muscle Energetics Perspective. Front. Physiol. 8, 433. doi:10.3389/fphys.2017.00433

PubMed Abstract | CrossRef Full Text | Google Scholar

Folland J. P., Allen S. J., Black M. I., Handsaker J. C., Forrester S. E. (2017). Running Technique Is an Important Component of Running Economy and Performance. Med. Sci. Sports Exerc. 49, 1412–1423. doi:10.1249/MSS.0000000000001245

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Pinillos F., García-Ramos A., Ramírez-Campillo R., Latorre-Román P. Á., Roche-Seruendo L. E. (2019a). How Do Spatiotemporal Parameters and Lower-Body Stiffness Change with Increased Running Velocity? A Comparison between Novice and Elite Level Runners. J. Hum. Kinet. 70, 25–38. doi:10.2478/hukin-2019-0036

CrossRef Full Text | Google Scholar

García-Pinillos F., Latorre-Román P. Á., Ramírez-Campillo R., Párraga-Montilla J. A., Roche-Seruendo L. E. (2019b). How Does the Slope Gradient Affect Spatiotemporal Parameters during Running? Influence of Athletic Level and Vertical and Leg Stiffness. Gait Posture 68, 72–77. doi:10.1016/j.gaitpost.2018.11.008

CrossRef Full Text | Google Scholar

Gómez-Molina J., Ogueta-Alday A., Stickley C., Cámara J., Cabrejas-Ugartondo J., García-López J. (2017). Differences in Spatiotemporal Parameters Between Trained Runners and Untrained Participants. J. Res. 31, 2169–2175. doi:10.1519/JSC.0000000000001679

CrossRef Full Text | Google Scholar

Gonosova Z., Linduska P., Bizovska L., Svoboda Z. (2018). Reliability of Ankle-Foot Complex Isokinetic Strength Assessment Using the Isomed 2000 Dynamometer. Medicina 54, 43. doi:10.3390/medicina54030043

PubMed Abstract | CrossRef Full Text | Google Scholar

Grimmer M., Seyfarth A. (2014). Mimicking Human-like Leg Function in Prosthetic Limbs. Dordrecht, Netherlands: Springer Netherlands, 111.

Google Scholar

Hamill J., Gruber A. H., Derrick T. R. (2014). Lower Extremity Joint Stiffness Characteristics during Running with Different Footfall Patterns. Eur. J. Sport Sci. 14, 130–136. doi:10.1080/17461391.2012.728249

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamner S. R., Delp S. L. (2013). Muscle Contributions to Fore-Aft and Vertical Body Mass Center Accelerations over a Range of Running Speeds. J. Biomechanics 46, 780–787. doi:10.1016/j.jbiomech.2012.11.024

CrossRef Full Text | Google Scholar

Hamner S. R., Seth A., Delp S. L. (2010). Muscle Contributions to Propulsion and Support during Running. J. Biomechanics 43, 2709–2716. doi:10.1016/j.jbiomech.2010.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasegawa H., Yamauchi T., Kraemer W. J. (2007). Foot Strike Patterns of Runners at the 15-km Point during an Elite-Level Half Marathon. J. Strength Cond. Res. 21, 888–893. doi:10.1519/R-22096.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Heise G. D., Martin P. E. (2001). Are Variations in Running Economy in Humans Associated with Ground Reaction Force Characteristics? Eur. J. Appl. Physiology 84, 438–442. doi:10.1007/s004210100394

CrossRef Full Text | Google Scholar

Heise G. D., Martin P. E. (1998). “Leg Spring” Characteristics and the Aerobic Demand of Running. Med. Sci. Sports Exerc. 30, 750–754. doi:10.1097/00005768-199805000-00017

PubMed Abstract | CrossRef Full Text | Google Scholar

Hopkins W. G., Marshall S. W., Batterham A. M., Hanin J. (2009). Progressive Statistics for Studies in Sports Medicine and Exercise Science. Med. Sci. Sports Exerc. 41, 3–12. doi:10.1249/MSS.0b013e31818cb278

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunter I., Smith G. A. (2007). Preferred and Optimal Stride Frequency, Stiffness and Economy: Changes with Fatigue during a 1-h High-Intensity Run. Eur. J. Appl. Physiol. 100, 653–661. doi:10.1007/s00421-007-0456-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kambič T., Lainščak M., Hadžić V. (2020). Reproducibility of Isokinetic Knee Testing Using the Novel Isokinetic SMM iMoment Dynamometer. PLoS One 15, e0237842. doi:10.1371/journal.pone.0237842

PubMed Abstract | CrossRef Full Text | Google Scholar

Komi P. V. (1990). Relevance of In Vivo Force Measurements to Human Biomechanics. J. Biomechanics 23 (Suppl. 1), 23–34. doi:10.1016/0021-9290(90)90038-5

CrossRef Full Text | Google Scholar

Kuitunen S., Komi P. V., Kyröläinen H. (2002). Knee and Ankle Joint Stiffness in Sprint Running. Med. Sci. Sports Exerc. 34, 166–173. doi:10.1097/00005768-200201000-00025

PubMed Abstract | CrossRef Full Text | Google Scholar

Kulmala J.-P., Kosonen J., Nurminen J., Avela J. (2018). Running in Highly Cushioned Shoes Increases Leg Stiffness and Amplifies Impact Loading. Sci. Rep. 8, 17496. doi:10.1038/s41598-018-35980-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kyr?L?Inen H., Belli A., Komi P. V. (2001). Biomechanical Factors Affecting Running Economy. Med. Sci. Sports Exerc. 33, 1330–1337. doi:10.1097/00005768-200108000-00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai A., Schache A. G., Lin Y.-C., Pandy M. G. (2014). Tendon Elastic Strain Energy in the Human Ankle Plantar-Flexors and its Role with Increased Running Speed. J. Exp. Biol. 217, 3159–3168. doi:10.1242/jeb.100826

PubMed Abstract | CrossRef Full Text | Google Scholar

Li F., Nassis G. P., Shi Y., Han G., Zhang X., Gao B., et al. (2021). Concurrent Complex and Endurance Training for Recreational Marathon Runners: Effects on Neuromuscular and Running Performance. Eur. J. Sport Sci. 21, 1243–1253. doi:10.1080/17461391.2020.1829080

PubMed Abstract | CrossRef Full Text | Google Scholar

Li F., Newton R. U., Shi Y., Sutton D., Ding H. (2019). Correlation of Eccentric Strength, Reactive Strength, and Leg Stiffness with Running Economy in Well-Trained Distance Runners. J. Strength Cond. Res. 35, 1491–1499. doi:10.1519/JSC.0000000000003446

PubMed Abstract | CrossRef Full Text | Google Scholar

Lima L. C. R., Blagrove R. (2020). Infographic. Strength Training-Induced Adaptations Associated with Improved Running Economy: Potential Mechanisms and Training Recommendations. Br. J. Sports Med. 54, 302–303. doi:10.1136/bjsports-2019-100840

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu M. Q., Anderson F. C., Schwartz M. H., Delp S. L. (2008). Muscle Contributions to Support and Progression over a Range of Walking Speeds. J. Biomechanics 41, 3243–3252. doi:10.1016/j.jbiomech.2008.07.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorimer A. V., Keogh J. W. L., Hume P. A. (2018). Using Stiffness to Assess Injury Risk: Comparison of Methods for Quantifying Stiffness and Their Reliability in Triathletes. PeerJ 6, e5845. doi:10.7717/peerj.5845

PubMed Abstract | CrossRef Full Text | Google Scholar

Lum D., Chua K., Rashid Aziz A. (2020). Isometric Mid-thigh Pull Force-Time Characteristics: a Good Indicator of Running Performance. J. Trainology. 9, 54–59. doi:10.17338/trainology.9.2_54

CrossRef Full Text | Google Scholar

Luna N. M., Alonso A. C., Brech G. C., Mochizuki L., Nakano E. Y., Greve J. M. (2012). Isokinetic Analysis of Ankle and Ground Reaction Forces in Runners and Triathletes. Clin. (Sao Paulo) 67, 1023–1028. doi:10.6061/clinics/2012(09)07

PubMed Abstract | CrossRef Full Text | Google Scholar

Maharaj J. N., Cresswell A. G., Lichtwark G. A. (2019). Tibialis Anterior Tendinous Tissue Plays a Key Role in Energy Absorption during Human Walking. J. Exp. Biol. 222, jeb191247. doi:10.1242/jeb.191247

PubMed Abstract | CrossRef Full Text | Google Scholar

Maloney S. J., Fletcher I. M. (2021). Lower Limb Stiffness Testing in Athletic Performance: a Critical Review. Sports Biomech. 20, 109–130. doi:10.1080/14763141.2018.1460395

PubMed Abstract | CrossRef Full Text | Google Scholar

Man H. S., Lam W. K., Lee J., Capio C. M., Leung A. K. L. (2016). Is Passive Metatarsophalangeal Joint Stiffness Related to Leg Stiffness, Vertical Stiffness and Running Economy during Sub-maximal Running? Gait Posture 49, 303–308. doi:10.1016/j.gaitpost.2016.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Mcguigan M. (2019). Testing Isometric Strength Testing and Evaluation of Strength and Power. Routledge Taylor Francis Group 52, 57–76. doi:10.4324/9780429028182-4

CrossRef Full Text | Google Scholar

Mo S., Huang M., Ng L., Cheung R. T. H. (2021). Footstrike Angle Cut-Off Values to Classify Footstrike Pattern in Runners. Res. Sports Med. 18, 1–11. doi:10.1080/15438627.2021.1954516

CrossRef Full Text | Google Scholar

Monte A., Baltzopoulos V., Maganaris C. N., Zamparo P. (2020). Gastrocnemius Medialis and Vastus Lateralis In Vivo Muscle‐tendon Behavior during Running at Increasing Speeds. Scand. J. Med. Sci. Sports 30, 1163–1176. doi:10.1111/sms.13662

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore I. S., Ashford K. J., Cross C., Hope J., Jones H. S. R., McCarthy-Ryan M. (2019). Humans Optimize Ground Contact Time and Leg Stiffness to Minimize the Metabolic Cost of Running. Front. Sports Act. Living 1, 53. doi:10.3389/fspor.2019.00053

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan D. W., Baldini F. D., Martin P. E., Kohrt W. M. (1989). Ten Kilometer Performance and Predicted Velocity at V02max Among Well-Trained Male Runners. Med. Sci. Sports Exerc. 21, 78–83. doi:10.1249/00005768-198902000-00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Morin J.-B., Dalleau G., Kyröläinen H., Jeannin T., Belli A. (2005). A Simple Method for Measuring Stiffness during Running. J. Appl. Biomech. 21, 167–180. doi:10.1123/jab.21.2.167

PubMed Abstract | CrossRef Full Text | Google Scholar

Novacheck T. F. (1998). The Biomechanics of Running. Gait Posture 7, 77–95. doi:10.1016/s0966-6362(97)00038-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pappas P., Paradisis G., Tsolakis C., Smirniotou A., Morin J.-B. (2014). Reliabilities of Leg and Vertical Stiffness during Treadmill Running. Sports Biomech. 13, 391–399. doi:10.1080/14763141.2014.981853

PubMed Abstract | CrossRef Full Text | Google Scholar

Piacentini M. F., De Ioannon G., Comotto S., Spedicato A., Vernillo G., La Torre A. (2013). Concurrent Strength and Endurance Training Effects on Running Economy in Master Endurance Runners. J. Strength Cond. Res. 27, 2295–2303. doi:10.1519/JSC.0b013e3182794485

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts T. J. (2002). The Integrated Function of Muscles and Tendons during Locomotion. Comp. Biochem. Physiology Part A Mol. Integr. Physiology 133, 1087–1099. doi:10.1016/s1095-6433(02)00244-1

CrossRef Full Text | Google Scholar

Saunders P. U., Pyne D. B., Telford R. D., Hawley J. A. (2004). Factors Affecting Running Economy in Trained Distance Runners. Sports Med. 34, 465–485. doi:10.2165/00007256-200434070-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva W. A., de Lira C. A. B., Vancini R. L., Andrade M. S. (2018). Hip Muscular Strength Balance Is Associated with Running Economy in Recreationally-Trained Endurance Runners. PeerJ 6, e5219. doi:10.7717/peerj.5219

PubMed Abstract | CrossRef Full Text | Google Scholar

Stensdotter A.-K., Hodges P. W., Mellor R., Sundelin G., H?Ger-Ross C. (2003). Quadriceps Activation in Closed and in Open Kinetic Chain Exercise. Med. Sci. Sports Exerc. 35, 2043–2047. doi:10.1249/01.MSS.0000099107.0370410.1249/01.mss.0000099107.03704.ae

PubMed Abstract | CrossRef Full Text | Google Scholar

Struzik A., Karamanidis K., Lorimer A., Keogh J. W. L., Gajewski J. (2021). Application of Leg, Vertical, and Joint Stiffness in Running Performance: a Literature Overview. Appl. Bionics Biomech. 2021, 9914278. doi:10.1155/2021/9914278

PubMed Abstract | CrossRef Full Text | Google Scholar

Sundby Ø. H., Gorelick M. L. S. (2014). Relationship Between Functional Hamstring. J. Strength Cond. Res. 28, 2214–2227. doi:10.1519/JSC.0000000000000376

PubMed Abstract | CrossRef Full Text | Google Scholar

Tam N., Tucker R., Santos-Concejero J., Prins D., Lamberts R. P. (2019). Running Economy: Neuromuscular and Joint-Stiffness Contributions in Trained Runners. Int. J. Sports. Physiol. Perform. 14, 16–22. doi:10.1123/ijspp.2018-0151

CrossRef Full Text | Google Scholar

Trowell D., Vicenzino B., Saunders N., Fox A., Bonacci J. (2020). Effect of Strength Training on Biomechanical and Neuromuscular Variables in Distance Runners: A Systematic Review and Meta-Analysis. Sports Med. 50, 133–150. doi:10.1007/s40279-019-01184-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsiokanos A., Kellis E., Jamurtas A., Kellis S., Κλλη. Σ. (2002). The Relationship between Jumping Performance and Isokinetic Strength of Hip and Knee Extensors and Ankle Plantar Flexors. Ies 10, 107–115. doi:10.1109/42.98123210.3233/ies-2002-0092

CrossRef Full Text | Google Scholar

Wang R., Hoffman J. R., Tanigawa S., Miramonti A. A., La Monica M. B., Beyer K. S., et al. (2016). Isometric Mid-thigh Pull Correlates with Strength, Sprint, and Agility Performance in Collegiate Rugby Union Players. J. Strength Cond. Res. 30, 3051–3056. doi:10.1519/JSC.0000000000001416

PubMed Abstract | CrossRef Full Text | Google Scholar

Westblad P., Svedenhag J., Rolf C. (1996). The Validity of Isokinetic Knee Extensor Endurance Measurements with Reference to Treadmill Running Capacities. Int. J. Sports Med. 17, 134–139. doi:10.1055/s-2007-972821

PubMed Abstract | CrossRef Full Text | Google Scholar

Zapparoli F. Y., Riberto M. (2017). Isokinetic Evaluation of the Hip Flexor and Extensor Muscles: A Systematic Review. J. Sport Rehabil. 26, 556–566. doi:10.1123/jsr.2016-0036

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: neuromuscular characteristics, isometric mid-thigh pull, running performance, eccentric strength, stiffness, recreational runner

Citation: Zhang Q, Nassis GP, Chen S, Shi Y and Li F (2022) Not Lower-Limb Joint Strength and Stiffness but Vertical Stiffness and Isometric Force-Time Characteristics Correlate With Running Economy in Recreational Male Runners. Front. Physiol. 13:940761. doi: 10.3389/fphys.2022.940761

Received: 10 May 2022; Accepted: 13 June 2022;
Published: 28 June 2022.

Edited by:

Tobias Siebert, University of Stuttgart, Germany

Reviewed by:

Thomas Muehlbauer, University of Duisburg-Essen, Germany
Danny Lum, Sport Singapore, Singapore

Copyright © 2022 Zhang, Nassis, Chen, Shi and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Fei Li, lifei@sus.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.