Correction: A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait.

[This corrects the article DOI: 10.1371/journal.pcbi.1009398.].

Morphology of proximal and distal human semitendinosus compartments and the effects of distal tendon harvesting for anterior cruciate ligament reconstruction.

The human semitendinosus muscle is characterized by a tendinous inscription separating proximal and distal neuromuscular compartments. As each compartment is innervated by separate nerve branches, potential exists for independent operation and control of compartments. However, the morphology and function of each compartment have not been thoroughly examined in an adult human population. Further, the distal semitendinosus tendon is typically harvested for use in anterior cruciate ligament reconstruction surgery, which induces long-term morphological changes to the semitendinosus muscle-tendon unit. It remains unknown if muscle morphological alterations following anterior cruciate ligament reconstruction are uniform between proximal and distal semitendinosus compartments. Here, we performed magnetic resonance imaging on 10 individuals who had undergone anterior cruciate ligament reconstruction involving an ipsilateral distal semitendinosus tendon graft 14 ± 4 months prior, extracting morphological parameters of the whole semitendinosus muscle and each individual compartment from both the (non-injured) contralateral and surgical legs. In the contralateral leg, volume and length of the proximal compartment were smaller than the distal compartment. No between-compartment differences in volume or length were found for anterior cruciate ligament reconstructed legs, likely due to greater shortening of the distal compared to the proximal compartment after anterior cruciate ligament reconstruction. The maximal anatomical cross-sectional area of both compartments was substantially smaller on the anterior cruciate ligament reconstructed leg but did not differ between compartments on either leg. The absolute and relative between-leg differences in proximal compartment morphology on the anterior cruciate ligament reconstructed leg were strongly correlated with the corresponding between-leg differences in distal compartment morphological parameters. Specifically, greater between-leg morphological differences in one compartment were highly correlated with large between-leg differences in the other compartment, and vice versa for smaller differences. These relationships indicate that despite the heterogeneity in compartment length and volume, compartment atrophy is not independent or random. Further, the tendinous inscription endpoints were generally positioned at the same proximodistal level as the compartment maximal anatomical cross-sectional areas, providing a wide area over which the tendinous inscription could mechanically interact with compartments. Overall, results suggest the two human semitendinosus compartments are not mechanically independent.

A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait.

Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.

Alterations in the Functional Knee Alignment Are Not an Effective Strategy to Modify the Mediolateral Distribution of Knee Forces During Closed Kinetic Chain Exercises.

Pain felt while performing rehabilitation exercises could be a reason for the low adherence of knee osteoarthritis patients to physical rehabilitation. Reducing compressive forces on the most affected knee regions may help to mitigate the pain. Knee frontal plane positioning with respect to pelvis and foot (functional knee alignment) has been shown to modify the mediolateral distribution of the tibiofemoral joint contact force in walking. Hence, different functional knee alignments could be potentially used to modify joint loading during rehabilitation exercises. The aim was to understand whether utilizing different alignments is an effective strategy to unload specific knee areas while performing rehabilitation exercises. Eight healthy volunteers performed 5 exercises with neutral, medial, and lateral knee alignment. A musculoskeletal model was modified for improved prediction of tibiofemoral contact forces and used to evaluate knee joint kinematics, moments, and contact forces. Functional knee alignment had only a small and inconsistent effect on the mediolateral distribution joint contact force. Moreover, the magnitude of tibiofemoral and patellofemoral contact forces, knee moments, and measured muscle activities was not significantly affected by the alignment. Our results suggest that altering the functional knee alignment is not an effective strategy to unload specific knee regions in physical rehabilitation.

Non-uniform displacement within ruptured Achilles tendon during isometric contraction.

The purpose of this study was investigate tendon displacement patterns in non-surgically treated patients 14 months after acute Achilles tendon rupture (ATR) and to classify patients into groups based on their Achilles tendon (AT) displacement patterns. Twenty patients were tested. Sagittal images of AT were acquired using B-mode ultrasonography during ramp contractions at a torque level corresponding to 30% of the maximal isometric plantarflexion torque of the uninjured limb. A speckle tracking algorithm was used to track proximal-distal movement of the tendon tissue at 6 antero-posterior locations. Two-way repeated measures ANOVA for peak tendon displacement was performed. K-means clustering was used to classify patients according to AT displacement patterns. The difference in peak relative displacement across locations was larger in the uninjured (1.29 ± 0.87 mm) than the injured limb (0.69 ± 0.68 mm), with a mean difference (95% CI) of 0.60 mm (0.14-1.05 mm, P < .001) between limbs. For the uninjured limb, cluster analysis formed 3 groups, while 2 groups were formed for the injured limb. The three distinct patterns of AT displacement during isometric plantarflexion in the uninjured limb may arise from subject-specific anatomical variations of AT sub-tendons, while the two patterns in the injured limb may reflect differential recovery after ATR with non-surgical treatment. Subject-specific tendon characteristics are a vital determinant of stress distribution across the tendon. Changes in stress distribution may lead to variation in the location and magnitude of peak displacement within the free AT. Quantifying internal tendon displacement patterns after ATR provides new insights into AT recovery.

Individual Leg Muscle Contributions to the Cost of Walking: Effects of Age and Walking Speed.

This study examined the contributions of individual muscles to changes in energetic cost of transport (COT) over seven walking speeds, and compared results between healthy young and elderly subjects. Twenty six participants (13 young aged 18-30; 13 old aged 70-80) were recruited. COT (O2/kg body mass/km) was calculated by standardizing the mean oxygen consumption recorded during steady state walking. Electromyography signals from 10 leg muscles were used to calculate the cumulative activity required to traverse a unit of distance (CMAPD) for each muscle at each speed. In the old group CMAPD was correlated with COT, presented higher and more variable values, and showed greater increases around optimal speed for all studied muscles. Soleus CMAPD was independent of speed in the young group, but this was not evident with aging. Greater energy cost of walking in older individuals seems to be attributable to increased energy cost of all lower limb muscles.

Association between osteocalcin and cognitive performance in healthy older adults.

INTRODUCTION: cognitive deterioration and reductions of bone health coincide with increasing age. We examine the relationship between bone composition and plasma markers of bone remodelling with measures of cognitive performance in healthy adults. METHODS: this cross-sectional study included 225 old (52% women, mean age: 74.4 ± 3.3 years) and 134 young (52% women, mean age: 23.4 ± 2.7 years) adult participants from the MyoAge project. Whole body bone mineral density was measured by dual-energy X-ray absorptiometry. Blood analyses included a panel of bone-related peptides (dickkopf-1, osteoprotegerin, osteocalcin (OC), osteopontin, sclerostin, parathyroid hormone and fibroblast growth factor 23), as well as serum calcium and 25-hydroxy vitamin D assays. A selection of cognitive domains (working memory capacity, episodic memory, executive functioning and global cognition) was assessed with a standardised neuropsychological test battery. RESULTS: adjusting for covariates and multiple testing revealed that plasma OC levels were positively associated with measures of executive functioning (ß = 0.444, P < 0.001) and global cognition (ß = 0.381, P = 0.001) in the older women. DISCUSSION: these correlative results demonstrate a positive association between OC, a factor known to regulate bone remodelling, with cognitive performance in older non-demented women. Further work should address possible mechanistic interpretations in humans.

Gait data from 51 healthy participants with motion capture, inertial measurement units, and computer vision.

We present a dataset comprising motion capture, inertial measurement unit data, and sagittal-plane video data from walking at three different instructed speeds (slow, comfortable, fast). The dataset contains 51 healthy participants with approximately 60 walking trials from each participant. Each walking trial contains data from motion capture, inertial measurement units, and computer vision. Motion capture data comprises ground reaction forces and moments from floor-embedded force plates and the 3D trajectories of subject-worn motion capture markers. Inertial measurement unit data comprises 3D accelerometer readings and 3D orientations from the lower limbs and pelvis. Computer vision data comprises 2D keypoint trajectories detected using the OpenPose human pose estimation algorithm from sagittal-plane video of the walking trial. Additionally, the dataset contains participant demographic and anthropometric information such as mass, height, sex, age, lower limb dimensions, and knee intercondylar distance measured from magnetic resonance images. The dataset can be used in musculoskeletal modelling and simulation to calculate kinematics and kinetics of motion and to compare data between motion capture, inertial measurement, and video capture.

Effect of patient specificity on predicting knee cartilage degeneration in obese adults: Musculoskeletal finite-element modeling of data from the CAROT trial.

Obesity is a known risk factor for development of osteoarthritis (OA). Numerical tools like finite-element (FE) models combined with degenerative algorithms have been developed to understand the interplay between OA and obesity. In this study, we aimed to predict knee cartilage degeneration in a cohort of obese adults to investigate the importance of patient-specific information on degeneration predictions. We used a validated FE modeling approach and three different age-dependent functions (step-wise, exponential, and linear) to simulate cartilage degradation under overloading in the knee joint. Gait motion analysis and magnetic resonance imaging data from 115 obese individuals with knee OA were used for musculoskeletal and FE modeling. Cartilage degeneration predictions were contrasted with Kellgren-Lawrence (KL) and Boston-Leeds Osteoarthritis Knee Score (BLOKS) grades. The findings show that overall, the similarities between numerical predictions and clinical measures were better for the medial (average area under the curve (AUC) = 0.62) compared to the lateral compartment (average AUC = 0.52) of the knee. Classification results for KL grades, full patient-specific models and patient-specific geometry with generic gait data showed higher AUC values (AUC = 0.71 and AUC = 0.68, respectively) compared to generic geometry and patient-specific gait (AUC = 0.48). For BLOKS grades, AUC values for both full patient-specific models and for patient-specific geometry with generic gait locomotion were higher (AUC = 0.66 and AUC = 0.64, respectively) compared to when the generic geometry and patient-specific gait were used (AUC = 0.53). In summary, our study highlights the importance of considering individual information in knee OA prediction. Nevertheless, our findings suggest that personalized gait play a smaller role in the OA prediction and classification capacity than personalized joint geometry.

Predicting Knee Joint Contact Force Peaks During Gait Using a Video Camera or Wearable Sensors.

PURPOSE: Estimating loading of the knee joint may be helpful in managing degenerative joint diseases. Contemporary methods to estimate loading involve calculating knee joint contact forces using musculoskeletal modeling and simulation from motion capture (MOCAP) data, which must be collected in a specialized environment and analyzed by a trained expert. To make the estimation of knee joint loading more accessible, simple input predictors should be used for predicting knee joint loading using artificial neural networks. METHODS: We trained feedforward artificial neural networks (ANNs) to predict knee joint loading peaks from the mass, height, age, sex, walking speed, and knee flexion angle (KFA) of subjects using their existing MOCAP data. We also collected an independent MOCAP dataset while recording walking with a video camera (VC) and inertial measurement units (IMUs). We quantified the prediction accuracy of the ANNs using walking speed and KFA estimates from (1) MOCAP data, (2) VC data, and (3) IMU data separately (i.e., we quantified three sets of prediction accuracy metrics). RESULTS: Using portable modalities, we achieved prediction accuracies between 0.13 and 0.37 root mean square error normalized to the mean of the musculoskeletal analysis-based reference values. The correlation between the predicted and reference loading peaks varied between 0.65 and 0.91. This was comparable to the prediction accuracies obtained when obtaining predictors from motion capture data. DISCUSSION: The prediction results show that both VCs and IMUs can be used to estimate predictors that can be used in estimating knee joint loading outside the motion laboratory. Future studies should investigate the usability of the methods in an out-of-laboratory setting.

Muscle Morphology Does Not Solely Determine Knee Flexion Weakness After Anterior Cruciate Ligament Reconstruction with a Semitendinosus Tendon Graft: A Combined Experimental and Computational Modeling Study.

The distal semitendinosus tendon is commonly harvested for anterior cruciate ligament reconstruction, inducing substantial morbidity at the knee. The aim of this study was to probe how morphological changes of the semitendinosus muscle after harvest of its distal tendon for anterior cruciate ligament reconstruction affects knee flexion strength and whether the knee flexor synergists can compensate for the knee flexion weakness. Ten participants 8-18 months after anterior cruciate ligament reconstruction with an ipsilateral distal semitendinosus tendon autograft performed isometric knee flexion strength testing (15°, 45°, 60°, and 90°; 0° = knee extension) positioned prone on an isokinetic dynamometer. Morphological parameters extracted from magnetic resonance images were used to inform a musculoskeletal model. Knee flexion moments estimated by the model were then compared with those measured experimentally at each knee angle position. A statistically significant between-leg difference in experimentally-measured maximal isometric strength was found at 60° and 90°, but not 15° or 45°, of knee flexion. The musculoskeletal model matched the between-leg differences observed in experimental knee flexion moments at 15° and 45° but did not well estimate between-leg differences with a more flexed knee, particularly at 90°. Further, the knee flexor synergists could not physiologically compensate for weakness in deep knee flexion. These results suggest additional factors other than knee flexor muscle morphology play a role in knee flexion weakness following anterior cruciate ligament reconstruction with a distal semitendinosus tendon graft and thus more work at neural and microscopic levels is required for informing treatment and rehabilitation in this demographic.

Physiological and functional evaluation of healthy young and older men and women: design of the European MyoAge study.

Within the European multi-centre MyoAge project, one workpackage was designed to investigate the contribution of age-related changes to muscle mass, contractile characteristics and neural control in relation to reductions in mobility in older age. The methodology has been described here. Test centres were located in Manchester, UK; Paris, France; Leiden, The Netherlands; Tartu, Estonia and Jyväskylä, Finland. In total, 182 young (18-30 years old, 52.2 % female) and 322 older adults (69-81 years old, 50 % female) have been examined. The participants were independent living, socially active and free from disease that impaired mobility levels. The older participants were selected based on physical activity levels, such that half exceeded current recommended physical activity levels and the other half had lower physical activity levels than is recommended to maintain health. Measurements consisted of blood pressure; anthropometry and body composition (dual-energy X-ray absorptiometry and magnetic resonance imaging); lung function; standing balance and cognitive function (CANTAB). Mobility was assessed using the Timed Up and Go, a 6 min walk, activity questionnaires and accelerometers to monitor habitual daily activities. Muscle strength, power, fatigue and neural activation were assessed using a combination of voluntary and electrically stimulated contractions. Fasting blood samples and skeletal muscle biopsies were collected for detailed examination of cell and molecular differences between young and older individuals. The results from this study will provide a detailed insight into "normal, healthy" ageing, linking whole-body function to the structure and function of the neuromuscular system and the molecular characteristics of skeletal muscle.

Circulating levels of adipokines and IGF-1 are associated with skeletal muscle strength of young and old healthy subjects.

It is known that adipose tissue mass increases with age, and that a number of hormones, collectively called adipokines, are produced by adipose tissue. For most of them it is not known whether their plasmatic levels change with age. Moreover, it is known that adipose tissue infiltration in skeletal muscle is related to sarcopenia and loss of muscle strength. In this study we investigated the age-related changes of representative adipokines and insulin-like growth factor (IGF)-1 and their effect on muscle strength. We studied the association between circulating levels of adiponectin, leptin, resistin and IGF-1 and muscle strength. This cross-sectional study included 412 subjects of different age (152 subjects aged 18-30 years and 260 subjects aged 69-81 years) recruited within the framework of the European research network project "Myoage". The levels of adiponectin (both in male and female subjects) and leptin (only in males) were significantly higher in old subjects compared to young, while those of IGF-1 were lower in old subjects. In old subjects adiponectin, resistin and the resistin/IGF-1 ratio (but not IGF-1 alone) were inversely associated with quadriceps torque, while only adiponectin was inversely associated with handgrip strength independently from percentage of fat mass, height, age, gender and geographical origin. The ratio of leptin to adiponectin was directly associated with handgrip strength in both young and old subjects. These results suggest that in humans the age-associated loss of strength is associated with the levels of representative adipokines and IGF-1.

Force-velocity profiling in ice hockey skating: reliability and validity of a simple, low-cost field method.

In recent years, a simple method for force-velocity (F-v) profiling, based on split times, has emerged as a potential tool to examine mechanical variables underlying running sprint performance in field conditions. In this study, the reliability and concurrent validity of F-v profiling based on split times were examined when used for ice hockey skating. It was also tested how a modification of the method, in which the start instant of the sprint is estimated based on optimisation (time shift method), affects the reliability and validity of the method. Both intra- and inter-rater reliability were markedly improved when using the time shift method (approximately 50% decrease in the standard error of measurement). Moreover, the results calculated using the time shift method highly correlated (r > 0.83 for all variables) with the results calculated from a continuously tracked movement of the athlete, which was considered here as the reference method. This study shows that a modification to the previously published simple method for F-v profiling improves intra- and inter-rater reliability of the method in ice hockey skating. The time shift method tested here can be used as a reliable tool to test a player's physical performance characteristic underlying sprint performance in ice hockey skating.

Open-source software library for real-time inertial measurement unit data-based inverse kinematics using OpenSim.

Background: Inertial measurements (IMUs) facilitate the measurement of human motion outside the motion laboratory. A commonly used open-source software for musculoskeletal simulation and analysis of human motion, OpenSim, includes a tool to enable kinematics analysis of IMU data. However, it only enables offline analysis, i.e., analysis after the data has been collected. Extending OpenSim's functionality to allow real-time kinematics analysis would allow real-time feedback for the subject during the measurement session and has uses in e.g., rehabilitation, robotics, and ergonomics. Methods: We developed an open-source software library for real-time inverse kinematics (IK) analysis of IMU data using OpenSim. The software library reads data from IMUs and uses multithreading for concurrent calculation of IK. Its operation delays and throughputs were measured with a varying number of IMUs and parallel computing IK threads using two different musculoskeletal models, one a lower-body and torso model and the other a full-body model. We published the code under an open-source license on GitHub. Results: A standard desktop computer calculated full-body inverse kinematics from treadmill walking at 1.5 m/s with data from 12 IMUs in real-time with a mean delay below 55 ms and reached a throughput of more than 90 samples per second. A laptop computer had similar delays and reached a throughput above 60 samples per second with treadmill walking. Minimal walking kinematics, motion of lower extremities and torso, were calculated from treadmill walking data in real-time with a throughput of 130 samples per second on the laptop and 180 samples per second on the desktop computer, with approximately half the delay of full-body kinematics. Conclusions: The software library enabled real-time inverse kinematical analysis with different numbers of IMUs and customizable musculoskeletal models. The performance results show that subject-specific full-body motion analysis is feasible in real-time, while a laptop computer and IMUs allowed the use of the method outside the motion laboratory.

Prediction of Knee Joint Compartmental Loading Maxima Utilizing Simple Subject Characteristics and Neural Networks.

Joint loading may affect the development of osteoarthritis, but patient-specific load estimation requires cumbersome motion laboratory equipment. This reliance could be eliminated using artificial neural networks (ANNs) to predict loading from simple input predictors. We used subject-specific musculoskeletal simulations to estimate knee joint contact forces for 290 subjects during over 5000 stance phases of walking and then extracted compartmental and total joint loading maxima from the first and second peaks of the stance phase. We then trained ANN models to predict the loading maxima from predictors that can be measured without motion laboratory equipment (subject mass, height, age, gender, knee abduction-adduction angle, and walking speed). When compared to the target data, our trained models had NRMSEs (RMSEs normalized to the mean of the response variable) between 0.14 and 0.42 and Pearson correlation coefficients between 0.42 and 0.84. The loading maxima were predicted most accurately using the models trained with all predictors. We demonstrated that prediction of knee joint loading maxima may be possible without laboratory-measured motion capture data. This is a promising step in facilitating knee joint loading predictions in simple environments, such as a physician's appointment. In future, the rapid measurement and analysis setup could be utilized to guide patients in rehabilitation to slow development of joint disorders, such as osteoarthritis.

Exploration of muscle-tendon biomechanics one year after Achilles tendon rupture and the compensatory role of flexor hallucis longus.

Achilles tendon (AT) rupture leads to long-term structural and functional impairments. Currently, the predictors of good recovery after rupture are poorly known. Thus, we aimed to explore the interconnections between structural, mechanical, and neuromuscular parameters and their associations with factors that could explain good recovery in patients with non-surgically treated AT rupture. A total of 35 patients with unilateral rupture (6 females) participated in this study. Muscle-tendon structural, mechanical, and neuromuscular parameters were measured 1-year after rupture. Interconnections between the inter-limb differences (Δ) were explored using partial correlations, followed by multivariable linear regression to find associations between the measured factors and the following markers that indicate good recovery: 1) tendon length, 2) tendon non-uniform displacement, and 3) flexor hallucis longus (FHL) normalized EMG amplitude difference between limbs. Δmedial gastrocnemius (MG) (ß = -0.12, p = 0.007) and Δlateral gastrocnemius (ß = -0.086, p = 0.030) subtendon lengths were associated with MG tendon Δstiffness. MG (ß = 11.56, p = 0.003) and soleus (ß = 2.18, p = 0.040) Δsubtendon lengths explained 48 % of variance in FHL EMG amplitude. Regression models for tendon length and non-uniform displacement were not significant. Smaller inter-limb differences in Achilles subtendon lengths were associated with smaller differences in the AT stiffness between limbs, and a smaller contribution of FHL muscle to the plantarflexion torque. In the injured limb, the increased contribution of FHL appears to partially counteract a smaller contribution from MG due to the elongated tendon, however the role of FHL should not be emphasized during rehabilitation to allow recovery of the TS muscles.

Sensitivity of simulated knee joint mechanics to selected human and bovine fibril-reinforced poroelastic material properties.

Fibril-reinforced poroviscoelastic material models are considered state-of-the-art in modeling articular cartilage biomechanics. Yet, cartilage material parameters are often based on bovine tissue properties in computational knee joint models, although bovine properties are distinctly different from those of humans. Thus, we aimed to investigate how cartilage mechanical responses are affected in the knee joint model during walking when fibril-reinforced poroviscoelastic properties of cartilage are based on human data instead of bovine. We constructed a finite element knee joint model in which tibial and femoral cartilages were modeled as fibril-reinforced poroviscoelastic material using either human or bovine data. Joint loading was based on subject-specific gait data. The resulting mechanical responses of knee cartilage were compared between the knee joint models with human or bovine fibril-reinforced poroviscoelastic cartilage properties. Furthermore, we conducted a sensitivity analysis to determine which fibril-reinforced poroviscoelastic material parameters have the greatest impact on cartilage mechanical responses in the knee joint during walking. In general, bovine cartilage properties yielded greater maximum principal stresses and fluid pressures (both up to 30%) when compared to the human cartilage properties during the loading response in both femoral and tibial cartilage sites. Cartilage mechanical responses were very sensitive to the collagen fibril-related material parameter variations during walking while they were unresponsive to proteoglycan matrix or fluid flow-related material parameter variations. Taken together, human cartilage material properties should be accounted for when the goal is to compare absolute mechanical responses of knee joint cartilage as bovine material parameters lead to substantially different cartilage mechanical responses.

Achilles tendon stiffness is unchanged one hour after a marathon.

Overuse-induced injuries have been proposed as a predisposing factor for Achilles tendon (AT) ruptures. If tendons can be overloaded, their mechanical properties should change during exercise. Because there data are lacking on the effects of a single bout of long-lasting exercise on AT mechanical properties, the present study measured AT stiffness before and after a marathon. AT stiffness was determined as the slope of the force-elongation curve between 10 and 80% of maximum voluntary force. AT force-elongation characteristics were measured in an ankle dynamometer using simultaneous motion-capture-assisted ultrasonography. Oxygen consumption and ankle kinematics were also measured on a treadmill at the marathon pace. All measurements were performed before and after the marathon. AT stiffness did not change significantly from the pre-race value of 197±62 N mm(-1) (mean ± s.d.) to the post-race value of 206±59 N mm(-1) (N=12, P=0.312). Oxygen consumption increased after the race by 7±10% (P<0.05) and ankle kinematic data revealed that in nine out of 12 subjects, the marathon induced a change in their foot strike technique. The AT of the physically active individuals seems to be able to resist mechanical changes under physiological stress. We therefore suggest that natural loading, like in running, may not overstress the AT or predispose it to injury. In addition, decreased running economy, as well as altered foot strike technique, was probably attributable to muscle fatigue.

In vivo localized gastrocnemius subtendon representation within the healthy and ruptured human Achilles tendon.

The Achilles tendon (AT) is composed of three distinct in-series elastic subtendons, arising from different muscles in the triceps surae. Independent activation of any of these muscles is thought to induce sliding between the adjacent AT subtendons. We aimed to investigate displacement patterns during voluntary contraction (VOL) and selective transcutaneous stimulation of medial (MGstim) and lateral (LGstim) gastrocnemius between ruptured and healthy tendons and to examine the representative areas of AT subtendons. Twenty-eight patients with unilateral AT rupture performed bilateral VOL at 30% of the maximal isometric uninjured plantarflexion torque. AT displacement was analyzed from sagittal B-mode ultrasonography images during VOL, MGstim, and LGstim. Three-way ANOVA revealed a significant two-way interaction of contraction type × location on the tendon displacement [F(10-815) = 3.72, P < 0.001]. The subsequent two-way analysis revealed a significant contraction type × location interaction for tendon displacement [F(10-410)=3.79, P < 0.001] in the uninjured limb only, where LGstim displacement pattern was significantly different from MGstim (P = 0.008) and VOL (P = 0.005). When comparing contraction types between limbs, there were no difference in the displacement patterns, but displacement amplitudes differed. There was no significant difference in the location of maximum or minimum displacement between limbs. The displacement pattern was not different in nonsurgically treated compared with uninjured tendons 1 yr after rupture. However, free tendon stiffness seems to be lower in the injured AT, leading to more displacement during electrically induced contractions compared with the uninjured. Our results suggest that near the calcaneus, LG subtendon is located in the most anterior region adjacent to medial gastrocnemius.NEW & NOTEWORTHY Using selective electrical stimulation, we report the distributions of medial and lateral gastrocnemius subtendon representations within the healthy and ruptured Achilles tendon. In the majority of our sample, lateral gastrocnemius subtendon was found in the most anterior region adjacent to medial gastrocnemius both in the healthy and ruptured, nonsurgically treated tendon. The tendon internal displacement pattern does not seem to differ, but displacement amplitude and nonuniformity differed between healthy and ruptured tendons 1 yr after rupture.