Skeletal muscle function underpins muscle spindle abundance.

Muscle spindle abundance is highly variable within and across species, but we currently lack any clear picture of the mechanistic causes or consequences of this variation. Previous use of spindle abundance as a correlate for muscle function implies a mechanical underpinning to this variation, but these ideas have not been tested. Herein, we use integrated medical imaging and subject-specific musculoskeletal models to investigate the relationship between spindle abundance, muscle architecture and in vivo muscle behaviour in the human locomotor system. These analyses indicate that muscle spindle number is tightly correlated with muscle fascicle length, absolute fascicle length change, velocity of fibre lengthening and active muscle forces during walking. Novel correlations between functional indices and spindle abundance are also recovered, where muscles with a high abundance predominantly function as springs, compared to those with a lower abundance mostly functioning as brakes during walking. These data demonstrate that muscle fibre length, lengthening velocity and fibre force are key physiological signals to the central nervous system and its modulation of locomotion, and that muscle spindle abundance may be tightly correlated to how a muscle generates work. These insights may be combined with neuromechanics and robotic studies of motor control to help further tease apart the functional drivers of muscle spindle composition.

Foot anatomy, walking energetics, and the evolution of human bipedalism.

Interspecies differences in locomotor efficiency have been extensively researched, but within-species variation in the metabolic cost of walking and its underlying causes have received much less attention. This is somewhat surprising given the importance of walking energetics to natural selection, and the fact that the mechanical efficiency of striding bipedalism in modern humans is thought to be related in some part to the unique morphology of the human foot. Previous studies of human running have linked specific anatomical traits in the foot to variations in locomotor energetics to provide insight into form-function relationships in human evolution. However, such studies are relatively rare, particularly for walking. In this study, relationships between a range of functional musculoskeletal traits in the human lower limb and the energetics of walking over compliant and noncompliant substrates are examined, with particular focus on the lower limb and foot. Twenty-nine young, healthy individuals walked across three surfaces-a noncompliant laboratory floor, and compliant 6 cm and 13 cm thick foams-at self-selected speeds while oxygen consumption was measured, from which the metabolic cost of transport was calculated. Lower limb lengths, calcaneus lengths, foot shape indices, and maximum isometric plantarflexion torques were also measured and subsequently tested for relationships with metabolic cost over these surfaces using linear regression. It was found that metabolic cost varied considerably between individuals within and across substrate types, but this variation was not statistically related to or explained by variations in musculoskeletal parameters considered to be adaptively important to efficient bipedal locomotion. This therefore provides no supportive evidence that variations in these gross anatomical parameters confer significant advantages to the efficiency of walking, and therefore suggest caution in the use of similar metrics to infer differences in walking energetics in closely related fossil species.

Predictions of Anterior Cruciate Ligament Dynamics From Subject-Specific Musculoskeletal Models and Dynamic Biplane Radiography.

In vivo knee ligament forces are important to consider for informing rehabilitation or clinical interventions. However, they are difficult to directly measure during functional activities. Musculoskeletal models and simulations have become the primary methods by which to estimate in vivo ligament loading. Previous estimates of anterior cruciate ligament (ACL) forces range widely, suggesting that individualized anatomy may have an impact on these predictions. Using ten subject-specific (SS) lower limb musculoskeletal models, which include individualized musculoskeletal geometry, muscle architecture, and six degree-of-freedom knee joint kinematics from dynamic biplane radiography (DBR), this study provides SS estimates of ACL force (anteromedial-aACL; and posterolateral-pACL bundles) during the full gait cycle of treadmill walking. These forces are compared to estimates from scaled-generic (SG) musculoskeletal models to assess the effect of musculoskeletal knee joint anatomy on predicted forces and the benefit of SS modeling in this context. On average, the SS models demonstrated a double force peak during stance (0.39-0.43 xBW per bundle), while only a single force peak during stance was observed in the SG aACL. No significant differences were observed between continuous SG and SS ACL forces; however, root mean-squared differences between SS and SG predictions ranged from 0.08 xBW to 0.27 xBW, suggesting SG models do not reliably reflect forces predicted by SS models. Force predictions were also found to be highly sensitive to ligament resting length, with ±10% variations resulting in force differences of up to 84%. Overall, this study demonstrates the sensitivity of ACL force predictions to SS anatomy, specifically musculoskeletal joint geometry and ligament resting lengths, as well as the feasibility for generating SS musculoskeletal models for a group of subjects to predict in vivo tissue loading during functional activities.

Subject-specific muscle properties from diffusion tensor imaging significantly improve the accuracy of musculoskeletal models.

Musculoskeletal modelling is an important platform on which to study the biomechanics of morphological structures in vertebrates and is widely used in clinical, zoological and palaeontological fields. The popularity of this approach stems from the potential to non-invasively quantify biologically important but difficult-to-measure functional parameters. However, while it is known that model predictions are highly sensitive to input values, it is standard practice to build models by combining musculoskeletal data from different sources resulting in 'generic' models for a given species. At present, there are little quantitative data on how merging disparate anatomical data in models impacts the accuracy of these functional predictions. This issue is addressed herein by quantifying the accuracy of both subject-specific human limb models containing individualised muscle force-generating properties and models built using generic properties from both elderly and young individuals, relative to experimental muscle torques obtained from an isokinetic dynamometer. The results show that subject-specific models predict isokinetic muscle torques to a greater degree of accuracy than generic models at the ankle (root-mean-squared error - 7.9% vs. 49.3% in elderly anatomy-based models), knee (13.2% vs. 57.3%) and hip (21.9% vs. 32.8%). These results have important implications for the choice of musculoskeletal properties in future modelling studies, and the relatively high level of accuracy achieved in the subject-specific models suggests that such models can potentially address questions about inter-subject variations of muscle functions. However, despite relatively high levels of overall accuracy, models built using averaged generic muscle architecture data from young, healthy individuals may lack the resolution and accuracy required to study such differences between individuals, at least in certain circumstances. The results do not wholly discourage the continued use of averaged generic data in musculoskeletal modelling studies but do emphasise the need for to maximise the accuracy of input values if studying intra-species form-function relationships in the musculoskeletal system.

Determining Subject-Specific Lower-Limb Muscle Architecture Data for Musculoskeletal Models Using Diffusion Tensor Imaging.

Accurate individualized muscle architecture data are crucial for generating subject-specific musculoskeletal models to investigate movement and dynamic muscle function. Diffusion tensor imaging (DTI) magnetic resonance (MR) imaging has emerged as a promising method of gathering muscle architecture data in vivo; however, its accuracy in estimating parameters such as muscle fiber lengths for creating subject-specific musculoskeletal models has not been tested. Here, we provide a validation of the method of using anatomical magnetic resonance imaging (MRI) and DTI to gather muscle architecture data in vivo by directly comparing those data obtained from MR scans of three human cadaveric lower limbs to those from dissections. DTI was used to measure fiber lengths and pennation angles, while the anatomical images were used to estimate muscle mass, which were used to calculate physiological cross-sectional area (PCSA). The same data were then obtained through dissections, where it was found that on average muscle masses and fiber lengths matched well between the two methods (4% and 1% differences, respectively), while PCSA values had slightly larger differences (6%). Overall, these results suggest that DTI is a promising technique to gather in vivo muscle architecture data, but further refinement and complementary imaging techniques may be needed to realize these goals.

Muscle moment arms and sensitivity analysis of a mouse hindlimb musculoskeletal model.

Musculoskeletal modelling has become a valuable tool with which to understand how neural, muscular, skeletal and other tissues are integrated to produce movement. Most musculoskeletal modelling work has to date focused on humans or their close relatives, with few examples of quadrupedal animal limb models. A musculoskeletal model of the mouse hindlimb could have broad utility for questions in medicine, genetics, locomotion and neuroscience. This is due to this species' position as a premier model of human disease, having an array of genetic tools for manipulation of the animal in vivo, and being a small quadruped, a category for which few models exist. Here, the methods used to develop the first three-dimensional (3D) model of a mouse hindlimb and pelvis are described. The model, which represents bones, joints and 39 musculotendon units, was created through a combination of previously gathered muscle architecture data from microdissections, contrast-enhanced micro-computed tomography (CT) scanning and digital segmentation. The model allowed muscle moment arms as well as muscle forces to be estimated for each musculotendon unit throughout a range of joint rotations. Moment arm analysis supported the reliability of musculotendon unit placement within the model, and comparison to a previously published rat hindlimb model further supported the model's reliability. A sensitivity analysis performed on both the force-generating parameters and muscle's attachment points of the model indicated that the maximal isometric muscle moment is generally most sensitive to changes in either tendon slack length or the coordinates of insertion, although the degree to which the moment is affected depends on several factors. This model represents the first step in the creation of a fully dynamic 3D computer model of the mouse hindlimb and pelvis that has application to neuromuscular disease, comparative biomechanics and the neuromechanical basis of movement. Capturing the morphology and dynamics of the limb, it enables future dissection of the complex interactions between the nervous and musculoskeletal systems as well as the environment.

Static versus dynamic muscle modelling in extinct species: a biomechanical case study of the Australopithecus afarensis pelvis and lower extremity.

The force a muscle generates is dependent on muscle structure, in which fibre length, pennation angle and tendon slack length all influence force production. Muscles are not preserved in the fossil record and these parameters must be estimated when constructing a musculoskeletal model. Here, we test the capability of digitally reconstructed muscles of the Australopithecus afarensis model (specimen AL 288-1) to maintain an upright, single-support limb posture. Our aim was to ascertain the influence that different architectural estimation methods have on muscle specialisation and on the subsequent inferences that can be extrapolated about limb function. Parameters were estimated for 36 muscles in the pelvis and lower limb and seven different musculoskeletal models of AL 288-1 were produced. These parameters represented either a 'static' Hill-type muscle model (n = 4 variants) which only incorporated force, or instead a 'dynamic' Hill-type muscle model with an elastic tendon and fibres that could vary force-length-velocity properties (n = 3 variants). Each muscle's fibre length, pennation angle, tendon slack length and maximal isometric force were calculated based upon different input variables. Static (inverse) simulations were computed in which the vertical and mediolateral ground reaction forces (GRF) were incrementally increased until limb collapse (simulation failure). All AL 288-1 variants produced somewhat similar simulated muscle activation patterns, but the maximum vertical GRF that could be exerted on a single limb was not consistent between models. Three of the four static-muscle models were unable to support >1.8 times body weight and produced models that under-performed. The dynamic-muscle models were stronger. Comparative results with a human model imply that similar muscle group activations between species are needed to sustain single-limb support at maximally applied GRFs in terms of the simplified static simulations (e.g., same walking pose) used here. This approach demonstrated the range of outputs that can be generated for a model of an extinct individual. Despite mostly comparable outputs, the models diverged mostly in terms of strength.

The impacts of muscle-specific force-velocity properties on predictions of mouse muscle function during locomotion.

Introduction: The accuracy of musculoskeletal models and simulations as methods for predicting muscle functional outputs is always improving. However, even the most complex models contain various assumptions and simplifications in how muscle force generation is simulated. One common example is the application of a generalised ("generic") force-velocity relationship, derived from a limited data set to each muscle within a model, uniformly across all muscles irrespective of whether those muscles have "fast" or "slow" contractile properties. Methods: Using a previously built and validated musculoskeletal model and simulation of trotting in the mouse hindlimb, this work examines the predicted functional impact of applying muscle-specific force-velocity properties to typically fast (extensor digitorum longus; EDL) and slow-contracting (soleus; SOL) muscles. Results: Using "real" data led to EDL producing more positive work and acting significantly more spring-like, and soleus producing more negative work and acting more brake-like in function compared to muscles modelled using "generic" force-velocity data. Extrapolating these force-velocity properties to other muscles considered "fast" or "slow" also substantially impacted their predicted function. Importantly, this also further impacted EDL and SOL function beyond that seen when changing only their properties alone, to a point where they show an improved match to ex vivo experimental data. Discussion: These data suggest that further improvements to how musculoskeletal models and simulations predict muscle function should include the use of different values defining their force-velocity relationship depending on their fibre-type composition.

The Functional and Anatomical Impacts of Healthy Muscle Ageing.

Even "healthy" muscle ageing is often associated with substantial changes in muscle form and function and can lead to increased injury risks and significant negative impacts on quality of life. However, the impacts of healthy muscle ageing on the fibre architecture and microstructure of different muscles and muscle groups throughout the lower limb, and how these are related to their functional capabilities, are not fully understood. Here, a previously established framework of magnetic resonance and diffusion tensor imaging was used to measure the muscle volumes, intramuscular fat, fibre lengths and physiological cross-sectional areas of 12 lower limb muscles in a cohort of healthily aged individuals, which were compared to the same data from a young population. Maximum muscle forces were also measured from an isokinetic dynamometer. The more substantial interpopulation differences in architecture and functional performance were located within the knee extensor muscles, while the aged muscles were also more heterogeneous in muscle fibre type and atrophy. The relationships between architecture and muscle strength were also more significant in the knee extensors compared to other functional groups. These data highlight the importance of the knee extensors as a potential focus for interventions to negate the impacts of muscle ageing.

The association between muscle architecture and muscle spindle abundance.

Across the human body, skeletal muscles have a broad range of biomechanical roles that employ complex proprioceptive control strategies to successfully execute a desired movement. This information is derived from peripherally located sensory apparatus, the muscle spindle and Golgi tendon organs. The abundance of these sensory organs, particularly muscle spindles, is known to differ considerably across individual muscles. Here we present a comprehensive data set of 119 muscles across the human body including architectural properties (muscle fibre length, mass, pennation angle and physiological cross-sectional area) and statistically test their relationships with absolute spindle number and relative spindle abundance (the residual value of the linear regression of the log-transformed spindle number and muscle mass). These data highlight a significant positive relationship between muscle spindle number and fibre length, emphasising the importance of fibre length as an input into the central nervous system. However, there appears to be no relationship between muscles architecturally optimised to function as displacement specialists and their provision of muscle spindles. Additionally, while there appears to be regional differences in muscle spindle abundance, independent of muscle mass and fibre length, our data provide no support for the hypothesis that muscle spindle abundance is related to anatomical specialisation.

In vivo human lower limb muscle architecture dataset obtained using diffusion tensor imaging.

'Gold standard' reference sets of human muscle architecture are based on elderly cadaveric specimens, which are unlikely to be representative of a large proportion of the human population. This is important for musculoskeletal modeling, where the muscle force-generating properties of generic models are defined by these data but may not be valid when applied to models of young, healthy individuals. Obtaining individualized muscle architecture data in vivo is difficult, however diffusion tensor magnetic resonance imaging (DTI) has recently emerged as a valid method of achieving this. DTI was used here to provide an architecture data set of 20 lower limb muscles from 10 healthy adults, including muscle fiber lengths, which are important inputs for Hill-type muscle models commonly used in musculoskeletal modeling. Maximum isometric force and muscle fiber lengths were found not to scale with subject anthropometry, suggesting that these factors may be difficult to predict using scaling or optimization algorithms. These data also highlight the high level of anatomical variation that exists between individuals in terms of lower limb muscle architecture, which supports the need of incorporating subject-specific force-generating properties into musculoskeletal models to optimize their accuracy for clinical evaluation.

A Dynamic Simulation of Musculoskeletal Function in the Mouse Hindlimb During Trotting Locomotion.

Mice are often used as animal models of various human neuromuscular diseases, and analysis of these models often requires detailed gait analysis. However, little is known of the dynamics of the mouse musculoskeletal system during locomotion. In this study, we used computer optimization procedures to create a simulation of trotting in a mouse, using a previously developed mouse hindlimb musculoskeletal model in conjunction with new experimental data, allowing muscle forces, activation patterns, and levels of mechanical work to be estimated. Analyzing musculotendon unit (MTU) mechanical work throughout the stride allowed a deeper understanding of their respective functions, with the rectus femoris MTU dominating the generation of positive and negative mechanical work during the swing and stance phases. This analysis also tested previous functional inferences of the mouse hindlimb made from anatomical data alone, such as the existence of a proximo-distal gradient of muscle function, thought to reflect adaptations for energy-efficient locomotion. The results do not strongly support the presence of this gradient within the mouse musculoskeletal system, particularly given relatively high negative net work output from the ankle plantarflexor MTUs, although more detailed simulations could test this further. This modeling analysis lays a foundation for future studies of the control of vertebrate movement through the development of neuromechanical simulations.

Are mice good models for human neuromuscular disease? Comparing muscle excursions in walking between mice and humans.

BACKGROUND: The mouse is one of the most widely used animal models to study neuromuscular diseases and test new therapeutic strategies. However, findings from successful pre-clinical studies using mouse models frequently fail to translate to humans due to various factors. Differences in muscle function between the two species could be crucial but often have been overlooked. The purpose of this study was to evaluate and compare muscle excursions in walking between mice and humans. METHODS: Recently published musculoskeletal models of the mouse hindlimb and human lower limb were used to simulate muscle-tendon dynamics during mouse and human walking, a key daily activity. Muscle fiber length changes (fiber excursions) of 25 muscle homologs in the two species were calculated from these simulations and then compared. To understand potential causes of differences in fiber excursions in walking, joint excursions and muscle moment arms were also compared across one gait cycle. RESULTS: Most muscles (19 out of 25 muscles) of the mouse hindlimb had much smaller fiber excursions as compared to human lower limb muscles during walking. For these muscles, fiber excursions in mice were only 48 ± 19% of those in humans. The differences in fiber excursion between the two species were primarily due to the reduced joint excursions and smaller muscle moment arms in mice as compared to humans. CONCLUSIONS: Since progressive neuromuscular diseases, such as Duchenne muscular dystrophy, are known to be accelerated by damage accumulated from active muscle lengthening, these results suggest that biomechanical differences in muscle function during walking between mice and humans may impede the translations of knowledge gained from mouse models to humans. This knowledge would add a fresh perspective on how pre-clinical studies on mice might be better designed to improve translation to human clinical trials.

Musculoskeletal Geometry, Muscle Architecture and Functional Specialisations of the Mouse Hindlimb.

Mice are one of the most commonly used laboratory animals, with an extensive array of disease models in existence, including for many neuromuscular diseases. The hindlimb is of particular interest due to several close muscle analogues/homologues to humans and other species. A detailed anatomical study describing the adult morphology is lacking, however. This study describes in detail the musculoskeletal geometry and skeletal muscle architecture of the mouse hindlimb and pelvis, determining the extent to which the muscles are adapted for their function, as inferred from their architecture. Using I2KI enhanced microCT scanning and digital segmentation, it was possible to identify 39 distinct muscles of the hindlimb and pelvis belonging to nine functional groups. The architecture of each of these muscles was determined through microdissections, revealing strong architectural specialisations between the functional groups. The hip extensors and hip adductors showed significantly stronger adaptations towards high contraction velocities and joint control relative to the distal functional groups, which exhibited larger physiological cross sectional areas and longer tendons, adaptations for high force output and elastic energy savings. These results suggest that a proximo-distal gradient in muscle architecture exists in the mouse hindlimb. Such a gradient has been purported to function in aiding locomotor stability and efficiency. The data presented here will be especially valuable to any research with a focus on the architecture or gross anatomy of the mouse hindlimb and pelvis musculature, but also of use to anyone interested in the functional significance of muscle design in relation to quadrupedal locomotion.