Skip to main content

Unilateral symptomatic Achilles tendinopathy has limited effects on bilateral lower limb ground reaction force asymmetries and muscular synergy attributes when walking at natural and fast speeds

Abstract

Background

Achilles tendinopathy (AT) may affect ground reaction force (GRF) and muscle synergy (MS) during walking due to pain, biological integrity changes in the tendon and neuroplastic adaptations. The objective of this study was to compare GRF asymmetries and MS attributes between symptomatic and asymptomatic lower limbs (LL) during walking at natural and fast speeds in adults with unilateral AT.

Methods

A convenience sample consisting of twenty-eight participants walked on an instrumented treadmill at natural (1.3 m/s) and fast (1.6 m/s) speeds. Peak GRF were measured in mediolateral, anteroposterior and vertical directions. Individualized electromyography (EMG) activation profiles were time- and amplitude-normalized for three consecutive gait cycles and MS were extracted using non-negative matrix factorization algorithms. MS were characterized by the number, composition (i.e., weighting of each muscle) and temporal profiles (i.e., duration and amplitude) of the MS extracted during walking. Paired Student’s t-tests assessed peak GRF and MS muscle weighting differences between sides whereas Pearson correlation coefficients characterized the similarities of the individualized EMG and MS activation temporal profiles within sides.

Results

AT had limited effects on peak GRF asymmetries and the number, composition and temporal profiles of MS between symptomatic and asymptomatic LL while walking on a level treadmill at natural and fast speeds. In most participants, four MS with a specific set of predominantly activated muscles were extracted across natural (71 and 61%) and fast (54 and 50%) walking speeds for the symptomatic and asymptomatic side respectively. Individualized EMG activation profiles were relatively similar between sides (r = 0.970 to 0.999). As for MS attributes, relatively similar temporal activation profiles (r = 0.988 to 0.998) and muscle weightings (p < 0.05) were found between sides for all four MS and the most solicited muscles. Although the faster walking speed increased the number of merged MS for both sides, it did not significantly alter MS symmetry.

Conclusion

Faster walking speed increased peak GRF values but had limited effects on GRF symmetries and MS attribute differences between the LL. Corticospinal neuroplastic adaptations associated with chronic unilateral AT may explain the preserved quasi-symmetric LL motor control strategy observed during natural and fast walking among adults with chronic unilateral AT.

Peer Review reports

Introduction

The Achilles tendon is the longest and most powerful tendon in the human body [1]. Tensile forces are transmitted through the tendon following contraction of the triceps surae muscles and enable ankle plantar flexion and related movements. Active plantar flexion is required for propulsion of the foot when generating functional movements, such as those required during walking and running [2]. To optimize the transmission of force during these functional movements, the Achilles tendon uses its potential to store and release elastic energy during the stretching and shortening phases, respectively [3]. For optimal function, a healthy tendon is crucial to withstand high tensile forces and to protect the triceps surae muscles from injury [2]. These high, repetitive and rapidly rising tensile forces transiting through the AT increase the risk of altering the biological integrity of the Achilles tendon [4]. In many cases, failure to achieve adaptive and restorative healing responses over time following Achilles tendon injury causes symptomatic and chronic Achilles tendinopathy (AT) [5,6,7].

The presence of localized pain, changes in the biological integrity of the Achilles tendon, and central adaptations typically characterize a symptomatic and chronic AT that can affect lower limb (LL) movement strategies during walking and running [8,9,10,11]. Runners with AT showed altered amplitude and duration of EMG activity at the lower limbs muscles during the different phases of running [12,13,14,15,16]. These muscular recruitment changes may also interact with changes in LL ground reaction forces (GRF) measured under the feet [11]. However, no significant differences in magnitude or timing of GRF were previously reported between adults with symptomatic AT and healthy counterparts [12, 16]. Overall, these movement strategy changes may consolidate over time and prompt cortical alterations related to central nervous system (CNS) motor planning, such as recruitment of a protective strategy with altered movements or reorganization of cortical regions [17,18,19].

To gain insight into the potential relationship between the above-mentioned changes in movement strategies and the CNS, a well-established hypothesis must be considered, namely that the CNS does not control each muscle individually, but instead adopts strategies that simplify the control of complex movements [20]. This hypothesis stipulates that the CNS has an organizational structure that synchronizes the amplitude, timing and duration of muscle activity to support coordinated movements [21]. In fact, the CNS initiates motor commands to select specific muscles to be activated at different intensities and in a coordinated way (i.e., motor modules or “muscle synergies” (MS)) by activating specific groups of motor neurons. These MS are thought to ease the potential complexity of the distinct neuronal activation of several individual muscles during coordinated movements.

Characterizing MS in terms of the number of motor modules, composition (i.e., weighting of each muscle per motor module) and temporal profile (i.e., duration and amplitude) during different functional activities is feasible, using non-negative matrix factorization (NNMF) algorithms [22, 23]. Such a methodological approach has allowed researchers to identify four MS typically observed during walking in healthy adults [24,25,26]. Each of these MS is activated at distinct phases of the gait cycle, ensuring a specific biomechanical function (Fig. 1): During weight acceptance, the gluteus medius, vastus medialis and rectus femoris are activated for leg stabilization (MS1). During pushoff, the soleus and medial gastrocnemius are activated for forward propulsion (MS2). During the early swing phase, the tibialis anterior and rectus femoris are activated for swing initiation (MS3). During the terminal swing phase, the semitendinosus and biceps femoris are activated for leg deceleration (MS4).

Fig. 1
figure 1

Muscular synergies identified by NNMF during walking among healthy adults

Some MS attributes (i.e., number, composition and temporal profile) are now considered to be potential determinants and predictors of intermuscular coordination and gait performance among adults with neurological [27,28,29] or musculoskeletal impairments [30,31,32]. In adults with neurological impairments (i.e., stroke or spinal cord injury), some authors have observed that fewer MS are required to account for muscle activation during walking, suggesting a reduction in overall motor activation complexity. In adults with musculoskeletal impairments, authors have observed that adults with gluteal tendinopathy had changes in MS composition with the gluteus minimus, gluteus medius, and tensor fascia latae, notably contributing to the MS related to single leg support [33]. Authors have observed that adults with unilateral anterior cruciate ligament deficiency have altered MS activation profiles when stabilizing the knee joint with a longer activation period of the tibialis anterior, quadriceps, and hamstring muscles during stance [32]. To our knowledge, no study has characterized the effect of chronic pain and changes in the Achilles tendon on GRF and MS during walking in adults with unilateral AT.

Hence, the primary objective of this study was to compare GRF symmetries (i.e., mediolateral, posteroanterior and vertical components) and MS attribute differences (i.e., number, composition and temporal profile) between symptomatic and asymptomatic LL during walking in adults with unilateral symptomatic AT. The secondary objective was to verify if increased walking speed would reveal additional or larger differences in terms of GRF asymmetries and MS attributes between LL. For the primary objective, it was hypothesized that peak GRF in the mediolateral, posteroanterior and vertical directions would be reduced at the symptomatic LL compared to the asymptomatic LL [10]. It was also anticipated that the number of synergies would be reduced at the symptomatic LL compared to the asymptomatic LL to reflect changes in motor control complexity due to AT [34]. In addition, a change in the motor recruitment strategy of the hip (gluteus medius) and knee (vastus medialis) stabilizer muscles was anticipated and expected to translate into MS1 composition and temporal profile differences with the asymptomatic LL. For the secondary objective, it was hypothesized that peak GRF asymmetries would be greater at fast speed compared to natural speed whereas MS attribute differences between the asymptomatic and symptomatic LL previously described would be amplified at fast speed compared to natural speed [26].

Methodology

Participants

A convenience sample consisting of twenty-eight (n = 28) adults with unilateral, symptomatic AT of the Achilles tendon participated in this cross-sectional study. The convenience sampling technique was selected based on previous LL studies evaluating changes in muscles synergies following induced or acquired pain [30, 32, 33]. To be included, potential participants had to report unilateral, localized pain at the insertion or midportion of the Achilles tendon for more than 3 months [7], experience pain on palpation at the enthesis or midportion of the Achilles tendon, and obtain a score lower than 90 out of 100 on the Victorian Institute of Sport Assessment-Achilles Questionnaire (VISA-A) [35]. Presence of AT was confirmed by both the presence of pain upon palpation of the tendon and the observation of tendon alteration (e.g., increased tendon thickness, reduced echogenicity) during musculoskeletal ultrasound imaging. Potential participants who reported bilateral pain during sport activities, had a body mass index (BMI) greater than 30.0 kg/m2, reported a history of Achilles tendon rupture or repair, were diagnosed with a metabolic, neurologic or systemic inflammatory disease, or had received any type of injection at the Achilles tendon in the past year were excluded.

Questionnaires

Basic sociodemographic data (e.g., age, sex), anthropometric data (e.g., height, mass) and AT-related information (e.g., affected side, time since onset of symptoms, etc.) were initially collected. Participants then completed two patient-reported outcome measures questionnaires: the VISA-A and the Lower Extremity Functional Scale (LEFS). The VISA-A includes eight questions targeting three dimensions: localized pain at the Achilles tendon, function in daily life and participation in sports activities [35,36,37]. The VISA-A questionnaire scores range from 0 to 100, with a low score indicating greater severity. The questionnaire is available in English and French and is valid, reliable, and sensitive to change with a minimal clinically important difference between 6 and 20 points [35, 38]. The LEFS questionnaire assesses function in adults with musculoskeletal disorders affecting the LL, including AT [39]. It includes 20 questions measuring the level of difficulty encountered when performing activities of daily living and sports. The LEFS has a maximum score of 80, with a high score indicating a high functional level. The LEFS, which is also available in English and French, is valid, reliable, and sensitive to change with a minimal clinically important difference of 9–12 points [38]. The VISA-A and LEFS were completed by participants either on paper or electronically via Lime Survey® in the language of their choice (i.e., French or English).

Musculoskeletal ultrasound imaging

One trained evaluator (ML), who is an experienced physical therapist, performed all components of both the clinical and imaging evaluations. For the latest, all ultrasound images of the Achilles tendon were recorded with a HD 11XE 1.0.6 ultrasonography system (Phillips Medical Systems, Bothell, WA), using the brightness mode and a 5–12 MHz linear array transducer with a 5-cm wide footprint. A previously described standardized protocol [40] was used for image acquisition. All parameters affecting image quality remained constant for all participants during each data collection session (i.e., probe frequency set at 12 MHz; depth = 2 cm; gain = 85; unique focus zone adjusted to 0.5 cm at the Achilles tendon; neutral time gain compensation). The most painful area along the symptomatic tendon was first located by palpation before being marked on the skin and mirrored on the asymptomatic side. Three (n = 3) images centered on the marked site were acquired in the longitudinal plane per side.

All images were analyzed using a custom program developed using MATLAB’s Image Processing Toolbox™ (MathWorks Inc., Natick, Ma, USA) [40, 41] to extract a standardized dataset of ultrasound biomarkers (i.e., mean thickness, echogenicity, skewness, mean homogeneity, and homogeneity at 0° and 90°) [42]. Images of symptomatic tendons were expected to show increased thickness, decreased echogenicity, increased skewness, increased mean homogeneity, and increased homogeneity at 90° compared to asymptomatic tendons [42].

Ground reaction forces

GRF were recorded in mediolateral (GRFML), posteroanterior (GRFPA) and vertical (GRFV) directions by a fully instrumented dual-belt treadmill with separate embedded force plates (Bertec, Columbus, USA; TM-09). Selected peak GRFML, GRFPA and GRFV were based on curves previously described among healthy adults [43] (Fig. 2). The GRFML curve for normal walking contains one lateral and two medial peaks: the lateral thrust force (LTF) at heel contact, followed by the medial impact force (MIF) during weight acceptance, and lastly the medial propulsive force (MPF) during terminal stance. The GRFPA curve for normal walking contains two peaks: the horizontal braking force (HBF) upon weight acceptance and the horizontal propulsive force (HPF) during the pushoff phase. The GRFV curve for normal walking contains two peaks and a trough: the vertical impact force (VIF) during weight acceptance, the minimal vertical peak (MVP) during midstance and the vertical propulsive force (VPF) during pushoff.

Fig. 2
figure 2

Selected ground reaction forces (GRF) in the a medio-lateral (GRFML), b postero-anterior (GRFPA) and c vertical (GRFV) directions

Surface electromyography

Bilateral surface EMG of the rectus femoris, vastus medialis, tibialis anterior, medial gastrocnemius, soleus, semitendinosus, biceps femoris and gluteus medius was recorded using a portable telemetric system (NORAXON USA Inc., Scottsdale, Arizona; Telemyo 900) with a frequency of 1200 Hertz (Hz). Self-adhesive surface electrodes (Ag/AgCl; Ambu_BlueSensor M) were placed in a bipolar configuration with a 1-cm inter-electrode distance over the muscle belly, perpendicular to muscle fibers orientation. Each skin site was previously shaved and cleaned with alcohol [44] in accordance with the SENIAM recommendations (refer to www.seniam.org). EMG signals were visually inspected during static voluntary contractions performed against manual resistance according to a standardized protocol [45]. The recorded EMG was filtered using a fourth-order, Butterworth, zero-lag bandpass filter with cut-off frequencies set at 20 Hz and 400 Hz. Subsequently, EMG signals were root mean squared with a centered 250 msec moving window to finally generate linear envelopes.

Gait assessment

Participants walked on the above-described instrumented treadmill. After a 3-minute familiarization walking period on the treadmill during which participants were asked to walk at a self-selected comfortable speed, participants walked at a natural speed (1.3 m/s) for 60 s and thereafter at a fast speed (1.6 m/s) for 60 s. Since the mean self-selected comfortable speed (1.3 ± 0.1 m/s) of participants was similar to the prescribed natural speed (1.3 m/s), only the latest natural speed was considered for the analyses. A rest period was provided if required between the tasks. Pain was monitored for the symptomatic and asymptomatic tendons while walking using a 10-point visual analog scale with the ethical perspective of not increasing pain to a level greater than 5/10 in accordance with the principle of nonmaleficence [46]. EMG and GRF data were collected for the last 40 s of each trial, and an average of 3 consecutive gait cycles (two consecutive heel strikes for the same leg) were selected according to the minimal EMG variation coefficient and then used for analysis [26]. Each GRF value and EMG envelope was time-normalized (0 to 100% in 1% increments) relative to each full gait cycle, then averaged together. Each GRF value was then normalized for each participant’s body weight (BW). The amplitude of each muscle was also normalized from its peak value (0 to 1), resulting in an 8 X 101 experimental EMG (EMGexp) matrix. The mean and standard deviation (SD) of the BW-normalized GRF and amplitude-normalized EMG during the time-normalized gait cycle were calculated across subjects.

Muscle synergies analysis

MS during walking were extracted by a NNMF algorithm using the Gamma model based on J divergence, a reliable method found to generate high coefficients of correlation and confidence levels [22] and used in previous studies to extract MS during various walking conditions [26, 47]. For each participant, the EMGexp data matrix was subjected to the NNMF algorithm. The NNMF algorithm broke down the EMGexp into two matrices and a reconstruction error (ε) (Eq. 1). The first matrix represents the muscle weighting (W), consisting of the contribution of each muscle (m) within each muscle synergy (s) (Eq. 2). The second matrix represents the activation timing profile (H), consisting of the muscle synergy (s) during the gait cycle for each time-normalized point (t) (Eq. 3). Agreement between EMGrec and EMGexp was then evaluated using the variance accounted for (VAF) criterion (Eq. 4).

$${EMG}_{exp}=W\times H+\varepsilon$$
(1)
$$\mathrm{W}=m\times s$$
(2)
$$\mathrm{H}=s\times t$$
(3)
$$VAF=1-\left(\sum {\left({EMG}_{exp}-{EMG}_{rec}\right)}^2/\sum {EMG_{exp}}^2\right)$$
(4)

The number of MS was determined by the least number of synergies that could explain a VAF for each muscle (VAFm) greater than 90% and a global VAF (VAFg) greater than 80% [48]. Whenever these criteria were reached, the reconstruction was deemed valid and the computation stopped. When the absolute difference of the coefficient of determination between the current and last pass was lower than 1 × e− 8 for 20 consecutive passes, or after 500 passes were run without convergence, the algorithm stopped. This procedure was done 20 times, and the result of the lowest reconstruction error with the lowest number of synergy modules within the validation criteria were considered adequate.

Muscle synergies were sorted out based on similarities in muscle weighting across participants and walking speeds using cosine similarity [49]. The inner product of the compared muscle synergy vectors were calculated, and the cosine angle between those synergies were measured. Sorting was performed by grouping muscle synergies based on the values of cosine similarity against a healthy reference [24,25,26]. The sequence of two muscle’s weightings values to be compared are viewed as vectors and the dot product of the vectors are divided by the product of their length providing an index of similarity with a maximal unitary value. Whenever the cosine similarities of W between the reference and other MS were over 0.80, MS were considered similar [50, 51]. Likewise, whenever two MS at the same walking speed were classified into the same MS group, these two synergies were considered to be merged together. The MS with the lowest cosine similarity was considered to be merged to the main MS presenting the highest correlation value.

Statistical analysis

Descriptive statistics (i.e., mean, standard deviation (SD), proportion, range) synthesized the sociodemographics, anthropometrics, questionnaires and ultrasound-related outcomes. To describe AT-related changes at the Achilles tendon, ultrasound biomarkers (i.e., mean thickness, echogenicity, skewness and homogeneity) were compared between the asymptomatic and symptomatic LL using paired Student’s t-tests. To test the first hypothesis, GRF measures were compared using paired Student’s t-tests between the asymptomatic and symptomatic LL. Effect sizes were also computed using Hedges’ g [52] to determine the absolute magnitude of the estimates. Thus, an effect size > 0.2 was considered small, > 0.5 was considered medium, and > 0.8 was considered large [53]. Thereafter, the relative difference between the symptomatic and asymptomatic LL was computed (Eq. 5) and used to test the second hypothesis as described hereunder.

$$Difference\ \left(\%\right)=\frac{\left( Symptomatic- Asymptomatic\right)\ }{(Asymptomatic)}\times 100$$
(5)

To test the first hypothesis further, cosine similarities were compared between LL using paired Student’s t-tests. In addition, the weighting of each muscle within a MS was also compared using Student’s t-tests. To verify the extent to which the MS activation timing profile and EMG activation profiles were similar between LL at natural and fast speed, Pearson product–moment correlation coefficients (r) were calculated. Strength of correlation coefficients were considered negligible between 0.00 and 0.30, low between 0.30 and 0.50, moderate between 0.50 and 0.70, high between 0.70 and 0.90, and very high between 0.90 and 1.00 [54]. To test the secondary hypothesis, the difference in GRF measures between walking speed (i.e., Fast – Natural) as a percentage of BW was compared using paired t-tests between LL. To assess the change in MS attributes when going from natural to fast speed, Pearson product-moment correlation coefficients (r) of MS activation timing profiles between natural and fast speeds for both LL were calculated. Statistical analyses were carried out with SPSS v25 software and the statistical significance threshold was set at 0.05.

Results

Characteristics of participants

A summary of the characteristics of the participants and scores achieved on the VISA-A and LEFS questionnaires are presented in Table 1.

Table 1 Mean (standard deviation) characteristics of participants

Ultrasound biomarkers

A summary of Achilles tendon ultrasound biomarkers is presented in Table 2. The mean thickness revealed a significant (p < 0.001) and large between-side difference (g = 1.21), reaching + 29.7% for the symptomatic tendon when compared to the asymptomatic one. The echogenicity revealed a significant (p < 0.001) and large between-side difference (g = − 0.81), reaching − 13.9% for the symptomatic tendon when compared to the asymptomatic one. Skewness revealed a significant (p = 0.012) and medium between-side difference (g = 0.67), reaching a difference of + 85.5% for the symptomatic tendon compared to the asymptomatic one. The mean global homogeneity and perpendicular homogeneity (90°) respectively revealed a significant (p < 0.001 and p < 0.001) and only small between-side differences (g = 0.23 and 0.28), reaching + 2.2% to 2.9% for the symptomatic tendon compared to the asymptomatic one.

Table 2 Mean (standard deviation) of ultrasound variables in the longitudinal plane

Ground reaction forces

During walking at natural speed (1.3 m/s), the mean stride length was 1.39 ± 0.08 m and cadence was 56.1 ± 3.1 strides/minute for both the symptomatic and asymptomatic sides. During walking at fast speed (1.6 m/s), the mean stride length was 1.59 ± 0.09 m and cadence was 60.5 ± 3.5 strides/minute for both the symptomatic and asymptomatic sides. A summary of the mean and standard deviation (SD) of GRFML, GRFPA and GRFV during the gait cycle and the selected peak GRF for the asymptomatic and symptomatic LL during natural and fast speed walking are presented in Table 3 and Fig. 3, respectively. Among the peak GRF, only the medial propulsive force during natural and fast walking speeds respectively revealed a significant (p < 0.001 and p = 0.022) and small between-side difference (g = 0.23 and g = 0.23), reaching − 5.46% and − 5.49% difference for the symptomatic LL compared to the asymptomatic one. Peak GRF differences between natural and fast speed were similar for the symptomatic tendon versus the asymptomatic one.

Table 3 Mean (standard deviation) kinetic variables
Fig. 3
figure 3

Group average (full line) and average + 1SD (dotted lines) of GRF in the a medio-lateral, b postero-anterior and c vertical directions for the symptomatic (Sympt) and Asymptomatic (Asympt) sides at natural and fast walking speeds

Number of muscle synergies and cosine similarity

Two to five MS were extracted to characterize gait using NNMF for all eight muscles for each LL (Table 4 & Fig. 4). Overall, four MS with a specific set of predominantly activated muscles (Fig. 1) were extracted among the majority of participants during natural (71% and 61%) and fast (54% and 50%) walking speeds respectively, for the symptomatic and asymptomatic LL. Specifically, MS1 was identified in most participants on the symptomatic LL at natural speed (93%) but in fewer participants at fast speed (86%). MS1 was also identified in most participants on the asymptomatic LL at natural speed (79%) but in fewer participants at fast speed (75%). MS2 was identified in all participants (100%) across all conditions. MS3 was identified in most participants on the symptomatic LL at natural speed (82%) but in fewer participants at fast speed (75%). MS3 was also identified in most participants on the asymptomatic LL at natural speed (93%) but in fewer participants at fast speed (75%). MS4 was identified in all participants on the symptomatic LL at natural speed (100%) but in fewer participants at fast speed (93%). MS4 was also identified in most participants on the asymptomatic LL at natural speed (89%) and at fast speed (89%). Hence, the number of merged synergies was higher with the asymptomatic LL compared to the symptomatic LL, with this phenomenon being accentuated further by an increase in walking speed. Additional analysis revealed that MS1 and MS3 merged more often into other synergies than MS2 and MS4. This is also confirmed by the cosine similarity values in which MS1 and MS3 have generally lower values than MS2 and MS4 (Fig. 4b). Also, the cosine similarity values were comparable (p > 0.05) between LL during natural and fast walking speeds for each MS (Fig. 4b).

Table 4 Muscle synergies detected and merged among walking conditions
Fig. 4
figure 4

a Number of synergies is determined by the VAF (variance accounted for) criterion for each walking condition. b Cosine similarity values (r) for the weightings of each muscle synergy

Muscle synergies activation timing profile

Activation timing profiles, demonstrating how activation of each MS (i.e., MS1 to MS4) varies over the gait cycle, are shown in Fig. 5a. For each of the four main MS, the global VAF and the Pearson coefficient of correlation (r) between the symptomatic and asymptomatic LL for the natural and fast speed were marked on each graph. In general, the activation timing profiles for all conditions were remarkably similar for all four MS. MS1 confirms that the activation of the vastus medialis and gluteus medius predominantly contributes to the synergistic pattern observed between 0% and 20% of the gait cycle. MS2 confirms the activation of the soleus and medial gastrocnemius between 30% and 50% of the gait cycle. MS3 confirms the activation of the tibialis anterior, vastus medialis and rectus femoris with two peaks of activity occurring at the beginning and at approximately 60% to 75% of the gait cycle. MS4 confirms the activation of the tibialis anterior, semitendinosus and biceps femoris between 85% and 100% of the gait cycle. The global VAS was equal to 0.978, 0.966, 0.991 and 0.973 for MS1, MS2, MS3 and MS4, respectively. Correlation between the activation profiles of the asymptomatic and symptomatic LL were very high (r > 0.973) in all comparisons.

Fig. 5
figure 5

Group average (n = 28) for each of the four muscle synergies at natural and fast walking speeds. a Activation timing profiles for each synergy over the gait cycle with corresponding global VAF (variance accounted for). Group average (full line) and average + 1SD (dotted lines). b Muscle weightings average and SD of the four synergies. RF = rectus femoris, VM = vastus medialis, TA = tibialis anterior, MG = medial gastrocnemius, SO = soleus, ST = semitendinosus, BF = biceps femoris, GM = gluteus medius

Muscle weightings

For each MS selected, muscle weightings were calculated to indicate the strength of representation of each muscle among each LL and condition (Fig. 5b). When comparing the muscle weighting between the LL, only two comparisons were found to be significantly different. For MS1, the muscle weighting of the biceps femoris was found to be significantly lower for the symptomatic LL (p = 0.008) compared to the asymptomatic LL. For MS4, the muscle weighting of the soleus was found to be lower for the symptomatic side (p = 0.03) when compared to the asymptomatic LL.

Individualized EMG activation profiles

Individualized experimental EMG activation profiles were illustrated to demonstrate the muscle activation pattern of each muscle over the gait cycle (Fig. 6). For each individualized EMG profile, the correlation between the symptomatic and asymptomatic LL for natural and fast speeds were marked on the graph. The correlation between the asymptomatic and symptomatic LL was found to be very high (r > 0.970) for all comparisons.

Fig. 6
figure 6

Individual EMG activation profiles over the gait cycle. Each muscle activity was normalized by maximum activation across each walking speed. Group average (full line) and average + 1SD (dotted lines)

Discussion

The effects of unilateral AT pain and Achilles tendon integrity changes on GRF asymmetries and MS attributes during walking were investigated. Overall, the results revealed that the presence of AT had limited effects on bilateral GRF asymmetries, muscle weightings and MS activation profiles between the symptomatic and asymptomatic LL during walking at natural and fast speeds. In most participants, four MS (i.e., MS1 to MS4) were sufficient to explain the majority of the VAF (i.e., > 90% of the VAF for each muscle and > 80% of the global VAF). For all conditions, MS demonstrated relatively similar temporal activation profiles and muscle weightings to those reported in previous studies among healthy adults [25, 34, 55], with very few key differences. Increasing the walking speed from a natural pace (1.3 m/s) to a fast pace (1.6 m/s) increased some peak GRF values and increased the number of merged MS for both LL, but neither significantly altered MS temporal activation profiles nor muscle weighting between the symptomatic and asymptomatic LL.

Limited effects on bilateral LL GRF asymmetries and MS attributes

The primary objective of this study was to compare GRF symmetries and MS attributes between the asymptomatic and symptomatic LL during walking in adults with unilateral symptomatic AT. The results partially supported the first hypothesis, which stated that peak GRF would be reduced at the symptomatic LL compared to the asymptomatic LL. First, only the two GRFML peaks were significantly different between the symptomatic and asymptomatic LL at natural and fast speeds. An increase in the GRFML peak during midstance on the symptomatic side may be explained, in part, by an increase in vastus medialis and medial gastrocnemius muscle activity to promote medial braking force [56]. A decrease in peak GRFML during terminal stance on the symptomatic side may also be explained, in part, by a reduction in hip adductor and medial gastrocnemius contractions to lower medial propulsive forces [56]. Because the absolute medial propulsive force differences between LL were ultimately very small (< 0.5% BW), it is also plausible that these potential muscle contraction differences during walking could not be detected or that they predominantly related to measurement errors associated with the treadmill or EMG recording system or both. Peak GRF in the postero-anterior and vertical directions were similar between LL at natural and fast walking speeds, which can be explained by the similar muscle activation patterns observed among all eight muscles analyzed [57].

Contrary to our hypothesis for MS, no unilateral change in the motor recruitment strategy of the hip or knee muscle stabilizers was observed during the support phase (i.e., MS1) for the symptomatic LL. Likewise, no motor recruitment strategy difference was observed at the ankle during the pushoff phase (i.e., MS2) for the symptomatic LL. Such an adaptation was anticipated to reduce the tensile force transiting through the symptomatic Achilles tendon and ultimately decrease the peak GRFV. Still contrary to our hypothesis, the number of merged synergies was slightly higher at the asymptomatic LL than the symptomatic one. Considering the bilateral nature of walking and the required adjustments of the CNS, these results may highlight potential motor control adaptations of the asymptomatic LL over time (i.e., chronic AT) to preserve symmetry during walking. As such, it was previously suggested that following a unilateral chronic musculoskeletal injury, neuroplastic adaptations and their effects on the CNS can explain changes in sensory and motor cortical representation, resulting in bilateral perceptual changes of body image and motor control [18]. Similarly, a relatively recent systematic review confirms the presence of sensory and motor alterations on the non-injured side of adults with unilateral tendon pain and related disability [58]. These central adaptations add to the peripheral adaptations highlighted in previous studies that have shown histological structure changes in the asymptomatic contralateral tendon [16, 59, 60]. Overall, the results of the present study align with evidence supporting the bilateral nature of tendinopathy and the involvement of CNS mechanisms.

MS attributes difference with respect to healthy adults

During walking at natural and fast speeds, the twenty-eight participants in our study had relatively similar individualized EMG activation profiles and MS attributes (i.e., activation timing profiles and muscle weightings) (MS1 to MS4) compared to those previously reported in the literature [26, 32, 34, 61,62,63]. However, some MS attribute differences in comparison with healthy adults warrant attention. For MS1, the gluteus medius and vastus medialis weightings were lower in this study. For MS4, the biceps femoris and semitendinosus weightings were lower, whereas the tibialis anterior weighting was greater. These differences could be explained bilaterally by a decrease in hip and knee stability during stance and an increase in tibialis anterior co-contractions during terminal swing, which has been previously described among adults with unilateral acute pain in the calf muscles [30]. Such changes could be interpreted as an adaptation of the muscle activation pattern in response to the injured tendon or aetiological causative factor of AT.

Increasing the walking speed did not alter GRF asymmetries and MS attributes differences between LL

The secondary objective was to verify if increasing the walking speed revealed additional differences in terms of GRF asymmetries and MS attributes between LL. The results only partially supported the hypothesis. Thus, increasing the walking speed from 1.3 m/s to 1.6 m/s resulted in higher peak GRF values The same phenomenon has been observed in healthy adults in previous studies [49, 56, 64, 65]. However, increased walking speed had limited effects on GRF and MS differences between LL. Also, the number of synergies merged to a greater extent at fast speed compared to natural speed, yet merged equally on both LL. An increase in some peak GRF values during fast speed is more likely associated with adjustments in the mechanical output of the MS rather than differences in the profile and activation timings.

Translating evidence into tendon rehabilitation clinical practice

Current AT rehabilitation in clinical practice focuses predominantly on tendon loading, as it stimulates the physiological adaptation of the muscle-tendon complex [4, 7]. In clinical practice, eccentric plantarflexor strengthening exercises and heavy-slow plantarflexor resistance training has been shown to reduce pain and improve function among adults with AT [66,67,68]. Based on the high recurrence rate of AT and potential persistent motor changes following AT, current clinical practice may fail to adequately address contributing peripheral or central factors that may impact motor control. On one hand, CNS cortical representation, interhemispheric inhibition and motor cortex excitability may warrant additional attention during rehabilitation [69, 70]. Accordingly, tendon loading exercises that would be externally paced could be an interesting complement to rehabilitation protocols to best modulate tendon pain and CNS motor control [71]. On the other hand, focusing solely on the symptomatic AT during rehabilitation may attenuate the potential beneficial effects given potential tendon changes on the asymptomatic tendon. Rehabilitation protocols may therefore also need to involve the asymptomatic LL to prompt beneficial tendon integrity and CNS motor control adaptations over the long term [4, 69].

Limitations

Some limitations associated with the present study require discussion. The relatively modest sample of participants (n = 28) with AT and the heterogeneity of the participants in terms of symptom duration and pain intensity may have reduced the statistical power and limited subgroup analysis. Also, the activity-specific muscle EMG amplitude normalization approach used may have limited the ability to compare peak EMG activation profile differences between conditions. As previously mentioned, the fact that the asymptomatic tendon was used as the comparator to assess the impact of AT may also have attenuated the magnitude of MS changes since both LL may have been affected in adults with unilateral AT [12, 58, 59]. Investigating adults with chronic AT and comparing them to healthy counterparts may have provided additional insights. As a mitigation strategy, supplementary materials present a cosine similarity analysis of the muscle weightings between our data and previously published normative data of healthy adults [34]. Caution remains advised if attempting to generalize the current findings with other tendons (e.g., patellar tendon, supraspinatus tendon), other walking speeds (e.g., slow speed, self-selected speed) [72], or other walking conditions (e.g., overground walking) [73]. Personal characteristics also known to influence walking (e.g., gender, height, leg length, body weight) would also need additional attention if within-subject comparisons were to be performed in the future. Finally, higher impact activity soliciting a greater amount of force or generating a higher pain level compared to walking, such as running and jumping, might have detected greater differences in peak GRF and MS attributes given the higher loading at the Achilles tendon [4, 11, 74].

Conclusion

The presence of AT had limited effects on peak GRF and MS number, composition and temporal profiles between the symptomatic and asymptomatic LL for level treadmill walking at natural and fast speeds. Peripheral (i.e., changes in the asymptomatic tendon) or central adaptations (i.e., corticospinal neuroplastic changes) related to chronic unilateral AT may explain the preserved quasi-symmetric LL motor control during natural and fast walking among adults with chronic unilateral AT. Increasing LL muscular demand further (e.g., running, jumping) may have altered the ability of adults with chronic AT to modulate excitatory and inhibitory control of their LL muscles. The paradigm shift in current tendon-focused rehabilitation strategies deserves continued attention to best address corticospinal neuroplasticity adaptations.

Availability of data and materials

The datasets used and analyzed in this study are available from the corresponding author upon reasonable request.

Abbreviations

AT:

Achilles tendinopathy

GRF:

Ground reaction forces

MS:

Muscle synergies

EMG:

Electromyography

LL:

Lower limb

CNS:

Central nervous system

NNMF:

Non-negative matrix factorization

VISA-A:

Victorian Institute of Sport Assessment-Achilles

LEFS:

Lower Extremity Functional Scale

SD:

Standard deviation

LTF:

Lateral thrust force

MIF:

Medial impact force

HBF:

Horizontal braking force

HPF:

Horizontal propulsive force

VIF:

Vertical impact force

MVP:

Minimal vertical peak

VPF:

Vertical propulsive force

VAF:

Variance accounted for

RF:

Rectus femoris

VM:

Vastus medialis

TA:

Tibialis anterior

MG:

Medial gastrocnemius

SO:

Soleus

ST:

Semitendinosus

BF:

Biceps femoris

GM:

Gluteus medius

References

  1. O'Brien M. Functional anatomy and physiology of tendons. Clin Sports Med. 1992;11(3):505–20.

    Article  CAS  Google Scholar 

  2. Winnicki K, Ochała-Kłos A, Rutowicz B, Pękala PA, Tomaszewski KA. Functional anatomy, histology and biomechanics of the human Achilles tendon — a comprehensive review. Ann Anatomy Anatomischer Anzeiger. 2020;229:151461. https://doi.org/10.1016/j.aanat.2020.151461.

    Article  PubMed  Google Scholar 

  3. Wilson A, Lichtwark G. The anatomical arrangement of muscle and tendon enhances limb versatility and locomotor performance. Philos Trans R Soc Lond B Biol Sci. 2011;366(1570):1540–53.

    Article  PubMed Central  Google Scholar 

  4. Martin RL, Chimenti R, Cuddeford T, Houck J, Matheson JW, McDonough CM, et al. Achilles pain, stiffness, and muscle power deficits: Midportion Achilles tendinopathy revision 2018. J Orthop Sports Phys Ther. 2018;48(5):A1–A38. https://doi.org/10.2519/jospt.2018.0302.

    Article  PubMed  Google Scholar 

  5. Leadbetter WB, Buckwalter JA, Gordon SL. Sports-induced inflammation: clinical and basic science concepts. Amer Acad Orthop. 1990. p. 189–204.

  6. Cook J, Purdam C. Is compressive load a factor in the development of tendinopathy? Br J Sports Med. 2012;46(3):163–8.

    Article  CAS  Google Scholar 

  7. Scott A, Backman LJ, Speed C. Tendinopathy: update on pathophysiology. J Orthop Sports Phys Ther. 2015;45(11):833–41. https://doi.org/10.2519/jospt.2015.5884.

    Article  PubMed  Google Scholar 

  8. Järvinen M, Jozsa L, Kannus P, Järvinen T, Kvist M, Leadbetter W. Histopathological findings in chronic tendon disorders. Scand J Med Sci Sports. 1997;7(2):86–95. https://doi.org/10.1111/j.1600-0838.1997.tb00124.x.

    Article  PubMed  Google Scholar 

  9. Maffulli N, Sharma P, Luscombe KL. Achilles tendinopathy: aetiology and management. J R Soc Med. 2004;97(10):472–6.

    Article  PubMed Central  Google Scholar 

  10. Ogbonmwan I, Kumar BD, Paton B. New lower-limb gait biomechanical characteristics in individuals with Achilles tendinopathy: A systematic review update. Gait Posture. 2018;62:146-56. https://doi.org/10.1016/j.gaitpost.2018.03.010.

  11. Sancho I, Malliaras P, Barton C, Willy RW, Morrissey D. Biomechanical alterations in individuals with Achilles tendinopathy during running and hopping: a systematic review with meta-analysis. Gait Posture. 2019;73:189–201. https://doi.org/10.1016/j.gaitpost.2019.07.121.

    Article  PubMed  Google Scholar 

  12. Baur H, Divert C, Hirschmüller A, Müller S, Belli A, Mayer F. Analysis of gait differences in healthy runners and runners with chronic Achilles tendon complaints. Isokinet Exerc Sci. 2004;12(2):111–6.

    Article  Google Scholar 

  13. Baur H, Müller S, Hirschmüller A, Cassel M, Weber J, Mayer F. Comparison in lower leg neuromuscular activity between runners with unilateral mid-portion Achilles tendinopathy and healthy individuals. J Electromyogr Kinesiol. 2011;21(3):499–505.

    Article  Google Scholar 

  14. Wyndow N, Cowan S, Wrigley T, Crossley K. Triceps surae activation is altered in male runners with Achilles tendinopathy. J Electromyogr Kinesiol. 2013;23(1):166–72.

    Article  CAS  Google Scholar 

  15. Franettovich Smith MM, Honeywill C, Wyndow N, Crossley KM, Creaby MW. Neuromotor control of gluteal muscles in runners with achilles tendinopathy. Med Sci Sports Exerc. 2014;46(3):594–9. https://doi.org/10.1249/mss.0000000000000133.

    Article  PubMed  Google Scholar 

  16. Azevedo LB, Lambert MI, Vaughan CL, O'Connor CM, Schwellnus MP. Biomechanical variables associated with achilles tendinopathy in runners. Br J Sports Med. 2008. https://doi.org/10.1136/bjsm.2008.053421.

  17. Hodges PW, Tucker K. Moving differently in pain: a new theory to explain the adaptation to pain. Pain. 2011;152(3):S90–S8.

    Article  Google Scholar 

  18. Pelletier R, Higgins J, Bourbonnais D. Is neuroplasticity in the central nervous system the missing link to our understanding of chronic musculoskeletal disorders? BMC Musculoskelet Disord. 2015;16(1):25. https://doi.org/10.1186/s12891-015-0480-y.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Tompra N, van Dieën JH, Coppieters MW. Central pain processing is altered in people with Achilles tendinopathy. Br J Sports Med. 2016;50(16):1004–7. https://doi.org/10.1136/bjsports-2015-095476.

    Article  PubMed  Google Scholar 

  20. Bernstein NA. The co-ordination and regulation of movements. London: Pergamon; 1967.

    Google Scholar 

  21. Cheung VCK, Seki K. Approaches to revealing the neural basis of muscle synergies: a review and a critique. J Neurophysiol. 2021;125(5):1580–97. https://doi.org/10.1152/jn.00625.2019.

    Article  PubMed  Google Scholar 

  22. Devarajan K, Cheung VC. On nonnegative matrix factorization algorithms for signal-dependent noise with application to electromyography data. Neural Comput. 2014;26(6):1128–68.

    Article  PubMed Central  Google Scholar 

  23. Lee DD, Seung HS. Algorithms for non-negative matrix factorization. In Advances in Neural Information Processing Systems 13 - Proceedings of the 2000 Conference, NIPS 2000. Neural information processing systems foundation. 2001.

  24. Neptune RR, Clark DJ, Kautz SA. Modular control of human walking: a simulation study. J Biomech. 2009;42(9):1282–7.

    Article  PubMed Central  Google Scholar 

  25. Lacquaniti F, Ivanenko YP, Zago M. Patterned control of human locomotion. J Physiol. 2012;590(10):2189–99.

    Article  CAS  PubMed Central  Google Scholar 

  26. Escalona MJ, Bourbonnais D, Goyette M, Le Flem D, Duclos C, Gagnon DH. Effects of varying Overground walking speeds on lower-extremity muscle synergies in healthy individuals. Mot Control. 2021;25(2):234–51. https://doi.org/10.1123/mc.2020-0008.

    Article  Google Scholar 

  27. Bowden MG, Clark DJ, Kautz SA. Evaluation of abnormal synergy patterns poststroke: relationship of the Fugl-Meyer assessment to hemiparetic locomotion. Neurorehabil Neural Repair. 2010;24(4):328–37.

    Article  Google Scholar 

  28. Cheung VC, Turolla A, Agostini M, Silvoni S, Bennis C, Kasi P, et al. Muscle synergy patterns as physiological markers of motor cortical damage. Proc Natl Acad Sci. 2012;109(36):14652–6.

    Article  CAS  PubMed Central  Google Scholar 

  29. Pérez-Nombela S, Barroso F, Torricelli D, de Los R-GA, Del-Ama A, Gómez-Soriano J, et al. Modular control of gait after incomplete spinal cord injury: differences between sides. Spinal Cord. 2017;55(1):79–86.

    Article  Google Scholar 

  30. van den Hoorn W, Hodges PW, van Dieën JH, Hug F. Effect of acute noxious stimulation to the leg or back on muscle synergies during walking. J Neurophysiol. 2015;113(1):244–54.

    Article  Google Scholar 

  31. Henriksen M, Aaboe J, Graven-Nielsen T, Bliddal H, Langberg H. Motor responses to experimental Achilles tendon pain. Br J Sports Med. 2011;45(5):393–8.

    Article  Google Scholar 

  32. Serrancolí G, Monllau JC, Font-Llagunes JM. Analysis of muscle synergies and activation–deactivation patterns in subjects with anterior cruciate ligament deficiency during walking. Clin Biomech. 2016;31:65–73.

    Article  Google Scholar 

  33. Allison K, Salomoni SE, Bennell KL, Wrigley TV, Hug F, Vicenzino B, et al. Hip abductor muscle activity during walking in individuals with gluteal tendinopathy. Scand J Med Sci Sports. 2018;28(2):686–95.

    Article  CAS  Google Scholar 

  34. Clark DJ, Ting LH, Zajac FE, Neptune RR, Kautz SA. Merging of healthy motor modules predicts reduced locomotor performance and muscle coordination complexity post-stroke. J Neurophysiol. 2010;103(2):844–57. https://doi.org/10.1152/jn.00825.2009.

    Article  PubMed  Google Scholar 

  35. Vestergrd Iversen J, Bartels EM, Langberg H. The Victorian Institute of Sports Assessment-Achilles questionnaire (VISA-A)-reliable tool for measuring achilles tendinopathy. Int J Sports Phys Ther. 2012;7(1):76–84.

    PubMed  Google Scholar 

  36. Robinson J, Cook JL, Purdam C, Visentini P, Ross J, Maffulli N, et al. The VISA-A questionnaire: a valid and reliable index of the clinical severity of Achilles tendinopathy. Br J Sports Med. 2001;35(5):335–41. https://doi.org/10.1136/bjsm.35.5.335.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kaux J-F, Delvaux F, Oppong-Kyei J, Dardenne N, Beaudart C, Buckinx F, et al. Validity and reliability of the French translation of the VISA-A questionnaire for Achilles tendinopathy. Disabil Rehabil. 2016:1–7. https://doi.org/10.3109/09638288.2016.1138553.

  38. McCormack J, Underwood F, Slaven E, Cappaert T. The minimum clinically important difference on the VISA-A and LEFS for patients with insertional achilles tendinopathy. Int J Sports Phys Ther. 2015;10(5):639.

    PubMed  PubMed Central  Google Scholar 

  39. Binkley JM, Stratford PW, Lott SA, Riddle DL. The lower extremity functional scale (LEFS): scale development, measurement properties, and clinical application. Phys Ther. 1999;79(4):371–83. https://doi.org/10.1093/ptj/79.4.371.

    Article  CAS  PubMed  Google Scholar 

  40. Nadeau M-J, Desrochers A, Lamontagne M, Larivière C, Gagnon DH. Quantitative ultrasound imaging of Achilles tendon integrity in symptomatic and asymptomatic individuals: reliability and minimal detectable change. J Foot Ankle Res. 2016;9(1):30. https://doi.org/10.1186/s13047-016-0164-3.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lalumiere M, Perrino S, Nadeau M-J, Larivière C, Lamontagne M, Desmeules F, et al. To what extent do musculoskeletal ultrasound biomarkers relate to pain, flexibility, strength, and function in individuals with chronic symptomatic Achilles tendinopathy? Front Rehabil Sci. 2021;2(19). https://doi.org/10.3389/fresc.2021.726313.

  42. Lalumiere M, Larivière C, Nadeau MJ, Paquette P, Lamontagne M, Desmeules F, et al. Proposing a minimal data set of musculoskeletal ultrasound imaging biomarkers to inform clinical practice: an analysis founded on the Achilles tendon. Ultrasound Med Biol. 2020;46(9):2222–35. https://doi.org/10.1016/j.ultrasmedbio.2020.04.024.

    Article  PubMed  Google Scholar 

  43. Roberts M, Mongeon D, Prince F. Biomechanical parameters for gait analysis: a systematic review of healthy human gait. Phys Ther Rehabil. 2017;4(1):6. https://doi.org/10.7243/2055-2386-4-6.

  44. Perotto AO. Anatomical guide for the electromyographer: the limbs and trunk: Charles C Thomas Publisher; 2011.

    Google Scholar 

  45. Hislop H, Avers D, Brown M. Daniels and Worthingham’s muscle testing-E-book: techniques of manual examination and performance testing: Elsevier Health Sciences; 2013.

    Google Scholar 

  46. Hawker G, Mian S, Kendzerska T, French M. Measures of adult pain: visual analog scale for pain (VAS pain), numeric rating scale for pain (NRS pain), McGill pain questionnaire (MPQ), short-form McGill pain questionnaire (SF-MPQ), chronic pain grade scale (CPGS), short Form-36 bodily pain scale (SF-36 BPS), and measure of intermittent and constant osteoarthritis pain (ICOAP). Arthritis Care Res. 2011;63(S11):S240–S52. https://doi.org/10.1002/acr.20543.

    Article  Google Scholar 

  47. Escalona MJ, Bourbonnais D, Goyette M, Duclos C, Gagnon DH. Wearable exoskeleton control modes selected during overground walking affect muscle synergies in adults with a chronic incomplete spinal cord injury. Spinal Cord Ser Cases. 2020;6(1):26. https://doi.org/10.1038/s41394-020-0269-6.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Torres-Oviedo G, Macpherson JM, Ting LH. Muscle synergy organization is robust across a variety of postural perturbations. J Neurophysiol. 2006;96(3):1530–46. https://doi.org/10.1152/jn.00810.2005.

    Article  PubMed  Google Scholar 

  49. Kibushi B, Hagio S, Moritani T, Kouzaki M. Speed-Dependent Modulation of Muscle Activity Based on Muscle Synergies during Treadmill Walking. Front Hum Neurosci. 2018;12:4. https://doi.org/10.3389/fnhum.2018.00004.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Booth ATC, van der Krogt MM, Harlaar J, Dominici N, Buizer AI. Muscle synergies in response to biofeedback-driven gait adaptations in children with cerebral palsy. Front Physiol. 2019;10(1208). https://doi.org/10.3389/fphys.2019.01208.

  51. Oliveira AS, Gizzi L, Farina D, Kersting UG. Motor modules of human locomotion: influence of EMG averaging, concatenation, and number of step cycles. Front Hum Neurosci. 2014;8:335.

    Article  PubMed Central  Google Scholar 

  52. Hedges LV. Distribution theory for Glass's estimator of effect size and related estimators. J Educ Stat. 1981;6(2):107–28. https://doi.org/10.3102/10769986006002107.

    Article  Google Scholar 

  53. Sawilowsky SS. New effect size rules of thumb. J Mod Appl Stat Methods. 2009;8(2):26. https://doi.org/10.22237/jmasm/1257035100.

    Article  Google Scholar 

  54. Mukaka MM. A guide to appropriate use of correlation coefficient in medical research. Malawi Med J. 2012;24(3):69–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Routson RL, Clark DJ, Bowden MG, Kautz SA, Neptune RR. The influence of locomotor rehabilitation on module quality and post-stroke hemiparetic walking performance. Gait Posture. 2013;38(3):511–7. https://doi.org/10.1016/j.gaitpost.2013.01.020.

    Article  PubMed  PubMed Central  Google Scholar 

  56. John CT, Seth A, Schwartz MH, Delp SL. Contributions of muscles to mediolateral ground reaction force over a range of walking speeds. J Biomech. 2012;45(14):2438–43. https://doi.org/10.1016/j.jbiomech.2012.06.037.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Liu MQ, Anderson FC, Schwartz MH, Delp SL. Muscle contributions to support and progression over a range of walking speeds. J Biomech. 2008;41(15):3243–52. https://doi.org/10.1016/j.jbiomech.2008.07.031.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Heales LJ, Lim ECW, Hodges PW, Vicenzino B. Sensory and motor deficits exist on the non-injured side of patients with unilateral tendon pain and disability—implications for central nervous system involvement: a systematic review with meta-analysis. Br J Sports Med. 2014;48(19):1400–6. https://doi.org/10.1136/bjsports-2013-092535.

    Article  CAS  PubMed  Google Scholar 

  59. Docking SI, Rosengarten SD, Daffy J, Cook J. Structural integrity is decreased in both Achilles tendons in people with unilateral Achilles tendinopathy. J Sci Med Sport. 2015;18(4):383–7. https://doi.org/10.1016/j.jsams.2014.06.004.

    Article  PubMed  Google Scholar 

  60. Andersson G, Forsgren S, Scott A, Gaida JE, Stjernfeldt JE, Lorentzon R, et al. Tenocyte hypercellularity and vascular proliferation in a rabbit model of tendinopathy: contralateral effects suggest the involvement of central neuronal mechanisms. Br J Sports Med. 2011;45(5):399–406. https://doi.org/10.1136/bjsm.2009.068122.

    Article  PubMed  Google Scholar 

  61. Winter DA, Yack HJ. EMG profiles during normal human walking: stride-to-stride and inter-subject variability. Electroencephalogr Clin Neurophysiol. 1987;67(5):402–11. https://doi.org/10.1016/0013-4694(87)90003-4.

    Article  CAS  PubMed  Google Scholar 

  62. Simonsen EB. Contributions to the understanding of gait control: University of Copenhagen; 2014.

    Google Scholar 

  63. Israely S, Leisman G, Carmeli E. Neuromuscular synergies in motor control in normal and poststroke individuals. Rev Neurosci. 2018;29(6):593–612. https://doi.org/10.1515/revneuro-2017-0058.

    Article  PubMed  Google Scholar 

  64. Burnett DR, Campbell-Kyureghyan NH, Cerrito PB, Quesada PM. Symmetry of ground reaction forces and muscle activity in asymptomatic subjects during walking, sit-to-stand, and stand-to-sit tasks. J Electromyogr Kinesiol. 2011;21(4):610–5. https://doi.org/10.1016/j.jelekin.2011.03.006.

    Article  PubMed  Google Scholar 

  65. Vaverka F, Elfmark M, Svoboda Z, Janura M. System of gait analysis based on ground reaction force assessment. Acta Gymnica. 2015;45(4):187–93.

    Article  Google Scholar 

  66. Murphy MC, Travers MJ, Chivers P, Debenham JR, Docking SI, Rio EK, et al. Efficacy of heavy eccentric calf training for treating mid-portion Achilles tendinopathy: a systematic review and meta-analysis. Br J Sports Med. 2019. https://doi.org/10.1136/bjsports-2018-099934.

  67. Wilson F, Walshe M, O'Dwyer T, Bennett K, Mockler D, Bleakley C. Exercise, orthoses and splinting for treating Achilles tendinopathy: a systematic review with meta-analysis. Br J Sports Med. 2018. https://doi.org/10.1136/bjsports-2017-098913.

  68. Van Der Vlist AC, Winters M, Weir A, Ardern CL, Welton NJ, Caldwell DM, et al. Which treatment is most effective for patients with Achilles tendinopathy? A living systematic review with network meta-analysis of 29 randomised controlled trials. Br J Sports Med. 2021;55(5):249–56.

    Article  Google Scholar 

  69. Paavola M, Kannus P, Paakkala T, Pasanen M, Järvinen M. Long-term prognosis of patients with Achilles tendinopathy an observational 8-year follow-up study. Am J Sports Med. 2000;28(5):634–42.

    Article  CAS  Google Scholar 

  70. Rio E, Kidgell D, Moseley GL, Gaida J, Docking S, Purdam C, et al. Tendon neuroplastic training: changing the way we think about tendon rehabilitation: a narrative review. Br J Sports Med. 2016;50(4):209–15. https://doi.org/10.1136/bjsports-2015-095215.

    Article  PubMed  Google Scholar 

  71. Ackerley SJ, Stinear CM, Byblow WD. Promoting use-dependent plasticity with externally-paced training. Clin Neurophysiol. 2011;122(12):2462–8. https://doi.org/10.1016/j.clinph.2011.05.011.

    Article  PubMed  Google Scholar 

  72. Astephen Wilson JL. Challenges in dealing with walking speed in knee osteoarthritis gait analyses. Clin Biomech (Bristol, Avon). 2012;27(3):210–2. https://doi.org/10.1016/j.clinbiomech.2011.09.009.

    Article  Google Scholar 

  73. Meyer C, Killeen T, Easthope CS, Curt A, Bolliger M, Linnebank M, et al. Familiarization with treadmill walking: How much is enough? Sci Rep. 2019;9(1):1–10.

    Article  Google Scholar 

  74. Gu Y, Li J, Lake M, Ren X, Zeng Y. The mechanical response of Achilles tendon during different kinds of sports. Commun Numer Methods Eng. 2008;24(12):2077–85.

    Article  Google Scholar 

Download references

Acknowledgments

The authors would like to acknowledge Philippe Gourdou for his assistance during data recording and processing.

Funding

This project was funded through a Clinical Research Partnership in Physiotherapy supported by the Ordre Professionnel de la Physiothérapie du Québec (OPPQ) and the Quebec Rehabilitation Research Network (REPAR). ML holds a doctoral training scholarship from the Fonds de la recherche en santé du Québec (FRSQ). FD and DHG hold a Research Career Salary Award from the FRSQ and DHG co-chairs the Initiative for the Development of New Technologies and Practices in Rehabilitation (INSPIRE) funded by the LRH Foundation. The equipment and material needed for work completed at the Pathokinesiology Laboratory were funded in part by the Canada Foundation for Innovation (CFI).

Author information

Authors and Affiliations

Authors

Contributions

ML and DHG contributed to the design of the project, development of the methodology, ethics approval process, collection of data, validation, analysis and interpretation, writing of the original draft of the manuscript, final review and editing of the manuscript, project administration, and acquisition of funding. SP and MG contributed to the acquisition of data, validation and analysis, the writing of the original draft of the manuscript, and the final review and editing of the manuscript. DB and FD contributed to the design of the project, development of the methodology, data analysis and interpretation, the final review and editing of the manuscript, and acquisition of funding. All authors contributed to the article and approved the submitted version.

Corresponding author

Correspondence to Dany H. Gagnon.

Ethics declarations

Ethics approval and consent to participate

All research on human subjects reported in this manuscript was prospectively approved by the Centre for Interdisciplinary Research in Rehabilitation of Greater Montreal (CRIR) Research Ethics Committees (certificate: CRIR-1224-0317). Participants reviewed and signed an informed consent form before entering the study. For all participant-related information and data (including ultrasound images), all direct personal identifiers were replaced with a project-specific code for the purposes of protecting their identity. However, the original identifiers were retained in such a way that they could be traced back to the source if needed only by the corresponding author.

Consent for publication

All subjects and authors have provided consent for this publication.

Competing interests

The authors have no conflicts of interest to report.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lalumiere, M., Bourbonnais, D., Goyette, M. et al. Unilateral symptomatic Achilles tendinopathy has limited effects on bilateral lower limb ground reaction force asymmetries and muscular synergy attributes when walking at natural and fast speeds. J Foot Ankle Res 15, 66 (2022). https://doi.org/10.1186/s13047-022-00570-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13047-022-00570-3

Keywords