Validation of a personalized ligament-constraining discrete element framework for computing ankle joint contact mechanics.

BACKGROUND AND OBJECTIVE: Computer simulations of joint contact mechanics have great merit to improve our current understanding of articular ankle pathology. Owed to its computational simplicity, discrete element analysis (DEA) is an encouraging alternative to finite element analysis (FEA). However, previous DEA models lack subject-specific anatomy and may oversimplify the biomechanics of the ankle. The objective of this study was to develop and validate a personalized DEA framework that permits movement of the fibula and incorporates personalized cartilage thickness as well as ligamentous constraints. METHODS: A linear and non-linear DEA framework, representing cartilage as compressive springs, was established, verified, and validated. Three-dimensional (3D) bony ankle models were constructed from cadaveric lower limb CT scans imaged during application of weight (85 kg) and/or torque (10 Nm). These 3D models were used to generate cartilage thickness and ligament insertion sites based on a previously validated statistical shape model. Ligaments were modelled as non-linear tension-only springs. Validation of contact stress prediction was performed using a simple, axially constrained tibiotalar DEA model against an equivalent FEA model. Validation of ligamentous constraints compared the final position of the ankle mortise to that of the cadaver after application of torque and sequential ligament sectioning. Finally, a combined ligamentous-constraining DEA model was validated for predicted contact stress against an equivalent ligament-constraining FEA model. RESULTS: The linear and non-linear DEA model reproduced a mean articular contact stress within 0.36 MPa and 0.39 MPa of the FEA calculated stress, respectively. With respect to the ligamentous validation, the DEA ligament-balancing algorithm could reproduce the position of the distal fibula within the ankle mortise to within 0.97 mm of the experimental observed distal fibula. When combining the ligament-constraining and contact stress algorithm, DEA was able to reproduce a mean articular contact stress to within 0.50 MPa of the FEA calculated contact stress. CONCLUSION: The DEA framework presented herein offers a computationally efficient alternative to FEA for the prediction of contact stress in the ankle joint, manifesting its potential to enhance the mechanical understanding of articular ankle pathologies on both a patient-specific and population-wide level. The novelty of this model lies in its personalized nature, inclusion of the distal tibiofibular joint and the use of non-linear ligament balancing to maintain the physiological ankle joint articulation.

Principal polynomial shape analysis: A non-linear tool for statistical shape modeling.

BACKGROUND AND OBJECTIVES: The most widespread statistical modeling technique is based on Principal Component Analysis (PCA). Although this approach has several appealing features, it remains hampered by its linearity. Principal Polynomial Analysis (PPA) can capture non-linearity in a sequential algorithm, while maintaining the interesting properties of PCA. PPA is, however, computationally expensive in handling shape surface data. To this end, we propose Principal Polynomial Shape Analysis (PPSA) as an adjusted approach for non-linear shape analyzes. The aim of this study was to assess PPSA's features, its model boundaries and its general applicability. METHODS: PCA and PPSA-based shape models were investigated on one verification and three model evaluation experiments. In the verification experiment, the estimated mean of the PCA and PPSA model on a data set of synthetic lower limbs of different lengths in different poses were compared to the real mean. Further, the PCA-based and PPSA shape models were tested for three challenging cases namely for shape model creation of gait marker data, for regression analysis on aging faces and for modeling pose variation in full body scans. For the latter, additionally a Fundamental Coordinate Model (FCM) and a PPSA model on Fundamental Coordinate(FC) space was created. The performances were evaluated based on model-based accuracy, generalization, compactness and specificity. RESULTS: In the verification experiment, the scaling error reduced from 75% to below 1% when employing PPSA instead of PCA for a training set with 180° angular variation. For the model evaluation experiments, the PPSA models described the data as accurate and generalized as the PCA-based shape models. The PPSA models were slightly more compact and specific (up to 30%) than the PCA-based models. In regression, PCA and PPSA-based parameterizations explained a similar amount of variation. Lastly, for the full body scans, applying PPSA to parameterizations improved the compactness and accuracy. CONCLUSIONS: PPSA describes the non-linear relationships between principal variations in a parameterized space. Compared to standard PCA-based shape models, capturing the non-linearity reduced the nonsense information in the shape components and improved the description of the data mean.

Personalised statistical modelling of soft tissue structures in the ankle.

BACKGROUND AND OBJECTIVE: Revealing the complexity behind subject-specific ankle joint mechanics requires simultaneous analysis of three-dimensional bony and soft-tissue structures. 3D musculoskeletal models have become pivotal in orthopedic treatment planning and biomechanical research. Since manual segmentation of these models is time-consuming and subject to manual errors, (semi-) automatic methods could improve the accuracy and enlarge the sample size of personalised 'in silico' biomechanical experiments and computer-assisted treatment planning. Therefore, our aim was to automatically predict ligament paths, cartilage topography and thickness in the ankle joint based on statistical shape modelling. METHODS: A personalised cartilage and ligamentous prediction algorithm was established using geometric morphometrics, based on an 'in-house' generated lower limb skeletal model (N = 542), tibiotalar cartilage (N = 60) and ankle ligament segmentations (N = 10). For cartilage, a population-averaged thickness map was determined by use of partial least-squares regression. Ligaments were wrapped around bony contours based on iterative shortest path calculation. Accuracy of ligament path and cartilage thickness prediction was quantified using leave-one-out experiments. The novel personalised thickness prediction was compared with a constant cartilage thickness of 1.50 mm by use of a paired sample T-test. RESULTS: Mean distance error of cartilage and ligament prediction was 0.12 mm (SD 0.04 mm) and 0.54 mm (SD 0.05 mm), respectively. No significant differences were found between the personalised thickness cartilage and segmented cartilage of the tibia (p = 0.73, CI [-1.60 .10-17, 1.13 .10-17]) and talus (p = 0.95, CI[ -1.35 .10-17, 1.28 .10-17]). For the constant thickness cartilage, a statistically significant difference was found in 89% and 92% of the tibial (p < 0.001, CI [0.51, 0.58]) and talar (p < 0.001, CI [0.33, 0.40]) cartilage area. CONCLUSIONS: In this study, we described a personalised prediction algorithm of cartilage and ligaments in the ankle joint. We were able to predict cartilage and main ankle ligaments with submillimeter accuracy. The proposed method has a high potential for generating large (virtual) sample sizes in biomechanical research and mitigates technological advances in computer-assisted orthopaedic surgery.

Separating positional noise from neutral alignment in multicomponent statistical shape models.

Given sufficient training samples, statistical shape models can provide detailed population representations for use in anthropological and computational genetic studies, injury biomechanics, musculoskeletal disease models or implant design optimization. While the technique has become extremely popular for the description of isolated anatomical structures, it suffers from positional interference when applied to coupled or articulated input data. In the present manuscript we describe and validate a novel approach to extract positional noise from such coupled data. The technique was first validated and then implemented in a multicomponent model of the lower limb. The impact of noise on the model itself as well as on the description of sexual dimorphism was evaluated. The novelty of our methodology lies in the fact that no rigid transformations are calculated or imposed on the data by means of idealized joint definitions and by extension the models obtained from them.

Ischiofemoral impingement: the evolutionary cost of pelvic obstetric adaptation.

The risk for ischiofemoral impingement has been mainly related to a reduced ischiofemoral distance and morphological variance of the femur. From an evolutionary perspective, however, there are strong arguments that the condition may also be related to sexual dimorphism of the pelvis. We, therefore, investigated the impact of gender-specific differences in anatomy of the ischiofemoral space on the ischiofemoral clearance, during static and dynamic conditions. A random sampling Monte-Carlo experiment was performed to investigate ischiofemoral clearance during stance and gait in a large (n = 40 000) virtual study population, while using gender-specific kinematics. Subsequently, a validated gender-specific geometric morphometric analysis of the hip was performed and correlations between overall hip morphology (statistical shape analysis) and standard discrete measures (conventional metric approach) with the ischiofemoral distance were evaluated. The available ischiofemoral space is indeed highly sexually dimorphic and related primarily to differences in the pelvic anatomy. The mean ischiofemoral distance was 22.2 ± 4.3 mm in the females and 29.1 ± 4.1 mm in the males and this difference was statistically significant (P < 0.001). Additionally, the ischiofemoral distance was observed to be a dynamic measure, and smallest during femoral extension, and this in turn explains the clinical sign of pain in extension during long stride walking. In conclusion, the presence of a reduced ischiofemroal distance and related risk to develop a clinical syndrome of ischiofemoral impingement is strongly dominated by evolutionary effects in sexual dimorphism of the pelvis. This should be considered when female patients present with posterior thigh/buttock pain, particularly if worsened by extension. Controlled laboratory study.

Statistical Shape Modeling of Skeletal Anatomy for Sex Discrimination: Their Training Size, Sexual Dimorphism, and Asymmetry.

Purpose: Statistical shape modeling provides a powerful tool for describing and analyzing human anatomy. By linearly combining the variance of the shape of a population of a given anatomical entity, statistical shape models (SSMs) identify its main modes of variation and may approximate the total variance of that population to a selected threshold, while reducing its dimensionality. Even though SSMs have been used for over two decades, they lack in characterization of their goodness of prediction, in particular when defining whether these models are actually representative for a given population. Methods: The current paper presents, to the authors' knowledge, the most extent lower limb anatomy shape model considering the pelvis, femur, patella, tibia, fibula, talus, and calcaneum to date. The present study includes the segmented training shapes (n = 542) obtained from 271 lower limb CT scans. The different models were evaluated in terms of accuracy, compactness, generalizability as well as specificity. Results: The size of training samples needed in each model so that it can be considered population covering was estimated to approximate around 200 samples, based on the generalizability properties of the different models. Simultaneously differences in gender and patterns in left-right asymmetry were identified and characterized. Size was found to be the most pronounced sexual discriminator whereas intra-individual variations in asymmetry were most pronounced at the insertion site of muscles. Conclusion: For models aimed at population covering descriptive studies, the number of training samples required should amount a sizeable 200 samples. The geometric morphometric method for sex discrimination scored excellent, however, it did not largely outperformed traditional methods based on discrete measures.

A discrete element model to predict anatomy of the psoas muscle and path of the tendon: Design implications for total hip arthroplasty.

BACKGROUND: The accurate estimation of a muscle's line of action is a fundamental requirement in computational modelling. We present a novel anatomical muscle wrapping technique and demonstrate its clinical use on the evaluation of the Psoas muscle mechanics in hip arthroplasty. METHODS: A volume preserving, spring model to parameterize muscle anatomy changes during motion is presented. Validation was performed by a CT scan of a cadaver model in multiple positions. The predicted psoas musculotendinous path was compared with the actual imaging findings. In a second stage, psoas kinetics were compared between a conventional versus a resurfacing hip arthroplasty during gait. FINDINGS: Anatomy prediction error was found to be 2.12 mm on average (SD 1.34 mm). When applied to psoas mechanics during walking, the muscle was found to wrap predominantly around the femoral head providing a biomechanically efficient and nearly constant moment arm for flexion during the entire gait cycle. However, this advantage was found to be lost in small diameter hip arthroplasty designs resulting in an important mechanical disadvantage. The moment arm for flexion, was on average 36% (SD 0.03%) lower in the small diameter conventional hip arthroplasty as compared to the large diameter head of the hip resurfacing and this difference was highly significant. (p < 0.001). INTERPRETATION: Despite the shortcomings of an "in silico" and cadaveric study, our findings are in accordance with previous clinical and gait studies. Furthermore, the findings are strongly in favour of large diameter implant designs, warranting their further development and optimisation.

Computer-based estimation of the hip joint reaction force and hip flexion angle in three different sitting configurations.

Sitting is part of our daily work and leisure activities and can be performed in different configurations. To date, the impact of different sitting configurations on hip joint loading has not been studied. We therefore evaluated the hip joint reaction force (HJRF) and hip flexion angle in a virtual representative male Caucasian population by means of musculoskeletal modelling of three distinct sitting configurations: a simple chair, a car seat and a kneeling chair configuration. The observed median HJRF in relation to body weight and hip flexion angle, respectively, was 22.3% body weight (%BW) and 63° for the simple chair, 22.5%BW and 79° for the car seat and 8.7%BW and 50° for the kneeling chair. Even though the absolute values of HJRF are low compared to the forces generated during dynamic activities, a relative reduction of over 50% in HJRF was observed in the kneeling chair configuration. Second, the hip flexion angles were both in the kneeling chair (-29°) and simple chair configuration (-16°) lower compared to the car seat and, as such, did not reach the threshold value for femoroacetabular conflict. In conclusion, the kneeling chair appears to hold the greatest potential as an ergonomic sitting configuration for the hip joint.

Accuracy of magnetic resonance studies in the detection of chondral and labral lesions in femoroacetabular impingement: systematic review and meta-analysis.

BACKGROUND: Several types of Magnetic resonance imaging (MRI) are commonly used in imaging of femoroacetabular impingement (FAI), however till now there are no clear protocols and recommendations for each type. The aim of this meta-analysis is to detect the accuracy of conventional magnetic resonance imaging (cMRI), direct magnetic resonance arthrography (dMRA) and indirect magnetic resonance arthrography (iMRA) in the diagnosis of chondral and labral lesions in femoroacetabular impingement (FAI). METHODS: A literature search was finalized on the 17th of May 2016 to collect all studies identifying the accuracy of cMRI, dMRA and iMRA in diagnosing chondral and labral lesions associated with FAI using surgical results (arthroscopic or open) as a reference test. Pooled sensitivity and specificity with 95% confidence intervals using a random-effects meta-analysis for MRI, dMRA and iMRA were calculated also area under receiver operating characteristic (ROC) curve (AUC) was retrieved whenever possible where AUC is equivocal to diagnostic accuracy. RESULTS: The search yielded 192 publications which were reviewed according inclusion and exclusion criteria then 21 studies fulfilled the eligibility criteria for the qualitative analysis with a total number of 828 cases, lastly 12 studies were included in the quantitative meta-analysis. Meta-analysis showed that as regard labral lesions the pooled sensitivity, specificity and AUC for cMRI were 0.864, 0.833 and 0.88 and for dMRA were 0.91, 0.58 and 0.92. While in chondral lesions the pooled sensitivity, specificity and AUC for cMRI were 0.76, 0.72 and 0.75 and for dMRA were 0.75, 0.79 and 0.83, while for iMRA were sensitivity of 0.722 and specificity of 0.917. CONCLUSIONS: The present meta-analysis showed that the diagnostic test accuracy was superior for dMRA when compared with cMRI for detection of labral and chondral lesions. The diagnostic test accuracy was superior for labral lesions when compared with chondral lesions in both cMRI and dMRA. Promising results are obtained concerning iMRA but further studies still needed to fully assess its diagnostic accuracy.

The pelvifemoral rhythm in cam-type femoroacetabular impingement.

BACKGROUND: There is growing evidence that femoroacetabular impingement is a potentially important risk factor for the development of early idiopathic osteoarthritis in the nondysplastic hip. Understanding of affected joint kinematics is a basic prerequisite in the evaluation of mechanical disorders in a clinical and research oriented setting. The aim of the present study was to compare pelvifemoral kinematics between subjects diagnosed with femoroacetabular impingement and healthy controls. METHODS: The authors collected motion data of the femur and pelvis on a total of 43 hips - 19 cam impingement hips and 24 healthy controls - using a validated electromagnetic tracking device. The pelvifemoral rhythm in supine position was defined during both active and passive hip flexion and statistically compared between both groups. FINDINGS: A significant increase in posterior pelvic rotation was observed during active hip flexion in the femoroacetabular impingement group compared with the control group (P<0.001). During passive hip flexion, however, posterior pelvic rotation between the impingement group and the controls did not differ significantly (P=0.628). INTERPRETATION: Posterior pelvic rotation during active high-end hip flexion is increased in femoroacetabular impingement, indicating the presence of an active compensational mechanism that decreases the extent of harmful joint conflict during high-flexion activities.

A statistical shape model of trochlear dysplasia of the knee.

BACKGROUND: Trochlear dysplasia is known as the primary predisposing factor for patellar dislocation. Current methods to describe trochlear dysplasia are mainly qualitative or based on a limited number of discrete measurements. The purpose of this study is to apply statistical shape analysis to take the full geometrical complexity of trochlear dysplasia into account. METHODS: Statistical shape analysis was applied to 20 normal and 20 trochlear dysplastic distal femur models, including the cartilage. RESULTS: This study showed that the trochlea was anteriorized, proximalized and lateralized and that the mediolateral width and the notch width were decreased in the trochlear dysplastic femur compared to the normal femur. The first three principal components of the trochlear dysplastic femurs, accounting for 79.7% of the total variation, were size, sulcus angle and notch width. Automated classification of the trochlear dysplastic and normal femora achieved a sensitivity of 85% and a specificity of 95%. CONCLUSIONS: This study shows that shape analysis is an outstanding method to visualise the location and magnitude of shape abnormalities. Improvement of automated classification and subtyping within the trochlear dysplastic group are expected when larger training sets are used. CLINICAL RELEVANCE: Classification of trochlear dysplasia, especially borderline cases may be facilitated by automated classification. Furthermore, the identification of a decreased notch width in association with an increased sulcus angle can also contribute to the diagnosis of trochlear dysplasia.

First insights into human acetabular labrum cell metabolism.

OBJECTIVE: Hip labrum pathology has only begun to emerge as a significant source of groin pain in the last decade since the development of hip arthroscopy. Few data are available on the anatomy, histology and function of this structure. Moreover, no metabolic data exist at cellular level. The aim of this study was to characterize extracellular matrix (ECM) genes and pro-inflammatory mediators expressed by these cells. METHODS: Isolated human acetabular labrum cells were cultured in alginate beads for 10 days and additionally stimulated with interleukin (IL)-1 for 24 h. Gene expression levels and secretion of different ECM genes, enzymes and cytokines were examined by quantitative polymerase chain reaction (qPCR) and enzyme-linked immunosorbent assay (ELISA) to assess the metabolic characteristics of labrum cells. Articular chondrocytes and meniscus cells served as controls. RESULTS: Labrum cells expressed high levels of COL1A1 and low levels of COL2A1, aggrecan and SOX-9 compared to chondrocytes. However, COL2A1 was more expressed by labrum cells than by meniscus cells. The expression of matrix metalloproteinase (MMP)-1/-2/-9, ADAMTS-4 and IL-6 was significantly higher in labrum cells than in chondrocytes. IL-1 suppressed the ECM gene expression levels of labrum cells, but increased the expression levels and release of MMP-1/-3/-9/-13 and ADAMTS-4 and IL-6 by these cells. Remarkably, MMP-9 was only significantly upregulated in acetabular labrum cells. CONCLUSIONS: The findings in this study demonstrated that the acetabular labrum is populated with unique highly active fibrochondrocyte-like cells. These cells are capable of expressing and releasing pro-inflammatory enzymes and cytokines and react to a pro-inflammatory stimulus. In this way, they contribute obviously to disturbed tissue function in hip labrum pathology.

Extra-articular hip endoscopy: A review of the literature.

The aim of this review is to evaluate the current available literature evidencing on peri-articular hip endoscopy (the third compartment). A comprehensive approach has been set on reports dealing with endoscopic surgery for recalcitrant trochanteric bursitis, snapping hip (or coxa-saltans; external and internal), gluteus medius and minimus tears and endoscopy (or arthroscopy) after total hip arthroplasty. This information can be used to trigger further research, innovation and education in extra-articular hip endoscopy.

In vitro validation and reliability study of electromagnetic skin sensors for evaluation of end range of motion positions of the hip.

There is growing evidence that femoroacetabular impingement (FAI) is a probable risk factor for the development of early osteoarthritis in the nondysplastic hip. As FAI arises with end range of motion activities, measurement errors related to skin movement might be higher than anticipated when using previously reported methods for kinematic evaluation of the hip. We performed an in vitro validation and reliability study of a noninvasive method to define pelvic and femur positions in end range of motion activities of the hip using an electromagnetic tracking device. Motion data, collected from sensors attached to the bone and skin of 11 cadaver hips, were simultaneously obtained and compared in a global reference frame. Motion data were then transposed in the hip joint local coordinate systems. Observer-related variability in locating the anatomical landmarks required to define the local coordinate system and variability of determining the hip joint center was evaluated. Angular root mean square (RMS) differences between the bony and skin sensors averaged 3.2° (SD 3.5°) and 1.8° (SD 2.3°) in the global reference frame for the femur and pelvic sensors, respectively. Angular RMS differences between the bony and skin sensors in the hip joint local coordinate systems ranged at end range of motion and dependent on the motion under investigation from 1.91 to 5.81°. The presented protocol for evaluation of hip motion seems to be suited for the 3-D description of motion relevant to the experimental and clinical evaluation of femoroacetabular impingement.

Shoulder prostheses treating cuff tear arthropathy: a comparative biomechanical study.

UNLABELLED: Painful cuff tear arthropathy (CTA) affects the independence of the elderly. Surgical treatment often consists of joint replacement, the functional outcome of which remains variable. Knowledge of the biomechanical properties of the different prosthetic designs can guide the orthopaedic surgeon in the choice of implant to predict its clinical result. A 3-D computer model of the glenohumeral joint is used to analyse the moment of the deltoid muscle in the scapular plane. A geometrical 3-D ball-and-socket model of the shoulder joint was used to calculate (1) the angle-force relationships, (2) the moment arm of the deltoid muscle and (3) the moment of the deltoid muscle components, for increasing degrees of arm elevation in the scapular plane. In this 3-D model, a clinical thoraco-scapular rhythm analysis was implemented, based on measurements in normal subjects, patients treated with an anatomical prosthesis and patients treated with an inversed delta III prosthesis. These data were compared for 10 different prosthetic treatment options. RESULTS: Muscle angle-force curves show a favourable slope in non-anatomical prosthetic designs, where the centre of rotation of the glenohumeral joint is medialized, the deltoid muscle is elongated and the humeral shaft is lateralized. On the contrary, anatomical prosthetic designs do not perform well in this computer analysis. CONCLUSIONS: From a biomechanical point of view, a shoulder prosthesis which medializes the centre of rotation, lengthens the deltoid muscle and increases the deltoid lever arm, results in a significantly more powerful abduction of the shoulder, despite complete loss of rotator cuff function. RELEVANCE: This study explains why a successful functional outcome can be expected in CTA with a reversed prosthesis.