Objective: Herewith, we report the development of Orthopedic Digital Image Analysis (ODIA) software that is developed to obtain quantitative measurements of knee osteoarthritis (OA) radiographs automatically. Manual segmentation and measurement of OA parameters currently hamper large-cohort analyses, and therefore, automated and reproducible methods are a valuable addition in OA research. This study aims to test the automated ODIA measurements and compare them with available manual Knee Imaging Digital Analysis (KIDA) measurements as comparison. Design: This study included data from the CHECK (Cohort Hip and Cohort Knee) initiative, a prospective multicentre cohort study in the Netherlands with 1,002 participants. Knee radiographs obtained at baseline of the CHECK cohort were included and mean medial/lateral joint space width (JSW), minimal JSW, joint line convergence angle (JLCA), eminence heights, and subchondral bone intensities were compared between ODIA and KIDA. Results: Of the potential 2,004 radiographs, 1,743 were included for analyses. Poor intraclass correlation coefficients (ICCs) were reported for the JLCA (0.422) and minimal JSW (0.299). The mean medial and lateral JSW, eminence height, and subchondral bone intensities reported a moderate to good ICC (0.7 or higher). Discrepancies in JLCA and minimal JSW between the 2 methods were mostly a problem in the lateral tibia plateau. Conclusions: The current ODIA tool provides important measurements of OA parameters in an automated manner from standard radiographs of the knee. Given the automated and computerized methodology that has very high reproducibility, ODIA is suitable for large epidemiological cohorts with various follow-up time points to investigate structural progression, such as CHECK or the Osteoarthritis Initiative (OAI).

Introduction: To date the aetiology of femoroacetabular impingement (FAI) is still not completely understood. There are mechanical theories that suggest symptomatic FAI is linked to sagittal pelvic morphology and spinopelvic-femoral dynamics. The aim of this study is to evaluate the relation of sagittal pelvic morphology and orientation to radiographic signs of FAI. Additionally, we test whether the relation between FAI and spinopelvic parameters differs in osteoarthritic hips. Methods: From a prospective, observational cohort study, 1002 patients between 45 and 65 years old with a first episode of knee or hip pain were followed for 8 years. All patients who had lateral lumbar radiographs and clinical and radiographic follow-up of the hips were included in the present study. Range of internal rotation of the hip as well as radiographic signs of FAI (alpha and Wiberg angle) and presence of hip osteoarthritis (Kellgren and Lawrence) were systematically measured at baseline. Pelvic incidence (PI), pelvic tilt (PT), sacral slope (SS)) were measured at 8-year follow-up. Associations between PI, PT, SS and FAI parameters were tested using generalised estimating equations. Results: 421 subjects, 842 hips, were included. No significant relations between PI, PT or SS and alpha or Wiberg angle were found. Comparison of hips with and without radiological sign(s) of FAI showed no differences in PI, PT or SS. There was no relation between range of internal rotation of the hip and spinopelvic parameters. Conclusion: Sagittal pelvic morphology and orientation are not related to the presence of radiological signs of FAI in this study population.

Magnetic resonance Imaging is the gold standard for assessment of soft tissues; however, X-ray-based techniques are required for evaluating bone-related pathologies. This study evaluated the performance of synthetic computed tomography (sCT), a novel MRI-based bone visualization technique, compared with CT, for the scoring of knee osteoarthritis. sCT images were generated from the 3T T1-weighted gradient-echo MR images using a trained machine learning algorithm. Two readers scored the severity of osteoarthritis in tibiofemoral and patellofemoral joints according to OACT, which enables the evaluation of osteoarthritis, from its characteristics of joint space narrowing, osteophytes, cysts and sclerosis in CT (and sCT) images. Cohen's κ was used to assess the interreader agreement for each modality, and intermodality agreement of CT- and sCT-based scores for each reader. We also compared the confidence level of readers for grading CT and sCT images using confidence scores collected during grading. Inter-reader agreement for tibiofemoral and patellofemoral joints were almost-perfect for both modalities (κ = 0.83–0.88). The intermodality agreement of osteoarthritis scores between CT and sCT was substantial to almost-perfect for tibiofemoral (κ = 0.63 and 0.84 for the two readers) and patellofemoral joints (κ = 0.78 and 0.81 for the two readers). The analysis of diagnosis confidence scores showed comparable visual quality of the two modalities, where both are showing acceptable confidence levels for scoring OA. In conclusion, in this single-center study, sCT and CT were comparable for the scoring of knee OA.

Femoral component orientation plays a key role in implant stability and therefore the success rate of total hip arthroplasty. To date, this topic has been studied using various definitions and a variety of imaging modalities and protocols. The aim of this study is a proof of concept that a new algorithm can be used to describe the femoral component's 3D orientation on the three orthogonal anatomical planes and relative to its mechanical axis using input from two orthogonal planes. CT scans of 18 patients with a total of 22 hip arthroplasties were collected. From these, orthogonal coronal and sagittal projections of the complete femur were acquired in the scanning position (MIPs) and relative to the femoral mechanical axis (corrected MIPs). On these images, the orientation of the neck of the femoral component in space and relative to the femoral axis, respectively, was measured by coronal inclination (CIF), sagittal inclination (SIF) and transverse version (TVF). With the algorithm, TVF was also calculated based on CIF and SIF. Differences between measured and calculated TVF and intra- and inter-observer reliability were evaluated using intra-class correlation coefficients (ICC). The error of non-orthogonal imaging (85° angle between the sagittal and coronal reconstructions) was tested on a third series of MIPs. The ICC between the calculated TVF and manually measured TVF, in space and relative to the femoral axis, was 0.98 for both with median absolute differences of 1.3 and 1.5°. For non-orthogonal images this was 0.70 with a median absolute difference of 5°. ICCs for intra-observer and inter-observer reliability for the calculated TVF values were 0.98 and 0.88, respectively. With this algorithm the transverse orientation of the neck of the femoral component can be assessed in space and relative to the mechanical femoral axis by combining its sagittal and coronal orientation. As long as the imaging visualizes two orthogonal planes, the orientation of an implant can be assessed in 3-D, regardless of the imaging modality.

BACKGROUND: The introduction of stem cell transplantation has improved life expectancy and cognitive outcome in patients with mucopolysaccharidosis I, but this condition remains associated with substantial residual disease in several parts of the body. Many patients have hip dysplasia with progressive medial flattening of the femoral head. Quantitative evidence on the effect of surgery on remodeling to sphericity of flattened femoral heads is lacking. In the present study, we used statistical shape modeling to quantify the effect of hip surgery on the sphericity of the femoral head in patients with mucopolysaccharidosis I. METHODS: We performed a retrospective case control study involving a series of 23 patients with hip dysplasia due to mucopolysaccharidosis I. Surgery was not offered to the first 11 children (control group). Following a change in treatment protocol, the next 12 children underwent bilateral proximal femoral varus derotation osteotomy and Pemberton osteotomy for the treatment of acetabular dysplasia with progressive femoral head flattening (surgery group). The surgery and control groups were compared with a reference group of patients with normal hips. Statistical shape modeling was used to quantify the shape of the femoral head (i.e., flattening and/or roundness of the epiphysis). RESULTS: The mean age at the time of stem cell transplantation in the surgery and control groups was comparable (1.2 years). The mean age at the time of surgical intervention was 5.5 years, and mean duration of postoperative follow-up was 3.3 years. Statistical shape modeling showed variations within the total group in terms of medial indentation, width, height, and sphericity of the femoral heads. In contrast to the progressive femoral head flattening in the control group, the surgery group showed improvement of the sphericity of the femoral head after surgery. The overall shape characteristics of the femoral head in the surgery group were similar to those of the reference group of patients with normal hips. CONCLUSIONS: To our knowledge, this is the first study in patients with mucopolysaccharidosis I that has shown quantitative remodeling of the dysplastic, flattened femoral head to normal sphericity after increasing containment of the femoral head. LEVEL OF EVIDENCE: Therapeutic Level III. See Instructions for Authors for a complete description of levels of evidence.

AIMS: The aetiologies of common degenerative spine, hip, and knee pathologies are still not completely understood. Mechanical theories have suggested that those diseases are related to sagittal pelvic morphology and spinopelvic-femoral dynamics. The link between the most widely used parameter for sagittal pelvic morphology, pelvic incidence (PI), and the onset of degenerative lumbar, hip, and knee pathologies has not been studied in a large-scale setting. METHODS: A total of 421 patients from the Cohort Hip and Cohort Knee (CHECK) database, a population-based observational cohort, with hip and knee complaints < 6 months, aged between 45 and 65 years old, and with lateral lumbar, hip, and knee radiographs available, were included. Sagittal spinopelvic parameters and pathologies (spondylolisthesis and degenerative disc disease (DDD)) were measured at eight-year follow-up and characteristics of hip and knee osteoarthritis (OA) at baseline and eight-year follow-up. Epidemiology of the degenerative disorders and clinical outcome scores (hip and knee pain and Western Ontario and McMaster Universities Osteoarthritis Index) were compared between low PI (< 50°), normal PI (50° to 60°), and high PI (> 60°) using generalized estimating equations. RESULTS: Demographic details were not different between the different PI groups. L4 to L5 and L5 to S1 spondylolisthesis were more frequently present in subjects with high PI compared to low PI (L4 to L5, OR 3.717; p = 0.024 vs L5 to S1 OR 7.751; p = 0.001). L5 to S1 DDD occurred more in patients with low PI compared to high PI (OR 1.889; p = 0.010), whereas there were no differences in L4 to L5 DDD among individuals with a different PI. The incidence of hip OA was higher in participants with low PI compared to normal (OR 1.262; p = 0.414) or high PI (OR 1.337; p = 0.274), but not statistically different. The incidence of knee OA was higher in individuals with a high PI compared to low PI (OR 1.620; p = 0.034). CONCLUSION: High PI is a risk factor for development of spondylolisthesis and knee OA. Low pelvic incidence is related to DDD, and may be linked to OA of the hip. Level of Evidence: 1b Cite this article: Bone Joint J 2020;102-B(9):1261-1267.

Objective: To design an automated workflow for hip radiographs focused on joint shape and tests its prognostic value for future hip osteoarthritis. Design: We used baseline and 8-year follow-up data from 1,002 participants of the CHECK-study. The primary outcome was definite radiographic hip osteoarthritis (rHOA) (Kellgren–Lawrence grade ≥2 or joint replacement) at 8-year follow-up. We designed a method to automatically segment the hip joint from radiographs. Subsequently, we applied machine learning algorithms (elastic net with automated parameter optimization) to provide the Shape-Score, a single value describing the risk for future rHOA based solely on joint shape. We built and internally validated prediction models using baseline demographics, physical examination, and radiologists scores and tested the added prognostic value of the Shape-Score using Area-Under-the-Curve (AUC). Missing data was imputed by multiple imputation by chained equations. Only hips with pain in the corresponding leg were included. Results: 84% were female, mean age was 56 (±5.1) years, mean BMI 26.3 (±4.2). Of 1,044 hips with pain at baseline and complete follow-up, 143 showed radiographic osteoarthritis and 42 were replaced. 91.5% of the hips had follow-up data available. The Shape-Score was a significant predictor of rHOA (odds ratio per decimal increase 5.21, 95%-CI (3.74–7.24)). The prediction model using demographics, physical examination, and radiologists scores demonstrated an AUC of 0.795, 95%-CI (0.757–0.834). After addition of the Shape-Score the AUC rose to 0.864, 95%-CI (0.833–0.895). Conclusions: Our Shape-Score, automatically derived from radiographs using a novel machine learning workflow, may strongly improve risk prediction in hip osteoarthritis.

Background and purpose — Being able to predict the hip–knee–ankle angle (HKAA) from standard knee radiographs allows studies on malalignment in cohorts lacking full-limb radiography. We aimed to develop an automated image analysis pipeline to measure the femoro-tibial angle (FTA) from standard knee radiographs and test various FTA definitions to predict the HKAA. Patients and methods — We included 110 pairs of standard knee and full-limb radiographs. Automatic search algorithms found anatomic landmarks on standard knee radiographs. Based on these landmarks, the FTA was automatically calculated according to 9 different definitions (6 described in the literature and 3 newly developed). Pearson and intra-class correlation coefficient [ICC]) were determined between the FTA and HKAA as measured on full-limb radiographs. Subsequently, the top 4 FTA definitions were used to predict the HKAA in a 5-fold cross-validation setting. Results — Across all pairs of images, the Pearson correlations between FTA and HKAA ranged between 0.83 and 0.90. The ICC values from 0.83 to 0.90. In the cross-validation experiments to predict the HKAA, these values decreased only minimally. The mean absolute error for the best method to predict the HKAA from standard knee radiographs was 1.8° (SD 1.3). Interpretation — We showed that the HKAA can be automatically predicted from standard knee radiographs with fair accuracy and high correlation compared with the true HKAA. Therefore, this method enables research of the relationship between malalignment and knee pathology in large (epidemiological) studies lacking full-limb radiography.

Objective: To assess the ability of radiography-based bone texture variables in proximal femur and acetabulum to predict incident radiographic hip osteoarthritis (rHOA) over a 10 years period. Design: Pelvic radiographs from CHECK at baseline (987 hips) were analyzed for bone texture using fractal signature analysis (FSA) in proximal femur and acetabulum. Elastic net (machine learning) was used to predict the incidence of rHOA (including Kellgren–Lawrence grade (KL) ≥ 2 or total hip replacement (THR)), joint space narrowing score (JSN, range 0–3), and osteophyte score (OST, range 0–3) after 10 years. Performance of prediction models was assessed using the area under the receiver operating characteristic curve (ROC AUC). Results: Of the 987 hips without rHOA at baseline, 435 (44%) had rHOA at 10-year follow-up. Of the 667 hips with JSN grade 0 at baseline, 471 (71%) had JSN grade ≥ 1 at 10-year follow-up. Of the 613 hips with OST grade 0 at baseline, 526 (86%) had OST grade ≥ 1 at 10-year follow-up. AUCs for the models including age, gender, and body mass index (BMI) to predict incident rHOA, JSN, and OST were 0.59, 0.54, and 0.51, respectively. The inclusion of bone texture variables in the models improved the prediction of incident rHOA (ROC AUC 0.68 and 0.71 when baseline KL was also included in the model) and JSN (ROC AUC 0.62), but not incident OST (ROC AUC 0.52). Conclusion: Bone texture analysis provides additional information for predicting incident rHOA or THR over 10 years.

The bones forming the talocrural joint (TCJ) and subtalar joint (STJ) are often assumed to be bilaterally symmetric. Therefore, the contralateral limb (i.e. the fibula, tibia, calcaneus and talus) is used as a template or an intra-subject control in clinical and research practice. However, the validity of the symmetry assumption is controversial, because insufficient information is available on the shape variations and bilateral (a)symmetry of the fibula, tibia, calcaneus and talus. Using three-dimensional spatially dense sampled representations of bone shapes extracted from bilateral computed tomography scans of 66 individuals (55 male, mean age: 61 ± 10 years; 11 female, mean age: 53 ± 15 years), we analyzed whether: (i) similar shape patterns exist in the left and right bones of the same type; (ii) gender has an effect on bone shape variations; (iii) intra-subject shape variation is smaller than that of inter-subject for a given shape variance direction. For the first set of analyses, all left and right instances of the same type of bone were considered as two separate groups, and statistically compared with each other on multiple aspects including group location (central tendency), variance-covariance scale (dispersion) and orientation (covariance structure) using distance-based permutational tests. For the second and third sets of analyses, all left and right bones of the same type were pooled into one group, and shape variations in the TCJ and STJ bones were extracted using principal component analysis. The effects of gender on age-adjusted bone shape differences were assessed using an analysis of covariance. Moreover, intra-class correlation was employed to evaluate intra- and inter-subject bone shape variations. For each bone type, both sides had similar shape patterns (P_{permutational}-values > 0.05). After Bonferroni adjustment, gender led to shape differences, which were mainly in the lateral and medial condyles of the tibia (P = 0.003), the length and height of the calcaneus (P < 0.001), the posterior and anterior talar articular surfaces of the calcaneus (P = 0.001), and in the posterior aspect of the talus (P = 0.001). Intra-subject shape variations in the tibial tuberosity together with the diameter of the tibia, and the curvature of the fibula shaft and the diameter of the fibula were as high as those of inter-subject. This result suggests that the shape symmetry assumption could be violated for some specific shape variations in the fibula and tibia.

Objectives: There is sparse evidence for a relationship between cardiovascular disease (CVD) and osteoarthritis (OA). We investigated the association between incidence of arterial calcifications and incidence of radiographic knee and/or hip OA. Design: We used baseline and 8-year follow-up data of Cohort Hip and Cohort Knee (CHECK). Knees and hips were either Kellgren-Lawrence (KL) grade 0 or 1 at baseline. Arterial calcifications were scored on hip and knee radiographs using a four-grade scale. Scores were summed for patient-level analyses. To investigate incidence, participants with arterial calcifications at baseline or missing follow-up were excluded. Incident OA was defined per joint as KL ≥ 2 or prosthesis at year eight. The association between incidenct of arterial calcifications and incident OA was studied using mixed-effects logistic regression. Results: Of 763 participants included, 623 (82%) were women. Mean (sd) age was 56 (5.1) years, mean (sd) body mass index (BMI) 26.2 (4.1) kg/m². Arterial calcifications developed in 174 participants (283 joints). OA developed in 456 participants (778 joints). Sex modified the association between arterial calcification and OA. In women, incident arterial calcification around a joint was positively associated with incident OA in that joint (adjusted OR 2.51 (95% CI 1.57-4.03)). In men, no association was observed on joint-level, but at patient-level the arterial calcification sum score was negatively associated with incident OA (adjusted OR per point increase 0.70 (95% CI 0.54-0.90)) indicating a systemic effect. Conclusions: We observed sex-dependent associations between incident arterial calcification and incident radiographic knee and/or hip OA, which differs between joint- and patient-level.