Anthropometric Scaling of Anatomical Datasets for Subject-Specific Musculoskeletal Modelling of the Shoulder

Quantifying changes in shoulder orientation between the prone and supine positions from magnetic resonance imaging

BACKGROUND

Predicting breast tissue motion using biomechanical models can provide navigational guidance during breast cancer treatment procedures. These models typically do not account for changes in posture between procedures. Difference in shoulder position can alter the shape of the pectoral muscles and breast. A greater understanding of the differences in the shoulder orientation between prone and supine could improve the accuracy of breast biomechanical models.

METHODS

19 landmarks were placed on the sternum, clavicle, scapula, and humerus of the shoulder girdle in prone and supine breast MRIs (N = 10). These landmarks were used in an optimization framework to fit subject-specific skeletal models and compare joint angles of the shoulder girdle between these positions.

FINDINGS

The mean Euclidean distance between joint locations from the fitted skeletal model and the manually identified joint locations was 15.7 mm ± 2.7 mm. Significant differences were observed between prone and supine. Compared to supine position, the shoulder girdle in the prone position had the lateral end of the clavicle in more anterior translation (i.e., scapula more protracted) (P < 0.05), the scapula in more protraction (P < 0.01), the scapula in more upward rotation (associated with humerus elevation) (P < 0.05); and the humerus more elevated (P < 0.05) for both the left and right sides.

INTERPRETATION

Shoulder girdle orientation was found to be different between prone and supine. These differences would affect the shape of multiple pectoral muscles, which would affect breast shape and the accuracy of biomechanical models.

The effect of muscle atrophy in people with unilateral transtibial amputation for three activities: Gait alone does not tell the whole story

Amputation imposes significant challenges in locomotion to millions of people with limb loss worldwide. The decline in the use of the residual limb results in muscle atrophy that affects musculoskeletal dynamics in daily activities. The aim of this study was to quantify the lower limb muscle volume discrepancy based on magnetic resonance (MR) imaging and to combine this with motion analysis and musculoskeletal modelling to quantify the effects in the dynamics of key activities of daily living. Eight male participants with traumatic unilateral transtibial amputation were recruited who were at least six months after receiving their definitive prostheses. The muscle volume discrepancies were found to be largest at the knee extensors (35 %, p = 0.008), followed by the hip abductors (17 %, p = 0.008). Daily activities (level walking, standing up from a chair and ascending one step) were measured in a motion analysis laboratory and muscle and joint forces quantified using a detailed musculoskeletal model for people with unilateral transtibial amputation which was calibrated in terms of the muscle volume discrepancies post-amputation at a subject-specific level. Knee extensor muscle forces were lower at the residual limb than the intact limb for all activities (p ≤ 0.008); residual limb muscle forces of the hip abductors (p ≤ 0.031) and adductors (p ≤ 0.031) were lower for standing-up and ascending one step. While the reduced knee extensor force has been reported by other studies, our results suggest a new biomechanically-based mitigation strategy to improve functional mobility, which could be achieved through strengthening of the hip abd/adductor muscles.

Biceps Tenodesis cannot be used as primary treatment option in baseball pitchers with intact rotator cuff muscles

Background Surgeons remain hesitant to perform biceps tenodesis in athletes with type II superior labrum anterior-to-posterior tears due to the lack of reported clinical outcomes for individual overhead throwing sports and associated concerns that this may predispose the joint to instability. This study aimed to assess the effect of biceps tenodesis on shoulder stability for major overhead throwing sports to aid sport-specific surgical decision-making for athletes with type II superior labrum anterior-to-posterior tears.

METHODS

This is a combined modelling and experimental study. Motion data and external forces were measured from 13 participants performing five overhead throwing motions. These data served as input into a musculoskeletal shoulder model that quantifies shoulder stability and muscle loading.

FINDINGS

The loading of the long head of the biceps brachii decreases significantly following biceps tenodesis in three overhead throwing motions (p = 0.02). The loss in joint stability following biceps tenodesis is compensated by a non-significant increase in rotator cuff muscle force which maintains shoulder stability across all overhead throwing motions, except baseball pitching (p = 0.01). The presence of a full-thickness supraspinatus tear post biceps tenodesis further decreases shoulder stability in four of the five overhead throwing motions (p = 0.01).

INTERPRETATION

The study findings demonstrate that an increase in rotator cuff muscle force maintains joint stability for all overhead throwing motions post biceps tenodesis, except baseball pitching. As the presence of a full-thickness tear of the supraspinatus significantly reduces joint stability, biceps tenodesis may be used as a primary treatment in overhead throwing athletes with intact rotator cuff muscles, except baseball pitchers.

LEVEL OF EVIDENCE

Controlled Laboratory Study; Level of Evidence 3.

The effects of anatomical errors on shoulder kinematics computed using multi-body models

Joint motion calculated using multi-body models and inverse kinematics presents many advantages over direct marker-based calculations. However, the sensitivity of the computed kinematics is known to be partly caused by the model and could also be influenced by the participants' anthropometry and sex. This study aimed to compare kinematics computed from an anatomical shoulder model based on medical images against a scaled-generic model and quantify the effects of anatomical errors and participants' anthropometry on the calculated joint angles. Twelve participants have had planar shoulder movements experimentally captured in a motion lab, and their shoulder anatomy imaged using an MRI scanner. A shoulder multi-body dynamics model was developed for each participant, using both an image-based approach and a scaled-generic approach. Inverse kinematics have been performed using the two different modelling procedures and the three different experimental motions. Results have been compared using Bland-Altman analysis of agreement and further analysed using multi-linear regressions. Kinematics computed via an anatomical and a scaled-generic shoulder models differed in average from 3.2 to 5.4 degrees depending on the task. The MRI-based model presented smaller limits of agreement to direct kinematics than the scaled-generic model. Finally, the regression model predictors, including anatomical errors, sex, and BMI of the participant, explained from 41 to 80% of the kinematic variability between model types with respect to the task. This study highlighted the consequences of modelling precision, quantified the effects of anatomical errors on the shoulder kinematics, and showed that participants' anthropometry and sex could indirectly affect kinematic outcomes.

Statistical Quantification of the Effects of Marker Misplacement and Soft-Tissue Artifact on Shoulder Kinematics and Kinetics

The assessment of shoulder kinematics and kinetics are commonly undertaken biomechanically and clinically by using rigid-body models and experimental skin-marker trajectories. However, the accuracy of these trajectories is plagued by inherent skin-based marker errors due to marker misplacements (offset) and soft-tissue artifacts (STA). This paper aimed to assess the individual contribution of each of these errors to kinematic and kinetic shoulder outcomes computed using a shoulder rigid-body model. Baseline experimental data of three shoulder planar motions in a young healthy adult were collected. The baseline marker trajectories were then perturbed by simulating typically observed population-based offset and/or STA using a probabilistic Monte-Carlo approach. The perturbed trajectories were then used together with a shoulder rigid-body model to compute shoulder angles and moments and study their accuracy and variability against baseline. Each type of error was studied individually, as well as in combination. On average, shoulder kinematics varied by 3%, 6% and 7% due to offset, STA or combined errors, respectively. Shoulder kinetics varied by 11%, 27% and 28% due to offset, STA or combined errors, respectively. In conclusion, to reduce shoulder kinematic and kinetic errors, one should prioritise reducing STA as they have the largest error contribution compared to marker misplacements.

Comparing optical and electromagnetic tracking systems to facilitate compatibility in sports kinematics data

Electromagnetic (EM) tracking has been used to quantify biomechanical parameters of the lower limb and lumbar spine during ergometer rowing to improve performance and reduce injury. Optical motion capture (OMC) is potentially better suited to measure comprehensive whole-body dynamics in rowing. This study compared accuracy and precision of EM and OMC displacements by simultaneously recording kinematics during rowing trials at low, middle, and high rates on an instrumented ergometer (n=12). Trajectories calculated from OMC and EM sensors attached to the pelvis, lumbar spine, and right leg were highly correlated, but EM tracking lagged behind ergometer and OMC tracking by approximately 6%, yielding large RMS errors. When this phase-lag was corrected by least squares minimization, agreement between systems improved. Both systems demonstrated an ability to adequately track large dynamic compound movements in the sagittal plane but struggled at times to precisely track small displacements and narrow angular ranges in medial/lateral and superior/inferior directions. An OMC based tracking methodology can obtain equivalence with a previously validated EM system, for spine and lower limb metrics. Improvements in speed and consistency of data acquisition with OMC are beneficial for dynamic motion studies. Compatibility ensures continuity by maintaining the ability to compare to prior work.

In Vitro Kinematic Analysis of Single Axis Radius Posterior-Substituting Total Knee Arthroplasty

This is an experimental study. As current posterior-substituting (PS) total knee arthroplasties have been reported to incompletely restore intrinsic joint biomechanics of the healthy knee, the recently designed single axis radius PS knee system was introduced to increase posterior femoral translation and promote ligament isometry. As there is a paucity of data available regarding its ability to replicate healthy knee biomechanics, this study aimed to assess joint and articular contact kinematics as well as ligament isometry of the contemporary single axis radius PS knee system. Implant kinematics were measured from 11 cadaveric knees using an in vitro robotic testing system. In addition, medial collateral ligament (MCL) and lateral collateral ligament (LCL) forces were quantified under simulated functional loads during knee flexion for the contemporary PS knee system. Posterior femoral translation between the intact knee and the single axis radius PS knee system differed significantly (p < 0.05) at 60, 90, and 120 degrees of flexion. The LCL force at 60 degrees (9.06 ± 2.81 N) was significantly lower (p < 0.05) than those at 30, 90, and 120 degrees of flexion, while MCL forces did not differ significantly throughout the range of tested flexion angles. The results from this study suggest that although the contemporary single axis radius PS knee system has the potential to mimic the intact knee kinematics under muscle loading during flexion extension due to its design features, single axis radius PS knee system did not fully replicate posterior femoral translation and ligament isometry of the healthy knee during knee flexion.

An Anatomical Atlas-Based Scaling Study for Quantifying Muscle and Hip Joint Contact Forces in Above and Through-Knee Amputees Using Validated Musculoskeletal Modelling

OBJECTIVE

Customisation of musculoskeletal modelling using magnetic resonance imaging (MRI) significantly improves the model accuracy, but the process is time consuming and computationally intensive. This study hypothesizes that linear scaling to a lower limb amputee model with anthropometric similarity can accurately predict muscle and joint contact forces.

METHODS

An MRI-based anatomical atlas, comprising 18 trans-femoral and through-knee traumatic lower limb amputee models, is developed. Gait data, using a 10-camera motion capture system with two force plates, and surface electromyography (EMG) data were collected. Muscle and hip joint contact forces were quantified using musculoskeletal modelling. The predicted muscle activations from the subject-specific models were validated using EMG recordings. Anthropometry based multiple linear regression models, which minimize errors in force predictions, are presented.

RESULTS

All predictions showed excellent (error interval c = 0-0.15), very good (c = 0.15-0.30) or good (c = 0.30-0.45) similarity to the EMG data, demonstrating accurate computation of muscle activations. The primary predictors of discrepancies in force predictions were differences in pelvis width (p < 0.001), body mass index (BMI, p < 0.001) and stump length to pelvis width ratio (p < 0.001) between the respective individual and underlying dataset.

CONCLUSION

Linear scaling to a model with the most similar pelvis width, BMI and stump length to pelvis width ratio results in modelling outcomes with minimal errors.

SIGNIFICANCE

This study provides robust tools to perform accurate analyses of musculoskeletal mechanics for high-functioning lower limb military amputees, thus facilitating the further understanding and improvement of the amputee's function. The atlas is available in an open source repository.

Review of musculoskeletal modelling in a clinical setting: Current use in rehabilitation design, surgical decision making and healthcare interventions

BACKGROUND

Musculoskeletal modelling is a common means by which to non-invasively analyse movement. Such models have largely been used to observe function in both healthy and patient populations. However, utility in a clinical environment is largely unknown. The aim of this review was to explore existing uses of musculoskeletal models as a clinical intervention, or decision-making, tool.

METHODS

A literature search was performed using PubMed and Scopus to find articles published since 2010 and relating to musculoskeletal modelling and joint and muscle forces.

FINDINGS

4662 abstracts were found, of which 39 relevant articles were reviewed. Journal articles were categorised into 5 distinct groups: non-surgical treatment, orthoses assessment, surgical decision making, surgical intervention assessment and rehabilitation regime assessment. All reviewed articles were authored by collaborations between clinicians and engineers/modellers. Current uses included insight into the development of osteoarthritis, identifying candidates for hamstring lengthening surgery, and the assessment of exercise programmes to reduce joint damage.

INTERPRETATION

There is little evidence showing the use of musculoskeletal modelling as a tool for patient care, despite the ability to assess long-term joint loading and muscle overuse during functional activities, as well as clinical decision making to avoid unfavourable treatment outcomes. Continued collaboration between model developers should aim to create clinically-friendly models which can be used with minimal input and experience by healthcare professionals to determine surgical necessity and suitability for rehabilitation regimes, and in the assessment of orthotic devices.

An open-source musculoskeletal model of the lumbar spine and lower limbs: a validation for movements of the lumbar spine

Musculoskeletal models of the lumbar spine have been developed with varying levels of detail for a wide range of clinical applications. Providing consistency is ensured throughout the modelling approach, these models can be combined with other computational models and be used in predictive modelling studies to investigate bone health deterioration and the associated fracture risk. To provide precise physiological loading conditions for such predictive modelling studies, a new full-body musculoskeletal model including a detailed and consistent representation of the lower limbs and the lumbar spine was developed. The model was assessed against in vivo measurements from the literature for a range of spine movements representative of daily living activities. Comparison between model estimations and electromyography recordings was also made for a range of lifting tasks. This new musculoskeletal model will provide a comprehensive physiological mechanical environment for future predictive finite element modelling studies on bone structural adaptation. It is freely available on https://simtk.org/projects/llsm/.

A biomechanical confirmation of the relationship between critical shoulder angle (CSA) and articular joint loading

BACKGROUND

The critical shoulder angle (CSA) has been shown to be correlated with shoulder disease states. The biomechanical hypothesis to explain this correlation is that the CSA changes the shear and compressive forces on the shoulder. The objective of this study is to test this hypothesis by use of a validated computational shoulder model. Specifically, this study assesses the impact on glenohumeral biomechanics of modifying the CSA.

METHODS

An inverse dynamics 3-dimensional musculoskeletal model of the shoulder was used to quantify muscle forces and glenohumeral joint forces. The CSA was changed by altering the attachment point of the middle deltoid into a normal CSA (33°), a reduced CSA of 28°, and an increased CSA of 38°. Subject-specific kinematics of slow and fast speed abduction in the scapular plane and slow and fast forward flexion measured by a 3-dimensional motion capture system were used to quantify joint reaction shear and compressive forces.

RESULTS

Increasing the CSA results in increased superior-inferior forces (shearing forces; integrated over the range of motion; P < .05). Reducing CSA results in increased lateromedial (compressive) forces for both the maximum and integrated sum of the forces over the whole motion (P < .01).

DISCUSSION/CONCLUSION

Changes in the CSA modify glenohumeral joint biomechanics with increasing CSA producing higher shear forces that could contribute to rotator cuff overuse, whereas reducing the CSA results in higher compressive forces that contribute to joint wear.

Accuracy and kinematics consistency of marker-based scaling approaches on a lower limb model: a comparative study with imagery data

Medical images are not typically included in protocol of motion laboratories. Thus, accurate scaling of musculoskeletal models from optoelectronic data are important for any biomechanical analysis. The aim of the current study was to identify a scaling method based on optoelectronic data, inspired from literature, which could offer the best trade-off between accurate geometrical parameters (segment lengths, orientation of joint axes, marker coordinates) and consistent inverse kinematics outputs (kinematic error, joint angles). The methods were applied on 26 subjects and assessed with medical imagery building EOS-based models, considered as a reference. The main contribution of this paper is to show that the marker-based scaling followed by an optimisation of orientation joint axes and markers local coordinates, gives the most consistent scaling and joint angles with EOS-based models. Thus, when a non-invasive mean with an optoelectronic system is considered, a marker-based scaling is preliminary needed to get accurate segment lengths and to optimise joint axes and marker local coordinates to reduce kinematic errors.AbbrevationsAJCAnkle joint centreCKEcumulative kinematic errorDoFdegree of freedomEBEOS-basedHBheight-basedHJChip joint centreKJCknee joint centreMBmarker-basedMSMmusculoskeletal modelsSPMstatistical parametric mappingSTAsoft tissue artifactEB_a.m∗EOS-based with optimised joint axes, and all model markers coordinatesMB_a.m∗marker-based with optimised joint axes, and all model markers coordinatesMB_l.a.mmarker-based with optimised segment lengths, joint axes, and selected model markers coordinatesASISanterior superior illiac spinePSISposterior superior illiac spine.