Does Characterizing Global Running Pattern Help to Prescribe Individualized Strength Training in Recreational Runners?

This study aimed to determine if concurrent endurance and strength training that matches the global running pattern would be more effective in increasing running economy (RE) than non-matched training. The global running pattern of 37 recreational runners was determined using the Volodalen® method as being aerial (AER) or terrestrial (TER). Strength training consisted of endurance running training and either plyometric (PLY) or dynamic weight training (DWT). Runners were randomly assigned to a matched (n = 18; DWT for TER, PLY for AER) or non-matched (n = 19; DWT for AER, PLY for TER) 8 weeks concurrent training program. RE, maximal oxygen uptake V̇O2max) and peak treadmill speed at V̇O2max (PTS) were measured before and after the training intervention. None of the tested performance related variables depicted a significant group effect or interaction effect between training and grouping (p ≥ 0.436). However, a significant increase in RE, V̇O2max, and PTS (p ≤ 0.003) was found after the training intervention. No difference in number of responders between matched and non-matched groups was observed for any of the performance related variables (p ≥ 0.248). In recreational runners, prescribing PLT or DWT according to the global running pattern of individuals, in addition to endurance training, did not lead to greater improvements in RE.

This study aimed to determine if concurrent endurance and strength training that matches the global running pattern would be more effective in increasing running economy (RE) than non-matched training. The global running pattern of 37 recreational runners was determined using the Volodalen R method as being aerial (AER) or terrestrial (TER). Strength training consisted of endurance running training and either plyometric (PLY) or dynamic weight training (DWT). Runners were randomly assigned to a matched (n = 18; DWT for TER, PLY for AER) or non-matched (n = 19; DWT for AER, PLY for TER) 8 weeks concurrent training program. RE, maximal oxygen uptakeVO 2 max) and peak treadmill speed atVO 2 max (PTS) were measured before and after the training intervention. None of the tested performance related variables depicted a significant group effect or interaction effect between training and grouping (p ≥ 0.436). However, a significant increase in RE,VO 2 max, and PTS (p ≤ 0.003) was found after the training intervention. No difference in number of responders between matched and nonmatched groups was observed for any of the performance related variables (p ≥ 0.248). In recreational runners, prescribing PLT or DWT according to the global running pattern of individuals, in addition to endurance training, did not lead to greater improvements in RE.

INTRODUCTION
Running economy (RE), which refers to the steady-state of oxygen consumption at a given running speed, is a critical factor of running performance (Conley and Krahenbuhl, 1980). RE improves after years of endurance running training, and especially if high volume, high intensity interval, or uphill running training are undertaken (Barnes and Kilding, 2014). Beyond running, different training strategies have been shown to potentially improve RE (Mikkola et al., 2007;Taipale et al., 2013;Barnes and Kilding, 2014). Among them, concurrent training, i.e., the use of strength training such as plyometric training (PLT) or dynamic weight training (DWT) in parallel with endurance running training, has been shown to further benefit RE (Barnes and Kilding, 2014). For instance, several studies that used this concurrent training method reported an improvement of RE ranging from 0 to 4.7% (Pellegrino et al., 2016;Meszler et al., 2019). However, the exact mechanisms leading to an improvement of RE after PLT or DWT remained unclear (Trowell et al., 2020).
PLT involves eccentric-concentric contraction cycles to allow the muscle-tendon unit to efficiently store and release elastic energy. During such cycle, there is a focus on a short ground contact time (t c ) and a high leg stiffness (k leg ) (Anderson, 1996). Common PLT exercises for runners are repeated rebound jumps or drop jumps. On the other hand, DWT involves a greater focus on concentric contractions aiming to produce a maximal power output, which is a compromise between speed and force generation (Kawamori and Haff, 2004). Common DWT exercises are, e.g., squats jumps and dynamic lunges. From a kinematic point of view, PLT implies shorter t c than DWT, de facto theoretically more in line with the mechanical demands of running. Indeed, running is characterized by a short contact phase (dependent on the running speed but generally smaller than 400 ms) followed by a flight phase. Therefore, the running pattern is a succession of plyometric contractions showing a spring like behavior, as suggested by the spring mass model (Blickhan, 1989). However, all runners do not share a running pattern that equally resembles to a spring. Indeed, some runners were shown to exhibit a more asymmetric contact phase (i.e., a longer duration of the propulsion phase than the braking phase) and less vertical oscillation of their center of mass during the flight time (t f ) than would be predicted by the spring-mass model . Thus, these running patterns are less accurately modeled by the spring-mass model. Following such ideas, it has been show that individuals could be classified into two categories termed aerial (AER) and terrestrial (TER) using the subjective Volodalen R scale (Gindre et al., 2015). Shorter t c and greater k leg are exhibited in AER than TER, while greater leg compression during stance is observed in TER compared to AER (Gindre et al., 2015). These kinematic differences might indicate that theoretically certain training modalities such as PLT or DWT might better suit AER or TER, respectively.
It is well established that individual differences exist in response to training, where high responders show large responses whereas low responders show small responses or no responses at all (Mann et al., 2014). Interestingly, Hautala et al. (2006) reported that low responders to an endurance based training program could become high responders to a strength based training program. Unfortunately, this variability in training responsiveness is not well understood and might be attributable to various factors including the absence of definition for high and low responders in the scientific literature and a one size fits all approach to exercise prescription (Mann et al., 2014). It has been purported that a more "personalized approach" to exercise prescription based on factors such as genotype, baseline phenotype, pre-training autonomic activity, individual homeostatic stress responses, recovery, and nutrition should improve training responsiveness (Mann et al., 2014). However, more research is still needed to clarify and quantify the role of these parameters. Also, for coaches, these factors are often hard to assess. In line with this view, Gindre et al. (2015) made the assumption that AER and TER could respond preferentially to different types of training interventions to improve RE. In other words, the knowledge of the global running pattern might provide useful indications for the prescription of training modalities toward an improvement in RE.
Hence, the purpose of the present study was to verify the effectiveness (i.e., mean increase in RE) and responsiveness (i.e., number of participants with a significant increase in RE) of two strength training modalities (i.e., PLT and DWT) on top of a standard endurance running training program to improve RE in runners with different global running patterns (i.e., AER and TER). We hypothesized that a training program that matches the underlying kinematics of the global running pattern (i.e., PLT for AER and DWT for TER) would be more efficient and thus would trigger a greater increase in RE and a lower rate of low-responders than a non-matched training program.

Participants
The study has been conducted over a 3 months period, which permitted to test 37 recreational and regular runners among which there were 5 females (age: 29.0 ± 9.0 years, height: 168 ± 6 cm, body mass: 59.3 ± 3.0 kg, weekly training hours: 2.0 ± 1.2 h) and 32 males (age: 29.4 ± 9.3 years, height: 177 ± 8 cm, body mass: 73.4 ± 12.4 kg, weekly training hours: 2.6 ± 1.3 h). For study inclusion, voluntary participants were required to be in good self-reported general health with no current or recent (<3 months) musculoskeletal injuries, and to have not previously undertaken any structured PLT or DWT. Two groups of runners were set up. The matched group consisted of AER following PLT and TER following DWT (n = 18). The non-matched group consisted of AER following DWT and TER following PLT (n = 19). As assessed by two-tailed non-matched t-tests, there were no significant differences in age, height, body mass, and weekly training hours between both groups ( Table 1). Participants were informed of the benefits and risks of the investigation prior to signing an institutionally approved informed consent document to participate in the study. They were informed that the data and results were confidential, and that they could withdraw at any time during the study, that was approved by an Institutional Review Board of the University of Bourgogne, Franche-Comté (CPP: 2014-A00336-41) and adhered to the latest Declaration of Helsinki of the World Medical Association (World Medical Association, 2013).

Experimental Approach to the Problem
After providing written informed consent, participants performed an initial baseline experimental session including a series of tests. These tests consisted of the assessment of the global running pattern using the Volodalen R scale to classify a runner as AER or TER, jump tests to evaluate the explosive concentric capacity and plyometric characteristics of the lower limb, a submaximal running test to determine RE, and a maximal incremental running test to determine peak treadmill speed (PTS) and maximal oxygen uptake (VO 2 max) (Gindre et al., 2015). Tests were interspersed by a 5-min passive recovery in a seated position. After that, each participant was randomly assigned to one of two 8 week concurrent training modalities, i.e., a standard endurance running training program combined with either PLT or DWT. After this assignment, participants were regrouped for statistical analysis based on whether their running pattern (AER or TER) was matched or non-matched with their prescribed strength training (DWT or PLT).

Global Running Pattern Assessment
During the warm-up of the initial baseline experimental session (5 min on a treadmill at 9 km·h −1 ), two running coaches with more than 3 years of experience using the Volodalen R method (CG and TL) paid attention to five key elements: vertical oscillation of the head, antero-posterior motion of the elbows, pelvis position at ground contact, foot position at ground contact, and foot strike pattern. Each element was scored from one to five, leading to a global subjective score (V R score) that represents the global running pattern of participants. This score ultimately allows the classification of runners into the two different categories (i.e., AER if V R score > 15 and TER otherwise). The Volodalen R method was fully described and studied elsewhere (Gindre et al., 2015) and was shown to be a reliable method to assess running pattern . The two coaches disagreed in their assessment of 3 individuals (8.1%). In these cases, the two coaches adopted a consensus following a discussion.

Endurance Running Training
All participants followed a basic endurance running training in line with what they were used to do before the study. Noteworthy, participants were not following a proper periodization training before starting the given training, i.e., they were not in a specific phase of a global periodization training. Training was divided into three different intensities based on their PTS: below 80%, between 80 and 95%, and between 95 and 105% of PTS. These percentages were chosen as to represent an aerobic, threshold and high intensity zone, respectively. The prescribed time in each of these three training zones during the 8 weeks training is described in Table 2. Main training volume (70-80%) was spent at running speeds below 80% of PTS.
Basic endurance sessions consisted of continuous running for 45-75 min, predominantly at a running speed below 80% of PTS with some unstructured bouts of faster running at 80-95% of PTS between 10 and 25 min per session. Interval sessions consisted of a 15 min easy warm-up at a running speed below 80% of PTS and involved repeated interval bouts ranging from 30 s to 2 min at 95-105% of PTS for an accumulated total of 6-12 min of fast running per session. In the beginning of the 8 weeks training plan, an example interval session consisted of 2 times 6 min of (30 s at 100% of PTS-30 s below 80% of PTS) with 2 min recovery between each 6 min block while at the end of the 8 weeks training plan, an example interval session consisted of 3 blocks of 2 repetitions of (2 min at 100% of PTS-1 min 30 s below 80% of PTS) with 5 min recovery between each block.

Plyometric or Dynamic Weight Training
Participants were asked to perform a predetermined circuit training composed of six exercises and designed as PLT or DWT (Figure 1). Details of the 8 weeks training are given in Table 2. Participants performed the same circuit training during the entire protocol but with progressive changes in the number of cycles and the exercise/rest ratio. Hence, despite different exercises, the total training load (total duration of effort and resting periods) was aimed to be equivalent between groups. Also, as the participants had no previous experience in resistance training, only body weight was used.

Jump Test
The squat jump test (SJ) was used to evaluate the explosive concentric capacity of the lower limbs (Bosco et al., 1983). Participants were required to jump vertically as high as possible from a static squat position and to start the landing with knees straight and ankles plantar-flexed. The depth of the squat was selfselected. Participants had to maintain the static squat position for two seconds prior to the jump. Squat jump height (SJ-h, in cm) was calculated from flight time (t f ) (Eq. 1) as measured by an optical measurement system (Optojump Next R , MicroGate Timing and Sport, Bolzano, Italy) sampling at 1,000 Hz. (1) Following the SJ, a five-repetition rebound jump test (5RJ) was used to evaluate the plyometric characteristics of the participants' lower limbs (Bosco et al., 1983;Dalleau et al., 2004). Participants were required to jump vertically as high as possible while minimizing ground contact time (t c ) and maximizing flight time (t f ). Participants were also instructed to minimize knee actions (i.e., flexion and extension) during the test. t f and t c  were measured by the Optojump Next R system. The average mechanical power during the positive (concentric) work per body mass (5RJ-P, in W·kg −1 ) was then calculated on the basis of the methods described by Bosco et al. (1983) using the following formula (Eq. 2).
All jumps were performed with hands placed on the hips and participants were wearing their habitual running shoes. After five practice trials of each jump, three repetitions of each jump test were performed with a 30 s rest between repetitions and a 2 min rest between the SJ and 5RJ tests. The best repetition of the SJ (based on the longest t f ) and 5RJ (based on the highest average mechanical power) was used for statistical analysis.

Submaximal Running Test
Participants ran for 5 min on a treadmill at 12 km·h −1 . Gas exchange was measured breath-by-breath using a gas analyser (Cortex Metamax 3B, Cortex Biophysik, Leipzig, Germany) and subsequently averaged over 10 s intervals throughout the test. Before each test, the gas analyzer was calibrated following the manufacturer's recommendations using ambient air (O 2 : 20.93% and CO 2 : 0.03%) and a gas mixture of known composition (O 2 : 15.00% and CO 2 : 5.00%). The spirometer was calibrated using a 3 L syringe. Respiratory exchange ratio (RER), oxygen uptake (VO 2 ), and carbon dioxide output (VO 2 ) were averaged over the last minute of the 5 min running trial. RER had to remain below 1.0 during the trials for the data to be included in the analysis, otherwise the corresponding data were excluded as deemed to not represent a submaximal effort. In such case, the selected submaximal speed was lowered iteratively by 1 km·h −1 until an RER below 1.0 was achieved. This resulted in submaximal testing speeds of 9 (n = 1), 10 (n = 6), 11 (n = 5), and 12 km·h −1 (n = 25). These speeds were kept the same for the post testing. RE was calculated from the running velocity divided by netVO 2 normalized to individual body mass (m·ml −1 ·kg −1 ) where neṫ VO 2 =VO 2 −restVO 2 , with restVO 2 given by the average over the last minute of a 5 min upright stance measure prior to the submaximal running test. This choice of units for RE have a conceptual advantage that numerical values are directly related to RE (i.e., the larger the numerical value, the better the RE) (Turner et al., 2003).

Maximal Incremental Test
Following the submaximal test, participants performed a maximal incremental running test on the treadmill. Starting at 8 km·h −1 , the treadmill speed was increased by 0.5 km·h −1 every minute until volitional exhaustion. The participants received strong verbal encouragement to ensure attainment of maximal values during the test.VO 2 max, averaged over 30 s, was said attained when two or more of the following criteria were met: an increase inVO 2 less than 2.1 ml·kg −1 ·min −1 between two consecutive stages, an RER greater than 1.1, and a heart rate (RS810, Polar Electro Oy, Kempele, Finland) of ±10 beats per minute of the predicted maximal heart rate value (i.e., 220 − age), as done by Howley et al. (1995). PTS is defined as the running speed of the last fully completed increment (MAS) plus the fraction of time spent in the following uncompleted increment (α) multiplied by the running speed increment ( s = 0.5 km·h −1 ) (Kuipers et al., 2003): PTS MAS + α s.

Statistical Analyses
Assuming a medium effect size [partial eta squared (η 2 p ) = 0.06] in RE improvement between matched and non-matched training groups, an α error of 0.05, and a power of 0.8, sample size calculations resulted in the requirement of 34 participants (Faul et al., 2007). However, the 37 participants were kept to slightly increase statistical power. Test-retest reliability coefficients ranged from r = 0.805 to 0.954 (p < 0.001) indicating good to excellent reliability of measurements. Descriptive statistics are presented using mean ± standard deviation. Effect sizes are reported as η 2 p values. The normality of the data and homogeneity of variances were verified using Shapiro-Wilk (p range: 0.163-0.943) and Levene's test (p range: 0.162-0.880), respectively. Unpaired two-sided Student's t-tests were used to compare participant characteristics between matched and nonmatched groups at baseline. Statistical analysis was performed using Jamovi [version 1.0.8 (Computer Software), retrieved from https://www.jamovi.org] with a level of significance set at p ≤ 0.05.

Effectiveness of the Training Interventions
A pre-post experimental design was used with two training groups (matched vs. non-matched). Effectiveness of the training protocol on performance parameters (primary criteria RE; secondary criteriaVO 2 max, PTS, SJ-h, and 5RJ-P) was assessed by repeated measures ANOVA (RM-ANOVA) with pre vs. post testing as within-subject factor and matched vs. non-matched grouping as between subject factor, and employing Bonferroni procedures for pair-wise post-hoc comparisons.

High Responders vs. Low Responders
Participants were all labeled as a responder or non-responder for the three performance variables that were significantly influenced (significant pre-post effect reported by the RM-ANOVA) by the protocol (i.e., RE,VO 2 max, and PTS) based on set % changes derived from the literature. A participant was determined as responder when RE increased by more than 2.6% ( Barnes and Kilding, 2015),VO 2 max by more than 5.9% (Dalleck et al., 2016), and PTS by more than 4% (arbitrary cut-off). Chi-squared analyses (χ 2 ) were performed on the number of responders and non-responders to assess if there was a difference of responsiveness for any of the three performance variables within matched and non-matched groups.

Effectiveness
No significant group effect or interaction effect were found between the training (pre-post) and grouping (matched vs. nonmatched) for any of the tested performance related variables. These results indicate that the effect of the applied training intervention on the performance related variables did not significantly differ between the matched and non-matched groups (p ≥ 0.436; Table 3). However, we found a significant increase in RE, PTS, andVO 2 max after the training intervention (p ≤ 0.003; Table 3). Noteworthy, post-hoc comparisons were not investigated as no interaction effect was reported.

Responsiveness to Training
No statistical difference in the responsiveness to training intervention were found between matched and non-matched groups for any of the performance related variables that were significantly influenced (significant pre-post effect reported by the RM-ANOVA; Table 3) by the protocol (p ≥ 0.248; Table 4). Individual responses are shown in Figure 2.

DISCUSSION
This study aimed at determining the effectiveness and responsiveness to two strength training modalities (i.e., PLT and DWT) combined with standard endurance training to improve RE in recreational runners. Identifying the global running pattern (i.e., TER or AER) and matching it to a PLT or DWT type of strength training prescription resulted in similar RE improvements and response to training than if no matching had been performed. As such, the results of the present study could not support our hypotheses. The following discussion is elaborating on possible explanations for the rejection of our hypothesis and directions for future research.
When following a certain training intervention, some individuals show a large positive response while others a small 3 | Mean ± SD for running economy (RE), maximal oxygen uptake (VO 2 max), peak treadmill speed (PTS), squat jump height (SJ-h), and average mechanical power during the positive (concentric) work per body mass of five repeated rebound jumps (5RJ-P) per training group, pre and post the training intervention as well as main effects (pre-post and group) and interaction effect (pre-post × group) for these five performance related variables.  or even no response (Mann et al., 2014). Moreover, individual responsiveness to training can vary by training mode (Hautala et al., 2006). Likewise, in the present study we aimed at improving RE and found a significant increase of 2.3% at groups level, but with individual effect ranging from a 4.8% decrease to a 13.8% increase. This was concomitant with an average PTS increase of 4.3% (ranging from a 3.8% decrease to a 12.5% increase) and aVO 2 max increase of 3.8% (ranging from a 10.6% decrease to a 14.8% increase). This wide range of individual responses highlights the importance of taking into account individual responses and not only the training effect on a group level (see Figure 2). It has been reported that average group level increases in RE ranged from 0 to 4.7% using heavy weight strength training (Johnston et al., 1997) and 0-4% using explosive training (Pellegrino et al., 2016;Meszler et al., 2019), which is in line with the current findings. Nevertheless, several studies reported no significant increase in RE following concurrent endurance and strength training in recreational athletes despite improvements in muscle strength (Ferrauti et al., 2010;Mikkola et al., 2011;Taipale et al., 2013;Damasceno et al., 2015). A possible explanation for the relatively small improvements in RE in our study can be found in the duration of the training period. A review by Denadai et al. (2017) highlighted that longer training periods (>8-21 weeks) are likely to result in greater RE improvements than shorter programs (6-8 weeks) due to the time-course of neuromuscular adaptations following concurrent training. A review on this topic by Rønnestad and Mujika (2014) stated that the possible mechanisms of how concurrent endurance and strength training can improve RE are related to a delayed activation of less efficient type II muscle fibers, an improved neuromuscular efficiency, the conversion of fast-twitch type IIx fibers into more fatigue resisted type IIa fibers, and an improved musculo-tendinous stiffness. During a short to medium term training period (up to 8 weeks), as was the case in the current study, the expected neuromuscular adaptations are an increased neural activation and a smaller relative proportion of type IIx than type IIa fibers, while an optimized musculo-tendinous stiffness is only achieved after longer training periods (>8-21 weeks). The absence of significant effects on any of the jumping performance parameters (SJ-h and 5RJ-P) in this study confirms a possible lack of musculo-tendinous adaptations after the 8 weeks training period. Therefore, when aiming at improving RE by strength training, longer training periods (>8-21 weeks) are advised.

Pre
Another possible explanation for unachieved RE improvement through concurrent training in some studies is that because the exact mechanisms behind RE improvement are still unknown, or the appropriate stimulus is not used. Some studies were able to induce neuromuscular adaptations (i.e., increase in maximal and explosive strength) through heavy weight training, while the intended effects on RE remained absent (Taipale et al., 2013;Vikmoen et al., 2016). A recent review on the effect of strength training on biomechanical and neuromuscular adaptations concluded that evidences that neuromuscular effects obtained by strength training transfer to running biomechanics are lacking (Trowell et al., 2020). In this study we tried to match the concurrent strength training to the global running pattern in order to obtain individualized neuromuscular adaptation in an attempt to maximize RE improvements. However, as we did not obtain better results, it seems that more insights are needed regarding the interaction between neuromuscular stimuli, its resulting adaptations, their transfer to the running pattern and their impact on RE.
A possible limitation to the current study might be that executed training sessions were unsupervised. Indeed, after an initial supervised strength training session, athletes were instructed to perform the training sessions on their own. Moreover, as runners were novice to strength training, the load of the training sessions was kept submaximal to avoid injuries. Higher training load might lead to greater RE improvements. Therefore, the adherence, intensity, order, and organization of the training sessions was not strictly controlled and could partly explain the low mean training responses as well as the large interindividual differences. On the other hand, it represents real-life conditions. In addition, even though runners were used to do interval trainings before starting the endurance training, the fact that they were now following a structured endurance training instead of their own "unstructured" one might partly explain the increase ofVO 2 max, PTS, and RE. Finally, RE was assessed at a fixed running speed (12 km·h −1 ), which was obviously not individualized for each participant. An alternative could have been to determine RE at different running speeds, as long as they fall below the respiratory compensation point and that a steady-state of oxygen consumption was reached within 3-15 min (Barnes and Kilding, 2015). These speeds could have been chosen to correspond to theoretical optimal running speeds to run 5, 10, 21, and 42 km races for each individual or to participants personal best on these distances. However, obtaining these speeds would have required to perform the maximal incremental running test before the submaximal one and to perform several submaximal tests, which would have increase the duration of the overall testing.
Future research should continue to focus on how neuromuscular adaptations, induced by individualized concurrent strength and endurance training, relate to changes in running biomechanics and improvement in performance (Trowell et al., 2020). Specific attention should go to the time-course of these adaptations to reveal if improvements are mainly made during the initial training period or if more long-term progress can be made, depending on the intervention. As well, a specific evaluation of the individually different responses to strength training, including PLY and DWT, should be done. Understanding these mechanisms should help predetermining which runner (i.e., global running pattern) needs which additional strength training to optimize performance. Such an individualized approach remains the ultimate goal for coaches and athletes.
As practical guidelines we can conclude that inter individual differences in training response to concurrent training are substantial (Figure 2). We encourage coaches and athletes to regularly evaluate the effectiveness of the prescribed training program and to keep looking for ways to individualize and optimize training responses. An initial assessment of the global running pattern as being rather AER or TER using the Volodalen method can be used as a way to identify the runners' preferences and be a guideline for training individualization.

CONCLUSION
As a conclusion, prescribing PLT or DWT strength training based on global running pattern, in addition to regular endurance training, did not lead to greater improvements in RE for recreational runners. In order to be able to optimize strength training prescription and its individualization in endurance runners, future research should aim to understand the exact mechanisms relating strength training to the resulting neuromuscular and biomechanical adaptations while running.

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

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Institutional Review Board of the University of Bourgogne, Franche-Comté (CPP: 2014-A00336-41). The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
CG, LM, and TL: conceptualization, methodology, and supervision. AT and TL: investigation. AP, BB, and TL: formal analysis and writing-original draft preparation. AP, BB, AT, LM, CG, and TL: writing-review and editing. All authors contributed to the article and approved the submitted version.

FUNDING
This study was supported by the University of Franche-Comté, Tomsk Polytechnic University CE Program, and Volodalen Swiss Sport Lab (CIFRE thesis 2013/0425).