ABSTRACT
The Achilles is a long tendon varying in cross-sectional area (CSA) considerably along its length. For the same force, a smaller CSA would experience higher tendon stress and we hypothesised that these areas would therefore undergo larger transverse deformations. A novel magnetic resonance imaging-based approach was implemented to quantify changes in tendon CSA from rest along the length of the Achilles tendon under load conditions corresponding to 10%, 20% and 30% of isometric plantar flexor maximum voluntary contraction (MVC). Reductions in tendon CSA occurring during contraction from the resting condition were assumed to be proportional to the longitudinal elongations within those regions (Poisson's ratio). Rather than tendon regions of smallest CSA undergoing the greatest deformations, the outcome was region specific, with the proximal (gastrocnemius) tendon portion showing larger transverse deformations upon loading compared with the distal portion of the Achilles (P<0.01). Transverse tendon deformation only occurred in selected regions of the distal Achilles tendon at 20% and 30% of MVC, but in contrast occurred throughout the proximal portion of the Achilles at all contraction levels (10%, 20% and 30% of MVC; P<0.01). Calculations showed that force on the proximal tendon portion was ∼60% lower, stress ∼70% lower, stiffness ∼30% lower and Poisson's ratio 6-fold higher compared with those for the distal portion of the Achilles tendon. These marked regional differences in mechanical properties may allow the proximal portion to function as a mechanical buffer to protect the stiffer, more highly stressed, distal portion of the Achilles tendon from injury.
INTRODUCTION
The technique for testing human tendon mechanical properties in vivo has developed over the last ∼15 years based upon the general principles of tendon tensile tests conducted in vitro (Butler et al., 1978; Cuming et al., 1978; Rigby, 1964). In vivo approaches have typically involved measuring longitudinal tendon elongations during isometric contraction using ultrasound imaging and relating these elongations to tendon force estimations derived from dynamometry measurements of joint torque (Maganaris and Paul, 1999). A necessary simplification with this approach is the measurement of tendon elongations from a single anatomical point. Typically for ‘free’ tendons, this measurement point is the myo-tendinous junction (MTJ), or an osteo-tendinous interface, as these are clear ‘reference’ points that can be identified and tracked with ultrasound imaging. Intra-muscular tendon–aponeurosis sites have also been used for tracking tendon elongation at the site where fascicles intersect the aponeurosis. Although likely to be an oversimplification, the use of a single ‘reference’ point is a necessary step in obtaining measurements to permit calculation of tendon stiffness and modulus, with many current variations on the in vivo approach. It also allows specific questions to be answered relating to the variations in tendon properties between different populations (Hansen et al., 2003; Karamanidis and Arampatzis, 2006; Maganaris et al., 2006), or changes in tendon properties as a result of interventions (Arampatzis et al., 2007; Kubo et al., 2002; Reeves et al., 2005; Wiesinger et al., 2015). However, the implicit assumption with this approach is that tendon elongations are homogeneous, or at least representative of the entire tendon, which may not always be the case.
For relatively long tendons, in vitro studies have shown that variations in elongation occur along the length of the tendon (Wren et al., 2001; Zernicke et al., 1984). In vivo measurements also confirm the non-uniformity of tendon elongations along the length of the free tendon or tendon–aponeurosis in the human lower limb (Finni et al., 2003; Maganaris and Paul, 2000; Magnusson et al., 2003). Using an ultrasound speckle tracking approach, the human Achilles tendon has been shown to have differences in elongation between superficial and deeper regions (Arndt et al., 2012; Slane and Thelen, 2015). It might be hypothesised that elongation inhomogeneities should be particularly evident in long tendons where the cross-sectional area (CSA) varies considerably along the tendon length. This is because regions of the tendon with smaller CSAs are expected to experience higher stress (stress=force/area) and therefore are likely to experience greater transverse and longitudinal elongations. The CSA of the Achilles tendon varies considerably along its length (Finni et al., 2003; Magnusson and Kjaer, 2003), and we might hypothesise that because of this, the tendon will experience different elongations along its length, with smaller CSA regions experiencing larger transverse and longitudinal elongations.
The Achilles tendon is composed of two tendon components, with the more proximal gastrocnemius tendon ‘fusing’ with the soleus tendon just distal to the soleus muscle (Cummins and Anson, 1946). The ‘fusing’ of these two tendon components in a spiral manner (Cummins and Anson, 1946) also highlights the potential for shear to occur within the Achilles tendon, which may further contribute to non-uniform elongations along its length and in part explain changes seen with speckle tracking (Arndt et al., 2012; Slane and Thelen, 2015) and other ultrasound-based approaches (Bojsen-Moller et al., 2004; Magnusson et al., 2003). Larger elongations occurring at certain regions along the tendon may indicate a propensity for tendon strain injuries. If these regions coincide with portions of the tendon where the CSA is smaller, this may facilitate understanding of the mechanisms of tendon injuries and ruptures. The most frequent site (in 85% of cases) for complete rupture of the Achilles tendon is reported to be the region 3–5 cm proximal to the calcaneus (Józsa et al., 1989). This also coincides with the region of the Achilles tendon where the CSA is at its smallest (Magnusson and Kjaer, 2003), supporting the hypothesis for a link between regions of higher stress, greater transverse/longitudinal elongations and propensity for tendon injury.
- Ao
average unloaded tendon section cross-sectional area (cm2)
- Aave
average loaded tendon section cross-sectional area (cm2)
- CSA
cross-sectional area (cm2)
- E
average Young's modulus of tendon section (MPa)
- EMG
electromyographic
- F
force on the tendon section (N)
- K
average stiffness of tendon section (N mm−1)
- Lo
original tendon section length (mm)
- MRI
magnetic resonance imaging
- MTJ
myo-tendinous junction
- MVC
maximal voluntary contraction
- PCSA
physiological cross-sectional area
- Vo
original volume of tendon section
- ΔL
change in tendon section length (mm)
- ΔV
change in volume of tendon section
- Δx
extension of tendon section
- εx
average strain of tendon section
- ν
Poisson's ratio for tendon section
- σ
average engineering stress in tendon section (MPa)
As current conventional in vivo techniques cannot distinguish transverse length-dependent deformations, we adopted a novel magnetic resonance imaging (MRI)-based approach to assess tendon deformation occurring along the length of the Achilles tendon, with the hypothesis that larger transverse deformation would be observed at regions of the tendon where the CSA is smaller. This experimental approach involved measuring tendon CSA changes in the transverse plane with the assumption that tensile loading causes a reduction in CSA proportional to the longitudinal elongation (Poisson's ratio); essentially, as the tendon is stretched it becomes proportionally thinner.
MATERIALS AND METHODS
Participants
Nine male participants (mean±s.d. age: 25±6 years, body mass: 80±10 kg and height: 1.81±0.05 m) gave informed consent to take part, after the study had received institutional ethics committee approval. Exclusion criteria included: contraindication to MRI scanning, prior lower limb surgery, tendinopathy affecting the Achilles or tendons of the lower limbs, and serious injury to the Achilles tendon or lower limb.
Determining the statistical power with the results obtained tested the adequacy of the sample size against the optimal of 80% power recommended by Cohen (1988). In the proximal (gastrocnemius) tendon region, statistical power was 100% and ranged between 40–79% in the region of smallest CSA.
Dynamometry
Participants lay prone on the chair of an isokinetic dynamometer (Cybex Norm, NY, USA) with the right foot fixed into the footplate at a neutral ankle position [i.e. 1.57 rad (90 deg) between foot plate and tibia] and the knee in full extension. Participants first performed isometric maximal voluntary contraction (MVC) of the plantar flexors on the dynamometer. Two isometric MVCs were performed (an additional contraction was performed if the two MVC torque values were not within 5% of each other) and the mean torque value was calculated. Joint torque values corresponding to 10%, 20% and 30% of a participant's MVC were then calculated and visually highlighted on a screen for the participant, displaying the torque trace in real time, linked to an acquisition system (Biopac Systems Inc., CA, USA). Using visual feedback, participants were then asked to perform a plantarflexion contraction, developing the torque value corresponding to 10%, 20% and 30% of their MVC (in a randomised order) and maintain this value for 3 min at each percentage of their MVC (with 3 min rest between MVC levels). Simultaneously during these plantarflexion contractions, electromyographic (EMG) activity was measured from the tibialis anterior muscle to assess the level of co-contraction (see ‘Antagonist muscle co-contraction’, below) and the MTJ of the gastrocnemius muscle–tendon unit was tracked using B-mode ultrasound imaging to assess tendon elongation (see ‘Measurement of tendon elongation and assessment of tendon creep’, below).
Antagonist muscle co-contraction
To estimate the level of dorsiflexor co-activation during plantarflexion efforts, the EMG activity was measured from a representative dorsiflexor muscle (tibialis anterior). The root mean square EMG activity was calculated from the raw signal and related to that measured from the same muscle (tibialis anterior) when acting as an agonist during dorsiflexion contractions at various percentages of the dorsiflexion MVC. The level of dorsiflexor co-activation and the corresponding co-activation torque during plantarflexion efforts at 10%, 20% and 30% of MVC were then calculated in line with previously described methods (Maganaris et al., 1998b).
Measurement of longitudinal tendon elongation and assessment of tendon creep
During dynamometry testing, the medial head of the gastrocnemius MTJ was scanned at 25 Hz using ultrasound imaging (Mylab 70, Esaote, Italy), in line with previously described methods (Maganaris and Paul, 2002; Reeves et al., 2005), to represent the longitudinal elongations of the whole Achilles tendon. It was previously reported that there are no significant differences in Achilles tendon stiffness and strain between scanning at the MTJ of the medial or lateral head of the gastrocnemius (Morrison et al., 2015). The ultrasound probe was secured in position using a custom-made holder to prevent any movement relative to the scanned structure, and an external marker placed onto the skin casting a line on the image and confirming that no movement took place. The longitudinal displacement of the gastrocnemius MTJ was tracked continuously during the 3 min plantarflexion contraction performed at each contraction level (10%, 20% and 30% of MVC). Scans were then digitised offline using Image J (National Institutes of Health, Bethesda, MD, USA) at 1 min intervals (i.e. tendon elongation was measured at 0, 1, 2 and 3 min time intervals). These scans were acquired to assess the degree of overall longitudinal Achilles tendon elongation and to identify whether any tendon creep was evident during these sustained contractions.
MRI scanning
Participants were positioned supine, with the knees fully extended within a 0.25 T MRI scanner (G-Scan, Esaote, Italy). Axial plane scans of participant's right lower leg were acquired using a Spin-Echo Fast Fourier sequence with the following scanning parameters: 1:59 min scanning time, 18 ms echo time, 1020 ms repetition time, 1 acquisition, 180×170 mm field of view, 256×256 pixels, 7 mm slice thickness, 1 mm inter-slice gap. These scans were acquired (foot fixed in the neutral position) with the participant at rest and also while developing an isometric plantarflexion contraction (with the knee fully extended), with a torque equivalent to 10%, 20% and 30% of their MVC (determined from the dynamometry measurements described above). Development of these prescribed isometric joint torque levels (10%, 20% and 30% of MVC) was achieved by using a custom-made, MRI-compatible lever system. The participant's external lever arm (defined as the distance from the ankle joint centre of rotation to the point of force application at the first metatarsal head) and the device lever arm were carefully measured using a ruler incremented in millimetres. Loads that would achieve the required plantarflexion torque levels were then calculated and positioned on one end of the lever system. The other end of the lever system was in contact with the participant's foot (the first metatarsal head region), where the application of force acted to secure the system in place around a pivot point on the MRI bed. When the scan began, an experimenter positioned inside the MRI released the load. Participants then needed to develop a plantarflexion joint torque (isometric contraction) equivalent to the prescribed level that would balance the load within a small, tightly controlled range, verified visually by the experimenter. This process was repeated for each level of joint torque (10%, 20% and 30% of MVC). Axial plane scans were acquired starting from the calcaneus and continuing ∼17 cm proximal from this point, capturing the entire length of the distal portion of the Achilles tendon and the vast majority of the proximal (gastrocnemius) tendon portion (Fig. 1). Moving proximally from the calcaneus, the first appearance of the Achilles tendon was designated as ‘scan 0’. Using digitising software (OsiriX, Pixmeo, Geneva, Switzerland), the CSA of the Achilles tendon was then measured on scans 2–21 (Figs 1 and 2). The distal portion of the Achilles tendon was defined as that composed of both soleus and gastrocnemius tendons. Moving proximally from the calcaneus, the first appearance of the soleus muscle (i.e. the soleus MTJ) therefore delineated where the distal portion of the Achilles tendon ended and the proximal portion (gastrocnemius) of the Achilles tendon began (Fig. 1).
Sagittal plane MRI scans were acquired for the purpose of quantifying the Achilles tendon moment arm length, which was used for calculating tendon force by dividing the measured plantarflexion joint torque by the Achilles tendon moment arm length. Scans were taken with the foot in the neutral ankle position, in 0.17 rad (10 deg) of dorsiflexion and 0.17 rad (10 deg) of plantarflexion. During the scans, participants were asked to perform a plantarflexion contraction to the level of ∼20% MVC so that the Achilles tendon moment arm was measured under load, because the Achilles tendon moment arm length is known to increase when a tensile load is applied compared with the resting state (Maganaris et al., 1998a). The Achilles tendon moment arm length was calculated for a neutral ankle position using the Reuleux method as previously described in detail (Maganaris et al., 1998a). Briefly, scans taken in dorsiflexion and plantarflexion were used to identify the instant centre of rotation on the talus bone in the neutral ankle scan. The Achilles tendon moment arm length was then measured as the perpendicular distance between the Achilles tendon action line and the joint centre of rotation on the talus in the neutral ankle scan.
Estimation of mechanical properties of the Achilles tendon proximal and distal components
Mechanical and material properties were calculated separately for the distal and proximal portions of the Achilles tendon as described below. Overall Achilles tendon length was estimated using the participant's height and multiplying by a factor of 0.00792 (calculated from tendon data from Reeves et al., 2005). MRI scans were used to measure the length of the distal portion of the Achilles tendon by summing the number of slices from the calcaneus insertion to the soleus MTJ (7 mm slice thickness+1 mm gap=8 mm resolution). The length of the proximal portion of the Achilles tendon was estimated by subtracting the length of the distal portion from the overall Achilles tendon length. The overall longitudinal tendon elongation was measured using ultrasound at the gastrocnemius medial MTJ and longitudinal elongation for each section, ΔL, was estimated by assuming that elongation at the soleus MTJ was 30% less than elongation of the overall Achilles tendon (Bojsen-Moller et al., 2004).
Statistical analysis
A repeated measures analysis of variance (ANOVA) with a Newman–Keuls multiple comparison post hoc test was applied to test for differences between contraction conditions (rest, 10%, 20% and 30% MVC). This statistical approach was used for the parameters of tendon CSA, tendon longitudinal elongation, plantar flexor joint torque and dorsiflexor muscle co-activation. A paired samples Student's t-test was used to test for differences in mechanical properties between the distal and proximal portions of the Achilles tendon. Values presented are means±s.d.
RESULTS
Regional Achilles tendon CSA changes upon loading
At 10% of MVC, CSAs were significantly smaller throughout the proximal Achilles tendon component compared with those at rest (P<0.01), but there was no significant difference in the distal Achilles tendon CSAs (Fig. 2A). At 20% and 30% of MVC, CSAs throughout the proximal Achilles tendon component remained significantly smaller compared with those at rest (P<0.01; Fig. 2B,C). At 20% of MVC, the most superior CSA of the distal Achilles tendon component (scan 9) was significantly smaller compared with that at rest (P<0.01), with no significant difference in any other CSAs throughout the distal Achilles tendon component (Fig. 2B). At 30% of MVC, CSAs in the superior region of the distal Achilles tendon (scans 4–9) were significantly smaller compared with those at rest (P<0.01; Fig. 2C), with no significant difference at the most inferior CSAs of the distal Achilles tendon component (scans 2 and 3).
Longitudinal Achilles tendon elongations
Longitudinal elongations of the whole Achilles tendon measured using ultrasound are shown in Fig. 3 (see below). Elongations of the whole Achilles tendon measured at 0, 1, 2 and 3 min during the sustained isometric contractions were not significantly different between any of these time points, for any of the contraction levels (Fig. 3).
Plantarflexion joint torque and Achilles tendon moment arm
Table 1 shows how participants were able to be very accurate with matching the target plantarflexion torque over the duration of the 3 min contraction on the dynamometer, with little variation around the mean. Table 2 shows the Achilles tendon moment arm lengths for each individual participant as well as the group mean and s.d.
Dorsiflexor co-activation
The level of dorsiflexor coactivation during plantar flexor efforts was <1% expressed as a proportion of the maximum agonist dorsiflexor EMG and <1 N m when the torque contribution from this level of co-activation was calculated. There were no significant differences in co-activation between contraction levels (% MVC) or over time during the sustained isometric contraction.
Mechanical properties of the distal and proximal components of the Achilles tendon
Table 3 shows the mechanical properties for the distal and proximal components of the Achilles tendon at the three different loading levels (10%, 20% and 30% MVC). The proximal component of the Achilles tendon was significantly longer (Lo), experienced ∼60% lower force (F), ∼70% lower stress (σ), ∼30% lower stiffness (K) and 6-fold higher Poisson's ratio (ν) compared with the distal Achilles tendon component (Table 3).
DISCUSSION
In this study, we adopted a novel MRI-based approach to study transverse tendon deformations occurring throughout the length of the Achilles tendon, with the initial hypothesis that deformation would be greater where tendon CSA is smaller. This experimental approach involved measuring tendon CSA changes in the transverse plane with the assumption that tensile loading causes a reduction in CSA proportional to the longitudinal elongation (Poisson's ratio); essentially, as the tendon is stretched it becomes proportionally thinner. Our initial hypothesis, that the tendon would undergo larger deformations in regions of smaller CSA, was based upon the assumption that material properties and forces were similar throughout the tendon's length. Therefore, areas where the tendon CSA was smaller would experience higher stress and undergo larger transverse and longitudinal deformations. The present results do not support this hypothesis, but instead support the intriguing notion that Achilles tendon deformation is region specific and closely related to differences in mechanical properties between the proximal (gastrocnemius) and distal components of the Achilles tendon (Figs 1 and 2, Table 3).
In the present study, we applied tensile loading via isometric muscle contraction at 10%, 20% and 30% of MVC. At 10% of MVC, we found reductions in CSA throughout the length of the proximal (gastrocnemius) Achilles tendon portion, but no changes in the distal portion of the Achilles tendon (Fig. 2A). Selected superior regions of the distal portion of the Achilles tendon only began to show reductions in CSA upon loading at 20% and 30% of MVC, with no changes to inferior regions of the distal tendon portion (Fig. 2B,C). These results show that rather than the region of smallest tendon CSA undergoing larger transverse deformations due to tensile loading, there is a marked difference between proximal and distal regions of the Achilles tendon. Specifically, in the proximal region of the Achilles tendon, force was ∼60% lower, stress was ∼70% lower, stiffness was ∼30% lower and Poisson's ratio was 6-fold higher compared with the distal portion of the Achilles tendon (Fig. 2, Table 3). Substantially greater transverse strain in the proximal compared with the distal region of the Achilles tendon explains the 6-fold higher Poisson's ratio value for this proximal Achilles tendon region.
The substantially lower forces acting on the proximal portion of the Achilles tendon (∼60% lower) can be explained because only forces generated proximal to the gastrocnemius MTJ by the gastrocnemius muscle act to deform this proximal Achilles tendon region. In contrast, the more distal portion of the Achilles tendon experiences much higher forces generated by both proximally located muscles: the gastrocnemius and soleus (Table 3). Substantially lower forces in the proximal Achilles tendon region contribute to explaining the markedly lower stiffness (∼30% lower) compared with the distal tendon portion (Table 3). The lower stiffness of the proximal Achilles tendon region may allow it to function as a ‘buffer’, absorbing high strains induced by high force eccentric contractions and thereby protecting the stiffer, more highly stressed, distal portion of the Achilles tendon. The capacity of the tendon to act as a buffer and protect the muscle fascicles from injury by attenuating the force and slowing the rate of lengthening for muscle fascicles has previously been identified (Konow and Roberts, 2015; Konow et al., 2012; Roberts and Konow, 2013). We speculate that the proximal portion of the Achilles tendon may function as a buffer not only to protect the muscle fascicles but also to protect the stiffer, more highly stressed, distal portion of the Achilles tendon.
A tendon's modulus reflects the intrinsic ‘material’ properties of the tendon and is given by the stiffness normalised to the dimensions of the tendon (Maganaris and Paul, 1999; Reeves et al., 2003). As the CSAs of these two tendon regions (proximal and distal Achilles tendon regions) are fairly equivalent (see Fig. 2 and Table 3) and our findings show the stiffness of the proximal tendon region to be lower, these results might indicate a lower modulus of the proximal Achilles tendon; however, this was not found to be the case (Table 3). The stiffness values in the present study ranged between 51 and 132 N mm−1, which agrees well with a number of previously reported values for the Achilles tendon ranging between 26 and 187 N mm−1 (Arampatzis et al., 2007; Karamanidis and Arampatzis, 2006; Kubo et al., 2002; Maganaris and Paul, 2002; Reeves et al., 2005), but is lower than some other reports for the same tendon, ranging between 486 and 2000 N mm−1 (Couppé et al., 2016; Kongsgaard et al., 2011; Magnusson et al., 2003, 2001). These apparent discrepancies across the literature can probably be explained by a number of factors including the proportion of force attributed to the Achilles tendon from the measured joint torque, differing sections of the force–elongation curve examined and variance in the anatomical structures tracked to measure longitudinal elongations.
Consistent with our findings in the ‘free’ Achilles tendon, a greater elongation of the gastrocnemius compared with the soleus tendon–aponeurosis has been observed using ultrasound imaging during maximal and sub-maximal contractions with the knee in full extension (Bojsen-Moller et al., 2004). Through the insertion and tracking of a needle into the ‘free’ Achilles tendon, it has been reported that the elongation and strain of the free Achilles tendon (corresponding to the distal portion of the Achilles tendon in the present study) was greater than an intra-muscular point tracked on the medial gastrocnemius tendon–aponeurosis (Magnusson et al., 2003). These previous findings are broadly in line with those of the present study when considering the two separate regions of the Achilles tendon; the present study found greater elongation and longitudinal strains in the distal Achilles tendon (equivalent to the ‘free’ tendon of the previous study) compared with the proximal Achilles tendon region (Table 3). Previous work using a 3D ultrasound approach found 3 mm elongation in the Achilles tendon and 3 mm elongation in the proximal gastrocnemius component of the Achilles tendon during contraction at 50% of MVC (Farris et al., 2013). Similar elongations between these two components of the Achilles tendon with markedly lower forces acting on the gastrocnemius tendon component are in line with our assertions of a lower stiffness for the proximal (gastrocnemius) compared with the distal Achilles tendon. Because of the greater length of the proximal Achilles tendon (gastrocnemius) compared with the distal region, these elongations correspond to smaller longitudinal strains in the proximal tendon component, consistent with MRI-based findings at 30% and 60% of MVC (Iwanuma et al., 2011). These previous reports, however, reflect overall ‘end-to-end’ length changes of these tendon components rather than more detailed region specificity as examined in the present study. An increased width measured at the gastrocnemius MTJ was noted upon contraction in a previous report (Farris et al., 2013). In another report, an increased width of the proximal (gastrocnemius) Achilles tendon component but a decreased width of the distal Achilles tendon component during contractions at 30% and 60% of MVC was noted (Iwanuma et al., 2011). Whilst girth measurements at a single site such as the MTJ may result from muscle bulging, tensile deformation of a tendon should only result in thinning of its overall CSA (Poisson’s ratio) and therefore consideration of only one axis in the transverse plane may not reveal the true nature of deformation during tendon loading. In the present study, the 6-fold higher Poisson's ratio of the proximal compared with the distal Achilles tendon region reflects the substantial ‘thinning’ of the proximal tendon region during tensile loading.
We initially hypothesised that larger CSA reductions, inferring larger longitudinal elongations, would occur in areas where the tendon CSA was smaller. In sharp contrast, we found larger CSA reductions, inferring larger longitudinal elongations, in the proximal Achilles tendon region where tendon CSA was actually greatest (Fig. 2). This finding raises at least two main possibilities. (1) Forces may not be distributed equally throughout the length of the Achilles tendon and areas of smaller CSA may experience lower forces and stresses. (2) There may be differences in the Achilles tendon's material properties along its length that could result from changes in the volume fraction, with collagen fibre-to-matrix ratios increasing in areas where the tendon CSA is smaller.
Although possible, option 1 seems unlikely as the forces experienced in the smaller CSA regions would need to be substantially lower than those in other regions to compensate for the smaller CSA. Although not implausible, it is difficult to understand how the forces transmitted through an in-series structure could vary so drastically. The density and area fraction of collagen fibrils have been shown to vary between different regions of the rabbit patellar tendon (Williams et al., 2008), raising the possibility that changes to the collagen fibre-to-matrix ratio (option 2) may occur along the length of the human Achilles tendon. If this collagen fibre-to-matrix ratio increases in regions of smaller tendon CSA, it may increase the modulus in these Achilles tendon regions and contribute towards explaining the current findings.
The finding of differential deformations along the length of the proximal and distal Achilles tendon components is in line with previous in vivo reports of non-uniform longitudinal elongations in the soleus tendon–aponeurosis and tibialis anterior muscle–tendon unit (Finni et al., 2003; Maganaris and Paul, 2000). Differences in elongation have also been reported between the gastrocnemius tendon–aponeurosis and the Achilles tendon (Bojsen-Moller et al., 2004; Magnusson et al., 2003). Our findings are also in line with in vitro reports of non-uniform elongations along the length of the human Achilles and other long tendons (Wren et al., 2001; Zernicke et al., 1984). As the gastrocnemius and soleus tendons ‘fuse’ to form the Achilles tendon distal to the soleus muscle, there is the potential for intra-tendinous shear to occur within the Achilles tendon. Indeed, the propensity for shear within the Achilles tendon has been indicated from in vivo (Bojsen-Moller et al., 2004; Magnusson et al., 2003) and in vitro (Lersch et al., 2012) human studies.
Three muscles (gastrocnemius medial and lateral heads and soleus) generate force applied to the Achilles tendon and their individual force contributions will vary according to their PCSA (Fukunaga et al., 1992). These three muscles may also be activated to different relative levels at any given joint torque. Indeed, activation of the lateral gastrocnemius muscle was shown to be relatively low at 30% of maximum voluntary plantarflexion contraction compared with that of the medial gastrocnemius (Masood et al., 2014). This highlights the complexity of loading on the Achilles tendon with contributions from three independent muscles, which have been suggested to constitute mechanically separate tendon compartments within the Achilles (Bojsen-Moller and Magnusson, 2015). Indeed, ultrasound speckle tracking methods have shown greater elongation within the deeper region of the Achilles tendon compared with the superficial layer during loading (Arndt et al., 2012; Franz et al., 2015; Slane and Thelen, 2015). This may suggest the presence of inter-fascicle sliding and a relative independence of gastrocnemius and soleus tendon components functioning within the Achilles tendon. The fibres of the Achilles tendon may spiral by up to 1.57 rad as they descend towards the attachment site on the calcaneus, with the degree of rotation varying according to the nature and extent of the fusion between gastrocnemius and soleus tendon components within the Achilles tendon (Cummins and Anson, 1946) and also varying considerably between individuals (Bojsen-Moller and Magnusson, 2015). This complex 3D micro-structure will have implications for the degree of deformation occurring along the different regions of the Achilles tendon, probably determining areas of shear and torsional stress concentration. This is reflected by markedly different values for Poisson's ratio between the proximal and distal Achilles tendon regions (Table 3). It might be speculated that rather than differences in tendon CSA as we initially hypothesised, the complex 3D micro-structure and the potential for shear within the tendon could be factors contributing to the high rupture rate of the Achilles tendon 3–5 cm proximal to the calcaneus (Józsa et al., 1989). To determine the effects of different structural distributions within and along the Achilles tendon, construction of a finite element model is required similar to that previously performed for the patellar tendon (Lavagnino et al., 2008).
The approach followed in the present study involved isometric contractions being held for 2 min duration during MRI scanning. During prolonged constant tensile loading, there is the potential for tendon creep to occur and we therefore tested for this possibility using ultrasound scanning. Our results showed no significant increase in longitudinal tendon elongation during this constant contraction over a 3 min period (i.e. longer than our actual MRI scanning period), indicating that creep was unlikely to be present in our MRI measurements for the range of loading levels examined (Fig. 3). The absence of creep may be explained by the lower level of contractions elicited here (up to 30% MVC) compared with other studies where creep has been observed (Cohen et al., 1976; Maganaris et al., 2002).
Co-activation of dorsiflexor muscles was also unlikely to play any role in our MRI measurements at the range of voluntary forces examined, as the estimated dorsiflexor co-activation torque was <1%. Our participants demonstrated on the dynamometer that they were capable of maintaining the target torque very accurately during plantarflexion contractions with visual feedback at the elicited torque levels (Table 1). We are therefore confident that participants maintained the torque level constant in the MRI scanner during the acquisition of the scans. Additionally, our plantarflexion lever system was very sensitive and the experimenter present in the MRI scanner could identify whether any deviations occurred from the balanced situation.
Certain limitations and assumptions of the current work should be noted. In the present study, we examined plantarflexion contractions up to the level of 30% MVC because of the constraints associated with the length of time contractions needed to be maintained in the MRI scanner (2 min); therefore, caution should be emphasised in extrapolating the current findings to higher forces approaching MVC. Although it might be argued that these force levels (10%, 20% and 30% MVC) are at the relatively low end in relation to MVC forces, the largest elongations occur within the initial low-force region of the tendon force–elongation curve (Butler et al., 1978). In fact, at force levels corresponding to 30% of MVC, tendon elongation is ∼50–60% of the elongation measured at 100% of MVC (Malliaras et al., 2013; Reeves et al., 2003, 2005). Hence, by examining tendon deformations at loads up to 30% of MVC, we actually cover a force range accounting for 50–60% of its entire elongation. Because of constraints of the MRI field of view, we were not able to scan the entire length of the proximal part of the Achilles tendon, which should be acknowledged as a limitation. For calculation of the mechanical properties of this proximal part of the tendon, it was therefore necessary to estimate its length using an anthropometric ratio from published data. When estimating the overall elongation of the distal part of the Achilles tendon in calculating mechanical properties, we applied a ratio from published data to the elongation directly measured higher up at the gastrocnemius MTJ. In calculating stiffness of the distal and proximal portions of the Achilles tendon, although the overall force–elongation relationship was curvilinear, a necessary simplification made in the present study (and all previous studies of a similar nature) was to assume linearity of specific but small sections of this force–elongation curve.
In conclusion, using a novel MRI-based approach, the present study shows marked differences along the length of the Achilles tendon, with larger transverse tendon deformations upon tensile loading within the proximal compared with the distal region of the Achilles tendon. These marked differences were reflected by the force on the proximal tendon region being ∼60% lower, stress ∼70% lower, stiffness ∼30% lower and Poisson's ratio 6-fold higher compared with the distal Achilles tendon component. These results suggest that the proximal component of the Achilles tendon may act as a mechanical buffer to protect the stiffer, more highly stressed, distal component of the Achilles tendon from injury.
Acknowledgements
Many thanks to Diane Chrapkowski for her role in the data collection for this study.
Footnotes
Author contributions
N.D.R. conceived the project idea. N.D.R. and G.C. were both involved in collecting the data for the study. N.D.R. conducted all image analysis. N.D.R. primarily drafted the manuscript with both N.D.R. and G.C. reviewing manuscript drafts and providing final approval.
Funding
This work was supported by funding from the Engineering and Physical Sciences Research Council (EPSRC), via the ‘Nano-info-bio’ project.
References
Competing interests
The authors declare no competing or financial interests.