Cluster Analysis of Finite Element Analysis and Bone Microarchitectural Parameters Identifies Phenotypes with High Fracture Risk

High-resolution peripheral quantitative computed tomography (HRpQCT) is increasingly used for exploring associations between bone microarchitectural and finite element analysis (FEA) parameters and fracture. We hypothesised that combining bone microarchitectural parameters, geometry, BMD and FEA estimates of bone strength from HRpQCT may improve discrimination of fragility fractures. The analysis sample comprised of 359 participants (aged 72-81 years) from the Hertfordshire Cohort Study. Fracture history was determined by self-report and vertebral fracture assessment. Participants underwent HRpQCT scans of the distal radius and DXA scans of the proximal femur and lateral spine. Poisson regression with robust variance estimation was used to derive relative risks for the relationship between individual bone microarchitectural and FEA parameters and previous fracture. Cluster analysis of these parameters was then performed to identify phenotypes associated with fracture prevalence. Receiver operating characteristic analysis suggested that bone microarchitectural parameters improved fracture discrimination compared to aBMD alone, whereas further inclusion of FEA parameters resulted in minimal improvements. Cluster analysis (k-means) identified four clusters. The first had lower Young modulus, cortical thickness, cortical volumetric density and Von Mises stresses compared to the wider sample; fracture rates were only significantly greater among women (relative risk [95%CI] compared to lowest risk cluster: 2.55 [1.28, 5.07], p = 0.008). The second cluster in women had greater trabecular separation, lower trabecular volumetric density and lower trabecular load with an increase in fracture rate compared to lowest risk cluster (1.93 [0.98, 3.78], p = 0.057). These findings may help inform intervention strategies for the prevention and management of osteoporosis.