Abstract
Purpose
This study aimed to evaluate the acute changes in the knee extensors maximum voluntary isometric contraction force (MVIC), rate of force development (RFD), and rate of EMG rise (RER) following a bout of downhill running.
Methods
MVIC and RFD at 0–50, 50–100, 100–200, and 0–200 ms were determined in thirteen men (22 ± 2 yr) before and after 30 min of downhill running (speed: 10 km h−1; slope: − 20%). Vastus lateralis maximum EMG (EMGmax) and RER at 0–30, 0–50, and 0–75 ms were also recorded.
Results
MVIC, RFD0–200, and EMGmax decreased by ~ 25% [Cohen’s d = − 1.09 (95% confidence interval: − 1.88/− 0.24)], ~ 15% [d = − 0.50 (− 1.26/0.30)], and ~ 22% [d = − 0.37 (− 1.13/0.42)] (all P < 0.05), respectively. RFD100–200 was also reduced [− 25%; d = − 0.70 (− 1.47/0.11); P < 0.001]. No change was observed at 0–50 ms and 50–100 ms (P ≥ 0.05). RER values were similar at each time interval (all P > 0.05).
Conclusion
Downhill running impairs the muscle capacity to produce maximum force and the overall ability to rapidly develop force. No change was observed for the early phase of the RFD and the absolute RER, suggesting no alterations in the neural mechanisms underlying RFD. RFD100–200 reduction suggests that impairments in the rapid force-generating capacity are located within the skeletal muscle, likely due to a reduction in muscle–tendon stiffness and/or impairments in the muscle contractile apparatus. These findings may help explain evidence of neuromuscular alterations in trail runners and following prolonged duration races wherein cumulative eccentric loading is high.
Similar content being viewed by others
Introduction
Trail running is an endurance-oriented discipline consisting of prolonged running on rough terrain with an alternation of uphill and downhill segments (Scheer et al. 2020). The popularity of these events has seen an exponential increase worldwide (Hoffman et al. 2010), and the combination of long distance and large sections of downhill running has been hypothesized to induce deleterious alterations to the performance of the neuromuscular system (Millet et al. 2011; Saugy et al. 2013; Temesi et al. 2014; Vernillo et al. 2019). Particularly during downhill running, the considerable role of the knee extensor muscles in energy absorption and dissipation (Khassetarash et al. 2020) results in high force eccentric muscle actions (Vernillo et al. 2017, 2020a), inducing muscle fatigability (Giandolini et al. 2016b; Bontemps et al. 2020). Accordingly, the research investigating the trail running performance identifies the downhill segments as particularly deleterious for the muscle function (Martin et al. 2004; Giandolini et al. 2016a; Ehrström et al. 2018).
Endurance running, both on roads and on trails, involves the repetitive activation of large muscle masses of the lower limbs for prolonged periods of time (Novacheck 1998). In addition, each step requires the ability to rapidly develop force (Nicol et al. 2006), a key aspect for trail running where technical segments and elevation changes could require a rapid muscle force generation. Much of the research investigating the influence of downhill running-induced fatigability has focused on the decline in maximum voluntary isometric contraction force (MVIC) (Giandolini et al. 2016b; Bontemps et al. 2020). Nonetheless, the capacity to develop a high amount of force in a very short amount of time is considered functionally more relevant than MVIC (Maffiuletti et al. 2016) and was recently proposed as an indicator of neuromuscular fatigue (D’Emanuele et al. 2021). Such a capacity is defined as the rate of force development (RFD) and can be assessed by measuring the slope of the force–time curve during short (~ 1 s) isometric impulsive contractions over a window of 200 ms, as suggested by Maffiuletti et al. (2016). Two phases, the early and late phase, have been identified to better understand all the neural and structural factors that could contribute to the RFD (Maffiuletti et al. 2016). Specifically, the early phase (the firsts 50 ms from the onset of the contraction) has been associated with the initial motor unit recruitment and firing rates (Klass et al. 2008), whereas maximum strength and other structural factors, such as muscle contractile properties (Aagaard et al. 2002) and muscle architecture (Coratella et al. 2020), influence the late phase (last 100 ms of the 200 ms time-window).
Accordingly, the analysis of RFD may be used to detect the fatigue-induced changes after an exercise, and an examination of the early and late phases may reflect the neural and structural fatigue-induced changes, respectively (Aagaard et al. 2002; Boccia et al. 2017). Furthermore, the RFD analysis can be complimented by the analysis of the rate of surface electromyography (EMG) signal rise (RER) to further detail the contribution of neural factors to RFD as a non-invasive surrogate of the capacity to deliver high-frequency motor unit bursts (i.e., lower maximal discharge frequency and incidence of doublet discharge during fast contractions) (Klass et al. 2008; Del Vecchio et al. 2019). Thus far, it has been shown that the RFD is impaired after level running events [e.g., half-marathon (Boccia et al. 2018)], whereas a limited number of studies have investigated the ability of rapidly develop force after downhill running. Specifically, Lima et al. (2020) examined the changes in RFD in the knee extensors immediately after a downhill run, reporting a ~ 20% decrease in peak RFD following a 30-min downhill running at 10 km h−1 and − 15% incline. However, neither the different phases of the force–time curve of RFD nor the RER were investigated, limiting the understanding of possible factors negatively affecting RFD following downhill running.
Therefore, the present study aimed to investigate the changes in the knee extensors RFD after a downhill run. If an effect of downhill running was observed on maximal force and RFD (200 ms window), the various phases of RFD, accompanied by the EMG rise, were investigated to determine the etiology (neural vs. muscular) of the fatigue-induced changes on the functional capacity of the neuromuscular system to develop a high amount of force in a very short period of time.
Methods
Participants
Thirteen healthy men (age: 22 ± 2 yr, height: 179 ± 7 cm, body mass: 76 ± 9 kg) volunteered to participate in the present study. Using sample size calculations provided by G*Power, the sample size was calculated a priori to show an effect size of around 1.2 [Cohen’s d (Cohen 1988)] for the effect of downhill running on the knee extensors MVIC force (Martin et al. 2004), with an α at 0.05 and a β at 0.2. Prior to their inclusion in this study, the participants were screened for the following exclusion criteria: smoking, current medication or drug consumption, and presence of apparent cardiovascular, metabolic, neurologic, or musculoskeletal disease. Furthermore, participants were excluded if they practiced regular physical activity involving substantial eccentric contractions within a six-month period prior to the start date of this study. Recreational running was tolerated except for prolonged (> 10 min) and repeated (more than 2 times per week) downhill running practice. Participants were also excluded if they were already familiarized with other types of eccentric exercises (e.g., strength training) prior to this study. They were also instructed to avoid (i) aspirin, ibuprofen, or other anti-inflammatory drugs; (ii) the consumption of caffeine and/or alcohol on the day of the experiment; and (iii) any strenuous exercise during the 48 h prior to testing. All procedures were approved by the local ethics committee (RCB number ID-RCB: 2019-A03012-55) and this study conformed to the standards set by the Declaration of Helsinki, except for registration in a database.
Experimental procedures
A pre-post cross-sectional design was used. Each participant attended the laboratory on two separate occasions. Room temperature and humidity were similar across sessions (21 °C, 45% RH). They first came for a familiarization session, and again one week later for the experimental session. During the familiarization session, the participants performed maximal, as well as impulsive, knee extensors voluntary isometric contractions and a single 5-min downhill run, replicating the experimental settings on the same treadmill used for the experimental session. During the experimental session, the participants performed a warm-up consisting of 10 min of light pedaling on a cycle ergometer (Monark 818E, Stockholm, Sweden) and 20 × 1-to-2-s knee extensors voluntary isometric contractions (separated by 10 s each) starting from a self-selected force and progressively increasing until the maximal volitional force was exerted (Coratella et al. 2020). After 5 min, MVIC and RFD were tested before performing a 30-min downhill run at a running speed of 10 km h−1 (2.78 m s−1) and a slope of − 20% (11.3°) on a motorized treadmill (Medic 2855, Genin Medical, La Roque-d'Anthéron, France). The protocol was chosen to induce significant neuromuscular fatigue, as shown by Martin et al. (2004) and (Khassetarash et al. 2021). Treadmill downhill running was selected since it was illustrated that overground graded running is well replicated on a treadmill (Firminger et al. 2018). After 90 s from the cessation of the downhill running bout, the participants performed the post-test evaluation of MVIC and RFD. This time window of 90 s was kept consistent for all participants, and it was chosen after pilot tests as it allowed enough time to move and install the participant on the isometric chair after the downhill running bout.
Neuromuscular evaluation
The testing protocol for MVIC and RFD was identical before (PRE) and after (POST). The downhill running session consisted in a 5-s MVIC followed by a series of eight impulsive contractions. Contractions were separated by 5 s of rest. At PRE, two MVICs were performed after having completed the warm-up protocol, and the peak value was retained for further analysis. If the difference between the two MVICs was > 5%, further trials were performed until the difference between two consecutive trials was < 5%.
The instruction given for all MVICs was to contract “as hard as possible,” while for impulsive contractions the participants were asked to extend the knee “as fast as possible” without any countermovement (Varesco et al. 2019). In case of a countermovement (determined by a force drop of 2 N below the baseline right before the impulsive contraction) or pre-tension (determined by a force level ≥ 2 N above the baseline right before the impulsive contraction) the contraction was repeated. The impulsive contractions were also repeated if the force level was < 70% of the MVIC that preceded the series of impulsive contractions (Varesco et al. 2019). All testing protocols were performed with real-time visual feedback.
Experimental setting
The participants were seated upright with hip and knee angles set at 90° of flexion on an isometric chair, equipped with a calibrated force transducer (Legcontrol, Mtraining, Ecole Valentin, France). The lever arm of the dynamometer was adjusted to firmly attach the leg 3 cm above the medial malleoli with two non-compliant belts. A belt strapped over the waist was used to minimize extraneous movements of the upper body. Passive resting force was subtracted from the signal so that the baseline was set at 0 N.
The EMG activity was recorded from the vastus lateralis muscle of the right leg as a surrogate for knee extensors activity (Place et al. 2010; Coratella et al. 2018). Self-adhesive Ag/AgCl surface electrodes (recording diameter = 10 mm, Kendall MediTrace foam electrode, MA) were placed in bipolar configuration over the muscle belly with an interelectrode (center to center) distance of 30 mm (Hermens et al. 2000). A reference (ground) electrode was placed on the right patella. Prior to electrode placement, the skin was shaved, abraded, and gently cleaned with an alcohol swab to lower the impedance. The EMG and force signals were collected simultaneously without analogical filters at a frequency of 2 kHz by PowerLab System (16/30-ML880/P, AD Instruments, Bellavista, Australia) and transferred on a personal computer using Labchart 8 software (ADInstruments) interface. The EMG signals were amplified (gain = 500) using an octal bio-amplifier (ML138, AD Instruments) with a bandwidth frequency range of 10–500 Hz obtained through a digital band-pass filter. EMG data were rectified and smoothed using a moving root mean square (RMS) with a time constant interval of 50 ms (Aagaard et al. 2002; Varesco et al. 2019).
Additional experiment
To confirm that changes in MVIC and RFD at POST results from the downhill running intervention, eight additional participants were tested as control group (age: 24 ± 5 yr, height: 176 ± 9 cm, body mass: 69 ± 8 kg). The setup and procedures were the same as the testing of the main study. However, these participants performed the neuromuscular evaluations PRE and POST with a resting period of equivalent length to the downhill running intervention.
Data analysis
For each maximal contraction, MVIC force and the associated RMS EMG signal (EMGmax) were calculated over the highest 1-s window. For each impulsive contraction, the peak RFD was determined as the maximum value on the force–time derivative curve over a 10-ms period to identify the best five impulsive contractions over all those performed (Varesco et al. 2019). The five contractions retained were averaged for further analysis, while the others were discarded. The onset of the force development was automatically defined as the point at which force exceeded the average resting baseline by ~ 2 N. The onset was also checked visually by an experienced investigator blinded to the condition (PRE or POST). Force was then measured at 25, 50, 100, and 200 ms after the onset (Fig. 1A). RFD was calculated at three time intervals of 0–50 ms (RFD0–50), 50–100 ms (RFD50–100), and 100–200 ms (RFD100–200) from the mean linear slope of the force–time curve (Aagaard et al. 2002; Coratella et al. 2020). The RFD (RFD0–200) was also calculated from contraction onset to 200 ms as an overall index of rapid force production ability at PRE and POST. The RER was calculated as the mean linear slope of the RMS EMG–time curve and analyzed at time intervals of 0–30 ms (RER0–30), 0–50 ms (RER0–50), and 0–75 ms (RER0–75) (Fig. 1B). The onset was visually identified on the RMS EMG–time curve as a rise of the RMS EMG signal over the baseline (Varesco et al. 2019). As the RMS EMG signal was nearly silent before the contraction, it was not necessary to set a threshold for noise (Fig. 1B). The long interval for the RER was measured on the first 75 ms of contraction as suggested by (Aagaard et al. 2002). This is also because, for our dataset, a decrease in RMS EMG signal amplitude occurred after ~ 80–100 ms from the onset (Fig. 1B). All RFD parameters are presented in both absolute and normalized to MVIC values. All RER intervals are presented in both absolute and normalized to EMGmax values. All data were analyzed offline using Labchart 8 Software (ADInstruments).
Statistical analysis
Data are reported as mean ± standard deviation. Assumption of normality (Shapiro–Wilk test) was tested for all variables. Generalized estimating equations were employed to account for the unbalanced participant numbers on data relative to baseline (n = 13 for the downhill running intervention and n = 8 for the control group) (Liang and Zeger 1986). If significant time (PRE vs. POST) × group (downhill running vs. resting) interactions were observed on MVIC and RFD0–200, the etiology was investigated using paired samples t-tests for normally distributed variables (MVIC, EMGmax, absolute and normalized RFD windows, normalized RER variables) to compare differences between PRE and POST. For non-normally distributed variables (absolute RER variables), Wilcoxon Signed-Rank Tests were used to compare differences between PRE and POST. As a measure of effect size for parametric tests, Cohen’s d was calculated with 95% confidence intervals (CI). Values of ± 0.20, ± 0.50, and over ± 0.80 were considered as small, medium, and large effect, respectively (Cohen 1988). As a measure of the effect size for non-parametric tests, eta squared (η2) was calculated with 95% CI. Values of ± 0.01, ± 0.06, and over ± 0.14 were considered as small, medium, and large effect, respectively (Cohen 1988). The criterion α-level was set to 0.05. The statistical analysis was conducted using IBM™ SPSS™ Statistics (Version 28.0, IBM Corp., Somers, New York, NY).
Results
Significant time × group interactions were observed for MVIC (χ2 (1) = 42.817, P < 0.001), normalized peak RFD (χ2 (1) = 24.909, P < 0.001), RFD0–200 (χ2 (1) = 5.847, P = 0.016), and normalized RFD0–200 (χ2 (1) = 5.700, P = 0.017). For all these parameters, the two groups were similar at PRE (all P ≥ 0.375). At POST, the two groups differed with the group of participants who performed the downhill running intervention that presented lower values than the group of participants who rested for equivalent length to the downhill running intervention (all P ≤ 0.018) (Table 1). Therefore, only changes in MVIC and RFD at POST for the downhill running group are presented.
For the group of participants who performed the downhill running intervention, at POST, MVIC and RFD0–200 decreased by ~ 25% [d = − 1.09 (95% CI − 1.88/− 0.24); P < 0.001; Fig. 2A] and ~ 15% [d = − 0.50 (− 1.26/0.30); P = 0.001; Fig. 2B], respectively. When normalized to the corresponding MVIC, RFD0–200 was ~ 11% greater at POST compared to PRE [d = 0.71 (− 0.10/1.48); P = 0.024; Fig. 2C]. Finally, no significant differences were found from PRE to POST for peak RFD (PRE = 7866 ± 3283 N s−1; POST = 7284 ± 2655 N s−1; d = 0.43 (− 0.15/0.10); P = 0.072) and for time to peak RFD (PRE = 0.054 ± 0.010 s; POST = 0.060 ± 0.015 s; d = − 0.39 (− 0.95/0.18); P = 0.965). For the group of participants who rested for equivalent length to the downhill running intervention, no changes were observed for MVIC, RFD0–200, and RFD0–200 normalized to MVIC (all P > 0.05; Fig. 2).
For the group of participants who performed the downhill running intervention, EMGmax showed a significant decrease of ~ 22% at POST [PRE: 0.28 ± 0.18 mV; POST: 0.22 ± 0.14 mV; d = − 0.37 (− 1.13/0.42); P = 0.038]. Absolute RER of the vastus lateralis muscle did not change at any time interval: 0–30 ms [PRE: 4.93 ± 2.84 mV; POST: 5.35 ± 2.64 mV; η2 = 0.086 (− 0.69/0.85); P = 0.311], 0–50 ms [PRE: 5.47 ± 3.64 mV; POST: 5.34 ± 2.84 mV; η2 = 0.017 (− 0.75/0.79); P = 0.650], or 0–75 ms [PRE: 4.03 ± 3.12 mV; POST: 3.89 ± 1.90 mV; η2 = 0.023 (− 0.55/0.99); P = 0.600]. Similar results were observed when RER was expressed as %EMGmax: RER0–30 [PRE: 2203 ± 1138%; POST: 3311 ± 2618%; d = 0.55 (− 0.25/1.31); P = 0.090], RER0–50 [PRE: 2278 ± 1147%; POST: 3052 ± 2037%; d = 0.47 (− 0.32/1.23); P = 0.153], and RER0–75 [PRE: 1504 ± 791%; POST: 2213 ± 1344%; d = 0.64 (− 0.17/1.41); P = 0.079].
Figure 3A and B shows the force values at different time points from the onset of the contraction, respectively. Only for the group of participants who performed the downhill running intervention, force at 200 ms decreased by ~ 16% at POST [d = − 0.54 (− 1.30/0.26); P = 0.008], while no differences were observed at 25 ms [d = − 0.42 (− 1.18/0.37); P = 0.227], 50 ms [d = − 0.38 (− 1.14/0.41); P = 0.203] or 100 ms [d = − 0.21 (− 0.97/0.57); P = 0.254]. When normalized to MVIC (Fig. 3C and D), the force at 100 ms increased by ~ 20% at POST [d = 0.83 (0.005/1.60), P = 0.016]. No differences were observed for normalized data at 25 ms [d = 0.37 (− 0.42/1.13); P = 0.267], 50 ms [d = 0.49 (− 0.30/1.25); P = 0.117] or 200 ms [d = 0.61 (− 0.19/1.38); P = 0.088]. No changes were observed for the group of participants who rested for equivalent length to the downhill running intervention (all P > 0.05).
Figure 4A shows the RFD at different time windows from the onset of the contraction. No changes were observed for the group of participants who rested for equivalent length to the downhill running intervention (all P > 0.05). For the group of participants who performed the downhill running intervention, RFD100–200 was ~ 25% lower at POST compared to PRE [d = − 0.70 (− 1.47/0.11); P < 0.001]. No differences were observed for RFD0–50 [d = − 0.40 (− 1.16/0.39); P = 0.125] or RFD50–100 [d = 0.09 (− 0.68/0.86), P = 0.599]. When normalized to MVIC (Fig. 4B), no differences were observed for the group of participants who rested for equivalent length to the downhill running intervention (all P > 0.05). Conversely, for the group of participants who performed the downhill running intervention, only RFD50–100 increased by ~ 26% from PRE to POST [d = 0.61 (− 0.19/1.38); P = 0.030]. No differences were observed for RFD0-50 [d = 0.35 (− 1.12/0.43); P = 0.179] or RFD100–200 [d = 0.07 (− 0.84/0.70); P = 0.763].
Discussion
The current study examined for the first time the capacity of the knee extensors to rapidly produce force after a 30-min eccentric-biased downhill run. The downhill running bout induced a drop in MVIC and impaired the overall ability to rapidly develop force, as indicated by the reduction in RFD0–200. When analyzing the different RFD time windows, we observed a reduction in the 100–200 ms phase, while the early phase of the RFD remains unchanged. This latter was accompanied by a lack of change in the absolute RER. Together, these data suggest that a fatiguing downhill run affects the force-generating capacity and the RFD during isometric impulsive contractions through structural rather than neural impairments.
Maximal isometric force
After the 30-min downhill run, MVIC decreased by ~ 25%, which was consistent with the values reported in the literature (Martin et al. 2004; Giandolini et al. 2016a; Ehrström et al. 2018; Lemire et al. 2020; Khassetarash et al. 2021). The force loss after downhill running has been related to impairments in both neural and contractile components (Giandolini et al. 2016a; Ehrström et al. 2018; Khassetarash et al. 2021). Specifically, the former seems attributable to a deficit in muscle voluntary activation caused by several factors, such as neurobiological alterations in the brain and/or changes in the intrinsic properties of the motoneuron pool (Gandevia 2001). In contrast, the latter appear primarily due to impairments in one or more steps involved in the excitation–contraction coupling, such as decoupling at the T-tubule–sarcoplasmic reticulum interface and/or a decreased Ca2+ release from the sarcoplasmic reticulum (Giandolini et al. 2016b; Bontemps et al. 2020). Furthermore, EMGmax was lower at POST compared to PRE, and this has been associated with a reduced firing rate of high threshold motor units (Balshaw et al. 2017).
Maximal rate of force development
RFD0–200 decreased by ~ 15% after the 30-min downhill run, which is less than the ~ 25% drop observed in MVIC. Consequently, normalized RFD (%MVIC) was greater, indicating that downhill running-induced fatigue affected the maximal force capacity to a greater extent compared to RFD. Previously, Maeo et al. (2017) evaluated the knee extensors MVIC and RFD responses following 45 min of downhill running at ~ 10 km h−1 and − 15% slope. Twenty-four hours after the exercise, the authors observed a decrease of ~ 16% and ~ 20% in MVIC and RFD0–200, respectively. However, neither the different phases of the force–time curve of RFD nor the RER were investigated. Furthermore, RFD was calculated on the rising in force during MVIC and attempting to measure RFD during an MVIC may have resulted in an underestimation of RFD (Maffiuletti et al. 2016). These limitations impede direct comparisons with the current study. Recently, using a similar protocol, Khassetarash et al. (2021) found a 16% decrease in MVIC, while we found much larger impairments (− 25%). Those authors did not evaluated RFD; however, they studied the possible etiology for the MVIC force loss using the interpolated twitch technique with electrical stimulations. They found impairments in both voluntary activation (~ 8.6%) and artificially evoked force responses (~ 4.9%). Therefore, it is possible that, in our study, the downhill running bout induced impairments in both maximal voluntary activation and contractile function. During a 5-s MVIC, the time is sufficient to recruit most of the motor unit pool (De Luca 1985; Maffiuletti et al. 2016). It is then possible that impairment in maximal voluntary activation or contractile function would manifest to a larger extent during MVIC than during RFD. While during RFD contraction it would be possible to detect the rapid activation of the muscle rather than the maximal one (Maffiuletti et al. 2016; D’Emanuele et al. 2021).
Early and late phases of the rate of force development
No changes were observed for both the early phase of the RFD and the absolute RER, suggesting no detectable alterations in the generation and/or transmission of firing bursts along the neuromuscular pathway (Klass et al. 2008; Del Vecchio et al. 2019).
In the present study, the downhill running-induced fatigue altered the late phase of the RFD, as indicated by a decrease in RFD100–200. Several possible mechanisms may have contributed to the decrease in the late phase of the RFD. First, the responsiveness of a muscle to high firing rates from the motoneuron pool depends on both the Ca2+ sensitivity and the number of actin-myosin bridges available for contractions (Allen et al. 2008). The repetitive eccentric muscle actions that occur during the braking phase of each landing step of a downhill run (Vernillo et al. 2017) were shown to cause reductions in the Ca2+ sensitivity, actin–myosin coupling (Giandolini et al. 2016b; Bontemps et al. 2020), and contiguity of the Z-lines of the sarcomere (Féasson et al. 2002). These structural alterations within the muscle likely affected the interface sarcoplasmic reticulum/T-tubule (Giandolini et al. 2016b), negatively influencing the late phase of the RFD after downhill running. Second, the decrease in the late phase of the RFD could have also been modulated by an alteration in tendon stiffness. Indeed, repetitive eccentric muscle actions were shown to induce considerable stress on the passive structures of the muscle (i.e., muscle connective tissue and tendons) (Reeves and Narici 2003). During a downhill run this could have led to repetitive elongations of the knee extensors muscle–tendon unit, reducing stiffness and consequently increasing compliance (Guilhem et al. 2011). The reduced stiffness of the muscle–tendon unit could then explain the observed decreased late-phase RFD after the downhill running with an altered force transmission (Hannah and Folland 2015). Third, downhill running may have induced changes in the muscle architecture. Fatiguing tasks may result in a decreased fascicle length (Csapo et al. 2011), and considering that fascicle length is a possible predictor of the RFD late phase (Coratella et al. 2020), it would then be possible that downhill running may have affected the fascicle length. However, it is acknowledged that this remains to be proven in future investigations. Our results suggest that the ability to rapidly activate the muscle after downhill running was preserved, as shown by the absence of significant change in RER and RFD0–50 from PRE to POST. Although speculative, it is possible that downhill running would only cause structural alterations (or a reduction in the total number of muscle fibers activated during the rapid contractions), thereby impairing RFD ability.
Limitations
Some possible limitations could have influenced the outcomes assessed in the present study. First, our POST measurements were collected after 90 s of recovery. As recovery kinetics at the level of the central nervous system is rapid (Kennedy et al. 2016; Mira et al. 2017; Vernillo et al. 2018; Krüger et al. 2019; Koral et al. 2020; Vernillo et al. 2020b; Ducrocq et al. 2021), we cannot rule out that the capacity to deliver maximal discharge frequency and incidence of doublet discharge during fast contractions (i.e., the neural factors influencing RFD) has been underestimated in the present study. Second, RFD contractions also depend on the ability of the individual to perform impulsive contractions. We limited its possible effect by performing a familiarization session where participants were instructed by an expert investigator to perform impulsive contractions and strongly encouraged during the experiments. Finally, trains of high frequency evoked contractions would have provided additional information independent of any volitional input (to better discern voluntary vs. involuntary muscle rapid activation capacity). However, these procedures are painful and could increase the risk of injury (Millet et al. 2012).
Conclusions
The present findings showed that RFD is altered after a single bout of downhill running. Specifically, we observed that only the late phase of the force–time curve during impulsive contractions was altered. It seems that a bout of downhill running influenced more the structural than the neuromuscular factors, as indicated by the drop in the late but not early phase of RFD. Since downhill running is an essential component of trail running performance, these findings may help explain evidence of neuromuscular alterations in trail runners and following prolonged duration races wherein cumulative eccentric loading is high.
Availability of data and materials
Derived data supporting the findings of this study are available on request from the corresponding author.
Code availability
Not applicable.
Abbreviations
- CI:
-
95% Confidence intervals
- EMG:
-
Electromyography
- MVIC:
-
Maximum voluntary isometric contraction
- POST:
-
Neuromuscular evaluation performed after the fatiguing exercise
- PRE:
-
Neuromuscular evaluation performed before the fatiguing exercise
- RER:
-
Rate of electromyography signal rise
- RFD:
-
Rate of force development
- RFD0–50 :
-
Rate of force development calculated at the time interval of 0–50 ms
- RFD50–100 :
-
Rate of force development calculated at the time interval of 50–100 ms
- RFD100–200 :
-
Rate of force development calculated at the time interval of 100–200 ms
- RMS:
-
Root mean square
References
Aagaard P, Simonsen EB, Andersen JL, Magnusson P, Dyhre-Poulsen P (2002) Increased rate of force development and neural drive of human skeletal muscle following resistance training. J Appl Physiol 93(4):1318–1326. https://doi.org/10.1152/japplphysiol.00283.2002
Allen DG, Lamb GD, Westerblad H (2008) Skeletal muscle fatigue: cellular mechanisms. Physiol Rev 88(1):287–332. https://doi.org/10.1152/physrev.00015.2007
Balshaw TG, Pahar M, Chesham R, Macgregor LJ, Hunter AM (2017) Reduced firing rates of high threshold motor units in response to eccentric overload. Physiol Rep. https://doi.org/10.14814/phy2.13111
Boccia G, Dardanello D, Zoppirolli C, Bortolan L, Cescon C, Schneebeli A, Vernillo G, Schena F, Rainoldi A, Pellegrini B (2017) Central and peripheral fatigue in knee and elbow extensor muscles after a long-distance cross-country ski race. Scand J Med Sci Sports 27(9):945–955. https://doi.org/10.1111/sms.12718
Boccia G, Dardanello D, Brustio PR, Tarperi C, Festa L, Zoppirolli C, Pellegrini B, Schena F, Rainoldi A (2018) Neuromuscular fatigue does not impair the rate of force development in ballistic contractions of submaximal amplitudes. Front Physiol 9:1503. https://doi.org/10.3389/fphys.2018.01503
Bontemps B, Vercruyssen F, Gruet M, Louis J (2020) Downhill running: what are the effects and how can we adapt? A narrative review. Sports Med 50(12):2083–2110. https://doi.org/10.1007/s40279-020-01355-z
Cohen J (1988) Statistical power analysis for the behavioral sciences. Lawrence Erlbaum, Mahwah
Coratella G, Grospretre S, Gimenez P, Mourot L (2018) Greater fatigability in knee-flexors vs. knee-extensors after a standardized fatiguing protocol. Eur J Sport Sci 18(8):1110–1118. https://doi.org/10.1080/17461391.2018.1469674
Coratella G, Longo S, Borrelli M, Doria C, Ce E, Esposito F (2020) Vastus intermedius muscle architecture predicts the late phase of the knee extension rate of force development in recreationally resistance-trained men. J Sci Med Sport 23(11):1100–1104. https://doi.org/10.1016/j.jsams.2020.04.006
Csapo R, Alegre LM, Baron R (2011) Time kinetics of acute changes in muscle architecture in response to resistance exercise. J Sci Med Sport 14(3):270–274. https://doi.org/10.1016/j.jsams.2011.02.003
D’Emanuele S, Maffiuletti NA, Tarperi C, Rainoldi A, Schena F, Boccia G (2021) Rate of force development as an indicator of neuromuscular fatigue: a scoping review. Front Hum Neurosci 15:1
De Luca CJ (1985) Control properties of motor units. J Exp Biol 115:125–136
Del Vecchio A, Negro F, Holobar A, Casolo A, Folland JP, Felici F, Farina D (2019) You are as fast as your motor neurons: speed of recruitment and maximal discharge of motor neurons determine the maximal rate of force development in humans. J Physiol 597(9):2445–2456. https://doi.org/10.1113/JP277396
Ducrocq GP, Hureau TJ, Bogseth T, Meste O, Blain GM (2021) Recovery from fatigue after cycling time trials in elite endurance athletes. Med Sci Sports Exerc 53(5):904–917. https://doi.org/10.1249/MSS.0000000000002557
Ehrström S, Gruet M, Giandolini M, Chapuis S, Morin JB, Vercruyssen F (2018) Acute and delayed neuromuscular alterations induced by downhill running in trained trail runners: beneficial effects of high-pressure compression garments. Front Physiol 9:1627. https://doi.org/10.3389/fphys.2018.01627
Féasson L, Stockholm D, Freyssenet D, Richard I, Duguez S, Beckmann JS, Denis C (2002) Molecular adaptations ofneuromuscular disease-associated proteins in response to eccentric exercise in human skeletal muscle. J Physiol 543(Pt1):297–306. https://doi.org/10.1113/jphysiol.2002.018689
Firminger CR, Vernillo G, Savoldelli A, Stefanyshyn DJ, Millet GY, Edwards WB (2018) Joint kinematics and ground reaction forces in overground versus treadmill graded running. Gait Posture 63:109–113. https://doi.org/10.1016/j.gaitpost.2018.04.042
Gandevia SC (2001) Spinal and supraspinal factors in human muscle fatigue. Physiol Rev 81(4):1725–1789. https://doi.org/10.1152/physrev.2001.81.4.1725
Giandolini M, Horvais N, Rossi J, Millet GY, Morin JB, Samozino P (2016a) Acute and delayed peripheral and central neuromuscular alterations induced by a short and intense downhill trail run. Scand J Med Sci Sports 26(11):1321–1333. https://doi.org/10.1111/sms.12583
Giandolini M, Vernillo G, Samozino P, Horvais N, Edwards WB, Morin JB, Millet GY (2016b) Fatigue associated with prolonged graded running. Eur J Appl Physiol 116(10):1859–1873. https://doi.org/10.1007/s00421-016-3437-4
Guilhem G, Cornu C, Guevel A (2011) Muscle architecture and EMG activity changes during isotonic and isokinetic eccentric exercises. Eur J Appl Physiol 111(11):2723–2733. https://doi.org/10.1007/s00421-011-1894-3
Hannah R, Folland JP (2015) Muscle-tendon unit stiffness does not independently affect voluntary explosive force production or muscle intrinsic contractile properties. Appl Physiol Nutr Metab 40(1):87–95. https://doi.org/10.1139/apnm-2014-0343
Hermens HJ, Freriks B, Disselhorst-Klug C, Rau G (2000) Development of recommendations for SEMG sensors and sensor placement procedures. J Electromyogr Kinesiol 10(5):361–374
Hoffman MD, Ong JC, Wang G (2010) Historical analysis of participation in 161 km ultramarathons in North America. Int J Hist Sport 27(11):1877–1891. https://doi.org/10.1080/09523367.2010.494385
Kennedy DS, McNeil CJ, Gandevia SC, Taylor JL (2016) Effects of fatigue on corticospinal excitability of the human knee extensors. Exp Physiol 101(12):1552–1564. https://doi.org/10.1113/EP085753
Khassetarash A, Vernillo G, Martinez A, Baggaley M, Giandolini M, Horvais N, Millet GY, Edwards WB (2020) Biomechanics of graded running: Part II-Joint kinematics and kinetics. Scand J Med Sci Sports 30(9):1642–1654. https://doi.org/10.1111/sms.13735
Khassetarash A, Vernillo G, Kruger RL, Edwards WB, Millet GY (2021) Neuromuscular, biomechanical, and energetic adjustments following repeated bouts of downhill running. J Sport Health Sci. https://doi.org/10.1016/j.jshs.2021.06.001
Klass M, Baudry S, Duchateau J (2008) Age-related decline in rate of torque development is accompanied by lower maximal motor unit discharge frequency during fast contractions. J Appl Physiol 104(3):739–746. https://doi.org/10.1152/japplphysiol.00550.2007
Koral J, Oranchuk DJ, Wrightson JG, Twomey R, Millet GY (2020) Mechanisms of neuromuscular fatigue and recovery in unilateral versus bilateral maximal voluntary contractions. J Appl Physiol (1985) 128(4):785–794. https://doi.org/10.1152/japplphysiol.00651.2019
Krüger RL, Aboodarda SJ, Jaimes LM, Samozino P, Millet GY (2019) Cycling performed on an innovative ergometer at different intensities-durations in men: neuromuscular fatigue and recovery kinetics. Appl Physiol Nutr Metab 44(12):1320–1328. https://doi.org/10.1139/apnm-2018-0858
Lemire M, Remetter R, Hureau TJ, Kouassi BYL, Lonsdorfer E, Geny B, Isner-Horobeti ME, Favret F, Dufour SP (2020) High-intensity downhill running exacerbates heart rate and muscular fatigue in trail runners. J Sports Sci. https://doi.org/10.1080/02640414.2020.1847502
Liang KY, Zeger SL (1986) Longitudinal data analysis using generalized linear models. Biometrika 73(1):13–22. https://doi.org/10.1093/biomet/73.1.13
Lima LCR, Nosaka K, Chen TC, Pinto RS, Greco CC, Denadai BS (2020) Decreased running economy is not associated with decreased force production capacity following downhill running in untrained, young men. Eur J Sport Sci. https://doi.org/10.1080/17461391.2020.1727570
Maeo S, Ando Y, Kanehisa H, Kawakami Y (2017) Localization of damage in the human leg muscles induced by downhill running. Sci Rep 7(1):5769. https://doi.org/10.1038/s41598-017-06129-8
Maffiuletti NA, Aagaard P, Blazevich AJ, Folland J, Tillin N, Duchateau J (2016) Rate of force development: physiological and methodological considerations. Eur J Appl Physiol 116(6):1091–1116. https://doi.org/10.1007/s00421-016-3346-6
Martin V, Millet GY, Martin A, Deley G, Lattier G (2004) Assessment of low-frequency fatigue with two methods of electrical stimulation. J Appl Physiol 97(5):1923–1929. https://doi.org/10.1152/japplphysiol.00376.2004
Millet GY, Tomazin K, Verges S, Vincent C, Bonnefoy R, Boisson RC, Gergele L, Feasson L, Martin V (2011) Neuromuscular consequences of an extreme mountain ultra-marathon. PLoS ONE 6(2):e17059. https://doi.org/10.1371/journal.pone.0017059
Millet GY, Bachasson D, Temesi J, Wuyam B, Feasson L, Verges S, Levy P (2012) Potential interests and limits of magnetic and electrical stimulation techniques to assess neuromuscular fatigue. Neuromuscul Disord 22(Suppl 3):S181-186. https://doi.org/10.1016/j.nmd.2012.10.007
Mira J, Lapole T, Souron R, Messonnier L, Millet GY, Rupp T (2017) Cortical voluntary activation testing methodology impacts central fatigue. Eur J Appl Physiol 117(9):1845–1857. https://doi.org/10.1007/s00421-017-3678-x
Nicol C, Avela J, Komi PV (2006) The stretch-shortening cycle: a model to study naturally occurring neuromuscular fatigue. Sports Med 36(11):977–999. https://doi.org/10.2165/00007256-200636110-00004
Novacheck TF (1998) The biomechanics of running. Gait Posture 7(1):77–95. https://doi.org/10.1016/s0966-6362(97)00038-6
Place N, Casartelli N, Glatthorn JF, Maffiuletti NA (2010) Comparison of quadriceps inactivation between nerve and muscle stimulation. Muscle Nerve 42(6):894–900. https://doi.org/10.1002/mus.21776
Reeves ND, Narici MV (2003) Behavior of human muscle fascicles during shortening and lengthening contractions in vivo. J Appl Physiol 95(3):1090–1096. https://doi.org/10.1152/japplphysiol.01046.2002
Saugy J, Place N, Millet GY, Degache F, Schena F, Millet GP (2013) Alterations of neuromuscular function after the World’s most challenging mountain ultra-marathon. PLoS ONE 8(6):e65596. https://doi.org/10.1371/journal.pone.0065596
Scheer V, Basset P, Giovanelli N, Vernillo G, Millet GP, Costa RJS (2020) Defining off-road running: a position statement from the Ultra Sports Science Foundation. Int J Sports Med 41(5):275–284. https://doi.org/10.1055/a-1096-0980
Temesi J, Rupp T, Martin V, Arnal PJ, Feasson L, Verges S, Millet GY (2014) Central fatigue assessed by transcranial magnetic stimulation in ultratrail running. Med Sci Sports Exerc 46(6):1166–1175. https://doi.org/10.1249/MSS.0000000000000207
Varesco G, Espeit L, Feasson L, Lapole T, Rozand V (2019) Rate of force development and rapid muscle activation characteristics of knee extensors in very old men. Exp Gerontol 124:110640. https://doi.org/10.1016/j.exger.2019.110640
Vernillo G, Giandolini M, Edwards WB, Morin JB, Samozino P, Horvais N, Millet GY (2017) Biomechanics and physiology of uphill and downhill running. Sports Med 47(4):615–629. https://doi.org/10.1007/s40279-016-0605-y
Vernillo G, Temesi J, Martin M, Millet GY (2018) Mechanisms of fatigue and recovery in upper versus lower limbs in men. Med Sci Sports Exerc 50(2):334–343. https://doi.org/10.1249/MSS.0000000000001445
Vernillo G, Aguiar M, Savoldelli A, Martinez A, Giandolini M, Horvais N, Edwards WB, Millet GY (2019) Regular changes in foot strike pattern during prolonged downhill running do not influence neuromuscular, energetics, or biomechanical parameters. Eur J Sport Sci. https://doi.org/10.1080/17461391.2019.1645212
Vernillo G, Martinez A, Baggaley M, Khassetarash A, Giandolini M, Horvais N, Edwards WB, Millet GY (2020a) Biomechanics of graded running: Part I—Stride parameters, external forces, muscle activations. Scand J Med Sci Sports 30(9):1632–1641. https://doi.org/10.1111/sms.13708
Vernillo G, Temesi J, Martin M, Kruger RL, Millet GY (2020b) Spinal contribution to neuromuscular recovery differs between elbow-flexor and knee-extensor muscles after a maximal sustained fatiguing task. J Neurophysiol 124(3):763–773. https://doi.org/10.1152/jn.00273.2020
Acknowledgements
The authors thank our participants for their time and effort. They also thank Ingrid Giuliana Ulloa Gutiérrez and Léo Vanini for their valuable technical assistance during data acquisition; Dr. Guillaume Y. Millet and Dr. Vincent Martin for sharing the raw data of their work (Martin et al. 2004) that was used for sample size calculations; and Hayley L. Patmore for her assistance in English proof reading. This research was supported by grants from the Région Bourgogne Franche Comté (2019-Y-09131 to LM), the French Ministry of National Education, Research and Technology (#EA3920 to LM), the Tomsk Polytechnic University Competitiveness Enhancement Program (#BBE-BCUT-108/2017-TPU CEP-HSTI-108/2017 to LM), and the Università degli Studi di Milano (#PSR2019_VERNILLO to GV).
Funding
Open access funding provided by Università degli Studi di Milano within the CRUI-CARE Agreement. This research was supported by grants from: the Région Bourgogne Franche Comté (2019-Y-09131 to LM), the French Ministry of National Education, Research and Technology (#EA3920 to LM), the Tomsk Polytechnic University Development Program (#TPU to LM), and the Università degli Studi di Milano (#PSR2019_VERNILLO to GV).
Author information
Authors and Affiliations
Contributions
GVa, GC, GL, LM, and GVe conceived and designed research; GVa and BC performed experiments; GVa, GC, and GVe interpreted results of experiments and drafted manuscript; GVa analyzed data and prepared figures; GVa, GC, VR, BC, GL, LM, and GVe edited and revised manuscript; GVa, GC, VR, BC, GL, LM, and GVe approved final version of manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors have no conflicts of interest to declare.
Ethics approval
All procedures were approved by the local ethics committee (RCB number ID-RCB: 2019-A03012-55).
Consent to participate
Informed consent was obtained from all participants in this study.
Consent for publication
Not applicable.
Additional information
Communicated by Nicolas Place.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Varesco, G., Coratella, G., Rozand, V. et al. Downhill running affects the late but not the early phase of the rate of force development. Eur J Appl Physiol 122, 2049–2059 (2022). https://doi.org/10.1007/s00421-022-04990-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00421-022-04990-8