Predicting OA progression to total hip replacement: can we do better than risk factors alone using active shape modelling as an imaging biomarker?

Development of Local Software for Automatic Measurement of Geometric Parameters in the Proximal Femur Using a Combination of a Deep Learning Approach and an Active Shape Model on X-ray Images

Proximal femur geometry is an important risk factor for diagnosing and predicting hip and femur injuries. Hence, the development of an automated approach for measuring these parameters could help physicians with the early identification of hip and femur ailments. This paper presents a technique that combines the active shape model (ASM) and deep learning methodologies. First, the femur boundary is extracted by a deep learning neural network. Then, the femur's anatomical landmarks are fitted to the extracted border using the ASM method. Finally, the geometric parameters of the proximal femur, including femur neck axis length (FNAL), femur head diameter (FHD), femur neck width (FNW), shaft width (SW), neck shaft angle (NSA), and alpha angle (AA), are calculated by measuring the distances and angles between the landmarks. The dataset of hip radiographic images consisted of 428 images, with 208 men and 220 women. These images were split into training and testing sets for analysis. The deep learning network and ASM were subsequently trained on the training dataset. In the testing dataset, the automatic measurement of FNAL, FHD, FNW, SW, NSA, and AA parameters resulted in mean errors of 1.19%, 1.46%, 2.28%, 2.43%, 1.95%, and 4.53%, respectively.

Association of Bone Shape and Alignment Analyzed Using Statistical Shape Modeling With Severity of First Metatarsophalangeal Joint Osteoarthritis

OBJECTIVE

We aimed to explore the relationship between bone shape and radiographic severity in individuals with first metatarsophalangeal joint osteoarthritis (first MTP joint OA).

METHODS

Weightbearing lateral and dorsoplantar radiographs were obtained for the symptomatic foot of 185 participants (105 females, aged 22 to 85 years) with clinically diagnosed first MTP joint OA. Participants were classified into none/mild, moderate, or severe categories using a standardized atlas. An 80-point model for lateral radiographs and 77-point model for dorsoplantar radiographs was used to define independent modes of variation using statistical shape modeling software. Odds ratios adjusted for confounders were calculated using ordinal regression to determine the association between radiographic severity and mode scores.

RESULTS

After assessment and grading of radiographs, 35 participants (18.9%) were included in the none/mild first MTP joint OA severity category, 69 (37.2%) in the moderate severity category, and 81 (43.7%) in the severe category. For lateral-view radiographs, 16 modes of variation were included, which collectively represented 83.2% of total shape variance. Of these, four modes were associated with radiographic severity. For dorsoplantar-view radiographs, 15 modes of variation were included, representing 82.6% of total shape variance. Of these, six modes were associated with radiographic severity.

CONCLUSIONS

Variations in the shape and alignment of the medial cuneiform, first metatarsal, and proximal and distal phalanx of the hallux are significantly associated with radiographic severity of first MTP joint OA. Prospective studies are required to determine whether bone shape characteristics are associated with the development and/or progression of this condition.

Analysis of cam location characteristics in FAI syndrome patients from 3D MR images demonstrates sex-specific differences

Cam femoroacetabular impingement (FAI) syndrome is associated with hip osteoarthritis (OA) development. Hip shape features, derived from statistical shape modeling (SSM), are predictive for OA incidence, progression, and arthroplasty. Currently, no three-dimensional (3D) SSM studies have investigated whether there are cam shape differences between male and female patients, which may be of potential clinical relevance for FAI syndrome assessments. This study analyzed sex-specific cam location and shape in FAI syndrome patients from clinical magnetic resonance examinations (M:F 56:41, age: 16-63 years) using 3D focused shape modeling-based segmentation (CamMorph) and partial least squares regression to obtain shape features (latent variables [LVs]) of cam morphology. Two-way analysis of variance tests were used to assess cam LV data for sex and cam volume severity differences. There was no significant interaction between sex and cam volume severity for the LV data. A sex main effect was significant for LV 1 (cam size) and LV 2 (cam location) with medium to large effect sizes (p < 0.001, d > 0.75). Mean results revealed males presented with a superior-focused cam, whereas females presented with an anterior-focused cam. When stratified by cam volume, cam morphologies were located superiorly in male and anteriorly in female FAI syndrome patients with negligible, mild, or moderate cam volumes. Both male and female FAI syndrome patients with major cam volumes had a global cam distribution. In conclusion, sex-specific cam location differences are present in FAI syndrome patients with negligible, mild, and moderate cam volumes, whereas major cam volumes were globally distributed in both male and female patients.

Intergenic Interactions of SBNO1, NFAT5 and GLT8D1 Determine the Susceptibility to Knee Osteoarthritis among Europeans of Russia

This study was conducted to examine the associations between genome-wide association studies (GWAS)-important single nucleotide polymorphisms (SNPs) and knee osteoarthritis (KOA) among Europeans of Russia. The present replicative study ("patient-control" design has been used) was carried out on 1000 DNA samples from KOA (n = 500) and KOA-free (n = 500) participants. Ten GWAS-important for KOA SNPs of eight candidate genes (LYPLAL1, GNL3, GLT8D1, SBNO1, WWP2, NFAT5, TGFA, GDF5) were studied. To assess the link between SNPs and KOA susceptibility, logistic regression (to establish independent SNP effects) and MB-MDR (to identify SNP-SNP interactions) were used. As a result of this genetic analysis, the associations of individual SNPs with KOA have not been proven. Eight loci out of ten tested SNPs interacted with each other (within twelve genetic models) and determined susceptibility to KOA. The greatest contribution to the disease development were made by three polymorphisms/genes such as rs6976 (C>T) GLT8D1, rs56116847 (G>A) SBNO1, rs6499244 (T>A) NFAT5 (each was included in 2/3 [8 out 12] KOA-responsible genetic interaction models). A two-locus epistatic interaction of rs56116847 (G >A) SBNO1 × rs6499244 (T>A) NFAT5 determined the maximum percentage (0.86%) of KOA entropy. KOA-associated SNPs are regulatory polymorphisms that affect the expression/splicing level, epigenetic modification of 72 genes in KOA-pathogenetically significant organs such as skeletal muscles, tibial arteries/nerves, thyroid, adipose tissue, etc. These putative KOA-effector genes are mainly involved in the organization/activity of the exoribonuclease complex and antigen processing/presentation pathways. In conclusion, KOA susceptibility among Europeans of Russia is mediated by intergenic interactions (but not the main effects) of GWAS-important SNPs.

Prediction of femoral head coverage from articulated statistical shape models of patients with developmental dysplasia of the hip

Developmental dysplasia of the hip (DDH) is commonly described as reduced femoral head coverage due to anterolateral acetabular deficiency. Although reduced coverage is the defining trait of DDH, more subtle and localized anatomic features of the joint are also thought to contribute to symptom development and degeneration. These features are challenging to identify using conventional approaches. Herein, we assessed the morphology of the full femur and hemi-pelvis using an articulated statistical shape model (SSM). The model determined the morphological and pose-based variations associated with DDH in a population of Japanese females and established which of these variations predict coverage. Computed tomography (CT) images of 83 hips from 47 patients were segmented for input into a correspondence-based SSM. The dominant modes of variation in the model initially represented scale and pose. After removal of these factors through individual bone alignment, femoral version and neck-shaft angle, pelvic curvature, and acetabular version dominated the observed variation. Femoral head oblateness and prominence of the acetabular rim and various muscle attachment sites of the femur and hemi-pelvis were found to predict 3D CT-based coverage measurements (R2  = 0.5-0.7 for the full bones, R2  = 0.9 for the joint). Statement of Clinical Significance: Currently, clinical measurements of DDH only consider the morphology of the acetabulum. However, the results of this study demonstrated that variability in femoral head shape and several muscle attachment sites were predictive of femoral head coverage. These morphological differences may provide insight into improved clinical diagnosis and surgical planning based on functional adaptations of patients with DDH.

A prospective cohort study on cam morphology and its role in progression of osteoarthritis

Background

Cam morphology contributes to the development of hip osteoarthritis (OA) but is less studied in the general population. This study describes its associations with clinical and imaging features of hip OA.

Methods

Anteroposterior hip radiographs of 1019 participants from the Tasmanian Older Adult Cohort (TASOAC) were scored at baseline for α angle (cam morphology) in both hips. Using the Altman's atlas, radiographic hip OA (ROA) was assessed at baseline. Hip pain and right hip structural changes were assessed on a subset of 245 magnetic resonance images (MRI) at 5 years. Joint registry data for total hip replacement (THR) was acquired 14 years from baseline.

Results

Of 1906 images, cam morphology was assessed in 1016 right and 890 left hips. Cross‐sectionally, cam morphology modestly associated with age (prevalence ratio [PR]: 1.02 P = .03) and body mass index (BMI) (PR: 1.03‐1.07, P = .03) and strongly related to male gender (PR: 2.96, P < .001). Radiographically, cam morphology was prevalent in those with decreased joint space (PR: 1.30 P = .03) and osteophytes (PR: 1.47, P = .03). Longitudinally, participants with right cam and high BMI had more hip pain (PR: 17.9, P = .02). At the end of 5 years of follow‐up these participants were also more likely to have structural changes such as bone marrow lesions (BMLs) (PR: 1.90 P = .04), cartilage defects (PR: 1.26, P = .04) and effusion‐synovitis at multiple sites (PR: 1.25 P = .02). Cam morphology at baseline in either hip predicted up to threefold risk of THR (PR: 3.19, P = .003) at the end of 14 years.

Conclusion

At baseline, cam morphology was linked with age, higher weight, male gender, early signs of radiographic OA such as joint space narrowing (JSN) and osteophytes (OST). At follow‐up, cam predicted development of hip BMLs, hip effusion‐synovitis, cartilage damage and THR. These findings suggest that cam morphology plays a significant role in early OA and can be a precursor or contribute to hip OA in later life.

Statistical Shape Modeling of US Images to Predict Hip Dysplasia Development in Infants

Background The current widely applied Graf classification used on US images for developmental dysplasia of the hip in infants does not enable prediction of the development and outcome of well-centered stable dysplastic hips (Graf type II). Purpose To use statistical shape modeling on US images to identify acetabular shape characteristics of Graf type II hips, which enable prediction of the development of Graf type II hips, and to identify which hips benefit from Pavlik harness treatment. Materials and Methods In this secondary analysis of a prospective multicenter randomized trial on treatment of 104 infants aged 3-4 months with Graf type IIb or IIc hip dysplasia conducted between 2009 and 2015, a statistical shape model was developed on baseline US images. With multivariable logistic regression adjusted for infant sex and treatment (Pavlik harness treatment vs active observation), shape modes were correlated with the outcomes of persistent hip dysplasia on US images (α angle <60°) after 12-week follow-up and residual hip dysplasia on pelvic radiographs (Tönnis classification: acetabular index greater than 2 standard deviations) around 1 year of age. An interaction term (treatment with mode) was used to investigate if this result depended on treatment. Results Baseline US images were available in 97 infants (mean age, 3.37 years ± 0.43 [standard deviation]; 89 [92%] girls; 90 cases of Graf type IIb hip dysplasia; 52 cases treated with Pavlik harness). Shape modes 2 and 3 of the statistical shape modeling were associated with persistent hip dysplasia on US images (odds ratio [OR] = 0.43; P = .007 and OR = 2.39; P = .02, respectively). Mode 2 was also associated with residual hip dysplasia on pelvic radiographs (OR = 0.09; P = .002). The interaction term remained significant after multivariable analysis, indicating that Pavlik harness treatment was beneficial in patients with negative mode 2 values (OR = 12.46; P = .01). Conclusion Statistical shape modeling of US images of infants with Graf type II dysplastic hips predicted which hips developed to normal or remained dysplastic and identified hips that benefited from Pavlik harness treatment. © RSNA, 2022.

Predictors of total hip replacement in community based older adults: a cohort study

OBJECTIVE

The purpose of this study is to describe predictors of total hip replacement (THR) in community dwelling older adults. A better understanding of predictors of THR can aid in triaging patients and researching preventative strategies.

DESIGN

At baseline, participants had assessment of radiographic OA and cam morphology (from pelvic radiographs), shape mode scores and hip bone mineral density (BMD; from dual energy X-ray absorptiometry (DXA)). After 2.6 and 5 years, participants reported hip pain using WOMAC (Western Ontario and McMaster Universities Osteoarthritis Index), and had hip structural changes assessed using magnetic resonance imaging (MRI). Risk of THR was analysed using mixed-effect Poisson regression.

RESULTS

Incidence of THR for OA over 14 years was 4.6% (37/801). As expected, WOMAC hip pain and hip radiographic OA both predicted risk of THR. Additionally, shape mode 2 score (decreasing acetabular coverage) (RR 1.83/SD; 95% CI 1.1-3.04), shape mode 4 score (non-spherical femoral head) (RR 0.59/SD; 95% CI 0.36-0.96), cam morphology (α > 60°) (RR 2.2/SD; 95% CI 1.33-3.36), neck of femur BMD (RR 2.09/SD, 95% CI 1.48-2.94) and bone marrow lesions (BMLs) increased risk of THR (RR 7.10/unit; 95% CI 1.09-46.29).

CONCLUSION

In addition to hip pain and radiographic hip OA, measures of hip shape, cam morphology, BMD and BMLs independently predict risk of THR. This supports the role of hip bone geometry and structure in the pathogenesis of end stage hip OA and has identified factors that can be used to improve prediction models for THR.

Statistical shape modeling of the hip and the association with hip osteoarthritis: a systematic review

OBJECTIVE

To summarize available evidence on the association between hip shape as quantified by statistical shape modeling (SSM) and the incidence or progression of hip osteoarthritis.

DESIGN

We conducted a systematic search of five electronic databases, based on a registered protocol (available: PROSPERO CRD42020145411). Articles presenting original data on the longitudinal relationship between radiographic hip shape (quantified by SSM) and hip OA were eligible. Quantitative meta-analysis was precluded because of the use of different SSM models across studies. We used the Newcastle-Ottawa Scale (NOS) for risk of bias assessment.

RESULTS

Nine studies (6,483 hips analyzed with SSM) were included in this review. The SSM models used to describe hip shape ranged from 16 points on the femoral head to 85 points on the proximal femur and hemipelvis. Multiple hip shape features and combinations thereof were associated with incident or progressive hip OA. Shape variants that seemed to be consistently associated with hip OA across studies were acetabular dysplasia, cam morphology, and deviations in acetabular version (either excessive anteversion or retroversion).

CONCLUSIONS

Various radiographic, SSM-defined hip shape features are associated with hip OA. Some hip shape features only seem to increase the risk for hip OA when combined together. The heterogeneity of the used SSM models across studies precludes the estimation of pooled effect sizes. Further studies using the same SSM model and definition of hip OA are needed to allow for the comparison of outcomes across studies, and to validate the found associations.

Statistical shape modelling provides a responsive measure of morphological change in knee osteoarthritis over 12 months

OBJECTIVES

Responsive biomarkers are needed to assess the progression of OA and their lack has hampered previous clinical trials. Statistical shape modelling (SSM) from radiographic images identifies those at greatest risk of fast-progression or joint replacement, but its sensitivity to change has not previously been measured. This study evaluates the responsiveness of SSM in knee OA in a 12-month observational study.

METHODS

A total of 109 people were recruited who had undergone knee radiographs in the previous 12 months, and were grouped based on severity of radiographic OA (Kellgren-Lawrence grading). An SSM was built from three dual-energy X-ray absorptiometry scans at 6-month intervals. Change-over-time and OA were assessed using generalized estimating equations, standardized response means (SRM) and reliable change indices.

RESULTS

Mode 1 showed typical features of radiographic OA and had a strong link with Kellgren-Lawrence grading but did not change significantly during the study. Mode 3 showed asymmetrical changes consistent with medial cartilage loss, osteophytes and joint malalignment, and was responsive to change, with a 12-month SRM of 0.63. The greatest change was observed in the moderate radiographic OA group (SRM 0.92) compared with the controls (SRM 0.21), and the reliable change index identified 14% of this group whose progression was clinically significant.

CONCLUSION

Shape changes linked the progression of osteophytosis with increasing malalignment within the joint. Modelling of the whole joint enabled quantification of change beyond the point where bone-to-bone contact has been made. The knee SSM is, therefore, a responsive biomarker for radiographic change in knees over 12 months.

Quantitative 3D imaging parameters improve prediction of hip osteoarthritis outcome

Osteoarthritis is an increasingly important health problem for which the main treatment remains joint replacement. Therapy developments have been hampered by a lack of biomarkers that can reliably predict disease, while 2D radiographs interpreted by human observers are still the gold standard for clinical trial imaging assessment. We propose a 3D approach using computed tomography—a fast, readily available clinical technique—that can be applied in the assessment of osteoarthritis using a new quantitative 3D analysis technique called joint space mapping (JSM). We demonstrate the application of JSM at the hip in 263 healthy older adults from the AGES-Reykjavík cohort, examining relationships between 3D joint space width, 3D joint shape, and future joint replacement. Using JSM, statistical shape modelling, and statistical parametric mapping, we show an 18% improvement in prediction of joint replacement using 3D metrics combined with radiographic Kellgren & Lawrence grade (AUC 0.86) over the existing 2D FDA-approved gold standard of minimum 2D joint space width (AUC 0.73). We also show that assessment of joint asymmetry can reveal significant differences between individuals destined for joint replacement versus controls at regions of the joint that are not captured by radiographs. This technique is immediately implementable with standard imaging technologies.

An automated workflow based on hip shape improves personalized risk prediction for hip osteoarthritis in the CHECK study

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.

Prognostic factors for progression of osteoarthritis of the hip: a systematic review

Background

Predicting which patients with hip osteoarthritis are more likely to show disease progression is important for healthcare professionals. Therefore, the aim of this review was to assess which factors are predictive of progression in patients with hip osteoarthritis.

Methods

A literature search was made up until 14 March 2019. Included were cohort and case-control studies evaluating the association between factors and progression (either clinical, radiological, or THR). Excluded were studies with a follow-up < 1 year or specific underlying pathologies of osteoarthritis. Risk of bias was assessed using the QUIPS tool. A best-evidence synthesis was conducted.

Results

We included 57 articles describing 154 different factors. Of these, a best-evidence synthesis was possible for 103 factors, separately for clinical and radiological progression, and progression to total hip replacement. We found strong evidence for more clinical progression in patients with comorbidity and more progression to total hip replacement for a higher Kellgren and Lawrence grade, superior or (supero) lateral femoral head migration, and subchondral sclerosis. Strong evidence for no association was found regarding clinical progression for gender, social support, pain medication, quality of life, and limited range of motion of internal rotation or external rotation. Also, strong evidence for no association was found regarding radiological progression for the markers CTX-I, COMP, NTX-I, PINP, and PIIINP and regarding progression to total hip replacement for body mass index.

Conclusion

Strong evidence suggested that 4 factors were predictive of progression of hip osteoarthritis, whereas 12 factors were not predictive of progression. Evidence for most of the reported factors was either limited or conflicting.

Protocol registration

PROSPERO, CRD42015010757

Electronic supplementary material

The online version of this article (10.1186/s13075-019-1969-9) contains supplementary material, which is available to authorized users.

Perinatal factors associate with vertebral size and shape but not lumbar lordosis in 10-year-old children

The intrauterine environment is known to influence foetal development and future health. Low birthweight has been linked to smaller vertebral canals in children and decreased adulthood spine bone mineral content. Perinatal factors affecting lumbar spine curvature have not yet been considered but could be important for adult spinal health, as lumbar movement during lifting, a risk factor for backpain, is associated with lordosis. To investigate this, lumbar spine magnetic resonance images at age 10 years and perinatal and maternal data (birthweight, placental weight, gestation length, crown-heel length, maternal age, height, weight and smoking status) from 161 children born in Aberdeen in 1988-1989 were acquired. Statistical shape modelling, using principal component analysis, quantified variations in lumbar spine shape and resulting modes of variation were assessed in combination with perinatal data using correlations and analyses of covariance, adjusted for potential confounders. Spine modes 1-3 (SM1-SM3) captured 75% of the variation in lumbar spine shape. The first and third modes described the total amount (SM1) and evenness of curvature distribution (SM3). SM2 accounted for variations in antero-posterior vertebral diameter relative to vertebral height, increasing positive scores representing a larger relative diameter. Adjusting for gestation length and sex, SM2 positively correlated with birthweight (r = 0.25, P < 0.01), placental weight (r = 0.20, P = 0.04), crown-heel length (r = 0.36, P < 0.001) and maternal weight (r = 0.19, P = 0.04), and negatively with maternal age (r = -0.22, P = 0.02). SM2 scores were lower in girls (P < 0.001) and in the low birthweight group (P = 0.02). There were no significant differences in SM1 and SM3 scores between birthweight groups, boys and girls or children of smokers (31%) and non-smokers (69%). In conclusion, some perinatal factors were associated with vertebral body morphology but had little effect on lumbar curvature.

Age at Onset of Walking in Infancy Is Associated With Hip Shape in Early Old Age

Bones' shapes and structures adapt to the muscle and reaction forces they experience during everyday movements. Onset of independent walking, at approximately 12 months, represents the first postnatal exposure of the lower limbs to the large forces associated with bipedal movements; accordingly, earlier walking is associated with greater bone strength. However, associations between early life loading and joint shape have not been explored. We therefore examined associations between walking age and hip shape at age 60 to 64 years in 1423 individuals (740 women) from the MRC National Survey of Health and Development, a nationally representative British birth cohort. Walking age in months was obtained from maternal interview at age 2 years. Ten modes of variation in hip shape (HM1 to HM10), described by statistical shape models, were ascertained from DXA images. In sex-adjusted analyses, earlier walking age was associated with higher HM1 and HM7 scores; these associations were maintained after further adjustment for height, body composition, and socioeconomic position. Earlier walking was also associated with lower HM2 scores in women only, and lower HM4 scores in men only. Taken together, this suggests that earlier walkers have proportionately larger (HM4) and flatter (HM1, HM4) femoral heads, wider (HM1, HM4, HM7) and flatter (HM1, HM7) femoral necks, a smaller neck-shaft angle (HM1, HM4), anteversion (HM2, HM7), and early development of osteophytes (HM1). These results suggest that age at onset of walking in infancy is associated with variations in hip shape in older age. Early walkers have a larger femoral head and neck and smaller neck-shaft angle; these features are associated with reduced hip fracture risk, but also represent an osteoarthritic-like phenotype. Unlike results of previous studies of walking age and bone mass, associations in this study were not affected by adjustment for lean mass, suggesting that associations may relate directly to skeletal loading in early life when joint shape changes rapidly. © 2018 American Society for Bone and Mineral Research.

Significant morphological change in osteoarthritic hips identified over 6-12 months using statistical shape modelling

OBJECTIVE

Predicting who will develop osteoarthritis, assessing how rapidly their disease will progress and monitoring early responses to treatment are key to the development of therapeutic agents able to treat this crippling disease and to their future clinical use. Statistical Shape Modelling (SSM) enables quantification of variations in multiple geometric measures describing the whole hip joint to be considered in concert. This prospective study evaluates the responsiveness of SSM to changes in hip-shape within 1 year.

METHODS

Sixty-two people, mean age 67.1 yrs, were recruited. Dual-energy X-ray Absorptiometry images were taken at three timepoints (baseline, 6 and 12 months). Based on Kellgren-Lawrence grading (KLG) of their baseline images, subjects were classified into control/doubtful OA: KLG < 1 in both hips; moderate OA: KLG = 2; and severe OA: KLG ≥ 3 in their most severe hip. Morphology was quantified using SSM and changes in shape were assessed using generalised estimating equations. Standardized response means (SRMs) were calculated for the first and second 6 month periods, then the full 12 months.

RESULTS

Disease severity ranged from KLG0-KLG4 in the 124 hips assessed at baseline. Three SSM modes (Modes 1, 3 and 4) were associated with OA severity. Across the whole cohort, SRM magnitudes ranged from 0.16 to 0.63. The greatest subgroup SRM (magnitude 0.91) was observed over 12 months in those subjects with moderate OA (KLG2).

CONCLUSIONS

We have demonstrated that SSM can capture changes in hip shape over 6 and 12 months across the entire hip joint providing a sensitive measure of hip OA progression.

3-dimensional metrics of proximal femoral shape deformities in Legg-Calvé-Perthes disease and slipped capital femoral epiphysis

Legg-Calvé-Perthes disease (LCPD) and slipped capital femoral epiphysis (SCFE) are two common pediatric hip disorders that affect the 3-dimensional shape and function of the proximal femur. This study applied the principles of continuum mechanics to statistical shape modeling (SSM) and determined 3-D metrics for the evaluation of shape deformations in normal growth, LCPD, and SCFE. CT scans were obtained from 32 patients with asymptomatic, LCPD, and SCFE hips ((0.5-0.9 mm)² in-plane resolution, 0.63 mm slice thickness). SSM was performed on segmented proximal femoral surfaces, and shape deformations were described by surface displacement, strain, and growth plate angle metrics. Asymptomatic normal femurs underwent coordinated, growth-associated surface displacements and anisotropic strains that were site-specific and highest at the greater trochanter. After size- and age-based shape adjustment, LCPD femurs exhibited large displacements and surface strains in the femoral head and neck, with associated changes in femoral head growth plate angles. Mild SCFE femurs had contracted femoral neck surfaces, and surface displacements in all regions tended to increase with severity of slip. The results of this paper provide new 3-D metrics for characterizing the shape and biomechanics of the proximal femur. Statement of Clinical Significance: Quantitative 3-D metrics of shape may be useful for understanding and monitoring disease progression, identifying target regions for shape modulation therapies, and objectively evaluating the success of such therapies. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1526-1535, 2018.

DXA-derived hip shape is related to osteoarthritis: findings from in the MrOS cohort

OBJECTIVE

Statistical shape modelling (SSM) of radiographs has been used to explore relationships between altered joint shape and hip osteoarthritis (OA). We aimed to apply SSM to Dual-energy X-ray Absorptiometry (DXA) hip scans, and examine associations between resultant hip shape modes (HSMs), radiographic hip OA (RHOA), and hip pain, in a large population based cohort.

METHOD

SSM was performed on baseline hip DXA scans from the Osteoporotic Fractures in Men (MrOS) Study. Associations between the top ten HSMs, and prevalent RHOA from pelvic radiographs obtained 4.6 years later, were analysed in 4100 participants. RHOA was defined as Croft score ≥2. Hip pain was based on pain on walking, hip pain on examination, and Western Ontario and McMaster Universities Arthritis Index (WOMAC).

RESULTS

The five HSMs associated with RHOA showed features of either pincer- or cam-type deformities. HSM 1 (increased pincer-type deformity) was positively associated with RHOA [1.23 (1.09, 1.39)] [odds ratio (OR) and 95% CI]. HSM 8 (reduced pincer-type deformity) was inversely associated with RHOA [0.79 (0.70, 0.89)]. HSM 10 (increased cam-type deformity) was positively associated with RHOA [1.21 (1.07, 1.37)]. HSM 3 and HSM 4 (reduced cam-type deformity) were inversely associated with RHOA [0.73 (0.65, 0.83) and 0.82 (0.73, 0.93), respectively]. HSM 3 was inversely related to pain on examination [0.84 (0.76, 0.92)] and walking [0.88, (0.81, 0.95)], and to WOMAC score [0.87 (0.80, 0.93)].

CONCLUSIONS

DXA-derived measures of hip shape are associated with RHOA, and to a lesser extent hip pain, possibly reflecting their role in the pathogenesis of hip OA.

Associations between body mass index across adult life and hip shapes at age 60 to 64: Evidence from the 1946 British birth cohort

OBJECTIVE

To examine the associations of body mass index (BMI) across adulthood with hip shapes at age 60-64years.

METHODS

Up to 1633 men and women from the MRC National Survey of Health and Development with repeat measures of BMI across adulthood and posterior-anterior dual-energy X-ray absorptiometry bone mineral density images of the proximal femur recorded at age 60-64 were included in analyses. Statistical shape modelling was applied to quantify independent variations in hip mode (HM), of which the first 6 were examined in relation to: i) BMI at each age of assessment; ii) BMI gain during different phases of adulthood; iii) age first overweight.

RESULTS

Higher BMI at all ages (i.e. 15 to 60-64) and greater gains in BMI were associated with higher HM2 scores in both sexes (with positive HM2 values representing a shorter femoral neck and a wider and flatter femoral head). Similarly, younger age first overweight was associated with higher HM2 scores but only in men once current BMI was accounted for. In men, higher BMI at all ages was also associated with lower HM4 scores (with negative HM4 values representing a flatter femoral head, a wider neck and smaller neck shaft angle) but no associations with BMI gain or prolonged exposure to high BMI were found. Less consistent evidence of associations was found between BMI and the other four HMs.

CONCLUSION

These results suggest that BMI across adulthood may be associated with specific variations in hip shapes in early old age.

Statistical shape modelling of hip and lumbar spine morphology and their relationship in the MRC National Survey of Health and Development

The anatomical shape of bones and joints is important for their proper function but quantifying this, and detecting pathological variations, is difficult to do. Numerical descriptions would also enable correlations between joint shapes to be explored. Statistical shape modelling (SSM) is a method of image analysis employing pattern recognition statistics to describe and quantify such shapes from images; it uses principal components analysis to generate modes of variation describing each image in terms of a set of numerical scores after removing global size variation. We used SSM to quantify the shapes of the hip and the lumbar spine in dual-energy x-ray absorptiometry (DXA) images from 1511 individuals in the MRC National Survey of Health and Development at ages 60-64 years. We compared shapes of both joints in men and women and hypothesised that hip and spine shape would be strongly correlated. We also investigated associations with height, weight, body mass index (BMI) and local (hip or lumber spine) bone mineral density. In the hip, all except one of the first 10 modes differed between men and women. Men had a wider femoral neck, smaller neck-shaft angle, increased presence of osteophytes and a loss of the femoral head/neck curvature compared with women. Women presented with a flattening of the femoral head and greater acetabular coverage of the femoral head. Greater weight was associated with a shorter, wider femoral neck and larger greater and lesser trochanters. Taller height was accompanied by a flattening of the curve between superior head and neck and a larger lesser trochanter. Four of the first eight modes describing lumbar spine shape differed between men and women. Women tended to have a more lordotic spine than men with relatively smaller but caudally increasing anterior-posterior (a-p) vertebral diameters. Men were more likely to have a straighter spine with larger vertebral a-p diameters relative to vertebral height than women, increasing cranially. A weak correlation was found between body weight and a-p vertebral diameter. No correlations were found between shape modes and height in men, whereas in women there was a weak positive correlation between height and evenness of spinal curvature. Linear relationships between hip and spine shapes were weak and inconsistent in both sexes, thereby offering little support for our hypothesis. In conclusion, men and women entering their seventh decade have small but statistically significant differences in the shapes of their hips and their spines. Associations with height, weight, BMI and BMD are small and correspond to subtle variations whose anatomical significance is not yet clear. Correlations between hip and spine shapes are small.

Cross-sectional associations between variations in ankle shape by statistical shape modeling, injury history, and race: the Johnston County Osteoarthritis Project

Background

Injury is an important risk factor for osteoarthritis (OA), a highly prevalent and disabling joint disease. Joint shape is linked to OA, but the interplay of injury and joint shape and their combined role in OA, particularly at the ankle, is not well known. Therefore, we explored cross-sectional associations between ankle shape and injury in a large community-based cohort.

Methods

Ankles without radiographic OA were selected from the current data collection of the Johnston County OA Project. Ankles with self-reported prior injury were included as injury cases (n = 108) along with 1:1 randomly selected non-injured ankles. To define ankle shape, a 68 point model on weight-bearing lateral ankle radiographs was entered into a statistical shape model, producing a mean shape and a set of continuous variables (modes) representing variation in that shape. Nineteen modes, explaining 80% of shape variance, were simultaneously included in a logistic regression model with injury status as the dependent variable, adjusted for intra-person correlation, sex, race, body mass index (BMI), baseline OA radiographic grade, and baseline symptoms.

Results

A total of 194 participants (213 ankles) were included; mean age 71 years, BMI 30 kg/m², 67% white and 71% women. Injured ankles were more often symptomatic and from whites. In a model adjusted only for intra-person correlation, associations were seen between injury status and modes 1, 6, 13, and 19. In a fully adjusted model, race strongly affected the estimate for mode 1 (which was no longer statistically significant).

Conclusions

This study showed variations in ankle shape and history of injury as well as with race. These novel findings may indicate a change in ankle morphology following injury, or that ankle morphology predisposes to injury, and suggest that ankle shape is a potentially important factor in the development of ankle OA.

Automated assessment of bone changes in cross-sectional micro-CT studies of murine experimental osteoarthritis

OBJECTIVE

The degradation of articular cartilage, which characterises osteoarthritis (OA), is usually paired with excessive bone remodelling, including subchondral bone sclerosis, cysts, and osteophyte formation. Experimental models of OA are widely used to investigate pathogenesis, yet few validated methodologies for assessing periarticular bone morphology exist and quantitative measurements are limited by manual segmentation of micro-CT scans. The aim of this work was to chart the temporal changes in periarticular bone in murine OA by novel, automated micro-CT methods.

METHODS

OA was induced by destabilisation of the medial meniscus (DMM) in 10-week old male mice and disease assessed cross-sectionally from 1- to 20-weeks post-surgery. A novel approach was developed to automatically segment subchondral bone compartments into plate and trabecular bone in micro-CT scans of tibial epiphyses. Osteophyte volume, as assessed by shape differences using 3D image registration, and by measuring total epiphyseal volume was performed.

RESULTS

Significant linear and volumetric structural modifications in subchondral bone compartments and osteophytes were measured from 4-weeks post-surgery and showed progressive changes at all time points; by 20 weeks, medial subchondral bone plate thickness increased by 160±19.5 μm and the medial osteophyte grew by 0.124±0.028 μm3. Excellent agreement was found when automated measurements were compared with manual assessments.

CONCLUSION

Our automated methods for assessing bone changes in murine periarticular bone are rapid, quantitative, and highly accurate, and promise to be a useful tool in future preclinical studies of OA progression and treatment. The current approaches were developed specifically for cross-sectional micro-CT studies but could be applied to longitudinal studies.

Three dimensional measurement of minimum joint space width in the knee from stereo radiographs using statistical shape models

Objectives

An important measure for the diagnosis and monitoring of knee osteoarthritis is the minimum joint space width (mJSW). This requires accurate alignment of the x-ray beam with the tibial plateau, which may not be accomplished in practice. We investigate the feasibility of a new mJSW measurement method from stereo radiographs using 3D statistical shape models (SSM) and evaluate its sensitivity to changes in the mJSW and its robustness to variations in patient positioning and bone geometry.

Materials and Methods

A validation study was performed using five cadaver specimens. The actual mJSW was varied and images were acquired with variation in the cadaver positioning. For comparison purposes, the mJSW was also assessed from plain radiographs. To study the influence of SSM model accuracy, the 3D mJSW measurement was repeated with models from the actual bones, obtained from CT scans.

Results

The SSM-based measurement method was more robust (consistent output for a wide range of input data/consistent output under varying measurement circumstances) than the conventional 2D method, showing that the 3D reconstruction indeed reduces the influence of patient positioning. However, the SSM-based method showed comparable sensitivity to changes in the mJSW with respect to the conventional method. The CT-based measurement was more accurate than the SSM-based measurement (smallest detectable differences 0.55 mm versus 0. 82 mm, respectively).

Conclusion

The proposed measurement method is not a substitute for the conventional 2D measurement due to limitations in the SSM model accuracy. However, further improvement of the model accuracy and optimisation technique can be obtained. Combined with the promising options for applications using quantitative information on bone morphology, SSM based 3D reconstructions of natural knees are attractive for further development.

Cite this article: E. A. van IJsseldijk, E. R. Valstar, B. C. Stoel, R. G. H. H. Nelissen, N. Baka, R. van’t Klooster, B. L. Kaptein. Three dimensional measurement of minimum joint space width in the knee from stereo radiographs using statistical shape models. Bone Joint Res 2016;320–327. DOI: 10.1302/2046-3758.58.2000626.

Can radiographs of hip fractures predict subsequent hip fractures? A shape modelling analysis

INTRODUCTION

The geometrical shape of the proximal femur has previously been shown to predict primary hip fractures. Hip fractures are routinely diagnosed on plain radiographs of the pelvis, and these have both hips viewable. We have investigated if statistical shape modelling of the uninvolved hip on plain radiographs, at the time of the first hip fracture episode, could predict a subsequent 'second fracture' on that (uninvolved) side.

MATERIALS AND METHODS

60 radiographs taken at the time of the index hip fracture were blinded and separated into two arms; patients sustaining one hip fracture only (n=30), and those who went on to sustain a second fracture (n=30), over the three-year follow-up period. Two separate shape models were used for these groups and compared using t-tests or Mann-Whitney U-tests, along with Cohen's d to measure the effect size of each measure.

RESULTS

We found no statistically significant difference in the shape of the femur between the first fracture and second fracture group (p>0.05) and no results reached a "medium" effect size (Cohen's d <0.5).

CONCLUSIONS

Shape modelling is feasible and can be applied in the routine clinical setting. However, we were unable to elucidate any predictive value in this relatively small sample. A reliable radiograph-based method of identifying patients at risk of second fracture would be of value in planning prevention, service provision, and cost analysis. Further work is required and a study with more patients might exclude the type 2 error in our work.

Novel statistical methodology reveals that hip shape is associated with incident radiographic hip osteoarthritis among African American women

INTRODUCTION

Hip shape is a risk factor for the development of hip osteoarthritis (OA), and current methods to assess hip shape from radiographs are limited; therefore this study explored current and novel methods to assess hip shape.

METHODS

Data from a prior case-control study nested in the Johnston County OA Project were used, including 382 hips (from 342 individuals). Hips were classified by radiographic hip OA (RHOA) status as RHOA cases (baseline Kellgren Lawrence grade [KLG] 0 or 1, follow-up [mean 6 years] KLG ≥ 2) or controls (KLG = 0 or 1 at both baseline and follow-up). Proximal femur shape was assessed using a 60-point model as previously described. The current analysis explored commonly used principal component analysis (PCA), as well as novel statistical methodologies suited to high dimension low sample size settings (Distance Weighted Discrimination [DWD] and Distance Projection Permutation [DiProPerm] hypothesis testing) to assess differences between cases and controls.

RESULTS

Using these novel methodologies, we were able to better characterize morphologic differences by sex and race. In particular, the proximal femurs of African American women demonstrated significantly different shapes between cases and controls, implying an important role for sex and race in the development of RHOA. Notably, discrimination was improved with the use of DWD and DiProPerm compared to PCA.

CONCLUSIONS

DWD with DiProPerm significance testing provides improved discrimination of variation in hip morphology between groups, and enables subgroup analyses even under small sample sizes.

Variations in Hip Shape Are Associated with Radiographic Knee Osteoarthritis: Cross-sectional and Longitudinal Analyses of the Johnston County Osteoarthritis Project

OBJECTIVE

Hip shape by statistical shape modeling (SSM) is associated with hip radiographic osteoarthritis (rOA). We examined associations between hip shape and knee rOA given the biomechanical interrelationships between these joints.

METHODS

Bilateral baseline hip shape assessments [for those with at least 1 hip with a Kellgren-Lawrence arthritis grading scale (KL) 0 or 1] from the Johnston County Osteoarthritis Project were available. Proximal femur shape was defined on baseline pelvis radiographs and evaluated by SSM, producing mean shape and continuous variables representing independent modes of variation (14 modes = 95% of shape variance). Outcomes included prevalent [baseline KL ≥ 2 or total knee replacement (TKR)], incident (baseline KL 0/1 with followup ≥ 2), and progressive knee rOA (KL increase of ≥ 1 or TKR). Limb-based logistic regression models for ipsilateral and contralateral comparisons were adjusted for age, sex, race, body mass index (BMI), and hip rOA, accounting for intraperson correlations.

RESULTS

We evaluated 681 hips and 682 knees from 342 individuals (61% women, 83% white, mean age 62 yrs, BMI 29 kg/m(2)). Ninety-nine knees (15%) had prevalent rOA (4 knees with TKR). Lower modes 2 and 3 scores were associated with ipsilateral prevalent knee rOA, and only lower mode 3 scores were associated with contralateral prevalent knee rOA. No statistically significant associations were seen for incident or progressive knee rOA.

CONCLUSION

Variations in hip shape were associated with prevalent, but not incident or progressive, knee rOA in this cohort, and may reflect biomechanical differences between limbs, genetic influences, or common factors related to both hip shape and knee rOA.

Bone morphology of the femur and tibia captured by statistical shape modelling predicts rapid bone loss in acute spinal cord injury patients

After spinal cord injury (SCI), bone loss in the paralysed limbs progresses at variable rates. Decreases in bone mineral density (BMD) in the first year range from 1% (slow) to 40% (rapid). In chronic SCI, fragility fractures commonly occur around the knee, with significant associated morbidity. Osteoporosis treatments await full evaluation in SCI, but should be initiated early and targeted towards patients exhibiting rapid bone loss. The potential to predict rapid bone loss from a single bone scan within weeks of a SCI was investigated using statistical shape modelling (SSM) of bone morphology, hypothesis: baseline bone shape predicts bone loss at 12-months post-injury at fracture-prone sites. In this retrospective cohort study 25 SCI patients (median age, 33 years) were scanned at the distal femur and proximal tibia using peripheral Quantitative Computed Tomography at <5 weeks (baseline), 4, 8 and 12 months post-injury. An SSM was made for each bone. Links between the baseline shape-modes and 12-month total and trabecular BMD loss were analysed using multiple linear regression. One mode from each SSM significantly predicted bone loss (age-adjusted P<0.05 R(2)=0.37-0.61) at baseline. An elongated intercondylar femoral notch (femur mode 4, +1 SD from the mean) was associated with 8.2% additional loss of femoral trabecular BMD at 12-months. A more concave posterior tibial fossa (tibia mode 3, +1 SD) was associated with 9.4% additional 12-month tibial trabecular BMD loss. Baseline bone shape determined from a single bone scan is a valid imaging biomarker for the prediction of 12-month bone loss in SCI patients.

A systematic review of the relationship between subchondral bone features, pain and structural pathology in peripheral joint osteoarthritis

INTRODUCTION

Bone is an integral part of the osteoarthritis (OA) process. We conducted a systematic literature review in order to understand the relationship between non-conventional radiographic imaging of subchondral bone, pain, structural pathology and joint replacement in peripheral joint OA.

METHODS

A search of the Medline, EMBASE and Cochrane library databases was performed for original articles reporting association between non-conventional radiographic imaging-assessed subchondral bone pathologies and joint replacement, pain or structural progression in knee, hip, hand, ankle and foot OA. Each association was qualitatively characterised by a synthesis of the data from each analysis based upon study design, adequacy of covariate adjustment and quality scoring.

RESULTS

In total 2456 abstracts were screened and 139 papers were included (70 cross-sectional, 71 longitudinal analyses; 116 knee, 15 hip, six hand, two ankle and involved 113 MRI, eight DXA, four CT, eight scintigraphic and eight 2D shape analyses). BMLs, osteophytes and bone shape were independently associated with structural progression or joint replacement. BMLs and bone shape were independently associated with longitudinal change in pain and incident frequent knee pain respectively.

CONCLUSION

Subchondral bone features have independent associations with structural progression, pain and joint replacement in peripheral OA in the hip and hand but especially in the knee. For peripheral OA sites other than the knee, there are fewer associations and independent associations of bone pathologies with these important OA outcomes which may reflect fewer studies; for example the foot and ankle were poorly studied. Subchondral OA bone appears to be a relevant therapeutic target.

SYSTEMATIC REVIEW

PROSPERO registration number: CRD 42013005009.

Validation of statistical shape modelling to predict hip osteoarthritis in females: data from two prospective cohort studies (Cohort Hip and Cohort Knee and Chingford)

OBJECTIVES

To prospectively investigate whether hip shape variants at baseline are associated with the need for future total hip replacement (THR) in women and to validate the resulting associated shape variants of the Cohort Hip and Cohort Knee (CHECK) cohort and the Chingford cohort.

METHODS

Female participants from the CHECK cohort without radiographic OA (Kellgren-Lawrence score <2) at baseline were included (1100 hips); 22 hips had a THR within 5 years of follow-up. For the Chingford cohort, with only female participants, hips without radiographic OA at baseline were selected and a nested case-control design was used, with 19 THR cases within 19 years of follow-up and 95 controls matched 5 to 1 for age and BMI. Hip shape on baseline anteroposterior pelvic radiographs was assessed by statistical shape modelling (SSM) using the same model for both cohorts.

RESULTS

In the CHECK and Chingford cohorts, the respective mean age was 55.8 (s.d. 5.1) and 53.6 (s.d. 5.4) and the BMI was 26.14 (s.d. 4.3) and 25.7 (s.d. 3.3), respectively. Multiple shape variants of the hip were significantly (P < 0.05) associated with future THR in both the CHECK (modes 4, 11, 15, 17 and 22) and Chingford (modes 2 and 17) cohorts. Mode 17 [odds ratio (OR) 0.51 (95% CI 0.33, 0.80) in the CHECK cohort], representing a flattened head-neck junction and flat greater trochanter, could be confirmed in the Chingford cohort [OR 0.41 (95% CI 0.23, 0.82)]. Modes 4 and 15 of the CHECK cohort also showed non-significant trends in the Chingford cohort.

CONCLUSION

Several baseline shape variants are associated with the future need for THR within a cohort. Despite differences in participant characteristics, radiographic protocol and follow-up time, we could validate at least one shape variant, suggesting that SSM is reasonably transferable between cohorts.

Investigation of association between hip osteoarthritis susceptibility loci and radiographic proximal femur shape

OBJECTIVE

To test whether previously reported hip morphology or osteoarthritis (OA) susceptibility loci are associated with proximal femur shape as represented by statistical shape model (SSM) modes and as univariate or multivariate quantitative traits.

METHODS

We used pelvic radiographs and genotype data from 929 subjects with unilateral hip OA who had been recruited previously for the Arthritis Research UK Osteoarthritis Genetics Consortium genome-wide association study. We built 3 SSMs capturing the shape variation of the OA-unaffected proximal femur in the entire mixed-sex cohort and for male/female-stratified cohorts. We selected 41 candidate single-nucleotide polymorphisms (SNPs) previously reported as being associated with hip morphology (for replication analysis) or OA (for discovery analysis) and for which genotype data were available. We performed 2 types of analysis for genotype-phenotype associations between these SNPs and the modes of the SSMs: 1) a univariate analysis using individual SSM modes and 2) a multivariate analysis using combinations of SSM modes.

RESULTS

The univariate analysis identified association between rs4836732 (within the ASTN2 gene) and mode 5 of the female SSM (P = 0.0016) and between rs6976 (within the GLT8D1 gene) and mode 7 of the mixed-sex SSM (P = 0.0003). The multivariate analysis identified association between rs5009270 (near the IFRD1 gene) and a combination of modes 3, 4, and 9 of the mixed-sex SSM (P = 0.0004). Evidence of associations remained significant following adjustment for multiple testing. All 3 SNPs had previously been associated with hip OA.

CONCLUSION

These de novo findings suggest that rs4836732, rs6976, and rs5009270 may contribute to hip OA susceptibility by altering proximal femur shape.

Proximal femur shape differs between subjects with lateral and medial knee osteoarthritis and controls: the Osteoarthritis Initiative

OBJECTIVE

This study examined the association of proximal femur shape with ipsilateral medial and lateral compartment knee osteoarthritis (OA).

DESIGN

Radiographs were obtained from the NIH-funded Osteoarthritis Initiative (OAI). Cases of isolated radiographic lateral compartment knee OA were defined on baseline radiographs as Kellgren/Lawrence (K/L) Grade ≥ 2 and joint space narrowing (JSN) > 0 in the lateral compartment and JSN = 0 in the medial compartment; isolated medial compartment knee OA had K/L ≥ 2 and JSN > 0 medially with JSN = 0 in the lateral compartment. Controls had K/L < 2 and JSN = 0 in both compartments. Controls were frequency matched to cases by sex and 10-year age intervals. We characterized the shape of the proximal femurs on radiographs using Active Shape Modeling (ASM) and determined the association of proximal femur shape with knee OA using logistic regression.

RESULTS

There were 168 lateral compartment knee OA cases (mean body mass index (BMI) 29.72 ± 5.26), 169 medial compartment knee OA cases (mean BMI 29.68 ± 4.83) and 168 controls (mean BMI 26.87 ± 4.2). Thirteen modes were derived for femur shape which described 95.5% of the total variance in proximal femur shape in the population. Modes 6, 8 and 12 were associated with prevalent lateral compartment knee OA. Medial compartment knee OA was associated with proximal femur modes 1, 5, 8, and 12.

CONCLUSIONS

Prevalent lateral and medial compartment knee OA are associated with different ipsilateral proximal femur shapes. Additional studies are needed to better define how the shape of the proximal femur influences compartment-specific knee OA.

Increasing shape modelling accuracy by adjusting for subject positioning: an application to the analysis of radiographic proximal femur symmetry using data from the Osteoarthritis Initiative

In total hip arthroplasty, the shape of the contra-lateral femur frequently serves as a template for preoperative planning. Previous research on contra-lateral femoral symmetry has been based on conventional hip geometric measurements (which reduce shape to a series of linear measurements) and did not take the effect of subject positioning on radiographic femur shape into account. The aim of this study was to analyse proximal femur symmetry based on statistical shape models (SSMs) which quantify global femoral shape while also adjusting for differences in subject positioning during image acquisition. We applied our recently developed fully automatic shape model matching (FASMM) system to automatically segment the proximal femur from AP pelvic radiographs to generate SSMs of the proximal femurs of 1258 Caucasian females (mean age: 61.3 SD=9.0). We used a combined SSM (capturing the left and right femurs) to identify and adjust for shape variation attributable to subject positioning as well as a single SSM (including all femurs as left femurs) to analyse proximal femur symmetry. We also calculated conventional hip geometric measurements (head diameter, neck width, shaft width and neck-shaft angle) using the output of the FASMM system. The combined SSM revealed two modes that were clearly attributable to subject positioning. The average difference (mean point-to-curve distance) between left and right femur shape was 1.0mm before and 0.8mm after adjusting for these two modes. The automatic calculation of conventional hip geometric measurements after adjustment gave an average absolute percent asymmetry of within 3.1% and an average absolute difference of within 1.1mm or 2.9° for all measurements. We conclude that (i) for Caucasian females the global shape of the right and left proximal femurs is symmetric without isolated locations of asymmetry; (ii) a combined left-right SSM can be used to adjust for radiographic shape variation due to subject positioning; and (iii) adjusting for subject positioning increases the accuracy of predicting the shape of the contra-lateral hip.

Statistical shape and appearance models of bones

When applied to bones, statistical shape models (SSM) and statistical appearance models (SAM) respectively describe the mean shape and mean density distribution of bones within a certain population as well as the main modes of variations of shape and density distribution from their mean values. The availability of this quantitative information regarding the detailed anatomy of bones provides new opportunities for diagnosis, evaluation, and treatment of skeletal diseases. The potential of SSM and SAM has been recently recognized within the bone research community. For example, these models have been applied for studying the effects of bone shape on the etiology of osteoarthritis, improving the accuracy of clinical osteoporotic fracture prediction techniques, design of orthopedic implants, and surgery planning. This paper reviews the main concepts, methods, and applications of SSM and SAM as applied to bone.

Osteochondral tissue engineering with biphasic scaffold: current strategies and techniques

The management of osteoarthritis (OA) remains challenging and controversial. Although several clinical options exist for the treatment of OA, regeneration of the damaged articular cartilage has proved difficult due to the limited healing capacity. With the advancements in tissue engineering and cell-based technologies over the past decade, new therapeutic options for patients with osteochondral lesions potentially exist. This review will focus on the feasibility of tissue-engineered biphasic scaffolds, which can mimic the native osteochondral complex, for osteochondral repair and highlight the recent development of these techniques toward tissue regeneration. Moreover, basic anatomy, strategy for osteochondral repair, the design and fabrication methods of scaffolds, as well as the choice of cells, growth factor, and materials will be discussed. Specifically, we focus on the latest preclinical animal studies using large animals and clinical trials with high clinical relevance. In turn, this will facilitate an understanding of the latest trends in osteochondral repair and contribute to the future application of such clinical therapies in patients with OA.

Quantifying the characteristics of the acetabulum and proximal femur using a semi-automated hip morphology software programme (SHIPS)

BACKGROUND

Hip-joint shape is an important factor that affects an individual's risk of developing osteoarthritis (OA). However, validated tools to quantify these morphological characteristics on clinical plain radiographs are few. We developed a Matlab-based Semi-automated HIP Software programme (SHIPS) that measures 10 morphologic characteristics of the hip that are risk factors for OA using a plain digitised antero-posterior pelvic radiograph. In this study we validated the accuracy and repeatability of this freeware-tool.

METHODS

Software accuracy was assessed using a test pelvic radiograph, and by repeated measurements of an AP-pelvic radiograph digitally recreated from pelvis computed-tomography images reformatted to create images rotated in 2-dimensions (2.5° increments, range -15° to +15°). Intra- and inter-observer repeatability was assessed using pelvic radiographs from 30 subjects analysed twice using the software by two readers, and expressed as coefficient of variation (CV). Clinical-repeatability was assessed by measuring sequential pelvic radiographs taken on the same day after re-positioning in 23 subjects.

RESULTS

Software accuracy was within 0.1% for linear-ratios and 0.4° for angular-measurements. Changes in pelvic inclination and rotation of ±15° resulted in <14% change in linear-measurement ratios and <7° change in angular-measurements. The intra-observer CV was between 0.3 to 4.1%, and inter-
observer CV 0.7 to 9.7% with the exception of horizontal-toit-externa (HTE, 14.6 and 24.2% respectively). Short-term clinical-repeatability varied from 0.4 to 6.1%, with the exception of HTE (37.4%).

CONCLUSION

The software showed good accuracy and repeatability for measurement of most hip-joint morphologic risk factors for OA apart from HTE. This tool has particular value in studying large or 
retrospective datasets where cross-sectional imaging is not feasible or available.

Osteoarthritis year 2013 in review: clinical

Some major themes over the past year in clinical research of osteoarthritis (OA) include obesity, muscle strength, pain mechanisms, novel disease modifying drugs, and risk factors for poor outcomes of joint replacement surgery. A systematic literature search was performed using PubMed from January 2012 to December 2012. The articles selected for this review represent topics the authors thought best highlight recent clinical OA research.

Development of a fully automatic shape model matching (FASMM) system to derive statistical shape models from radiographs: application to the accurate capture and global representation of proximal femur shape

OBJECTIVE

To evaluate the accuracy and sensitivity of a fully automatic shape model matching (FASMM) system to derive statistical shape models (SSMs) of the proximal femur from non-standardised anteroposterior (AP) pelvic radiographs.

DESIGN

AP pelvic radiographs obtained with informed consent and appropriate ethical approval were available for 1105 subjects with unilateral hip osteoarthritis (OA) who had been recruited previously for The arcOGEN Study. The FASMM system was applied to capture the shape of the unaffected (i.e., without signs of radiographic OA) proximal femur from these radiographs. The accuracy and sensitivity of the FASMM system in calculating geometric measurements of the proximal femur and in shape representation were evaluated relative to validated manual methods.

RESULTS

De novo application of the FASMM system had a mean point-to-curve error of less than 0.9 mm in 99% of images (n = 266). Geometric measurements generated by the FASMM system were as accurate as those obtained manually. The analysis of the SSMs generated by the FASMM system for male and female subject groups identified more significant differences (in five of 17 SSM modes after Bonferroni adjustment) in their global proximal femur shape than those obtained from the analysis of conventional geometric measurements. Multivariate gender-classification accuracy was higher when using SSM mode values (76.3%) than when using conventional hip geometric measurements (71.8%).

CONCLUSIONS

The FASMM system rapidly and accurately generates a global SSM of the proximal femur from radiographs of varying quality and resolution. This system will facilitate complex morphometric analysis of global shape variation across large datasets. The FASMM system could be adapted to generate SSMs from the radiographs of other skeletal structures such as the hand, knee or pelvis.

Total hip replacement but not clinical osteoarthritis can be predicted by the shape of the hip: a prospective cohort study (CHECK)

OBJECTIVE

To investigate the association between baseline hip shape and both clinical hip osteoarthritis (OA) and total hip replacement (THR) at 5-year follow-up.

DESIGN

Individuals from the Cohort Hip and Cohort Knee (CHECK) study, with early symptomatic OA, having standardized anteroposterior pelvic radiographs at baseline and 5-year follow-up (n = 723) were included. Hip shape on the radiographs was assessed using statistical shape modeling (SSM). Hips fulfilling the American College of Rheumatology (ACR) criteria at follow-up were classified as clinical OA. The association between each mode of shape variation and both outcome measures was calculated by Generalized Estimating Equations (GEE).

RESULTS

The included individuals comprised 575 females and 148 males (mean age 55.9 ± 5.2 years). At baseline, 8% fulfilled the ACR criteria, 76% had no radiographic hip OA [Kellgren & Lawrence (K&L) = 0] and 24% had doubtful OA (K&L = 1). At follow-up, 147 hips (10.4%) fulfilled the ACR criteria and 35 hips (2.5%) had received THR. Five shape variants (modes) at baseline associated significantly with THR within 5 years. When combined in one GEE model, these shape variants resulted in a predictive power indicated by an area under the curve of 0.81. No shape variants associated with the presence of clinical OA at follow-up.

CONCLUSION

The shape of the hip as quantified by an SSM has a good predictive value for THR, whereas variation in shape cannot predict clinical OA. Minor shape variants may be used as a radiographic biomarker to predict the future risk of THR.

Osteoarthritis year 2012 in review: imaging

PURPOSE

This article reviews original publications related to imaging in osteoarthritis (OA) published in English from September 2011 through March 2012. In vitro data and animal studies are not covered.

METHODS

To extract relevant studies, an extensive PubMed database search was performed using the query terms "osteoarthritis" in combination with "MRI", "imaging", "radiography", "ultrasound", "computed tomography" and "nuclear medicine". Publications were sorted according to relevance based on potential impact to the OA research community with the over all goal of a balanced overview of all aspects of imaging. Focus was on publications in high-impact special-interest journals. The literature will be presented by topics covering radiography, morphologic magnetic resonance imaging (MRI), compositional and high-field MRI, quantitative MRI, ultrasound, other joints and systematic reviews. Original research that was presented as a podium or poster presentation at osteoarthritis research society international (OARSI) 2012 will not be included.

RESULTS AND CONCLUSIONS

For the search topics "MRI" and "osteoarthritis" a decrease in overall publications was observed over the 6 months following September 2011 when compared to the previous 6 months (-38.1%). For the terms "radiography" and "osteoarthritis" a decrease of 56.9% was noted. The 6 months since the last OARSI conference were characterized by several MRI-based studies dealing with epidemiologic and methodologic aspects of disease. Other modalities such as radiography or ultrasound received much less attention. Most imaging research is still concentrated on the knee although interest in other sites, especially the hand, has increased since the last OARSI meeting.