# Downhill running increases markers of muscle damage and impairs the maximal voluntary force production as well as the late phase of the rate of voluntary force development

Giuseppe Coratella, Giorgio Varesco, Vianney Rozand, Benjamin Cuinet, Veronica Sansoni, Giovanni Lombardi, Gianluca Vernillo, Laurent Mourot

► To cite this version:

Giuseppe Coratella, Giorgio Varesco, Vianney Rozand, Benjamin Cuinet, Veronica Sansoni, et al. Downhill running increases markers of muscle damage and impairs the maximal voluntary force production as well as the late phase of the rate of voluntary force development. European Journal of Applied Physiology, 2024, 124, pp.1875-1883. 10.1007/s00421-023-05412-z. hal-04596352

HAL Id: hal-04596352

https://hal.science/hal-04596352

Submitted on 31 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Downhill running increases markers of muscle damage and impairs the maximal voluntary force production as well as the late phase of the rate of voluntary force development

Giuseppe Coratella<fnref n="1" /> ID · Giorgio Varesco<fnref n="2" />,<fnref n="3" /> ID · Vianney Rozand<fnref n="2" /> ID · Benjamin Cuinet<fnref n="4" /> · Veronica Sansoni<fnref n="5" /> ID ·

Giovanni Lombardi<fnref n="5" />,<fnref n="6" /> ID · Gianluca Vernillo<fnref n="1" /> ID · Laurent Mourot<fnref n="4" /> ID

Received: 27 July 2023 / Accepted: 31 December 2023 / Published online: 10 January 2024

© The Author(s) 2024

## Abstract

**Purpose** To examined the time-course of the early and late phase of the rate of voluntary force development (RVFD) and muscle damage markers after downhill running.

**Methods** Ten recreational runners performed a 30-min downhill run at 10 km h⁻<fnref n="1" /> and -20% (-11.3°) on a motorized treadmill. At baseline and each day up to 4 days RVFD, knee extensors maximum voluntary isometric force (MVIC), serum creatine kinase (CK) concentration, quadriceps swelling, and soreness were assessed. The early (0–50 ms) and late (100–200 ms) phase of the RVFD, as well as the force developed at 50 and 200 ms, were also determined.

**Results** MVIC showed moderate decrements (p < 0.05) and recovered after 4 days (p > 0.05). Force at 50 ms and the early phase were not impaired (p > 0.05). Conversely, force at 200 ms and the late phase showed moderate decrements (p < 0.05) and recovered after 3 and 4 days, respectively (p > 0.05). CK concentration, quadriceps swelling, and soreness increased (p < 0.05) were overall fully resolved after 4 days (p > 0.05).

**Conclusion** Downhill running affected the knee extensors RVFD late but not early phase. The RVFD late phase may be used as an additional marker of muscle damage in trail running.

**Keywords** Creatine kinase · Delayed onset muscle soreness · Eccentric · Maximum voluntary contraction · Muscle

swelling · Musculoskeletal · Quadriceps · Trail running

## Abbreviations

CK Creatine kinase

MVIC Maximum voluntary isometric contraction

POST24 Neuromuscular evaluation performed 24 h after the fatiguing exercise

POST48 Neuromuscular evaluation performed 48 h after the fatiguing exercise

POST72 Neuromuscular evaluation performed 72 h after the fatiguing exercise

POST96 Neuromuscular evaluation performed 96 h after the fatiguing exercise

PRE Neuromuscular evaluation performed before the fatiguing exercise

RVFD Rate of voluntary force development

Communicated by Toshio Moritani.

Giuseppe Coratella and Giorgio Varesco have contributed equally to this work.

<figure><img src="image_1.png" /></figure>

Gianluca Vernillo
gianluca.vernillo@gmail.com

<footnote n="1">¹ Department of Biomedical Sciences for Health, Università Degli Studi Di Milano, Milan, Italy</footnote>

<footnote n="2">² Université Jean Monnet Saint-Etienne, Inter-University Laboratory of Human Movement Biology, 42023 Saint-Etienne, France</footnote>

<footnote n="3">³ Laboratory Movement—Interactions—Performance, MIP Lab, UR 4334, Nantes Université, F-44000 Nantes, France</footnote>

<footnote n="4">⁴ Prognostic Factors and Regulatory Factors of Cardiac and Vascular Pathologies (EA3920), Exercise Performance Health Innovation (EPHI) Platform, University of Franche-Comté, Besançon, France</footnote>

<footnote n="5">⁵ Laboratory of Experimental Biochemistry and Molecular Biology, IRCCS Instituto Ortopedico Galeazzi, 20157 Milan, Italy</footnote>

<footnote n="6">⁶ Department of Athletics, Strength and Conditioning, Poznań University of Physical Education, 61-871 Poznań, Poland</footnote>
# Downhill running increases markers of muscle damage and impairs the maximal voluntary force production as well as the late phase of the rate of voluntary force development

Giuseppe Coratella¹ ID · Giorgio Varesco²,³ ID · Vianney Rozand² ID · Benjamin Cuinet⁴ · Veronica Sansoni⁵ ID ·

Giovanni Lombardi⁵,⁶ ID · Gianluca Vernillo¹ ID · Laurent Mourot⁴ ID

Received: 27 July 2023 / Accepted: 31 December 2023 / Published online: 10 January 2024

© The Author(s) 2024

## Abstract

**Purpose** To examined the time-course of the early and late phase of the rate of voluntary force development (RVFD) and muscle damage markers after downhill running.

**Methods** Ten recreational runners performed a 30-min downhill run at 10 km h⁻¹ and -20% (-11.3°) on a motorized treadmill. At baseline and each day up to 4 days RVFD, knee extensors maximum voluntary isometric force (MVIC), serum creatine kinase (CK) concentration, quadriceps swelling, and soreness were assessed. The early (0–50 ms) and late (100–200 ms) phase of the RVFD, as well as the force developed at 50 and 200 ms, were also determined.

**Results** MVIC showed moderate decrements (p < 0.05) and recovered after 4 days (p > 0.05). Force at 50 ms and the early phase were not impaired (p > 0.05). Conversely, force at 200 ms and the late phase showed moderate decrements (p < 0.05) and recovered after 3 and 4 days, respectively (p > 0.05). CK concentration, quadriceps swelling, and soreness increased (p < 0.05) were overall fully resolved after 4 days (p > 0.05).

**Conclusion** Downhill running affected the knee extensors RVFD late but not early phase. The RVFD late phase may be used as an additional marker of muscle damage in trail running.

**Keywords** Creatine kinase · Delayed onset muscle soreness · Eccentric · Maximum voluntary contraction · Muscle

swelling · Musculoskeletal · Quadriceps · Trail running

## Abbreviations

CK Creatine kinase

MVIC Maximum voluntary isometric contraction

POST24 Neuromuscular evaluation performed 24 h after the fatiguing exercise

POST48 Neuromuscular evaluation performed 48 h after the fatiguing exercise

POST72 Neuromuscular evaluation performed 72 h after the fatiguing exercise

POST96 Neuromuscular evaluation performed 96 h after the fatiguing exercise

PRE Neuromuscular evaluation performed before the fatiguing exercise

RVFD Rate of voluntary force development

Communicated by Toshio Moritani.

Giuseppe Coratella and Giorgio Varesco have contributed equally to this work.

<figure><img src="image_2.png" /></figure>

Gianluca Vernillo
gianluca.vernillo@gmail.com

¹ Department of Biomedical Sciences for Health, Università Degli Studi Di Milano, Milan, Italy

² Université Jean Monnet Saint-Etienne, Inter-University Laboratory of Human Movement Biology, 42023 Saint-Etienne, France

³ Laboratory Movement—Interactions—Performance, MIP Lab, UR 4334, Nantes Université, F-44000 Nantes, France

⁴ Prognostic Factors and Regulatory Factors of Cardiac and Vascular Pathologies (EA3920), Exercise Performance Health Innovation (EPHI) Platform, University of Franche-Comté, Besançon, France

⁵ Laboratory of Experimental Biochemistry and Molecular Biology, IRCCS Instituto Ortopedico Galeazzi, 20157 Milan, Italy

⁶ Department of Athletics, Strength and Conditioning, Poznań University of Physical Education, 61-871 Poznań, Poland
RVFD<sub>0–50</sub> Rate of voluntary force development calculated at the 0–50 ms time interval

RVFD<sub>100–200</sub> Rate of voluntary force development calculated at the 100–200 ms time interval

## Introduction

Performing eccentric exercise is known to induce symptoms of muscle damage (Clarkson and Hubal 2002). Decreased force, as well as increased muscle soreness, swelling and/or muscle protein activity are among the most common indicators used to monitor the reparation process taking place the days following the eccentric session (Nosaka and Clarkson 1995). However, they manifest with an unsynchronized time-course after the eccentric exercise. Interestingly, eccentric exercise can be performed in different modalities, and muscle damage has been investigated using dynamic constant external load or isokinetic device (Coratella and Bertinato 2015), flywheel instruments (Coratella et al. 2016), eccentric cycling (Peñailillo et al. 2015; Mavropalias et al. 2020), and downhill running (Byrnes et al. 1985; Khassetarash et al. 2021). Particularly, downhill running places unique biomechanical, neuromuscular, and physiological challenges on the human body (Giandolini et al. 2016; Vernillo et al. 2017, 2020; Bontemps et al. 2020; Khassetarash et al. 2020, 2021). Therefore, it has raised a scientific interest given the increasing participation in trail and ultratrail running (Hoffman et al. 2010), whose downhill running seems to be a crucial aspect of the performance (Bontemps et al. 2020). Using downhill running as an eccentric exercise model could (i) inform aspects of training preparation in trail and ultratrail events, and (ii) minimize the risk of overuse injuries associated with greater impact loading and muscle damage.

While maximal force is generally considered as the gold standard to evaluate global changes in neuromuscular function, it was recently suggested that this evaluation may be incomplete, and the rise in muscle force should also be considered for a more comprehensive evaluation (Maffiuletti et al. 2016). Although instantaneous force values at specific time points could provide further insights, the rate of voluntary force development (RVFD) has been more recently adopted for a more comprehensive evaluation of a task-induced neuromuscular changes (D’Emanuele et al. 2021). The RVFD denotes the ability of the neuromuscular system to rapidly increase force over time (e.g., usually in the first 200–250 ms) (Maffiuletti et al. 2016). Specifically, the slope from the first 50–100 ms of the force–time curve has been described as the early phase, and mainly depends on neural aspects such as the firing of the motor units for the muscles involved in the task (Del Vecchio et al. 2019). In contrast, the slope from 100 to 250 ms has been described as the late phase, and depends on a series of structural factors such as

muscle morphology (Aagaard et al. 2002), muscle architecture (Coratella et al. 2020), and tendon stiffness (Maffiuletti et al. 2016).

Muscle damage is initiated by a disruption of the sarcomeres and an impairment in the excitation–contraction coupling (Proske and Morgan 2001). Both mechanisms may have implications on the RVFD, specifically on its late phase since they refer to the muscle structure instead of the neural component. After an eccentric cycling bout for 30 min at 60% of maximal concentric power output, Peñailillo et al. (2015) showed that RVFD can be considered an indirect marker of eccentric exercise-induced muscle damage. However, several physiological and biomechanical differences between running and cycling exist (Millet et al. 2009), so whether or not the RVFD may present the same behaviour after downhill running remains to be investigated. In addition, Peñailillo et al. (2015) stopped the examination at 48 h after the eccentric protocol, while it is known that muscle damage may last much further (Clarkson and Hubal 2002). Therefore, the present study aimed to examine the time-course alterations of the early and late phase of RVFD induced by a session of downhill running, together with other traditional muscle damage markers. This information could be of practical significance in trail running and following prolonged duration races wherein cumulative downhill eccentric loading is high.

## Materials and methods

### Participants

Ten healthy male recreational runners volunteered to participate in the present study (age: 22 ± 3 years, height: 181 ± 8 cm, body mass: 79 ± 8 kg) (Peñailillo et al. 2015). Participants were considered “recreational runners” if they ran less than 4 days per week (Kuru 2016; Clermont et al. 2020). 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 eccentric-based resistance training within a 6-month period prior to the start date of this study, as well as prolonged (> 20 min) and repeated (more than two 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 for the entire duration of the investigation. All procedures were approved

<figure><img src="image_2.png" /></figure>
RVFD<sub>0–50</sub> Rate of voluntary force development calculated at the 0–50 ms time interval

RVFD<sub>100–200</sub> Rate of voluntary force development calculated at the 100–200 ms time interval

## Introduction

Performing eccentric exercise is known to induce symptoms of muscle damage (Clarkson and Hubal 2002). Decreased force, as well as increased muscle soreness, swelling and/or muscle protein activity are among the most common indicators used to monitor the reparation process taking place the days following the eccentric session (Nosaka and Clarkson 1995). However, they manifest with an unsynchronized time-course after the eccentric exercise. Interestingly, eccentric exercise can be performed in different modalities, and muscle damage has been investigated using dynamic constant external load or isokinetic device (Coratella and Bertinato 2015), flywheel instruments (Coratella et al. 2016), eccentric cycling (Peñailillo et al. 2015; Mavropalias et al. 2020), and downhill running (Byrnes et al. 1985; Khassetarash et al. 2021). Particularly, downhill running places unique biomechanical, neuromuscular, and physiological challenges on the human body (Giandolini et al. 2016; Vernillo et al. 2017, 2020; Bontemps et al. 2020; Khassetarash et al. 2020, 2021). Therefore, it has raised a scientific interest given the increasing participation in trail and ultratrail running (Hoffman et al. 2010), whose downhill running seems to be a crucial aspect of the performance (Bontemps et al. 2020). Using downhill running as an eccentric exercise model could (i) inform aspects of training preparation in trail and ultratrail events, and (ii) minimize the risk of overuse injuries associated with greater impact loading and muscle damage.

While maximal force is generally considered as the gold standard to evaluate global changes in neuromuscular function, it was recently suggested that this evaluation may be incomplete, and the rise in muscle force should also be considered for a more comprehensive evaluation (Maffiuletti et al. 2016). Although instantaneous force values at specific time points could provide further insights, the rate of voluntary force development (RVFD) has been more recently adopted for a more comprehensive evaluation of a task-induced neuromuscular changes (D’Emanuele et al. 2021). The RVFD denotes the ability of the neuromuscular system to rapidly increase force over time (e.g., usually in the first 200–250 ms) (Maffiuletti et al. 2016). Specifically, the slope from the first 50–100 ms of the force–time curve has been described as the early phase, and mainly depends on neural aspects such as the firing of the motor units for the muscles involved in the task (Del Vecchio et al. 2019). In contrast, the slope from 100 to 250 ms has been described as the late phase, and depends on a series of structural factors such as

muscle morphology (Aagaard et al. 2002), muscle architecture (Coratella et al. 2020), and tendon stiffness (Maffiuletti et al. 2016).

Muscle damage is initiated by a disruption of the sarcomeres and an impairment in the excitation–contraction coupling (Proske and Morgan 2001). Both mechanisms may have implications on the RVFD, specifically on its late phase since they refer to the muscle structure instead of the neural component. After an eccentric cycling bout for 30 min at 60% of maximal concentric power output, Peñailillo et al. (2015) showed that RVFD can be considered an indirect marker of eccentric exercise-induced muscle damage. However, several physiological and biomechanical differences between running and cycling exist (Millet et al. 2009), so whether or not the RVFD may present the same behaviour after downhill running remains to be investigated. In addition, Peñailillo et al. (2015) stopped the examination at 48 h after the eccentric protocol, while it is known that muscle damage may last much further (Clarkson and Hubal 2002). Therefore, the present study aimed to examine the time-course alterations of the early and late phase of RVFD induced by a session of downhill running, together with other traditional muscle damage markers. This information could be of practical significance in trail running and following prolonged duration races wherein cumulative downhill eccentric loading is high.

## Materials and methods

### Participants

Ten healthy male recreational runners volunteered to participate in the present study (age: 22 ± 3 years, height: 181 ± 8 cm, body mass: 79 ± 8 kg) (Peñailillo et al. 2015). Participants were considered “recreational runners” if they ran less than 4 days per week (Kuru 2016; Clermont et al. 2020). 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 eccentric-based resistance training within a 6-month period prior to the start date of this study, as well as prolonged (> 20 min) and repeated (more than two 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 for the entire duration of the investigation. All procedures were approved

<figure><img src="image_3.png" /></figure>
by the local ethics committee and this study conformed to
the standards set by the Declaration of Helsinki (1965 and
further modifications), except for registration in a database.

Experimental protocol

To investigate the duration of the downhill-induced muscle
damage, the participants were tested at baseline (PRE) and
24 (POST24), 48 (POST48), 72 (POST72) and 96 (POST96)
hours after the downhill bout. Participants first came for two
familiarization sessions (separated by between 3 and 7 days),
and again 1 week after the second familiarisation session
for the experimental session. During the first familiarisation
session, participants were familiarised with the procedures
to identify the maximum voluntary isometric contraction
(MVIC) and the RVFD. During the second familiarisation
session, participants performed a single 5-min downhill run,
replicating the experimental settings on the same treadmill
used for the experimental session. Both the baseline assess-
ments and the downhill running session were conducted on
the same days. Room temperature and humidity were similar
across sessions (21 °C, 45% RH). For each assessment ses-
sion, the order was fixed as follows: venous blood samples,
muscle soreness, muscle swelling, MVIC and RVFD (Chap-
man et al. 2008).

Maximum voluntary isometric contraction

Unilateral right limb knee extensors MVIC was assessed
with the participants seated with the hip and knee angles
at 90° on a customized seat with a calibrated force trans-
ducer (Legcontrol, Mtraining, Ecole Valentin, France). In
line with a previous investigation (Varesco et al. 2022), the
lever arm of the dynamometer was adjusted to firmly attach
the leg 3 cm above the medial malleoli with two non-com-
pliant belts and a belt strapped over the waist was used to
minimize extraneous movements of the upper body. Pas-
sive resting force was subtracted from the signal so that the
baseline was set at 0 N (Varesco et al. 2022). Force signals
were collected without analogical filters at a frequency of
2 kHz by PowerLab System (16/30-ML880/P, AD Instru-
ments, Bellavista, Australia). The MVIC assessment was
preceded by a standardized warm up consisting of 20×2-s
knee extensions separated by 10 s each, and the partici-
pants were instructed to gradually increase the force up to
the maximum volitional force (Coratella et al. 2020). After
5 min, two 4-s MVIC attempts were performed, separated by
5 min of passive rest. At baseline, if the difference between
the two MVICs was > 5%, further trials were performed until
the difference between two consecutive trials was < 5%. The
participants were instructed to “push as hard as possible”
and strong standardized encouragements were provided by

the operators at each attempt. The maximum value was then
retained for further analysis.

Rate of voluntary force development

The RVFD was assessed after 10 min of passive rest on the
same customized seat and with the participants in the same
position as for the MVIC. The participants were requested to
“push as fast as possible” during eight impulsive unilateral
right limb isometric knee extensions (Varesco et al. 2019).
In case of a countermovement (determined by a force drop
of 2 N below the baseline right before the impulsive con-
traction) 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. 2022). All testing protocols were performed with real-
time visual feedback. For each impulsive contraction, we
computed the first derivative of the force–time curve, and
retrieved the peak value as the point of maximal accelera-
tion during the contraction. We then extracted and average
value from data 5 ms before to 5 ms after this data point
to compute peak RVFD value (Varesco et al. 2019). The
onset of the voluntary force development was automatically
defined as the point at which force exceeded the average rest-
ing baseline by ~ 2 N. The onset was also checked visually by
an experienced operator blinded to the condition. The force
was measured at 50 and 200 ms and the RVFD was calcu-
lated fitting a linear model over the 0–50 ms (RVFD<sub>0–50</sub>)
and 100–200 ms (RVFD<sub>100–200</sub>) time-windows data points
and extracting the slope of the model (Fig. 1). All data were
analysed offline using Labchart 8 Software (ADInstruments,
Bella Vista, Australia) and then retained for further analysis.

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Representative participant's typical traces of the rate of voluntary force development (RVFD) at different time intervals from the onset before (PRE) and 24 h after (POST24) the downhill bout</figcaption></figure>
by the local ethics committee and this study conformed to
the standards set by the Declaration of Helsinki (1965 and
further modifications), except for registration in a database.

Experimental protocol

To investigate the duration of the downhill-induced muscle
damage, the participants were tested at baseline (PRE) and
24 (POST24), 48 (POST48), 72 (POST72) and 96 (POST96)
hours after the downhill bout. Participants first came for two
familiarization sessions (separated by between 3 and 7 days),
and again 1 week after the second familiarisation session
for the experimental session. During the first familiarisation
session, participants were familiarised with the procedures
to identify the maximum voluntary isometric contraction
(MVIC) and the RVFD. During the second familiarisation
session, participants performed a single 5-min downhill run,
replicating the experimental settings on the same treadmill
used for the experimental session. Both the baseline assess-
ments and the downhill running session were conducted on
the same days. Room temperature and humidity were similar
across sessions (21 °C, 45% RH). For each assessment ses-
sion, the order was fixed as follows: venous blood samples,
muscle soreness, muscle swelling, MVIC and RVFD (Chap-
man et al. 2008).

Maximum voluntary isometric contraction

Unilateral right limb knee extensors MVIC was assessed
with the participants seated with the hip and knee angles
at 90° on a customized seat with a calibrated force trans-
ducer (Legcontrol, Mtraining, Ecole Valentin, France). In
line with a previous investigation (Varesco et al. 2022), the
lever arm of the dynamometer was adjusted to firmly attach
the leg 3 cm above the medial malleoli with two non-com-
pliant belts and a belt strapped over the waist was used to
minimize extraneous movements of the upper body. Pas-
sive resting force was subtracted from the signal so that the
baseline was set at 0 N (Varesco et al. 2022). Force signals
were collected without analogical filters at a frequency of
2 kHz by PowerLab System (16/30-ML880/P, AD Instru-
ments, Bellavista, Australia). The MVIC assessment was
preceded by a standardized warm up consisting of 20×2-s
knee extensions separated by 10 s each, and the partici-
pants were instructed to gradually increase the force up to
the maximum volitional force (Coratella et al. 2020). After
5 min, two 4-s MVIC attempts were performed, separated by
5 min of passive rest. At baseline, if the difference between
the two MVICs was > 5%, further trials were performed until
the difference between two consecutive trials was < 5%. The
participants were instructed to “push as hard as possible”
and strong standardized encouragements were provided by

the operators at each attempt. The maximum value was then
retained for further analysis.

Rate of voluntary force development

The RVFD was assessed after 10 min of passive rest on the
same customized seat and with the participants in the same
position as for the MVIC. The participants were requested to
“push as fast as possible” during eight impulsive unilateral
right limb isometric knee extensions (Varesco et al. 2019).
In case of a countermovement (determined by a force drop
of 2 N below the baseline right before the impulsive con-
traction) 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. 2022). All testing protocols were performed with real-
time visual feedback. For each impulsive contraction, we
computed the first derivative of the force–time curve, and
retrieved the peak value as the point of maximal accelera-
tion during the contraction. We then extracted and average
value from data 5 ms before to 5 ms after this data point
to compute peak RVFD value (Varesco et al. 2019). The
onset of the voluntary force development was automatically
defined as the point at which force exceeded the average rest-
ing baseline by ~ 2 N. The onset was also checked visually by
an experienced operator blinded to the condition. The force
was measured at 50 and 200 ms and the RVFD was calcu-
lated fitting a linear model over the 0–50 ms (RVFD<sub>0–50</sub>)
and 100–200 ms (RVFD<sub>100–200</sub>) time-windows data points
and extracting the slope of the model (Fig. 1). All data were
analysed offline using Labchart 8 Software (ADInstruments,
Bella Vista, Australia) and then retained for further analysis.

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Representative participant’s typical traces of the rate of voluntary force development (RVFD) at different time intervals from the onset before (PRE) and 24 h after (POST24) the downhill bout</figcaption></figure>
## Serum creatine kinase

Venous blood samples (~5 mL) were collected by standard venipuncture of the antecubital vein in SST II Advance Vacutainer® and in K₂EDTA tubes (Becton, Dickinson & Co., Franklin Lakes, NJ, USA). Then it was immediately centrifuged at 3000×g (15 min, 4 °C) and serum and plasma was aliquoted and stored at –80 °C until assayed. Serum concentration of muscle creatine kinase isoform (CKM), including skeletal (CK-MM) and cardiac (CK-MB) muscle isoenzymes, was measured by a SimpleStep ELISA® (Abcam, Cambridge, UK). Sensitivity, intra-assay (CVw), and inter-assay (CVb) variations were 270 pg/mL 2.2% and 7.4%, respectively.

## Muscle swelling

To evaluate muscle swelling, the right thigh circumference was measured at 50% of the distance between the iliac spine and the patella using a seamstress meter with the participant seated, knee at 90° flexion (with 0° being full extension) and the lower limbs relaxed. To ensure consistency, the site was marked with a semipermanent ink. Three measurements were taken, and the mean was then retained for further analysis.

## Muscle soreness

Muscle soreness was investigated as pressure pain threshold using an algometer (Pain Diagnostics and Thermography, New York, USA) on the right vastus medialis at 50% of the distance between the anterior superior iliac spine and the patella. The probe was placed perpendicular to the site and the investigator gradually applied force until the participant reported pain (Fischer 1987). Three measurements were performed, and the mean was then retained for further analysis.

## Downhill running

The participants performed a 30-min downhill run at 10 km h⁻¹ and with a slope of –20% (–11.3°) on a motorized treadmill (Medic 2855, Genin Medical, La Roque-d'Anthéron, France) (Varesco et al. 2022). We used treadmill running given the similarity in running biomechanics observed with the downhill overground running (Firminger et al. 2018).

## Statistical analysis

Results are given as means ± standard deviations. The normality of distribution was verified using the Shapiro–Wilk normality test. To test differences between PRE and POST24/48/72/96, a longitudinal analysis was

performed using generalized estimating equations (GEE; i.e., GEE under ‘generalized linear model’ procedure in SPSS v. 28) to take into account the correlated nature of observations within each participant (i.e., within-participant measurements) (Liang and Zeger 1986). If a significant main effect for time was observed, Bonferroni’s test was used for post hoc analysis. As a measure of effect size, Cohen’s d (d) was calculated and interpreted as follows: <0.19=trivial, 0.20–0.59=small, 0.60–1.19=moderate, 1.20–1.99=larg e, >2.00=very large (Hopkins et al. 2009). The statistical analyses were conducted using IBM™ SPSS™ Statistics (version 28.0.0; IBM Corp., Somers, New York, NY) with the criterion α level set to 0.05.

## Results

### Rate of voluntary force development

Force at 50 ms did not show a time effect [$_{w}$χ²(4)=2.0, p=0.732] (Fig. 2A). Force at 200 ms showed a time effect [$_{w}$χ²(4)=16.7, p=0.002] (Fig. 2B). Force at 200 ms was only lower at POST48 compared to PRE (467±132 N vs 375±113 N, 83±20% of PRE values, p=0.045, d=0.75). RVFD₀⁻₅₀ showed a time effect [$_{w}$χ²(4)=9.8, p=0.044]. However, when PRE was compared to the other time points, no difference was detected (all p=1.000) (Fig. 2C). RVFD₁₀₀⁻₂₀₀ showed a time effect [$_{w}$χ²(4)=42.2, p<0.001]. Compared to PRE (1339±418 N s⁻¹), RFD100–200 was lower at POST24 (894±559 N s⁻¹, 63±28% of PRE values, p<0.001, d=0.90), and remained lower until POST72 (993±568 N s⁻¹, 74±31% of PRE values, p=0.030, d=0.69) (Fig. 2D).

### Maximum voluntary isometric contraction

MVIC force showed a time effect [$_{w}$χ²(4)=17.7, p=0.001]. Compared to PRE (608±156 N), MVIC force was lower at POST24 (500±109 N, 84±13% of PRE values, p=0.002, d=0.80), and remained lower until POST72 (534±118 N, 90±13% of PRE values, p=0.026, d=0.53) (Fig. 3A).

## Serum creatine kinase

Creatine kinase (CK) showed a time effect [$_{w}$χ²(4)=46.2, p<0.001]. Compared to PRE (32,923±23,874 pg mL⁻¹), CK was higher at POST24 (147,273±55,447 pg mL⁻¹, 613±394% of PRE values, p<0.001, d=2.68), and POST48 (69,495±41,670 pg mL⁻¹, 232±118% of PRE values, p<0.001, d=1.08) (Fig. 3B).
## Serum creatine kinase

Venous blood samples (~5 mL) were collected by standard venipuncture of the antecubital vein in SST II Advance Vacutainer® and in K₂EDTA tubes (Becton, Dickinson & Co., Franklin Lakes, NJ, USA). Then it was immediately centrifuged at 3000×g (15 min, 4 °C) and serum and plasma was aliquoted and stored at –80 °C until assayed. Serum concentration of muscle creatine kinase isoform (CKM), including skeletal (CK-MM) and cardiac (CK-MB) muscle isoenzymes, was measured by a SimpleStep ELISA® (Abcam, Cambridge, UK). Sensitivity, intra-assay (CVw), and inter-assay (CVb) variations were 270 pg/mL 2.2% and 7.4%, respectively.

## Muscle swelling

To evaluate muscle swelling, the right thigh circumference was measured at 50% of the distance between the iliac spine and the patella using a seamstress meter with the participant seated, knee at 90° flexion (with 0° being full extension) and the lower limbs relaxed. To ensure consistency, the site was marked with a semipermanent ink. Three measurements were taken, and the mean was then retained for further analysis.

## Muscle soreness

Muscle soreness was investigated as pressure pain threshold using an algometer (Pain Diagnostics and Thermography, New York, USA) on the right vastus medialis at 50% of the distance between the anterior superior iliac spine and the patella. The probe was placed perpendicular to the site and the investigator gradually applied force until the participant reported pain (Fischer 1987). Three measurements were performed, and the mean was then retained for further analysis.

## Downhill running

The participants performed a 30-min downhill run at 10 km h⁻¹ and with a slope of –20% (–11.3°) on a motorized treadmill (Medic 2855, Genin Medical, La Roque-d'Anthéron, France) (Varesco et al. 2022). We used treadmill running given the similarity in running biomechanics observed with the downhill overground running (Firminger et al. 2018).

## Statistical analysis

Results are given as means ± standard deviations. The normality of distribution was verified using the Shapiro–Wilk normality test. To test differences between PRE and POST24/48/72/96, a longitudinal analysis was

performed using generalized estimating equations (GEE; i.e., GEE under ‘generalized linear model’ procedure in SPSS v. 28) to take into account the correlated nature of observations within each participant (i.e., within-participant measurements) (Liang and Zeger 1986). If a significant main effect for time was observed, Bonferroni’s test was used for post hoc analysis. As a measure of effect size, Cohen’s d (d) was calculated and interpreted as follows: <0.19=trivial, 0.20–0.59=small, 0.60–1.19=moderate, 1.20–1.99=larg e, >2.00=very large (Hopkins et al. 2009). The statistical analyses were conducted using IBM™ SPSS™ Statistics (version 28.0.0; IBM Corp., Somers, New York, NY) with the criterion α level set to 0.05.

## Results

### Rate of voluntary force development

Force at 50 ms did not show a time effect [$_{w}$χ²(4)=2.0, p=0.732] (Fig. 2A). Force at 200 ms showed a time effect [$_{w}$χ²(4)=16.7, p=0.002] (Fig. 2B). Force at 200 ms was only lower at POST48 compared to PRE (467±132 N vs 375±113 N, 83±20% of PRE values, p=0.045, d=0.75). RVFD₀⁻₅₀ showed a time effect [$_{w}$χ²(4)=9.8, p=0.044]. However, when PRE was compared to the other time points, no difference was detected (all p=1.000) (Fig. 2C). RVFD₁₀₀⁻₂₀₀ showed a time effect [$_{w}$χ²(4)=42.2, p<0.001]. Compared to PRE (1339±418 N s⁻¹), RFD100–200 was lower at POST24 (894±559 N s⁻¹, 63±28% of PRE values, p<0.001, d=0.90), and remained lower until POST72 (993±568 N s⁻¹, 74±31% of PRE values, p=0.030, d=0.69) (Fig. 2D).

### Maximum voluntary isometric contraction

MVIC force showed a time effect [$_{w}$χ²(4)=17.7, p=0.001]. Compared to PRE (608±156 N), MVIC force was lower at POST24 (500±109 N, 84±13% of PRE values, p=0.002, d=0.80), and remained lower until POST72 (534±118 N, 90±13% of PRE values, p=0.026, d=0.53) (Fig. 3A).

## Serum creatine kinase

Creatine kinase (CK) showed a time effect [$_{w}$χ²(4)=46.2, p<0.001]. Compared to PRE (32,923±23,874 pg mL⁻¹), CK was higher at POST24 (147,273±55,447 pg mL⁻¹, 613±394% of PRE values, p<0.001, d=2.68), and POST48 (69,495±41,670 pg mL⁻¹, 232±118% of PRE values, p<0.001, d=1.08) (Fig. 3B).
<figure><img src="image_4.png" /><figcaption>**Fig. 2** The time-course of the force exerted at 50 ms (A) and 200 ms (B), and of the early (C) and the late phase (D) of the knee extensors RVFD. \*p<0.05 vs PRE</figcaption></figure>

## Muscle swelling

Muscle swelling showed a time effect [$\omega\chi^2$ (4) = 55.0, p<0.001]. Compared to PRE (51 ± 3 cm), muscle swelling was higher at POST24 (52 ± 2 cm, 102 ± 1% of PRE values, p<0.001, d=0.43), and remained higher until POST72 (51 ± 3 cm, 102 ± 2% of PRE values, p=0.005, d=0.32) (Fig. 3C).

## Muscle soreness

Muscle soreness showed a time effect [$\omega\chi^2$ (4) = 19.0, p<0.001]. Compared to PRE (7.7 ± 2.5 kPa), muscle soreness decreased at POST48 (5.5 ± 1.7 kPa, 77 ± 27% of PRE

values, p=0.007, d=1.03), and POST72 (6.0 ± 1.8 kPa, 81 ± 18% of PRE values, p=0.005, d=0.78) (Fig. 3D).

## Discussion

The present study was conceived to determine whether RVFD is sensitive to the muscle damage following downhill running. The effectiveness of the downhill running bout in inducing muscle damage was confirmed by a reduction in the MVIC, as well as by an increase in muscle swelling, CK concentration, and muscle soreness. Overall, our results showed that a downhill running session impaired the RVFD<sub>100-200</sub> for 3 days after the bout, while the RVFD<sub>0-50</sub>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** The time-course of the force exerted at 50 ms (A) and 200 ms (B), and of the early (C) and the late phase (D) of the knee extensors RVFD. \*p<0.05 vs PRE</figcaption></figure>

## Muscle swelling

Muscle swelling showed a time effect [$\omega\chi^2$ (4) = 55.0, p<0.001]. Compared to PRE (51 ± 3 cm), muscle swelling was higher at POST24 (52 ± 2 cm, 102 ± 1% of PRE values, p<0.001, d=0.43), and remained higher until POST72 (51 ± 3 cm, 102 ± 2% of PRE values, p=0.005, d=0.32) (Fig. 3C).

## Muscle soreness

Muscle soreness showed a time effect [$\omega\chi^2$ (4) = 19.0, p<0.001]. Compared to PRE (7.7 ± 2.5 kPa), muscle soreness decreased at POST48 (5.5 ± 1.7 kPa, 77 ± 27% of PRE

values, p=0.007, d=1.03), and POST72 (6.0 ± 1.8 kPa, 81 ± 18% of PRE values, p=0.005, d=0.78) (Fig. 3D).

## Discussion

The present study was conceived to determine whether RVFD is sensitive to the muscle damage following downhill running. The effectiveness of the downhill running bout in inducing muscle damage was confirmed by a reduction in the MVIC, as well as by an increase in muscle swelling, CK concentration, and muscle soreness. Overall, our results showed that a downhill running session impaired the RVFD<sub>100-200</sub> for 3 days after the bout, while the RVFD<sub>0-50</sub>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** The time-course of the maximum voluntary isometric contraction (MVIC) (A), serum creatine kinase (CK) concentration (B), muscle swelling (C) and soreness (algometry) (D). \*p<0.05 vs PRE</figcaption></figure>

was not affected. Taken together, the present results suggest that the late phase of RVFD appears sensitive to the muscle damage induced after an eccentric-biased downhill running, while the early phase does not.

We found moderate decrements for the RVFD<sub>100-200</sub>, while the RVFD<sub>0-50</sub> was not affected by a bout of downhill running. While the early phase is mostly determined by neural factors, the late phase of RVFD is mainly based on muscular determinants (Maffiuletti et al. 2016). Among others, the type-II fibres have been pointed as more advantageous for producing faster torque than type-I fibres because of a greater amount of Ca<sup>2+</sup> released per action potential, faster muscle proteins isoforms and cross-bridge formation rates (Maffiuletti et al. 2016). This had implications on the skeletal muscle capacity to produce rapid torque (Harridge et al. 1996). On the other hand, type-II fibres are more sensitive to muscle damage as compared to type-I fibres (Lieber

and Friden 1988). Indeed, greater muscle damage was observed when comparing the same relative eccentric exercise performed by elbow flexors vs leg extensors (Jamurtas et al. 2005), with the former showing greater type-II fibre prevalence than the latter (Saltin and Gollnick 2011). To our knowledge, the only study that has previously assessed RVFD following a downhill running bout did not distinguish the early from the late phase, while it reported the force exerted at 30 ms, 50 ms, 100 ms and 200 ms (Maeo et al. 2017). Notwithstanding, this was investigated in both the knee extensors and the plantar flexors, opening for interesting comparisons (Maeo et al. 2017). For instance, the authors observed that the force exerted at 200 ms was impaired in the knee extensors but not in the plantar flexors, in line with the possible relationship between the type-II prevalence and the late phase (Maeo et al. 2017). However, the authors also found impairments in the force exerted at

<figure><img src="image_6.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** The time-course of the maximum voluntary isometric contraction (MVIC) (A), serum creatine kinase (CK) concentration (B), muscle swelling (C) and soreness (algometry) (D). \*p<0.05 vs PRE</figcaption></figure>

was not affected. Taken together, the present results suggest that the late phase of RVFD appears sensitive to the muscle damage induced after an eccentric-biased downhill running, while the early phase does not.

We found moderate decrements for the RVFD<sub>100-200</sub>, while the RVFD<sub>0-50</sub> was not affected by a bout of downhill running. While the early phase is mostly determined by neural factors, the late phase of RVFD is mainly based on muscular determinants (Maffiuletti et al. 2016). Among others, the type-II fibres have been pointed as more advantageous for producing faster torque than type-I fibres because of a greater amount of Ca<sup>2+</sup> released per action potential, faster muscle proteins isoforms and cross-bridge formation rates (Maffiuletti et al. 2016). This had implications on the skeletal muscle capacity to produce rapid torque (Harridge et al. 1996). On the other hand, type-II fibres are more sensitive to muscle damage as compared to type-I fibres (Lieber

and Friden 1988). Indeed, greater muscle damage was observed when comparing the same relative eccentric exercise performed by elbow flexors vs leg extensors (Jamurtas et al. 2005), with the former showing greater type-II fibre prevalence than the latter (Saltin and Gollnick 2011). To our knowledge, the only study that has previously assessed RVFD following a downhill running bout did not distinguish the early from the late phase, while it reported the force exerted at 30 ms, 50 ms, 100 ms and 200 ms (Maeo et al. 2017). Notwithstanding, this was investigated in both the knee extensors and the plantar flexors, opening for interesting comparisons (Maeo et al. 2017). For instance, the authors observed that the force exerted at 200 ms was impaired in the knee extensors but not in the plantar flexors, in line with the possible relationship between the type-II prevalence and the late phase (Maeo et al. 2017). However, the authors also found impairments in the force exerted at

<figure><img src="image_7.png" /></figure>
50 ms (Maeo et al. 2017), which contrasts with what we observed in the present study. However, besides the different downhill running protocols [45 min at -15% (Maeo et al. 2017) vs 30 min at -20% in the present study), our participants were recreational runners, while the others (Maeo et al. 2017) were described as young adults with no specific experience in running. Since the extent of muscle damage is greater in untrained vs trained individuals (Newton et al. 2008), it is possible that those participants may have experienced greater muscle damage, also affecting the neural aspects as already seen after eccentric exercise (Prasartwuth et al. 2005). Summarizing, the muscle damage following the downhill running bout may have affected the muscle structure and especially the type-II fibres, whose characteristics are important determinants of the RVFD late phase (Peñailillo et al. 2015).

The exposure to downhill running is followed by alterations in muscle damage markers. We observed decrements in knee extensors MVIC that recovered after 4 days, coupled with increments in CK blood concentration, muscle swelling and delayed soreness. The decrement in maximum force derives from both neural and structural factors, while the increment in CK blood concentration is an index of the rupture of the myonuclei, and the increment in muscle swelling and soreness is mainly due to the inflammation process following the eccentric exercise (Proske and Morgan 2001).

## Limitations

First, the present results refer to downhill running, and other eccentric exercise modalities may produce different effects. Second, different durations and/or slopes, as well as different populations may result in dissimilar outcomes. Third, the sample size was small and should be increased for more definitive conclusions. Fourth, further muscle damage markers such as the echo-intensity or the shift in optimum angle could have deepened the discussion. Last, we only investigated the knee extensors, and different muscles in combination with different eccentric exercises may differ in behaviour.

## Conclusions

A downhill running bout impaired the late but not the early phase of the knee extensors RVFD for the following 72 h, with concomitant decrements in knee extensors maximum force and increments in serum CK concentration, muscle swelling and soreness. Overall, all muscle damage markers recovered in 4 days. The rate of voluntary force development, and especially its late phase, appears sensitive to muscle damage and could be used as an additional muscle

damage marker in trail running and following prolonged duration races wherein cumulative eccentric loading is high.

**Acknowledgements** The authors thank the participants for their time and effort. They also thank Léo Vanini for his valuable technical assistance during data acquisition, and Hayley L. Patmore for her assistance in English proof reading.

**Authors' contributions** GC, GVa, GL, GVe, and LM conceived and designed research; GVa, and BC performed experiments; GC, GVa, VS, GL, and GVe interpreted results of experiments; GVa, VS, GL, analyzed data; GC, and GVe prepared figures; GC drafted manuscript; GC, GVa, VR, BC, VS, GL, GVe, and LM edited and revised manuscript; GC, GVa, VR, BC, VS, GL, GVe, and LM approved final version of manuscript.

**Funding** Open access funding provided by Università degli Studi di Milano within the CRUI-CARE Agreement. This research was supported by the Ultra-Endurance Sports Science & Medicine Grant Program.

**Data availability** Derived data supporting the findings of this study are available on request from the corresponding author.

**Code availability** Not applicable.

## Declarations

**Conflicts 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 the study.

**Consent for publication** Not applicable.

**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/.

## 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

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
50 ms (Maeo et al. 2017), which contrasts with what we observed in the present study. However, besides the different downhill running protocols [45 min at -15% (Maeo et al. 2017) vs 30 min at -20% in the present study), our participants were recreational runners, while the others (Maeo et al. 2017) were described as young adults with no specific experience in running. Since the extent of muscle damage is greater in untrained vs trained individuals (Newton et al. 2008), it is possible that those participants may have experienced greater muscle damage, also affecting the neural aspects as already seen after eccentric exercise (Prasartwuth et al. 2005). Summarizing, the muscle damage following the downhill running bout may have affected the muscle structure and especially the type-II fibres, whose characteristics are important determinants of the RVFD late phase (Peñailillo et al. 2015).

The exposure to downhill running is followed by alterations in muscle damage markers. We observed decrements in knee extensors MVIC that recovered after 4 days, coupled with increments in CK blood concentration, muscle swelling and delayed soreness. The decrement in maximum force derives from both neural and structural factors, while the increment in CK blood concentration is an index of the rupture of the myonuclei, and the increment in muscle swelling and soreness is mainly due to the inflammation process following the eccentric exercise (Proske and Morgan 2001).

## Limitations

First, the present results refer to downhill running, and other eccentric exercise modalities may produce different effects. Second, different durations and/or slopes, as well as different populations may result in dissimilar outcomes. Third, the sample size was small and should be increased for more definitive conclusions. Fourth, further muscle damage markers such as the echo-intensity or the shift in optimum angle could have deepened the discussion. Last, we only investigated the knee extensors, and different muscles in combination with different eccentric exercises may differ in behaviour.

## Conclusions

A downhill running bout impaired the late but not the early phase of the knee extensors RVFD for the following 72 h, with concomitant decrements in knee extensors maximum force and increments in serum CK concentration, muscle swelling and soreness. Overall, all muscle damage markers recovered in 4 days. The rate of voluntary force development, and especially its late phase, appears sensitive to muscle damage and could be used as an additional muscle

damage marker in trail running and following prolonged duration races wherein cumulative eccentric loading is high.

**Acknowledgements** The authors thank the participants for their time and effort. They also thank Léo Vanini for his valuable technical assistance during data acquisition, and Hayley L. Patmore for her assistance in English proof reading.

**Authors' contributions** GC, GVa, GL, GVe, and LM conceived and designed research; GVa, and BC performed experiments; GC, GVa, VS, GL, and GVe interpreted results of experiments; GVa, VS, GL, analyzed data; GC, and GVe prepared figures; GC drafted manuscript; GC, GVa, VR, BC, VS, GL, GVe, and LM edited and revised manuscript; GC, GVa, VR, BC, VS, GL, GVe, and LM approved final version of manuscript.

**Funding** Open access funding provided by Università degli Studi di Milano within the CRUI-CARE Agreement. This research was supported by the Ultra-Endurance Sports Science & Medicine Grant Program.

**Data availability** Derived data supporting the findings of this study are available on request from the corresponding author.

**Code availability** Not applicable.

## Declarations

**Conflicts 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 the study.

**Consent for publication** Not applicable.

**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/.

## 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

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
Byrnes WC, Clarkson PM, White JS, Hsieh SS, Frykman PN, Maughan RJ (1985) Delayed onset muscle soreness following repeated bouts of downhill running. J Appl Physiol (1985) 59(3):710–715. https://doi.org/10.1152/jappl.1985.59.3.710

Chapman DW, Newton M, McGuigan MR, Nosaka K (2008) Comparison between old and young men for responses to fast velocity maximal lengthening contractions of the elbow flexors. Eur J Appl Physiol 104(3):531–539. https://doi.org/10.1007/s00421-008-0806-7

Clarkson PM, Hubal MJ (2002) Exercise-induced muscle damage in humans. Am J Phys Med Rehabil 81(11 Suppl):S52–S69. https://doi.org/10.1097/00002060-200211001-00007

Clermont CA, Duffett-Leger L, Hettinga BA, Ferber R (2020) Runners’ perspectives on ‘smart’ wearable technology and its use for preventing injury. Int J Hum Comput Int 36(1):31–40. https://doi.org/10.1080/10447318.2019.1597575

Coratella G, Bertinato L (2015) Isoload vs isokinetic eccentric exercise: a direct comparison of exercise-induced muscle damage and repeated bout effect. Sport Sci Health 11(1):87–96. https://doi.org/10.1007/s11332-014-0213-x

Coratella G, Chemello A, Schena F (2016) Muscle damage and repeated bout effect induced by enhanced eccentric squats. J Sports Med Phys Fitness 56(12):1540–1546

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

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

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

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

Fischer AA (1987) Pressure algometry over normal muscles. Standard values, validity and reproducibility of pressure threshold. Pain 30(1):115–126. https://doi.org/10.1016/0304-3959(87)90089-3

Giandolini M, Vernillo G, Samozino P, Horvais N, Edwards WB, Morin JB, Millet GY (2016) Fatigue associated with prolonged graded running. Eur J Appl Physiol 116(10):1859–1873. https://doi.org/10.1007/s00421-016-3437-4

Harridge SD, Bottinelli R, Canepari M, Pellegrino MA, Reggiani C, Esb-jornsson M, Saltin B (1996) Whole-muscle and single-fibre contractile properties and myosin heavy chain isoforms in humans. Pflugers Arch 432(5):913–920. https://doi.org/10.1007/s004240050215

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

Hopkins WG, Marshall SW, Batterham AM, Hanin J (2009) Progressive statistics for studies in sports medicine and exercise science. Med Sci Sports Exerc 41(1):3–13. https://doi.org/10.1249/MSS.0b013e31818cb278

Jamurtas AZ, Theocharis V, Tofas T, Tsiokanos A, Yfanti C, Paschalis V, Koutedakis Y, Nosaka K (2005) Comparison between leg and arm eccentric exercises of the same relative intensity on indices of muscle damage. Eur J Appl Physiol 95(2–3):179–185. https://doi.org/10.1007/s00421-005-1345-0

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

Kuru A (2016) Exploring experience of runners with sports tracking technology. Int J Hum Comput Int 32(11):847–860. https://doi.org/10.1080/10447318.2016.120246

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

Lieber RL, Friden J (1988) Selective damage of fast glycolytic muscle fibres with eccentric contraction of the rabbit tibialis anterior. Acta Physiol Scand 133(4):587–588. https://doi.org/10.1111/j.1748-1716.1988.tb08446.x

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

Mavropalias G, Koeda T, Barley OR, Poon WCK, Fisher AJ, Blazevich AJ, Nosaka K (2020) Comparison between high- and low-intensity eccentric cycling of equal mechanical work for muscle damage and the repeated bout effect. Eur J Appl Physiol 120(5):1015–1025. https://doi.org/10.1007/s00421-020-04341-5

Millet GP, Vleck VE, Bentley DJ (2009) Physiological differences between cycling and running: lessons from triathletes. Sports Med 39(3):179–206. https://doi.org/10.2165/00007256-200993030-00002

Newton MJ, Morgan GT, Sacco P, Chapman DW, Nosaka K (2008) Comparison of responses to strenuous eccentric exercise of the elbow flexors between resistance-trained and untrained men. J Strength Cond Res 22(2):597–607. https://doi.org/10.1519/JSC.0b013e3181660003

Nosaka K, Clarkson PM (1995) Muscle damage following repeated bouts of high force eccentric exercise. Med Sci Sports Exerc 27(9):1263–1269

Peñailillo L, Blazevich A, Numazawa H, Nosaka K (2015) Rate of force development as a measure of muscle damage. Scand J Med Sci Sports 25(3):417–427. https://doi.org/10.1111/sms.12241

Prasartwuth O, Taylor JL, Gandevia SC (2005) Maximal force, voluntary activation and muscle soreness after eccentric damage to human elbow flexor muscles. J Physiol 567(Pt 1):337–348. https://doi.org/10.1113/jphysiol.2005.087767

Prospe U, Morgan DL (2001) Muscle damage from eccentric exercise: mechanism, mechanical signs, adaptation and clinical applications. J Physiol 537(Pt 2):333–345. https://doi.org/10.1111/j.1469-7793.2001.00333.x

Saltin B, Gollnick PD (2011) Skeletal muscle adaptability: significance for metabolism and performance. Compr Physiol 555–631

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

Varesco G, Coratella G, Rozand V, Cuinet B, Lombardi G, Mourot L, Vernillo G (2022) Downhill running affects the late but not the early phase of the rate of force development. Eur J Appl Physiol 122(9):2049–2059. https://doi.org/10.1007/s00421-022-04990-8

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
Byrnes WC, Clarkson PM, White JS, Hsieh SS, Frykman PN, Maughan RJ (1985) Delayed onset muscle soreness following repeated bouts of downhill running. J Appl Physiol (1985) 59(3):710–715. https://doi.org/10.1152/jappl.1985.59.3.710

Chapman DW, Newton M, McGuigan MR, Nosaka K (2008) Comparison between old and young men for responses to fast velocity maximal lengthening contractions of the elbow flexors. Eur J Appl Physiol 104(3):531–539. https://doi.org/10.1007/s00421-008-0806-7

Clarkson PM, Hubal MJ (2002) Exercise-induced muscle damage in humans. Am J Phys Med Rehabil 81(11 Suppl):S52–S69. https://doi.org/10.1097/00002060-200211001-00007

Clermont CA, Duffett-Leger L, Hettinga BA, Ferber R (2020) Runners’ perspectives on ‘smart’ wearable technology and its use for preventing injury. Int J Hum Comput Int 36(1):31–40. https://doi.org/10.1080/10447318.2019.1597575

Coratella G, Bertinato L (2015) Isoload vs isokinetic eccentric exercise: a direct comparison of exercise-induced muscle damage and repeated bout effect. Sport Sci Health 11(1):87–96. https://doi.org/10.1007/s11332-014-0213-x

Coratella G, Chemello A, Schena F (2016) Muscle damage and repeated bout effect induced by enhanced eccentric squats. J Sports Med Phys Fitness 56(12):1540–1546

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

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

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

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

Fischer AA (1987) Pressure algometry over normal muscles. Standard values, validity and reproducibility of pressure threshold. Pain 30(1):115–126. https://doi.org/10.1016/0304-3959(87)90089-3

Giandolini M, Vernillo G, Samozino P, Horvais N, Edwards WB, Morin JB, Millet GY (2016) Fatigue associated with prolonged graded running. Eur J Appl Physiol 116(10):1859–1873. https://doi.org/10.1007/s00421-016-3437-4

Harridge SD, Bottinelli R, Canepari M, Pellegrino MA, Reggiani C, Esb-jornsson M, Saltin B (1996) Whole-muscle and single-fibre contractile properties and myosin heavy chain isoforms in humans. Pflugers Arch 432(5):913–920. https://doi.org/10.1007/s004240050215

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

Hopkins WG, Marshall SW, Batterham AM, Hanin J (2009) Progressive statistics for studies in sports medicine and exercise science. Med Sci Sports Exerc 41(1):3–13. https://doi.org/10.1249/MSS.0b013e31818cb278

Jamurtas AZ, Theocharis V, Tofas T, Tsiokanos A, Yfanti C, Paschalis V, Koutedakis Y, Nosaka K (2005) Comparison between leg and arm eccentric exercises of the same relative intensity on indices of muscle damage. Eur J Appl Physiol 95(2–3):179–185. https://doi.org/10.1007/s00421-005-1345-0

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

Kuru A (2016) Exploring experience of runners with sports tracking technology. Int J Hum Comput Int 32(11):847–860. https://doi.org/10.1080/10447318.2016.120246

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

Lieber RL, Friden J (1988) Selective damage of fast glycolytic muscle fibres with eccentric contraction of the rabbit tibialis anterior. Acta Physiol Scand 133(4):587–588. https://doi.org/10.1111/j.1748-1716.1988.tb08446.x

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

Mavropalias G, Koeda T, Barley OR, Poon WCK, Fisher AJ, Blazevich AJ, Nosaka K (2020) Comparison between high- and low-intensity eccentric cycling of equal mechanical work for muscle damage and the repeated bout effect. Eur J Appl Physiol 120(5):1015–1025. https://doi.org/10.1007/s00421-020-04341-5

Millet GP, Vleck VE, Bentley DJ (2009) Physiological differences between cycling and running: lessons from triathletes. Sports Med 39(3):179–206. https://doi.org/10.2165/00007256-200993030-00002

Newton MJ, Morgan GT, Sacco P, Chapman DW, Nosaka K (2008) Comparison of responses to strenuous eccentric exercise of the elbow flexors between resistance-trained and untrained men. J Strength Cond Res 22(2):597–607. https://doi.org/10.1519/JSC.0b013e3181660003

Nosaka K, Clarkson PM (1995) Muscle damage following repeated bouts of high force eccentric exercise. Med Sci Sports Exerc 27(9):1263–1269

Peñailillo L, Blazevich A, Numazawa H, Nosaka K (2015) Rate of force development as a measure of muscle damage. Scand J Med Sci Sports 25(3):417–427. https://doi.org/10.1111/sms.12241

Prasartwuth O, Taylor JL, Gandevia SC (2005) Maximal force, voluntary activation and muscle soreness after eccentric damage to human elbow flexor muscles. J Physiol 567(Pt 1):337–348. https://doi.org/10.1113/jphysiol.2005.087767

Prospe U, Morgan DL (2001) Muscle damage from eccentric exercise: mechanism, mechanical signs, adaptation and clinical applications. J Physiol 537(Pt 2):333–345. https://doi.org/10.1111/j.1469-7793.2001.00333.x

Saltin B, Gollnick PD (2011) Skeletal muscle adaptability: significance for metabolism and performance. Compr Physiol 555–631

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

Varesco G, Coratella G, Rozand V, Cuinet B, Lombardi G, Mourot L, Vernillo G (2022) Downhill running affects the late but not the early phase of the rate of force development. Eur J Appl Physiol 122(9):2049–2059. https://doi.org/10.1007/s00421-022-04990-8

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, Martinez A, Baggaley M, Khassetarash A, Giandolini M, Horvais N, Edwards WB, Millet GY (2020) 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

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Vernillo G, Martinez A, Baggaley M, Khassetarash A, Giandolini M, Horvais N, Edwards WB, Millet GY (2020) 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

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.