The effects of knee support on the sagittal lower-body joint kinematics and kinetics of deep squats

Little work has been done to examine the deep squat position (>130° sagittal knee flexion). In baseball and softball, catchers perform this squat an average of 146 times per nine-inning game. To alleviate some of the stress on their knees caused by this repetitive loading, some catchers wear foam knee supports.ObjectivesThis work quantifies the effects of knee support on lower-body joint kinematics and kinetics in the deep squat position.MethodsSubjects in this study performed the deep squat with no support, foam support, and instrumented support. In order to measure the force through the knee support, instrumented knee supports were designed and fabricated. We then developed an inverse dynamic model to incorporate the support loads. From the model, joint angles and moments were calculated for the three conditions.ResultsWith support there is a significant reduction in the sagittal moment at the knee of 43% on the dominant side and 63% on the non-dominant side compared to without support. These reductions are a result of the foam supports carrying approximately 20% of body weight on each side.ConclusionKnee support reduces the moment necessary to generate the deep squat position common to baseball catchers. Given the short moment arm of the patella femoral tendon, even small changes in moment can have a large effect in the tibial-femoral contact forces, particularly at deep squat angles. Reducing knee forces may be effective in decreasing incidence of osteochondritis dissecans.     

Use of distraction loading to estimate subject-specific knee ligament slack lengths

Knee ligaments guide and restrain joint motion, and their properties influence joint mechanics. Inverse modeling schemes have been used to estimate specimen-specific ligament properties, where external joint forces are assumed to balance with internal ligament and contact forces. This study simplifies this assumption by adjusting experimental loads to remove internal contact forces. The purpose of this study was to use novel experimental loading in an inverse modeling scheme to estimate ligament slack lengths, perform validation using additional loading scenarios, and evaluate sensitivity to the applied loading. Joint kinematics and kinetics were experimentally measured for a set of load cases. An optimization scheme used a specimen-specific forward kinematics model to estimate ligament slack lengths by minimizing the residual between model and experimentally measured kinetics. The calibrated model was used for a form of validation by evaluating non-optimized load cases. Additionally, uncertainty analysis related kinetic errors to previously reported kinematic errors. The six DOF tibial reactions realized RMS errors less than 23 N and 0.75 Nm for optimized load cases, and 33 N and 2.25 Nm for the non-optimized load cases. The uncertainty analysis, which was performed using the optimized load cases, showed average kinetic RMS errors less than 26 N and 0.45 Nm. The model’s recruitment patterns were similar to those found in clinical and cadaveric studies. This study demonstrated that experimental distraction loading can be used in an inverse modeling scheme to estimate ligament slack lengths with a forward kinematics model.     

Computational stability of human knee joint at early stance in Gait: Effects of muscle coactivity and anterior cruciate ligament deficiency

As one of the most complex and vulnerable structures of body, the human knee joint should maintain dynamic equilibrium and stability in occupational and recreational activities. The evaluation of its stability and factors affecting it is vital in performance evaluation/enhancement, injury prevention and treatment managements. Knee stability often manifests itself by pain, hypermobility and giving-way sensations and is usually assessed by the passive joint laxity tests. Mechanical stability of both the human knee joint and the lower extremity at early stance periods of gait (0% and 5%) were quantified here for the first time using a hybrid musculoskeletal model of the lower extremity. The roles of muscle coactivity, simulated by setting minimum muscle activation at 0–10% levels and ACL deficiency, simulated by reducing ACL resistance by up to 85%, on the stability margin as well as joint biomechanics (contact/muscle/ligament forces) were investigated. Dynamic stability was analyzed using both linear buckling and perturbation approaches at the final deformed configurations in gait. The knee joint was much more stable at 0% stance than at 5% due to smaller ground reaction and contact forces. Muscle coactivity, when at lower intensities (<3% of its maximum active force), increased dynamic stability margin. Greater minimum activation levels, however, acted as an ineffective strategy to enhance stability. Coactivation also substantially increased muscle forces, joint loads and ACL force and hence the risk of further injury and degeneration. A deficiency in ACL decreases total ACL force (by 31% at 85% reduced stiffness) and the stability margin of the knee joint at the heel strike. It also markedly diminishes forces in lateral hamstrings (by up to 39%) and contact forces on the lateral plateau (by up to 17%). Current work emphasizes the need for quantification of the lower extremity stability margin in gait.     

Elevated gastrocnemius forces compensate for decreased hamstrings forces during the weight-acceptance phase of single-leg jump landing: implications for anterior cruciate ligament injury risk

Approximately 320,000 anterior cruciate ligament (ACL) injuries in the United States each year are non-contact injuries, with many occurring during a single-leg jump landing. To reduce ACL injury risk, one option is to improve muscle strength and/or the activation of muscles crossing the knee under elevated external loading. This study?s purpose was to characterize the relative force production of the muscles supporting the knee during the weight-acceptance (WA) phase of single-leg jump landing and investigate the gastrocnemii forces compared to the hamstrings forces. Amateur male Western Australian Rules Football players completed a single-leg jump landing protocol and six participants were randomly chosen for further modeling and simulation. A three-dimensional, 14-segment, 37 degree-of-freedom, 92 muscle-tendon actuated model was created for each participant in OpenSim. Computed muscle control was used to generate 12 muscle-driven simulations, 2 trials per participant, of the WA phase of single-leg jump landing. A one-way ANOVA and Tukey post-hoc analysis showed both the quadriceps and gastrocnemii muscle force estimates were significantly greater than the hamstrings (p<0.001). Elevated gastrocnemii forces corresponded with increased joint compression and lower ACL forces. The elevated quadriceps and gastrocnemii forces during landing may represent a generalized muscle strategy to increase knee joint stiffness, protecting the knee and ACL from external knee loading and injury risk. These results contribute to our understanding of how muscle?s function during single-leg jump landing and should serve as the foundation for novel muscle-targeted training intervention programs aimed to reduce ACL injuries in sport.     

OpenSim as a preliminary kinematic testing platform for the development of total knee arthroplasty implants

The design of a total knee replacement implant needs to take account the complex surfaces of the knee which it is replacing. Ensuring design performance of the implant requires in vitro testing of the implant. A considerable amount of time is required to produce components and evaluate them inside an experimental setting. Numerous adjustments in the design of an implant and testing each individual design can be time consuming and expensive.Our solution is to use the OpenSim simulation software to rapidly test multiple design configurations of implants. This study modeled a testing rig which characterized the motion and laxity of knee implants. Three different knee implant designs were used to test and validate the accuracy of the simulation: symmetrical, asymmetric, and anatomic. Kinematics were described as distances measured from the center of each femoral condyle to a plane intersecting the most posterior points of the tibial condyles between 0 and 135° of flexion with 15° increments. Excluding the initial flexion measurement (∼0°) results, the absolute differences between all experimental and simulation results (neutral path, anterior-posterior shear, internal-external torque) for the symmetric, asymmetric, and anatomical designs were 1.98 mm ± 1.15, 1.17 mm ± 0.89, and 1.24 mm ± 0.97, respectively. Considering all designs, the accuracy of the simulation across all tests was 1.46 mm ± 1.07. It was concluded that the results of the simulation were an acceptable representation of the testing rig and hence applicable as a design tool for new total knees.     

The influence of ligament modelling strategies on the predictive capability of finite element models of the human knee joint

In finite element (FE) models knee ligaments can represented either by a group of one-dimensional springs, or by three-dimensional continuum elements based on segmentations. Continuum models closer approximate the anatomy, and facilitate ligament wrapping, while spring models are computationally less expensive. The mechanical properties of ligaments can be based on literature, or adjusted specifically for the subject. In the current study we investigated the effect of ligament modelling strategy on the predictive capability of FE models of the human knee joint. The effect of literature-based versus specimen-specific optimized material parameters was evaluated. Experiments were performed on three human cadaver knees, which were modelled in FE models with ligaments represented either using springs, or using continuum representations. In spring representation collateral ligaments were each modelled with three and cruciate ligaments with two single-element bundles. Stiffness parameters and pre-strains were optimized based on laxity tests for both approaches. Validation experiments were conducted to evaluate the outcomes of the FE models.Models (both spring and continuum) with subject-specific properties improved the predicted kinematics and contact outcome parameters. Models incorporating literature-based parameters, and particularly the spring models (with the representations implemented in this study), led to relatively high errors in kinematics and contact pressures. Using a continuum modelling approach resulted in more accurate contact outcome variables than the spring representation with two (cruciate ligaments) and three (collateral ligaments) single-element-bundle representations. However, when the prediction of joint kinematics is of main interest, spring ligament models provide a faster option with acceptable outcome.     

A subject-specific finite element musculoskeletal framework for mechanics analysis of a total knee replacement

Concurrent use of finite element (FE) and musculoskeletal (MS) modeling techniques is capable of considering the interactions between prosthetic mechanics and subject dynamics after a total knee replacement (TKR) surgery is performed. However, it still has not been performed in terms of favorable prediction accuracy and systematic experimental validation. In this study, we presented a methodology to develop a subject-specific FE-MS model of a human right lower extremity including the interactions among the subject-specific MS model, the knee joint model with ligament bundles, and the deformable FE prosthesis model. In order to evaluate its accuracy, the FE-MS model was compared with a traditional hinge-constraint MS model and experimentally verified over a gait cycle. Both models achieved good temporal agreement between the predicted muscle force and the electromyography results, though the magnitude on models is different. A higher predicted accuracy, quantified by the root-mean-square error (RMSE) and the squared Pearson correlation coefficient (r²), was found in the FE-MS model (RMSE = 177.2 N, r² = 0.90) when compared with the MS model (RMSE = 224.1 N, r² = 0.81) on the total tibiofemoral contact force. The contact mechanics, including the contact area, pressure, and stress were synchronously simulated, and the maximum contact pressure, 22.06 MPa, occurred on the medial side of the tibial insert without exceeding the yield strength of the ultra-high-molecular-weight polyethylene, 24.79 MPa. The approach outlines an accurate knee joint biomechanics analysis and provides an effective method of applying individualized prosthesis design and verification in TKR.     

Effect of knee musculature on anterior cruciate ligament strain in vivo

Squatting is a commonly prescribed exercise following reconstruction of the anterior cruciate ligament (ACL). The objective of this paper was to measure the in vivo strain patterns of the normal ACL and the load at the knee for the simple squat and for squatting with a “sport cord”. A sport cord is a large elastic rubber tube used for added resistance. Strain patterns were deduced using displacement data from a Hall Effect Strain Transducer (HEST), while joint loads were determined by a mathematical model with inputs from a force plate and electrogoniometers. ACL strain for the free squat in one subject had a maximum of 2% at a knee angle of 10° and was slack for knee angles >17°. In squatting with a sport cord, peak strain was 1% at 10° and was slack at knee angles >14°. Since these peak strains are low, squatting appears to be a safe exercise for conservative rehabilitation of ACL reconstruction patients. In addition, the sport cord is a recommended augmentation to the activity. We believe that the decrease in strain with the sport cord results from added joint stiffness due to greater compressive forces at the tibiofemoral joint. This greater compressive force results from the approximately 10% increase in quadriceps activity. From shear force data predicted by the mathematical model, the maximum anterior drawer force for free squatting (50 N) was considerably less than for sport cord squatting (430 N). Therefore, the value of shear force at the tibiofemoral joint only partially determines the load placed on the ACL.     

Comparative assessment of knee joint models used in multi-body kinematics optimisation for soft tissue artefact compensation

Estimating joint kinematics from skin-marker trajectories recorded using stereophotogrammetry is complicated by soft tissue artefact (STA), an inexorable source of error. One solution is to use a bone pose estimator based on multi-body kinematics optimisation (MKO) embedding joint constraints to compensate for STA. However, there is some debate over the effectiveness of this method. The present study aimed to quantitatively assess the degree of agreement between reference (i.e., artefact-free) knee joint kinematics and the same kinematics estimated using MKO embedding six different knee joint models. The following motor tasks were assessed: level walking, hopping, cutting, running, sit-to-stand, and step-up. Reference knee kinematics was taken from pin-marker or biplane fluoroscopic data acquired concurrently with skin-marker data, made available by the respective authors. For each motor task, Bland-Altman analysis revealed that the performance of MKO varied according to the joint model used, with a wide discrepancy in results across degrees of freedom (DoFs), models and motor tasks (with a bias between −10.2° and 13.2° and between −10.2 mm and 7.2 mm, and with a confidence interval up to ±14.8° and ±11.1 mm, for rotation and displacement, respectively). It can be concluded that, while MKO might occasionally improve kinematics estimation, as implemented to date it does not represent a reliable solution to the STA issue.     

Effects of fatigue on lower limb,pelvis and trunk kinematics and muscle activation: Gender differences

Background: Muscle fatigue is associated with biomechanical changes that may lead to anterior cruciate ligament (ACL) injuries. Alterations in trunk and pelvis kinematics may also be involved in ACL injury. Although some studies have compared the effects of muscle fatigue on lower limb kinematics between men and women, little is known about its effects on pelvis and trunk kinematics. The aim of the study was to compare the effects of fatigue on lower limb, pelvis and trunk kinematics and muscle activation between men and women during landing. Methods: The participants included forty healthy subjects. We performed kinematic analysis of the trunk, pelvis, hip and knee and muscle activation analysis of the gluteal muscles, vastus lateralis and biceps femoris, during a single-leg landing before and after fatigue. Results: Men had greater trunk flexion than women after fatigue. After fatigue, a decrease in peak knee flexion and an increase in Gmax and BF activation were observed. Conclusion: The increase in the trunk flexion can decrease the anterior tibiofemoral shear force resulted from the lower knee flexion angle, thereby decreasing the stress on the ACL.     

Sensitivity of medial and lateral knee contact force predictions to frontal plane alignment and contact locations

Musculoskeletal models are increasingly used to estimate medial and lateral knee contact forces, which are difficult to measure in vivo. The sensitivity of contact force predictions to modeling parameters is important to the interpretation and implication of results generated by the model. The purpose of this study was to quantify the sensitivity of knee contact force predictions to simultaneous errors in frontal plane knee alignment and contact locations under different dynamic conditions. We scaled a generic musculoskeletal model for N = 23 subjects’ stature and radiographic knee alignment, then perturbed frontal plane alignment and mediolateral contact locations within experimentally-possible ranges of 10° to −10° and 10 to −10 mm, respectively. The sensitivity of first peak, second peak, and mean medial and lateral knee contact forces to knee adduction angle and contact locations was modeled using linear regression. Medial loads increased, and lateral loads decreased, by between 3% and 6% bodyweight for each degree of varus perturbation. Shifting the medial contact point medially increased medial loads and decreased lateral loads by between 1% and 4% bodyweight per millimeter. This study demonstrates that realistic measurement errors of 5 mm (contact distance) or 5° (frontal plane alignment) could result in a combined 50% BW error in subject specific contact force estimates. We also show that model sensitivity varies between subjects as a result of differences in gait dynamics. These results demonstrate that predicted knee joint contact forces should be considered as a range of possible values determined by model uncertainty.     

Development of a novel MATLAB-based framework for implementing mechanical joint stability constraints within OpenSim musculoskeletal models

The Static Optimization (SO) solver in OpenSim estimates muscle activations and forces that only equilibrate applied moments. In this study, SO was enhanced through an open-access MATLAB interface, where calculated muscle activations can additionally satisfy crucial mechanical stability requirements. This Stability-Constrained SO (SCSO) is applicable to many OpenSim models and can potentially produce more biofidelic results than SO alone, especially when antagonistic muscle co-contraction is required to stabilize body joints. This hypothesis was tested using existing models and experimental data in the literature. Muscle activations were calculated by SO and SCSO for a spine model during two series of static trials (i.e. simulation 1 and 2), and also for a lower limb model (supplementary material 2). In simulation 1, symmetric and asymmetric flexion postures were compared, while in simulation 2, various external load heights were compared, where increases in load height did not change the external lumbar flexion moment, but necessitated higher EMG activations. During the tasks in simulation 1, the predicted muscle activations by SCSO demonstrated less average deviation from the EMG data (6.8% −7.5%) compared to those from SO (10.2%). In simulation 2, SO predicts constant muscle activations and forces, while SCSO predicts increases in the average activations of back and abdominal muscles that better match experimental data. Although the SCSO results are sensitive to some parameters (e.g. musculotendon stiffness), when considering the strategy of the central nervous system in distributing muscle forces and in activating antagonistic muscles, the assigned activations by SCSO are more biofidelic than SO.     

Contact pressure models for spiral phyllotaxis and their computer simulation

We present a simple, biologically motivated model for the creation of phyllotactical patterns in capitula and their computer simulation. An in-depth investigation of Ridley's contact pressure model is performed and a refinement of contact pressure between primordia on the basis of local centroidal Voronoi relaxation is presented. Using this method in combination with Hofmeisters rule for placing new primordia creates stable patterns with Fibonacci spirals of high degree for a large range of initial conditions.     

A comparison of optimisation methods and knee joint degrees of freedom on muscle force predictions during single-leg hop landings

The aim of this paper was to compare the effect of different optimisation methods and different knee joint degrees of freedom (DOF) on muscle force predictions during a single legged hop. Nineteen subjects performed single-legged hopping manoeuvres and subject-specific musculoskeletal models were developed to predict muscle forces during the movement. Muscle forces were predicted using static optimisation (SO) and computed muscle control (CMC) methods using either 1 or 3 DOF knee joint models. All sagittal and transverse plane joint angles calculated using inverse kinematics or CMC in a 1 DOF or 3 DOF knee were well-matched (RMS error<3°). Biarticular muscles (hamstrings, rectus femoris and gastrocnemius) showed more differences in muscle force profiles when comparing between the different muscle prediction approaches where these muscles showed larger time delays for many of the comparisons. The muscle force magnitudes of vasti, gluteus maximus and gluteus medius were not greatly influenced by the choice of muscle force prediction method with low normalised root mean squared errors (<48%) observed in most comparisons. We conclude that SO and CMC can be used to predict lower-limb muscle co-contraction during hopping movements. However, care must be taken in interpreting the magnitude of force predicted in the biarticular muscles and the soleus, especially when using a 1 DOF knee. Despite this limitation, given that SO is a more robust and computationally efficient method for predicting muscle forces than CMC, we suggest that SO can be used in conjunction with musculoskeletal models that have a 1 or 3 DOF knee joint to study the relative differences and the role of muscles during hopping activities in future studies.     

Effect of posterior cruciate ligament creep on muscular co-activation around knee: A pilot study

The effect of posterior cruciate ligament (PCL) on muscle co-activation (MCO) is not known though MCO has been extensively studied. The purpose of the study was to investigate the effect of PCL creep on MCO and on joint moment around the knee. Twelve males and twelve females volunteered for this study. PCL creep was estimated via tibial posterior displacement which was elicited by a 20 kg dumbbell hanged on horizontal shank near patella for 10 min. Electromyography activity from both rectus femoris and biceps femoris as well as muscle strength on the right thigh was recorded synchronically during knee isokinetic flexion–extension performance in speed of 60 deg/s as well as 120 deg/s on a dynamometer before and after PCL creep. A one-way ANOVA with repeated measures was used to evaluate the effect of creep, gender and speed. The results showed that significant tibial posterior displacement was found (p = 0.01) in both male and female groups. No significant increase of joint moment was found in flexion as well as in extension phase in both female and male groups. There was a significant effect of speed (p = 0.036) on joint moment in extension phase. Co-activation index (CI) decreased significantly (p = 0.049) in extension phase with a significant effect of gender (p ⩽ 0.001). It was concluded that creep developed in PCL due to static posterior load on the proximal tibia could significantly elicit the increase of the activation of agonist muscles but with no compensation from the antagonist in flexion as well as in extension phase. The creep significantly elicited the decrease of the antagonist–agonist CI in extension phase. MCO in females was reduced significantly in extension phase. It was suggested that PCL creep might be one of risk factors to the knee injury in sports activity.     

Robotic simulation of identical athletic-task kinematics on cadaveric limbs exhibits a lack of differences in knee mechanics between contralateral pairs

Limb asymmetry is a known factor for increased ACL injury risk. These asymmetries are normally observed during in vivo testing. Prior studies have developed in vitro testing methodologies driven by in vivo kinematics to investigate knee mechanics relative to ACL injury. The objective of this study was to determine if mechanical side-to-side asymmetries persist in contralateral pairs during in vitro simulation testing. In vivo kinematics were recorded for male and female drop vertical jump and sidestep cutting tasks. The recorded kinematics were used to robotically simulate the motions on 7 contralateral pairs of cadaveric lower extremities specimens. ACL and MCL force, torque, and strains were recorded and analyzed for differences between contralateral pairs. There was a general lack of mechanical differences between limb sides. Adduction peak torque for the male sidestep cut movement was significantly different between limb sides (p=0.04). However, this is consistent with ACL injury mechanics in that movement in the frontal plane (abduction/adduction) increases injury risk and it is possible loading differences in this plane may have resulted from tolerances within the setup process. The findings of this study indicate that contralateral knee joints were representative of each other during biomechanical in vitro tests. In future cadaveric robotic simulations, contralateral limbs can be used interchangeably. In addition, direct comparisons of the structural behaviors of isolated conditions for contralateral knee joints can be performed.     

Automatic registration of MRI-based joint models to high-speed biplanar radiographs for precise quantification of in vivo anterior cruciate ligament deformation during gait

Understanding in vivo joint mechanics during dynamic activity is crucial for revealing mechanisms of injury and disease development. To this end, laboratories have utilized computed tomography (CT) to create 3-dimensional (3D) models of bone, which are then registered to high-speed biplanar radiographic data captured during movement in order to measure in vivo joint kinematics. In the present study, we describe a system for measuring dynamic joint mechanics using 3D surface models of the joint created from magnetic resonance imaging (MRI) registered to high-speed biplanar radiographs using a novel automatic registration algorithm. The use of MRI allows for modeling of both bony and soft tissue structures. Specifically, the attachment site footprints of the anterior cruciate ligament (ACL) on the femur and tibia can be modeled, allowing for measurement of dynamic ACL deformation. In the present study, we demonstrate the precision of this system by tracking the motion of a cadaveric porcine knee joint. We then utilize this system to quantify in vivo ACL deformation during gait in four healthy volunteers.     

A method for measurement of joint kinematics in vivo by registration of 3-D geometric models with cine phase contrast magnetic resonance imaging data

A new method is presented for measuring joint kinematics by optimally matching modeled trajectories of geometric surface models of bones with cine phase contrast (cine-PC) magnetic resonance imaging data. The incorporation of the geometric bone models (GBMs) allows computation of kinematics based on coordinate systems placed relative to full 3-D anatomy, as well as quantification of changes in articular contact locations and relative velocities during dynamic motion. These capabilities are additional to those of cine-PC based techniques that have been used previously to measure joint kinematics during activity. Cine-PC magnitude and velocity data are collected on a fixed image plane prescribed through a repetitively moved skeletal joint. The intersection of each GBM with a simulated image plane is calculated as the model moves along a computed trajectory, and cine-PC velocity data are sampled from the regions of the velocity images within the area of this intersection. From the sampled velocity data, the instantaneous linear and angular velocities of a coordinate system fixed to the GBM are estimated, and integration of the linear and angular velocities is used to predict updated trajectories. A moving validation phantom that produces motions and velocity data similar to those observed in an experiment on human knee kinematics was designed. This phantom was used to assess cine-PC rigid body tracking performance by comparing the kinematics of the phantom measured by this method to similar measurements made using a magnetic tracking system. Average differences between the two methods were measured as 2.82 mm rms for anterior/posterior tibial position, and 2.63 deg rms for axial rotation. An intertrial repeatability study of human knee kinematics using the new method produced rms differences in anterior/posterior tibial position and axial rotation of 1.44 mm and 2.35 deg. The performance of the method is concluded to be sufficient for the effective study of kinematic changes caused to knees by soft tissue injuries.     

The association between patellar alignment on magnetic resonance imaging and radiographic manifestations of knee osteoarthritis

The aim of our study was to evaluate the association between patellar alignment by using magnetic resonance imaging images and radiographic manifestations of patello-femoral osteoarthritis (OA). Subjects were recruited to participate in a natural history study of symptomatic knee OA. We examined the relation of patellar alignment in the sagittal plane (patellar length ratio (PLR)) and the transverse plane (sulcus angle (SA), lateral patellar tilt angle (LPTA), and bisect offset (BO)) to radiographic features of patello-femoral OA, namely joint space narrowing and patellar osteophytes, using a proportional odds logistic regression model while adjusting for age, sex, and bone mass index (BMI). The study sample consisted of 126 males (average age 68.0 years, BMI 31.2) and 87 females (average age 64.7 years, BMI 31.6), 75% of whom had tibiofemoral OA (a Kellgren-Lawrence score of 2 or more). PLR showed a statistically significant association with joint space narrowing and osteophytosis in the lateral compartment. SA showed significant association with medial joint space narrowing and with lateral and medial patellar osteophytosis. LPTA and BO showed significant association with both radiographic indices of the lateral compartment. Clear linear trends were found in association between PLR, LPTA and BO, and with outcomes associated with lateral patello-femoral OA. SA, LPTA, and BO showed linear trends of association with medial joint space narrowing. Results of our study clearly suggest the association between indices of patellar alignment and such features of patello-femoral OA as osteophytosis and joint space narrowing. Additional studies will be required to establish the normal and abnormal ranges of patellar alignment indices and their longitudinal relation to patello-femoral OA.