Dual-energy X-ray absorptiometry derived knee shape may provide a useful imaging biomarker for predicting total knee replacement: Findings from a study of 37,843 people in UK Biobank

Objective

We aimed to create an imaging biomarker for knee shape using knee dual-energy x-ray absorptiometry (DXA) scans and investigate its potential association with subsequent total knee replacement (TKR), independently of radiographic features of knee osteoarthritis and established risk factors.

Methods

Using a 129-point statistical shape model, knee shape (expressed as a B-score) and minimum joint space width (mJSW) of the medial joint compartment (binarized as above or below the first quartile) were derived. Osteophytes were manually graded in a subset of images and an overall score was assigned. Cox proportional hazards models were used to examine the associations of B-score, mJSW and osteophyte score with TKR risk, adjusting for age, sex, height and weight.

Results

The analysis included 37,843 individuals (mean age 63.7 years). In adjusted models, B-score was associated with TKR: each unit increase in B-score, reflecting one standard deviation from the mean healthy shape, corresponded to a hazard ratio (HR) of 2.25 (2.08, 2.43), while a lower mJSW had a HR of 2.28 (1.88, 2.77). Among the 6719 images scored for osteophytes, mJSW was replaced by osteophyte score in the most strongly predictive model for TKR. In ROC analyses, a model combining B-score, osteophyte score, and demographics outperformed a model including demographics alone (AUC ​= ​0.87 vs 0.73).

Conclusions

Using statistical shape modelling, we derived a DXA-based imaging biomarker for knee shape that was associated with kOA progression. When combined with osteophytes and demographic data, this biomarker may help identify individuals at high risk of TKR, facilitating targeted interventions.

Femoral neck width genetic risk score is a novel independent risk factor for hip fractures

Femoral neck width (FNW) derived from DXA scans may provide a useful adjunct to hip fracture prediction. Therefore, we investigated whether FNW is related to hip fracture risk independently of femoral neck bone mineral density (FN-BMD), using a genetic approach. FNW was derived from points automatically placed on the proximal femur using hip DXA scans from 38 150 individuals (mean age 63.8 yr, 48.0% males) in UK Biobank (UKB). Genome-wide association study (GWAS) identified 71 independent genome-wide significant FNW SNPs, comprising genes involved in cartilage differentiation, hedgehog, skeletal development, in contrast to SNPs identified by FN-BMD GWAS which primarily comprised runx1/Wnt signaling genes (MAGMA gene set analyses). FNW and FN-BMD SNPs were used to generate genetic instruments for multivariable Mendelian randomization. Greater genetically determined FNW increased risk of all hip fractures (odds ratio [OR] 1.53; 95% CI, 1.29-1.82 per SD increase) and femoral neck fractures (OR 1.58;1.30-1.92), but not trochanteric or forearm fractures. In contrast, greater genetically determined FN-BMD decreased fracture risk at all 4 sites. FNW and FN-BMD SNPs were also used to generate genetic risk scores (GRSs), which were examined in relation to incident hip fracture in UKB (excluding the FNW GWAS population; n = 338 742, 3222 cases) using a Cox proportional hazards model. FNW GRS was associated with increased risk of all incident hip fractures (HR 1.08;1.05-1.12) and femoral neck fractures (hazard ratio [HR] 1.10;1.06-1.15), but not trochanteric fractures, whereas FN-BMD GRS was associated with reduced risk of all hip fracture types. We conclude that the underlying biology regulating FNW and FN-BMD differs, and that DXA-derived FNW is causally related to hip fractures independently of FN-BMD, adding information beyond FN-BMD for hip fracture prediction. Hence, FNW derived from DXA analyses or a FNW GRS may contribute clinically useful information beyond FN-BMD for hip fracture prediction.

Associations between life course longitudinal growth and hip shapes at ages 60-64 years: evidence from the MRC National Survey of Health and Development

OBJECTIVE

We sought to examine associations between height gain across childhood and adolescence with hip shape in individuals aged 60-64 years from the Medical Research Council National Survey of Health and Development, a nationally representative British birth cohort.

METHODS

Height was measured at ages 2, 4, 6, 7, 11 and 15 years, and self-reported at age 20 years. 10 modes of variation in hip shape (HM1-10), described by statistical shape models, were previously ascertained from DXA images taken at ages 60-64 years. Associations between (1) height at each age; (2) Super-Imposition by Translation And Rotation (SITAR) growth curve variables of height size, tempo and velocity; and (3) height gain during specific periods of childhood and adolescence, and HM1-10 were tested.

RESULTS

Faster growth velocity was associated with a wider, flatter femoral head and neck, as described by positive scores for HM6 (regression coefficient 0.014; 95% CI 0.08 to 0.019; p<0.001) and HM7 (regression coefficient 0.07; 95% CI 0.002 to 0.013; p=0.009), and negative scores for HM10 (regression coefficient -0.006; 95% CI -0.011 to 0.00, p=0.04) and HM2 (males only, regression coefficient -0.017; 95% CI -0.026 to -0.09; p<0.001). Similar associations were observed with greater height size and later height tempo. Examination of height gains during specific periods of childhood and adolescence identified those during the adolescence period as being most consistently associated.

CONCLUSION

Our analyses suggest that individual growth patterns, particularly in the adolescent period, are associated with modest variations in hip shape at 60-64 years, which are consistent with features seen in osteoarthritis.

Hip geometric parameters are associated with radiographic and clinical hip osteoarthritis: Findings from a cross-sectional study in UK Biobank

OBJECTIVE

To examine the extent to which geometric parameters derived from dual-energy x-ray absorptiometry (DXA) scans in the UK Biobank study are related to hip osteoarthritis (HOA) independently of sex, age and body size.

DESIGN

Femoral neck width (FNW), diameter of the femoral head (DFH) and hip axis length (HAL) were derived automatically from left hip DXA scans in UK Biobank using outline points placed around the hip by a machine-learning program. Correlations were calculated between geometric parameters, age, height, and weight. Logistic regression was used to examine the relationship of geometric parameters with radiographic HOA, hospital diagnosed HOA (HESOA), and Cox proportional hazards model to evaluate the relationship with total hip replacement (THR). Analyses were adjusted for sex, age, height, weight, and geometric parameters.

RESULTS

The study consisted of 40,312 participants. In age and sex-adjusted analyses, FNW, HAL and DFH were related to increased risk of radiographic HOA. In a model adjusted for age, sex, height, weight and other geometric parameters, both FNW and HAL retained independent relationships with radiographic HOA [FNW: odds ratios 2.38 (2.18-2.59), HAL: 1.25 (1.15-1.36)], while DFH was now protective [0.55 (0.50-0.61)]. Only FNW was independently related to HESOA [2.20 (1.80-2.68)] and THR [hazard ratios 2.51 (1.89-3.32)].

CONCLUSION

Greater FNW and HAL were independently related to an increased risk of radiographic HOA, whereas greater DFH appeared to be protective. Greater FNW was independently associated with HESOA and THR. These results suggest that DXA-derived geometric parameters, particularly FNW, could help determine HOA and THR risk.

Comparison between UK Biobank and Shanghai Changfeng suggests distinct hip morphology may contribute to ethnic differences in the prevalence of hip osteoarthritis

OBJECTS

Joint morphology is a risk factor for hip osteoarthritis (HOA) and could explain ethnic differences in HOA prevalence. Therefore, we aimed to compare the prevalence of radiographic HOA (rHOA) and hip morphology between the predominantly White UK Biobank (UKB) and exclusively Chinese Shanghai Changfeng (SC) cohorts.

METHODS

Left hip iDXA scans were used to quantify rHOA, from a combination of osteophytes (grade ≥1) and joint space narrowing (grade ≥1), and hip morphology. Using an 85-point Statistical Shape Model (SSM) we evaluated cam (alpha angle ≥60°) and pincer (lateral centre-edge angle (LCEA) ≥45°) morphology and acetabular dysplasia (LCEA <25°). Diameter of femoral head (DFH), femoral neck width (FNW), and hip axis length (HAL) were also obtained from these points. Results were adjusted for differences in age, height, and weight and stratified by sex.

RESULTS

Complete data were available for 5924 SC and 39,020 White UKB participants with mean ages of 63.4 and 63.7 years old. rHOA prevalence was considerably lower in female (2.2% versus 13.1%) and male (12.0% and 25.1%) SC compared to UKB participants. Cam morphology, rarely seen in females, was less common in SC compared with UKB males (6.3% versus 16.5%). Composite SSM modes, scaled to the same overall size, revealed SC participants to have a wider femoral head compared to UKB participants. FNW and HAL were smaller in SC compared to UKB, whereas DFH/FNW ratio was higher in SC.

CONCLUSIONS

rHOA prevalence is lower in Chinese compared with White individuals. Several differences in hip shape were observed, including frequency of cam morphology, FNW, and DFH/FNW ratio. These characteristics have previously been identified as risk factors for HOA and may contribute to observed ethnic differences in HOA prevalence.

The identification of distinct protective and susceptibility mechanisms for hip osteoarthritis: findings from a genome-wide association study meta-analysis of minimum joint space width and Mendelian randomisation cluster analyses

BACKGROUND

Hip minimum joint space width (mJSW) provides a proxy for cartilage thickness. This study aimed to conduct a genome-wide association study (GWAS) of mJSW to (i) identify new genetic determinants of mJSW and (ii) identify which mJSW loci convey hip osteoarthritis (HOA) risk and would therefore be of therapeutic interest.

METHODS

GWAS meta-analysis of hip mJSW derived from plain X-rays and DXA was performed, stratified by sex and adjusted for age and ancestry principal components. Mendelian randomisation (MR) and cluster analyses were used to examine causal effect of mJSW on HOA.

FINDINGS

50,745 individuals were included in the meta-analysis. 42 SNPs, which mapped to 39 loci, were identified. Mendelian randomisation (MR) revealed little evidence of a causal effect of mJSW on HOA (ORIVW 0.98 [95% CI 0.82-1.18]). However, MR-Clust analysis suggested the null MR estimates reflected the net effect of two distinct causal mechanisms cancelling each other out, one of which was protective, whereas the other increased HOA susceptibility. For the latter mechanism, all loci were positively associated with height, suggesting mechanisms leading to greater height and mJSW increase the risk of HOA in later life.

INTERPRETATIONS

One group of mJSW loci reduce HOA risk via increased mJSW, suggesting possible utility as targets for chondroprotective therapies. The second group of mJSW loci increased HOA risk, despite increasing mJSW, but were also positively related to height, suggesting they contribute to HOA risk via a growth-related mechanism.

FUNDING

Primarily funded by the Medical Research Council and Wellcome Trust.

A Genome-Wide Association Study Meta-Analysis of Alpha Angle Suggests Cam-Type Morphology May Be a Specific Feature of Hip Osteoarthritis in Older Adults

OBJECTIVE

To examine the genetic architecture of cam morphology using alpha angle (AA) as a proxy measure and conduct an AA genome-wide association study (GWAS) followed by Mendelian randomization (MR) to evaluate its causal relationship with hip osteoarthritis (OA).

METHODS

Observational analyses examined associations between AA measurements derived from hip dual x-ray absorptiometry (DXA) scans from the UK Biobank study and radiographic hip OA outcomes and subsequent total hip replacement. Following these analyses, an AA GWAS meta-analysis was performed (N = 44,214) using AA measurements previously derived in the Rotterdam Study. Linkage disequilibrium score regression assessed the genetic correlation between AA and hip OA. Genetic associations considered significant (P < 5 × 10-8 ) were used as AA genetic instrument for 2-sample MR analysis.

RESULTS

DXA-derived AA showed expected associations between AA and radiographic hip OA (adjusted odds ratio [OR] 1.63 [95% confidence interval (95% CI) 1.58, 1.67]) and between AA and total hip replacement (adjusted hazard ratio 1.45 [95% CI 1.33, 1.59]) in the UK Biobank study cohort. The heritability of AA was 10%, and AA had a moderate genetic correlation with hip OA (rg  = 0.26 [95% CI 0.10, 0.43]). Eight independent genetic signals were associated with AA. Two-sample MR provided weak evidence of causal effects of AA on hip OA risk (inverse variance weighted OR 1.84 [95% CI 1.14, 2.96], P = 0.01). In contrast, genetic predisposition for hip OA had stronger evidence of a causal effect on increased AA (inverse variance weighted β = 0.09 [95% CI 0.04, 0.13], P = 4.58 × 10-5 ).

CONCLUSION

Expected observational associations between AA and related clinical outcomes provided face validity for the DXA-derived AA measurements. Evidence of bidirectional associations between AA and hip OA, particularly for risk of hip OA on AA, suggests that hip shape modeling secondary to a genetic predisposition to hip OA contributes to the well-established relationship between hip OA and cam morphology in older adults.

Machine Learning-Derived Acetabular Dysplasia and Cam Morphology Are Features of Severe Hip Osteoarthritis: Findings From UK Biobank

The contribution of shape changes to hip osteoarthritis (HOA) remains unclear, as is the extent to which these vary according to HOA severity. In the present study, we used statistical shape modeling (SSM) to evaluate relationships between hip shape and HOA of different severities using UK Biobank DXA images. We performed a cross-sectional study in individuals with left hip dual-energy X-ray absorptiometry (DXA) scans. Statistical shape modeling (SSM) was used to quantify hip shape. Radiographic HOA (rHOA) was classified using osteophyte size and number and joint space narrowing. HOA outcomes ranged in severity from moderate (grade 2) to severe (grade ≥3) rHOA, hospital-diagnosed HOA, and subsequent total hip replacement (THR). Confounder-adjusted logistic regression between the top 10 hip shape modes (HSMs) and OA outcomes was performed. Further models adjusted for alpha angle (AA) and lateral center-edge angle (LCEA), reflecting acetabular dysplasia and cam morphology, respectively. Composite HSM figures were produced combining HSMs associated with separate OA outcomes. A total of 40,311 individuals were included (mean 63.7 years, 47.8% male), of whom 5.7% had grade 2 rHOA, 1.7% grade ≥3 rHOA, 1.3% hospital-diagnosed HOA, and 0.6% underwent THR. Composite HSM figures for grade 2 rHOA revealed femoral neck widening, increased acetabular coverage, and enlarged lesser and greater trochanters. In contrast, grade ≥3 rHOA, hospital-diagnosed HOA, and THR were suggestive of cam morphology and reduced acetabular coverage. Associations between HSMs depicting cam morphology and reduced acetabular coverage and more severe HOA were attenuated by AA and LCEA adjustment, respectively. Relationships between hip shape and HOA differed according to severity. Notably, cam morphology and acetabular dysplasia were features of severe HOA, but unrelated to moderate disease, suggesting possible prognostic utility. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

A novel semi-automated classifier of hip osteoarthritis on DXA images shows expected relationships with clinical outcomes in UK Biobank

OBJECTIVE

Conventional scoring methods for radiographic hip OA (rHOA) are subjective and show inconsistent relationships with clinical outcomes. To provide a more objective rHOA scoring method, we aimed to develop a semi-automated classifier based on DXA images and confirm its relationships with clinical outcomes.

METHODS

Hip DXAs in UK Biobank (UKB) were marked up for osteophyte area from which acetabular, superior and inferior femoral head osteophyte grades were derived. Joint space narrowing (JSN) grade was obtained automatically from minimum joint space width (mJSW) measures. Clinical outcomes related to rHOA comprised hip pain, hospital diagnosed OA (HES OA) and total hip replacement. Logistic regression and Cox proportional hazard modelling were used to examine associations between overall rHOA grade (0-4; derived from combining osteophyte and JSN grades) and the clinical outcomes.

RESULTS

A toal of 40 340 individuals were included in the study (mean age 63.7), of whom 81.2% had no evidence of rHOA, while 18.8% had grade ≥1 rHOA. Grade ≥1 osteophytes at each location and JSN were associated with hip pain, HES OA and total hip replacement. Associations with all three clinical outcomes increased progressively according to rHOA grade, with grade 4 rHOA and total hip replacement showing the strongest association [57.70 (38.08-87.44)].

CONCLUSIONS

Our novel semi-automated tool provides a useful means for classifying rHOA on hip DXAs, given its strong and progressive relationships with clinical outcomes. These findings suggest DXA scanning can be used to classify rHOA in large DXA-based cohort studies supporting further research, with the future potential for population-based screening.

Walking on water: revisiting the role of water in articular cartilage biomechanics in relation to tissue engineering and regenerative medicine

The importance, and the difficulty, of generating biosynthetic articular cartilage is widely recognized. Problems arise from obtaining sufficient stiffness, toughness and longevity in the material and integration of new material into existing cartilage and bone. Much work has been done on chondrocytes and tissue macromolecular components while water, which comprises the bulk of the tissue, is largely seen as a passive component; the ‘solid matrix’ is believed to be the main load-bearing element most of the time. Water is commonly seen as an inert filler whose restricted flow through the tissue is believed to be sufficient to generate the properties measured. We propose that this model should be turned on its head. Water comprises 70–80% of the matrix and has a bulk modulus considerably greater than that of cartilage. We suggest that the macromolecular components structure the water to support the loads applied. Here, we shall examine the structure and organization of the main macromolecules, collagen, aggrecan and hyaluronan, and explore how water interacts with their polyelectrolyte nature. This may inform the biosynthetic process by identifying starting points to enable developing tissue properties to guide the cells into producing the appropriate macromolecular composition and structure.

OA20 Radiographic hip osteoarthritis classified semi-automatically on dual-energy x-ray absorptiometry scans is strongly predictive of total hip replacement: findings from UK Biobank

Background/Aims

Radiographic hip osteoarthritis (rHOA) is traditionally defined on hip x-rays, using subjective methods such as Kellgren-Lawrence scoring. Associations between subjective rHOA measures and symptoms are inconsistent. Applying digital tools to high-resolution dual-energy X-ray absorptiometry (DXA) scans, we aimed to develop a novel semi-automated classifier for rHOA and evaluate the face validity of the classifier based on relationships with hip pain, hospital diagnosed OA (HES OA), and risk of total hip replacement (THR).

Methods

Using hip DXAs in UK Biobank, osteophyte grades 0-3 were assigned based on manually measured osteophyte area. Minimum joint space width (mJSW) was automatically measured using outline points placed by a machine learning-based algorithm and then used to categorise individuals into joint space narrowing (JSN) grades 0-3. Osteophyte and JSN grades were combined, using a novel system giving greater to weight to osteophytes, to categorise individuals into rHOA grades 0-4. Logistic regression giving odds ratios (OR) was used to examine associations between rHOA grade and hip pain, and HES OA. Cox proportional hazard models giving hazard ratios (HR) were used to examine associations between rHOA grade and subsequent THR. Our adjusted model included age, sex, height and weight as covariates.

Results

40,340 individuals were included in the study (mean age 63.7 [range 44-82], 19294/21046 male/female). 32758 (81.2%) had rHOA grade 0, 4565 (11.3%) grade 1, 2317 (5.7%) grade 2, 543 (1.3%) grade 3, 157 (0.4%) grade 4, with all features of rHOA being more common in males than females. rHOA grades ≥2 were associated with all three clinical outcomes in both unadjusted and adjusted models, a clear dose-response relationship was seen with each increase in grade showing a large rise in OR/HRs (Table 1). Grade 4 rHOA was strongly predictive of THR (HR 57.70 [95%CI 38.08-87.44]).

Conclusion

We successfully applied a novel semi-automated classifier to over 40,000 individuals from UKB. The validity of our classifier was supported by the strong and progressive relationships observed between rHOA and hip pain, and HES OA and risk of THR. We conclude that hip DXAs provide a promising means of defining rHOA, with potential screening applications in the clinic.

Disclosure

B.G. Faber: None. R. Ebsim: None. F.R. Saunders: None. M. Frysz: None. C. Lindner: None. J.S. Gregory: None. R.M. Aspden: None. N.C. Harvey: None. G. Davey Smith: None. T. Cootes: None. J.H. Tobias: None.

Analysis of hydration and subchondral bone density on the viscoelastic properties of bovine articular cartilage

Background

Articular cartilage is known to be a viscoelastic material, however little research has explored the impact of cartilage water content and bone density on its viscoelasticity. This study aimed to isolate subchondral bone density and hydration of articular cartilage and analyse their effects on the viscoelastic properties of articular cartilage.

Methods

Dynamic mechanical analysis was used to test samples at frequencies of 1, 8, 12, 29, 49, 71, and 88 Hz. Synthetic bone material with densities of 663.7 kg/m³ and 156.8 kg/m³ were used to mimic the bone mineral density (BMD). Dehydration occurred in a stepwise manner at relative humidity (RH) levels of 100%, 30%, and 1%. These relative humidity levels led to water contents of approximately 76%, 8.5%, and ≈ 0% by mass, respectively.

Results

Samples from eight bovine femoral heads were tested under a sinusoidal load. Storage stiffness was lower on the lower substrate density. Storage stiffness, though, increased as cartilage samples were dehydrated from a water content of 76% to 8.5%; decreasing again as the water content was further reduced. Loss stiffness was lower on a lower density substrate and decreased as the water content decreased.

Conclusions

In conclusions, a decrease in hydration decreases the loss stiffness, but a non-linear relationship between hydration and storage stiffness may exist. Additionally, higher BMD values led to greater storage and loss stiffnesses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05169-0.

Osteophyte size and location on hip DXA scans are associated with hip pain: Findings from a cross sectional study in UK Biobank

OBJECTIVE

It remains unclear how the different features of radiographic hip osteoarthritis (rHOA) contribute to hip pain. We examined the relationship between rHOA, including its individual components, and hip pain using a novel dual-energy x-ray absorptiometry (DXA)-based method.

METHODS

Hip DXAs were obtained from UK Biobank. A novel automated method obtained minimum joint space width (mJSW) from points placed around the femoral head and acetabulum. Osteophyte areas at the lateral acetabulum, superior and inferior femoral head were derived manually. Semi-quantitative measures of osteophytes and joint space narrowing (JSN) were combined to define rHOA. Logistic regression was used to examine the relationships between these variables and hip pain, obtained via questionnaires.

RESULTS

6807 hip DXAs were examined. rHOA was present in 353 (5.2%) individuals and was associated with hip pain [OR 2.42 (1.78-3.29)] and hospital diagnosed OA [6.01 (2.98-12.16)]. Total osteophyte area but not mJSW was associated with hip pain in mutually adjusted models [1.31 (1.23-1.39), 0.95 (0.87-1.04) respectively]. On the other hand, JSN as a categorical variable showed weak associations between grade ≥ 1 and grade ≥ 2 JSN with hip pain [1.30 (1.06-1.60), 1.80 (1.34-2.42) respectively]. Acetabular, superior and inferior femoral osteophyte areas were all independently associated with hip pain [1.13 (1.06-1.20), 1.13 (1.05-1.24), 1.10 (1.03-1.17) respectively].

CONCLUSION

In this cohort, the relationship between rHOA and prevalent hip pain was explained by 2-dimensional osteophyte area, but not by the apparent mJSW. Osteophytes at different locations showed important, potentially independent, associations with hip pain, possibly reflecting the contribution of distinct biomechanical pathways.

Cam morphology but neither acetabular dysplasia nor pincer morphology is associated with osteophytosis throughout the hip: findings from a cross-sectional study in UK Biobank

OBJECTIVES

To examine whether acetabular dysplasia (AD), cam and/or pincer morphology are associated with radiographic hip osteoarthritis (rHOA) and hip pain in UK Biobank (UKB) and, if so, what distribution of osteophytes is observed.

DESIGN

Participants from UKB with a left hip dual-energy X-ray absorptiometry (DXA) scan had alpha angle (AA), lateral centre-edge angle (LCEA) and joint space narrowing (JSN) derived automatically. Cam and pincer morphology, and AD were defined using AA and LCEA. Osteophytes were measured manually and rHOA grades were calculated from JSN and osteophyte measures. Logistic regression was used to examine the relationships between these hip morphologies and rHOA, osteophytes, JSN, and hip pain.

RESULTS

6,807 individuals were selected (mean age: 62.7; 3382/3425 males/females). Cam morphology was more prevalent in males than females (15.4% and 1.8% respectively). In males, cam morphology was associated with rHOA [OR 3.20 (95% CI 2.41-4.25)], JSN [1.53 (1.24-1.88)], and acetabular [1.87 (1.48-2.36)], superior [1.94 (1.45-2.57)] and inferior [4.75 (3.44-6.57)] femoral osteophytes, and hip pain [1.48 (1.05-2.09)]. Broadly similar associations were seen in females, but with weaker statistical evidence. Neither pincer morphology nor AD showed any associations with rHOA or hip pain.

CONCLUSIONS

Cam morphology was predominantly seen in males in whom it was associated with rHOA and hip pain. In males and females, cam morphology was associated with inferior femoral head osteophytes more strongly than those at the superior femoral head and acetabulum. Further studies are justified to characterise the biomechanical disturbances associated with cam morphology, underlying the observed osteophyte distribution.

The influence of adult hip shape genetic variants on adolescent hip shape: Findings from a population-based DXA study

OBJECTIVE

Hip shape is a well-recognized risk factor for hip osteoarthritis (OA) and hip fracture. We aimed to investigate whether the genetic variants known to be associated with adult hip shape were also associated with adolescent hip shape.

METHODS

Hip DXA scans, obtained in offspring from the Avon Longitudinal Study of Parents and Children (ALSPAC) at two time points (mean ages 13.8 and 17.8 years), were used to quantify hip morphology using a 53-point Statistical Shape Model (SSM). Principal component analysis was used to generate hip shape modes (HSMs). Genetic variants which had previously shown genome-wide significant association with specific HSMs in adults were tested for association with the same HSMs in adolescents (at each timepoint separately) using SNPTEST v2.

RESULTS

Complete genotypic and phenotypic data were available for 3550 and 3175 individuals at 14 and 18 years, respectively. The strongest evidence for association with adolescent hip shape was for a variant located near SOX9 (rs2158915) with consistent effects across both time points for HSM1 (age 14: beta -0.05, p = 9.9 × 10-8; age 18: -0.05, p = 3.3 × 10-6) and HSM5 (age 14: beta -0.07, p = 1.6 × 10-4; age 18: -0.1, p = 2.7 × 10-6). There was also strong evidence of association between rs10743612 (near PTHLH) and HSM1 (age 14: 0.05, p = 1.1 × 10-5; age 18: 0.04, p = 0.003) and between rs6537291 (near HHIP) and HSM2 (age 14: -0.06, p = 0.001; age 18: -0.07, p = 0.001) across both time points. The genes with the strongest associations with hip shape in adolescents, (SOX9, PTHLH and HHIP) are known to be involved in endochondral bone formation. HSM1 indicates narrower aspect ratio of the upper femur, whereas both HSM2 and HSM5 reflect variation in the femoral head size and femoral neck width, features previously found to be related to the risk of OA in later life. The SOX9 locus has previously been found to associate with increased risk of hip fracture.

CONCLUSION

In conclusion, variants implicated in endochondral bone formation appear to consistently influence hip shape between adolescence and adulthood, including those aspects related to risk of hip OA and/or fracture in later life.

Motor development in infancy and spine shape in early old age: Findings from a British birth cohort study

Spine shape changes dramatically in early life, influenced by attainment of developmental milestones such as independent walking. Whether these associations persist across life is unknown. Therefore, we investigated associations between developmental milestones and spine shape, as determined using statistical shape models (SSMs) of lumbar spine from dual-energy X-ray absorptiometry scans in 1327 individuals (688 female) at 60 to 64 years in the MRC National Survey of Health and Development. Lumbar lordosis angle (L4 inferior endplate to T12 superior endplate) was measured using the two-line Cobb method. In analyses adjusted for sex, height, lean and fat mass, socioeconomic position, and birthweight, later walking age was associated with greater lordosis described by SSM1 (regression coefficient, 0.023; 95% CI, 0.000-0.047; P = .05) and direct angle measurement. Modest associations between walking age and less variation in anterior-posterior vertebral size caudally (SSM6) were also observed (0.021; 95% CI, -0.002 to 0.044; P = .07). Sex interactions showed that later walking was associated with larger relative vertebral anterior-posterior dimensions in men (SSM3; -0.043; 95% CI, -0.075 to 0.01; P = .01) but not women (0.018; 95% CI, -0.0007 to 0.043; P = .17). Similar associations were observed between age at independent standing and SSMs but there was little evidence of association between sitting age and spine shape. Unadjusted associations between walking age and SSMs 1 and 6 remained similar after adjustment for potential confounders and mediators. This suggests that these associations may be explained by altered mechanical loading of the spine during childhood growth, although other factors could contribute. Early life motor development, particularly walking, may have a lasting effect on the features of spine morphology with clinical significance.

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.

Subregional statistical shape modelling identifies lesser trochanter size as a possible risk factor for radiographic hip osteoarthritis, a cross-sectional analysis from the Osteoporotic Fractures in Men Study

OBJECTIVE

Statistical shape modelling (SSM) of hip dual-energy X-ray absorptiometry (DXA) scans has identified relationships between hip shape and radiographic hip OA (rHOA). We aimed to further elucidate shape characteristics related to rHOA by focusing on subregions identified from whole-hip shape models.

METHOD

SSM was applied to hip DXAs obtained in the Osteoporotic Fractures in Men Study. Whole-hip shape modes (HSMs) associated with rHOA were combined to form a composite at-risk-shape. Subsequently, subregional HSMs (cam-type and lesser trochanter modes) were built, and associations with rHOA were examined by logistic regression. Subregional HSMs were further characterised, by examining associations with 3D-HSMs derived from concurrent hip CT scans.

RESULTS

4,098 participants were identified with hip DXAs and radiographs. Composite shapes from whole-hip HSMs revealed that lesser trochanter size and cam-type femoral head are related to rHOA. From sub-regional models, lesser trochanter mode (LTM)1 [OR 0.74; 95%CI 0.63.0.87] and cam-type mode (CTM)3 [OR 1.27; 1.13.1.42] were associated with rHOA, associations being similar to those for whole hip HSMs. 515 MrOS participants had hip DXAs and 3D-HSMs derived from hip CT scans. LTM1 was associated with 3D-HSMs that also represented a larger lesser trochanter [3D-HSM7 (beta (β)-0.23;-0.33,-0.14) and 3D-HSM9 (β0.36; 0.27.0.45)], and CTM3 with 3D-HSMs describing cam morphology [3D-HSM3 (β-0.16;-0.25,-0.07) and 3D-HSM6 (β 0.19; 0.10.0.28)].

CONCLUSION

Subregional SSM of hip DXA scans suggested larger lesser trochanter and cam morphology underlie associations between overall hip shape and rHOA. 3D hip modelling suggests our subregional SSMs represent true anatomical variations in hip shape, warranting further investigation.

Is intrinsic lumbar spine shape associated with lumbar disc degeneration? An exploratory study

Background

Lumbar disc degeneration (LDD) is a condition associated with recurrent low back pain (LBP). Knowledge regarding effective management is limited. As a step towards the identification of risk, prognostic or potentially modifiable factors in LDD patients, the aim of this study was to explore the hypothesis that intrinsic lumbar spine shape is associated with LDD and clinical outcomes in symptomatic adults.

Methods

3 T MRI was used to acquire T2-weighted sagittal images (L1-S1) from 70 healthy controls and LDD patients (mean age 49 years, SD 11, range 31–71 years). Statistical Shape Modelling (SSM) was used to describe lumbar spine shape. SSM identified variations in lumbar shape as ‘modes’ of variation and quantified deviation from the mean. Intrinsic shape differences were determined between LDD groups using analysis of variance with post-hoc comparisons. The relationship between intrinsic shape and self-reported function, mental health and quality of life were also examined.

Results

The first 7 modes of variation explained 91% of variance in lumbar shape. Higher LDD sum scores correlated with a larger lumbar lordosis (Mode 1 (55% variance), P = 0.02), even lumbar curve distribution (Mode 2 (12% variance), P = 0.05), larger anterior-posterior (A-P) vertebral diameter (Mode 3 (10% variance), P = 0.007) and smaller L4-S1 disc spaces (Mode 7 (2% variance), P ≤ 0.001). In the presence of recurrent LBP, LDD was associated with a larger A-P vertebral diameter (Mode 3) and a more even lumbar curvature with smaller L5/S1 disc spaces (Mode 4), which was significantly associated with patient quality of life (P = 0.002–0.04, r_p = 0.43–0.61)).

Conclusions

This exploratory study provides new evidence that intrinsic shape phenotypes are associated with LDD and quality of life in patients. Longitudinal studies are required to establish the potential role of these risk or prognostic shape phenotypes.

The effect of pubertal timing, as reflected by height tempo, on proximal femur shape: Findings from a population-based study in adolescents

OBJECTIVE

To examine the relationship between pubertal timing (using measures of height tempo) and proximal femur shape in a large adolescent cohort.

METHODS

Hip DXA scans were obtained in offspring from the Avon Longitudinal Study of Parents and Children. To quantify hip morphology, the images were analyzed using Shape software based on a 53-point statistical shape model and independent modes of variation (hip shape mode (HSM) scores) for each image were generated. Height tempo (which corresponds to age at peak height velocity (aPHV)) was estimated from serial height measurements collected between age 5-20 years. Multivariable linear regression was used to examine cross-sectional associations between height tempo and the top ten HSMs at age 14 and 18, adjusting for sex and fat mass index (FMI).

RESULTS

Complete outcome and covariate data were available from 3827 and 3507 participants at age 14 and 18 years, respectively. Mean aPHV was 13.5 and 11.8 years for males and females, respectively. At age 14, height tempo was associated with a majority of modes, except for HSM4 and there was strong evidence of interaction by sex. In males, all modes showed evidence of an association with tempo, independent of FMI, with the strongest observed for HSM8 (adjusted β 0.38 (0.33, 0.43) p = 4.1 × 10^-50). Compared with males, the associations were generally weaker in females, with the strongest effect observed for HSM8 (adjusted β 0.10 (0.05, 0.14) p = 1.6 × 10^-5). The overall effect of later pubertal timing on proximal femur shape in males was a narrower femoral neck and larger superolateral head, whereas in females these changes were hard to discern. When assessed at age 18, there was little relationship between tempo and proximal femur shape in either sex.

CONCLUSION

Our results indicate that significant changes in hip shape occur during puberty, including aspects of shape which may be related to future risk of hip OA and/or fracture. However, puberty timing per se does not appear to exert long lasting effects on proximal femur shape.

Associations between prenatal indicators of mechanical loading and proximal femur shape: findings from a population-based study in ALSPAC offspring

OBJECTIVES

Hip development is influenced by mechanical loading, but associations between prenatal loading and hip shape in later life remain unexplored.

METHODS

We examined associations between prenatal loading indicators (gestation length, oligohydramnios (OH) and breech) obtained from obstetric records and hip shape modes (HSMs) generated using dual-energy X-ray absorptiometry images taken at age 14- and 18-years in participants from the UK Avon Longitudinal Study of Parents and Children (ALSPAC). These associations were examined in 2453 (30 OH, 105 breech) and 2330 (27 OH, 95 breech) participants with complete data at age 14- and 18-years respectively using confounder-adjusted models.

RESULTS

At 14 years HSM2 was 0.59SD lower in OH males, and HSM5 (-0.31SD) and HSM9 (-0.32SD) were lower in OH in both sexes. At 18 years HSM1 (-0.44SD) and HSM2 (-0.71SD) were lower and HSM6 (0.61SD) and HSM8 (1.06SD) were higher in OH males, whilst HSM5 was lower in OH in both sexes. OH appeared to be associated with a wider femoral neck and head, and larger lesser/greater trochanters. Only weak associations were observed between gestation length/breech and HSMs.

CONCLUSIONS

These results suggest that prenatal skeletal loading, in particular oligohydramnios, may influence adolescent joint shape with associations generally stronger in males.

Describing the application of statistical shape modelling to DXA images to quantify the shape of the proximal femur at ages 14 and 18 years in the Avon Longitudinal Study of Parents and Children

Bones are complex objects with considerable variation in the shape and structure often attributed to anatomical, environmental or genetic differences. In addition, bone shape has been of interest in relation to its associations with disease processes. Hip shape is an important determinant of hip osteoarthritis and osteoporotic hip fracture; however, its quantification is difficult. While previous studies largely focused on individual geometrical indices of hip geometry such as neck-shaft angle or femoral neck width, statistical shape modelling offers the means to quantify the entire contour of the proximal femur, including lesser trochanter and acetabular eyebrow. We describe the derivation of independent modes of variation (hip shape mode scores) to characterise variation in hip shape from dual-energy X-ray absorptiometry (DXA) images in the Avon Longitudinal Study of Parents and Children (ALSPAC) offspring, using statistical shape modelling. ALSPAC is a rich source of phenotypic and genotypic data which provides a unique opportunity to investigate the environmental and genetic influences on hip shape in adolescence, as well as comparison with adult hip shape.

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.

Using statistical shape modelling of DXA images to quantify the shape of the proximal femur at ages 14 and 18 years in the Avon Longitudinal Study of Parents and Children

Hip shape is an important determinant of hip osteoarthritis and osteoporotic hip fracture; however, little is known about its development in childhood and adolescence. While previous studies largely focused on individual geometrical indices of hip geometry such as neck-shaft angle or femoral neck width, statistical shape modelling offers the means to quantify the entire contour of the proximal femur, including lesser trochanter and acetabular eyebrow. We describe the derivation of independent modes of variation (hip shape mode scores) to characterise variation in hip shape from dual-energy X-ray absorptiometry (DXA) images in the Avon Longitudinal Study of Parents and Children (ALSPAC) offspring, using statistical shape modelling. ALSPAC is a rich source of phenotypic and genotypic data which provides a unique opportunity to investigate the environmental and genetic influences on hip shape in adolescence, as well as comparison with adult hip shape.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Osteoarthritis as an organ disease: from the cradle to the grave

Considered for decades as a cartilage disease, recent studies of osteoarthritis (OA) take us back to the concepts discussed at the naming of the disorder as "bone-joint-inflammation". By describing the joint as an organ, can OA be called an organ disease - similar to heart disease? Is there a systemic (which system?) involvement? Would this help with diagnosis or therapy? Hyperplasia of the joint tissues is one of the most notable early features of the disease: articular cartilage thickens, chondrocytes proliferate and increase matrix biosynthesis, but not its incorporation; the subchondral bone densifies but is hypomineralised and there is an increase in bone marrow fat content. Associations between OA and hypertension, hypercholesterolaemia and blood glucose suggest systemic and metabolic components are involved. The source of pain is still unknown but here is evidence for peripheral and central sensitisation. Joint deformity is difficult to quantify, but statistical shape modelling provides a tool to use as an imaging biomarker. A genome-wide association study meta-analysis has identified novel genes associated with hip shape with many genes related to tissue growth and development. There are associations between hip shapes and age of first walking as well as with obesity through adulthood. These life-course events and a recapitulation in old age of developmental processes suggest that the cradle may affect our path to the grave. These observations suggest that tissue regeneration approaches, treating only the cartilage in OA joints, may only be of limited benefit.

Raman Microscopy and Bone

Raman microscopy is a nondestructive technique requiring minimal sample preparation that can be used to measure the chemical properties of the mineral and collagen parts of bone simultaneously. Modern Raman instruments contain the necessary components and software to acquire the standard information required in most bone studies. The spatial resolution of the technique is about a micron. As it is nondestructive and small samples can be used, it forms a useful part of a bone characterization toolbox.

Mechanical Properties of Bone Ex Vivo

The primary functions of bone are to provide support and protection-mechanical functions. The aim of this chapter is to set out some of the methods that can be used to measure these properties in cortical and cancelleous bone from large (e.g., human or bovine) and small (e.g., mouse) animals. The difference between the properties of the sample (extrinsic properties) and the properties of the material (intrinsic properties) is introduced and techniques for measuring them suggested. The addition of other tests to give a complete characterization of a bone sample is presented.

Variation in lifting kinematics related to individual intrinsic lumbar curvature: an investigation in healthy adults

Objective

Lifting postures are frequently implicated in back pain. We previously related responses to a static load with intrinsic spine shape, and here we investigate the role of lumbar spine shape in lifting kinematics.

Methods

Thirty healthy adults (18-65 years) performed freestyle, stoop and squat lifts with a weighted box (6-15 kg, self-selected) while being recorded by Vicon motion capture. Internal spine shape was characterised using statistical shape modelling (SSM) from standing mid-sagittal MRIs. Associations were investigated between spine shapes quantified by SSM and peak flexion angles.

Results

Two SSM modes described variations in overall lumbar curvature (mode 1 (M1), 55% variance) and the evenness of curvature distribution (mode 2 (M2), 12% variance). M1 was associated with greater peak pelvis (r=0.38, p=0.04) and smaller knee flexion (r=-0.40, p=0.03) angles; individuals with greater curviness preferred to lift with a stooped lifting posture. This was confirmed by analysis of those individuals with very curvy or very straight spines (|M1|>1 SD). There were no associations between peak flexion angles and mode scores in stoop or squat trials (p>0.05). Peak flexion angles were positively correlated between freestyle and squat trials but not between freestyle and stoop or squat and stoop, indicating that individuals adjusted knee flexion while maintaining their preferred range of lumbar flexion and that 'squatters' adapted better to different techniques than 'stoopers'.

Conclusion

Spinal curvature affects preferred lifting styles, and individuals with curvier spines adapt more easily to different lifting techniques. Lifting tasks may need to be tailored to an individual's lumbar spine shape.

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.

Sex Differences in the Associations between L-Arginine Pathway Metabolites, Skeletal Muscle Mass and Function, and their Responses to Resistance Exercise, in Old Age

OBJECTIVES

The current study was designed to explore the associations between L-arginine metabolites and muscle mass and function in old age, which are largely unknown.

DESIGN

The study used a randomised, double-blind, placebo-controlled design.

SETTING

The study was carried out in a laboratory setting.

PARTICIPANTS

50 healthy older adults [median age 70 years (IQR 67-73); 27 males].

INTERVENTION

Participants undertook an 18-week resistance exercise program, and a nutritional intervention (fish oil vs. placebo).

MEASUREMENTS

Serum homoarginine, ornithine, citrulline, asymmetric dimethylarginine (ADMA), NG-monomethyl-L-arginine (L-NMMA), and symmetric dimethylarginine (SDMA), maximal voluntary contraction (MVC) and isokinetic torque of the knee extensors at 30° s-1 (MIT), muscle cross sectional area (MCSA) and quality (MQ) were measured at baseline and after the intervention.

RESULTS

No significant exercise-induced changes were observed in metabolite concentrations. There were significant sex differences in the associations between metabolites and muscle parameters. After adjusting for age, glomerular filtration rate and fish oil intervention, citrulline (P=0.002) and ornithine (P=0.022) were negatively associated with MCSA at baseline in males but not females. However, baseline citrulline was negatively correlated with exercise-induced changes in MVC (P=0.043) and MQ (P=0.026) amongst females. Furthermore, amongst males, baseline homoarginine was positively associated with exercise-induced changes in MVC (P=0.026), ADMA was negatively associated with changes in MIT (P=0.026), L-NMMA (p=0.048) and ornithine (P<0.001) were both positively associated with changes in MCSA, and ornithine was negatively associated with changes in MQ (P=0.039).

CONCLUSION

Therefore, barring citrulline, there are significant sex differences in the associations between L-arginine metabolites and muscle mass and function in healthy older adults. These metabolites might enhance sarcopenia risk stratification, and the success of exercise programs, in old age.

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.

Age-related changes in the effects of strength training on lower leg muscles in healthy individuals measured using MRI

Background

We previously measured the rate of regaining muscle strength during rehabilitation of lower leg muscles in patients following lower leg casting. Our primary aim in this study was to measure the rate of gain of strength in healthy individuals undergoing a similar training regime. Our secondary aim was to test the ability of MRI to provide a biomarker for muscle function.

Methods

Men and women were recruited in three age groups: 20–30, 50–65 and over 70 years. Their response to resistance training of the right lower leg twice a week for 8 weeks was monitored using a dynamometer and MRI of tibialis anterior, soleus and gastrocnemius muscles at 2 weekly intervals to measure muscle size (anatomical cross-sectional area (ACSA)) and quality (T₂ relaxation). Forty-four volunteers completed the study.

Results

Baseline strength declined with age. Training had no effect in middle-aged females or in elderly men in dorsiflexion. Other groups significantly increased both plantarflexion and dorsiflexion strength at rates up to 5.5 N m week^-1 in young females in plantarflexion and 1.25 N m week^-1 in young males in dorsiflexion. No changes were observed in ACSA or T₂ in any age group in any muscle.

Conclusion

Exercise training improves muscle strength in males at all ages except the elderly in dorsiflexion. Responses in females were less clear with variation across age and muscle groups. These results were not reflected in simple MRI measures that do not, therefore, provide a good biomarker for muscle atrophy or the efficacy of rehabilitation.

Sex differences in the effect of fish-oil supplementation on the adaptive response to resistance exercise training in older people: a randomized controlled trial

BACKGROUND

Resistance exercise increases muscle mass and function in older adults, but responses are attenuated compared with younger people. Data suggest that long-chain n-3 polyunsaturated fatty acids (PUFAs) may enhance adaptations to resistance exercise in older women. To our knowledge, this possibility has not been investigated in men.

OBJECTIVE

We sought to determine the effects of long-chain n-3 PUFA supplementation on resistance exercise training-induced increases in muscle mass and function and whether these effects differ between older men and women.

DESIGN

Fifty men and women [men: n = 27, mean ± SD age: 70.6 ± 4.5 y, mean ± SD body mass index (BMI; in kg/m2): 25.6 ± 4.2; women: n = 23, mean ± SD age: 70.7 ± 3.3 y, mean ± SD BMI: 25.3 ± 4.7] were randomly assigned to either long-chain n-3 PUFA (n = 23; 3 g fish oil/d) or placebo (n = 27; 3 g safflower oil/d) and participated in lower-limb resistance exercise training twice weekly for 18 wk. Muscle size, strength, and quality (strength per unit muscle area), functional abilities, and circulating metabolic and inflammatory markers were measured before and after the intervention.

RESULTS

Maximal isometric torque increased after exercise training to a greater (P < 0.05) extent in the long-chain n-3 PUFA group than in the placebo group in women, with no differences (P > 0.05) between groups in men. In both sexes, the effect of exercise training on maximal isokinetic torque at 30, 90, and 240° s-1, 4-m walk time, chair-rise time, muscle anatomic cross-sectional area, and muscle fat did not differ (P > 0.05) between groups. There was a greater (P < 0.05) increase in muscle quality in women after exercise training in the long-chain n-3 PUFA group than in the placebo group, with no such differences in men (P > 0.05). Long-chain n-3 PUFAs resulted in a greater decrease (P < 0.05) than the placebo in plasma triglyceride concentrations in both sexes, with no differences (P > 0.05) in glucose, insulin, or inflammatory markers.

CONCLUSION

Long-chain n-3 PUFA supplementation augments increases in muscle function and quality in older women but not in older men after resistance exercise training. This trial was registered at clinicaltrials.gov as NCT02843009.

Fibre Reinforcement and Mechanical Stability in Articular Cartilage

The gel phase of articular cartilage is reinforced by collagen fibrils. These fibrils have low flexural and torsional stiffness, but are able to provide reinforcement if deformation of the tissue increases their tensile stress. An estimate suggests that the lengths of collagen fibrils in articular cartilage are at least of the same order as their critical length so that tensile stress in the tissue will increase the stress in the fibrils rather than simply pull them out of the gel. In the surface zone the collagen fibrils are oriented so that the efficiency of reinforcement, η, is about 0.6 tangential to the surface; tension in the fibrils is thus able to withstand swelling pressure within the tissue whose condition for stability resembles that of a pressure vessel. Swelling pressure allows the tissue to support applied pressure. An intermediate zone has a roughly isotropic η value of about 0.2, while in the deep zone collagen fibrils appear to tie the cartilage to the subchondral bone; in this deep zone η has a value of about 0.6 perpendicular to the surface direction. There is also some preferred orientation of collagen fibrils in the plane of the articular surface within the surface zone; in patellar cartilage the preferred orientations can be related to the direction of stress which could be generated by movement of the joint.

Viscoelastic properties of human and bovine articular cartilage: a comparison of frequency-dependent trends

Background

The purpose of this study was to compare the frequency-dependent viscoelastic properties of human and bovine cartilage.

Methods

Full-depth cartilage specimens were extracted from bovine and human femoral heads. Using dynamic mechanical analysis, the viscoelastic properties of eight bovine and six human specimens were measured over the frequency range 1 Hz to 88 Hz. Significant differences between bovine and human cartilage viscoelastic properties were assessed using a Mann–Whitney test (p < 0.05).

Results

Throughout the range of frequencies tested and for both species, the storage modulus was greater than the loss modulus and both were frequency-dependent. The storage and loss moduli of all human and bovine cartilage specimens presented a logarithmic relationship with respect to frequency. The mean human storage modulus ranged from 31.9 MPa to 43.3 MPa, while the mean bovine storage modulus ranged from 54.0 MPa to 80.5 MPa; bovine storage moduli were 1.7 to 1.9 times greater than the human modulus. Similarly, the loss modulus of bovine cartilage was 2.0 to 2.1 times greater than human. The mean human loss modulus ranged from 5.3 MPa to 8.5 MPa while bovine moduli ranged from 10.6 MPa to 18.1 MPa.

Conclusion

Frequency-dependent viscoelastic trends of bovine articular cartilage were consistent with those of human articular cartilage; this includes a similar frequency dependency and high-frequency plateau. Bovine cartilage was, however, ‘stiffer’ than human by a factor of approximately 2. With these provisos, bovine articular cartilage may be a suitable dynamic model for human articular cartilage.

Sex differences in the response to resistance exercise training in older people

Resistance exercise training is known to be effective in increasing muscle mass in older people. Acute measurement of protein metabolism data has indicated that the magnitude of response may differ between sexes. We compared adaptive responses in muscle mass and function to 18 weeks resistance exercise training in a cohort of older (>65 years) men and women. Resistance exercise training improved knee extensor maximal torque, 4 m walk time, time to complete five chair rises, muscle anatomical cross‐sectional area (ACSA) and muscle quality with no effect on muscle fat/water ratio or plasma glucose, insulin, triacylglycerol, IL‐6, and TNF‐α. Differences between sexes were observed for knee extensor maximal torque and muscle quality with greater increases observed in men versus women (P < 0.05). Maximal torque increased by 15.8 ± 10.6% in women and 41.7 ± 25.5% in men, whereas muscle quality increased by 8.8 ± 17.5% in women and by 33.7 ± 25.6% in men. In conclusion, this study has demonstrated a difference in the magnitude of adaptation, of some of the outcome measures employed, in response to 18 weeks of resistance exercise training between men and women. The mechanisms underlying this observation remain to be established.

Genome-wide association of multiple complex traits in outbred mice by ultra-low-coverage sequencing

Two bottlenecks impeding the genetic analysis of complex traits in rodents are access to mapping populations able to deliver gene-level mapping resolution, and the need for population specific genotyping arrays and haplotype reference panels. Here we combine low coverage sequencing (0.15X) with a novel method to impute the ancestral haplotype space in 1,887 commercially available outbred mice. We mapped 156 unique quantitative trait loci for 92 phenotypes at 5% false discovery rate. Gene-level mapping resolution was achieved at about a fifth of loci, implicating Unc13c and Pgc1-alpha at loci for the quality of sleep, Adarb2 for home cage activity, Rtkn2 for intensity of reaction to startle, Bmp2 for wound healing, Il15 and Id2 for several T-cell measures and Prkca for bone mineral content. These findings have implications for diverse areas of mammalian biology and demonstrate how GWAS can be extended via low-coverage sequencing to species with highly recombinant outbred populations.

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.

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.