# Subject specific muscle synergies and mechanical output during cycling with arms or legs

Théo Cartier, Laurent Vigouroux, Elke Viehweger, Guillaume Rao

► To cite this version:

Théo Cartier, Laurent Vigouroux, Elke Viehweger, Guillaume Rao. Subject specific muscle synergies and mechanical output during cycling with arms or legs. PeerJ, 2022, 10.7717/peerj.13155. hal-03644422

HAL Id: hal-03644422

https://amu.hal.science/hal-03644422

Submitted on 19 Apr 2022

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
# Subject specific muscle synergies and mechanical output during cycling with arms or legs

Théo Cartier<fnref n="1" />, Laurent Vigouroux<fnref n="1" />, Elke Viehweger<fnref n="2" />,<fnref n="3" /> and Guillaume Rao<fnref n="1" />

## ABSTRACT

**Background.** Upper (UL) and lower limb (LL) cycling is extensively used for several applications, especially for rehabilitation for which neuromuscular interactions between UL and LL have been shown. Nevertheless, the knowledge on the muscular coordination modality for UL is poorly investigated and it is still not known whether those mechanisms are similar or different to those of LL. The aim of this study was thus to put in evidence common coordination mechanism between UL and LL during cycling by investigating the mechanical output and the underlying muscle coordination using synergy analysis.

**Methods.** Twenty-five revolutions were analyzed for six non-experts' participants during sub-maximal cycling with UL or LL. Crank torque and muscle activity of eleven muscles UL or LL were recorded. Muscle synergies were extracted using nonnegative matrix factorization (NNMF) and group- and subject-specific analysis were conducted.

**Results.** Four synergies were extracted for both UL and LL. UL muscle coordination was organized around several mechanical functions (pushing, downing, and pulling) with a proportion of propulsive torque almost 80% of the total revolution while LL muscle coordination was organized around a main function (pushing) during the first half of the cycling revolution. LL muscle coordination was robust between participants while UL presented higher interindividual variability.

**Discussion.** We showed that a same principle of muscle coordination exists for UL during cycling but with more complex mechanical implications. This study also brings further results suggesting each individual has unique muscle signature.

**Subjects** Anatomy and Physiology, Neurology, Orthopedics, Biomechanics

**Keywords** Arm cranking, Upper limb, Lower limb, Cycling, Muscle synergies, Muscle coordination

## INTRODUCTION

Cycling with the upper (UL) or lower limbs (LL) both consists in producing torque in a circular movement with fixed cranks length around a fixed rotation axis. Cycling is extensively used for improving health, training, physiological evaluation, or rehabilitation with LL (Faria, Parker & Faria, 2005; Fujiwara PhD et al., 2005; Oja et al., 2011; Sitko et al., 2020) and to a lesser extent with UL (Hübner-Wozniak et al., 2004;

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

<footnote n="1">¹ Aix Marseille Univ, CNRS, ISM, Marseille, France</footnote>

<footnote n="2">² Department of Biomedical Engineering, University of Basel, Allschwil, Switzerland</footnote>

<footnote n="3">³ Department of Orthopedics and Gait Laboratory, University Children's Hospital of Both Basel, Basel, Switzerland</footnote>
Elmer, Danvind & Holmberg, 2013; Zinner et al., 2016; Zhou et al., 2017). An increasing number of studies highlight the interest of UL cycling also for LL rehabilitation: for example, Kaupp et al. (2018) showed functional and neurological improvements of walking capacities after a rehabilitation protocol involving UL cycling. Moreover, modulations of cutaneous reflex and muscle activities of the LL was also shown after UL cycling training (Zehr & Chua, 2000; Zehr, 2005; Barzi & Zehr, 2008; Mezzarane, Nakajima & Zehr, 2014). These examples suggest similarities in movement execution and interactions in the production of movement with the UL and LL that can be observed at muscle activity level. This latter point is particularly interesting from a clinical point of view since intervention at the UL could benefit the LL functional aspects.

To deeply understand the movement of cycling, it remains crucial to understand the underlying muscle coordination. Muscle coordination leads to the question of understanding how the large number of available muscles are coordinated to execute a movement (Bernstein, 1967; Latash, 2012). The muscle synergy analysis is frequently used to answer this question (d'Avella, Saltiel & Bizzi, 2003; Turpin, Uriac & Dalleau, 2021). This analysis assumes that the Central Nervous System (CNS) simplifies the abundance of muscles by using smaller numbers of functional elements activated in time. Muscles synergies thus allow a global understanding of the muscle coordination in comparison to the production of a mechanical output. The existence of muscle synergies on the LL during cycling has been widely demonstrated even when manipulating various mechanical constraints such as posture or output torque intensity (Hug et al., 2011). During LL cycling, it has been shown that muscle synergies are shared across participant despite interindividual variability of EMG pattern (Hug et al., 2010; De Marchis et al., 2013). Further studies showed that this movement is executed, from a kinematic and kinetics point of view, in the same way by any healthy person, expert or not (Mornieux et al., 2008).

In comparison to LL cycling, UL cycling muscular coordination remains poorly understood despite the interests for training and clinical applications. A single recent study aimed to extract muscle synergies during UL cycling (Botzheim et al., 2021). Using a set of four muscles for each arm, these authors showed that bimanual muscle coordination can be described by muscle synergies. However, no simultaneous analysis of the mechanical output was performed and only a few numbers of representative muscles were used to fully understand the functional synergies during UL cycling. More generally, it remains unknown whether the UL synergies are similarly organized than in LL with respect to the mechanical output. Unfortunately, this lack of knowledge limits the potential for these exercises to be objectively used for training and clinical applications.

Complementary, several studies highlight the necessity to take into account individual differences (Bartlett, Wheat & Robins, 2007) and a recent work of Hug et al. (2019) put in evidence that individuals have unique muscle signatures. Each individual muscle signature could result from numerous individual factors such as development, experiences, or diseases discrepancies. Contrary to LL, UL musculoskeletal system offers the possibility to pull or push equally, allowing several possible strategies in movement execution. Individual differences in synergies can thus be differently managed in the UL compared to LL.
The first objective of this study was thus to fill the knowledge gap by investigating the simultaneous temporal evolutions of the mechanical output and of the EMG and muscle synergies of the Upper Limb in comparison to the Lower Limb during cycling. To individualize our analysis, our second objective was to compare interindividual variability between UL and LL at kinematics, kinetics, and muscle coordination levels.

# MATERIALS & METHODS

## Participants

Six young healthy adults (age: 27.2 ± 2.9 years; weight: 69 ± 7.5 kg; height: 174.2 ± 5.1 cm) volunteered to participate in the study and suffered of no injuries in the past 6 months. To limit the effect of participants' level, expertise in either UL or LL cyclic exercise constituted an exclusion criterion. Participants gave their written consent and were informed of the possible risks and discomfort associated with the experiment. The experimental protocol was approved by the national ethic committee (CERSTAPS, IRB00012476-2020-30-11-75) and the study was conducted in accordance with the Helsinki declaration.

## Experimental procedure

The experimental session consisted in two phases performed successively with the UL or the LL in random order. Before each phase, participants performed a 5-minutes standardized warm up with the tested limb at self-selected cadence (SSC: 85.0 ± 10.4 vs 86.3 ± 9.8 rpm, *p* > 0.05, for arms and legs respectively). The self-selected cadence was then used and controlled for the entire protocol. To estimate maximal power (MP), participants performed a single 10 s isokinetic all-out with strong verbal encouragements. MP was computed as the mean peak power during the 10 s. The same relative amount of power was set at 30% of MP to reach a significant intensity without any excessive fatigue and to compare both limbs at the same relative mechanical demands. Participants then performed a submaximal 2-minute test at 30% of MP at self-selected cadence (102.8 ± 14.7 vs 171.7 ± 34.4 W, *p* < 0.05, for arms and legs respectively). The two experimental phases were separated by at least 35 min rests to limit the fatigue effect and change the experimental equipment (EMG, markers etc.) from one limb to another.

## Material and data collection

### Ergometer

The task was performed on a cycle ergometer (Excalibur for LL and Brachumera for UL, Lode, Groningen, The Netherlands) equipped with standard cranks, flat pedals for feet and standard horizontal pegs for hands. The type of pedals interface was chosen to study the most ecological situation for each limb. For LL cycling, participants were asked to seat comfortably on the saddle and maintain their position during the entire experimental procedure. For UL cycling, participants were seated in front of the ergometer in pronated grip (Fig. 1). The crank axis height was positioned at the glenohumeral level. Participants were preferentially seated on a chair and were instructed not to grab the chair legs with their own legs. Torque and pedalling frequency were recorded each 2° by the Lode Ergometry Manager software. Ergometers were calibrated before each evaluation.
<figure><img src="image_3.png" /><figcaption>**Figure 1** **Experimental set-up.** Left column: Example of experimental set-up for UL (top) and LL (bottom) and schematic representation of the cycle revolution for each. The yellow dot represent the SH to FH transition and the blue dot the FH to SH transition. FH, first half, SH, second half.</figcaption></figure>

Full-size

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

DOI: 10.7717/peerj.13155/fig-1

## Kinematics

Kinematics was recorded using 5 Qualisys cameras (Oqus 7) at 200 Hz with Qualisys Track Manager software (Qualisys, Sweden). Trajectories of five reflective markers placed on distal extremity of the 5th *metacarpal*, *ulnar styloid*, *lateral epicondyle*, *acromion* and *iliac crest* and distal extremity of the 5th *metatarsal*, *lateral malleolus*, *lateral knee epicondyle*, *greater trochanter* and *iliac crest* were used to compute the 2D angular trajectories of the ankle, knee, hip and wrist, elbow, shoulder joints.

## EMG

The electromyography of 22 muscles were recorded using EMG system (Wireless System, DELSYS, USA) sampled at 2000 Hz with the Delsys Trigno Acquisition software. Before electrode application, the skin was shaved and cleaned with an alcoholic solution to improve the electrode/skin impedance. SENIAM recommendations were followed to ensure correct electrode placement and avoid crosstalk between EMG signals. A functional test was carried out for each muscle to ensure correct electrode placement. The following LL muscles were recorded: *Tibialis Anterior* (TA), *Soleus* (Sol), *Gastrocnemius Medialis/Lateralis* (GM, GL), *Vastus Medialis/Lateralis* (VM, VL), *Biceps Femoris* (BF), *Semimembranosus* (SM),
Rectus Femoris (RF), Tensor Fasciae Latae (TF) and Gluteus Maximum (GMAX) (Hug & Dorel, 2009; Hug et al., 2010; Hug et al., 2011; De Marchis et al., 2013). The following UL muscles were recorded: Extensor Carpi Radialis Longus (ECRL), Flexor Carpi Radialis (FCR), BrachioRadialis (BR), Biceps Brachii (BB), Triceps Brachis Longus/Shortus (TBL, TBS), Deltoid Anterior/Posterior (DA, DP); (Zehr & Chua, 2000; Barzi & Zehr, 2008; Zehr, Loadman & Hundza, 2012; Chaytor et al., 2020), Pectoralis Major (PM), Latissimus Dorsi (LD), Flexors Digitorum Superficialis (FDS). EMG, torque, and kinematic recording were synchronized via a trigger impulse.

## Data processing

### Torque and kinematics

Marker trajectories and torque data were low pass filtered (4th order, zero-time lag Butterworth 6 Hz). Absolute 2D angle of the LL (Hip, Knee, Ankle) and UL (Shoulder, Elbow Wrist) were computed using vector coordinates of each segment in a sagittal plane of the ergometer. Pedalling revolutions were extracted using the position of the distal marker. For LL, 0° was defined as the higher position of the pedal and 180° as the lower position. For UL, 0° was defined as the nearest position of the pedal, when the crank is horizontally directed towards the participant and far 180° the farthest when the crank is horizontally directed to the opposite of the participant. Two phases are considered in the movement, the First Half (FH) from 0° to 180° and the Second Half (SH) from 180° to 360°. The FH to SH transition was thus when the limb is fully extended and the SH to FH transition when the limb is grouped. To normalize torque between participants, absolute values were divided by body mass.

### EMG pre-processing

EMG data were first filtered using a bandpass 4th order Butterworth zero-time lag filter (20–400 Hz) then full wave rectified and low pass filtered (4th order Butterworth, zero-time lag 4 Hz) to obtain EMG envelops. For amplitude analysis, the EMG data were normalized by the maximal value recorded during the 10 s all-out. A dedicated detection algorithm has been used to select the 25 most representative cycles over the entire dataset (Sangeux & Polak, 2015).

After revolution extraction and pre-processing, torque, kinematics, and EMG data were resampled to 200 points by revolution.

### Synergy extraction

Since the aim of our study was to quantify muscle coordination, EMG data were normalized by the maximal value across the cycles before synergy extraction and synchronous muscle synergy model was used to extracted synergies Eq. (1):

$$V_O(t) = V_R(t) + \varepsilon(t) \quad \text{with} \quad V_R(t) = \sum_{s=1}^{K} H_s(t) W_s \qquad (1)$$

where $V_O(t)$ indicates the initial EMG values of all muscles at time instant t, $V_R(t)$ is the reconstructed signal resulting from the linear combination of $W_s$ (the muscle synergy vector for the s-th synergy) and $H_s(t)$ (the value of the synergy activation coefficient for
<figure><img src="image_5.png" /><figcaption>**Figure 2** Reconstructed (red) vs original (blue) EMG pattern. Example of reconstructed vs original EMG pattern for representative participant for lower limb muscles (left) and upper limb muscles (right).</figcaption></figure>

Full-size DOI: 10.7717/peerj.13155/fig-2

the s-th synergy at time instant t). K is the number of synergies necessary to reconstruct the initial EMG dataset, and ε(t) represents noise. This model is a compromise that allows the analysis of both spatial (through W) and temporal (through H) structures of the motor outputs (Turpin, Uriac & Dalleau, 2021).

Nonnegative Matrix Factorization (NNMF) algorithm was applied to a matrix containing the envelop of the 11 muscles across 25 cycles. The algorithm approximates the initial signal $V_O$ with a reconstructed matrix $V_R = W \times H$ by minimizing the Frobenius norm $||V_O - V_R||$ where W is a K × 11 matrix of the synergy vectors and H is a K × 5000 (25 cycles * 200 time points). The algorithm was run for K in range from 1 to 10 and 10 times for each K synergies combinations to avoid local minima and keep the best solution for further analysis.

The number of synergy sufficient to reconstruct the original matrix $V_O$ was set by the calculation of the Variance Accounted For (VAF) by $V_R$ on the entire matrix (VAF$_{tot}$) and for each muscle individually (VAF$_m$) with following equations:

$$VAF_{tot} = 1 - \frac{\sum (VO_{tot} - VR_{tot})^2}{\sum VO_{tot}^2} \qquad (2)$$

$$VAF_m = 1 - \frac{\sum (VO_m - VR_m)^2}{\sum VO_m^2}. \qquad (3)$$

The minimum number of synergies that explain at least 90% of VAF$_{tot}$ and 75% of VAF$_m$ was chosen. An example of reconstructed vs original EMG pattern was shown on Fig. 2. When both VAF constraints are satisfied (i.e., total VAF and by muscle VAF), 4 synergies were extracted for each participant to allow an analysis and a functional comparison between participants.
## Statistical analysis

For each temporal data (Angle, Torque, EMG and Hsyn), the inter- or intra-individual variability was calculated using a Pearson correlation coefficient r (r<sub>inter</sub> or r<sub>intra</sub>) and a Variance Ratio parameter (VR<sub>inter</sub> or VR<sub>intra</sub>).

For interindividual variability, the input matrix X was composed of a mean curve across the 25 cycles for each participant (final size = 6 × 200). For intra-individual variability, the input matrix X was composed of the 25 cycles for each participant, (final size = 25 × 200).

The VR was computed as follow:

$$VR_a = \frac{\sum_{i=1}^{k} \sum_{j=1}^{n} (X_{ij} - \bar{X}_i)^2 / k(n-1)}{\sum_{i=1}^{k} \sum_{j=1}^{n} (X_{ij} - \bar{X})^2 / (kn-1)} \quad \text{with} \quad \bar{X} = \frac{1}{k} \sum_{i=1}^{k} \bar{X}_i. \quad (4)$$

where a stand for inter or intra, k is the number of samples (k = 200), n is the number of participants/cycles (i.e., 6/25), X<sub>ij</sub> is the value at the ith sample for the jth participant/cycle, X̄<sub>i</sub> is the mean value at the ith sample over the n participants/cycles.

The Variance Ratio parameter has been used in other studies (Burden, Trew & Baltzopoulos, 2003; Hug et al., 2010; De Marchis et al., 2013) as an indicator of the overall variation of the data with respect to a mean profile: the higher the VR, the higher the inter or intra participants variability in the dataset. The Pearson correlation coefficient r was computed for inter individual variability between each pair of participants excluding autocorrelation and a mean ± sd value is presented (Tables 1–2). For intra-individual variability, r was computed between each pair of cycle for each participant excluding autocorrelation and a mean ± sd value across all participants is presented (Tables 3–4). The Pearson correlation coefficient r indicates the similarity in the shape profile of the data without being impacted by amplitude differences. We also used this coefficient r to compute similarities across participant in muscle synergy vectors (W<sub>s</sub>).

Because VR and r do not take account a possible time shift between two identical curves, the normalized cross-correlation function was computed and allowed to extract the percentage of shift (t<sub>max</sub>) needed to maximize the cross-correlation function (r<sub>max</sub>). r<sub>max</sub> and t<sub>max</sub> were only computed for interindividual variability, intra-individual variability (through VR and r) being low for each participant. Even though these four variables (VR, r, t<sub>max</sub>, and r<sub>max</sub>) are redundant, they allow to nuance the interpretation of the variability according to differences in time, amplitude, or patterns. Furthermore, these different data allow an easier comparison of our results with the literature.

We performed correlations between the variables presenting high t<sub>max</sub> (H<sub>s</sub>, Torque) and the RPM to further analyze the relationships between the Realization variables and the Performance variables (see Discussion section for the definition). The time-shift of each synergy activation coefficient (H<sub>s</sub>) was computed by extracting the normalized time instants of max(H<sub>s</sub>) and Pearson r correlations were performed for all participants between each time-shift and pedaling cadence.
<table><caption>**Table 1** Interindividual variability of upper limb variables.</caption><thead><tr><th></th><th>VR</th><th>r<br>(mean ± sd)</th><th>r<sub>max</sub><br>(mean ± sd)</th><th>t<sub>max</sub> %<br>(mean ± sd)</th></tr></thead><tbody><tr><td colspan="5"><strong>KINEMATICS & DYNAMICS</strong></td></tr><tr><td>Torque</td><td>0.26</td><td>0.72 ± 0.12</td><td>0.84 ± 0.10</td><td>2.37 ± 2.46</td></tr><tr><td>Shoulder</td><td>0.12</td><td>0.93 ± 0.03</td><td>0.99 ± 0.01</td><td>0.87 ± 1.01</td></tr><tr><td>Elbow</td><td>0.27</td><td>0.84 ± 0.12</td><td>0.97 ± 0.02</td><td>0.00 ± 0.00</td></tr><tr><td>Wrist</td><td>0.79</td><td>0.25 ± 0.26</td><td>0.81 ± 0.15</td><td>11.77 ± 9.46</td></tr><tr><td colspan="5"><strong>EMG</strong></td></tr><tr><td>FDS</td><td>0.56</td><td>0.48 ± 0.42</td><td>0.81 ± 0.12</td><td>15.40 ± 12 ± 69</td></tr><tr><td>FCR</td><td>0.45</td><td>0.57 ± 0.40</td><td>0.77 ± 0.13</td><td>8.93 ± 8.93</td></tr><tr><td>BR</td><td>0.73</td><td>0.28 ± 0.56</td><td>0.88 ± 0.08</td><td>16.00 ± 11.41</td></tr><tr><td>ECRL</td><td>0.49</td><td>0.57 ± 0.33</td><td>0.91 ± 0.04</td><td>11.06 ± 7.05</td></tr><tr><td>DA</td><td>0.67</td><td>0.36 ± 0.49</td><td>0.86 ± 0.10</td><td>13.47 ± 10.09</td></tr><tr><td>DP</td><td>1.00</td><td>0.10 ± 0.50</td><td>0.71 ± 0.17</td><td>15.70 ± 10.37</td></tr><tr><td>LD</td><td>0.67</td><td>0.32 ± 0.49</td><td>0.87 ± 0.07</td><td>12.50 ± 8.12</td></tr><tr><td>PM</td><td>0.91</td><td>0.22 ± 0.40</td><td>0.73 ± 0.10</td><td>19.93 ± 14.14</td></tr><tr><td>TBL</td><td>0.46</td><td>0.53 ± 0.44</td><td>0.80 ± 0.16</td><td>23.47 ± 31.18</td></tr><tr><td>BB</td><td>0.61</td><td>0.38 ± 0.46</td><td>0.88 ± 0.07</td><td>13.87 ± 9.15</td></tr><tr><td>TBS</td><td>0.89</td><td>0.16 ± 0.52</td><td>0.65 ± 0.12</td><td>15.37 ± 14.06</td></tr><tr><td colspan="5"><strong>Synergies</strong></td></tr><tr><td></td><th>W<sub>s</sub></th><th></th><th>H<sub>s</sub></th><th></th></tr><tr><td>S<sub>0</sub></td><td>0.75 ± 0.17</td><td>0.61</td><td>0.40 ± 0.48</td><td>0.80 ± 0.15</td></tr><tr><td>S<sub>1</sub></td><td>0.71 ± 0.11</td><td>0.41</td><td>0.46 ± 0.42</td><td>0.90 ± 0.06</td></tr><tr><td>S<sub>2</sub></td><td>0.75 ± 0.12</td><td>0.70</td><td>0.30 ± 0.52</td><td>0.89 ± 0.08</td></tr><tr><td>S<sub>3</sub></td><td>0.75 ± 0.17</td><td>0.64</td><td>0.38 ± 0.46</td><td>0.88 ± 0.08</td></tr></tbody></table>

**Notes.**

Variance Ratio (VR), r-Pearson (r), maximum value of cross-correlation function (rmax) and time lag (tmax) of Kinematics, torque, EMG patterns, synergies vectors (Ws) and synergy activation coefficients (Hs).

# RESULTS

## Torque

Figure 3 shows the averaged torque applied on the right pedal for UL (top) and LL (bottom). Maximal peak torque was higher for LL than UL (0.72 ± 0.12 Nm/kg vs 0.18 ± 0.30 Nm/kg). For UL, the output torque was positive during most of the cycle (85 ± 9%) excepted at the end of SH and maximal torque was produced at the beginning of SH (around 199 ± 42°). For LL, the torque was positive during FH (52 ± 2% of the total cycle) and negative during SH. Maximal torque was produced about mid-FH (87 ± 8°).

## Joint kinematics

Figure 4 shows the kinematics of proximal, intermediate, and distal joints for LL (left) and UL (right). Hip and knee angular displacements result from extension during FH and flexion during SH. Ankle angular displacement result from plantar flexion in FH then dorsiflexion during SH. Shoulder angular displacement results from antepulsion until mid of FH then retropulsion until almost the end of crank revolution. Elbow angular
<table><caption>**Table 2** Interindividual variability of lower limbs variables.</caption><thead><tr><th></th><th>VR</th><th>r<br>(mean ± sd)</th><th>r<sub>max</sub><br>(mean ± sd)</th><th>t<sub>max</sub> %<br>(mean ± sd)</th></tr></thead><tbody><tr><td colspan="5"><strong>KINEMATICS & DYNAMICS</strong></td></tr><tr><td>Torque</td><td>0.04</td><td>0.97 ± 0.01</td><td>0.99 ± 0.01</td><td>0.81 ± 0.78</td></tr><tr><td>Hip</td><td>0.15</td><td>0.98 ± 0.01</td><td>0.99 ± 0.00</td><td>0.00 ± 0.00</td></tr><tr><td>Knee</td><td>0.05</td><td>0.99 ± 0.00</td><td>1.00 ± 0.00</td><td>0.00 ± 0.00</td></tr><tr><td>Ankle</td><td>0.50</td><td>0.76 ± 0.13</td><td>0.89 ± 0.08</td><td>2.37 ± 2.77</td></tr><tr><td colspan="5"><strong>EMG</strong></td></tr><tr><td>SOL</td><td>0.20</td><td>0.84 ± 0.15</td><td>0.93 ± 0.04</td><td>3.40 ± 3.72</td></tr><tr><td>GL</td><td>0.37</td><td>0.68 ± 0.26</td><td>0.91 ± 0.05</td><td>9.87 ± 6.10</td></tr><tr><td>GM</td><td>0.19</td><td>0.82 ± 0.17</td><td>0.97 ± 0.03</td><td>5.73 ± 3.50</td></tr><tr><td>TA</td><td>0.57</td><td>0.37 ± 0.45</td><td>0.84 ± 0.10</td><td>24.93 ± 27.13</td></tr><tr><td>GMAX</td><td>0.07</td><td>0.97 ± 0.02</td><td>0.97 ± 0.02</td><td>0.00 ± 0.00</td></tr><tr><td>RF</td><td>0.35</td><td>0.80 ± 0.12</td><td>0.80 ± 0.12</td><td>0.43 ± 1.05</td></tr><tr><td>VL</td><td>0.05</td><td>0.96 ± 0.03</td><td>0.96 ± 0.03</td><td>0.10 ± 0.39</td></tr><tr><td>VM</td><td>0.08</td><td>0.94 ± 0.05</td><td>0.94 ± 0.04</td><td>0.33 ± 0.84</td></tr><tr><td>BF</td><td>0.45</td><td>0.56 ± 0.32</td><td>0.77 ± 0.14</td><td>8.90 ± 6.84</td></tr><tr><td>SEM</td><td>0.45</td><td>0.64 ± 0.28</td><td>0.88 ± 0.06</td><td>8.70 ± 6.10</td></tr><tr><td>TF</td><td>0.75</td><td>0.25 ± 0.58</td><td>0.68 ± 0.20</td><td>25.53 ± 19.91</td></tr><tr><td colspan="5"><strong>Synergies</strong></td></tr><tr><td></td><th>W<sub>s</sub></th><th></th><th>H<sub>s</sub></th><th></th></tr><tr><td>S<sub>1</sub></td><td>0.75 ± 0.13</td><td>0.06</td><td>0.95 ± 0.04</td><td>0.16 ± 0.52</td></tr><tr><td>S<sub>2</sub></td><td>0.73 ± 0.15</td><td>0.21</td><td>0.73 ± 0.15</td><td>5.73 ± 3.53</td></tr><tr><td>S<sub>3</sub></td><td>0.39 ± 0.37</td><td>0.48</td><td>0.54 ± 0.36</td><td>9.07 ± 5.74</td></tr><tr><td>S<sub>4</sub></td><td>0.87 ± 0.11</td><td>0.30</td><td>0.75 ± 0.24</td><td>6.27 ± 4.90</td></tr></tbody></table>

**Notes.**

Variance Ratio (VR), r-Pearson (r), maximum value of cross-correlation function (rmax) and time lag (tmax) of Kinematics, torque, EMG patterns, synergies vectors (Ws) and synergy activation coefficients (Hs).

displacement results from extension during FH then flexion during SH. Wrist angular displacement results from extension during FH then flexion during SH.

## Synergies description

### UL

Four synergies were extracted for each participant (Fig. 5-UL) and explained 96 ± 0.01% of the total variance of the initial EMG dataset (VAF<sub>tot</sub>) and 95 ± 0.05% of each initial individual EMG dataset (VAF<sub>m</sub>). The **first synergy** (S<sub>1</sub>) included the TBS and TBL during the transition from SH to FH and the first part of FH. The **second synergy** (S<sub>2</sub>) was composed by the activations of the hand flexors (FDS and FCR) and LD during end of FH and the transition from FH to SH (90–180°). The **third synergy** (S<sub>3</sub>) involved hand flexors (FDS/FCR) and extensors (ECRL), elbow flexors (BB and BR), elbow extensors (TBS) and shoulder extensor (DP) during SH (180–270°). The **last synergy** (S<sub>4</sub>) involved DA and PM in anticipation to the SH to FH transition (270–0°).
<table><caption>**Table 3** Upper limbs intra individual variability.</caption><thead><tr><th></th><th>VR</th><th>r (mean ± sd)</th></tr></thead><tbody><tr><td><strong>KINEMATICS AND DYNAMICS</strong></td><td></td><td></td></tr><tr><td>Torque</td><td>0.19</td><td>0.89 ± 0.13</td></tr><tr><td>Shoulder</td><td>0.06</td><td>0.98 ± 0.03</td></tr><tr><td>Elbow</td><td>0.15</td><td>0.94 ± 0.13</td></tr><tr><td>Wrist</td><td>0.58</td><td>0.70 ± 0.33</td></tr><tr><td><strong>EMG</strong></td><td></td><td></td></tr><tr><td>FDS</td><td>0.21</td><td>0.90 ± 0.09</td></tr><tr><td>FCR</td><td>0.17</td><td>0.91 ± 0.08</td></tr><tr><td>BR</td><td>0.14</td><td>0.95 ± 0.07</td></tr><tr><td>ECRL</td><td>0.24</td><td>0.89 ± 0.12</td></tr><tr><td>DA</td><td>0.18</td><td>0.93 ± 0.08</td></tr><tr><td>DP</td><td>0.43</td><td>0.78 ± 0.24</td></tr><tr><td>LD</td><td>0.32</td><td>0.81 ± 0.26</td></tr><tr><td>PM</td><td>0.37</td><td>0.82 ± 0.20</td></tr><tr><td>TBL</td><td>0.11</td><td>0.95 ± 0.05</td></tr><tr><td>BB</td><td>0.13</td><td>0.96 ± 0.06</td></tr><tr><td>TBS</td><td>0.21</td><td>0.87 ± 0.15</td></tr><tr><td><strong>SYNERGIES</strong></td><td></td><td></td></tr><tr><td>H<sub>1</sub></td><td>0.11</td><td>0.95 ± 0.06</td></tr><tr><td>H<sub>2</sub></td><td>0.12</td><td>0.94 ± 0.08</td></tr><tr><td>H<sub>3</sub></td><td>0.07</td><td>0.98 ± 0.02</td></tr><tr><td>H<sub>4</sub></td><td>0.13</td><td>0.95 ± 0.05</td></tr></tbody></table>

**Notes.**

Variance Ratio (VR) and r-Pearson (r) of Kinematics, torque, EMG patterns, synergies vectors ($W_s$) and synergy activation coefficients ($H_s$).

## LL

Four synergies were extracted for each participant (Fig. 5-LL) and explained 96 ± 0.02% of $VAF_{tot}$ and 95 ± 0.05% of $VAF_m$. The **first synergy ($S_1$)** involved hip flexors (RF, TF) and extensors (GMAX, BF, SEM), knee extensors (RF, VM, VL) and SOL and was active around the upper transition and the first part of FH (270-90°). The **second synergy ($S_2$)** involved ankle plantar-flexors (SOL, GL, GM), hip extensors (BF, SEM) and occurred during at mid-FH (around 90°). The **third synergy ($S_3$)** located around the bottom transition (180°), involved ankle plantar-flexors (GL) and hip flexors (TF). The **fourth synergy ($S_4$)** involved principally the TA and to a lesser extent the RF and TF during SH (180-0°).

## Interindividual variability

The results concerning the interindividual variability are described in Tables 1 and 2.

### Joint kinematics

UL showed slightly more angular kinematics variability than LL. The distal joint showed more variability than the proximal and intermediate joints for both UL and LL. Negligible time shifts in kinematic patterns were found for both limbs (Tables 1–2 - Kinematics & dynamics).
<table><caption>**Table 4** Lower limbs intra individual variability.</caption><thead><tr><th></th><th>VR<sub>intra</sub></th><th>r<sub>intra</sub><br>(mean ± sd)</th></tr></thead><tbody><tr><td colspan="3"><strong>KINEMATICS AND DYNAMICS</strong></td></tr><tr><td>Torque</td><td>0.02</td><td>1.00 ± 0.00</td></tr><tr><td>Hip</td><td>0.01</td><td>1.00 ± 0.00</td></tr><tr><td>Knee</td><td>0.00</td><td>1.00 ± 0.00</td></tr><tr><td>Ankle</td><td>0.14</td><td>0.96 ± 0.05</td></tr><tr><td colspan="3"><strong>EMG</strong></td></tr><tr><td>SOL</td><td>0.06</td><td>0.98 ± 0.02</td></tr><tr><td>GL</td><td>0.09</td><td>0.97 ± 0.04</td></tr><tr><td>GM</td><td>0.07</td><td>0.82 ± 0.04</td></tr><tr><td>TA</td><td>0.32</td><td>0.95 ± 0.30</td></tr><tr><td>GMAX</td><td>0.14</td><td>0.90 ± 0.06</td></tr><tr><td>RF</td><td>0.23</td><td>0.98 ± 0.13</td></tr><tr><td>VL</td><td>0.05</td><td>0.98 ± 0.03</td></tr><tr><td>VM</td><td>0.06</td><td>0.97 ± 0.03</td></tr><tr><td>BF</td><td>0.15</td><td>0.93 ± 0.07</td></tr><tr><td>SEM</td><td>0.22</td><td>0.89 ± 0.15</td></tr><tr><td>TF</td><td>0.21</td><td>0.87 ± 0.18</td></tr><tr><td colspan="3"><strong>SYNERGIES</strong></td></tr><tr><td>S<sub>1</sub></td><td>0.04</td><td>0.98 ± 0.02</td></tr><tr><td>S<sub>2</sub></td><td>0.05</td><td>0.98 ± 0.02</td></tr><tr><td>S<sub>3</sub></td><td>0.10</td><td>0.95 ± 0.05</td></tr><tr><td>S<sub>4</sub></td><td>0.26</td><td>0.88 ± 0.19</td></tr></tbody></table>

**Notes.**

Variance Ratio (VR) and r-Pearson (r) of Kinematics, torque, EMG patterns, synergies vectors (W<sub>s</sub>) and synergy activation coefficients (H<sub>s</sub>).

### Torque

UL showed more interindividual variability in torque production than LL (VR<sub>UL</sub> = 0.26, r<sub>UL</sub> = 0.72 ± 12 vs VR<sub>LL</sub> = 0.04, r<sub>LL</sub> = 0.97 ± 0.01 respectively, *p* < 0.05). A low time shift was found at the maximum of the cross-correlation for both limbs (t<sub>max-UL</sub> = 2.37 ± 2.46; t<sub>max-LL</sub> = 0.81 ± 0.78). These results indicates that interindividual variability in torque pattern mainly lies in differences in the pattern that is produced and not in a time shifting of similar patterns. Nevertheless, considering the amount of variability, even with a higher value for UL, a similar pattern of torque is produced by participants (Table 1, Fig. 3).

### EMG

Compared to LL, UL EMG patterns showed high interindividual variability for all muscles (VR<sub>EMG-UL</sub> = 0.68 ± 0.22 > VR<sub>EMG-LL</sub> = 0.32 ± 0.22, r<sub>EMG-UL</sub> = 0.36 ± 0.17 < r<sub>EMG-LL</sub> = 0.71 ± 0.24, *p* < 0.05, Tables 2–3-EMG). For LL, only GL, TA, BF, SEM and TF muscles showed high variability (VR > 0.37, r > 0.64 ± 0.28), these muscles being mainly involved in S<sub>3</sub> and S<sub>4</sub> and not associated to the production of propulsive torque. An important time shift was found at the maximum of cross-correlation function for UL (t<sub>max-EMG-UL</sub> = 15,01 ± 4% of the total cycle, Table 1-EMG) compared to LL (t<sub>max-EMG-LL</sub> = 7.9 ± 9.3%,
<figure><img src="image_6.png" /><figcaption>**Figure 3** Crank torque across the cycling revolution. Mean ± sd right crank torque of each subject for upper limbs (top) and lower limbs (bottom). A schematic representation of the cycling revolution. The blue point illustrates the FH to SH transition, and the yellow point illustrates the SH to FH transition or upper transition for upper and lower limbs respectively. FH, First Half; SH, Second Half.</figcaption></figure>

Full-size 🗠 DOI: 10.7717/peerj.13155/fig-3

$p < 0.05$, excepted for TA and TF, Table 2-EMG). See Figs. S3, and S4 for mean ± sd representation of each individual muscle pattern for UL and LL respectively.

## Synergies

The interindividual variability observed at the EMG level is reflected in the synergy activation coefficients, with greater interindividual variability for UL than LL ($VR_{Hs-UL} = 0.59 \pm 0.13 > VR_{Hs-LL} = 0,26 \pm 0.18$, $r_{Hs-UL} = 0.38 \pm 0.07 < r_{Hs-LL} = 0.74 \pm 0.17$, $p < 0.05$). A larger time shift was found for all UL synergies ($t_{max-Hs-UL} = 12.96 \pm 1.17 > t_{max-Hs-LL} = 5.31 \pm 3.73\%$, $p < 0.05$) despite a good index of similarity at the maximum of $r_{max}$ function ($r_{max-Hs-UL}: = 0.87 \pm 0.05$ and $r_{max-Hs-LL} = 0.93 \pm 0.06$). UL and LL $W_s$ showed similar correlation coefficients ($0.74 \pm 0.02$ vs $0.69 \pm 0.21$ for UL and LL respectively) excepted for $S_3$ of LL ($r_{S3} = 0.39 \pm 0.37$).

## Correlation in time shifts

The correlations in time shifting of the synergies were all significant for UL (Fig. S1, $r_{UL} > 0.81$, $p_{UL} < 0.00$) and only for $S_2$ and $S_3$ on LL (Fig. S2, $r = 0.70$, $p < 0.00$). This result indicates that all synergy activation coefficients are time-shifted by a similar amount for each participant and thus, despite the high interindividual variability in EMG patterns, participants produced the same coordination pattern (i.e., 4 similar synergies, produced at similar relative interval) but with delayed timings.

A weak negative correlation is observed between all UL synergy activation coefficients and RPM ($r_{RPM-H1} = -0.58$, $r_{RPM-H2} = -0.50$, $r_{RPM-H3} = -0.50$, $r_{RPM-H4} = -0.49$, $p < 0.00$) and to a lesser extent for LL synergy activation coefficients and RPM ($r_{RPM-H1} = -0.28$, $r_{RPM-H2} = -0.35$, $r_{RPM-H4} = -0.55$, $p < 0.00$). No correlation was found between $H_s$ and torque for UL.
<figure><img src="image_7.png" /><figcaption>**Figure 4** Kinematic pattern across the cycling revolution. Mean ± sd angle (degree) of proximal, in-termediate, and distal articulations of LL (left) and UL (right) torque of each subject for UL (top) and LL (bottom). The blue point illustrates the FH to SH transition, and the yellow point illustrates the SH to FH transition or upper transition for upper and lower limbs respectively. FH, First Half; SH, Second Half; UL, upper limbs; LL, lower limbs.</figcaption></figure>

Full-size

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

DOI: 10.7717/peerj.13155/fig-4

## Intraindividual variability

The results concerning the interindividual variability are described in Tables 3 and 4. Intraindividual variability was inferior to interindividual variability and showed a good repeatability of kinematics, torque, EMG, and synergy activation coefficients for each participant.

## DISCUSSION

The first objective of this study was to investigate whether common muscular coordination processes related to mechanical output were observable during upper limbs (UL) and lower limbs (LL) cycling. To this aim, muscles coordination was analyzed through synergies in relation to the torque produced during cycling for both limbs. Even though the movement is strongly constrained, and the relative task intensity is similar for UL and LL, the torque outputs produced by both limbs showed strong differences. The results showed that cycling with the LL resulted in a reciprocal contribution of one leg then the other and a production of propulsive torque during half the cycle only. Cycling with the UL resulted in a different torque pattern arising from a simultaneous contribution of both arms, each producing
<figure><img src="image_9.png" /><figcaption>**Figure 5** Muscle's synergy parameters and variability. Muscle synergies for LL (left) and UL (right). (A) Muscle's synergies vectors ($W_s$) for each muscle (mean ± sd). (B) Mean synergy activation coefficient ($H_s$) for each subject. (C) Boxplot of activation coefficient peak occurrences across the cycle revolution during 25 cycles for each subject (each row). The boxplot showed interindividual variability in the occurrence of the peaks for each activation coefficient: the participants 1 and 2 showed synchronous synergy occurrence, the participant 3 showed an early peak for each synergy and the participants 4, 5 and 6 showed latter peaks for each synergy. The detail of $r_{inter-H}$ showed high correlation between participants 1 and 2 ($r_{inter-1/2}$ = 0.95) and between participants 4, 5 and 6 ($r_{inter-4/5/6}$ > 0.94). The third participant showed weak or negative correlation with all the others.</figcaption></figure>

Full-size 🗠 DOI: 10.7717/peerj.13155/fig-5

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

propulsive torque during most of the cycle. Despite these differences in mechanical output, the same number of muscle synergies (4) were extracted to account for the variations in muscle activations of 11 muscles for each limb. The simultaneous analysis of the synergies related to the mechanical output showed different functional implications in the movement production between UL and LL.

Our second objective was to compare the interindividual and intraindividual variability of the UL and LL torque, joint kinematics, EMG, and synergy-related variables. The results showed a higher variability of the UL in kinematics, torque, EMG and in muscle coordination compared to LL. More precisely, the variability in EMG patterns was mainly reflected in the temporal components of muscle synergies: each synergy activation coefficient ($H_s$) of UL showed important time shifts between participants. Nevertheless,
despite this significant variability, each participant, when analyzed individually, produced
reproducible synergy activation coefficients.

Four synergies were extracted for both UL and LL, despite a large interindividual variability in EMG patterns for UL. Our results on muscle synergy vectors showed close values of correlation between participants for both limbs. The *r*-value obtained are lower than those obtained by *De Marchis et al. (2013)* or *Hug et al. (2010)*, which can certainly be explained by differences in the methodologies used (average EMG patterns vs 25 consecutive cycles) or by the type of population chosen (expert vs non-expert). It is thus reasonable to consider that the spatial constitution of the synergies is shared between individuals. This indicates that the structure of synergies was robust between participants and support the assumption of a spatially fixed synergies.

Concerning the temporal aspects of the synergies, our results showed a strong interindividual variability in the EMG patterns which is reflected in the activation coefficients of the synergies, as constated by several authors (Torres-Oviedo & Ting, 2007; Safavynia & Ting, 2012; De Marchis et al., 2013). More specifically, each UL temporal activation coefficients showed large interindividual variability, mainly resulting from differences in the onset in apparition of each synergy pattern between participants rather than differences in patterns produced (Table 1). These results showed that same four synergy activation patterns were produced between participants, with a constant interval between each synergy (Fig. S1), but the expression of these four synergies was highly variable in the timing of apparition during the cycle. In comparison, LL synergies, particularly Syn1 and Syn2 showed low inter-individual variability which is an important difference between UL and LL.

The four synergies extracted during LL cycling agree with the synergies described in the literature in novice participants (De Marchis et al., 2013) or professional cyclists (Hug et al., 2010). S₁ and S₂ were associated with the production of propulsive torque during the power phase (from top dead center to bottom dead center). These two synergies strongly mobilized the hip and knee extensor muscles and were therefore both involved in the same mechanical function: producing propulsive torque by pushing the pedal with a leg extension (Fig. 4 - left). S₃ and S₄ were expressed during the resistive torque production phase and may be involved to optimize the lower transition and minimize the resistive torque induced by the legs weight during the flexion of the limb (Fig. 4 - left). Importantly, a merging of S₂ and S₃ is frequently observed in literature with expertise. S₃ is associated with tangential application of the force on the pedal at the end of the push. A merging of these synergies is currently depicted as better orientation of the force on the crank/pedal at the end of the propulsive phase. It is therefore not surprising that both showed variability and correlation in timing of apparition. Being related, if S₂ appears earlier, S₃ logically also appears earlier. The same explanation stands for the variability observed at the level of synergy vector (W₃).

To our knowledge, a single study looked at UL muscle synergies (Botzheim et al., 2021) and observed two distinct synergies and attributed these synergies to the two phases in the movement: pushing and pulling. If that was enough to describe bimanual coordination, it is impossible to rely on this observation to describe the mechanical output. Indeed,
these authors did not used crank torque measurement nor indicates for the output power.
Thanks to our methodology our study distinguishes a total of 4 synergies, three (S<sub>1</sub> to S<sub>3</sub>)
corresponded to propulsive torque and one (S<sub>4</sub>) to the reduction of resistive torques. S<sub>1</sub>
primarily involved the shoulder flexors and elbow extensor muscles meaning that the first
increase in torque is related to a pushing action on the pedals. This synergy is simultaneous
to shoulder antepulsion and elbow extension (Fig. 4 -right). S<sub>3</sub> involved the elbow flexors,
the wrist/finger flexors, and the shoulder extensor. This synergy is simultaneous with
the peak of propulsive torque production and indicates that pulling is primarily used to
produce torque. This synergy is simultaneous to reduction of shoulder antepulsion and
elbow flexion (i.e., reducing limb length). Our results showed that a third synergy (S<sub>2</sub>) was
involved in the production (or conservation) of propulsive torque during the transition
from FH to SH. This synergy involved the wrist/finger flexors and the *latissimus dorsi*, the
muscle responsible for the lowering of the arm. It seems that this synergy makes it possible
to produced/conserved propulsive torque by lowering the arm in the transition from FH
to SH. The simultaneous reduction of shoulder antepulsion, while elbow is still extending
is in line with this explanation. Furthermore, this synergy probably takes advantage of the
gravitational force to optimize the lowering action by using the weight of the arm. However,
given the high relative intensity of the cycling movement, the gravitational influence is
probably minimal. As for the synergies observed on LL, the synergy S<sub>4</sub> may be responsible
for minimizing the resistive torque around the transition from SH to FH.

These results demonstrated that for a similar task demand (i.e., moving the crank at an identical relative intensity), the muscular coordination of the UL is organized around several mechanical functions (i.e., pushing, lowering, and pulling) to produce propulsive torque while the LL produce propulsive torque using a single main function (i.e., push). These differences lead to a higher mechanical efficiency for the UL pedaling (i.e., lower maximal torque and longer propulsive phase). Such differences may result from the daily use of UL in tasks of reaching, gripping, and manipulating in 3D space. The musculoskeletal system could be optimized to both extend and flex the arm. The LL are used mainly in locomotion tasks thanks to a musculoskeletal system optimized for antigravity functions (therefore pushing upward and propelling the upper body). Consequently, with torque production over a greater proportion of the pedaling cycle, which is synonymous with greater mechanical efficiency in performing the movement, UL demonstrates an important adaptation capacity to perform this movement for which our subjects and most of the population are non-expert.

A main difference between UL and LL is the level of variability observed. The results
showed a higher variability for the UL, both for torque (*performance* variability) and for
EMG and muscles coordination (*realization* variability). Nevertheless, when considering
the variability for the UL, our results showed less variability in *performance* in comparison
with *realization*. Thus, a similar torque pattern is produced between participants, with
minor differences in shape and amplitude. On the contrary, our results showed that the
variability of *realization* is large and lies specifically in the temporal aspects. This observation
is surprising since one would have expected that differences in the timing onsets of the
synergies would imply similar differences in the timing onsets in torques production.
<table><caption>**Table 5** Cadence of each participant for UL (right) and LL (right).</caption><thead><tr><th></th><th>UL<br>(RPM, mean ± sd)</th><th>LL<br>(RPM, mean ± sd)</th></tr></thead><tbody><tr><td>P01</td><td>90.6 ± 0.6</td><td>82.5 ± 0.3</td></tr><tr><td>P02</td><td>78.1 ± 1.3</td><td>89.6 ± 2.0</td></tr><tr><td>P03</td><td>99.6 ± 1.3</td><td>64.1 ± 0.7</td></tr><tr><td>P04</td><td>71.4 ± 0.5</td><td>82.1 ± 1.8</td></tr><tr><td>P05</td><td>92.6 ± 0.9</td><td>102.0 ± 0.8</td></tr><tr><td>P06</td><td>81.7 ± 0.9</td><td>80.0 ± 1.1</td></tr></tbody></table>

**Notes.**

LL showed more interindividual differences in cadence (ranging from 64.1 ± 0.7 to 102.0 ± 0.8) than UL (ranging from 71.4 ± 0.5 to 99.6 ± 1.3). UL: upper limbs, LL, lower limbs.

This suggests that there was no direct link between the variability observed in muscle coordination and the variability of the resulting torque.

Several hypotheses can highlight these discrepancies: first, the difference in timing can be explained by differences in self-selected cadence as shown by the correlation results between the appearance of synergies and cadence. Nevertheless, the differences in cadence cannot totally explain the temporal variability of the synergy activation patterns, especially since on LL the results showed a wider range of cadence between participants and no phase shift in synergies (Table 5). Second, greater level of variability in the UL can also be explained by differences in expertise between the UL and LL cycling. Although the participants were chosen novices for both tasks, LL cycling is associated to greater expertise since it has been shown that cycling shares common neural circuits with walking (Barroso et al., 2014), which ultimately induces an expertise in the realization of cycling since all the participants knew how to walk. However, in the literature, it has been shown during a bench press exercise that a greater inter-individual variability in synergy activation patterns is synonym of expertise instead of non-experts (Kristiansen et al., 2015). Thus, a difference in expertise does not seem sufficient to explain the differences in realization variability observed between UL and LL.

Before that conclusion can be made on the topic of expertise, further studies are needed to understand whether expertise (multiple sessions of UL cycling) tends to reduce the realization variability between individuals or whether individual strategies persist. Nevertheless, from our results for “non-experts” (as most patients and athletes), a flexibility in the use of synergies is adopted when performing the movement with UL unlike LL for which robust coordination is adopted by everyone.

## CONCLUSION

This study compared the muscles coordination during cycling with Upper Limbs (UL) and Lower Limbs (LL). Torque profiles differed between UL and LL and UL EMG patterns showed a larger interindividual variability than LL. Nevertheless, as already reported for LL, a functional organization of the UL muscle activations do exist during UL cycling and is shared among participants. UL cycling is achieved with a higher flexibly in the temporal onsets of the synergies compared to LL cycling besides this higher flexibility in
realization was not reflected in the mechanical output. The overall results of this study could be considered when combining UL and LL cycling for an assessment/rehabilitation process. In particular for UL, our results showed that an individualized approach should be recommended to evaluate the coordination state of the patient and to adapt the exercise parameter.

**ADDITIONAL INFORMATION AND DECLARATIONS**

**Funding**

The authors received no funding for this work.

**Competing Interests**

The authors declare there are no competing interests.

**Author Contributions**

* Théo Cartier conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.

* Laurent Vigouroux, Elke Viehweger and Guillaume Rao conceived and designed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

**Human Ethics**

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

CERSTAPS approved the study (IRB00012476-2020-30-11-75).

**Data Availability**

The following information was supplied regarding data availability:

The data is available at FigShare: Cartier, Théo (2022): Data_2_Share.zip. figshare.

Dataset. https://doi.org/10.6084/m9 figshare.19099754.v1.

**Supplemental Information**

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/peerj.13155#supplemental-information.

**REFERENCES**

Barroso FO, Torricelli D, Moreno JC, Taylor J, Gomez-Soriano J, Bravo-Esteban E, Piazza S, Santos C, Pons JL. 2014. Shared muscle synergies in human walking and cycling. *Journal of Neurophysiology* 112:1984–1998 DOI 10.1152/jn.00220.2014.

Bartlett R, Wheat J, Robins M. 2007. Is movement variability important for sports biomechanists? *Sports Biomechanics* 6:224–243 DOI 10.1080/14763140701322994.

Barzi Y, Zehr EP. 2008. Rhythmic arm cycling suppresses hyperactive soleus H-reflex amplitude after stroke. *Clinical Neurophysiology* 119:1443–1452 DOI 10.1016/j.clinph.2008.02.016.
Bernstein N. 1967. *The co-ordination and regulation of movements*. Oxford: Pergamon Press.

Botzheim L, Laczko J, Torricelli D, Mravcsik M, Pons JL, Oliveira Barroso F. 2021. Effects of gravity and kinematic constraints on muscle synergies in arm cycling. *Journal of Neurophysiology* 125:1367–1381 DOI 10.1152/jn.00415.2020.

Burden AM, Trew M, Baltzopoulos V. 2003. Normalisation of gait EMGs: a re-examination. *Journal of Electromyography and Kinesiology* 13:519–532 DOI 10.1016/S1050-6411(03)00082-8.

Chaytor CP, Forman D, Byrne J, Loucks-Atkinson A, Power KE. 2020. Changes in muscle activity during the flexion and extension phases of arm cycling as an effect of power output are muscle-specific. *PeerJ* 8:e9759 DOI 10.7717/peerj.9759.

d'Avella A, Saltiel P, Bizzi E. 2003. Combinations of muscle synergies in the construction of a natural motor behavior. *Nature Neuroscience* 6:300–308 DOI 10.1038/nn1010.

De Marchis C, Schmid M, Bibbo D, Bernabucci I. 2013. Inter-individual variability of forces and modular muscle coordination in cycling: a study on untrained subjects. *Human Movement Science* 32:1480–1494 DOI 10.1016/j.humov.2013.07.018.

Elmer SJ, Danvind J, Holmberg H-C. 2013. Development of a novel eccentric arm cycle ergometer for training the upper body. *Medicine & Science in Sports & Exercise* 45:206–211 DOI 10.1249/MSS.0b013e318269c79c.

Faria EW, Parker DL, Faria IE. 2005. The science of cycling: physiology and training. Part 1. *Sports Medicine* 35:285–312 DOI 10.2165/00007256-200535040-00002.

Fujiwara T, Liu M, Tanuma A, Hase K, Tsuji T. 2005. Pedaling exercise for neuromuscular re-education: a review. *Critical Reviews in Physical and Rehabilitation Medicine* 17:163–178 DOI 10.1615/CritRevPhysRehabilMed.v17.i3.10.

Hübner-Wozniak E, Kosmol A, Lutoslawska G, Bem E. 2004. Anaerobic performance of arms and legs in male and female free style wrestlers. *Journal of Science and Medicine in Sport* 7:473–480 DOI 10.1016/S1440-2440(04)80266-4.

Hug F, Dorel S. 2009. Electromyographic analysis of pedaling: a review. *Journal of Electromyography and Kinesiology* 19:182–198 DOI 10.1016/j.jelekin.2007.10.010.

Hug F, Turpin NA, Couturier A, Dorel S. 2011. Consistency of muscle synergies during pedaling across different mechanical constraints. *Journal of Neurophysiology* 106:91–103 DOI 10.1152/jn.01096.2010.

Hug F, Turpin NA, Guével A, Dorel S. 2010. Is interindividual variability of EMG patterns in trained cyclists related to different muscle synergies? *Journal of Applied Physiology* 108:1727–1736 DOI 10.1152/japplphysiol.01305.2009.

Hug F, Vogel C, Tucker K, Dorel S, Deschamps T, Le Carpentier É, Lacour-paille L. 2019. Individuals have unique muscle activation signatures as revealed during gait and pedaling. *Journal of Applied Physiology* 127:1165–1174 DOI 10.1152/japplphysiol.01101.2018.

Kaupp C, Pearcey GEP, Klarner T, Sun Y, Cullen H, Barss TS, Zehr EP. 2018. Rhythmic arm cycling training improves walking and neurophysiological integrity in chronic stroke: the arms can give legs a helping hand in rehabilitation. *Journal of Neurophysiology* 119:1095–1112 DOI 10.1152/jn.00570.2017.
**Kristiansen M, Madeleine P, Hansen EA, Samani A. 2015.** Inter-subject variability of muscle synergies during bench press in power lifters and untrained individuals: muscle coordination during bench press. *Scandinavian Journal of Medicine & Science in Sports* **25**:89-97 DOI 10.1111/sms.12167.

**Latash ML. 2012.** The bliss (not the problem) of motor abundance (not redundancy). *Experimental Brain Research* **217**:1-5 DOI 10.1007/s00221-012-3000-4.

**Mezzarane RA, Nakajima T, Zehr EP. 2014.** After stroke bidirectional modulation of soleus stretch reflex amplitude emerges during rhythmic arm cycling. *Frontiers in Human Neuroscience* **8**:136 DOI 10.3389/fnhum.2014.00136.

**Mornieux G, Stapelfeldt B, Gollhofer A, Belli A. 2008.** Effects of Pedal Type and Pull-Up Action during Cycling. *International Journal of Sports Medicine* **29**:817-822 DOI 10.1055/s-2008-1038374.

**Oja P, Titze S, Bauman A, de Geus B, Krenn P, Reger-Nash B, Kohlberger T. 2011.** Health benefits of cycling: a systematic review: cycling and health. *Scandinavian Journal of Medicine & Science in Sports* **21**:496-509 DOI 10.1111/j.1600-0838.2011.01299.x.

**Safavynia SA, Ting LH. 2012.** Task-level feedback can explain temporal recruitment of spatially fixed muscle synergies throughout postural perturbations. *Journal of Neurophysiology* **107**:159-177 DOI 10.1152/jn.00653.2011.

**Sangeux M, Polak J. 2015.** A simple method to choose the most representative stride and detect outliers. *Gait & Posture* **41**:726-730 DOI 10.1016/j.gaitpost.2014.12.004.

**Sitko S, Cirer-Sastre R, Corbi F, López-Laval I. 2020.** Power assessment in road cycling: a narrative review. *Sustainability* **12**:5216 DOI 10.3390/su12125216.

**Torres-Oviedo G, Ting LH. 2007.** Muscle synergies characterizing human postural responses. *Journal of Neurophysiology* **98**:2144-2156 DOI 10.1152/jn.01360.2006.

**Turpin NA, Uriac S, Dalleau G. 2021.** How to improve the muscle synergy analysis methodology? *European Journal of Applied Physiology* **121**:1009-1025 DOI 10.1007/s00421-021-04604-9.

**Zehr EP. 2005.** Neural control of rhythmic human movement: the common core hypothesis. *Exercise and Sport Sciences Reviews* **33**(7):54-60.

**Zehr EP, Chua R. 2000.** Modulation of human cutaneous reflexes during rhythmic cyclical arm movement. *Experimental Brain Research* **135**:241-250 DOI 10.1007/s002210000515.

**Zehr EP, Loadman PM, Hundza SR. 2012.** Neural control of rhythmic arm cycling after stroke. *Journal of Neurophysiology* **108**:891-905 DOI 10.1152/jn.01152.2011.

**Zhou R, Alvarado L, Kim S, Chong SL, Mushahwar VK. 2017.** Modulation of corticospinal input to the legs by arm and leg cycling in people with incomplete spinal cord injury. *Journal of Neurophysiology* **118**:2507-2519 DOI 10.1152/jn.00663.2016.

**Zinner C, Morales-Alamo D, Ørtenblad N, Larsen FJ, Schiffer TA, Willis SJ, Gelabert-Rebato M, Perez-Valera M, Boushel R, Calbet JAL, Holmberg H-C. 2016.** The physiological mechanisms of performance enhancement with sprint interval training differ between the upper and lower extremities in humans. *Frontiers in Physiology* **7**:426 DOI 10.3389/fphys.2016.00426.