The relationship between tibiofemoral geometry and musculoskeletal function during normal activity

Body shape from birth to adulthood is associated with skeletal development: A Mendelian randomization study

BACKGROUND

Observational studies have shown that childhood obesity is associated with adult bone health but yield inconsistent results. We aimed to explore the potential causal association between body shape and skeletal development.

METHODS

We used two-sample Mendelian randomization (MR) to estimate causal relationships between body shape from birth to adulthood and skeletal phenotypes, with exposures including placental weight, birth weight, childhood obesity, BMI, lean mass, fat mass, waist circumference, and hip circumference. Independent genetic instruments associated with the exposures at the genome-wide significance level (P < 5 × 10-8) were selected from corresponding large-scale genome-wide association studies. The inverse-variance weighted analysis was chosen as the primary method, and complementary MR analyses included the weighted median, MR-Egger, weighted mode, and simple mode.

RESULTS

The MR analysis shows strong evidence that childhood (β = -1.29 × 10-3, P = 8.61 × 10-5) and adulthood BMI (β = -1.28 × 10-3, P = 1.45 × 10-10) were associated with humerus length. Tibiofemoral angle was negatively associated with childhood BMI (β = -3.60 × 10-1, P = 3.00 × 10-5) and adolescent BMI (β = -3.62 × 10-1, P = 2.68 × 10-3). In addition, genetically predicted levels of appendicular lean mass (β = 1.16 × 10-3, P = 1.49 × 10-13), whole body fat mass (β = 1.66 × 10-3, P = 1.35 × 10-9), waist circumference (β = 1.72 × 10-3, P = 6.93 × 10-8) and hip circumference (β =1.28 × 10-3, P = 4.34 × 10-6) were all associated with tibia length. However, we found no causal association between placental weight, birth weight and bone length/width.

CONCLUSIONS

This large-scale MR analysis explores changes in growth patterns in the length/width of major bone sites, highlighting the important role of childhood body shape in bone development and providing insights into factors that may drive bone maturation.

Scan-Free and Fully Automatic Tracking of Native Knee Anatomy from Dynamic Stereo-Radiography with Statistical Shape and Intensity Models

Kinematic tracking of native anatomy from stereo-radiography provides a quantitative basis for evaluating human movement. Conventional tracking procedures require significant manual effort and call for acquisition and annotation of subject-specific volumetric medical images. The current work introduces a framework for fully automatic tracking of native knee anatomy from dynamic stereo-radiography which forgoes reliance on volumetric scans. The method consists of three computational steps. First, captured radiographs are annotated with segmentation maps and anatomic landmarks using a convolutional neural network. Next, a non-convex polynomial optimization problem formulated from annotated landmarks is solved to acquire preliminary anatomy and pose estimates. Finally, a global optimization routine is performed for concurrent refinement of anatomy and pose. An objective function is maximized which quantifies similarities between masked radiographs and digitally reconstructed radiographs produced from statistical shape and intensity models. The proposed framework was evaluated against manually tracked trials comprising dynamic activities, and additional frames capturing a static knee phantom. Experiments revealed anatomic surface errors routinely below 1.0 mm in both evaluation cohorts. Median absolute errors of individual bone pose estimates were below 1.0∘ or mm for 15 out of 18 degrees of freedom in both evaluation cohorts. Results indicate that accurate pose estimation of native anatomy from stereo-radiography may be performed with significantly reduced manual effort, and without reliance on volumetric scans.

Fully automatic tracking of native knee kinematics from stereo-radiography with digitally reconstructed radiographs

Precise measurement of joint-level motion from stereo-radiography facilitates understanding of human movement. Conventional procedures for kinematic tracking require significant manual effort and are time intensive. The current work introduces a method for fully automatic tracking of native knee kinematics from stereo-radiography sequences. The framework consists of three computational steps. First, biplanar radiograph frames are annotated with segmentation maps and key points using a convolutional neural network. Next, initial bone pose estimates are acquired by solving a polynomial optimization problem constructed from annotated key points and anatomic landmarks from digitized models. A semidefinite relaxation is formulated to realize the global minimum of the non-convex problem. Pose estimates are then refined by registering computed tomography-based digitally reconstructed radiographs to masked radiographs. A novel rendering method is also introduced which enables generating digitally reconstructed radiographs from computed tomography scans with inconsistent slice widths. The automatic tracking framework was evaluated with stereo-radiography trials manually tracked with model-image registration, and with frames which capture a synthetic leg phantom. The tracking method produced pose estimates which were consistently similar to manually tracked values; and demonstrated pose errors below 1.0 degree or millimeter for all femur and tibia degrees of freedom in phantom trials. Results indicate the described framework may benefit orthopaedics and biomechanics applications through acceleration of kinematic tracking.

A frame orientation optimisation method for consistent interpretation of kinematic signals

In clinical movement biomechanics, kinematic data are often depicted as waveforms (i.e. signals), characterising the motion of articulating joints. Clinically meaningful interpretations of the underlying joint kinematics, however, require an objective understanding of whether two different kinematic signals actually represent two different underlying physical movement patterns of the joint or not. Previously, the accuracy of IMU-based knee joint angles was assessed using a six-degrees-of-freedom joint simulator guided by fluoroscopy-based signals. Despite implementation of sensor-to-segment corrections, observed errors were clearly indicative of cross-talk, and thus inconsistent reference frame orientations. Here, we address these limitations by exploring how minimisation of dedicated cost functions can harmonise differences in frame orientations, ultimately facilitating consistent interpretation of articulating joint kinematic signals. In this study, we present and investigate a frame orientation optimisation method (FOOM) that aligns reference frames and corrects for cross-talk errors, hence yielding a consistent interpretation of the underlying movement patterns. By executing optimised rotational sequences, thus producing angular corrections around each axis, we enable a reproducible frame definition and hence an approach for reliable comparison of kinematic data. Using this approach, root-mean-square errors between the previously collected (1) IMU-based data using functional joint axes, and (2) simulated fluoroscopy-based data relying on geometrical axes were almost entirely eliminated from an initial range of 0.7°-5.1° to a mere 0.1°-0.8°. Our results confirm that different local segment frames can yield different kinematic patterns, despite following the same rotation convention, and that appropriate alignment of reference frame orientation can successfully enable consistent kinematic interpretation.

Can a Musculoskeletal Model Adapted to Knee Implant Geometry Improve Prediction of 3D Contact Forces and Moments?

Tibiofemoral contact loads are crucial parameters in the onset and progression of osteoarthrosis. While contact loads are frequently estimated from musculoskeletal models, their customization is often limited to scaling musculoskeletal geometry or adapting muscle lines. Moreover, studies have usually focused on superior-inferior contact force without investigating three-dimensional contact loads. Using experimental data from six patients with instrumented total knee arthroplasty (TKA), this study customized a lower limb musculoskeletal model to consider the positioning and the geometry of the implant at knee level. Static optimization was performed to estimate tibiofemoral contact forces and contact moments as well as musculotendinous forces. Predictions from both a generic and a customized model were compared to the instrumented implant measurements. Both models accurately predict superior-inferior (SI) force and abduction-adduction (AA) moment. Notably, the customization improves prediction of medial-lateral (ML) force and flexion-extension (FE) moments. However, there is subject-dependent variability in the prediction of anterior-posterior (AP) force. The customized models presented here predict loads on all joint axes and in most cases improve prediction. Unexpectedly, this improvement was more limited for patients with more rotated implants, suggesting a need for further model adaptations such as muscle wrapping or redefinition of hip and ankle joint centers and axes.

Kinetics influence of multibody kinematics optimisation for soft tissue artefact compensation

Soft tissue artefact (STA) remains a major source of error in human movement analysis. The multibody kinematics optimisation (MKO) approach is widely stated as a solution to reduce the effects of STA. This study aimed at assessing the influence of the MKO STA-compensation on the errors of estimation of the knee intersegment moments. Experimental data were issued from the CAMS-Knee dataset where six participants with instrumented total knee arthroplasty performed five activities of daily living: gait, downhill walking, stair descent, squat, and sit-to-stand. Kinematics was measured both on the basis of skin markers and a mobile mono-plane fluoroscope, used to obtain the STA-free bone movement. For four different lower limb models and one corresponding to a single-body kinematics optimization (SKO), knee intersegmental moments (estimated using model-derived kinematics and ground reaction force) were compared with an estimate based on the fluoroscope. Considering all participants and activities, mean root mean square differences were the largest along the adduction/abduction axis: of 3.22Nm with a SKO approach, 3.49Nm with the three-DoF knee model, and 7.66Nm, 8.52Nm, and 8.54Nm with the one-DoF knee models. Results showed that adding joint kinematics constraints can increase the estimation errors of the intersegmental moment. These errors came directly from the errors in the estimation of the position of the knee joint centre induced by the constraints. When using a MKO approach, we recommend to analyse carefully joint centre position estimates that do not remain close to the one obtained with a SKO approach.

Determining the relationship between tibiofemoral geometry and passive motion with partial least squares regression

Tibiofemoral geometry influences knee passive motion and understanding their relationship can provide insight into knee function and mechanisms of injury. However, the complexity of the geometric constraints has made characterizing the relationship challenging. The aim of this study was to determine the tibiofemoral bone geometries that explain the variation in passive motion using a partial least squares regression (PLSR) model. The PLSR model was developed for 29 healthy cadaver specimens (10 female, 19 male) with femur and tibia geometries retrieved from MRI images and six degree-of-freedom tibiofemoral kinematics determined during a flexion cycle with minimal medial pressure. The first 13 partial least squares (PLS) components explained 90% of the variation in the kinematics and accounted for 89% of the variation in geometry. The first three PLS components which shared geometric changes to particular surface congruencies of the tibial and femoral condyles explained the most amount of variation in the kinematics, primarily in anterior-posterior translation. Meanwhile, variations in femoral condyle width and the intercondylar space, tibia plateau size and conformity, and tibia eminences heights in PLS 2 and 4 explained the greatest amount of variation in internal-external rotation. PLS 4 exhibiting variation in overall size of the knee accounted for greatest amount of variation in geometry (50%) and had the greatest influence on the abduction-adduction motion and some on internal-external rotation but, overall, explained only a small proportion of the kinematics (10%). Elucidating the complex relationship between tibiofemoral bone geometry and passive kinematics may help personalize treatments for improved functional outcomes in patients.

Knee Joint Contact Forces during High-Risk Dynamic Tasks: 90° Change of Direction and Deceleration Movements

Pivoting sports expose athletes to a high risk of knee injuries, mainly due to mechanical overloading of the joint which shatters overall tissue integrity. The present study explored the magnitude of tibiofemoral contact forces (TFCF) in high-risk dynamic tasks. A novel musculoskeletal model with modifiable frontal plane knee alignment was developed to estimate the total, medial, and lateral TFCF developed during vigorous activities. Thirty-one competitive soccer players performing deceleration and 90° sidestepping tasks were assessed via 3D motion analysis by using a marker-based optoelectronic system and TFCF were assessed via OpenSim software. Statistical parametric mapping was used to investigate the effect of frontal plane alignment, compartment laterality, and varus-valgus genu on TFCF. Further, in consideration of specific risk factors, sex influence was also assessed. A strong correlation (R = 0.71 ÷ 0.98, p < 0.001) was found between modification of compartmental forces and changes in frontal plane alignment. Medial and lateral TFCF were similar throughout most of the tasks with the exception of the initial phase, where the lateral compartment had to withstand to higher loadings (1.5 ÷ 3 BW higher, p = 0.010). Significant sex differences emerged in the late phase of the deceleration task. A comprehensive view of factors influencing the mediolateral distribution of TFCF would benefit knee injury prevention and rehabilitation in sport activities.

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.

Altered Muscle Contributions are Required to Support the Stance Limb During Voluntary Toe-Walking

Toe-walking characterizes several neuromuscular conditions and is associated with a reduction in gait stability and efficiency, as well as in life quality. The optimal choice of treatment depends on a correct understanding of the underlying pathology and on the individual biomechanics of walking. The objective of this study was to describe gait deviations occurring in a cohort of healthy adult subjects when mimicking a unilateral toe-walking pattern compared to their normal heel-to-toe gait pattern. The focus was to characterize the functional adaptations of the major lower-limb muscles which are required in order to toe walk. Musculoskeletal modeling was used to estimate the required muscle contributions to the joint sagittal moments. The support moment, defined as the sum of the sagittal extensive moments at the ankle, knee, and hip joints, was used to evaluate the overall muscular effort necessary to maintain stance limb stability and prevent the collapse of the knee. Compared to a normal heel-to-toe gait pattern, toe-walking was characterized by significantly different lower-limb kinematics and kinetics. The altered kinetic demands at each joint translated into different necessary moment contributions from most muscles. In particular, an earlier and prolonged ankle plantarflexion contribution was required from the soleus and gastrocnemius during most of the stance phase. The hip extensors had to provide a higher extensive moment during loading response, while a significantly higher knee extension contribution from the vasti was necessary during mid-stance. Compensatory muscular activations are therefore functionally required at every joint level in order to toe walk. A higher support moment during toe-walking indicates an overall higher muscular effort necessary to maintain stance limb stability and prevent the collapse of the knee. Higher muscular demands during gait may lead to fatigue, pain, and reduced quality of life. Toe-walking is indeed associated with significantly larger muscle forces exerted by the quadriceps to the patella and prolonged force transmission through the Achilles tendon during stance phase. Optimal treatment options should therefore account for muscular demands and potential overloads associated with specific compensatory mechanisms.

Modeling Musculoskeletal Dynamics during Gait: Evaluating the Best Personalization Strategy through Model Anatomical Consistency

No consensus exists on how to model human articulations within MSK models for the analysis of gait dynamics. We propose a method to evaluate joint models and we apply it to three models with different levels of personalization. The method evaluates the joint model’s adherence to the MSK hypothesis of negligible joint work by quantifying ligament and cartilage deformations resulting from joint motion; to be anatomically consistent, these deformations should be minimum. The contrary would require considerable external work to move the joint, violating a strong working hypothesis and raising concerns about the credibility of the MSK outputs. Gait analysis and medical resonance imaging (MRI) from ten participants were combined to build lower limb subject-specific MSK models. MRI-reconstructed anatomy enabled three levels of personalization using different ankle joint models, in which motion corresponded to different ligament elongation and cartilage co-penetration. To estimate the impact of anatomical inconsistency in MSK outputs, joint internal forces resulting from tissue deformations were computed for each joint model and MSK simulations were performed ignoring or considering their contribution. The three models differed considerably for maximum ligament elongation and cartilage co-penetration (between 5.94 and 50.69% and between −0.53 and −5.36 mm, respectively). However, the model dynamic output from the gait simulations were similar. When accounting for the internal forces associated with tissue deformation, outputs changed considerably, the higher the personalization level the smaller the changes. Anatomical consistency provides a solid method to compare different joint models. Results suggest that consistency grows with personalization, which should be tailored according to the research question. A high level of anatomical consistency is recommended when individual specificity and the behavior of articular structures is under investigation.

Statistical-Shape Prediction of Lower Limb Kinematics During Cycling, Squatting, Lunging, and Stepping-Are Bone Geometry Predictors Helpful?

Purpose: Statistical shape methods have proven to be useful tools in providing statistical predications of several clinical and biomechanical features as to analyze and describe the possible link with them. In the present study, we aimed to explore and quantify the relationship between biometric features derived from imaging data and model-derived kinematics. Methods: Fifty-seven healthy males were gathered under strict exclusion criteria to ensure a sample representative of normal physiological conditions. MRI-based bone geometry was established and subject-specific musculoskeletal simulations in the Anybody Modeling System enabled us to derive personalized kinematics. Kinematic and shape findings were parameterized using principal component analysis. Partial least squares regression and canonical correlation analysis were then performed with the goal of predicting motion and exploring the possible association, respectively, with the given bone geometry. The relationship of hip flexion, abduction, and rotation, knee flexion, and ankle flexion with a subset of biometric features (age, length, and weight) was also investigated. Results: In the statistical kinematic models, mean accuracy errors ranged from 1.60° (race cycling) up to 3.10° (lunge). When imposing averaged kinematic waveforms, the reconstruction errors varied between 4.59° (step up) and 6.61° (lunge). A weak, yet clinical irrelevant, correlation between the modes describing bone geometry and kinematics was observed. Partial least square regression led to a minimal error reduction up to 0.42° compared to imposing gender-specific reference curves. The relationship between motion and the subject characteristics was even less pronounced with an error reduction up to 0.21°. Conclusion: The contribution of bone shape to model-derived joint kinematics appears to be relatively small and lack in clinical relevance.

Mechanical, biochemical, and morphological topography of ovine knee cartilage

The knee is the most common site for translational cartilage research in sheep, though topographic features of articular cartilage across surfaces are unspecified. We aimed to characterize the mechanical, morphological, and biochemical properties of articular cartilage across ovine knee surfaces and document variations between and within surface locations. Regions of interest (ROIs) were delineated across surfaces of 10 healthy ovine knees. Articular cartilage at each ROI was measured for creep indentation, thickness, and glycosaminoglycan (GAG) and collagen content. Variables were compared between surface locations (trochlea, and lateral [LFC] and medial [MFC] femoral condyles) and between ROIs within each surface location. Correlations between variables were also assessed. Articular surface location had a significant effect on creep (P < .0001), thickness (P < .0001), and collagen (P = .0007), but not GAG (P = .28). Significant differences in percent creep between ROIs were found within the LFC (P < .0001), MFC (P < .0001), and trochlea (P = .0002). Cartilage thickness was different between ROIs within the LFC, MFC, and trochlea (all P < .0001). The LFC (P = .002) and trochlea (P = .01) each had significant differences in GAG between ROIs. Collagen content between ROIs was different within the LFC (P = .0003), MFC (P = .0005), and trochlea (P < .0001). Collagen content was correlated with thickness (r = -.55), percent creep (r = .47), and GAG (r = -.21). Percent creep was correlated with thickness (r = -.64) and GAG (r = -.19). Topographic variations in mechanical, morphological, and biochemical properties exist across knee cartilage surfaces in sheep. Recognition of this variability is important to optimize study protocols and improve accuracy of results.

Quantification of the errors associated with marker occlusion in stereophotogrammetric systems and implications on gait analysis

Optoelectronic stereophotogrammetric systems (OSSs) represent the standard for gait analysis. Despite widespread, their reported accuracy in nominal working conditions shows a variability of several orders of magnitude, ranging from few microns to several millimetres. No clear explanation for this variability has been provided yet. We hypothesized that this reflects an error affecting OSS outcomes when some of the tracked markers are totally or partially occluded. The aim of this paper is to quantify this error in static and dynamic conditions, also distinguishing between total and partial marker occlusion. A Vicon system featuring 8 cameras is employed in this study. Two camera distributions, one designed to maximize OSS accuracy and another one representative of a typical gait setup, are investigated. For both the setups, static and dynamic tests are performed, evaluating the different impact of partial and total marker occlusions. Marker occlusions significantly affected the system performances. The maximum measure variation reached 1.86 mm and 7.20 mm in static and dynamic conditions, respectively, both obtained in the case of partial occlusion. This systematic source of error is likely to affect gait measures: markers placed on the patient body are often visible only by half of the cameras, with swinging arms and legs providing moving occlusions. The maximum error observed in this study can potentially affect the kinematics outcomes of conventional gait models, particularly on frontal and coronal plane, and consequently the peak muscle forces estimated with musculoskeletal models.