Unique shape variations of hind and midfoot bones in flatfoot subjects-A statistical shape modeling approach.

Flatfoot deformity is a prevalent hind- and midfoot disorder. Given its complexity, single-plane radiological measurements omit case-specific joint interaction and bone shape variations. Three-dimensional medical imaging assessment using statistical shape models provides a complete approach in characterizing bone shape variations unique to flatfoot condition. This study used statistical shape models to define specific bone shape variations of the subtalar, talonavicular, and calcaneocuboid joints that characterize flatfoot deformity, that differentiate them from healthy controls. Bones of the aforementioned joints were segmented from computed tomography scans of 40 feet. The three-dimensional hindfoot alignment angle categorized the population into 18 flatfoot subjects (≥7° valgus) and 22 controls. Statistical shape models for each joint were defined using the entire study cohort. For each joint, an average weighted shape parameter was calculated for each mode of variation, and then compared between flatfoot and controls. Significance was set at p < 0.05, with values between 0.05 ≤ p < 0.1 considered trending towards significance. The flatfoot population showed a more adducted talar head, inferiorly inclined talar neck, and posteriorly orientated medial subtalar articulation compare to controls, coupled with more navicular eversion, shallower navicular cup, and more prominent navicular tuberosity. The calcaneocuboid joint presented trends of a more adducted calcaneus, more abducted cuboid, narrower calcaneal roof, and less prominent cuboid beak compared to controls. Statistical shape model analysis identified unique shape variations which may enhance understanding and computer-aided models of the intricacies of flatfoot, leading to better diagnosis and, ultimately, surgical treatment.

Patellar cartilage increase following ACL reconstruction with and without meniscal pathology: a two-year prospective MRI morphological study.

BACKGROUND: Anterior cruciate ligament reconstruction (ACLR) together with concomitant meniscal injury are risk factors for the development of tibiofemoral (TF) osteoarthritis (OA), but the potential effect on the patellofemoral (PF) joint is unclear. The aim of this study was to: (i) investigate change in patellar cartilage morphology in individuals 2.5 to 4.5 years after ACLR with or without concomitant meniscal pathology and in healthy controls, and (ii) examine the association between baseline patellar cartilage defects and patellar cartilage volume change. METHODS: Thirty two isolated ACLR participants, 25 ACLR participants with combined meniscal pathology and nine healthy controls underwent knee magnetic resonance imaging (MRI) with 2-year intervals (baseline = 2.5 years post-ACLR). Patellar cartilage volume and cartilage defects were assessed from MRI using validated methods. RESULTS: Both ACLR groups showed patellar cartilage volume increased over 2 years (p < 0.05), and isolated ACLR group had greater annual percentage cartilage volume increase compared with controls (mean difference 3.6, 95% confidence interval (CI) 1.0, 6.3%, p = 0.008) and combined ACLR group (mean difference 2.2, 95% CI 0.2, 4.2%, p = 0.028). Patellar cartilage defects regressed in the isolated ACLR group over 2 years (p = 0.02; Z = - 2.33; r = 0.3). Baseline patellar cartilage defect score was positively associated with annual percentage cartilage volume increase (Regression coefficient B = 0.014; 95% CI 0.001, 0.027; p = 0.03) in the pooled ACLR participants. CONCLUSIONS: Hypertrophic response was evident in the patellar cartilage of ACLR participants with and without meniscal pathology. Surprisingly, the increase in patellar cartilage volume was more pronounced in those with isolated ACLR. Although cartilage defects stabilised in the majority of ACLR participants, the severity of patellar cartilage defects at baseline influenced the magnitude of the cartilage hypertrophic response over the subsequent ~ 2 years.

In-silico techniques to inform and improve the personalized prescription of shoe insoles.

Background: Corrective shoe insoles are prescribed for a range of foot deformities and are typically designed based on a subjective assessment limiting personalization and potentially leading to sub optimal treatment outcomes. The incorporation of in silico techniques in the design and customization of insoles may improve personalized correction and hence insole efficiency. Methods: We developed an in silico workflow for insole design and customization using a combination of measured motion capture, inverse musculoskeletal modelling as well as forward simulation approaches to predict the kinematic response to specific insole designs. The developed workflow was tested on twenty-seven participants containing a combination of healthy participants (7) and patients with flatfoot deformity (20). Results: Average error between measured and simulated kinematics were 4.7 ± 3.1, 4.5 ± 3.1, 2.3 ± 2.3, and 2.3 ± 2.7° for the chopart obliquity, chopart anterior-posterior axis, tarsometatarsal first ray, and tarsometatarsal fifth ray joints respectively. Discussion: The developed workflow offers distinct advantages to previous modeling workflows such as speed of use, use of more accessible data, use of only open-source software, and is highly automated. It provides a solid basis for future work on improving predictive accuracy by adapting the currently implemented insole model and incorporating additional data such as plantar pressure.

Population-based in silico modeling of anatomical shape variation of the knee and its impact on joint loading in knee osteoarthritis.

Anatomical knee joint features and osteoarthritis (OA) severity are associated, however confirming causals link to altered knee loading is challenging. This study leverages statistical shape models (SSM) to investigate the relationship between joint shape/alignment and knee loading during gait in knee OA (KOA) patients to understand their contribution to elevated medial knee loading in OA. Musculoskeletal (MSK) models were created for the mean as well as the first eight SSM principal modes of variation (-3,-2,-1, +1, +2, +3 standard deviations for each mode) and used as input to a MSK modeling framework. Using an identical KOA gait pattern (i.e., joint kinematics and ground reaction forces), we ran simulations for each MSK model and evaluated medial compartment loading magnitude and contact distribution at the instant of first and second peak of knee joint loading. An increase in external rotation, posterior tibia translation and a decrease in medial joint space and medial femoral condylar size predisposed the medial compartment knee joint to overloading during gait. This was coupled with an anterior and medial shift in contact location with increasing external rotated tibial position and increasing posterior tibial translation with respect to the femur. Next, results also highlighted a posterior shift of the medial compartment loading location with decreasing medial joint space. This study provides important population-based insights on how knee shape and alignment predispose individuals with KOA to elevated medial compartmental knee loading. This information can be crucial in assessing the risk for medial KOA development and progression.

Hindfoot kinematics and kinetics - A combined in vivo and in silico analysis approach.

BACKGROUND: The complex anatomical structure of the foot-ankle imposes challenges to accurately quantify detailed hindfoot kinematics and estimate musculoskeletal loading parameters. Most systems used to capture or estimate dynamic joint function oversimplify the anatomical structure by reducing its complexity. RESEARCH QUESTION: Can four dimensional computed tomography (4D CT) imaging in combination with an innovative foot manipulator capture in vivo hindfoot kinematics during a simulated stance phase of walking and can talocrural and subtalar articular joint mechanics be estimated based on a detailed in silico musculoskeletal foot-ankle model. METHODS: A foot manipulator imposed plantar/dorsiflexion and inversion/eversion representing a healthy stance phase of gait in 12 healthy participants while simultaneously acquiring 4D CT images. Participant-specific 3D hindfoot rotations and translations were calculated based on bone-specific anatomical coordinate systems. Articular cartilage contact area and contact pressure of the talocrural and subtalar joints were estimated using an extended foot-ankle model updated with an elastic foundation contact model upon prescribing the participant-specific rotations measured in the 4D CT measurement. RESULTS: Plantar/dorsiflexion predominantly occurred at the talocrural joint (RoM 15.9±3.9°), while inversion/eversion (RoM 5.9±3.9°) occurred mostly at the subtalar joint, with the contact area being larger at the subtalar than at the talocrural joint. Contact pressure was evenly distributed between the talocrural and subtalar joint at the beginning of the simulated stance phase but was then redistributed from the talocrural to the subtalar joint with increasing dorsiflexion. SIGNIFICANCE: In a clinical case study, the healthy participants were compared with four patients after surgically treaded intra-articular calcaneal fracture. The proposed workflow was able to detect small but meaningful differences in hindfoot kinematics and kinetics, indicative of remaining hindfoot pathomechanics that may influence the onset and progression of degenerative joint diseases.

Signatures of disease progression in knee osteoarthritis: insights from an integrated multi-scale modeling approach, a proof of concept.

Introduction: Knee osteoarthritis (KOA) is characterized by articular cartilage degeneration. It has been widely accepted that the mechanical joint environment plays a significant role in the onset and progression of this disease. In silico models have been used to study the interplay between mechanical loading and cartilage degeneration, hereby relying mainly on two key mechanoregulatory factors indicative of collagen degradation and proteoglycans depletion. These factors are the strain in collagen fibril direction (SFD) and maximum shear strain (MSS) respectively. Methods: In this study, a multi-scale in silico modeling approach was used based on a synergy between musculoskeletal and finite element modeling to evaluate the SFD and MSS. These strains were evaluated during gait based on subject-specific gait analysis data collected at baseline (before a 2-year follow-up) for a healthy and progressive early-stage KOA subject with similar demographics. Results: The results show that both SFD and MSS factors allowed distinguishing between a healthy subject and a KOA subject, showing progression at 2 years follow-up, at the instance of peak contact force as well as during the stance phase of the gait cycle. At the peak of the stance phase, the SFD were found to be more elevated in the KOA patient with the median being 0.82% higher in the lateral and 0.4% higher in the medial compartment of the tibial cartilage compared to the healthy subject. Similarly, for the MSS, the median strains were found to be 3.6% higher in the lateral and 0.7% higher in the medial tibial compartment of the KOA patient compared to the healthy subject. Based on these intersubject SFD and MSS differences, we were additionally able to identify that the tibial compartment of the KOA subject at risk of progression. Conclusion/discussion: We confirmed the mechanoregulatory factors as potential biomarkers to discriminate patients at risk of disease progression. Future studies should evaluate the sensitivity of the mechanoregulatory factors calculated based on this multi-scale modeling workflow in larger patient and control cohorts.

Hamstring harvest results in significantly reduced knee muscular protection during side-step cutting two years after anterior cruciate ligament reconstruction.

The purpose of this study was to determine the effect of donor muscle morphology following tendon harvest in anterior cruciate ligament (ACL) reconstruction on muscular support of the tibiofemoral joint during sidestep cutting. Magnetic resonance imaging (MRI) was used to measure peak cross-sectional area (CSA) and volume of the semitendinosus (ST) and gracilis (GR) muscles and tendons (bilaterally) in 18 individuals following ACL reconstruction. Participants performed sidestep cutting tasks in a biomechanics laboratory during which lower-limb electromyography, ground reaction loads, whole-body motions were recorded. An EMG driven neuro-musculoskeletal model was subsequently used to determine force from 34 musculotendinous units of the lower limb and the contribution of the ST and GR to muscular support of the tibiofemoral joint based on a normal muscle-tendon model (Standard model). Then, differences in peak CSA and volume between the ipsilateral/contralateral ST and GR were used to adjust their muscle-tendon parameters in the model followed by a recalibration to determine muscle force for 34 musculotendinous units (Adjusted model). The combined contribution of the donor muscles to muscular support about the medial and lateral compartments were reduced by 52% and 42%, respectively, in the adjusted compared to standard model. While the semimembranosus (SM) increased its contribution to muscular stabilisation about the medial and lateral compartment by 23% and 30%, respectively. This computer simulation study demonstrated the muscles harvested for ACL reconstruction reduced their support of the tibiofemoral joint during sidestep cutting, while the SM may have the potential to partially offset these reductions. This suggests donor muscle impairment could be a factor that contributes to ipsilateral re-injury rates to the ACL following return to sport.

Multi-level personalization of neuromusculoskeletal models to estimate physiologically plausible knee joint contact forces in children.

Neuromusculoskeletal models are a powerful tool to investigate the internal biomechanics of an individual. However, commonly used neuromusculoskeletal models are generated via linear scaling of generic templates derived from elderly adult anatomies and poorly represent a child, let alone children with a neuromuscular disorder whose musculoskeletal structures and muscle activation patterns are profoundly altered. Model personalization can capture abnormalities and appropriately describe the underlying (altered) biomechanics of an individual. In this work, we explored the effect of six different levels of neuromusculoskeletal model personalization on estimates of muscle forces and knee joint contact forces to tease out the importance of model personalization for normal and abnormal musculoskeletal structures and muscle activation patterns. For six children, with and without cerebral palsy, generic scaled models were developed and progressively personalized by (1) tuning and calibrating musculotendon units' parameters, (2) implementing an electromyogram-assisted approach to synthesize muscle activations, and (3) replacing generic anatomies with image-based bony geometries, and physiologically and physically plausible muscle kinematics. Biomechanical simulations of gait were performed in the OpenSim and CEINMS software on ten overground walking trials per participant. A mixed-ANOVA test, with Bonferroni corrections, was conducted to compare all models' estimates. The model with the highest level of personalization produced the most physiologically plausible estimates. Model personalization is crucial to produce physiologically plausible estimates of internal biomechanical quantities. In particular, personalization of musculoskeletal anatomy and muscle activation patterns had the largest effect overall. Increased research efforts are needed to ease the creation of personalized neuromusculoskeletal models.

Torsion Tool: An automated tool for personalising femoral and tibial geometries in OpenSim musculoskeletal models.

Common practice in musculoskeletal modelling is to use scaled musculoskeletal models based on a healthy adult, but this does not consider subject-specific geometry, such as tibial torsion and femoral neck-shaft and anteversion angles (NSA and AVA). The aims of this study were to (1) develop an automated tool for creating OpenSim models with subject-specific tibial torsion and femoral NSA and AVA, (2) evaluate the femoral component, and (3) release the tool open-source. The Torsion Tool (https://simtk.org/projects/torsiontool) is a MATLAB-based tool that requires an individual's tibial torsion, NSA and AVA estimates as input and rotates corresponding bones and associated muscle points of a generic musculoskeletal model. Performance of the Torsion Tool was evaluated comparing femur bones as personalised with the Torsion Tool and scaled generic femurs with manually segmented bones as golden standard for six typically developing children and thirteen children with cerebral palsy. The tool generated femur geometries closer to the segmentations, with lower maximum (-19%) and root mean square (-18%) errors and higher Jaccard indices (+9%) compared to generic femurs. Furthermore, the tool resulted in larger improvements for participants with higher NSA and AVA deviations. The Torsion Tool allows an automatic, fast, and user-friendly way of personalising femoral and tibial geometry in an OpenSim musculoskeletal model. Personalisation is expected to be particularly relevant in pathological populations, as will be further investigated by evaluating the effects on simulation outcomes.

ESB Clinical Biomechanics Award 2020: Pelvis and hip movement strategies discriminate typical and pathological femoral growth - Insights gained from a multi-scale mechanobiological modelling framework.

BACKGROUND: Many children with cerebral palsy (CP) develop skeletal deformities during childhood. So far, it is unknown why some children with CP develop bony deformities whereas others do not. The aims of this study were to (i) investigate what loading characteristics lead to typical and pathological femoral growth, and (ii) evaluate why some children with CP develop femoral deformities whereas other do not. METHODS: A multi-scale mechanobiological modelling workflow was used to simulate femoral growth based on three-dimensional motion capture data of six typically developing children and 16 children with CP. Based on the growth results, the participants with CP were divided into two groups: typical growth group and pathological growth group. Gait kinematics and femoral loading were compared between simulations resulting in typical growth and those resulting in pathologic growth. FINDINGS: Hip joint contact forces were less posteriorly-oriented in the pathological growth simulations compared to the typical ones. Compared to the typically developing participants, the CP group with pathological femoral growth presented increased knee flexion and no hip extension. The CP group with simulated typical growth presented similar sagittal plane joint kinematics but differed in the frontal plane pelvic and hip movement strategy, which normalized the hip joint contact force and therefore contributed to typical femoral growth trends. INTERPRETATION: Our simulation results identified specific gait features, which may contribute to pathological femoral growth. Furthermore, the hip joint contact force orientation in the sagittal plane seems to be the dominant factor for determining femoral growth simulations.