- 1Faculty of Medical Engineering and Technomathematics, Aachen University of Applied Sciences, Aachen, Germany
- 2Institute of Movement and Neurosciences, German Sport University Cologne, Cologne, Germany
- 3Department of Physical Performance, Norwegian School of Sport Sciences, Oslo, Norway
- 4Orthopädie am Gürzenich, Cologne, Germany
- 5Institute of Biomechanics and Orthopedics, German Sport University Cologne, Cologne, Germany
Achilles tendon rupture (ATR) remains associated with functional limitations years after injury. Architectural remodeling of the gastrocnemius medialis (GM) muscle is typically observed in the affected leg and may compensate force deficits caused by a longer tendon. Yet patients seem to retain functional limitations during—low-force—walking gait. To explore the potential limits imposed by the remodeled GM muscle-tendon unit (MTU) on walking gait, we examined the contractile behavior of muscle fascicles during the stance phase. In a cross-sectional design, we studied nine former patients (males; age: 45 ± 9 years; height: 180 ± 7 cm; weight: 83 ± 6 kg) with a history of complete unilateral ATR, approximately 4 years post-surgery. Using ultrasonography, GM tendon morphology, muscle architecture at rest, and fascicular behavior were assessed during walking at 1.5 m⋅s–1 on a treadmill. Walking patterns were recorded with a motion capture system. The unaffected leg served as control. Lower limbs kinematics were largely similar between legs during walking. Typical features of ATR-related MTU remodeling were observed during the stance sub-phases corresponding to series elastic element (SEE) lengthening (energy storage) and SEE shortening (energy release), with shorter GM fascicles (36 and 36%, respectively) and greater pennation angles (8° and 12°, respectively). However, relative to the optimal fascicle length for force production, fascicles operated at comparable length in both legs. Similarly, when expressed relative to optimal fascicle length, fascicle contraction velocity was not different between sides, except at the time-point of peak series elastic element (SEE) length, where it was 39 ± 49% lower in the affected leg. Concomitantly, fascicles rotation during contraction was greater in the affected leg during the whole stance-phase, and architectural gear ratios (AGR) was larger during SEE lengthening. Under the present testing conditions, former ATR patients had recovered a relatively symmetrical walking gait pattern. Differences in seen AGR seem to accommodate the profound changes in MTU architecture, limiting the required fascicle shortening velocity. Overall, the contractile behavior of the GM fascicles does not restrict length- or velocity-dependent force potentials during this locomotor task.
Introduction
Recovery from Achilles tendon rupture (ATR) is a tedious process that results in permanent functional deficits in the majority of patients (Nilsson-Helander et al., 2010; Heikkinen et al., 2017). Such deficits are characterized by a weakness in end-range plantarflexion (Mullaney et al., 2006) and limited heel raise height (Silbernagel et al., 2012), and are often reflected by deficiencies during locomotor tasks (Willy et al., 2017). Gait deficiencies appear to be commensurate with movement velocity (Willy et al., 2017; Jandacka et al., 2018) and may therefore be related to the force-velocity relation of the affected muscles. Nonetheless, slow movements, — which require lower force levels — seem also impaired in most (Tengman and Riad, 2013; Agres et al., 2015; Willy et al., 2017; Speedtsberg et al., 2019) but not all investigated cases (Jandacka et al., 2017). The links between such deficits and the remodeling of the muscle-tendon unit (MTU) following ATR are, to date, poorly understood.
Long after recovery and regardless of treatment strategy, the MTU of ATR patients is characterized by a longer tendon (Silbernagel et al., 2012; Peng et al., 2019; Svensson et al., 2019). Recent studies suggest that the increased tendon stiffness (Agres et al., 2015; Stäudle et al., 2021) and shorter gastrocnemius medialis (GM) muscle fascicles (Baxter et al., 2018; Peng et al., 2019; Svensson et al., 2019) typically seen in ATR patients may compensate for their longer tendons, albeit incompletely (Stäudle et al., 2021). Using a musculoskeletal model to simulate maximum isometric contractions at various joint angles, we have shown that the shorter GM fascicles in the affected leg enables sarcomeres to operate close to their optimal length, but at the expense of a narrowed range for active force generation (Stäudle et al., 2021). The insights gained from these findings are, however, insufficient for predicting triceps surae mechanics in dynamic situations, where force-velocity conditions may set additional constraints.
As sarcomeres operate close to their optimal length during walking (Ishikawa et al., 2007), length-dependent deficits in muscle strength are expected to be rather small in ATR patients due to corresponding shorter fascicle lengths throughout the walking stance phase. Whereas, when considering the force-velocity relation of a muscle, shorter fascicles are expected to produce less force than longer ones at the same velocity because of their lower number of in-series sarcomeres (Hill, 1953; Bahler et al., 1968). This point may be critical in the case of ATR patients during walking, because of their shorter GM fascicle length and because of the force-limiting role of contractile velocity in the walking gait (Neptune and Sasaki, 2005; Farris and Sawicki, 2012). As the behavior of the remodeled muscle fascicles of ATR patients during walking gait has not yet been investigated and its impact on force generation is unclear, the purpose of this study was to investigate the hypothesis of altered contractile behavior of the GM causing a velocity-based deficit during walking in ATR patients. We expected the shorter GM fascicles of the affected leg to operate at a comparable length range but at a higher contractile velocity, relative to their optimal length, than the fascicles of the unaffected leg. Using combined ultrasonography and motion capture methods, we measured GM muscle mechanics in former ATR patients (more than 2 years post-surgery) walking at 1.5 m s–1.
Materials and Methods
Subjects
Male patients (20–60 years) were recruited for this study if they had suffered a complete ATR that had been treated surgically within 7 days after injury and were at least 2 years post-surgery. Subjects were excluded if they had a concomitant soleus muscle tear, sural nerve injury, or recurrent or bilateral ATR. The “Physical Activity Readiness Questionnaire” (Thomas et al., 1992) was used to exclude volunteers with cardiovascular or musculoskeletal disorders. The institutional review board of the German Sport University Cologne approved the study (approval number: 12/72), and all subjects provided written informed consent prior to voluntary participation.
This study is part of a comprehensive investigation on ATR patients’ functional deficits. For this purpose, sample size calculations were based on ATR patients’ strength deficits, as described previously (Stäudle et al., 2021). The a priori power analysis suggested a minimum sample size of 10 subjects.
Furthermore, GM tendon length, fascicle length, pennation angle and muscle thickness with a resting muscle had already been included in the previous study (Stäudle et al., 2021).
Study Design and Experimental Protocol
A cross-sectional design was used for this study. Data acquisition took place according to a pseudo randomized order between the affected and unaffected leg, while the unaffected leg served as control for matched comparison. During the first of two testing sessions (Figure 1), GM tendon length, fascicle length, pennation angle and muscle thickness were examined using ultrasonography with the subjects lying prone, with a resting muscle and ankle and knee joint angles in anatomical position (0°), as described previously (Stäudle et al., 2021). During the second session, subjects were familiarized to the treadmill by walking for about 5 min at 1.5 m⋅s–1 (h/p cosmos pulsar 4.0, 2.2 kW, Traunstein, Germany) using their own running shoes. Kinematic and ultrasonographic data were thereafter collected from each leg during six consecutive stance phases.
Figure 1. Experimental design. (A) Setup to determine gastrocnemius medialis (GM) tendon length and muscle architecture, showing the scanning location S1 and S2 of the ultrasound probe, visualizing the calcaneal insertion and the myotendinous junction of the GM tendon (red dots), respectively, defining GM tendon resting length (lt,rest). Scanning location 3 (S3, GM mid-belly) visualizes GM muscle architecture in resting muscle (A) and during walking trials (B), detailed by the ultrasound image. The treadmill setup (B) includes the motion capture system and the high-speed cameras. t0: time-point of rupture, t1: time-point of surgery, t2: time-point of measurements, φf: pennation angle.
Motion Capture of the Legs
A motion capture system (eight infrared cameras) sampling at 100 Hz (Vicon, Vicon Motion Systems Ltd., Oxford, United Kingdom) was used to capture knee and ankle joint kinematics. Anatomical landmarks were labeled via reflective markers representing the subjects’ greater trochanter, lateral femur condyle, lateral malleolus, calcaneus, and second metatarsal head. A static reference was captured with the ankle and knee joints in anatomical position to define 0°. For gait event detection (touch-down, toe-off), two additional high-speed cameras (Basler, 100 Hz, Ahrensburg, Germany) were positioned anteriorly and posteriorly to the treadmill belt.
Measurements of Muscle Fascicle Behavior
During walking, B-mode ultrasound (Prosound α7, ALOKA, Tokyo, Japan) image sequences (73 Hz) were recorded using a t-shaped 6 cm linear array transducer (UST-5713T, 13 MHz) fixed in a custom-made cast to the mid-belly of the GM muscle via self-adhesive bandages. For time-synchronization, a rectangular voltage pulse was generated and sent to all data capturing devices.
Ultrasonography is a frequently used method to quantify muscle architecture under dynamic conditions (Cronin and Lichtwark, 2013) and shows good reliability within session for fascicle length and pennation angle in the present dataset (intraclass correlation coefficient: 0.99). In previous studies using similar methods, our group measured good inter-rater reliability with an intraclass correlation coefficient of 0.97 for pennation angle (Pohle-Fröhlich et al., 2020) and good reliability between days represented by a coefficient of multiple correlation of 0.93 for fascicle length and 0.87 for pennation angle (Aggeloussis et al., 2010).
A semi-automatic tracking algorithm (UltraTrack Software, version 4.2) was used to quantify GM muscle architecture. A dominant fascicle was drawn over a visible fascicle fragment and tracked across all frames. Additionally, superficial and deep aponeuroses were segmented and tracked. Fascicle length was defined as the distance between the insertions of the fascicles on the superficial and deep aponeuroses (Farris et al., 2016; Farris and Lichtwark, 2016).
In the rare exception where the transducer’s 6 cm width field of view failed to display the entire tracked fascicle (Figure 1), the missing portion was manually extrapolated, assuming that the fascicles and aponeurosis extend linearly. Linear extrapolation is associated with an error of less than 6% during maximal contractions of the GM muscle (Muramatsu et al., 2002). The pennation angle was defined as the angle between the muscle fascicle and the deep aponeurosis. Changes in series elastic element (SEE) length were estimated by subtracting muscle shortening amplitude from changes in MTU length. To this end, muscle shortening patterns were obtained from the length of the geometric projection of fascicles onto the axis of the deeper aponeurosis (Fukunaga et al., 2001). The length of the GM MTU was determined via a multiple linear regression equation using normative data based on joint angles and shank length (Hawkins and Hull, 1990).
Data Processing
A custom-made script (MATLAB R2020b, The MathWorks, Inc., Natick, MA, United States) was used to analyze the data. Fascicle length and pennation angle were smoothed with a 5th order Butterworth lowpass filter at a 10 Hz cut-off frequency. Data were then time-normalized by being resampled to 101 data points and then averaged per stance. Muscle fascicle velocities were calculated as the time derivative of the respective lengths using the central difference method (Robertson et al., 2013). Marker trajectories of the kinematic measurement were smoothed with a Woltring filter (Generalized Cross Validation, smoothing: 10) (Vicon Nexus 2.2.2, Vicon Motion Systems Ltd., Oxford, United Kingdom).
All outcome parameters were obtained at the time-point of peak SEE length and calculated as average values during stance sub-phases of SEE lengthening and shortening. SEE lengthening was defined as the duration between initial ground contact and peak SEE length, the latter presumably indicating maximal SEE loading. SEE shortening was defined as the remaining duration to toe-off.
Kinematic parameters were also examined at additional time-points, to obtain an exhaustive characterization of potential changes in gait pattern. Thus, ankle and knee joint angles were analyzed at touch-down and toe-off, peak ankle joint dorsiflexion, and the first knee joint flexion angle local maximum. The range of motion of each joint was defined as the difference between the angles’ minima and maxima.
Variables that characterize fascicle behavior were also studied, including operating fascicle length, fascicle velocity, and pennation angle. In addition, fascicle and muscle shortening amplitudes, changes in pennation angle and architectural gear ratio (AGR) were analyzed for each of the two sub-phases. A modified version of the AGR (Brainerd and Azizi, 2005) was calculated as the ratio between muscle shortening amplitude along the axis of the deeper aponeurosis [calculated as the product of fascicle length by the cosine of the pennation angle (Fukunaga et al., 2001)] and fascicle shortening amplitude during the SEE sub-phases of walking (Werkhausen et al., 2019).
In addition to absolute values, the average operating fascicle length was expressed relative to optimal fascicle length and termed normalized operating fascicle length. Average fascicle velocity was expressed relative to the velocity of one optimal fascicle length per second and termed normalized fascicle velocity. Optimal fascicle length was estimated from resting fascicle length and normative data of sarcomere length using the following equation:
where lf,rest is the fascicle length in resting condition (anatomically neutral position of the knee and ankle joint, 0°), ls,rest the sarcomere length at the identical joint angles reported by Sanchez et al. (2015) (3.09 μm), and ls,o the optimal sarcomere length of 2.725 μm defined by the mean value of the plateau region of the human sarcomere force-length relation (2.64–2.81 μm) (Herzog et al., 1990; Ward et al., 2009).
To illustrate the fraction of maximal GM force produced during walking, force potentials were estimated relative to operating length (length-dependent force potential) or velocity (velocity-dependent force potential). The length-dependent force potential was obtained using the default active-force-length curve of OpenSim based on quintic Bezier splines (Millard et al., 2013; Seth et al., 2018). As outlined in the OpenSim API guide (Millard, 2021), default parameters were chosen that the curve approximated the theoretical active-force-length curve of human sarcomeres (Nigg and Herzog, 1994) with the descending limb adapted from in vitro human fiber data (Gollapudi and Lin, 2009).
A dimensionless velocity-dependent force potential was obtained using the Hill-equation (Hill, 1938) for concentric contractions, relative values of force, velocity, and Hill’s constant a (Zajac, 1989) as follows:
vmax is the GM muscle’s maximum shortening velocity, estimated as multiples of optimal fascicle length per second and considering the GM fraction of fast twitch fibers (Winters and Stark, 1988) as follows:
where the fraction of fast twitch (FT) fibers for the GM muscle was assumed to be 49.2% (Johnson et al., 1973). The normalized Hill constant arel was calculated from 0.1 + 0.4FT (Winters and Stark, 1988; Bohm et al., 2019), yielding 0.297.
Statistical Analysis
Two tailed paired t-tests or Wilcoxon signed-rank tests were used to identify differences between the legs using Prism (version 7.04, GraphPad Software, Inc., San Diego, CA, United States). Prior to these tests, data were checked for normal distribution (Shapiro–Wilk normality test). Values are reported as means ± standard deviations (M ± SD) and the level of statistical significance was set to α ≤ 0.05. A statistical software package [G*Power version 3.1.9.6 (Faul et al., 2007)] was used to calculate absolute effect sizes, which were defined as small (0.2), moderate (0.5), or large (0.8), as appropriate (Cohen, 1988). An equivalent to effect sizes for signed rank tests was obtained by dividing the sum of the signed ranks by the total rank sum (matched pairs rank-biserial correlation) (Cureton, 1956; Kerby, 2014).
Results
Out of 59 former ATR patients who responded to advertisements at a medical center (n = 45) or in media announcements (n = 14), 14 were recruited for this study. Of the remainder, 23 patients did not meet inclusion criteria, eight had moved, and 14 declined to participate. Nine of the 14 recruited subjects performed all tests and were included in the present data. The data from the five remaining subjects could not be obtained because of incomplete ultrasound scans (n = 2) or schedule conflicts (n = 3).
The subjects (age: 44.7 ± 9.1 years; height: 180 ± 7 cm; weight: 82.8 ± 5.9 kg) had been operated with either a modified Bunnel, a Kessler end-to-end or a Dresdner Instrument (minimal invasive) 3.4 ± 1.8 days after rupture and, 4.4 ± 2.1 years prior to the study, on average. Structural differences in muscles and tendons were observed in the affected legs, with shorter GM muscle fascicles [30.8 ± 10.8%, t(8) = 7.28, p < 0.001, dz = 2.42], greater pennation angles [5 ± 3°, t(8) = 5.08, p = 0.001, dz = 1.69], reduced muscle thickness [9.8 ± 12.4%, t(8) = 2.40, p = 0.043, dz = 0.80], and longer GM tendons [13.9 ± 11.7%, t(8) = 3.73, p = 0.006, dz = 1.25] (Supplementary Table 1).
Kinematics
Figure 2 represents the knee and ankle joint angles and their corresponding angular velocities during stance. Except for a less plantarflexed averaged ankle joint angle during SEE shortening [3 ± 4°, t(8) = 2.34, p = 0.048, dz = 0.78] and a slower ankle angular velocity at peak SEE length [21.5 ± 25.2%, t(8) = 2.41, p = 0.043, dz = 0.80], the analyzed kinematic or spatio-temporal parameters did not show any significant differences (Figure 2, Table 1, and Supplementary Table 2). On average, the step frequency was 119 ± 5 steps⋅min–1.
Figure 2. Ankle (A) and knee (C) joint angles and corresponding angular velocities (B,D) during walking stance. The solid black and red lines, with their corresponding shaded areas [standard deviation (SD)], represent mean traces for the unaffected and affected leg, respectively. The dashed vertical lines with their corresponding shaded areas (SD) represent the mean time points of peak series elastic element (SEE) length (peak lSEE) separating the SEE lengthening and shortening sub-phases. Statistical testing was run on average values during the SEE lengthening and shortening sub-phases, and on single data points at peak SEE length. *p < 0.05, significant side-to-side difference, ns: non-significant.
Gastrocnemius Medialis Muscle Fascicle Length and Pennation Angle During Walking
The timing of maximum SEE elongation was used to separate the stance phase into SEE lengthening and shortening sub-phases. This time-point did not differ between the affected and unaffected leg [80.9 ± 2.2% vs. 82.3 ± 1.1%, t(8) = 1.76, p = 0.116, dz = 0.59].
During SEE lengthening and shortening sub-phases, the average operating fascicle length was shorter in the affected leg (36.2 ± 8.5% and 36.3 ± 10.0%, respectively), which was consistent with a shorter operating fascicle length at peak SEE length in the same leg (36.2 ± 10.6%) (Table 2 and Figure 3A). Shortening amplitudes did not differ significantly during the lengthening sub-phase, whereas smaller shortening amplitudes were observed during the shortening sub-phase in the affected leg (26.0 ± 32.2%) (Table 3). However, after normalization to optimal length, the average operating fascicle length did not differ between legs in either of the stance sub-phases or peak SEE length (Table 2 and Figure 4A), nor did the average length-dependent force potential (Table 2 and Figure 5A).
Figure 3. Gastrocnemius medialis operating fascicle length (A), pennation angle (B) and fascicle velocity (C) during walking stance. The solid black and red lines, with their corresponding shaded areas [standard deviation (SD)], represent mean traces for the unaffected and affected leg, respectively. The dashed vertical lines with their corresponding shaded areas (SD) represent the mean time points of peak series elastic element (SEE) length (peak lSEE) separating the SEE lengthening and shortening sub-phases. Statistical testing was run on average values during SEE lengthening and shortening sub-phases, and on single data points at peak SEE length. **p < 0.01, significant side-to-side difference; ns: non-significant.
Table 3. Changes in muscle architecture of the gastrocnemius medialis muscle during the series elastic element (SEE) lengthening and shortening sub-phases.
Figure 4. Differences between the affected and unaffected legs in gastrocnemius medialis normalized fascicle length (A) and normalized fascicle velocity (B). Statistical testing was run on average values during the series elastic element (SEE) lengthening and shortening sub-phases, and on single data points at peak SEE length (lSEE). The lower and upper parts of the box plots represent the first and third quartile, respectively. The length of the whisker delineates the minimum and maximum values. The horizontal line in the box represents the median of the sample; +, sample mean; *p < 0.05, significant side-to-side difference; ns: non-significant.
Figure 5. Differences between the affected and unaffected legs in gastrocnemius medialis length-dependent (A) and velocity-dependent [ (B) force potential. Statistical testing was run on average values during the series elastic element (SEE) lengthening and shortening sub-phases, and on single data points at peak SEE length (lSEE). The lower and upper parts of the box plots represent the first and third quartile, respectively. The length of the whisker delineates the minimum and maximum values. The horizontal line in the box represents the median of the sample; +, sample mean; *p < 0.05, significant side-to-side difference; ns: non-significant.
Average fascicle velocity did not differ between legs during the lengthening sub-phase, but differences were present during the SEE shortening phase (31.6 ± 28.3% lower in the affected leg, Table 2 and Figure 3C). To the time-point of peak SEE length, fascicle velocity was lower (56.5 ± 28.0%) (Table 2 and Figure 3C). No side-to-side differences were detected when velocities were expressed relative to optimal length in either sub-phase (Table 2 and Figure 4B), but being lower in the effected leg to the time-point of peak SEE length (39.2 ± 49.4%) (Table 2 and Figure 4B).
Likewise, the velocity-dependent force potential did not differ between legs during either sub-phase, while a greater velocity-dependent force potential was observed for the affected leg at the discrete time-point of peak SEE length (14.8 ± 17.7%) (Table 2 and Figure 5B).
The average pennation angle was greater in the affected leg during SEE lengthening and shortening (8 ± 3° and 12 ± 4°, respectively), which was consistent with a greater pennation angle at peak SEE length time-point in the same leg (10 ± 4°) (Table 2 and Figure 3B). Greater changes in pennation angle were also found in the affected leg during SEE lengthening (4 ± 4°) and shortening (3 ± 3°) sub-phases (Table 3). Representative ultrasound images of the affected and unaffected GM at different time-points during ground contact are presented in the Supplementary Figure 1.
Muscle shortening amplitude did not differ between the legs in either sub-phase, while the AGR was greater in the affected leg during the SEE lengthening (8.7 ± 4.9%), but not during the SEE shortening sub-phase (Table 3).
Discussion
This study aimed to understand how the behavior of the MTU is altered in former ATR patients during walking, to accommodate the longer Achilles tendon and shorter GM fascicles caused by the injury. Despite the drastic changes in the MTU architecture of the affected leg, gross walking parameters were not found to be dissimilar between the legs. The GM fascicle behavior was different between the affected and unaffected legs, however, unfavorable contractile conditions for force production during stance were not observed. Relative to their shorter length, the operating length of the affected fascicles was preserved, and the shortening velocity did not differ or was actually lower than in the healthy leg at the time of peak force production (as estimated from the timing of peak SEE length). These results indicate that under the tested walking conditions, the MTU recovery in former ATR patients may be sufficient to enable symmetrical walking gait. In addition, they suggest that long after recovery, GM fascicle behavior does not limit force generation potential and that remaining functional deficits during walking are attributable to other factors.
Leg Kinematics During Walking
The similar ankle and knee joint kinematics between the affected and unaffected legs indicate that the ATR patients had recovered a symmetrical gait pattern during treadmill walking. Two exceptions were found with a less plantarflexed average ankle joint angle during the push-off sub-phase and a lower angular velocity of that joint in the affected leg at the instant of peak SEE length. Since no side-to-side difference in angular excursion was detected during SEE shortening, the functional significance of a more dorsiflexed angle during the push-off of the affected leg is elusive. These differences may be congruent with several, but not all (Jandacka et al., 2018), reports documenting small kinematic differences in patients with similar characteristics. Several authors investigating leg kinematics during overground walking found greater dorsiflexion angles in the affected leg than in the unaffected leg at touch-down (∼1°) (Tengman and Riad, 2013; Speedtsberg et al., 2019), at toe-off (∼4°) (Willy et al., 2017), and at peak joint flexion (∼1–4°) (Don et al., 2007; Tengman and Riad, 2013; Agres et al., 2015; Speedtsberg et al., 2019). Such alterations in ankle joint kinematics, particularly greater dorsiflexion angles, could be related to the lengthened tendon caused by ATR (Agres et al., 2015). For this reason, side-to-side differences in joint kinematics may be proportional to the degree of MTU remodeling—and slackness reduction—characterizing individual recovery. Of note, the contrasting experimental modalities (walking on a treadmill or overground) between studies may have set slight differences in gait requirements, arguably less challenging in the present study. Walking on a treadmill seems indeed to reduce ground reaction forces compared to walking overground (Riley et al., 2007; Parvataneni et al., 2009), which may have lowered the GM force requirements for our patients. Regardless, the present results suggest that the remaining side-to-side differences in kinematics are small and that the gross walking pattern is relatively symmetrical. They also suggest that any difference in fascicle behavior is likely not related to dissimilarities in joint configuration.
The Gastrocnemius Medialis Fascicle Behavior Does Not Seem to Limit Force Production During Walking Stance
In support of our hypothesis, the operating GM fascicle length of the affected leg did not seem to limit potential force production. If fascicle length normalized to optimal length reflects sarcomere length, the sarcomeres of the affected leg appeared to work in a region of their force-length relation similar to that in the unaffected leg. These observations reinforce the notion that fascicle remodeling in ATR patients is achieved via a reduction in sarcomeres in series, as previously shown in animal studies on muscles immobilized at short lengths (Williams and Goldspink, 1973).
Unlike their operating length, we expected the fascicles of the affected leg to operate in an unfavorable region of their force-velocity relationship, explaining some of the functional limitations previously observed during walking. In fact, our data demonstrate that for most of the stance phase, fascicles do not operate at different normalized velocities between the legs. Conversely, affected fascicles operate in a more favorable region of their force-velocity relationship at peak SEE length, which is supported by a 15 ± 18% higher velocity-dependent force potential.
The observation of a slower contractile velocity also seems consistent with the slower ankle joint angular velocity at the time of peak SEE length. Such a reduction was likely not caused by an altered amount of fascicle shortening, since we did not find any difference in this parameter during the lengthening sub-phase. Therefore, it is difficult to elucidate the mechanisms responsible for the unchanged or reduced velocity from the present data.
One architectural factor may explain how the shorter fascicles of the affected leg could avoid contracting at a higher velocity than in the healthy leg: their rotation about their insertion point (i.e., change in pennation angle). This parameter being larger (4 ± 4°) in the affected leg may have limited the required fascicle contraction amplitude for a given muscle shortening, thereby providing a greater AGR. The present data support this interpretation, with a 9 ± 5% greater AGR in the affected leg during SEE lengthening. Thus, the contribution of the fascicle rotation to total muscle shortening during the SEE lengthening sub-phase accounted 16 ± 4% (1.6 ± 0.4 mm) on the affected leg and just 9 ± 4% (1.0 ± 0.5 mm) in the unaffected leg.
Variable AGR in general provides a mechanism to modulate performance during mechanically diverse functions (Azizi et al., 2008). The greater pennation angle after MTU remodeling of ATR patients may also promote fascicle rotation (Brainerd and Azizi, 2005), supporting architectural gearing. Architectural gearing may, therefore, constitute another adaptive mechanism in ATR patients, allowing maintaining the required muscle shortening velocity with a reduced fascicle velocity.
In the same line of thought, the altered mechanical properties of the GM tendon observed in ATR patients (Agres et al., 2015; Geremia et al., 2015; Stäudle et al., 2021) may also affect the contractile behavior of the muscle fascicles, via the effect of in series compliance on fascicle operating length and contraction velocity. We hypothesize that an increased stiffness counteracts the larger strain expected in the longer tendons of the affected side. This trade-off may thus limit the shortening of muscle fibers, while more compliant tendons would promote the opposite behavior. Since fascicular velocity was maintained or decreased on the affected side despite longer tendons, tendon stiffness was likely higher in the present patients, which could represent another adaptive mechanism in ATR patients to maintain function.
Functional Implications
The present findings indicate that affected GM sarcomeres still operate within a favorable length range and at favorable velocities for force production during treadmill walking. This is particularly interesting, as, contrary to our hypothesis, fascicle behavior does not appear to explain the functional limitations observed in ATR patients during walking (Tengman and Riad, 2013; Willy et al., 2017; Speedtsberg et al., 2019).
However, the question remains whether ATR patients can fully regain the plantar flexion force required for walking. Substantial atrophy (Heikkinen et al., 2017), as shown here by a 10 ± 12% reduction in GM thickness, is usually measured in the transversal plane of the affected fascicles, reflecting a reduced capacity for maximal force production compared to the unaffected leg. In the present experiment, the remodeled—smaller—GM muscle may have produced marginally less force during most of the stance phase. However, we contend that the force requirements of our walking experimental conditions were sufficiently low and met by compensatory mechanisms without substantially affecting gait symmetry. Although this study was not designed to investigate such mechanisms, an increase in muscle activation such as that previously reported (Wenning et al., 2021) likely occurred here.
It follows that locomotor deficits may be force- and gait-dependent in ATR patients. This hypothesis is consistent with the data of Jandacka et al. (2018), who did not find functional deficit in ATR patients during overground walking, but observed such deficits in the affected leg during running. Future studies should systematically investigate the effects of gait velocity and ground reaction forces on fascicle behavior and force production capacity to describe muscle function in a more comprehensive manner.
Limitations
A few aspects deserve further consideration for the interpretation of the present data. The final sample size (n = 9) was lower than the sample size (n ≥ 10) suggested after our a priori power analysis. Dropouts and exclusions of data of insufficient quality could not be prevented. However, the effect sizes of the differences discussed here are satisfactory; nevertheless, for completeness, the results should be verified in a future study against a matched control population with a greater sample size.
Furthermore, we used mean values of variables describing fascicular behavior. Although this approach omits testing differences at discrete time points, this choice was justified by two criteria. Firstly, it was our intention to capture substantial differences relative to tissue loading and unloading. Secondly, there was no obvious rationale to support the choice of discrete time points. Additionally, statistical parametric mapping analyses done a posteriori (Supplementary Figure 2; Bankosz and Winiarski, 2020) did not suggest that significant differences may have been masked by the present approach.
In general, the length- and velocity-dependent force potentials are based on the assumptions of biological consistency and best-fit approximations. Consequently, no model is exempt from inaccuracies. Our models could have been enhanced with additional data collection as previously described for the force-length relationship (de Brito Fontana and Herzog, 2016; Nikolaidou et al., 2017; Bohm et al., 2019, 2021) and the force-velocity relationship (Hauraix et al., 2015). Unfortunately, the large project encompassing this study did not allow additional measurements.
Similarly, for the calculations of the velocity-dependent force potential, maximum shortening velocity and normalized Hill constant (arel) were estimated based on normative data for fiber type distribution in GM muscles (type I fibers: 50.8%, type II fibers: 49.2%) (Johnson et al., 1973). As we cannot discard the possible influence of this parameter on side-to-side differences, speculating about the potential role of fiber type distribution in the present data should be avoided.
Conclusion
The results of this study indicate that sarcomere operating length and lower limbs kinematics are relatively preserved during walking in ATR patients more than 2 years post-surgery. Despite atrophy in the longitudinal and transversal planes of the affected GM muscle, fascicles were found to operate at contractile velocities comparable to or lower than those of the non-operated leg during the stance phase. An increased fascicle rotation during stance may contribute to preserving fascicle contractile velocity. Collectively, these results suggest that the contractile behavior of the GM muscle does not limit force production under the present walking condition and in ATR patients with a similar level of recovery.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the institutional review board of the German Sport University Cologne. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
G-PB and KA conceptualized the research. BS, GL, G-PB, and KA designed the research. BS acquired the data. BS and KA analyzed the data. BS, OS, and KA interpreted the data and drafted the manuscript. BS, OS, GL, G-PB, and KA revised the manuscript. All authors approved the final manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.792576/full#supplementary-material
Abbreviations
AGR, architectural gear ratio; ATR, Achilles tendon rupture; FT, fast twitch; ICC, intra-class correlation coefficient; GM, gastrocnemius medialis; MTU, muscle-tendon unit; SEE, series elastic element.
References
Aggeloussis, N., Giannakou, E., Albracht, K., and Arampatzis, A. (2010). Reproducibility of fascicle length and pennation angle of gastrocnemius medialis in human gait in vivo. Gait Posture 31, 73–77. doi: 10.1016/j.gaitpost.2009.08.249
Agres, A. N., Duda, G. N., Gehlen, T. J., Arampatzis, A., Taylor, W. R., and Manegold, S. (2015). Increased unilateral tendon stiffness and its effect on gait 2-6 years after Achilles tendon rupture. Scand. J. Med. Sci. Sports 25, 860–867. doi: 10.1111/sms.12456
Azizi, E., Brainerd, E. L., and Roberts, T. J. (2008). Variable gearing in pennate muscles. Proc. Natl. Acad. Sci. U.S.A. 105, 1745–1750. doi: 10.1073/pnas.0709212105
Bahler, A. S., Fales, J. T., and Zierler, K. L. (1968). The dynamic properties of mammalian skeletal muscle. J. Gen. Physiol. 51, 369–384. doi: 10.1085/jgp.51.3.369
Bankosz, Z., and Winiarski, S. (2020). Statistical parametric mapping reveals subtle gender differences in angular movements in table tennis topspin backhand. Int. J. Environ. Res. Public Health 17:6996. doi: 10.3390/ijerph17196996
Baxter, J. R., Hullfish, T. J., and Chao, W. (2018). Functional deficits may be explained by plantarflexor remodeling following Achilles tendon rupture repair: preliminary findings. J. Biomech. 79, 238–242. doi: 10.1016/j.jbiomech.2018.08.016
Bohm, S., Mersmann, F., Santuz, A., and Arampatzis, A. (2019). The force-length-velocity potential of the human soleus muscle is related to the energetic cost of running. Proc. Biol. Sci. 286:20192560. doi: 10.1098/rspb.2019.2560
Bohm, S., Mersmann, F., Santuz, A., and Arampatzis, A. (2021). Enthalpy efficiency of the soleus muscle contributes to improvements in running economy. Proc. Biol. Sci. 288:20202784. doi: 10.1098/rspb.2020.2784
Brainerd, E. L., and Azizi, E. (2005). Muscle fiber angle, segment bulging and architectural gear ratio in segmented musculature. J. Exp. Biol. 208, 3249–3261. doi: 10.1242/jeb.01770
Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences. Amsterdam: Elsevier Science.
Cronin, N. J., and Lichtwark, G. (2013). The use of ultrasound to study muscle-tendon function in human posture and locomotion. Gait Posture 37, 305–312. doi: 10.1016/j.gaitpost.2012.07.024
Cureton, E. E. (1956). Rank-biserial correlation. Psychometrika 21, 287–290. doi: 10.1007/bf02289138
de Brito Fontana, H., and Herzog, W. (2016). Vastus lateralis maximum force-generating potential occurs at optimal fascicle length regardless of activation level. Eur. J. Appl. Physiol. 116, 1267–1277. doi: 10.1007/s00421-016-3381-3
Don, R., Ranavolo, A., Cacchio, A., Serrao, M., Costabile, F., Iachelli, M., et al. (2007). Relationship between recovery of calf-muscle biomechanical properties and gait pattern following surgery for achilles tendon rupture. Clin. Biomech. (Bristol, Avon) 22, 211–220. doi: 10.1016/j.clinbiomech.2006.10.001
Farris, D. J., and Lichtwark, G. A. (2016). UltraTrack: software for semi-automated tracking of muscle fascicles in sequences of B-mode ultrasound images. Comput. Methods Programs Biomed. 128, 111–118. doi: 10.1016/j.cmpb.2016.02.016
Farris, D. J., and Sawicki, G. S. (2012). Human medial gastrocnemius force-velocity behavior shifts with locomotion speed and gait. Proc. Natl. Acad. Sci. U.S.A. 109, 977–982. doi: 10.1073/pnas.1107972109
Farris, D. J., Lichtwark, G. A., Brown, N. A., and Cresswell, A. G. (2016). The role of human ankle plantar flexor muscle-tendon interaction and architecture in maximal vertical jumping examined in vivo. J. Exp. Biol. 219, 528–534. doi: 10.1242/jeb.126854
Faul, F., Erdfelder, E., Lang, A. G., and Buchner, A. (2007). G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav. Res. Methods 39, 175–191. doi: 10.3758/bf03193146
Fukunaga, T., Kubo, K., Kawakami, Y., Fukashiro, S., Kanehisa, H., and Maganaris, C. N. (2001). In vivo behaviour of human muscle tendon during walking. Proc. Biol. Sci. 268, 229–233. doi: 10.1098/rspb.2000.1361
Geremia, J. M., Bobbert, M. F., Casa Nova, M., Ott, R. D., Lemos Fde, A., Lupion Rde, O., et al. (2015). The structural and mechanical properties of the Achilles tendon 2 years after surgical repair. Clin. Biomech. (Bristol, Avon) 30, 485–492. doi: 10.1016/j.clinbiomech.2015.03.005
Gollapudi, S. K., and Lin, D. C. (2009). Experimental determination of sarcomere force-length relationship in type-I human skeletal muscle fibers. J. Biomech. 42, 2011–2016. doi: 10.1016/j.jbiomech.2009.06.013
Hauraix, H., Nordez, A., Guilhem, G., Rabita, G., and Dorel, S. (2015). In vivo maximal fascicle-shortening velocity during plantar flexion in humans. J. Appl. Physiol. (1985) 119, 1262–1271. doi: 10.1152/japplphysiol.00542.2015
Hawkins, D., and Hull, M. L. (1990). A method for determining lower extremity muscle-tendon lengths during flexion/extension movements. J. Biomech. 23, 487–494. doi: 10.1016/0021-9290(90)90304-l
Heikkinen, J., Lantto, I., Piilonen, J., Flinkkila, T., Ohtonen, P., Siira, P., et al. (2017). Tendon length, calf muscle atrophy, and strength deficit after acute achilles tendon rupture: long-term follow-up of patients in a previous study. J. Bone Joint Surg. Am. 99, 1509–1515. doi: 10.2106/JBJS.16.01491
Herzog, W., Abrahamse, S. K., and Ter Keurs, H. E. (1990). Theoretical determination of force-length relations of intact human skeletal muscles using the cross-bridge model. Pflugers Arch. 416, 113–119. doi: 10.1007/BF00370231
Hill, A. V. (1938). The heat of shortening and the dynamic constants of muscle. Proc. R. Soc. Lond. Ser. B Biol. Sci. 126, 136–195. doi: 10.1152/advan.00072.2005
Ishikawa, M., Pakaslahti, J., and Komi, P. V. (2007). Medial gastrocnemius muscle behavior during human running and walking. Gait Posture 25, 380–384. doi: 10.1016/j.gaitpost.2006.05.002
Jandacka, D., Plesek, J., Skypala, J., Uchytil, J., Silvernail, J. F., and Hamill, J. (2018). Knee joint kinematics and kinetics during walking and running after surgical achilles tendon repair. Orthop. J. Sports Med. 6:2325967118779862. doi: 10.1177/2325967118779862
Jandacka, D., Silvernail, J. F., Uchytil, J., Zahradnik, D., Farana, R., and Hamill, J. (2017). Do athletes alter their running mechanics after an Achilles tendon rupture? J. Foot Ankle Res. 10:53. doi: 10.1186/s13047-017-0235-0
Johnson, M. A., Polgar, J., Weightman, D., and Appleton, D. (1973). Data on the distribution of fibre types in thirty-six human muscles. An autopsy study. J. Neurol. Sci. 18, 111–129. doi: 10.1016/0022-510x(73)90023-3
Kerby, D. S. (2014). The simple difference formula: an approach to teaching nonparametric correlation. Compr. Psychol. 3, 2165–2228.
Millard, M. (2021). OpenSim:: ActiveForceLengthCurve Class Reference. In API Guide for C++ developers. Available online at: https://simtk.org/api_docs/opensim/api_docs/classOpenSim_1_1ActiveForceLengthCurve.html (accessed September 15, 2021).
Millard, M., Uchida, T., Seth, A., and Delp, S. L. (2013). Flexing computational muscle: modeling and simulation of musculotendon dynamics. J. Biomech. Eng. 135:021005. doi: 10.1115/1.4023390
Mullaney, M. J., Mchugh, M. P., Tyler, T. F., Nicholas, S. J., and Lee, S. J. (2006). Weakness in end-range plantar flexion after Achilles tendon repair. Am. J. Sports Med. 34, 1120–1125. doi: 10.1177/0363546505284186
Muramatsu, T., Muraoka, T., Kawakami, Y., Shibayama, A., and Fukunaga, T. (2002). In vivo determination of fascicle curvature in contracting human skeletal muscles. J. Appl. Physiol. (1985) 92, 129–134. doi: 10.1152/jappl.2002.92.1.129
Neptune, R. R., and Sasaki, K. (2005). Ankle plantar flexor force production is an important determinant of the preferred walk-to-run transition speed. J. Exp. Biol. 208, 799–808. doi: 10.1242/jeb.01435
Nigg, B. M., and Herzog, W. (1994). Biomechanics Of The Musculo-Skeletal System. Hoboken, NJ: Wiley.
Nikolaidou, M. E., Marzilger, R., Bohm, S., Mersmann, F., and Arampatzis, A. (2017). Operating length and velocity of human M. vastus lateralis fascicles during vertical jumping. R. Soc. Open Sci. 4:170185. doi: 10.1098/rsos.170185
Nilsson-Helander, K., Silbernagel, K. G., Thomee, R., Faxen, E., Olsson, N., Eriksson, B. I., et al. (2010). Acute achilles tendon rupture: a randomized, controlled study comparing surgical and nonsurgical treatments using validated outcome measures. Am. J. Sports Med. 38, 2186–2193. doi: 10.1177/0363546510376052
Parvataneni, K., Ploeg, L., Olney, S. J., and Brouwer, B. (2009). Kinematic, kinetic and metabolic parameters of treadmill versus overground walking in healthy older adults. Clin. Biomech. (Bristol, Avon) 24, 95–100. doi: 10.1016/j.clinbiomech.2008.07.002
Peng, W. C., Chao, Y. H., Fu, A. S. N., Fong, S. S. M., Rolf, C., Chiang, H., et al. (2019). Muscular morphomechanical characteristics after an achilles repair. Foot Ankle Int. 40, 568–577. doi: 10.1177/1071100718822537
Pohle-Fröhlich, R., Dalitz, C., Richter, C., Hahnen, T., Stäudle, B., and Albracht, K. (2020). “Estimation of muscle fascicle orientation in ultrasonic images,” in Proceedings of the 15th International Joint conference on Computer Vision, Imaging and computer Graphics Theory and Applications, Valletta, 79–86.
Riley, P. O., Paolini, G., Della Croce, U., Paylo, K. W., and Kerrigan, D. C. (2007). A kinematic and kinetic comparison of overground and treadmill walking in healthy subjects. Gait Posture 26, 17–24. doi: 10.1016/j.gaitpost.2006.07.003
Robertson, D. G. E., Caldwell, G. E., Hamill, J., Kamen, G., and Whittlesey, S. (2013). Research Methods in Biomechanics. Champaign, IL: Human Kinetics.
Sanchez, G. N., Sinha, S., Liske, H., Chen, X., Nguyen, V., Delp, S. L., et al. (2015). In vivo imaging of human sarcomere twitch dynamics in individual motor units. Neuron 88, 1109–1120. doi: 10.1016/j.neuron.2015.11.022
Seth, A., Hicks, J. L., Uchida, T. K., Habib, A., Dembia, C. L., Dunne, J. J., et al. (2018). OpenSim: simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. PLoS Comput. Biol. 14:e1006223. doi: 10.1371/journal.pcbi.1006223
Silbernagel, K. G., Steele, R., and Manal, K. (2012). Deficits in heel-rise height and achilles tendon elongation occur in patients recovering from an Achilles tendon rupture. Am. J. Sports Med. 40, 1564–1571. doi: 10.1177/0363546512447926
Speedtsberg, M. B., Kastoft, R., Barfod, K. W., Penny, J. O., and Bencke, J. (2019). Gait function and postural control 4.5 years after nonoperative dynamic treatment of acute achilles tendon ruptures. Orthop. J. Sports Med. 7:2325967119854324. doi: 10.1177/2325967119854324
Stäudle, B., Seynnes, O., Laps, G., Goll, F., Bruggemann, G. P., and Albracht, K. (2021). Recovery from achilles tendon repair: a combination of postsurgery outcomes and insufficient remodeling of muscle and tendon. Med. Sci. Sports Exerc. 53, 1356–1366. doi: 10.1249/MSS.0000000000002592
Svensson, R. B., Couppé, C., Agergaard, A.-S., Josefsen, C. O., Jensen, M. H., Weisskirchner Barfod, K., et al. (2019). Persistent functional loss following ruptured Achilles tendon is associated with reduced gastrocnemius muscle fascicle length, elongated gastrocnemius and soleus tendon, and reduced muscle cross-sectional area. Transl. Sports Med. 2, 316–324. doi: 10.1002/tsm2.103
Tengman, T., and Riad, J. (2013). Three-dimensional gait analysis following achilles tendon rupture with nonsurgical treatment reveals long-term deficiencies in muscle strength and function. Orthop. J. Sports Med. 1:2325967113504734. doi: 10.1177/2325967113504734
Thomas, S., Reading, J., and Shephard, R. J. (1992). Revision of the physical activity readiness questionnaire (PAR-Q). Can. J. Sport Sci. 17, 338–345.
Ward, S. R., Eng, C. M., Smallwood, L. H., and Lieber, R. L. (2009). Are current measurements of lower extremity muscle architecture accurate? Clin. Orthop. Relat. Res. 467, 1074–1082. doi: 10.1007/s11999-008-0594-8
Wenning, M., Mauch, M., Heitner, A., Lienhard, J., Ritzmann, R., and Paul, J. (2021). Neuromechanical activation of triceps surae muscle remains altered at 3.5 years following open surgical repair of acute Achilles tendon rupture. Knee Surg. Sports Traumatol. Arthrosc. 29, 2517–2527. doi: 10.1007/s00167-021-06512-z
Werkhausen, A., Cronin, N. J., Albracht, K., Paulsen, G., Larsen, A. V., Bojsen-Moller, J., et al. (2019). Training-induced increase in Achilles tendon stiffness affects tendon strain pattern during running. PeerJ 7:e6764. doi: 10.7717/peerj.6764
Williams, P. E., and Goldspink, G. (1973). The effect of immobilization on the longitudinal growth of striated muscle fibres. J. Anat. 116, 45–55.
Willy, R. W., Brorsson, A., Powell, H. C., Willson, J. D., Tranberg, R., and Gravare Silbernagel, K. (2017). Elevated knee joint kinetics and reduced ankle kinetics are present during jogging and hopping after achilles tendon ruptures. Am. J. Sports Med. 45, 1124–1133. doi: 10.1177/0363546516685055
Winters, J. M., and Stark, L. (1988). Estimated mechanical properties of synergistic muscles involved in movements of a variety of human joints. J. Biomech. 21, 1027–1041. doi: 10.1016/0021-9290(88)90249-7
Keywords: tendon rupture, muscle fascicle behavior, walking gait, force generation, ultrasound imaging
Citation: Stäudle B, Seynnes O, Laps G, Brüggemann G-P and Albracht K (2022) Altered Gastrocnemius Contractile Behavior in Former Achilles Tendon Rupture Patients During Walking. Front. Physiol. 13:792576. doi: 10.3389/fphys.2022.792576
Received: 10 October 2021; Accepted: 01 February 2022;
Published: 01 March 2022.
Edited by:
Christoph Centner, University of Freiburg, GermanyReviewed by:
Fábio Juner Lanferdini, Federal University of Santa Catarina, BrazilMartino V. Franchi, University of Padua, Italy
Markus Tilp, University of Graz, Austria
Jeam Marcel Geremia, Federal University of Rio Grande do Sul, Brazil
Copyright © 2022 Stäudle, Seynnes, Laps, Brüggemann and Albracht. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Benjamin Stäudle, staeudle@fh-aachen.de